Plasma Physics (physics.plasm-ph)

Abstract: We study with a 3D PIC simulation discontinuities between an electron-positron pair plasma and magnetized electrons and protons. A pair plasma is injected at one simulation boundary with a speed 0.6$c$ along its normal. It expands into an electron-proton plasma and a magnetic field that points orthogonally to the injection direction. Diamagnetic currents expel the magnetic field from within the pair plasma and pile it up in front of it. It pushes electrons, which induces an electric field pulse ahead of the magnetic one. This initial electromagnetic pulse (EMP) confines the pair plasma magnetically and accelerates protons electrically. The fast flow of the injected pair plasma across the protons behind the initial EMP triggers the filamentation instability. Some electrons and positrons cross the injection boundary and build up a second EMP. Electron-cyclotron drift instabilities perturb the plasma ahead of both EMPs seeding a Rayleigh-Taylor-type instability. Despite equally strong perturbations ahead of both EMPs, the second EMP is much more stable than the initial one. We attribute the rapid collapse of the initial EMP to the filamentation instability, which perturbed the plasma behind it. The Rayleigh-Taylor-type instability transforms the planar EMPs into transition layers, in which magnetic flux ropes and electrostatic forces due to uneven numbers of electrons and positrons slow down and compress the pair plasma and accelerate protons. In our simulation, the expansion speed of the pair cloud decreased by about an order of magnitude and its density increased by the same factor. Its small thickness implies that it is capable of separating a relativistic pair outflow from an electron-proton plasma, which is essential for collimating relativistic jets of pair plasma in collisionless astrophysical plasma.

Abstract: The properties of hydrogen at warm dense matter (WDM) conditions are of high importance for the understanding of astrophysical objects and technological applications such as inertial confinement fusion. In this work, we present extensive new \emph{ab initio} path integral Monte Carlo (PIMC) results for the electronic properties in the Coulomb potential of a fixed ionic configuration. This gives us new insights into the complex interplay between the electronic localization around the protons with their density response to an external harmonic perturbation. We find qualitative agreement between our simulation data and a heuristic model based on the assumption of a local uniform electron gas model, but important trends are not captured by this simplification. In addition to being interesting in their own right, we are convinced that our results will be of high value for future projects, such as the rigorous benchmarking of approximate theories for the simulation of WDM, most notably density functional theory.

Abstract: Magnetic holes (MHs) are coherent structures associated with strong magnetic field depressions in magnetized plasmas. They are observed in many astrophysical environments at a wide range of scales but their origin is still under debate. In this work we investigate the formation of sub-ion scale MHs using a fully kinetic 2D simulation of plasma turbulence initialized with parameters typical of the Earth's magnetosheath. Our analysis shows that the turbulence is capable of generating sub-ion scale MHs from large scale fluctuations via the following mechanism: first, the nonlinear large scale dynamics spontaneously leads to the development of thin and elongated electron velocity shears; these structures then become unstable to the electron Kelvin-Helmholtz instability and break up into small scale electron vortices; the electric current carried by these vortices locally reduces the magnetic field, inducing the formation of sub-ion scale MHs. The MHs thus produced exhibit features consistent with satellite observations and with previous numerical studies. We finally discuss the kinetic properties of the observed sub-ion scale MHs, showing that they are characterized by complex non-Maxwellian electron velocity distributions exhibiting anisotropic and agyrotropic features.

Abstract: A novel quantum algorithm for solving the Boltzmann-Maxwell equations of the 6D collisionless plasma is proposed. The equation describes the kinetic behavior of plasma particles in electromagnetic fields and is known for the classical first-principles equations in various domains, from space to laboratory plasmas. We have constructed a quantum algorithm for a future large-scale quantum computer to accelerate its costly computation. This algorithm consists mainly of two routines: the Boltzmann solver and the Maxwell solver. Quantum algorithms undertake these dual procedures, while classical algorithms facilitate their interplay. Each solver has a similar structure consisting of three steps: Encoding, Propagation, and Integration. We conducted a preliminary implementation of the quantum algorithm and performed a parallel validation against a comparable classical approach. IBM Qiskit was used to implement all quantum circuits.