Multiphase equation of state for carbon addressing high pressures and temperatures
Name
PhysRevB.89.224109.pdf
Size
1.72 MB
Format
Adobe PDF
Checksum (MD5)
65f31836b16d406112c158e476d171fd
Author(s)
Benedict, Lorin X.
Driver, Kevin P.
Hamel, Sebastien
Militzer, Burkhard
Qi, Tingting
Correa, Alfredo A.
Schwegler, Eric
Saul, Alberto Andres
Date Issued
June 2014
Journal
Physical Review B
Publisher
American Physical Society
Citation
Benedict, Lorin X., Kevin P. Driver, Sebastien Hamel, Burkhard Militzer, Tingting Qi, Alfredo A. Correa, A. Saul, and Eric Schwegler. “Multiphase Equation of State for Carbon Addressing High Pressures and Temperatures.” Phys. Rev. B 89, no. 22 (June 2014). © 2014 American Physical Society
Version
Final published version
Abstract
We present a 5-phase equation of state for elemental carbon which addresses a wide range of density and temperature conditions: 3g/cc < ρ < 20g/cc, 0K < T < ∞. The phases considered are diamond, BC8, simple cubic, simple hexagonal, and the liquid/plasma state. The solid phase free energies are constrained by density functional theory (DFT) calculations. Vibrational contributions to the free energy of each solid phase are treated within the quasiharmonic framework. The liquid free energy model is constrained by fitting to a combination of DFT molecular dynamics performed over the range 10 000K < T < 100 000K, and path integral quantum Monte Carlo calculations for T > 100 000K (both for ρ between 3 and 12 g/cc, with select higher-ρ DFT calculations as well). The liquid free energy model includes an atom-in-jellium approach to account for the effects of ionization due to temperature and pressure in the plasma state, and an ion-thermal model which includes the approach to the ideal gas limit. The precise manner in which the ideal gas limit is reached is greatly constrained by both the highest-temperature DFT data and the path integral data, forcing us to discard an ion-thermal model we had used previously in favor of a new one. Predictions are made for the principal Hugoniot and the room-temperature isotherm, and comparisons are made to recent experimental results.
MIT Department
Massachusetts Institute of Technology. Department of Civil and Environmental Engineering
Terms of Use
Article is made available in accordance with the publisher's policy and may be subject to US copyright law. Please refer to the publisher's site for terms of use.
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1103/PhysRevB.89.224109
