@article{WilhelmTelezhinskyDwarkadasetal.2020, author = {Wilhelm, Alina and Telezhinsky, Igor and Dwarkadas, Vikram V. and Pohl, Martin}, title = {Stochastic re-acceleration and magnetic-field damping in Tycho's supernova remnant}, series = {Astronomy and astrophysics}, volume = {639}, journal = {Astronomy and astrophysics}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201936079}, pages = {14}, year = {2020}, abstract = {Context. Tycho's supernova remnant (SNR) is associated with the historical supernova (SN) event SN 1572 of Type Ia. The explosion occurred in a relatively clean environment, and was visually observed, providing an age estimate. This SNR therefore represents an ideal astrophysical test-bed for the study of cosmic-ray acceleration and related phenomena. A number of studies suggest that shock acceleration with particle feedback and very efficient magnetic-field amplification combined with Alfvenic drift are needed to explain the rather soft radio spectrum and the narrow rims observed in X-rays. Aims. We show that the broadband spectrum of Tycho's SNR can alternatively be well explained when accounting for stochastic acceleration as a secondary process. The re-acceleration of particles in the turbulent region immediately downstream of the shock should be efficient enough to impact particle spectra over several decades in energy. The so-called Alfvenic drift and particle feedback on the shock structure are not required in this scenario. Additionally, we investigate whether synchrotron losses or magnetic-field damping play a more profound role in the formation of the non-thermal filaments. Methods. We solved the full particle transport equation in test-particle mode using hydrodynamic simulations of the SNR plasma flow. The background magnetic field was either computed from the induction equation or follows analytic profiles, depending on the model considered. Fast-mode waves in the downstream region provide the diffusion of particles in momentum space. Results. We show that the broadband spectrum of Tycho can be well explained if magnetic-field damping and stochastic re-acceleration of particles are taken into account. Although not as efficient as standard diffusive shock acceleration, stochastic acceleration leaves its imprint on the particle spectra, which is especially notable in the emission at radio wavelengths. We find a lower limit for the post-shock magnetic-field strength similar to 330 mu G, implying efficient amplification even for the magnetic-field damping scenario. Magnetic-field damping is necessary for the formation of the filaments in the radio range, while the X-ray filaments are shaped by both the synchrotron losses and magnetic-field damping.}, language = {en} } @article{WarrenEllisonBarkovetal.2017, author = {Warren, Donald C. and Ellison, Donald C. and Barkov, Maxim V. and Nagataki, Shigehiro}, title = {Nonlinear Particle Acceleration and Thermal Particles in GRB Afterglows}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {835}, 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.3847/1538-4357/aa56c3}, pages = {15}, year = {2017}, abstract = {The standard model for GRB afterglow emission treats the accelerated electron population as a simple power law, N(E) proportional to E-p for p greater than or similar to 2. However, in standard Fermi shock acceleration, a substantial fraction of the swept-up particles do not enter the acceleration process at all. Additionally, if acceleration is efficient, then the nonlinear back-reaction of accelerated particles on the shock structure modifies the shape of the nonthermal tail of the particle spectra. Both of these modifications to the standard synchrotron afterglow impact the luminosity, spectra, and temporal variation of the afterglow. To examine the effects of including thermal particles and nonlinear particle acceleration on afterglow emission, we follow a hydrodynamical model for an afterglow jet and simulate acceleration at numerous points during the evolution. When thermal particles are included, we find that the electron population is at no time well fitted by a single power law, though the highest-energy electrons are; if the acceleration is efficient, then the power-law region is even smaller. Our model predicts hard-soft-hard spectral evolution at X-ray energies, as well as an uncoupled X-ray and optical light curve. Additionally, we show that including emission from thermal particles has drastic effects (increases by factors of 100 and 30, respectively) on the observed flux at optical and GeV energies. This enhancement of GeV emission makes afterglow detections by future gamma-ray observatories, such as CTA, very likely.}, language = {en} } @article{tAliuArchambaultArlenetal.2013, author = {tAliu, E. and Archambault, S. and Arlen, T. and Aune, T. and Beilicke, M. and Benbow, W. and Bird, R. and Bouvier, A. and Bradbury, S. M. and Buckley, J. H. and Bugaev, V. and Byrum, K. and Cannon, A. and Cesarini, A. and Ciupik, L. and Collins-Hughes, E. and Connolly, M. P. and Cui, W. and Dickherber, R. and Duke, C. and Dumm, J. and Dwarkadas, Vikram V. and Errando, M. and Falcone, A. and Federici, Simone and Feng, Q. and Finley, J. P. and Finnegan, G. and Fortson, L. and Furniss, A. and Galante, N. and Gall, D. and Gillanders, G. H. and Godambe, S. and Gotthelf, E. V. and Griffin, S. and Grube, J. and Gyuk, G. and Hanna, D. and Holder, J. and Huan, H. and Hughes, G. and Humensky, T. B. and Kaaret, P. and Karlsson, N. and Kertzman, M. and Khassen, Y. and Kieda, D. and Krawczynski, H. and Krennrich, F. and Lang, M. J. and Lee, K. and Madhavan, A. S. and Maier, G. and Majumdar, P. and McArthur, S. and McCann, A. and Millis, J. and Moriarty, P. and Mukherjee, R. and Nelson, T. and de Bhroithe, A. O'Faolain and Ong, R. A. and Orr, M. and Otte, A. N. and Pandel, D. and Park, N. and Perkins, J. S. and Pohl, Martin and Popkow, A. and Prokoph, H. and Quinn, J. and Ragan, K. and Reyes, L. C. and Reynolds, P. T. and Roache, E. and Rose, H. J. and Ruppel, Jens and Saxon, D. B. and Schroedter, M. and Sembroski, G. H. and Sentuerk, G. D. and Skole, C. and Telezhinsky, Igor O. and Tesic, G. and Theiling, M. and Thibadeau, S. and Tsurusaki, K. and Tyler, J. and Varlotta, A. and Vassiliev, V. V. and Vincent, S. and Wakely, S. P. and Ward, J. E. and Weekes, T. C. and Weinstein, A. and Weisgarber, T. and Welsing, R. and Williams, D. A. and Zitzer, B.}, title = {Discovery of TeV Gamma-Ray emission toward supernova remnant SNR G78.2+2.1}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {770}, 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/770/2/93}, pages = {7}, year = {2013}, abstract = {We report the discovery of an unidentified, extended source of very-high-energy gamma-ray emission, VER J2019+407, within the radio shell of the supernova remnant SNR G78.2+2.1, using 21.4 hr of data taken by the VERITAS gamma-ray observatory in 2009. These data confirm the preliminary indications of gamma-ray emission previously seen in a two-year (2007-2009) blind survey of the Cygnus region by VERITAS. VER J2019+407, which is detected at a post-trials significance of 7.5 standard deviations in the 2009 data, is localized to the northwestern rim of the remnant in a region of enhanced radio and X-ray emission. It has an intrinsic extent of 0 degrees.23 +/- 0 degrees.03(stat-0 degrees.02sys)(+0 degrees.04) and its spectrum is well-characterized by a differential power law (dN/dE = N-0 x (E/TeV)-Gamma) with a photon index of Gamma = 2.37 +/- 0.14(stat) +/- 0.20(sys) and a flux normalization of N-0 = 1.5 +/- 0.2(stat) +/- 0.4(sys) x 10(-12) photon TeV-1 cm(-2) s(-1). This yields an integral flux of 5.2 +/- 0.8(stat) +/- 1.4(sys) x 10(-12) photon cm(-2) s(-1) above 320 GeV, corresponding to 3.7\% of the Crab Nebula flux. We consider the relationship of the TeV gamma-ray emission with the GeV gamma-ray emission seen from SNR G78.2+2.1 as well as that seen from a nearby cocoon of freshly accelerated cosmic rays. Multiple scenarios are considered as possible origins for the TeV gamma-ray emission, including hadronic particle acceleration at the SNR shock.}, language = {en} } @article{SushchBrosePohl2018, author = {Sushch, Iurii and Brose, Robert and Pohl, Martin}, title = {Modeling of the spatially resolved nonthermal emission from the Vela Jr. supernova remnant}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {618}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201832879}, pages = {11}, year = {2018}, abstract = {Vela Jr. (RX J0852.0-4622) is one of just a few known supernova remnants (SNRs) with a resolved shell across the whole electromagnetic spectrum from radio to very-high-energy (>100 GeV; VHE) gamma-rays. Its proximity and large size allow for detailed spatially resolved observations of the source, making Vela Jr. one of the primary sources used for the study of particle acceleration and emission mechanisms in SNRs. High-resolution X-ray observations reveal a steepening of the spectrum toward the interior of the remnant. In this study we aim for a self-consistent radiation model of Vela Jr. which at the same time would explain the broadband emission from the source and its intensity distribution. We solve the full particle transport equation combined with the high-resolution one-dimensional (1D) hydrodynamic simulations (using Pluto code) and subsequently calculate the radiation from the remnant. The equations are solved in the test particle regime. We test two models for the magnetic field profile downstream of the shock: damped magnetic field, which accounts for the damping of strong magnetic turbulence downstream, and transported magnetic field. Neither of these scenarios can fully explain the observed radial dependence of the X-ray spectrum under spherical symmetry. We show, however, that the softening of the spectrum and the X-ray intensity profile can be explained under the assumption that the emission is enhanced within a cone.}, language = {en} } @article{RettigPohl2012, author = {Rettig, R. and Pohl, M.}, title = {The properties of non-thermal X-ray filaments in young supernova remnants}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {545}, journal = {Astronomy and astrophysics : an international weekly journal}, number = {5}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201219409}, pages = {9}, year = {2012}, abstract = {Context. Young supernova remnants (SNRs) exhibit narrow filaments of non-thermal X-ray emission whose widths can be limited either by electron energy losses or damping of the magnetic field. Aims. We want to investigate whether or not different models of these filaments can be observationally tested. Methods. Using observational parameters of four historical remnants, we calculated the filament profiles and compared the spectra of the filaments with those of the total non-thermal emission. For that purpose, we solved a one-dimensional stationary transport equation for the isotropic differential number density of the electrons. Results. We find that the difference between the spectra of filament and total non-thermal emission above 1 keV is more pronounced in the damping model than in the energy-loss model. Conclusions. A considerable damping of the magnetic field can result in an observable difference between the spectra of filament and total non-thermal emission, thus potentially permitting an observational discrimination between the energy-loss model and the damping model of the X-ray filaments.}, language = {en} } @article{PohlWilhelmTelezhinsky2015, author = {Pohl, Manuela and Wilhelm, Alina and Telezhinsky, Igor O.}, title = {Reacceleration of electrons in supernova remnants}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {574}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201425027}, pages = {11}, year = {2015}, abstract = {Context. radio spectra of many shell-type supernova remnants show deviations from those expected on theoretical grounds. Aims. In this paper we determine the effect of stochastic reacceleration on the spectra of electrons in the GeV band and at lower energies, and we investigate whether reacceleration can explain the observed variation in radio spectral indices. Methods. We explicitely calculated the momentum diffusion coefficient for 3 types of turbulence expected downstream of the forward shock: fast-mode waves, small-scale non-resonant modes, and large-scale modes arising from turbulent dynamo activity. After noting that low-energy particles are efficiently coupled to the quasi-thermal plasma, a simplified cosmic-ray transport equation can be formulated and is numerically solved. Results. Only fast-mode waves can provide momentum diffusion fast enough to significantly modify the spectra of particles. Using a synchrotron emissivity that accurately reflects a highly turbulent magnetic field, we calculated the radio spectral index and find that soft spectra with index a alpha less than or similar to -0.6 can be maintained over more than 2 decades in radio frequency, even if the electrons experience reacceleration for only one acceleration time. A spectral hardening is possible but considerably more frequency-dependent. The spectral modification imposed by stochastic reacceleration downstream of the forward shock depends only weakly on the initial spectrum provided by, e.g., diffusive shock acceleration at the shock itself.}, 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{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{NiemiecPohlBretetal.2012, author = {Niemiec, Jacek and Pohl, Martin and Bret, Antoine and Wieland, Volkmar}, title = {Nonrelativistic parallel shocks in unmagnetized and weakly magnetized plasmas}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {759}, 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/759/1/73}, pages = {20}, year = {2012}, abstract = {We present results of 2D3V particle-in-cell simulations of nonrelativistic plasma collisions with absent or parallel large-scale magnetic field for parameters applicable to the conditions at young supernova remnants. We study the collision of plasma slabs of different density, leading to two different shocks and a contact discontinuity. Electron dynamics play an important role in the development of the system. While nonrelativistic shocks in both unmagnetized and magnetized plasmas can be mediated by Weibel-type instabilities, the efficiency of shock-formation processes is higher when a large-scale magnetic field is present. The electron distributions downstream of the forward and reverse shocks are generally isotropic, whereas that is not always the case for the ions. We do not see any significant evidence of pre-acceleration, neither in the electron population nor in the ion distribution.}, language = {en} } @article{KobzarNiemiecPohletal.2017, author = {Kobzar, Oleh and Niemiec, Jacek and Pohl, Martin and Bohdan, Artem}, title = {Spatio-temporal evolution of the non-resonant instability in shock precursors of young supernova remnants}, series = {Monthly notices of the Royal Astronomical Society}, volume = {469}, journal = {Monthly notices of the Royal Astronomical Society}, publisher = {Oxford Univ. Press}, address = {Oxford}, organization = {ANTARES Collaboration;H E S S Collaboration}, issn = {0035-8711}, doi = {10.1093/mnras/stx1201}, pages = {4985 -- 4998}, year = {2017}, abstract = {A non-resonant cosmic ray (CR) current-driven instability may operate in the shock precursors of young supernova remnants and be responsible for magnetic-field amplification, plasma heating and turbulence. Earlier simulations demonstrated magnetic-field amplification, and in kinetic studies a reduction of the relative drift between CRs and thermal plasma was observed as backreaction. However, all published simulations used periodic boundary conditions, which do not account for mass conservation in decelerating flows and only allow the temporal development to be studied. Here we report results of fully kinetic particle-in-cell simulations with open boundaries that permit inflow of plasma on one side of the simulation box and outflow at the other end, hence allowing an investigation of both the temporal and the spatial development of the instability. Magnetic-field amplification proceeds as in studies with periodic boundaries and, observed here for the first time, the reduction of relative drifts causes the formation of a shock-like compression structure at which a fraction of the plasma ions are reflected. Turbulent electric field generated by the non-resonant instability inelastically scatters CRs, modifying and anisotropizing their energy distribution. Spatial CR scattering is compatible with Bohm diffusion. Electromagnetic turbulence leads to significant non-adiabatic heating of the background plasma maintaining bulk equipartition between ions and electrons. The highest temperatures are reached at sites of large-amplitude electrostatic fields. Ion spectra show supra-thermal tails resulting from stochastic scattering in the turbulent electric field. Together, these modifications in the plasma flow will affect the properties of the shock and particle acceleration there.}, language = {en} } @article{FraschettiPohl2017, author = {Fraschetti, F. and Pohl, Martin}, title = {Particle acceleration model for the broad-band baseline spectrum of the Crab nebula}, series = {Monthly notices of the Royal Astronomical Society}, volume = {471}, journal = {Monthly notices of the Royal Astronomical Society}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stx1833}, pages = {4866 -- 4874}, year = {2017}, abstract = {We develop a simple one-zone model of the steady-state Crab nebula spectrum encompassing both the radio/soft X-ray and the GeV/multi-TeV observations. By solving the transport equation for GeV-TeV electrons injected at the wind termination shock as a log-parabola momentum distribution and evolved via energy losses, we determine analytically the resulting differential energy spectrum of photons. We find an impressive agreement with the observed spectrum of synchrotron emission, and the synchrotron self-Compton component reproduces the previously unexplained broad 200-GeV peak that matches the Fermi/Large Area Telescope (LAT) data beyond 1 GeV with the Major Atmospheric Gamma Imaging Cherenkov (MAGIC) data. We determine the parameters of the single log-parabola electron injection distribution, in contrast with multiple broken power-law electron spectra proposed in the literature. The resulting photon differential spectrum provides a natural interpretation of the deviation from power law customarily fitted with empirical multiple broken power laws. Our model can be applied to the radio-to-multi-TeV spectrum of a variety of astrophysical outflows, including pulsar wind nebulae and supernova remnants, as well as to interplanetary shocks.}, language = {en} } @article{DeBeckerdelValleRomeroetal.2017, author = {De Becker, M. and del Valle, Maria Victoria and Romero, G. E. and Peri, C. S. and Benaglia, P.}, title = {X- ray study of bow shocks in runaway stars}, series = {Monthly notices of the Royal Astronomical Society}, volume = {471}, journal = {Monthly notices of the Royal Astronomical Society}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stx1826}, pages = {4452 -- 4464}, year = {2017}, abstract = {Massive runaway stars produce bow shocks through the interaction of their winds with the interstellar medium, with the prospect for particle acceleration by the shocks. These objects are consequently candidates for non-thermal emission. Our aim is to investigate the X-ray emission from these sources. We observed with XMM-Newton a sample of five bow shock runaways, which constitutes a significant improvement of the sample of bow shock runaways studied in X-rays so far. A careful analysis of the data did not reveal any X-ray emission related to the bow shocks. However, X-ray emission from the stars is detected, in agreement with the expected thermal emission from stellar winds. On the basis of background measurements we derive conservative upper limits between 0.3 and 10 keV on the bow shocks emission. Using a simple radiation model, these limits together with radio upper limits allow us to constrain some of the main physical quantities involved in the non-thermal emission processes, such as the magnetic field strength and the amount of incident infrared photons. The reasons likely responsible for the non-detection of non-thermal radiation are discussed. Finally, using energy budget arguments, we investigate the detectability of inverse Compton X-rays in a more extended sample of catalogued runaway star bow shocks. From our analysis we conclude that a clear identification of non-thermal X-rays from massive runaway bow shocks requires one order of magnitude (or higher) sensitivity improvement with respect to present observatories.}, language = {en} } @article{ChenPohlBoettcher2015, author = {Chen, Xuhui and Pohl, Martin and B{\"o}ttcher, Markus}, title = {Particle diffusion and localized acceleration in inhomogeneous AGN jets - I. Steady-state spectra}, series = {Monthly notices of the Royal Astronomical Society}, volume = {447}, journal = {Monthly notices of the Royal Astronomical Society}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stu2438}, pages = {530 -- 544}, year = {2015}, abstract = {We study the acceleration, transport, and emission of particles in relativistic jets. Localized stochastic particle acceleration, spatial diffusion, and synchrotron as well as synchrotron self-Compton (SSC) emission are considered in a leptonic model. To account for inhomogeneity, we use a 2D axisymmetric cylindrical geometry for both relativistic electrons and magnetic field. In this first phase of our work, we focus on steady-state spectra that develop from a time-dependent model. We demonstrate that small isolated acceleration region in a much larger emission volume are sufficient to accelerate particles to high energy. Diffusive escape from these small regions provides a natural explanation for the spectral form of the jet emission. The location of the acceleration regions within the jet is found to affect the cooling break of the spectrum in this diffusive model. Diffusion-caused energy-dependent inhomogeneity in the jets predicts that the SSC spectrum is harder than the synchrotron spectrum. There can also be a spectral hardening towards the high-energy section of the synchrotron spectrum, if particle escape is relatively slow. These two spectral hardening effects indicate that the jet inhomogeneity might be a natural explanation for the unexpected hard. gamma-ray spectra observed in some blazars.}, language = {en} } @article{ChenPohlBottcheretal.2016, author = {Chen, Xuhui and Pohl, Martin and Bottcher, Markus and Gao, Shan}, title = {Particle diffusion and localized acceleration in inhomogeneous AGN jets - II. Stochastic variation}, series = {Monthly notices of the Royal Astronomical Society}, volume = {458}, journal = {Monthly notices of the Royal Astronomical Society}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stw528}, pages = {3260 -- 3271}, year = {2016}, abstract = {We study the stochastic variation of blazar emission under a 2D spatially resolved leptonic jet model we previously developed. Random events of particle acceleration and injection in small zones within the emission region are assumed to be responsible for flux variations. In addition to producing spectral energy distributions that describe the observed flux of Mrk 421, we further analyse the timing properties of the simulated light curves, such as the power spectral density (PSD) at different bands, flux-flux correlations, aswell as the cross-correlation function between X-rays and TeV gamma-rays. We find spectral breaks in the PSD at a time-scale comparable to the dominant characteristic time-scale in the system, which is usually the predefined decay time-scale of an acceleration event. Cooling imposes a delay, and so PSDs taken at lower energy bands in each emission component (synchrotron or inverse Compton) generally break at longer time-scales. The flux-flux correlation between X-rays and TeV gamma-rays can be either quadratic or linear, depending on whether or not there are large variation of the injection into the particle acceleration process. When the relationship is quadratic, the TeV flares lag the X-ray flares, and the optical and GeV flares are large enough to be comparable to the ones in X-ray. When the relationship is linear, the lags are insignificant, and the optical and GeV flares are small.}, language = {en} } @article{BroseTelezhinskyPohl2016, author = {Brose, Robert and Telezhinsky, Igor O. and Pohl, Martin}, title = {Transport of magnetic turbulence in supernova remnants}, series = {Physical review letters}, volume = {593}, journal = {Physical review letters}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201527345}, pages = {8}, year = {2016}, abstract = {Context. Supernova remnants are known as sources of Galactic cosmic rays for their nonthermal emission of radio waves, X-rays, and gamma rays. However, the observed soft broken power-law spectra are hard to reproduce within standard acceleration theory based on the assumption of Bohm diffusion and steady-state calculations. Aims. We point out that a time-dependent treatment of the acceleration process together with a self-consistent treatment of the scattering turbulence amplification is necessary. Methods. We numerically solve the coupled system of transport equations for cosmic rays and isotropic Alfvenic turbulence. The equations are coupled through the growth rate of turbulence determined by the cosmic-ray gradient and the spatial diffusion coefficient of cosmic rays determined by the energy density of the turbulence. The system is solved on a comoving expanding grid extending upstream for dozens of shock radii, allowing for the self-consistent study of cosmic-ray diffusion in the vicinity of their acceleration site. The transport equation for cosmic rays is solved in a test-particle approach. Results. We demonstrate that the system is typically not in a steady state. In fact, even after several thousand years of evolution, no equilibrium situation is reached. The resulting time-dependent particle spectra strongly differ from those derived assuming a steady state and Bohm diffusion. Our results indicate that proper accounting for the evolution of the scattering turbulence and hence the particle diffusion coefficient is crucial for the formation of the observed soft spectra. In any case, the need to continuously develop magnetic turbulence upstream of the shock introduces nonlinearity in addition to that imposed by cosmic-ray feedback.}, language = {en} } @article{BroseSushchPohletal.2019, author = {Brose, Robert and Sushch, Iuri and Pohl, Martin and Luken, K. J. and Filipovic, M. D. and Lin, R.}, title = {Nonthermal emission from the reverse shock of the youngest galactic supernova remnant G1.9+0.3}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {627}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {1432-0746}, doi = {10.1051/0004-6361/201834430}, pages = {9}, year = {2019}, abstract = {Context. The youngest Galactic supernova remnant G1.9+0.3 is an interesting target for next-generation gamma-ray observatories. So far, the remnant is only detected in the radio and the X-ray bands, but its young age of approximate to 100 yr and inferred shock speed of approximate to 14 000 km s(-1) could make it an efficient particle accelerator. Aims. We aim to model the observed radio and X-ray spectra together with the morphology of the remnant. At the same time, we aim to estimate the gamma-ray flux from the source and evaluate the prospects of its detection with future gamma-ray experiments. Methods. We performed spherical symmetric 1D simulations with the RATPaC code, in which we simultaneously solved the transport equation for cosmic rays, the transport equation for magnetic turbulence, and the hydro-dynamical equations for the gas flow. Separately computed distributions of the particles accelerated at the forward and the reverse shock were then used to calculate the spectra of synchrotron, inverse Compton, and pion-decay radiation from the source. Results. The emission from G1.9+0.3 can be self-consistently explained within the test-particle limit. We find that the X-ray flux is dominated by emission from the forward shock while most of the radio emission originates near the reverse shock, which makes G1.9+0.3 the first remnant with nonthermal radiation detected from the reverse shock. The flux of very-high-energy gamma-ray emission from G1.9+0.3 is expected to be close to the sensitivity threshold of the Cherenkov Telescope Array. The limited time available to grow large-scale turbulence limits the maximum energy of particles to values below 100 TeV, hence G1.9+0.3 is not a PeVatron.}, language = {en} } @article{BohdanNiemiecPohletal.2019, author = {Bohdan, Artem and Niemiec, Jacek and Pohl, Martin and Matsumoto, Yosuke and Amano, Takanobu and Hoshino, Masahiro}, title = {Kinetic Simulations of Nonrelativistic Perpendicular Shocks of Young Supernova Remnants}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {878}, 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/ab1b6d}, pages = {11}, year = {2019}, abstract = {Electron injection at high Mach number nonrelativistic perpendicular shocks is studied here for parameters that are applicable to young SNR shocks. Using high-resolution large-scale two-dimensional fully kinetic particle-in-cell simulations and tracing individual particles, we in detail analyze the shock-surfing acceleration (SSA) of electrons at the leading edge of the shock foot. The central question is to what degree the process can be captured in 2D3V simulations. We find that the energy gain in SSA always arises from the electrostatic field of a Buneman wave. Electron energization is more efficient in the out-of-plane orientation of the large-scale magnetic field because both the phase speed and the amplitude of the waves are higher than for the in-plane scenario. Also, a larger number of electrons is trapped by the waves compared to the in-plane configuration. We conclude that significant modifications of the simulation parameters are needed to reach the same level of SSA efficiency as in simulations with out-of-plane magnetic field or 3D simulations.}, language = {en} } @article{BohdanNiemiecKobzaretal.2017, author = {Bohdan, Artem and Niemiec, Jacek and Kobzar, Oleh and Pohl, Martin}, title = {Electron Pre-acceleration at Nonrelativistic High-Mach-number Perpendicular Shocks}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {847}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, issn = {0004-637X}, doi = {10.3847/1538-4357/aa872a}, pages = {17}, year = {2017}, abstract = {We perform particle-in-cell simulations of perpendicular nonrelativistic collisionless shocks to study electron heating and pre-acceleration for parameters that permit the extrapolation to the conditions at young supernova remnants. Our high-resolution large-scale numerical experiments sample a representative portion of the shock surface and demonstrate that the efficiency of electron injection is strongly modulated with the phase of the shock reformation. For plasmas with low and moderate temperature (plasma beta beta p =5.10(-4) and 0.5 beta p =), we explore the nonlinear shock structure and electron pre-acceleration for various orientations of the large-scale magnetic field with respect to the simulation plane, while keeping it at 90 degrees to the shock normal. Ion reflection off of the shock leads to the formation of magnetic filaments in the shock ramp, resulting from Weibel-type instabilities, and electrostatic Buneman modes in the shock foot. In all of the cases under study, the latter provides first-stage electron energization through the shock-surfing acceleration mechanism. The subsequent energization strongly depends on the field orientation and proceeds through adiabatic or second-order Fermi acceleration processes for configurations with the out-of-plane and in-plane field components, respectively. For strictly out-of-plane field, the fraction of suprathermal electrons is much higher than for other configurations, because only in this case are the Buneman modes fully captured by the 2D simulation grid. Shocks in plasma with moderate bp provide more efficient pre-acceleration. The relevance of our results to the physics of fully 3D systems is discussed.}, language = {en} } @article{ArchambaultBeilickeBenbowetal.2013, author = {Archambault, S. and Beilicke, M. and Benbow, W. and Berger, K. and Bird, R. and Bouvier, A. and Buckley, J. H. and Bugaev, V. and Byrum, K. and Cerruti, M. and Chen, Xuhui and Ciupik, L. and Connolly, M. P. and Cui, W. and Duke, C. and Dumm, J. and Errando, M. and Falcone, A. and Federici, Simone and Feng, Q. and Finley, J. 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 Krawczynski, H. and Lang, M. J. and Madhavan, A. S. and Maier, G. and Majumdar, P. and McArthur, S. and McCann, A. and Moriarty, P. and Mukherjee, R. and Nieto, D. and de Bhroithe, A. O'Faolain and Ong, R. A. and Otte, A. N. and Pandel, D. 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 Reyes, L. C. and Reynolds, P. T. and Richards, G. T. and Roache, E. and Sembroski, G. H. and Sheidaei, F. and Smith, A. W. and Staszak, D. and Telezhinsky, Igor O. and Theiling, M. and Tucci, J. V. and Tyler, J. and Varlotta, A. and Vincent, S. and Wakely, S. P. and Weekes, T. C. and Weinstein, A. and Williams, D. A. and Zitzer, B. and McCollough, M. L.}, title = {Veritas observatons of the microquasar cygnus X-3}, 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 = {2}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, organization = {VERITAS Collaboration, Smithsonian Astrophys Observ}, issn = {0004-637X}, doi = {10.1088/0004-637X/779/2/150}, pages = {10}, year = {2013}, abstract = {We report results from TeV gamma-ray observations of the microquasar Cygnus X-3. The observations were made with the Very Energetic Radiation Imaging Telescope Array System (VERITAS) over a time period from 2007 June 11 to 2011 November 28. VERITAS is most sensitive to gamma rays at energies between 85 GeV and 30 TeV. The effective exposure time amounts to a total of about 44 hr, with the observations covering six distinct radio/X-ray states of the object. No significant TeV gamma-ray emission was detected in any of the states, nor with all observations combined. The lack of a positive signal, especially in the states where GeV gamma rays were detected, places constraints on TeV gamma-ray production in Cygnus X-3. We discuss the implications of the results.}, 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} }