@article{MizunoPohlNiemiecetal.2011,
  author    = {Mizuno, Yosuke and Pohl, Martin and Niemiec, Jacek and Zhang, Bing and Nishikawa, Ken-Ichi and Hardee, Philip E.},
  title     = {Magnetic-field amplification by turbulence in a relativistic shockpropagating through an inhomogeneous medium},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {726},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {2},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.1088/0004-637X/726/2/62},
  pages     = {11},
  year      = {2011},
  abstract  = {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.},
  language  = {en}
}
@article{AliuArchambaultBeheraetal.2013,
  author    = {Aliu, E. and Archambault, S. and Behera, B. and Berger, K. and Beilicke, M. and Benbow, W. and Bird, R. and Bouvier, A. and Bugaev, V. and Cerruti, M. and Chen, Xuhui and Ciupik, L. and Connolly, M. P. and Cui, W. and Dumm, J. and Falcone, A. and Federici, Simone and Feng, Q. and Finley, J. P. and Fortin, P. and Fortson, L. and Furniss, A. and Galante, N. and Gillanders, G. H. and Griffin, S. and Griffiths, S. T. and Grube, J. and Gyuk, G. and Hanna, D. and Holder, J. and Hughes, G. and Humensky, T. B. and Kaaret, P. and Kertzman, M. and Khassen, Y. and Kieda, D. and Krennrich, F. and Lang, M. J. and Maier, G. and Majumdar, P. and McArthur, S. and McCann, A. and Moriarty, P. and Mukherjee, R. and de Bhroithe, A. O'Faolain and Ong, R. A. and Otte, A. N. and Park, N. and Perkins, J. S. and Pohl, Martin and Popkow, A. and Prokoph, H. and Quinn, J. and Ragan, K. and Rajotte, J. and Ratliff, G. and Reynolds, P. T. and Richards, G. T. and Roache, E. and Sembroski, G. H. and Sheidaei, F. and Skole, C. and Smith, A. W. and Staszak, D. and Telezhinsky, Igor O. and Tyler, J. and Varlotta, A. and Vincent, S. and Wakely, S. P. and Weekes, T. C. and Weinstein, A. and Welsing, R. and Zajczyk, A. and Zitzer, B.},
  title     = {Multiwavelemght oservatons of the TeV binary LS I+61 degrees 303 with veritas, fermi-lat, and swift/XRT during a TeV outburst},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {779},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {1},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.1088/0004-637X/779/1/88},
  pages     = {7},
  year      = {2013},
  abstract  = {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.},
  language  = {en}
}
@article{MizunoPohlNiemiecetal.2014,
  author    = {Mizuno, Yosuke and Pohl, Martin and Niemiec, Jacek and Zhang, Bing and Nishikawa, Ken-Ichi and Hardee, Philip E.},
  title     = {Magnetic field amplification and saturation in turbulence behind a relativistic shock},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {439},
  journal   = {Monthly notices of the Royal Astronomical Society},
  number    = {4},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stu196},
  pages     = {3490 -- 3503},
  year      = {2014},
  abstract  = {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.},
  language  = {en}
}
@article{NishikawaHardeeDutanetal.2014,
  author    = {Nishikawa, Ken-Ichi and Hardee, P. E. and Dutan, I. and Niemiec, J. and Medvedev, M. and Mizuno, Y. and Meli, A. and Sol, H. and Zhang, B. and Pohl, Martin and Hartmann, D. H.},
  title     = {Magnetic agnetic field generation in core-sheath jets via the kinetic Kelvin-Helmholtz instability},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {793},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {1},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.1088/0004-637X/793/1/60},
  pages     = {16},
  year      = {2014},
  abstract  = {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.},
  language  = {en}
}
@article{AbeysekaraArchambaultArcheretal.2016,
  author    = {Abeysekara, A. U. and Archambault, S. and Archer, A. and Benbow, W. and Bird, R. and Biteau, Jonathan and Buchovecky, M. and Buckley, J. H. and Bugaev, V. and Byrum, K. and Cardenzana, J. V. and Cerruti, M. and Chen, Xuhui and Christiansen, J. L. and Ciupik, L. and Connolly, M. P. and Cui, W. and Dickinson, H. J. and Dumm, J. and Eisch, J. D. and Errando, M. and Falcone, A. and Feng, Q. and Finley, J. P. and Fleischhack, H. and Flinders, A. and Fortin, P. and Fortson, L. and Furniss, A. and Gillanders, G. H. and Griffin, S. and Grube, J. and Gyuk, G. and Huetten, M. and Hanna, D. and Holder, J. and Humensky, T. B. and Johnson, C. A. and Kaaret, P. and Kar, P. and Kelley-Hoskins, N. and Kertzman, M. and Kieda, D. and Krause, M. and Krennrich, F. and Lang, M. J. and Maier, G. and McArthur, S. and McCann, A. and Meagher, K. and Moriarty, P. and Mukherjee, R. and Nieto, D. and Ong, R. A. and Otte, A. N. and Park, N. and Pelassa, V. and Petrashyk, A. and Petry, D. and Pohl, Martin and Popkow, A. and Pueschel, Elisa and Quinn, J. and Ragan, K. and Ratliff, G. and Reyes, L. C. and Reynolds, P. T. and Reynolds, K. and Richards, G. T. and Roache, E. and Rulten, C. and Santander, M. and Sembroski, G. H. and Shahinyan, K. and Smith, A. W. and Staszak, D. and Telezhinsky, Igor O. and Tucci, J. V. and Tyler, J. and Vincent, S. and Wakely, S. P. and Weiner, O. M. and Weinstein, A. and Wilhelm, Alina and Williams, D. A. and Zitzer, B.},
  title     = {VERITAS and multiwavelength observations of the BL Lacertae object 1ES 1741+196},
  series = {Monthly notices of the Royal Astronomical Society},
  volume    = {459},
  journal   = {Monthly notices of the Royal Astronomical Society},
  publisher = {Oxford Univ. Press},
  address   = {Oxford},
  issn      = {0035-8711},
  doi       = {10.1093/mnras/stw664},
  pages     = {2550 -- 2557},
  year      = {2016},
  abstract  = {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.},
  language  = {en}
}
@article{NishikawaFrederiksenNordlundetal.2016,
  author    = {Nishikawa, Ken-Ichi and Frederiksen, J. T. and Nordlund, A. and Mizuno, Y. and Hardee, P. E. and Niemiec, J. and Gomez, J. L. and Dutan, I. and Meli, A. and Sol, H. and Pohl, Martin and Hartmann, D. H.},
  title     = {EVOLUTION OF GLOBAL RELATIVISTIC JETS: COLLIMATIONS AND EXPANSION WITH kKHI AND THE WEIBEL INSTABILITY},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {820},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.3847/0004-637X/820/2/94},
  pages     = {14},
  year      = {2016},
  abstract  = {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.},
  language  = {en}
}
@article{VafinRafighiPohletal.2018,
  author    = {Vafin, Sergei and Rafighi, Iman and Pohl, Martin and Niemiec, Jacek},
  title     = {The Electrostatic Instability for Realistic Pair Distributions in Blazar/EBL Cascades},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {857},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {1},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.3847/1538-4357/aab552},
  pages     = {12},
  year      = {2018},
  abstract  = {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.},
  language  = {en}
}
@article{VafinDekaPohletal.2019,
  author    = {Vafin, Sergei and Deka, Pranab Jyoti and Pohl, Martin and Bohdan, Artem},
  title     = {Revisit of Nonlinear Landau Damping for Electrostatic Instability Driven by Blazar-induced Pair Beams},
  series = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  volume    = {873},
  journal   = {The astrophysical journal : an international review of spectroscopy and astronomical physics},
  number    = {1},
  publisher = {IOP Publ. Ltd.},
  address   = {Bristol},
  issn      = {0004-637X},
  doi       = {10.3847/1538-4357/ab017b},
  pages     = {12},
  year      = {2019},
  abstract  = {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.},
  language  = {en}
}