Showing 1–50 of 261 results for author: Kluźniak, W

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

Unveiling extended gamma-ray emission around HESS J1813-178

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

Abstract: HESS J1813$-$178 is a very-high-energy $γ$-ray source spatially coincident with the young and energetic pulsar PSR J1813$-$1749 and thought to be associated with its pulsar wind nebula (PWN). Recently, evidence for extended high-energy emission in the vicinity of the pulsar has been revealed in the Fermi Large Area Telescope (LAT) data. This motivates revisiting the HESS J1813$-$178 region, taking… ▽ More HESS J1813$-$178 is a very-high-energy $γ$-ray source spatially coincident with the young and energetic pulsar PSR J1813$-$1749 and thought to be associated with its pulsar wind nebula (PWN). Recently, evidence for extended high-energy emission in the vicinity of the pulsar has been revealed in the Fermi Large Area Telescope (LAT) data. This motivates revisiting the HESS J1813$-$178 region, taking advantage of improved analysis methods and an extended data set. Using data taken by the High Energy Stereoscopic System (H.E.S.S.) experiment and the Fermi-LAT, we aim to describe the $γ$-ray emission in the region with a consistent model, to provide insights into its origin. We performed a likelihood-based analysis on 32 hours of H.E.S.S. data and 12 years of Fermi-LAT data and fit a spectro-morphological model to the combined datasets. These results allowed us to develop a physical model for the origin of the observed $γ$-ray emission in the region. In addition to the compact very-high-energy $γ$-ray emission centered on the pulsar, we find a significant yet previously undetected component along the Galactic plane. With Fermi-LAT data, we confirm extended high-energy emission consistent with the position and elongation of the extended emission observed with H.E.S.S. These results establish a consistent description of the emission in the region from GeV energies to several tens of TeV. This study suggests that HESS J1813$-$178 is associated with a $γ$-ray PWN powered by PSR J1813$-$1749. A possible origin of the extended emission component is inverse Compton emission from electrons and positrons that have escaped the confines of the pulsar and form a halo around the PWN. △ Less

Submitted 25 March, 2024; originally announced March 2024.

Comments: 13+5 pages, 13+11 figures. Accepted for publication in A&A. Corresponding authors: T.Wach, A.Mitchell, V.Joshi, P.Chambéry

Journal ref: A&A 686, A149 (2024)

General Relativistic Magnetohydrodynamic Simulations of Accreting Tori: Resolution Study

Authors: Angelos Karakonstantakis, Debora Lančová, Włodek Kluźniak, Miljenko Čemeljić

Abstract: We present two-dimensional general relativistic radiative magnetohydrodynamical simulations of accretion disks around non-rotating stellar-mass black hole. We study the evolution of an equilibrium accreting torus in different grid resolutions to determine an adequate resolution to produce a stable turbulent disk driven by magneto-rotational instability. We evaluate the quality parameter, $Q_θ$, fr… ▽ More We present two-dimensional general relativistic radiative magnetohydrodynamical simulations of accretion disks around non-rotating stellar-mass black hole. We study the evolution of an equilibrium accreting torus in different grid resolutions to determine an adequate resolution to produce a stable turbulent disk driven by magneto-rotational instability. We evaluate the quality parameter, $Q_θ$, from the ratio of MRI wavelength to the grid zone size and examine the effect of resolution in various quantitative values such as the accretion rate, magnetisation, fluxes of physical quantities and disk scale-height. We also analyse how the resolution affects the formation of plasmoids produced in the magnetic reconnection events. △ Less

Submitted 24 March, 2024; originally announced March 2024.

Comments: Published in the Proceedings of RAGtime 23-25, https://proceedings.physics.cz/proceedings-of-ragtime-23-25/

Spectrum and extension of the inverse-Compton emission of the Crab Nebula from a combined Fermi-LAT and H.E.S.S. analysis

Authors: F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, A. Baktash, 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, J. Borowska, F. Bradascio, M. Breuhaus, R. Brose, A. Brown, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin , et al. (137 additional authors not shown)

Abstract: The Crab Nebula is a unique laboratory for studying the acceleration of electrons and positrons through their non-thermal radiation. Observations of very-high-energy $γ$ rays from the Crab Nebula have provided important constraints for modelling its broadband emission. We present the first fully self-consistent analysis of the Crab Nebula's $γ$-ray emission between 1 GeV and $\sim$100 TeV, that is… ▽ More The Crab Nebula is a unique laboratory for studying the acceleration of electrons and positrons through their non-thermal radiation. Observations of very-high-energy $γ$ rays from the Crab Nebula have provided important constraints for modelling its broadband emission. We present the first fully self-consistent analysis of the Crab Nebula's $γ$-ray emission between 1 GeV and $\sim$100 TeV, that is, over five orders of magnitude in energy. Using the open-source software package Gammapy, we combined 11.4 yr of data from the Fermi Large Area Telescope and 80 h of High Energy Stereoscopic System (H.E.S.S.) data at the event level and provide a measurement of the spatial extension of the nebula and its energy spectrum. We find evidence for a shrinking of the nebula with increasing $γ$-ray energy. Furthermore, we fitted several phenomenological models to the measured data, finding that none of them can fully describe the spatial extension and the spectral energy distribution at the same time. Especially the extension measured at TeV energies appears too large when compared to the X-ray emission. Our measurements probe the structure of the magnetic field between the pulsar wind termination shock and the dust torus, and we conclude that the magnetic field strength decreases with increasing distance from the pulsar. We complement our study with a careful assessment of systematic uncertainties. △ Less

Submitted 21 March, 2024; v1 submitted 19 March, 2024; originally announced March 2024.

Comments: 18+6 pages, 15+2 figures. Accepted for publication in A&A. Corresponding authors: M. Meyer, L. Mohrmann, T. Unbehaun. v2: after A&A language editing

Journal ref: A&A 686, A308 (2024)

Fires in the deep: The luminosity distribution of early-time gamma-ray-burst afterglows in light of the Gamow Explorer sensitivity requirements

Authors: D. A. Kann, N. E. White, G. Ghirlanda, S. R. Oates, A. Melandri, M. Jelinek, A. de Ugarte Postigo, A. J. Levan, A. Martin-Carrillo, G. S. -H. Paek, L. Izzo, M. Blazek, C. Thone, J. F. Agui Fernandez, R. Salvaterra, N. R. Tanvir, T. -C. Chang, P. O'Brien, A. Rossi, D. A. Perley, M. Im, D. B. Malesani, A. Antonelli, S. Covino, C. Choi , et al. (36 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) are ideal probes of the Universe at high redshift (z > 5), pinpointing the locations of the earliest star-forming galaxies and providing bright backlights that can be used to spectrally fingerprint the intergalactic medium and host galaxy during the period of reionization. Future missions such as Gamow Explorer are being proposed to unlock this potential by increasing the r… ▽ More Gamma-ray bursts (GRBs) are ideal probes of the Universe at high redshift (z > 5), pinpointing the locations of the earliest star-forming galaxies and providing bright backlights that can be used to spectrally fingerprint the intergalactic medium and host galaxy during the period of reionization. Future missions such as Gamow Explorer are being proposed to unlock this potential by increasing the rate of identification of high-z GRBs to rapidly trigger observations from 6-10 m ground telescopes, JWST, and the Extremely Large Telescopes. Gamow was proposed to the NASA 2021 Medium-Class Explorer (MIDEX) program as a fast-slewing satellite featuring a wide-field lobster-eye X-ray telescope (LEXT) to detect and localize GRBs, and a 30 cm narrow-field multi-channel photo-z infrared telescope (PIRT) to measure their photometric redshifts using the Lyman-alpha dropout technique. To derive the PIRT sensitivity requirement we compiled a complete sample of GRB optical-near-infrared afterglows from 2008 to 2021, adding a total of 66 new afterglows to our earlier sample, including all known high-z GRB afterglows. We performed full light-curve and spectral-energy-distribution analyses of these afterglows to derive their true luminosity at very early times. For all the light curves, where possible, we determined the brightness at the time of the initial finding chart of Gamow, at different high redshifts and in different NIR bands. We then followed the evolution of the luminosity to predict requirements for ground and space-based follow-up. We find that a PIRT sensitivity of 15 micro-Jy (21 mag AB) in a 500 s exposure simultaneously in five NIR bands within 1000s of the GRB trigger will meet the Gamow mission requirement to recover > 80% of all redshifts at z > 5. △ Less

Submitted 29 February, 2024; originally announced March 2024.

Comments: 44 pages, 10 figures, 1 table. Accepted for publication in Astronomy and Astrophysics 15 Feb 2024. Abstract abridged for arXiv

Journal ref: A&A 686, A56 (2024)

Curvature in the very-high energy gamma-ray spectrum of M87

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, F. Bradascio, R. Brose, F. Brun, B. Bruno, T. Bulik C. Burger-Scheidlin, T. Bylund, S. Casanova, R. Cecil, J. Celic, M. Cerruti , et al. (110 additional authors not shown)

Abstract: The radio galaxy M87 is a variable very-high energy (VHE) gamma-ray source, exhibiting three major flares reported in 2005, 2008, and 2010. Despite extensive studies, the origin of the VHE gamma-ray emission is yet to be understood. In this study, we investigate the VHE gamma-ray spectrum of M87 during states of high gamma-ray activity, utilizing 20.2$\,$ hours the H.E.S.S. observations. Our findi… ▽ More The radio galaxy M87 is a variable very-high energy (VHE) gamma-ray source, exhibiting three major flares reported in 2005, 2008, and 2010. Despite extensive studies, the origin of the VHE gamma-ray emission is yet to be understood. In this study, we investigate the VHE gamma-ray spectrum of M87 during states of high gamma-ray activity, utilizing 20.2$\,$ hours the H.E.S.S. observations. Our findings indicate a preference for a curved spectrum, characterized by a log-parabola model with extra-galactic background light (EBL) model above 0.3$\,$TeV at the 4$σ$ level, compared to a power-law spectrum with EBL. We investigate the degeneracy between the absorption feature and the EBL normalization and derive upper limits on EBL models mainly sensitive in the wavelength range 12.4$\,$$μ$m - 40$\,$$μ$m. △ Less

Submitted 25 April, 2024; v1 submitted 20 February, 2024; originally announced February 2024.

Comments: 10 pages, 7 figures. Accepted for publication in A&A. Corresponding authors: Victor Barbosa Martins, Rahul Cecil, Iryna Lypova, Manuel Meyer, Perri Zilberman. Supplementary material: https://zenodo.org/records/10781524

Journal ref: A&A, 685, A96 (2024)

Acceleration and transport of relativistic electrons in the jets of the microquasar SS 433

Authors: F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, 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, J. Borowska, M. Bouyahiaou, M. Breuhau, R. Brose, A. M. Brown, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin, S. Caroff , et al. (140 additional authors not shown)

Abstract: SS 433 is a microquasar, a stellar binary system with collimated relativistic jets. We observed SS 433 in gamma rays using the High Energy Stereoscopic System (H.E.S.S.), finding an energy-dependent shift in the apparent position of the gamma-ray emission of the parsec-scale jets. These observations trace the energetic electron population and indicate the gamma rays are produced by inverse-Compton… ▽ More SS 433 is a microquasar, a stellar binary system with collimated relativistic jets. We observed SS 433 in gamma rays using the High Energy Stereoscopic System (H.E.S.S.), finding an energy-dependent shift in the apparent position of the gamma-ray emission of the parsec-scale jets. These observations trace the energetic electron population and indicate the gamma rays are produced by inverse-Compton scattering. Modelling of the energy-dependent gamma-ray morphology constrains the location of particle acceleration and requires an abrupt deceleration of the jet flow. We infer the presence of shocks on either side of the binary system at distances of 25 to 30 parsecs and conclude that self-collimation of the precessing jets forms the shocks, which then efficiently accelerate electrons. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: Submitted 20th Apr. 2023, published 25th January 2024 (accepted version)

Journal ref: Science383,402-406(2024)

Truncated, Tilted Discs as a Possible Source of Quasi-Periodic Oscillations

Authors: A. Deepika Bollimpalli, P. Chris Fragile, W. Janosz Dewberry, Włodek Kluźniak

Abstract: Many accreting black holes and neutron stars exhibit rapid variability in their X-ray light curves, termed quasi-periodic oscillations (QPOs). The most commonly observed type is the low-frequency ($\lesssim 10$ Hz), type-C QPO, while only a handful of sources exhibit high-frequency QPOs ($\gtrsim 60$ Hz). The leading model for the type-C QPO is Lense-Thirring precession of a hot, geometrically thi… ▽ More Many accreting black holes and neutron stars exhibit rapid variability in their X-ray light curves, termed quasi-periodic oscillations (QPOs). The most commonly observed type is the low-frequency ($\lesssim 10$ Hz), type-C QPO, while only a handful of sources exhibit high-frequency QPOs ($\gtrsim 60$ Hz). The leading model for the type-C QPO is Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with the black hole's spin axis. However, existing versions of this model have not taken into account the effects of a surrounding, geometrically thin disc on the precessing, inner, geometrically thick flow. In Bollimpalli et. al 2023, using a set of GRMHD simulations of tilted, truncated accretion discs, we confirmed that the outer thin disc slows down the precession rate of the precessing torus, which has direct observational implications for type-C QPOs. In this paper, we provide a detailed analysis of those simulations and compare them with an aligned truncated disc simulation. We find that the misalignment of the disc excites additional variability in the inner hot flow, which is absent in the comparable aligned-disc simulations. This suggests that the misalignment may be a crucial requirement for producing QPOs. We attribute this variability to global vertical oscillations of the inner torus at epicyclic frequencies corresponding to the transition radius. This explanation is consistent with current observations of higher frequency QPOs in black hole X-ray binary systems. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Comments: 17 pages, 15 figures, 1 table; accepted for publication in MNRAS

Energy dissipation in astrophysical simulations: results of the Orszag-Tang test problem

Authors: Fatemeh Kayanikhoo, Miljenko Cemeljic, Maciek Wielgus, Wlodek Kluzniak

Abstract: The magnetic field through the magnetic reconnection process affects the dynamics and structure of astrophysical systems. Numerical simulations are the tools to study the evolution of these systems. However, the resolution, dimensions, resistivity, and turbulence of the system are some important parameters to take into account in the simulations. In this paper, we investigate the evolution of magn… ▽ More The magnetic field through the magnetic reconnection process affects the dynamics and structure of astrophysical systems. Numerical simulations are the tools to study the evolution of these systems. However, the resolution, dimensions, resistivity, and turbulence of the system are some important parameters to take into account in the simulations. In this paper, we investigate the evolution of magnetic energy in astrophysical simulations by performing a standard test problem for MHD codes, Orszag-Tang. We estimate the numerical dissipation in the simulations using state-of-the-art numerical simulation code in astrophysics, PLUTO. The estimated numerical resistivity in 2D simulations corresponds to the Lundquist number $\approx 10^{4}$ in the resolution of $512\times512$ grid cells. It is also shown that the plasmoid unstable reconnection layer can be resolved with sufficient resolutions. Our analysis demonstrates that in non-relativistic magnetohydrodynamics simulations, magnetic and kinetic energies undergo conversion into internal energy, resulting in plasma heating. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 8 pages, 3 figures, Accepted in Proceedings of RAGtime 23-25

MSC Class: 85-08 ACM Class: J.2

Chasing Gravitational Waves with the Cherenkov Telescope Array

Authors: Jarred Gershon Green, Alessandro Carosi, Lara Nava, Barbara Patricelli, Fabian Schüssler, Monica Seglar-Arroyo, Cta Consortium, :, Kazuki Abe, Shotaro Abe, Atreya Acharyya, Remi Adam, Arnau Aguasca-Cabot, Ivan Agudo, Jorge Alfaro, Nuria Alvarez-Crespo, Rafael Alves Batista, Jean-Philippe Amans, Elena Amato, Filippo Ambrosino, Ekrem Oguzhan Angüner, Lucio Angelo Antonelli, Carla Aramo, Cornelia Arcaro, Luisa Arrabito , et al. (545 additional authors not shown)

Abstract: The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very… ▽ More The detection of gravitational waves from a binary neutron star merger by Advanced LIGO and Advanced Virgo (GW170817), along with the discovery of the electromagnetic counterparts of this gravitational wave event, ushered in a new era of multimessenger astronomy, providing the first direct evidence that BNS mergers are progenitors of short gamma-ray bursts (GRBs). Such events may also produce very-high-energy (VHE, > 100GeV) photons which have yet to be detected in coincidence with a gravitational wave signal. The Cherenkov Telescope Array (CTA) is a next-generation VHE observatory which aims to be indispensable in this search, with an unparalleled sensitivity and ability to slew anywhere on the sky within a few tens of seconds. New observing modes and follow-up strategies are being developed for CTA to rapidly cover localization areas of gravitational wave events that are typically larger than the CTA field of view. This work will evaluate and provide estimations on the expected number of of gravitational wave events that will be observable with CTA, considering both on- and off-axis emission. In addition, we will present and discuss the prospects of potential follow-up strategies with CTA. △ Less

Submitted 5 February, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: Presented at the 38th International Cosmic Ray Conference (ICRC 2023), 2023 (arXiv:2309.08219)

Report number: CTA-ICRC/2023/30

Discovery of a Radiation Component from the Vela Pulsar Reaching 20 Teraelectronvolts

Authors: The 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, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, J. Borowska, F. Bradascio, M. Breuhaus, R. Brose, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin , et al. (157 additional authors not shown)

Abstract: Gamma-ray observations have established energetic isolated pulsars as outstanding particle accelerators and antimatter factories in the Galaxy. There is, however, no consensus regarding the acceleration mechanisms and the radiative processes at play, nor the locations where these take place. The spectra of all observed gamma-ray pulsars to date show strong cutoffs or a break above energies of a fe… ▽ More Gamma-ray observations have established energetic isolated pulsars as outstanding particle accelerators and antimatter factories in the Galaxy. There is, however, no consensus regarding the acceleration mechanisms and the radiative processes at play, nor the locations where these take place. The spectra of all observed gamma-ray pulsars to date show strong cutoffs or a break above energies of a few gigaelectronvolt (GeV). Using the H.E.S.S. array of Cherenkov telescopes, we discovered a novel radiation component emerging beyond this generic GeV cutoff in the Vela pulsar's broadband spectrum. The extension of gamma-ray pulsation energies up to at least 20 teraelectronvolts (TeV) shows that Vela pulsar can accelerate particles to Lorentz factors higher than $4\times10^7$. This is an order of magnitude larger than in the case of the Crab pulsar, the only other pulsar detected in the TeV energy range. Our results challenge the state-of-the-art models for high-energy emission of pulsars while providing a new probe, i.e. the energetic multi-TeV component, for constraining the acceleration and emission processes in their extreme energy limit. △ Less

Submitted 9 October, 2023; originally announced October 2023.

Comments: 38 pages, 6 figures. This preprint has not undergone peer review or any post-submission improvements or corrections. The Version of Record of this article is published in Nature Astronomy, Nat Astron (2023), and is available online at https://doi.org/10.1038/s41550-023-02052-3

Prospects for $γ$-ray observations of the Perseus galaxy cluster with the Cherenkov Telescope Array

Authors: The Cherenkov Telescope Array Consortium, :, K. Abe, S. Abe, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Alfaro, N. Alvarez-Crespo, R. Alves Batista, J. -P. Amans, E. Amato, E. O. Angüner, L. A. Antonelli, C. Aramo, M. Araya, C. Arcaro, L. Arrabito, K. Asano, Y. Ascasíbar, J. Aschersleben , et al. (542 additional authors not shown)

Abstract: Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster med… ▽ More Galaxy clusters are expected to be dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay at gamma-ray energies and are predicted to be sources of large-scale gamma-ray emission due to hadronic interactions in the intracluster medium. We estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuse gamma-ray emission from the Perseus galaxy cluster. We perform a detailed spatial and spectral modelling of the expected signal for the DM and the CRp components. For each, we compute the expected CTA sensitivity. The observing strategy of Perseus is also discussed. In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratio within the radius $R_{500}$ down to about $X_{500}<3\times 10^{-3}$, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index $α_{\rm CRp}=2.3$. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure $α_{\rm CRp}$ down to about $Δα_{\rm CRp}\simeq 0.1$ and the CRp spatial distribution with 10% precision. Regarding DM, CTA should improve the current ground-based gamma-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to $\sim 5$, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models with $τ_χ>10^{27}$s for DM masses above 1 TeV. These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. △ Less

Submitted 7 September, 2023; originally announced September 2023.

Comments: 93 pages (including author list, appendix and references), 143 figures. Submitted to JCAP

Energy distribution and substructure formation in astrophysical MHD simulations

Authors: Fatemeh Kayanikhoo, Miljenko Cemeljic, Maciek Wielgus, Wlodek Kluzniak

Abstract: During substructure formation in magnetized astrophysical plasma, dissipation of magnetic energy facilitated by magnetic reconnection affects the system dynamics by heating and accelerating the ejected plasmoids. Numerical simulations are a crucial tool for investigating such systems. In astrophysical simulations, the energy dissipation, reconnection rate and substructure formation critically depe… ▽ More During substructure formation in magnetized astrophysical plasma, dissipation of magnetic energy facilitated by magnetic reconnection affects the system dynamics by heating and accelerating the ejected plasmoids. Numerical simulations are a crucial tool for investigating such systems. In astrophysical simulations, the energy dissipation, reconnection rate and substructure formation critically depend on the onset of reconnection of numerical or physical origin. In this paper, we hope to assess the reliability of the state-of-the-art numerical codes, PLUTO and KORAL by quantifying and discussing the impact of dimensionality, resolution, and code accuracy on magnetic energy dissipation, reconnection rate, and substructure formation. We quantitatively compare results obtained with relativistic and non-relativistic, resistive and non-resistive, as well as two- and three-dimensional setups performing the Orszag-Tang test problem. We find the sufficient resolution in each model, for which numerical error is negligible and the resolution does not significantly affect the magnetic energy dissipation and reconnection rate. The non-relativistic simulations show that at sufficient resolution, magnetic and kinetic energies convert to internal energy and heat up the plasma. The results show that in the relativistic system, energy components undergo mutual conversion during the simulation time, which leads to a substantial increase in magnetic energy at 20\% and 90\% of the total simulation time of $10$ light-crossing times -- the magnetic field is amplified by a factor of five due to relativistic shocks. We also show that the reconnection rate in all our simulations is higher than $0.1$, indicating plasmoid-mediated regime. It is shown that in KORAL simulations magnetic energy is slightly larger and more substructures are captured than in PLUTO simulations. △ Less

Submitted 15 December, 2023; v1 submitted 30 August, 2023; originally announced August 2023.

Comments: 17 pages, 20 Figures,

MSC Class: 85-08 ACM Class: J.2

Magnetically threaded accretion disks in resistive magnetohydrodynamic simulations and asymptotic expansion

Authors: Miljenko Čemeljić, Włodek Kluźniak, Varadarajan Parthasarathy

Abstract: Aims. A realistic model of magnetic linkage between a central object and its accretion disk is a prerequisite for understanding the spin history of stars and stellar remnants. To this end, we aim to provide an analytic model in agreement with magnetohydrodynamic (MHD) simulations. Methods. For the first time, we wrote a full set of stationary asymptotic expansion equations of a thin magnetic acc… ▽ More Aims. A realistic model of magnetic linkage between a central object and its accretion disk is a prerequisite for understanding the spin history of stars and stellar remnants. To this end, we aim to provide an analytic model in agreement with magnetohydrodynamic (MHD) simulations. Methods. For the first time, we wrote a full set of stationary asymptotic expansion equations of a thin magnetic accretion disk, including the induction and energy equations. We also performed a resistive MHD simulation of an accretion disk around a star endowed with a magnetic dipole, using the publicly available code PLUTO. We compared the analytical results with the numerical solutions, and discussed the results in the context of previous solutions of the induction equation describing the star-disk magnetospheric interaction. Results. We found that the magnetic field threading the disk is suppressed by orders of magnitude inside thin disks, so the presence of the stellar magnetic field does not strongly affect the velocity field, nor the density profile inside the disk. Density and velocity fields found in the MHD simulations match the radial and vertical profiles of the analytic solution. Qualitatively, the MHD simulations result in an internal magnetic field similar to the solutions previously obtained by solving the induction equation in the disk alone. However, the magnetic field configuration is quantitatively affected by magnetic field inflation outside the disk; this is reflected in the net torque. The torque on the star is an order of magnitude larger in the magnetic than in the non-magnetic case. Spin-up of the star occurs on a timescale comparable to the accretion timescale in the MHD case, and is an order of magnitude slower in the absence of a stellar magnetic field. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 17 pages, accepted in Astronomy and Astrophysics. arXiv admin note: text overlap with arXiv:1907.12592

Journal ref: A&A 678, A57 (2023)

Multiwavelength Observations of the Blazar PKS 0735+178 in Spatial and Temporal Coincidence with an Astrophysical Neutrino Candidate IceCube-211208A

Authors: A. Acharyya, C. B. Adams, A. Archer, P. Bangale, J. T. Bartkoske, P. Batista, W. Benbow, A. Brill, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. Errando, A. Falcone, Q. Feng, G. M. Foote, L. Fortson, A. Furniss, G. Gallagher, W. Hanlon, D. Hanna, O. Hervet, C. E. Hinrichs, J. Hoang, J. Holder, T. B. Humensky , et al. (185 additional authors not shown)

Abstract: We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ra… ▽ More We report on multiwavelength target-of-opportunity observations of the blazar PKS 0735+178, located 2.2$^\circ$ away from the best-fit position of the IceCube neutrino event IceCube-211208A detected on December 8, 2021. The source was in a high-flux state in the optical, ultraviolet, X-ray, and GeV gamma-ray bands around the time of the neutrino event, exhibiting daily variability in the soft X-ray flux. The X-ray data from Swift-XRT and NuSTAR characterize the transition between the low-energy and high-energy components of the broadband spectral energy distribution (SED), and the gamma-ray data from Fermi -LAT, VERITAS, and H.E.S.S. require a spectral cut-off near 100 GeV. Both X-ray and gamma-ray measurements provide strong constraints on the leptonic and hadronic models. We analytically explore a synchrotron self-Compton model, an external Compton model, and a lepto-hadronic model. Models that are entirely based on internal photon fields face serious difficulties in matching the observed SED. The existence of an external photon field in the source would instead explain the observed gamma-ray spectral cut-off in both leptonic and lepto-hadronic models and allow a proton jet power that marginally agrees with the Eddington limit in the lepto-hadronic model. We show a numerical lepto-hadronic model with external target photons that reproduces the observed SED and is reasonably consistent with the neutrino event despite requiring a high jet power. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: 21 pages, 3 figures, accepted by ApJ

Constraints on the intergalactic magnetic field using Fermi-LAT and H.E.S.S. blazar observations

Authors: H. E. S. S., Fermi-LAT Collaborations, :, F. Aharonian, J. Aschersleben, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, B. Bi, M. Bouyahiaoui, M. Breuhaus, R. Brose, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin, T. Bylund, S. Caroff, S. Casanova, J. Celic, M. Cerruti, T. Chand, S. Chandra , et al. (113 additional authors not shown)

Abstract: Magnetic fields in galaxies and galaxy clusters are believed to be the result of the amplification of intergalactic seed fields during the formation of large-scale structures in the universe. However, the origin, strength, and morphology of this intergalactic magnetic field (IGMF) remain unknown. Lower limits on (or indirect detection of) the IGMF can be obtained from observations of high-energy g… ▽ More Magnetic fields in galaxies and galaxy clusters are believed to be the result of the amplification of intergalactic seed fields during the formation of large-scale structures in the universe. However, the origin, strength, and morphology of this intergalactic magnetic field (IGMF) remain unknown. Lower limits on (or indirect detection of) the IGMF can be obtained from observations of high-energy gamma rays from distant blazars. Gamma rays interact with the extragalactic background light to produce electron-positron pairs, which can subsequently initiate electromagnetic cascades. The $γ$-ray signature of the cascade depends on the IGMF since it deflects the pairs. Here we report on a new search for this cascade emission using a combined data set from the Fermi Large Area Telescope and the High Energy Stereoscopic System. Using state-of-the-art Monte Carlo predictions for the cascade signal, our results place a lower limit on the IGMF of $B > 7.1\times10^{-16}$ G for a coherence length of 1 Mpc even when blazar duty cycles as short as 10 yr are assumed. This improves on previous lower limits by a factor of 2. For longer duty cycles of $10^4$ ($10^7$) yr, IGMF strengths below $1.8\times10^{-14}$ G ($3.9\times10^{-14}$ G) are excluded, which rules out specific models for IGMF generation in the early universe. △ Less

Submitted 8 June, 2023; originally announced June 2023.

Comments: 20 pages, 7 figures, 4 tables. Accepted for publication in ApJ Letters. Auxiliary data is provided in electronic format at https://zenodo.org/record/8014311

Journal ref: ApJ Letters 2023, Volume 950, Number 2 950, L16

Constraining the cosmic-ray pressure in the inner Virgo Cluster using H.E.S.S. observations of M 87

Authors: H. E. S. S. Collaboration, :, F. Aharonian, F. Ait Benkhali, C. Arcaro, J. Aschersleben, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlöhr, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, J. Borowska, F. Bradascio, M. Breuhaus, R. Brose, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin, T. Bylund , et al. (139 additional authors not shown)

Abstract: The origin of the gamma-ray emission from M87 is currently a matter of debate. This work aims to localize the VHE (100 GeV-100 TeV) gamma-ray emission from M87 and probe a potential extended hadronic emission component in the inner Virgo Cluster. The search for a steady and extended gamma-ray signal around M87 can constrain the cosmic-ray energy density and the pressure exerted by the cosmic rays… ▽ More The origin of the gamma-ray emission from M87 is currently a matter of debate. This work aims to localize the VHE (100 GeV-100 TeV) gamma-ray emission from M87 and probe a potential extended hadronic emission component in the inner Virgo Cluster. The search for a steady and extended gamma-ray signal around M87 can constrain the cosmic-ray energy density and the pressure exerted by the cosmic rays onto the intra-cluster medium, and allow us to investigate the role of the cosmic rays in the active galactic nucleus feedback as a heating mechanism in the Virgo Cluster. H.E.S.S. telescopes are sensitive to VHE gamma rays and have been utilized to observe M87 since 2004. We utilized a Bayesian block analysis to identify M87 emission states with H.E.S.S. observations from 2004 until 2021, dividing them into low, intermediate, and high states. Because of the causality argument, an extended ($\gtrsim$kpc) signal is allowed only in steady emission states. Hence, we fitted the morphology of the 120h low state data and found no significant gamma-ray extension. Therefore, we derived for the low state an upper limit of 58"(corresponding to $\approx$4.6kpc) in the extension of a single-component morphological model described by a rotationally symmetric 2D Gaussian model at 99.7% confidence level. Our results exclude the radio lobes ($\approx$30 kpc) as the principal component of the VHE gamma-ray emission from the low state of M87. The gamma-ray emission is compatible with a single emission region at the radio core of M87. These results, with the help of two multiple-component models, constrain the maximum cosmic-ray to thermal pressure ratio $X_{CR,max.}$$\lesssim$$0.32$ and the total energy in cosmic-ray protons (CRp) to $U_{CR}$$\lesssim$5$\times10^{58}$ erg in the inner 20kpc of the Virgo Cluster for an assumed CRp power-law distribution in momentum with spectral index $α_{p}$=2.1. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: 15 pages, 7 figures. Accepted for publication in A&A. Corresponding authors: Victor Barbosa Martins, Stefan Ohm, Cornelia Arcaro, Natalia Żywucka, Mathieu de Naurois

Journal ref: A&A 675, A138 (2023)

Astrophysical cloaking of a naked singularity

Authors: Ronaldo S. S. Vieira, Włodek Kluźniak

Abstract: A massive naked singularity would be cloaked by accreted matter, and thus may appear to a distant observer as an opaque \mbox{(quasi-)}spherical surface of a fluid, not unlike that of a star or planet. We present here analytical solutions for levitating atmospheres around a wide class of spherically symmetric naked singularities. Such an atmosphere can be constructed in every spacetime which posse… ▽ More A massive naked singularity would be cloaked by accreted matter, and thus may appear to a distant observer as an opaque \mbox{(quasi-)}spherical surface of a fluid, not unlike that of a star or planet. We present here analytical solutions for levitating atmospheres around a wide class of spherically symmetric naked singularities. Such an atmosphere can be constructed in every spacetime which possesses a zero-gravity radius and which is a solution of a (modified-)gravity theory possessing the usual conservation laws for matter. Its density peaks at the zero-gravity radius and the atmospheric fluid is supported against infall onto the singularity by gravity alone. In an astrophysical context, an opaque atmosphere would be formed in a very short time by accretion of ambient matter onto the singularity -- in a millisecond for an X-ray binary, in a thousand seconds for a singularity traversing interstellar space, and a thousand years for a singularity that is the central engine of an AGN. △ Less

Submitted 28 June, 2023; v1 submitted 12 April, 2023; originally announced April 2023.

Comments: 9 pages, 4 figures. Accepted in MNRAS

Journal ref: MNRAS 523, 4615--4623 (2023)

Detection of extended gamma-ray emission around the Geminga pulsar with H.E.S.S

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, B. Bi, M. Böttcher, C. Boisson, J. Bolmont, J. Borowska, M. Bouyahiaoui, F. Bradascio, R. Brose, F. Brun, B. Bruno, T. Bulik, C. Burger Scheidlin, F. Cangemi , et al. (143 additional authors not shown)

Abstract: Geminga is an enigmatic radio-quiet gamma-ray pulsar located at a mere 250 pc distance from Earth. Extended very-high-energy gamma-ray emission around the pulsar was discovered by Milagro and later confirmed by HAWC, which are both water Cherenkov detector-based experiments. However, evidence for the Geminga pulsar wind nebula in gamma rays has long evaded detection by imaging atmospheric Cherenko… ▽ More Geminga is an enigmatic radio-quiet gamma-ray pulsar located at a mere 250 pc distance from Earth. Extended very-high-energy gamma-ray emission around the pulsar was discovered by Milagro and later confirmed by HAWC, which are both water Cherenkov detector-based experiments. However, evidence for the Geminga pulsar wind nebula in gamma rays has long evaded detection by imaging atmospheric Cherenkov telescopes (IACTs) despite targeted observations. The detection of gamma-ray emission on angular scales > 2 deg poses a considerable challenge for the background estimation in IACT data analysis. With recent developments in understanding the complementary background estimation techniques of water Cherenkov and atmospheric Cherenkov instruments, the H.E.S.S. IACT array can now confirm the detection of highly extended gamma-ray emission around the Geminga pulsar with a radius of at least 3 deg in the energy range 0.5-40 TeV. We find no indications for statistically significant asymmetries or energy-dependent morphology. A flux normalisation of $(2.8\pm0.7)\times10^{-12}$ cm$^{-2}$s$^{-1}$TeV$^{-1}$ at 1 TeV is obtained within a 1 deg radius region around the pulsar. To investigate the particle transport within the halo of energetic leptons around the pulsar, we fitted an electron diffusion model to the data. The normalisation of the diffusion coefficient obtained of $D_0 = 7.6^{+1.5}_{-1.2} \times 10^{27}$ cm$^2$s$^{-1}$, at an electron energy of 100 TeV, is compatible with values previously reported for the pulsar halo around Geminga, which is considerably below the Galactic average. △ Less

Submitted 5 April, 2023; originally announced April 2023.

Comments: 16 pages, 15 figures, 7 tables. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 673, A148 (2023)

Sensitivity of the Cherenkov Telescope Array to spectral signatures of hadronic PeVatrons with application to Galactic Supernova Remnants

Authors: The Cherenkov Telescope Array Consortium, F. Acero, A. Acharyya, R. Adam, A. Aguasca-Cabot, I. Agudo, A. Aguirre-Santaella, J. Alfaro, R. Aloisio, N. Álvarez Crespo, R. Alves Batista, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, C. Aramo, C. Arcaro, T. Armstrong, K. Asano, Y. Ascasibar, J. Aschersleben, M. Backes, A. Baktash, C. Balazs, M. Balbo , et al. (334 additional authors not shown)

Abstract: The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The pote… ▽ More The local Cosmic Ray (CR) energy spectrum exhibits a spectral softening at energies around 3~PeV. Sources which are capable of accelerating hadrons to such energies are called hadronic PeVatrons. However, hadronic PeVatrons have not yet been firmly identified within the Galaxy. Several source classes, including Galactic Supernova Remnants (SNRs), have been proposed as PeVatron candidates. The potential to search for hadronic PeVatrons with the Cherenkov Telescope Array (CTA) is assessed. The focus is on the usage of very high energy $γ$-ray spectral signatures for the identification of PeVatrons. Assuming that SNRs can accelerate CRs up to knee energies, the number of Galactic SNRs which can be identified as PeVatrons with CTA is estimated within a model for the evolution of SNRs. Additionally, the potential of a follow-up observation strategy under moonlight conditions for PeVatron searches is investigated. Statistical methods for the identification of PeVatrons are introduced, and realistic Monte--Carlo simulations of the response of the CTA observatory to the emission spectra from hadronic PeVatrons are performed. Based on simulations of a simplified model for the evolution for SNRs, the detection of a $γ$-ray signal from in average 9 Galactic PeVatron SNRs is expected to result from the scan of the Galactic plane with CTA after 10 hours of exposure. CTA is also shown to have excellent potential to confirm these sources as PeVatrons in deep observations with $\mathcal{O}(100)$ hours of exposure per source. △ Less

Submitted 27 March, 2023; originally announced March 2023.

Comments: 34 pages, 16 figures, Accepted for publication in Astroparticle Physics

Search for the evaporation of primordial black holes with H.E.S.S

Authors: H. E. S. S. collaboration, :, F. Aharonian, F. Ait Benkhali, J. Aschersleben, M. Boettcher, M. Backes, V. Barbosa Martins, R. Batzo, Y. Becherini, D. Berge, B. Bi, C. Boisson, J. Bolmont, M. de Bony de Lavergne, J. Borowska, F. Bradascio, R. Brose, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin, S. Caro, S. Casanova, J. Celic , et al. (124 additional authors not shown)

Abstract: Primordial Black Holes (PBHs) are hypothetical black holes predicted to have been formed from density fluctuations in the early Universe. PBHs with an initial mass around $10^{14}-10^{15}$g are expected to end their evaporation at present times in a burst of particles and very-high-energy (VHE) gamma rays. Those gamma rays may be detectable by the High Energy Stereoscopic System (H.E.S.S.), an arr… ▽ More Primordial Black Holes (PBHs) are hypothetical black holes predicted to have been formed from density fluctuations in the early Universe. PBHs with an initial mass around $10^{14}-10^{15}$g are expected to end their evaporation at present times in a burst of particles and very-high-energy (VHE) gamma rays. Those gamma rays may be detectable by the High Energy Stereoscopic System (H.E.S.S.), an array of imaging atmospheric Cherenkov telescopes. This paper reports on the search for evaporation bursts of VHE gamma rays with H.E.S.S., ranging from 10 to 120 seconds, as expected from the final stage of PBH evaporation and using a total of 4816 hours of observations. The most constraining upper limit on the burst rate of local PBHs is $2000$ pc$^{-3}$ yr$^{-1}$ for a burst interval of 120 seconds, at the 95\% confidence level. The implication of these measurements for PBH dark matter are also discussed. △ Less

Submitted 22 March, 2023; originally announced March 2023.

Comments: Accepted for publication in JCAP, corresponding authors: F.Brun, J-F. Glicenstein, V. Marandon, T. Tavernier

H.E.S.S. follow-up observations of GRB221009A

Authors: H. E. S. S. Collaboration, :, F. Aharonian, F. Ait Benkhali, J. Aschersleben, H. Ashkar, M. Backes, A. Baktash, 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, J. Borowska, M. Bouyahiaoui, F. Bradascio, M. Breuhaus, R. Brose, F. Brun, B. Bruno , et al. (138 additional authors not shown)

Abstract: GRB221009A is the brightest gamma-ray burst ever detected. To probe the very-high-energy (VHE, $>$\!100 GeV) emission, the High Energy Stereoscopic System (H.E.S.S.) began observations 53 hours after the triggering event, when the brightness of the moonlight no longer precluded observations. We derive differential and integral upper limits using H.E.S.S. data from the third, fourth, and ninth nigh… ▽ More GRB221009A is the brightest gamma-ray burst ever detected. To probe the very-high-energy (VHE, $>$\!100 GeV) emission, the High Energy Stereoscopic System (H.E.S.S.) began observations 53 hours after the triggering event, when the brightness of the moonlight no longer precluded observations. We derive differential and integral upper limits using H.E.S.S. data from the third, fourth, and ninth nights after the initial GRB detection, after applying atmospheric corrections. The combined observations yield an integral energy flux upper limit of $Φ_\mathrm{UL}^{95\%} = 9.7 \times 10^{-12}~\mathrm{erg\,cm^{-2}\,s^{-1}}$ above $E_\mathrm{thr} = 650$ GeV. The constraints derived from the H.E.S.S. observations complement the available multiwavelength data. The radio to X-ray data are consistent with synchrotron emission from a single electron population, with the peak in the SED occurring above the X-ray band. Compared to the VHE-bright GRB190829A, the upper limits for GRB221009A imply a smaller gamma-ray to X-ray flux ratio in the afterglow. Even in the absence of a detection, the H.E.S.S. upper limits thus contribute to the multiwavelength picture of GRB221009A, effectively ruling out an IC dominated scenario. △ Less

Submitted 18 March, 2023; originally announced March 2023.

Comments: 10 pages, 4 figures. Accepted for publication in APJL. Corresponding authors: J. Damascene Mbarubucyeye, H. Ashkar, S. J. Zhu, B. Reville, F. Schüssler

HESS J1809$-$193: a halo of escaped electrons around a pulsar wind nebula?

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, M. Böttcher, C. Boisson, J. Bolmont, J. Borowska, M. Bouyahiaoui, F. Bradascio, M. Breuhaus, R. Brose, F. Brun, B. Bruno, T. Bulik, C. Burger-Scheidlin, T. Bylund, S. Caroff , et al. (130 additional authors not shown)

Abstract: Context. HESS J1809$-$193 is an unassociated very-high-energy $γ$-ray source located on the Galactic plane. While it has been connected to the nebula of the energetic pulsar PSR J1809$-$1917, supernova remnants and molecular clouds present in the vicinity also constitute possible associations. Recently, the detection of $γ$-ray emission up to energies of $\sim$100 TeV with the HAWC observatory has… ▽ More Context. HESS J1809$-$193 is an unassociated very-high-energy $γ$-ray source located on the Galactic plane. While it has been connected to the nebula of the energetic pulsar PSR J1809$-$1917, supernova remnants and molecular clouds present in the vicinity also constitute possible associations. Recently, the detection of $γ$-ray emission up to energies of $\sim$100 TeV with the HAWC observatory has led to renewed interest in HESS J1809$-$193. Aims. We aim to understand the origin of the $γ$-ray emission of HESS J1809$-$193. Methods. We analysed 93.2 h of data taken on HESS J1809$-$193 above 0.27 TeV with the High Energy Stereoscopic System (H.E.S.S.), using a multi-component, three-dimensional likelihood analysis. In addition, we provide a new analysis of 12.5 yr of Fermi-LAT data above 1 GeV within the region of HESS J1809$-$193. The obtained results are interpreted in a time-dependent modelling framework. Results. For the first time, we were able to resolve the emission detected with H.E.S.S. into two components: an extended component that exhibits a spectral cut-off at $\sim$13 TeV, and a compact component that is located close to PSR J1809$-$1917 and shows no clear spectral cut-off. The Fermi-LAT analysis also revealed extended $γ$-ray emission, on scales similar to that of the extended H.E.S.S. component. Conclusions. Our modelling indicates that based on its spectrum and spatial extent, the extended H.E.S.S. component is likely caused by inverse Compton emission from old electrons that form a halo around the pulsar wind nebula. The compact component could be connected to either the pulsar wind nebula or the supernova remnant and molecular clouds. Due to its comparatively steep spectrum, modelling the Fermi-LAT emission together with the H.E.S.S. components is not straightforward. (abridged) △ Less

Submitted 27 February, 2023; originally announced February 2023.

Comments: 14 pages, 10 figures. Accepted for publication in A&A. Corresponding authors: Vikas Joshi, Lars Mohrmann

Journal ref: A&A 672, A103 (2023)

Effect of geometrically thin discs on precessing, thick flows: Relevance to type-C QPOs

Authors: D. A. Bollimpalli, P. C. Fragile, W. Kluźniak

Abstract: Type-C quasi-periodic oscillations (QPOs) are the low-frequency QPOs most commonly observed during the hard spectral state of X-ray binary systems. The leading model for these QPOs is the Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with respect to the black hole spin axis. However, none of the work done to date has accounted for the effects of a surrou… ▽ More Type-C quasi-periodic oscillations (QPOs) are the low-frequency QPOs most commonly observed during the hard spectral state of X-ray binary systems. The leading model for these QPOs is the Lense-Thirring precession of a hot, geometrically thick accretion flow that is misaligned with respect to the black hole spin axis. However, none of the work done to date has accounted for the effects of a surrounding, geometrically thin disc on this precession, as would be the case in the truncated disc picture of the hard state. To address this, we perform a set of GRMHD simulations of truncated discs misaligned with the spin axes of their central black holes. Our results confirm that the inner-hot flow still undergoes precession, though at a rate that is only 5 percent of what is predicted for an isolated, precessing torus. We find that the exchange of angular momentum between the outer, thin and the inner, thick disc causes this slow-down in the precession rate and discuss its relevance to type-C QPOs. △ Less

Submitted 6 October, 2022; originally announced October 2022.

Comments: Submitted to MNRAS Letters

Gamma-ray observations of MAXI J1820+070 during the 2018 outburst

Authors: H. Abe, S. Abe, V. A. Acciari, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari , et al. (418 additional authors not shown)

Abstract: MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS expe… ▽ More MAXI J1820+070 is a low-mass X-ray binary with a black hole as a compact object. This binary underwent an exceptionally bright X-ray outburst from March to October 2018, showing evidence of a non-thermal particle population through its radio emission during this whole period. The combined results of 59.5 hours of observations of the MAXI J1820+070 outburst with the H.E.S.S., MAGIC and VERITAS experiments at energies above 200 GeV are presented, together with Fermi-LAT data between 0.1 and 500 GeV, and multiwavelength observations from radio to X-rays. Gamma-ray emission is not detected from MAXI J1820+070, but the obtained upper limits and the multiwavelength data allow us to put meaningful constraints on the source properties under reasonable assumptions regarding the non-thermal particle population and the jet synchrotron spectrum. In particular, it is possible to show that, if a high-energy gamma-ray emitting region is present during the hard state of the source, its predicted flux should be at most a factor of 20 below the obtained Fermi-LAT upper limits, and closer to them for magnetic fields significantly below equipartition. During the state transitions, under the plausible assumption that electrons are accelerated up to ~ 500 GeV, the multiwavelength data and the gamma-ray upper limits lead consistently to the conclusion that a potential high-energy and very-high-energy gamma-ray emitting region should be located at a distance from the black hole ranging between 10^11 and 10^13 cm. Similar outbursts from low-mass X-ray binaries might be detectable in the near future with upcoming instruments such as CTA. △ Less

Submitted 6 October, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: 15 pages, 5 figures. Accepted for publication in MNRAS

Accretion disc backflow in resistive MHD simulations

Authors: R. Mishra, M. Čemeljić, W. Kluźniak

Abstract: We investigate accretion onto a central star, with the size, rotation rate, and magnetic dipole of a young stellar object, to study the flow pattern (velocity and density) of the fluid within and outside of the disc. We perform resistive MHD simulations of thin $α$-discs, varying the parameters such as the stellar rotation rate and magnetic field, and (anomalous) coefficients of viscosity and resi… ▽ More We investigate accretion onto a central star, with the size, rotation rate, and magnetic dipole of a young stellar object, to study the flow pattern (velocity and density) of the fluid within and outside of the disc. We perform resistive MHD simulations of thin $α$-discs, varying the parameters such as the stellar rotation rate and magnetic field, and (anomalous) coefficients of viscosity and resistivity in the disc. To provide a benchmark for the results and to compare with known analytic results, we also perform purely hydrodynamic simulations (HD) for the same problem. Although obtained for a different situation with differing inner boundary condition, the disc structure in the HD simulations closely follows the analytic solution of Kluźniak and Kita (2000) -- in particular a region of "midplane" backflow exists in the right range of radii, depending on the viscosity parameter. In the MHD solutions, whenever the magnetic Prandtl number does not exceed a certain critical value, the midplane backflow exists throughout the accretion disc, extending all the way down to the inner transition zone where the disc transitions to a magnetic funnel flow. For values of the magnetic Prandtl number close to the critical value the backflow and the inner disc undergo a quasiperiodic radial oscillation, otherwise the backflow is steady, as is the disc solution. From our results, supplemented by our reading of the literature, we conclude that midplane backflow is a real feature of at least some accretion discs, whether HD $α$-discs or MHD discs, including ones driven by MRI turbulence. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Comments: 11 pages, 11 figures, 1 table Submitted to Monthly Notices of the Royal Astronomical Society (MNRAS)

A deep spectromorphological study of the $γ$-ray emission surrounding the young massive stellar cluster Westerlund 1

Authors: F. Aharonian, H. Ashkar, M. Backes, V. Barbosa Martins, Y. Becherini, D. Berge, B. Bi, M. Böttcher, M. de Bony de Lavergne, F. Bradascio, R. Brose, F. Brun, T. Bulik, C. Burger-Scheidlin, F. Cangemi, S. Caroff, S. Casanova, M. Cerruti, T. Chand, S. Chandra, A. Chen, O. Chibueze, P. Cristofari, J. Damascene Mbarubucyeye, A. Djannati-Ataï , et al. (134 additional authors not shown)

Abstract: Young massive stellar clusters are extreme environments and potentially provide the means for efficient particle acceleration. Indeed, they are increasingly considered as being responsible for a significant fraction of cosmic rays (CRs) accelerated within the Milky Way. Westerlund 1, the most massive known young stellar cluster in our Galaxy is a prime candidate for studying this hypothesis. While… ▽ More Young massive stellar clusters are extreme environments and potentially provide the means for efficient particle acceleration. Indeed, they are increasingly considered as being responsible for a significant fraction of cosmic rays (CRs) accelerated within the Milky Way. Westerlund 1, the most massive known young stellar cluster in our Galaxy is a prime candidate for studying this hypothesis. While the very-high-energy $γ$-ray source HESS J1646-458 has been detected in the vicinity of Westerlund 1 in the past, its association could not be firmly identified. We aim to identify the physical processes responsible for the $γ$-ray emission around Westerlund 1 and thus to better understand the role of massive stellar clusters in the acceleration of Galactic CRs. Using 164 hours of data recorded with the High Energy Stereoscopic System (H.E.S.S.), we carried out a deep spectromorphological study of the $γ$-ray emission of HESS J1646-458. We furthermore employed H I and CO observations of the region to infer the presence of gas that could serve as target material for interactions of accelerated CRs. We detected large-scale ($\sim 2^\circ$ diameter) $γ$-ray emission with a complex morphology, exhibiting a shell-like structure and showing no significant variation with $γ$-ray energy. The combined energy spectrum of the emission extends to several tens of TeV, and is uniform across the entire source region. We did not find a clear correlation of the $γ$-ray emission with gas clouds as identified through H I and CO observations. We conclude that, of the known objects within the region, only Westerlund 1 can explain the bulk of the $γ$-ray emission. Several CR acceleration sites and mechanisms are conceivable, and discussed in detail. (abridged) △ Less

Submitted 10 November, 2022; v1 submitted 22 July, 2022; originally announced July 2022.

Comments: 15 pages, 9 figures. Corresponding authors: L. Mohrmann, S. Ohm, R. Rauth, A. Specovius. v2: corrected affiliation of M. Vecchi (still incorrect in journal version)

Journal ref: A&A 666, A124 (2022)

Search for dark matter annihilation signals in the H.E.S.S. Inner Galaxy Survey

Authors: H. E. S. S. Collaboration, H. Abdalla, F. Aharonian, F. Ait Benkhali, E. O. Anguner, C. Armand, H. Ashkar, M. Backes, V. Baghmanyan, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernlohr, B. Bi, M. Bottcher, J. Bolmont, M. de Bony de Lavergne, R. Brose, F. Brun, F. Cangemi, S. Caroff, M. Cerruti, T. Chand, A. Chen , et al. (116 additional authors not shown)

Abstract: The central region of the Milky Way is one of the foremost locations to look for dark matter (DM) signatures. We report the first results on a search for DM particle annihilation signals using new observations from an unprecedented gamma-ray survey of the Galactic Center (GC) region, ${\it i.e.}$, the Inner Galaxy Survey, at very high energies ($\gtrsim$ 100 GeV) performed with the H.E.S.S. array… ▽ More The central region of the Milky Way is one of the foremost locations to look for dark matter (DM) signatures. We report the first results on a search for DM particle annihilation signals using new observations from an unprecedented gamma-ray survey of the Galactic Center (GC) region, ${\it i.e.}$, the Inner Galaxy Survey, at very high energies ($\gtrsim$ 100 GeV) performed with the H.E.S.S. array of five ground-based Cherenkov telescopes. No significant gamma-ray excess is found in the search region of the 2014-2020 dataset and a profile likelihood ratio analysis is carried out to set exclusion limits on the annihilation cross section $\langle σv\rangle$. Assuming Einasto and Navarro-Frenk-White (NFW) DM density profiles at the GC, these constraints are the strongest obtained so far in the TeV DM mass range. For the Einasto profile, the constraints reach $\langle σv\rangle$ values of $\rm 3.7\times10^{-26} cm^3s^{-1}$ for 1.5 TeV DM mass in the $W^+W^-$ annihilation channel, and $\rm 1.2 \times 10^{-26} cm^3s^{-1}$ for 0.7 TeV DM mass in the $τ^+τ^-$ annihilation channel. With the H.E.S.S. Inner Galaxy Survey, ground-based $γ$-ray observations thus probe $\langle σv\rangle$ values expected from thermal-relic annihilating TeV DM particles. △ Less

Submitted 21 July, 2022; originally announced July 2022.

Comments: Accepted in Phys. Rev. Lett., includes Supplemental materials. 5+10 pages, 2+8 figures, 4 tables

Journal ref: Phys. Rev. Lett. 129, 111101 (2022)

Observational properties of puffy disks: radiative GRMHD spectra of mildly sub-Eddington accretion

Authors: Maciek Wielgus, Debora Lancova, Odele Straub, Wlodek Kluzniak, Ramesh Narayan, David Abarca, Agata Rozanska, Frederic Vincent, Gabriel Torok, Marek Abramowicz

Abstract: Numerical general relativistic radiative magnetohydrodynamic simulations of accretion disks around a stellar mass black hole with a luminosity above 0.5 of the Eddington value reveal their stratified, elevated vertical structure. We refer to these thermally stable numerical solutions as puffy disks. Above a dense and geometrically thin core of dimensionless thickness $h/r \sim 0.1$, crudely resemb… ▽ More Numerical general relativistic radiative magnetohydrodynamic simulations of accretion disks around a stellar mass black hole with a luminosity above 0.5 of the Eddington value reveal their stratified, elevated vertical structure. We refer to these thermally stable numerical solutions as puffy disks. Above a dense and geometrically thin core of dimensionless thickness $h/r \sim 0.1$, crudely resembling a classic thin accretion disk, a puffed-up, geometrically thick layer of lower density and $h/r \sim 1.0$ is formed. We discuss the observational properties of puffy disks, in particular the geometrical obscuration of the inner disk by the elevated puffy region at higher observing inclinations, and collimation of the radiation along the accretion disk spin axis, which may explain the apparent super-Eddington luminosity of some X-ray objects. We also present synthetic spectra of puffy disks, and show that they are qualitatively similar to those of a Comptonized thin disk. We demonstrate that the existing xspec spectral fitting models provide good fits to synthetic observations of puffy disks, but cannot correctly recover the input black hole spin. The puffy region remains optically thick to scattering; in its spectral properties the puffy disk roughly resembles that of a warm corona sandwiching the disk core. We suggest that puffy disks may correspond to X-ray binary systems of luminosities above 0.3 of the Eddington luminosity in the intermediate spectral states. △ Less

Submitted 14 July, 2022; v1 submitted 17 February, 2022; originally announced February 2022.

Journal ref: MNRAS, Volume 514, Issue 1, pp.780-789 (2022)

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)

Evidence for gamma-ray emission from the remnant of Kepler's supernova based on deep H.E.S.S. observations

Authors: H. E. S. S. Collaboration, F. Aharonian, F. Ait Benkhali, E. O. Anguner, H. Ashkar, M. Backes, V. Barbosa Martins, R. Batzofin, Y. Becherini, D. Berge, K. Bernloehr, M. Boettcher, C. Boisson, J. Bolmont, M. de Bony de Lavergne, M. Breuhaus, R. Brose, F. Brun, T. Bulik, T. Bylund, F. Cangemi, S. Caroff, S. Casanova, M. Cerruti, T. Chand , et al. (136 additional authors not shown)

Abstract: Observations with imaging atmospheric Cherenkov telescopes (IACTs) have enhanced our knowledge of nearby supernova (SN) remnants with ages younger than 500 years by establishing Cassiopeia A and the remnant of Tycho's SN as very-high-energy (VHE) gamma-ray sources. The remnant of Kepler's SN, which is the product of the most recent naked-eye supernova in our Galaxy, is comparable in age to the oth… ▽ More Observations with imaging atmospheric Cherenkov telescopes (IACTs) have enhanced our knowledge of nearby supernova (SN) remnants with ages younger than 500 years by establishing Cassiopeia A and the remnant of Tycho's SN as very-high-energy (VHE) gamma-ray sources. The remnant of Kepler's SN, which is the product of the most recent naked-eye supernova in our Galaxy, is comparable in age to the other two, but is significantly more distant. If the gamma-ray luminosities of the remnants of Tycho's and Kepler's SNe are similar, then the latter is expected to be one of the faintest gamma-ray sources within reach of the current generation IACT arrays. Here we report evidence at a statistical level of 4.6 sigma for a VHE signal from the remnant of Kepler's SN based on deep observations by the High Energy Stereoscopic System (H.E.S.S.) with an exposure of 152 hours. The measured integral flux above an energy of 226 GeV is ~0.3% of the flux of the Crab Nebula. The spectral energy distribution (SED) reveals a gamma-ray emitting component connecting the VHE emission observed with H.E.S.S. to the emission observed at GeV energies with Fermi-LAT. The overall SED is similar to that of the remnant of Tycho's SN, possibly indicating the same non-thermal emission processes acting in both these young remnants of thermonuclear SNe. △ Less

Submitted 23 March, 2024; v1 submitted 15 January, 2022; originally announced January 2022.

Comments: Published in Astronomy and Astrophysics. A typo in the normalization constant in Sect. 4 was corrected

Journal ref: A&A 662, A65 (2022)

A MeerKAT, e-MERLIN, H.E.S.S. and Swift search for persistent and transient emission associated with three localised FRBs

Authors: James O. Chibueze, M. Caleb, L. Spitler, H. Ashkar, F. Schussler, B. W. Stappers, C. Venter, I. Heywood, A. M. S. Richards, D. R. A. Williams, M. Kramer, R. Beswick, M. C. Bezuidenhout, R. P. Breton, L. N. Driessen, F. Jankowski, E. F. Keane, M. Malenta, M. Mickaliger, V. Morello, H. Qiu, K. Rajwade, S. Sanidas, M. Surnis, T. W. Scragg , et al. (134 additional authors not shown)

Abstract: We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, an… ▽ More We report on a search for persistent radio emission from the one-off Fast Radio Burst (FRB) 20190714A, as well as from two repeating FRBs, 20190711A and 20171019A, using the MeerKAT radio telescope. For FRB 20171019A we also conducted simultaneous observations with the High Energy Stereoscopic System (H.E.S.S.) in very high energy gamma rays and searched for signals in the ultraviolet, optical, and X-ray bands. For this FRB, we obtain a UV flux upper limit of 1.39x10^-16 erg/cm^-2/s/Amstrong, X-ray limit of ~ 6.6x10^-14 erg/cm^-2/s and a limit on the very-high-energy gamma-ray flux (Phi) (E > 120 GeV) < 1.7 x 10^-12 erg/cm^-2/s. We obtain a radio upper limit of ~15 microJy/beam for persistent emission at the locations of both FRBs 20190711A and 20171019A, but detect diffuse radio emission with a peak brightness of ~53 microJy/beam associated with FRB 20190714A at z = 0.2365. This represents the first detection of the radio continuum emission potentially associated with the host (galaxy) of FRB 20190714A, and is only the third known FRB to have such an association. Given the possible association of a faint persistent source, FRB 20190714A may potentially be a repeating FRB whose age lies between that of FRB 20121102A and FRB 20180916A. A parallel search for repeat bursts from these FRBs revealed no new detections down to a fluence of 0.08 Jy ms for a 1 ms duration burst. △ Less

Submitted 31 December, 2021; originally announced January 2022.

Comments: 15 pages, 9 figures

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

Levitating atmospheres around naked singularities

Authors: Ronaldo S. S. Vieira, Włodek Kluźniak

Abstract: For a wide class of spherically symmetric naked singularities there is a sphere within which gravity is effectively repulsive. In such spacetimes accreting matter cannot reach the singularity and will instead form a levitating atmosphere, which is kept suspended by gravity alone. The density of the atmosphere has a maximum at a definite radius. In its qualitative properties the atmosphere is analo… ▽ More For a wide class of spherically symmetric naked singularities there is a sphere within which gravity is effectively repulsive. In such spacetimes accreting matter cannot reach the singularity and will instead form a levitating atmosphere, which is kept suspended by gravity alone. The density of the atmosphere has a maximum at a definite radius. In its qualitative properties the atmosphere is analogous to the recently discussed atmospheres that are supported by radiation pressure above luminous neutron stars, however for the levitating atmospheres around a naked singularity no radiation need be present. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: 5 pages, no figures

Journal ref: Proceedings of RAGtime 20-22: Workshops on black holes and neutron stars, 15-19 Oct., 16-20 Sept., 19-23 Oct. 2018/2019/2020, pp. 351-355, ed. by Z. Stuchlik, G. Torok and V. Karas, Opava, Czech Republic

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)

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.

Beamed emission from a neutron-star ULX in a GRRMHD simulation

Authors: David Abarca, Kyle Parfrey, Włodek Kluźniak

Abstract: We perform a global 2.5D general-relativistic radiation magnetohydrodynamic simulation of super-critical accretion onto a neutron star with a $2\times 10^{10}$ G dipolar magnetic field, as a model of a neutron-star-powered ultraluminous X-ray source (ULX). We compute a lower limit on the total luminosity of $\sim 2.5\,L_\mathrm{Edd}$, and find the radiation to be highly beamed by the accretion dis… ▽ More We perform a global 2.5D general-relativistic radiation magnetohydrodynamic simulation of super-critical accretion onto a neutron star with a $2\times 10^{10}$ G dipolar magnetic field, as a model of a neutron-star-powered ultraluminous X-ray source (ULX). We compute a lower limit on the total luminosity of $\sim 2.5\,L_\mathrm{Edd}$, and find the radiation to be highly beamed by the accretion disk outflows. The apparent isotropic luminosity, which is a function of the viewing angle, reaches a maximum above $100\,L_\mathrm{Edd}$, consistent with the luminosities observed in ULXs. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: 9 pages, 4 figures, Submitted to ApJ Letters

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)

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)

Backflow in simulated MHD accretion disks

Authors: Ruchi Mishra, Miljenko Čemeljić, Wlodek Kluźniak

Abstract: We perform resistive MHD simulations of accretion disk with alpha-viscosity, accreting onto a rotating star endowed with a magnetic dipole. We find backflow in the presence of strong magnetic field and large resistivity, and probe for the dependence on Prandtl number. We find that in the magnetic case the distance from the star at which backflow begins, the stagnation radius, is different than in… ▽ More We perform resistive MHD simulations of accretion disk with alpha-viscosity, accreting onto a rotating star endowed with a magnetic dipole. We find backflow in the presence of strong magnetic field and large resistivity, and probe for the dependence on Prandtl number. We find that in the magnetic case the distance from the star at which backflow begins, the stagnation radius, is different than in the hydrodynamic case, and the backflow shows non-stationary behavior. We compare the results with hydrodynamics simulations. △ Less

Submitted 24 December, 2020; originally announced December 2020.

Comments: 6 pages, 2 figures; to appear in RAGtime22 conference Proceedings

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)

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)

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)

Asymmetric Jet launching

Authors: Aleksandra Kotek, Miljenko Čemeljić, Włodek Kluźniak, Deepika A. Bollimpalli

Abstract: In resistive and viscous magnetohydrodynamical (MHD) simulations we obtain axial jets launched from the innermost magnetosphere of a star-disk system. We found that in a part of the parameter space continuous asymmetric jets, which are propagating in opposite directions, are launched. We compare the speed of propagation and rotation of obtained jets with recent observational results. In resistive and viscous magnetohydrodynamical (MHD) simulations we obtain axial jets launched from the innermost magnetosphere of a star-disk system. We found that in a part of the parameter space continuous asymmetric jets, which are propagating in opposite directions, are launched. We compare the speed of propagation and rotation of obtained jets with recent observational results. △ Less

Submitted 7 June, 2020; originally announced June 2020.

Comments: Proceedings of Polish Astronomical Society from Annual PAS meeting in Olsztyn

Backflow in Accretion Disk

Authors: Ruchi Mishra, Miljenko Čemeljić, Wlodek Kluźniak

Abstract: Analytical solution for a thin hydrodynamic accretion disk shows that for some values of the viscosity parameter part of the accretion flow in the disk is not towards the star, but in the opposite direction. We study thin disks by performing hydrodynamic simulations and compare the numerical with analytical results. We confirm that for viscosity coefficient smaller than a critical value, there is… ▽ More Analytical solution for a thin hydrodynamic accretion disk shows that for some values of the viscosity parameter part of the accretion flow in the disk is not towards the star, but in the opposite direction. We study thin disks by performing hydrodynamic simulations and compare the numerical with analytical results. We confirm that for viscosity coefficient smaller than a critical value, there is a backflow in simulations near the disk mid-plane. The distance from the star to the starting point of backflow is increasing with viscosity, as predicted by the analytical solution. When the viscosity coefficient is larger than critical, there is no backflow in the disk. △ Less

Submitted 6 June, 2020; v1 submitted 2 June, 2020; originally announced June 2020.

Comments: "submitted as a conference contribution (39. PTA meeting; Olsztyn, September 9-12, 2019) to the Polish Astronomical Society Proceedings"

Transport of dust grain particles in the accretion disk

Authors: Robert Jaros, Miljenko Čemeljić, Włodek Kluźniak, Dejan Vinković, Cezary Turski

Abstract: Entrainment of dust particles in the flow inside and outside of the proto-planetary disk has implications for the disk evolution and composition of planets. Using quasi-stationary solutions in our star-disk simulations as a background, we add dust particles of different radii in post-processing of the results, using our Python tool DUSTER. The distribution and motion of particles in the disk is fo… ▽ More Entrainment of dust particles in the flow inside and outside of the proto-planetary disk has implications for the disk evolution and composition of planets. Using quasi-stationary solutions in our star-disk simulations as a background, we add dust particles of different radii in post-processing of the results, using our Python tool DUSTER. The distribution and motion of particles in the disk is followed in the cases with and without the backflow in the disk. We also compare the results with and without the radiation pressure included in the computation. △ Less

Submitted 1 June, 2020; originally announced June 2020.

Looking for the underlying cause of black hole X-ray variability in GRMHD simulations

Authors: D. A. Bollimpalli, R. Mahmoud, C. Done, P. C. Fragile, W. Kluźniak, R. Narayan, C. J. White

Abstract: Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrin… ▽ More Long-term observations have shown that black hole X-ray binaries exhibit strong, aperiodic variability on time-scales of a few milliseconds to seconds. The observed light curves display various characteristic features like a log-normal distribution of flux and a linear rms-flux relation, which indicate that the underlying variability process is stochastic in nature. It is also thought to be intrinsic to accretion. This variability has been modelled as inward propagating fluctuations of mass accretion rate, although the physical process driving the fluctuations remains puzzling. In this work, we analyse five exceptionally long duration general relativistic magnetohydrodynamic (GRMHD) simulations of optically thin, geometrically thick, black hole accretion flows to look for hints of propagating fluctuations in the simulation data. We find that the accretion profiles from these simulations do show evidence for inward propagating fluctuations below the viscous frequency by featuring strong radial coherence and positive time lags when comparing smaller to larger radii, although these time lags are generally shorter than the viscous time-scale and frequency independent. Our simulations also support the notion that the fluctuations in $\dot{M}$ build up in a multiplicative manner, as the simulations exhibit linear rms-mass flux relations, as well as log-normal distributions of their mass fluxes. When combining the mass fluxes from the simulations with an assumed emissivity profile, we additionally find broad agreement with observed power spectra and time lags, including a recovery of the frequency dependency of the time lags. △ Less

Submitted 27 July, 2020; v1 submitted 1 June, 2020; originally announced June 2020.

Comments: 22 pages, 20 figures, 4 tables

Journal ref: MNRAS 496 (2020), 3808-3828

Towards pseudo-Newtonian black hole jets: comparison of forces

Authors: Maciej Jabłoński, Aleksandra Kotek, Miljenko Čemeljić, Włodek Kluźniak

Abstract: We perform resistive and viscous MHD simulations of axial jet launching from a magnetized thin accretion disk around a supermassive black hole. We compare the forces in our pseudo-Newtonian simulations for the Paczyński-Wiita and Kluźniak-Lee potentials. The results will help to find magnetic field configurations and physical parameters conducive to jet launching. We perform resistive and viscous MHD simulations of axial jet launching from a magnetized thin accretion disk around a supermassive black hole. We compare the forces in our pseudo-Newtonian simulations for the Paczyński-Wiita and Kluźniak-Lee potentials. The results will help to find magnetic field configurations and physical parameters conducive to jet launching. △ Less

Submitted 11 May, 2020; originally announced May 2020.

Comments: 3 pages, 5 figures, contribution to the Proceedings of the XXXIX Polish Astronomical Society Meeting (09-12 September 2019, Olsztyn, Poland)