Melting of (Mg, Fe)SiO 3 -Perovskite to 625 Kilobars: Indication of a High Melting Temperature in the Lower Mantle

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.

Temperature dependence of nitrogen solubility in bridgmanite and evolution of nitrogen storage capacity in the lower mantle

Relative nitrogen abundance normalized by carbonaceous chondrites in the bulk silicate Earth appears to be depleted compared to other volatile elements. Especially, nitrogen behavior in the deep part of the Earth such as the lower mantle is not clearly understood. Here, we experimentally investigated the temperature dependence of nitrogen solubility in bridgmanite which occupies 75 wt.% of the lower mantle. The experimental temperature ranged from 1400 to 1700 °C at 28 GPa in the redox state corresponding to the shallow lower mantle. The maximum nitrogen solubility in bridgmanite (MgSiO₃) increased from 1.8 ± 0.4 to 5.7 ± 0.8 ppm with increasing temperature from 1400 to 1700 °C. The nitrogen storage capacity of Mg-endmember bridgmanite under the current temperature conditions is 3.4 PAN (PAN: mass of present atmospheric nitrogen). Furthermore, the nitrogen solubility of bridgmanite increased with increasing temperature, in contrast to the nitrogen solubility of metallic iron. Thus, the nitrogen storage capacity of bridgmanite can be larger than that of metallic iron during the solidification of the magma ocean. Such a "hidden" nitrogen reservoir formed by bridgmanite in the lower mantle may have depleted the apparent nitrogen abundance ratio in the bulk silicate Earth.

New frontiers in extreme conditions science at synchrotrons and free electron lasers

Synchrotrons and free electron lasers are unique facilities to probe the atomic structure and electronic properties of matter at extreme thermodynamical conditions. In this context, 'matter at extreme pressures and temperatures' was one of the science drivers for the construction of low emittance 4th generation synchrotron sources such as the Extremely Brilliant Source of the European Synchrotron Radiation Facility and hard x-ray free electron lasers, such as the European x-ray free electron laser. These new user facilities combine static high pressure and dynamic shock compression experiments to outstanding high brilliance and submicron beams. This combination not only increases the data-quality but also enlarges tremendously the accessible pressure, temperature and density space. At the same time, the large spectrum of available complementary x-ray diagnostics for static and shock compression studies opens unprecedented insights into the state of matter at extremes. The article aims at highlighting a new horizon of scientific opportunities based on the synergy between extremely brilliant synchrotrons and hard x-ray free electron lasers.

Melting and density of MgSiO3 determined by shock compression of bridgmanite to 1254GPa

The essential data for interior and thermal evolution models of the Earth and super-Earths are the density and melting of mantle silicate under extreme conditions. Here, we report an unprecedently high melting temperature of MgSiO3 at 500 GPa by direct shockwave loading of pre-synthesized dense MgSiO3 (bridgmanite) using the Z Pulsed Power Facility. We also present the first high-precision density data of crystalline MgSiO3 to 422 GPa and 7200 K and of silicate melt to 1254 GPa. The experimental density measurements support our density functional theory based molecular dynamics calculations, providing benchmarks for theoretical calculations under extreme conditions. The excellent agreement between experiment and theory provides a reliable reference density profile for super-Earth mantles. Furthermore, the observed upper bound of melting temperature, 9430 K at 500 GPa, provides a critical constraint on the accretion energy required to melt the mantle and the prospect of driving a dynamo in massive rocky planets.

A CO2 laser heating system for in situ high pressure-temperature experiments at HPCAT

We present a CO₂ laser heating setup for synchrotron x-ray diffraction inside a diamond anvil cell, situated at HPCAT (Sector 16, Advanced Photon Source, Argonne National Lab, Illinois, USA), which is modular and portable between the HPCAT experiment hutches. The system allows direct laser heating of wide bandgap insulating materials to thousands of degrees at static high pressures up to the Mbar regime. Alignment of the focused CO₂ laser spot is performed using a mid-infrared microscope, which addressed past difficulties with aligning the invisible radiation. The implementation of the mid-infrared microscope alongside a mirror pinhole spatial filter system allows precise alignment of the heating laser spot and optical pyrometry measurement location to the x-ray probe. A comparatively large heating spot (∼50 μm) relative to the x-ray beam (<10 μm) reduces the risk of temperature gradients across the probed area. Each component of the heating system and its diagnostics have been designed with portability in mind and compatibility with the various experimental hutches at the HPCAT beamlines. We present measurements on ZrO₂ at 5.5 GPa which demonstrate the improved room-temperature diffraction data quality afforded by annealing with the CO₂ laser. We also present in situ measurements at 5.5 GPa up to 2800 K in which we do not observe the postulated fluorite ZrO₂ structure, in agreement with recent findings.

Viscosity jump in the lower mantle inferred from melting curves of ferropericlase

Convection provides the mechanism behind plate tectonics, which allows oceanic lithosphere to be subducted into the mantle as "slabs" and new rock to be generated by volcanism. Stagnation of subducting slabs and deflection of rising plumes in Earth's shallow lower mantle have been suggested to result from a viscosity increase at those depths. However, the mechanism for this increase remains elusive. Here, we examine the melting behavior in the MgO-FeO binary system at high pressures using the laser-heated diamond-anvil cell and show that the liquidus and solidus of (Mg _x Fe_1-x )O ferropericlase (x = ~0.52-0.98), exhibit a local maximum at ~40 GPa, likely caused by the spin transition of iron. We calculate the relative viscosity profiles of ferropericlase using homologous temperature scaling and find that viscosity increases 10-100 times from ~750 km to ~1000-1250 km, with a smaller decrease at deeper depths, pointing to a single mechanism for slab stagnation and plume deflection.

Melting temperatures of MgO under high pressure by micro-texture analysis

Periclase (MgO) is the second most abundant mineral after bridgmanite in the Earth's lower mantle, and its melting behaviour under pressure is important to constrain rheological properties and melting behaviours of the lower mantle materials. Significant discrepancies exist between the melting temperatures of MgO determined by laser-heated diamond anvil cell (LHDAC) and those based on dynamic compressions and theoretical predictions. Here we show the melting temperatures in earlier LHDAC experiments are underestimated due to misjudgment of melting, based on micro-texture observations of the quenched samples. The high melting temperatures of MgO suggest that the subducted cold slabs should have higher viscosities than previously thought, suggesting that the inter-connecting textural feature of MgO would not play important roles for the slab stagnation in the lower mantle. The present results also predict that the ultra-deep magmas produced in the lower mantle are peridotitic, which are stabilized near the core–mantle boundary.

Melting behaviour of MgO under pressure remains unclear despite the importance of constraining the rheology and composition of the Earth's mantle. Here, the authors show that melting temperatures in earlier static experiments were underestimated based on micro-texture analysis of the quenched samples.

Modeling the melting of multicomponent systems: the case of MgSiO3 perovskite under lower mantle conditions

Knowledge of the melting properties of materials, especially at extreme pressure conditions, represents a long-standing scientific challenge. For instance, there is currently considerable uncertainty over the melting temperatures of the high-pressure mantle mineral, bridgmanite (MgSiO3-perovskite), with current estimates of the melting T at the base of the mantle ranging from 4800 K to 8000 K. The difficulty with experimentally measuring high pressure melting temperatures has motivated the use of ab initio methods, however, melting is a complex multi-scale phenomenon and the timescale for melting can be prohibitively long. Here we show that a combination of empirical and ab-initio molecular dynamics calculations can be used to successfully predict the melting point of multicomponent systems, such as MgSiO3 perovskite. We predict the correct low-pressure melting T, and at high-pressure we show that the melting temperature is only 5000 K at 120 GPa, a value lower than nearly all previous estimates. In addition, we believe that this strategy is of general applicability and therefore suitable for any system under physical conditions where simpler models fail.

Mapping temperatures and temperature gradients during flash heating in a diamond-anvil cell

Here, we couple two-dimensional, 4-color multi-wavelength imaging radiometry with laser flash heating to determine temperature profiles and melting temperatures under high pressures in a diamond-anvil cell. This technique combines the attributes of flash heating (e.g., minimal chemical reactions, thermal runaway, and sample instability), with those of multi-wavelength imaging radiometry (e.g., 2D temperature mapping and reduction of chromatic aberrations). Using this new technique in conjunction with electron microscopy makes a powerful tool to determine melting temperatures at high pressures generated by a diamond-anvil cell.

Demixing instability in dense molten MgSiO3 and the phase diagram of MgO

The phase diagrams of MgSiO3 and MgO are studied from first-principles theory for pressures and temperatures up to 600 GPa and 20,000 K. Through the evaluation of finite-temperature Gibbs free energies, using density-functional theory within the generalized gradient approximation as well as with hybrid exchange-correlation functionals, we find evidence for a vast pressure-temperature regime where molten MgSiO3 decomposes into liquid SiO2 and solid MgO, with a volume change of approximately 1.2%. The demixing transition is driven by the crystallization of MgO--the reaction only occurs below the high-pressure MgO melting curve. The predicted transition pressure at 10,000 K is in close proximity to an anomaly reported in recent laser-driven shock experiments of MgSiO3. We also present new results for the high-pressure melting curve of MgO and its B1-B2 solid phase transition, with a triple point at 364 GPa and 12,000 K.

Evidence for a phase transition in silicate melt at extreme pressure and temperature conditions

Laser-driven shock compression experiments reveal the presence of a phase transition in MgSiO(3) over the pressure-temperature range 300-400 GPa and 10 000-16 000 K, with a positive Clapeyron slope and a volume change of ∼6.3 (±2.0) percent. The observations are most readily interpreted as an abrupt liquid-liquid transition in a silicate composition representative of terrestrial planetary mantles, implying potentially significant consequences for the thermal-chemical evolution of extrasolar planetary interiors. In addition, the present results extend the Hugoniot equation of state of MgSiO(3) single crystal and glass to 950 GPa.

Melting of peridotite to 140 gigapascals

Interrogating physical processes that occur within the lowermost mantle is a key to understanding Earth's evolution and present-day inner composition. Among such processes, partial melting has been proposed to explain mantle regions with ultralow seismic velocities near the core-mantle boundary, but experimental validation at the appropriate temperature and pressure regimes remains challenging. Using laser-heated diamond anvil cells, we constructed the solidus curve of a natural fertile peridotite between 36 and 140 gigapascals. Melting at core-mantle boundary pressures occurs at 4180 ± 150 kelvin, which is a value that matches estimated mantle geotherms. Molten regions may therefore exist at the base of the present-day mantle. Melting phase relations and element partitioning data also show that these liquids could host many incompatible elements at the base of the mantle.

Mineral phases of the Earth´s mantle

Our knowledge of the structure of the Earth´s interior has been obtained by analysing seismic waves that travel in the Earth, and the reference Earth global models used by geophysicists are essentially seismological. Depth profiles of the seismic waves velocities reveal that the deep Earth is divided in several shells, separated by velocity and density discontinuities. The main discontinuity located at a depth of 2900 km corresponds to the transition between the mantle and the core. The Earth´s mantle can be further divided into the upper mantle and the lower mantle, with a transition zone characterised by two prominent increases in velocities observed at 410- and 660-km depths. This article will be focused on the mineral phases of the Earth´s mantle. The interpretation of seismological models in terms of chemical composition and temperature relies on the knowledge of the nature, structure and elastic properties of the candidate materials. We will describe to what extent recent advances in experimental mineral physics and X-ray diffraction have yielded essential knowledge on the structure and high-pressure high-temperature behaviour of pertinent materials, and major improvements in our understanding of the chemical and mineralogical composition of the Earth´s mantle.

X-ray diffraction measurements of Mo melting to 119 GPa and the high pressure phase diagram

In this paper, we report angle-dispersive X-ray diffraction data of molybdenum melting, measured in a double-sided laser-heated diamond-anvil cell up to a pressure of 119 GPa and temperatures up to 3400 K. The new melting temperatures are in excellent agreement with earlier measurements up to 90 GPa that relied on optical observations of melting and in strong contrast to most theoretical estimates. The X-ray measurements show that the solid melts from the bcc structure throughout the reported pressure range and provide no evidence for a high temperature transition from bcc to a close-packed structure, or to any other crystalline structure. This observation contradicts earlier interpretations of shock data arguing for such a transition. Instead, the values for the Poisson ratios of shock compressed Mo, obtained from the sound speed measurements, and the present X-ray evidence of loss of long-range order suggest that the 210 GPa (approximately 4100 K) transition in the shock experiment is from the bcc structure to a new, highly viscous, structured melt.

High-pressure chemistry of nitride-based materials

Besides temperature at one atmosphere, the applied pressure is another important parameter for influencing and controlling reaction pathways and final reaction products. This is relevant not only for the genesis of natural minerals, but also for synthetic chemical products and technological materials. The present critical review (316 references) highlights recent developments that utilise high pressures and high-temperatures for the synthesis of new materials with unique properties, such as high hardness, or interesting magnetic or optoelectronic features. Novel metal nitrides, oxonitrides as well as the new class of nitride-diazenide compounds, all formed under high-pressure conditions, are highlighted. Pure oxides and carbides are not considered here. Moreover, syntheses under high-pressure conditions require special equipment and preparation techniques, completely different from those used for conventional synthetic approaches at ambient pressure. Therefore, we also summarize the high-pressure techniques used for the synthesis of new materials on a laboratory scale. In particular, our attention is focused on reactive gas pressure devices with pressures between 1.2 and 600 MPa, multi-anvil apparatus at P < 25 GPa and the diamond anvil cell, which allows work at pressures of 100 GPa and higher. For example, some of these techniques have been successfully upgraded to an industrial scale for the synthesis of diamond and cubic boron nitride.

Structure and freezing of MgSiO3 liquid in Earth's lower mantle

First-principles molecular-dynamics simulations show that over the pressure regime of Earth's mantle the mean silicon-oxygen coordination number of magnesium metasilicate liquid changes nearly linearly from 4 to 6. The density contrast between liquid and crystal decreases by a factor of nearly 5 over the mantle pressure regime and is 4% at the core-mantle boundary. The ab initio melting curve, obtained by integration of the Clausius-Clapeyron equation, yields a melting temperature at the core-mantle boundary of 5400 +/- 600 kelvins.

Melting curve of MgO from first-principles simulations

First-principles calculations based on density functional theory, both with the local density approximation (LDA) and with generalized gradient corrections (GGA), have been used to simulate solid and liquid MgO in direct coexistence in the range of pressure 0 < or = p < or = 135 GPa. The calculated LDA zero pressure melting temperature is T(LDA)m = 3110 +/- 50 K, in good agreement with the experimental data. The GGA zero pressure melting temperature T(GGA)m = 2575 +/- 100 K is significantly lower than the LDA one, but the difference between the GGA and the LDA is greatly reduced at high pressure. The LDA zero pressure melting slope is dT/dp approximately 100 K/GPa, which is more than 3 times higher than the currently available experimental one from Zerr and Boehler [Nature (London) 371, 506 (1994)]. At the core mantle boundary pressure of 135 GPa MgO melts at Tm = 8140 +/- 150 K.

High-pressure melting of MgSiO3

The melting curve of MgSiO(3) perovskite has been determined by means of ab initio molecular dynamics complemented by effective pair potentials, and a new phenomenological model of melting. Using first principles ground state calculations, we find that the MgSiO(3) perovskite phase transforms into post perovskite at pressures above 100 GPa, in agreement with recent theoretical and experimental studies. We find that the melting curve of MgSiO(3), being very steep at pressures below 60 GPa, rapidly flattens on increasing pressure. The experimental controversy on the melting of the MgSiO(3) perovskite at high pressures is resolved, confirming the data by Zerr and Boehler.

Iron partitioning in Earth's mantle: toward a deep lower mantle discontinuity

We measured the spin state of iron in ferropericlase (Mg0.83Fe0.17)O at high pressure and found a high-spin to low-spin transition occurring in the 60- to 70-gigapascal pressure range, corresponding to depths of 2000 kilometers in Earth's lower mantle. This transition implies that the partition coefficient of iron between ferropericlase and magnesium silicate perovskite, the two main constituents of the lower mantle, may increase by several orders of magnitude, depleting the perovskite phase of its iron. The lower mantle may then be composed of two different layers. The upper layer would consist of a phase mixture with about equal partitioning of iron between magnesium silicate perovskite and ferropericlase, whereas the lower layer would consist of almost iron-free perovskite and iron-rich ferropericlase. This stratification is likely to have profound implications for the transport properties of Earth's lowermost mantle.

Solidus of Earth's deep mantle

The solidus of a pyrolite-like composition, approximating that of the lower mantle, was measured up to 59 gigapascals by using CO2 laser heating in a diamond anvil cell. The solidus temperatures are at least 700 kelvin below the melting temperatures of magnesiowustite, which in the deep mantle has the lowest melting temperatures of the three major components-magnesiowustite, Mg-Si-perovskite, and Ca-Si-perovskite. The solidus in the deep mantle is more than 1500 kelvin above the average present-day geotherm, but at the core-mantle boundary it is near the core temperature. Thus, partial melting of the mantle is possible at the core-mantle boundary.

Multivariable dependence of Fe-Mg partitioning in the lower mantle

High-pressure diamond-cell experiments indicate that the iron-magnesium partitioning between (Fe,Mg)SiO3-perovskite and magnesiowustite in Earth's lower mantle depends on the pressure, temperature, bulk iron/magnesium ratio, and ferric iron content. The perovskite stability field expands with increasing pressure and temperature. The ferric iron component preferentially dissolves in perovskite and raises the apparent total iron content but had little effect on the partitioning of the ferrous iron. The ferrous iron depletes in perovskite at the top of the lower mantle and gradually increases at greater depth. These changes in iron-magnesium composition should affect geochemical and geophysical properties of the deep interior.

Melting of (Mg,Fe)2SiO4 at the Core-Mantle Boundary of the Earth

The lower mantle of the Earth is believed to be largely composed of (Mg,Fe)O (magnesiowustite) and (Mg,Fe)SiO3 (perovskite). Radiative temperatures of single-crystal olivine [(Mg0.9,Fe0.1)2SiO4] decreased abruptly from 7040 +/- 315 to 4300 +/- 270 kelvin upon shock compression above 80 gigapascals. The data indicate that an upper bound to the solidus of the magnesiowustite and perovskite assemblage at 4300 +/- 270 kelvin is 130 +/- 3 gigapascals. These conditions correspond to those for partial melting at the base of the mantle, as has been suggested occurs within the ultralow-velocity zone beneath the central Pacific.

Discussion comment on melting criteria and imaging spectroradiometry in laser-heated diamond-cell experiments (by R. Jeanloz & A. Kavner)

Jeanloz & Kavner have a very timely contribution which raises some very important issues concerning the measurement of temperature using spectroradiometry. In this discussion of the paper, we intend to show that while all the issues raised by the authors should be of real concern to all workers, some of the issues have indeed been resolved with improved technique. However, we agree with the authorn that there are several issues that must be resolved and caution is necessary in interpreting the observations.

Perovskite Temperature Profile

In the News article "Finding ;sustainable' ways to prevent parasitic diseases" by Rebecca Kolberg (24 June, p. 1859), mention is made of nylon gauze that is being distributed worldwide free of charge to filter Guinea worm larvae from drinking water. The cloth is donated by E. I. Dupont de Nemours & Company, Inc., and Precision Fabrics Group for distribution by Global 2000, Inc.

Implications for Mantle Dynamics from the High Melting Temperature of Perovskite

Recent studies have implied that (Mg, Fe)SiO(3)-perovskite, a likely dominant mineral phase in the lower mantle, may have a high melting temperature. The implications of these findings for the dynamics of the lower mantle were investigated with the use of numerical convection models. The results showed that low homologous temperatures (0.3 to 0.5) would prevail in the modeled lower mantle, regardless of the effective Rayleigh number and internal heating rates. High-temperature ductile creep is possible under relatively cold conditions. In models with low rates of internal heating, local maxima of viscosity developed in the mid-lower mantle that were similar to those obtained from inversion of geoid, topography, and plate velocities.

Temperatures in Earth's Core Based on Melting and Phase Transformation Experiments on Iron

Experiments on melting and phase transformations on iron in a laser-heated, diamond-anvil cell to a pressure of 150 gigapascals (approximately 1.5 million atmospheres) show that iron melts at the central core pressure of 363.85 gigapascals at 6350 +/- 350 kelvin. The central core temperature corresponding to the upper temperature of iron melting is 6150 kelvin. The pressure dependence of iron melting temperature is such that a simple model can be used to explain the inner solid core and the outer liquid core. The inner core is nearly isothermal (6150 kelvin at the center to 6130 kelvin at the inner core-outer core boundary), is made of hexagonal closest-packed iron, and is about 1 percent solid (MgSiO(3) + MgO). By the inclusion of less than 2 percent of solid impurities with iron, the outer core densities along a thermal gradient (6130 kelvin at the base of the outer core and 4000 kelvin at the top) can be matched with the average seismic densities of the core.