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

A Method for Predicting the Melting Temperature of Ionic Compounds

Predicting the melting temperature of materials has always been a topic of great concern. This article proposes an alternative model for determining the melting temperature of materials based on the main idea of the Lindemann melting criterion combined with the first-principles calculations of density functional theory. To verify the accuracy of the melting model, this article selected typical ionic crystals of MgO and 10 alkali metal halides as the validation objects. The calculation results indicate that the melting temperature of the MgO crystals and I-VII compounds is in good agreement with the experimental results.

Stability of hydrides in sub-Neptune exoplanets with thick hydrogen-rich atmospheres

Many sub-Neptune exoplanets have been believed to be composed of a thick hydrogen-dominated atmosphere and a high-temperature heavier-element-dominant core. From an assumption that there is no chemical reaction between hydrogen and silicates/metals at the atmosphere-interior boundary, the cores of sub-Neptunes have been modeled with molten silicates and metals (magma) in previous studies. In large sub-Neptunes, pressure at the atmosphere-magma boundary can reach tens of gigapascals where hydrogen is a dense liquid. A recent experiment showed that hydrogen can induce the reduction of Fe[Formula: see text] in (Mg,Fe)O to Fe[Formula: see text] metal at the pressure-temperature conditions relevant to the atmosphere-interior boundary. However, it is unclear whether Mg, one of the abundant heavy elements in the planetary interiors, remains oxidized or can be reduced by H. Our experiments in the laser-heated diamond-anvil cell found that heating of MgO + Fe to 3,500 to 4,900 K (close to or above their melting temperatures) in an H medium leads to the formation of Mg[Formula: see text]FeH[Formula: see text] and H[Formula: see text]O at 8 to 13 GPa. At 26 to 29 GPa, the behavior of the system changes, and Mg-H in an H fluid and H[Formula: see text]O were detected with separate FeH[Formula: see text]. The observations indicate the dissociation of the Mg-O bond by H and subsequent production of hydride and water. Therefore, the atmosphere-magma interaction can lead to a fundamentally different mineralogy for sub-Neptune exoplanets compared with rocky planets. The change in the chemical reaction at the higher pressures can also affect the size demographics (i.e., "radius cliff") and the atmosphere chemistry of sub-Neptune exoplanets.

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.

Melting and defect transitions in FeO up to pressures of Earth's core-mantle boundary

The high-pressure melting curve of FeO controls key aspects of Earth's deep interior and the evolution of rocky planets more broadly. However, existing melting studies on wüstite were conducted across a limited pressure range and exhibit substantial disagreement. Here we use an in-situ dual-technique approach that combines a suite of >1000 x-ray diffraction and synchrotron Mössbauer measurements to report the melting curve for Fe1-xO wüstite to pressures of Earth's lowermost mantle. We further observe features in the data suggesting an order-disorder transition in the iron defect structure several hundred kelvin below melting. This solid-solid transition, suggested by decades of ambient pressure research, is detected across the full pressure range of the study (30 to 140 GPa). At 136 GPa, our results constrain a relatively high melting temperature of 4140 ± 110 K, which falls above recent temperature estimates for Earth's present-day core-mantle boundary and supports the viability of solid FeO-rich structures at the roots of mantle plumes. The coincidence of the defect order-disorder transition with pressure-temperature conditions of Earth's mantle base raises broad questions about its possible influence on key physical properties of the region, including rheology and conductivity.

Machine-Learning Accelerated First-Principles Accurate Modeling of the Solid-Liquid Phase Transition in MgO under Mantle Conditions

While accurate measurements of MgO under extreme high-pressure conditions are needed to understand and model planetary behavior, these studies are challenging from both experimental and computational modeling perspectives. Herein, we accelerate density functional theory (DFT) accurate calculations using deep neural network potentials (DNPs) trained over multiple phases and study the melting behavior of MgO via the two-phase coexistence (TPC) approach at 0-300 GPa and ≤9600 K. The resulting DNP-TPC melting curve is in excellent agreement with existing experimental studies. We show that the mitigation of finite-size effects that typically skew the predicted melting temperatures in DFT-TPC simulations in excess of several hundred kelvin requires models with ∼16 000 atoms and >100 ps molecular dynamics trajectories. In addition, the DNP can successfully describe MgO metallization well at increased pressures that are captured by DFT but missed by classical interatomic potentials.

Toward better understanding of the high-pressure structural transformation in beryllium by the statistical moment method

Beryllium is a vital alkaline-earth metal for plasma physics, space science, and nuclear technology. Unfortunately, its accurate phase diagram is clouded by many controversial results, even though solid beryllium can only exist with hcp or bcc crystalline structures. Herein, we offer a simple quantum-statistical solution to the above problem. Our core idea is to develop the moment expansion technique to determine the Helmholtz free energy under extreme conditions. This strategy helps elucidate the underlying correlation among symmetric characteristics, vibrational excitations, and physical stabilities. In particular, our analyses reveal that the appearance of anharmonic effects forcefully straightens up the hcp-bcc boundary. This phenomenon explains why it has been difficult to detect bcc signatures via diamond-anvil-cell measurements. Besides, we modify the work-heat equivalence principle to quickly obtain the high-pressure melting profile from the room-temperature equation of state. The hcp-bcc-liquid triple point of beryllium is found at 165 GPa and 4559 K. Our theoretical findings agree excellently with cutting-edge ab initio simulations adopting the phonon quasiparticle method and the thermodynamic integration. Finally, we consider the principal Hugoniot curve and its secondary branches to explore the behaviors of beryllium under shock compression. Our predictions would be advantageous for designing inertial-confinement-fusion experiments.

P–V–T Equation of State of Iridium Up to 80 GPa and 3100 K

In the present study, the high-pressure high-temperature equation of the state of iridium has been determined through a combination of in situ synchrotron X-ray diffraction experiments using laser-heating diamond-anvil cells (up to 48 GPa and 3100 K) and density-functional theory calculations (up to 80 GPa and 3000 K). The melting temperature of iridium at 40 GPa was also determined experimentally as being 4260 (200) K. The results obtained with the two different methods are fully consistent and agree with previous thermal expansion studies performed at ambient pressure. The resulting thermal equation of state can be described using a third-order Birch–Murnaghan formalism with a Berman thermal-expansion model. The present equation of the state of iridium can be used as a reliable primary pressure standard for static experiments up to 80 GPa and 3100 K. A comparison with gold, copper, platinum, niobium, rhenium, tantalum, and osmium is also presented. On top of that, the radial-distribution function of liquid iridium has been determined from experiments and calculations.

A Review of the Melting Curves of Transition Metals at High Pressures Using Static Compression Techniques

The accurate determination of melting curves for transition metals is an intense topic within high pressure research, both because of the technical challenges included as well as the controversial data obtained from various experiments. This review presents the main static techniques that are used for melting studies, with a strong focus on the diamond anvil cell; it also explores the state of the art of melting detection methods and analyzes the major reasons for discrepancies in the determination of the melting curves of transition metals. The physics of the melting transition is also discussed.

Fluid-like elastic response of superionic NH3 in Uranus and Neptune

Nondipolar magnetic fields exhibited at Uranus and Neptune may be derived from a unique geometry of their icy mantle with a thin convective layer on top of a stratified nonconvective layer. The presence of superionic H₂O and NH₃ has been thought as an explanation to stabilize such nonconvective regions. However, a lack of experimental data on the physical properties of those superionic phases has prevented the clarification of this matter. Here, our Brillouin measurements for NH₃ show a two-stage reduction in longitudinal wave velocity (V _p) by ∼9% and ∼20% relative to the molecular solid in the temperature range of 1,500 K and 2,000 K above 47 GPa. While the first V _p reduction observed at the boundary to the superionic α phase was most likely due to the onset of the hydrogen diffusion, the further one was likely attributed to the transition to another superionic phase, denoted γ phase, exhibiting the higher diffusivity. The reduction rate of V _p in the superionic γ phase, comparable to that of the liquid, implies that this phase elastically behaves almost like a liquid. Our measurements show that superionic NH₃ becomes convective and cannot contribute to the internal stratification.

Optical and electronic solutions for power stabilization of CO2 lasers

High pressure-temperature conditions can be readily achieved through the laser-heated diamond anvil cell (LH-DAC). A stable laser source is required for reliable in situ measurements of the sample, as the sample is small with a thermal time constant of the order of microseconds. Here, we show that the power instabilities typical of CO₂ gas lasers used in LH-DAC's are ±5% at the second timescale and ∼±50% at the microsecond timescale. We also demonstrate that the pointing instability of the laser requires either a diffuser or an integrating sphere for reliable total power measurements with small sized detectors. We present a simple solution for stabilizing the power of a CO₂ gas laser on the second timescale by the direct modulation of the current across the tube and another solution that stabilizes the power to the microsecond timescale by externally modulating the CO₂ laser beam. Both solutions can achieve a ±0.3% power stability.

High-pressure/high-temperature phase diagram of zinc

The phase diagram of zinc (Zn) has been explored up to 140 GPa and 6000 K, by combining optical observations, x-ray diffraction, and ab initio calculations. In the pressure range covered by this study, Zn is found to retain a hexagonal close-packed (hcp) crystal symmetry up to the melting temperature. The known decrease of the axial ratio (c/a) of the hcp phase of Zn under compression is observed in x-ray diffraction experiments from 300 K up to the melting temperature. The pressure at which c/a reaches [Formula: see text] (≈10 GPa) is slightly affected by temperature. When this axial ratio is reached, we observed that single crystals of Zn, formed at high temperature, break into multiple poly-crystals. In addition, a noticeable change in the pressure dependence of c/a takes place at the same pressure. Both phenomena could be caused by an isomorphic second-order phase transition induced by pressure in Zn. The reported melt curve extends previous results from 24 to 135 GPa. The pressure dependence obtained for the melting temperature is accurately described up to 135 GPa by using a Simon-Glatzel equation: [Formula: see text], where P is the pressure in GPa. The determined melt curve agrees with previous low-pressure studies and with shock-wave experiments, with a melting temperature of 5060(30) K at 135 GPa. Finally, a thermal equation of state is reported, which at room-temperature agrees with the literature.

High-Pressure Synthesis: A New Frontier in the Search for Next-Generation Intermetallic Compounds

The application of high pressure adds an additional dimension to chemical phase space, opening up an unexplored expanse bearing tremendous potential for discovery. Our continuing mission is to explore this new frontier, to seek out new intermetallic compounds and new solid-state bonding. Simple binary elemental systems, in particular those composed of pairs of elements that do not form compounds under ambient pressures, can yield novel crystalline phases under compression. Thus, high-pressure synthesis can provide access to solid-state compounds that cannot be formed with traditional thermodynamic methods. An emerging approach for the rapid exploration of composition-pressure-temperature phase space is the use of hand-held high-pressure devices known as diamond anvil cells (DACs). These devices were originally developed by geologists as a way to study minerals under conditions relevant to the earth's interior, but they possess a host of capabilities that make them ideal for high-pressure solid-state synthesis. Of particular importance, they offer the capability for in situ spectroscopic and diffraction measurements, thereby enabling continuous reaction monitoring-a powerful capability for solid-state synthesis. In this Account, we provide an overview of this approach in the context of research we have performed in the pursuit of new intermetallic compounds. We start with a discussion of pressure as a fundamental experimental variable that enables the formation of intermetallic compounds that cannot be isolated under ambient conditions. We then introduce the DAC apparatus and explain how it can be repurposed for use as a synthetic vessel with which to explore this phase space, going to extremes of pressure where no chemist has gone before. The remainder of the Account is devoted to discussions of recent experiments we have performed with this approach that have led to the discovery of novel intermetallic compounds in the Fe-Bi, Cu-Bi, and Ni-Bi systems, with a focus on the cutting-edge methods that made these experiments possible. We review the use of in situ laser heating at high pressure, which led to the discovery of FeBi₂, the first binary intermetallic compound in the Fe-Bi system. Our work in the Cu-Bi system is described in the context of in situ experiments carried out in the DAC to map its high-pressure phase space, which revealed two intermetallic phases (Cu₁₁Bi₇ and CuBi). Finally, we review the discovery of β-NiBi, a novel high-pressure phase in the Ni-Bi system. We hope that this Account will inspire the next generation of solid-state chemists to boldly explore high-pressure phase space.

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.

Synthesis of boron-doped diamond and its application as a heating material in a multi-anvil high-pressure apparatus

We developed methods to use synthesized boron-doped diamond (BDD) as a heater in a multi-anvil high-pressure apparatus. The synthesized BDD heater could stably generate an ultra-high temperature without the issues (anomalous melt, pressure drop, and instability of heating) arising from oxidation of boron into boron oxide and graphite-diamond conversion. We synthesized BDD blocks and tubes with boron contents of 0.5-3.0 wt. % from a mixture of graphite and amorphous boron at 15 GPa and 2000 °C. The electrical conductivity of BDD increased with increasing boron content. The stability of the heater and heating reproducibility were confirmed through repeated cycles of heating and cooling. Temperatures as high as ∼3700 °C were successfully generated at higher than 10 GPa using the BDD heater. The effect of the BDD heater on the pressure-generation efficiency was evaluated using MgO pressure scale by in situ X-ray diffraction study at the SPring-8 synchrotron. The pressure-generation efficiency was lower than that using a graphite-boron composite heater up to 1500 tons. The achievement of stable temperature generation above 3000 °C enables melting experiments of silicates and determination of some physical properties (such as viscosity) of silicate melts under the Earth's lower mantle conditions.