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 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.
more »
« less
Melting Curve of Potassium Chloride from in situ Ionic Conduction Measurements
We report experimental constraints on the melting curve of potassium chloride (KCl) between 3.2 and 9 GPa from in situ ionic conduction measurements using a multi-anvil apparatus. On the basis of concurrent measurements of KCl and sodium chloride (NaCl) at 1 bar using the differential thermal analysis (DTA) method and Pt sphere marker, we show that the peak rate of increase in ionic current with temperature upon heating coincides with latent heat ledge and fall of Pt sphere, thus establishing the criterion for melting detection from ionic conduction measurements. Applying this criterion to high pressures, we found that the melting point of KCl rose steeply with increasing pressure to exceed 2443 ± 100 K at 9 GPa. Fitting the results of this study together with existing data at pressures below 4 GPa and above 20 GPa, we obtained the Simon’s melting equation for KCl in the simple cubic B2 structure between 1.8 and 50 GPa: T m = 1323 ( P − 1.87 2.2 ( 1 ) + 1 ) 1 2.7 ( 1 ) , where T is in K and P is in GPa. Starting at 1 bar, the melting point of KCl increases at an average rate of ~150 K/GPa to cross that of Pt near 9 GPa. The highly refractory nature of KCl makes it a sensitive pressure calibrant for the large-volume pressure at moderate pressures and a potential sample container for experiments at moderate pressures and very high temperatures.
more »
« less
- Award ID(s):
- 1763189
- NSF-PAR ID:
- 10142546
- Date Published:
- Journal Name:
- Minerals
- Volume:
- 10
- Issue:
- 3
- ISSN:
- 2075-163X
- Page Range / eLocation ID:
- 250
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Experimental techniques, such as cryo-electron microscopy, require biological samples to be recovered at cryogenic temperatures ( T ≈ 100 K) with water being in an amorphous ice state. However, (bulk) water can exist in two amorphous ices at P < 1 GPa, low-density amorphous (LDA) ice at low pressures and high-density amorphous ice (HDA) at high pressures; HDA is ≈20–25% denser than LDA. While fast/plunge cooling at 1 bar brings the sample into LDA, high-pressure cooling (HPC), at sufficiently high pressure, produces HDA. HDA can also be produced by isothermal compression of LDA at cryogenic temperatures. Here, we perform classical molecular dynamics simulations to study the effects of LDA, HDA, and the LDA–HDA transformation on the structure and hydration of a small peptide, polyalanine. We follow thermodynamic paths corresponding to (i) fast/plunge cooling at 1 bar, (ii) HPC at P = 400 MPa, and (iii) compression/decompression cycles at T = 80 K. While process (i) produced LDA in the system, path (iii) produces HDA. Interestingly, the amorphous ice produced in process (ii) is an intermediate amorphous ice (IA) with properties that fall in-between those of LDA and HDA. Remarkably, the structural changes in polyalanine are negligible at all conditions studied (0–2000 MPa, 80–300 K) even when water changes among the low and high-density liquid states as well as the amorphous solids LDA, IA, and HDA. The similarities and differences in the hydration of polyalanine vitrified in LDA, IA, and HDA are described. Since the studied thermodynamic paths are suitable for the cryopreservation of biomolecules, we also study the structure and hydration of polyalanine along isobaric and isochoric heating paths, which can be followed experimentally for the recovery of cryopreserved samples. Upon heating, the structure of polyalanine remains practically unchanged. We conclude with a brief discussion of the practical advantages of (a) using HDA and IA as a cryoprotectant environment (as opposed to LDA), and (b) the use of isochoric heating as a recovery process (as opposed to isobaric heating).more » « less
-
Understanding the viscosity of mantle-derived magmas is needed to model their migration mechanisms and ascent rate from the source rock to the surface. High pressure–temperature experimental data are now available on the viscosity of synthetic melts, pure carbonatitic to carbonate–silicate compositions, anhydrous basalts, dacites and rhyolites. However, the viscosity of volatile-bearing melilititic melts, among the most plausible carriers of deep carbon, has not been investigated. In this study, we experimentally determined the viscosity of synthetic liquids with ~31 and ~39 wt% SiO2, 1.60 and 1.42 wt% CO2 and 5.7 and 1 wt% H2O, respectively, at pressures from 1 to 4.7 GPa and temperatures between 1265 and 1755 °C, using the falling-sphere technique combined with in situ X-ray radiography. Our results show viscosities between 0.1044 and 2.1221 Pa·s, with a clear dependence on temperature and SiO2 content. The atomic structure of both melt compositions was also determined at high pressure and temperature, using in situ multi-angle energy-dispersive X-ray diffraction supported by ex situ microFTIR and microRaman spectroscopic measurements. Our results yield evidence that the T–T and T–O (T = Si,Al) interatomic distances of ultrabasic melts are higher than those for basaltic melts known from similar recent studies. Based on our experimental data, melilititic melts are expected to migrate at a rate ~from 2 to 57 km·yr−1 in the present-day or the Archaean mantle, respectively.more » « less
-
more » « less
Abstract The goal of classifying shock metamorphic features in meteorites is to estimate the corresponding shock pressure conditions. However, the temperature variability of shock metamorphism is equally important and can result in a diverse and heterogeneous set of shock features in samples with a common overall shock pressure. In particular, high-pressure (HP) minerals, which were previously used as a solid indicator of high shock pressure in meteorites, require complex pressure–temperature–time (
P–T–t ) histories to form and survive. First, parts of the sample must be heated to melting temperatures, at high pressure, to enable rapid formation of HP minerals before pressure release. Second, the HP minerals must be rapidly cooled to below a critical temperature, before the pressure returns to ambient conditions, to avoid retrograde transformation to their low-pressure polymorphs. These two constraints require the sample to contain large temperature heterogeneities, e.g. melt veins in a cooler groundmass, during shock. In this study, we calculated shock temperatures and possibleP–T paths of chondritic and differentiated mafic–ultramafic rocks for various shock pressures. TheseP–T conditions and paths, combined with observations from shocked meteorites, are used to constrain shock conditions andP–T –t histories of HP-mineral bearing samples. The need for rapid thermal quench of HP phases requires a relatively low bulk-shock temperature and therefore moderate shock pressures below ~ 30 GPa, which matches the stabilities of these HP minerals. The low-temperature moderate-pressure host rock generally shows moderate shock-deformation features consistent with S4 and, less commonly, S5 shock stages. Shock pressures in excess of 50 GPa in meteorites result in melt breccias with high overall post-shock temperatures that anneal out HP-mineral signatures. The presence of ringwoodite, which is commonly considered an indicator of the S6 shock stage, is inconsistent with pressures in excess of 30 GPa and does not represent shock conditions different from S4 shock conditions. Indeed, ringwoodite and coexisting HP minerals should be considered as robust evidence for moderate shock pressures (S4) rather than extreme shock (S6) near whole-rock melting. -
more » « less
Abstract CeOs4Sb12, a member of the skutterudite family, has an unusual semimetallic low-temperature
H —temperatureT phase diagram. We have conducted measurements of electrical transport and megahertz conductivity on CeOs4Sb12single crystals under pressures of up to 3 GPa and in high magnetic fields of up to 41 T to investigate the influence of pressure on the differentH –T phase boundaries. While the high-temperature valence transition between the metallicT by pressures of the order of 1 GPa, we observed only a marginal suppression of the
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.3390/min10030250
