Deep Mantle Conditions Research Articles

Pressure stabilizes ferrous iron in bridgmanite under hydrous deep lower mantle conditions

Earth’s lower mantle is a potential water reservoir. The physical and chemical properties of the region are in part controlled by the Fe3+/ΣFe ratio and total iron content in bridgmanite. However, the water effect on the chemistry of bridgmanite remains unclear. We carry out laser-heated diamond anvil cell experiments under hydrous conditions and observe dominant Fe2+ in bridgmanite (Mg, Fe)SiO3 above 105 GPa under the normal geotherm conditions corresponding to depth > 2300 km, whereas Fe3+-rich bridgmanite is obtained at lower pressures. We further observe FeO in coexistence with hydrous NiAs-type SiO2 under similar conditions, indicating that the stability of ferrous iron is a combined result of H2O effect and high pressure. The stability of ferrous iron in bridgmanite under hydrous conditions would provide an explanation for the nature of the low-shear-velocity anomalies in the deep lower mantle. In addition, entrainment from a hydrous dense layer may influence mantle plume dynamics and contribute to variations in the redox conditions of the mantle.

Phase Relation and Equation of State of Iron‐Titanium Oxyhydroxides With α‐PbO2 Type Crystal Structure at Deep Mantle Conditions

AbstractPhase relation and equation of state of iron‐titanium oxyhydroxide with α‐PbO2‐type crystal structure (hydrous α‐phase) was investigated at pressure from middle upper mantle to mantle transition zone. The experiments were performed using multi‐anvil apparatus with in situ synchrotron X‐ray diffraction. The metal‐diamond sample container was used to maintain a closed system with respect to water while allowing X‐ray transmittance. Starting materials were mixture of reagent grade goethite and anatase with Fe:Ti = 1:1 and 1:3. The hydrous α‐phase was synthesized at 11–14 GPa and 800–900°C. The X‐ray diffraction data were obtained over a wide range of pressures from 5 to 22 GPa and temperatures from room temperature to 930°C. Using pressure‐volume‐temperature data of the hydrous α‐phase with Fe:Ti = 1:1 collected at room temperature to 800°C and at 8.9–20.8 GPa, we determined the isothermal bulk modulus (KT0 = 183 (1) GPa) and the thermal expansivity (α0 = 3.29 (1) × 10−5K−1) at ambient condition. The stability field of the hydrous α‐phase and phase relation of FeOOH‐TiO2 system at 900°C was well constrained. It was found that the hydrous α‐phase decomposes into baddeleyite‐type TiO2 + ε‐FeOOH at pressure of approximately 20–21 GPa, and into ilmenite + rutile at 5–6 GPa. This stability field equivalent to depth of 180–600 km in the subduction zone. Our results suggest that the hydrous α‐phase is an important water reservoir in the middle upper mantle to mantle transition zone.

Sound velocities and single-crystal elasticity of hydrous Fo90 olivine to 12 GPa

Nominally anhydrous minerals (NAMs) may contain significant amounts of water and constitute an important reservoir for mantle hydrogen. The colloquial term ‘water’ in NAMs is related to the presence of hydroxyl-bearing (OH−) point defects in their crystal structure, where hydrogen is bonded to lattice oxygen and is charge-balanced by cation vacancies. This hydrous component may therefore have substantial effects on the thermoelastic parameters of NAMs, comparable to other major crystal-chemical substitutions (e.g., Fe, Al). Assessment of water concentrations in natural minerals from mantle xenoliths indicates that olivine commonly stores ∼0–200 ppm of water. However, the lack of samples originating from depths exceeding ∼250 km coupled with the rapid diffusion of hydrogen in olivine at magmatic temperatures makes the determination of the olivine water content in the upper mantle challenging. On the other hand, numerous experimental data show that, at pressures and temperatures corresponding to deep upper mantle conditions, the water storage capacity of olivine increases to 0.2–0.5 wt%. Therefore, determining the elastic properties of olivine samples with more realistic water contents for deep upper mantle conditions may help in interpreting both seismic velocity anomalies in potentially hydrous regions of Earth's mantle as well as the observed seismic velocity and density contrasts across the 410-km discontinuity.Here, we report simultaneous single-crystal X-ray diffraction and Brillouin scattering experiments at room temperature up to 11.96(2) GPa on hydrous [0.20(3) wt% H2O] Fo90 olivine to assess its full elastic tensor, and complement these results with a careful re-analysis of all the available single-crystal elasticity data from the literature for anhydrous Fo90 olivine. While the bulk (K) and shear (G) moduli of hydrous Fo90 olivine are virtually identical to those of the corresponding anhydrous phase, their pressure derivatives K′ and G′ are slightly larger, although consistent within mutual uncertainties. We then defined linear relations between the water concentration in Fo90 olivine, the elastic moduli and their pressure derivatives, which were then used to compute the sound velocities of Fo90 olivine with higher degrees of hydration. Even for water concentrations as high as 0.5 wt%, the sound wave velocities of hydrous and anhydrous olivines were found to be identical within uncertainties at pressures corresponding to the base of the upper mantle. Contrary to previous claims, our data suggest that water in olivine is not seismically detectable, at least for contents consistent with deep upper mantle conditions. In addition to that, our data reveal that the hydration of olivine is unlikely to be a key factor in reconciling seismic velocity and density contrasts across the 410-km discontinuity with a pyrolitic mantle.

The Post‐Perovskite Transition in Fe‐ and Al‐Bearing Bridgmanite: Effects on Seismic Observables

AbstractThe primary phase of the Earth’s lower mantle, (Al, Fe)‐bearing bridgmanite, transitions to the post‐perovskite (PPv) phase at Earth’s deep mantle conditions. Despite extensive experimental and ab initio investigations, there are still important aspects of this transformation that need clarification. Here, we address this transition in (Al3+, Fe3+)‐, (Al3+)‐, (Fe2+)‐, and (Fe3+)‐bearing bridgmanite using ab initio calculations and validate our results against experiments on similar compositions. Consistent with experiments, our results show that the onset transition pressure and the width of the two‐phase region depend distinctly on the chemical composition: (a) Fe3+‐, Al3+‐, or (Al3+, Fe3+)‐alloying increases the transition pressure, while Fe2+‐alloying has the opposite effect; (b) in the absence of coexisting phases, the pressure‐depth range of the Pv‐PPv transition is likely too broad to cause a sharp D” discontinuity (<30 km); (c) the average Clapeyron slope of the two‐phase regions are consistent with previous measurements, calculations in MgSiO3, and inferences from seismic data. In addition, (d) we observe a softening of the bulk modulus in the two‐phase region. The consistency between our results and experiments gives us the confidence to proceed and examine this transition in aggregates with different compositions computationally, which will be fundamental for resolving the most likely chemical composition of the D" region by analyses of tomographic images.

Melting curve and transport properties of ammonia ice up to the deep mantle conditions of Uranus and Neptune

Motivated by the poor knowledge of the implication of planetary ice on the internal convection and anomalous magnetic field of ice giants, we perform extensive first-principles molecular dynamics simulation on ammonia ice to investigate its thermophysical behaviors at high temperature ($T$) and pressure ($P$) conditions relevant to the deep mantle of Neptune and Uranus. The melting curve, transport properties, and sound velocity of ammonia up to 350 GPa and 5500 K, spanning from fluid mixtures to a highly compressed superionic regime, are determined. A first-order phase transition of ammonia from the plastic state to the superionic state is observed explicitly, which is associated with the reduction of entropy. In contrast to previous predictions, the melting curve elucidates the existence of superionic ammonia in deep planetary interiors. Inspection of the reduction for transport properties and sound velocity along the isentrope of ice giants further evidences that superionic ammonia can sufficiently sustain the planetary dynamo mechanism and may contribute to the internal stratification which is responsible for the generation of the multipolar magnetic field. Finally, a comprehensive phase diagram of ammonia extended to a higher T-P regime is constructed, which provides a clear picture for studying the fundamental behavior and phase transition of ammonia in deep planets, and enhances our understanding of the interior structure and thermal convection of these planetary systems.

Ab initio thermal expansion and thermoelastic properties of ringwoodite (<i>γ</i>-Mg<sub>2</sub>SiO<sub>4</sub>) at mantle transition zone conditions

Abstract. Thermal convection in the Earth's mantle is driven by lateral variations in temperature and density, which are substantially controlled by the local volume thermal expansion of the constituent mineral phases. Ringwoodite is a major component of the lower mantle transition zone, but its thermal expansivity and thermoelastic properties are still affected by large uncertainties. Ambient thermal expansion coefficient (αV0​​​​​​​), for instance, can vary as much as 100 % according to different experimental investigations available from the literature. In this work, we perform ab initio density functional theory calculations of vibrational properties of spinel-structured Mg2SiO4 ringwoodite in order to provide reliable thermophysical data up to mantle transition zone conditions. Temperature- and pressure-dependent thermal expansivity has been obtained by phonon dispersion calculations in the framework of quasi-harmonic approximation (QHA) up to 25 GPa and 2000 K. Theoretical analysis of vibrational spectra reveals that accurate prediction of IR and silent modes, along with their relative mode Grüneisen parameters, is crucial to define thermal expansivity. A six-parameter analytical function is able to reproduce ab initio values fairly well in the whole investigated P–T range, i.e., αV(P,T)=(1.6033×10-5+8.839×10-9T+11.586×10-3T-1-6.055T-2+804.31T-3) ×exp⁡(-2.52×10-2P), with temperature in kelvin and pressure in gigapascal. Ab initio static and isothermal bulk moduli have been derived for ringwoodite along with their P, T and cross derivatives, i.e., K0 = 184.3 GPa, KT,300 K = 176.6 GPa, K0′ = 4.13, KT,300K′ = 4.16, ∂KT∂TP = −0.0233 GPa K−1 and ∂2KT∂P∂T0=1.0×10-4 K−1. Computed thermal expansivity and thermoelastic properties support the evidence that QHA performs remarkably well for Mg2SiO4 ringwoodite up to mantle transition zone temperatures. Since volume thermal expansion of ringwoodite is strongly pressure-dependent and its pressure dependence becomes more marked with the increasing temperature, internally consistent assessments and empirical extrapolation of thermoelastic data to deep mantle conditions should be taken with care to avoid inaccurate or spurious predictions in phase equilibrium and mantle convection numerical modeling.

Garnet stability in the deep lunar mantle: Constraints on the physics and chemistry of the interior of the Moon

High-pressure, high-temperature experiments have been conducted at deep lunar mantle conditions to constrain the garnet stability field. Using the Taylor Whole Moon composition, garnet is found to be stable at pressures above 3 GPa and temperatures below 1700°C, yielding a smaller stability field than previously suggested on the basis of thermodynamic calculations. Experimental data are used to model equilibrium crystallization in a ‘two-stage’ model of lunar magma ocean (LMO) crystallization starting from a fully molten Moon. In the first stage, isothermal (1600°C) equilibrium crystallization of the LMO would produce garnet-bearing lherzolite cumulates (containing up to ∼20 wt.% garnet) in the lowermost lunar mantle. Garnet in the deep lunar mantle would significantly decrease the Al2O3 content of the residual LMO and impact HREE/MREE fractionation. Numerical modeling of the second stage (residual LMO fractional crystallization) shows a delay in plagioclase saturation compared to models of single-stage fractional crystallization of a whole-Moon LMO of the Taylor Whole Moon composition, thinning the anorthositic crust from 95 km to 75 km. To reach the upper limit of current estimates of the average lunar crustal thickness (∼45 km), the two-stage scenario needs to be accompanied by a total of 10% liquid trapped in cumulates and 70% efficiency of plagioclase flotation. We also conduct trace element evolution modeling and reproduce a REE pattern identical to high K, REE, and P (KREEP) compositions after 99.8% solidification, when starting with a CI chondritic REE abundance. The density of a garnet-bearing deep lunar mantle is significantly higher than the density of olivine/orthopyroxene mixtures without garnet. The present-day lowest mantle in the Moon could therefore be characterized by chemical interactions between the earliest (garnet-bearing) and latest (ilmenite-bearing) products of LMO crystallization.

Investigating Magma Ocean Solidification on Earth Through Laser-Heated Diamond Anvil Cell Experiments.

We carried out a series of silicate fractional crystallization experiments at lower mantle pressures using the laser‐heated diamond anvil cell. Phase relations and the compositional evolution of the cotectic melt and equilibrium solids along the liquid line of descent were determined and used to assemble the melting phase diagram. In a pyrolitic magma ocean, the first mineral to crystallize in the deep mantle is iron‐depleted calcium‐bearing bridgmanite. From the phase diagram, we estimate that the initial 33%–36% of the magma ocean will crystallize to form such a buoyant bridgmanite. Substantial calcium solubility in bridgmanite is observed up to 129 GPa, and significantly delays the crystallization of the calcium silicate perovskite phase during magma ocean solidification. Residual melts are strongly iron‐enriched as crystallization proceeds, making them denser than any of the coexisting solids at deep mantle conditions, thus supporting the terrestrial basal magma ocean hypothesis (Labrosse et al., 2007).

The influence of \u03b4-(Al,Fe)OOH on seismic heterogeneities in Earth\u2019s lower mantle

The high-pressure phases of oxyhydroxides (δ-AlOOH, ε-FeOOH, and their solid solution), candidate components of subducted slabs, have wide stability fields, thus potentially influencing volatile circulation and dynamics in the Earth’s lower mantle. Here, we report the elastic wave velocities of δ-(Al,Fe)OOH (Fe/(Al + Fe) = 0.13, δ-Fe13) to 79 GPa, determined by nuclear resonant inelastic X-ray scattering. At pressures below 20 GPa, a softening of the phonon spectra is observed. With increasing pressure up to the Fe3+ spin crossover (~ 45 GPa), the Debye sound velocity (vD) increases. At higher pressures, the low spin δ-Fe13 is characterized by a pressure-invariant vD. Using the equation of state for the same sample, the shear-, compressional-, and bulk-velocities (vS, vP, and vΦ) are calculated and extrapolated to deep mantle conditions. The obtained velocity data show that δ-(Al,Fe)OOH may cause low-vΦ and low-vP anomalies in the shallow lower mantle. At deeper depths, we find that this hydrous phase reproduces the anti-correlation between vS and vΦ reported for the large low seismic velocity provinces, thus serving as a potential seismic signature of hydrous circulation in the lower mantle.

Behavior and properties of water in silicate melts under deep mantle conditions

Water (H2O) as one of the most abundant fluids present in Earth plays crucial role in the generation and transport of magmas in the interior. Though hydrous silicate melts have been studied extensively, the experimental data are confined to relatively low pressures and the computational results are still rare. Moreover, these studies imply large differences in the way water influences the physical properties of silicate magmas, such as density and electrical conductivity. Here, we investigate the equation of state, speciation, and transport properties of water dissolved in Mg1−xFexSiO3 and Mg2(1−x)Fe2xSiO4 melts (for x = 0 and 0.25) as well as in its bulk (pure) fluid state over the entire mantle pressure regime at 2000–4000 K using first-principles molecular dynamics. The simulation results allow us to constrain the partial molar volume of the water component in melts along with the molar volume of pure water. The predicted volume of silicate melt + water solution is negative at low pressures and becomes almost zero above 15 GPa. Consequently, the hydrous component tends to lower the melt density to similar extent over much of the mantle pressure regime irrespective of composition. Our results also show that hydrogen diffuses fast in silicate melts and enhances the melt electrical conductivity in a way that differs from electrical conduction in the bulk water. The speciation of the water component varies considerably from the bulk water structure as well. Water is dissolved in melts mostly as hydroxyls at low pressure and as –O–H–O–, –O–H–O–H– and other extended species with increasing pressure. On the other hand, the pure water behaves as a molecular fluid below 15 GPa, gradually becoming a dissociated fluid with further compression. On the basis of modeled density and conductivity results, we suggest that partial melts containing a few percent of water may be gravitationally trapped both above and below the upper mantle-transition region. Moreover, such hydrous melts can give rise to detectable electrical conductance by means of electromagnetic sounding observations.

Reversal of carbonate-silicate cation exchange in cold slabs in Earth\u2019s lower mantle

The stable forms of carbon in Earth’s deep interior control storage and fluxes of carbon through the planet over geologic time, impacting the surface climate as well as carrying records of geologic processes in the form of diamond inclusions. However, current estimates of the distribution of carbon in Earth’s mantle are uncertain, due in part to limited understanding of the fate of carbonates through subduction, the main mechanism that transports carbon from Earth’s surface to its interior. Oxidized carbon carried by subduction has been found to reside in MgCO3 throughout much of the mantle. Experiments in this study demonstrate that at deep mantle conditions MgCO3 reacts with silicates to form CaCO3. In combination with previous work indicating that CaCO3 is more stable than MgCO3 under reducing conditions of Earth’s lowermost mantle, these observations allow us to predict that the signature of surface carbon reaching Earth’s lowermost mantle may include CaCO3.

Recycled calcium carbonate is an efficient oxidation agent under deep upper mantle conditions

Observations of high ferric iron content in diamond garnet inclusions and mantle plume melts suggest a highly heterogeneous distribution of ferric iron in the mantle. Recycling of oxidized materials such as carbonates from Earth’s surface by subduction could explain the observed variations. Here we present high-pressure high-temperature multi-anvil experiments to determine the redox reactions between calcium-, magnesium-, or iron-carbonate and ferrous iron-bearing silicate mineral (garnet or fayalite) at conditions representative of subduction zones with intermediate thermal structures. We show that both garnet and fayalite can be oxidized to ferric iron-rich garnets accompanied by reduction of calcium carbonate to form graphite. The ferric iron content in the synthetic garnets increases with increasing pressure, and is correlated with the Ca content in the garnets. We suggest that recycled sedimentary calcium carbonate could influence the evolution of the mantle oxidation state by efficiently increasing the ferric iron content in the deep upper mantle.

Spectroscopic evidence for the Fe3+ spin transition in iron-bearing δ-AlOOH at high pressure

Abstractδ-AlOOH has emerged as a promising candidate for water storage in the lower mantle and could have delivered water into the bottom of the mantle. To date, it still remains unclear how the presence of iron affects its elastic, rheological, vibrational, and transport properties, especially across the spin crossover. Here, we conducted high-pressure X-ray emission spectroscopy experiments on a δ-(Al0.85Fe0.15) OOH sample up to 53 GPa using silicone oil as the pressure transmitting medium in a diamond-anvil cell. We also carried out laser Raman measurements on δ-(Al0.85Fe0.15)OOH and δ-(Al0.52Fe0.48)OOH up to 57 and 62 GPa, respectively, using neon as the pressure-transmitting medium. Evolution of Raman spectra of δ-(Al0.85Fe0.15)OOH with pressure shows two new bands at 226 and 632 cm−1 at 6.0 GPa, in agreement with the transition from an ordered (P21nm) to a disordered hydrogen bonding structure (Pnnm) for δ-AlOOH. Similarly, the two new Raman bands at 155 and 539 cm−1 appear in δ-(Al0.52Fe0.48)OOH between 8.5 and 15.8 GPa, indicating that the incorporation of 48 mol% FeOOH could postpone the order-disorder transition upon compression. On the other hand, the satellite peak (Kβ′) intensity of δ-(Al0.85Fe0.15)OOH starts to decrease at ~30 GPa and it disappears completely at 42 GPa. That is, δ-(Al0.85Fe0.15)OOH undergoes a gradual electronic spin-pairing transition at 30–42 GPa. Furthermore, the pressure dependence of Raman shifts of δ-(Al0.85Fe0.15)OOH discontinuously decreases at 32–37 GPa, suggesting that the improved hydrostaticity by the use of neon pressure medium could lead to a relatively narrow spin crossover. Notably, the pressure dependence of Raman shifts and optical color of δ-(Al0.52Fe0.48)OOH dramatically change at 41–45 GPa, suggesting that it probably undergoes a relatively sharp spin transition in the neon pressure medium. Together with literature data on the solid solutions between δ-AlOOH and ε-FeOOH, we found that the onset pressure of the spin transition in δ-(Al,Fe)OOH increases with increasing FeOOH content. These results shed new insights into the effects of iron on the structural evolution and vibrational properties of δ-AlOOH. The presence of FeOOH in δ-AlOOH can substantially influence its high-pressure behavior and stability at the deep mantle conditions and play an important role in the deep-water cycle.

The Electrical Conductivity of Liebermannite: Implications for Water Transport Into the Earth's Lower Mantle

AbstractLiebermannite (KAlSi3O8) is a principal mineral phase expected to be thermodynamically stable in deeply subducted continental and oceanic crusts. The crystal structure of liebermannite exhibits tunnels that are formed between the assemblies of double chains of edge‐sharing (Si, Al) O6 octahedral units, which act as a repository for large incompatible alkali ions. In this study, we investigate the electrical conductivity of liebermannite at 12, 15, and 24 GPa and temperature of 1500 K to track subduction pathways of continental sediments into the Earth's lower mantle. Further, we looked at whether liebermannite could sequester incompatible H2O at deep mantle conditions. We observe that the superionic conductivity of liebermannite due to the thermally activated hopping of K+ ions results in high electrical conductivity of more than 1 S/m. Infrared spectral features of hydrous liebermannite indicate the presence of both molecular H2O and hydroxyl (OH−) groups in its crystal structure. The observed high electrical conductivity in the mantle transition zone beneath Northeastern China and the lower mantle beneath the Philippine Sea can be attributed to the subduction pathways of continental sediments deep into the Earth's mantle. While major mineral phases in pyrolitic compositions are almost devoid of H2O under lower mantle conditions, our study demonstrates that liebermannite could be an important host of H2O in these conditions. We propose that the relatively high H2O contents of ocean island basalts derived from deep mantle plumes are primarily related to deeply subducted continental sediments, in which liebermannite is the principal H2O carrier.

Genesis of the Ray-Iz chromitite, Polar Urals: Inferences to mantle conditions and recycling processes

Ophiolitic chromitites in the eastern sector of the Alpine-Himalayan and Polar Ural orogenic belts preserve evidence of formation under ultra-high-pressure (UHP), deep mantle conditions. The textural and compositional characteristics of microscopic and submicroscopic inclusions in magnesiochromite from the Ray-Iz ophiolites are considered as clues of their complex evolutionary history. Bulk-rock platinum-group-element (PGE) compositions of refractory and metallurgical chromitite ores (102–309 ppb) with generally IPGE>PPGE are consistent with supra-subduction zone (SSZ) metasomatism. The abundant unusual UHP phases (e.g., micro-diamonds, moissanite, coesite…etc.), together with clinopyroxene lamellae, and globular silicate inclusions in the Cr-spinel are seen as indications of deep mantle recycling of the mantle section of the Ray-Iz ophiolites. The calculated fO2 values for dunite and chromitite (+1.17 to +4.16, generally above the FMQ buffer) are in line with interaction of mid-ocean ridge (MOR) or back-arc ophiolites with oxidized, Mg-rich silicic (boninitic) melts in a SSZ environment.

Fate of Hydrous Fe‐Rich Silicate Melt in Earth's Deep Mantle

AbstractDensity of silicate melt dictates melt migration and establishes the gross structure of Earth's interior. However, due to technical challenges, the melt density of relevant compositions is poorly known at deep mantle conditions. Particularly, water may be dissolved in such melts in large amounts and can potentially affect their density at extreme pressure and temperature conditions. Here we perform first‐principles molecular dynamics simulations to evaluate the density of Fe‐rich, eutectic‐like silicate melt (E melt) with varying water content up to about 12 wt %. Our results show that water mixes nearly ideally with the nonvolatile component in silicate melt and can decrease the melt density significantly. They also suggest that hydrous melts can be gravitationally stable in the lowermost mantle given its likely high iron content, providing a mechanism to explain seismically slow and dense layers near the core‐mantle boundary.

Metastable silica high pressure polymorphs as structural proxies of deep Earth silicate melts

Modelling of processes involving deep Earth liquids requires information on their structures and compression mechanisms. However, knowledge of the local structures of silicates and silica (SiO2) melts at deep mantle conditions and of their densification mechanisms is still limited. Here we report the synthesis and characterization of metastable high-pressure silica phases, coesite-IV and coesite-V, using in situ single-crystal X-ray diffraction and ab initio simulations. Their crystal structures are drastically different from any previously considered models, but explain well features of pair-distribution functions of highly densified silica glass and molten basalt at high pressure. Built of four, five-, and six-coordinated silicon, coesite-IV and coesite-V contain SiO6 octahedra, which, at odds with 3rd Pauling’s rule, are connected through common faces. Our results suggest that possible silicate liquids in Earth’s lower mantle may have complex structures making them more compressible than previously supposed.

Pressure stability field of Mg-perovskite under deep mantle conditions: A topological approach based on Bader's analysis coupled with catastrophe theory

The pressure stability field of the Mg-perovskite phase was investigated by characterizing the evolution of the electron arrangement in the crystal. Ab initio calculations of the perovskite structures in the range 0–185 GPa were performed at the HF/DFT (Hartree-Fock/Density Functional Theory) exchange–correlation terms level. The electron densities, calculated throughout the ab-initio wave functions, were analysed by means of the Bader's theory, coupled with Thom's catastrophe theory. To the best of our knowledge the approach is used for the first time. The topological results show the occurrence of two topological anomalies at P~20 GPa and P~110 GPa which delineate the pressure range where Mg-perovskite is stable.The paper accomplishes the twofold objectives of providing a contribution in shading light into the behaviour of the dominant component of the Earth's lower mantle across the D’’ layer and of proposing a novel approach in predicting the stability of a compound at extreme conditions.

Synthesis and crystal structure of Mg-bearing Fe9O11: New insight in the complexity of Fe-Mg oxides at conditions of the deep upper mantle

A novel Mg-bearing iron oxide Mg0.87(1)Fe4.13(1)2+Fe43+O11 was synthesized at 12 GPa and 1300 °C using a large volume press. Rietveld structural analysis was conducted with a laboratory X‑ray diffraction pattern obtained at ambient conditions. The crystal structure, which has one oxygen trigonal prism site and four octahedral sites for the cations, was found to be isostructural with Ca2Fe7O11. The unit-cell lattice parameters are a = 9.8441(5) Å, b = 2.8920(1) Å, c = 14.1760(6) Å, β = 99.956(4)°, V = 397.50(3) Å3, and Z = 2 (monoclinic, C2/m). Mg and Fe cations are disordered on the trigonal prism site and on two of the four octahedral sites, and the remaining Fe is accommodated at the other two octahedral sites. The present structure is closely related to the other recently discovered Fe oxide structures, e.g., Fe4O5 and Fe5O6, by distortion derived either from incorporation (Fe4O5) or removal (Fe5O6) of an edge-shared FeO6 single octahedral chain in their structures. The present synthesis at deep upper mantle conditions and the structural relationships observed between various novel Mg-Fe oxides indicate that a series of different phases become stable above 10 GPa and that their relative stabilities (Fe2+/Fe3+) must be controlled by oxygen fugacity.

Single crystal growth, characterization and high-pressure Raman spectroscopy of impurity-free magnesite (MgCO3)

The understanding of the physical and chemical properties of magnesite (MgCO3) under deep-mantle conditions is highly important to capture the essence of deep-carbon storage in Earth’s interior. To develop standard rating scales, the impurity-free magnesite single crystal, paying particular attention to the case of avoiding adverse impacts of Ca2+, Fe2+, and Mn2+ impurities in natural magnesite, is undoubtedly necessary for all research of magnesite, including crystalline structural phase transitions, anisotropic elasticity and conductivity, and equation of state (EoS). Thus, a high-quality single crystal of impurity-free magnesite was grown successfully for the first time using the self-flux method under high pressure–temperature conditions. The size of the magnesite single crystal, observed in a plane-polarized microscope, exceeds 200 μm, and the crystal exhibits a rhombohedral structure to cleave along the (101) plane. In addition, its composition of Mg0.999 ± 0.001CO3 was quantified through electron probing analysis. The structural property was investigated by means of single crystal X-ray diffraction and the unit cell dimensions obtained in the rhombohedral symmetry of the $$R\bar {3}c$$ space group are a = 4.6255 (3) and c = 14.987 (2), and the final R = 0.0243 for 718 reflections. High-pressure Raman spectroscopy of the magnesite single crystal was performed up to 27 GPa at ambient temperature. All Raman active bands, ν i, without any splitting increased almost linearly with increasing pressure. In combination with the high-pressure Raman results $$\frac{{{\text{d}}{\nu _i}}}{{{\text{d}}P}}$$ and the bulk modulus K T (103 GPa) reported from magnesite EoS studies, the mode Gruneisen parameters (1.49, 1.40, 0.26, and 0.27) of each vibration (T, L, ν 4, and ν 1) were calculated.