Showing 1–50 of 222 results for author: González, J B

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.

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 23 July, 2024; v1 submitted 7 March, 2024; originally announced March 2024.

Comments: 44 pages JCAP style (excluding author list and references), 19 figures; minor changes to match published version

Journal ref: JCAP 07 (2024) 047

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

Multiwavelength Variability Analysis of Fermi-LAT Blazars

Authors: P. Peñil, J. Otero-Santos, M. Ajello, S. Buson, A. Domínguez, L. Marcotulli, N. Torres-Albà, J. Becerra González, J. A. Acosta-Pulido

Abstract: Blazars present highly variable $γ$-ray emission. This variability, which can range from a few minutes to several years, is also observed at other wavelengths across the entire electromagnetic spectrum. We make use of the first 12 years of data from the Fermi Large Area Telescope (LAT), complemented with multiwavelength (MWL) archival data from different observatories and facilities in radio, infr… ▽ More Blazars present highly variable $γ$-ray emission. This variability, which can range from a few minutes to several years, is also observed at other wavelengths across the entire electromagnetic spectrum. We make use of the first 12 years of data from the Fermi Large Area Telescope (LAT), complemented with multiwavelength (MWL) archival data from different observatories and facilities in radio, infrared and optical bands, to study the possible periodic emission from 19 blazars previously claimed as periodic candidates. A periodicity analysis is performed with a pipeline for periodicity searches. Moreover, we study the cross-correlations between the $γ$-ray and MWL light curves. Additionally, we use the fractional variability and the structure function to evaluate the variability timescales. We find five blazars showing hints of periodic modulation with $\geq$3.0$σ$ ($\approx$0$σ$ post-trials), with periods ranging from 1.2 to 4 years, both in their $γ$-ray and MWL emission. The results provide clues for understanding the physical mechanisms generating the observed periodicity. △ Less

Submitted 3 March, 2024; originally announced March 2024.

Comments: 21 pages, 5 figures

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

Optical spectroscopy of blazars for the Cherenkov Telescope Array -- III

Authors: F. D'Ammando, P. Goldoni, W. Max-Moerbeck, J. Becerra Gonzalez, E. Kasai, D. A. Williams, N. Alvarez-Crespo, M. Backes, U. Barres de Almeida, C. Boisson, G. Cotter, V. Fallah Ramazani, O. Hervet, E. Lindfors, D. Mukhi-Nilo, S. Pita, M. Splettstoesser, B. van Soelen

Abstract: Due to their almost featureless optical/UV spectra, it is challenging to measure the redshifts of BL Lacs. As a result, about 50% of gamma-ray BL Lacs lack a firm measurement of this property, which is fundamental for population studies, indirect estimates of the EBL, and fundamental physics probes. This paper is the third in a series of papers aimed at determining the redshift of a sample of blaz… ▽ More Due to their almost featureless optical/UV spectra, it is challenging to measure the redshifts of BL Lacs. As a result, about 50% of gamma-ray BL Lacs lack a firm measurement of this property, which is fundamental for population studies, indirect estimates of the EBL, and fundamental physics probes. This paper is the third in a series of papers aimed at determining the redshift of a sample of blazars selected as prime targets for future observations with the next generation, ground-based VHE gamma-ray astronomy observatory, Cherenkov Telescope Array Observatory (CTAO). The accurate determination of the redshift of these objects is an important aid in source selection and planning of future CTAO observations. The selected targets were expected to be detectable with CTAO in observations of 30 hours or less. We performed deep spectroscopic observations of 41 of these blazars using the Keck II, Lick, SALT, GTC, and ESO/VLT telescopes. We carefully searched for spectral lines in the spectra and whenever features of the host galaxy were detected, we attempted to model the properties of the host galaxy. The magnitudes of the targets at the time of the observations were also compared to their long-term light curves. Spectra from 24 objects display spectral features or a high S/N. From these, 12 spectroscopic redshifts were determined, ranging from 0.2223 to 0.7018. Furthermore, 1 tentative redshift (0.6622) and 2 redshift lower limits at z > 0.6185 and z > 0.6347 were obtained. The other 9 BL Lacs showed featureless spectra, despite the high S/N (> 100) observations. Our comparisons with long-term optical light curves tentatively suggest that redshift measurements are more straightforward during an optical low state of the AGN. Overall, we have determined 37 redshifts and 6 spectroscopic lower limits as part of our programme thus far. △ Less

Submitted 15 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Astronomy and Astrophysics. 17 pages, 4 Figures, 10 Tables

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

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

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)

Variability studies of active galactic nuclei from the long-term monitoring program with the Cherenkov Telescope Array

Authors: Guillaume Grolleron, Josefa Becerra González, Jonathan Biteau, Matteo Cerruti, Roger Grau, Lucas Gréaux, Talvikki Hovatta, Jean-Philippe Lenain, Elina Lindfors, Walter Max-Moerbeck, Davide Miceli, Abelardo Moralejo, Kari Nilsson, Elisa Pueschel, Arkadipta Sarkar, Sofia Suutarinen

Abstract: Blazars are active galactic nuclei (AGN) with a relativistic jet oriented toward the observer. This jet is composed of accelerated particles which can display emission over the entire electromagnetic spectrum. Spectral variability has been observed on short- and long-time scales in AGN, with a power spectral density (PSD) that can show a break at frequencies below the well-known red-noise process.… ▽ More Blazars are active galactic nuclei (AGN) with a relativistic jet oriented toward the observer. This jet is composed of accelerated particles which can display emission over the entire electromagnetic spectrum. Spectral variability has been observed on short- and long-time scales in AGN, with a power spectral density (PSD) that can show a break at frequencies below the well-known red-noise process. This break frequency in the PSD has been observed in X-rays to scale with the accretion regime and the mass of the central black hole. It is expected that a break could also be seen in the very-high-energy gamma rays, but constraining the shape of the PSD in these wavelengths has not been possible with the current instruments. The Cherenkov Telescope Array (CTA) will be more sensitive by a factor of five to ten depending on energy than the current generation of imaging atmospheric Cherenkov telescopes, therefore it will be possible with CTA to reconstruct the PSD with a high accuracy, bringing new information about AGN variability. In this work, we focus on the AGN long-term monitoring program planned with CTA. The program is proposed to begin with early-start observing campaigns with CTA precursors. This would allow us to probe longer time scales on the AGN PSD. △ Less

Submitted 21 September, 2023; originally announced September 2023.

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

Report number: CTA-ICRC/2023/19

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

TeV pion bumps in the gamma-ray spectra of flaring blazars

Authors: M. Petropoulou, A. Mastichiadis, G. Vasilopoulos, D. Paneque, J. Becerra González, F. Zanias

Abstract: (Abridged) Very high-energy (VHE, $E>100$ GeV) observations of the blazar Mrk 501 with MAGIC in 2014 provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component. Motivated by this rare observation, we prop… ▽ More (Abridged) Very high-energy (VHE, $E>100$ GeV) observations of the blazar Mrk 501 with MAGIC in 2014 provided evidence for an unusual narrow spectral feature at about 3 TeV during an extreme X-ray flaring activity. The one-zone synchrotron-self Compton scenario, widely used in blazar broadband spectral modeling, fails to explain the narrow TeV component. Motivated by this rare observation, we propose an alternative model for the production of narrow features in the VHE spectra of flaring blazars. These spectral features may result from the decay of neutral pions ($π^0$ bumps) that are in turn produced via interactions of protons (of tens of TeV energy) with energetic photons, whose density increases during hard X-ray flares. We explored the conditions needed for the emergence of narrow $π^0$ bumps in VHE blazar spectra during X-ray flares reaching synchrotron energies $\sim 100$ keV using time-dependent radiative transfer calculations. We focused on high-synchrotron peaked (HSP) blazars, which comprise the majority of VHE-detected extragalactic sources. We find that synchrotron-dominated flares with peak energies $\gtrsim100$ keV can be ideal periods for the search of $π^0$ bumps in the VHE spectra of HSP blazars. The flaring region is optically thin to photopion production, its energy content is dominated by the relativistic proton population, and the inferred jet power is highly super-Eddington. Application of the model to the spectral energy distribution of Mrk 501 on MJD 56857.98 shows that the VHE spectrum of the flare is described well by the sum of a synchrotron self-Compton (SSC) component and a distinct $π^0$ bump centered at 3 TeV. Spectral fitting of simulated SSC+$π^0$ spectra for the Cherenkov Telescope Array (CTA) show that a $π^0$ bump could be detected at a 5$σ$ significance level with a 30-min exposure. △ Less

Submitted 5 March, 2024; v1 submitted 27 August, 2023; originally announced August 2023.

Comments: 10 pages, 7 figures, accepted version in A&A

A re-emerging bright soft-X-ray state of the changing-look Active Galactic Nucleus 1ES~1927+654: a multi-wavelength view

Authors: Ritesh Ghosh, Sibasish Laha, Eileen Meyer, Agniva Roychowdhury, Xiaolong Yang, J. A. Acosta Pulido, Suvendu Rakshit, Shivangi Pandey, Josefa Becerra Gonzalez, Ehud Behar, Luigi C. Gallo, Francesca Panessa, Stefano Bianchi, Fabio La Franca, Nicolas Scepi, Mitchell C. Begelman, Anna Lia Longinotti, Elisabeta Lusso, Samantha Oates, Matt Nicholl, S. Bradley Cenko, Brendan O Connor, Erica Hammerstein, Jincen Jose, Krisztina Eva Gabanyi , et al. (2 additional authors not shown)

Abstract: 1ES1927+654 is a nearby active galactic nucleus that has shown an enigmatic outburst in optical/UV followed by X-rays, exhibiting strange variability patterns at timescales of months-years. Here we report the unusual X-ray, UV, and radio variability of the source in its post-flare state (Jan 2022- May 2023). Firstly, we detect an increase in the soft X-ray (0.3-2 keV) flux from May 2022- May 2023… ▽ More 1ES1927+654 is a nearby active galactic nucleus that has shown an enigmatic outburst in optical/UV followed by X-rays, exhibiting strange variability patterns at timescales of months-years. Here we report the unusual X-ray, UV, and radio variability of the source in its post-flare state (Jan 2022- May 2023). Firstly, we detect an increase in the soft X-ray (0.3-2 keV) flux from May 2022- May 2023 by almost a factor of five, which we call the bright-soft-state. The hard X-ray 2-10 keV flux increased by a factor of two, while the UV flux density did not show any significant changes ($\le 30\%$) in the same period. The integrated energy pumped into the soft and hard X-ray during this period of eleven months is $\sim 3.57\times 10^{50}$ erg and $5.9\times 10^{49}$ erg, respectively. From the energetics, it is evident that whatever is producing the soft excess (SE) is pumping out more energy than either the UV or hard X-ray source. Since the energy source presumably is ultimately the accretion of matter onto the SMBH, the SE emitting region must be receiving the majority of this energy. In addition, the source does not follow the typical disc-corona relation found in AGNs, neither in the initial flare (in 2017-2019) nor in the current bright soft state (2022-2023). We found that the core (<1 pc) radio emission at 5 GHz gradually increased till March 2022 but showed a dip in August 2022. The Güdel Benz relation ($L_{\rm radio}/L_{\rm X-ray}\sim 10^{-5}$), however, is still within the expected range for radio-quiet AGN and further follow-up radio observations are currently being undertaken. △ Less

Submitted 11 August, 2023; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: Accepted for publication in ApJ

The Third Fermi Large Area Telescope Catalog of Gamma-ray Pulsars

Authors: David A. Smith, Philippe Bruel, Colin J. Clark, Lucas Guillemot, Matthew T. Kerr, Paul Ray, Soheila Abdollahi, Marco Ajello, Luca Baldini, Jean Ballet, Matthew Baring, Cees Bassa, Josefa Becerra Gonzalez, Ronaldo Bellazzini, Alessandra Berretta, Bhaswati Bhattacharyya, Elisabetta Bissaldi, Raffaella Bonino, Eugenio Bottacini, Johan Bregeon, Marta Burgay, Toby Burnett, Rob Cameron, Fernando Camilo, Regina Caputo , et al. (134 additional authors not shown)

Abstract: We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray M… ▽ More We present 294 pulsars found in GeV data from the Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope. Another 33 millisecond pulsars (MSPs) discovered in deep radio searches of LAT sources will likely reveal pulsations once phase-connected rotation ephemerides are achieved. A further dozen optical and/or X-ray binary systems co-located with LAT sources also likely harbor gamma-ray MSPs. This catalog thus reports roughly 340 gamma-ray pulsars and candidates, 10% of all known pulsars, compared to $\leq 11$ known before Fermi. Half of the gamma-ray pulsars are young. Of these, the half that are undetected in radio have a broader Galactic latitude distribution than the young radio-loud pulsars. The others are MSPs, with 6 undetected in radio. Overall, >235 are bright enough above 50 MeV to fit the pulse profile, the energy spectrum, or both. For the common two-peaked profiles, the gamma-ray peak closest to the magnetic pole crossing generally has a softer spectrum. The spectral energy distributions tend to narrow as the spindown power $\dot E$ decreases to its observed minimum near $10^{33}$ erg s$^{-1}$, approaching the shape for synchrotron radiation from monoenergetic electrons. We calculate gamma-ray luminosities when distances are available. Our all-sky gamma-ray sensitivity map is useful for population syntheses. The electronic catalog version provides gamma-ray pulsar ephemerides, properties and fit results to guide and be compared with modeling results. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 142 pages. Accepted by the Astrophysical Journal Supplement

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

Variability and evolution of the optical polarization of a sample of gamma-ray blazars

Authors: J. Otero-Santos, J. A. Acosta-Pulido, J. Becerra González, C. M. Raiteri, M. I. Carnerero, N. Castro Segura, O. González-Martín, A. Luashvili

Abstract: We present a polarization variability analysis of a sample of 26 $γ$-ray blazars monitored by the Steward Observatory between 2008 and 2018 in the optical band. We investigate the properties and long-term variability of their optical polarization, searching for differences between blazar types. We observe that BL Lac objects are typically less polarized and less variable than flat spectrum radio q… ▽ More We present a polarization variability analysis of a sample of 26 $γ$-ray blazars monitored by the Steward Observatory between 2008 and 2018 in the optical band. We investigate the properties and long-term variability of their optical polarization, searching for differences between blazar types. We observe that BL Lac objects are typically less polarized and less variable than flat spectrum radio quasars (FSRQs). Moreover, BL Lacs display a distribution of their polarization angle typically oriented in a preferential direction, contrary to the rather random distribution of FSRQs. For the latter blazar type, as well as those sources showing a bright stellar emission, we take into account the depolarizing effect introduced by the broad line region and the host galaxy on the measured polarization degree. In this sample we also observe that BL Lacs present an uncorrelated evolution of the flux and the polarization. Contrary, FSRQs show a correlation before the depolarization correction, that is lost however after considering this effect. In addition, we study the behaviour of the polarization angle, searching for angle rotations in its long-term evolution. We derive that the FSRQs studied here show rotations more frequently than BL Lac objects by a factor $\sim$1.5. During these periods we also observe a systematic decrease of the polarization fraction, as well as a marginal flux increase, not significant however to connect rotations with optical flares. We interpret these results within the extended shock-in-jet scenario, able to explain the overall features observed here for the polarization of the blazar sample. △ Less

Submitted 6 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS

A Contribution of the HAWC Observatory to the TeV era in the High Energy Gamma-Ray Astrophysics: The case of the TeV-Halos

Authors: Ramiro Torres-Escobedo, Hao Zhou, Eduardo de la Fuente, A. U. Abeysekara, A. Albert, R. Alfaro, C. Alvarez, J. D. Álvarez, J. R. Angeles Camacho, J. C. Arteaga-Velázquez, K. P. Arunbabu, D. Avila Rojas, H. A. Ayala Solares, R. Babu, V. Baghmanyan, A. S. Barber, J. Becerra Gonzalez, E. Belmont-Moreno, S. Y. BenZvi, D. Berley, C. Brisbois, K. S. Caballero-Mora, T. Capistrán, A. Carramiñana, S. Casanova , et al. (108 additional authors not shown)

Abstract: We present a short overview of the TeV-Halos objects as a discovery and a relevant contribution of the High Altitude Water Čerenkov (HAWC) observatory to TeV astrophysics. We discuss history, discovery, knowledge, and the next step through a new and more detailed analysis than the original study in 2017. TeV-Halos will contribute to resolving the problem of the local positron excess observed on th… ▽ More We present a short overview of the TeV-Halos objects as a discovery and a relevant contribution of the High Altitude Water Čerenkov (HAWC) observatory to TeV astrophysics. We discuss history, discovery, knowledge, and the next step through a new and more detailed analysis than the original study in 2017. TeV-Halos will contribute to resolving the problem of the local positron excess observed on the Earth. To clarify the latter, understanding the diffusion process is mandatory. △ Less

Submitted 13 April, 2023; originally announced April 2023.

Comments: Work presented in the 21st International Symposium on Very High Energy Cosmic Ray Interactions(ISVHECRI 2022) as part of the Ph. D. Thesis of Ramiro Torres-Escobedo (SJTU, Shanghai, China). Accepted for publication in SciPost Physics Proceedings (ISSN 2666-4003). 11 pages, 3 Figures. Short overview of HAWC and TeV Halos objects until 2022

Fermi-GBM Discovery of GRB 221009A: An Extraordinarily Bright GRB from Onset to Afterglow

Authors: S. Lesage, P. Veres, M. S. Briggs, A. Goldstein, D. Kocevski, E. Burns, C. A. Wilson-Hodge, P. N. Bhat, D. Huppenkothen, C. L. Fryer, R. Hamburg, J. Racusin, E. Bissaldi, W. H. Cleveland, S. Dalessi, C. Fletcher, M. M. Giles, B. A. Hristov, C. M. Hui, B. Mailyan, C. Malacaria, S. Poolakkil, O. J. Roberts, A. von Kienlin, J. Wood , et al. (115 additional authors not shown)

Abstract: We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing ana… ▽ More We report the discovery of GRB 221009A, the highest flux gamma-ray burst ever observed by the Fermi Gamma-ray Burst Monitor (GBM). This GRB has continuous prompt emission lasting more than 600 seconds which smoothly transitions to afterglow visible in the GBM energy range (8 keV--40 MeV), and total energetics higher than any other burst in the GBM sample. By using a variety of new and existing analysis techniques we probe the spectral and temporal evolution of GRB 221009A. We find no emission prior to the GBM trigger time (t0; 2022 October 9 at 13:16:59.99 UTC), indicating that this is the time of prompt emission onset. The triggering pulse exhibits distinct spectral and temporal properties suggestive of the thermal, photospheric emission of shock-breakout, with significant emission up to $\sim$15 MeV. We characterize the onset of external shock at t0+600 s and find evidence of a plateau region in the early-afterglow phase which transitions to a slope consistent with Swift-XRT afterglow measurements. We place the total energetics of GRB 221009A in context with the rest of the GBM sample and find that this GRB has the highest total isotropic-equivalent energy ($\textrm{E}_{γ,\textrm{iso}}=1.0\times10^{55}$ erg) and second highest isotropic-equivalent luminosity ($\textrm{L}_{γ,\textrm{iso}}=9.9\times10^{53}$ erg/s) based on redshift of z = 0.151. These extreme energetics are what allowed us to observe the continuously emitting central engine of GBM from the beginning of the prompt emission phase through the onset of early afterglow. △ Less

Submitted 12 July, 2023; v1 submitted 24 March, 2023; originally announced March 2023.

Comments: 26 pages 7 figures - accepted for publication in ApJL

Follow-up survey for the binary black hole merger GW200224_222234 using Subaru/HSC and GTC/OSIRIS

Authors: Takayuki Ohgami, Josefa Becerra Gonzalez, Nozomu Tominaga, Tomoki Morokuma, Yousuke Utsumi, Yuu Niino, Masaomi Tanaka, Smaranika Banerjee, Frederick Poidevin, Jose Antonio Acosta-Pulido, Ismael Perez-Fournon, Teo Munoz-Darias, Hiroshi Akitaya, Kenshi Yanagisawa, Mahito Sasada, Michitoshi Yoshida, Mirko Simunovic, Ryou Ohsawa, Ichi Tanaka, Tsuyoshi Terai, Yuhei Takagi, The J-GEM collaboration

Abstract: The LIGO/Virgo detected a gravitational wave (GW) event, named GW200224_222234 (a.k.a. S200224ca) and classified as a binary-black-hole coalescence, on February 24, 2020. Given its relatively small localization skymap (71 deg$^2$ for a 90% credible region; revised to 50 deg$^2$ in GWTC-3), we performed target-of-opportunity observations using the Subaru/Hyper Suprime-Cam (HSC) in the $r2$- and… ▽ More The LIGO/Virgo detected a gravitational wave (GW) event, named GW200224_222234 (a.k.a. S200224ca) and classified as a binary-black-hole coalescence, on February 24, 2020. Given its relatively small localization skymap (71 deg$^2$ for a 90% credible region; revised to 50 deg$^2$ in GWTC-3), we performed target-of-opportunity observations using the Subaru/Hyper Suprime-Cam (HSC) in the $r2$- and $z$-bands. Observations were conducted on February 25 and 28 and March 23, 2020, with the first epoch beginning 12.3 h after the GW detection. The survey covered the highest probability sky area of 56.6 deg$^2$, corresponding to a 91% probability. This was the first deep follow-up ($m_{r}\gtrsim24, m_{z}\gtrsim23$) for a binary-black-hole merger covering $>$90% of the localization. By performing image subtraction and candidate screening including light curve fitting with transient templates and examples, we found 22 off-nucleus transients that were not ruled out as the counterparts of GW200224_222234 with only our Subaru/HSC data. We also performed GTC/OSIRIS spectroscopy of the probable host galaxies for five candidates; two are likely to be located within the 3D skymap, whereas the others are not. In conclusion, 19 transients remain as possible optical counterparts of GW200224_222234; however, we could not identify a unique promising counterpart. If there are no counterparts in the remaining candidates, the upper limits of optical luminosity are $νL_ν < 5.2^{+2.4}_{-1.9}\times 10^{41}$ erg s$^{-1}$ and $νL_ν < 1.8^{+0.8}_{-0.6}\times 10^{42}$ erg s$^{-1}$ in the $r2$- and $z$-bands, respectively, at $\sim$12 h after GW detection. We also discuss improvements in the strategies of optical follow-ups for future GW events. △ Less

Submitted 18 February, 2023; originally announced February 2023.

Comments: 22 pages, 15 figures. Accepted for publication in the Astrophysical Journal

The Fermi-LAT Light Curve Repository

Authors: S. Abdollahi, M. Ajello, L. Baldini, J. Ballet, D. Bastieri, J. Becerra Gonzalez, R. Bellazzini, A. Berretta, E. Bissaldi, R. Bonino, A. Brill, P. Bruel, E. Burns, S. Buson, A. Cameron, R. Caputo, P. A. Caraveo, N. Cibrario, S. Ciprini, P. Cristarella Orestano, M. Crnogorcevic, S. Cutini, F. D'Ammando, S. De Gaetano, S. W. Digel , et al. (88 additional authors not shown)

Abstract: The Fermi Large Area Telescope (LAT) light curve repository (LCR) is a publicly available, continually updated library of gamma-ray light curves of variable Fermi-LAT sources generated over multiple timescales. The Fermi-LAT LCR aims to provide publication-quality light curves binned on timescales of 3 days, 7 days, and 30 days for 1525 sources deemed variable in the source catalog of the first 10… ▽ More The Fermi Large Area Telescope (LAT) light curve repository (LCR) is a publicly available, continually updated library of gamma-ray light curves of variable Fermi-LAT sources generated over multiple timescales. The Fermi-LAT LCR aims to provide publication-quality light curves binned on timescales of 3 days, 7 days, and 30 days for 1525 sources deemed variable in the source catalog of the first 10 years of Fermi-LAT observations. The repository consists of light curves generated through full likelihood analyses that model the sources and the surrounding region, providing fluxes and photon indices for each time bin. The LCR is intended as a resource for the time-domain and multi-messenger communities by allowing users to quickly search LAT data to identify correlated variability and flaring emission episodes from gamma-ray sources. We describe the sample selection and analysis employed by the LCR and provide an overview of the associated data access portal. △ Less

Submitted 14 February, 2023; v1 submitted 4 January, 2023; originally announced January 2023.

Comments: Accepted for publication in ApJ Supplement Series

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)

Optical Spectroscopy of Blazars for the Cherenkov Telescope Array -- II

Authors: E. Kasai, P. Goldoni, S. Pita, D. A. Williams, W. Max-Moerbeck, O. Hervet, G. Cotter, M. Backes, C. Boisson, J. Becerra González, U. Barres de Almeida, F. D'Ammando, V. Fallah Ramazani, E. Lindfors

Abstract: Active galactic nuclei (AGNs) make up about 35 per cent of the more than 250 sources detected in very-high-energy (VHE) gamma rays to date with Imaging Atmospheric Cherenkov Telescopes. Apart from four nearby radio galaxies and two AGNs of unknown type, all known VHE AGNs are blazars. Knowledge of the cosmological redshift of gamma-ray blazars is key to enabling the study of their intrinsic emissi… ▽ More Active galactic nuclei (AGNs) make up about 35 per cent of the more than 250 sources detected in very-high-energy (VHE) gamma rays to date with Imaging Atmospheric Cherenkov Telescopes. Apart from four nearby radio galaxies and two AGNs of unknown type, all known VHE AGNs are blazars. Knowledge of the cosmological redshift of gamma-ray blazars is key to enabling the study of their intrinsic emission properties, as the interaction between gamma rays and the extragalactic background light (EBL) results in a spectral softening. Therefore, the redshift determination exercise is crucial to indirectly placing tight constraints on the EBL density and to studying blazar population evolution across cosmic time. Due to the powerful relativistic jets in blazars, most of their host galaxies' spectral features are outshined, and dedicated high signal-to-noise spectroscopic observations are required. Deep medium- to high-resolution spectroscopy of 33 gamma-ray blazar optical counterparts was performed with the European Southern Observatory New Technology Telescope, Keck II telescope, Shane 3-meter telescope and the Southern African Large Telescope. From the sample, spectra from 25 objects display spectral features or are featureless and have high signal-to-noise. The other eight objects have low quality featureless spectra. We systematically searched for absorption and emission features and estimated, when possible, the fractional host galaxy flux in the measured total flux. Our measurements yielded 14 firm spectroscopic redshifts, ranging from 0.0838 to 0.8125, one tentative redshift, and two lower limits: one at z > 0.382 and the other at z > 0.629. △ Less

Submitted 9 November, 2022; originally announced November 2022.

Comments: 19 pages, 5 figures This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society (https://doi.org/10.1093/mnras/stac3167) Published by Oxford University Press on behalf of the Royal Astronomical Society

Multiwavelength periodicity search in a sample of $γ$-ray bright blazars

Authors: J. Otero-Santos, P. Peñil, J. A. Acosta-Pulido, J. Becerra González, C. M. Raiteri, M. I. Carnerero, M. Villata

Abstract: We present the results of a long-term periodicity search in a sample of $γ$-ray blazars within a multiwavelength context. These blazars have been selected from the Steward Observatory sample as part of its optical monitoring program between 2008 and 2018. We study 15 sources with a temporal coverage in their optical total and polarized emission sufficiently large ($>9$ years) to perform a reliable… ▽ More We present the results of a long-term periodicity search in a sample of $γ$-ray blazars within a multiwavelength context. These blazars have been selected from the Steward Observatory sample as part of its optical monitoring program between 2008 and 2018. We study 15 sources with a temporal coverage in their optical total and polarized emission sufficiently large ($>9$ years) to perform a reliable long-term periodicity analysis. We collect data from several observatories to extend the coverage, enabling the search of longer periods. In addition, data are also gathered in the high-energy ($E>100$ MeV) $γ$-ray band from the $\textit{Fermi}$ Large Area Telescope; and in the 15-GHz radio band from the Owens Valley Radio Observatory. We identify 5 promising candidates to host quasi-periodic emission, AO 0235+164, PKS 1222+216, Mrk 501, BL Lacertae and 1ES 2344+514 with periods in one or more bands and statistical significances $\sim$3$σ$ after trial factor correction. AO 0235+164 shows a period of $\sim$8.2 years in the R band; PKS 1222+216 has a quasi-periodic modulation in its total and polarized optical emission of $\sim$1.6 years; Mrk 501 displays a $\sim$5-year quasi-periodicity in optical and radio wavelengths; BL Lacertae presents a period of $\sim$1.8 years in its polarized emission; and 1ES 2344+514 shows a hint of a $\sim$5.5-year period in its optical R band. We interpret these results in the framework of the most common models and scenarios, namely the presence of a binary supermassive black hole system; or geometrical effects like helical or precessing jets. △ Less

Submitted 11 November, 2022; v1 submitted 28 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS

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)

The Fourth Catalog of Active Galactic Nuclei Detected by the Fermi Large Area Telescope -- Data Release 3

Authors: The Fermi-LAT collaboration, :, Marco Ajello, Luca Baldini, Jean Ballet, Denis Bastieri, Josefa Becerra Gonzalez, Ronaldo Bellazzini, Alessandra Berretta, Elisabetta Bissaldi, Raffaella Bonino, Ari Brill, Philippe Bruel, Sara Buson, Regina Caputo, Patrizia Caraveo, Teddy Cheung, Graziano Chiaro, Nicolo Cibrario, Stefano Ciprini, Milena Crnogorcevic, Sara Cutini, Filippo D'Ammando, Salvatore De Gaetano, Niccolo Di Lalla , et al. (79 additional authors not shown)

Abstract: An incremental version of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi-Large Area Telescope is presented. This version (4LAC-DR3) derives from the third data release of the 4FGL catalog based on 12 years of E>50 MeV gamma-ray data, where the spectral parameters, spectral energy distributions (SEDs), yearly light curves, and associations have been updated for all source… ▽ More An incremental version of the fourth catalog of active galactic nuclei (AGNs) detected by the Fermi-Large Area Telescope is presented. This version (4LAC-DR3) derives from the third data release of the 4FGL catalog based on 12 years of E>50 MeV gamma-ray data, where the spectral parameters, spectral energy distributions (SEDs), yearly light curves, and associations have been updated for all sources. The new reported AGNs include 587 blazar candidates and four radio galaxies. We describe the properties of the new sample and outline changes affecting the previously published one. We also introduce two new parameters in this release, namely the peak energy of the SED high-energy component and the corresponding flux. These parameters allow an assessment of the Compton dominance, the ratio of the Inverse-Compton to the synchrotron peak luminosities, without relying on X-ray data. △ Less

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

Comments: 13 pages, 7 figures, accepted for publication in The Astrophysical Journal Supplement Series. Fits files are available at https://fermi.gsfc.nasa.gov/ssc/data/access/lat/4LACDR3/ or alternatively at https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-h.fits and https://www.ssdc.asi.it/fermi4lac-DR3/table-4LAC-DR3-l.fits

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

ASTRI Mini-Array Core Science at the Observatorio del Teide

Authors: S. Vercellone, C. Bigongiari, A. Burtovoi, M. Cardillo, O. Catalano, A. Franceschini, S. Lombardi, L. Nava, F. Pintore, A. Stamerra, F. Tavecchio, L. Zampieri, R. Alves Batista, E. Amato, L. A. Antonelli, C. Arcaro, J. Becerra Gonzalez, G. Bonnoli, M. Bottcher, G. Brunetti, A. A. Compagnino, S. Crestan, A. D Ai, M. Fiori, G. Galanti , et al. (62 additional authors not shown)

Abstract: The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Project led by the Italian National Institute for Astrophysics (INAF) is developing and will deploy at the Observatorio del Teide a mini-array (ASTRI Mini-Array) composed of nine telescopes similar to the small-size dual-mirror Schwarzschild-Couder telescope (ASTRI-Horn) currently operating on the slopes of Mt. Etna in Sicily.… ▽ More The ASTRI (Astrofisica con Specchi a Tecnologia Replicante Italiana) Project led by the Italian National Institute for Astrophysics (INAF) is developing and will deploy at the Observatorio del Teide a mini-array (ASTRI Mini-Array) composed of nine telescopes similar to the small-size dual-mirror Schwarzschild-Couder telescope (ASTRI-Horn) currently operating on the slopes of Mt. Etna in Sicily. The ASTRI Mini-Array will surpass the current Cherenkov telescope array differential sensitivity above a few tera-electronvolt (TeV), extending the energy band well above hundreds of TeV. This will allow us to explore a new window of the electromagnetic spectrum, by convolving the sensitivity performance with excellent angular and energy resolution figures. In this paper we describe the Core Science that we will address during the first four years of operation, providing examples of the breakthrough results that we will obtain when dealing with current open questions, such as the acceleration of cosmic rays, cosmology and fundamental physics and the new window, for the TeV energy band, of the time-domain astrophysics. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: Published in Journal of High Energy Astrophysics. 46 Figures, 7 Tables

Journal ref: Journal of High Energy Astrophysics, Volume 35, August 2022, Pages 1-42

Extragalactic Observatory Science with the ASTRI Mini-Array at the Observatorio del Teide

Authors: F. G. Saturni, C. H. E. Arcaro, B. Balmaverde, J. Becerra González, A. Caccianiga, M. Capalbi, A. Lamastra, S. Lombardi, F. Lucarelli, R. Alves Batista, L. A. Antonelli, E. M. de Gouveia Dal Pino, R. Della Ceca, J. G. Green, A. Pagliaro, C. Righi, F. Tavecchio, S. Vercellone, A. Wolter, E. Amato, C. Bigongiari, M. Böttcher, G. Brunetti, P. Bruno, A. Bulgarelli , et al. (25 additional authors not shown)

Abstract: The ASTRI Mini-Array is a next-generation system of nine imaging atmospheric Cherenkov telescopes that is going to be built at the Observatorio del Teide site. After a first phase, in which the instrument will be operated as an experiment prioritizing a schedule of primary science cases, an observatory phase is foreseen in which other significant targets will be pointed. We focus on the observatio… ▽ More The ASTRI Mini-Array is a next-generation system of nine imaging atmospheric Cherenkov telescopes that is going to be built at the Observatorio del Teide site. After a first phase, in which the instrument will be operated as an experiment prioritizing a schedule of primary science cases, an observatory phase is foreseen in which other significant targets will be pointed. We focus on the observational feasibility of extragalactic sources and on astrophysical processes that best complement and expand the ASTRI Mini-Array core science, presenting the most relevant examples that are at reach of detection over long-term time scales and whose observation can provide breakthrough achievements in the very-high energy extragalactic science. Such examples cover a wide range of $γ$-ray emitters, including the study of AGN low states in the multi-TeV energy range, the possible detection of Seyfert galaxies with long exposures and the searches of dark matter lines above 10 TeV. Simulations of the presented objects show that the instrument performance will be competitive at multi-TeV energies with respect to current arrays of Cherenkov telescopes. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: 27 pages, 10 figures, 7 tables, published on JHEAp

Journal ref: Journal of High Energy Astrophysics 35 (2022), 91-111

The relation between optical and $γ$-ray emission in BL Lac sources

Authors: G. La Mura, J. Becerra Gonzalez, G. Chiaro, S. Ciroi, J. Otero-Santos

Abstract: The relativistic jets produced by some Active Galactic Nuclei (AGNs) are among the most efficient persistent sources of non-thermal radiation and represent an ideal laboratory for studying high-energy interactions. In particular, when the relativistic jet propagates along the observer's line of sight, the beaming effect produces dominant signatures in the observed spectral energy distribution (SED… ▽ More The relativistic jets produced by some Active Galactic Nuclei (AGNs) are among the most efficient persistent sources of non-thermal radiation and represent an ideal laboratory for studying high-energy interactions. In particular, when the relativistic jet propagates along the observer's line of sight, the beaming effect produces dominant signatures in the observed spectral energy distribution (SED), from the radio domain up to the highest energies, with the further possibility of resulting in radiation-particle multi-messenger associations. In this work, we investigate the relationships between the emission of $γ$ rays and the optical spectra of a sample of AGN, selected from BL Lac sources detected by the $Fermi$ Large Area Telescope ($Fermi$-LAT). We find that there is a close relationship between the optical and $γ$-ray spectral indices. Despite all the limitations due to the non-simultaneity of the data, this observation strongly supports a substantial role of Synchrotron-Self Compton (SSC) radiation in a single zone leptonic scenario for most sources. This result simplifies the application of theoretical models to explore the physical parameters of the jets in this type of sources. △ Less

Submitted 17 April, 2023; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: 20 pages, 8 figures, published on MNRAS. Fixed error in target coordinates list

Spectral variability of the born-again ejecta in A 58

Authors: Borja Montoro Molina, Martín A. Guerrero, Jesús A. Toalá, Janis B. Rodríguez González

Abstract: Born-again planetary nebulae (PNe) allow investigating stellar evolution, dust production, and nebular shocks in human timescales. Here we present an analysis of multi-epoch optical spectroscopic observations of the born-again PN A 58 around V605 Aql, which experienced a very late thermal pulse (VLTP) about a century ago. The H-deficient ejecta has experienced a considerable brightening in the tim… ▽ More Born-again planetary nebulae (PNe) allow investigating stellar evolution, dust production, and nebular shocks in human timescales. Here we present an analysis of multi-epoch optical spectroscopic observations of the born-again PN A 58 around V605 Aql, which experienced a very late thermal pulse (VLTP) about a century ago. The H-deficient ejecta has experienced a considerable brightening in the time period considered, from 1996 to 2021, with notable changes also in many emission line ratios. Neither the reduction of the extinction caused by the dilution of the ejecta nor the increase of the ionizing photon flux from the central star seem capable to produce these spectral changes, which are instead attributed to shocks in the bipolar H-poor outflow dissociating molecular material and propagating through the outer nebula. △ Less

Submitted 21 June, 2022; originally announced June 2022.

Comments: 9 pages, 5 figures

A deep survey of short GRB host galaxies over $z\sim0-2$: implications for offsets, redshifts, and environments

Authors: B. O'Connor, E. Troja, S. Dichiara, P. Beniamini, S. B. Cenko, C. Kouveliotou, J. Becerra Gonzalez, J. Durbak, P. Gatkine, A. Kutyrev, T. Sakamoto, R. Sanchez-Ramirez, S. Veilleux

Abstract: A significant fraction (30\%) of well-localized short gamma-ray bursts (sGRBs) lack a coincident host galaxy. This leads to two main scenarios: \textit{i}) that the progenitor system merged outside of the visible light of its host, or \textit{ii}) that the sGRB resided within a faint and distant galaxy that was not detected by follow-up observations. Discriminating between these scenarios has impo… ▽ More A significant fraction (30\%) of well-localized short gamma-ray bursts (sGRBs) lack a coincident host galaxy. This leads to two main scenarios: \textit{i}) that the progenitor system merged outside of the visible light of its host, or \textit{ii}) that the sGRB resided within a faint and distant galaxy that was not detected by follow-up observations. Discriminating between these scenarios has important implications for constraining the formation channels of neutron star mergers, the rate and environments of gravitational wave sources, and the production of heavy elements in the Universe. In this work, we present the results of our observing campaign targeted at 31 sGRBs that lack a putative host galaxy. Our study effectively doubles the sample of well-studied sGRB host galaxies, now totaling 72 events of which $28\%$ lack a coincident host to deep limits ($r$\,$\gtrsim$\,$26$ or $F110W$\,$\gtrsim$\,$27$ AB mag), and represents the largest homogeneously selected catalog of sGRB offsets to date. We find that 70\% of sub-arcsecond localized sGRBs occur within 10 kpc of their host's nucleus, with a median projected physical offset of $5.6$ kpc. Using this larger population, we discover an apparent redshift evolution in their locations: bursts at low-$z$ occur at $2\times$ larger offsets compared to those at $z$\,$>$\,$0.5$. This evolution could be due to a physical evolution of the host galaxies themselves or a bias against faint high-$z$ galaxies. Furthermore, we discover a sample of hostless sGRBs at $z$\,$\gtrsim$\,$1$ that are indicative of a larger high-$z$ population, constraining the redshift distribution and disfavoring log-normal delay time models. △ Less

Submitted 11 July, 2022; v1 submitted 19 April, 2022; originally announced April 2022.

Comments: Accepted by MNRAS

Limits on the hard X-ray emission from the periodic fast radio burst FRB 180916.J0158+65

Authors: Sibasish Laha, Zorawar Wadiasingh, Tyler Parsotan, Amy Lien, George Younes, Bing Zhang, S. Bradley Cenko, Eleonora Troja, Samantha Oates, Matt Nicholl, Eileen Meyer, Josefa Becerra González, Ritesh Ghosh, Noel Klingler

Abstract: FRB 180916.J0158+65 is one of the nearest, periodically repeating, and actively bursting fast radio burst (FRB) which has been localized to the outskirts of a spiral galaxy. In this work we study the FRB with the hard X-ray $14-195$ keV data from the Burst Alert Telescope (BAT) on board The Neil Gehrels Swift Observatory. BAT uses coded mask technology giving a localization of $\lesssim 3$ arc-min… ▽ More FRB 180916.J0158+65 is one of the nearest, periodically repeating, and actively bursting fast radio burst (FRB) which has been localized to the outskirts of a spiral galaxy. In this work we study the FRB with the hard X-ray $14-195$ keV data from the Burst Alert Telescope (BAT) on board The Neil Gehrels Swift Observatory. BAT uses coded mask technology giving a localization of $\lesssim 3$ arc-minute in the hard X-ray band, along with an accurate background estimation. BAT has been observing the source location in survey mode since February 2020. The survey mode observations involves background subtracted spectra, integrated over a time span ranging $300-2000$ seconds, at the source location (from Feb 2020-Jan 2022). We analyzed all the $\sim 230$ survey mode observations from BAT and checked for any signal in any of the observations. We did not detect any signal at $>5σ$ confidence level in any of the observations. We could estimate a $5σ$ upper limit on the $14-195$ keV flux, which ranged between $4.5\times 10^{-10} - 7.6\times 10^{-9}\, \rm erg\, cm^{-2}\, s^{-1}$. At the source distance this relates to a $5σ$ upper limit on luminosity of $5.08\times 10^{44}- 8.5\times 10^{45} \rm erg\, s^{-1}$. With this estimate, we could rule out any persistent X-ray emission, at the source location for these snapshots of BAT observations. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: Accepted for publication in ApJ

Simultaneous view of the FRB~180301 with FAST and NICER during a bursting phase

Authors: Sibasish Laha, George Younes, Zorawar Wadiasingh, Bo-Jun Wang, Ke-Jia Lee, Noel Klingler, Bing Zhang, Heng Xu, Chin-Feng Zhang, Wei-Wei Zhu, Ritesh Ghosh, Amy Lien, Eleonora Troja, S. Bradley Cenko, Samantha Oates, Matt Nicholl, Josefa Becerra González, Eileen Meyer, Tyler Parsotan

Abstract: FRB180301 is one of the most actively repeating fast radio bursts (FRBs) which has shown polarization angle changes in its radio burst emission, an indication for their likely origin in the magnetosphere of a highly-magnetized neutron star. We carried out a multi-wavelength campaign with the FAST radio telescope and NICER X-ray observatory to investigate any possible X-ray emission temporally coin… ▽ More FRB180301 is one of the most actively repeating fast radio bursts (FRBs) which has shown polarization angle changes in its radio burst emission, an indication for their likely origin in the magnetosphere of a highly-magnetized neutron star. We carried out a multi-wavelength campaign with the FAST radio telescope and NICER X-ray observatory to investigate any possible X-ray emission temporally coincident with the bright radio bursts. The observations took place on 2021 March 4, 9 and 19. We detected five bright radio bursts with FAST, four of which were strictly simultaneous with the NICER observations. The peak flux-density of the radio bursts ranged between $28-105$ mJy, the burst fluence between $27-170$ mJy-ms, and the burst durations between $1.7-12.3$ ms. The radio bursts from FRB~180301 exhibited complex time domain structure, and sub-pulses were detected in individual bursts, with no significant circular polarisation. The linear degree of polarisation in L-band reduced significantly compared to the 2019 observations. We do not detect any X-ray emission in excess of the background during the 5ms, 10ms, 100ms, 1sec and 100sec time intervals at/around the radio-burst barycenter-corrected arrival times, at a $>5σ$ confidence level. The $5σ$ upper limits on the X-ray a) persistent flux is $<7.64\times 10^{-12}\, \rm erg\, cm^{-2}\, s^{-1}$ , equivalent to $L_{\rm X}<2.50 \times 10^{45} \rm erg\, s^{-1}$ and b) 5 ms fluence is $<2\times 10^{-11} \rm erg\, cm^{-2}$, at the radio burst regions. Using the $5$ ms X-ray fluence upper limit, we can estimate the radio efficiency $η_{R/X} \equiv L_{\rm Radio}/L_{\rm X-ray} \gtrsim 10^{-8}$. The derived upper limit on $η_{R/X}$ is consistent with both magnetospheric models and synchrotron maser models involving relativistic shocks. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: Resubmitted to ApJ following referee's minor comments

A radio, optical, UV and X-ray view of the enigmatic changing look Active Galactic Nucleus 1ES~1927+654 from its pre- to post-flare states

Authors: Sibasish Laha, Eileen Meyer, Agniva Roychowdhury, Josefa Becerra González, J. A. Acosta-Pulido, Aditya Thapa, Ritesh Ghosh, Ehud Behar, Luigi C. Gallo, Gerard A. Kriss, Francesca Panessa, Stefano Bianchi, Fabio La Franca, Nicolas Scepi, Mitchell C. Begelman, Anna Lia Longinotti, Elisabeta Lusso, Samantha Oates, Matt Nicholl, S. Bradley Cenko

Abstract: The nearby type-II AGN 1ES1927+654 went through a violent changing-look (CL) event beginning December 2017 during which the optical and UV fluxes increased by four magnitudes over a few months, and broad emission lines newly appeared in the optical/UV. By July 2018 the X-ray coronal emission had completely vanished, only to reappear a few months later. In this work we report the evolution of the r… ▽ More The nearby type-II AGN 1ES1927+654 went through a violent changing-look (CL) event beginning December 2017 during which the optical and UV fluxes increased by four magnitudes over a few months, and broad emission lines newly appeared in the optical/UV. By July 2018 the X-ray coronal emission had completely vanished, only to reappear a few months later. In this work we report the evolution of the radio, optical, UV and X-rays from the pre-flare state through mid-2021 with new and archival data from the Very Long Baseline Array (VLBA), the European VLBI Network, the Very Large Array (VLA), the Telescopio Nazionale Galileo (TNG), Gran Telescopio Canarias (GTC), The Neil Gehrels Swift observatory and XMM-Newton. The main results from our work are: (1) The source has returned to its pre-CL state in optical, UV, and X-ray; the disk-corona relation has been re-established as has been in the pre-CL state, with an $α_{\rm OX}\sim 1.02$. The optical spectra are dominated by narrow emission lines. (2) The UV light curve follows a shallower slope of $\propto t^{-0.91\pm 0.04}$ compared to that predicted by a tidal disruption event. We conjecture that a magnetic flux inversion event is the possible cause for this enigmatic event. (3) The compact radio emission which we tracked in the pre-CL (2014), during CL (2018) and post-CL(2021) at spatial scales $<1$ pc was at its lowest level during the changing look event in 2018, nearly contemporaneous with a low $2-10$ keV emission. The radio to X-ray ratio of the compact source $L_{\rm Radio}/L_{\rm X-ray}\sim 10^{-5.5}$, follows the Gudel-Benz relation, typically found in coronally active stars, and several AGN. (4) We do not detect any presence of nascent jets at the spatial scales of $\sim 5-10$ pc. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Comments: Resubmitted to ApJ following minor comments from the referee

A statistical study of the optical spectral variability in gamma-ray blazars

Authors: J. Otero-Santos, J. A. Acosta-Pulido, J. Becerra González, A. Luashvili, N. Castro Segura, O. González-Martín, C. M. Raiteri, M. I. Carnerero

Abstract: Blazars optical emission is generally dominated by relativistic jets, although the host galaxy, accretion disk and broad line region (BLR) may also contribute significantly. Disentangling their contributions has been challenging for years due to the dominance of the jet. To quantify the contributions to the spectral variability, we use the statistical technique for dimensionality reduction Non-Neg… ▽ More Blazars optical emission is generally dominated by relativistic jets, although the host galaxy, accretion disk and broad line region (BLR) may also contribute significantly. Disentangling their contributions has been challenging for years due to the dominance of the jet. To quantify the contributions to the spectral variability, we use the statistical technique for dimensionality reduction Non-Negative Matrix Factorization on a spectroscopic data set of 26 $γ$-ray blazars. This technique allows to model large numbers of spectra in terms of a reduced number of components.We use a priori knowledge to obtain components associated to meaningful physical processes. The sources are classified according to their optical spectrum as host-galaxy dominated BL Lac objects (BL Lacs), BL Lacs, or Flat Spectrum Radio Quasars (FSRQs). Host-galaxy sources show less variability, as expected, and bluer-when-brighter trends, as the other BL Lacs. For FSRQs, more complicated colour-flux behaviours are observed: redder-when-brighter for low states saturating above a certain level and, in some cases, turning to bluer-when-brighter. We are able to reproduce the variability observed during 10 years using only 2 to 4 components, depending on the type. The simplest scenario corresponds to host-galaxy blazars, whose spectra are reconstructed using the stellar population and a power law for the jet. BL Lac spectra are reproduced using from 2 to 4 power laws. Different components can be associated to acceleration/cooling processes taking place in the jet. The reconstruction of FSRQs also incorporates a QSO-like component to account for the BLR, plus a very steep power law, associated to the accretion disk. △ Less

Submitted 17 February, 2022; originally announced February 2022.

Comments: Accepted for publication in MNRAS

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

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

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

Gamma-ray signatures from pair cascades in recombination-line radiation fields

Authors: Christoph Wendel, Josefa Becerra González, Amit Shukla, David Paneque, Karl Mannheim

Abstract: Beams of ultra-relativistic electrons in blazar jets develop pair cascades interacting with ambient soft photons. Employing coupled kinetic equations with escape terms, we model the unsaturated pair cascade spectrum. We assume that the gamma rays predominantly scatter off recombination-line photons from clouds photoionised by the irradiation from the accretion disk and the jet. The cascade spectru… ▽ More Beams of ultra-relativistic electrons in blazar jets develop pair cascades interacting with ambient soft photons. Employing coupled kinetic equations with escape terms, we model the unsaturated pair cascade spectrum. We assume that the gamma rays predominantly scatter off recombination-line photons from clouds photoionised by the irradiation from the accretion disk and the jet. The cascade spectrum is rather insensitive to the injection of hard electron spectra associated with the short-time variability of blazars. Adopting physical parameters representative of Markarian 501 and 3C 279, respectively, we numerically obtain spectral energy distributions showing distinct features imprinted by the recombination-line photons. The hints for a peculiar feature at ~3 TeV in the spectrum of Markarian 501, detected with the MAGIC telescopes during a strong X-ray flux activity in 2014 July, can be explained in this scenario as a result of the up-scattering of line photons by beam electrons and the low pair-creation optical depth. Inspecting a high-fidelity Fermi-LAT spectrum of 3C 279 from January 2018 reveals troughs in the spectrum that coincide with the threshold energies for gamma rays producing pairs in collisions with recombination-line photons, and the absence of exponential attenuation. Our finding implies that the gamma rays in 3C 279 escape from the edge of the broad emission line region. △ Less

Submitted 5 August, 2021; originally announced August 2021.

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