Search

Energy dissipation in collisionless shocks is a key mechanism in various astrophysical environments. Its non-linear nature complicates analytical understanding and necessitate Particle-in-Cell (PIC) simulations. This study examines the impact of reducing the ion-to-electron mass ratio ($m_r$), to decrease computational cost, on energy partitioning in 1D3V (one spatial and three velocity-space dimensions) PIC simulations of strong, non-relativistic, parallel electron-ion collisionless shocks using the SHARP code. We compare simulations with a reduced mass ratio ($m_r = 100$) to those with a realistic mass ratio ($m_r = 1836$) for shocks with high ($\mathcal{M}_A = 21.3$) and low ($\mathcal{M}_A = 5.3$) Alfv$\acute{\text{e}}$n Mach numbers. Our findings show that the mass ratio significantly affects particle acceleration and thermal energy dissipation. At high $\mathcal{M}_A$, a reduced mass ratio leads to more efficient electron acceleration and an unrealistically high ion flux at higher momentum. At low $\mathcal{M}_A$, it causes complete suppression of electron acceleration, whereas the realistic mass ratio enables efficient electron acceleration. The reduced mass ratio also results in excessive electron heating and lower heating in downstream ions at both Mach numbers, with slightly more magnetic field amplification at low $\mathcal{M}_A$. Consequently, the electron-to-ion temperature ratio is high at low $\mathcal{M}_A$ due to reduced ion heating and remains high at high $\mathcal{M}_A$ due to increased electron heating. In contrast, simulations with the realistic $m_r$ show that the ion-to-electron temperature ratio is independent of the upstream magnetic field, a result not observed in reduced $m_r$ simulations., Comment: 9 pages, 6 figures; accepted for publication in ApJ Letters (submitted August 3, 2024)

Cosmic ray (CR) feedback is critical for galaxy formation as CRs drive galactic winds, regularize star formation in galaxies, and escape from active galactic nuclei to heat the cooling cores of galaxy clusters. The feedback strength of CRs depends on their coupling to the background plasma and, as such, on the effective CR transport speed. Traditionally, this has been hypothesized to depend on the balance between wave growth of CR-driven instabilities and their damping. Here, we study the physics of CR-driven instabilities from first principles, starting from a gyrotropic distribution of CR ions that stream along a background magnetic field. We develop a theory of the underlying processes that organize the particles' orbits and in particular their gyrophases, which provides an intuitive physical picture of (i) wave growth as the CR gyrophases start to bunch up lopsidedly towards the local wave magnetic field, (ii) instability saturation as a result of CRs overtaking the wave and damping its amplitude without isotropizing CRs in the wave frame, and (iii) CR back-reaction onto the unstable plasma waves as the CR gyrophases follow a pendulum motion around the wave magnetic field. Using our novel fluid-particle-in-cell code fluid-SHARP, we validate our theory on the evolution and excitation of individual unstable modes, such as forward and backward propagating Alfv\'en and whistler waves. We show that these kinetic simulations support our theoretical considerations, thus potentially foreshadowing a fundamental revision of the theory of CR transport in galaxies and galaxy clusters., Comment: Comments are welcome

We study the underlying physics of cosmic-ray (CR) driven instabilities that play a crucial role for CR transport across a wide range of scales, from interstellar to galaxy cluster environments. By examining the linear dispersion relation of CR-driven instabilities in a magnetised electron-ion background plasma, we establish that both, the intermediate and gyroscale instabilities have a resonant origin and show that these resonances can be understood via a simple graphical interpretation. These instabilities destabilise wave modes parallel to the large-scale background magnetic field at significantly distinct scales and with very different phase speeds. Furthermore, we show that approximating the electron-ion background plasma with either magnetohydrodynamics (MHD) or Hall-MHD fails to capture the fastest growing instability in the linear regime, namely the intermediate-scale instability. This finding highlights the importance of accurately characterising the background plasma for resolving the most unstable wave modes. Finally, we discuss the implications of the different phase speeds of unstable modes on particle-wave scattering. Further work is needed to investigate the relative importance of these two instabilities in the non-linear, saturated regime and to develop a physical understanding of the effective CR transport coefficients in large-scale CR hydrodynamics theories., Comment: 13 pages + references, 4 figures. To appear in JPP Letters

The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run-times and often to setups in less than three spatial dimensions. Traditionally, this is circumvented either via hybrid-PIC methods (adopting massless electrons) or via magneto-hydrodynamic-PIC methods (modelling the background plasma as a single charge-neutral magneto-hydrodynamical fluid). Because both methods preclude modelling important plasma-kinetic effects, we introduce a new fluid-PIC code that couples a fully explicit and charge-conservative multi-fluid solver to the PIC code SHARP through a current-coupling scheme and solve the full set of Maxwell's equations. This avoids simplifications typically adopted for Ohm's Law and enables us to fully resolve the electron temporal and spatial scales while retaining the versatility of initializing any number of ion, electron, or neutral species with arbitrary velocity distributions. The fluid solver includes closures emulating Landau damping so that we can account for this important kinetic process in our fluid species. Our fluid-PIC code is second-order accurate in space and time. The code is successfully validated against several test problems, including the stability and accuracy of shocks and the dispersion relation and damping rates of waves in unmagnetized and magnetized plasmas. It also matches growth rates and saturation levels of the gyro-scale and intermediate-scale instabilities driven by drifting charged particles in magnetized thermal background plasmas in comparison to linear theory and PIC simulations. This new fluid-SHARP code is specially designed for studying high-energy cosmic rays interacting with thermal plasmas over macroscopic timescales., Comment: 35 pages, 11 figures, submitted to JPP. Comments are welcome

Thermal electrons cannot directly participate in the process of diffusive acceleration at electron-ion shocks because their Larmor radii are smaller than the shock transition width: this is the well-known electron injection problem of diffusive shock acceleration. Instead, an efficient pre-acceleration process must exist that scatters electrons off of electromagnetic fluctuations on scales much shorter than the ion gyro radius. The recently found intermediate-scale instability provides a natural way to produce such fluctuations in parallel shocks. The instability drives comoving (with the upstream plasma) ion-cyclotron waves at the shock front and only operates when the drift speed is smaller than half of the electron Alfven speed. Here, we perform particle-in-cell simulations with the SHARP code to study the impact of this instability on electron acceleration at parallel non-relativistic, electron-ion shocks. To this end, we compare a shock simulation in which the intermediate-scale instability is expected to grow to simulations where it is suppressed. In particular, the simulation with an Alfvenic Mach number large enough to quench the intermediate instability shows a great reduction (by two orders of magnitude) of the electron acceleration efficiency. Moreover, the simulation with a reduced ion-to-electron mass ratio (where the intermediate instability is also suppressed) not only artificially precludes electron acceleration but also results in erroneous electron and ion heating in the downstream and shock transition regions. This finding opens up a promising route for a plasma physical understanding of diffusive shock acceleration of electrons, which necessarily requires realistic mass ratios in simulations of collisionless electron-ion shocks., Comment: 12 pages, 8 figures, 1 table, Accepted to ApJ

Cosmic ray (CR)-driven instabilities play a decisive role during particle acceleration at shocks and CR propagation in galaxies and galaxy clusters. These instabilities amplify magnetic fields and modulate CR transport so that the intrinsically collisionless CR population is tightly coupled to the thermal plasma and provides dynamical feedback. Here, we show that CRs with a finite pitch angle drive electromagnetic waves (along the background magnetic field) unstable on intermediate scales between the gyro-radii of CR ions and electrons as long as CRs are drifting with a velocity less than half of the Alfv\'en speed of electrons. By solving the linear dispersion relation, we show that this new instability typically grows faster by more than an order of magnitude in comparison to the commonly discussed resonant instability at the ion gyroscale. We find the growth rate for this intermediate-scale instability and identify the growing modes as background ion-cyclotron modes in the frame that is comoving with the CRs. We confirm the theoretical growth rate with a particle-in-cell (PIC) simulation and study the non-linear saturation of this instability. We identify three important astro-physical applications of this intermediate-scale instability, which is expected to 1. modulate CR transport and strengthen CR feedback in galaxies and galaxy clusters, 2. enable electron injection into the diffusive shock acceleration process, and 3. decelerate CR escape from the sites of particle acceleration which would generate gamma-ray halos surrounding CR sources such as supernova remnants., Comment: 16 pages, 13 figures, 2 tables. Accepted to ApJ. Comments are welcome. An animation for the growth of the new instability can be found at http://www.youtube.com/watch?v=yP2dBCPOQYc

We study the longitudinal stability of beam-plasma systems in the presence of a density inhomogeneity in the background plasma. Previous works have focused on the non-relativistic regime where hydrodynamical models are used to evolve pre-existing Langmuir waves within inhomogeneous background plasmas. Here, for the first time we study the problem with kinetic equations in a fully-relativistic way. We do not assume the existence of Langmuir waves, and we focus on the rate and the mechanism by which waves are excited in such systems from an initial perturbation. We derive the structure of the unstable modes and compute an analytical approximation for their growth rates. Our computation is limited to dilute and cold beams, and shows an excellent agreement with particle-in-cell simulations performed using the SHARP code. We show that, due to such an inhomogeneity, the virulent beam-plasma instabilities in the intergalactic medium are not suppressed but their counterparts in the solar wind can be suppressed as evidenced by propagating type-III solar radio bursts., Comment: 23 pages, 5 figures, accepted for publication in the Journal of Plasma Physics; comments are welcome

An intergalactic magnetic field stronger than $3\times10^{-13}$~G would explain the lack of a bright, extended degree-scale, GeV-energy inverse Compton component in the gamma-ray spectra of TeV-blazars. A robustly predicted consequence of the presence of such a field is the existence of degree-scale GeV-energy gamma-ray halos -- gamma-ray bow ties -- about TeV-bright active galactic nuclei, corresponding to more than half of all radio galaxies. However, the emitting regions of these halos are confined to and aligned with the direction of the relativistic jets associated with gamma-ray sources. Based on the orientation of radio jets, we align and stack corresponding degree-scale gamma-ray images of isolated Fanaroff-Riley class I and II objects and exclude the existence of these halos at overwhelming confidence, limiting the intergalactic field strength to $<10^{-15}$~G for large-scale fields and progressively larger in the diffusive regime when the correlation length of the field becomes small in comparison to 1 Mpc. When combined with prior limits on the strength of the intergalactic magnetic field, this excludes a purely magnetic explanation for the absence of halos. Thus, it requires the existence of novel physical processes that preempt the creation of halos, e.g., the presence of beam-plasma instabilities in the intergalactic medium or a drastic cutoff of the very high energy spectrum of these sources., Comment: Accepted to ApJ, 20 pages, 15 figures

We explore how inhomogeneity in the background plasma number density alters the growth of electrostatic unstable wavemodes of beam plasma systems. This is particularly interesting for blazar-driven beam-plasma instabilities, which may be suppressed by inhomogeneities in the intergalactic medium as was recently claimed in the literature. Using high resolution Particle-In-Cell simulations with the SHARP code, we show that the growth of the instability is local, i.e., regions with almost homogeneous background density will support the growth of the Langmuir waves, even when they are separated by strongly inhomogeneous regions, resulting in an overall slower growth of the instability. We also show that if the background density is continuously varying, the growth rate of the instability is lower; though in all cases, the system remains within the linear regime longer and the instability is not extinguished. In all cases, the beam loses approximately the same fraction of its initial kinetic energy in comparison to the uniform case at non-linear saturation. Thus, inhomogeneities in the intergalactic medium are unlikely to suppress the growth of blazar-driven beam-plasma instabilities., Comment: 10 pages, 6 figures, Accepted by ApJ, comments welcome

Many astrophysical plasmas are prone to beam-plasma instabilities. For relativistic and dilute beams, the {\it spectral} support of the beam-plasma instabilities is narrow, i.e., the linearly unstable modes that grow with rates comparable to the maximum growth rate occupy a narrow range of wave numbers. This places stringent requirements on the box-sizes when simulating the evolution of the instabilities. We identify the implied lower limits on the box size imposed by the longitudinal beam plasma instability, i.e., typically the most stringent condition required to correctly capture the linear evolution of the instabilities in multidimensional simulations. We find that sizes many orders of magnitude larger than the resonant wavelength are typically required. Using one-dimensional particle-in-cell simulations, we show that the failure to sufficiently resolve the spectral support of the longitudinal instability yields slower growth and lower levels of saturation, potentially leading to erroneous physical conclusion., Comment: 7 pages, 9 figures, accepted by ApJ

Numerical heating in particle-in-cell (PIC) codes currently precludes the accurate simulation of cold, relativistic plasma over long periods, severely limiting their applications in astrophysical environments. We present a spatially higher-order accurate relativistic PIC algorithm in one spatial dimension, which conserves charge and momentum exactly. We utilize the smoothness implied by the usage of higher-order interpolation functions to achieve a spatially higher-order accurate algorithm (up to fifth order). We validate our algorithm against several test problems -- thermal stability of stationary plasma, stability of linear plasma waves, and two-stream instability in the relativistic and non-relativistic regimes. Comparing our simulations to exact solutions of the dispersion relations, we demonstrate that SHARP can quantitatively reproduce important kinetic features of the linear regime. Our simulations have a superior ability to control energy non-conservation and avoid numerical heating in comparison to common second-order schemes. We provide a natural definition for convergence of a general PIC algorithm: the complement of physical modes captured by the simulation, i.e., those that lie above the Poisson noise, must grow commensurately with the resolution. This implies that it is necessary to simultaneously increase the number of particles per cell and decrease the cell size. We demonstrate that traditional ways for testing for convergence fail, leading to plateauing of the energy error. This new PIC code enables us to faithfully study the long-term evolution of plasma problems that require absolute control of the energy and momentum conservation., Comment: 26 pages, 19 figures, discussion about performance is added, published in ApJ

Many-degree-scale gamma-ray halos are expected to surround extragalactic high-energy gamma ray sources. These arise from the inverse Compton emission of an intergalactic population of relativistic electron/positron pairs generated by the annihilation of >100 GeV gamma rays on the extragalactic background light. These are typically anisotropic due to the jetted structure from which they originate or the presence of intergalactic magnetic fields. Here we propose a novel method for detecting these inverse-Compton gamma-ray halos based upon this anisotropic structure. Specifically, we show that by stacking suitably defined angular power spectra instead of images it is possible to robustly detect gamma-ray halos with existing Fermi Large Area Telescope (LAT) observations for a broad class of intergalactic magnetic fields. Importantly, these are largely insensitive to systematic uncertainties within the LAT instrumental response or associated with contaminating astronomical sources., Comment: Submitted to ApJ, 22 pages, 27 figures, 9 pages in appendix

Pair creation on the cosmic infrared background and subsequent inverse-Compton scattering on the CMB potentially reprocesses the TeV emission of blazars into faint GeV halos with structures sensitive to intergalactic magnetic fields (IGMF). Previous work has shown that these halos are then highly-anisotropic and extended. If the coherence length of the IGMF is greater than the inverse-Compton cooling length of the pairs, then the orientation of the gamma-ray halo will be correlated with the direction of the magnetic field which is unknown and expected to change for each source. In order to constructively add each source we then use angular power spectra which are insensitive to the jet orientation. By looking at known GeV blazars detected by Fermi, we exclude the existence of an IGMF with coherence lengths $>100$ Mpc at greater than $3.9\sigma$ with current-day strengths in the range $10^{-16}$ to $10^{-15}$ G, and at 2$\sigma$ from $10^{-17}$ to $10^{-14}$ G. This provides a direct measurement of the non-existence of gamma-ray halos, providing an important check on previous results., Comment: 12 pages, 14 figures, 1 appendix. Accepted to ApJ

The annihilation of TeV photons from extragalactic TeV sources and the extragalactic background light produces ultrarelativistic $e^{\pm}$ beams, which are subject to powerful plasma instabilities that sap their kinetic energy. Here we study the linear phase of the plasma instabilities that these pair beams drive. To this end, we calculate the linear growth rate of the beam plasma and oblique instability in the electrostatic approximation in both the reactive and kinetic regimes, assuming a Maxwell-J{\"u}ttner distribution for the pair beam. We reproduce the well-known reactive and kinetic growth rates for both the beam plasma and oblique mode. We demonstrate for the oblique instability that there is a broad spectrum of unstable modes that grow at the maximum rate for a wide range of beam temperatures and wave vector orientations relative to the beam. We also delineate the conditions for applicability for the reactive and kinetic regimes and find that the beam plasma mode transitions to the reactive regime at a lower Lorentz factor than the oblique mode due to a combination of their different scalings and the anisotropy of the velocity dispersions. Applying these results to the ultrarelativistic $e^{\pm}$ beams from TeV blazars, we confirm that these beams are unstable to both, the kinetic oblique mode and the reactive beam-plasma mode. These results are important in understanding how powerful plasma instabilities may sap the energy of the ultrarelativistic $e^{\pm}$ beams as they propagate through intergalactic space., Comment: 13 pages, 4 figure, accepted to ApJ

Extended inverse Compton halos are generally anticipated around extragalactic sources of gamma rays with energies above 100 GeV. These result from inverse Compton scattered cosmic microwave background photons by a population of high-energy electron/positron pairs produced by the annihilation of the high-energy gamma rays on the infrared background. Despite the observed attenuation of the high-energy gamma rays, the halo emission has yet to be directly detected. Here, we demonstrate that in most cases these halos are expected to be highly anisotropic, distributing the up-scattered gamma rays along axes defined either by the radio jets of the sources or oriented perpendicular to a global magnetic field. We present a pedagogical derivation of the angular structure in the inverse Compton halo and provide an analytic formalism that facilitates the generation of mock images. We discuss exploiting this fact for the purpose of detecting gamma-ray halos in a set of companion papers., Comment: Accepted by ApJ, 18 pages, 17 figures, 13 pages in appendix

TeV-blazars potentially heat the intergalactic medium (IGM) as their gamma rays interact with photons of the extragalactic background light to produce electron-positron pairs, which lose their kinetic energy to the surrounding medium through plasma instabilities. This results in a heating mechanism that is only weakly sensitive to the local density, and therefore approximately spatially uniform, naturally producing an inverted temperature-density relation in underdense regions. In this paper we go beyond the approximation of uniform heating and quantify the heating rate fluctuations due to the clustered distribution of blazars and how this impacts on the thermal history of the IGM. We analytically compute a filtering function that relates the heating rate fluctuations to the underlying dark matter density field. We implement it in the cosmological code GADGET-3 and perform large scale simulations to determine the impact of inhomogeneous heating. We show that, because of blazar clustering, blazar heating is inhomogeneous for z>= 2. At high redshift, the temperature-density relation shows an important scatter and presents a low temperature envelope of unheated regions, in particular at low densities and within voids. However, the median temperature of the IGM is close to that in the uniform case, albeit slightly lower at low redshift. We find that blazar heating is more complex than initially assumed and that the temperature-density relation is not unique. Our analytic model for the heating rate fluctuations couples well with large scale simulations and provides a cost-effective alternative to subgrid models.

Very-high energy gamma-rays from extragalactic sources pair-produce off of the extragalactic background light, yielding an electron-positron pair beam. This pair beam is unstable to various plasma instabilities, especially the "oblique" instability, which can be the dominant cooling mechanism for the beam. However, recently, it has been claimed that nonlinear Landau damping renders it physically irrelevant by reducing the effective damping rate to a low level. Here, we show with numerical calculations that the effective damping rate is $8\times 10^{-4}$ of the growth rate of the linear instability, which is sufficient for the "oblique" instability to be the dominant cooling mechanism of these pair beams. In particular, we show that previous estimates of this rate ignored the exponential cutoff in the scattering amplitude at large wavenumber and assumed that the damping of scattered waves entirely depends on collisions, ignoring collisionless processes. We find that the total wave energy eventually grows to approximate equipartition with the beam by increasingly depositing energy into long wavelength modes. As we have not included the effect of nonlinear wave-wave interactions on these long wavelength modes, this scenario represents the "worst-case" scenario for the oblique instability. As it continues to drain energy from the beam at a faster rate than other processes, we conclude that the "oblique" instability is sufficiently strong to make it the physically dominant cooling mechanism for high-energy pair beams in the intergalactic medium., Comment: 6 pages, 3 figures, accepted to ApJ

We discuss a model for a universe with discrete matter content instead of the continuous perfect fluid taken in FRW models. We show how the redshift in such a universe deviates from the corresponding one in an FRW cosmology. This illustrates the fact that averaging the matter content in a universe and then evolving it in time, is not the same as evolving a universe with discrete matter content. The main reason for such deviation is the fact that the photons in such a universe mainly travel in an empty space rather than the continuous perfect fluid in FRW geometry., Comment: submitted as AIIMS diploma essay at http://archive.aims.ac.za/2010