Refine
Has Fulltext
- no (8)
Document Type
- Article (8)
Language
- English (8) (remove)
Is part of the Bibliography
- yes (8) (remove)
Keywords
- relativistic processes (8) (remove)
Institute
This work revisits the electrostatic instability for blazar-induced pair beams propagating through the intergalactic medium (IGM) using linear analysis and PIC simulations. We study the impact of the realistic distribution function of pairs resulting from the interaction of high-energy gamma-rays with the extragalactic background light. We present analytical and numerical calculations of the linear growth rate of the instability for the arbitrary orientation of wave vectors. Our results explicitly demonstrate that the finite angular spread of the beam dramatically affects the growth rate of the waves, leading to the fastest growth for wave vectors quasi-parallel to the beam direction and a growth rate at oblique directions that is only a factor of 2-4 smaller compared to the maximum. To study the nonlinear beam relaxation, we performed PIC simulations that take into account a realistic wide-energy distribution of beam particles. The parameters of the simulated beam-plasma system provide an adequate physical picture that can be extrapolated to realistic blazar-induced pairs. In our simulations, the beam looses only 1% of its energy, and we analytically estimate that the beam would lose its total energy over about 100 simulation times. An analytical scaling is then used to extrapolate the parameters of realistic blazar-induced pair beams. We find that they can dissipate their energy slightly faster by the electrostatic instability than through inverse-Compton scattering. The uncertainties arising from, e.g., details of the primary gamma-ray spectrum are too large to make firm statements for individual blazars, and an analysis based on their specific properties is required.
We revisit the effect of nonlinear Landau (NL) damping on the electrostatic instability of blazar-induced pair beams, using a realistic pair-beam distribution. We employ a simplified 2D model in k-space to study the evolution of the electric-field spectrum and to calculate the relaxation time of the beam. We demonstrate that the 2D model is an adequate representation of the 3D physics. We find that nonlinear Landau damping, once it operates efficiently, transports essentially the entire wave energy to small wave numbers where wave driving is weak or absent. The relaxation time also strongly depends on the intergalactic medium temperature, T-IGM, and for T-IGM << 10 eV, and in the absence of any other damping mechanism, the relaxation time of the pair beam is longer than the inverse Compton (IC) scattering time. The weak late-time beam energy losses arise from the accumulation of wave energy at small k, that nonlinearly drains the wave energy at the resonant k of the pair-beam instability. Any other dissipation process operating at small k would reduce that wave-energy drain and hence lead to stronger pair-beam energy losses. As an example, collisions reduce the relaxation time by an order of magnitude, although their rate is very small. Other nonlinear processes, such as the modulation instability, could provide additional damping of the nonresonant waves and dramatically reduce the relaxation time of the pair beam. An accurate description of the spectral evolution of the electrostatic waves is crucial for calculating the relaxation time of the pair beam.
Magnetic agnetic field generation in core-sheath jets via the kinetic Kelvin-Helmholtz instability
(2014)
We have investigated magnetic field generation in velocity shears via the kinetic Kelvin-Helmholtz instability (kKHI) using a relativistic plasma jet core and stationary plasma sheath. Our three-dimensional particle-in-cell simulations consider plasma jet cores with Lorentz factors of 1.5, 5, and 15 for both electron-proton and electron-positron plasmas. For electron-proton plasmas, we find generation of strong large-scale DC currents and magnetic fields that extend over the entire shear surface and reach thicknesses of a few tens of electron skin depths. For electron-positron plasmas, we find generation of alternating currents and magnetic fields. Jet and sheath plasmas are accelerated across the shear surface in the strong magnetic fields generated by the kKHI. The mixing of jet and sheath plasmas generates a transverse structure similar to that produced by the Weibel instability.
In the study of relativistic jets one of the key open questions is their interaction with the environment. Here. we study the initial evolution of both electron-proton (e(-) - p(+)) and electron-positron (e(+/-)) relativistic jets, focusing on their lateral interaction with ambient plasma. We follow the evolution of toroidal magnetic fields generated by both the kinetic Kelvin-Helmholtz and Mushroom instabilities. For an e(-) - p(+) jet, the induced magnetic field collimates the jet and electrons are perpendicularly accelerated. As the instabilities saturate and subsequently weaken, the magnetic polarity switches from clockwise to counterclockwise in the middle of the jet. For an e(+/-) jet, we find strong mixing of electrons and positrons with the ambient plasma, resulting in the creation of a bow shock. The merging of current filaments generates density inhomogeneities that. initiate a forward shock. Strong jet-ambient plasma mixing prevents a full development of the jet (on the scale studied), revealing evidence for both jet collimation and particle acceleration in the forming bow shock. Differences in the magnetic field structure generated by e(-) - p(+) and e(+/-) jets may contribute to the polarization properties of the observed emission in AGN jets and gamma-ray bursts.
We perform two-dimensional relativistic magnetohydrodynamic simulations of a mildly relativistic shock propagating through an inhomogeneous medium. We show that the postshock region becomes turbulent owing to preshock density inhomogeneity, and the magnetic field is strongly amplified due to the stretching and folding of field lines in the turbulent velocity field. The amplified magnetic field evolves into a filamentary structure in two-dimensional simulations. The magnetic energy spectrum is flatter than the Kolmogorov spectrum and indicates that a so-called small-scale dynamo is occurring in the postshock region. We also find that the amount of magnetic-field amplification depends on the direction of the mean preshock magnetic field, and the timescale of magnetic-field growth depends on the shock strength.
We have investigated via 2D relativistic magnetohydrodynamic simulations the long-term evolution of turbulence created by a relativistic shock propagating through an inhomogeneous medium. In the post-shock region, magnetic field is strongly amplified by turbulent motions triggered by pre-shock density inhomogeneities. Using a long-simulation box we have followed the magnetic field amplification until it is fully developed and saturated. The turbulent velocity is subrelativistic even for a strong shock. Magnetic field amplification is controlled by the turbulent motion and saturation occurs when the magnetic energy is comparable to the turbulent kinetic energy. Magnetic field amplification and saturation depend on the initial strength and direction of the magnetic field in the pre-shock medium, and on the shock strength. If the initial magnetic field is perpendicular to the shock normal, the magnetic field is first compressed at the shock and then can be amplified by turbulent motion in the post-shock region. Saturation occurs when the magnetic energy becomes comparable to the turbulent kinetic energy in the post-shock region. If the initial magnetic field in the pre-shock medium is strong, the post-shock region becomes turbulent but significant field amplification does not occur. If the magnetic energy after shock compression is larger than the turbulent kinetic energy in the post-shock region, significant field amplification does not occur. We discuss possible applications of our results to gamma-ray bursts and active galactic nuclei.
We present the results of a multiwavelength observational campaign on the TeV binary system LS I +61 degrees 303 with the VERITAS telescope array (>200 GeV), Fermi-LAT (0.3-300 GeV), and Swift/XRT (2-10 keV). The data were taken from 2011 December through 2012 January and show a strong detection in all three wavebands. During this period VERITAS obtained 24.9 hr of quality selected livetime data in which LS I +61 degrees 303 was detected at a statistical significance of 11.9 sigma. These TeV observations show evidence for nightly variability in the TeV regime at a post-trial significance of 3.6 sigma. The combination of the simultaneously obtained TeV and X-ray fluxes do not demonstrate any evidence for a correlation between emission in the two bands. For the first time since the launch of the Fermi satellite in 2008, this TeV detection allows the construction of a detailed MeV-TeV spectral energy distribution from LS I +61 degrees 303. This spectrum shows a distinct cutoff in emission near 4 GeV, with emission seen by the VERITAS observations following a simple power-law above 200 GeV. This feature in the spectrum of LS I +61 degrees 303, obtained from overlapping observations with Fermi-LAT and VERITAS, may indicate that there are two distinct populations of accelerated particles producing the GeV and TeV emission.
We present results from multiwavelength observations of the BL Lacertae object 1ES 1741 + 196, including results in the very high energy gamma-ray regime using the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The VERITAS time-averaged spectrum, measured above 180 GeV, is well modelled by a power law with a spectral index of 2.7 +/- 0.7(stat) +/- 0.2(syst). The integral flux above 180 GeV is (3.9 +/- 0.8(stat) +/- 1.0(syst)) x 10(-8) m(-2) s(-1), corresponding to 1.6 per cent of the Crab nebula flux on average. The multiwavelength spectral energy distribution of the source suggests that 1ES 1741+196 is an extreme-high-frequency-peaked BL Lacertae object. The observations analysed in this paper extend over a period of six years, during which time no strong flares were observed in any band. This analysis is therefore one of the few characterizations of a blazar in a non-flaring state.