Equation of state of warm-dense boron nitride combining computation, modeling, and experiment

Abstract

The equation of state (EOS) of materials at warm dense conditions poses significant challenges to both theory and experiment. We report a combined computational, modeling, and experimental investigation leveraging new theoretical and experimental capabilities to investigate warm-dense boron nitride (BN). The simulation methodologies include path integral Monte Carlo (PIMC), several density functional theory (DFT) molecular dynamics methods [plane-wave pseudopotential, Fermi operator expansion (FOE), and spectral quadrature (SQ)], activity expansion (ACTEX), and all-electron Green's function Korringa-Kohn-Rostoker (MECCA), and compute the pressure and internal energy of BN over a broad range of densities ($\rho$) and temperatures ($T$). Our experiments were conducted at the Omega laser facility and measured the Hugoniot of BN to unprecedented pressures (12--30 Mbar). The EOSs computed using different methods cross validate one another, and the experimental Hugoniot are in good agreement with our theoretical predictions. We assess that the largest discrepancies between theoretical predictions are $<$4% in pressure and $<$3% in energy and occur at $10^6$ K. We find remarkable consistency between the EOS from DFT calculations performed on different platforms and using different exchange-correlation functionals and those from PIMC using free-particle nodes. This provides strong evidence for the accuracy of both PIMC and DFT in the warm-dense regime. Moreover, SQ and FOE data have significantly smaller error bars than PIMC, and so represent significant advances for efficient computation at high $T$. We also construct tabular EOS models and clarify the ionic and electronic structure of BN over a broad $T-\rho$ range and quantify their roles in the EOS. The tabular models may be utilized for future simulations of laser-driven experiments that include BN as a candidate ablator material.

arXiv comment19 pages, 14 figures, 4 tables

References in corpus (34)

• Optimized norm-conserving Vanderbilt pseudopotentials
• Libxc: a library of exchange and correlation functionals for density functional theory
• Nearsightedness of Electronic Matter
• O(N) methods in electronic structure calculations
• Path Integral Monte Carlo Calculation of the Deuterium Hugoniot
• Molecular Dissociation in Hot, Dense Hydrogen
• FPEOS: A First-Principles Equation of State Table of Deuterium for Inertial Confinement Fusion Applications
• Anisotropic thermal transport in bulk hexagonal boron nitride
• A multiphase equation of state for carbon addressing high pressures and temperatures
• The Decay Properties of the Finite Temperature Density Matrix in Metals
• A real-time approach to the optical properties of solids and nano-structures: the time-dependent Bethe-Salpeter equation
• All-Electron Path Integral Monte Carlo Simulations of Warm Dense Matter: Application to Water and Carbon Plasmas
• Path Integral Monte Carlo Simulation of the Low-Density Hydrogen Plasma
• First Principles Calculations of Shock Compressed Fluid Helium
• Calculation of a Deuterium Double Shock Hugoniot from Ab initio Simulations
• Equation of state of cubic boron nitride at high pressures and temperatures
• Analysis of laser shock experiments on precompressed samples using a quartz reference and application to warm dense hydrogen and helium
• First-Principles Equation of State Calculations of Warm Dense Nitrogen
• Accurate optical properties from first principles: a Quasiparticle Self consistent GW plus Bethe-Salpeter Equation approach
• SQDFT: Spectral Quadrature method for large-scale parallel $\mathcal{O}(N)$ Kohn-Sham calculations at high temperature
• Quantum Monte Carlo calculations in solids with downfolded Hamiltonians
• First-Principles Equation of State and Electronic Properties of Warm Dense Oxygen
• First-principles Equation of State and Shock Compression Predictions of Warm Dense Hydrocarbons
• Simulations of Dense Atomic Hydrogen in the Wigner Crystal Phase
• First-Principles Prediction of the Softening of the Silicon Shock Hugoniot Curve
• Path integral Monte Carlo simulations of dense carbon-hydrogen plasmas
• Fundamental high pressure calibration from all-electron quantum Monte Carlo calculations
• Characterization of the State of Hydrogen
• Equation of state and shock compression of warm dense sodium - a first-principles study
• First-Principles Simulations of Warm Dense Lithium Fluoride
• Hydrogen-Helium Mixtures at High Pressure
• Theoretical and experimental investigation of the equation of state of boron plasmas
• On nearsightedness in metallic systems for $\mathcal{O}(N)$ Density Functional Theory calculations: A case study on Aluminum
• Effect of nanostructuration on compressibility of cubic BN

Cited by in corpus (14)

• First-Principles Equation of State Database for Warm Dense Matter Computation
• Roadmap on Electronic Structure Codes in the Exascale Era
• Path Integral Monte Carlo and Density Functional Molecular Dynamics Simulations of Warm, Dense MgSiO${}_3$
• Properties of carbon up to 10 million kelvin from Kohn-Sham density functional theory molecular dynamics
• Magnesium Oxide at Extreme Temperatures and Pressures Studied with First-Principles Simulations
• Real-space formulation of the stress tensor for $\mathcal{O}(N)$ density functional theory: application to high temperature calculations
• Benchmarking boron carbide equation of state using computation and experiment
• Thermal and Pressure Ionization in Warm, Dense MgSiO${}_3$ Studied with First-Principles Computer Simulations
• Comparison of ablators for the polar direct drive exploding pusher platform
• Phase transformation in boron under shock compression
• Nature of the bonded-to-atomic transition in liquid silica to TPa pressures
• Real-space density kernel method for Kohn-Sham density functional theory calculations at high temperature
• Simultaneous compression and opacity data from time-series radiography with a Lagrangian marker
• First-principles Equations of State and Structures of Liquid Metals in Multi-megabar Conditions