Fate of MgSiO3 melts at core-mantle boundary conditions

Pressure dependent structure of amorphous magnesium aluminosilicates: The effect of replacing magnesia by alumina at the enstatite composition

The effect of replacing magnesia by alumina on the pressure-dependent structure of amorphous enstatite was investigated by applying in situ high-pressure neutron diffraction with magnesium isotope substitution to glassy (MgO)0.375(Al2O3)0.125(SiO2)0.5. The replacement leads to a factor of 2.4 increase in the rate-of-change of the Mg-O coordination number with pressure, which increases from 4.76(4) at ambient pressure to 6.51(4) at 8.2 GPa, and accompanies a larger probability of magnesium finding bridging oxygen atoms as nearest-neighbors. The Al-O coordination number increases from 4.17(7) to 5.24(8) over the same pressure interval at a rate that increases when the pressure is above ∼3.5 GPa. On recovering the glass to ambient conditions, the Mg-O and Al-O coordination numbers reduce to 5.32(4) and 4.42(6), respectively. The Al-O value is in accordance with the results from solid-state 27Al nuclear magnetic resonance spectroscopy, which show the presence of six-coordinated aluminum species that are absent in the uncompressed material. These findings explain the appearance of distinct pressure-dependent structural transformation regimes in the preparation of permanently densified magnesium aluminosilicate glasses. They also indicate an anomalous minimum in the pressure dependence of the bulk modulus with an onset that suggests a pressure-dependent threshold for transitioning between scratch-resistant and crack-resistant material properties.

Davemaoite as the mantle mineral with the highest melting temperature

Knowledge of high-pressure melting curves of silicate minerals is critical for modeling the thermal-chemical evolution of rocky planets. However, the melting temperature of davemaoite, the third most abundant mineral in Earth's lower mantle, is still controversial. Here, we investigate the melting curves of two minerals, MgSiO3 bridgmanite and CaSiO3 davemaoite, under their stability field in the mantle by performing first-principles molecular dynamics simulations based on the density functional theory. The melting curve of bridgmanite is in excellent agreement with previous studies, confirming a general consensus on its melting temperature. However, we predict a much higher melting curve of davemaoite than almost all previous estimates. Melting temperature of davemaoite at the pressure of core-mantle boundary (~136 gigapascals) is about 7700(150) K, which is approximately 2000 K higher than that of bridgmanite. The ultrarefractory nature of davemaoite is critical to reconsider many models in the deep planetary interior, for instance, solidification of early magma ocean and geodynamical behavior of mantle rocks.

Ultrafast x-ray detection of low-spin iron in molten silicate under deep planetary interior conditions

The spin state of Fe can alter the key physical properties of silicate melts, affecting the early differentiation and the dynamic stability of the melts in the deep rocky planets. The low-spin state of Fe can increase the affinity of Fe for the melt over the solid phases and the electrical conductivity of melt at high pressures. However, the spin state of Fe has never been measured in dense silicate melts due to experimental challenges. We report detection of dominantly low-spin Fe in dynamically compressed olivine melt at 150 to 256 gigapascals and 3000 to 6000 kelvin using laser-driven shock wave compression combined with femtosecond x-ray diffraction and x-ray emission spectroscopy using an x-ray free electron laser. The observation of dominantly low-spin Fe supports gravitationally stable melt in the deep mantle and generation of a dynamo from the silicate melt portion of rocky planets.

Local structural investigation of non-crystalline materials at high pressure: the case of GeO2glass

Local structures play a crucial role in the structural polyamorphism and novel electronic properties of amorphous materials, but their accurate measurement at high pressure remains a formidable challenge. In this article, we use the local structure of network-forming GeO₂glass as an example, to present our recent approaches and advances in high-energy x-ray diffraction, high-pressure x-ray absorption fine structure, andab initiofirst-principles density functional theory calculations and simulations. Although GeO₂glass is one of the best studied materials in the field of high pressure research due to its importance in glass theory and geophysical significance, there are still some long-standing puzzles, such as the existence of appreciable distinct fivefold^[5]Ge coordination at low pressure and the sixfold-plus^[6+]Ge coordination at ultrahigh pressure. Our work sheds light on the origin of pressure-induced polyamorphism of GeO₂glass, and the^[5]Ge polyhedral units may be the dominant species in the densification mechanism of network-forming glasses from tetrahedral to octahedral amorphous structures.

Electronic, Structural, and Mechanical Properties of SiO_{2} Glass at High Pressure Inferred from its Refractive Index

We report the first direct measurements of the refractive index of silica glass up to 145 GPa that allowed quantifying its density, bulk modulus, Lorenz-Lorentz polarizability, and band gap. These properties show two major anomalies at ∼10 and ∼40  GPa. The anomaly at ∼10  GPa signals the onset of the increase in Si coordination, and the anomaly at ∼40  GPa corresponds to a nearly complete vanishing of fourfold Si. More generally, we show that the compressibility and density of noncrystalline solids can be accurately measured in simple optical experiments up to at least 110 GPa.

Liquid structure under extreme conditions: high-pressure x-ray diffraction studies

Under extreme conditions of high pressure and temperature, liquids can undergo substantial structural transformations as their atoms rearrange to minimise energy within a more confined volume. Understanding the structural response of liquids under extreme conditions is important across a variety of disciplines, from fundamental physics and exotic chemistry to materials and planetary science.In situexperiments and atomistic simulations can provide crucial insight into the nature of liquid-liquid phase transitions and the complex phase diagrams and melting relations of high-pressure materials. Structural changes in natural magmas at the high-pressures experienced in deep planetary interiors can have a profound impact on their physical properties, knowledge of which is important to inform geochemical models of magmatic processes. Generating the extreme conditions required to melt samples at high-pressure, whilst simultaneously measuring their liquid structure, is a considerable challenge. The measurement, analysis, and interpretation of structural data is further complicated by the inherent disordered nature of liquids at the atomic-scale. However, recent advances in high-pressure technology mean that liquid diffraction measurements are becoming more routinely feasible at synchrotron facilities around the world. This topical review examines methods for high pressure synchrotron x-ray diffraction of liquids and the wide variety of systems which have been studied by them, from simple liquid metals and their remarkable complex behaviour at high-pressure, to molecular-polymeric liquid-liquid transitions in pnicogen and chalcogen liquids, and density-driven structural transformations in water and silicate melts.

Isothermal equation of state of crystalline and glassy materials from optical measurements in diamond anvil cells

Here, we present a method to study the equation of state of opaque amorphous and crystalline materials in diamond anvil cells. The approach is based on measurements of sample dimensions using high-resolution optical microscopy. Data on the volumetric strain as a function of pressure allow deriving the isothermal equation of state of the studied material. The analysis of optical images is fully automatized and allows measuring the sample dimensions with the precision of about 60 nm. The methodology was validated by studying isothermal compression of ω-Ti up to 30 GPa in a Ne pressure transmitting medium. Within the accuracy of the measurements, the bulk modulus of ω-Ti determined using optical microscopy was similar to that obtained from x-ray diffraction. For glassy carbon compressed to ∼30 GPa, the previously unknown bulk modulus was found to be equal to K₀ = 28 (2) GPa [K' = 5.5(5)].

In situ X-ray diffraction of silicate liquids and glasses under dynamic and static compression to megabar pressures

Properties of liquid silicates under high-pressure and high-temperature conditions are critical for modeling the dynamics and solidification mechanisms of the magma ocean in the early Earth, as well as for constraining entrainment of melts in the mantle and in the present-day core-mantle boundary. Here we present in situ structural measurements by X-ray diffraction of selected amorphous silicates compressed statically in diamond anvil cells (up to 157 GPa at room temperature) or dynamically by laser-generated shock compression (up to 130 GPa and 6,000 K along the MgSiO3 glass Hugoniot). The X-ray diffraction patterns of silicate glasses and liquids reveal similar characteristics over a wide pressure and temperature range. Beyond the increase in Si coordination observed at 20 GPa, we find no evidence for major structural changes occurring in the silicate melts studied up to pressures and temperatures exceeding Earth's core mantle boundary conditions. This result is supported by molecular dynamics calculations. Our findings reinforce the widely used assumption that the silicate glasses studies are appropriate structural analogs for understanding the atomic arrangement of silicate liquids at these high pressures.

Intraplate volcanism originating from upwelling hydrous mantle transition zone

Most magmatism occurring on Earth is conventionally attributed to passive mantle upwelling at mid-ocean ridges, to slab devolatilization at subduction zones, or to mantle plumes. However, the widespread Cenozoic intraplate volcanism in northeast China^1-3 and the young petit-spot volcanoes^4-7 offshore of the Japan Trench cannot readily be associated with any of these mechanisms. In addition, the mantle beneath these types of volcanism is characterized by zones of anomalously low seismic velocity above and below the transition zone^8-12 (a mantle level located at depths between 410 and 660 kilometres). A comprehensive interpretation of these phenomena is lacking. Here we show that most (or possibly all) of the intraplate and petit-spot volcanism and low-velocity zones around the Japanese subduction zone can be explained by the Cenozoic interaction of the subducting Pacific slab with a hydrous mantle transition zone. Numerical modelling indicates that 0.2 to 0.3 weight per cent of water dissolved in mantle minerals that are driven out from the transition zone in response to subduction and retreat of a tectonic plate is sufficient to reproduce the observations. This suggests that a critical amount of water may have accumulated in the transition zone around this subduction zone, as well as in others of the Tethyan tectonic belt¹³ that are characterized by intraplate or petit-spot volcanism and low-velocity zones in the underlying mantle.

Density of NaAlSi2O6 Melt at High Pressure and Temperature Measured by In-Situ X-ray Microtomography

In this study, the volumetric compression of jadeite (NaAlSi2O6) melt at high pressures was determined by three-dimensional volume imaging using the synchrotron-based X-ray microtomography technique in a rotation-anvil device. Combined with the sample mass, measured using a high-precision analytical balance prior to the high-pressure experiment, the density of jadeite melt was obtained at high pressures and high temperatures up to 4.8 GPa and 1955 K. The density data were fitted to a third-order Birch-Murnaghan equation of state, resulting in a best-fit isothermal bulk modulus K T 0 of 10.8 − 5.3 + 1.9 GPa and its pressure derivative K T 0 ′ of 3.4 − 0.4 + 6.6 . Comparison with data for silicate melts of various compositions from the literature shows that alkali-rich, polymerized melts are generally more compressible than alkali-poor, depolymerized ones. The high compressibility of jadeite melt at high pressures implies that polymerized sodium aluminosilicate melts, if generated by low-degree partial melting of mantle peridotite at ~250–400 km depth in the deep upper mantle, are likely denser than surrounding mantle materials, and thus gravitationally stable.

Formation of bridgmanite-enriched layer at the top lower-mantle during magma ocean solidification

Thermochemical heterogeneities detected today in the Earth’s mantle could arise from ongoing partial melting in different mantle regions. A major open question, however, is the level of chemical stratification inherited from an early magma-ocean (MO) solidification. Here we show that the MO crystallized homogeneously in the deep mantle, but with chemical fractionation at depths around 1000 km and in the upper mantle. Our arguments are based on accurate measurements of the viscosity of melts with forsterite, enstatite and diopside compositions up to ~30 GPa and more than 3000 K at synchrotron X-ray facilities. Fractional solidification would induce the formation of a bridgmanite-enriched layer at ~1000 km depth. This layer may have resisted to mantle mixing by convection and cause the reported viscosity peak and anomalous dynamic impedance. On the other hand, fractional solidification in the upper mantle would have favored the formation of the first crust.

Following the impact of the protoplanet Theia, planet Earth likely transformed into a magma ocean. New high temperature and pressure experiments by Xie et al. suggest that a layer enriched in bridgmanite formed during the magma ocean phase of Earth–remnants of this ancient layer today may be responsible for the viscosity peak between 660 and 1500 km in present solid mantle.

Structure and thermodynamics of aqueous urea solutions from ambient to kilobar pressures: From thermodynamic modeling, experiments, and first principles simulations to an accurate force field description

Molecular simulations based on classical force fields are a powerful method for shedding light on the complex behavior of biomolecules in solution. When cosolutes are present in addition to water and biomolecules, subtle balances of weak intermolecular forces have to be accounted for. This imposes high demands on the quality of the underlying force fields, and therefore force field development for small cosolutes is still an active field. Here, we present the development of a new urea force field from studies of urea solutions at ambient and elevated hydrostatic pressures based on a combination of experimental and theoretical approaches. Experimental densities and solvation shell properties from ab initio molecular dynamics simulations at ambient conditions served as the target properties for the force field optimization. Since urea is present in many marine life forms, elevated hydrostatic pressure was rigorously addressed: densities at high pressure were measured by vibrating tube densitometry up to 500 bar and by X-ray absorption up to 5 kbar. Densities were determined by the perturbed-chain statistical associating fluid theory equation of state. Solvation properties were determined by embedded cluster integral equation theory and ab initio molecular dynamics. Our new force field is able to capture the properties of urea solutions at high pressures without further high-pressure adaption, unlike trimethylamine-N-oxide, for which a high-pressure adaption is necessary.

Density variations of TMAO solutions in the kilobar range: Experiments, PC-SAFT predictions, and molecular dynamics simulations

We present measurements, molecular dynamics (MD) simulations, and predictions using Perturbed-Chain Statistical Associating Fluid Theory (PC-SAFT) of the density of aqueous solutions in a pressure range from 1 bar to 5000 bar, a pressure regime that is highly relevant for both biochemical applications and the fundamental understanding of solvation. The accurate determination of density data of pressurized solutions remains challenging. We determined relative density changes from the variations in X-ray absorption through the sample and developed a new water parameter set for PC-SAFT modeling that is appropriate for high pressure conditions in the kilobar regime. As a showcase, we studied trimethylamine N-oxide (TMAO) solutions and demonstrated that their compressibility decreases with the TMAO content. This result is linked to the stabilizing effect of TMAO on the local H-bond network of water. Experiments and calculations, which represent two independent methods, are in very good agreement and are in accordance with results of force field molecular dynamics simulations of the same systems.

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-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.

Iron diapirs entrain silicates to the core and initiate thermochemical plumes

Segregation of the iron core from rocky silicates is a massive evolutionary event in planetary accretion, yet the process of metal segregation remains obscure, due to obstacles in simulating the extreme physical properties of liquid iron and silicates at finite length scales. We present new experimental results studying gravitational instability of an emulsified liquid gallium layer, initially at rest at the interface between two glucose solutions. Metal settling coats liquid metal drops with a film of low density material. The emulsified metal pond descends as a coherent Rayleigh-Taylor instability with a trailing fluid-filled conduit. Scaling to planetary interiors and high pressure mineral experiments indicates that molten silicates and volatiles are entrained toward the iron core and initiate buoyant thermochemical plumes that later oxidize and hydrate the upper mantle. Surface volcanism from thermochemical plumes releases oxygen and volatiles linking atmospheric growth to the Earth's mantle and core processes.

SiO_{2} Glass Density to Lower-Mantle Pressures

The convection or settling of matter in the deep Earth's interior is mostly constrained by density variations between the different reservoirs. Knowledge of the density contrast between solid and molten silicates is thus of prime importance to understand and model the dynamic behavior of the past and present Earth. SiO_{2} is the main constituent of Earth's mantle and is the reference model system for the behavior of silicate melts at high pressure. Here, we apply our recently developed x-ray absorption technique to the density of SiO_{2} glass up to 110 GPa, doubling the pressure range for such measurements. Our density data validate recent molecular dynamics simulations and are in good agreement with previous experimental studies conducted at lower pressure. Silica glass rapidly densifies up to 40 GPa, but the density trend then flattens to become asymptotic to the density of SiO_{2} minerals above 60 GPa. The density data present two discontinuities at ∼17 and ∼60  GPa that can be related to a silicon coordination increase from 4 to a mixed 5/6 coordination and from 5/6 to sixfold, respectively. SiO_{2} glass becomes denser than MgSiO_{3} glass at ∼40  GPa, and its density becomes identical to that of MgSiO_{3} glass above 80 GPa. Our results on SiO_{2} glass may suggest that a variation of SiO_{2} content in a basaltic or pyrolitic melt with pressure has at most a minor effect on the final melt density, and iron partitioning between the melts and residual solids is the predominant factor that controls melt buoyancy in the lowermost mantle.

Miniature diamond anvils for X-ray Raman scattering spectroscopy experiments at high pressure

X-ray Raman scattering (XRS) spectroscopy is an inelastic scattering method that uses hard X-rays of the order of 10 keV to measure energy-loss spectra at absorption edges of light elements (Si, Mg, O etc.), with an energy resolution below 1 eV. The high-energy X-rays employed with this technique can penetrate thick or dense sample containers such as the diamond anvils employed in high-pressure cells. Here, we describe the use of custom-made conical miniature diamond anvils of less than 500 µm thickness which allow pressure generation of up to 70 GPa. This set-up overcomes the limitations of the XRS technique in very high-pressure measurements (>10 GPa) by drastically improving the signal-to-noise ratio. The conical shape of the base of the diamonds gives a 70° opening angle, enabling measurements in both low- and high-angle scattering geometry. This reduction of the diamond thickness to one-third of the classical diamond anvils considerably lowers the attenuation of the incoming and the scattered beams and thus enhances the signal-to-noise ratio significantly. A further improvement of the signal-to-background ratio is obtained by a recess of ∼20 µm that is milled in the culet of the miniature anvils. This recess increases the sample scattering volume by a factor of three at a pressure of 60 GPa. Examples of X-ray Raman spectra collected at the O K-edge and Si L-edge in SiO₂ glass at high pressures up to 47 GPa demonstrate the significant improvement and potential for spectroscopic studies of low-Z elements at high pressure.

Direct tomography imaging for inelastic X-ray scattering experiments at high pressure

A method to separate the non-resonant inelastic X-ray scattering signal of a micro-metric sample contained inside a diamond anvil cell (DAC) from the signal originating from the high-pressure sample environment is described. Especially for high-pressure experiments, the parasitic signal originating from the diamond anvils, the gasket and/or the pressure medium can easily obscure the sample signal or even render the experiment impossible. Another severe complication for high-pressure non-resonant inelastic X-ray measurements, such as X-ray Raman scattering spectroscopy, can be the proximity of the desired sample edge energy to an absorption edge energy of elements constituting the DAC. It is shown that recording the scattered signal in a spatially resolved manner allows these problems to be overcome by separating the sample signal from the spurious scattering of the DAC without constraints on the solid angle of detection. Furthermore, simple machine learning algorithms facilitate finding the corresponding detector pixels that record the sample signal. The outlined experimental technique and data analysis approach are demonstrated by presenting spectra of the Si L_2,3-edge and O K-edge of compressed α-quartz. The spectra are of unprecedented quality and both the O K-edge and the Si L_2,3-edge clearly show the existence of a pressure-induced phase transition between 10 and 24 GPa.

Solid-liquid density and spin crossovers in (Mg, Fe)O system at deep mantle conditions

The low/ultralow-velocity zones in the Earth's mantle can be explained by the presence of partial melting, critically depending on density contrast between the melt and surrounding solid mantle. Here, first-principles molecular dynamics simulations of (Mg, Fe) O ferropericlase in the solid and liquid states show that their densities increasingly approach each other as pressure increases. The isochemical density difference between them diminishes from 0.78 (±0.7) g/cm3 at zero pressure (3000 K) to 0.16 (±0.04) g/cm3 at 135 GPa (4000 K) for pure and alloyed compositions containing up to 25% iron. The simulations also predict a high-spin to low-spin transition of iron in the liquid ferropericlase gradually occurring over a pressure interval centered at 55 GPa (4000 K) accompanied by a density increase of 0.14 (±0.02) g/cm3. Temperature tends to widen the transition to higher pressure. The estimated iron partition coefficient between the solid and liquid ferropericlase varies from 0.3 to 0.6 over the pressure range of 23 to 135 GPa. Based on these results, an excess of as low as 5% iron dissolved in the liquid could cause the solid-liquid density crossover at conditions of the lowermost mantle.