Laboratory astroparticle physics

Abstract:

We present numerical simulations used to interpret laser-driven plasma experiments at the GSI Helmholtz Centre for Heavy Ion Research. The mechanisms by which non-thermal particles are accelerated, in astrophysical environments e.g., the solar wind, supernova remnants, and gamma ray bursts, is a topic of intense study. When shocks are present the primary acceleration mechanism is believed to be first-order Fermi, which accelerates particles as they cross a shock. Second-order Fermi acceleration can also contribute, utilizing magnetic mirrors for particle energization. Despite this mechanism being less efficient, the ubiquity of magnetized turbulence in the universe necessitates its consideration. Another acceleration mechanism is the lower-hybrid drift instability, arising from gradients of both density and magnetic field, which produce lower-hybrid waves with an electric field which energizes particles as they cross these waves. With the combination of high-powered laser systems and particle accelerators it is possible to study the mechanisms behind cosmic-ray acceleration in the laboratory. In this work, we combine experimental results and high-fidelity threedimensional simulations to estimate the efficiency of ion acceleration in a weakly magnetized interaction region. We validate the FLASH MHD code with experimental results and use OSIRIS particle-in-cell (PIC) code to verify the initial formation of the interaction region, showing good agreement between codes and experimental results. We find that the plasma conditions in the experiment are conducive to the lower-hybrid drift instability, yielding an increase in energy ∆E of ∼ 264 keV for 242 MeV calcium ions.

Abstract:

We present two methods for computing the dynamic structure factor for warm dense hydrogen without invoking either the Born-Oppenheimer approximation or the Chihara decomposition, by employing a wave-packet description that resolves the electron dynamics during ion evolution. First, a semiclassical method is discussed, which is corrected based on known quantum constraints, and second, a direct computation of the density response function within the molecular dynamics. The wave-packet models are compared to PIMC and DFT-MD for the static and low-frequency behavior. For the high-frequency behavior the models recover the expected behavior in the limits of small and large momentum transfers and show the characteristic flattening of the plasmon dispersion for intermediate momentum transfers due to interactions, in agreement with commonly used models for x-ray Thomson scattering. By modeling the electrons and ions on an equal footing, both the ion and free electron part of the spectrum can now be treated within a single framework where we simultaneously resolve the ion-acoustic and plasmon mode, with a self-consistent description of collisions and screening.

Abstract:

We extract electron transport properties from atomistic simulations of a two-component plasma by mapping the long-wavelength behaviour to a two-fluid model. The mapping procedure is performed via Markov Chain Monte Carlo sampling over multiple spectra simultaneously. The free-electron dynamic structure factor and its properties have been investigated in the hydrodynamic formulation to justify its application to the long-wavelength behaviour of warm dense matter. We have applied this method to warm dense hydrogen modelled with wave packet molecular dynamics and showed that the inferred electron transport properties are in agreement with a variety of reference calculations, except for the electron viscosity, where a substantive decrease is observed when compared to classical models.