We report a theoretical equation of state (EOS) table for boron across a wide range of temperatures (5.1×10^{4}-5.2×10^{8} K) and densities (0.25-49 g/cm^{3}) and experimental shock Hugoniot data at ...unprecedented high pressures (5608±118 GPa). The calculations are performed with first-principles methods combining path-integral Monte Carlo (PIMC) at high temperatures and density-functional-theory molecular-dynamics (DFT-MD) methods at lower temperatures. PIMC and DFT-MD cross-validate each other by providing coherent EOS (difference <1.5 Hartree/boron in energy and <5% in pressure) at 5.1×10^{5} K. The Hugoniot measurement is conducted at the National Ignition Facility using a planar shock platform. The pressure-density relation found in our shock experiment is on top of the shock Hugoniot profile predicted with our first-principles EOS and a semiempirical EOS table (LEOS 50). We investigate the self-diffusivity and the effect of thermal and pressure-driven ionization on the EOS and shock compression behavior in high-pressure and -temperature conditions. We also study the sensitivity of a polar direct-drive exploding pusher platform to pressure variations based on applying pressure multipliers to LEOS 50 and by utilizing a new EOS model based on our ab initio simulations via one-dimensional radiation-hydrodynamic calculations. The results are valuable for future theoretical and experimental studies and engineering design in high-energy density research.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
Over the last six years many experiments have been done at the National Ignition Facility
to measure the Hugoniot of materials, such as CH plastic at extreme pressures, up to 800
Mbar. The “Gbar” ...design employs a strong spherically converging shock launched through a
solid ball of material using a hohlraum radiation drive. The shock front conditions are
characterized using x-ray radiography. In this paper we examine the role of radiation in
heating the unshocked material in front of the shock to understand the impact it has on
equation of state measurements and how it drives the measured data off the theoretical
Hugoniot curve. In particular, the two main sources of radiation heating are the
preheating of the unshocked material by the high-energy kilo-electron-volt x-rays in the
hohlraum and the heating of the material in front of the shock, as the shocked material
becomes hot enough to radiate significantly. Using our model, we estimate that preheating
can reach 4 eV in unshocked material, and that radiation heating can begin to drive data
off the Hugoniot significantly, as pressures reach above 400 Mb.
Using first-principles molecular dynamics, we calculated the equation of state and shock Hugoniot of various boron phases. We find a large mismatch between Hugoniots based on existing knowledge of ...the equilibrium phase diagram and those measured by shock experiments, which could be reconciled if the α-B12/β phases are significantly over-pressurized in boron under shock compression. Our results also indicate that there exist an anomaly and negative Clapeyron slope along the melting curve of boron at 100 GPa and 1500–3000 K. These results enable an in-depth understanding of matter under shock compression, in particular the significance of compression-rate dependence of phase transitions and kinetic effects in experimental measurements.
Using first-principles molecular dynamics, we calculated the equation of state and shock Hugoniot of various boron phases. We find a large mismatch between Hugoniots based on existing knowledge of the equilibrium phase diagram and those measured by shock experiments, which could be reconciled if the α-B12/β phases are significantly over-pressurized in boron under shock compression. Our results also indicate that there exist an anomaly and negative Clapeyron slope along the melting curve of boron at 100 GPa and 1500–3000 K. These results enable an in-depth understanding of matter under shock compression, in particular the significance of compression-rate dependence of phase transitions and kinetic effects in experimental measurements. Display omitted
•Using boron as an example, we show that materials phase diagrams are different depending on the compression technique.•We find an anomaly and a negative Clapeyron slope in the melting curve of boron.•We provide first-principles equations of state and shock Hugoniot curves for boron in multiple phases at up to 1000 GPa.•We predict the possible existence of new boron phases at 100 GPa and 2000 K, and of post-α-Ga phases at megabar pressures.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
X-ray absorption spectroscopy and ab initio modeling of the experimental spectra have been used to investigate the effects of surface passivation on the unoccupied electronic states of CdSe quantum ...dots (QDs). Significant differences are observed in the unoccupied electronic structure of the CdSe QDs, which are shown to arise from variations in specific ligand-surface bonding interactions.
Full text
Available for:
IJS, KILJ, NUK, PNG, UL, UM
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. 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×10^{6} 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.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
Warm dense carbon is generated at 0.3–2.0 g/cc and 1–7 eV by proton heating. The release equation of state (EOS) after heating and thermal conductivity of warm dense carbon are studied experimentally ...in this regime using a Au/C dual-layer target to initiate a temperature gradient and two picosecond time-resolved diagnostics to probe the surface expansion and heat flow. Comparison between the data and simulations using various EOSs and thermal conductivity models is quantified with a statistical χ2 analysis. In conclusion, out of seven EOS tables and five thermal conductivity models, only L9061 with the Lee-More model provides a probability above 50% to match all data.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
We study the problem of electron-ion temperature equilibration in plasmas. We consider pure H at various densities and temperatures and Ar-doped H at temperatures high enough so that the Ar is fully ...ionized. Two theoretical approaches are used: classical molecular dynamics (MD) with statistical two-body potentials and a generalized Lenard-Balescu (GLB) theory capable of treating multicomponent weakly coupled plasmas. The GLB is used in two modes: (1) with the quantum dielectric response in the random-phase approximation (RPA) together with the pure Coulomb interaction and (2) with the classical (ℏ→0) dielectric response (both with and without local-field corrections) together with the statistical potentials. We find that the MD results are described very well by classical GLB including the statistical potentials and without local-field corrections (RPA only); worse agreement is found when static local-field effects are included, in contradiction to the classical pure-Coulomb case with like charges. The results of the various approaches are all in excellent agreement with pure-Coulomb quantum GLB when the temperature is high enough. In addition, we show that classical calculations with statistical potentials derived from the exact quantum two-body density matrix produce results in far better agreement with pure-Coulomb quantum GLB than classical calculations performed with older existing statistical potentials.
Full text
Available for:
CMK, CTK, FMFMET, IJS, NUK, PNG, UM
The ground and excited state properties of small helium clusters, 4HeN, containing nanoscale (approx 3-10 A) planar aromatic molecules have been studied with quantum Monte Carlo methods. Ground state ...structures and energies are obtained from importance-sampled, rigid-body diffusion Monte Carlo. Excited state energies due to helium vibrational motion are evaluated using the projection operator, imaginary time spectral evolution technique. We examine the adsorption of N helium atoms (N < or = 24) on a series of planar aromatic molecules (benzene, naphthalene, anthracene, tetracene, phthalo-cyanine). The first layer of helium atoms is well-localized on the molecule surface, and we find well-defined localized excitations due to in-plane vibrational motion of helium on the molecule surface. We discuss the implications of these confined excitations for the molecule spectroscopy.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ