@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{TelezhinskyDwarkadasPohl2012, author = {Telezhinsky, Igor O. and Dwarkadas, Vikram V. and Pohl, Martin}, title = {Time-dependent escape of cosmic rays from supernova remnants, and their interaction with dense media}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {541}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201118639}, pages = {11}, year = {2012}, abstract = {Context. Supernova remnants (SNRs) are thought to be the main source of Galactic cosmic rays (CRs) up to the "knee" in CR spectrum. During the evolution of a SNR, the bulk of the CRs are confined inside the SNR shell. The highest-energy particles leave the system continuously, while the remaining adiabatically cooled particles are released when the SNR has expanded sufficiently and decelerated so that the magnetic field at the shock is no longer able to confine them. Particles escaping from the parent system may interact with nearby molecular clouds, producing.-rays in the process via pion decay. The soft gamma-ray spectra observed for a number of SNRs interacting with molecular clouds, however, challenge current theories of non-linear particle acceleration that predict harder spectra. Aims. We study how the spectrum of escaped particles depends on the time-dependent acceleration history in both Type Ia and core-collapse SNRs, as well as on different assumptions about the diffusion coefficient in the vicinity of the SNR. Methods. We solve the CR transport equation in a test-particle approach combined with numerical simulations of SNR evolution. Results. We extend our method for calculating the CR acceleration in SNRs to trace the escaped particles in a large volume around SNRs. We calculate the evolution of the spectra of CRs that have escaped from a SNR into a molecular cloud or dense shell for two diffusion models. We find a strong confinement of CRs in a close region around the SNR, and a strong dilution effect for CRs that were able to propagate out as far as a few SNR radii.}, 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{PetrukKuzyoOrlandoetal.2018, author = {Petruk, Oleh and Kuzyo, T. and Orlando, S. and Pohl, Martin and Miceli, M. and Bocchino, F. and Beshley, V. and Brose, Robert}, title = {Post-adiabatic supernova remnants in an interstellar magnetic field}, series = {Monthly notices of the Royal Astronomical Society}, volume = {479}, journal = {Monthly notices of the Royal Astronomical Society}, number = {3}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/sty1750}, pages = {4253 -- 4270}, year = {2018}, abstract = {We present very-high-resolution 1D MHD simulations of the late-stage supernova remnants (SNRs). In the post-adiabatic stage, the magnetic field has an important and significant dynamical effect on the shock dynamics, the flow structure, and hence the acceleration and emission of cosmic rays. We find that the tangential component of the magnetic field provides pressure support that to a fair degree prevents the collapse of the radiative shell and thus limits the total compression ratio of the partially or fully radiative forward shock. A consequence is that the spectra of cosmic rays would not be as hard as in hydrodynamic simulations. We also investigated the effect on the flow profiles of the magnetic-field inclination and a large-scale gradient in the gas density and/or the magnetic field. A positive density gradient shortens the evolutionary stages, whereas a shock obliquity lowers the shock compression. The compression of the tangential component of the magnetic field leads to its dominance in the downstream region of post-adiabatic shocks for a wide range of orientation of the upstream field, which may explain why one preferentially observes tangential radio polarization in old SNRs. As most cosmic rays are produced at late stages of SNR evolution, the post-adiabatic phase and the influence of the magnetic field during it are most important for modeling the cosmic-ray acceleration at old SNRs and the gamma-ray emission from late-stage SNRs interacting with clouds.}, 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{MeyerVelazquezPetruketal.2022, author = {Meyer, Dominique M.-A. and Velazquez, Pablo F. and Petruk, Oleh and Chiotellis, Alexandros and Pohl, Martin and Camps-Farina, Artemi and Petrov, Miroslav and Reynoso, Estela M. and Toledo-Roy, Juan C. and Schneiter, E. Matias and Castellanos-Ramirez, Antonio and Esquivel, Alejandro}, title = {Rectangular core-collapse supernova remnants}, series = {Monthly notices of the Royal Astronomical Society}, volume = {515}, journal = {Monthly notices of the Royal Astronomical Society}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/stac1832}, pages = {594 -- 605}, year = {2022}, abstract = {Core-collapse supernova remnants are the gaseous nebulae of galactic interstellar media (ISM) formed after the explosive death of massive stars. Their morphology and emission properties depend both on the surrounding circumstellar structure shaped by the stellar wind-ISM interaction of the progenitor star and on the local conditions of the ambient medium. In the warm phase of the Galactic plane (n approximate to 1 cm(-3), T approximate to 8000 K), an organized magnetic field of strength 7 mu G has profound consequences on the morphology of the wind bubble of massive stars at rest. In this paper, we show through 2.5D magnetohydrodynamical simulations, in the context of a Wolf-Rayet-evolving 35 M 0 star, that it affects the development of its supernova remnant. When the supernova remnant reaches its middle age (15-20 kyr), it adopts a tubular shape that results from the interaction between the isotropic supernova ejecta and the anisotropic, magnetized, shocked stellar progenitor bubble into which the supernova blast wave expands. Our calculations for non-thermal emission, i.e. radio synchrotron and inverse-Compton radiation, reveal that such supernova remnants can, due to projection effects, appear as rectangular objects in certain cases. This mechanism for shaping a supernova remnant is similar to the bipolar and elliptical planetary nebula production by wind-wind interaction in the low-mass regime of stellar evolution. If such a rectangular core-collapse supernova remnant is created, the progenitor star must not have been a runaway star. We propose that such a mechanism is at work in the shaping of the asymmetric core-collapse supernova remnant Puppis A.}, 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} }