A role for subducting clays in the water transportation into the Earth's lower mantle

Subducting sedimentary layer typically contains water and hydrated clay minerals. The stability of clay minerals under such hydrous subduction environment would therefore constraint the lithology and physical properties of the subducting slab interface. Here we show that pyrophyllite (Al2Si4O10(OH)2), one of the representative clay minerals in the alumina-silica-water (Al2O3-SiO2-H2O, ASH) system, breakdowns to contain further hydrated minerals, gibbsite (Al(OH)3) and diaspore (AlO(OH)), when subducts along a water-saturated cold subduction geotherm. Such a hydration breakdown occurs at a depth of ~135 km to uptake water by ~1.8 wt%. Subsequently, dehydration breakdown occurs at ~185 km depth to release back the same amount of water, after which the net crystalline water content is preserved down to ~660 km depth, delivering a net amount of ~5.0 wt% H2O in a phase assemblage containing δ-AlOOH and phase Egg (AlSiO3(OH)). Our results thus demonstrate the importance of subducting clays to account the delivery of ~22% of water down to the lower mantle.

Structure Responsible for the Superconducting State in La3Ni2O7 at High-Pressure and Low-Temperature Conditions

Very recently, a new superconductor with Tc = 80 K has been reported in nickelate (La3Ni2O7) at around 15-40 GPa conditions (Nature, 621, 493, 2023), which is the second type of unconventional superconductor, besides cuprates, with Tc above liquid nitrogen temperature. However, the phase diagram plotted in this report was mostly based on the transport measurement under low-temperature and high-pressure conditions, and the assumed corresponding X-ray diffraction (XRD) results were carried out at room temperature. This encouraged us to carry out in situ high-pressure and low-temperature synchrotron XRD experiments to determine which phase is responsible for the high Tc state. In addition to the phase transition from the orthorhombic Amam structure to the orthorhombic Fmmm structure, a tetragonal phase with the space group of I4/mmm was discovered when the sample was compressed to around 19 GPa at 40 K where the superconductivity takes place in La3Ni2O7. The calculations based on this tetragonal structure reveal that the electronic states that approached the Fermi energy were mainly dominated by the eg orbitals (3dz2 and 3dx2-y2) of Ni atoms, which are located in the oxygen octahedral crystal field. The correlation between Tc and this structural evolution, especially Ni-O octahedra regularity and the in-plane Ni-O-Ni bonding angles, is analyzed. This work sheds new light to identify what is the most likely phase responsible for superconductivity in double-layered nickelate.

Submillisecond in situ X-ray diffraction measurement system with changing temperature and pressure using diamond anvil cells at BL10XU/SPring-8

Recently, there has been a high demand for elucidating kinetics and visualizing reaction processes under extreme dynamic conditions, such as chemical reactions under meteorite impact conditions, structural changes under nonequilibrium conditions, and in situ observations of dynamic changes. To accelerate material science studies and Earth science fields under dynamic conditions, a submillisecond in situ X-ray diffraction measurement system has been developed using a diamond anvil cell to observe reaction processes under rapidly changing pressure and temperature conditions replicating extreme dynamic conditions. The development and measurements were performed at the high-pressure beamline BL10XU/SPring-8 by synchronizing a high-speed hybrid pixel array detector, laser heating and temperature measurement system, and gas-pressure control system that enables remote and rapid pressure changes using the diamond anvil cell. The synchronized system enabled momentary heating and rapid cooling experiments up to 5000 K via laser heating as well as the visualization of structural changes in high-pressure samples under extreme dynamic conditions during high-speed pressure changes.

High-Pressure Synthesis of Light Lanthanide Dodecaborides (PrB12 and CeB12): Effects of Valence Fluctuation on Volume and Formation Pressure

Light lanthanide dodecaborides, RB₁₂ (R = Pr and Ce), were synthesized from a stoichiometric mixture of hexaborides and boron using a laser-heated diamond anvil cell under high-pressure and high-temperature conditions. Contrary to the expectation that lighter lanthanide elements require higher pressure to crystallize RB₁₂, in situ X-ray diffraction experiments reveal that cerium dodecaboride crystallizes at 26 GPa, which is significantly lower than that required to form the heavier praseodymium dodecaboride (35 GPa). In addition to the lower formation pressure, an anomalous volume reduction is also observed in CeB₁₂, which can be explained by a valence fluctuation between Ce³⁺ and Ce⁴⁺ indicated by X-ray absorption near-edge structure measurements. A polyhedral coordination change from a truncated cube in RB₆ to a truncated octahedron in RB₁₂ and associated shortening of the R-B bond length result in an increase in bulk modulus and hardness.

Viscosity of hcp iron at Earth's inner core conditions from density functional theory

The inner core, extending to 1,221 km above the Earth’s center at pressures between 329 and 364 GPa, is primarily composed of solid iron. Its rheological properties influence both the Earth’s rotation and deformation of the inner core which is a potential source of the observed seismic anisotropy. However, the rheology of the inner core is poorly understood. We propose a mineral physics approach based on the density functional theory to infer the viscosity of hexagonal close packed (hcp) iron at the inner core pressure (P) and temperature (T). As plastic deformation is rate-limited by atomic diffusion under the extreme conditions of the Earth’s center, we quantify self-diffusion in iron non-empirically. The results are applied to model steady-state creep of hcp iron. Here, we show that dislocation creep is a key mechanism driving deformation of hcp iron at inner core conditions. The associated viscosity agrees well with the estimates from geophysical observations supporting that the inner core is significantly less viscous than the Earth’s mantle. Such low viscosity rules out inner core translation, with melting on one side and solidification on the opposite, but allows for the occurrence of the seismically observed fluctuations in inner core differential rotation.

CO3+1 network formation in ultra-high pressure carbonate liquids

Carbonate liquids are an important class of molten salts, not just for industrial applications, but also in geological processes. Carbonates are generally expected to be simple liquids, in terms of ionic interactions between the molecular carbonate anions and metal cations, and therefore relatively structureless compared to more “polymerized” silicate melts. But there is increasing evidence from phase relations, metal solubility, glass spectroscopy and simulations to suggest the emergence of carbonate “networks” at length scales longer than the component molecular anions. The stability of these emergent structures are known to be sensitive to temperature, but are also predicted to be favoured by pressure. This is important as a recent study suggests that subducted surface carbonate may melt near the Earth’s transition zone (~44 km), representing a barrier to the deep carbon cycle depending on the buoyancy and viscosity of these liquids. In this study we demonstrate a major advance in our understanding of carbonate liquids by combining simulations and high pressure measurements on a carbonate glass, (K₂CO₃-MgCO₃) to pressures in excess of 40 GPa, far higher than any previous in situ study. We show the clear formation of extended low-dimensional carbonate networks of close CO₃²⁻ pairs and the emergence of a “three plus one” local coordination environment, producing an unexpected increase in viscosity with pressure. Although carbonate melts may still be buoyant in the lower mantle, an increased viscosity by at least three orders of magnitude will restrict the upward mobility, possibly resulting in entrainment by the down-going slab.

Compressed glassy carbon maintaining graphite-like structure with linkage formation between graphene layers

Amorphous diamond, formed by high-pressure compression of glassy carbon, is of interests for new carbon materials with unique properties such as high compressive strength. Previous studies attributed the ultrahigh strength of the compressed glassy carbon to structural transformation from graphite-like sp²-bonded structure to diamond-like sp³-bonded structure. However, there is no direct experimental determination of the bond structure of the compressed glassy carbon, because of experimental challenges. Here we succeeded to experimentally determine pair distribution functions of a glassy carbon at ultrahigh pressures up to 49.0 GPa by utilizing our recently developed double-stage large volume cell. Our results show that the C-C-C bond angle in the glassy carbon remains close to 120°, which is the ideal angle for the sp²-bonded honey-comb structure, up to 49.0 GPa. Our data clearly indicate that the glassy carbon maintains graphite-like structure up to 49.0 GPa. In contrast, graphene interlayer distance decreases sharply with increasing pressure, approaching values of the second neighbor C-C distance above 31.4 GPa. Linkages between the graphene layers may be formed with such a short distance, but not in the form of tetrahedral sp³ bond. The unique structure of the compressed glassy carbon may be the key to the ultrahigh strength.

The thermal expansion of gold: point defect concentrations and pre-melting in a face-centred cubic metal

On the basis of ab initio computer simulations, pre-melting phenomena have been suggested to occur in the elastic properties of hexagonal close-packed iron under the conditions of the Earth's inner core just before melting. The extent to which these pre-melting effects might also occur in the physical properties of face-centred cubic metals has been investigated here under more experimentally accessible conditions for gold, allowing for comparison with future computer simulations of this material. The thermal expansion of gold has been determined by X-ray powder diffraction from 40 K up to the melting point (1337 K). For the entire temperature range investigated, the unit-cell volume can be represented in the following way: a second-order Grüneisen approximation to the zero-pressure volumetric equation of state, with the internal energy calculated via a Debye model, is used to represent the thermal expansion of the 'perfect crystal'. Gold shows a nonlinear increase in thermal expansion that departs from this Grüneisen-Debye model prior to melting, which is probably a result of the generation of point defects over a large range of temperatures, beginning at T/Tm > 0.75 (a similar homologous T to where softening has been observed in the elastic moduli of Au). Therefore, the thermodynamic theory of point defects was used to include the additional volume of the vacancies at high temperatures ('real crystal'), resulting in the following fitted parameters: Q = (V0K0)/γ = 4.04 (1) × 10-18 J, V0 = 67.1671 (3) Å3, b = (K0' - 1)/2 = 3.84 (9), θD = 182 (2) K, (vf/Ω)exp(sf/kB) = 1.8 (23) and hf = 0.9 (2) eV, where V0 is the unit-cell volume at 0 K, K0 and K0' are the isothermal incompressibility and its first derivative with respect to pressure (evaluated at zero pressure), γ is a Grüneisen parameter, θD is the Debye temperature, vf, hf and sf are the vacancy formation volume, enthalpy and entropy, respectively, Ω is the average volume per atom, and kB is Boltzmann's constant.

Pressure-induced structural change in MgSiO3 glass at pressures near the Earth's core-mantle boundary

Knowledge of the structure and properties of silicate magma under extreme pressure plays an important role in understanding the nature and evolution of Earth's deep interior. Here we report the structure of MgSiO₃ glass, considered an analog of silicate melts, up to 111 GPa. The first (r1) and second (r2) neighbor distances in the pair distribution function change rapidly, with r1 increasing and r2 decreasing with pressure. At 53-62 GPa, the observed r1 and r2 distances are similar to the Si-O and Si-Si distances, respectively, of crystalline MgSiO₃ akimotoite with edge-sharing SiO₆ structural motifs. Above 62 GPa, r1 decreases, and r2 remains constant, with increasing pressure until 88 GPa. Above this pressure, r1 remains more or less constant, and r2 begins decreasing again. These observations suggest an ultrahigh-pressure structural change around 88 GPa. The structure above 88 GPa is interpreted as having the closest edge-shared SiO₆ structural motifs similar to those of the crystalline postperovskite, with densely packed oxygen atoms. The pressure of the structural change is broadly consistent with or slightly lower than that of the bridgmanite-to-postperovskite transition in crystalline MgSiO₃ These results suggest that a structural change may occur in MgSiO₃ melt under pressure conditions corresponding to the deep lower mantle.

Ultrahigh-pressure polyamorphism in GeO2 glass with coordination number >6

Knowledge of pressure-induced structural changes in glasses is important in various scientific fields as well as in engineering and industry. However, polyamorphism in glasses under high pressure remains poorly understood because of experimental challenges. Here we report new experimental findings of ultrahigh-pressure polyamorphism in GeO2 glass, investigated using a newly developed double-stage large-volume cell. The Ge-O coordination number (CN) is found to remain constant at ∼6 between 22.6 and 37.9 GPa. At higher pressures, CN begins to increase rapidly and reaches 7.4 at 91.7 GPa. This transformation begins when the oxygen-packing fraction in GeO2 glass is close to the maximal dense-packing state (the Kepler conjecture = ∼0.74), which provides new insights into structural changes in network-forming glasses and liquids with CN higher than 6 at ultrahigh-pressure conditions.

Pressure-induced stiffness of Au nanoparticles to 71 GPa under quasi-hydrostatic loading

The compressibility of nanocrystalline gold (n-Au, 20 nm) has been studied by x-ray total scattering using high-energy monochromatic x-rays in the diamond anvil cell under quasi-hydrostatic conditions up to 71 GPa. The bulk modulus, K0, of the n-Au obtained from fitting to a Vinet equation of state is ~196(3) GPa, which is about 17% higher than for the corresponding bulk materials (K0: 167 GPa). At low pressures (<7 GPa), the compression behavior of n-Au shows little difference from that of bulk Au. With increasing pressure, the compressive behavior of n-Au gradually deviates from the equation of state (EOS) of bulk gold. Analysis of the pair distribution function, peak broadening and Rietveld refinement reveals that the microstructure of n-Au is nearly a single-grain/domain at ambient conditions, but undergoes substantial pressure-induced reduction in grain size until 10 GPa. The results indicate that the nature of the internal microstructure in n-Au is associated with the observed EOS difference from bulk Au at high pressure. Full-pattern analysis confirms that significant changes in grain size, stacking faults, grain orientation and texture occur in n-Au at high pressure. We have observed direct experimental evidence of a transition in compressional mechanism for n-Au at ~20 GPa, i.e. from a deformation dominated by nucleation and motion of lattice dislocations (dislocation-mediated) to a prominent grain boundary mediated response to external pressure. The internal microstructure inside the nanoparticle (nanocrystallinity) plays a critical role for the macro-mechanical properties of nano-Au.

Multiple seismic reflectors in Earth's lowermost mantle

The modern view of Earth's lowermost mantle considers a D″ region of enhanced (seismologically inferred) heterogeneity bounded by the core-mantle boundary and an interface some 150-300 km above it, with the latter often attributed to the postperovskite phase transition (in MgSiO3). Seismic exploration of Earth's deep interior suggests, however, that this view needs modification. So-called ScS and SKKS waves, which probe the lowermost mantle from above and below, respectively, reveal multiple reflectors beneath Central America and East Asia, two areas known for subduction of oceanic plates deep into Earth's mantle. This observation is inconsistent with expectations from a thermal response of a single isochemical postperovskite transition, but some of the newly observed structures can be explained with postperovskite transitions in differentiated slab materials. Our results imply that the lowermost mantle is more complex than hitherto thought and that interfaces and compositional heterogeneity occur beyond the D″ region sensu stricto.

Implementation of micro-ball nanodiamond anvils for high-pressure studies above 6 Mbar

Since invention of the diamond anvil cell technique in the late 1950s for studying materials at extreme conditions, the maximum static pressure generated so far at room temperature was reported to be about 400 GPa. Here we show that use of micro-semi-balls made of nanodiamond as second-stage anvils in conventional diamond anvil cells drastically extends the achievable pressure range in static compression experiments to above 600 GPa. Micro-anvils (10–50 μm in diameter) of superhard nanodiamond (with a grain size below ∼50 nm) were synthesized in a large volume press using a newly developed technique. In our pilot experiments on rhenium and gold we have studied the equation of state of rhenium at pressures up to 640 GPa and demonstrated the feasibility and crucial necessity of the in situ ultra high-pressure measurements for accurate determination of material properties at extreme conditions.

The study of materials at high pressure has been limited by the conditions achievable using single-crystal diamond anvils. The use of anvils that incorporate a second stage consisting of two hemispherical nanocrystalline diamond micro-balls, extends the range of static pressures that can be generated in the lab.

Mineralogical effects on the detectability of the postperovskite boundary

The discovery of a phase transition in Mg-silicate perovskite (Pv) to postperovskite (pPv) at lowermost mantle pressure-temperature (P - T) conditions may provide an explanation for the discontinuous increase in shear wave velocity found in some regions at a depth range of 200 to 400 km above the core-mantle boundary, hereafter the D('') discontinuity. However, recent studies on binary and ternary systems showed that reasonable contents of Fe(2+) and Al for pyrolite increase the thickness (width of the mixed phase region) of the Pv - pPv boundary (400-600 km) to much larger than the D('') discontinuity (≤ 70 km). These results challenge the assignment of the D('') discontinuity to the Pv - pPv boundary in pyrolite (homogenized mantle composition). Furthermore, the mineralogy and composition of rocks that can host a detectable Pv → pPv boundary are still unknown. Here we report in situ measurements of the depths and thicknesses of the Pv → pPv transition in multiphase systems (San Carlos olivine, pyrolitic, and midocean ridge basaltic compositions) at the P - T conditions of the lowermost mantle, searching for candidate rocks with a sharp Pv - pPv discontinuity. Whereas the pyrolitic mantle may not have a seismologically detectable Pv → pPv transition due to the effect of Al, harzburgitic compositions have detectable transitions due to low Al content. In contrast, Al-rich basaltic compositions may have a detectable Pv - pPv boundary due to their distinct mineralogy. Therefore, the observation of the D('') discontinuity may be related to the Pv → pPv transition in the differentiated oceanic lithosphere materials transported to the lowermost mantle by subducting slabs.

Crystal structure and thermoelastic properties of (Mg0.91Fe0.09)SiO3 postperovskite up to 135 GPa and 2,700 K

Intriguing seismic observations have been made for the bottom 400 km of Earth's mantle (the D'' region) over the past few decades, yet the origin of these seismic structures has not been well understood. Recent theoretical calculations have predicted many unusual changes in physical properties across the postperovskite transition, perovskite (Pv) --> postperovskite (PPv), that may provide explanations for the seismic observations. Here, we report measurements of the crystal structure of (Mg(0.91)Fe(0.09))SiO(3)-PPv under quasi-hydrostatic conditions up to the pressure (P)-temperature (T) conditions expected for the core-mantle boundary (CMB). The measured crystal structure is in excellent agreement with the first-principles calculations. We found that bulk sound speed (V(Phi)) decreases by 2.4 +/- 1.4% across the PPv transition. Combined with the predicted shear-wave velocity (V(S)) increase, our measurements indicate that lateral variations in mineralogy between Pv and PPv may result in the anticorrelation between the V(Phi) and V(S) anomalies at the D'' region. Also, density increases by 1.6 +/- 0.4% and Grüneisen parameter decreases by 21 +/- 15% across the PPv transition, which will dynamically stabilize the PPv lenses observed in recent seismic studies.

Toward an internally consistent pressure scale

Our ability to interpret seismic observations including the seismic discontinuities and the density and velocity profiles in the earth's interior is critically dependent on the accuracy of pressure measurements up to 364 GPa at high temperature. Pressure scales based on the reduced shock-wave equations of state alone may predict pressure variations up to 7% in the megabar pressure range at room temperature and even higher percentage at high temperature, leading to large uncertainties in understanding the nature of the seismic discontinuities and chemical composition of the earth's interior. Here, we report compression data of gold (Au), platinum (Pt), the NaCl-B2 phase, and solid neon (Ne) at 300 K and high temperatures up to megabar pressures. Combined with existing experimental data, the compression data were used to establish internally consistent thermal equations of state of Au, Pt, NaCl-B2, and solid Ne. The internally consistent pressure scales provide a tractable, accurate baseline for comparing high pressure-temperature experimental data with theoretical calculations and the seismic observations, thereby advancing our understanding fundamental high-pressure phenomena and the chemistry and physics of the earth's interior.

Simultaneous determination of mean pressure and deviatoric stress based on numerical tensor analysis: a case study for polycrystalline x-ray diffraction of gold enclosed in a methanol-ethanol mixture

It is known that the {100} and {111} planes of cubic crystals subjected to uniaxial deviatoric stress conditions have strain responses that are free from the effect of lattice preferred orientation. By utilizing this special character, one can unambiguously and simultaneously determine the mean pressure and deviatoric stress from polycrystalline diffraction data of the cubic sample. Here we introduce a numerical tensor calculation method based on the generalized Hooke's law to simultaneously determine the hydrostatic component of the stress (mean pressure) and deviatoric stress in the sample. The feasibility of this method has been tested by examining the experimental data of the Au pressure marker enclosed in a diamond anvil cell using a pressure medium of methanol-ethanol mixture. The results demonstrated that the magnitude of the deviatoric stress is ∼0.07 GPa at the mean pressure of 10.5 GPa, which is consistent with previous results of Au strength under high pressure. Our results also showed that even a small deviatoric stress (∼0.07 GPa) could yield a ∼0.3 GPa mean pressure error at ∼10 GPa.