Properties of warm dense polystyrene plasmas along the principal Hugoniot

Measuring the principal Hugoniot of inertial-confinement-fusion-relevant TMPTA plastic foams

Wetted-foam layers are of significant interest for inertial-confinement-fusion capsules, due to the control they provide over the convergence ratio of the implosion and the opportunity this affords to minimize hydrodynamic instability growth. However, the equation of state for fusion-relevant foams are not well characterized, and many simulations rely on modeling such foams as a homogeneous medium with the foam average density. To address this issue, an experiment was performed using the VULCAN Nd:glass laser at the Central Laser Facility. The aim was to measure the principal Hugoniot of TMPTA plastic foams at 260mg/cm^{3}, corresponding to the density of liquid DT-wetted-foam layers, and their "hydrodynamic equivalent" capsules. A VISAR was used to obtain the shock velocity of both the foam and an α-quartz reference layer, while streaked optical pyrometry provided the temperature of the shocked material. The measurements confirm that, for the 20-120 GPa pressure range accessed, this material can indeed be well described using the equation of state of the homogeneous medium at the foam density.

First-principles equation of state of CHON resin for inertial confinement fusion applications

A wide-range (0 to 1044.0 g/cm^{3} and 0 to 10^{9} K) equation-of-state (EOS) table for a CH_{1.72}O_{0.37}N_{0.086} quaternary compound has been constructed based on density-functional theory (DFT) molecular-dynamics (MD) calculations using a combination of Kohn-Sham DFT MD, orbital-free DFT MD, and numerical extrapolation. The first-principles EOS data are compared with predictions of simple models, including the fully ionized ideal gas and the Fermi-degenerate electron gas models, to chart their temperature-density conditions of applicability. The shock Hugoniot, thermodynamic properties, and bulk sound velocities are predicted based on the EOS table and compared to those of C-H compounds. The Hugoniot results show the maximum compression ratio of the C-H-O-N resin is larger than that of CH polystyrene due to the existence of oxygen and nitrogen; while the other properties are similar between CHON and CH. Radiation hydrodynamic simulations have been performed using the table for inertial confinement fusion targets with a CHON ablator and compared with a similar design with CH. The simulations show CHON outperforms CH as the ablator for laser-direct-drive target designs.

Species Separation and Hydrogen Streaming upon Shock Release from Polystyrene under Inertial Confinement Fusion Conditions

Shock release from inertial confinement fusion (ICF) shells poses a great challenge to single-fluid hydrodynamic equations, especially for describing materials composed of different ion species. This has been evidenced by a recent experiment [Haberberger et al., Phys. Rev. Lett. 123, 235001 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.235001], in which low-density plasmas (10^{19} to 10^{20}  cm^{-3}) are measured to move far ahead of what radiation-hydrodynamic simulations predict. To understand such experimental observations, we have performed large-scale nonequilibrium molecular-dynamics simulations of shock release in polystyrene (CH) at experimental conditions. These simulations revealed that upon shock releasing from the back surface of a CH foil, hydrogen can stream out of the bulk of the foil due to its mass being lighter than carbon. This released hydrogen, exhibiting a much broader velocity distribution than carbon, forms low-density plasmas moving in nearly constant velocities ahead of the in-flight shell, which is in quantitative agreement with the experimental measurements. Such kinetic effect of species separation is currently missing in single-fluid radiation-hydrodynamics codes for ICF simulations.

Direct-drive double-shell implosion: A platform for burning-plasma physics studies

Double-shell ignition designs have been studied with the indirect-drive inertial confinement fusion (ICF) scheme in both simulations and experiments in which the inner-shell kinetic energy was limited to ∼10-15 kJ, even driven by megajoule-class lasers such as the National Ignition Facility. Since direct-drive ICF can couple more energy to the imploding shells, we have performed a detailed study on direct-drive double-shell (D^{3}S) implosions with state-of-the-art physics models implemented in radiation-hydrodynamic codes (lilac and draco), including nonlocal thermal transport, cross-beam energy transfer (CBET), and first-principles-based material properties. To mitigate classical unstable interfaces, we have proposed the use of a tungsten-beryllium-mixed inner shell with gradient-density layers that can be made by magnetron sputtering. In our D^{3}S designs, a 70-μm-thick beryllium outer shell is driven symmetrically by a high-adiabat (α≥10), 1.9-MJ laser pulse to a peak velocity of ∼240 km/s. Upon spherical impact, the outer shell transfers ∼30-40 kJ of kinetic energy to the inner shell filled with deuterium-tritium gas or liquid, giving neutron-yield energies of ∼6 MJ in one-dimensional simulations. Two-dimensional high-mode draco simulations indicated that such high-adiabat D^{3}S implosions are not susceptible to laser imprint, but the long-wavelength perturbations from the laser port configuration along with CBET can be detrimental to the target performance. Nevertheless, neutron yields of ∼0.3-1.0-MJ energies can still be obtained from our high-mode draco simulations. The robust α-particle bootstrap is readily reached, which could provide a viable platform for burning-plasma physics studies. Once CBET mitigation and/or more laser energy becomes available, we anticipate that break-even or moderate energy gain might be feasible with the proposed D^{3}S scheme.

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

We have developed an experimental platform for the National Ignition Facility that uses spherically converging shock waves for absolute equation-of-state (EOS) measurements along the principal Hugoniot. In this Letter, we present one indirect-drive implosion experiment with a polystyrene sample that employs radiographic compression measurements over a range of shock pressures reaching up to 60 Mbar (6 TPa). This significantly exceeds previously published results obtained on the Nova laser [R. Cauble et al., Phys. Rev. Lett. 80, 1248 (1998)PRLTAO0031-900710.1103/PhysRevLett.80.1248] at a strongly improved precision, allowing us to discriminate between different EOS models. We find excellent agreement with Kohn-Sham density-functional-theory-based molecular dynamics simulations.

First-principles equation of state and shock compression predictions of warm dense hydrocarbons

We use path integral Monte Carlo and density functional molecular dynamics to construct a coherent set of equations of state (EOS) for a series of hydrocarbon materials with various C:H ratios (2:1, 1:1, 2:3, 1:2, and 1:4) over the range of 0.07-22.4gcm^{-3} and 6.7×10^{3}-1.29×10^{8}K. The shock Hugoniot curve derived for each material displays a single compression maximum corresponding to K-shell ionization. For C:H = 1:1, the compression maximum occurs at 4.7-fold of the initial density and we show radiation effects significantly increase the shock compression ratio above 2 Gbar, surpassing relativistic effects. The single-peaked structure of the Hugoniot curves contrasts with previous work on higher-Z plasmas, which exhibit a two-peak structure corresponding to both K- and L-shell ionization. Analysis of the electronic density of states reveals that the change in Hugoniot structure is due to merging of the L-shell eigenstates in carbon, while they remain distinct for higher-Z elements. Finally, we show that the isobaric-isothermal linear mixing rule for carbon and hydrogen EOS is a reasonable approximation with errors better than 1% for stellar-core conditions.

First-principles equation-of-state table of beryllium based on density-functional theory calculations

Beryllium has been considered a superior ablator material for inertial confinement fusion (ICF) target designs. An accurate equation-of-state (EOS) of beryllium under extreme conditions is essential for reliable ICF designs. Based on density-functional theory (DFT) calculations, we have established a wide-range beryllium EOS table of density ρ = 0.001 to 500 g/cm3 and temperature T = 2000 to 108 K. Our first-principle equation-of-state (FPEOS) table is in better agreement with the widely used SESAME EOS table (SESAME 2023) than the average-atom INFERNO and Purgatorio models. For the principal Hugoniot, our FPEOS prediction shows ∼10% stiffer than the last two models in the maximum compression. Although the existing experimental data (only up to 17 Mbar) cannot distinguish these EOS models, we anticipate that high-pressure experiments at the maximum compression region should differentiate our FPEOS from INFERNO and Purgatorio models. Comparisons between FPEOS and SESAME EOS for off-Hugoniot conditions show that the differences in the pressure and internal energy are within ∼20%. By implementing the FPEOS table into the 1-D radiation-hydrodynamic code LILAC, we studied the EOS effects on beryllium-shell-target implosions. The FPEOS simulation predicts higher neutron yield (∼15%) compared to the simulation using the SESAME 2023 EOS table.

Dissociation along the principal Hugoniot of the Laser Mégajoule ablator material

Glow discharge polymer hydrocarbon (GDP-CH) is used as the ablator material in inertial confinement fusion (ICF) capsules for the Laser Mégajoule and National Ignition Facility. Due to its fabrication process, GDP-CH chemical composition and structure differ from commercially available plastics and detailed knowledge of its properties in the warm dense matter regime is needed to achieve accurate design of ICF capsules. First-principles ab initio simulations of the GDP-CH principal Hugoniot up to 8 Mbar were performed using the quantum molecular dynamics (QMD) code abinit and showed that atomic bond dissociation has an effect on the compressibility. Results from these simulations are used to parametrize a quantum semiempirical model in order to generate a tabulated equation of state that includes dissociation. Hugoniot measurements obtained from an experiment conducted at the LULI2000 laser facility confirm QMD simulations as well as EOS modeling. We conclude by showing the EOS model influence on shock timing in a hydrodynamic simulation.

Experimental and ab initio investigations of microscopic properties of laser-shocked Ge-doped ablator

Plastic materials (CH) doped with mid-Z elements are used as ablators in inertial confinement fusion (ICF) capsules and in their surrogates. Hugoniot equation of state (EOS) and electronic properties of CH doped with germanium (at 2.5% and 13% dopant fractions) are investigated experimentally up to 7 Mbar using velocity and reflectivity measurements of shock fronts on the GEKKO laser at Osaka University. Reflectivity and temperature measurements were updated using a quartz standard. Shocked quartz reflectivity was measured at 532 and 1064 nm. Theoretical investigation of shock pressure and reflectivity was then carried out by ab initio simulations using the quantum molecular dynamics (QMD) code abinit and compared with tabulated average atom EOS models. We find that shock states calculated by QMD are in better agreement with experimental data than EOS models because of a more accurate description of ionic structure. We finally discuss electronic properties by comparing reflectivity data to a semiconductor gap closure model and to QMD simulations.

Impact of oxygen on the 300-K isotherm of Laser Megajoule ablator using ab initio simulation

The ablator material for inertial confinement fusion (ICF) capsules on the Laser Mégajoule is a glow-discharge polymer (GDP) plastic. Its equation of state (EOS) is of primary importance for the design of such capsules, since it has direct consequences on shock timing and is essential to mitigate hydrodynamic instabilities. Using ab initio molecular dynamics (AIMD), we have investigated the 300-K isotherm of amorphous CH(1.37)O(0.08) plastic, whose structure is close to GDP plastic. The 300-K isotherm, which is often used as a cold curve within tabular EOS, is an important contribution of the EOS in the multimegabar pressure range. AIMD results are compared to analytic models within tabular EOS, pointing out large discrepancies. In addition, we show that the effect of oxygen decreases 300-K isotherm pressure by 10%-15%. The implication of these observations is the ability to improve ICF target performance, which is essential to achieve fusion ignition.

First-principles equation of state of polystyrene and its effect on inertial confinement fusion implosions

Obtaining an accurate equation of state (EOS) of polystyrene (CH) is crucial to reliably design inertial confinement fusion (ICF) capsules using CH/CH-based ablators. With first-principles calculations, we have investigated the extended EOS of CH over a wide range of plasma conditions (ρ=0.1to100g/cm(3) and T=1000 to 4,000,000 K). When compared with the widely used SESAME-EOS table, the first-principles equation of state (FPEOS) of CH has shown significant differences in the low-temperature regime, in which strong coupling and electron degeneracy play an essential role in determining plasma properties. Hydrodynamic simulations of cryogenic target implosions on OMEGA using the FPEOS table of CH have predicted ∼30% decrease in neutron yield in comparison with the usual SESAME simulations. This is attributed to the ∼5% reduction in implosion velocity that is caused by the ∼10% lower mass ablation rate of CH predicted by FPEOS. Simulations using CH-FPEOS show better agreement with measurements of Hugoniot temperature and scattered light from ICF implosions.

Transport properties and equation of state for HCNO mixtures in and beyond the warm dense matter regime

We present simulations of a four-component mixture of HCNO with orbital free molecular dynamics (OFMD). These simulations were conducted for 5-200 eV with densities ranging between 0.184 and 36.8 g/cm3. We extract the equation of state from the simulations and compare to average atom models. We found that we only need to add a cold curve model to find excellent agreement. Additionally, we studied mass transport properties. We present fits to the self-diffusion and shear viscosity that are able to reproduce the transport properties over the parameter range studied. We compare these OFMD results to models based on the Coulomb coupling parameter and one-component plasmas.

Equation of state of a dense plasma by orbital-free and quantum molecular dynamics: examination of two isothermal-isobaric mixing rules

We test two isothermal-isobaric mixing rules, respectively based on excess-pressure and total-pressure equilibration, applied to the equation of state of a dense plasma. While the equation of state is generally known for pure species, that of arbitrary mixtures is not available so that the validation of accurate mixing rules, that implies resorting to first-principles simulations, is very useful. Here we consider the case of a plastic with composition C(2)H(3) and we implement two complementary ab initio approaches adapted to the dense plasma domain: quantum molecular dynamics, limited to low temperature by its computational cost, and orbital-free molecular dynamics, that can be implemented at high temperature. The temperature and density range considered is 1-10 eV and 0.6-10 g/cm(3) for quantum molecular dynamics, and 5-1000 eV and 1-10 g/cm(3) for orbital-free molecular dynamics. Simulations for the full C(2)H(3) mixture are the benchmark against which to assess the mixing rules, and both pressure and internal energy are compared. We find that the mixing rule based on excess-pressure equilibration is overall more accurate than that based on total-pressure equilibration; except for quantum molecular dynamics and a thermodynamic domain characterized by very low or negative excess pressures, it gives pressures which are generally within statistical error or within 1% of the exact ones. Besides, its superiority is amplified in the calculation of a principal Hugoniot.