Search

We present a study of spin-unpolarized and spin-polarized two-dimensional uniform electron liquids using variational and diffusion quantum Monte Carlo (VMC and DMC) methods with Slater-Jastrow-backflow trial wave functions. Ground-state VMC and DMC energies are obtained in the density range $1 \leq r_\text{s} \leq 40$. Single-particle and many-body finite-size errors are corrected using canonical-ensemble twist-averaged boundary conditions and extrapolation of twist-averaged energies to the thermodynamic limit of infinite system size. System-size-dependent errors in Slater-Jastrow-backflow DMC energies caused by partially converged VMC energy minimization calculations are discussed. We find that, for $1 \leq r_\text{s} \leq 5$, optimizing the backflow function at each twist lowers the twist-averaged DMC energy at finite system size. However, nonsystematic system-size-dependent effects remain in the DMC energies, which can be partially removed by extrapolation from multiple finite system sizes to infinite system size. We attribute these nonsystematic effects to the close competition between fluid and defected crystal phases at different system sizes at low density. The DMC energies in the thermodynamic limit are used to parameterize a local spin density approximation correlation functional for inhomogeneous electron systems. Our zero-temperature phase diagram shows a single transition from a paramagnetic fluid to a hexagonal Wigner crystal at $r_\text{s}=35(1)$, with no region of stability for a ferromagnetic fluid.

We present in-situ x-ray diffraction and velocity measurements of Fe$_2$O$_3$ under laser shock compression at pressures between 38-116 GPa. None of the phases reported by static compression studies were observed. Instead, we observed an isostructural phase transition from $\alpha$-Fe$_2$O$_3$ to a new $\alpha^\prime$-Fe$_2$O$_3$ phase at a pressure of 50-62 GPa. The $\alpha^\prime$-Fe$_2$O$_3$ phase differs from $\alpha$-Fe$_2$O$_3$ by an 11% volume drop and a different unit cell compressibility. We further observed a two-wave structure in the velocity profile, which can be related to an intermediate regime where both $\alpha$ and $\alpha^\prime$ phases coexist. Density functional theory calculations with a Hubbard parameter indicate that the observed unit cell volume drop can be associated with a spin transition following a magnetic collapse., Comment: 7 pages, 5 figures

We report on experiments where solid-density Mg plasmas are created by heating with the focused output of the Linac Coherent Light Source x-ray free-electron-laser. We study the K-shell emission from the Helium and Lithium-like ions using Bragg crystal spectroscopy. Observation of the dielectronic satellites in Lithium-like ions confirms that the M-shell electrons appear bound for these high charge states. An analysis of the intensity of these satellites indicates that when modelled with an atomic-kinetics code, the ionisation potential depression model employed needs to produce depressions for these ions which lie between those predicted by the well known Stewart-Pyatt and Ecker-Kroll models. These results are largely consistent with recent Density Functional Theory calculations.

The variational and diffusion quantum Monte Carlo methods are used to calculate the correlation energy of the paramagnetic three-dimensional homogeneous electron gas at intermediate to high density. Ground state energies in finite cells are determined using Slater-Jastrow-backflow trial wave functions, and finite-size errors are removed using twist-averaged boundary conditions and extrapolation of the energy per particle to the thermodynamic limit of infinite system size. Our correlation energies in the thermodynamic limit are lower (i.e., more negative, and therefore more accurate according to the variational principle) than previous results, and can be used for the parameterization of density functionals to be applied to high-density systems.

The design of inertial confinement fusion (ICF) ignition targets requires radiation-hydrodynamics simulations with accurate models of the fundamental material properties (i.e., equation of state, opacity, and conductivity). Validation of these models are required via experimentation. A feasibility study of using spatially-integrated, spectrally-resolved, X-ray Thomson scattering (XRTS) measurements to diagnose the temperature, density, and ionization of the compressed DT shell and hot spot of a laser direct-drive implosion at two-thirds convergence was conducted. Synthetic scattering spectra were generated using 1-D implosion simulations from the LILAC code that were post processed with the X-ray Scattering (XRS) model which is incorporated within SPECT3D. Analysis of two extreme adiabat capsule conditions showed that the plasma conditions for both compressed DT shells could be resolved., Comment: 11 pages, 8 figures

We present picosecond-resolution streaked K-shell spectra from 400 nm-diameter nickel nanowire arrays, demonstrating the ability to generate large volumes of high energy density plasma when combined with the longer pulses typical of the largest short pulse lasers. After irradiating the wire array with 100 J, 600 fs ultra-high-contrast laser pulses focussed to $>10^{20}$ W/cm$^{2}$ at the Orion laser facility, we combine atomic kinetics modeling of the streaked spectra with 2D collisional particle-in-cell simulations to describe the evolution of material conditions within these samples for the first time. We observe a three-fold enhancement of helium-like emission compared to a flat foil in a near-solid-density plasma sustaining keV temperatures for tens of picoseconds, the result of strong electric return currents heating the wires and causing them to explode and collide., Comment: 5 pages, 5 figures

The free-free opacity in plasmas is fundamental to our understanding of energy transport in stellar interiors and for inertial confinement fusion research. However, theoretical predictions in the challenging dense plasma regime are conflicting and there is a dearth of accurate experimental data to allow for direct model validation. Here we present time-resolved transmission measurements in solid-density Al heated by an XUV free-electron laser. We use a novel functional optimization approach to extract the temperature-dependent absorption coefficient directly from an oversampled pool of single-shot measurements, and find a pronounced enhancement of the opacity as the plasma is heated to temperatures of order the Fermi energy. Plasma heating and opacity-enhancement is observed on ultrafast time scales, within the duration of the femtosecond XUV pulse. We attribute further rises in the opacity on ps timescales to melt and the formation of warm-dense matter.

Computer simulations are invaluable tools for scientific discovery. However, accurate simulations are often slow to execute, which limits their applicability to extensive parameter exploration, large-scale data analysis, and uncertainty quantification. A promising route to accelerate simulations by building fast emulators with machine learning requires large training datasets, which can be prohibitively expensive to obtain with slow simulations. Here we present a method based on neural architecture search to build accurate emulators even with a limited number of training data. The method successfully accelerates simulations by up to 2 billion times in 10 scientific cases including astrophysics, climate science, biogeochemistry, high energy density physics, fusion energy, and seismology, using the same super-architecture, algorithm, and hyperparameters. Our approach also inherently provides emulator uncertainty estimation, adding further confidence in their use. We anticipate this work will accelerate research involving expensive simulations, allow more extensive parameters exploration, and enable new, previously unfeasible computational discovery.

The development of high-brightness free-electron lasers (FEL) has revolutionised our ability to create and study matter in the high-energy-density (HED) regime. Current diagnostic techniques have been very successful in yielding information on fundamental thermodynamic plasma properties, but provide only limited or indirect information on the detailed quantum structure of these systems, and on how it is affected by ionization dynamics. Here we show how the electronic structure of solid-density nickel, heated to temperatures of 10's of eV on femtosecond timescales, can be studied by resonant (Raman) inelastic x-ray scattering (RIXS) using the Linac Coherent Light Source FEL. We present single-shot measurements of the valence density of states in the x-ray-heated transient system, and extract simultaneously electron temperatures, ionization, and ionization potential energies. The RIXS spectrum provides a wealth of information on the valence structure of the HED system that goes beyond what can be extracted from x-ray absorption or emission spectroscopy alone.

Proton radiography is a technique in high energy density science to diagnose magnetic and/or electric fields in a plasma by firing a proton beam and detecting its modulated intensity profile on a screen. Current approaches to retrieve the integrated field from the modulated intensity profile require the unmodulated beam intensity profile before the interaction, which is rarely available experimentally due to shot-to-shot variability. In this paper, we present a statistical method to retrieve the integrated field without needing to know the exact source profile. We apply our method to experimental data, showing the robustness of our approach. Our proposed technique allows not only for the retrieval of the path-integrated fields, but also of the statistical properties of the fields.

The free-free opacity in dense systems is a property that both tests our fundamental understanding of correlated many-body systems, and is needed to understand the radiative properties of high energy-density plasmas. Despite its importance, predictive calculations of the free-free opacity remain challenging even in the condensed matter phase for simple metals. Here we show how the free-free opacity can be modelled at finite-temperatures via time-dependent density functional theory, and illustrate the importance of including local field corrections, core polarization and self-energy corrections. Our calculations for ground-state Al are shown to agree well with experimental opacity measurements performed on the Artemis laser facility across a wide range of x-ray to ultraviolet wavelengths. We extend our calculations across the melt to the warm-dense matter regime, and find good agreement with advanced plasma models based on inverse bremsstrahlung at temperatures above 10 eV.

Our understanding of physical systems generally depends on our ability to match complex computational modelling with measured experimental outcomes. However, simulations with large parameter spaces suffer from inverse problem instabilities, where similar simulated outputs can map back to very different sets of input parameters. While of fundamental importance, such instabilities are seldom resolved due to the intractably large number of simulations required to comprehensively explore parameter space. Here we show how Bayesian machine learning can be used to address inverse problem instabilities, and apply it to two popular experimental diagnostics in plasma physics. We find that the extraction of information from measurements simply on the basis of agreement with simulations is unreliable, and leads to a significant underestimation of uncertainties. We describe how to statistically quantify the effect of unstable inverse models, and describe an approach to experimental design that mitigates its impact.

X-ray emission spectroscopy is a well-established technique used to study continuum lowering in dense plasmas. It relies on accurate atomic physics models to robustly reproduce high-resolution emission spectra, and depends on our ability to identify spectroscopic signatures such as emission lines or ionization edges of individual charge states within the plasma. Here we describe a method that forgoes these requirements, enabling the validation of different continuum lowering models based solely on the total intensity of plasma emission in systems driven by narrow-bandwidth x-ray pulses across a range of wavelengths. The method is tested on published Al spectroscopy data and applied to the new case of solid-density partially-ionized Fe plasmas, where extracting ionization edges directly is precluded by the significant overlap of emission from a wide range of charge states.

We report new measurements of the acoustic excitation of an Al5056 superconductive bar when hit by an electron beam, in a previously unexplored temperature range, down to 0.35 K. These data, analyzed together with previous results of the RAP experiment obtained for T > 0.54 K, show a vibrational response enhanced by a factor 4.9 with respect to that measured in the normal state. This enhancement explains the anomalous large signals due to cosmic rays previously detected in the NAUTILUS gravitational wave detector., Comment: 28 pages, 13 figures

We present calculations of the free-free XUV opacity of warm, solid-density aluminum at photon energies between the plasma frequency at 15 eV and the L-edge at 73 eV, using both density functional theory combined with molecular dynamics and a semi-analytical model in the RPA framework with the inclusion of local field corrections. As the temperature is increased from room temperature to 10 eV, with the ion and electron temperatures equal, we calculate an increase in the opacity in the range over which the degree of ionization is constant. The effect is less pronounced if only the electron temperature is allowed to increase. The physical significance of these increases is discussed in terms of intense XUV-laser matter interactions on both femtosecond and picosecond time-scales., Comment: 4 pages, 3 figures

Peak amplitude measurements of the fundamental mode of oscillation of a suspended aluminum alloy bar hit by an electron beam show that the amplitude is enhanced by a factor ~3.5 when the material is in the superconducting state. This result is consistent with the cosmic ray observations made by the resonant gravitational wave detector NAUTILUS, made of the same alloy, when operated in the superconducting state. A comparison of the experimental data with the predictions of the model describing the underlying physical process is also presented.