Spin crossover in ferropericlase and velocity heterogeneities in the lower mantle

Full-waveform tomography reveals iron spin crossover in Earth's lower mantle

Three-dimensional models of Earth's seismic structure can be used to identify temperature-dependent phenomena, including mineralogical phase and spin transformations, that are obscured in 1-D spherical averages. Full-waveform tomography maps seismic wave-speeds inside the Earth in three dimensions, at a higher resolution than classical methods. By providing absolute wave speeds (rather than perturbations) and simultaneously constraining bulk and shear wave speeds over the same frequency range, it becomes feasible to distinguish variations in temperature from changes in composition or spin state. We present a quantitative joint interpretation of bulk and shear wave speeds in the lower mantle, using a recently published full-waveform tomography model. At all depths the diversity of wave speeds cannot be explained by an isochemical mantle. Between 1000 and 2500 km depth, hypothetical mantle models containing an electronic spin crossover in ferropericlase provide a significantly better fit to the wave-speed distributions, as well as more realistic temperatures and silica contents, than models without a spin crossover. Below 2500 km, wave speed distributions are explained by an enrichment in silica towards the core-mantle boundary. This silica enrichment may represent the fractionated remains of an ancient basal magma ocean.

Low-spin state of Fe in Fe-doped NiOOH electrocatalysts

Doping with Fe boosts the electrocatalytic performance of NiOOH for the oxygen evolution reaction (OER). To understand this effect, we have employed state-of-the-art electronic structure calculations and thermodynamic modeling. Our study reveals that at low concentrations Fe exists in a low-spin state. Only this spin state explains the large solubility limit of Fe and similarity of Fe-O and Ni-O bond lengths measured in the Fe-doped NiOOH phase. The low-spin state renders the surface Fe sites highly active for the OER. The low-to-high spin transition at the Fe concentration of ~ 25% is consistent with the experimentally determined solubility limit of Fe in NiOOH. The thermodynamic overpotentials computed for doped and pure materials, η = 0.42 V and 0.77 V, agree well with the measured values. Our results indicate a key role of the low-spin state of Fe for the OER activity of Fe-doped NiOOH electrocatalysts.

Structural transition and re-emergence of iron's total electron spin in (Mg,Fe)O at ultrahigh pressure

Fe-bearing MgO [(Mg_1−xFe_x)O] is considered a major constituent of terrestrial exoplanets. Crystallizing in the B1 structure in the Earth’s lower mantle, (Mg_1−xFe_x)O undergoes a high-spin (S = 2) to low-spin (S = 0) transition at ∼45 GPa, accompanied by anomalous changes of this mineral’s physical properties, while the intermediate-spin (S = 1) state has not been observed. In this work, we investigate (Mg_1−xFe_x)O (x ≤ 0.25) up to 1.8 TPa via first-principles calculations. Our calculations indicate that (Mg_1−xFe_x)O undergoes a simultaneous structural and spin transition at ∼0.6 TPa, from the B1 phase low-spin state to the B2 phase intermediate-spin state, with Fe’s total electron spin S re-emerging from 0 to 1 at ultrahigh pressure. Upon further compression, an intermediate-to-low spin transition occurs in the B2 phase. Depending on the Fe concentration (x), metal–insulator transition and rhombohedral distortions can also occur in the B2 phase. These results suggest that Fe and spin transition may affect planetary interiors over a vast pressure range.

Iron spin transition occurs at ultrahigh pressure. The total electron spin increases from 0 to 1 as the structural transition of (Mg,Fe)O occurs (~0.6 TPa) and drops back to 0 at higher pressure. Its effects on exoplanet interiors are anticipated.

Seismological expression of the iron spin crossover in ferropericlase in the Earth's lower mantle

The two most abundant minerals in the Earth’s lower mantle are bridgmanite and ferropericlase. The bulk modulus of ferropericlase (Fp) softens as iron d-electrons transition from a high-spin to low-spin state, affecting the seismic compressional velocity but not the shear velocity. Here, we identify a seismological expression of the iron spin crossover in fast regions associated with cold Fp-rich subducted oceanic lithosphere: the relative abundance of fast velocities in P- and S-wave tomography models diverges in the ~1,400-2,000 km depth range. This is consistent with a reduced temperature sensitivity of P-waves throughout the iron spin crossover. A similar signal is also found in seismically slow regions below ~1,800 km, consistent with broadening and deepening of the crossover at higher temperatures. The corresponding inflection in P-wave velocity is not yet observed in 1-D seismic profiles, suggesting that the lower mantle is composed of non-uniformly distributed thermochemical heterogeneities which dampen the global signature of the Fp spin crossover.

This study identifies the predicted seismic expression of the high-to-low iron spin crossover in the deep Earth mineral ferropericlase. A depth-dependent signal is detected in the fastest and slowest regions, related to lateral temperature variations, of several global seismic tomography models.

Mineralogy of the deep lower mantle in the presence of H2O

Understanding the mineralogy of the Earth's interior is a prerequisite for unravelling the evolution and dynamics of our planet. Here, we conducted high pressure-temperature experiments mimicking the conditions of the deep lower mantle (DLM, 1800-2890 km in depth) and observed surprising mineralogical transformations in the presence of water. Ferropericlase, (Mg, Fe)O, which is the most abundant oxide mineral in Earth, reacts with H2O to form a previously unknown (Mg, Fe)O2H x (x ≤ 1) phase. The (Mg, Fe)O2H x has a pyrite structure and it coexists with the dominant silicate phases, bridgmanite and post-perovskite. Depending on Mg content and geotherm temperatures, the transformation may occur at 1800 km for (Mg0.6Fe0.4)O or beyond 2300 km for (Mg0.7Fe0.3)O. The (Mg, Fe)O2H x is an oxygen excess phase that stores an excessive amount of oxygen beyond the charge balance of maximum cation valences (Mg2+, Fe3+ and H+). This important phase has a number of far-reaching implications including extreme redox inhomogeneity, deep-oxygen reservoirs in the DLM and an internal source for modulating oxygen in the atmosphere.

Experimental evidence for silica-enriched Earth's lower mantle with ferrous iron dominant bridgmanite

Determination of the chemical composition of the Earth's mantle is of prime importance to understand the evolution, dynamics, and origin of the Earth. However, there is a lack of experimental data on sound velocity of iron-bearing Bridgmanite (Brd) under relevant high-pressure conditions of the whole mantle, which prevents constraints on the mineralogical model of the lower mantle. To uncover these issues, we have conducted sound-velocity measurement of iron-bearing Brd in a diamond-anvil cell (DAC) up to 124 GPa using Brillouin scattering spectroscopy. Here we show that the sound velocities of iron-bearing Brd throughout the whole pressure range of lower mantle exhibit an apparent linear reduction with the iron content. Our data fit remarkably with the seismic structure throughout the lower mantle with Fe²⁺-enriched Brd, indicating that the greater part of the lower mantle could be occupied by Fe²⁺-enriched Brd. Our lower-mantle model shows a distinctive Si-enriched composition with Mg/Si of 1.14 relative to the upper mantle (Mg/Si = 1.25), which implies that the mantle convection has been inefficient enough to chemically homogenize the Earth's whole mantle.

Velocity and density characteristics of subducted oceanic crust and the origin of lower-mantle heterogeneities

Seismic heterogeneities detected in the lower mantle were proposed to be related to subducted oceanic crust. However, the velocity and density of subducted oceanic crust at lower-mantle conditions remain unknown. Here, we report ab initio results for the elastic properties of calcium ferrite-type phases and determine the velocities and density of oceanic crust along different mantle geotherms. We find that the subducted oceanic crust shows a large negative shear velocity anomaly at the phase boundary between stishovite and CaCl₂-type silica, which is highly consistent with the feature of mid-mantle scatterers. After this phase transition in silica, subducted oceanic crust will be visible as high-velocity heterogeneities as imaged by seismic tomography. This study suggests that the presence of subducted oceanic crust could provide good explanations for some lower-mantle seismic heterogeneities with different length scales except large low shear velocity provinces (LLSVPs).

Evaluating the Role of Iron-Rich (Mg,Fe)O in Ultralow Velocity Zones

The composition of ultralow velocity zones (ULVZs) remains an open question, despite advances in both seismology and experimental work. We investigate the hypothesis of iron-rich (Mg,Fe)O (magnesiowüstite) as a cause of ULVZ seismic signatures. We report new quasi-hydrostatic X-ray diffraction measurements to constrain the equation of state of (Mg0.06Fe0.94)O with fit parameters V0 = 9.860 ± 0.007 Å3, K0T = 155.3 ± 2.2 GPa, K’0T = 3.79 ± 0.11, as well as synchrotron Mössbauer spectroscopy measurements to characterize the high-pressure magnetic and spin state of magnesiowüstite. We combine these results with information from previous studies to calculate the elastic behavior at core–mantle boundary conditions of magnesiowüstite, as well as coexisting bridgmanite and calcium silicate perovskite. Forward models of aggregate elastic properties are computed, and from these, we construct an inverse model to determine the proportions of magnesiowüstite that best reproduce ULVZ observations within estimated mutual uncertainties. We find that the presence of magnesiowüstite can explain ULVZ observations exhibiting 1:2 VP:VS reduction ratios relative to the Preliminary Reference Earth Model (PREM), as well as certain 1:3 VP:VS reductions within estimated uncertainty bounds. Our work quantifies the viability of compositionally distinct ULVZs containing magnesiowüstite and contributes to developing a framework for a methodical approach to evaluating ULVZ hypotheses.

Effects of vacancy defects on Fe properties incorporated in MgO

Distributions of Fe in MgO containing Mg vacancy, O vacancy, and Schottky defect are investigated based on the density functional theory (DFT). Our results show that since Mg vacancy will remove electrons from MgO, Fe tends to get close to Mg vacancy but far from O vacancy. The Mg vacancy can decrease the magnetic moment of iron and change its valence state from 2+  to 3+, which leads to ~5% decrease of Fe-O bond length comparable to the effect of 30 GPa external pressure. Furthermore, iron incorporation can increase the Schottky defect concentration of MgO especially in the environment of the Earth's lower mantle, where ~20 mol% Fe-bearing MgO locates at extreme high temperature conditions.

Abnormal Elasticity of Single-Crystal Magnesiosiderite across the Spin Transition in Earth's Lower Mantle

Brillouin light scattering and impulsive stimulated light scattering have been used to determine the full elastic constants of magnesiosiderite [(Mg_{0.35}Fe_{0.65})CO_{3}] up to 70 GPa at room temperature in a diamond-anvil cell. Drastic softening in C_{11}, C_{33}, C_{12}, and C_{13} elastic moduli associated with the compressive stress component and stiffening in C_{44} and C_{14} moduli associated with the shear stress component are observed to occur within the spin transition between ∼42.4 and ∼46.5  GPa. Negative values of C_{12} and C_{13} are also observed within the spin transition region. The Born criteria constants for the crystal remain positive within the spin transition, indicating that the mixed-spin state remains mechanically stable. Significant auxeticity can be related to the electronic spin transition-induced elastic anomalies based on the analysis of Poisson's ratio. These elastic anomalies are explained using a thermoelastic model for the rhombohedral system. Finally, we conclude that mixed-spin state ferromagnesite, which is potentially a major deep-carbon carrier, is expected to exhibit abnormal elasticity, including a negative Poisson's ratio of -0.6 and drastically reduced V_{P} by 10%, in Earth's midlower mantle.

Elasticity of Ferropericlase across the Spin Crossover in the Earth's Lower Mantle

Knowing the elasticity of ferropericlase across the spin transition can help explain seismic and mineralogical models of the lower-mantle including the origin of seismic heterogeneities in the middle to lowermost parts of the lower mantle1234. However, the effects of spin transition on full elastic constants of ferropericlase remain experimentally controversial due to technical challenges in directly measuring sound velocities under lower-mantle conditions12345. Here we have reliably measured both V_P and V_S of a single-crystal ferropericlase ((Mg_0.92,Fe_0.08)O) using complementary Brillouin Light Scattering and Impulsive Stimulated Light Scattering coupled with a diamond anvil cell up to 96 GPa. The derived elastic constants show drastically softened C₁₁ and C₁₂ within the spin transition at 40–60 GPa while C₄₄ is not affected. The spin transition is associated with a significant reduction of the aggregate V_P/V_S via the aggregate V_P softening because V_S softening does not visibly occur within the transition. Based on thermoelastic modelling along an expected geotherm, the spin crossover in ferropericlase can contribute to 2% reduction in V_P/V_S in a pyrolite mineralogical model in mid lower-mantle. Our results imply that the middle to lowermost parts of the lower-mantle would exhibit enhanced seismic heterogeneities due to the occurrence of the mixed-spin and low-spin ferropericlase.