Structural and elastic properties of DNA and chromatin

We investigate a chain consisting of two coupled worm-like chains withconstant distance between the strands. The effects due todouble-strandedness of the chain are studied. In a previous analyticalstudy of this system an intrinsic twist-stretch coupling and atendency of kinking is predicted. Even though a local twist structureis observed the predicted features are not recovered. A new model for DNA at the base-pair level is presented. Thebase-pairs are treated as flat rigid ellipsoids and thesugar-phosphate backbones are represented as stiff harmonic springs.The base-pair stacking interaction is modeled by a variant of theGay-Berne potential. It is shown by systematic coarse-graininghow the elastic constants of a worm-like chain are related to thelocal fluctuations of the base-pair step parameters. Even though a lotof microscopic details of the base-pair geometry is neglected themodel can be optimized to obtain a B-DNA conformation as ground stateand reasonable elastic properties. Moreover the model allows tosimulate much larger length scales than it is possible with atomisticsimulations due to the simplification of the force-field and inparticular due to the possibility of non-local Monte-Carlo moves. Asa first application the behavior under stretching is investigated. Inagreement with micromanipulation experiments on single DNA moleculesone observes a force-plateau in the force-extension curvescorresponding to an overstretching transition from B-DNA to aso-called S-DNA state. The model suggests a structure for S-DNA withhighly inclined base-pairs in order to enable at least partialbase-pair stacking. Finally a simple model for chromatin is introduced to study itsstructural and elastic properties. The underlying geometry of themodeled fiber is based on a crossed-linker model. The chromatosomesare treated as disk-like objects. Excluded volume and short rangenucleosomal interaction are taken into account by a variant of theGay-Berne potential. It is found that the bending rigidity and thestretching modulus of the fiber increase with more compact fibers. Fora reasonable parameterization of the fiber for physiologicalconditions and sufficiently high attraction between the nucleosomes aforce-extension curve is found similar to stretching experiments onsingle chromatin fibers. For very small stretching forces a kinkedfiber forming a loop is observed. If larger forces are applied theloop formation is stretched out and a decondensation of the fibertakes place.

Modeling, design and experimental characterization of Micro-Electro-Mechanical-Systems for gas chromatographic applications

Design parameters, process flows, electro-thermal-fluidic simulations and experimental characterizations of Micro-Electro-Mechanical-Systems (MEMS) suited for gas-chromatographic (GC) applications are presented and thoroughly described in this thesis, whose topic belongs to the research activities the Institute for Microelectronics and Microsystems (IMM)-Bologna is involved since several years, i.e. the development of micro-systems for chemical analysis, based on silicon micro-machining techniques and able to perform analysis of complex gaseous mixtures, especially in the field of environmental monitoring. In this regard, attention has been focused on the development of micro-fabricated devices to be employed in a portable mini-GC system for the analysis of aromatic Volatile Organic Compounds (VOC) like Benzene, Toluene, Ethyl-benzene and Xylene (BTEX), i.e. chemical compounds which can significantly affect environment and human health because of their demonstrated carcinogenicity (benzene) or toxicity (toluene, xylene) even at parts per billion (ppb) concentrations. The most significant results achieved through the laboratory functional characterization of the mini-GC system have been reported, together with in-field analysis results carried out in a station of the Bologna air monitoring network and compared with those provided by a commercial GC system. The development of more advanced prototypes of micro-fabricated devices specifically suited for FAST-GC have been also presented (silicon capillary columns, Ultra-Low-Power (ULP) Metal OXide (MOX) sensor, Thermal Conductivity Detector (TCD)), together with the technological processes for their fabrication. The experimentally demonstrated very high sensitivity of ULP-MOX sensors to VOCs, coupled with the extremely low power consumption, makes the developed ULP-MOX sensor the most performing metal oxide sensor reported up to now in literature, while preliminary test results proved that the developed silicon capillary columns are capable of performances comparable to those of the best fused silica capillary columns. Finally, the development and the validation of a coupled electro-thermal Finite Element Model suited for both steady-state and transient analysis of the micro-devices has been described, and subsequently implemented with a fluidic part to investigate devices behaviour in presence of a gas flowing with certain volumetric flow rates.

Synthese und Photochemie von Dendrimeren mit stilbenoiden Chromophoren

Dendrimere spielen als strukturtreue Nanopartikel eine herausragende Rolle. Ziel dieser Arbeit war, Dendrimere mit einer hohen Dichte an photoaktiven Chromophoren herzu-stellen und zu untersuchen. Dazu wurden die terminalen Aminogruppen von Poly(propylenimin)dendrimeren 1. und 2. Generation, Astramol DAB-Am-4R und DAB-Am-8R, mit Stilbenen und Styrylstilbenen als Chromophor verknüpft. Mittels Wittig-Horner- und Heck-Reaktion wurden (E)-Stilbene aufgebaut, die auf der einen Seite drei Propoxygruppen zur Verbesserung der Löslichkeit und auf der anderen Seite eine passende Funktionalität zur Verknüpfung mit dem dendritischen Core tragen. Als Verknüpfungsmethoden wurden die Verknüpfung als Amid (PSDA), Schiffsche Base (PSDS) und Harnstoff (PSDH) getestet. Die Schiffschen Basen wurden außerdem zur Erhöhung der Hydrolysestabilität zum sekundären Amin reduziert (PSDR und PQDR). Durch die Verknüpfung mit dem Core werden die stilbenoiden Chromophore sehr stark photoaktiviert. Das beruht auf einem Singulett-Energietransfer (Förster-Mechanismus) von Chromophor zu Chromophor. Dieser Prozeß konkurriert zu den Deaktivierungsprozessen, verlängert die mittlere S1-Lebensdauer und erhöht somit die Chancen der Photochemie. Der Styrylstilben-Chromophor hat darüber hinaus einen erheblichen Teil seiner UV-Absorbtion bereits im Tageslicht und photopolymerisiert daher bereits im Tageslicht. Vor allem bei den Dendrimeren 2. Generation stellte sich die Frage nach der vollständigen, d.h. achtfachen Umsetzung; das Core sollte als Knäuel vorliegen, die Arme zum Teil nach innen gefaltet und somit dem Reaktand nur bedingt zugänglich. Auch dort konnten unter optimierten Reaktionsbedingungen alle Aminogruppen umgesetzt werden. Die vollständige Umsetzung der Dendrimere wurde mittels NMR und massenspektroskopischen Methoden untersucht. Bei den Absorptionsspektren der Dendrimere 1. Generation ändert sich die Lage der Maxima je nach Art der Verknüpfung der Chromophore mit dem Core. Die Verlängerung des Chromophors um eine Styryleinheit bedingt eine beträchtliche Rotverschiebung. Die Lage der Emissionsmaxima differiert stärker als die Lage der Absorptionsmaxima. Den geringsten Stokes-Shift weist der Harnstoff auf, dann folgt das sekundäre Amin, dann die Schiffsche Base. Dies weist auf unterschiedlich relaxierte S1-Geometrien hin. Die Verbindungen PSDS1, PSDR1 und PSDH1 aus 3,4,5-Tripropoxystilbeneinheit und Astramol-Core 1. Generation DAB-Am-4 wurden in einer Konzentration von 10-5 mol/L belichtet. Der vollständige Photoabbau durch Belichtung in Chloroform mit einer Xenon-Lampe erfolgte ohne jeglichen Filter innerhalb von zehn Minuten (PSDH1), 20 Minuten (PSDR1) und einer Stunde (PSDS1). Allen drei Verbindungen gemeinsam ist das Entstehen eines intermediären neuen Maximums geringer Intensität, das um etwa 100 nm bathochrom verschoben ist. Das Harnstoffsystem weist außerdem ein weiteres intermediäres Maximum bei 614 nm auf. Diese Maxima können (laut früherer Untersuchungen) durch Oxidation entstandenen chinoiden Strukturen zugeordnet werden, deren Lebensdauer (im Sekundenbereich) zu kurz für eine NMR-Charakterisierung ist. PSDR1 wurde außerdem bei höheren Konzentrationen (10-4 und 10-3 mol/L) mit einer Quecksilberlampe mit Pyrex-Filter (lambda > 300 nm) belichtet. Dabei wird, wie erwartet, eine Verbreiterung der NMR-Signale beobachtet. Es bildet sich zunächst cis-Stilben. Außerdem läßt sich bei 4.3 ppm ein Signal beobachten, das von inter- oder intramolekular gebildeten Methinprotonen herrührt. Auch wenn laut MOPAC- und Kraftfeldrechnung die Doppelbindungen ungünstig für eine [pi2s + pi2s]-Cyclodimerisierung zueinander stehen, kann im photochemisch angeregten Zustand eine Geometrie vorherrschen, die die intramolekulare Kopf-Kopf-Cyclobutanbildung ermöglicht. Die massenspektrometrischen Untersuchungen der Belichtungsprodukte (FD, ESI, MALDI-TOF) zeigen als höchste Masse lediglich das Monomer. Allerdings kann dadurch nicht auf eine rein intramolekulare Reaktion geschlossen werden. Die fortschreitende statistische CC-Verknüpfung kann schnell zu vernetzten Nanopartikeln führen, die im Massenspektrometer nicht fliegen. Die NMR-Spektren der mit zunehmender Vernetzung immer schlechter löslich werdenden Teilchen belegen die Oligomerisierung.

Atomistic simulations of cosolvent/water mixtures

In this thesis, atomistic simulations are performed to investigate hydrophobic solvation and hydrophobic interactions in cosolvent/water binary mixtures. Many cosolvent/water binary mixtures exhibit non-ideal behavior caused by aggregation at the molecular scale level although they are stable and homogenous at the macroscopic scale. Force-field based atomistic simulations provide routes to relate atomistic-scale structure and interactions to thermodynamic solution properties. The predicted solution properties are however sensitive to the parameters used to describe the molecular interactions. In this thesis, a force field for tertiary butanol (TBA) and water mixtures is parameterized by making use of the Kirkwood-Buff theory of solution. The new force field is capable of describing the alcohol-alcohol, water-water and alcohol-water clustering in the solution as well as the solution components’ chemical potential derivatives in agreement with experimental data. With the new force field, the preferential solvation and the solvation thermodynamics of a hydrophobic solute in TBA/water mixtures have been studied. First, methane solvation at various TBA/water concentrations is discussed in terms of solvation free energy-, enthalpy- and entropy- changes, which have been compared to experimental data. We observed that the methane solvation free energy varies smoothly with the alcohol/water composition while the solvation enthalpies and entropies vary nonmonotonically. The latter occurs due to structural solvent reorganization contributions which are not present in the free energy change due to exact enthalpy-entropy compensation. It is therefore concluded that the enthalpy and entropy of solvation provide more detailed information on the reorganization of solvent molecules around the inserted solute. Hydrophobic interactions in binary urea/water mixtures are next discussed. This system is particularly relevant in biology (protein folding/unfolding), however, changes in the hydrophobic interaction induced by urea molecules are not well understood. In this thesis, this interaction has been studied by calculating the free energy (potential of mean force), enthalpy and entropy changes as a function of the solute-solute distance in water and in aqueous urea (6.9 M) solution. In chapter 5, the potential of mean force in both solution systems is analyzed in terms of its enthalpic and entropic contributions. In particular, contributions of solvent reorganization in the enthalpy and entropy changes are studied separately to better understand what are the changes in interactions in the system that contribute to the free energy of association of the nonpolar solutes. We observe that in aqueous urea the association between nonpolar solutes remains thermodynamically favorable (i.e., as it is the case in pure water). This observation contrasts a long-standing belief that clusters of nonpolar molecules dissolve completely in the presence of urea molecules. The consequences of our observations for the stability of proteins in concentrated urea solutions are discussed in the chapter 6 of the thesis.

Characterization of geothermal reservoirs’ parameters by inverse problem resolution and geostatistical simulations

BTES (borehole thermal energy storage)systems exchange thermal energy by conduction with the surrounding ground through borehole materials. The spatial variability of the geological properties and the space-time variability of hydrogeological conditions affect the real power rate of heat exchangers and, consequently, the amount of energy extracted from / injected into the ground. For this reason, it is not an easy task to identify the underground thermal properties to use when designing. At the current state of technology, Thermal Response Test (TRT) is the in situ test for the characterization of ground thermal properties with the higher degree of accuracy, but it doesn’t fully solve the problem of characterizing the thermal properties of a shallow geothermal reservoir, simply because it characterizes only the neighborhood of the heat exchanger at hand and only for the test duration. Different analytical and numerical models exist for the characterization of shallow geothermal reservoir, but they are still inadequate and not exhaustive: more sophisticated models must be taken into account and a geostatistical approach is needed to tackle natural variability and estimates uncertainty. The approach adopted for reservoir characterization is the “inverse problem”, typical of oil&gas field analysis. Similarly, we create different realizations of thermal properties by direct sequential simulation and we find the best one fitting real production data (fluid temperature along time). The software used to develop heat production simulation is FEFLOW 5.4 (Finite Element subsurface FLOW system). A geostatistical reservoir model has been set up based on literature thermal properties data and spatial variability hypotheses, and a real TRT has been tested. Then we analyzed and used as well two other codes (SA-Geotherm and FV-Geotherm) which are two implementation of the same numerical model of FEFLOW (Al-Khoury model).

Molecular organization of n-cyanobiphenyl liquid crystals on a molybdenite surface

The alignement and anchoring of liquid crystals on solid surfaces is a key problem for modern device technology that until now has been treated empirically, but that can now be tackled by atomistic computer simulations. Molecular dynamics (MD) simulations were used in this thesis work to study two films of 7 and 8 n-alkyl-4’cyanobiphenyl (7CB and 8CB) liquid crystals , with a thickness of 15 nm, confined between two (001) surfaces of MoS2 (molybdenite). The isotropic and nematic phases of both liquid crystals were simulated, and the resulting structures characterized structurally. A new force field was designed to model the interactions between the liquid crystal (LC) molecules and the surface of molybdenite, while an accurate force field developed previously was used to model the 7CB and 8CB molecules. The results show that the (001) molybdenite surface induces a planar orientation in both the liquid crystals. For the nematic phase of 8CB, one of the two solid/LC interfaces is composed of a first layer of molecules aligned parallel to the surface, followed by a second layer of molecules aligned perpendicular to the surface (also called, homeotropic). The effect of the surface appears to be local in nature as it is confined to the first 15 Angström of the LC film. Conversely, for the nematic phase of 7CB, a planar ordering is established into the LC film. The LC molecules at the interface with the molybdenite appear to align preferentially their alkyl chains toward the solid substrate. The resulting tilt angle of molecules was found to be in good agreement with experimental measurements available in literature. Despite the fact that the MD simulations spanned a time range of more than 100 ns, the nematic phases of both 7CB and 8CB were found not to be completely formed. In order to confirm the findings presented in this thesis, we propose to extend the current study.

Modeling and simulations of some anisotropic soft-matter systems: from biaxial to chiral materials

We have modeled various soft-matter systems with molecular dynamics (MD) simulations. The first topic concerns liquid crystal (LC) biaxial nematic (Nb) phases, that can be possibly used in fast displays. We have investigated the phase organization of biaxial Gay-Berne (GB) mesogens, considering the effects of the orientation, strength and position of a molecular dipole. We have observed that for systems with a central dipole, nematic biaxial phases disappear when increasing dipole strength, while for systems characterized by an offset dipole, the Nb phase is stabilized at very low temperatures. In a second project, in view of their increasing importance as nanomaterials in LC phases, we are developing a DNA coarse-grained (CG) model, in which sugar and phosphate groups are represented with Lennard-Jones spheres, while bases with GB ellipsoids. We have obtained shape, position and orientation parameters for each bead, to best reproduce the atomistic structure of a B-DNA helix. Starting from atomistic simulations results, we have completed a first parametrization of the force field terms, accounting for bonded (bonds, angles and dihedrals) and non-bonded interactions (H-bond and stacking). We are currently validating the model, by investigating stability and melting temperature of various sequences. Finally, in a third project, we aim to explain the mechanism of enantiomeric discrimination due to the presence of a chiral helix of poly(gamma-benzyl L-glutamate) (PBLG), in solution of dimethylformamide (DMF), interacting with chiral or pro-chiral molecules (in our case heptyl butyrate, HEP), after tuning properly an atomistic force field (AMBER). We have observed that DMF and HEP molecules solvate uniformly the PBLG helix, but the pro-chiral solute is on average found closer to the helix with respect to the DMF. The solvent presents a faster isotropic diffusion, twice as HEP, also indicating a stronger interaction of the solute with the helix.

Atomistic simulations of liquid crystals in the bulk and at their interfaces

In this thesis, we have dealt with several problems concerning liquid crystals (LC) phases, either in the bulk or at their interfaces, by the use of atomistic molecular dynamics (MD) simulations. We first focused our attention on simulating and characterizing the bulk smectic phase of 4-n-octyl-4'-cyanobiphenyl (8CB), allowing us to investigate the antiparallel molecular arrangement typical of SmAd smectic phases. A second topic of study was the characterization of the 8CB interface with vacuum by simulating freely suspended thin films, which allowed us to determine the influence of the interface on the orientational and positional order. Then we investigated the LC-water and LC-electrolyte water solution interface. This interface has recently found application in the development of sensors for several compounds, including biological molecules, and here we tried to understand the re-orientation mechanism of LC molecules at the interface which is behind the functioning of these sensors. The characterization of this peculiar interface has incidentally led us to develop a polarizable force field for the pentyl-cyanobiphenyl mesogen, whose process of parametrization and validation is reported here in detail. We have shown that this force field is a significant improvement over its previous, static charge non polarizable version in terms of density, orientational order parameter and translational diffusion.

Molekulardynamik-Simulationen supramolekularer Komplexe unter Einfluss einer externen Kraft

Supramolekulare Komplexe werden durch nichtkovalente Bindungen stabilisiert. Legt man eine externe Kraft an einen solchen Komplex an, ist es möglich, diese Bindungen zu öffnen. Anhand der dafür benötigten Kraft läßt sich die Stabilität des Komplexes bestimmen. Im Rahmen dieser Arbeit wurden zwei supramolekulare Komplexe, die unterschiedliche Arten von nichtkovalenten Bindungen enthalten, mit Hilfe von Molekulardynamik (MD) Simulationen untersucht. In beiden Fällen wurden die relevanten Bindungsstrukturen und deren Stabilität ermittelt.rnZum einen wurden zwei synthetische Calix[4]aren-Catenan-Dimersysteme betrachtet, in denen die beiden Monomere über Wasserstoffbrückenbindungen aneinander gebunden sind. Die Besonderheit dieser Komplexe ist, dass die Monomere aufgrund von verschlauften Alkylketten (Catenan-Struktur) nicht vollständig voneinander getrennt werden können. In Abhängigkeit der Länge derrnAlkylketten findet man für die beiden Komplexe eine unterschiedliche Zahl von relevanten Bindungsstrukturen (Zustände). Für ein System mit relativ kurzen Alkylketten findet man zwei Zustände, eine kompakte Struktur, die auch im Gleichgewicht beobachtet wird und eine gestreckte Struktur, die nur unter dem Einfluss der externen Kraft stabil ist. Verlängert man die Alkylketten,rnbeobachtet man einen weiteren Zustand, in dem das Dimer vollständig gestreckt ist und die Monomere eine größere Separation aufweisen.rnBei dem zweiten System, das untersucht wurde, handelte es sich um einen Carbohydrat-Kation-Carbohydrat Komplex, der für die Selbstadhäsion von Meeresschwämmen eine wichtige Rolle spielt. Experimentell ist bekannt, dass sich dieser Komplex zwar mit Calciumionen, nicht aber mit Magnesiumionen bildet. Im Rahmen dieser Arbeit wurde gezeigt, dass die wesentlichen Unterschiede der beiden Kationarten in Bezug auf die Komplexbildung auf den kleineren Ionenradius des Magnesiumions zurückzuführen sind. Aufgrund des kleineren Radius bindet ein solvatisiertes Magnesiumion die Hydrathülle stärker und die Komplexbindung wird kinetisch gehemmt. Zum anderen bindet im Magnesiumkomplex nur eines der beiden Carbohydrate direkt an das Kation.rnDas andere Carbohydrat bindet nur indirekt über ein Wassermolekül an das Kation. Da diese indirekte Bindung gegenüber einer direkten Bindung schwächer ist, weist der Magensiumkomplex eine geringere Stabilität auf als ein vergleichbarer Calciumkomplex.rnDes Weiteren wurde untersucht, inwieweit die Ergebnisse von MD Simulationen vom verwendeten Modell (Kraftfeld) abhängen. Allgemein ist bekannt, dass die Ergebnisse von Gleichgewichtssimulationen kraftfeldabhängig sind. Im Rahmen diese Arbeit konnte gezeigt werden, dass sich für Zugsimulationen, in denen eine externe Kraft an das System angelegt wird, eine ähnliche Kraftfeldabhängigkeit ergibt. Da sich die Unterschiede der Ergebnisse auf Unterschiede in den Gleichgewichtssimulationen zurückführen lassen, kann man annehmen, dass die externe Kraft keine zusätzliche Einschränkung in Bezug auf die Zuverlässigkeit der Kraftfelder darstellt.rnAbgesehen von den MD Simulationen wurde eine in der Literatur beschriebene Methode zur Analyse von Zweizustandssystemen unter dem Einfluss einer konstanten externen Kraft erweitert. Ein Komplex läßt sich als Zweizustandssystem beschreiben, wenn dieser zwei relevante Bindungsstrukturen aufweist. Wird an solch einen Komplex eine konstante Kraft angelegt, lassen sich Übergänge zwischen den beiden Zuständen beobachten. Ist das System weit entfernt vom Gleichgewicht, kann es problematisch sein, einen der beiden Übergänge vollständig aufzulösen. In diesen Fällen wird nun vorgeschlagen, die beiden Übergänge zu einem sogenannten Kreisübergang zusammen zu fassen und diese zu zählen. Bestimmt man die Zahl der Übergänge pro Zeit in Abhängigkeit der angelegten Kraft, können die Übergangsraten bestimmt werden. Um die Methode zu validieren wurden kinetische Monte-Carlo Simulationen durchgeführt. Es zeigt sich, dass schon mit relativ kleinen Datensätzen gute Ergebnisse erzielt werden können.

Cicero, Tullia, and Marcus: Gender-Specific Concerns for Family Tradition?

A fair number of Cicero's letters reveal his concern for his daughter Tullia and his son Marcus. Recent scholarship has read these letters as evidence for a ‘natural’ emotional attachment of a father to his children, in reaction to Philippe Ariès's opposite claim. This chapter considers whether Cicero's letters can be analysed only as expressions of paternal affection. The fact that the pater familias Cicero occupies a political position simultaneously in his nuclear family, his domus, and the Senate, results in a concern for his prestige within the social field of the aristocracy. And this concern is necessarily conferred upon his support of the education and the social and political career of his children. The chapter traces the gender-specific differences between Cicero's treatment of Tullia and Marcus, shows the social construction of parental affection, and contributes to a further understanding of the different functions of daughters and sons in the social force field of family memory.

Structure-activity relationships from in vitro efficacies of the thiazolide series against the intracellular apicomplexan protozoan Neospora caninum

Nitazoxanide (2-acetolyloxy-N-(5-nitro 2-thiazolyl) benzamide; NTZ) represents the parent compound of a novel class of broad-spectrum anti-parasitic compounds named thiazolides. NTZ is active against a wide variety of intestinal and tissue-dwelling helminths, protozoa, enteric bacteria and a number of viruses infecting animals and humans. While potent, this poses a problem in practice, since this obvious non-selectivity can lead to undesired side effects in both humans and animals. In this study, we used real time PCR to determine the in vitro activities of 29 different thiazolides (NTZ-derivatives), which carry distinct modifications on both the thiazole- and the benzene moieties, against the tachyzoite stage of the intracellular protozoan Neospora caninum. The goal was to identify a highly active compound lacking the undesirable nitro group, which would have a more specific applicability, such as in food animals. By applying self-organizing molecular field analysis (SOMFA), these data were used to develop a predictive model for future drug design. SOMFA performs self-alignment of the molecules, and takes into account the steric and electrostatic properties, in order to determine 3D-quantitative structure activity relationship models. The best model was obtained by overlay of the thiazole moieties. Plotting of predicted versus experimentally determined activity produced an r2 value of 0.8052 and cross-validation using the "leave one out" methodology resulted in a q2 value of 0.7987. A master grid map showed that large steric groups at the R2 position, the nitrogen of the amide bond and position Y could greatly reduce activity, and the presence of large steric groups placed at positions X, R4 and surrounding the oxygen atom of the amide bond, may increase the activity of thiazolides against Neospora caninum tachyzoites. The model obtained here will be an important predictive tool for future development of this important class of drugs.

Molecular modeling of EPON 862-DETDA polymer

EPON 862 is an epoxy resin which is cured with the hardening agent DETDA to form a crosslinked epoxy polymer and is used as a component in modern aircraft structures. These crosslinked polymers are often exposed to prolonged periods of temperatures below glass transition range which cause physical aging to occur. Because physical aging can compromise the performance of epoxies and their composites and because experimental techniques cannot provide all of the necessary physical insight that is needed to fully understand physical aging, efficient computational approaches to predict the effects of physical aging on thermo-mechanical properties are needed. In this study, Molecular Dynamics and Molecular Minimization simulations are being used to establish well-equilibrated, validated molecular models of the EPON 862-DETDA epoxy system with a range of crosslink densities using a united-atom force field. These simulations are subsequently used to predict the glass transition temperature, thermal expansion coefficients, and elastic properties of each of the crosslinked systems for validation of the modeling techniques. The results indicate that glass transition temperature and elastic properties increase with increasing levels of crosslink density and the thermal expansion coefficient decreases with crosslink density, both above and below the glass transition temperature. The results also indicate that there may be an upper limit to crosslink density that can be realistically achieved in epoxy systems. After evaluation of the thermo-mechanical properties, a method is developed to efficiently establish molecular models of epoxy resins that represent the corresponding real molecular structure at specific aging times. Although this approach does not model the physical aging process, it is useful in establishing a molecular model that resembles the physically-aged state for further use in predicting thermo-mechanical properties as a function of aging time. An equation has been predicted based on the results which directly correlate aging time to aged volume of the molecular model. This equation can be helpful for modelers who want to study properties of epoxy resins at different levels of aging but have little information about volume shrinkage occurring during physical aging.

DENSITY OF AMORPHOUS CARBON BY USING DENSITY FUNCTIONAL THEORY

Amorphous carbon has been investigated for a long time. Since it has the random orientation of carbon atoms, its density depends on the position of each carbon atom. It is important to know the density of amorphous carbon to use it for modeling advance carbon materials in the future. Two methods were used to create the initial structures of amorphous carbon. One is the random placement method by randomly locating 100 carbon atoms in a cubic lattice. Another method is the liquid-quench method by using reactive force field (ReaxFF) to rapidly decrease the system of 100 carbon atoms from the melting temperature. Density functional theory (DFT) was used to refine the position of each carbon atom and the dimensions of the boundaries to minimize the ground energy of the structure. The average densities of amorphous carbon structures created by the random placement method and the liquid-quench method are 2.59 and 2.44 g/cm3, respectively. Both densities have a good agreement with previous works. In addition, the final structure of amorphous carbon generated by the liquid-quench method has lower energy.

Preliminary results and documentation of sediment cores CRP2 and CRP-2A from the Ross Sea off Cape Roberts, Antarctica

The site for CRP-2, 14 km east of Cape Roberts (77.006°S; 163.719°E), was selected to overlap the early Miocene strata cored in nearby CRP-1, and to sample deeper into the east-dipping strata near the western margin ofe he Victoria Land Basin to investigate Palaeogene climatic and tectonic history. CRP-2 was cored from 5 to 57 mbsf (metres below the sea floor) (core recovery 91 %), with a deviation resulting in CRP-2A being cored at the same site. CRP-2A reached down to 624mbsf (recovery 95%), and to strata with an age of c. 33-35 Ma. Drilling took place from 16 October to 25 November 1998, on 2.0-2.2 m of sea ice and through 178 m of water. Core fractures and other physical properties, such as sonic velocity, density and magnetic susceptibility, were measured throughout the core. Down-hole logs for these and other properties were run from 63 to 167 mbsf and subsequently from 200 to 623 mbsf, although density and velocity data could be obtained only to 440 mbsf because of hole collapse. Sonic velocity averages c. 2.0 km S-1 for the upper part of the hole, but there is an sharp increase to c. 3.0 km s-1 and also a slight angular unconformity, at 306 mbsf, corresponding most likely to the early/late Oligocene boundary (c. 28-30 Ma). Velocity then increases irregularly to around 3.6 km s-1 at the bottom of the hole, which is estimated to lie 120 m above the V4/V5 boundary. The higher velocities below 306 mbsf probably reflect more extensive carbonate and common pyrite cementation, in patches, nodules, bedding-parallel masses and as vein infills. Dip of the strata also increases down-hole from 3° in the upper 300 in to over 10° at the bottom. Temperature gradient is 21° k-1. Over 2 000 fractures were logged through the hole. Borehole televiewer imagery was obtained for the interval from 200 to 440 mbsf to orient the fractures for stress field analysis. Lithostratigraphical descriptions on a scale of 1:20 are presented for the full length of the core, along with core box images, as a 200 page supplement to this issue. The hole initially passed through a layer of muddy gravel to 5.5 mbsf (Lithological Sub-Unit or LSU 1.1), and then into a Quaternary diatom-bearing clast-rich diamicton to 21 mbsf (LSU 2. l), with an interval of alternating compact diamicton and loose sand, and containing a rich Pliocene foraminiferal fauna, to 27 mbsf (LSU 2.2). The unit beneath this (LSU 3.1) has similar physical properties (sonic velocity, porosity, magnetic susceptibility) and includes diamictites of similar character to those of LSU 2.1 and 2.2, but an early Miocene (c. 19 Ma) diatom assemblage at 28 mbsf (top of LSU 3.1) shows that this sub-unit is part of the older section. The strata beneath 27 mbsf, primary target for the project, extend from early Miocene to perhaps latest Eocene age, and are largely cyclic glacimarine nearshore to offshore sediments. They are described as 41 lithological sub-units and interpreted in terms of 12 recurrent lithofacies. These are 1) mudstone, 2) inter-stratified mudstone and sandstone, 3) muddy very fine to coarse sandstone, 4) well-sorted stratified fine sandstone, 5) moderately to well-sorted, medium-grained sandstone, 6) stratified diamictite, 7) massive diamictite, 8) rhythmically inter-stratified sandstone and mudstone, 9) clast-supported conglomerate, 10) matrix-supported conglomerate, 11) mudstone breccia and 12) volcaniclastic sediment. Sequence stratigraphical analysis has identified 22 unconformity-bounded depositional sequences in pre- Pliocene strata. They typically comprise a four-part architecture involving, in ascending order, 1) a sharp-based coarse-grained unit (Facies 6,7,9 or 10), 2) a fining-upward succession of sandstones (Facies 3 and 4), 3) a mudstone interval (Facies l), in some cases coarsening upward to muddy sandstones (Facies 3), and 4) a sharp-based sandstone dominated succession (mainly Facies 4). The cyclicity recorded by the strata is interpreted in terms of a glacier ice margin retreating and advancing from land to the west, and of rises and falls in sea level. Analysis of sequence periodicity awaits afirmer chronology. However, apreliminary spectral analysis of magnetic susceptibility for a deepwater mudstone within one of the sequences (from 339 to 347 mbsf) reveals ratios between hierarchical levels that are similar to those of the three Milankovitch orbital forcing periodicities. The strata contain a wide range of fossils, the most abundant being marine diatoms. These commonly form up to 5% of the sediment, though in places the core is barren (notably between 300 and 412 mbsf). Fifty samples out of 250 reviewed were studied in detail. The assemblages define ten biostratigraphical zones, some of them based on local or as yet undescribed forms. The assemblages are neritic, and largely planktonic, suggesting that the sea floor was mostly below the photic zone throughout deposition of the corcd sequence. Calcareous nannofossils, representing incursions of ocean surface waters, are much less common (72 out of 183 samples examined) and restricted to mudstone intervals a few tens of metres thick, but are important for dating. Foraminifera are also sparse (73 out of 135 samples) and represented only by calcareous benthic species. Changing assemblages indicate a shift from inshore environments in the early Oligocenc to outer shelf in the late Oligocenc, returning to inshore in the early Miocene. Marine palynomorplis yielded large numbers of well-preserved forms from most of the 116 samples examined. The new in situ assemblagc found last year in CRP-1 is extended down into the late Oligocene and a further new assemblage is found in the early Oligoccnc. Many taxa are new, and cannot us yet contribute to an improved understanding of chronology or ecology. Marine invertebrate macrofossils, mostly molluscs and serpulid tubes, are scattered throughout the core. Preservation is good in mudstones but poor in other lithologies. Climate on land is reflected in the content of terrestrial palynomorphs, which are extremely scarce down to c. 300 mbsf. Some forms are reworked, and others represent a low growing sparse tundra with at least one species of Nothofagus. Beneath this level, a significantly greater diversity and abundance suggests a milder climate and a low diversity woody vegetation in the early Oligocene, but still far short of the richness found in known Eocene strata of the region. Sedimentary facies in the oldest strata also suggest a milder climate in the oldest strata cored, with indications of substantial glacial melt-water discharges, but are typical of a coldcr climate in late Oligocene and early Miocene times. Clast analyses from diamictites reveal weak to random fabrics, suggesting either lack of ice-contact deposition or post-depositional modification, but periods when ice grounded at the drill site are inferred from thin zones of in-situ brecciated rock and soft-sediment folding. These are more common above c. 300 mbsf, perhaps reflecting more extensive glacial advances during deposition of those strata. Erosion of the adjacent Transantarctic Mountains through Jurassic basalt and dolerite-intruded Beacon strata into basement rocks beneath is recorded by petrographical studies of clast and sand grain assemblages. Core below 310 mbsf contains a dominance of fine-grained Jurassic dolerite and basalt fragments along with Beacon-derived coal debris and rounded quartz grains, whereas the strata above this level have a much higher proportion of basement derived granitoids, implying that the large areas of the adjacent mountains had been eroded to basement by the end of the early Oligocene. There is little indication of rift-related volcanism below 310 mbsf. Above this, however, basaltic and trachytic tephras are common, especially from 280 to 200 mbsf, from 150 to 46 mbsf, and in Pliocene LSU 2.2 from 21 to 27 mbsf. The largest volcanic eruptions generated layers of coarse (up to 1 cm) trachytic pumice lapilli between 97 and 114 mbsf. The thickest of these (1.2 m at 112 mbsf) may have produced an eruptive column extending tens of km into the stratosphere. A source within a few tens of km of the drill site is considered most likely. Present age estimates for the pre-Pliocene sequence are based mainly on biostratigraphy (using mainly marine diatoms and to a lesser extent calcareous nannofossils), with the age of the tephra from 112 to 114 mbsf (21.44k0.05 Ma from 84 crystals by Ar-Ar) as a key reference point. Although there are varied and well-preserved microfossil assemblages through most of the sequence (notably of diatoms and marine palynomorphs), they comprise largely taxa either known only locally or as yet undescribed. In addition, sequence stratigraphical analysis and features in the core itself indicate numerous disconformities. The present estimate from diatom assemblages is that the interval from 27 to 130 mbsf is early Miocene in age (c. 19 to 23.5 Ma), consistent with the Ar-Ar age from 112 to 114 mbsf. Diatom assemblages also indicate that the late Oligocene epoch extends from c. 130 to 307 mbsf, which is supported by late Oligocene nannofossils from 130 to 185 mbsf. Strata from 307 to 412 mbsf have no age-diagnostic assemblages, but below this early Oligocene diatoms and nannofossils have been recovered. A nannoflora at the bottom of the hole is consistent with an earliest Oligocene or latest Eocene age. Magnetostratigraphical studies based on about 1000 samples, 700 of which have so far undergone demagnetisation treatment, have provided a polarity stratigraphy of 12 pre-Pliocene magnetozones. Samples above 270 mbsf are of consistently high quality. Below this, magnetic behaviour is more variable. A preliminary age-depth plot using the Magnetic Polarity Time Scale (MPTS) and constrained by biostratigraphical data suggests that episodes of relatively rapid sedimentation took place at CRP-2 during Oligocene times (c. 100 m/My), but that more than half of the record was lost in a few major and many minor disconformities. Age estimates from Sr isotopes in shell debris and further tephra dating are expected to lead to a better comparison with the MPTS. CRP-2/2A has recorded a history of subsidence of the Victoria Land Basin margin that is similar to that found in CIROS-170 km to the south, reflecting stability in both basin and the adjacent mountains in late Cenozoic times, but with slow net accumulation in the middle Cenozoic. The climatic indicators from both drill holes show a similar correspondence, indicating polar conditions for the Quaternary but with sub-polar conditions in the early Miocene-late Oligocene and indications of warmer conditions still in the early Oligocene. Correlation between the CRP-2A core and seismic records shows that seismic units V3 and V4, both widespread in the Victoria Land Basin, represent a period of fluctuating ice margins and glacimarine sedimentation. The next drill hole, CRP-3, is expected to core deep into V5 and extend this record of climate and tectonics still further back in time.