@article{SushchBrosePohletal.2022, author = {Sushch, Iurii and Brose, Robert and Pohl, Martin and Plotko, Pavlo and Das, Samata}, title = {Leptonic nonthermal emission from supernova remnants evolving in the circumstellar magnetic field}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {926}, 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/ac3cb8}, pages = {14}, year = {2022}, abstract = {The very-high-energy (VHE; E > 100 GeV) gamma-ray emission observed from a number of supernova remnants (SNRs) indicates particle acceleration to high energies at the shock of the remnants and a potentially significant contribution to Galactic cosmic rays. It is extremely difficult to determine whether protons (through hadronic interactions and subsequent pion decay) or electrons (through inverse Compton scattering on ambient photon fields) are responsible for this emission. For a successful diagnostic, a good understanding of the spatial and energy distribution of the underlying particle population is crucial. Most SNRs are created in core-collapse explosions and expand into the wind bubble of their progenitor stars. This circumstellar medium features a complex spatial distribution of gas and magnetic field which naturally strongly affects the resulting particle population. In this work, we conduct a detailed study of the spectro-spatial evolution of the electrons accelerated at the forward shock of core-collapse SNRs and their nonthermal radiation, using the RATPaC code that is designed for the time- and spatially dependent treatment of particle acceleration at SNR shocks. We focus on the impact of the spatially inhomogeneous magnetic field through the efficiency of diffusion and synchrotron cooling. It is demonstrated that the structure of the circumstellar magnetic field can leave strong signatures in the spectrum and morphology of the resulting nonthermal emission.}, 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{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} } @misc{PetrukKuzyoOrlandoetal.2019, 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 = {Erratum: Post-adiabatic supernova remnants in an interstellar magnetic field: oblique shocks and non-uniform environment. - (Monthly notices of the Royal Astronomical Society. - 479, (2018), pg. 4253 - 4270)}, series = {Monthly notices of the Royal Astronomical Society}, volume = {482}, journal = {Monthly notices of the Royal Astronomical Society}, number = {2}, publisher = {Oxford Univ. Press}, address = {Oxford}, issn = {0035-8711}, doi = {10.1093/mnras/sty2861}, pages = {1979 -- 1980}, year = {2019}, abstract = {This is a correction notice for 'Post-adiabatic supernova remnants in an interstellar magnetic field: oblique shocks and non-uniform environment' (DOI: https://doi.org/10.1093/mnras/sty1750), which was published in MNRAS 479, 4253-4270 (2018). The publisher regrets to inform that the colour was missing from the colour scales in Figs 8(a)-(d) and Figs 9(a) and (b). This has now been corrected online. The publisher apologizes for this error.}, 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} } @phdthesis{Brose2020, author = {Brose, Robert}, title = {From dawn till dusk}, doi = {10.25932/publishup-47086}, url = {http://nbn-resolving.de/urn:nbn:de:kobv:517-opus4-470865}, school = {Universit{\"a}t Potsdam}, pages = {xiii, 146}, year = {2020}, abstract = {Supernova remnants are believed to be the source of cosmic rays with energies up to 10^15 eV that are produced within our Galaxy. The acceleration mechanism associated with the collision-less shocks in supernova remnants - diffusive shock acceleration - predicts a spectral index of the accelerated non-thermal particles of s = 2. However, measurements of non-thermal emission in radio, X-rays and gamma-rays reveal significant deviations of the particles spectral index from the canonical value of s = 2. 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 ≈100 yrs and inferred shock speed of ≈ 14, 000 km/s could make it an efficient particle accelerator. I performed spherical symmetric 1D simulations with the RATPaC code, in which I 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. The emission from G1.9+0.3 can be self-consistently explained within the test-particle limit. I 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 non-thermal 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. Although there are many models for the acceleration of cosmic rays in Supernova remnants, the escape of cosmic rays from these sources is yet understudied. I use our time-dependent acceleration code RATPaC to study the acceleration of cosmic rays and their escape in post-adiabatic Supernova remnants and calculate the subsequent gamma-ray emission from inverse-Compton scattering and Pion decay. My simulations span 100,000 years, thus covering the free-expansion, the Sedov-Taylor, and the beginning of the post-adiabatic phase of the remnant's evolution. At later stages of the evolution cosmic rays over a wide range of energy can reside outside of the remnant, creating spectra that are softer than predicted by standard diffusive shock acceleration and feature breaks in the 10 - 100 GeV-range. The total spectrum of cosmic rays released into the interstellar medium has a spectral index of s ≈ 2.4 above roughly 10 GeV which is close to that required by Galactic propagation models. I further find the gamma-ray luminosity to peak around an age of 4,000 years for inverse-Compton-dominated high-energy emission. Remnants expanding in low-density media emit generally more inverse-Compton radiation matching the fact that the brightest known supernova remnants - RCW86, Vela Jr, HESSJ1721-347 and RXJ1713.7-3946 - are all expanding in low density environments. The importance of feedback from the cosmic-rays on the hydrodynamical evolution of the remnants is debated as a possibility to obtain soft cosmic-ray spectra at low energies. I performed spherically symmetric 1-D simulations with a modified version of the RATPaC code, in which I simultaneously solve the transport equation for cosmic rays and the hydrodynamical equations, including the back-reaction of the cosmic-ray pressure on the flow profiles. Besides the known modification of the flow profiles and the consequently curved cosmic-ray spectra, steady-state models for non-linear diffusive shock acceleration overpredict the total compression ratio that can be reached with cosmic-ray feedback, as there is limited time for building these modifications. Further, I find modifications to the downstream flow structure that change the evolutionary behavior of the remnant and trigger a cosmic-ray-induced instability close to the contact discontinuity, if and when the cosmic-ray pressure becomes dominant there.}, language = {en} } @article{BenbowBirdBrilletal.2019, author = {Benbow, W. and Bird, R. and Brill, A. and Brose, Robert and Chromey, A. J. and Daniel, M. K. and Feng, Q. and Finley, J. P. and Fortson, L. and Furniss, A. and Gillanders, G. H. and Giuri, C. and Gueta, O. and Hanna, D. and Halpern, J. P. and Hassan, Tarek and Holder, J. and Hughes, G. and Humensky, T. B. and Joyce, Amy M. and Kaaret, P. and Kar, P. and Kelley-Hoskins, N. and Kertzman, M. and Kieda, D. and Krause, M. and Lang, M. J. and Lin, T. T. Y. and Maier, Gernot and Matthews, N. and Moriarty, P. and Mukherjee, R. and Nieto, D. and Nievas-Rosillos, M. and Ong, R. A. and Park, N. and Petrashyk, A. and Pohl, Martin and Pueschel, Elisa and Quinn, John and Ragan, K. and Reynolds, P. T. and Richards, G. T. and Roache, E. and Rulten, C. and Sadeh, Iftach and Santander, M. and Sembroski, G. H. and Shahinyan, K. and Sushch, Iurii and Wakely, S. P. and Wells, R. M. and Wilcox, P. and Wilhelm, Alina and Williams, David A. and Williamson, T. J.}, title = {Direct measurement of stellar angular diameters by the VERITAS Cherenkov telescopes}, series = {Nature astronomy}, volume = {3}, journal = {Nature astronomy}, number = {6}, publisher = {Nature Publ. Group}, address = {London}, issn = {2397-3366}, doi = {10.1038/s41550-019-0741-z}, pages = {511 -- 516}, year = {2019}, abstract = {The angular size of a star is a critical factor in determining its basic properties1. Direct measurement of stellar angular diameters is difficult: at interstellar distances stars are generally too small to resolve by any individual imaging telescope. This fundamental limitation can be overcome by studying the diffraction pattern in the shadow cast when an asteroid occults a star2, but only when the photometric uncertainty is smaller than the noise added by atmospheric scintillation3. Atmospheric Cherenkov telescopes used for particle astrophysics observations have not generally been exploited for optical astronomy due to the modest optical quality of the mirror surface. However, their large mirror area makes them well suited for such high-time-resolution precision photometry measurements4. Here we report two occultations of stars observed by the Very Energetic Radiation Imaging Telescope Array System (VERITAS)5 Cherenkov telescopes with millisecond sampling, from which we are able to provide a direct measurement of the occulted stars' angular diameter at the ≤0.1 mas scale. This is a resolution never achieved before with optical measurements and represents an order of magnitude improvement over the equivalent lunar occultation method6. We compare the resulting stellar radius with empirically derived estimates from temperature and brightness measurements, confirming the latter can be biased for stars with ambiguous stellar classifications.}, language = {en} } @article{ArcherBenbowBirdetal.2018, author = {Archer, A. and Benbow, Wystan and Bird, Ralph and Brose, Robert and Buchovecky, M. and Bugaev, V and Cui, Wei and Danie, M. K. and Falcone, A. and Feng, Qi and Finley, John P. and Flinders, A. and Fortson, L. and Furniss, Amy and Gillanders, Gerard H. and Huttens, M. and Hanna, David and Hervet, O. and Holder, J. and Hughes, G. and Humensky, T. B. and Johnson, Caitlin A. and Kaaret, Philip and Kar, P. and Kelley-Hoskins, N. and Kieda, David and Krause, Maria and Krennrich, F. and Kumar, S. and Lang, M. J. and Lin, T. T. Y. and McArthur, S. and Moriarty, P. and Mukherjee, Reshmi and Nieto, Daniel and Ong, R. A. and Otte, A. N. and Park, Nahee and Petrashyk, A. and Pohl, Martin and Popkow, Alexis and Pueschel, Elisa and Quinn, J. and Ragan, K. and Reynold, P. T. and Richards, Gregory T. and Roache, E. and Rulten, C. and Sadeh, I and Sembroski, G. H. and Shahinyan, Karlen and Tyler, J. and Wakely, S. P. and Weiner, O. M. and Weinstein, A. and Wells, R. M. and Wilcox, P. and Wilhelm, Alina and Williams, David A. and Brisken, W. F. and Pontrelli, P.}, title = {HESS J1943+213}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {862}, journal = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, number = {1}, publisher = {IOP Publ. Ltd.}, address = {Bristol}, organization = {VERITAS Collaboration}, issn = {0004-637X}, doi = {10.3847/1538-4357/aacbd0}, pages = {15}, year = {2018}, abstract = {HESS J1943+213 is a very high energy (VHE; > 100 GeV) gamma-ray source in the direction of the Galactic plane. Studies exploring the classification of the source are converging toward its identification as an extreme synchrotron BL Lac object. Here we present 38 hr of VERITAS observations of HESS J1943+213 taken over 2 yr. The source is detected with a significance of similar to 20 standard deviations, showing a remarkably stable flux and spectrum in VHE gamma-rays. Multifrequency Very Long Baseline Array (VLBA) observations of the source confirm the extended, jet-like structure previously found in the 1.6 GHz band with the European VLBI Network and detect this component in the 4.6 and 7.3 GHz bands. The radio spectral indices of the core and the jet and the level of polarization derived from the VLBA observations are in a range typical for blazars. Data from VERITAS, Fermi-LAT, Swift-XRT, the FLWO 48 ' telescope, and archival infrared and hard X-ray observations are used to construct and model the spectral energy distribution (SED) of the source with a synchrotron self-Compton model. The well-measured gamma-ray peak of the SED with VERITAS and Fermi-LAT provides constraining upper limits on the source redshift. Possible contribution of secondary gamma-rays from ultra-high-energy cosmic-ray-initiated electromagnetic cascades to the gamma-ray emission is explored, finding that only a segment of the VHE spectrum can be accommodated with this process. A variability search is performed across X-ray and gamma-ray bands. No statistically significant flux or spectral variability is detected.}, language = {en} } @article{ArcherBenbowBirdetal.2019, author = {Archer, A. and Benbow, Wystan and Bird, Ralph and Brose, Robert and Buchovecky, M. and Buckley, J. H. and Chromey, A. J. and Cui, Wei and Falcone, A. and Feng, Qi and Finley, J. P. and Fortson, Lucy and Furniss, Amy and Gent, A. and Gueta, O. and Hanna, David and Hassan, T. and Hervet, Olivier and Holder, J. and Hughes, G. and Humensky, T. B. and Johnson, Caitlin A. and Kaaret, Philip and Kar, P. and Kelley-Hoskins, N. and Kertzman, M. and Kieda, David and Krennrich, F. and Kumar, S. and Lang, M. J. and Lin, T. T. Y. and McCann, A. and Moriarty, P. and Mukherjee, Reshmi and Ong, R. A. and Otte, Adam Nepomuk and Pandel, D. and Park, N. and Petrashyk, A. and Pohl, Martin and Pueschel, Elisa and Quinn, J. and Ragan, K. and Richards, Gregory T. and Roache, E. and Sadeh, I and Santander, Marcos and Scott, S. S. and Sembroski, G. H. and Shahinyan, Karlen and Sushch, Iurii and Tyler, J. and Wakely, S. P. and Weinstein, A. and Wells, R. M. and Wilcox, P. and Wilhelm, Alina and Williams, D. A. and Williamson, T. J. and Zitzer, B.}, title = {A Search for Pulsed Very High-energy Gamma-Rays from 13 Young Pulsars in Archival VERITAS Data}, series = {The astrophysical journal : an international review of spectroscopy and astronomical physics}, volume = {876}, 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/ab14f4}, pages = {14}, year = {2019}, abstract = {We conduct a search for periodic emission in the very high-energy (VHE) gamma-ray band (E > 100 GeV) from a total of 13 pulsars in an archival VERITAS data set with a total exposure of over 450 hr. The set of pulsars includes many of the brightest young gamma-ray pulsars visible in the Northern Hemisphere. The data analysis resulted in nondetections of pulsed VHE gamma-rays from each pulsar. Upper limits on a potential VHE gamma-ray flux are derived at the 95\% confidence level above three energy thresholds using two methods. These are the first such searches for pulsed VHE emission from each of the pulsars, and the obtained limits constrain a possible flux component manifesting at VHEs as is seen for the Crab pulsar.}, language = {en} }