skip to main content

DOE PAGESDOE PAGES

This content will become publicly available on June 25, 2019

Title: Path integral Monte Carlo simulations of warm dense aluminum

Here, we perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of $0.1-32.4$~g$$\,$$cm$$^{-3}$$ and $10^4-10^8$~K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multi-shock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multi-shock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and Z to reveal the evolution of the plasma structure and ionization behavior.
Authors:
 [1] ;  [2] ;  [3]
  1. Univ. of California, Berkeley, CA (United States). Dept. of Earth and Planetary Science; Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  2. Univ. of California, Berkeley, CA (United States). Dept. of Earth and Planetary Science; Univ. Claude Bernard Lyon, Lyon (France). Higher Normal School and Geology Lab. of Lyon
  3. Univ. of California, Berkeley, CA (United States). Dept. of Earth and Planetary Science and Dept. of Astronomy
Publication Date:
Report Number(s):
LLNL-JRNL-741093
Journal ID: ISSN 2470-0045; PLEEE8; SC0010517, AC02-05CH11231
Grant/Contract Number:
SC0016248; SC0010517; AC52-07NA27344; AC02-05CH11231
Type:
Accepted Manuscript
Journal Name:
Physical Review E
Additional Journal Information:
Journal Volume: 97; Journal Issue: 6; Related Information: https://journals.aps.org/pre/supplemental/10.1103/PhysRevE.97.063207; Journal ID: ISSN 2470-0045
Publisher:
American Physical Society (APS)
Research Org:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). National Energy Research Scientific Computing Center (NERSC)
Sponsoring Org:
USDOE Office of Science (SC), Fusion Energy Sciences (FES) (SC-24)
Contributing Orgs:
Univ. of Illinois at Urbana-Champaign, IL (United States). National Center for Supercomputing Applications (NCSA). Blue Waters
Country of Publication:
United States
Language:
English
Subject:
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Equation of State; Physics - Plasma physics
OSTI Identifier:
1457418
Alternate Identifier(s):
OSTI ID: 1457086; OSTI ID: 1458619

Driver, K. P., Soubiran, F., and Militzer, B.. Path integral Monte Carlo simulations of warm dense aluminum. United States: N. p., Web. doi:10.1103/PhysRevE.97.063207.
Driver, K. P., Soubiran, F., & Militzer, B.. Path integral Monte Carlo simulations of warm dense aluminum. United States. doi:10.1103/PhysRevE.97.063207.
Driver, K. P., Soubiran, F., and Militzer, B.. 2018. "Path integral Monte Carlo simulations of warm dense aluminum". United States. doi:10.1103/PhysRevE.97.063207.
@article{osti_1457418,
title = {Path integral Monte Carlo simulations of warm dense aluminum},
author = {Driver, K. P. and Soubiran, F. and Militzer, B.},
abstractNote = {Here, we perform first-principles path integral Monte Carlo (PIMC) and density functional theory molecular dynamics (DFT-MD) calculations to explore warm dense matter states of aluminum. Our equation of state (EOS) simulations cover a wide density-temperature range of $0.1-32.4$~g$\,$cm$^{-3}$ and $10^4-10^8$~K. Since PIMC and DFT-MD accurately treat effects of the atomic shell structure, we find two compression maxima along the principal Hugoniot curve attributed to K-shell and L-shell ionization. The results provide a benchmark for widely used EOS tables, such as SESAME, QEOS, and models based on Thomas-Fermi and average-atom techniques. A subsequent multi-shock analysis provides a quantitative assessment for how much heating occurs relative to an isentrope in multi-shock experiments. Finally, we compute heat capacity, pair-correlation functions, the electronic density of states, and Z to reveal the evolution of the plasma structure and ionization behavior.},
doi = {10.1103/PhysRevE.97.063207},
journal = {Physical Review E},
number = 6,
volume = 97,
place = {United States},
year = {2018},
month = {6}
}

Works referenced in this record:

Projector augmented-wave method
journal, December 1994

From ultrasoft pseudopotentials to the projector augmented-wave method
journal, January 1999

Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set
journal, October 1996