DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Benchmarking boron carbide equation of state using computation and experiment

Abstract

Boron carbide ( B 4 C ) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B 4 C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B 4 C over a broad range of temperatures ( ~ 6 × 10 3 5 × 10 8 K) and densities (0.025–50 g / cm 3 ). We assess that the largest discrepancies between theoretical predictions are 5 % near the compression maximum at 1– 2 × 106 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~ 18 % and occur at temperatures between 6 × 10 3 2 × 10 5 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.

Authors:
ORCiD logo [1]; ORCiD logo [2];  [2];  [2];  [3]; ORCiD logo [2];  [2];  [2];  [2];  [2];  [2]; ORCiD logo [2];  [2];  [2];  [2];  [2]
  1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States); Univ. of Rochester, NY (United States). Lab. for Laser Energetics
  2. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  3. Univ. of California, Berkeley, CA (United States)
Publication Date:
Research Org.:
Lawrence Livermore National Laboratory (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE National Nuclear Security Administration (NNSA)
OSTI Identifier:
1775037
Report Number(s):
LLNL-JRNL-812984
Journal ID: ISSN 2470-0045; 1020196; TRN: US2209296
Grant/Contract Number:  
AC52-07NA27344; SC0016248
Resource Type:
Accepted Manuscript
Journal Name:
Physical Review. E
Additional Journal Information:
Journal Volume: 102; Journal Issue: 5; Journal ID: ISSN 2470-0045
Publisher:
American Physical Society (APS)
Country of Publication:
United States
Language:
English
Subject:
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND SUPERFLUIDITY; equations of state; high-energy-density plasmas; plasma ionization; shock waves; hot dense plasma; warm-dense matter; density functional calculations; path-integral Monte Carlo

Citation Formats

Zhang, Shuai, Marshall, Michelle C., Yang, Lin H., Sterne, Philip A., Militzer, Burkhard, Däne, Markus, Gaffney, James A., Shamp, Andrew, Ogitsu, Tadashi, Caspersen, Kyle, Lazicki, Amy E., Erskine, David, London, Richard A., Celliers, Peter M., Nilsen, Joseph, and Whitley, Heather D. Benchmarking boron carbide equation of state using computation and experiment. United States: N. p., 2020. Web. doi:10.1103/physreve.102.053203.
Zhang, Shuai, Marshall, Michelle C., Yang, Lin H., Sterne, Philip A., Militzer, Burkhard, Däne, Markus, Gaffney, James A., Shamp, Andrew, Ogitsu, Tadashi, Caspersen, Kyle, Lazicki, Amy E., Erskine, David, London, Richard A., Celliers, Peter M., Nilsen, Joseph, & Whitley, Heather D. Benchmarking boron carbide equation of state using computation and experiment. United States. https://doi.org/10.1103/physreve.102.053203
Zhang, Shuai, Marshall, Michelle C., Yang, Lin H., Sterne, Philip A., Militzer, Burkhard, Däne, Markus, Gaffney, James A., Shamp, Andrew, Ogitsu, Tadashi, Caspersen, Kyle, Lazicki, Amy E., Erskine, David, London, Richard A., Celliers, Peter M., Nilsen, Joseph, and Whitley, Heather D. Tue . "Benchmarking boron carbide equation of state using computation and experiment". United States. https://doi.org/10.1103/physreve.102.053203. https://www.osti.gov/servlets/purl/1775037.
@article{osti_1775037,
title = {Benchmarking boron carbide equation of state using computation and experiment},
author = {Zhang, Shuai and Marshall, Michelle C. and Yang, Lin H. and Sterne, Philip A. and Militzer, Burkhard and Däne, Markus and Gaffney, James A. and Shamp, Andrew and Ogitsu, Tadashi and Caspersen, Kyle and Lazicki, Amy E. and Erskine, David and London, Richard A. and Celliers, Peter M. and Nilsen, Joseph and Whitley, Heather D.},
abstractNote = {Boron carbide (B4C) is of both fundamental scientific and practical interest due to its structural complexity and how it changes upon compression, as well as its many industrial uses and potential for use in inertial confinement fusion (ICF) and high-energy density physics experiments. Here, we report the results of a comprehensive computational study of the equation of state (EOS) of B4C in the liquid, warm dense matter, and plasma phases. Our calculations are cross-validated by comparisons with Hugoniot measurements up to 61 megabar from planar shock experiments performed at the National Ignition Facility (NIF). Our computational methods include path integral Monte Carlo, activity expansion, as well as all-electron Green's function Korringa-Kohn-Rostoker and molecular dynamics that are both based on density functional theory. We calculate the pressure-internal energy EOS of B4C over a broad range of temperatures (~6×103–5×108 K) and densities (0.025–50 g/cm3). We assess that the largest discrepancies between theoretical predictions are ≲5% near the compression maximum at 1–2×106 K. This is the warm-dense state in which the K shell significantly ionizes and has posed grand challenges to theory and experiment. By comparing with different EOS models, we find a Purgatorio model (LEOS 2122) that agrees with our calculations. The maximum discrepancies in pressure between our first-principles predictions and LEOS 2122 are ~18% and occur at temperatures between 6×103–2×105 K, which we believe originate from differences in the ion thermal term and the cold curve that are modeled in LEOS 2122 in comparison with our first-principles calculations. To account for potential differences in the ion thermal term, we have developed three new equation-of-state models that are consistent with theoretical calculations and experiment. We apply these new models to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.},
doi = {10.1103/physreve.102.053203},
journal = {Physical Review. E},
number = 5,
volume = 102,
place = {United States},
year = {Tue Nov 03 00:00:00 EST 2020},
month = {Tue Nov 03 00:00:00 EST 2020}
}

Works referenced in this record:

Generalized Gradient Approximation Made Simple
journal, October 1996

  • Perdew, John P.; Burke, Kieron; Ernzerhof, Matthias
  • Physical Review Letters, Vol. 77, Issue 18, p. 3865-3868
  • DOI: 10.1103/PhysRevLett.77.3865

Equation of state of boron nitride combining computation, modeling, and experiment
journal, April 2019


Line-imaging velocimeter for shock diagnostics at the OMEGA laser facility
journal, November 2004

  • Celliers, P. M.; Bradley, D. K.; Collins, G. W.
  • Review of Scientific Instruments, Vol. 75, Issue 11
  • DOI: 10.1063/1.1807008

Statistical mechanics of Coulomb gases of arbitrary charge
journal, December 1974


High-precision shock equation of state measurements for metallic fluid carbon between 15 and 20 Mbar
journal, October 2020

  • Millot, Marius; Sterne, Philip A.; Eggert, Jon H.
  • Physics of Plasmas, Vol. 27, Issue 10
  • DOI: 10.1063/5.0007304

High-precision measurements of the equation of state of hydrocarbons at 1–10 Mbar using laser-driven shock waves
journal, May 2010

  • Barrios, M. A.; Hicks, D. G.; Boehly, T. R.
  • Physics of Plasmas, Vol. 17, Issue 5
  • DOI: 10.1063/1.3358144

Raman spectroscopy of pressure-induced amorphous boron carbide
journal, March 2006

  • Yan, X. Q.; Li, W. J.; Goto, T.
  • Applied Physics Letters, Vol. 88, Issue 13
  • DOI: 10.1063/1.2189826

Nanotwins soften boron-rich boron carbide (B 13 C 2 )
journal, March 2017

  • An, Qi; Goddard, William A.
  • Applied Physics Letters, Vol. 110, Issue 11
  • DOI: 10.1063/1.4978644

Path integral Monte Carlo and density functional molecular dynamics simulations of warm dense MgSiO 3
journal, January 2020

  • González-Cataldo, Felipe; Soubiran, François; Peterson, Henry
  • Physical Review B, Vol. 101, Issue 2
  • DOI: 10.1103/PhysRevB.101.024107

First-principles equation of state and electronic properties of warm dense oxygen
journal, October 2015

  • Driver, K. P.; Soubiran, F.; Zhang, Shuai
  • The Journal of Chemical Physics, Vol. 143, Issue 16
  • DOI: 10.1063/1.4934348

High-density carbon ablator experiments on the National Ignition Facility
journal, May 2014

  • MacKinnon, A. J.; Meezan, N. B.; Ross, J. S.
  • Physics of Plasmas, Vol. 21, Issue 5
  • DOI: 10.1063/1.4876611

Preparation of a B4C hollow microsphere through gel-casting for an inertial confinement fusion (ICF) target
journal, January 2017


Large-scale quantum mechanical simulations of high-Z metals
journal, July 2007

  • Yang, L. H.; Hood, Randolph Q.; Pask, J. E.
  • Journal of Computer-Aided Materials Design, Vol. 14, Issue 3
  • DOI: 10.1007/s10820-007-9053-1

Ab initio calculations of boron and its carbides
conference, January 1991

  • Kleinman, Leonard
  • Boron-rich solids, AIP Conference Proceedings
  • DOI: 10.1063/1.40861

Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis
journal, August 1980

  • Vosko, S. H.; Wilk, L.; Nusair, M.
  • Canadian Journal of Physics, Vol. 58, Issue 8
  • DOI: 10.1139/p80-159

First-principles equation of state and shock compression predictions of warm dense hydrocarbons
journal, July 2017


Comparison of plastic, high density carbon, and beryllium as indirect drive NIF ablators
journal, May 2018

  • Kritcher, A. L.; Clark, D.; Haan, S.
  • Physics of Plasmas, Vol. 25, Issue 5
  • DOI: 10.1063/1.5018000

Depressurization Amorphization of Single-Crystal Boron Carbide
journal, February 2009


Optimized norm-conserving Vanderbilt pseudopotentials
journal, August 2013


\mathcal{O}(N) methods in electronic structure calculations
journal, February 2012


Making inertial confinement fusion models more predictive
journal, August 2019

  • Gaffney, Jim A.; Brandon, Scott T.; Humbird, Kelli D.
  • Physics of Plasmas, Vol. 26, Issue 8
  • DOI: 10.1063/1.5108667

Statistical Mechanics of Reacting Coulomb Gases
journal, August 1973


Nanoindentation and Raman spectroscopy studies of boron carbide single crystals
journal, November 2002

  • Domnich, Vladislav; Gogotsi, Yury; Trenary, Michael
  • Applied Physics Letters, Vol. 81, Issue 20
  • DOI: 10.1063/1.1521580

Hugoniot experiments with unsteady waves
journal, July 2014

  • Fratanduono, D. E.; Munro, D. H.; Celliers, P. M.
  • Journal of Applied Physics, Vol. 116, Issue 3
  • DOI: 10.1063/1.4890014

A measurement of the equation of state of carbon envelopes of white dwarfs
journal, August 2020


Equation of state, adiabatic sound speed, and Grüneisen coefficient of boron carbide along the principal Hugoniot to 700 GPa
journal, November 2016


Quantum and classical molecular dynamics simulation of boron carbide behavior under pressure
conference, January 2017

  • Korotaev, P.; Kuksin, A.; Pokatashkin, P.
  • SHOCK COMPRESSION OF CONDENSED MATTER - 2015: Proceedings of the Conference of the American Physical Society Topical Group on Shock Compression of Condensed Matter, AIP Conference Proceedings
  • DOI: 10.1063/1.4971602

Development of Path Integral Monte Carlo Simulations with Localized Nodal Surfaces for Second-Row Elements
journal, October 2015


The National Ignition Facility: Ushering in a new age for high energy density science
journal, April 2009

  • Moses, E. I.; Boyd, R. N.; Remington, B. A.
  • Physics of Plasmas, Vol. 16, Issue 4
  • DOI: 10.1063/1.3116505

Computation of the high temperature Coulomb density matrix in periodic boundary conditions
journal, July 2016


The Crystal Structure of Boron Carbide
journal, November 1943

  • Clark, H. K.; Hoard, J. L.
  • Journal of the American Chemical Society, Vol. 65, Issue 11
  • DOI: 10.1021/ja01251a026

Path integral Monte Carlo simulations of warm dense sodium
journal, December 2016


Three-dimensional simulations of low foot and high foot implosion experiments on the National Ignition Facility
journal, March 2016

  • Clark, D. S.; Weber, C. R.; Milovich, J. L.
  • Physics of Plasmas, Vol. 23, Issue 5
  • DOI: 10.1063/1.4943527

Equations of state and melting curve of boron carbide in the high-pressure range of shock compression
journal, March 2017

  • Molodets, A. M.; Golyshev, A. A.; Shakhrai, D. V.
  • Journal of Experimental and Theoretical Physics, Vol. 124, Issue 3
  • DOI: 10.1134/S1063776117030049

A unified formulation of the constant temperature molecular dynamics methods
journal, July 1984

  • Nosé, Shuichi
  • The Journal of Chemical Physics, Vol. 81, Issue 1
  • DOI: 10.1063/1.447334

Formation of composites in equilibrium plasmas
journal, January 1979


Boron carbide—A comprehensive review
journal, January 1990


Shock-Induced Localized Amorphization in Boron Carbide
journal, March 2003


Structural damage in boron carbide under contact loading
journal, August 2004


Equation of state of dense, partially degenerate, reacting plasmas
journal, September 1981


Absolute Hugoniot measurements from a spherically convergent shock using x-ray radiography
journal, May 2018

  • Swift, Damian C.; Kritcher, Andrea L.; Hawreliak, James A.
  • Review of Scientific Instruments, Vol. 89, Issue 5
  • DOI: 10.1063/1.5032142

Development and modeling of a polar-direct-drive exploding pusher platform at the National Ignition Facility
journal, July 2018

  • Ellison, C. Leland; Whitley, Heather D.; Brown, Colin R. D.
  • Physics of Plasmas, Vol. 25, Issue 7
  • DOI: 10.1063/1.5025724

First-principles equation of state calculations of warm dense nitrogen
journal, February 2016


Purgatorio—a new implementation of the Inferno algorithm
journal, May 2006

  • Wilson, B.; Sonnad, V.; Sterne, P.
  • Journal of Quantitative Spectroscopy and Radiative Transfer, Vol. 99, Issue 1-3
  • DOI: 10.1016/j.jqsrt.2005.05.053

Measuring the shock impedance mismatch between high-density carbon and deuterium at the National Ignition Facility
journal, April 2018


A new quotidian equation of state (QEOS) for hot dense matter
journal, January 1988

  • More, R. M.; Warren, K. H.; Young, D. A.
  • Physics of Fluids, Vol. 31, Issue 10
  • DOI: 10.1063/1.866963

All-Electron Path Integral Monte Carlo Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas
journal, March 2012


A new global equation of state model for hot, dense matter
journal, September 1995

  • Young, David A.; Corey, Ellen M.
  • Journal of Applied Physics, Vol. 78, Issue 6
  • DOI: 10.1063/1.359955

Magnesium oxide at extreme temperatures and pressures studied with first-principles simulations
journal, December 2019

  • Soubiran, François; González-Cataldo, Felipe; Driver, Kevin P.
  • The Journal of Chemical Physics, Vol. 151, Issue 21
  • DOI: 10.1063/1.5126624

Precision equation-of-state measurements on National Ignition Facility ablator materials from 1 to 12 Mbar using laser-driven shock waves
journal, May 2012

  • Barrios, M. A.; Boehly, T. R.; Hicks, D. G.
  • Journal of Applied Physics, Vol. 111, Issue 9
  • DOI: 10.1063/1.4712050

Special points for Brillouin-zone integrations
journal, June 1976

  • Monkhorst, Hendrik J.; Pack, James D.
  • Physical Review B, Vol. 13, Issue 12, p. 5188-5192
  • DOI: 10.1103/PhysRevB.13.5188

Effect of Mixing at the Fuel–Ablator Interface on the Burning of Inertial Confinement Fusion Targets Upon Direct Irradiation with a Megajoule Laser Pulse
journal, March 2017

  • Gus’kov, S. Yu.; Demchenko, N. N.; Zmitrenko, N. V.
  • Journal of Russian Laser Research, Vol. 38, Issue 2
  • DOI: 10.1007/s10946-017-9631-y

Properties of B 4 C in the shocked state for pressures up to 1.5 TPa
journal, May 2017


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


Shock compression behaviors of boron carbide (B[sub 4]C)
journal, January 2006

  • Zhang, Y.; Mashimo, T.; Uemura, Y.
  • Journal of Applied Physics, Vol. 100, Issue 11
  • DOI: 10.1063/1.2399334

Off-Hugoniot characterization of alternative inertial confinement fusion ablator materials.
journal, May 2016


High-pressure deformation and amorphization in boron carbide
journal, June 2019

  • Awasthi, Amnaya P.; Subhash, Ghatu
  • Journal of Applied Physics, Vol. 125, Issue 21
  • DOI: 10.1063/1.5091795

Theoretical phase diagram of boron carbide from ambient to high pressure and temperature
journal, May 2019

  • Jay, Antoine; Hardouin Duparc, Olivier; Sjakste, Jelena
  • Journal of Applied Physics, Vol. 125, Issue 18
  • DOI: 10.1063/1.5091000

Implosion performance of subscale beryllium capsules on the NIF
journal, May 2019

  • Zylstra, A. B.; MacLaren, S.; Yi, S. A.
  • Physics of Plasmas, Vol. 26, Issue 5
  • DOI: 10.1063/1.5098319

An Optimized Method for Treating Long-Range Potentials
journal, March 1995

  • Natoli, Vincent; Ceperley, David M.
  • Journal of Computational Physics, Vol. 117, Issue 1
  • DOI: 10.1006/jcph.1995.1054

Applications of nuclear physics
journal, January 2017


Equations of State for Ablator Materials in Inertial Confinement Fusion Simulations
journal, May 2016


New insights into the enigma of boron carbide inverse molecular behavior
journal, July 2014


Path integrals in the theory of condensed helium
journal, April 1995


Using neutrons to measure keV temperatures in highly compressed plastic at multi-Gbar pressures
journal, December 2016


Absolute Equation-of-State Measurement for Polystyrene from 25 to 60 Mbar Using a Spherically Converging Shock Wave
journal, July 2018


Theoretical and experimental investigation of the equation of state of boron plasmas
journal, August 2018


Hugoniot and release measurements in diamond shocked up to 26 Mbar
journal, April 2017


Dynamic behavior of boron carbide
journal, April 2004

  • Vogler, T. J.; Reinhart, W. D.; Chhabildas, L. C.
  • Journal of Applied Physics, Vol. 95, Issue 8
  • DOI: 10.1063/1.1686902

Equation of state and shock compression of warm dense sodium—A first-principles study
journal, February 2017

  • Zhang, Shuai; Driver, Kevin P.; Soubiran, François
  • The Journal of Chemical Physics, Vol. 146, Issue 7
  • DOI: 10.1063/1.4976559

X-ray diffraction study of B 4 C under high pressure
journal, March 2010


Developing quartz and molybdenum as impedance-matching standards in the 100-Mbar regime
journal, May 2019


Gigabar material properties experiments on nif and omega
conference, January 2012

  • Swift, Damian; Hawreliak, James; Braun, David
  • SHOCK COMPRESSION OF CONDENSED MATTER - 2011: Proceedings of the Conference of the American Physical Society Topical Group on Shock Compression of Condensed Matter, AIP Conference Proceedings
  • DOI: 10.1063/1.3686321

Shock Hugoniot measurements of CH at Gbar pressures at the NIF
journal, March 2016


Hugoniot equation of state and dynamic strength of boron carbide
journal, April 2015

  • Grady, Dennis E.
  • Journal of Applied Physics, Vol. 117, Issue 16
  • DOI: 10.1063/1.4918604

Improved tetrahedron method for Brillouin-zone integrations
journal, June 1994


Equation of state of hot, dense magnesium derived with first-principles computer simulations
journal, September 2020

  • González-Cataldo, Felipe; Soubiran, François; Militzer, Burkhard
  • Physics of Plasmas, Vol. 27, Issue 9
  • DOI: 10.1063/5.0017555

Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas
journal, March 2018

  • Zhang, Shuai; Militzer, Burkhard; Benedict, Lorin X.
  • The Journal of Chemical Physics, Vol. 148, Issue 10
  • DOI: 10.1063/1.5001208

Quantum and classical molecular dynamics simulation of boron carbide behavior under pressure
conference, January 2017

  • Korotaev, P.; Kuksin, A.; Pokatashkin, P.
  • SHOCK COMPRESSION OF CONDENSED MATTER - 2015: Proceedings of the Conference of the American Physical Society Topical Group on Shock Compression of Condensed Matter, AIP Conference Proceedings
  • DOI: 10.1063/1.4971602