skip to main content


Title: Thermal Pressure in the Thermal Equation of State for Solid and a Proposed Substitute
Abstract

The thermal equation of state (TEOS) for solids is a mathematic model among pressure, temperature and density, and is essential for geophysical, geochemical, and other high pressure–temperature (high P–T) researches. However, in the last few decades, there has been a growing concern about the accuracy of the pressure scales of the calibrants, and efforts have been made to improve it by either introducing a reference standard or building new thermal pressure models. The existing thermal equation of state,P(V,T) = P(V,T0) + Pth(V,T), consists of an isothermal compression and an isochoric heating, while the thermal pressure is the pressure change in the isochoric heating. In this paper, we demonstrate that, for solids in a soft pressure medium in a diamond anvil cell, the thermal pressure can neither be determined from a single heating process, nor from the thermal pressure of its calibrant. To avoid the thermal pressure, we propose to replace the thermal pressure with a well-known thermal expansion model, and integrate it with the isothermal compression model to yields a Birch–Murnaghan-expansion TEOS model, called VPT TEOS. The predicted pressure of MgO and Au at ambient pressure from Birch–Murnaghan-expansion VPT TEOS model matches the experimental pressure of zero (0) GPa very well, while the pressure prediction from the approximated Anderson PVT TEOS exhibit a big deviation and a wrong trend.

 
more » « less
Award ID(s):
2216805
NSF-PAR ID:
10372512
Author(s) / Creator(s):
;
Publisher / Repository:
Springer Science + Business Media
Date Published:
Journal Name:
International Journal of Thermophysics
Volume:
43
Issue:
11
ISSN:
0195-928X
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Thermoelastic behavior of transition metal boride Os2B3 was studied under quasi-hydrostatic and isothermal conditions in a Paris-Edinburgh cell to 5.4 GPa and 1273 K. In-situ Energy Dispersive X-ray diffraction was used to determine interplanar spacings of the hexagonal crystal structure and thus the volume and axial compression. P-V-T data were fitted to a 3rd Order Birch-Murnaghan equation of state with a temperature modification to determine thermal elastic constants. The bulk modulus was shown to be K0 = 402 ± 21 GPa when the first pressure derivative was held to K0’ = 4.0 from the room temperature P-V curve. Under a quadratic fit α=α_0+α_1 T-α_2 T^(-2), the thermal expansion coefficients were determined to be α_0=1.862×10^(-5) K-1, α_1=0.841×10^(-9) K-2, and α_2=-0.525 K. Density functional theory (DFT) with the quasi-harmonic approximation (QHA) were further employed to study Os2B3, including its P-V-T curves, phonon spectra, bulk modulus, specific heat, thermal expansion, and the Grüneisen parameter. A good agreement between the first-principle theory and experimental observations was achieved, highlighting the success of the Armiento-Mattsson 2005 generalized gradient approximation functional employed in this study and QHA for describing thermodynamic properties of Os2B3. 
    more » « less
  2. We report the P-V-T equation of state measurements of B4C to 50GPa and approximately 2500K in laser-heated diamond anvil cells. We obtain an ambient temperature, third-order Birch–Murnaghan fit to the P-V data that yields a bulk modulus K0 of 221(2) GPa and derivative, (dK/dP)0 of 3.3(1). These were used in fits with both a Mie–Grüneisen– Debye model and a temperature-dependent, Birch– Murnaghan equation of state that includes thermal pressure estimated by thermal expansion (α) and a temperature-dependent bulk modulus (dK0/dT). The ambient pressure thermal expansion coefficient (α0+α1T), Grüneisen γ (V)=γ 0(V/V0)q and volumedependent Debye temperature, were used as input parameters for these fits and found to be sufficient to describe the data in the whole P-T range of this study. 
    more » « less
  3. null (Ed.)
    Abstract Calcium carbonate (CaCO3) is one of the most abundant carbonates on Earth's surface and transports carbon to Earth's interior via subduction. Although some petrological observations support the preservation of CaCO3 in cold slabs to lower mantle depths, the geophysical properties and stability of CaCO3 at these depths are not known, due in part to complicated polymorphic phase transitions and lack of constraints on thermodynamic properties. Here we measured thermal equation of state of CaCO3-Pmmn, the stable polymorph of CaCO3 through much of the lower mantle, using synchrotron X-ray diffraction in a laser-heated diamond-anvil cell up to 75 GPa and 2200 K. The room-temperature compression data for CaCO3-Pmmn are fit with third-order Birch-Murnaghan equation of state, yielding KT0 = 146.7 (±1.9) GPa and K′0 = 3.4(±0.1) with V0 fixed to the value determined by ab initio calculation, 97.76 Å3. High-temperature compression data are consistent with zero-pressure thermal expansion αT = a0 + a1T with a0 = 4.3(±0.3)×10-5 K-1, a1 = 0.8(±0.2)×10-8 K-2, temperature derivative of the bulk modulus (∂KT/∂T)P = –0.021(±0.001) GPa/K; the Grüneisen parameter γ0 = 1.94(±0.02), and the volume independent constant q = 1.9(±0.3) at a fixed Debye temperature θ0 = 631 K predicted via ab initio calculation. Using these newly determined thermodynamic parameters, the density and bulk sound velocity of CaCO3-Pmmn and (Ca,Mg)-carbonate-bearing eclogite are quantitatively modeled from 30 to 80 GPa along a cold slab geotherm. With the assumption that carbonates are homogeneously mixed into the slab, the results indicate the presence of carbonates in the subducted slab is unlikely to be detected by seismic observations, and the buoyancy provided by carbonates has a negligible effect on slab dynamics. 
    more » « less
  4. Abstract

    The crystal structure and bonding environment of K2Ca(CO3)2bütschliite were probed under isothermal compression via Raman spectroscopy to 95 GPa and single crystal and powder X-ray diffraction to 12 and 68 GPa, respectively. A second order Birch-Murnaghan equation of state fit to the X-ray data yields a bulk modulus,$${K}_{0}=46.9$$K0=46.9GPa with an imposed value of$${K}_{0}^{\prime}= 4$$K0=4for the ambient pressure phase. Compression of bütschliite is highly anisotropic, with contraction along thec-axis accounting for most of the volume change. Bütschliite undergoes a phase transition to a monoclinicC2/mstructure at around 6 GPa, mirroring polymorphism within isostructural borates. A fit to the compression data of the monoclinic phase yields$${V}_{0}=322.2$$V0=322.2 Å3$$,$$,$${K}_{0}=24.8$$K0=24.8GPa and$${K}_{0}^{\prime}=4.0$$K0=4.0using a third order fit; the ability to access different compression mechanisms gives rise to a more compressible material than the low-pressure phase. In particular, compression of theC2/mphase involves interlayer displacement and twisting of the [CO3] units, and an increase in coordination number of the K+ion. Three more phase transitions, at ~ 28, 34, and 37 GPa occur based on the Raman spectra and powder diffraction data: these give rise to new [CO3] bonding environments within the structure.

     
    more » « less
  5. Abstract

    Synchrotron‐based high‐pressure/high‐temperature single‐crystal X‐ray diffraction experiments to ~24 GPa and 700 K were conducted on eclogitic garnets (low‐Fe: Prp28Alm38Grs33Sps1and high‐Fe: Prp14Alm62Grs19Adr3Sps2) and omphacites (low‐Fe: Quad57Jd42Ae1and high‐Fe: Quad53Jd27Ae20), using an externally heated diamond anvil cell. Fitting the pressure‐volume‐temperature data to a third‐order Birch‐Murnaghan equation of state yields the thermoelastic parameters including bulk modulus (KT0), its pressure derivative (KT0), temperature derivative ((∂KT/∂T)P), and thermal expansion coefficient (αT). The densities of the high‐Fe and low‐Fe eclogites were then modeled along typical geotherms of the normal mantle and the subducted oceanic crust to the transition zone depth (550 km). The metastable low‐Fe eclogite could be a reason for the stagnant slabs within the upper range of the transition zone. Eclogite would be responsible for density anomalies within 100–200 km in the upper mantle of Asia.

     
    more » « less