Melting and melt structure of MgO at high pressures

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.

Quantum critical point and spin fluctuations in lower-mantle ferropericlase

Ferropericlase [(Mg,Fe)O] is one of the most abundant minerals of the earth's lower mantle. The high-spin (HS) to low-spin (LS) transition in the Fe(2+) ions may dramatically alter the physical and chemical properties of (Mg,Fe)O in the deep mantle. To understand the effects of compression on the ground electronic state of iron, electronic and magnetic states of Fe(2+) in (Mg0.75Fe0.25)O have been investigated using transmission and synchrotron Mössbauer spectroscopy at high pressures and low temperatures (down to 5 K). Our results show that the ground electronic state of Fe(2+) at the critical pressure Pc of the spin transition close to T = 0 is governed by a quantum critical point (T = 0, P = P(c)) at which the energy required for the fluctuation between HS and LS states is zero. Analysis of the data gives P(c) = 55 GPa. Thermal excitation within the HS or LS states (T > 0 K) is expected to strongly influence the magnetic as well as physical properties of ferropericlase. Multielectron theoretical calculations show that the existence of the quantum critical point at temperatures approaching zero affects not only physical properties of ferropericlase at low temperatures but also its properties at P-T of the earth's lower mantle.

Rocky core solubility in Jupiter and giant exoplanets

Gas giants are believed to form by the accretion of hydrogen-helium gas around an initial protocore of rock and ice. The question of whether the rocky parts of the core dissolve into the fluid H-He layers following formation has significant implications for planetary structure and evolution. Here we use ab initio calculations to study rock solubility in fluid hydrogen, choosing MgO as a representative example of planetary rocky materials, and find MgO to be highly soluble in H for temperatures in excess of approximately 10,000 K, implying the potential for significant redistribution of rocky core material in Jupiter and larger exoplanets.

Infrared absorption of MgO at high pressures and temperatures: a molecular dynamic study

We calculate by molecular dynamics the optical functions of MgO in the far infrared region 100-1000 cm(-1), for pressures up to 40 GPa and temperatures up to 4000 K. An ab initio parametrized many-body force field is used to generate the trajectories. Infrared spectra are obtained from the time correlation of the polarization, and from Kramers-Kronig relations. The calculated spectra agree well with experimental data at ambient pressure. We find that the infrared absorption of MgO at CO(2) laser frequencies increases substantially with both pressure and temperature and we argue that this may explain the underestimation, with respect to theoretical calculations, of the high-pressure melting temperature of MgO determined in CO(2) laser-heated diamond-anvil cell experiments.

Effect of strain on the thermal conductivity of solids

We present a systematic study of the effect of strain (equivalent to uniform pressure) on the thermal conductivity of an insulating solid. Following a theoretical analysis that uncovers the dependence of the thermal conductivity on temperature and strain, we present classical molecular dynamics calculations of the thermal conductivity. We find that the molecular dynamics results closely match the theoretical result.

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.

New insights into the melting behavior of MgO from molecular dynamics simulations: the importance of premelting effects

We provide a plausible resolution of a long-standing controversy relevant to the geophysics community, namely, that the experimental slope of the melting curve Tm(P) of MgO at low pressures is about 3 times smaller than that obtained from computer simulation of the melting of the normal rock-salt-structured crystal. With increasing temperature at zero pressure, our simulations predict a solid-solid phase transition (from a rock salt to a wurtzite crystalline lattice) to occur just before melting. The coexistence of wurtzite and liquid phases at low pressures is found to be described by a Clapeyron slope which is in much better agreement with the experimental results of Zerr and Boehler [Nature (London) 371, 506 (1994)] than the calculated melting line for the rock salt structure. We also show that the existence of a certain concentration of lattice defects in the rock salt phase cannot provide an alternative explanation.