Benchmarking boron carbide equation of state using computation and experiment.

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.

Theoretical and experimental investigation of the equation of state of boron plasmas.

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.

Self-diffusivity and interdiffusivity of molten aluminum-copper alloys under pressure, derived from molecular dynamics.

We use molecular dynamics (MD) to simulate diffusion in molten aluminum-copper (AlCu) alloys. The self-diffusivities and Maxwell-Stefan diffusivities are calculated for AlCu mixtures using the Green-Kubo formulas at temperatures from 1000 to 4000 K and pressures from 0 to 25 GPa, along with additional points at higher temperatures and pressures. The diffusivities are corrected for finite-size effects. The Maxwell-Stefan diffusivity is compared to the diffusivity calculated from the self-diffusivities using a generalization of the Darken equation. We find that the effects of cross-correlation are small. Using the calculated self-diffusivities, we have assessed whether dilute hard-sphere and dilute Lennard-Jones models apply to the molten mixture. Neither of the two dilute gas diffusivities describes the diffusivity in molten Al and Cu. We report generalized analytic models for the self-diffusivities and interdiffusivity (mutual diffusivity) that fit the MD results well. The MD-derived transport coefficients are in good agreement with the available experimental data. We also report MD calculations of the viscosity and an analytic fit to those results. The ionic thermal conductivity is discussed briefly.

Phase separation in hydrogen-helium mixtures at Mbar pressures.

The properties of hydrogen-helium mixtures at Mbar pressures and intermediate temperatures (4000 to 10000 K) are calculated with first-principles molecular dynamics simulations. We determine the equation of state as a function of density, temperature, and composition and, using thermodynamic integration, we estimate the Gibbs free energy of mixing, thereby determining the temperature, at a given pressure, when helium becomes insoluble in dense metallic hydrogen. These results are directly relevant to models of the interior structure and evolution of Jovian planets. We find that the temperatures for the demixing of helium and hydrogen are sufficiently high to cross the planetary adiabat of Saturn at pressures approximately 5 Mbar; helium is partially miscible throughout a significant portion of the interior of Saturn, and to a lesser extent in Jupiter.

Energetics and kinetics of vacancy diffusion and aggregation in shocked aluminium via orbital-free density functional theory.

A possible mechanism for shock-induced failure in aluminium involves atomic vacancies diffusing through the crystal lattice and agglomerating to form voids, which continue to grow, ultimately resulting in ductile fracture. We employ orbital-free density functional theory, a linear-scaling first-principles quantum mechanics method, to study vacancy formation, diffusion, and aggregation in aluminium under shock loading conditions of compression and tension. We calculate vacancy formation and migration energies, and find that while nearest-neighbor vacancy pairs are unstable, next-nearest-neighbor vacancy pairs are stable. As the number of nearby vacancies increases, we predict that vacancy clusters preferentially grow through next-nearest-neighbor vacancies. The energetics are found to be greatly affected by expansion and compression, leading to insight as to how vacancies behave under shock conditions.

Finding transition states for crystalline solid-solid phase transformations.

We present a method to identify transition states and minimum energy paths for martensitic solid-solid phase transformations, thereby allowing quantification of the activation energies of such transformations. Our approach is a generalization of a previous method for identifying transition states for chemical reactions, namely the climbing image-nudged elastic band algorithm, where here the global deformation of the crystalline lattice (volume and shape fluctuations) becomes the reaction coordinate instead of atomic motion. We also introduce an analogue to the Born-Oppenheimer approximation that allows a decoupling of nuclear motion and lattice deformation, where the nuclear positions along the path are determined variationally according to current deformation state. We then apply this technique to characterize the energetics of elemental lithium phase transformations as a function of applied pressure, where we see a validation of the Born-Oppenheimer-like approximation, small energy barriers (expected for martensitic transformations), and a pronounced pressure dependence of various properties characterizing the phase transitions.

Importance of shear in the bcc-to-hcp transformation in iron.

Iron shows a pressure-induced martensitic phase transformation from the ground state ferromagnetic bcc phase to a nonmagnetic hcp phase at approximately 13 GPa. The exact transformation pressure (TP) and pathway are not known. Here we present a multiscale model containing a quantum-mechanics-based multiwell energy function accounting for the bcc and hcp phases of Fe and a construction of kinematically compatible and equilibrated mixed phases. This model suggests that shear stresses have a significant influence on the bcc<-->hcp transformation. In particular, the presence of modest shear accounts for the scatter in measured TPs. The formation of mixed phases also provides an explanation for the observed hysteresis in TP.