skip to main content

This content will become publicly available on February 1, 2023

Title: The phase diagrams of beryllium and magnesium oxide at megabar pressures
Abstract We perform ab initio simulations of beryllium (Be) and magnesium oxide (MgO) at megabar pressures and compare their structural and thermodynamic properties. We make a detailed comparison of our two recently derived phase diagrams of Be (Wu et al 2021 Phys. Rev. B 104 014103) and MgO (Soubiran and Militzer 2020 Phys. Rev. Lett. 125 175701) using the thermodynamic integration technique, as they exhibit striking similarities regarding their shape. We explore whether the Lindemann criterion can explain the melting temperatures of these materials through the calculation of the Debye temperature at high pressure. From our free energy calculations, we find that the melting line of both materials is well represented by the Simon–Glazel fit T m ( P ) = T 0 (1 + P / a ) 1/ c , where T 0 = 1564 K, a = 15.8037 GPa and c = 2.4154 for Be, while T 0 = 3010 K, a = 10.5797 GPa and c = 2.8683 for the MgO in the B1. For the B2 phase, we use the values a = 26.1163 GPa and c = 2.2426. Both materials exhibit negative Clapeyron slopes on the boundaries between the two solid phases that are strongly more » affected by anharmonic effects, which also influence the location of the solid–solid–liquid triple point. We find that the quasi-harmonic approximation underestimates the stability range of the low-pressure phases, namely hcp for Be and B1 for MgO. We also compute the phonon dispersion relations at low and high pressure for each of the phases of these materials, and also explore how the phonon density of states is modified by temperature. Finally, we derive secondary shock Hugoniot curves in addition to the principal Hugoniot curve for both materials, and study their offsets in pressure between solid and liquid branches. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Journal of Physics: Condensed Matter
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Rocksalt structure nitrides emerge as a promising class of semiconductors for high-temperature thermoelectric and plasmonic applications. Controlling the bandgap and strain is essential for the development of a wide variety of electronic devices. Here we use (Ti 0.5 Mg 0.5 ) 1−x Al x N as a model system to explore and demonstrate the tunability of both the bandgap and the strain state in rocksalt structure nitrides, employing a combined experimental and computational approach. (Ti 0.5 Mg 0.5 ) 1−x Al x N layers with x ≤ 0.44 deposited on MgO(001) substrates by reactive co-sputtering at 700 °C are epitaxial single crystals with a solid-solution B1 rocksalt structure. The lattice mismatch with the substrate decreases with increasing x , leading to a transition in the strain-state from partially relaxed (74% and 38% for x = 0 and 0.09) to fully strained for x ≥ 0.22. First-principles calculations employing 64-atom Special Quasirandom Structures (SQS) indicate that the lattice constant decreases linearly with x according to a = (4.308 − 0.234 x ) Å for 0 ≤ x ≤ 1. In contrast, the measured relaxed lattice parameter a o = (4.269 − 0.131 x ) Å is linear only for x ≤more »0.33, its composition dependence is less pronounced, and x > 0.44 leads to the nucleation of secondary phases. The fundamental (indirect) bandgap predicted using the same SQS supercells and the HSE06 functional increases from 1.0 to 2.6 eV for x = 0–0.75. In contrast, the onset of the measured optical absorption due to interband transitions increases only from 2.3 to 2.6 eV for x = 0–0.44, suggesting that the addition of Al in the solid solution relaxes the electron momentum conservation and causes a shift from direct to indirect gap transitions. The resistivity increases from 9.0 to 708 μΩ m at 77 K and from 6.8 to 89 μΩ m at 295 K with increasing x = 0–0.44, indicating an increasing carrier localization associated with a randomization of cation site occupation and the increasing bandgap which also causes a 33% reduction in the optical carrier concentration. The overall results demonstrate bandgap and strain engineering in rocksalt nitride semiconductors and show that, in contrast to conventional covalent semiconductors, the random cation site occupation strongly affects optical transitions.« less
  2. The mineral apatite, Ca10(PO4)6(F,OH,Cl)2, incorporates sulfur (S) during crystallization from S-bearing hydrothermal fluids and silicate melts. Our previous studies of natural and experimental apatite demonstrate that the oxidation state of S in apatite varies systematically as a function of oxygen fugacity (fO2). The S oxidation states –1 and –2 were quantitatively identified in apatite crystallized from reduced, S-bearing hydrothermal fluids and silicate melts by using sulfur K-edge X‑ray absorption near-edge structure spectroscopy (S-XANES) where S6+/ΣS in apatite increases from ~0 at FMQ-1 to ~1 at FMQ+2, where FMQ refers to the fayalite-magnetite-quartz fO2 buffer. In this study, we employ quantum-mechanical calculations to investigate the atomistic structure and energetics of S(-I) and S(-II) incorporated into apatite and elucidate incorporation mechanisms. One S(-I) species (disulfide, S22−) and two S(-II) species (bisulfide, HS−, and sulfide, S2−) are investigated as possible forms of reduced S species in apatite. In configuration models for the simulation, these reduced S species are positioned along the c-axis channel, originally occupied by the column anions F, Cl, and OH in the end-member apatites. In the lowest-energy configurations of S-incorporated apatite, disulfide prefers to be positioned halfway between the mirror planes at z = 1/4 and 3/4. In contrast, themore »energy-optimized bisulfide is located slightly away from the mirror planes by ~0.04 fractional units in the c direction. The energetic stability of these reduced S species as a function of position along the c-axis can be explained by the geometric and electrostatic constraints of the Ca and O planes that constitute the c-axis channel. The thermodynamics of incorporation of disulfide and bisulfide into apatite are evaluated by using solid-state reaction equations where the apatite host and a solid S-bearing source phase (pyrite and Na2S2(s) for disulfide; troilite and Na2S(s) for sulfide) are the reactants, and the S-incorporated apatite and an anion sink phase are the products. The Gibbs free energy (ΔG) is lower for incorporation with Na-bearing phases than with Fe-bearing phases, which is attributed to the higher energetic stability of the iron sulfide minerals as a source phase for S than the sodium sulfide phases. The thermodynamics of incorporation of reduced S are also evaluated by using reaction equations involving dissolved disulfide and sulfide species [HnS2(aq)(2–n) and HnS(aq)(2–n); n = 0, 1, and 2] as a source phase. The ΔG of S-incorporation increases for fluorapatite and chlorapatite and decreases for hydroxylapatite as these species are protonated (i.e., as n changes from 0 to 2). These thermodynamic results demonstrate that the presence of reduced S in apatite is primarily controlled by the chemistry of magmatic and hydrothermal systems where apatite forms (e.g., an abundance of Fe; solution pH). Ultimately, our methodology developed for evaluating the thermodynamics of S incorporation in apatite as a function of temperature, pH, and composition is highly applicable to predicting the trace and volatile element incorporation in minerals in a variety of geological systems. In addition to solid-solid and solid-liquid equilibria treated here at different temperatures and pH, the methodology can be easily extended also to different pressure conditions by just performing the quantum-mechanical calculations at elevated pressures.« less
  3. Abstract The properties of all materials at one atmosphere of pressure are controlled by the configurations of their valence electrons. At extreme pressures, neighboring atoms approach so close that core-electron orbitals overlap, and theory predicts the emergence of unusual quantum behavior. We ramp-compress monovalent elemental sodium, a prototypical metal at ambient conditions, to nearly 500 GPa (5 million atmospheres). The 7-fold increase of density brings the interatomic distance to 1.74 Å well within the initial 2.03 Å of the Na + ionic diameter, and squeezes the valence electrons into the interstitial voids suggesting the formation of an electride phase. The laser-driven compression results in pressure-driven melting and recrystallization in a billionth of a second. In situ x-ray diffraction reveals a series of unexpected phase transitions upon recrystallization, and optical reflectivity measurements show a precipitous decrease throughout the liquid and solid phases, where the liquid is predicted to have electronic localization. These data reveal the presence of a rich, temperature-driven polymorphism where core electron overlap is thought to stabilize the formation of peculiar electride states.
  4. The degree to which the Earth’s mantle stores and cycles water in excess of the storage capacity of nominally anhydrous minerals is dependent upon the stability of hydrous phases under mantle-relevant pressures, temperatures, and compositions. Two hydrous phases, phase D and phase H, are stable to the pressures and temperatures of the Earth’s lower mantle, suggesting that the Earth’s lower mantle may participate in the cycling of water. We build on our prior work of density functional theory calculations on phase H with the stability, structure, and bonding of hydrous phases D, and we predict the aluminum partitioning with H in the Al 2 O 3 -SiO 2 -MgO-H 2 O system. We address the solid solutions through a statistical sampling of site occupancy and calculation of the partition function from the grand canonical ensemble. We show that each phase has a wide solid solution series between MgSi 2 O 6 H 2 -Al 2 SiO 6 H 2 and MgSiO 4 H 2 -2 δ AlOOH + SiO 2 , in which phase H is more aluminum rich than phase D at a given bulk composition. We predict that the addition of Al to both phases D and Hmore »stabilizes each phase to higher temperatures through additional configurational entropy. While we have shown that phase H does not exhibit symmetric hydrogen bonding at high pressure, we report here that phase D undergoes a gradual increase in the number of symmetric H-bonds beginning at ∼30 GPa, and it is only ∼50% complete at 60 GPa.« less
  5. Abstract Studies of molecular mixtures containing hydrogen sulfide (H 2 S) could open up new routes towards hydrogen-rich high-temperature superconductors under pressure. H 2 S and ammonia (NH 3 ) form hydrogen-bonded molecular mixtures at ambient conditions, but their phase behavior and propensity towards mixing under pressure is not well understood. Here, we show stable phases in the H 2 S–NH 3 system under extreme pressure conditions to 4 Mbar from first-principles crystal structure prediction methods. We identify four stable compositions, two of which, (H 2 S) (NH 3 ) and (H 2 S) (NH 3 ) 4 , are stable in a sequence of structures to the Mbar regime. A re-entrant stabilization of (H 2 S) (NH 3 ) 4 above 300 GPa is driven by a marked reversal of sulfur-hydrogen chemistry. Several stable phases exhibit metallic character. Electron–phonon coupling calculations predict superconducting temperatures up to 50 K, in the Cmma phase of (H 2 S) (NH 3 ) at 150 GPa. The present findings shed light on how sulfur hydride bonding and superconductivity are affected in molecular mixtures. They also suggest a reservoir for hydrogen sulfide in the upper mantle regions of icy planets in a potentially metallic mixture, which couldmore »have implications for their magnetic field formation.« less