In this study we investigate jet interaction at a microscopic level in a cosmological environment, which responds to a key open question in the study of relativistic jets. Using small simulation systems during prior research, we initially studied the evolution of both electron-proton and electron-positron relativistic jets containing helical magnetic fields, by focusing on their interactions with an ambient plasma. Here, using larger jet radii, we have performed simulations of global jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). We found that the evolution of global jets strongly depends on the size of the jet radius. For example, phase bunching of jet electrons, in particular in the electron-proton jet, is mixed with larger jet radius due to the more complicated structures of magnetic fields with excited kinetic instabilities. In our simulation study these kinetic instabilities lead to new types of instabilities in global jets. In the electron-proton jet simulation a modified recollimation occurs and jet electrons are strongly perturbed. In the electron-positron jet simulation mixed kinetic instabilities occur at early times followed by a turbulence-like structure. Simulations using much larger (and longer) systems are further required in order to thoroughly investigate the evolution of global jets containing helical magnetic fields.

We study the interaction of relativistic jets with their environment, using 3-dimensional relativistic particle-in-cell simulations for two cases of jet composition: (i) electron-proton ($e^{-}-p^{+}$) and (ii) electron-positron ($e^{\pm}$) plasmas containing helical magnetic fields. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability and the Mushroom instability. We have found that these kinetic instabilities are suppressed and new types of instabilities can grow. For the $e^{-}-p^{+}$ jet, a recollimation-like instability occurs and jet electrons are strongly perturbed, whereas for the $e^{\pm}$ jet, a recollimation-like instability occurs at early times followed by kinetic instability and the general structure is similar to a simulation without a helical magnetic field. We plan to perform further simulations using much larger systems to confirm these new findings.

In the study of relativistic jets one of the key open questions is their interaction with the environment on the microscopic level. Here, we study the initial evolution of both electron$-$proton ($e^{-}-p^{+}$) and electron$-$positron ($e^{\pm}$) relativistic jets containing helical magnetic fields, focusing on their interaction with an ambient plasma. We have performed simulations of "global" jets containing helical magnetic fields in order to examine how helical magnetic fields affect kinetic instabilities such as the Weibel instability, the kinetic Kelvin-Helmholtz instability (kKHI) and the Mushroom instability (MI). In our initial simulation study these kinetic instabilities are suppressed and new types of instabilities can grow. In the $e^{-}-p^{+}$ jet simulation a recollimation-like instability occurs and jet electrons are strongly perturbed. In the $e^{\pm}$ jet simulation a recollimation-like instability occurs at early times followed by a kinetic instability and the general structure is similar to a simulation without helical magnetic field. Simulations using much larger systems are required in order to thoroughly follow the evolution of global jets containing helical magnetic fields.

One of the key open questions in the study of relativistic jets is their interaction with the environment. Here, we study the initial evolution of both electron-proton and electron-positron relativistic jets, focusing on their lateral interaction with the ambient plasma. We trace the generation and evolution of the toroidal magnetic field generated by both kinetic Kelvin-Helmholtz (kKH) and Mushroom instabilities (MI). This magnetic field collimates the jet. We show that in electron-proton jet, electrons are perpendicularly accelerated with jet collimation. The magnetic polarity switches from the clockwise to anti-clockwise in the middle of jet, as the instabilities weaken. For the electron-positron jet, we find strong mixture of electron-positron with the ambient plasma, that results in the creation of a bow shock. Merger of magnetic field current filaments generate density bumps which initiate a forward shock. The strong mixing between jet and ambient particles prevents full development of the jet on the studied scale. Our results therefore provide a direct evidence for both jet collimation and particle acceleration in the created bow shock. Differences in the magnetic field structures generated by electron-proton and electron-positron jets may contribute to observable differences in the polarized properties of emission by electrons.

We devise and explore an iterative optimization procedure for controlling particle populations in particle-in-cell (PIC) codes via merging and splitting of computational macro-particles. Our approach, is to compute an optimal representation of the global particle phase space structure while decreasing or increasing the entire particle population, based on k-means clustering of the data. In essence the procedure amounts to merging or splitting particles by statistical means, throughout the entire simulation volume in question, while minimizing a 6-dimensional total distance measure to preserve the physics. Particle merging is by far the most demanding procedure when considering conservation laws of physics; it amounts to lossy compression of particle phase space data. We demonstrate that our k-means approach conserves energy and momentum to high accuracy, even for high compression ratios, $\mathcal{R} \approx 3$ --- \emph{i.e.}, $N_{f} \lesssim 0.33N_{i}$. Interestingly, we find that an accurate particle splitting step can be performed using k-means as well; this from an argument of symmetry. The split solution, using k-means, places splitted particles optimally, to obtain maximal spanning on the phase space manifold. Implementation and testing is done using an electromagnetic PIC code, the \ppcode. Nonetheless, the k-means framework is general; it is not limited to Vlasov-Maxwell type PIC codes. We discuss advantages and drawbacks of this optimal phase space reconstruction.

We have investigated particle acceleration and emission from shocks and shear flows associated with an unmagnetized relativistic jet plasma propagating into an unmagnetized ambient plasma. Strong electro-magnetic fields are generated in the jet shock via the filamentation (Weibel) instability. Shock field strength and structure depend on plasma composition (($e^{\pm}$ or $e^-$- $p^+$ plasmas) and Lorentz factor. In the velocity shear between jet and ambient plasmas, strong AC ($e^{\pm}$ plasmas) or DC ($e^-$- $p^+$ plasmas) magnetic fields are generated via the kinetic Kelvin-Helmholtz instability (kKHI), and the magnetic field structure also depends on the jet Lorentz factor. We have calculated, self-consistently, the radiation from electrons accelerated in shock generated magnetic fields. The spectra depend on the jet's initial Lorentz factor and temperature via the resulting particle acceleration and magnetic field generation. Our ongoing "Global" jet simulations containing shocks and velocity shears will provide us with the ability to calculate and model the complex time evolution and/or spectral structure observed from gamma-ray bursts, AGN jets, and supernova remnants.

The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime.

We have investigated generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of KKHI of our jet-sheath configuration is slightly different even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field $E_{\rm z}$ and the magnetic field $B_{\rm y}$. After the $B_{\rm y}$ component is excited, an induced electric field $E_{\rm x}$ becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios $m_{\rm i}/m_{\rm e} = 1836$ and $m_{\rm i}/m_{\rm e} = 20$ are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case ($\gamma_{\rm j} = 1.5$) is larger than for a relativistic jet case ($\gamma_{\rm j} = 15$).

We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Strong magnetic fields generated in the trailing jet shock lead to transverse deflection and acceleration of the electrons. We have self-consistently calculated the radiation from the electrons accelerated in the turbulent magnetic fields. We find that the synthetic spectra depend on the bulk Lorentz factor of the jet, the jet temperature, and the strength of the magnetic fields generated in the shock. We have also begun study of electron acceleration in the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Our calculated spectra should lead to a better understanding of the complex time evolution and/or spectral structure from gamma-ray bursts, relativistic jets, and supernova remnants.

We present the Photon-Plasma code, a modern high order charge conserving particle-in-cell code for simulating relativistic plasmas. The code is using a high order implicit field solver and a novel high order charge conserving interpolation scheme for particle-to-cell interpolation and charge deposition. It includes powerful diagnostics tools with on-the-fly particle tracking, synthetic spectra integration, 2D volume slicing, and a new method to correctly account for radiative cooling in the simulations. A robust technique for imposing (time-dependent) particle and field fluxes on the boundaries is also presented. Using a hybrid OpenMP and MPI approach the code scales efficiently from 8 to more than 250.000 cores with almost linear weak scaling on a range of architectures. The code is tested with the classical benchmarks particle heating, cold beam instability, and two-stream instability. We also present particle-in-cell simulations of the Kelvin-Helmholtz instability, and new results on radiative collisionless shocks.

Spontaneous rapid growth of strong magnetic fields is rather ubiquitous in high-energy density environments ranging from astrophysical sources (e.g., gamma-ray bursts and relativistic shocks), to reconnection, to laser-plasma interaction laboratory experiments, where they are produced by kinetic streaming instabilities of the Weibel type. Relativistic electrons propagating through these sub-Larmor-scale magnetic fields radiate in the jitter regime, in which the anisotropy of the magnetic fields and the particle distribution have a strong effect on the produced radiation. Here we develop the general theory of jitter radiation, which includes (i) anisotropic magnetic fields and electron velocity distributions, (ii) the effects of trapped electrons and (iii) extends the description to large deflection angles of radiating particles thus establishing a cross-over between the classical jitter and synchrotron regimes. Our results are in remarkable agreement with the radiation spectra obtained from particle-in-cell simulations of the classical Weibel instability. Particularly interesting is the onset of the field growth, when the transient hard synchrotron-violating spectra are common as a result of the dominant role of the trapped population. This effect can serve as a distinct observational signature of the violent field growth in astrophysical sources and lab experiments. It is also interesting that a system with small-scale fields tends to evolve toward the small-angle jitter regime, which can, under certain conditions, dominate the overall emission of a source.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system. We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system generates a jitter/synchrotron spectrum.

Radiation from many astrophysical sources, e.g. gamma-ray bursts and active galactic nuclei, is believed to arise from relativistically shocked collisionless plasmas. Such sources often exhibit highly transient spectra evolving rapidly, compared with source lifetimes. Radiation emitted from these sources is typically associated with non-linear plasma physics, complex field topologies and non-thermal particle distributions. In such circumstances a standard synchrotron paradigm may fail to produce accurate conclusions regarding the underlying physics. Simulating spectral emission and spectral evolution numerically in various relativistic shock scenarios is then the only viable method to determine the detailed physical origin of the emitted spectra. In this Letter we present synthetic radiation spectra representing the early stage development of the filamentation (streaming) instability of an initially unmagnetized plasma, which is relevant for both collisionless shock formation and reconnection dynamics in relativistic astrophysical outflows, as well as for laboratory astrophysics experiments. Results were obtained using a highly efficient "in situ" diagnostics method, based on detailed particle-in-cell modeling of collisionless plasmas. The synthetic spectra obtained here are compared with those predicted by a semi-analytical model for jitter radiation from the filamentation instability, the latter including self-consistent generated field topologies and particle distributions obtained from the simulations reported upon here. Spectra exhibit dependence on the presence - or absence - of an inert plasma constituent, when comparing baryonic plasmas (i.e. containing protons) with pair plasmas. The results also illustrate that considerable care should be taken when using lower-dimensional models to obtain information about the astrophysical phenomena generating observed spectra.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. We have also performed simulations with electron-ion jets. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability for electron-positron jets and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks for both cases. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value for pair plasmas. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system with two different cases for Lorentz factors (15 and 100). We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we have investigated long-term particle acceleration associated with an relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations have been performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. The acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to the afterglow emission. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and microquasars commonly exhibit power-law emission spectra. Recent PIC simulations of relativistic electron-ion (or electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In collisionless, relativistic shocks, particle (electron, positron, and ion) acceleration is due to plasma waves and their associated instabilities (e.g., the Weibel (filamentation) instability) created in the shock region. The simulations show that the Weibel instability is responsible for generating and amplifying highly non-uniform, small-scale magnetic fields. These fields contribute to the electron's transverse deflection behind the jet head. The resulting "jitter" radiation from deflected electrons has different properties compared to synchrotron radiation, which assumes a uniform magnetic field. Jitter radiation may be important for understanding the complex time evolution and/or spectra in gamma-ray bursts, relativistic jets in general, and supernova remnants.

In continuation of previous work, numerical results are presented, concerning relativistically counter-streaming plasmas. Here, the relativistic mixed mode instability evolves through, and beyond, the linear saturation -- well into the nonlinear regime. Besides confirming earlier findings, that wave power initially peaks on the mixed mode branch, it is observed that, during late time evolution wave power is transferred to other wave numbers. It is argued that the isotropization of power in wavenumber space may be a consequence of weak turbulence. Further, some modifications to the ideal weak turbulence limit is observed. Development of almost isotropic predominantly electrostatic -- partially electromagnetic -- turbulent spectra holds relevance when considering the spectral emission signatures of the plasma, namely bremsstrahlung, respectively magneto-bremsstrahlung (synchrotron radiation and jitter radiation) from relativistic shocks in astrophysical jets and shocks from gamma-ray bursts and active galactic nuclei.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the presence of relativistic jets, instabilities such as the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability create collisionless shocks, which are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons in small-scale magnetic fields has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation, a case of diffusive synchrotron radiation, may be important to understand the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Modeling plasma physical processes in astrophysical context demands for both detailed kinetics and large scale development of the electromagnetic field densities. We present a new framework for modeling plasma physics of hot tenuous plasmas by a two-split scheme, in which the large scale fields are modeled by means of a particle-in-cell (PIC) code, and in which binary collision processes and single-particle processes are modeled through a Monte-Carlo approach. Our novel simulation tool -- the PhotonPlasma code -- is a unique hybrid model; it combines a highly parallelized (Vlasov) particle-in-cell approach with continuous weighting of particles and a sub-Debye Monte-Carlo binary particle interaction framework. As an illustration of the capabilities we present results from a numerical study of Gamma-Ray Burst - Circumburst Medium interaction and plasma preconditioning via Compton scattering. We argue that important microphysical processes can only viably be investigated by means of hybrid codes such as the PhotonPlasma code. Our first results from 3D simulations with this new simulation tool suggest that magnetic fields and plasma filaments are created in the wakefield of prompt gamma-ray bursts. Furthermore, the photon flux density gradient impacts on particle acceleration in the burst head and wakefield. We discuss some possible implications of the circumburst medium being preconditioned for a trailing afterglow shock front. We also discuss important improvements for future studies of GRB wakefields processes, using the PhotonPlasma code.

We present a numerical study of Gamma-Ray Burst - Circumburst Medium interaction and plasma preconditioning via Compton scattering. The simulation tool employed is a unique hybrid model; it combines a highly parallelized (Vlasov) particle-in-cell approach with continuous weighting of particles and a sub-Debye Monte-Carlo binary particle interaction framework. These first results from 3D simulations with this new simulation tool, the PhotonPlasma code, suggests that magnetic fields and plasma density filaments are created in the wakefield of prompt gamma-ray bursts, and that the photon flux density gradient has a significant impact on particle acceleration in the burst head and wakefield. We discuss some possible implications of the circumburst medium preconditioning for the trailing afterglow, and also discuss which additional processes will be needed to improve future studies within this unique and powerful simulation framework.

Gamma ray bursts are among the most energetic events in the known universe. A highly relativistic fireball is ejected. In most cases the burst itself is followed by an afterglow, emitted under deceleration as the fireball plunges through the circum-stellar media. To interpret the observations of the afterglow emission, two physical aspects need to be understood: 1) The origin and nature of the magnetic field in the fireball and 2) the particle velocity distribution function behind the shock. Both are necessary in existing afterglow models to account for what is believed to be synchrotron radiation. To answer these questions, we need to understand the microphysics at play in collisionless shocks. Using 3D particle-in-cell simulations we can gain insight in the microphysical processes that take place in such shocks. We discuss the results of such computer experiments. It is shown how a Weibel-like two-stream plasma instability is able to create a strong transverse intermittent magnetic field and points to a connected mechanism for in situ particle acceleration in the shock region.

The outflows from gamma ray bursts, active galactic nuclei and relativistic jets in general interact with the surrounding media through collisionless shocks. With three dimensional relativistic particle-in-cell simulations we investigate such shocks. The results from these experiments show that small--scale magnetic filaments with strengths of up to percents of equipartition are generated and that electrons are accelerated to power law distributions N(E)~E^{-p} in the vicinity of the filaments through a new acceleration mechanism. The acceleration is locally confined, instantaneous and differs from recursive acceleration processes such as Fermi acceleration. We find that the proposed acceleration mechanism competes with thermalization and becomes important at high Lorentz factors.

Collisionless plasma shock theory, which applies for example to the afterglow of gamma ray bursts, still contains key issues that are poorly understood. In this paper we study charged particle dynamics in a highly relativistic collisionless shock numerically using ~10^9 particles. We find a power law distribution of accelerated electrons, which upon detailed investigation turns out to originate from an acceleration mechanism that is decidedly different from Fermi acceleration. Electrons are accelerated by strong filamentation instabilities in the shocked interpenetrating plasmas and coincide spatially with the power law distributed current filamentary structures. These structures are an inevitable consequence of the now well established Weibel-like two-stream instability that operates in relativistic collisionless shocks. The electrons are accelerated and decelerated instantaneously and locally; a scenery that differs qualitatively from recursive acceleration mechanisms such as Fermi acceleration. The slopes of the electron distribution power laws are in concordance with the particle power law spectra inferred from observed afterglow synchrotron radiation in gamma ray bursts, and the mechanism can possibly explain more generally the origin of non-thermal radiation from shocked inter- and circum-stellar regions and from relativistic jets.

Using a three dimensional relativistic particle-in-cell code we have performed numerical experiments of plasma shells colliding at relativistic velocities. Such scenarios are found in many astrophysical objects e.g. the relativistic outflow from gamma ray bursts, active galactic nuclei jets and supernova remnants. We show how a Weibel-like two-stream instability is capable of generating small-scale magnetic filaments with strength up to percents of equipartition. Such field topology is ideal for the generation of jitter radiation as opposed to synchrotron radiation. We also explain how the field generating mechanism involves acceleration of electrons to power law distributions (N(E)~E^(-p))through a non-Fermi acceleration mechanism. The results add to our understanding of collisionless shocks.

We present results from three-dimensional particle simulations of collisionless shocks with relativistic counter-streaming ion-electron plasmas. Particles are followed over many skin depths downstream of the shock. Open boundaries allow the experiments to be continued for several particle crossing times. The experiments confirm the generation of strong magnetic and electric fields by a Weibel-like kinetic streaming instability, and demonstrate that the electromagnetic fields propagate far downstream of the shock. The magnetic fields are predominantly transversal, and are associated with merging ion current channels. The total magnetic energy grows as the ion channels merge, and as the magnetic field patterns propagate down stream. The electron populations are quickly thermalized, while the ion populations retain distinct bulk speeds in shielded ion channels and thermalize much more slowly. These results may help explain the origin of the magnetic fields responsible for afterglow synchrotron/jitter radiation from Gamma-Ray Bursts.

We present numerical results from plasma particle simulations of collisionless shocks and ultra-relativistic counter-streaming plasmas. We demonstrate how the field-particle interactions lead to particle acceleration behind the shock-front. Further, we demonstrate how ultra relativistic counter-streaming plasmas create large scale patchy magnetic field structures and that these field structures propagate down-stream of the shock front. These results may help explain the origin of the magnetic fields and accelerated electrons responsible for afterglow synchrotron radiation from gamma ray bursts.