Showing 1–50 of 142 results for author: Khangulyan, D

H.E.S.S. observations of the 2021 periastron passage of PSR B1259-63/LS 2883

Authors: H. E. S. S. Collaboration, F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, J. Borowska, M. Bouyahiaoui, R. Brose, A. Brown, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin, S. Caroff, S. Casanova , et al. (119 additional authors not shown)

Abstract: PSR B1259-63 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 year period, around an O9.5Ve star. At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of ~100 hours of data taken from $t_p-24$ day… ▽ More PSR B1259-63 is a gamma-ray binary system that hosts a pulsar in an eccentric orbit, with a 3.4 year period, around an O9.5Ve star. At orbital phases close to periastron passages, the system radiates bright and variable non-thermal emission. We report on an extensive VHE observation campaign conducted with the High Energy Stereoscopic System, comprised of ~100 hours of data taken from $t_p-24$ days to $t_p+127$ days around the system's 2021 periastron passage. We also present the timing and spectral analyses of the source. The VHE light curve in 2021 is consistent with the stacked light curve of all previous observations. Within the light curve, we report a VHE maximum at times coincident with the third X-ray peak first detected in the 2021 X-ray light curve. In the light curve -- although sparsely sampled in this time period -- we see no VHE enhancement during the second disc crossing. In addition, we see no correspondence to the 2021 GeV flare in the VHE light curve. The VHE spectrum obtained from the analysis of the 2021 dataset is best described by a power law of spectral index $Γ= 2.65 \pm 0.04_{\text{stat}}$ $\pm 0.04_{\text{sys}}$, a value consistent with the previous H.E.S.S. observations of the source. We report spectral variability with a difference of $ΔΓ= 0.56 ~\pm~ 0.18_{\text{stat}}$ $~\pm~0.10_{\text{sys}}$ at 95% c.l., between sub-periods of the 2021 dataset. We also find a linear correlation between contemporaneous flux values of X-ray and TeV datasets, detected mainly after $t_p+25$ days, suggesting a change in the available energy for non-thermal radiation processes. We detect no significant correlation between GeV and TeV flux points, within the uncertainties of the measurements, from $\sim t_p-23$ days to $\sim t_p+126$ days. This suggests that the GeV and TeV emission originate from different electron populations. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: accepted to A&A

Upper Limit on the Coronal Cosmic Ray Energy Budget in Seyfert Galaxies

Authors: Yoshiyuki Inoue, Shinsuke Takasao, Dmitry Khangulyan

Abstract: The IceCube collaboration has reported possible detections of high-energy neutrinos from nearby Seyfert galaxies. While central hot coronae are proposed as the primary neutrino production site, the exact coronal cosmic-ray energy budget has been loosely constrained. In this study, we propose a new stringent upper bound on the coronal cosmic-ray energy budget of Seyfert galaxies, considering both a… ▽ More The IceCube collaboration has reported possible detections of high-energy neutrinos from nearby Seyfert galaxies. While central hot coronae are proposed as the primary neutrino production site, the exact coronal cosmic-ray energy budget has been loosely constrained. In this study, we propose a new stringent upper bound on the coronal cosmic-ray energy budget of Seyfert galaxies, considering both accretion dynamics and observed properties of radio-quiet Seyfert galaxies. Notably, even under the calorimetric condition where cosmic rays lose all their energies, our limit indicates that the coronal neutrino flux of NGC~1068 is about an order of magnitude fainter than the observed levels. This discrepancy suggests the need for further theoretical and observational investigations on the IceCube signals from Seyfert galaxies. △ Less

Submitted 6 July, 2024; v1 submitted 15 January, 2024; originally announced January 2024.

Comments: 8 pages, no figures, accepted for publication in PASJ

Report number: RIKEN-iTHEMS-Report-24

Naked forward shock seen in the TeV afterglow data of GRB221009A

Authors: Dmitry Khangulyan, Felix Aharonian, Andrew M. Taylor

Abstract: We explore the implications of the light curve of the early TeV gamma-ray afterglow of GRB221009A reported by the LHAASO collaboration. We show that the reported offset of the reference time, $T_*$, allows the determination of the relativistic jet activation time, which occurs approximately $200\,\mathrm{s}$ after the GBM trigger time and closely precedes the moment at which GBM was saturated. We… ▽ More We explore the implications of the light curve of the early TeV gamma-ray afterglow of GRB221009A reported by the LHAASO collaboration. We show that the reported offset of the reference time, $T_*$, allows the determination of the relativistic jet activation time, which occurs approximately $200\,\mathrm{s}$ after the GBM trigger time and closely precedes the moment at which GBM was saturated. We find that while the LHAASO data do not exclude the homogeneous circumburst medium scenario, the progenitor wind scenario looks preferable, finding excellent agreement with the expected size of the stellar bubble. We conclude that the initial growth of the light curve is dominated by processes internal to the jet or by gamma-gamma attenuation on the photons emitted during the prompt phase. Namely, either the activation of the acceleration process or the decrease of internal gamma-gamma absorption can naturally explain the initial rapid flux increase. The subsequent slow flux growth phase observed up to $T_*+18\,\mathrm{s}$ is explained by the build-up of the synchrotron radiation -- the target for inverse Compton scattering, which is also supported by a softer TeV spectrum measured during this period. The duration of this phase allows an almost parameter-independent determination of the jet's initial Lorentz factor, $Γ_0\approx600$, and magnetic field strength, $B'\sim0.3\,\mathrm{G}$. These values appear to match well those previously revealed through spectral modeling of the GRB emission. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 20 pages, 6 figures, submitted

On the properties of inverse Compton spectra generated by up-scattering a power-law distribution of target photons

Authors: Dmitry Khangulyan, Felix Aharonian, Andrew M. Taylor

Abstract: Relativistic electrons are an essential component in many astrophysical sources, and their radiation may dominate the high-energy bands. Inverse Compton (IC) emission is the radiation mechanism that plays the most important role in these bands. The basic properties of IC, such as the total and differential cross sections, have long been studied; the properties of the IC emission depend strongly no… ▽ More Relativistic electrons are an essential component in many astrophysical sources, and their radiation may dominate the high-energy bands. Inverse Compton (IC) emission is the radiation mechanism that plays the most important role in these bands. The basic properties of IC, such as the total and differential cross sections, have long been studied; the properties of the IC emission depend strongly not only on the emitting electron distribution but also on the properties of the target photons. This complicates the phenomenological studies of sources, where target photons are supplied from a broad radiation component. We study the spectral properties of IC emission generated by a power-law distribution of electrons on a power-law distribution of target photons. We approximate the resulting spectrum by a broken-power-law distribution and show that there can be up to three physically motivated spectral breaks. If the target photon spectrum extends to sufficiently low energies, $\varepsilon_{\mathrm{min}}< m_e^2c^4/E_{\mathrm{max}}$ ($m_e$ and $c$ are electron mass and speed of light, respectively; $\varepsilon_{\mathrm{min}}$ and $E_{\mathrm{max}}$ are the minimum/maximum energies of target photons and electrons, respectively), then the high energy part of the IC component has a spectral slope typical for the Thomson regime with an abrupt cutoff close to $E_{\mathrm{max}}$. The spectra typical for the Klein-Nishina regime are formed above $m_e^2c^4/\varepsilon_{\mathrm{min}}$. If the spectrum of target photons features a cooling break, i.e., a change of the photon index by $0.5$ at $\varepsilon_{\mathrm{br}}$, then the transition to the Klein-Nishina regime proceeds through an intermediate change of the photon index by $0.5$ at $m_e^2c^4/\varepsilon_{\mathrm{br}}$. △ Less

Submitted 23 July, 2023; originally announced July 2023.

Comments: 12 pages, 7 figures

Gamma-Ray Emission in the Seyfert Galaxy NGC 4151: Investigating the Role of Jet and Coronal Activities

Authors: Yoshiyuki Inoue, Dmitry Khangulyan

Abstract: NGC 4151, a nearby Seyfert galaxy, has recently been reported to emit gamma rays in the GeV range, posing an intriguing astrophysical mystery. The star formation rate of NGC 4151 is too low to explain the observed GeV flux, but the galaxy is known for its coronal activity in X-ray and jet activity in radio. We propose that either the combination of these two activities or the jet activity alone ca… ▽ More NGC 4151, a nearby Seyfert galaxy, has recently been reported to emit gamma rays in the GeV range, posing an intriguing astrophysical mystery. The star formation rate of NGC 4151 is too low to explain the observed GeV flux, but the galaxy is known for its coronal activity in X-ray and jet activity in radio. We propose that either the combination of these two activities or the jet activity alone can account for the gamma-ray spectrum. An energy-dependent variability search will allow one to distinguish between the two scenarios, as the coronal component can only contribute at energies of $\lesssim1$ GeV. Our analysis also indicates that it might still be difficult to see coronal neutrinos from the apparently X-ray brightest Seyfert NGC 4151 with current-generation neutrino observatories. △ Less

Submitted 16 January, 2024; v1 submitted 8 April, 2023; originally announced April 2023.

Comments: 7 pages, 4 figures, published in PASJ

Report number: RIKEN-iTHEMS-Report-23

Unveiling properties of the non-thermal X-ray production in the gamma-ray binary LS 5039 using the long-term pattern of its fast X-ray variability

Authors: Hiroki Yoneda, Valenti Bosch-Ramon, Teruaki Enoto, Dmitry Khangulyan, Paul S. Ray, Tod Strohmayer, Toru Tamagawa, Zorawar Wadiasingh

Abstract: Gamma-ray binary systems, a subclass of high-mass X-ray binaries, show non-thermal emissions from radio to TeV. While efficient electron acceleration is considered to take place in them, the nature of the acceleration mechanism and the physical environments in these systems have been a long-standing question. In this work, we report on long-term recurrent patterns in the short-term variability of… ▽ More Gamma-ray binary systems, a subclass of high-mass X-ray binaries, show non-thermal emissions from radio to TeV. While efficient electron acceleration is considered to take place in them, the nature of the acceleration mechanism and the physical environments in these systems have been a long-standing question. In this work, we report on long-term recurrent patterns in the short-term variability of the soft X-ray emission of LS 5039, one of the brightest gamma-ray binary systems. The Neutron star Interior Composition Explorer (NICER) observed LS 5039 four times from 2018 to 2021. By comparing them with the previous Suzaku and NuSTAR long-exposure observations, we studied the long-term evolution of the orbital light curve in the soft X-ray band. Although the observations by NICER and Suzaku are separated by $\sim$14 years, i.e., more than 10^3 orbits, the orbital light curves show remarkable consistency after calculating their running averages with a window width 70 ks. Furthermore, all of the light curves show short-term variability with a time scale of $\sim$10 ks. Since the column density did not vary when the flux changed abruptly, such a short-term variability seems to be an intrinsic feature of the X-ray emission. We propose that the short-term variability is caused by clumps (or inhomogeneities) of the companion star wind impacting the X-ray production site. The observed time scale matches well with the lifetime of the clumps interacting with the pulsar wind and the dynamical time scale of the relativistic intrabinary shock in the pulsar wind scenario. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: 17 pages, 13 figures, accepted for publication in ApJ

The formation of hard VHE spectra from GRB afterglow via Two-Zone Synchrotron Self-Compton Emission

Authors: Dmitry Khangulyan, Andrew M. Taylor, Felix Aharonian

Abstract: Electron Compton scattering of target photons into the gamma-ray energy band (inverse Compton scattering --IC--) is commonly expected to dominate the very high energy spectra in gamma-ray bursts especially during the afterglow phase. For sufficiently large center-of-mass energies in these collisions, the effect of the electron recoil starts reducing the scattering cross section (the Klein-Nishina… ▽ More Electron Compton scattering of target photons into the gamma-ray energy band (inverse Compton scattering --IC--) is commonly expected to dominate the very high energy spectra in gamma-ray bursts especially during the afterglow phase. For sufficiently large center-of-mass energies in these collisions, the effect of the electron recoil starts reducing the scattering cross section (the Klein-Nishina regime). The IC spectra generated in the Klein-Nishina regime is softer and has a smaller flux level compared to the synchrotron spectra produced by the same electrons. The detection of afterglow emission from nearby GRB 190819A in the very high energy (VHE) domain with H.E.S.S. has revealed an unexpected feature: the slope of the VHE spectrum matches well the slope of the X-ray spectra, despite expectations that for the IC production process, the impact of the Klein-Nishina effect should be strong. The multi-wavelength spectral energy distribution appears to be inconsistent with predictions of one-zone synchrotron-self-Compton models. We study the possible impact of two-zone configuration on the properties of IC emission when the magnetic field strength differs considerably between the two zones. Synchrotron photons from the strong magnetic field zone provide the dominant target for cooling of the electrons in the weak magnetic field zone, which results in a formation of hard electron distribution and consequently of a hard IC emission. We show that the two-zone model can provide a good description of the X-ray XRT and VHE H.E.S.S. data. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 13 pages, 7 figures, ApJ submitted

ALMA detection of parsec-scale blobs at the head of kiloparsec-scale jet in the nearby Seyfert galaxy NGC 1068

Authors: Tomonari Michiyama, Yoshiyuki Inoue, Akihiro Doi, Dmitry Khangulyan

Abstract: We present Atacama Large Millimeter/submillimeter Array observations at $\approx100$ GHz with $0.05$ arcsec (3 pc) resolution of the kiloparsec-scale jet seen in the nearby Seyfert galaxy NGC 1068, and we report the presence of parsec-scale blobs at the head of the jet. The combination of the detected radio flux ($\approx0.8$ mJy), spectral index ($\approx0.5$), and the blob size ($\approx10$ pc)… ▽ More We present Atacama Large Millimeter/submillimeter Array observations at $\approx100$ GHz with $0.05$ arcsec (3 pc) resolution of the kiloparsec-scale jet seen in the nearby Seyfert galaxy NGC 1068, and we report the presence of parsec-scale blobs at the head of the jet. The combination of the detected radio flux ($\approx0.8$ mJy), spectral index ($\approx0.5$), and the blob size ($\approx10$ pc) suggests a strong magnetic field of $B\approx240\,μ$G. Such a strong magnetic field most likely implies magnetic field amplification by streaming cosmic rays. The estimated cosmic-ray power by the jet may exceed the limit set by the star formation activity in this galaxy. This result suggests that even modest-power jets can increase the galactic cosmic-ray content while propagating through the galactic bulge. △ Less

Submitted 17 August, 2022; originally announced August 2022.

Comments: 9 pages, 3 figures, accepted for publication in ApJL

Report number: RIKEN-iTHEMS-Report-22

Spatially resolved study of the SS 433/W50 west region with Chandra: X-ray structure and spectral variation of non-thermal emission

Authors: Kazuho Kayama, Takaaki Tanaka, Hiroyuki Uchida, Takeshi Go Tsuru, Takahiro Sudoh, Yoshiyuki Inoue, Dmitry Khangulyan, Naomi Tsuji, Hiroaki Yamamoto

Abstract: The X-ray binary SS 433, embedded in the W50 nebula (or supernova remnant W50), shows bipolar jets that are ejected with mildly relativistic velocities, and extend toward the east and west out to scales of tens of parsecs. Previous X-ray observations revealed twin lobes along the jet precession axis that contain compact bright knots dominated by synchrotron radiation, which provide evidence of ele… ▽ More The X-ray binary SS 433, embedded in the W50 nebula (or supernova remnant W50), shows bipolar jets that are ejected with mildly relativistic velocities, and extend toward the east and west out to scales of tens of parsecs. Previous X-ray observations revealed twin lobes along the jet precession axis that contain compact bright knots dominated by synchrotron radiation, which provide evidence of electron acceleration in this system. Particle acceleration in this system is substantiated by the recently detected gamma rays with energies up to at least 25 TeV. To further elucidate the origin of the knots and particle acceleration sites in SS 433/W50, we report here on detailed, spatially resolved X-ray spectroscopy of its western lobe with Chandra. We detect synchrotron emission along the jet precession axis, as well as optically thin thermal emission that is more spatially extended. Between the two previously known knots, w1 and w2, we discover another synchrotron knot, which we call w1.5. We find no significant synchrotron emission between SS 433 and the innermost X-ray knot (w1), suggesting that electrons only begin to be accelerated at w1. The X-ray spectra become gradually steeper from w1 to w2, and then rapidly so immediately outside of w2. Comparing with a model taking into account electron transport and cooling along the jet, this result indicates that the magnetic field in w2 is substantially enhanced, which also explains its brightness. We discuss possible origins of the enhanced magnetic field of w2 as well as scenarios to explain the other two knots. △ Less

Submitted 12 July, 2022; originally announced July 2022.

Comments: 14 pages, 11 figures, Accepted to PASJ

Report number: RIKEN-iTHEMS-Report-22

Time-resolved hadronic particle acceleration in the recurrent Nova RS Ophiuchi

Authors: H. E. S. S. Collaboration, F. Aharonian, F. Ait Benkhali, E. O. Angüner, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun, S. Caroff, S. Casanova, M. Cerruti, T. Chand, A. Chen , et al. (150 additional authors not shown)

Abstract: Recurrent Novae are repeating thermonuclear explosions in the outer layers of white dwarfs, due to the accretion of fresh material from a binary companion. The shock generated by ejected material slamming into the companion star's wind, accelerates particles to very-high-energies. We report very-high-energy (VHE, $\gtrsim100$\,GeV) gamma rays from the recurrent nova RS\,Ophiuchi up to a month afte… ▽ More Recurrent Novae are repeating thermonuclear explosions in the outer layers of white dwarfs, due to the accretion of fresh material from a binary companion. The shock generated by ejected material slamming into the companion star's wind, accelerates particles to very-high-energies. We report very-high-energy (VHE, $\gtrsim100$\,GeV) gamma rays from the recurrent nova RS\,Ophiuchi up to a month after its 2021 outburst, using the High Energy Stereoscopic System. The VHE emission has a similar temporal profile to lower-energy GeV emission, indicating a common origin, with a two-day delay in peak flux. These observations constrain models of time-dependent particle energization, favouring a hadronic emission scenario over the leptonic alternative. This confirms that shocks in dense winds provide favourable environments for efficient cosmic-ray acceleration to very-high-energies. △ Less

Submitted 28 March, 2022; v1 submitted 16 February, 2022; originally announced February 2022.

Comments: Submitted 1st Nov. 2021, first release 10th March 2022 (accepted version)

Searching for TeV gamma-ray emission from SGR\,1935+2154 during its 2020 X-ray and radio bursting phase

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Anguner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlohr, B. Bi, M. Bottcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose , et al. (230 additional authors not shown)

Abstract: Magnetar hyperflares are the most plausible explanation for fast radio bursts (FRB) -- enigmatic powerful radio pulses with durations of several milliseconds and high brightness temperatures. The first observational evidence for this scenario was obtained in 2020 April when a FRB was detected from the direction of the Galactic magnetar and soft gamma-ray repeater SGR\,1935+2154. The FRB was preced… ▽ More Magnetar hyperflares are the most plausible explanation for fast radio bursts (FRB) -- enigmatic powerful radio pulses with durations of several milliseconds and high brightness temperatures. The first observational evidence for this scenario was obtained in 2020 April when a FRB was detected from the direction of the Galactic magnetar and soft gamma-ray repeater SGR\,1935+2154. The FRB was preceded by two gamma-ray outburst alerts by the BAT instrument aboard the Swift satellite, which triggered follow-up observations by the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. has observed SGR\,1935+2154 for 2 hr on 2020 April 28. The observations are coincident with X-ray bursts from the magnetar detected by INTEGRAL and Fermi-GBM, thus providing the first very high energy (VHE) gamma-ray observations of a magnetar in a flaring state. High-quality data acquired during these follow-up observations allow us to perform a search for short-time transients. No significant signal at energies $E>0.6$~TeV is found and upper limits on the persistent and transient emission are derived. We here present the analysis of these observations and discuss the obtained results and prospects of the H.E.S.S. follow-up program for soft gamma-ray repeaters. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: Received 2021 June 13; accepted 2021 June 28; published 2021 September 29

LMC N132D: A mature supernova remnant with a power-law gamma-ray spectrum extending beyond 8 TeV

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun , et al. (212 additional authors not shown)

Abstract: We analyzed 252 hours of High Energy Stereoscopic System (H.E.S.S.) observations towards the supernova remnant (SNR) LMC N132D that were accumulated between December 2004 and March 2016 during a deep survey of the Large Magellanic Cloud, adding 104 hours of observations to the previously published data set to ensure a > 5 sigma detection. To broaden the gamma-ray spectral coverage required for mod… ▽ More We analyzed 252 hours of High Energy Stereoscopic System (H.E.S.S.) observations towards the supernova remnant (SNR) LMC N132D that were accumulated between December 2004 and March 2016 during a deep survey of the Large Magellanic Cloud, adding 104 hours of observations to the previously published data set to ensure a > 5 sigma detection. To broaden the gamma-ray spectral coverage required for modeling the spectral energy distribution, an analysis of Fermi-LAT Pass 8 data was also included. We unambiguously detect N132D at very high energies (VHE) with a significance of 5.7 sigma. We report the results of a detailed analysis of its spectrum and localization based on the extended H.E.S.S. data set. The joint analysis of the extended H.E.S.S and Fermi-LAT data results in a spectral energy distribution in the energy range from 1.7 GeV to 14.8 TeV, which suggests a high luminosity of N132D at GeV and TeV energies. We set a lower limit on a gamma-ray cutoff energy of 8 TeV with a confidence level of 95%. The new gamma-ray spectrum as well as multiwavelength observations of N132D when compared to physical models suggests a hadronic origin of the VHE gamma-ray emission. SNR N132D is a VHE gamma-ray source that shows a spectrum extending to the VHE domain without a spectral cutoff at a few TeV, unlike the younger oxygen-rich SNR Cassiopeia A. The gamma-ray properties of N132D may be affected by an interaction with a nearby molecular cloud that partially lies inside the 95% confidence region of the source position. [Abridged] △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 655, A7 (2021)

Reconstruction of multiple Compton scattering events in MeV gamma-ray Compton telescopes towards GRAMS: the physics-based probabilistic model

Authors: Hiroki Yoneda, Hirokazu Odaka, Yuto Ichinohe, Satoshi Takashima, Tsuguo Aramaki, Kazutaka Aoyama, Jonathan Asaadi, Lorenzo Fabris, Yoshiyuki Inoue, Georgia Karagiorgi, Dmitry Khangulyan, Masato Kimura, Jonathan Leyva, Reshmi Mukherjee, Taichi Nakasone, Kerstin Perez, Mayu Sakurai, William Seligman, Masashi Tanaka, Naomi Tsuji, Kohei Yorita, Jiancheng Zeng

Abstract: Aimed at progress in mega-electron volt (MeV) gamma-ray astronomy, which has not yet been well-explored, Compton telescope missions with a variety of detector concepts have been proposed so far. One of the key techniques for these future missions is an event reconstruction algorithm that is able to determine the scattering orders of multiple Compton scattering events and to identify events in whic… ▽ More Aimed at progress in mega-electron volt (MeV) gamma-ray astronomy, which has not yet been well-explored, Compton telescope missions with a variety of detector concepts have been proposed so far. One of the key techniques for these future missions is an event reconstruction algorithm that is able to determine the scattering orders of multiple Compton scattering events and to identify events in which gamma rays escape from the detectors before they deposit all of their energies. We revisit previous event reconstruction methods and propose a modified algorithm based on a probabilistic method. First, we present a general formalism of the probabilistic model of Compton scattering describing physical interactions inside the detector and measurement processes. Then, we also introduce several approximations in the calculation of the probability functions for efficient computation. For validation, the developed algorithm has been applied to simulation data of a Compton telescope using a liquid argon time projection chamber, which is a new type of Compton telescope proposed for the GRAMS project. We have confirmed that it works successfully for up to 8-hit events, including correction of incoming gamma-ray energies for escape events. The proposed algorithm can be used for next-generation MeV gamma-ray missions featured by large-volume detectors, e.g., GRAMS. △ Less

Submitted 10 August, 2022; v1 submitted 5 July, 2021; originally announced July 2021.

Comments: 28 pages, 26 figures, accepted for publication in Astroparticle Physics

TeV emission of Galactic plane sources with HAWC and H.E.S.S

Authors: H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun, P. Brun , et al. (299 additional authors not shown)

Abstract: The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy gamma-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their… ▽ More The High Altitude Water Cherenkov (HAWC) observatory and the High Energy Stereoscopic System (H.E.S.S.) are two leading instruments in the ground-based very-high-energy gamma-ray domain. HAWC employs the water Cherenkov detection (WCD) technique, while H.E.S.S. is an array of Imaging Atmospheric Cherenkov Telescopes (IACTs). The two facilities therefore differ in multiple aspects, including their observation strategy, the size of their field of view and their angular resolution, leading to different analysis approaches. Until now, it has been unclear if the results of observations by both types of instruments are consistent: several of the recently discovered HAWC sources have been followed up by IACTs, resulting in a confirmed detection only in a minority of cases. With this paper, we go further and try to resolve the tensions between previous results by performing a new analysis of the H.E.S.S. Galactic plane survey data, applying an analysis technique comparable between H.E.S.S. and HAWC. Events above 1 TeV are selected for both datasets, the point spread function of H.E.S.S. is broadened to approach that of HAWC, and a similar background estimation method is used. This is the first detailed comparison of the Galactic plane observed by both instruments. H.E.S.S. can confirm the gamma-ray emission of four HAWC sources among seven previously undetected by IACTs, while the three others have measured fluxes below the sensitivity of the H.E.S.S. dataset. Remaining differences in the overall gamma-ray flux can be explained by the systematic uncertainties. Therefore, we confirm a consistent view of the gamma-ray sky between WCD and IACT techniques. △ Less

Submitted 8 September, 2021; v1 submitted 3 July, 2021; originally announced July 2021.

Broad-Band High-Energy Emission of the Gamma-ray Binary System LS 5039: Spectral and Temporal Features using NuSTAR and Fermi Observations

Authors: Hiroki Yoneda, Dmitry Khangulyan, Teruaki Enoto, Kazuo Makishima, Kairi Mine, Tsunefumi Mizuno, Tadayuki Takahashi

Abstract: We report on detailed analysis of the hard X-ray and GeV gamma-ray spectra of LS 5039, one of the brightest gamma-ray binary system in the Galaxy. The NuSTAR observation covering its entire orbit in 2016 allowed us for the first time to study the orbital variability of the spectrum above 10 keV. The hard X-ray spectrum is well described with a single power-law component up to 78 keV. The X-ray flu… ▽ More We report on detailed analysis of the hard X-ray and GeV gamma-ray spectra of LS 5039, one of the brightest gamma-ray binary system in the Galaxy. The NuSTAR observation covering its entire orbit in 2016 allowed us for the first time to study the orbital variability of the spectrum above 10 keV. The hard X-ray spectrum is well described with a single power-law component up to 78 keV. The X-ray flux showed a slight deviation from those observed previously with Suzaku in 2007. The fast X-ray brightening observed with Suzaku, around the inferior conjunction, was not observed in this observation. We also analyzed 11 years of Fermi Large Area Telescope data of LS 5039. The GeV spectrum around the inferior conjunction was well described with two non-thermal components; a power law with a photon index of $\sim 3$ and a cut-off power law with a cutoff energy of $\sim 2$ GeV. The orbital flux variability also changed gradually around a few GeV. These results indicate that there are two emission components in the GeV band, and the dominant component above $\sim 1$ GeV does not depend on the orbital phase. By combining these results, we update the spectral energy distribution of LS 5039 with the highest available statistics. Theoretical models proposed so far cannot explain the obtained multi-wavelength spectrum, especially the emission from $\sim$ 1 MeV to $\sim$ 400 MeV, and we discuss a possibility that particle acceleration in LS 5039 is different from the shock acceleration. △ Less

Submitted 18 October, 2021; v1 submitted 24 June, 2021; originally announced June 2021.

Comments: 19 pages, 15 figures, 8 tables, accepted for publication in ApJ

Fast radio bursts by high-frequency synchrotron maser emission generated at the reverse shock of a powerful magnetar flare

Authors: Dmitry Khangulyan, Maxim V. Barkov, Sergey B. Popov

Abstract: We consider a magnetar flare model for fast radio bursts (FRBs). We show that millisecond burst of sufficient power can be generated by synchrotron maser emission ignited at the reverse shock propagating through the weakly magnetized material that forms the magnetar flare. If the maser emission is generated in an anisotropic regime (due to the geometry of the production region or presence of an in… ▽ More We consider a magnetar flare model for fast radio bursts (FRBs). We show that millisecond burst of sufficient power can be generated by synchrotron maser emission ignited at the reverse shock propagating through the weakly magnetized material that forms the magnetar flare. If the maser emission is generated in an anisotropic regime (due to the geometry of the production region or presence of an intense external source of stimulating photons) the duration of the maser flashes is similar to the magnetar flare duration even if the shock front radius is large. Our scenario allows relaxing the requirements for several key parameters: the magnetic field strength at the production site, luminosity of the flare, and the production site bulk Lorentz factor. To check the feasibility of this model, we study the statistical relation between powerful magnetar flares and the rate of FRBs. The expected ratio is derived by convoluting the redshift-dependent magnetar density with their flare luminosity function above the energy limit determined by the FRB detection threshold. We obtain that only a small fraction, $\sim10^{-5}$, of powerful magnetar flares trigger FRBs. This ratio agrees surprisingly well with our estimates: we obtained that $10\%$ of magnetars should be in the evolutionary phase suitable for the production of FRBs, and only $10^{-4}$ of all flares are expected to be weakly magnetized, which is a necessary condition for the high-frequency maser emission. △ Less

Submitted 16 February, 2022; v1 submitted 17 June, 2021; originally announced June 2021.

Comments: 16 pages, 2 figure, ApJ in press

Evidence of 100 TeV $γ$-ray emission from HESS J1702-420: A new PeVatron candidate

Authors: H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun, P. Brun, M. Bryan , et al. (211 additional authors not shown)

Abstract: The identification of PeVatrons, hadronic particle accelerators reaching the knee of the cosmic ray spectrum (few $10^{15}$ eV), is crucial to understand the origin of cosmic rays in the Galaxy. We provide an update on the unidentified source HESS J1702-420, a promising PeVatron candidate. We present new observations of HESS J1702-420 made with the High Energy Stereoscopic System (H.E.S.S.), and p… ▽ More The identification of PeVatrons, hadronic particle accelerators reaching the knee of the cosmic ray spectrum (few $10^{15}$ eV), is crucial to understand the origin of cosmic rays in the Galaxy. We provide an update on the unidentified source HESS J1702-420, a promising PeVatron candidate. We present new observations of HESS J1702-420 made with the High Energy Stereoscopic System (H.E.S.S.), and processed using improved analysis techniques. The analysis configuration was optimized to enhance the collection area at the highest energies. We applied a three-dimensional (3D) likelihood analysis to model the source region and adjust non thermal radiative spectral models to the $γ$-ray data. We also analyzed archival data from the Fermi Large Area Telescope (LAT) to constrain the source spectrum at $γ$-ray energies >10 GeV. We report the detection of a new source component called HESS J1702-420A, that was separated from the bulk of TeV emission at a $5.4σ$ confidence level. The power law $γ$-ray spectrum of HESS J1702-420A extends with an index of $Γ=1.53\pm0.19_\text{stat}\pm0.20_\text{sys}$ and without curvature up to the energy band 64-113 TeV, in which it was detected by H.E.S.S. at a $4.0σ$ confidence level. This brings evidence for the source emission up to $100\,\text{TeV}$, which makes HESS J1702-420A a compelling candidate site for the presence of extremely high energy cosmic rays. Remarkably, in a hadronic scenario, the cut-off energy of the proton distribution powering HESS J1702-420A is found to be higher than 0.5 PeV at a 95% confidence level. HESS J1702-420A becomes therefore one of the most solid PeVatron candidates detected so far in H.E.S.S. data, altough a leptonic origin of its emission could not be ruled out either. △ Less

Submitted 14 June, 2021; v1 submitted 11 June, 2021; originally announced June 2021.

Comments: Accepted for publication in the 2. Astrophysical processes section of Astronomy & Astrophysics

Search for dark matter annihilation signals from unidentified Fermi-LAT objects with H.E.S.S

Authors: H. E. S. S. Collaboration, H. Abdallah, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun , et al. (205 additional authors not shown)

Abstract: Cosmological $N$-body simulations show that Milky Way-sized galaxies harbor a population of unmerged dark matter subhalos. These subhalos could shine in gamma-rays and be eventually detected in gamma-ray surveys as unidentified sources. We performed a thorough selection among unidentified Fermi-LAT Objects (UFOs) to identify them as possible TeV-scale dark matter subhalo candidates. We search for… ▽ More Cosmological $N$-body simulations show that Milky Way-sized galaxies harbor a population of unmerged dark matter subhalos. These subhalos could shine in gamma-rays and be eventually detected in gamma-ray surveys as unidentified sources. We performed a thorough selection among unidentified Fermi-LAT Objects (UFOs) to identify them as possible TeV-scale dark matter subhalo candidates. We search for very-high-energy (E $\gtrsim$ 100 GeV) gamma-ray emissions using H.E.S.S. observations towards four selected UFOs. Since no significant very-high-energy gamma-ray emission is detected in any dataset of the four observed UFOs nor in the combined UFO dataset, strong constraints are derived on the product of the velocity-weighted annihilation cross section $\langle σv \rangle$ by the $J$-factor for the dark matter models. The 95% C.L. observed upper limits derived from combined H.E.S.S. observations reach $\langle σv \rangle J$ values of 3.7$\times$10$^{-5}$ and 8.1$\times$10$^{-6}$ GeV$^2$cm$^{-2}$s$^{-1}$ in the $W^+W^-$ and $τ^+τ^-$ channels, respectively, for a 1 TeV dark matter mass. Focusing on thermal WIMPs, the H.E.S.S. constraints restrict the $J$-factors to lie in the range 6.1$\times$10$^{19}$ - 2.0$\times$10$^{21}$ GeV$^2$cm$^{-5}$, and the masses to lie between 0.2 and 6 TeV in the $W^+W^-$ channel. For the $τ^+τ^-$ channel, the $J$-factors lie in the range 7.0$\times$10$^{19}$ - 7.1$\times$10$^{20}$ GeV$^2$cm$^{-5}$ and the masses lie between 0.2 and 0.5 TeV. Assuming model-dependent predictions from cosmological N-body simulations on the $J$-factor distribution for Milky Way-sized galaxies, the dark matter models with masses greater than 0.3 TeV for the UFO emissions can be ruled out at high confidence level. △ Less

Submitted 15 June, 2021; v1 submitted 1 June, 2021; originally announced June 2021.

Comments: 11 pages, 7 figures, matches accepted version in The Astrophysical Journal

Journal ref: Astrophys.J., 918, 17 (2021)

Gamma-ray and Neutrino Signals from Accretion Disk Coronae of Active Galactic Nuclei

Authors: Yoshiyuki Inoue, Dmitry Khangulyan, Akihiro Doi

Abstract: To explain X-ray spectra of active galactic nuclei (AGN), non-thermal activity in AGN coronae such as pair cascade models has been extensively discussed in the past literature. Although X-ray and gamma-ray observations in the 1990s disfavored such pair cascade models, recent millimeter-wave observations of nearby Seyferts establish the existence of weak non-thermal coronal activity. Besides, the I… ▽ More To explain X-ray spectra of active galactic nuclei (AGN), non-thermal activity in AGN coronae such as pair cascade models has been extensively discussed in the past literature. Although X-ray and gamma-ray observations in the 1990s disfavored such pair cascade models, recent millimeter-wave observations of nearby Seyferts establish the existence of weak non-thermal coronal activity. Besides, the IceCube collaboration reported NGC 1068, a nearby Seyfert, as the hottest spot in their 10-yr survey. These pieces of evidence are enough to investigate the non-thermal perspective of AGN coronae in depth again. This article summarizes our current observational understandings of AGN coronae and describes how AGN coronae generate high-energy particles. We also provide ways to test the AGN corona model with radio, X-ray, MeV gamma-ray, and high-energy neutrino observations. △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: Invited review in Galaxies special issue "Searching for Astrophysical Sources of Cosmic Rays, Gamma-Rays and Neutrinos: Real-Time Multimessenger Programs and Theoretical Models". This paper is based on our previous papers arXiv:1904.00554 & 1909.02239

Report number: RIKEN-iTHEMS-Report-21

Journal ref: Galaxies 2021, 9(2), 36

Search for dark matter annihilation in the dwarf irregular galaxy WLM with H.E.S.S

Authors: H. E. S. S. Collaboration, H. Abdallah, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, F. Brun , et al. (211 additional authors not shown)

Abstract: We search for an indirect signal of dark matter through very high-energy gamma rays from the Wolf-Lundmark-Melotte (WLM) dwarf irregular galaxy. The pair annihilation of dark matter particles would produce Standard Model particles in the final state such as gamma rays, which might be detected by ground-based Cherenkov telescopes. Dwarf irregular galaxies represent promising targets as they are dar… ▽ More We search for an indirect signal of dark matter through very high-energy gamma rays from the Wolf-Lundmark-Melotte (WLM) dwarf irregular galaxy. The pair annihilation of dark matter particles would produce Standard Model particles in the final state such as gamma rays, which might be detected by ground-based Cherenkov telescopes. Dwarf irregular galaxies represent promising targets as they are dark matter dominated objects with well measured kinematics and small uncertainties on their dark matter distribution profiles. In 2018, the H.E.S.S. five-telescope array observed the dwarf irregular galaxy WLM for 18 hours. We present the first analysis based on data obtained from an imaging atmospheric Cherenkov telescope for this subclass of dwarf galaxy. As we do not observe any significant excess in the direction of WLM, we interpret the result in terms of constraints on the velocity-weighted cross section for dark matter pair annihilation as a function of the dark matter particle mass for various continuum channels as well as the prompt gamma-gamma emission. For the $τ^+τ^-$ channel the limits reach a $\langle σv \rangle$ value of about $4\times 10^{-22}$ cm3s-1 for a dark matter particle mass of 1 TeV. For the prompt gamma-gamma channel, the upper limit reaches a $\langle σv \rangle$ value of about $5 \times10^{-24}$ cm3s-1 for a mass of 370 GeV. These limits represent an improvement of up to a factor 200 with respect to previous results for the dwarf irregular galaxies for TeV dark matter search. △ Less

Submitted 10 May, 2021; originally announced May 2021.

Comments: 13 pages, 5 figures

Journal ref: Phys. Rev. D 103, 102002 (2021)

Observation of a sudden cessation of a very-high-energy gamma-ray flare in PKS 1510-089 with H.E.S.S. and MAGIC in May 2016

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Arm, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, M. de Bony de Lavergne, J. Bregeon , et al. (409 additional authors not shown)

Abstract: The flat spectrum radio quasar (FSRQ) PKS 1510-089 is known for its complex multiwavelength behavior, and is one of only a few FSRQs detected at very high energy (VHE, $E>100\,$GeV) $γ$-rays. VHE $γ$-ray observations with H.E.S.S. and MAGIC during late May and early June 2016 resulted in the detection of an unprecedented flare, which reveals for the first time VHE $γ$-ray intranight variability in… ▽ More The flat spectrum radio quasar (FSRQ) PKS 1510-089 is known for its complex multiwavelength behavior, and is one of only a few FSRQs detected at very high energy (VHE, $E>100\,$GeV) $γ$-rays. VHE $γ$-ray observations with H.E.S.S. and MAGIC during late May and early June 2016 resulted in the detection of an unprecedented flare, which reveals for the first time VHE $γ$-ray intranight variability in this source. While a common variability timescale of $1.5\,$hr is found, there is a significant deviation near the end of the flare with a timescale of $\sim 20\,$min marking the cessation of the event. The peak flux is nearly two orders of magnitude above the low-level emission. For the first time, curvature is detected in the VHE $γ$-ray spectrum of PKS 1510-089, which is fully explained through absorption by the extragalactic background light. Optical R-band observations with ATOM reveal a counterpart of the $γ$-ray flare, even though the detailed flux evolution differs from the VHE ightcurve. Interestingly, a steep flux decrease is observed at the same time as the cessation of the VHE flare. In the high energy (HE, $E>100\,$MeV) $γ$-ray band only a moderate flux increase is observed with Fermi-LAT, while the HE $γ$-ray spectrum significantly hardens up to a photon index of 1.6. A search for broad-line region (BLR) absorption features in the $γ$-ray spectrum indicates that the emission region is located outside of the BLR. Radio VLBI observations reveal a fast moving knot interacting with a standing jet feature around the time of the flare. As the standing feature is located $\sim 50\,$pc from the black hole, the emission region of the flare may have been located at a significant distance from the black hole. If this correlation is indeed true, VHE $γ$ rays have been produced far down the jet where turbulent plasma crosses a standing shock. △ Less

Submitted 18 December, 2020; originally announced December 2020.

Comments: 25 pages, 17 figures. Accepted for publication by Astronomy & Astrophysics. Corresponding authors: M. Zacharias, J. Sitarek, D. Sanchez, T. Terzic

Journal ref: A&A 648, A23 (2021)

Systematic Study of Acceleration Efficiency in Young Supernova Remnants with Nonthermal X-ray Observations

Authors: Naomi Tsuji, Yasunobu Uchiyama, Dmitry Khangulyan, Felix Aharonian

Abstract: Cutoff energy in a synchrotron radiation spectrum of a supernova remnant (SNR) contains a key parameter of ongoing particle acceleration. We systematically analyze 11 young SNRs, including all historical SNRs, to measure the cutoff energy, thus shedding light on the nature of particle acceleration at the early stage of SNR evolution. The nonthermal (synchrotron) dominated spectra in filament-like… ▽ More Cutoff energy in a synchrotron radiation spectrum of a supernova remnant (SNR) contains a key parameter of ongoing particle acceleration. We systematically analyze 11 young SNRs, including all historical SNRs, to measure the cutoff energy, thus shedding light on the nature of particle acceleration at the early stage of SNR evolution. The nonthermal (synchrotron) dominated spectra in filament-like outer rims are selectively extracted and used for spectral fitting because our model assumes that accelerated electrons are concentrated in the vicinity of the shock front due to synchrotron cooling. The cutoff energy parameter ($\varepsilon_0$) and shock speed ($v_{\rm sh}$) are related as $ \varepsilon_0 \propto v_{\rm sh}^2 η^{-1}$ with a Bohm factor of $η$. Five SNRs provide us with spatially resolved $\varepsilon_0$-$v_{\rm sh}$ plots across the remnants, indicating a variety of particle acceleration. With all SNRs considered together, the systematic tendency of $η$ clarifies a correlation between $η$ and an age of $t$ (or an expansion parameter of $m$) as $η\propto t^{-0.4}$ ($η\propto m^{4}$). This might be interpreted as the magnetic field becomes more turbulent and self-generated, as particles are accelerated at a greater rate with time. The maximum energy achieved in SNRs can be higher if we consider the newly observed time dependence on $η$. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 26 pages, 7 figure, 6 tables; Accepted for publication in ApJ

Report number: RIKEN-iTHEMS-Report-20

Journal ref: The Astrophysical Journal, 907:117 (21pp), 2021 February 1

The Crab nebula variability at short timescales with the Cherenkov Telescope Array

Authors: Enrique Mestre, Emma de Ona Wilhelmi, Dmitry Khangulyan, Roberta Zanin, Fabio Acero, Diego F. Torres

Abstract: Since 2009, several rapid and bright flares have been observed at high energies (>100 MeV) from the direction of the Crab Nebula. Several hypotheses have been put forward to explain this phenomenon, but the origin is still unclear. The detection of counterparts at higher energies with the next generation of Cherenkov telescopes will be determinant to constrain the underlying emission mechanisms. W… ▽ More Since 2009, several rapid and bright flares have been observed at high energies (>100 MeV) from the direction of the Crab Nebula. Several hypotheses have been put forward to explain this phenomenon, but the origin is still unclear. The detection of counterparts at higher energies with the next generation of Cherenkov telescopes will be determinant to constrain the underlying emission mechanisms. We aim at studying the capability of the Cherenkov Telescope Array (CTA) to explore the physics behind the flares, by performing simulations of the Crab Nebula spectral energy distribution, both in flaring and steady state, for different parameters related to the physical conditions in the nebula. In particular, we explore the data recorded by Fermi during two particular flares that occurred in 2011 and 2013. The expected GeV and TeV gamma-ray emission is derived using different radiation models. The resulting emission is convoluted with the CTA response and tested for detection, obtaining an exclusion region for the space of parameters that rule the different flare emission models. Our simulations show different scenarios that may be favourable for achieving the detection of the flares in Crab with CTA, in different regimes of energy. In particular, we find that observations with low sub-100 GeV energy threshold telescopes could provide the most model-constraining results. △ Less

Submitted 9 October, 2022; v1 submitted 17 November, 2020; originally announced November 2020.

Comments: 11 pages, 9 figures; typos corrected

Journal ref: Volume 501, Issue 1, February 2021, Pages 337-346

An extreme particle accelerator in the Galactic plane: HESS J1826$-$130

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, P. Bordas, M. Breuhaus , et al. (215 additional authors not shown)

Abstract: The unidentified very-high-energy (VHE; E $>$ 0.1 TeV) $γ$-ray source, HESS J1826$-$130, was discovered with the High Energy Stereoscopic System (HESS) in the Galactic plane. The analysis of 215 h of HESS data has revealed a steady $γ$-ray flux from HESS J1826$-$130, which appears extended with a half-width of 0.21$^{\circ}$ $\pm$ 0.02$^{\circ}_{\text{stat}}$ $\pm$ 0.05$^{\circ}_{\text{sys}}$. The… ▽ More The unidentified very-high-energy (VHE; E $>$ 0.1 TeV) $γ$-ray source, HESS J1826$-$130, was discovered with the High Energy Stereoscopic System (HESS) in the Galactic plane. The analysis of 215 h of HESS data has revealed a steady $γ$-ray flux from HESS J1826$-$130, which appears extended with a half-width of 0.21$^{\circ}$ $\pm$ 0.02$^{\circ}_{\text{stat}}$ $\pm$ 0.05$^{\circ}_{\text{sys}}$. The source spectrum is best fit with either a power-law function with a spectral index $Γ$ = 1.78 $\pm$ 0.10$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$ and an exponential cut-off at 15.2$^{+5.5}_{-3.2}$ TeV, or a broken power-law with $Γ_{1}$ = 1.96 $\pm$ 0.06$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$, $Γ_{2}$ = 3.59 $\pm$ 0.69$_{\text{stat}}$ $\pm$ 0.20$_{\text{sys}}$ for energies below and above $E_{\rm{br}}$ = 11.2 $\pm$ 2.7 TeV, respectively. The VHE flux from HESS J1826$-$130 is contaminated by the extended emission of the bright, nearby pulsar wind nebula (PWN), HESS J1825$-$137, particularly at the low end of the energy spectrum. Leptonic scenarios for the origin of HESS J1826$-$130 VHE emission related to PSR J1826$-$1256 are confronted by our spectral and morphological analysis. In a hadronic framework, taking into account the properties of dense gas regions surrounding HESS J1826$-$130, the source spectrum would imply an astrophysical object capable of accelerating the parent particle population up to $\gtrsim$200 TeV. Our results are also discussed in a multiwavelength context, accounting for both the presence of nearby supernova remnants (SNRs), molecular clouds, and counterparts detected in radio, X-rays, and TeV energies. △ Less

Submitted 25 October, 2020; originally announced October 2020.

Comments: 9 Pages, 5 Figures

Journal ref: A&A 644, A112 (2020)

Sign of hard X-ray pulsation from the gamma-ray binary system LS 5039

Authors: Hiroki Yoneda, Kazuo Makishima, Teruaki Enoto, Dmitry Khangulyan, Takahiro Matsumoto, Tadayuki Takahashi

Abstract: To understand the nature of the brightest gamma-ray binary system LS 5039, hard X-ray data of the object, taken with the Suzaku and NuSTAR observatories in 2007 and 2016, respectively, were analyzed. The two data sets jointly gave tentative evidence for a hard X-ray periodicity, with a period of $\sim 9$ s and a period increase rate by $\sim 3 \times 10^{-10}$ s s$^{-1}$. Therefore, the compact ob… ▽ More To understand the nature of the brightest gamma-ray binary system LS 5039, hard X-ray data of the object, taken with the Suzaku and NuSTAR observatories in 2007 and 2016, respectively, were analyzed. The two data sets jointly gave tentative evidence for a hard X-ray periodicity, with a period of $\sim 9$ s and a period increase rate by $\sim 3 \times 10^{-10}$ s s$^{-1}$. Therefore, the compact object in LS 5039 is inferred to be a rotating neutron star, rather than a black hole. Furthermore, several lines of arguments suggest that this object has a magnetic field of several times $\sim 10^{10}$ T, two orders of magnitude higher than those of typical neutron stars. The object is hence suggested to be a magnetar, which would be the first to be found in a binary. The results also suggest that the highly efficient particle acceleration process, known to be operating in LS 5039, emerges through interactions between dense stellar winds from the massive primary star, and ultra-strong magnetic fields of the magnetar. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Comments: 8 pages, 3 figures, 1 table, accepted for publication in PRL

Search for dark matter signals towards a selection of recently-detected DES dwarf galaxy satellites of the Milky Way with H.E.S.S

Authors: H. E. S. S. Collaboration, :, H. Abdallah, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, M. Breuhaus , et al. (206 additional authors not shown)

Abstract: Dwarf spheroidal galaxy satellites of the Milky Way are prime targets for indirect detection of dark matter with gamma rays due to their proximity, high dark matter content and absence of non-thermal emission processes. Recently, the Dark Energy Survey (DES) revealed the existence of new ultra-faint dwarf spheroidal galaxies in the southern-hemisphere sky, therefore ideally located for ground-base… ▽ More Dwarf spheroidal galaxy satellites of the Milky Way are prime targets for indirect detection of dark matter with gamma rays due to their proximity, high dark matter content and absence of non-thermal emission processes. Recently, the Dark Energy Survey (DES) revealed the existence of new ultra-faint dwarf spheroidal galaxies in the southern-hemisphere sky, therefore ideally located for ground-based observations with the imaging atmospheric Cherenkov telescope array H.E.S.S. We present a search for very-high-energy ($E\gtrsim100$ GeV) gamma-ray emission using H.E.S.S. observations carried out recently towards Reticulum II, Tucana II, Tucana III, Tucana IV and Grus II satellites. No significant very-high-energy gamma-ray excess is found from the observations on any individual object nor in the combined analysis of all the datasets. Using the most recent modeling of the dark matter distribution in the dwarf galaxy halo, we compute for the first time on DES satellites individual and combined constraints from Cherenkov telescope observations on the annihilation cross section of dark matter particles in the form of Weakly Interacting Massive Particles. The combined 95% C.L. observed upper limits reach $\langle σv \rangle \simeq 1 \times 10^{-23}$ cm$^3$s$^{-1}$ in the $W^+W^-$ channel and $4 \times 10^{-26}$ cm$^3$s$^{-1}$ in the $γγ$ channels for a dark matter mass of 1.5 TeV. The H.E.S.S. constraints well complement the results from Fermi-LAT, HAWC, MAGIC and VERITAS and are currently the most stringent in the $γγ$ channels in the multi-GeV/multi-TeV mass range. △ Less

Submitted 4 September, 2020; v1 submitted 3 August, 2020; originally announced August 2020.

Comments: 20 pages, 12 figures

Journal ref: Phys. Rev. D 102, 062001 (2020)

Resolving acceleration to very high energies along the Jet of Centaurus A

Authors: The H. E. S. S. Collaboration, :, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, H. Ashkar, M. Backes, V. Barbosa Martins, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus, F. Brun , et al. (202 additional authors not shown)

Abstract: The nearby radio galaxy Centaurus A belongs to a class of Active Galaxies that are very luminous at radio wavelengths. The majority of these galaxies show collimated relativistic outflows known as jets, that extend over hundreds of thousands of parsecs for the most powerful sources. Accretion of matter onto the central super-massive black hole is believed to fuel these jets and power their emissio… ▽ More The nearby radio galaxy Centaurus A belongs to a class of Active Galaxies that are very luminous at radio wavelengths. The majority of these galaxies show collimated relativistic outflows known as jets, that extend over hundreds of thousands of parsecs for the most powerful sources. Accretion of matter onto the central super-massive black hole is believed to fuel these jets and power their emission, with the radio emission being related to the synchrotron radiation of relativistic electrons in magnetic fields. The origin of the extended X-ray emission seen in the kiloparsec-scale jets from these sources is still a matter of debate, although Cen A's X-ray emission has been suggested to originate in electron synchrotron processes. The other possible explanation is Inverse Compton (IC) scattering with CMB soft photons. Synchrotron radiation needs ultra-relativistic electrons ($\sim50$ TeV), and given their short cooling times, requires some continuous re-acceleration mechanism to be active. IC scattering, on the other hand, does not require very energetic electrons, but requires jets that stay highly relativistic on large scales ($\geq$1 Mpc) and that remain well-aligned with the line of sight. Some recent evidence disfavours inverse Compton-CMB models, although other evidence seems to be compatible with them. In principle, the detection of extended gamma-ray emission, directly probing the presence of ultra-relativistic electrons, could distinguish between these options, but instruments have hitherto been unable to resolve the relevant structures. At GeV energies there is also an unusual spectral hardening in Cen A, whose explanation is unclear. Here we report observations of Cen A at TeV energies that resolve its large-scale jet. We interpret the data as evidence for the acceleration of ultra-relativistic electrons in the jet, and favour the synchrotron explanation for the X-rays. △ Less

Submitted 9 July, 2020; originally announced July 2020.

Comments: 32 pages, 5 figure

Journal ref: Nature volume 582, pages 356-359 (2020)

Physical Conditions and Particle Acceleration in the Kiloparsec Jet of Centaurus A

Authors: Takahiro Sudoh, Dmitry Khangulyan, Yoshiyuki Inoue

Abstract: The non-thermal emission from the kiloparsec-scale jet of Centaurus A exhibits two notable features, bright diffuse emission and many compact knots, which have been intensively studied in X-ray and radio observations. H.E.S.S. recently reported that the very-high-energy gamma-ray emission from this object is extended along the jet direction beyond a kiloparsec from the core. Here, we combine these… ▽ More The non-thermal emission from the kiloparsec-scale jet of Centaurus A exhibits two notable features, bright diffuse emission and many compact knots, which have been intensively studied in X-ray and radio observations. H.E.S.S. recently reported that the very-high-energy gamma-ray emission from this object is extended along the jet direction beyond a kiloparsec from the core. Here, we combine these observations to constrain the physical conditions of the kpc-jet and study the origin of the non-thermal emission. We show that the diffuse jet is weakly magnetized ($η_B\sim10^{-2}$) and energetically dominated by thermal particles. We also show that knots are the sites of both amplified magnetic field and particle (re-)acceleration. To keep sufficient energy in thermal particles, the magnetic and non-thermal particle energy in the knot regions are tightly constrained. The most plausible condition is an energy equipartition between them, $η_B\simη_e\sim0.1$. Such weak magnetic energy implies that particles in the knots are in the slow cooling regime. We suggest that the entire kpc-scale diffuse emission could be powered by particles that are accelerated at and escaped from knots. △ Less

Submitted 2 September, 2020; v1 submitted 22 June, 2020; originally announced June 2020.

Comments: Main text 6 pages, 2 figures, accepted for publication in ApJ Letters

Report number: RIKEN-iTHEMS-Report-20

Detection of Small Flares from the Crab Nebula with Fermi-LAT

Authors: Masanori Arakawa, Masaaki Hayashida, Dmitry Khangulyan, Yasunobu Uchiyama

Abstract: Gamma radiation from the Crab pulsar wind nebula (PWN) shows significant variability at $\sim100$ MeV energies, recently revealed with spaceborne gamma-ray telescopes. Here we report the results of a systematic search for gamma-ray flares using a 7.4-year data set acquired with the Fermi Large Area Telescope. Analyzing the off-pulse phases of the Crab pulsar, we found seven previously unreported l… ▽ More Gamma radiation from the Crab pulsar wind nebula (PWN) shows significant variability at $\sim100$ MeV energies, recently revealed with spaceborne gamma-ray telescopes. Here we report the results of a systematic search for gamma-ray flares using a 7.4-year data set acquired with the Fermi Large Area Telescope. Analyzing the off-pulse phases of the Crab pulsar, we found seven previously unreported low-intensity flares ("small flares"). The "small flares'' originate from the variable synchrotron component of the Crab PWN and show clearly different features from the steady component of the Crab PWN emission. They are characterized by larger fluxes and harder photon indices, similar to previously reported flares. These flares show day-scale time variability and imply a strong magnetic field of $B_{\rm min}\approx 1~\mathrm{mG}$ at the site of the gamma-ray production. This result seems to be inconsistent with the typical values revealed with modeling of the non-thermal emission from the nebula. The detection of the "small flares'' gives a hint of production of gamma rays above $100$ MeV in a part of the nebula with properties which are different from the main emitters, e.g., due to bulk relativistic motion. △ Less

Submitted 16 May, 2020; originally announced May 2020.

Comments: 18 pages, 9 figures, accepted for publication in ApJ

Probing the magnetic field in the GW170817 outflow using H.E.S.S. observations

Authors: H. E. S. S. Collaboration, :, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa-Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy , et al. (209 additional authors not shown)

Abstract: The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short $γ$-ray bursts and BNS mergers. It also confirmed the forging of heavy elements in the ejecta (a so-called kilonova) via the r-process nucleosynthesis. The appearance of non-thermal radio and X-ray emission, as well as the brightening, which laste… ▽ More The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short $γ$-ray bursts and BNS mergers. It also confirmed the forging of heavy elements in the ejecta (a so-called kilonova) via the r-process nucleosynthesis. The appearance of non-thermal radio and X-ray emission, as well as the brightening, which lasted more than 100 days, were somewhat unexpected. Current theoretical models attempt to explain this temporal behavior as either originating from a relativistic off-axis jet or a kilonova-like outflow. In either scenario, there is some ambiguity regarding how much energy is transported in the non-thermal electrons versus the magnetic field of the emission region. Combining the VLA (radio) and Chandra (X-ray) measurements with observations in the GeV-TeV domain can help break this ambiguity, almost independently of the assumed origin of the emission. Here we report for the first time on deep H.E.S.S. observations of GW170817 / GRB 170817A between 124 and 272 days after the BNS merger with the full H.E.S.S. array of telescopes, as well as on an updated analysis of the prompt (<5 days) observations with the upgraded H.E.S.S. phase-I telescopes. We discuss implications of the H.E.S.S. measurement for the magnetic field in the context of different source scenarios. △ Less

Submitted 18 May, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: 7 pages, 2 figures, 1 table

Journal ref: The Astrophysical Journal Letters, 894, L16, 2020

Very high energy $γ$-ray emission from two blazars of unknown redshift and upper limits on their distance

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, A. Barnacka, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus , et al. (204 additional authors not shown)

Abstract: We report on the detection of very-high-energy (VHE; $E > 100$ GeV) $γ$-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multi-wavelength observations with Fermi/LAT, XRT and UVOT on board the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum t… ▽ More We report on the detection of very-high-energy (VHE; $E > 100$ GeV) $γ$-ray emission from the BL Lac objects KUV 00311-1938 and PKS 1440-389 with the High Energy Stereoscopic System (H.E.S.S.). H.E.S.S. observations were accompanied or preceded by multi-wavelength observations with Fermi/LAT, XRT and UVOT on board the Swift satellite, and ATOM. Based on an extrapolation of the Fermi/LAT spectrum towards the VHE $γ$-ray regime, we deduce a 95% confidence level upper limit on the unknown redshift of KUV 00311-1938 of z < 0.98, and of PKS 1440-389 of z < 0.53. When combined with previous spectroscopy results the redshift of KUV 00311-1938 is constrained to $0.51 \leq z < 0.98$ and for PKS 1440-389 to $0.14 \lessapprox z < 0.53$. △ Less

Submitted 20 April, 2020; v1 submitted 7 April, 2020; originally announced April 2020.

Comments: Accepted for publication in MNRAS. 13 pages, 4 figures

Extension of the synchrotron radiation of electrons to very high energies in clumpy environments

Authors: Dmitry Khangulyan, Felix Aharonian, Carlo Romoli, Andrew Taylor

Abstract: The synchrotron cooling of relativistic electrons is one of the most effective radiation mechanisms in astrophysics. It not only accompanies the process of particle acceleration but also has feedback on the formation of the energy distribution of the parent electrons. The radiative cooling time of electrons decreases with energy as $t_{\rm syn} \propto 1/E$; correspondingly the overall radiation e… ▽ More The synchrotron cooling of relativistic electrons is one of the most effective radiation mechanisms in astrophysics. It not only accompanies the process of particle acceleration but also has feedback on the formation of the energy distribution of the parent electrons. The radiative cooling time of electrons decreases with energy as $t_{\rm syn} \propto 1/E$; correspondingly the overall radiation efficiency increases with energy. On the other hand, this effect strictly limits the maximum energy of individual photons. Even in the so-called extreme accelerators, where the acceleration proceeds at the highest possible rate, $t_{\rm acc}^{-1} = eBc/E$, allowed in an ideal magnetohydrodynamic plasma, the synchrotron radiation cannot extend well beyond the characteristic energy determined by the electron mass and the fine-structure constant: $h ν^{\rm max} \sim m_e c^2/α\sim 70 \rm\,MeV$. In this paper, we propose a model in which the formation of synchrotron radiation takes place in compact magnetic blobs located inside the particle accelerator and develop a formalism for calculations of synchrotron radiation emerging from such systems. We demonstrate that for certain combinations of parameters characterizing the accelerator and the magnetic blobs, the synchrotron radiation can extend beyond this limit by a several orders of magnitude. This scenario requires a weak magnetization of the particle accelerator, and an efficient conversion of gas internal energy into magnetic energy in sufficiently small blobs. The required size of the blobs is constrained by the magnetic mirroring effect, that can prevent particle penetration into the regions of strong magnetic field under certain conditions. △ Less

Submitted 11 May, 2021; v1 submitted 2 March, 2020; originally announced March 2020.

Comments: 14 pages, 6 figures, accepted for publication in ApJ

Journal ref: ApJ 914:76 (2021)

Fast moving pulsars as probes of interstellar medium

Authors: Maxim Barkov, Maxim Lyutikov, Dmitry Khangulyan

Abstract: Pulsars moving through ISM produce bow shocks detected in hydrogen H$α$ line emission. The morphology of the bow shock nebulae allows one to probe the properties of ISM on scales $\sim 0.01$ pc and smaller. We performed 2D RMHD modeling of the pulsar bow shock and simulated the corresponding H$α$ emission morphology. We find that even a mild spatial inhomogeneity of ISM density, $δρ/ρ\sim 1$, lead… ▽ More Pulsars moving through ISM produce bow shocks detected in hydrogen H$α$ line emission. The morphology of the bow shock nebulae allows one to probe the properties of ISM on scales $\sim 0.01$ pc and smaller. We performed 2D RMHD modeling of the pulsar bow shock and simulated the corresponding H$α$ emission morphology. We find that even a mild spatial inhomogeneity of ISM density, $δρ/ρ\sim 1$, leads to significant variations of the shape of the shock seen in H$α$ line emission. We successfully reproduce the morphology of the Guitar Nebula. We infer quasi-periodic density variations in the warm component of ISM with a characteristic length of $\sim0.1$~pc. Structures of this scale might be also responsible for the formation of the fine features seen at the forward shock of Tycho SNR in X-rays. Formation of such short periodic density structures in the warm component of ISM is puzzling, and bow-shock nebulae provide unique probes to study this phenomenon. △ Less

Submitted 3 June, 2020; v1 submitted 27 February, 2020; originally announced February 2020.

Comments: 12 pages, 7 figures, accepted MNRAS

Detection of very-high-energy γ-ray emission from the colliding wind binary η Car with H.E.S.S

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, T. Armstrong, H. Ashkar, M. Backes, V. Barbosa Martins, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus, F. Brun , et al. (210 additional authors not shown)

Abstract: Aims. Colliding wind binary systems have long been suspected to be high-energy (HE; 100 MeV < E < 100 GeV) γ-ray emitters. η Car is the most prominent member of this object class and is confirmed to emit phase-locked HE γ rays from hundreds of MeV to ~100 GeV energies. This work aims to search for and characterise the very-high-energy (VHE; E >100 GeV) γ-ray emission from η Car around the last per… ▽ More Aims. Colliding wind binary systems have long been suspected to be high-energy (HE; 100 MeV < E < 100 GeV) γ-ray emitters. η Car is the most prominent member of this object class and is confirmed to emit phase-locked HE γ rays from hundreds of MeV to ~100 GeV energies. This work aims to search for and characterise the very-high-energy (VHE; E >100 GeV) γ-ray emission from η Car around the last periastron passage in 2014 with the ground-based High Energy Stereoscopic System (H.E.S.S.). Methods. The region around η Car was observed with H.E.S.S. between orbital phase p = 0.78 - 1.10, with a closer sampling at p {\approx} 0.95 and p {\approx} 1.10 (assuming a period of 2023 days). Optimised hardware settings as well as adjustments to the data reduction, reconstruction, and signal selection were needed to suppress and take into account the strong, extended, and inhomogeneous night sky background (NSB) in the η Car field of view. Tailored run-wise Monte-Carlo simulations (RWS) were required to accurately treat the additional noise from NSB photons in the instrument response functions. Results. H.E.S.S. detected VHE γ-ray emission from the direction of η Car shortly before and after the minimum in the X-ray light-curve close to periastron. Using the point spread function provided by RWS, the reconstructed signal is point-like and the spectrum is best described by a power law. The overall flux and spectral index in VHE γ rays agree within statistical and systematic errors before and after periastron. The γ-ray spectrum extends up to at least ~400 GeV. This implies a maximum magnetic field in a leptonic scenario in the emission region of 0.5 Gauss. No indication for phase-locked flux variations is detected in the H.E.S.S. data. △ Less

Submitted 6 February, 2020; originally announced February 2020.

Comments: 9 pages, 4 figures, 3 tables, in press with A&A

Journal ref: A&A 635, A167 (2020)

H.E.S.S. and Fermi-LAT observations of PSR B1259-63/LS 2883 during its 2014 and 2017 periastron passages

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, H. Ashkar, M. Backes, V. Barbosa Martins, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus, F. Brun, P. Brun , et al. (201 additional authors not shown)

Abstract: PSR B1259-63/LS 2883 is a gamma-ray binary system consisting of a pulsar in an eccentric orbit around a bright Oe stellar-type companion star that features a dense circumstellar disc. The high- and very-high-energy (HE, VHE) gamma-ray emission from PSR B1259-63/LS 2883 around the times of its periastron passage are characterised, in particular, at the time of the HE gamma-ray flares reported to ha… ▽ More PSR B1259-63/LS 2883 is a gamma-ray binary system consisting of a pulsar in an eccentric orbit around a bright Oe stellar-type companion star that features a dense circumstellar disc. The high- and very-high-energy (HE, VHE) gamma-ray emission from PSR B1259-63/LS 2883 around the times of its periastron passage are characterised, in particular, at the time of the HE gamma-ray flares reported to have occurred in 2011, 2014, and 2017. Spectra and light curves were derived from observations conducted with the H.E.S.S.-II array in 2014 and 2017. A local double-peak profile with asymmetric peaks in the VHE light curve is measured, with a flux minimum at the time of periastron $t_p$ and two peaks coinciding with the times at which the neutron star crosses the companion's circumstellar disc ($\sim t_p \pm 16$ d). A high VHE gamma-ray flux is also observed at the times of the HE gamma-ray flares ($\sim t_p + 30$ d) and at phases before the first disc crossing ($\sim t_p - 35$ d). PSR B1259-63/LS 2883 displays periodic flux variability at VHE gamma-rays without clear signatures of super-orbital modulation in the time span covered by H.E.S.S. observations. In contrast, the photon index of the measured power-law spectra remains unchanged within uncertainties for about 200 d around periastron. Lower limits on exponential cut-off energies up to $\sim 40$ TeV are placed. At HE gamma-rays, PSR B1259-63/LS 2883 has now been detected also before and after periastron, close to the disc crossing times. Repetitive flares with distinct variability patterns are detected in this energy range. Such outbursts are not observed at VHEs, although a relatively high emission level is measured. The spectra obtained in both energy regimes displays a similar slope, although a common physical origin either in terms of a related particle population, emission mechanism, or emitter location is ruled out. △ Less

Submitted 12 December, 2019; originally announced December 2019.

Comments: 16 pages, 7 figures Accepted for publication in A&A

Journal ref: A&A 633, A102 (2020)

A very-high-energy component deep in the Gamma-ray Burst afterglow

Authors: H. Abdalla R. Adam F. Aharonian F. Ait Benkhali E. O. Anguener M. Arakawa, C. Arcaro, C. Armand, H. Ashkar, M. Backes, V. Barbosa Martins, M. Barnard Y. Becherini, D. Berge, K. Bernloehr, E. Bissaldi, R. Blackwell, M. Boettcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus, F. Brun, P. Brun, M. Bryan, M. Buechele, T. Bulik, T. Bylund, M. Capasso, S. Caroff , et al. (197 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) are brief flashes of gamma rays, considered to be the most energetic explosive phenomena in the Universe. The emission from GRBs comprises a short (typically tens of seconds) and bright prompt emission, followed by a much longer afterglow phase. During the afterglow phase, the shocked outflow -- produced by the interaction between the ejected matter and the circumburst medi… ▽ More Gamma-ray bursts (GRBs) are brief flashes of gamma rays, considered to be the most energetic explosive phenomena in the Universe. The emission from GRBs comprises a short (typically tens of seconds) and bright prompt emission, followed by a much longer afterglow phase. During the afterglow phase, the shocked outflow -- produced by the interaction between the ejected matter and the circumburst medium -- slows down, and a gradual decrease in brightness is observed. GRBs typically emit most of their energy via gamma-rays with energies in the kiloelectronvolt-to-megaelectronvolt range, but a few photons with energies of tens of gigaelectronvolts have been detected by space-based instruments. However, the origins of such high-energy (above one gigaelectronvolt) photons and the presence of very-high-energy (more than 100 gigaelectronvolts) emission have remained elussive. Here we report observations of very-high-energy emission in the bright GRB 180720B deep in the GRB afterglow -ten hours after the end of the prompt emission phase, when the X-ray flux had already decayed by four orders of magnitude. Two possible explanations exist for the observed radiation: inverse Compton emission and synchrotron emission of ultrarelativistic electrons. Our observations show that the energy fluxes in the X-ray and gamma-ray range and their photon indices remain comparable to each other throughout the afterglow. This discovery places distinct constraints on the GRB environment for both emission mechanisms, with the inverse Compton explanation alleviating the particle energy requirements for the emission observed at late times. The late timing of this detection has consequences for the future observations of GRBs at the highest energies. △ Less

Submitted 20 November, 2019; originally announced November 2019.

Comments: Preprint version of Nature paper. Contacts: E.Ruiz-Velasco, F. Aharonian, E.Bissaldi, C.Hoischen, R.D Parsons, Q.Piel, A.Taylor, D.Khangulyan

Modelling the interaction between relativistic and non-relativistic winds in binary pulsar systems: strong magnetization of the pulsar wind

Authors: Sergey V. Bogovalov, Dmitry Khangulyan, Alexander V. Koldoba, Galina V. Ustyugova, Felix Aharonian

Abstract: We present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with the circumbinary environment. We focus on studying the impact of the high wind magnetization on the geometrical structure of the shocked flow. This work is an extension of our earlier studies that focused on a purely hydrodynamic interaction and weak wind magnetization.… ▽ More We present a numerical study of the properties of the flow produced by the collision of a magnetized anisotropic pulsar wind with the circumbinary environment. We focus on studying the impact of the high wind magnetization on the geometrical structure of the shocked flow. This work is an extension of our earlier studies that focused on a purely hydrodynamic interaction and weak wind magnetization. We consider the collision in the axisymmetric approximation, that is, the pulsar rotation axis is assumed to be oriented along the line between the pulsar and the optical star. The increase of the magnetization results in the expansion of the opening cone in which the shocked pulsar wind propagates. This effect is explained in the frameworks of the conventional theory of collimation of magnetized winds. This finding has a direct implication for scenarios that involve Doppler boosting as the primary mechanism behind the GeV flares detected with the Fermi/LAT from PSR B1259-63/LS2883. The maximum enhancement of the apparent emission is determined by the ratio of $4π$ to the solid in which the shocked pulsar wind propagates. Our simulations suggest that this enhancement factor is decreased by the impact of the magnetic field. △ Less

Submitted 18 November, 2019; originally announced November 2019.

Comments: 8 pages, 5 figures, MNRAS accepted

Journal ref: 2019MNRAS.490.3601B

Detection of UHE gamma rays from the Crab Nebula: Physical Implications

Authors: Dmitry Khangulyan, Masanori Arakawa, Felix Aharonian

Abstract: The Crab Nebula is an extreme particle accelerator boosting the energy of electrons up to a few PeV ($10^{15} \ \rm eV$), close to the maximum energy allowed by theory. The physical conditions in the acceleration site and the nature of the acceleration process itself remain highly uncertain. The key information about the highest energy accelerated particles is contained in the synchrotron and inve… ▽ More The Crab Nebula is an extreme particle accelerator boosting the energy of electrons up to a few PeV ($10^{15} \ \rm eV$), close to the maximum energy allowed by theory. The physical conditions in the acceleration site and the nature of the acceleration process itself remain highly uncertain. The key information about the highest energy accelerated particles is contained in the synchrotron and inverse Compton (IC) channels of radiation at energies above 1~MeV and 100~TeV, respectively. The recent report of detection of ultra-high energy gamma-ray signal from the Crab Nebula up to 300~TeV allows one to determine the energy distribution of the highest energy electrons and to derive the magnetic field strength in the acceleration region, $B\leq120\rm\,μG$, in a parameter-free way. This estimate brings new constraints on the properties of non-thermal particle distributions and puts important constraints on the MHD models for the Crab Nebula, in particular on the feasible magnetization and anisotropy of the pulsar wind. The calculations of synchrotron and IC emission show that future observations with instruments allowing detection of the Crab Nebula above 300~TeV and above 1~MeV will clarify the conditions allowing acceleration of electrons beyond PeV energies in the Crab Nebula. In particular, one will (1) verify the hypothetical multi-component composition of the electron energy distribution, and (2) determine the magnetic field strength in the regions responsible for the acceleration of PeV electrons. △ Less

Submitted 18 November, 2019; originally announced November 2019.

Comments: 8 pages, 5 figures, MNRAS accepted

H.E.S.S. detection of very-high-energy gamma-ray emission from the quasar PKS 0736+017

Authors: H. E. S. S. Collaboration, :, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, H. Ashkar, M. Backes, V. Barbosa Martins, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus, F. Brun , et al. (203 additional authors not shown)

Abstract: Flat-spectrum radio-quasars (FSRQs) are rarely detected at very-high-energies (VHE; E>100 GeV) due to their low-frequency-peaked SEDs. At present, only 6 FSRQs are known to emit VHE photons, representing only 7% of the VHE extragalactic catalog. Following the detection of MeV-GeV gamma-ray flaring activity from the FSRQ PKS 0736+017 (z=0.189) with Fermi, the H.E.S.S. array of Cherenkov telescopes… ▽ More Flat-spectrum radio-quasars (FSRQs) are rarely detected at very-high-energies (VHE; E>100 GeV) due to their low-frequency-peaked SEDs. At present, only 6 FSRQs are known to emit VHE photons, representing only 7% of the VHE extragalactic catalog. Following the detection of MeV-GeV gamma-ray flaring activity from the FSRQ PKS 0736+017 (z=0.189) with Fermi, the H.E.S.S. array of Cherenkov telescopes triggered ToO observations on February 18, 2015, with the goal of studying the gamma-ray emission in the VHE band. H.E.S.S. ToO observations were carried out during the nights of February 18, 19, 21, and 24, 2015. Together with Fermi-LAT, the multi-wavelength coverage of the flare includes Swift observations in soft-X-rays and optical/UV, and optical monitoring (photometry and spectro-polarimetry) by the Steward Observatory, the ATOM, the KAIT and the ASAS-SN telescope. VHE emission from PKS 0736+017 was detected with H.E.S.S. during the night of February 19, 2015, only. Fermi data indicate the presence of a gamma-ray flare, peaking at the time of the H.E.S.S. detection, with a flux doubling time-scale of around six hours. The gamma-ray flare was accompanied by at least a 1 mag brightening of the non-thermal optical continuum. No simultaneous observations at longer wavelengths are available for the night of the H.E.S.S. detection. The gamma-ray observations with H.E.S.S. and Fermi are used to put constraints on the location of the gamma-ray emitting region during the flare: it is constrained to be just outside the radius of the broad-line-region with a bulk Lorentz factor $\simeq 20$, or at the level of the radius of the dusty torus with Gamma > 60. PKS 0736+017 is the seventh FSRQ known to emit VHE photons and, at z=0.189, is the nearest so far. The location of the gamma-ray emitting region during the flare can be tightly constrained thanks to opacity, variability, and collimation arguments. △ Less

Submitted 12 November, 2019; originally announced November 2019.

Comments: In press in A&A. 13 pages, 5 figures. Abstract abridged

Journal ref: A&A 633, A162 (2020)

Multi-wavelength Emission from Galactic Jets: the Case of the Microquasar SS433

Authors: Takahiro Sudoh, Yoshiyuki Inoue, Dmitry Khangulyan

Abstract: SS433 is a Galactic microquasar with powerful jets, where very-high-energy particles are produced. We study particle acceleration in the jets of SS433 in the light of the recent multi-wavelength data from radio to TeV gamma ray. We first present a general framework for the particle acceleration, cooling, and transport in relativistic jets. We then apply this to two X-ray knots in the jets of SS433… ▽ More SS433 is a Galactic microquasar with powerful jets, where very-high-energy particles are produced. We study particle acceleration in the jets of SS433 in the light of the recent multi-wavelength data from radio to TeV gamma ray. We first present a general framework for the particle acceleration, cooling, and transport in relativistic jets. We then apply this to two X-ray knots in the jets of SS433, focusing on leptonic emission. Our detailed treatment of particle transport and evolution produces substantially different predictions from previous papers. For both regions, our model can account for the multi-wavelength data except for the GeV data. This suggests that GeV emission is mostly from different regions and/or mechanisms. We find that the acceleration process should be efficient, which could be realized by diffusive shock acceleration close to the Bohm limit. Provided that protons are accelerated at the same efficiency as electrons, our results imply that SS433 is a PeVatron, i.e., can accelerate protons beyond a PeV. Future hard X-ray and MeV gamma-ray observations can critically test our models by detecting the spectral turnover or cutoff. △ Less

Submitted 17 January, 2020; v1 submitted 31 October, 2019; originally announced November 2019.

Comments: Main text 11 pages, 3 figures. To appear in ApJ

Report number: RIKEN-iTHEMS-Report-19

Resolving the Crab pulsar wind nebula at teraelectronvolt energies

Authors: H. E. S. S. Collaboration, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Arm, M. Backes, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, P. Bordas, J. Bregeon, F. Brun, P. Brun, M. Bryan, M. Büchele , et al. (199 additional authors not shown)

Abstract: The Crab nebula is one of the most studied cosmic particle accelerators, shining brightly across the entire electromagnetic spectrum up to very high-energy gamma rays. It is known from radio to gamma-ray observations that the nebula is powered by a pulsar, which converts most of its rotational energy losses into a highly relativistic outflow. This outflow powers a pulsar wind nebula (PWN), a regio… ▽ More The Crab nebula is one of the most studied cosmic particle accelerators, shining brightly across the entire electromagnetic spectrum up to very high-energy gamma rays. It is known from radio to gamma-ray observations that the nebula is powered by a pulsar, which converts most of its rotational energy losses into a highly relativistic outflow. This outflow powers a pulsar wind nebula (PWN), a region of up to 10~light-years across, filled with relativistic electrons and positrons. These particles emit synchrotron photons in the ambient magnetic field and produce very high-energy gamma rays by Compton up-scattering of ambient low-energy photons. While the synchrotron morphology of the nebula is well established, it was up to now not known in which region the very high-energy gamma rays are emitted. Here we report that the Crab nebula has an angular extension at gamma-ray energies of 52 arcseconds (assuming a Gaussian source width), significantly larger than at X-ray energies. This result closes a gap in the multi-wavelength coverage of the nebula, revealing the emission region of the highest energy gamma rays. These gamma rays are a new probe of a previously inaccessible electron and positron energy range. We find that simulations of the electromagnetic emission reproduce our new measurement, providing a non-trivial test of our understanding of particle acceleration in the Crab nebula. △ Less

Submitted 23 September, 2019; v1 submitted 20 September, 2019; originally announced September 2019.

Comments: 23 pages, 5 figures

On the Origin of High Energy Neutrinos from NGC 1068: The Role of Non-Thermal Coronal Activity

Authors: Yoshiyuki Inoue, Dmitry Khangulyan, Akihiro Doi

Abstract: NGC 1068, a nearby type-2 Seyfert galaxy, is reported as the hottest neutrino spot in the 10-year survey data of IceCube. Although there are several different possibilities for the generation of high-energy neutrinos in astrophysical sources, feasible scenarios allowing such emission in NGC 1068 have not yet been firmly defined. We show that the flux level of GeV and neutrino emission observed fro… ▽ More NGC 1068, a nearby type-2 Seyfert galaxy, is reported as the hottest neutrino spot in the 10-year survey data of IceCube. Although there are several different possibilities for the generation of high-energy neutrinos in astrophysical sources, feasible scenarios allowing such emission in NGC 1068 have not yet been firmly defined. We show that the flux level of GeV and neutrino emission observed from NGC 1068 implies that the neutrino emission can be produced only in the vicinity of the supermassive black hole in the center of the galaxy. The coronal parameters, such as magnetic field strength and corona size, making this emission possible are consistent with the spectral excess registered in the millimeter range. The suggested model and relevant physical parameters are similar to those revealed for several nearby Seyferts. Due to the internal gamma-ray attenuation, the suggested scenario cannot be verified by observations of NGC 1068 in the GeV and TeV gamma-ray energy bands. However, the optical depth is expected to become negligible for MeV gamma rays, thus future observations in this band will be able to prove our model. △ Less

Submitted 18 February, 2020; v1 submitted 5 September, 2019; originally announced September 2019.

Comments: 7 pages, 2 figures, new ALMA data are added, accepted by ApJL

Report number: RIKEN-iTHEMS-Report-19

Constraints on the emission region of 3C 279 during strong flares in 2014 and 2015 through VHE gamma-ray observations with H.E.S.S

Authors: H. E. S. S. Collaboration, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, H. Ashkar, M. Backes, V. Barbosa Martins, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus, F. Brun, P. Brun , et al. (202 additional authors not shown)

Abstract: The flat spectrum radio quasar 3C 279 is known to exhibit pronounced variability in the high-energy ($100\,$MeV$<E<100\,$GeV) $γ$-ray band, which is continuously monitored with Fermi-LAT. During two periods of high activity in April 2014 and June 2015 Target-of-Opportunity observations were undertaken with H.E.S.S. in the very-high-energy (VHE, $E>100\,$GeV) $γ$-ray domain. While the observation i… ▽ More The flat spectrum radio quasar 3C 279 is known to exhibit pronounced variability in the high-energy ($100\,$MeV$<E<100\,$GeV) $γ$-ray band, which is continuously monitored with Fermi-LAT. During two periods of high activity in April 2014 and June 2015 Target-of-Opportunity observations were undertaken with H.E.S.S. in the very-high-energy (VHE, $E>100\,$GeV) $γ$-ray domain. While the observation in 2014 provides an upper limit, the observation in 2015 results in a signal with $8.7\,σ$ significance above an energy threshold of $66\,$GeV. No VHE variability has been detected during the 2015 observations. The VHE photon spectrum is soft and described by a power-law index of $4.2\pm 0.3$. The H.E.S.S. data along with a detailed and contemporaneous multiwavelength data set provide constraints on the physical parameters of the emission region. The minimum distance of the emission region from the central black hole is estimated using two plausible geometries of the broad-line region and three potential intrinsic spectra. The emission region is confidently placed at $r\gtrsim 1.7\times10^{17}\,$cm from the black hole, i.e., beyond the assumed distance of the broad-line region. Time-dependent leptonic and lepto-hadronic one-zone models are used to describe the evolution of the 2015 flare. Neither model can fully reproduce the observations, despite testing various parameter sets. Furthermore, the H.E.S.S. data are used to derive constraints on Lorentz invariance violation given the large redshift of 3C 279. △ Less

Submitted 12 June, 2019; originally announced June 2019.

Comments: 20 pages, 12 figures, accepted by A&A

Journal ref: A&A 627, A159 (2019)

Upper Limits on Very-High-Energy Gamma-ray Emission from Core-Collapse Supernovae Observed with H.E.S.S

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, H. Ashkar, M. Backes, V. Barbosa Martins, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus, F. Brun, P. Brun , et al. (203 additional authors not shown)

Abstract: Young core-collapse supernovae with dense-wind progenitors may be able to accelerate cosmic-ray hadrons beyond the knee of the cosmic-ray spectrum, and this may result in measurable gamma-ray emission. We searched for gamma-ray emission from ten supernovae observed with the High Energy Stereoscopic System (H.E.S.S.) within a year of the supernova event. Nine supernovae were observed serendipitousl… ▽ More Young core-collapse supernovae with dense-wind progenitors may be able to accelerate cosmic-ray hadrons beyond the knee of the cosmic-ray spectrum, and this may result in measurable gamma-ray emission. We searched for gamma-ray emission from ten supernovae observed with the High Energy Stereoscopic System (H.E.S.S.) within a year of the supernova event. Nine supernovae were observed serendipitously in the H.E.S.S. data collected between December 2003 and December 2014, with exposure times ranging from 1.4 hours to 53 hours. In addition we observed SN 2016adj as a target of opportunity in February 2016 for 13 hours. No significant gamma-ray emission has been detected for any of the objects, and upper limits on the $>1$ TeV gamma-ray flux of the order of $\sim$10$^{-13}$ cm$^{-2}$s$^{-1}$ are established, corresponding to upper limits on the luminosities in the range $\sim$2 $\times$ 10$^{39}$ erg s$^{-1}$ to $\sim$1 $\times$ 10$^{42}$ erg s$^{-1}$. These values are used to place model-dependent constraints on the mass-loss rates of the progenitor stars, implying upper limits between $\sim$2 $\times 10^{-5}$ and $\sim$2 $\times 10^{-3}$M$_{\odot}$yr$^{-1}$ under reasonable assumptions on the particle acceleration parameters. △ Less

Submitted 23 April, 2019; originally announced April 2019.

Comments: Accepted for publication in A&A

H.E.S.S. observations of the flaring gravitationally lensed galaxy PKS 1830-211

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Anguener, M. Arakawa, C. Arcaro, C. Armand, M. Arrieta, M. Backes, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, K. Bernloehr, R. Blackwell, M. Boettcher, C. Boisson, J. Bolmont, S. Bonnefoy, P. Bordas, J. Bregeon, F. Brun, P. Brun , et al. (202 additional authors not shown)

Abstract: PKS 1830-211 is a known macrolensed quasar located at a redshift of z=2.5. Its high-energy gamma-ray emission has been detected with the Fermi-LAT instrument and evidence for lensing was obtained by several authors from its high-energy data. Observations of PKS 1830-211 were taken with the H.E.S.S. array of Imaging Atmospheric Cherenkov Telescopes in August 2014, following a flare alert by the F… ▽ More PKS 1830-211 is a known macrolensed quasar located at a redshift of z=2.5. Its high-energy gamma-ray emission has been detected with the Fermi-LAT instrument and evidence for lensing was obtained by several authors from its high-energy data. Observations of PKS 1830-211 were taken with the H.E.S.S. array of Imaging Atmospheric Cherenkov Telescopes in August 2014, following a flare alert by the Fermi- LAT collaboration. The H.E.S.S observations were aimed at detecting a gamma-ray flare delayed by 20-27 days from the alert flare, as expected from observations at other wavelengths. More than twelve hours of good quality data were taken with an analysis threshold of $\sim67$ GeV. The significance of a potential signal is computed as a function of the date as well as the average significance over the whole period. Data are compared to simultaneous observations by Fermi-LAT. No photon excess or significant signal is detected. An upper limit on PKS 1830-211 flux above 67 GeV is computed and compared to the extrapolation of the Fermi-LAT flare spectrum. △ Less

Submitted 10 April, 2019; originally announced April 2019.

Comments: Accepted for publication in MNRAS

On high-energy particles in accretion disk coronae of supermassive black holes: implications for MeV gamma rays and high-energy neutrinos from AGN cores

Authors: Yoshiyuki Inoue, Dmitry Khangulyan, Susumu Inoue, Akihiro Doi

Abstract: Recent observations with ALMA have revealed evidence for non-thermal synchrotron emission from the core regions of two nearby Seyfert galaxies. This suggests that the coronae of accretion disks in active galactic nuclei (AGNs) can be conducive to the acceleration of non-thermal electrons, in addition to the hot, thermal electrons responsible for their X-ray emission through thermal Comptonization.… ▽ More Recent observations with ALMA have revealed evidence for non-thermal synchrotron emission from the core regions of two nearby Seyfert galaxies. This suggests that the coronae of accretion disks in active galactic nuclei (AGNs) can be conducive to the acceleration of non-thermal electrons, in addition to the hot, thermal electrons responsible for their X-ray emission through thermal Comptonization. Here we investigate the mechanism of such particle acceleration, based on observationally inferred parameters for AGN disk coronae. One possibility to account for the observed non-thermal electrons is diffusive shock acceleration, as long as the gyrofactor $η_g$ does not exceed $\sim10^6$. These non-thermal electrons can generate gamma rays via inverse Compton scattering of disk photons, which can appear in the MeV band, while those with energies above $\sim100$ MeV would be attenuated via internal $γγ$ pair production. The integrated emission from all AGNs with thermal and non-thermal Comptonization can reproduce the observed cosmic background radiation in X-rays as well as gamma-rays up to $\sim 10$ MeV. Furthermore, if protons are accelerated in the same conditions as electrons and $η_g\sim30$, our observationally motivated model is also able to account for the diffuse neutrino flux at energies below 100-300 TeV. The next generation of MeV gamma-ray and neutrino facilities can test these expectations by searching for signals from bright, nearby Seyfert galaxies such as NGC 4151 and IC 4329A. △ Less

Submitted 3 June, 2019; v1 submitted 1 April, 2019; originally announced April 2019.

Comments: 24 pages, 10 figures, accepted for publication in ApJ

Report number: RIKEN-iTHEMS-Report-19

The 2014 TeV Gamma-ray Flare of Mrk 501 Seen with H.E.S.S.: Temporal and Spectral Constraints on Lorentz Invariance Violation

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, M. Arrieta, M. Backes, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, P. Bordas, J. Bregeon, F. Brun , et al. (210 additional authors not shown)

Abstract: The blazar Mrk 501 (z=0.034) was observed at very-high-energy (VHE, $E\gtrsim 100$~GeV) gamma-ray wavelengths during a bright flare on the night of 2014 June 23-24 (MJD 56832) with the H.E.S.S. phase-II array of Cherenkov telescopes. Data taken that night by H.E.S.S. at large zenith angle reveal an exceptional number of gamma-ray photons at multi-TeV energies, with rapid flux variability and an en… ▽ More The blazar Mrk 501 (z=0.034) was observed at very-high-energy (VHE, $E\gtrsim 100$~GeV) gamma-ray wavelengths during a bright flare on the night of 2014 June 23-24 (MJD 56832) with the H.E.S.S. phase-II array of Cherenkov telescopes. Data taken that night by H.E.S.S. at large zenith angle reveal an exceptional number of gamma-ray photons at multi-TeV energies, with rapid flux variability and an energy coverage extending significantly up to 20 TeV. This data set is used to constrain Lorentz invariance violation (LIV) using two independent channels: a temporal approach considers the possibility of an energy dependence in the arrival time of gamma rays, whereas a spectral approach considers the possibility of modifications to the interaction of VHE gamma rays with extragalactic background light (EBL) photons. The non-detection of energy-dependent time delays and the non-observation of deviations between the measured spectrum and that of a supposed power-law intrinsic spectrum with standard EBL attenuation are used independently to derive strong constraints on the energy scale of LIV ($E_{\rm{QG}}$) in the subluminal scenario for linear and quadratic perturbations in the dispersion relation of photons. For the case of linear perturbations, the 95% confidence level limits obtained are $E_{\rm{QG},1} > 3.6 \times 10^{17} \ \rm{GeV} $ using the temporal approach and $E_{\rm{QG},1} > 2.6 \times 10^{19} \ \rm{GeV}$ using the spectral approach. For the case of quadratic perturbations, the limits obtained are $E_{\rm{QG},2} > 8.5 \times 10^{10} \ \rm{GeV} $ using the temporal approach and $E_{\rm{QG},2} > 7.8 \times 10^{11} \rm{ GeV}$ using the spectral approach. △ Less

Submitted 16 January, 2019; originally announced January 2019.

Comments: Published in ApJ

Journal ref: The Astrophysical Journal, 2019, Volume 870, Number 2

Particle Transport within the Pulsar Wind Nebula HESS J1825-137

Authors: H. E. S. S. Collaboration, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, M. Arrieta, M. Backes, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, P. Bordas, J. Bregeon, F. Brun, P. Brun, M. Bryan , et al. (203 additional authors not shown)

Abstract: Aims: We present a detailed view of the pulsar wind nebula (PWN) HESS J1825-137. We aim to constrain the mechanisms dominating the particle transport within the nebula, accounting for its anomalously large size and spectral characteristics. Methods: The nebula is studied using a deep exposure from over 12 years of H.E.S.S. I operation, together with data from H.E.S.S. II improving the low energy s… ▽ More Aims: We present a detailed view of the pulsar wind nebula (PWN) HESS J1825-137. We aim to constrain the mechanisms dominating the particle transport within the nebula, accounting for its anomalously large size and spectral characteristics. Methods: The nebula is studied using a deep exposure from over 12 years of H.E.S.S. I operation, together with data from H.E.S.S. II improving the low energy sensitivity. Enhanced energy-dependent morphological and spatially-resolved spectral analyses probe the Very High Energy (VHE, E > 0.1 TeV) gamma-ray properties of the nebula. Results: The nebula emission is revealed to extend out to 1.5 degrees from the pulsar, ~1.5 times further than previously seen, making HESS J1825--137, with an intrinsic diameter of ~100 pc, potentially the largest gamma-ray PWN currently known. Characterisation of the nebula's strongly energy-dependent morphology enables the particle transport mechanisms to be constrained. A dependence of the nebula extent with energy of R $\propto$ E^αwith α= -0.29 +/- 0.04 (stat) +/- 0.05 (sys) disfavours a pure diffusion scenario for particle transport within the nebula. The total gamma-ray flux of the nebula above 1~TeV is found to be (1.12 +/- 0.03 (stat) +/- 0.25 (sys)) $\times 10^{-11}$ cm$^{-2}$ s$^{-1}$, corresponding to ~64% of the flux of the Crab Nebula. Conclusions: HESS J1825-137 is a PWN with clear energy-dependent morphology at VHE gamma-ray energies. This source is used as a laboratory to investigate particle transport within middle-aged PWNe. Deep observations of this highly spatially-extended PWN enable a spectral map of the region to be produced, providing insights into the spectral variation within the nebula. △ Less

Submitted 23 November, 2018; v1 submitted 30 October, 2018; originally announced October 2018.

Comments: 20 pages, 10 figures, 5 tables, accepted for publication by Astronomy and Astrophysics

Journal ref: A&A 621, A116 (2019)

VHE $γ$-ray discovery and multi-wavelength study of the blazar 1ES 2322-409

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, M. Arrieta, M. Backes, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, P. Bordas, J. Bregeon, F. Brun , et al. (210 additional authors not shown)

Abstract: A hotspot at a position compatible with the BL Lac object 1ES 2322-409 was serendipitously detected with H.E.S.S. during observations performed in 2004 and 2006 on the blazar PKS 2316-423. Additional data on 1ES 2322-409 were taken in 2011 and 2012, leading to a total live-time of 22.3h. Point-like very-high-energy (VHE; E>100GeV) $γ$-ray emission is detected from a source centred on the 1ES 2322-… ▽ More A hotspot at a position compatible with the BL Lac object 1ES 2322-409 was serendipitously detected with H.E.S.S. during observations performed in 2004 and 2006 on the blazar PKS 2316-423. Additional data on 1ES 2322-409 were taken in 2011 and 2012, leading to a total live-time of 22.3h. Point-like very-high-energy (VHE; E>100GeV) $γ$-ray emission is detected from a source centred on the 1ES 2322-409 position, with an excess of 116.7 events at a significance of 6.0$σ$. The average VHE $γ$-ray spectrum is well described with a power law with a photon index $Γ=3.40\pm0.66_{\text{stat}}\pm0.20_{\text{sys}}$ and an integral flux $Φ(E>200GeV) = (3.11\pm0.71_{\rm stat}\pm0.62_{\rm sys})\times10^{-12} cm^{-2} s^{-1}$, which corresponds to 1.1$\%$ of the Crab nebula flux above 200 GeV. Multi-wavelength data obtained with Fermi LAT, Swift XRT and UVOT, RXTE PCA, ATOM, and additional data from WISE, GROND and Catalina, are also used to characterise the broad-band non-thermal emission of 1ES 2322-409. The multi-wavelength behaviour indicates day-scale variability. Swift UVOT and XRT data show strong variability at longer scales. A spectral energy distribution (SED) is built from contemporaneous observations obtained around a high state identified in Swift data. A modelling of the SED is performed with a stationary homogeneous one-zone synchrotron-self-Compton (SSC) leptonic model. The redshift of the source being unknown, two plausible values were tested for the modelling. A systematic scan of the model parameters space is performed, resulting in a well-constrained combination of values providing a good description of the broad-band behaviour of 1ES 2322-409. △ Less

Submitted 10 October, 2018; originally announced October 2018.

Comments: Accepted for publication in MNRAS. 13 pages, 8 figures

Searches for gamma-ray lines and `pure WIMP' spectra from Dark Matter annihilations in dwarf galaxies with H.E.S.S

Authors: H. E. S. S. Collaboration, :, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, M. Arrieta, M. Backes, M. Barnard, Y. Becherini, J. Becker Tjus, D. Berge, S. Bernhard, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, P. Bordas, J. Bregeon, F. Brun , et al. (212 additional authors not shown)

Abstract: Dwarf spheroidal galaxies are among the most promising targets for detecting signals of Dark Matter (DM) annihilations. The H.E.S.S. experiment has observed five of these systems for a total of about 130 hours. The data are re-analyzed here, and, in the absence of any detected signals, are interpreted in terms of limits on the DM annihilation cross section. Two scenarios are considered: i) DM anni… ▽ More Dwarf spheroidal galaxies are among the most promising targets for detecting signals of Dark Matter (DM) annihilations. The H.E.S.S. experiment has observed five of these systems for a total of about 130 hours. The data are re-analyzed here, and, in the absence of any detected signals, are interpreted in terms of limits on the DM annihilation cross section. Two scenarios are considered: i) DM annihilation into mono-energetic gamma-rays and ii) DM in the form of pure WIMP multiplets that, annihilating into all electroweak bosons, produce a distinctive gamma-ray spectral shape with a high-energy peak at the DM mass and a lower-energy continuum. For case i), upper limits at 95\% confidence level of about $\langle σv \rangle \lesssim 3 \times 10^{-25}$ cm$^3$ s$^{-1}$ are obtained in the mass range of 400 GeV to 1 TeV. For case ii), the full spectral shape of the models is used and several excluded regions are identified, but the thermal masses of the candidates are not robustly ruled out. △ Less

Submitted 1 October, 2018; originally announced October 2018.

Comments: Accepted for publication in JCAP. 18 pages, 7 figures

Report number: DESY 18-157 , CERN-TH-2018-201