skip to main content

This content will become publicly available on April 1, 2025

Title: Hydrogen Diffusion in the Lower Mantle Revealed by Machine Learning Potentials

Hydrogen may be incorporated into nominally anhydrous minerals including bridgmanite and post‐perovskite as defects, making the Earth's deep mantle a potentially significant water reservoir. The diffusion of hydrogen and its contribution to the electrical conductivity in the lower mantle are rarely explored and remain largely unconstrained. Here we calculate hydrogen diffusivity in hydrous bridgmanite and post‐perovskite, using molecular dynamics simulations driven by machine learning potentials of ab initio quality. Our findings reveal that hydrogen diffusivity significantly increases with increasing temperature and decreasing pressure, and is considerably sensitive to hydrogen incorporation mechanism. Among the four defect mechanisms examined, (Mg + 2H)Siand (Al + H)Sishow similar patterns and yield the highest hydrogen diffusivity. Hydrogen diffusion is generally faster in post‐perovskite than in bridgmanite, and these two phases exhibit distinct diffusion anisotropies. Overall, hydrogen diffusion is slow on geological time scales and may result in heterogeneous water distribution in the lower mantle. Additionally, the proton conductivity of bridgmanite for (Mg + 2H)Siand (Al + H)Sidefects aligns with the same order of magnitude of lower mantle conductivity, suggesting that the water distribution in the lower mantle may be inferred by examining the heterogeneity of electrical conductivity.

more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Subject(s) / Keyword(s):
["Hydrogen","Diffusion","Lower mantle","Bridgmanite"]
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    High pressure and temperature experiments were carried out on the oxide mixtures corresponding to the bridgmanite stoichiometry under the hydrous shallow lower mantle conditions (24–25 GPa and 1673–1873 K with 5–10 wt. % of water in the starting material). Oxide mixtures investigated correspond to MgSiO3, (Mg, Fe)SiO3, (Mg, Al, Si)O3, and (Mg, Fe, Al, Si)O3. Melting was observed in all runs. Partitioning of various elements, including Mg, Fe, Si, and H is investigated. Melting under hydrous lower mantle conditions leads to increased (Mg + Fe)O/SiO2in the melt compared to the residual solids. The residual solids often contain a large amount of stishovite, and the melt contains higher (Mg,Fe)O/SiO2ratio than the initial material. (Mg + Fe)O‐rich hydrous melt could explain the low‐velocity anomalies observed in the shallow lower mantle and a large amount of stishovite in the residual solid may be responsible for the scattering of seismic waves in the mid‐lower mantle and may explain the “stishovite paradox. Since stishovite‐rich materials are formed only when silica‐rich source rock (MORB) is melted (not a typical peridotitic rock [bulk silicate Earth]), seismic scattering in the lower mantle provides a clue on the circulation of subducted MORB materials. To estimate hydrogen content, we use a new method of estimating the water content of unquenchable melts, and also propose a new interpretation of the significance of superhydrous phase B inclusions in bridgmanite. The results provide revised values of water partitioning between solid minerals and hydrous melts that are substantially higher than previous estimates.

    more » « less
  2. Abstract

    Thermal conductivity plays a pivotal role in understanding the dynamics and evolution of Earth's interior. The Earth's lower mantle is dominated by MgSiO3polymorphs which may incorporate trace amounts of water. However, the thermal conductivity of MgSiO3‐H2O binary system remains poorly understood. Here, we calculate the thermal conductivity of water‐free and water‐bearing bridgmanite, post‐perovskite, and MgSiO3melt, using a combination of Green‐Kubo method with molecular dynamics simulations based on a machine learning potential of ab initio quality. The thermal conductivities of water‐free bridgmanite and post‐perovskite overall agree well with previous theoretical and experimental studies. The presence of water mildly reduces the thermal conductivity of the host minerals, significantly weakens the temperature dependence of the thermal conductivity, and reduces the thermal anisotropy of post‐perovskite. Overall, water reduces the thermal conductivity difference between bridgmanite and post‐perovskite, and thus may attenuate lateral heterogeneities of the core‐mantle boundary heat flux.

    more » « less
  3. Abstract

    High‐quality single‐crystals of (Al,Fe)‐bearing bridgmanite, Mg0.88Fe3+0.065Fe2+0.035Al0.14Si0.90O3, of hundreds of micrometer size were synthesized at 24 GPa and 1800 °C in a Kawai‐type apparatus from the starting hydrous melt containing ~6.7 wt% water. Analyses of synthesized bridgmanite using petrographic microscopy, scanning electron microscopy, and transmission electron microscopy show that the crystals are chemically homogeneous and inclusion free in micrometer‐ to nanometer‐spatial resolutions. Nanosecondary ion mass spectrometry (NanoSIMS) analyses on selected platelets show ~1,020(±70) ppm wt water (hydrogen). The high water concentration in the structure of bridgmanite was further confirmed using polarized and unpolarized Fourier‐transform infrared spectroscopy (FTIR) analyses with two pronounced OH‐stretching bands at ~3,230 and ~3,460 cm−1. Our results indicate that lower‐mantle bridgmanite can accommodate relatively high amount of water. Therefore, dehydration melting at the topmost lower mantle by downward flow of transition zone materials would require water content exceeding ~0.1 wt%.

    more » « less
  4. null (Ed.)
    Abstract Electronic states of iron in the lower mantle's dominant mineral, (Mg,Fe,Al)(Fe,Al,Si)O3 bridgmanite, control physical properties of the mantle including density, elasticity, and electrical and thermal conductivity. However, the determination of electronic states of iron has been controversial, in part due to different interpretations of Mössbauer spectroscopy results used to identify spin state, valence state, and site occupancy of iron. We applied energy-domain Mössbauer spectroscopy to a set of four bridgmanite samples spanning a wide range of compositions: 10–50% Fe/total cations, 0–25% Al/total cations, 12–100% Fe3+/total Fe. Measurements performed in the diamond-anvil cell at pressures up to 76 GPa below and above the high to low spin transition in Fe3+ provide a Mössbauer reference library for bridgmanite and demonstrate the effects of pressure and composition on electronic states of iron. Results indicate that although the spin transition in Fe3+ in the bridgmanite B-site occurs as predicted, it does not strongly affect the observed quadrupole splitting of 1.4 mm/s, and only decreases center shift for this site to 0 mm/s at ~70 GPa. Thus center shift can easily distinguish Fe3+ from Fe2+ at high pressure, which exhibits two distinct Mössbauer sites with center shift ~1 mm/s and quadrupole splitting 2.4–3.1 and 3.9 mm/s at ~70 GPa. Correct quantification of Fe3+/total Fe in bridgmanite is required to constrain the effects of composition and redox states in experimental measurements of seismic properties of bridgmanite. In Fe-rich, mixed-valence bridgmanite at deep-mantle-relevant pressures, up to ~20% of the Fe may be a Fe2.5+ charge transfer component, which should enhance electrical and thermal conductivity in Fe-rich heterogeneities at the base of Earth's mantle. 
    more » « less
  5. Abstract

    Fe‐Al‐bearing bridgmanite may be the dominant host for ferric iron in Earth's lower mantle. Here we report the synthesis of (Mg0.5Fe3+0.5)(Al0.5Si0.5)O3bridgmanite (FA50) with the highest Fe3+‐Al3+coupled substitution known to date. X‐ray diffraction measurements showed that at ambient conditions, the FA50 adopted the LiNbO3structure. Upon compression at room temperature to 18 GPa, it transformed back into the bridgmanite structure, which remained stable up to 102 GPa and 2,600 K. Fitting Birch‐Murnaghan equation of state of FA50 bridgmanite yieldsV0 = 172.1(4) Å3,K0 = 229(4) GPa withK0′ = 4(fixed). The calculated bulk sound velocity of the FA50 bridgmanite is ~7.7% lower than MgSiO3bridgmanite, mainly because the presence of ferric iron increases the unit‐cell mass by 15.5%. This difference likely represents the upper limit of sound velocity anomaly introduced by Fe3+‐Al3+substitution. X‐ray emission and synchrotron Mössbauer spectroscopy measurements showed that after laser annealing, ~6% of Fe3+cations exchanged with Al3+and underwent the high‐ to low‐spin transition at 59 GPa. The low‐spin proportion of Fe3+increased gradually with pressure and reached 17–31% at 80 GPa. Since the cation exchange and spin transition in this Fe3+‐Al3+‐enriched bridgmanite do not cause resolvable unit‐cell volume reduction, and the increase of low‐spin Fe3+fraction with pressure occurs gradually, the spin transition would not produce a distinct seismic signature in the lower mantle. However, it may influence iron partitioning and isotopic fractionation, thus introducing chemical heterogeneity in the lower mantle.

    more » « less