Showing 1–50 of 242 results for author: Dorner, D

Measurement of $CP$ asymmetries in $B^0 \to K^0_S π^0 γ$ decays at Belle II

Authors: Belle II Collaboration, I. Adachi, L. Aggarwal, H. Ahmed, H. Aihara, N. Akopov, A. Aloisio, N. Anh Ky, D. M. Asner, H. Atmacan, T. Aushev, V. Aushev, M. Aversano, R. Ayad, V. Babu, H. Bae, S. Bahinipati, P. Bambade, Sw. Banerjee, S. Bansal, M. Barrett, J. Baudot, A. Baur, A. Beaubien, F. Becherer , et al. (414 additional authors not shown)

Abstract: We report measurements of time-dependent $CP$ asymmetries in $B^0 \to K^0_S π^0 γ$ decays based on a data sample of $(388\pm6)\times10^6$ $B\bar{B}$ events collected at the $Υ(4S)$ resonance with the Belle II detector. The Belle II experiment operates at the SuperKEKB asymmetric-energy $e^+e^-$ collider. We measure decay-time distributions to determine $CP$-violating parameters $S$ and $C$. We det… ▽ More We report measurements of time-dependent $CP$ asymmetries in $B^0 \to K^0_S π^0 γ$ decays based on a data sample of $(388\pm6)\times10^6$ $B\bar{B}$ events collected at the $Υ(4S)$ resonance with the Belle II detector. The Belle II experiment operates at the SuperKEKB asymmetric-energy $e^+e^-$ collider. We measure decay-time distributions to determine $CP$-violating parameters $S$ and $C$. We determine these parameters for two ranges of $K^0_S π^0$ invariant mass: $m(K^0_S π^0)\in (0.8, 1.0)$ $GeV/c^2$, which is dominated by $B^0 \to K^{*0} (\to K^0_S π^0) γ$ decays, and a complementary region $m(K^0_S π^0)\in (0.6, 0.8)\cup(1.0, 1.8)$ $GeV/c^2$. Our results have improved precision as compared to previous measurements and are consistent with theory predictions. △ Less

Submitted 12 July, 2024; originally announced July 2024.

Comments: 10 pages, 4 figures

Report number: Belle II Preprint 2024-009, KEK Preprint 2024-1

A detailed study of the very-high-energy Crab pulsar emission with the LST-1

Authors: CTA-LST Project, :, K. Abe, S. Abe, A. Abhishek, F. Acero, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González , et al. (272 additional authors not shown)

Abstract: Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov… ▽ More Context: There are currently three pulsars firmly detected by imaging atmospheric Cherenkov telescopes (IACTs), two of them reaching TeV energies, challenging models of very-high-energy (VHE) emission in pulsars. More precise observations are needed to better characterize pulsar emission at these energies. The LST-1 is the prototype of the Large-Sized Telescope, that will be part of the Cherenkov Telescope Array Observatory (CTAO). Its improved performance over previous IACTs makes it well suited for studying pulsars. Aims: To study the Crab pulsar emission with the LST-1, improving and complementing the results from other telescopes. These observations can also be used to characterize the potential of the LST-1 to study other pulsars and detect new ones. Methods: We analyzed a total of $\sim$103 hours of gamma-ray observations of the Crab pulsar conducted with the LST-1 in the period from September 2020 to January 2023. The observations were carried out at zenith angles less than 50 degrees. A new analysis of the Fermi-LAT data was also performed, including $\sim$14 years of observations. Results: The Crab pulsar phaseogram, long-term light-curve, and phase-resolved spectra are reconstructed with the LST-1 from 20 GeV to 450 GeV for P1 and up to 700 GeV for P2. The pulsed emission is detected with a significance of 15.2$σ$. The two characteristic emission peaks of the Crab pulsar are clearly detected (>10$σ$), as well as the so-called bridge emission (5.7$σ$). We find that both peaks are well described by power laws, with spectral indices of $\sim$3.44 and $\sim$3.03 respectively. The joint analysis of Fermi-LAT and LST-1 data shows a good agreement between both instruments in the overlapping energy range. The detailed results obtained in the first observations of the Crab pulsar with LST-1 show the potential that CTAO will have to study this type of sources. △ Less

Submitted 2 July, 2024; originally announced July 2024.

Comments: Accepted by A&A

Constraints on Lorentz invariance violation from the extraordinary Mrk 421 flare of 2014 using a novel analysis method

Authors: MAGIC Collaboration, S. Abe, J. Abhir, A. Abhishek, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, J. Bernete , et al. (192 additional authors not shown)

Abstract: The Lorentz Invariance Violation (LIV), a proposed consequence of certain quantum gravity (QG) scenarios, could instigate an energy-dependent group velocity for ultra-relativistic particles. This energy dependence, although suppressed by the massive QG energy scale $E_\mathrm{QG}$, expected to be on the level of the Planck energy $1.22 \times 10^{19}$ GeV, is potentially detectable in astrophysica… ▽ More The Lorentz Invariance Violation (LIV), a proposed consequence of certain quantum gravity (QG) scenarios, could instigate an energy-dependent group velocity for ultra-relativistic particles. This energy dependence, although suppressed by the massive QG energy scale $E_\mathrm{QG}$, expected to be on the level of the Planck energy $1.22 \times 10^{19}$ GeV, is potentially detectable in astrophysical observations. In this scenario, the cosmological distances traversed by photons act as an amplifier for this effect. By leveraging the observation of a remarkable flare from the blazar Mrk\,421, recorded at energies above 100 GeV by the MAGIC telescopes on the night of April 25 to 26, 2014, we look for time delays scaling linearly and quadratically with the photon energies. Using for the first time in LIV studies a binned-likelihood approach we set constraints on the QG energy scale. For the linear scenario, we set $95\%$ lower limits $E_\mathrm{QG}>2.7\times10^{17}$ GeV for the subluminal case and $E_\mathrm{QG}> 3.6 \times10^{17}$ GeV for the superluminal case. For the quadratic scenario, the $95\%$ lower limits for the subluminal and superluminal cases are $E_\mathrm{QG}>2.6 \times10^{10}$ GeV and $E_\mathrm{QG}>2.5\times10^{10}$ GeV, respectively. △ Less

Submitted 11 June, 2024; originally announced June 2024.

Measurement of the branching fractions of $\bar{B}\to D^{(*)} K^- K^{(*)0}_{(S)}$ and $\bar{B}\to D^{(*)}D_s^{-}$ decays at Belle II

Authors: Belle II Collaboration, I. Adachi, L. Aggarwal, H. Aihara, N. Akopov, A. Aloisio, N. Althubiti, N. Anh Ky, D. M. Asner, H. Atmacan, T. Aushev, V. Aushev, M. Aversano, R. Ayad, V. Babu, H. Bae, S. Bahinipati, P. Bambade, Sw. Banerjee, S. Bansal, M. Barrett, J. Baudot, A. Baur, A. Beaubien, F. Becherer , et al. (382 additional authors not shown)

Abstract: We present measurements of the branching fractions of eight $\overline B{}^0\to D^{(*)+} K^- K^{(*)0}_{(S)}$, $B^{-}\to D^{(*)0} K^- K^{(*)0}_{(S)}$ decay channels. The results are based on data from SuperKEKB electron-positron collisions at the $Υ(4S)$ resonance collected with the Belle II detector, corresponding to an integrated luminosity of $362~\text{fb}^{-1}$. The event yields are extracted… ▽ More We present measurements of the branching fractions of eight $\overline B{}^0\to D^{(*)+} K^- K^{(*)0}_{(S)}$, $B^{-}\to D^{(*)0} K^- K^{(*)0}_{(S)}$ decay channels. The results are based on data from SuperKEKB electron-positron collisions at the $Υ(4S)$ resonance collected with the Belle II detector, corresponding to an integrated luminosity of $362~\text{fb}^{-1}$. The event yields are extracted from fits to the distributions of the difference between expected and observed $B$ meson energy, and are efficiency-corrected as a function of $m(K^-K^{(*)0}_{(S)})$ and $m(D^{(*)}K^{(*)0}_{(S)})$ in order to avoid dependence on the decay model. These results include the first observation of $\overline B{}^0\to D^+K^-K_S^0$, $B^-\to D^{*0}K^-K_S^0$, and $\overline B{}^0\to D^{*+}K^-K_S^0$ decays and a significant improvement in the precision of the other channels compared to previous measurements. The helicity-angle distributions and the invariant mass distributions of the $K^- K^{(*)0}_{(S)}$ systems are compatible with quasi-two-body decays via a resonant transition with spin-parity $J^P=1^-$ for the $K^-K_S^0$ systems and $J^P= 1^+$ for the $K^-K^{*0}$ systems. We also present measurements of the branching fractions of four $\overline B{}^0\to D^{(*)+} D_s^-$, $B^{-}\to D^{(*)0} D_s^- $ decay channels with a precision compatible to the current world averages. △ Less

Submitted 10 June, 2024; originally announced June 2024.

Comments: Prepared for submission to JHEP. 34 pages, 14 figures

Report number: Belle II Preprint: 2024-014, KEK Preprint: 2024-8

Search for the decay $B^{0}\toγγ$ using Belle and Belle II data

Authors: Belle, Belle II Collaborations, :, I. Adachi, L. Aggarwal, H. Aihara, N. Akopov, A. Aloisio, S. Al Said, N. Althubiti, N. Anh Ky, D. M. Asner, H. Atmacan, T. Aushev, V. Aushev, M. Aversano, R. Ayad, V. Babu, H. Bae, S. Bahinipati, P. Bambade, Sw. Banerjee, S. Bansal, M. Barrett, J. Baudot , et al. (385 additional authors not shown)

Abstract: We report the result of a search for the rare decay $B^{0} \to γγ$ using a combined dataset of $753\times10^{6}$ $B\bar{B}$ pairs collected by the Belle experiment and $387\times10^{6}$ $B\bar{B}$ pairs collected by the Belle II experiment from decays of the $\rm Υ(4S)$ resonance produced in $e^{+}e^{-}$ collisions. A simultaneous fit to the Belle and Belle II data sets yields… ▽ More We report the result of a search for the rare decay $B^{0} \to γγ$ using a combined dataset of $753\times10^{6}$ $B\bar{B}$ pairs collected by the Belle experiment and $387\times10^{6}$ $B\bar{B}$ pairs collected by the Belle II experiment from decays of the $\rm Υ(4S)$ resonance produced in $e^{+}e^{-}$ collisions. A simultaneous fit to the Belle and Belle II data sets yields $11.0^{+6.5}_{-5.5}$ signal events, corresponding to a 2.5$σ$ significance. We determine the branching fraction $\mathcal{B}(B^{0} \to γγ) = (3.7^{+2.2}_{-1.8}(\rm stat)\pm0.5(\rm syst))\times10^{-8}$ and set a 90% credibility level upper limit of $\mathcal{B}(B^{0} \to γγ) < 6.4\times10^{-8}$. △ Less

Submitted 30 May, 2024; originally announced May 2024.

Report number: Belle II Preprint: 2024-017, KEK Preprint: 2024-13

Measurement of the energy dependence of the $e^+e^- \to B\bar{B}$, $B\bar{B}{}^*$, and $B^*\bar{B}{}^*$ cross sections at Belle~II

Authors: Belle II Collaboration, I. Adachi, L. Aggarwal, H. Ahmed, H. Aihara, N. Akopov, A. Aloisio, N. Althubiti, N. Anh Ky, D. M. Asner, H. Atmacan, T. Aushev, V. Aushev, M. Aversano, R. Ayad, V. Babu, H. Bae, S. Bahinipati, P. Bambade, Sw. Banerjee, S. Bansal, M. Barrett, J. Baudot, M. Bauer, A. Baur , et al. (444 additional authors not shown)

Abstract: We report measurements of the $e^+e^- \to B\bar{B}$, $B\bar{B}{}^*$, and $B^*\bar{B}{}^*$ cross sections at four energies, 10653, 10701, 10746 and 10805 MeV, using data collected by the Belle~II experiment. We reconstruct one $B$ meson in a large number of hadronic final states and use its momentum to identify the production process. In the first $2-5$ MeV above $B^*\bar{B}{}^*$ threshold, the… ▽ More We report measurements of the $e^+e^- \to B\bar{B}$, $B\bar{B}{}^*$, and $B^*\bar{B}{}^*$ cross sections at four energies, 10653, 10701, 10746 and 10805 MeV, using data collected by the Belle~II experiment. We reconstruct one $B$ meson in a large number of hadronic final states and use its momentum to identify the production process. In the first $2-5$ MeV above $B^*\bar{B}{}^*$ threshold, the $e^+e^- \to B^*\bar{B}{}^*$ cross section increases rapidly. This may indicate the presence of a pole close to the threshold. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 30 pages, 15 figures, submitted to JHEP

Report number: Belle II Preprint 2024-016, KEK Preprint 2024-12

Test of light-lepton universality in $τ$ decays with the Belle II experiment

Authors: Belle II Collaboration, I. Adachi, K. Adamczyk, L. Aggarwal, H. Aihara, N. Akopov, A. Aloisio, N. Anh Ky, D. M. Asner, H. Atmacan, V. Aushev, M. Aversano, R. Ayad, V. Babu, H. Bae, S. Bahinipati, P. Bambade, Sw. Banerjee, S. Bansal, M. Barrett, J. Baudot, A. Baur, A. Beaubien, F. Becherer, J. Becker , et al. (406 additional authors not shown)

Abstract: We present a measurement of the ratio $R_μ= \mathcal{B}(τ^-\to μ^-\barν_μν_τ) / \mathcal{B}(τ^-\to e^-\barν_eν_τ)$ of branching fractions $\mathcal{B}$ of the $τ$ lepton decaying to muons or electrons using data collected with the Belle II detector at the SuperKEKB $e^+e^-$ collider. The sample has an integrated luminosity of 362 fb$^{-1}$ at a centre-of-mass energy of 10.58 GeV. Using an optimise… ▽ More We present a measurement of the ratio $R_μ= \mathcal{B}(τ^-\to μ^-\barν_μν_τ) / \mathcal{B}(τ^-\to e^-\barν_eν_τ)$ of branching fractions $\mathcal{B}$ of the $τ$ lepton decaying to muons or electrons using data collected with the Belle II detector at the SuperKEKB $e^+e^-$ collider. The sample has an integrated luminosity of 362 fb$^{-1}$ at a centre-of-mass energy of 10.58 GeV. Using an optimised event selection, a binned maximum likelihood fit is performed using the momentum spectra of the electron and muon candidates. The result, $R_μ= 0.9675 \pm 0.0007 \pm 0.0036$, where the first uncertainty is statistical and the second is systematic, is the most precise to date. It provides a stringent test of the light-lepton universality, translating to a ratio of the couplings of the muon and electron to the $W$ boson in $τ$ decays of $0.9974 \pm 0.0019$, in agreement with the standard model expectation of unity. △ Less

Submitted 23 May, 2024; originally announced May 2024.

Report number: Belle II Preprint 2024-002, KEK Preprint 2023-49

Search for lepton-flavor-violating $��^- \to μ^-μ^+μ^-$ decays at Belle II

Authors: Belle II Collaboration, I. Adachi, L. Aggarwal, H. Aihara, N. Akopov, A. Aloisio, N. Althubiti, N. Anh Ky, D. M. Asner, H. Atmacan, V. Aushev, M. Aversano, R. Ayad, V. Babu, H. Bae, S. Bahinipati, P. Bambade, Sw. Banerjee, S. Bansal, M. Barrett, J. Baudot, A. Baur, A. Beaubien, F. Becherer, J. Becker , et al. (407 additional authors not shown)

Abstract: We present the result of a search for the charged-lepton-flavor violating decay $τ^- \to μ^-μ^+μ^-$ using a $424fb^{-1}$ sample of data recorded by the Belle II experiment at the SuperKEKB $e^{-}e^{+}$ collider. The selection of $e^{-}e^{+}\toτ^+τ^-$ events is based on an inclusive reconstruction of the non-signal tau decay, and on a boosted decision tree to suppress background. We observe one sig… ▽ More We present the result of a search for the charged-lepton-flavor violating decay $τ^- \to μ^-μ^+μ^-$ using a $424fb^{-1}$ sample of data recorded by the Belle II experiment at the SuperKEKB $e^{-}e^{+}$ collider. The selection of $e^{-}e^{+}\toτ^+τ^-$ events is based on an inclusive reconstruction of the non-signal tau decay, and on a boosted decision tree to suppress background. We observe one signal candidate, which is compatible with the expectation from background processes. We set a $90\%$ confidence level upper limit of $1.9 \times 10^{-8}$ on the branching fraction of the \taumu decay, which is the most stringent bound to date. △ Less

Submitted 12 May, 2024; originally announced May 2024.

Report number: Belle II Preprint 2024-012 KEK Preprint 2024-6

Measurement of the branching fraction of the decay $B^- \to D^0 ρ(770)^-$ at Belle II

Authors: Belle II Collaboration, I. Adachi, L. Aggarwal, H. Aihara, N. Akopov, A. Aloisio, N. Anh Ky, D. M. Asner, H. Atmacan, V. Aushev, M. Aversano, R. Ayad, V. Babu, H. Bae, S. Bahinipati, P. Bambade, Sw. Banerjee, S. Bansal, M. Barrett, J. Baudot, A. Baur, A. Beaubien, F. Becherer, J. Becker, J. V. Bennett , et al. (367 additional authors not shown)

Abstract: We measure the branching fraction of the decay $B^- \to D^0 ρ(770)^-$ using data collected with the Belle II detector. The data contain 387 million $B\overline{B}$ pairs produced in $e^+e^-$ collisions at the $Υ(4S)$ resonance. We reconstruct $8360\pm 180$ decays from an analysis of the distributions of the $B^-$ energy and the $ρ(770)^-$ helicity angle. We determine the branching fraction to be… ▽ More We measure the branching fraction of the decay $B^- \to D^0 ρ(770)^-$ using data collected with the Belle II detector. The data contain 387 million $B\overline{B}$ pairs produced in $e^+e^-$ collisions at the $Υ(4S)$ resonance. We reconstruct $8360\pm 180$ decays from an analysis of the distributions of the $B^-$ energy and the $ρ(770)^-$ helicity angle. We determine the branching fraction to be $(0.939 \pm 0.021\mathrm{(stat)} \pm 0.050\mathrm{(syst)})\%$, in agreement with previous results. Our measurement improves the relative precision of the world average by more than a factor of two. △ Less

Submitted 27 June, 2024; v1 submitted 16 April, 2024; originally announced April 2024.

Report number: Belle II Preprint 2024-011, KEK Preprint 2024-4

Journal ref: PRD 109, 111103 (2024)

Dark Matter Line Searches with the Cherenkov Telescope Array

Authors: S. Abe, J. Abhir, A. Abhishek, 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, G. Ambrosi, L. Angel, C. Aramo, C. Arcaro, T. T. H. Arnesen, L. Arrabito, K. Asano, Y. Ascasibar, J. Aschersleben, H. Ashkar , et al. (540 additional authors not shown)

Abstract: Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of sele… ▽ More Monochromatic gamma-ray signals constitute a potential smoking gun signature for annihilating or decaying dark matter particles that could relatively easily be distinguished from astrophysical or instrumental backgrounds. We provide an updated assessment of the sensitivity of the Cherenkov Telescope Array (CTA) to such signals, based on observations of the Galactic centre region as well as of selected dwarf spheroidal galaxies. We find that current limits and detection prospects for dark matter masses above 300 GeV will be significantly improved, by up to an order of magnitude in the multi-TeV range. This demonstrates that CTA will set a new standard for gamma-ray astronomy also in this respect, as the world's largest and most sensitive high-energy gamma-ray observatory, in particular due to its exquisite energy resolution at TeV energies and the adopted observational strategy focussing on regions with large dark matter densities. Throughout our analysis, we use up-to-date instrument response functions, and we thoroughly model the effect of instrumental systematic uncertainties in our statistical treatment. We further present results for other potential signatures with sharp spectral features, e.g.~box-shaped spectra, that would likewise very clearly point to a particle dark matter origin. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 43 pages JCAP style (excluding author list and references), 19 figures

The variability patterns of the TeV blazar PG 1553+113 from a decade of MAGIC and multi-band observations

Authors: MAGIC Collaboration, H. Abe, S. Abe, J. Abhir, V. A. Acciari, I. Agudo, 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, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari , et al. (242 additional authors not shown)

Abstract: PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m a… ▽ More PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

Rapid Variability of Mrk 421 During Extreme Flaring as Seen Through the Eyes of XMM-Newton

Authors: A. Gokus, J. Wilms, M. Kadler, D. Dorner, M. A. Nowak, A. Kreikenbohm, K. Leiter, T. Bretz, B. Schleicher, A. G. Markowitz, K. Pottschmidt, K. Mannheim, I. Kreykenbohm, M. Langejahn, F. McBride, T. Beuchert, T. Dauser, M. Kreter, J. Abhir, D. Baack, M. Balbo, A. Biland, K. Brand, J. Buss, L. Eisenberger , et al. (9 additional authors not shown)

Abstract: By studying the variability of blazars across the electromagnetic spectrum, it is possible to resolve the underlying processes responsible for rapid flux increases, so-called flares. We report on an extremely bright X-ray flare in the high-peaked BL Lacertae object Mrk 421 that occurred simultaneously with enhanced $γ$-ray activity detected at very high energies (VHE) by FACT on 2019 June 9. We tr… ▽ More By studying the variability of blazars across the electromagnetic spectrum, it is possible to resolve the underlying processes responsible for rapid flux increases, so-called flares. We report on an extremely bright X-ray flare in the high-peaked BL Lacertae object Mrk 421 that occurred simultaneously with enhanced $γ$-ray activity detected at very high energies (VHE) by FACT on 2019 June 9. We triggered an observation with XMM-Newton, which observed the source quasi-continuously for 25 hours. We find that the source was in the brightest state ever observed using XMM-Newton, reaching a flux of $2.8\times10^{-9}$ erg cm$^{-2}$ s$^{-1}$ over an energy range of 0.3 - 10 keV. We perform a spectral and timing analysis to reveal the mechanisms of particle acceleration and to search for the shortest source-intrinsic time scales. Mrk 421 exhibits the typical harder-when-brighter behaviour throughout the observation and shows a clock-wise hysteresis pattern, which indicates that the cooling dominates over the acceleration process. While the X-ray emission in different sub-bands is highly correlated, we can exclude large time lags as the computed zDCFs are consistent with a zero lag. We find rapid variability on time scales of 1 ks for the 0.3 - 10 keV band and down to 300s in the hard X-ray band (4 - 10 keV). Taking these time scales into account, we discuss different models to explain the observed X-ray flare, and find that a plasmoid-dominated magnetic reconnection process is able to describe our observation best. △ Less

Submitted 29 February, 2024; originally announced February 2024.

Comments: 13 pages, 7 figures, accepted for publication in MNRAS

Performance and first measurements of the MAGIC Stellar Intensity Interferometer

Authors: MAGIC Collaboration, S. Abe, J. Abhir, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, E. Bernardini, M. Bernardos, J. Bernete, A. Berti , et al. (195 additional authors not shown)

Abstract: In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the sys… ▽ More In recent years, a new generation of optical intensity interferometers has emerged, leveraging the existing infrastructure of Imaging Atmospheric Cherenkov Telescopes (IACTs). The MAGIC telescopes host the MAGIC-SII system (Stellar Intensity Interferometer), implemented to investigate the feasibility and potential of this technique on IACTs. After the first successful measurements in 2019, the system was upgraded and now features a real-time, dead-time-free, 4-channel, GPU-based correlator. These hardware modifications allow seamless transitions between MAGIC's standard very-high-energy gamma-ray observations and optical interferometry measurements within seconds. We establish the feasibility and potential of employing IACTs as competitive optical Intensity Interferometers with minimal hardware adjustments. The measurement of a total of 22 stellar diameters are reported, 9 corresponding to reference stars with previous comparable measurements, and 13 with no prior measurements. A prospective implementation involving telescopes from the forthcoming Cherenkov Telescope Array Observatory's northern hemisphere array, such as the first prototype of its Large-Sized Telescopes, LST-1, is technically viable. This integration would significantly enhance the sensitivity of the current system and broaden the UV-plane coverage. This advancement would enable the system to achieve competitive sensitivity with the current generation of long-baseline optical interferometers over blue wavelengths. △ Less

Submitted 7 February, 2024; originally announced February 2024.

Comments: 18 pages, 13 figures, submitted to MNRAS

Insights into the broad-band emission of the TeV blazar Mrk 501 during the first X-ray polarization measurements

Authors: S. Abe, J. Abhir, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, A. Bautista, J. Baxter, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, J. Bernete, A. Berti, J. Besenrieder , et al. (239 additional authors not shown)

Abstract: We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and… ▽ More We present the first multi-wavelength study of Mrk 501 including very-high-energy (VHE) gamma-ray observations simultaneous to X-ray polarization measurements from the Imaging X-ray Polarimetry Explorer (IXPE). We use radio-to-VHE data from a multi-wavelength campaign organized between 2022-03-01 and 2022-07-19. The observations were performed by MAGIC, Fermi-LAT, NuSTAR, Swift (XRT and UVOT), and several instruments covering the optical and radio bands. During the IXPE pointings, the VHE state is close to the average behavior with a 0.2-1 TeV flux of 20%-50% the emission of the Crab Nebula. Despite the average VHE activity, an extreme X-ray behavior is measured for the first two IXPE pointings in March 2022 with a synchrotron peak frequency >1 keV. For the third IXPE pointing in July 2022, the synchrotron peak shifts towards lower energies and the optical/X-ray polarization degrees drop. The X-ray polarization is systematically higher than at lower energies, suggesting an energy-stratification of the jet. While during the IXPE epochs the polarization angle in the X-ray, optical and radio bands align well, we find a clear discrepancy in the optical and radio polarization angles in the middle of the campaign. We model the broad-band spectra simultaneous to the IXPE pointings assuming a compact zone dominating in the X-rays and VHE, and an extended zone stretching further downstream the jet dominating the emission at lower energies. NuSTAR data allow us to precisely constrain the synchrotron peak and therefore the underlying electron distribution. The change between the different states observed in the three IXPE pointings can be explained by a change of magnetization and/or emission region size, which directly connects the shift of the synchrotron peak to lower energies with the drop in polarization degree. △ Less

Submitted 16 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomy & Astrophysics. 19 pages, 9 figures. Corresponding authors: Lea Heckmann, Axel Arbet Engels, David Paneque

Constraints on axion-like particles with the Perseus Galaxy Cluster with MAGIC

Authors: MAGIC Collaboration, H. Abe, S. Abe, J. Abhir, V. A. Acciari, I. Agudo, 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, J. Bernete, A. Berti , et al. (189 additional authors not shown)

Abstract: Axion-like particles (ALPs) are pseudo-Nambu-Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 hrs of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the r… ▽ More Axion-like particles (ALPs) are pseudo-Nambu-Goldstone bosons that emerge in various theories beyond the standard model. These particles can interact with high-energy photons in external magnetic fields, influencing the observed gamma-ray spectrum. This study analyzes 41.3 hrs of observational data from the Perseus Galaxy Cluster collected with the MAGIC telescopes. We focused on the spectra the radio galaxy in the center of the cluster: NGC 1275. By modeling the magnetic field surrounding this target, we searched for spectral indications of ALP presence. Despite finding no statistical evidence of ALP signatures, we were able to exclude ALP models in the sub-micro electronvolt range. Our analysis improved upon previous work by calculating the full likelihood and statistical coverage for all considered models across the parameter space. Consequently, we achieved the most stringent limits to date for ALP masses around 50 neV, with cross sections down to $g_{aγ} = 3 \times 10^{-12}$ GeV$^{-1}$. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: 25 pages, 10 figures, accepted for publication in Physics of the Dark Universe

TELAMON: Effelsberg monitoring of AGN jets with very-high-energy astroparticle emission -- I. Program description and sample characterization

Authors: F. Eppel, M. Kadler, J. Heßdörfer, P. Benke, L. Debbrecht, J. Eich, A. Gokus, S. Hämmerich, D. Kirchner, G. F. Paraschos, F. Rösch, W. Schulga, J. Sinapius, P. Weber, U. Bach, D. Dorner, P. G. Edwards, M. Giroletti, A. Kraus, O. Hervet, S. Koyama, T. P. Krichbaum, K. Mannheim, E. Ros, M. Zacharias , et al. (1 additional authors not shown)

Abstract: Aims. We introduce the TELAMON program which is using the Effelsberg 100-m telescope to monitor the radio spectra of active galactic nuclei (AGN) under scrutiny in astroparticle physics, specifically TeV blazars and candidate neutrino-associated AGN. Here, we present and characterize our main sample of TeV-detected blazars. Methods. We analyze the data sample from the first ~2.5 years of observa… ▽ More Aims. We introduce the TELAMON program which is using the Effelsberg 100-m telescope to monitor the radio spectra of active galactic nuclei (AGN) under scrutiny in astroparticle physics, specifically TeV blazars and candidate neutrino-associated AGN. Here, we present and characterize our main sample of TeV-detected blazars. Methods. We analyze the data sample from the first ~2.5 years of observations between August 2020 and February 2023 in the range from 14 GHz to 45 GHz. During this pilot phase, we have observed all 59 TeV-detected blazars in the Northern Hemisphere (i.e., Dec. >0°) known at the time of observation. We discuss the basic data reduction and calibration procedures used for all TELAMON data and introduce a sub-band averaging method used to calculate average light curves for the sources in our sample. Results. The TeV-selected sources in our sample exhibit a median flux density of 0.12 Jy at 20 mm, 0.20 Jy at 14 mm and 0.60 Jy at 7 mm. The spectrum for most of the sources is consistent with a flat radio spectrum and we find a median spectral index ($S(ν)\proptoν^α$) of $α=-0.11$. Our results on flux density and spectral index are consistent with previous studies of TeV-selected blazars. Compared to the GeV-selected F-GAMMA sample, TELAMON sources are significantly fainter in the radio band. This is consistent with the double-humped spectrum of blazars being shifted towards higher frequencies for TeV-emitters (in particular for high-synchrotron peaked BL Lac type objects), which results in a lower radio flux density. The spectral index distribution of our TeV-selected blazar sample is not significantly different from the GeV-selected F-GAMMA sample. Moreover, we present a strategy to track the light curve evolution of sources in our sample for future variability and correlation analysis. △ Less

Submitted 11 January, 2024; originally announced January 2024.

Comments: 16 pages, 11 figures, 4 tables, accepted for publication in A&A

First characterization of the emission behavior of Mrk421 from radio to VHE gamma rays with simultaneous X-ray polarization measurements

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

Abstract: We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured… ▽ More We perform the first broadband study of Mrk421 from radio to TeV gamma rays with simultaneous measurements of the X-ray polarization from IXPE. The data were collected within an extensive multiwavelength campaign organized between May and June 2022 using MAGIC, Fermi-LAT, NuSTAR, XMM-Newton, Swift, and several optical and radio telescopes to complement IXPE. During the IXPE exposures, the measured 0.2-1 TeV flux is close to the quiescent state and ranges from 25% to 50% of the Crab Nebula without intra-night variability. Throughout the campaign, the VHE and X-ray emission are positively correlated at a $4σ$ significance level. The IXPE measurements unveil a X-ray polarization degree that is a factor of 2-5 higher than in the optical/radio bands; that implies an energy-stratified jet in which the VHE photons are emitted co-spatially with the X-rays, in the vicinity of a shock front. The June 2022 observations exhibit a rotation of the X-ray polarization angle. Despite no simultaneous VHE coverage being available during a large fraction of the swing, the Swift-XRT monitoring unveils an X-ray flux increase with a clear spectral hardening. It suggests that flares in high synchrotron peaked blazars can be accompanied by a polarization angle rotation, as observed in some flat spectrum radio quasars. Finally, during the polarization angle rotation, NuSTAR data reveal two contiguous spectral hysteresis loops in opposite directions (clockwise and counter-clockwise), implying important changes in the particle acceleration efficiency on $\sim$hour timescales. △ Less

Submitted 17 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics. 29 pages, 22 figures. Corresponding authors: Axel Arbet Engels, Felix Schmuckermaier, David Paneque

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

MAGIC detection of GRB 201216C at $z=1.1$

Authors: H. Abe, S. Abe, V. A. Acciari, I. Agudo, 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, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari , et al. (195 additional authors not shown)

Abstract: Gamma-ray bursts (GRBs) are explosive transient events occurring at cosmological distances, releasing a large amount of energy as electromagnetic radiation over several energy bands. We report the detection of the long GRB~201216C by the MAGIC telescopes. The source is located at $z=1.1$ and thus it is the farthest one detected at very high energies. The emission above \SI{70}{\GeV} of GRB~201216C… ▽ More Gamma-ray bursts (GRBs) are explosive transient events occurring at cosmological distances, releasing a large amount of energy as electromagnetic radiation over several energy bands. We report the detection of the long GRB~201216C by the MAGIC telescopes. The source is located at $z=1.1$ and thus it is the farthest one detected at very high energies. The emission above \SI{70}{\GeV} of GRB~201216C is modelled together with multi-wavelength data within a synchrotron and synchrotron-self Compton (SSC) scenario. We find that SSC can explain the broadband data well from the optical to the very-high-energy band. For the late-time radio data, a different component is needed to account for the observed emission. Differently from previous GRBs detected in the very-high-energy range, the model for GRB~201216C strongly favors a wind-like medium. The model parameters have values similar to those found in past studies of the afterglows of GRBs detected up to GeV energies. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: 13 pages, 6 figures, 2 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 3, January 2024, Pages 5856-5867

Multi-year characterisation of the broad-band emission from the intermittent extreme BL Lac 1ES~2344+514

Authors: H. Abe, S. Abe, V. A. Acciari, I. Agudo, 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, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari, A. Biland, O. Blanch , et al. (210 additional authors not shown)

Abstract: The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, f… ▽ More The BL Lac 1ES 2344+514 is known for temporary extreme properties (e.g., a shift of the synchrotron SED peak energy $ν_{synch,p}$ above 1keV). While those extreme states were so far observed only during high flux levels, additional multi-year observing campaigns are required to achieve a coherent picture. Here, we report the longest investigation of the source from radio to VHE performed so far, focusing on a systematic characterisation of the intermittent extreme states. While our results confirm that 1ES 2344+514 typically exhibits $ν_{synch,p}>$1keV during elevated flux periods, we also find periods where the extreme state coincides with low flux activity. A strong spectral variability thus happens in the quiescent state, and is likely caused by an increase of the electron acceleration efficiency without a change in the electron injection luminosity. We also report a strong X-ray flare (among the brightest for 1ES 2344+514) without a significant shift of $ν_{synch,p}$. During this particular flare, the X-ray spectrum is among the softest of the campaign. It unveils complexity in the spectral evolution, where the common harder-when-brighter trend observed in BL Lacs is violated. During a low and hard X-ray state, we find an excess of the UV flux with respect to an extrapolation of the X-ray spectrum to lower energies. This UV excess implies that at least two regions contribute significantly to the infrared/optical/ultraviolet/X-ray emission. Using the simultaneous MAGIC, XMM-Newton, NuSTAR, and AstroSat observations, we argue that a region possibly associated with the 10 GHz radio core may explain such an excess. Finally, we investigate a VHE flare, showing an absence of simultaneous variability in the 0.3-2keV band. Using a time-dependent leptonic modelling, we show that this behaviour, in contradiction to single-zone scenarios, can instead be explained by a two-component model. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 682, A114 (2024)

Performance of the joint LST-1 and MAGIC observations evaluated with Crab Nebula data

Authors: H. Abe, K. Abe, S. Abe, V. A. Acciari, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, T. Aniello, S. Ansoldi, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, D. Baack, A. Babić, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batković , et al. (344 additional authors not shown)

Abstract: Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both syste… ▽ More Aims. LST-1, the prototype of the Large-Sized Telescope for the upcoming Cherenkov Telescope Array Observatory, is concluding its commissioning in Observatorio del Roque de los Muchachos on the island of La Palma. The proximity of LST-1 (Large-Sized Telescope 1) to the two MAGIC (Major Atmospheric Gamma Imaging Cherenkov) telescopes permits observations of the same gamma-ray events with both systems. Methods. We describe the joint LST-1+MAGIC analysis pipeline and use simultaneous Crab Nebula observations and Monte Carlo simulations to assess the performance of the three-telescope system. The addition of the LST-1 telescope allows the recovery of events in which one of the MAGIC images is too dim to survive analysis quality cuts. Results. Thanks to the resulting increase in the collection area and stronger background rejection, we find a significant improvement in sensitivity, allowing the detection of 30% weaker fluxes in the energy range between 200 GeV and 3 TeV. The spectrum of the Crab Nebula, reconstructed in the energy range ~60 GeV to ~10 TeV, is in agreement with previous measurements. △ Less

Submitted 3 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 680, A66 (2023)

Joint searches by FACT, H.E.S.S., MAGIC and VERITAS for VHE gamma-ray emission associated with neutrinos detected by IceCube

Authors: Fabian Schüssler, Halim Ashkar, Elisa Bernardini, Alessio Berti, Federica Bradascio, Sara Buson, Daniela Dorner, Weidong Jin, Gasper Kukec Mezek, Marcos Santander, Konstancja Satalecka, Bernd Schleicher, Mohanraj Senniappan, Ilaria Viale

Abstract: The sources of the astrophysical flux of high-energy neutrinos detected by IceCube are still largely unknown, but searches for temporal and spatial correlation between neutrinos and electromagnetic radiation are a promising approach in this endeavor. All major imaging atmospheric Cherenkov telescopes (IACTs) - FACT, H.E.S.S., MAGIC, and VERITAS - operate an active follow-up program of target-of-op… ▽ More The sources of the astrophysical flux of high-energy neutrinos detected by IceCube are still largely unknown, but searches for temporal and spatial correlation between neutrinos and electromagnetic radiation are a promising approach in this endeavor. All major imaging atmospheric Cherenkov telescopes (IACTs) - FACT, H.E.S.S., MAGIC, and VERITAS - operate an active follow-up program of target-of-opportunity observations of neutrino alerts issued by IceCube. These programs use several complementary neutrino alert streams. A publicly distributed alert stream is formed by individual high-energy neutrino candidate events of potentially astrophysical origin, such as IceCube-170922A (which could be linked to the flaring blazar TXS\,0506+056). A privately distributed alert stream is formed by clusters of neutrino events in time and space around either pre-selected gamma-ray sources or anywhere in the sky. Here, we present joint searches for multi-wavelength emission associated with a set of IceCube alerts, both private and public, received through mid-January 2021. We will give an overview of the programs of the participating IACTs. We will showcase the various follow-up and data analysis strategies employed in response to the different alert types and various possible counterpart scenarios. Finally, we will present results from a combined analysis of the VHE gamma-ray observations obtained across all involved instruments, as well as relevant multi-wavelength data. △ Less

Submitted 27 September, 2023; originally announced September 2023.

Comments: Proceedings 38th International Cosmic Ray Conference (ICRC2023)

TELAMON: Effelsberg Monitoring of AGN Jets with Very-High-Energy Astroparticle Emissions -- Polarization properties

Authors: J. Heßdörfer, M. Kadler, P. Benke, L. Debbrecht, J. Eich, F. Eppel, A. Gokus, S. Hämmerich, D. Kirchner, G. F. Paraschos, F. Rösch, W. Schulga, J. Sinapius, P. Weber, U. Bach, D. Berge, S. Buson, D. Dorner, P. G. Edwards, C. M. Fromm, M. Giroletti, O. Hervet, A. Kappes, S. Koyama, A. Kraus , et al. (11 additional authors not shown)

Abstract: We present recent results of the TELAMON program, which is using the Effelsberg 100-m telescope to monitor the radio spectra of active galactic nuclei (AGN) under scrutiny in astroparticle physics, namely TeV blazars and neutrino-associated AGN. Our sample includes all known Northern TeV-emitting blazars as well as blazars positionally coincident with IceCube neutrino alerts. Polarization can give… ▽ More We present recent results of the TELAMON program, which is using the Effelsberg 100-m telescope to monitor the radio spectra of active galactic nuclei (AGN) under scrutiny in astroparticle physics, namely TeV blazars and neutrino-associated AGN. Our sample includes all known Northern TeV-emitting blazars as well as blazars positionally coincident with IceCube neutrino alerts. Polarization can give additional insight into the source properties, as the polarized emission is often found to vary on different timescales and amplitudes than the total intensity emission. Here, we present an overview of the polarization properties of the TeV-emitting TELAMON sources at four frequencies in the 20 mm and 7 mm bands. While at 7 mm roughly $82\,\%$ of all observed sources are found to be significantly polarized, for 20 mm the percentage is $\sim58\,\%$. We find that most of the sources exhibit mean fractional polarizations of $<5\%$, matching the expectations of rather low polarization levels in these sources from previous studies at lower radio frequencies. Nevertheless, we demonstrate examples of how the polarized emission can provide additional information over the total intensity. △ Less

Submitted 31 October, 2023; v1 submitted 26 September, 2023; originally announced September 2023.

Comments: 9 pages, 2 figures, Proceedings of the 38th International Cosmic Ray Conference (ICRC2023)

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

Observations of the Crab Nebula and Pulsar with the Large-Sized Telescope Prototype of the Cherenkov Telescope Array

Authors: CTA-LST Project, :, H. Abe, K. Abe, S. Abe, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, C. Arcaro, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González, E. Bernardini , et al. (467 additional authors not shown)

Abstract: CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing a… ▽ More CTA (Cherenkov Telescope Array) is the next generation ground-based observatory for gamma-ray astronomy at very-high energies. The Large-Sized Telescope prototype (LST-1) is located at the Northern site of CTA, on the Canary Island of La Palma. LSTs are designed to provide optimal performance in the lowest part of the energy range covered by CTA, down to $\simeq 20$ GeV. LST-1 started performing astronomical observations in November 2019, during its commissioning phase, and it has been taking data since then. We present the first LST-1 observations of the Crab Nebula, the standard candle of very-high energy gamma-ray astronomy, and use them, together with simulations, to assess the basic performance parameters of the telescope. The data sample consists of around 36 hours of observations at low zenith angles collected between November 2020 and March 2022. LST-1 has reached the expected performance during its commissioning period - only a minor adjustment of the preexisting simulations was needed to match the telescope behavior. The energy threshold at trigger level is estimated to be around 20 GeV, rising to $\simeq 30$ GeV after data analysis. Performance parameters depend strongly on energy, and on the strength of the gamma-ray selection cuts in the analysis: angular resolution ranges from 0.12 to 0.40 degrees, and energy resolution from 15 to 50%. Flux sensitivity is around 1.1% of the Crab Nebula flux above 250 GeV for a 50-h observation (12% for 30 minutes). The spectral energy distribution (in the 0.03 - 30 TeV range) and the light curve obtained for the Crab Nebula agree with previous measurements, considering statistical and systematic uncertainties. A clear periodic signal is also detected from the pulsar at the center of the Nebula. △ Less

Submitted 19 July, 2023; v1 submitted 22 June, 2023; originally announced June 2023.

Comments: Accepted in ApJ. v3: updated author list and acknowledgements, fixed typos and other minor issues

Probing AGN variability with the Cherenkov Telescope Array

Authors: F. Cangemi, T. Hovatta, E. Lindfors, M. Cerruti, J. Becerra-Gonzalez, J. Biteau, C. Boisson, M. Böttcher, E. de Gouveia Dal Pino, D. Dorner, G. Grolleron, J. -P. Lenain, M. Manganaro, W. Max-Moerbeck, P. Morris, K. Nilsson, L. Passos Reis, P. Romano, O. Sergijenko, F. Tavecchio, S. Vercellone, S. Wagner, M. Zacharias

Abstract: Relativistic jets launched by Active Galactic Nuclei are among the most powerful particle accelerators in the Universe. The emission over the entire electromagnetic spectrum of these relativistic jets can be extremely variable with scales of variability from less than few minutes up to several years. These variability patterns, which can be very complex, contain information about the acceleration… ▽ More Relativistic jets launched by Active Galactic Nuclei are among the most powerful particle accelerators in the Universe. The emission over the entire electromagnetic spectrum of these relativistic jets can be extremely variable with scales of variability from less than few minutes up to several years. These variability patterns, which can be very complex, contain information about the acceleration processes of the particles and the area(s) of emission. Thanks to its sensitivity, five-to twenty-times better than the current generation of Imaging Atmospheric Cherenkov Telescopes depending on energy, the Cherenkov Telescope Array will be able to follow the emission from these objects with a very accurate time sampling and over a wide spectral coverage from 20 GeV to > 20 TeV and thus reveal the nature of the acceleration processes at work in these objects. We will show the first results of our lightcurve simulations and long-term behavior of AGN as will be observed by CTA, based on state-of-art particle acceleration models. △ Less

Submitted 27 April, 2023; originally announced April 2023.

Comments: Proceedings for Gamma 2022, 8 pages, 6 figures

Correcting Imaging Atmospheric Cherenkov Telescope data with atmospheric profiles obtained with an elastic light detecting and ranging system

Authors: Felix Schmuckermaier, Markus Gaug, Christian Fruck, Abelardo Moralejo, Alexander Hahn, Dijana Dominis Prester, Daniela Dorner, Lluís Font, Saša Mićanović, Razmik Mirzoyan, David Paneque, Lovro Pavletić, Julian Sitarek, Martin Will

Abstract: Context. We are operating an elastic LIDAR for the monitoring of atmospheric conditions during regular observations of the MAGIC Telescopes. Aims. We present and evaluate methods to convert aerosol extinction profiles, obtained with the LIDAR, into corrections of the reconstructed gamma-ray event energy and Instrument Response Functions of Imaging Atmospheric Cherenkov Telescopes. Methods. We asse… ▽ More Context. We are operating an elastic LIDAR for the monitoring of atmospheric conditions during regular observations of the MAGIC Telescopes. Aims. We present and evaluate methods to convert aerosol extinction profiles, obtained with the LIDAR, into corrections of the reconstructed gamma-ray event energy and Instrument Response Functions of Imaging Atmospheric Cherenkov Telescopes. Methods. We assess the performance of these correction schemes with almost seven years of Crab Nebula data taken by the MAGIC Telescopes under various zenith angles and different aerosol extinction scenarios of Cherenkov light. Results. The methods enable the reconstruction of data taken under non-optimal atmospheric conditions with aerosol transmissions down to around 0.65 with systematic uncertainties comparable to those for data taken under optimal conditions. For the first time, the correction of data affected by clouds has been included in the assessment. The data can also be corrected when the transmission is lower than 0.65, but the results are less accurate and suffer from larger systematics. △ Less

Submitted 28 February, 2023; v1 submitted 23 February, 2023; originally announced February 2023.

Comments: 25 pages, 20 figures, 4 tables, accepted for publication in Astronomy & Astrophysics

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

Search for Gamma-ray Spectral Lines from Dark Matter Annihilation up to 100 TeV towards the Galactic Center with MAGIC

Authors: MAGIC Collaboration, 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 , et al. (188 additional authors not shown)

Abstract: Line-like features in TeV $γ$-rays constitute a ''smoking gun'' for TeV-scale particle dark matter and new physics. Probing the Galactic Center region with ground-based Cherenkov telescopes enables the search for TeV spectral features in immediate association with a dense dark matter reservoir at a sensitivity out of reach for satellite $γ$-ray detectors, and direct detection and collider experime… ▽ More Line-like features in TeV $γ$-rays constitute a ''smoking gun'' for TeV-scale particle dark matter and new physics. Probing the Galactic Center region with ground-based Cherenkov telescopes enables the search for TeV spectral features in immediate association with a dense dark matter reservoir at a sensitivity out of reach for satellite $γ$-ray detectors, and direct detection and collider experiments. We report on 223 hours of observations of the Galactic Center region with the MAGIC stereoscopic telescope system reaching $γ$-ray energies up to 100 TeV. We improved the sensitivity to spectral lines at high energies using large-zenith-angle observations and a novel background modeling method within a maximum-likelihood analysis in the energy domain. No line-like spectral feature is found in our analysis. Therefore, we constrain the cross section for dark matter annihilation into two photons to $\langle σv \rangle \lesssim 5 \times 10^{-28}\,\mathrm{cm^3\,s^{-1}}$ at 1 TeV and $\langle σv \rangle \lesssim 1 \times 10^{-25}\,\mathrm{cm^3\,s^{-1}}$ at 100 TeV, achieving the best limits to date for a dark matter mass above 20 TeV and a cuspy dark matter profile at the Galactic Center. Finally, we use the derived limits for both cuspy and cored dark matter profiles to constrain supersymmetric wino models. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: Accepted in Phys. Rev. Lett., includes Supplemental materials. 9+4 pages, 3+3 figures, 3 tables. Corresponding authors: T. Inada, D. Kerszberg, M. Hütten

Report number: KEK-TH-2487, KEK-Cosmo-0307

MAGIC observations provide compelling evidence of the hadronic multi-TeV emission from the putative PeVatron SNR G106.3+2.7

Authors: MAGIC Collaboration, H. Abe, S. Abe, V. A. Acciari, I. Agudo, 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 , et al. (192 additional authors not shown)

Abstract: The SNR G106.3+2.7, detected at 1--100 TeV energies by different $γ$-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape which can be divided into a head and a tail region with different physical conditions. However, it is not identified in which region the 100 TeV emission is produced due to the limited position accuracy and/or angular resolution of exi… ▽ More The SNR G106.3+2.7, detected at 1--100 TeV energies by different $γ$-ray facilities, is one of the most promising PeVatron candidates. This SNR has a cometary shape which can be divided into a head and a tail region with different physical conditions. However, it is not identified in which region the 100 TeV emission is produced due to the limited position accuracy and/or angular resolution of existing observational data. Additionally, it remains unclear whether the origin of the $γ$-ray emission is leptonic or hadronic. With the better angular resolution provided by these new MAGIC data compared to earlier $γ$-ray datasets, we aim to reveal the acceleration site of PeV particles and the emission mechanism by resolving the SNR G106.3+2.7 with 0.1$^\circ$ resolution at TeV energies. We detected extended $γ$-ray emission spatially coincident with the radio continuum emission at the head and tail of SNR G106.3+2.7. The fact that we detected a significant $γ$-ray emission with energies above 6.0 TeV from the tail region only suggests that the emissions above 10 TeV, detected with air shower experiments (Milagro, HAWC, Tibet AS$γ$ and LHAASO), are emitted only from the SNR tail. Under this assumption, the multi-wavelength spectrum of the head region can be explained with either hadronic or leptonic models, while the leptonic model for the tail region is in contradiction with the emission above 10 TeV and X-rays. In contrast, the hadronic model could reproduce the observed spectrum at the tail by assuming a proton spectrum with a cutoff energy of $\sim 1$ PeV for the tail region. Such a high energy emission in this middle-aged SNR (4--10 kyr) can be explained by considering the scenario that protons escaping from the SNR in the past interact with surrounding dense gases at present. △ Less

Submitted 28 November, 2022; originally announced November 2022.

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

Journal ref: A&A 671, A12 (2023)

Long-term multi-wavelength study of 1ES 0647+250

Authors: MAGIC Collaboration, 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. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch , et al. (195 additional authors not shown)

Abstract: The BL Lac object 1ES 0647+250 is one of the few distant $γ$-ray emitting blazars detected at very high energies (VHE, $\gtrsim$100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during its low activity in the years 2009-2011, as well as during three flaring activities in the years 2014, 2019 and 2020, with the highest VHE flux in the latter epoch. An extensive multi-in… ▽ More The BL Lac object 1ES 0647+250 is one of the few distant $γ$-ray emitting blazars detected at very high energies (VHE, $\gtrsim$100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during its low activity in the years 2009-2011, as well as during three flaring activities in the years 2014, 2019 and 2020, with the highest VHE flux in the latter epoch. An extensive multi-instrument data set was collected within several coordinated observing campaigns throughout these years. We aim to characterise the long-term multi-band flux variability of 1ES 0647+250, as well as its broadband spectral energy distribution (SED) during four distinct activity states selected in four different epochs, in order to constrain the physical parameters of the blazar emission region under certain assumptions. We evaluate the variability and correlation of the emission in the different energy bands with the fractional variability and the Z-transformed Discrete Correlation Function, as well as its spectral evolution in X-rays and $γ$ rays. Owing to the controversy in the redshift measurements of 1ES 0647+250 reported in the literature, we also estimate its distance in an indirect manner through the comparison of the GeV and TeV spectra from simultaneous observations with Fermi-LAT and MAGIC during the strongest flaring activity detected to date. Moreover, we interpret the SEDs from the four distinct activity states within the framework of one-component and two-component leptonic models, proposing specific scenarios that are able to reproduce the available multi-instrument data. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 20 pages, 7 figures. Accepted in A&A. Corresponding authors: Jorge Otero-Santos; Daniel Morcuende; Vandad Fallah Ramazani; Daniela Dorner; David Paneque (mailto: contact.magic@mpp.mpg.de)

Journal ref: A&A 670, A49 (2023)

A lower bound on intergalactic magnetic fields from time variability of 1ES 0229+200 from MAGIC and Fermi/LAT observations

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

Abstract: Extended and delayed emission around distant TeV sources induced by the effects of propagation of gamma rays through the intergalactic medium can be used for the measurement of the intergalactic magnetic field (IGMF). We search for delayed GeV emission from the hard-spectrum TeV blazar 1ES 0229+200 with the goal to detect or constrain the IGMF-dependent secondary flux generated during the propagat… ▽ More Extended and delayed emission around distant TeV sources induced by the effects of propagation of gamma rays through the intergalactic medium can be used for the measurement of the intergalactic magnetic field (IGMF). We search for delayed GeV emission from the hard-spectrum TeV blazar 1ES 0229+200 with the goal to detect or constrain the IGMF-dependent secondary flux generated during the propagation of TeV gamma rays through the intergalactic medium. We analyze the most recent MAGIC observations over a 5 year time span and complement them with historic data of the H.E.S.S. and VERITAS telescopes along with a 12-year long exposure of the Fermi/LAT telescope. We use them to trace source evolution in the GeV-TeV band over one-and-a-half decade in time. We use Monte Carlo simulations to predict the delayed secondary gamma-ray flux, modulated by the source variability, as revealed by TeV-band observations. We then compare these predictions for various assumed IGMF strengths to all available measurements of the gamma-ray flux evolution. We find that the source flux in the energy range above 200 GeV experiences variations around its average on the 14 years time span of observations. No evidence for the flux variability is found in 1-100 GeV energy range accessible to Fermi/LAT. Non-detection of variability due to delayed emission from electromagnetic cascade developing in the intergalactic medium imposes a lower bound of B>1.8e-17 G for long correlation length IGMF and B>1e-14 G for an IGMF of the cosmological origin. Though weaker than the one previously derived from the analysis of Fermi/LAT data, this bound is more robust, being based on a conservative intrinsic source spectrum estimate and accounting for the details of source variability in the TeV energy band. We discuss implications of this bound for cosmological magnetic fields which might explain the baryon asymmetry of the Universe. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 10 pages, 5 figures, accepted to A&A. Corresponding authors: Ievgen Vovk, Paolo Da Vela (mailto:contact.magic@mpp.mpg.de) and Andrii Neronov (mailto:Andrii.Neronov@unige.ch)

Journal ref: A&A 670, A145 (2023)

Multi-messenger characterization of Mrk 501 during historically low X-ray and $γ$-ray activity

Authors: MAGIC collaboration, H. Abe, S. Abe, V. A. Acciari, I. Agudo, 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 , et al. (300 additional authors not shown)

Abstract: We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the… ▽ More We study the broadband emission of Mrk 501 using multi-wavelength observations from 2017 to 2020 performed with a multitude of instruments, involving, among others, MAGIC, Fermi-LAT, NuSTAR, Swift, GASP-WEBT, and OVRO. Mrk 501 showed an extremely low broadband activity, which may help to unravel its baseline emission. Nonetheless, significant flux variations are detected at all wavebands, with the highest occurring at X-rays and very-high-energy (VHE) $γ$-rays. A significant correlation ($>$3$σ$) between X-rays and VHE $γ$-rays is measured, supporting leptonic scenarios to explain the variable parts of the emission, also during low activity. This is further supported when we extend our data from 2008 to 2020, and identify, for the first time, significant correlations between Swift-XRT and Fermi-LAT. We additionally find correlations between high-energy $γ$-rays and radio, with the radio lagging by more than 100 days, placing the $γ$-ray emission zone upstream of the radio-bright regions in the jet. Furthermore, Mrk 501 showed a historically low activity in X-rays and VHE $γ$-rays from mid-2017 to mid-2019 with a stable VHE flux ($>$0.2 TeV) of 5% the emission of the Crab Nebula. The broadband spectral energy distribution (SED) of this 2-year-long low-state, the potential baseline emission of Mrk 501, can be characterized with one-zone leptonic models, and with (lepto)-hadronic models fulfilling neutrino flux constraints from IceCube. We explore the time evolution of the SED towards the low-state, revealing that the stable baseline emission may be ascribed to a standing shock, and the variable emission to an additional expanding or traveling shock. △ Less

Submitted 5 March, 2023; v1 submitted 5 October, 2022; originally announced October 2022.

Comments: 55 pages, 30 figures, 14 tables, accepted by APJS. Corresponding authors are L. Heckmann, D. Paneque, S. Gasparyan, M. Cerruti, and N. Sahakyan

Journal ref: ApJS 266 37 (2023)

Multi-wavelength study of the galactic PeVatron candidate LHAASO J2108+5157

Authors: S. Abe, A. Aguasca-Cabot, I. Agudo, N. Alvarez Crespo, L. A. Antonelli, C. Aramo, A. Arbet-Engels, M. Artero, K. Asano, P. Aubert, A. Baktash, A. Bamba, A. Baquero Larriva, L. Baroncelli, U. Barres de Almeida, J. A. Barrio, I. Batkovic, J. Baxter, J. Becerra González, E. Bernardini, M. I. Bernardos, J. Bernete Medrano, A. Berti, P. Bhattacharjee, N. Biederbeck , et al. (245 additional authors not shown)

Abstract: LHAASO J2108+5157 is one of the few known unidentified Ultra-High-Energy (UHE) gamma-ray sources with no Very-High-Energy (VHE) counterpart, recently discovered by the LHAASO collaboration. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good quality data. In… ▽ More LHAASO J2108+5157 is one of the few known unidentified Ultra-High-Energy (UHE) gamma-ray sources with no Very-High-Energy (VHE) counterpart, recently discovered by the LHAASO collaboration. We observed LHAASO J2108+5157 in the X-ray band with XMM-Newton in 2021 for a total of 3.8 hours and at TeV energies with the Large-Sized Telescope prototype (LST-1), yielding 49 hours of good quality data. In addition, we analyzed 12 years of Fermi-LAT data, to better constrain emission of its High-Energy (HE) counterpart 4FGL J2108.0+5155. We found an excess (3.7 sigma) in the LST-1 data at energies E > 3 TeV. Further analysis in the whole LST-1 energy range assuming a point-like source, resulted in a hint (2.2 sigma) of hard emission which can be described with a single power law with photon index Gamma = 1.6 +- 0.2 between 0.3 - 100 TeV. We did not find any significant extended emission which could be related to a Supernova Remnant (SNR) or Pulsar Wind Nebula (PWN) in the XMM-Newton data, which puts strong constraints on possible synchrotron emission of relativistic electrons. The LST-1 and LHAASO observations can be explained as inverse Compton-dominated leptonic emission of relativistic electrons with a cutoff energy of $100^{+70}_{-30}$ TeV. The low magnetic field in the source imposed by the X-ray upper limits on synchrotron emission is compatible with a hypothesis of a PWN or a TeV halo. The lack of a pulsar in the neighborhood of the UHE source is a challenge to the PWN/TeV-halo scenario. The UHE gamma rays can also be explained as $π^0$ decay-dominated hadronic emission due to interaction of relativistic protons with one of the two known molecular clouds in the direction of the source. The hard spectrum in the LST-1 band is compatible with protons escaping a shock around a middle-aged SNR because of their high low-energy cut-off. △ Less

Submitted 16 March, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 26 pages, 9 figures, 5 tables; Astronomy and Astrophysics (accepted on 6 March 2023)

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

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

Rapid X-ray Variability in Mkn 421 during a Multiwavelength Campaign

Authors: Alex G. Markowitz, Krzysztof Nalewajko, Gopal Bhatta, Gulab C. Dewangan, Sunil Chandra, Daniela Dorner, Bernd Schleicher, Urszula Pajdosz-Smierciak, Lukasz Stawarz, Staszek Zola, Michal Ostrowski, Daniele Carosati, Saikruba Krishnan, Rumen Bachev, Erika Benitez, Kosmas Gazeas, David Hiriart, Shao-Ming Hu, Valeri Larionov, Alessandro Marchini, Katsura Matsumoto, A. A. Nikiforova, Tapio Pursimo, Claudia M. Raiteri, Daniel E. Reichart , et al. (25 additional authors not shown)

Abstract: The study of short-term variability properties in AGN jets has the potential to shed light on their particle acceleration and emission mechanisms. We report results from a four-day coordinated multi-wavelength campaign on the highly-peaked blazar (HBL) Mkn 421 in 2019 January. We obtained X-ray data from AstroSAT, BVRI photometry with the Whole Earth Blazar Telescope (WEBT), and TeV data from FACT… ▽ More The study of short-term variability properties in AGN jets has the potential to shed light on their particle acceleration and emission mechanisms. We report results from a four-day coordinated multi-wavelength campaign on the highly-peaked blazar (HBL) Mkn 421 in 2019 January. We obtained X-ray data from AstroSAT, BVRI photometry with the Whole Earth Blazar Telescope (WEBT), and TeV data from FACT to explore short-term multi-wavelength variability in this HBL. The X-ray continuum is rapidly variable on time-scales of tens of ks. Fractional variability amplitude increases with energy across the synchrotron hump, consistent with previous studies; we interpret this observation in the context of a model with multiple cells whose emission spectra contain cutoffs that follow a power-law distribution. We also performed time-averaged and time-resolved (time-scales of 6 ks) spectral fits; a broken power-law model fits all spectra well; time-resolved spectral fitting reveals the usual hardening when brightening behaviour. Intra-X-ray cross correlations yield evidence for the 0.6-0.8 keV band to likely lead the other bands by an average of 4.6 +- 2.6 ks, but only during the first half of the observation. The source displayed minimal night-to-night variability at all wavebands thus precluding significant interband correlations during our campaign. The broadband SED is modeled well with a standard one-zone leptonic model, yielding jet parameters consistent with those obtained from previous SEDs of this source. △ Less

Submitted 8 June, 2022; originally announced June 2022.

Comments: This article has been accepted for publication in The Monthly Notices of the Royal Astronomical Society (2022), and is published in Volume 513, Issue 2, pp.1662-1679. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved. 19 pages; 12 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, 2022, Volume 513, Issue 2, pp.1662-1679

Characterizing the aerosol atmosphere above the Observatorio del Roque de los Muchachos by analyzing seven years of data taken with an GaAsP HPD-readout, absolutely calibrated elastic LIDAR

Authors: Christian Fruck, Markus Gaug, Alexander Hahn, Victor Acciari, Jürgen Besenrieder, Dijana Dominis Prester, Daniela Dorner, David Fink, Lluís Font, Saša Mićanović, Razmik Mirzoyan, Dominik Müller, Lovro Pavletić, Felix Schmuckermaier, Martin Will

Abstract: We present a new elastic LIDAR concept, based on a bi-axially mounted Nd:YAG laser and a telescope with HPD readout, combined with fast FADC signal digitization and offline pulse analysis. The LIDAR return signals have been extensively quality checked and absolutely calibrated. We analyze seven years of quasi-continuous LIDAR data taken during those nights when the MAGIC telescopes were operating.… ▽ More We present a new elastic LIDAR concept, based on a bi-axially mounted Nd:YAG laser and a telescope with HPD readout, combined with fast FADC signal digitization and offline pulse analysis. The LIDAR return signals have been extensively quality checked and absolutely calibrated. We analyze seven years of quasi-continuous LIDAR data taken during those nights when the MAGIC telescopes were operating. Characterization of the nocturnal ground layer yields zenith and azimuth angle dependent aerosol extinction scale heights for clear nights. We derive aerosol transmission statistics for light emitted from various altitudes throughout the year and separated by seasons. We find further seasonal dependencies of cloud base and top altitudes, but none for the LIDAR ratios of clouds. Finally, the night sky background light is characterized using the LIDAR photon backgrounds. abstract.txt △ Less

Submitted 19 February, 2022; originally announced February 2022.

Comments: 40 pages, submitted to MNRAS

ACM Class: J.2.3; J.2.5

Journal ref: Monthly Notices of the Royal Academic Society, Volume 515, Issue 3, September 2022, Pages 4520-4550

Proton acceleration in thermonuclear nova explosions revealed by gamma rays

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, H. Bökenkamp , et al. (186 additional authors not shown)

Abstract: Classical novae are cataclysmic binary star systems in which the matter of a companion star is accreted on a white dwarf (WD). Accumulation of hydrogen in a layer eventually causes a thermonuclear explosion on the surface of the WD, brightening the WD to ~10^5 solar luminosities and triggering ejection of the accumulated matter.They provide extreme conditions required to accelerate particles, elec… ▽ More Classical novae are cataclysmic binary star systems in which the matter of a companion star is accreted on a white dwarf (WD). Accumulation of hydrogen in a layer eventually causes a thermonuclear explosion on the surface of the WD, brightening the WD to ~10^5 solar luminosities and triggering ejection of the accumulated matter.They provide extreme conditions required to accelerate particles, electrons or protons, to high energies. Here we present the detection of gamma rays by the MAGIC telescopes from the 2021 outburst of RS Ophiuchi (RS Oph), a recurrent nova with a red giant (RG) companion, that allowed us, for the first time, to accurately characterize the emission from a nova in the 60 GeV to 250 GeV energy range. The theoretical interpretation of the combined Fermi-LAT and MAGIC data suggests that protons are accelerated to hundreds of GeV in the nova shock. Such protons should create bubbles of enhanced Cosmic Ray density, on the order of 10 pc, from the recurrent novae. △ Less

Submitted 10 November, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: THIS VERSION IS AFTER THE JOURNAL REVIEW

Journal ref: Nature Astronomy, Volume 6, p. 689-697 (2022)

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

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

Abstract: The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the Ic… ▽ More The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $γ$ rays (Fermi/LAT) and the very-high-energy (VHE) $γ$ rays (MAGIC). In $γ$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $γ$ rays. According to the model presented here, the December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments. △ Less

Submitted 1 May, 2022; v1 submitted 5 February, 2022; originally announced February 2022.

Comments: 18 pages, 6 figures; V2 to match published version, typo fixed

Combined searches for dark matter in dwarf spheroidal galaxies observed with the MAGIC telescopes, including new data from Coma Berenices and Draco

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, H. Bökenkamp , et al. (169 additional authors not shown)

Abstract: Milky Way dwarf spheroidal galaxies (dSphs) are among the best candidates to search for signals of dark matter annihilation with Imaging Atmospheric Cherenkov Telescopes, given their high mass-to-light ratios and the fact that they are free of astrophysical gamma-ray emitting sources. Since 2011, MAGIC has performed a multi-year observation program in search for Weakly Interacting Massive Particle… ▽ More Milky Way dwarf spheroidal galaxies (dSphs) are among the best candidates to search for signals of dark matter annihilation with Imaging Atmospheric Cherenkov Telescopes, given their high mass-to-light ratios and the fact that they are free of astrophysical gamma-ray emitting sources. Since 2011, MAGIC has performed a multi-year observation program in search for Weakly Interacting Massive Particles (WIMPs) in dSphs. Results on the observations of Segue 1 and Ursa Major II dSphs have already been published and include some of the most stringent upper limits (ULs) on the velocity-averaged cross-section $\langle σ_{\mathrm{ann}} v \rangle$ of WIMP annihilation from observations of dSphs. In this work, we report on the analyses of 52.1 h of data of Draco dSph and 49.5 h of Coma Berenices dSph observed with the MAGIC telescopes in 2018 and in 2019 respectively. No hint of a signal has been detected from either of these targets and new constraints on the $\langle σ_{\mathrm{ann}} v \rangle$ of WIMP candidates have been derived. In order to improve the sensitivity of the search and reduce the effect of the systematic uncertainties due to the $J$-factor estimates, we have combined the data of all dSphs observed with the MAGIC telescopes. Using 354.3 h of dSphs good quality data, 95 % CL ULs on $\langle σ_{\mathrm{ann}} v \rangle$ have been obtained for 9 annihilation channels. For most of the channels, these results reach values of the order of $10^{-24} $cm$^3$/s at ${\sim}1$ TeV and are the most stringent limits obtained with the MAGIC telescopes so far. △ Less

Submitted 29 November, 2021; originally announced November 2021.

Comments: 17 pages, 5 figures

Journal ref: Physics of the Dark Universe, Volume 35, March 2022, 100912

Multiwavelength study of the gravitationally lensed blazar QSO B0218+357 between 2016 and 2020

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli , et al. (186 additional authors not shown)

Abstract: We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We const… ▽ More We report multiwavelength observations of the gravitationally lensed blazar QSO B0218+357 in 2016-2020. Optical, X-ray and GeV flares were detected. The contemporaneous MAGIC observations do not show significant very-high-energy (VHE, >= 100 GeV) gamma-ray emission. The lack of enhancement in radio emission measured by OVRO indicates the multi-zone nature of the emission from this object. We constrain the VHE duty cycle of the source to be < 16 2014-like flares per year (95% confidence). For the first time for this source, a broadband low-state SED is constructed with a deep exposure up to the VHE range. A flux upper limit on the low-state VHE gamma-ray emission of an order of magnitude below that of the 2014 flare is determined. The X-ray data are used to fit the column density of (8.10 +- 0.93 stat ) x 10^21 cm^-2 of the dust in the lensing galaxy. VLBI observations show a clear radio core and jet components in both lensed images, yet no significant movement of the components is seen. The radio measurements are used to model the source-lens-observer geometry and determine the magnifications and time delays for both components. The quiescent emission is modeled with the high-energy bump explained as a combination of synchrotron-self-Compton and external Compton emission from a region located outside of the broad line region. The bulk of the low-energy emission is explained as originating from a tens-of-parsecs scale jet. △ Less

Submitted 25 November, 2021; originally announced November 2021.

Comments: 21 pages, 12 figures, 5 tables, accepted for publication in MNRAS

Statistical properties of flux variations in blazar light curves at GeV and TeV energies

Authors: Sarah M. Wagner, Paul R. Burd, Daniela Dorner, Karl Mannheim, Sara Buson, Andrea Gokus, Greg Madejski, Jeffrey D. Scargle, Axel Arbet-Engels, Dominik Baack, Matteo Balbo, Adrian Biland, Thomas Bretz, Jens Buss, Laura Eisenberger, Dominik Elsaesser, Dorothee Hildebrand, Roman Iotov, Adelina Kalenski, Dominik Neise, Maximilian Noethe, Aleksander Paravac, Wolfgang Rhode, Bernd Schleicher, Vitalii Sliusar , et al. (1 additional authors not shown)

Abstract: Despite numerous detections of individual flares, the physical origin of the rapid variability observed from blazars remains uncertain. Using Bayesian blocks and the Eisenstein-Hut HOP algorithm, we characterize flux variations of high significance in the $γ$-ray light curves of two samples of blazars. Daily binned long-term light curves of TeV-bright blazars observed with the First G-APD Cherenko… ▽ More Despite numerous detections of individual flares, the physical origin of the rapid variability observed from blazars remains uncertain. Using Bayesian blocks and the Eisenstein-Hut HOP algorithm, we characterize flux variations of high significance in the $γ$-ray light curves of two samples of blazars. Daily binned long-term light curves of TeV-bright blazars observed with the First G-APD Cherenkov Telescope (FACT) are compared to those of GeV-bright blazars observed with the Large Area Telescope on board the $Fermi$ Gamma-ray Space Telescope ($Fermi$-LAT). We find no evidence for systematic asymmetry of the flux variations based on the derived rise and decay time scales. Additionally, we show that the daily-binned blazar light curves can be described by an exponential stochastic Ornstein-Uhlenbeck (OU) process with parameters depending on energy. Our analysis suggests that the flux variability in both samples is a superposition of faster fluctuations. This is, for instance, challenging to explain by shock-acceleration but expected for magnetic reconnection. △ Less

Submitted 8 November, 2021; v1 submitted 27 October, 2021; originally announced October 2021.

Comments: for associated HOP algorithm, see https://github.com/swagner-astro/lightcurves and for associated OU algorithm, see https://github.com/PRBurd/astro-wue

Observation of the gamma-ray binary HESS J0632+057 with the H.E.S.S., MAGIC, and VERITAS telescopes

Authors: C. B. Adams, W. Benbow, A. Brill, J. H. Buckley, M. Capasso, A. J. Chromey, M. Errando, A. Falcone, K. A. Farrell, Q. Feng, J P. Finley, G. Foote, L. Fortson, A. Furniss, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, D. Hanna, T. Hassan, O. Hervet, J. Holder, B. Hona, T. B. Humensky, W. Jin , et al. (387 additional authors not shown)

Abstract: The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these obs… ▽ More The results of gamma-ray observations of the binary system HESS J0632+057 collected during 450 hours over 15 years, between 2004 and 2019, are presented. Data taken with the atmospheric Cherenkov telescopes H.E.S.S., MAGIC, and VERITAS at energies above 350 GeV were used together with observations at X-ray energies obtained with Swift-XRT, Chandra, XMM-Newton, NuSTAR, and Suzaku. Some of these observations were accompanied by measurements of the Hα emission line. A significant detection of the modulation of the VHE gamma-ray fluxes with a period of 316.7+-4.4 days is reported, consistent with the period of 317.3+-0.7 days obtained with a refined analysis of X-ray data. The analysis of data of four orbital cycles with dense observational coverage reveals short timescale variability, with flux-decay timescales of less than 20 days at very high energies. Flux variations observed over the time scale of several years indicate orbit-to-orbit variability. The analysis confirms the previously reported correlation of X-ray and gamma-ray emission from the system at very high significance, but can not find any correlation of optical Hα parameters with X-ray or gamma-ray energy fluxes in simultaneous observations. The key finding is that the emission of HESS J0632+057 in the X-ray and gamma-ray energy bands is highly variable on different time scales. The ratio of gamma-ray to X-ray flux shows the equality or even dominance of the gamma-ray energy range. This wealth of new data is interpreted taking into account the insufficient knowledge of the ephemeris of the system, and discussed in the context of results reported on other gamma-ray binary systems. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Comments: Accepted for publication by The Astrophysical Journal

Searching for VHE gamma-ray emission associated with IceCube neutrino alerts using FACT, H.E.S.S., MAGIC, and VERITAS

Authors: Konstancja Satalecka, Elisa Bernardini, Daniela Dorner, Gašper Kukec Mezek, Weidong Jin

Abstract: The realtime follow-up of neutrino events is a promising approach to search for astrophysical neutrino sources. It has so far provided compelling evidence for a neutrino point source: the flaring gamma-ray blazar TXS 0506+056 observed in coincidence with the high-energy neutrino IceCube-170922A detected by IceCube. The detection of very-high-energy gamma rays (VHE,… ▽ More The realtime follow-up of neutrino events is a promising approach to search for astrophysical neutrino sources. It has so far provided compelling evidence for a neutrino point source: the flaring gamma-ray blazar TXS 0506+056 observed in coincidence with the high-energy neutrino IceCube-170922A detected by IceCube. The detection of very-high-energy gamma rays (VHE, $\mathrm{E} > 100\,\mathrm{GeV}$) from this source helped establish the coincidence and constrained the modeling of the blazar emission at the time of the IceCube event. The four major imaging atmospheric Cherenkov telescope arrays (IACTs) - FACT, H.E.S.S., MAGIC, and VERITAS - operate an active follow-up program of target-of-opportunity observations of neutrino alerts sent by IceCube. This program has two main components. One are the observations of known gamma-ray sources around which a cluster of candidate neutrino events has been identified by IceCube (Gamma-ray Follow-Up, GFU). Second one is the follow-up of single high-energy neutrino candidate events of potential astrophysical origin such as IceCube-170922A. GFU has been recently upgraded by IceCube in collaboration with the IACT groups. We present here recent results from the IACT follow-up programs of IceCube neutrino alerts and a description of the upgraded IceCube GFU system. △ Less

Submitted 9 September, 2021; originally announced September 2021.

Comments: Presented at the 37th International Cosmic Ray Conference (ICRC 2021). See arXiv:2107.06966 for all IceCube contributions. See arXiv:2108.05257 for all H.E.S.S. contributions

Report number: PoS-ICRC2021-960

Journal ref: PoS(ICRC2021)960

Long-term multi-band photometric monitoring of Mrk 501

Authors: Axel Arbet-Engels, Dominik Baack, Matteo Balbo, Adrian Biland, Thomas Bretz, Jens Buss, Daniela Dorner, Laura Eisenberger, Dominik Elsaesser, Dorothee Hildebrand, Roman Iotov, Adelina Kalenski, Karl Mannheim, Alison Mitchell, Dominik Neise, Maximilian Noethe, Aleksander Paravac, Wolfgang Rhode, Bernd Schleicher, Vitalii Sliusar, Roland Walter

Abstract: Radio-to-TeV observations of the bright nearby (z=0.034) blazar Markarian 501 (Mrk 501), performed from December 2012 to April 2018, are used to study the emission mechanisms in its relativistic jet. We examined the multi-wavelength variability and the correlations of the light curves obtained by eight different instruments, including the First G-APD Cherenkov Telescope (FACT), observing Mrk 501 i… ▽ More Radio-to-TeV observations of the bright nearby (z=0.034) blazar Markarian 501 (Mrk 501), performed from December 2012 to April 2018, are used to study the emission mechanisms in its relativistic jet. We examined the multi-wavelength variability and the correlations of the light curves obtained by eight different instruments, including the First G-APD Cherenkov Telescope (FACT), observing Mrk 501 in very high-energy (VHE) gamma-rays at TeV energies. We identified individual TeV and X-ray flares and found a sub-day lag between variability in these two bands. Simultaneous TeV and X-ray variations with almost zero lag are consistent with synchrotron self-Compton (SSC) emission, where TeV photons are produced through inverse Compton scattering. The characteristic time interval of 5-25 days between TeV flares is consistent with them being driven by Lense-Thirring precession. △ Less

Submitted 30 November, 2021; v1 submitted 7 September, 2021; originally announced September 2021.

Comments: 10 pages, 11 figures, Astronomy and astrophysics pre-print. arXiv admin note: text overlap with arXiv:2101.10651

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

Search for Very High-Energy Emission from the millisecond pulsar PSR J0218+4232

Authors: V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Ž. Bošnjak , et al. (176 additional authors not shown)

Abstract: PSR J0218+4232 is one of the most energetic millisecond pulsars known and has long been considered as one of the best candidates for very high-energy (VHE; >100 GeV) gamma-ray emission. Using 11.5 years of Fermi Large Area Telescope (LAT) data between 100 MeV and 870 GeV, and ~90 hours of MAGIC observations in the 20 GeV to 20 TeV range, we have searched for the highest energy gamma-ray emission f… ▽ More PSR J0218+4232 is one of the most energetic millisecond pulsars known and has long been considered as one of the best candidates for very high-energy (VHE; >100 GeV) gamma-ray emission. Using 11.5 years of Fermi Large Area Telescope (LAT) data between 100 MeV and 870 GeV, and ~90 hours of MAGIC observations in the 20 GeV to 20 TeV range, we have searched for the highest energy gamma-ray emission from PSR J0218+4232. Based on the analysis of the LAT data, we find evidence for pulsed emission above 25 GeV, but see no evidence for emission above 100 GeV (VHE) with MAGIC. We present the results of searches for gamma-ray emission, along with theoretical modeling, to interpret the lack of VHE emission. We conclude that, based on the experimental observations and theoretical modeling, it will remain extremely challenging to detect VHE emission from PSR J0218+4232 with the current generation of Imaging Atmospheric Cherenkov Telescopes (IACTs), and maybe even with future ones, such as the Cherenkov Telescope Array (CTA). △ Less

Submitted 25 August, 2021; originally announced August 2021.

Comments: 19 pages, 9 figures, Accepted to ApJ

Relentless multi-wavelength variability of Mrk 421 and Mrk 501

Authors: Vitalii Sliusar, Axel Arbet-Engels, Dominik Baack, Matteo Balbo, Adrian Biland, Thomas Bretz, Jens Buss, Daniela Dorner, Laura Eisenberger, Dominik Elsaesser, Dorothee Hildebrand, Roman Iotov, Adelina Kalenski, Karl Mannheim, Alison Mitchell, Dominik Neise, Maximilian Noethe, Aleksander Paravac, Wolfgang Rhode, Bernd Schleicher, Roland Walter

Abstract: Mrk 421 and Mrk 501 are two close, bright and well-studied high-synchrotron-peaked blazars, which feature bright and persistent GeV and TeV emission. We use the longest and densest dataset of unbiased observations of these two sources, obtained at TeV and GeV energies during five years with FACT and Fermi-LAT. To characterize the variability and derive constraints on the emission mechanism, we aug… ▽ More Mrk 421 and Mrk 501 are two close, bright and well-studied high-synchrotron-peaked blazars, which feature bright and persistent GeV and TeV emission. We use the longest and densest dataset of unbiased observations of these two sources, obtained at TeV and GeV energies during five years with FACT and Fermi-LAT. To characterize the variability and derive constraints on the emission mechanism, we augment the dataset with contemporaneous multi-wavelength observations from radio to X-rays. We correlate the light curves, identify individual flares in TeV energies and X-rays, and look for inter-band connections, which are expected from the shock propagations within the jet. For Mrk 421, we find that the X-rays and TeV energies are well correlated with close to zero lag, supporting the SSC emission scenario. The timing between the TeV, X-ray flares in Mrk 421 is consistent with periods expected in the case of Lense-Thirring precession of the accretion disc. The variability of Mrk 501 on long-term periods is also consistent with SSC, with a sub-day lag between X-rays and TeV energies. Fractional variability for both blazars shows a two bump structure with the highest variability in the X-ray and TeV bands. △ Less

Submitted 18 August, 2021; originally announced August 2021.

Comments: 8 pages, 3 figures

Multi-wavelength study of Mrk 421 during a TeV flare

Authors: A. Gokus, A. Kreikenbohm, K. Leiter, T. Bretz, T. Dauser, D. Dorner, D. Elsaesser, F. Eppel, J. Hessdoerfer, M. Kadler, A. Kraus, M. Kreter, I. Kreykenbohm, M. Langejahn, K. Mannheim, P. Thalhammer, J. Wilms, A. Arbet-Engels, D. Baack, M. Balbo, A. Biland, J. Buss, L. Eisenberger, D. Hildebrand, R. Iotov , et al. (9 additional authors not shown)

Abstract: The blazar Mrk 421 shows frequent, short flares in the TeV energy regime. Due to the fast nature of such episodes, we often fail to obtain sufficient simultaneous information about flux variations in several energy bands. To overcome this lack of multi-wavelength (MWL) coverage, especially for the pre- and post-flare periods, we have set up a monitoring program with the FACT telescope (TeV energie… ▽ More The blazar Mrk 421 shows frequent, short flares in the TeV energy regime. Due to the fast nature of such episodes, we often fail to obtain sufficient simultaneous information about flux variations in several energy bands. To overcome this lack of multi-wavelength (MWL) coverage, especially for the pre- and post-flare periods, we have set up a monitoring program with the FACT telescope (TeV energies) and the Neil Gehrels Swift Observatory (X-rays). On 2019 June 9, Mrk 421 showed a TeV outburst reaching a flux level of more than two times the flux of the Crab Nebula at TeV energies. We acquired simultaneous data in the X-rays with additional observations by XMM-Newton and INTEGRAL. For the first time, we can study a TeV blazar in outburst taking advantage of highly sensitive X-ray data from XMM-Newton and INTEGRAL combined. Our dataset is complemented by pointed radio observations by Effelsberg at GHz frequencies. We present our first results, including the γ-ray and X-ray light curves, a timing analysis of the X-ray data obtained with XMM-Newton , as well as the radio spectra before, during and after the flare. △ Less

Submitted 4 August, 2021; originally announced August 2021.

Comments: 8 pages, 4 figures; Proceedings of the 37th International Cosmic Ray Conference (ICRC2021). 12-23 July, 2021. Berlin, Germany. Online at https://pos.sissa.it/395/

TELAMON: Effelsberg Monitoring of AGN Jets with Very-High-Energy Astroparticle Emissions

Authors: M. Kadler, U. Bach, D. Berge, S. Buson, D. Dorner, P. G. Edwards, F. Eppel, M. Giroletti, A. Gokus, O. Hervet, J. Heßdörfer, S. Koyama, A. Kraus, T. P. Krichbaum, E. Lindfors, K. Mannheim, R. de Menezes, R. Ojha, G. F. Paraschos, E. Pueschel, F. Rösch, E. Ros, B. Schleicher, J. Sinapius, J. Sitarek , et al. (2 additional authors not shown)

Abstract: We introduce the TELAMON program, which is using the Effelsberg 100-m telescope to monitor the radio spectra of active galactic nuclei (AGN) under scrutiny in astroparticle physics, namely TeV blazars and candidate neutrino-associated AGN. Thanks to its large dish aperture and sensitive instrumentation, the Effelsberg telescope can yield radio data superior over other programs in the low flux-dens… ▽ More We introduce the TELAMON program, which is using the Effelsberg 100-m telescope to monitor the radio spectra of active galactic nuclei (AGN) under scrutiny in astroparticle physics, namely TeV blazars and candidate neutrino-associated AGN. Thanks to its large dish aperture and sensitive instrumentation, the Effelsberg telescope can yield radio data superior over other programs in the low flux-density regime down to several 10 mJy. This is a particular strength in the case of TeV-emitting blazars, which are often comparatively faint radio sources of the high-synchrotron peaked type. We perform high-cadence high-frequency observations every 2-4 weeks at multiple frequencies up to 44 GHz. This setup is well suited to trace dynamical processes in the compact parsec-scale jets of blazars related to high-energy flares or neutrino detections. Our sample currently covers about 40 sources and puts its focus on AGN with very-high-energy astroparticle emission, i.e., TeV blazars and neutrino-associated AGN. Here, we introduce the TELAMON program characteristics and present first results obtained since fall 2020. △ Less

Submitted 1 August, 2021; originally announced August 2021.

Comments: 8 pages, 5 figures, Proceedings of the 37th International Cosmic Ray Conference (ICRC2021). 12-23 July, 2021. Berlin, Germany. Online at https://pos.sissa.it/395/

First detection of VHE gamma-ray emission from TXS~1515--273, study of its X-ray variability and spectral energy distribution

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, Ž. Bošnjak , et al. (174 additional authors not shown)

Abstract: We report here on the first multi-wavelength (MWL) campaign on the blazar TXS 1515-273, undertaken in 2019 and extending from radio to very-high-energy gamma rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotronpeaked source. MAGIC observations result… ▽ More We report here on the first multi-wavelength (MWL) campaign on the blazar TXS 1515-273, undertaken in 2019 and extending from radio to very-high-energy gamma rays (VHE). Up until now, this blazar had not been the subject of any detailed MWL observations. It has a rather hard photon index at GeV energies and was considered a candidate extreme high-synchrotronpeaked source. MAGIC observations resulted in the first-time detection of the source in VHE with a statistical significance of 7.6$σ$. The average integral VHE flux of the source is 6 $\pm$ 1% of the Crab nebula flux above 400 GeV. X-ray coverage was provided by Swift-XRT, XMMNewton, and NuSTAR. The long continuous X-ray observations were separated by $\sim$ 9 h, both showing clear hour scale flares. In the XMM-Newton data, both the rise and decay timescales are longer in the soft X-ray than in the hard X-ray band, indicating the presence of a particle cooling regime. The X-ray variability timescales were used to constrain the size of the emission region and the strength of the magnetic field. The data allowed us to determine the synchrotron peak frequency and classify the source as a flaring high, but not extreme, synchrotron peaked object. Considering the constraints and variability patterns from the X-ray data, we model the broad-band spectral energy distribution. We applied a simple one-zone model, which could not reproduce the radio emission and the shape of the optical emission, and a two-component leptonic model with two interacting components, enabling us to reproduce the emission from radio to VHE band. △ Less

Submitted 20 July, 2021; originally announced July 2021.

Investigation of the correlation patterns and the Compton dominance variability of Mrk 421 in 2017

Authors: MAGIC Collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, L. Bellizzi, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, O. Blanch, Ž. Bošnjak, G. Busetto, R. Carosi , et al. (263 additional authors not shown)

Abstract: We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with special focus on the multi-band flux correlations. The dataset has been collected through an extensive multiwavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicin… ▽ More We present a detailed characterisation and theoretical interpretation of the broadband emission of the paradigmatic TeV blazar Mrk 421, with special focus on the multi-band flux correlations. The dataset has been collected through an extensive multiwavelength campaign organised between 2016 December and 2017 June. The instruments involved are MAGIC, FACT, Fermi-LAT, Swift, GASP-WEBT, OVRO, Medicina and Metsähovi. Additionally, four deep exposures (several hours long) with simultaneous MAGIC and NuSTAR observations allowed a precise measurement of the falling segments of the two spectral components. The very-high-energy (VHE; E > 100 GeV) gamma rays and X-rays are positively correlated at zero time lag, but the strength and characteristics of the correlation change substantially across the various energy bands probed. The VHE versus X-ray fluxes follow different patterns, partly due to substantial changes in the Compton dominance during a few days without a simultaneous increase in the X-ray flux (i.e. orphan gamma-ray activity). Studying the broadband spectral energy distribution (SED) during the days including NuSTAR observations, we show that these changes can be explained within a one-zone leptonic model with a blob that increases its size over time. Our multi-band correlation study also hints at an anti-correlation between UV/optical and X-ray at a significance higher than 3 sigmas. A VHE flare observed on 2017 February 4 shows gamma-ray variability on multi-hour timescales, with a factor 10 increase in the TeV flux but only a moderate increase in the keV flux. The related broadband SED is better described by a two-zone leptonic scenario rather than by a one-zone scenario. We find that the flare can be produced by the appearance of a compact second blob populated by high energetic electrons spanning a narrow range of Lorentz factors. △ Less

Submitted 10 June, 2021; originally announced June 2021.

Comments: 38 pages, 23 figures, accepted for publication in Astronomy and Astrophysics. Corresponding authors: Axel Arbet-Engels, David Paneque

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