-
Shock Hugoniot of diamond from 3 to 80 TPa
Authors:
Damian C. Swift,
Andrea L. Kritcher,
Amy Lazicki,
James A. Hawreliak,
Tilo Doeppner,
Heather D. Whitley,
Joseph Nilsen,
Benjamin Bachmann,
Michael MacDonald,
Brian Maddox,
Natalie Kostinski,
Gilbert W. Collins,
Siegfried Glenzer,
Stephen D. Rothman,
Dominik Kraus,
Roger W. Falcone
Abstract:
The principal Hugoniot of carbon, initially diamond, was measured from 3 to 80 TPa (30 to 800 million atmospheres), the highest pressure ever achieved, using radiography of spherically-converging shocks. The shocks were generated by ablation of a plastic coating by soft x-rays in a laser-heated hohlraum at the National Ignition Facility (NIF). Experiments were performed with low and high drive pow…
▽ More
The principal Hugoniot of carbon, initially diamond, was measured from 3 to 80 TPa (30 to 800 million atmospheres), the highest pressure ever achieved, using radiography of spherically-converging shocks. The shocks were generated by ablation of a plastic coating by soft x-rays in a laser-heated hohlraum at the National Ignition Facility (NIF). Experiments were performed with low and high drive powers, spanning different but overlapping pressure ranges. The radius-time history of the shock, and the profile of mass density behind, were determined by profile-matching from a time-resolved x-ray radiograph across the diameter of the sphere. Above ~50 TPa, the heating induced by the shock was great enough to ionize a significant fraction of K-shell electrons, reducing the opacity to the 10.2 keV probe x-rays. The opacity and mass density were deduced simultaneously using the constraint that the total mass of the sample was constant. The Hugoniot and opacity were consistent with density functional theory calculations of the electronic states and equation of state (EOS), and varied significantly from theoretical Hugoniots based on Thomas-Fermi theory. Theoretical models used to predict the compressibility of diamond ablator experiments at the NIF, producing the highest neutron yields so far from inertial confinement fusion experiments, are qualitatively consistent with our EOS measurements but appear to overpredict the compressibility slightly. These measurements help to evaluate theoretical techniques and constrain wide-range EOS models applicable to white dwarf stars, which are the ultimate evolutionary form of at least 97% of stars in the galaxy.
△ Less
Submitted 30 March, 2022;
originally announced March 2022.
-
Improved analysis of converging shock radiographs
Authors:
Damian C. Swift,
Andrea L. Kritcher,
Amy Lazicki,
Natalie Kostinski,
Brian R. Maddox,
Madison E. Martin,
Tilo Doeppner,
Joseph Nilsen,
Heather D. Whitley
Abstract:
We previously reported an experimental platform to induce a spherically-convergent shock in a sample using laser-driven ablation, probed with time-resolved x-ray radiography, and an analysis method to deduce states along the principal shock Hugoniot simultaneously with the x-ray opacity. We have now developed a modified method of analysis that is numerically better-conditioned and faster, and usua…
▽ More
We previously reported an experimental platform to induce a spherically-convergent shock in a sample using laser-driven ablation, probed with time-resolved x-ray radiography, and an analysis method to deduce states along the principal shock Hugoniot simultaneously with the x-ray opacity. We have now developed a modified method of analysis that is numerically better-conditioned and faster, and usually provides a better representation of the radiograph with correspondingly lower uncertainties. The previous approach was based on optimizing parameters in a model of the density distribution as a function of radius and time, warped to follow loci such as the shock and the outside of the sample. The converging shock configuration can be described more efficiently in terms of the shocked density and sound speed, expressed as functions of the shock speed Studies of the Hugoniot from various theoretical equations of state (EOS) indicate that, in the typical range of states explored by these experiments, these functions can be described by low-order polynomials. Similarly, few-parameter functions were found suitable for representing the variation of x-ray opacity with shock pressure. This approach was found to perform better in most cases than an alternative method based on parameterization of the EOS.
△ Less
Submitted 16 March, 2022;
originally announced March 2022.
-
Development of slurry targets for high repetition-rate XFEL experiments
Authors:
Raymond F. Smith,
Vinay Rastogi,
Amy E. Lazicki,
Martin G. Gorman,
Richard Briggs,
Amy L. Coleman,
Carol Davis,
Saransh Singh,
David McGonegle,
Samantha M. Clarke,
Travis Volz,
Trevor Hutchinson,
Christopher McGuire,
Dayne E. Fratanduono,
Damian C. Swift,
Eric Folsom,
Cynthia A. Bolme,
Arianna E. Gleason,
Federica Coppari,
Hae Ja Lee,
Bob Nagler,
Eric Cunningham,
Eduardo Granados,
Phil Heimann,
Richard G. Kraus
, et al. (4 additional authors not shown)
Abstract:
Combining an x-ray free electron laser (XFEL) with high power laser drivers enables the study of phase transitions, equation-of-state, grain growth, strength, and transformation pathways as a function of pressure to 100s GPa along different thermodynamic compression paths. Future high-repetition rate laser operation will enable data to be accumulated at >1 Hz which poses a number of experimental c…
▽ More
Combining an x-ray free electron laser (XFEL) with high power laser drivers enables the study of phase transitions, equation-of-state, grain growth, strength, and transformation pathways as a function of pressure to 100s GPa along different thermodynamic compression paths. Future high-repetition rate laser operation will enable data to be accumulated at >1 Hz which poses a number of experimental challenges including the need to rapidly replenish the target. Here, we present a combined shock-compression and X-ray diffraction study on vol% epoxy(50)-crystalline grains(50) (slurry) targets, which can be fashioned into extruded ribbons for high repetition-rate operation. For shock-loaded NaCl-slurry samples, we observe pressure, density and temperature states within the embedded NaCl grains consistent with observations for shock-compressed single-crystal NaCl.
△ Less
Submitted 11 January, 2022;
originally announced January 2022.
-
Benchmarking boron carbide equation of state using computation and experiment
Authors:
Shuai Zhang,
Michelle C. Marshall,
Lin H. Yang,
Philip A. Sterne,
Burkhard Militzer,
Markus Daene,
James A. Gaffney,
Andrew Shamp,
Tadashi Ogitsu,
Kyle Caspersen,
Amy E. Lazicki,
David Erskine,
Richard A. London,
Peter M. Celliers,
Joseph Nilsen,
Heather D. Whitley
Abstract:
Boron carbide (B$_4$C) is of both fundamental scientific and practical interest 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 t…
▽ More
Boron carbide (B$_4$C) is of both fundamental scientific and practical interest 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 ($\sim$6$\times$10$^3$--5$\times$10$^8$ K) and densities (0.025--50 g/cm$^{3}$). We assess that the largest discrepancies between theoretical predictions are $\lesssim$5% near the compression maximum at 1--2$\times10^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 $\sim$18% and occur at temperatures between 6$\times$10$^3$--2$\times$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. In addition, we have developed three new equation of state models and applied them to 1D hydrodynamic simulations of a polar direct-drive NIF implosion, demonstrating that these new models are now available for future ICF design studies.
△ Less
Submitted 3 August, 2020;
originally announced August 2020.
-
Comparison of ablators for the polar direct drive exploding pusher platform
Authors:
Heather D. Whitley,
G. Elijah Kemp,
Charles Yeamans,
Zachary Walters,
Brent E. Blue,
Warren Garbett,
Marilyn Schneider,
R. Stephen Craxton,
Emma M. Garcia,
Patrick W. McKenty,
Maria Gatu-Johnson,
Kyle Caspersen,
John I. Castor,
Markus Däne,
C. Leland Ellison,
James Gaffney,
Frank R. Graziani,
John Klepeis,
Natalie Kostinski,
Andrea Kritcher,
Brandon Lahmann,
Amy E. Lazicki,
Hai P. Le,
Richard A. London,
Brian Maddox
, et al. (14 additional authors not shown)
Abstract:
We examine the performance of pure boron, boron carbide, high density carbon, and boron nitride ablators in the polar direct drive exploding pusher (PDXP) platform. The platform uses the polar direct drive configuration at the National Ignition Facility to drive high ion temperatures in a room temperature capsule and has potential applications for plasma physics studies and as a neutron source. Th…
▽ More
We examine the performance of pure boron, boron carbide, high density carbon, and boron nitride ablators in the polar direct drive exploding pusher (PDXP) platform. The platform uses the polar direct drive configuration at the National Ignition Facility to drive high ion temperatures in a room temperature capsule and has potential applications for plasma physics studies and as a neutron source. The higher tensile strength of these materials compared to plastic enables a thinner ablator to support higher gas pressures, which could help optimize its performance for plasma physics experiments, while ablators containing boron enable the possiblity of collecting addtional data to constrain models of the platform. Applying recently developed and experimentally validated equation of state models for the boron materials, we examine the performance of these materials as ablators in 2D simulations, with particular focus on changes to the ablator and gas areal density, as well as the predicted symmetry of the inherently 2D implosion.
△ Less
Submitted 30 December, 2020; v1 submitted 28 June, 2020;
originally announced June 2020.
-
Simultaneous compression and opacity data from time-series radiography with a Lagrangian marker
Authors:
Damian C. Swift,
Andrea L. Kritcher,
James A. Hawreliak,
James Gaffney,
Amy Lazicki,
Andrew MacPhee,
Benjamin Bachmann,
Tilo Doeppner,
Joseph Nilsen,
Heather D. Whitley,
Gilbert W. Collins,
Siegfried Glenzer,
Stephen D. Rothman,
Dominik Kraus,
Roger W. Falcone
Abstract:
Time-resolved radiography can be used to obtain absolute shock Hugoniot states by simultaneously measuring at least two mechanical parameters of the shock, and this technique is particularly suitable for one-dimensional converging shocks where a single experiment probes a range of pressures as the converging shock strengthens. However, at sufficiently high pressures, the shocked material becomes h…
▽ More
Time-resolved radiography can be used to obtain absolute shock Hugoniot states by simultaneously measuring at least two mechanical parameters of the shock, and this technique is particularly suitable for one-dimensional converging shocks where a single experiment probes a range of pressures as the converging shock strengthens. However, at sufficiently high pressures, the shocked material becomes hot enough that the x-ray opacity falls significantly. If the system includes a Lagrangian marker, such that the mass within the marker is known, this additional information can be used to constrain the opacity as well as the Hugoniot state. In the limit that the opacity changes only on shock heating, and not significantly on subsequent isentropic compression, the opacity of shocked material can be determined uniquely. More generally, it is necessary to assume the form of the variation of opacity with isentropic compression, or to introduce multiple marker layers. Alternatively, assuming either the equation of state or the opacity, the presence of a marker layer in such experiments enables the non-assumed property to be deduced more accurately than from the radiographic density reconstruction alone. An example analysis is shown for measurements of a converging shock wave in polystyrene, at the National Ignition Facility.
△ Less
Submitted 15 October, 2020; v1 submitted 23 June, 2020;
originally announced June 2020.
-
Equation of state and strength of diamond in high pressure ramp loading
Authors:
Damian C. Swift,
Olivier Heuze,
Amy Lazicki,
Sebastien Hamel,
Lorin X. Benedict,
Raymond F. Smith,
James M. McNaney,
Graeme J. Ackland
Abstract:
Diamond is used extensively as a component in high energy density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first principles theoretical EOS models, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density f…
▽ More
Diamond is used extensively as a component in high energy density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first principles theoretical EOS models, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Greuneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. The resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.
△ Less
Submitted 31 August, 2021; v1 submitted 6 April, 2020;
originally announced April 2020.
-
Equation of state of warm-dense boron nitride combining computation, modeling, and experiment
Authors:
Shuai Zhang,
Amy Lazicki,
Burkhard Militzer,
Lin H. Yang,
Kyle Caspersen,
Jim A. Gaffney,
Markus W. Däne,
John E. Pask,
Walter R. Johnson,
Abhiraj Sharma,
Phanish Suryanarayana,
Duane D. Johnson,
Andrey V. Smirnov,
Philip A. Sterne,
David Erskine,
Richard A. London,
Federica Coppari,
Damian Swift,
Joseph Nilsen,
Art J. Nelson,
Heather D. Whitley
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 functio…
▽ More
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 ($ρ$) 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-ρ$ 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.
△ Less
Submitted 2 February, 2019;
originally announced February 2019.
-
Theoretical and experimental investigation of the equation of state of boron plasmas
Authors:
Shuai Zhang,
Burkhard Militzer,
Michelle C. Gregor,
Kyle Caspersen,
Lin H. Yang,
Tadashi Ogitsu,
Damian Swift,
Amy Lazicki,
D. Erskine,
Richard A. London,
P. M. Celliers,
Joseph Nilsen,
Philip A. Sterne,
Heather D. Whitley
Abstract:
We report a theoretical equation of state (EOS) table for boron across a wide range of temperatures (5.1$\times$10$^4$-5.2$\times$10$^8$ K) and densities (0.25-49 g/cm$^3$), and experimental shock Hugoniot data at unprecedented high pressures (5608$\pm$118 GPa). The calculations are performed with full, first-principles methods combining path integral Monte Carlo (PIMC) at high temperatures and de…
▽ More
We report a theoretical equation of state (EOS) table for boron across a wide range of temperatures (5.1$\times$10$^4$-5.2$\times$10$^8$ K) and densities (0.25-49 g/cm$^3$), and experimental shock Hugoniot data at unprecedented high pressures (5608$\pm$118 GPa). The calculations are performed with full, 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$\times$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 semi-empirical 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 study the performance sensitivity of a polar direct-drive exploding pusher platform to pressure variations based on comparison of the first-principles calculations with LEOS 50 via 1D hydrodynamic simulations. The results are valuable for future theoretical and experimental studies and engineering design in high energy density research. (LLNL-JRNL-748227)
△ Less
Submitted 30 April, 2018;
originally announced April 2018.
-
Analysis of laser shock experiments on precompressed samples using a quartz reference and application to warm dense hydrogen and helium
Authors:
Stephanie Brygoo,
Marius Millot,
Paul Loubeyre,
Amy E. Lazicki,
Sebastien Hamel,
Tingting Qi,
Peter M. Celliers,
Federica Coppari,
Jon H. Eggert,
Dayne E. Fratanduono,
Damien G. Hicks,
J. Ryan Rygg,
Raymond F. Smith,
Damian C. Swift,
Gilbert W. Collins,
Raymond Jeanloz
Abstract:
Megabar (1 Mbar = 100 GPa) laser shocks on precompressed samples allow reaching unprecedented high densities and moderately high 10000-100000K temperatures. We describe here a complete analysis framework for the velocimetry (VISAR) and pyrometry (SOP) data produced in these experiments. Since the precompression increases the initial density of both the sample of interest and the quartz reference f…
▽ More
Megabar (1 Mbar = 100 GPa) laser shocks on precompressed samples allow reaching unprecedented high densities and moderately high 10000-100000K temperatures. We describe here a complete analysis framework for the velocimetry (VISAR) and pyrometry (SOP) data produced in these experiments. Since the precompression increases the initial density of both the sample of interest and the quartz reference for pressure-density, reflectivity and temperature measurements, we describe analytical corrections based on available experimental data on warm dense silica and density-functional-theory based molecular dynamics computer simulations. Using our improved analysis framework we report a re-analysis of previously published data on warm dense hydrogen and helium, compare the newly inferred pressure, density and temperature data with most advanced equation of state models and provide updated reflectivity values.
△ Less
Submitted 12 October, 2015;
originally announced October 2015.
-
Superconductivity and Lattice Instability in Compressed Lithium from Fermi Surface Hot Spots
Authors:
Deepa Kasinathan,
J. Kunes,
A. Lazicki,
H. Rosner,
C. S. Yoo,
R. T. Scalettar,
W. E. Pickett
Abstract:
The highest superconducting temperature T$_c$ observed in any elemental metal (Li with T$_c$ ~ 20 K at pressure P ~ 40 GPa) is shown to arise from critical (formally divergent) electron-phonon coupling to the transverse T$_1$ phonon branch along intersections of Kohn anomaly surfaces with the Fermi surface. First principles linear response calculations of the phonon spectrum and spectral functio…
▽ More
The highest superconducting temperature T$_c$ observed in any elemental metal (Li with T$_c$ ~ 20 K at pressure P ~ 40 GPa) is shown to arise from critical (formally divergent) electron-phonon coupling to the transverse T$_1$ phonon branch along intersections of Kohn anomaly surfaces with the Fermi surface. First principles linear response calculations of the phonon spectrum and spectral function $α^2 F(ω)$ reveal (harmonic) instability already at 25 GPa. Our results imply that the fcc phase is anharmonically stabilized in the 25-38 GPa range.
△ Less
Submitted 20 October, 2005;
originally announced October 2005.