@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} } @article{HollerSchoeckEgeretal.2012, author = {Holler, M. and Schoeck, F. M. and Eger, P. and Kiessling, D. and Valerius, K. and Stegmann, Christian}, title = {Spatially resolved X-ray spectroscopy and modeling of the nonthermal emission of the pulsar wind nebula in G0.9+0.1}, series = {ASTRONOMY \& ASTROPHYSICS}, volume = {539}, journal = {ASTRONOMY \& ASTROPHYSICS}, number = {2}, publisher = {EDP SCIENCES S A}, address = {LES ULIS CEDEX A}, issn = {0004-6361}, doi = {10.1051/0004-6361/201118121}, pages = {8}, year = {2012}, abstract = {Aims. We performed a spatially resolved spectral X-ray study of the pulsar wind nebula ( PWN) in the supernova remnant G0.9+ 0.1. Furthermore, we modeled its nonthermal emission in the X-ray and very high-energy (VHE, E > 100 GeV) gamma-ray regime. Methods. Using Chandra ACIS-S3 data, we investigated the east-west dependence of the spectral properties of G0.9+ 0.1 by calculating hardness ratios. We analyzed the EPIC-MOS and EPIC-pn data of two on-axis observations of the XMM-Newton telescope and extracted spectra of four annulus-shaped regions, centered on the region of brightest emission of the source. A radially symmetric leptonic model was applied in order to reproduce the observed X-ray emission of the inner part of the PWN. Using the optimized model parameter values obtained from the X-ray analysis, we then compared the modeled inverse Compton (IC) radiation with the published H.E.S.S. gamma-ray data. Results. The spectral index within the four annuli increases with growing distance to the pulsar, whereas the surface brightness drops. With the adopted model we are able to reproduce the characteristics of the X-ray spectra. The model results for the VHE. radiation, however, strongly deviate from the H.E.S.S. data.}, language = {en} } @article{HenaultBrunetOskinovaGuerreroetal.2012, author = {Henault-Brunet, V. and Oskinova, Lida and Guerrero, Mart{\´i}n A. and Sun, W. and Chu, Y.-H. and Evans, C. J. and Gallagher, J. S. and Gruendl, R. A. and Reyes-Iturbide, J.}, title = {Discovery of a Be/X-ray pulsar binary and associated supernova remnant in the wing of the small magellanic cloud}, series = {Monthly notices of the Royal Astronomical Society}, volume = {420}, journal = {Monthly notices of the Royal Astronomical Society}, number = {1}, publisher = {Wiley-Blackwell}, address = {Malden}, issn = {0035-8711}, doi = {10.1111/j.1745-3933.2011.01183.x}, pages = {L13 -- L17}, year = {2012}, abstract = {We report on a new Be/X-ray pulsar binary located in the Wing of the Small Magellanic Cloud (SMC). The strong pulsed X-ray source was discovered with the Chandra and XMM-Newton X-ray observatories. The X-ray pulse period of 1062 s is consistently determined from both Chandra and XMM-Newton observations, revealing one of the slowest rotating X-ray pulsars known in the SMC. The optical counterpart of the X-ray source is the emission-line star 2dFS 3831. Its B0-0.5(III)e+ spectral type is determined from VLT-FLAMES and 2dF optical spectroscopy, establishing the system as a Be/X-ray binary (Be-XRB). The hard X-ray spectrum is well fitted by a power law with additional thermal and blackbody components, the latter reminiscent of persistent Be-XRBs. This system is the first evidence of a recent supernova in the low-density surroundings of NGC602. We detect a shell nebula around 2dFS 3831 in H alpha and [OIII] images and conclude that it is most likely a supernova remnant. If it is linked to the supernova explosion that created this new X-ray pulsar, its kinematic age of (2-4) x 10(4) yr provides a constraint on the age of the pulsar.}, language = {en} } @article{GvaramadzeKniazevOskinova2019, author = {Gvaramadze, V. V. and Kniazev, Alexei Y. and Oskinova, Lida}, title = {Discovery of a putative supernova remnant around the long-period X-ray pulsar SXP 1323 in the Small Magellanic Cloud}, series = {Monthly notices of the Royal Astronomical Society}, volume = {485}, journal = {Monthly notices of the Royal Astronomical Society}, number = {1}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnrasl/slz018}, pages = {L6 -- L10}, year = {2019}, abstract = {We report the discovery of a circular shell centred on the Be X-ray binary (BeXB) SXP 1323 in the Small Magellanic Cloud. The shell was detected in an H α image obtained with the Very Large Telescope. Follow-up spectroscopy with the Southern African Large Telescope showed that the shell expands with a velocity of ≈100kms-1 and that its emission is due to shock excitation. We suggest that this shell is a remnant of the supernova explosion that led to the formation of SXP 1323's neutron star ≈40000 yr ago. SXP 1323 represents the second known case of a BeXB associated with a supernova remnant (the first one is SXP 1062). Interestingly, both of these BeXBs harbour long-period pulsars and are located in a low-metallicity galaxy.}, 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{FedericiPohlTelezhinskyetal.2015, author = {Federici, Simone and Pohl, Martin and Telezhinsky, Igor O. and Wilhelm, Alina and Dwarkadas, Vikram V.}, title = {Analysis of GeV-band gamma-ray emission from supernova remnant RX J1713.7-3946}, series = {Astronomy and astrophysics : an international weekly journal}, volume = {577}, journal = {Astronomy and astrophysics : an international weekly journal}, publisher = {EDP Sciences}, address = {Les Ulis}, issn = {0004-6361}, doi = {10.1051/0004-6361/201424947}, pages = {9}, year = {2015}, abstract = {Context. RX J1713.7-3946 is the brightest shell-type supernova remnant (SNR) of the TeV gamma-ray sky. Earlier Fermi-LAT results on low energy gamma-ray emission suggested that, despite large uncertainties in the background determination, the spectrum is inconsistent with a hadronic origin. Aims. We update the GeV-band spectra using improved estimates for the diffuse Galactic gamma-ray emission and more than double the volume of data. We further investigate the viability of hadronic emission models for RX J1713.7-3946. Methods. We produced a high-resolution map of the diffuse Galactic gamma-ray background corrected for the HI self-absorption and used it in the analysis of more than five years worth of Fermi-LAT data. We used hydrodynamic scaling relations and a kinetic transport equation to calculate the acceleration and propagation of cosmic rays in SNR. We then determined spectra of hadronic gamma-ray emission from RX J1713.7-3946, separately for the SNR interior and the cosmic-ray precursor region of the forward shock, and computed flux variations that would allow us to test the model with observations. Results. We find that RX J1713.7-3946 is now detected by Fermi-LAT with very high statistical significance, and the source morphology is best described by that seen in the TeV band. The measured spectrum of RX J1713.7-3946 is hard with index gamma = 1.53 +/- 0.07, and the integral flux above 500 MeV is F = (5 : 5 +/- 1 : 1) x 10(-9) photons cm(-2) s(-1). We demonstrate that scenarios based on hadronic emission from the cosmic-ray precursor region are acceptable for RX J1713.7-3946, and we predict a secular flux increase at a few hundred GeV at the level of around 15\% over ten years, which may be detectable with the upcoming Cherenkov Telescope Array (CTA) observatory.}, 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{AliuArchambaultAuneetal.2014, author = {Aliu, E. and Archambault, S. and Aune, T. and Behera, B. and Beilicke, M. and Benbow, W. and Berger, K. and Bird, R. and Buckley, J. H. and Bugaev, V. and Cardenzana, J. V. and Cerruti, M. and Chen, Xuhui and Ciupik, L. and Collins-Hughes, E. and Connolly, M. P. and Cui, W. and Dumm, J. and Dwarkadas, Vikram V. and Errando, M. and Falcone, A. and Federici, Simone and Feng, Q. and Finley, J. P. and Fleischhack, H. and Fortin, P. and Fortson, L. and Furniss, A. and Galante, N. and Gall, D. 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 Kumar, S. and Lang, M. J. and Madhavan, A. S. and Maier, G. and McCann, A. J. and Meagher, K. and Millis, J. and Moriarty, P. and Mukherjee, R. and Nieto, D. and Ong, R. A. and Otte, A. N. and Pandel, D. and Park, N. and Pohl, Martin and Popkow, A. and Prokoph, H. and Quinn, J. and Ragan, K. and Rajotte, J. and Ratliff, G. and Reyes, L. C. and Reynolds, P. T. and Richards, G. T. and Roache, E. and Rousselle, J. and Sembroski, G. H. and Shahinyan, K. and Sheidaei, F. and Smith, A. W. and Staszak, D. and Telezhinsky, Igor O. and Tsurusaki, K. and Tucci, J. V. and Tyler, J. and Varlotta, A. and Vassiliev, V. V. and Vincent, S. and Wakely, S. P. and Ward, J. E. and Weinstein, A. and Welsing, R. and Wilhelm, Alina}, title = {Investigating the TeV morpholoy of MGRO J1908+06 with veritas}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {787}, 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/787/2/166}, pages = {7}, year = {2014}, abstract = {We report on deep observations of the extended TeV gamma-ray source MGRO J1908+06 made with the VERITAS very high energy gamma-ray observatory. Previously, the TeV emission has been attributed to the pulsar wind nebula (PWN) of the Fermi-LAT pulsar PSR J1907+0602. We detect MGRO J1908+06 at a significance level of 14 standard deviations (14 sigma) and measure a photon index of 2.20 +/- 0.10(stat) +/- 0.20(sys). The TeV emission is extended, covering the region near PSR J1907+0602 and also extending toward SNR G40.5-0.5. When fitted with a two-dimensional Gaussian, the intrinsic extension has a standard deviation of sigma(src) = 0 degrees.44 +/- 0 degrees.02. In contrast to other TeV PWNe of similar age in which the TeV spectrum softens with distance from the pulsar, the TeV spectrum measured near the pulsar location is consistent with that measured at a position near the rim of G40.5-0.5, 0 degrees.33 away.}, language = {en} }