PIC Simulations of the Temperature Anisotropy-Driven Weibel Instability: Analyzing the perpendicular mode

An instability driven by the thermal anisotropy of a single electron species is investigated in a 2D particle-in-cell (PIC) simulation. This instability is the one considered by Weibel and it differs from the beam driven filamentation instability. A comparison of the simulation results with analytic theory provides similar exponential growth rates of the magnetic field during the linear growth phase of the instability. We observe in accordance with previous works the growth of electric fields during the saturation phase of the instability. Some components of this electric field are not accounted for by the linearized theory. A single-fluid-based theory is used to determine the source of this nonlinear electric field. It is demonstrated that the magnetic stress tensor, which vanishes in a 1D geometry, is more important in this 2-dimensional model used here. The electric field grows to an amplitude, which yields a force on the electrons that is comparable to the magnetic one. The peak energy density of each magnetic field component in the simulation plane agrees with previous estimates. Eddy currents develop, which let the amplitude of the third magnetic field component grow, which is not observed in a 1D simulation.Comment: accepted by Plasma Physics and Controlled Fusio

Electromagnetic field generation in the downstream of electrostatic shocks due to electron trapping

A new magnetic field generation mechanism in electrostatic shocks is found, which can produce fields with magnetic energy density as high as 0.01 of the kinetic energy density of the flows on time scales $\tilde \, 10^4 \, {\omega}_{pe}^{-1}$. Electron trapping during the shock formation process creates a strong temperature anisotropy in the distribution function, giving rise to the pure Weibel instability. The generated magnetic field is well-confined to the downstream region of the electrostatic shock. The shock formation process is not modified and the features of the shock front responsible for ion acceleration, which are currently probed in laser-plasma laboratory experiments, are maintained. However, such a strong magnetic field determines the particle trajectories downstream and has the potential to modify the signatures of the collisionless shock

Physics of collisionless shocks - theory and simulation

Collisionless shocks occur in various fields of physics. In the context of space and astrophysics they have been investigated for many decades. However, a thorough understanding of shock formation and particle acceleration is still missing. Collisionless shocks can be distinguished into electromagnetic and electrostatic shocks. Electromagnetic shocks are of importance mainly in astrophysical environments and they are mediated by the Weibel or filamentation instability. In such shocks, charged particles gain energy by diffusive shock acceleration. Electrostatic shocks are characterized by a strong electrostatic field, which leads to electron trapping. Ions are accelerated by reflection from the electrostatic potential. Shock formation and particle acceleration will be discussed in theory and simulations

The impact of kinetic effects on the properties of relativistic electron-positron shocks

We assess the impact of non-thermally shock-accelerated particles on the magnetohydrodynamic (MHD) jump conditions of relativistic shocks. The adiabatic constant is calculated directly from first principle particle-in-cell simulation data, enabling a semi-kinetic approach to improve the standard fluid model and allowing for an identification of the key parameters that define the shock structure. We find that the evolving upstream parameters have a stronger impact than the corrections due to non-thermal particles. We find that the decrease of the upstream bulk speed yields deviations from the standard MHD model up to 10%. Furthermore, we obtain a quantitative definition of the shock transition region from our analysis. For Weibel-mediated shocks the inclusion of a magnetic field in the MHD conservation equations is addressed for the first time

Exploring the nature of collisionless shocks under laboratory conditions

Collisionless shocks are pervasive in astrophysics and they are critical to understand cosmic ray acceleration. Laboratory experiments with intense lasers are now opening the way to explore and characterise the underlying microphysics, which determine the acceleration process of collisionless shocks. We determine the shock character - electrostatic or electromagnetic - based on the stability of electrostatic shocks to transverse electromagnetic fluctuations as a function of the electron temperature and flow velocity of the plasma components, and we compare the analytical model with particle-in-cell simulations. By making the connection with the laser parameters driving the plasma flows, we demonstrate that shocks with different and distinct underlying microphysics can be explored in the laboratory with state-of-the-art laser systems

Spatial-temporal evolution of the current filamentation instability

The spatial-temporal evolution of the purely transverse current filamentation instability is analyzed by deriving a single partial differential equation for the instability and obtaining the analytical solutions for the spatially and temporally growing current filament mode. When the beam front always encounters fresh plasma, our analysis shows that the instability grows spatially from the beam front to the back up to a certain critical beam length; then the instability acquires a purely temporal growth. This critical beam length increases linearly with time and in the non-relativistic regime it is proportional to the beam velocity. In the relativistic regime the critical length is inversely proportional to the cube of the beam Lorentz factor $\gamma_{0b}$. Thus, in the ultra-relativistic regime the instability immediately acquires a purely temporal growth all over the beam. The analytical results are in good agreement with multidimensional particle-in-cell simulations performed with OSIRIS. Relevance of current study to recent and future experiments on fireball beams is also addressed

Shock formation in electron-ion plasmas: mechanism and timing

We analyse the full shock formation process in electron-ion plasmas in theory and simulations. It is accepted that electromagnetic shocks in initially unmagnetised relativistic plasmas are triggered by the filamentation instability. However, the transition from the first unstable phase to the quasi-steady shock is still missing. We derive a theoretical model for the shock formation time, taking into account the filament merging in the non-linear phase of the filamentation instability. This process is much slower than in electron-positron pair shocks, so that the shock formation is longer by a factor proportional to sqrt(m_i/m_e) ln(m_i/m_e)

The filamentation instability driven by warm electron beams: Statistics and electric field generation

The filamentation instability of counterpropagating symmetric beams of electrons is examined with 1D and 2D particle-in-cell (PIC) simulations, which are oriented orthogonally to the beam velocity vector. The beams are uniform, warm and their relative speed is mildly relativistic. The dynamics of the filaments is examined in 2D and it is confirmed that their characteristic size increases linearly in time. Currents orthogonal to the beam velocity vector are driven through the magnetic and electric fields in the simulation plane. The fields are tied to the filament boundaries and the scale size of the flow-aligned and the perpendicular currents are thus equal. It is confirmed that the electrostatic and the magnetic forces are equally important, when the filamentation instability saturates in 1D. Their balance is apparently the saturation mechanism of the filamentation instability for our initial conditions. The electric force is relatively weaker but not negligible in the 2D simulation, where the electron temperature is set higher to reduce the computational cost. The magnetic pressure gradient is the principal source of the electrostatic field, when and after the instability saturates in the 1D simulation and in the 2D simulation.Comment: 10 pages, 6 figures, accepted by the Plasma Physics and Controlled Fusion (Special Issue EPS 2009

Acceleration in perpendicular relativistic shocks for plasmas consisting of leptons and hadrons

We investigate the acceleration of light particles in perpendicular shocks for plasmas consisting of a mixture of leptonic and hadronic particles. Starting from the full set of conservation equations for the mixed plasma constituents, we generalize the magneto-hydrodynamical jump conditions for a multi-component plasma, including information about the specific adiabatic constants for the different species. The impact of deviations from the standard model of an ideal gas is compared in theory and particle-in-cell simulations, showing that the standard-MHD model is a good approximation. The simulations of shocks in electron-positron-ion plasmas are for the first time multi-dimensional, transverse effects are small in this configuration and 1D simulations are a good representation if the initial magnetization is chosen high. 1D runs with a mass ratio of 1836 are performed, which identify the Larmor frequency \omega_{ci} as the dominant frequency that determines the shock physics in mixed component plasmas. The maximum energy in the non-thermal tail of the particle spectra evolves in time according to a power-law proportional to t^\alpha with \alpha in the range 1/3 < \alpha < 1, depending on the initial parameters. A connection is made with transport theoretical models by Drury (1983) and Gargate & Spitkovsky (2011), which predict an acceleration time proportional to \gamma and the theory for small wavelength scattering by Kirk & Reville (2010), which predicts a behavior rather as proportional to \gamma^2. Furthermore, we compare different magnetic field orientations with B_0 inside and out of the plane, observing qualitatively different particle spectra than in pure electron-ion shocks