Cost-Effective Force Field Tailored for Solid-Phase Simulations of OLED Materials

A united atom force field is empirically derived by minimizing the difference between experimental and simulated crystal cells and melting temperatures for eight compounds representative of organic electronic materials used in OLEDs and other devices: biphenyl, carbazole, fluorene, 9,9′-(1,3-phenylene)bis(9H-carbazole)-1,3-bis(N-carbazolyl)benzene (mCP), 4,4′-bis(N-carbazolyl)-1,1′-biphenyl (pCBP), phenazine, phenylcarbazole, and triphenylamine. The force field is verified against dispersion-corrected DFT calculations and shown to also successfully reproduce the crystal structure for two larger compounds employed as hosts in phosphorescent and thermally activated delayed fluorescence OLEDs: N,N′-di(1-naphthyl)-N,N′-diphenyl-(1,1′-biphenyl)-4,4′-diamine (NPD), and 1,3,5-tri(1-phenyl-1H-benzo[d]imidazol-2-yl)phenyl (TPBI). The good performances of the force field coupled to the large computational savings granted by the united atom approximation make it an ideal choice for the simulation of the morphology of emissive layers for OLED materials in crystalline or glassy phases.

Parallel simulation system for earthquake generation: Fault analysis modules and parallel coupling analysis

Solid earth simulations have recently been developed to address issues such as natural disasters, global environmental destruction and the conservation of natural resources. The simulation of solid earth phenomena involves the analysis of complex structures including strata, faults, and heterogeneous material properties. Simulation of the generation and cycle of earthquakes is particularly important, but such simulations require the analysis of complex fault dynamics. GeoFEM is a parallel finite-element analysis system intended for solid earth field phenomena problems. This paper describes recent development in the GeoFEM project for the simulation of earthquake generation and cycles.

A two-phase composite in simple shear: Effective mechanical anisotropy development and localization potential

We present a combined shape and mechanical anisotropy evolution model for a two-phase inclusion-bearing rock subject to large deformation. A single elliptical inclusion embedded in a homogeneous but anisotropic matrix is used to represent a simplified shape evolution enforced on all inclusions. The mechanical anisotropy develops due to the alignment of elongated inclusions. The effective anisotropy is quantified using the differential effective medium (DEM) approach. The model can be run for any deformation path and an arbitrary viscosity ratio between the inclusion and host phase. We focus on the case of simple shear and weak inclusions. The shape evolution of the representative inclusion is largely insensitive to the anisotropy development and to parameter variations in the studied range. An initial hardening stage is observed up to a shear strain of gamma = 1 irrespective of the inclusion fraction. The hardening is followed by a softening stage related to the developing anisotropy and its progressive rotation toward the shear direction. The traction needed to maintain a constant shear rate exhibits a fivefold drop at gamma = 5 in the limiting case of an inviscid inclusion. Numerical simulations show that our analytical model provides a good approximation to the actual evolution of a two-phase inclusion-host composite. However, the inclusions develop complex sigmoidal shapes resulting in the formation of an S-C fabric. We attribute the observed drop in the effective normal viscosity to this structural development. We study the localization potential in a rock column bearing varying fraction of inclusions. In the inviscid inclusion case, a strain jump from gamma = 3 to gamma = 100 is observed for a change of the inclusion fraction from 20% to 33%.

Atlantic 231Pa/230Th and biogenic opal compilation of the Holocene and the LGM

The strength and geometry of the Atlantic meridional overturning circulation is tightly coupled to climate on glacial-interglacial and millennial timescales, but has proved difficult to reconstruct, particularly for the Last Glacial Maximum. Today, the return flow from the northern North Atlantic to lower latitudes associated with the Atlantic meridional overturning circulation reaches down to approximately 4,000 m. In contrast, during the Last Glacial Maximum this return flow is thought to have occurred primarily at shallower depths. Measurements of sedimentary 231Pa/230Th have been used to reconstruct the strength of circulation in the North Atlantic Ocean, but the effects of biogenic silica on 231Pa/230Th-based estimates remain controversial. Here we use measurements of 231Pa/230Th ratios and biogenic silica in Holocene-aged Atlantic sediments and simulations with a two-dimensional scavenging model to demonstrate that the geometry and strength of the Atlantic meridional overturning circulation are the primary controls of 231Pa/230Th ratios in modern Atlantic sediments. For the glacial maximum, a simulation of Atlantic overturning with a shallow, but vigorous circulation and bulk water transport at around 2,000 m depth best matched observed glacial Atlantic 231Pa/230Th values. We estimate that the transport of intermediate water during the Last Glacial Maximum was at least as strong as deep water transport today.

Resonance studies of artificial earth satellites

Orbit determination from artificial satellite observations is a key process in obtaining information about the Earth and its environment. A study of the perturbations experienced by these satellites enables knowledge to be gained of the upper atmosphere, the gravity field, ocean tides, solid-Earth tides and solar radiation. The gravity field is expressed as a double infinite series of associated Legendre functions (tesseral harmonics). In contemporary global gravity field models the overall geoid is well determined. An independent check on these gravity field harmonics of a particular order may be made by analysis of satellites that pass through resonance of that order. For such satellites the perturbations of the orbital elements close to resonance are analysed to derive lumped harmonic coefficients. The orbital parameters of 1984-106A have been determined at 43 epochs, during which time the satellite was close to 14th order resonance. Analysis of the inclination and eccentricity yielded 6 lumped harmonic coefficients of order 14 whilst analysis of the mean motion yielded additional pairs of lumped harmonics of orders 14, 28 and 42, with the 14th order harmonics superseding those obtained from analysis of the inclination. This thesis concentrates in detail on the theoretical changes of a near-circular satellite orbit perturbed by the Earth's gravity field under the influence of minimal air-drag whilst in resonance with the Earth. The satellite 1984-106A experienced the interesting property of being temporarily trapped with respect to a secondary resonance parameter due to the low air-drag in 1987. This prompted the theoretical investigation of such a phenomenon. Expressions obtained for the resonance parameter led to the determination of 8 lumped harmonic coefficients, coincidental to those already obtained. All the derived lumped harmonic values arc used to test the accuracy of contemporary gravity field models and the underlying theory in this thesis.

Thick dyke emplacement and internal flow: A structural and magnetic fabric study of the deep-seated dolerite dyke of Foum Zguid (southern Morocco)

Knowledge on forced magma injection and magma flow in dykes is crucial for the understanding of how magmas migrate through the crust to the Earth's surface. Because many questions still persist, we used the long, thick, and deep-seated Foum Zguid dyke (Morocco) to investigate dyke emplacement and internal flow by means of magnetic methods, structural analysis, petrography, and scanning electron microscopy. We also investigated how the host rocks accommodated the intrusion. Regarding internal flow: 1. Important variations of the rock magnetic properties and magnetic fabric occur with distance from dyke wall; 2. anisotropy of anhysteretic remanent magnetization reveals that anisotropy of magnetic susceptibility (AMS) results mainly from the superposition of subfabrics with distinct coercivities and that the imbrication between magnetic foliation and dyke plane is more reliable to deduce flow than the orientation of the AMS maximum principal axis; and 3. a dominant upward flow near the margins can be inferred. The magnetic fabric closest to the dyke wall likely records magma flow best due to fast cooling, whereas in the core the magnetic properties have been affected by high-temperature exsolution and metasomatic effects due to slow cooling. Regarding dyke emplacement, this study shows that the thick forceful intrusion induced deformation by homogeneous flattening and/or folding of the host sedimentary strata. Dewatering related to heat, as recorded by thick quartz veins bordering the dyke in some localities, may have also helped accommodating dyke intrusion. The spatial arrangement of quartz veins and their geometrical relationship with the dyke indicate a preintrusive to synintrusive sinistral component of strike slip.

Deep structure of the Santos basin-São Paulo plateau system, SE Brazil

The structure and nature of the crust underlying the Santos Basin-São Paulo Plateau System (SSPS), in the SE Brazilian margin, are discussed based on five wide-angle seismic profiles acquired during the Santos Basin (SanBa) experiment in 2011. Velocity models allow us to precisely divide the SSPS in six domains from unthinned continental crust (Domain CC) to normal oceanic crust (Domain OC). A seventh domain (Domain D), a triangular shape region in the SE of the SSPS, is discussed by Klingelhoefer et al. (2014). Beneath the continental shelf, a similar to 100km wide necking zone (Domain N) is imaged where the continental crust thins abruptly from similar to 40km to less than 15km. Toward the ocean, most of the SSPS (Domains A and C) shows velocity ranges, velocity gradients, and a Moho interface characteristic of the thinned continental crust. The central domain (Domain B) has, however, a very heterogeneous structure. While its southwestern part still exhibits extremely thinned (7km) continental crust, its northeastern part depicts a 2-4km thick upper layer (6.0-6.5km/s) overlying an anomalous velocity layer (7.0-7.8km/s) and no evidence of a Moho interface. This structure is interpreted as atypical oceanic crust, exhumed lower crust, or upper continental crust intruded by mafic material, overlying either altered mantle in the first two cases or intruded lower continental crust in the last case. The deep structure and v-shaped segmentation of the SSPS confirm that an initial episode of rifting occurred there obliquely to the general opening direction of the South Atlantic Central Segment.

Magma flow pattern in dykes of the Azores revealed by ansiotropy of magnetic susceptibility

The localization of magma melting areas at the lithosphere bottom in extensional volcanic domains is poorly understood. Large polygenetic volcanoes of long duration and their associated magma chambers suggest that melting at depth may be focused at specific points within the mantle. To validate the hypothesis that the magma feeding a mafic crust, comes from permanent localized crustal reservoirs, it is necessary to map the fossilized magma flow within the crustal planar intrusions. Using the AMS, we obtain magmatic flow vectors from 34 alkaline basaltic dykes from São Jorge, São Miguel and Santa Maria islands in the Azores Archipelago, a hot-spot related triple junction. The dykes contain titanomagnetite showing a wide spectrum of solid solution ranging from Ti-rich to Ti-poor compositions with vestiges of maghemitization. Most of the dykes exhibit a normal magnetic fabric. The orientation of the magnetic lineation k1 axis is more variable than that of the k3 axis, which is generally well grouped. The dykes of São Jorge and São Miguel show a predominance of subhorizontal magmatic flows. In Santa Maria the deduced flow pattern is less systematic changing from subhorizontal in the southern part of the island to oblique in north. These results suggest that the ascent of magma beneath the islands of Azores is predominantly over localized melting sources and then collected within shallow magma chambers. According to this concept, dykes in the upper levels of the crust propagate laterally away from these magma chambers thus feeding the lava flows observed at the surface.

On the origin of the ultradeep East Barents Sea basin

Very large subsidence, with up to 20 km thick sediment layers, is observed in the East Barents Sea basin. Subsidence started in early Paleozoic, accelerated in Permo-Triassic times, finished during the middle Cretaceous, and was followed by moderate uplift in Cenozoic times. The observed gravity signal suggests that the East Barents Sea is at present in isostatic balance and indicates that a mass excess is required in the lithosphere to produce the observed large subsidence. Several origins have been proposed for the mass excess. We use 1-D thermokinematic modeling and 2-D isostatic density models of continental lithosphere to evaluate these competing hypotheses. The crustal density in 2-D thermokinematic models resulting from pressure-, temperature-, and composition-dependent phase change models is computed along transects crossing the East Barents Sea. The results indicate the following. (1) Extension can only explain the observed subsidence provided that a 10 km thick serpentinized mantle lens beneath the basin center is present. We conclude that this is unlikely given that this highly serpentinized layer should be formed below a sedimentary basin with more than 10 km of sediments and crust at least 10 km thick. (2) Phase changes in a compositionally homogeneous crust do not provide enough mass excess to explain the present-day basin geometry. (3) Phase change induced densification of a preexisting lower crustal gabbroic body, interpreted as a mafic magmatic underplate, can explain the basin geometry and observed gravity anomalies. The following model is proposed for the formation of the East Barents Sea basin: (1) Devonian rifting and extension related magmatism resulted in moderate thinning of the crust and a mafic underplate below the central basin area explaining initial late Paleozoic subsidence. (2) East-west shortening during the Permian and Triassic resulted in densification of the previously emplaced mafic underplated body and enhanced subsidence dramatically, explaining the present-day deep basin geometry.

Thermomechanical analysis of strain localization in a ductile detachment zone

At the latitude of the Thor-Odin dome (British Columbia) the Columbia River Detachment defines the eastern margin of the Shuswap metamorphic core complex and localizes in a 1 km thick muscovite-bearing quartzite mylonite. We present a combined Ar-40/Ar-39, (micro) structural, and oxygen isotope study of the deformation history in the detachment and evaluate the spatial and temporal relationships between microstructure formation and localization of strain. High-precision Ar-40/Ar-39 geochronology from different levels in the mylonite delineates a pattern of increasingly younger (49.0 to 47.9 Ma) deformation ages in deeper levels of the mylonitic footwall. The correlation of Ar-40/Ar-39 ages with decreasing deformation temperatures (similar to 550 degrees-400 degrees C) in the top 200 m of the mylonite indicates that deformation migrated downward from the contact with the hanging wall. Strain localization was diachronous in progressively deeper levels of the footwall and was likely controlled by fluid-assisted strain hardening due to advective heat removal and contemporaneous reaction weakening due to dissolution-reprecipitation of white mica. The observed constant high-stress microstructures across the entire detachment indicate that flow stress was buffered by the interplay of strain rate and temperature, where high strain rates at elevated temperature produced the same microstructure as lower strain rates under decreasing temperature conditions. The combined data suggest that the complex interplay among temporally nonuniform rates of footwall exhumation, heat advection, and embrittlement by meteoric fluids strongly determines the thermomechanical behavior of extensional detachments.

Multiple plume events in the genesis of the peri-Caribbean Cretaceous oceanic plateau province

The oceanic crust fragments exposed in central America, in north-western South America, and in the Caribbean islands have been considered to represent accreted remnants of the Caribbean-Colombian Oceanic Plateau (CCOP). On the basis of trace element and Nd, Sr, and Pb isotopic compositions we infer that cumulate rocks, basalts, and diabases from coastal Ecuador have a different source than the basalts from the Dominican Republic. The latter suite includes the 86 Ma basalts of the Duarte Complex which are light rare earth element (REE) -enriched and display (relative to normal mid-ocean ridge basalts, NMORB) moderate enrichments in large ion lithophile elements, together with high Nb, Ta, Pb, and low Th contents. Moreover, they exhibit a rather restricted range of Nd and Pb isotopic ratios consistent with their derivation from an ocean island-type mantle source, the composition of which includes the HIMU (high U-238/Pb-204) component characteristic of the Galapagos hotspot. In contrast, the 123 Ma Ecuadorian oceanic rocks have flat REE patterns and (relative to NMORB) are depleted in Zr, Hf, Th, and U. Moreover, they show a wide range of Nd and Pb isotopic ratios intermediate between those of ocean island basalts and NMORB. It is unlikely, on geochemical grounds, that the plume source of the Ecuadorian fragments was similar to that of the Galapagos. In addition, because of the NNE motion of the Farallon plate during the Early Cretaceous, the Ecuadorian oceanic plateau fragments could not have been derived from the Galapagos hotspot but were likely formed at a ridge-centered or near-ridge hotspot somewhere in the SE Pacific.