Modeling viscosity of (Mg,Fe)O at lowermost mantle conditions

The viscosity of the lower mantle results from the rheological behavior of its two main constituent minerals, aluminous (Mg,Fe)SiO_3 bridgmanite and (Mg,Fe)O ferropericlase. Understanding the transport properties of lower mantle aggregates is of primary importance in geophysics and it is a challenging task, due to the extreme time-varying conditions to which such aggregates are subjected. In particular, viscosity is a crucial transport property that can vary over several orders of magnitude. It thus has a first-order control on the structure and dynamics of the mantle. Here we focus on the creep behavior of (Mg,Fe)O at the bottom of the lower mantle, where the presence of thermo-chemical anomalies such as ultralow-velocity zones (ULVZ) may significantly alter the viscosity contrast characterizing this region. Two different iron concentrations of (Mg_(1–x)Fe_x)O are considered: one mirroring the average composition of ferropericlase throughout most of the lower mantle (x = 0.20) and another representing a candidate magnesiowüstite component of ULVZs near the base of the mantle (x = 0.84). The investigated pressure-temperature conditions span from 120 GPa and 2800 K, corresponding to the average geotherm at this depth, to core-mantle boundary conditions of 135 GPa and 3800 K. In this study, dislocation creep of (Mg,Fe)O is investigated by dislocation dynamics (DD) simulations, a modeling tool which considers the collective motion and interactions of dislocations. To model their behavior, a 2.5 dimensional dislocation dynamics approach is employed. Within this method, both glide and climb mechanisms can be taken into account, and the interplay of these features results in a steady-state condition. This allows the retrieval of the creep strain rates at different temperatures, pressures, applied stresses and iron concentrations across the (Mg,Fe)O solid solution, providing information on the viscosity for these materials. A particularly low viscosity is obtained for magnesiowüstite with respect to ferropericlase, the difference being around 10 orders of magnitude. Thus, the final section of this work is devoted to the assessment of the dynamic implications of such a weak phase within ULVZs, in terms of the viscosity contrast with respect to the surrounding lowermost mantle

Static Safety for an Actor Dedicated Process Calculus by Abstract Interpretation

The actor model eases the definition of concurrent programs with non uniform behaviors. Static analysis of such a model was previously done in a data-flow oriented way, with type systems. This approach was based on constraint set resolution and was not able to deal with precise properties for communications of behaviors. We present here a new approach, control-flow oriented, based on the abstract interpretation framework, able to deal with communication of behaviors. Within our new analyses, we are able to verify most of the previous properties we observed as well as new ones, principally based on occurrence counting

Interpenetrated polymer networks of Poly(β-cyclodextrin) and Polyvinylpyrrolidone with synergistic and selective sorption capacities

Interpenetrating polymer network (IPN) hydrogels were synthesised using β-cyclodextrin (β-CD) and N-vynil-2-pyrrolidone (NVP) crosslinked with epichlorohydrin and divinylbenzene, respectively, and prepared by four different procedures: simultaneous, sequential, hybrid and a novel one named hybrid-sequential. The IPNs prepared have been characterised by infrared spectroscopy and thermal analysis. The equilibrium swelling in water and the sorption of model substances into the IPNs have also been studied. The model sorbates (1-naphthol, 2-acetylnaphthalene and tannic acid) were selected according to the affinities towards each one of the two constituent polymers. Our studies reveal that these IPNs can be applied for the sorption of substances that can interact with the network by two mechanisms, i.e. inclusion within cyclodextrin cavities and/or via specific interactions with the functional groups present. Besides, due to the complementary character of their constituent polymers, these networks could also serve to retain two substances of different nature such as cetirizine and pseudoephedrine

Dust amorphization in protoplanetary disks

High-energy irradiation of the circumstellar material might impact the structure and the composition of a protoplanetary disk and hence the process of planet formation. In this paper, we present a study on the possible influence of the stellar irradiation, indicated by X-ray emission, on the crystalline structure of the circumstellar dust. The dust crystallinity is measured for 42 class II T Tauri stars in the Taurus star-forming region using a decomposition fit of the 10 micron silicate feature, measured with the Spitzer IRS instrument. Since the sample includes objects with disks of various evolutionary stages, we further confine the target selection, using the age of the objects as a selection parameter. We correlate the X-ray luminosity and the X-ray hardness of the central object with the crystalline mass fraction of the circumstellar dust and find a significant anti-correlation for 20 objects within an age range of approx. 1 to 4.5 Myr. We postulate that X-rays represent the stellar activity and consequently the energetic ions of the stellar winds which interact with the circumstellar disk. We show that the fluxes around 1 AU and ion energies of the present solar wind are sufficient to amorphize the upper layer of dust grains very efficiently, leading to an observable reduction of the crystalline mass fraction of the circumstellar, sub-micron sized dust. This effect could also erase other relations between crystallinity and disk/star parameters such as age or spectral type.Comment: accepted for publication by A&

In-situ measurement of texture development rate in CaIrO₃ post-perovskite

The rate of crystallographic preferred orientation (CPO) development during deformation of post-perovskite is crucial in interpreting seismic anisotropy in the lowermost mantle but the stability field of MgSiO3 post-perovskite prevents high-strain deformation experiments being performed on it. Therefore, to constrain the rate of CPO development in post-perovskite, we deformed CaIrO3, a low-pressure analogue of MgSiO3 post-perovskite, in simple shear at 3.2GPa and 400○C to a shear strain (γ) of 0.81. From X-ray diffraction patterns acquired during deformation, we invert for CPO as a function of strain. By comparing the CPO that develops with visco-plastic self-consistent (VPSC) models we constrain the critical resolved shear stresses (CRSS) of the non-primary slip-systems in CaIrO3 to be of order 6 times stronger than the primary [100](010) slip system. This value is significantly less than has been assumed by previous studies and if applicable to MgSiO3 implies that seismic anisotropy in the D′ layer develops slower than has previously been assumed

Silicate dust in the environment of RS Ophiuchi following the 2006 eruption

We present further Spitzer Space Telescope observations of the recurrent nova RS Ophiuchi, obtained over the period 208-430 days after the 2006 eruption. The later Spitzer IRS data show that the line emission and free-free continuum emission reported earlier is declining, revealing incontrovertible evidence for the presence of silicate emission features at 9.7 and 18microns. We conclude that the silicate dust survives the hard radiation impulse and shock blast wave from the eruption. The existence of the extant dust may have significant implications for understanding the propagation of shocks through the red giant wind and likely wind geometry.Comment: 12 pages, 4 figures, accepted for publication in ApJ (Letters

Modeling viscosity of (Mg,Fe)O at lowermost mantle conditions

The viscosity of the lower mantle results from the rheological behavior of its two main constituent minerals, aluminous (Mg,Fe)SiO_3 bridgmanite and (Mg,Fe)O ferropericlase. Understanding the transport properties of lower mantle aggregates is of primary importance in geophysics and it is a challenging task, due to the extreme time-varying conditions to which such aggregates are subjected. In particular, viscosity is a crucial transport property that can vary over several orders of magnitude. It thus has a first-order control on the structure and dynamics of the mantle. Here we focus on the creep behavior of (Mg,Fe)O at the bottom of the lower mantle, where the presence of thermo-chemical anomalies such as ultralow-velocity zones (ULVZ) may significantly alter the viscosity contrast characterizing this region. Two different iron concentrations of (Mg_(1–x)Fe_x)O are considered: one mirroring the average composition of ferropericlase throughout most of the lower mantle (x = 0.20) and another representing a candidate magnesiowüstite component of ULVZs near the base of the mantle (x = 0.84). The investigated pressure-temperature conditions span from 120 GPa and 2800 K, corresponding to the average geotherm at this depth, to core-mantle boundary conditions of 135 GPa and 3800 K. In this study, dislocation creep of (Mg,Fe)O is investigated by dislocation dynamics (DD) simulations, a modeling tool which considers the collective motion and interactions of dislocations. To model their behavior, a 2.5 dimensional dislocation dynamics approach is employed. Within this method, both glide and climb mechanisms can be taken into account, and the interplay of these features results in a steady-state condition. This allows the retrieval of the creep strain rates at different temperatures, pressures, applied stresses and iron concentrations across the (Mg,Fe)O solid solution, providing information on the viscosity for these materials. A particularly low viscosity is obtained for magnesiowüstite with respect to ferropericlase, the difference being around 10 orders of magnitude. Thus, the final section of this work is devoted to the assessment of the dynamic implications of such a weak phase within ULVZs, in terms of the viscosity contrast with respect to the surrounding lowermost mantle

Upper- and mid-mantle interaction between the Samoan plume and the Tonga-Kermadec slabs

Mantle plumes are thought to play a key role in transferring heat from the core\u2013mantle boundary to the lithosphere, where it can significantly influence plate tectonics. On impinging on the lithosphere at spreading ridges or in intra-plate settings, mantle plumes may generate hotspots, large igneous provinces and hence considerable dynamic topography. However, the active role of mantle plumes on subducting slabs remains poorly understood. Here we show that the stagnation at 660 km and fastest trench retreat of the Tonga slab in Southwestern Pacific are consistent with an interaction with the Samoan plume and the Hikurangi plateau. Our findings are based on comparisons between 3D anisotropic tomography images and 3D petrological-thermo-mechanical models, which self-consistently explain several unique features of the Fiji\u2013Tonga region. We identify four possible slip systems of bridgmanite in the lower mantle that reconcile the observed seismic anisotropy beneath the Tonga slab (VSH4VSV) with thermo-mechanical calculations

The role of diffusion-driven pure climb creep on the rheology of bridgmanite under lower mantle conditions

The viscosity of Earth’s lower mantle is poorly constrained due to the lack of knowledge on some fundamental variables that affect the deformation behaviour of its main mineral phases. This study focuses on bridgmanite, the main lower mantle constituent, and assesses its rheology by developing an approach based on mineral physics. Following and revising the recent advances in this field, pure climb creep controlled by diffusion is identified as the key mechanism driving deformation in bridgmanite. The strain rates of this phase under lower mantle pressures, temperatures and stresses are thus calculated by constraining diffusion and implementing a creep theoretical model. The viscosity of MgSiO_3 bridgmanite resulting from pure climb creep is consequently evaluated and compared with the viscosity profiles available from the literature. We show that the inferred variability of viscosity in these profiles can be fully accounted for with the chosen variables of our calculation, e.g., diffusion coefficients, vacancy concentrations and applied stresses. A refinement of these variables is advocated in order to further constrain viscosity and match the observables