Search

Thermal conductivity is one of the most crucial physical properties of matter when it comes to understanding heat transport, hydrodynamic evolution, and energy balance in systems ranging from astrophysical objects to fusion plasmas. In the warm dense matter regime, experimental data are very scarce so that many theoretical models remain untested. Here we present the first thermal conductivity measurements of aluminum at 0.5-2.7 g/cc and 2-10 eV, using a recently developed platform of differential heating. A temperature gradient is induced in a Au/Al dual-layer target by proton heating, and subsequent heat flow from the hotter Au to the Al rear surface is detected by two simultaneous time-resolved diagnostics. A systematic data set allows for constraining both thermal conductivity and equation-of-state models. Simulations using Purgatorio model or Sesame S27314 for Al thermal conductivity and LEOS for Au/Al release equation-of-state show good agreement with data after 15 ps. Discrepancy still exists at early time 0-15 ps, likely due to non-equilibrium conditions.

The gravitational pressure in many astrophysical objects exceeds one gigabar (one billion atmospheres) 1-3 , creating extreme conditions where the distance between nuclei approaches the size of the K shell. This close proximity modifies these tightly bound states and, above a certain pressure, drives them into a delocalized state 4 . Both processes substantially affect the equation of state and radiation transport and, therefore, the structure and evolution of these objects. Still, our understanding of this transition is far from satisfactory and experimental data are sparse. Here we report on experiments that create and diagnose matter at pressures exceeding three gigabars at the National Ignition Facility 5 where 184 laser beams imploded a beryllium shell. Bright X-ray flashes enable precision radiography and X-ray Thomson scattering that reveal both the macroscopic conditions and the microscopic states. The data show clear signs of quantum-degenerate electrons in states reaching 30 times compression, and a temperature of around two million kelvins. At the most extreme conditions, we observe strongly reduced elastic scattering, which mainly originates from K-shell electrons. We attribute this reduction to the onset of delocalization of the remaining K-shell electron. With this interpretation, the ion charge inferred from the scattering data agrees well with ab initio simulations, but it is significantly higher than widely used analytical models predict 6 ., (© 2023. The Author(s), under exclusive licence to Springer Nature Limited.)

We report on the study of electron kinetics induced by intense femtosecond (fs) laser excitation of electrons in the 5d band of Au. Changes in the electron system are observed from the temporal evolution of ac conductivity and conduction electron density. The results reveal an increase of electron thermalization time with excitation energy density, contrary to the Fermi-liquid behavior of the decrease of thermalization time associated with the heating of conduction electrons. This is attributed to the severe mitigation of photoexcitation by Auger decay. The study also uncovers the shortening of 5d hole lifetime with the increase of photoexcitation rates. These unique findings provide valuable insights for understanding electron kinetics under extreme nonequilibrium conditions.

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.

Oxygen, the third most abundant element in the universe, plays a key role in the chemistry of condensed matter and biological systems. Here, we report evidence for a hitherto unexplored Auger transition in oxides, where a valence band electron fills a vacancy in the 2s state of oxygen, transferring sufficient energy to allow electron emission. We used a beam of positrons with kinetic energies of [Formula: see text] eV to create O 2s holes via matter-antimatter annihilation. This made possible the elimination of the large secondary electron background that has precluded definitive measurements of the low-energy electrons emitted through this process. Our experiments indicate that low-energy electron emission following the Auger decay of O 2s holes from adsorbed oxygen and oxide surfaces are very efficient. Specifically, our results indicate that the low energy electron emission following the Auger decay of O 2s hole is nearly as efficient as electron emission following the relaxation of O 1s holes in [Formula: see text]. This has important implications for the understanding of Auger-stimulated ion desorption, Coulombic decay, photodynamic cancer therapies, and may yield important insights into the radiation-induced reactive sites for corrosion and catalysis.

White dwarfs represent the final state of evolution for most stars 1-3 . Certain classes of white dwarfs pulsate 4,5 , leading to observable brightness variations, and analysis of these variations with theoretical stellar models probes their internal structure. Modelling of these pulsating stars provides stringent tests of white dwarf models and a detailed picture of the outcome of the late stages of stellar evolution 6 . However, the high-energy-density states that exist in white dwarfs are extremely difficult to reach and to measure in the laboratory, so theoretical predictions are largely untested at these conditions. Here we report measurements of the relationship between pressure and density along the principal shock Hugoniot (equations describing the state of the sample material before and after the passage of the shock derived from conservation laws) of hydrocarbon to within five per cent. The observed maximum compressibility is consistent with theoretical models that include detailed electronic structure. This is relevant for the equation of state of matter at pressures ranging from 100 million to 450 million atmospheres, where the understanding of white dwarf physics is sensitive to the equation of state and where models differ considerably. The measurements test these equation-of-state relations that are used in the modelling of white dwarfs and inertial confinement fusion experiments 7,8 , and we predict an increase in compressibility due to ionization of the inner-core orbitals of carbon. We also find that a detailed treatment of the electronic structure and the electron degeneracy pressure is required to capture the measured shape of the pressure-density evolution for hydrocarbon before peak compression. Our results illuminate the equation of state of the white dwarf envelope (the region surrounding the stellar core that contains partially ionized and partially degenerate non-ideal plasmas), which is a weak link in the constitutive physics informing the structure and evolution of white dwarf stars 9 .

Atom-in-jellium calculations of the Einstein frequency were used to calculate the mean displacement of an ion over a wide range of compression and temperature. Expressed as a fraction of the Wigner-Seitz radius, the displacement is a measure of the asymptotic freedom of the ion at high temperature, and thus of the change in heat capacity from six to three quadratic degrees of freedom per atom. A functional form for free energy was proposed based on the Maxwell-Boltzmann distribution as a correction to the Debye free energy, with a single free parameter representing the effective density of potential modes to be saturated. This parameter was investigated using molecular dynamics simulations, and found to be ∼0.2 per atom. In this way, the ion-thermal contribution can be calculated for a wide-range equation of state (EOS) without requiring a large number of molecular dynamics simulations. Example calculations were performed for carbon, including the sensitivity of key EOS loci to ionic freedom.

We report measurements of K-shell fluorescence lines induced by fast electrons in ramp-compressed Co targets. The fluorescence emission was stimulated by fast electrons generated through short-pulse laser-solid interaction with an Al target layer. Compression up to 2.1× solid density was achieved while maintaining temperatures well below the Fermi energy, effectively removing the thermal effects from consideration. We observed small but unambiguous redshifts in the Kβ fluorescence line relative to unshifted Cu Kα. Redshifts up to 2.6 eV were found to increase with compression and to be consistent with predictions from self-consistent models based on density-functional theory.

A notorious challenge in high-pressure science is to develop an equation of state (EOS) that explicitly treats chemical reactions. For instance, many materials tend to dissociate at high pressures and temperatures where the chemical bonds that hold them together break down. We present an EOS for carbon dioxide (CO 2 ) that allows for dissociation and captures the key material behavior in a wide range of pressure-temperature conditions. Carbon dioxide is an ideal prototype for the development of a wide-ranging EOS that allows for chemical-dissociation equilibria since it is one of the simplest polyatomic systems and because it is of great interest in planetary science and in the study of detonations. Here, we show that taking dissociation into account significantly improves the accuracy of the resulting EOS compared to other EOSs that either neglect chemistry completely or treat CO 2 dissociation in a more rudimentary way.

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. 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.

We present laser-driven shock compression experiments on cryogenic liquid deuterium to 550 GPa along the principal Hugoniot and reflected-shock data up to 1 TPa. High-precision interferometric Doppler velocimetry and impedance-matching analysis were used to determine the compression accurately enough to reveal a significant difference as compared to state-of-the-art ab initio calculations and thus, no single equation of state model fully matches the principal Hugoniot of deuterium over the observed pressure range. In the molecular-to-atomic transition pressure range, models based on density functional theory calculations predict the maximum compression accurately. However, beyond 250 GPa along the principal Hugoniot, first-principles models exhibit a stiffer response than the experimental data. Similarly, above 500 GPa the reflected shock data show 5%-7% higher compression than predicted by all current models.

Recent path-integral Monte Carlo and quantum molecular dynamics simulations have shown that computationally efficient average-atom models can predict thermodynamic states in warm dense matter to within a few percent. One such atom-in-jellium model has typically been used to predict the electron-thermal behavior only, although it was previously developed to predict the entire equation of state (EOS). We report completely atom-in-jellium EOS calculations for Be, Al, Si, Fe, and Mo, as elements representative of a range of atomic number and low-pressure electronic structure. Comparing the more recent method of pseudoatom molecular dynamics, atom-in-jellium results were similar: sometimes less accurate, sometimes more. All these techniques exhibited pronounced effects of electronic shell structure in the shock Hugoniot which are not captured by Thomas-Fermi based EOS. These results demonstrate the value of a hierarchical approach to EOS construction, using average-atom techniques with shell structure to populate a wide-range EOS surface efficiently, complemented by more rigorous three-dimensional multiatom calculations to validate and adjust the EOS.

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.

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.

Carbon-hydrogen plasmas and hydrocarbon materials are of broad interest to laser shock experimentalists, high energy density physicists, and astrophysicists. Accurate equations of state (EOSs) of hydrocarbons are valuable for various studies from inertial confinement fusion to planetary science. By combining path integral Monte Carlo (PIMC) results at high temperatures and density functional theory molecular dynamics results at lower temperatures, we compute the EOSs for hydrocarbons from simulations performed at 1473 separate (ρ, T)-points distributed over a range of compositions. These methods accurately treat electronic excitation effects with neither adjustable parameter nor experimental input. PIMC is also an accurate simulation method that is capable of treating many-body interaction and nuclear quantum effects at finite temperatures. These methods therefore provide a benchmark-quality EOS that surpasses that of semi-empirical and Thomas-Fermi-based methods in the warm dense matter regime. By comparing our first-principles EOS to the LEOS 5112 model for CH, we validate the specific heat assumptions in this model but suggest that the Grüneisen parameter is too large at low temperatures. Based on our first-principles EOSs, we predict the principal Hugoniot curve of polystyrene to be 2%-5% softer at maximum shock compression than that predicted by orbital-free density functional theory and SESAME 7593. By investigating the atomic structure and chemical bonding of hydrocarbons, we show a drastic decrease in the lifetime of chemical bonds in the pressure interval from 0.4 to 4 megabar. We find the assumption of linear mixing to be valid for describing the EOS and the shock Hugoniot curve of hydrocarbons in the regime of partially ionized atomic liquids. We make predictions of the shock compression of glow-discharge polymers and investigate the effects of oxygen content and C:H ratio on its Hugoniot curve. Our full suite of first-principles simulation results may be used to benchmark future theoretical investigations pertaining to hydrocarbon EOSs and should be helpful in guiding the design of future experiments on hydrocarbons in the gigabar regime.

The first gamma spectra associated with the annihilation of positrons with individual core levels (Cu 3p and Ag 4p) are presented. The spectra were obtained by measuring the energy of gamma rays time coincident with Auger electrons emitted as a result of positrons annihilating with a selected core level. Relativistic calculations show good agreement with experiment over a limited range of momenta. However, statistically significant differences indicate that the measurements can provide an impetus to new calculations of many body effects in positron-core electrons annihilation.

We report the first observation of electron momentum contributions to the Doppler broadening of radiation produced by in-flight two-photon annihilation in solids. In these experiments an approximately 2.5 MeV positron beam impinged on thin polyethylene, aluminum, and gold targets. Since energetic positrons easily penetrate the nuclear Coulomb potential and do not cause a strong charge polarization, the experimental annihilation line shapes agree well with calculations based on a simple independent-particle model. Moreover, annihilations with the deepest core electrons are greatly enhanced.