Predictive design of REHEDS experiments with radiation-hydrodynamic simulations requires knowledge of material properties (e.g. equations of state (EOS), transport coefficients, and radiation physics). Interpreting experimental results requires accurate models of diagnostic observables (e.g. detailed emission, absorption, and scattering spectra). In conditions of Local Thermodynamic Equilibrium (LTE), these material properties and observables can be pre-computed with relatively high accuracy and subsequently tabulated on simple temperature-density grids for fast look-up by simulations. When radiation and electron temperatures fall out of equilibrium, however, non-LTE effects can profoundly change material properties and diagnostic signatures. Accurately and efficiently incorporating these non-LTE effects has been a longstanding challenge for simulations. At present, most simulations include non-LTE effects by invoking highly simplified inline models. These inline non-LTE models are both much slower than table look-up and significantly less accurate than the detailed models used to populate LTE tables and diagnose experimental data through post-processing or inversion. Because inline non-LTE models are slow, designers avoid them whenever possible, which leads to known inaccuracies from using tabular LTE. Because inline models are simple, they are inconsistent with tabular data from detailed models, leading to ill-known inaccuracies, and they cannot generate detailed synthetic diagnostics suitable for direct comparisons with experimental data. This project addresses the challenge of generating and utilizing efficient, accurate, and consistent non-equilibrium material data along three complementary but relatively independent research lines. First, we have developed a relatively fast and accurate non-LTE average-atom model based on density functional theory (DFT) that provides a complete set of EOS, transport, and radiative data, and have rigorously tested it against more sophisticated first-principles multi-atom DFT models, including time-dependent DFT. Next, we have developed a tabular scheme and interpolation methods that compactly capture non-LTE effects for use in simulations and have implemented these tables in the GORGON magneto-hydrodynamic (MHD) code. Finally, we have developed post-processing tools that use detailed tabulated non-LTE data to directly predict experimental observables from simulation output.
We present a new analysis methodology that allows for the self-consistent integration of multiple diagnostics including nuclear measurements, x-ray imaging, and x-ray power detectors to determine the primary stagnation parameters, such as temperature, pressure, stagnation volume, and mix fraction in magnetized liner inertial fusion (MagLIF) experiments. The analysis uses a simplified model of the stagnation plasma in conjunction with a Bayesian inference framework to determine the most probable configuration that describes the experimental observations while simultaneously revealing the principal uncertainties in the analysis. We validate the approach by using a range of tests including analytic and three-dimensional MHD models. An ensemble of MagLIF experiments is analyzed, and the generalized Lawson criterion χ is estimated for all experiments.
Laboratory experiments typically test opacity models by measuring spectrally resolved transmission of a sample using bright backlight radiation. A potential problem is that any unaccounted background signal contaminating the spectrum will artificially reduce the inferred opacity. Methods developed to measure background signals in opacity experiments at the Sandia Z facility are discussed. Preliminary measurements indicate that backgrounds are 9%–11% of the backlight signal at wavelengths less than 10 Å. Background is thus a relatively modest correction for all Z opacity data published to date. Future work will determine how important background is at longer wavelengths.
White Dwarf (WD) stars are the most common stellar remnant in the universe. WDs usually have a hydrogen or helium atmosphere, and helium WD (called DB) spectra can be used to solve outstanding problems in stellar and galactic evolution. DB origins, which are still a mystery, must be known to solve these problems. DB masses are crucial for discriminating between different proposed DB evolutionary hypotheses. Current DB mass determination methods deliver conflicting results. The spectroscopic mass determination method relies on line broadening models that have not been validated at DB atmosphere conditions. We performed helium benchmark experiments using the White Dwarf Photosphere Experiment (WDPE) platform at Sandia National Laboratories' Z-machine that aims to study He line broadening at DB conditions. Using hydrogen/helium mixture plasmas allows investigating the importance of He Stark and van der Waals broadening simultaneously. Accurate experimental data reduction methods are essential to test these line-broadening theories. In this paper, we present data calibration methods for these benchmark He line shape experiments. We give a detailed account of data processing, spectral power calibrations, and instrument broadening measurements. Uncertainties for each data calibration step are also derived. We demonstrate that our experiments meet all benchmark experiment accuracy requirements: WDPE wavelength uncertainties are <1 Å, spectral powers can be determined to within 15%, densities are accurate at the 20% level, and instrumental broadening can be measured with 20% accuracy. Fulfilling these stringent requirements enables WDPE experimental data to provide physically meaningful conclusions about line broadening at DB conditions.
I am a staff scientist at Sandia National Laboratories (SNL). I work on multiple fundamental-science projects and lead modeling/data-analysis for the stellar opacity experiments. I also work with theorists to refine plasma material-property calculations to make NNSA simulations more predictable. These are challenging and important problems for national security. In fact, our experiments raised questions about “opacity”, the property of matter that controls energy transport inside stars. This work was published by Nature due to its serious implication over broad applications, and our team received NNSA Defense Program Award of Excellence, both in 2015. Work at SNL is collegial and full of learning through interactions with the world’s finest experimentalists and theorists. When I started my Ph.D. program, working at a national lab was a dream position. Today, I collaborate with top scientists to solve mission-critical problems. The Stewardship Science Academic Programs (SSAP) and National Laser Users’ Facility (NLUF) played a vital role in the efficient learning and seamless transition to SNL.
The magnetized liner inertial fusion (MagLIF) scheme relies on coupling laser energy into an underdense fuel raising the fuel adiabat at the start of the implosion. To deposit energy into the fuel, the laser must first penetrate a laser entrance hole (LEH) foil which can be a significant energy sink and introduce mix. In this paper, we report on experiments investigating laser energy coupling into MagLIF-relevant gas cell targets with LEH foil thicknesses varying from 0.5 μm to 3 μm. Two-dimensional (2D) axisymmetric simulations match the experimental results well for 0.5 μm and 1 μm thick LEH foils but exhibit whole-beam self-focusing and excessive penetration of the laser into the gas for 2 μm and 3 μm thick LEH foils. Better agreement for the 2 μm-thick foil is found when using a different thermal conductivity model in 2D simulations, while only 3D Cartesian simulations come close to matching the 3 μm-thick foil experiments. The study suggests that simulations may over-predict the tendency for the laser to self-focus during MagLIF preheat when thicker LEH foils are used. This effect is pronounced with 2D simulations where the azimuthally symmetric density channel effectively self-focuses the rays that are forced to traverse the center of the plasma. The extra degree of freedom in 3D simulations significantly reduces this effect. The experiments and simulations also suggest that, in this study, the amount of energy coupled into the gas is highly correlated with the laser propagation length regardless of the LEH foil thickness.
The spectroscopic method relies on hydrogen Balmer absorption lines to infer white dwarf (WD) masses. These masses depend on the choice of atmosphere model, hydrogen atomic line shape calculation, and which Balmer series members are included in the spectral fit. In addition to those variables, spectroscopic masses disagree with those derived using other methods. Here we present laboratory experiments aimed at investigating the main component of the spectroscopic method: hydrogen line shape calculations. These experiments use X-rays from Sandia National Laboratories' Z-machine to create a uniform ∼15 cm3 hydrogen plasma and a ∼4 eV backlighter that enables recording high-quality absorption spectra. The large plasma, volumetric X-ray heating that fosters plasma uniformity, and the ability to collect absorption spectra at WD photosphere conditions are improvements over past laboratory experiments. Analysis of the experimental absorption spectra reveals that electron density (ne ) values derived from the Hγ line are ∼34% ± 7.3% lower than from Hβ. Two potential systematic errors that may contribute to this difference were investigated. A detailed evaluation of self-emission and plasma gradients shows that these phenomena are unlikely to produce any measurable Hβ-Hγ ne difference. WD masses inferred with the spectroscopic method are proportional to the photosphere density. Hence, the measured Hβ-Hγ n rm e difference is qualitatively consistent with the trend that WD masses inferred from their Hβ line are higher than that resulting from the analysis of Hβ and Hγ. This evidence may suggest that current hydrogen line shape calculations are not sufficiently accurate to capture the intricacies of the Balmer series.
This report will describe an improved computer code for two-photon opacity. The new code incorporates many recent advances and is ready to start to face the experiments. It incorporates the difficult mathematical techniques for handling free states and free-free matrix elements.
The burning core of an inertial confinement fusion (ICF) plasma produces bright x-rays at stagnation that can directly diagnose core conditions essential for comparison to simulations and understanding fusion yields. These x-rays also backlight the surrounding shell of warm, dense matter, whose properties are critical to understanding the efficacy of the inertial confinement and global morphology. We show that the absorption and fluorescence spectra of mid-Z impurities or dopants in the warm dense shell can reveal the optical depth, temperature, and density of the shell and help constrain models of warm, dense matter. This is illustrated by the example of a high-resolution spectrum collected from an ICF plasma with a beryllium shell containing native iron impurities. Analysis of the iron K-edge provides model-independent diagnostics of the shell density (2.3 × 1024 e/cm3) and temperature (10 eV), while a 12-eV red shift in Kβ and 5-eV blue shift in the K-edge discriminate among models of warm dense matter: Both shifts are well described by a self-consistent field model based on density functional theory but are not fully consistent with isolated-atom models using ad-hoc density effects.
The MagLIF campaign operated by Sandia conducted a total of four shot days in FY17 (one on OMEGA and three on OMEGA-EP) aimed at characterizing the laser heating of underdense plasmas (D2, Ar) at parameters that are relevant to the Magnetized Liner Inertial Fusion (MagLIF) ICF scheme being pursued at Sandia National Laboratories [1] [2]. MagLIF combines fuel preheat, magnetization and pulsed power implosion to significantly relax the implosion velocity and pR required for self-heating. Effective fuel preheat requires coupling several kJ of laser energy into the 10 mm long, underdense (typically ne/nc<0.1) fusion fuel without introducing significant mix. Barriers to achieving this include the presence laser plasma instabilities (LPI) as laser energy is coupled to the initially cold fuel, and the presence of a thin, polyimide laser entrance hole (LEH) foil that the laser must pass through and that can be a significant perturbation.
The recent iron opacity measurements performed at Sandia National Laboratory by Bailey and collaborators have raised questions about the completeness of the physical models normally used to understand partially ionized hot dense plasmas. We describe calculations of two-photon absorption, which is a candidate for the observed extra opacity. Our calculations do not yet match the experiments but show that the two-photon absorption process is strong enough to require careful consideration.
The interpretation of x-ray spectra emerging from x-ray binaries and active galactic nuclei accreted plasmas relies on complex physical models for radiation generation and transport in photoionized plasmas. These models have not been sufficiently experimentally validated. We have developed a highly reproducible benchmark experiment to study spectrum formation from a photoionized silicon plasma in a regime comparable to astrophysical plasmas. Ionization predictions are higher than inferred from measured absorption spectra. Self-emission measured at adjustable column densities tests radiation transport effects, demonstrating that the resonant Auger destruction assumption used to interpret black hole accretion spectra is inaccurate.
Iron opacity calculations presently disagree with measurements at an electron temperature of ∼180-195 eV and an electron density of (2-4)×1022cm-3, conditions similar to those at the base of the solar convection zone. The measurements use x rays to volumetrically heat a thin iron sample that is tamped with low-Z materials. The opacity is inferred from spectrally resolved x-ray transmission measurements. Plasma self-emission, tamper attenuation, and temporal and spatial gradients can all potentially cause systematic errors in the measured opacity spectra. In this article we quantitatively evaluate these potential errors with numerical investigations. The analysis exploits computer simulations that were previously found to reproduce the experimentally measured plasma conditions. The simulations, combined with a spectral synthesis model, enable evaluations of individual and combined potential errors in order to estimate their potential effects on the opacity measurement. The results show that the errors considered here do not account for the previously observed model-data discrepancies.
The MAGnetized Liner Inertial Fusion (MagLIF) scheme has achieved thermonuclear fusion yields on the Z Facility by imploding a cylindrical liner filled with D 2 fuel that is preheated with a multi-kJ laser and pre-magnetized with an axial Bz=10 T magnetic field. Preheating (T e = 100- 200 eV) and pre-magnetizing (10-30 T) the fuel serves to reduce the implosion velocity required to achieve multi-keV fusion-relevant temperatures at stagnation with a modest radial convergence. The challenge of fuel preheat in MagLIF is to deposit multiple-kJ of energy into the underdense (n e /n c %3C0.1) fuel over %7E10 mm target length efficiently and without introducing contaminants. Once the fuel is heated the applied axial magnetic field (o ce t e %7E 10) needs to suppress electron thermal conduction sufficiently to prevent unacceptable heat losses to the liner walls. In this LDRD we investigated laser energy deposition at two facilities: The OMEGA-EP laser at the Laboratory for Laser Energetics and the Z-beamlet laser at Sandia National Labs utilizing the PECOS chamber. Multiple experiments were carried out investigating laser transmission through LEH foils, laser heating of underdense gasses and the effects of magnetization on laser preheat. The studies find that magneto-hydrodynamic simulations are able to reproduce energy deposition at MagLIF-like conditions but that at the intensities currently used to preheat MagLIF significant laser plasma instabilities (LPI) occur which partly explain the inability of codes to reproduce previous MagLIF preheat studies. The experiments find that reducing the intensity and smoothing the beam dramatically reduces the amount of stimulated Brillouin backscatter and produces deposition profiles more similar to those produced in simulations. The experiments have provided a large and varied dataset that can be compared to simulations. As part of the LDRD new experimental capabilities have also been developed that will be used to design future MagLIF integrated experiments and investigate fuel magnetization.
This powerpoint presentation goes over the Fe opacity measurement platform, including how the experiment works, what can be gathered from the measurements, what can be gathered from the simulations, and the limitations of the experiment.
We present a platform on the OMEGA EP Laser Facility that creates and diagnoses the conditions present during the preheat stage of the MAGnetized Liner Inertial Fusion (MagLIF) concept. Experiments were conducted using 9 kJ of 3ω (355 nm) light to heat an underdense deuterium gas (electron density: 2.5×1020 cm-3=0.025 of critical density) magnetized with a 10 T axial field. Results show that the deuterium plasma reached a peak electron temperature of 670 ± 140 eV, diagnosed using streaked spectroscopy of an argon dopant. The results demonstrate that plasmas relevant to the preheat stage of MagLIF can be produced at multiple laser facilities, thereby enabling more rapid progress in understanding magnetized preheat. Results are compared with magneto-radiation-hydrodynamics simulations, and plans for future experiments are described.
We spectroscopically measure multiple hydrogen Balmer line profiles from laboratory plasmas to investigate the theoretical line profiles used in white dwarf (WD) atmosphere models. X-ray radiation produced at the Z Pulsed Power Facility at Sandia National Laboratories initiates plasma formation in a hydrogen-filled gas cell, replicating WD photospheric conditions. Here we present time-resolved measurements of Hβ and fit this line using different theoretical line profiles to diagnose electron density, ne, and n = 2 level population, n2. Aided by synthetic tests, we characterize the validity of our diagnostic method for this experimental platform. During a single experiment, we infer a continuous range of electron densities increasing from ne ∼ 4 to ∼30 × 1016 cm-3 throughout a 120-ns evolution of our plasma. Also, we observe n2 to be initially elevated with respect to local thermodynamic equilibrium (LTE); it then equilibrates within ∼55 ns to become consistent with LTE. This supports our electron-temperature determination of Te ∼ 1.3 eV (∼15,000 K) after this time. At ne 1017 cm-3, we find that computer-simulation-based line-profile calculations provide better fits (lower reduced χ2) than the line profiles currently used in the WD astronomy community. The inferred conditions, however, are in good quantitative agreement. This work establishes an experimental foundation for the future investigation of relative shapes and strengths between different hydrogen Balmer lines.