Showing 1–50 of 202 results for author: Neronov, A

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 Ultra Heavy Dark Matter Properties from Dwarf Spheroidal Galaxies with LHAASO Observations

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes… ▽ More In this work we try to search for signals generated by ultra-heavy dark matter at the Large High Altitude Air Shower Observatory (LHAASO) data. We look for possible gamma-ray by dark matter annihilation or decay from 16 dwarf spheroidal galaxies in the field of view of LHAASO. Dwarf spheroidal galaxies are among the most promising targets for indirect detection of dark matter which have low fluxes of astrophysical $γ$-ray background while large amount of dark matter. By analyzing more than 700 days observational data at LHAASO, no significant dark matter signal from 1 TeV to 1 EeV is detected. Accordingly we derive the most stringent constraints on the ultra-heavy dark matter annihilation cross-section up to EeV. The constraints on the lifetime of dark matter in decay mode are also derived. △ Less

Submitted 12 June, 2024; originally announced June 2024.

Comments: 17 pages, 12 figures, accepted by PRL

Data quality control system and long-term performance monitor of the LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, W. Bian, A. V. Bukevich, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, H. X. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. Chen , et al. (263 additional authors not shown)

Abstract: The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To… ▽ More The KM2A is the largest sub-array of the Large High Altitude Air Shower Observatory (LHAASO). It consists of 5216 electromagnetic particle detectors (EDs) and 1188 muon detectors (MDs). The data recorded by the EDs and MDs are used to reconstruct primary information of cosmic ray and gamma-ray showers. This information is used for physical analysis in gamma-ray astronomy and cosmic ray physics. To ensure the reliability of the LHAASO-KM2A data, a three-level quality control system has been established. It is used to monitor the status of detector units, stability of reconstructed parameters and the performance of the array based on observations of the Crab Nebula and Moon shadow. This paper will introduce the control system and its application on the LHAASO-KM2A data collected from August 2021 to July 2023. During this period, the pointing and angular resolution of the array were stable. From the observations of the Moon shadow and Crab Nebula, the results achieved using the two methods are consistent with each other. According to the observation of the Crab Nebula at energies from 25 TeV to 100 TeV, the time averaged pointing errors are estimated to be $-0.003^{\circ} \pm 0.005^{\circ}$ and $0.001^{\circ} \pm 0.006^{\circ}$ in the R.A. and Dec directions, respectively. △ Less

Submitted 13 June, 2024; v1 submitted 20 May, 2024; originally announced May 2024.

Comments: 15 pages, 9 figures

Pulsar timing array constraints on the cosmological magnetic field from quark confinement epoch

Authors: T. Boyer, A. Neronov

Abstract: Recent evidence for the stochastic gravitational wave backgorund reported by the pulsar timing arrays (PTA) can be interpreted as a signal from the cosmological phase transition. We use up-to-date models of the gravitational wave power spectra to compare constraints on the parameters of the phase transition for the three different available PTA measurements and to work out a refined estimate of th… ▽ More Recent evidence for the stochastic gravitational wave backgorund reported by the pulsar timing arrays (PTA) can be interpreted as a signal from the cosmological phase transition. We use up-to-date models of the gravitational wave power spectra to compare constraints on the parameters of the phase transition for the three different available PTA measurements and to work out a refined estimate of the cosmological magnetic field that should result from this transition. We find that the PTA data, combined with a constraint from the abundance of primordial black holes, are consistent with a possibility of a moderate strength first-order phase transition during quark confinement and yield a rather precise prediction for the initial parameters of the magnetic field, with the magnetic field energy density in near equipartition with photon energy density and correlation length close to one co-moving parsec. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 8 pages, 5 figures

Discovery of Very-high-energy Gamma-ray Emissions from the Low Luminosity AGN NGC 4278 by LHAASO

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) i… ▽ More The first source catalog of Large High Altitude Air Shower Observatory reported the detection of a very-high-energy gamma ray source, 1LHAASO J1219+2915. In this paper a further detailed study of the spectral and temporal behavior of this point-like source have been carried. The best-fit position of the TeV source ($\rm{RA}=185.05^{\circ}\pm0.04^{\circ}$, $\rm{Dec}=29.25^{\circ}\pm0.03^{\circ}$) is compatible with NGC 4278 within $\sim0.03$ degree. Variation analysis shows an indication of the variability at a few months level in the TeV band, which is consistent with low frequency observations. Based on these observations, we report the detection of TeV $γ$-ray emissions from this low-luminosity AGN NGC 4278. The observations by LHAASO-WCDA during active period has a significance level of 8.8\,$σ$ with best-fit photon spectral index $\varGamma=2.56\pm0.14$ and a flux $f_{1-10\,\rm{TeV}}=(7.0\pm1.1_{\rm{sta}}\pm0.35_{\rm{syst}})\times10^{-13}\,\rm{photons\,cm^{-2}\,s^{-1}}$, or approximately $5\%$ of the Crab Nebula. The discovery of VHE from NGC 4278 indicates that the compact, weak radio jet can efficiently accelerate particles and emit TeV photons. △ Less

Submitted 13 May, 2024; originally announced May 2024.

Comments: 11 pages, 5 figures

LHAASO-KM2A detector simulation using Geant4

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (254 additional authors not shown)

Abstract: KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with… ▽ More KM2A is one of the main sub-arrays of LHAASO, working on gamma ray astronomy and cosmic ray physics at energies above 10 TeV. Detector simulation is the important foundation for estimating detector performance and data analysis. It is a big challenge to simulate the KM2A detector in the framework of Geant4 due to the need to track numerous photons from a large number of detector units (>6000) with large altitude difference (30 m) and huge coverage (1.3 km^2). In this paper, the design of the KM2A simulation code G4KM2A based on Geant4 is introduced. The process of G4KM2A is optimized mainly in memory consumption to avoid memory overffow. Some simpliffcations are used to signiffcantly speed up the execution of G4KM2A. The running time is reduced by at least 30 times compared to full detector simulation. The particle distributions and the core/angle resolution comparison between simulation and experimental data of the full KM2A array are also presented, which show good agreement. △ Less

Submitted 7 April, 2024; originally announced April 2024.

Constraint on magnetized galactic outflows from LOFAR rotation measure data

Authors: J. Blunier, A. Neronov

Abstract: Outflows from galaxies, driven by active galactic nuclei and star formation activity, spread magnetic fields into the intergalactic medium. Importance of this process can be assessed using cosmological magneto-hydrodynamic numerical modelling of the baryonic feedback on the Large Scale Structure, like that of Illustris-TNG simulations. We use the Faraday Rotation Measure data of LOFAR Two-Meter Sk… ▽ More Outflows from galaxies, driven by active galactic nuclei and star formation activity, spread magnetic fields into the intergalactic medium. Importance of this process can be assessed using cosmological magneto-hydrodynamic numerical modelling of the baryonic feedback on the Large Scale Structure, like that of Illustris-TNG simulations. We use the Faraday Rotation Measure data of LOFAR Two-Meter Sky Survey (LOTSS) to test the Illustris-TNG baryonic feedback model. We show that the Illustris-TNG over-predicts the root-mean-square of the residual Rotation Measure in LOTSS. This suggests that the "pollution" of the intergalactic medium by the magnetized outflows from galaxies is less important than the estimate from the Illustris-TNG. This fact provides a support to the hypothesis that the volume-filling large scale magnetic fields found in the voids of the Large Scale Structure are likely of cosmological origin. △ Less

Submitted 20 March, 2024; originally announced March 2024.

Comments: 5 pages, 5 figures

Measurements of All-Particle Energy Spectrum and Mean Logarithmic Mass of Cosmic Rays from 0.3 to 30 PeV with LHAASO-KM2A

Authors: The LHAASO Collaboration, Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen , et al. (256 additional authors not shown)

Abstract: We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at… ▽ More We present the measurements of all-particle energy spectrum and mean logarithmic mass of cosmic rays in the energy range of 0.3-30 PeV using data collected from LHAASO-KM2A between September 2021 and December 2022, which is based on a nearly composition-independent energy reconstruction method, achieving unprecedented accuracy. Our analysis reveals the position of the knee at $3.67 \pm 0.05 \pm 0.15$ PeV. Below the knee, the spectral index is found to be -$2.7413 \pm 0.0004 \pm 0.0050$, while above the knee, it is -$3.128 \pm 0.005 \pm 0.027$, with the sharpness of the transition measured with a statistical error of 2%. The mean logarithmic mass of cosmic rays is almost heavier than helium in the whole measured energy range. It decreases from 1.7 at 0.3 PeV to 1.3 at 3 PeV, representing a 24% decline following a power law with an index of -$0.1200 \pm 0.0003 \pm 0.0341$. This is equivalent to an increase in abundance of light components. Above the knee, the mean logarithmic mass exhibits a power law trend towards heavier components, which is reversal to the behavior observed in the all-particle energy spectrum. Additionally, the knee position and the change in power-law index are approximately the same. These findings suggest that the knee observed in the all-particle spectrum corresponds to the knee of the light component, rather than the medium-heavy components. △ Less

Submitted 26 March, 2024; v1 submitted 15 March, 2024; originally announced March 2024.

Comments: 8 pages, 3 figures

Journal ref: Physical Review Letters 132, 131002 (2024)

Resistively controlled primordial magnetic turbulence decay

Authors: A. Brandenburg, A. Neronov, F. Vazza

Abstract: Magnetic fields generated in the early Universe undergo turbulent decay during the radiation-dominated era. The decay is governed by a decay exponent and a decay time. It has been argued that the latter is prolonged by magnetic reconnection, which depends on the microphysical resistivity and viscosity. Turbulence, on the other hand, is not usually expected to be sensitive to microphysical dissipat… ▽ More Magnetic fields generated in the early Universe undergo turbulent decay during the radiation-dominated era. The decay is governed by a decay exponent and a decay time. It has been argued that the latter is prolonged by magnetic reconnection, which depends on the microphysical resistivity and viscosity. Turbulence, on the other hand, is not usually expected to be sensitive to microphysical dissipation, which affects only very small scales. We want to test and quantify the reconnection hypothesis in decaying hydromagnetic turbulence. We performed high-resolution numerical simulations with zero net magnetic helicity using the Pencil Code with up to $2048^3$ mesh points and relate the decay time to the Alfvén time for different resistivities and viscosities. The decay time is found to be longer than the Alfvén time by a factor that increases with increasing Lundquist number to the 1/4 power. The decay exponent is as expected from the conservation of the Hosking integral, but a timescale dependence on resistivity is unusual for developed turbulence and not found for hydrodynamic turbulence. In two dimensions, the Lundquist number dependence is shown to be leveling off above values of $\approx25,000$, independently of the value of the viscosity. Our numerical results suggest that resistivity effects have been overestimated in earlier work. Instead of reconnection, it may be the magnetic helicity density in smaller patches that is responsible for the resistively slow decay. The leveling off at large Lundquist number cannot currently be confirmed in three dimensions. △ Less

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

Comments: 15 pages, 17 figures, 3 tables; A&A (in press)

Report number: NORDITA-2024-001

Journal ref: Astron. Astrophys. 687, A186 (2024)

EUSO-SPB1 Mission and Science

Authors: JEM-EUSO Collaboration, :, G. Abdellaoui, S. Abe, J. H. Adams. Jr., D. Allard, G. Alonso, L. Anchordoqui, A. Anzalone, E. Arnone, K. Asano, R. Attallah, H. Attoui, M. Ave Pernas, R. Bachmann, S. Bacholle, M. Bagheri, M. Bakiri, J. Baláz, D. Barghini, S. Bartocci, M. Battisti, J. Bayer, B. Beldjilali, T. Belenguer , et al. (271 additional authors not shown)

Abstract: The Extreme Universe Space Observatory on a Super Pressure Balloon 1 (EUSO-SPB1) was launched in 2017 April from Wanaka, New Zealand. The plan of this mission of opportunity on a NASA super pressure balloon test flight was to circle the southern hemisphere. The primary scientific goal was to make the first observations of ultra-high-energy cosmic-ray extensive air showers (EASs) by looking down on… ▽ More The Extreme Universe Space Observatory on a Super Pressure Balloon 1 (EUSO-SPB1) was launched in 2017 April from Wanaka, New Zealand. The plan of this mission of opportunity on a NASA super pressure balloon test flight was to circle the southern hemisphere. The primary scientific goal was to make the first observations of ultra-high-energy cosmic-ray extensive air showers (EASs) by looking down on the atmosphere with an ultraviolet (UV) fluorescence telescope from suborbital altitude (33~km). After 12~days and 4~hours aloft, the flight was terminated prematurely in the Pacific Ocean. Before the flight, the instrument was tested extensively in the West Desert of Utah, USA, with UV point sources and lasers. The test results indicated that the instrument had sensitivity to EASs of approximately 3 EeV. Simulations of the telescope system, telescope on time, and realized flight trajectory predicted an observation of about 1 event assuming clear sky conditions. The effects of high clouds were estimated to reduce this value by approximately a factor of 2. A manual search and a machine-learning-based search did not find any EAS signals in these data. Here we review the EUSO-SPB1 instrument and flight and the EAS search. △ Less

Submitted 12 January, 2024; originally announced January 2024.

Comments: 18 pages, 19 figures

Journal ref: Astropart Phys 154 (2024) 102891

JEM-EUSO Collaboration contributions to the 38th International Cosmic Ray Conference

Authors: S. Abe, J. H. Adams Jr., D. Allard, P. Alldredge, R. Aloisio, L. Anchordoqui, A. Anzalone, E. Arnone, M. Bagheri, B. Baret, D. Barghini, M. Battisti, R. Bellotti, A. A. Belov, M. Bertaina, P. F. Bertone, M. Bianciotto, F. Bisconti, C. Blaksley, S. Blin-Bondil, K. Bolmgren, S. Briz, J. Burton, F. Cafagna, G. Cambiè , et al. (133 additional authors not shown)

Abstract: This is a collection of papers presented by the JEM-EUSO Collaboration at the 38th International Cosmic Ray Conference (Nagoya, Japan, July 26-August 3, 2023) This is a collection of papers presented by the JEM-EUSO Collaboration at the 38th International Cosmic Ray Conference (Nagoya, Japan, July 26-August 3, 2023) △ Less

Submitted 13 December, 2023; originally announced December 2023.

Neutrino signal from Cygnus region of the Milky Way

Authors: A. Neronov, D. Semikoz, D. Savchenko

Abstract: Interactions of cosmic ray protons and nuclei in their sources and in the interstellar medium produce "hadronic" gamma-ray emission. Gamma-rays can also be of "leptonic" origin, i.e. originating from high-energy electrons accelerated together with protons. It is difficult to distinguish between hadronic and leptonic emission mechanisms based on gamma-ray data alone. This can be done via detection… ▽ More Interactions of cosmic ray protons and nuclei in their sources and in the interstellar medium produce "hadronic" gamma-ray emission. Gamma-rays can also be of "leptonic" origin, i.e. originating from high-energy electrons accelerated together with protons. It is difficult to distinguish between hadronic and leptonic emission mechanisms based on gamma-ray data alone. This can be done via detection of neutrinos, because only hadronic processes lead to neutrino production. We use publicly available ten-year IceCube neutrino telescope dataset to demonstrate the hadronic nature of high-energy emission from the direction of Cygnus region of the Milky Way. We find a 3-sigma excess of neutrino events from an extended Cygnus Cocoon, with the flux comparable to the flux of gamma-rays in the multi-TeV energy range seen by HAWC and LHAASO telescopes. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 4 pages 3 figures

Developments and results in the context of the JEM-EUSO program obtained with the ESAF Simulation and Analysis Framework

Authors: S. Abe, J. H. Adams Jr., D. Allard, P. Alldredge, L. Anchordoqui, A. Anzalone, E. Arnone, B. Baret, D. Barghini, M. Battisti, J. Bayer, R. Bellotti, A. A. Belov, M. Bertaina, P. F. Bertone, M. Bianciotto, P. L. Biermann, F. Bisconti, C. Blaksley, S. Blin-Bondil, P. Bobik, K. Bolmgren, S. Briz, J. Burton, F. Cafagna , et al. (150 additional authors not shown)

Abstract: JEM--EUSO is an international program for the development of space-based Ultra-High Energy Cosmic Ray observatories. The program consists of a series of missions which are either under development or in the data analysis phase. All instruments are based on a wide-field-of-view telescope, which operates in the near-UV range, designed to detect the fluorescence light emitted by extensive air showers… ▽ More JEM--EUSO is an international program for the development of space-based Ultra-High Energy Cosmic Ray observatories. The program consists of a series of missions which are either under development or in the data analysis phase. All instruments are based on a wide-field-of-view telescope, which operates in the near-UV range, designed to detect the fluorescence light emitted by extensive air showers in the atmosphere. We describe the simulation software ESAFin the framework of the JEM--EUSO program and explain the physical assumptions used. We present here the implementation of the JEM--EUSO, POEMMA, K--EUSO, TUS, Mini--EUSO, EUSO--SPB1 and EUSO--TA configurations in ESAF. For the first time ESAF simulation outputs are compared with experimental data. △ Less

Submitted 21 November, 2023; originally announced November 2023.

Journal ref: Eur. Phys. J. C 83, 1028 (2023)

Does or did the supernova remnant Cassiopeia A operate as a PeVatron?

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE;… ▽ More For decades, supernova remnants (SNRs) have been considered the prime sources of Galactic Cosmic rays (CRs). But whether SNRs can accelerate CR protons to PeV energies and thus dominate CR flux up to the knee is currently under intensive theoretical and phenomenological debate. The direct test of the ability of SNRs to operate as CR PeVatrons can be provided by ultrahigh-energy (UHE; $E_γ\geq 100$~TeV) $γ$-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising target for UHE observations. This paper presents the observation of Cas A and its vicinity by the LHAASO KM2A detector. The exceptional sensitivity of LHAASO KM2A in the UHE band, combined with the young age of Cas A, enabled us to derive stringent model-independent limits on the energy budget of UHE protons and nuclei accelerated by Cas A at any epoch after the explosion. The results challenge the prevailing paradigm that Cas A-type SNRs are major suppliers of PeV CRs in the Milky Way. △ Less

Submitted 25 October, 2023; originally announced October 2023.

Comments: 11 pages, 3 figures, Accepted by the APJL

Very high energy gamma-ray emission beyond 10 TeV from GRB 221009A

Authors: Zhen Cao, F. Aharonian, Q. An, A. Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the t… ▽ More The highest energy gamma-rays from gamma-ray bursts (GRBs) have important implications for their radiation mechanism. Here we report for the first time the detection of gamma-rays up to 13 TeV from the brightest GRB 221009A by the Large High Altitude Air-shower Observatory (LHAASO). The LHAASO-KM2A detector registered more than 140 gamma-rays with energies above 3 TeV during 230$-$900s after the trigger. The intrinsic energy spectrum of gamma-rays can be described by a power-law after correcting for extragalactic background light (EBL) absorption. Such a hard spectrum challenges the synchrotron self-Compton (SSC) scenario of relativistic electrons for the afterglow emission above several TeV. Observations of gamma-rays up to 13 TeV from a source with a measured redshift of z=0.151 hints more transparency in intergalactic space than previously expected. Alternatively, one may invoke new physics such as Lorentz Invariance Violation (LIV) or an axion origin of very high energy (VHE) signals. △ Less

Submitted 22 November, 2023; v1 submitted 13 October, 2023; originally announced October 2023.

Comments: 49pages, 11figures

Journal ref: Science Advances, 9, eadj2778 (2023) 15 November 2023

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)

LISA and $γ$-ray telescopes as multi-messenger probes of a first-order cosmological phase transition

Authors: A. Roper Pol, A. Neronov, C. Caprini, T. Boyer, D. Semikoz

Abstract: We study two possible cosmological consequences of a first-order phase transition in the temperature range of $1$ GeV to $10^3$ TeV: the generation of a stochastic gravitational wave background (SGWB) within the sensitivity of the Laser Interferometer Space Antenna (LISA) and, simultaneously, primordial magnetic fields that would evolve through the Universe's history and could be compatible with t… ▽ More We study two possible cosmological consequences of a first-order phase transition in the temperature range of $1$ GeV to $10^3$ TeV: the generation of a stochastic gravitational wave background (SGWB) within the sensitivity of the Laser Interferometer Space Antenna (LISA) and, simultaneously, primordial magnetic fields that would evolve through the Universe's history and could be compatible with the lower bound from $γ$-ray telescopes on intergalactic magnetic fields (IGMF) at present time. We find that, if even a small fraction of the kinetic energy in sound waves is converted into MHD turbulence, a first order phase transition occurring at temperature between $1$ and $10^6$ GeV can give rise to an observable SGWB signal in LISA and, at the same time, an IGMF compatible with the lower bound from the $γ$-ray telescope MAGIC, for all proposed evolutionary paths of the magnetic fields throughout the radiation dominated era. For two values of the fraction of energy density converted into turbulence, $\varepsilon_{\rm turb}=0.1$ and $1$, we provide the range of first-order phase transition parameters (strength $α$, duration $β^{-1}$, bubbles wall speed $v_w$, and temperature $T_*$), together with the corresponding range of magnetic field strength $B$ and correlation length $λ$, that would lead to the SGWB and IGMF observable with LISA and MAGIC. The resulting magnetic field strength at recombination can also correspond to the one that has been proposed to induce baryon clumping, previously suggested as a possible way to ease the Hubble tension. In the limiting case $\varepsilon_{\rm turb} \ll 1$, the SGWB is only sourced by sound waves, however, an IGMF is still generated. We find that values as small as $\varepsilon_{\rm turb} \sim O(10^{-13})$ (helical) and $O (10^{-9})$ (non-helical) can provide IGMF compatible with MAGIC's lower bound. △ Less

Submitted 20 July, 2023; originally announced July 2023.

Comments: 10 pages, 4 figures

Hadronic nature of high-energy emission from the Galactic Ridge

Authors: A. Neronov, D. Semikoz, J. Aublin, M. Lamoureux, A. Kouchner

Abstract: We show that the IceCube observation of the Galactic neutrino flux component confirms the hint of detection of neutrinos from the Galactic Ridge (the inner part of the Milky Way disk within the Galactic longitude |l|<30 degrees), previously reported by the ANTARES collaboration. This confirmation indicates that the bulk of the high-energy flux from the Galactic Ridge in multi-TeV band is produced… ▽ More We show that the IceCube observation of the Galactic neutrino flux component confirms the hint of detection of neutrinos from the Galactic Ridge (the inner part of the Milky Way disk within the Galactic longitude |l|<30 degrees), previously reported by the ANTARES collaboration. This confirmation indicates that the bulk of the high-energy flux from the Galactic Ridge in multi-TeV band is produced by interactions of high-energy protons and atomic nuclei, rather than electrons. We show that both ANTARES and IceCube measurements agree with the Fermi-LAT telescope measurements of the gamma-ray emission from the Ridge. The multi-messenger (neutrino plus gamma-ray) spectrum of the Ridge over a broad energy range from 10 GeV to 10 TeV is consistent with a model of pion decay emission produced by a power-law distribution of protons with a slope Gamma~2.5, harder than that of the locally observed cosmic ray spectrum. This provides for the first time an unambiguous multi-messenger demonstration of the variability of the spectrum of cosmic rays across the Galactic disk. △ Less

Submitted 16 July, 2023; originally announced July 2023.

Comments: 4 pages, 1 figure

The second data release from the European Pulsar Timing Array: IV. Implications for massive black holes, dark matter and the early Universe

Authors: J. Antoniadis, P. Arumugam, S. Arumugam, P. Auclair, S. Babak, M. Bagchi, A. -S. Bak Nielsen, E. Barausse, C. G. Bassa, A. Bathula, A. Berthereau, M. Bonetti, E. Bortolas, P. R. Brook, M. Burgay, R. N. Caballero, C. Caprini, A. Chalumeau, D. J. Champion, S. Chanlaridis, S. Chen, I. Cognard, M. Crisostomi, S. Dandapat, D. Deb , et al. (89 additional authors not shown)

Abstract: The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling sup… ▽ More The European Pulsar Timing Array (EPTA) and Indian Pulsar Timing Array (InPTA) collaborations have measured a low-frequency common signal in the combination of their second and first data releases respectively, with the correlation properties of a gravitational wave background (GWB). Such signal may have its origin in a number of physical processes including a cosmic population of inspiralling supermassive black hole binaries (SMBHBs); inflation, phase transitions, cosmic strings and tensor mode generation by non-linear evolution of scalar perturbations in the early Universe; oscillations of the Galactic potential in the presence of ultra-light dark matter (ULDM). At the current stage of emerging evidence, it is impossible to discriminate among the different origins. Therefore, in this paper, we consider each process separately, and investigate the implications of the signal under the hypothesis that it is generated by that specific process. We find that the signal is consistent with a cosmic population of inspiralling SMBHBs, and its relatively high amplitude can be used to place constraints on binary merger timescales and the SMBH-host galaxy scaling relations. If this origin is confirmed, this is the first direct evidence that SMBHBs merge in nature, adding an important observational piece to the puzzle of structure formation and galaxy evolution. As for early Universe processes, the measurement would place tight constraints on the cosmic string tension and on the level of turbulence developed by first-order phase transitions. Other processes would require non-standard scenarios, such as a blue-tilted inflationary spectrum or an excess in the primordial spectrum of scalar perturbations at large wavenumbers. Finally, a ULDM origin of the detected signal is disfavoured, which leads to direct constraints on the abundance of ULDM in our Galaxy. △ Less

Submitted 15 May, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 23 figures, replaced to match the version published in Astronomy & Astrophysics, note the change in the numbering order in the series (now paper IV)

Measurement of size of gamma-ray source in blazar B0218+357 from microlensing at 100 GeV energy

Authors: Ie. Vovk, A. Neronov

Abstract: Observations of the effect of microlensing in gravitationally lensed quasars can be used to study the structure of active galactic nuclei on distance scales down to the sizes of the supermassive black holes powering source activity. We search for the microlensing in the signal from a gravitationally lensed blazar B0218+357 in very-high-energy gamma-ray band. We combine observations of a bright fla… ▽ More Observations of the effect of microlensing in gravitationally lensed quasars can be used to study the structure of active galactic nuclei on distance scales down to the sizes of the supermassive black holes powering source activity. We search for the microlensing in the signal from a gravitationally lensed blazar B0218+357 in very-high-energy gamma-ray band. We combine observations of a bright flare of the source in 2014 with Fermi/LAT and MAGIC telescopes in 0.1-100 GeV energy range. Using the time-delayed leading and trailing signals from two gravitationally lensed images of the source, we measure magnification factor at the moment of the flare. We use the scaling of the maximal magnification factor with the source size to constrain the size of gamma-ray emission region in the entire 0.1-100 GeV energy range. The magnification factor in the very-high-energy band derived from a comparison of Fermi/LAT and MAGIC data is $μ_{VHE} = 36^{+40}_{-26}$, substantially larger than that in the radio band. This suggests that one of the source images is strongly affected by microlensing at the moment of the flare. Assuming that the microlensing is produced by a stellar mass object in the lens galaxy, we constrain the size of the emission region in the $E>100$ GeV band to be $\mathrm{R_{VHE} = 2.2^{+27}_{-1.7} \times 10^{13}~cm}$. We note that the spectrum of the microlensed source was unusually hard at the moment of the flare and speculate that this hardening may be due to the energy dependent microlensing effect. This interpretation suggests that the source size decreases with energy in 0.1-100 GeV energy range studied. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: 4 pages, 2 figures

Neutrino signal from Seyfert galaxies

Authors: A. Neronov, D. Savchenko, D. V. Semikoz

Abstract: IceCube collaboration has previously reported an evidence for neutrino signal from a Seyfert galaxy NGC 1068. This may suggest that all Seyfert galaxies emit neutrinos. To test this hypothesis, we identify the best candidate neutrino sources among nearby Seyfert galaxies, based on their hard X-ray properties. Only two other sources, NGC 4151 and NGC 3079 are expected to be detectable in 10 years o… ▽ More IceCube collaboration has previously reported an evidence for neutrino signal from a Seyfert galaxy NGC 1068. This may suggest that all Seyfert galaxies emit neutrinos. To test this hypothesis, we identify the best candidate neutrino sources among nearby Seyfert galaxies, based on their hard X-ray properties. Only two other sources, NGC 4151 and NGC 3079 are expected to be detectable in 10 years of IceCube data. We find an evidence for neutrino signal from both sources in publicly available ten-year IceCube dataset. The chance coincidence probability to find the observed neutrino count excesses in the directions of the two out of two expected sources, in addition to the previously reported brightest source, is p<2.6e-7. This corresponds to the detection of Seyfert galaxies as a neutrino source class. △ Less

Submitted 15 June, 2023; originally announced June 2023.

Comments: 4 pages

Constraint on intergalactic magnetic field from Fermi/LAT observations of the "pair echo" of GRB 221009A

Authors: Ie. Vovk, A. Korochkin, A. Neronov, D. Semikoz

Abstract: Delayed "pair echo" signal from interactions of very-high-energy gamma rays in the intergalactic medium can be used for detection of the inter-galactic magnetic field (IGMF). We use the data of Fermi/LAT telescope coupled with LHAASO observatory measurements to confirm the presence of IGMF along the line of sight to the gamma-ray burst GRB221009A. Comparing the Fermi/LAT measurements with the expe… ▽ More Delayed "pair echo" signal from interactions of very-high-energy gamma rays in the intergalactic medium can be used for detection of the inter-galactic magnetic field (IGMF). We use the data of Fermi/LAT telescope coupled with LHAASO observatory measurements to confirm the presence of IGMF along the line of sight to the gamma-ray burst GRB221009A. Comparing the Fermi/LAT measurements with the expected level of the pair echo flux, set by the multi-TeV LHAASO detection, we derive a lower bound $10^{-19}$ G on the IGMF with correlation length $l$ larger than 1 Mpc, improving as $l^{-1/2}$ for shorter correlation lengths. This provides an independent verification of existence of a lower bound on IGMF in the voids of the Large Scale Structure, previously derived from the observations of active galactic nuclei. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 4 pages, 2 figues

Energy-dependent periodicities of LS I +61$^\circ$ 303 in the GeV band

Authors: M. Chernyakova, D. Malyshev, A. Neronov, D. Savchenko

Abstract: LS I +61$^\circ$303 is a rare representative of the gamma-ray binaries with a compact object known to be a pulsar. We report on the periodicity and spectral analysis of this source performed with more than 14 years of Fermi/LAT data. The periodicity of LS I +61$^\circ$303 is strongly energy dependent. Two periods $P_1 = 26.932\pm 0.004 (stat)\pm 0.008 (syst)$ and… ▽ More LS I +61$^\circ$303 is a rare representative of the gamma-ray binaries with a compact object known to be a pulsar. We report on the periodicity and spectral analysis of this source performed with more than 14 years of Fermi/LAT data. The periodicity of LS I +61$^\circ$303 is strongly energy dependent. Two periods $P_1 = 26.932\pm 0.004 (stat)\pm 0.008 (syst)$ and $P_2 = 26.485 \pm 0.004 (stat)\pm 0.007 (syst)$ are detected only at $E>1$ GeV and at $E<0.3$ GeV correspondingly. Within $1σ$ (stat+syst) the periods are consistent with orbital ($P_2$) and beat orbital/superorbital ($P_1$) periods. We present the orbital light curves of the system in several energy bands and the results of the spectral analysis. We discuss the possible origin of the change in the variability pattern between 0.1 and 1 GeV energy. △ Less

Submitted 1 August, 2023; v1 submitted 9 June, 2023; originally announced June 2023.

Comments: accepted to MNRAS

The First LHAASO Catalog of Gamma-Ray Sources

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022.… ▽ More We present the first catalog of very-high energy and ultra-high energy gamma-ray sources detected by the Large High Altitude Air Shower Observatory (LHAASO). The catalog was compiled using 508 days of data collected by the Water Cherenkov Detector Array (WCDA) from March 2021 to September 2022 and 933 days of data recorded by the Kilometer Squared Array (KM2A) from January 2020 to September 2022. This catalog represents the main result from the most sensitive large coverage gamma-ray survey of the sky above 1 TeV, covering declination from $-$20$^{\circ}$ to 80$^{\circ}$. In total, the catalog contains 90 sources with an extended size smaller than $2^\circ$ and a significance of detection at $> 5σ$. Based on our source association criteria, 32 new TeV sources are proposed in this study. Among the 90 sources, 43 sources are detected with ultra-high energy ($E > 100$ TeV) emission at $> 4σ$ significance level. We provide the position, extension, and spectral characteristics of all the sources in this catalog. △ Less

Submitted 27 November, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 40 pages, 13 figures, 4 tables

Journal ref: The Astrophysical Journal Supplement Series, 271 (2024) 25

Measurement of ultra-high-energy diffuse gamma-ray emission of the Galactic plane from 10 TeV to 1 PeV with LHAASO-KM2A

Authors: Zhen Cao, F. Aharonian, Q. An, Axikegu, Y. X. Bai, Y. W. Bao, D. Bastieri, X. J. Bi, Y. J. Bi, J. T. Cai, Q. Cao, W. Y. Cao, Zhe Cao, J. Chang, J. F. Chang, A. M. Chen, E. S. Chen, Liang Chen, Lin Chen, Long Chen, M. J. Chen, M. L. Chen, Q. H. Chen, S. H. Chen, S. Z. Chen , et al. (255 additional authors not shown)

Abstract: The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer ar… ▽ More The diffuse Galactic $γ$-ray emission, mainly produced via interactions between cosmic rays and the interstellar medium and/or radiation field, is a very important probe of the distribution, propagation, and interaction of cosmic rays in the Milky Way. In this work we report the measurements of diffuse $γ$-rays from the Galactic plane between 10 TeV and 1 PeV energies, with the square kilometer array of the Large High Altitude Air Shower Observatory (LHAASO). Diffuse emissions from the inner ($15^{\circ}<l<125^{\circ}$, $|b|<5^{\circ}$) and outer ($125^{\circ}<l<235^{\circ}$, $|b|<5^{\circ}$) Galactic plane are detected with $29.1σ$ and $12.7σ$ significance, respectively. The outer Galactic plane diffuse emission is detected for the first time in the very- to ultra-high-energy domain ($E>10$~TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of $-2.99\pm0.04$, which is different from the curved spectrum as expected from hadronic interactions between locally measured cosmic rays and the line-of-sight integrated gas content. Furthermore, the measured flux is higher by a factor of $\sim3$ than the prediction. A similar spectrum with an index of $-2.99\pm0.07$ is found in the outer Galaxy region, and the absolute flux for $10\lesssim E\lesssim60$ TeV is again higher than the prediction for hadronic cosmic ray interactions. The latitude distributions of the diffuse emission are consistent with the gas distribution, while the longitude distributions show clear deviation from the gas distribution. The LHAASO measurements imply that either additional emission sources exist or cosmic ray intensities have spatial variations. △ Less

Submitted 19 August, 2023; v1 submitted 9 May, 2023; originally announced May 2023.

Comments: 12 pages, 8 figures, 5 tables; accepted for publication in Physical Review Letters; source mask file provided as ancillary file

Journal ref: Phys. Rev. Lett. 131, 151001 (2023)

Hint for a TeV neutrino emission from the Galactic Ridge with ANTARES

Authors: A. Albert, S. Alves, M. André, M. Ardid, S. Ardid, J. -J. Aubert, J. Aublin, B. Baret, S. Basa, Y. Becherini, B. Belhorma, M. Bendahman, F. Benfenati, V. Bertin, S. Biagi, M. Bissinger, J. Boumaaza, M. Bouta, M. C. Bouwhuis, H. Brânzaş, R. Bruijn, J. Brunner, J. Busto, B. Caiffi, D. Calvo , et al. (129 additional authors not shown)

Abstract: Interactions of cosmic ray protons, atomic nuclei, and electrons in the interstellar medium in the inner part of the Milky Way produce a $γ$-ray flux from the Galactic Ridge. If the $γ$-ray emission is dominated by proton and nuclei interactions, a neutrino flux comparable to the $γ$-ray flux is expected from the same sky region. Data collected by the ANTARES neutrino telescope are used to constra… ▽ More Interactions of cosmic ray protons, atomic nuclei, and electrons in the interstellar medium in the inner part of the Milky Way produce a $γ$-ray flux from the Galactic Ridge. If the $γ$-ray emission is dominated by proton and nuclei interactions, a neutrino flux comparable to the $γ$-ray flux is expected from the same sky region. Data collected by the ANTARES neutrino telescope are used to constrain the neutrino flux from the Galactic Ridge in the 1-100 TeV energy range. Neutrino events reconstructed both as tracks and showers are considered in the analysis and the selection is optimized for the search of an excess in the region $|l| < 30°$, $|b| < 2°$. The expected background in the search region is estimated using an off-zone region with similar sky coverage. Neutrino signal originating from a power-law spectrum with spectral index ranging from $Γ_ν=1$ to $4$ is simulated in both channels. The observed energy distributions are fitted to constrain the neutrino emission from the Ridge. The energy distributions in the signal region are inconsistent with the background expectation at $\sim 96\%$ confidence level. The mild excess over the background is consistent with a neutrino flux with a power law with a spectral index $2.45^{+0.22}_{-0.34}$ and a flux normalization $dN_ν/dE_ν= 4.0^{+2.7}_{-2.0} \times 10^{-16} \text{GeV}^{-1} \text{cm}^{-2} \text{s}^{-1} \text{sr}^{-1}$ at 40 TeV reference energy. Such flux is consistent with the expected neutrino signal if the bulk of the observed $γ$-ray flux from the Galactic Ridge originates from interactions of cosmic ray protons and nuclei with a power-law spectrum extending well into the PeV energy range. △ Less

Submitted 12 May, 2023; v1 submitted 22 December, 2022; originally announced December 2022.

Comments: 8 pages, 5 figures

Journal ref: Physics Letters B, Volume 841, 2023, 137951

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)

The contribution of magnetized galactic outflows to extragalactic Faraday rotation

Authors: Andres Aramburo-Garcia, Kyrylo Bondarenko, Alexey Boyarsky, Andrii Neronov, Anna Scaife, Anastasia Sokolenko

Abstract: Galactic outflows driven by star formation and active galactic nuclei blow bubbles into their local environments, causing galactic magnetic fields to be carried into intergalactic space. We explore the redshift-dependent effect of these magnetized bubbles on the Faraday Rotation Measure (RM) of extragalactic radio sources. Using the IllustrisTNG cosmological simulations, we separate the contributi… ▽ More Galactic outflows driven by star formation and active galactic nuclei blow bubbles into their local environments, causing galactic magnetic fields to be carried into intergalactic space. We explore the redshift-dependent effect of these magnetized bubbles on the Faraday Rotation Measure (RM) of extragalactic radio sources. Using the IllustrisTNG cosmological simulations, we separate the contribution from magnetic bubbles from that of the volume-filling magnetic component expected to be due to the seed field originating in the Early Universe. We use this separation to extract the redshift dependence of each component and to compare TNG model predictions with observation measurements of the NRAO VLA Sky Survey (NVSS). We find that magnetized bubbles provide a sizeable contribution to the extragalactic RM, with redshift-independent $\langle |{\rm RM}| \rangle \simeq 13$ rad/m$^2$ for sources at redshifts $z\ge 2$. This is close to the mean residual RM of $16$ rad/m$^2$ found from NVSS data in this redshift range. Using the IllustrisTNG simulations, we also evaluate a simple model for the contribution to residual RM from individual host galaxies and show that this contribution is negligible at high-redshift. While the contribution from magnetic bubbles in the IllustrisTNG model is currently compatible with observational measurements of residual RM, the next-generation RM sky surveys, which will be free from the wrapping uncertainty, have larger statistics and better sensitivity should be able to observe predicted flat contribution from magnetic bubbles at large redshifts. This should allow to experimentally probe magnetic bubbles and check models of galaxy feedback in cosmological simulations. △ Less

Submitted 13 April, 2022; originally announced April 2022.

Comments: 6 pages, 3 figures

Revision of Faraday rotation measure constraints on the primordial magnetic field using the IllustrisTNG simulation

Authors: Andres Aramburo-Garcia, Kyrylo Bondarenko, Alexey Boyarsky, Andrii Neronov, Anna Scaife, Anastasia Sokolenko

Abstract: Previously derived Faraday rotation constraints on the volume-filling intergalactic magnetic field (IGMF) have used analytic models that made a range of simplifying assumptions about magnetic field evolution in the intergalactic medium and did not consider the effect of baryonic feedback on large-scale structures. In this work we revise existing Faraday rotation constraints on the IGMF using a num… ▽ More Previously derived Faraday rotation constraints on the volume-filling intergalactic magnetic field (IGMF) have used analytic models that made a range of simplifying assumptions about magnetic field evolution in the intergalactic medium and did not consider the effect of baryonic feedback on large-scale structures. In this work we revise existing Faraday rotation constraints on the IGMF using a numerical model of the intergalactic medium from the IllustrisTNG cosmological simulation that includes a sophisticated model of the baryonic feedback. We use the IllustrisTNG model to calculate the rotation measure and compare the resulting mean and median of the absolute value of the rotation measure with data from the NRAO VLA Sky Survey (NVSS). The numerical model of the intergalactic medium includes a full magneto-hydrodynamic model of the compressed primordial magnetic field as well as a model of the regions where the magnetic field is not primordial, but is rather produced by the process of baryonic feedback. Separating these two types of regions, we are able to assess the influence of the primordial magnetic field on the Faraday rotation signal. We find that by correcting for regions of compressed primordial field and accounting for the fact that part of the intergalactic medium is occupied by magnetic fields spread by baryonic feedback processes rather than by the primordial field relaxes the Faraday rotation bound by a factor of $\simeq 3$. This results in $B_0<1.8\times10^{-9}$ G for large correlation length IGMFs. △ Less

Submitted 7 September, 2022; v1 submitted 12 April, 2022; originally announced April 2022.

Comments: Replaced with the published version

Report number: FERMILAB-PUB-22-299-V

Signatures of anisotropic diffusion around PeVatrons in 100 TeV gamma-ray data

Authors: G. Giacinti, T. Abounnasr, A. Neronov, D. Semikoz

Abstract: The Tibet AS$γ$ collaboration has reported a diffuse $γ$-ray emission signal from the Galactic Plane. We consider the possibility that the diffuse emission from the outer Galactic Plane at the highest energies is produced by cosmic rays spreading from a single supernova-type source either in the Local or Perseus arm of the Milky Way. We show that anisotropic diffusion of multi-PeV cosmic rays alon… ▽ More The Tibet AS$γ$ collaboration has reported a diffuse $γ$-ray emission signal from the Galactic Plane. We consider the possibility that the diffuse emission from the outer Galactic Plane at the highest energies is produced by cosmic rays spreading from a single supernova-type source either in the Local or Perseus arm of the Milky Way. We show that anisotropic diffusion of multi-PeV cosmic rays along the Galactic magnetic field can produce an extended source spanning ten(s) of degrees on the sky, with a flux-per-unit-solid-angle consistent with Tibet AS$γ$ measurements. Observations of this new type of very extended sources, and measurements of their morphologies, can be used to characterize the anisotropic diffusion of PeV cosmic rays in the Galactic magnetic field, and to constrain the locations and properties of past PeVatrons. △ Less

Submitted 21 March, 2022; originally announced March 2022.

Comments: 7 pages, 3 figures

Journal ref: Phys. Rev. D 106, 123029 (2022)

Status of the K-EUSO Orbital Detector of Ultra-high Energy Cosmic Rays

Authors: P. Klimov, M. Battisti, A. Belov, M. Bertaina, M. Bianciotto, S. Blin-Bondil, M. Casolino, T. Ebisuzaki, F. Fenu, C. Fuglesang, W. Marszał, A. Neronov, E. Parizot, P. Picozza, Z. Plebaniak, G. Prévôt, M. Przybylak N. Sakaki, S. Sharakin, K. Shinozaki, J. Szabelski, Y. Takizawa, D. Trofimov, I. Yashin, M. Zotov

Abstract: K-EUSO (KLYPVE-EUSO) is a planned orbital mission aimed at studying ultra-high energy cosmic rays (UHECRs) by detecting fluorescence and Cherenkov light emitted by extensive air showers in the nocturnal atmosphere of Earth in the ultraviolet (UV) range. The observatory is being developed within the JEM-EUSO collaboration and is planned to be deployed on the International Space Station after 2025 a… ▽ More K-EUSO (KLYPVE-EUSO) is a planned orbital mission aimed at studying ultra-high energy cosmic rays (UHECRs) by detecting fluorescence and Cherenkov light emitted by extensive air showers in the nocturnal atmosphere of Earth in the ultraviolet (UV) range. The observatory is being developed within the JEM-EUSO collaboration and is planned to be deployed on the International Space Station after 2025 and operated for at least two years. The telescope, consisting of $\sim10^{5}$ independent pixels, will allow a spatial resolution of $\sim0.6$ km on the ground, and, from a 400 km altitude, it will achieve a large and full sky exposure to sample the highest energy range of the UHECR spectrum. We provide a comprehensive review of the current status of the development of the K-EUSO experiment, paying special attention to its hardware parts and expected performance. We demonstrate how results of the K-EUSO mission can complement the achievements of the existing ground-based experiments and push forward the intriguing studies of ultra-high energy cosmic rays, as well as bring new knowledge about other phenomena manifesting themselves in the atmosphere in the UV range. △ Less

Submitted 30 January, 2022; originally announced January 2022.

Comments: 26 pages

Journal ref: Universe 8(2), 88 (2022)

JEM-EUSO Collaboration contributions to the 37th International Cosmic Ray Conference

Authors: G. Abdellaoui, S. Abe, J. H. Adams Jr., D. Allard, G. Alonso, L. Anchordoqui, A. Anzalone, E. Arnone, K. Asano, R. Attallah, H. Attoui, M. Ave Pernas, M. Bagheri, J. Baláz, M. Bakiri, D. Barghini, S. Bartocci, M. Battisti, J. Bayer, B. Beldjilali, T. Belenguer, N. Belkhalfa, R. Bellotti, A. A. Belov, K. Benmessai , et al. (267 additional authors not shown)

Abstract: Compilation of papers presented by the JEM-EUSO Collaboration at the 37th International Cosmic Ray Conference (ICRC), held on July 12-23, 2021 (online) in Berlin, Germany. Compilation of papers presented by the JEM-EUSO Collaboration at the 37th International Cosmic Ray Conference (ICRC), held on July 12-23, 2021 (online) in Berlin, Germany. △ Less

Submitted 28 January, 2022; originally announced January 2022.

Comments: html page with links to the JEM-EUSO Collaboration papers presented at ICRC-2021, Berlin, Germany

Gravitational wave signal from primordial magnetic fields in the Pulsar Timing Array frequency band

Authors: Alberto Roper Pol, Chiara Caprini, Andrii Neronov, Dmitri Semikoz

Abstract: The NANOGrav, Parkes, European, and International Pulsar Timing Array (PTA) Collaborations have reported evidence for a common-spectrum process that can potentially correspond to a stochastic gravitational wave background (SGWB) in the 1--100 nHz frequency range. We consider the scenario in which this signal is produced by magnetohydrodynamic (MHD) turbulence in the early Universe, induced by a no… ▽ More The NANOGrav, Parkes, European, and International Pulsar Timing Array (PTA) Collaborations have reported evidence for a common-spectrum process that can potentially correspond to a stochastic gravitational wave background (SGWB) in the 1--100 nHz frequency range. We consider the scenario in which this signal is produced by magnetohydrodynamic (MHD) turbulence in the early Universe, induced by a nonhelical primordial magnetic field at the energy scale corresponding to the quark confinement phase transition. We perform MHD simulations to study the dynamical evolution of the magnetic field and compute the resulting SGWB. We show that the SGWB output from the simulations can be very well approximated by assuming that the magnetic anisotropic stress is constant in time, over a time interval related to the eddy turnover time. The analytical spectrum that we derive under this assumption features a change of slope at a frequency corresponding to the GW source duration that we confirm with the numerical simulations. We compare the SGWB signal with the PTA data to constrain the temperature scale at which the SGWB is sourced, as well as the amplitude and characteristic scale of the initial magnetic field. We find that the generation temperature is constrained to be in the 1--200 MeV range, the magnetic field amplitude must be $>1$\% of the radiation energy density at that time, and the magnetic field characteristic scale is constrained to be $>10$\% of the horizon scale. We show that the turbulent decay of this magnetic field will lead to a field at recombination that can help to alleviate the Hubble tension and can be tested by measurements in the voids of the Large Scale Structure with gamma-ray telescopes like the Cherenkov Telescope Array. △ Less

Submitted 11 May, 2022; v1 submitted 14 January, 2022; originally announced January 2022.

Comments: 28 pages, 9 figures, 1 table, accepted for publication in Phys. Rev. D, files and code available in https://github.com/AlbertoRoper/GW_turbulence

Modelling the propagation of very-high-energy gamma rays with the CRbeam code: Comparison with CRPropa and ELMAG codes

Authors: O. Kalashev, A. Korochkin, A. Neronov, D. Semikoz

Abstract: Very-high-energy gamma rays produce electron positron pairs in interactions with low-energy photons of extragalactic background light during propagation through the intergalactic medium. The electron-positron pairs generate secondary gamma rays detectable by gamma-ray telescopes. This secondary emission can be used to detect intergalactic magnetic fields (IGMF) in the voids of large-scale structur… ▽ More Very-high-energy gamma rays produce electron positron pairs in interactions with low-energy photons of extragalactic background light during propagation through the intergalactic medium. The electron-positron pairs generate secondary gamma rays detectable by gamma-ray telescopes. This secondary emission can be used to detect intergalactic magnetic fields (IGMF) in the voids of large-scale structure. A new gamma-ray observatory, namely, Cherenkov Telescope Array (CTA), will provide an increase in sensitivity for detections of these secondary gamma-ray emission and enable the measurement of its properties for sources at cosmological distances. The interpretation of the CTA data, including detection of IGMF and study of its properties and origins, will require precision modeling of the primary and secondary gamma-ray fluxes. We asses the precision of the modeling of the secondary gamma-ray emission using model calculations with publicly available Monte-Carlo codes CRPropa and ELMAG and compare their predictions with theoretical expectations and with model calculations of a newly developed CRbeam code. We find that model predictions of different codes differ by up to 50% for low-redshift sources, with discrepancies increasing up to order-of-magnitude level with the increasing source redshifts. We identify the origin of these discrepancies and demonstrate that after eliminating the inaccuracies found, the discrepancies between the three codes are reduced to 10% when modeling nearby sources with z~0.1. We argue that the new CRbeam code provides reliable predictions for spectral, timing and imaging properties of the secondary gamma-ray signal for both nearby and distant sources with z~1. Thus, it can be used to study gamma-ray sources and IGMF with a level of precision that is appropriate for the prospective CTA study of the effects of gamma-ray propagation through the intergalactic medium. △ Less

Submitted 12 April, 2023; v1 submitted 11 January, 2022; originally announced January 2022.

Comments: V2: 15 pages, 12 figures, version accepted by A&A. While the paper was under review, new versions of the codes CRPropa3-3.2 and ELMAG 3.03 were released in which many of the comments presented in this paper were taken into account (see text for the details)

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

Limit on intergalactic magnetic field from ultra-high-energy cosmic ray hotspot in Perseus-Pisces region

Authors: Andrii Neronov, Dmitri Semikoz, Oleg Kalashev

Abstract: Telescope Array collaboration has reported an evidence for existence of a source of ultra-high-energy cosmic ray (UHECR) events in Perseus-Pisces supercluster. We show that the mere existence of such a source imposes an upper bound on the strength of intergalactic magnetic field (IGMF) in the Taurus void lying between the Perseus-Pisces supercluster and the Milky Way galaxy. This limit is at the l… ▽ More Telescope Array collaboration has reported an evidence for existence of a source of ultra-high-energy cosmic ray (UHECR) events in Perseus-Pisces supercluster. We show that the mere existence of such a source imposes an upper bound on the strength of intergalactic magnetic field (IGMF) in the Taurus void lying between the Perseus-Pisces supercluster and the Milky Way galaxy. This limit is at the level of 10^{-10} G for a field with correlation length larger than the distance of the super-cluster ~ 70 Mpc. This bound is an order-of-magnitude stronger that the previously known bound on IGMF from radio Faraday rotation measurements and it is the first upper bound on magnetic field in the voids of the Large Scale Structure. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 5 pages, 3 figures

Report number: INR-TH-2021-027

Detectability of large correlation length inflationary magnetic field with Cherenkov telescopes

Authors: Alexander Korochkin, Andrii Neronov, Guilhem Lavaux, Marius Ramsoy, Dmitri Semikoz

Abstract: Magnetic fields occupying the voids of the large scale structure may be a relic from the Early Universe originating from either Inflation or from cosmological phase transitions. We explore the possibility of identifying the inflationary origin of the void magnetic fields and measuring its parameters with gamma-ray astronomy methods. The large correlation length inflationary field is expected to im… ▽ More Magnetic fields occupying the voids of the large scale structure may be a relic from the Early Universe originating from either Inflation or from cosmological phase transitions. We explore the possibility of identifying the inflationary origin of the void magnetic fields and measuring its parameters with gamma-ray astronomy methods. The large correlation length inflationary field is expected to impose a characteristic asymmetry of extended gamma-ray emission that is correlated between different sources on the sky. We show that a set of nearby blazars for which the extended emission is observable in the 0.1-1 TeV band with CTA can be used for the test of inflationary origin of the void magnetic fields. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: 8 pages, 6 figures, accepted to JETP

Energy spectra of secondaries in proton-proton interactions

Authors: S. Koldobskiy, M. Kachelrieß, A. Lskavyan, A. Neronov, S. Ostapchenko, D. V. Semikoz

Abstract: We compare the predictions of AAfrag for the spectra of secondary photons, neutrinos, electrons, and positrons produced in proton-proton collisions to those of the parameterisations of Kamae et al., Kelner et al., and Kafexhiu et al. We find that the differences in the normalisation of the photon energy spectra reach 20$-$50% at intermediate values of the transferred energy fraction $x$, growing u… ▽ More We compare the predictions of AAfrag for the spectra of secondary photons, neutrinos, electrons, and positrons produced in proton-proton collisions to those of the parameterisations of Kamae et al., Kelner et al., and Kafexhiu et al. We find that the differences in the normalisation of the photon energy spectra reach 20$-$50% at intermediate values of the transferred energy fraction $x$, growing up to a factor of two for $x\rightarrow\,$1, while the differences in the neutrino spectra are even larger. We argue that LHCf results on the forward production of photons favor the use of the QGSJET-II-04m model on which AAfrag is based. The differences in the normalisation have important implications in the context of multi-messenger astronomy, in particular, for the prediction of neutrino fluxes based on gamma-ray flux measurements, or regarding the inference of the cosmic ray spectrum, based on gamma-ray data. We note also that the positron-electron ratio from hadronic interactions increases with energy towards the cutoff, an effect which is missed using the average electron-positron spectrum from Kelner et al. Finally, we describe the publicly available python package aafragpy, which provides the secondary spectra of photons, neutrinos, electrons, and positrons. This package complements the AAfrag results for protons with energies above 4 GeV with previous analytical parameterizations of particle spectra for lower energy protons. △ Less

Submitted 17 December, 2021; v1 submitted 1 October, 2021; originally announced October 2021.

Comments: 8 pages, 10 figures

Journal ref: Phys. Rev. D 104, 123027 (2021)

New limit on high Galactic latitude PeV gamma-ray flux from Tibet AS-gamma data

Authors: A. Neronov, D. Semikoz, Ie. Vovk

Abstract: Tibet AS-gamma collaboration has recently reported detection of gamma-rays with energies up to Peta-electronvolt from parts of the Galactic plane. We notice that the analysis of gamma-ray flux by the Tibet AS-gamma experiment also implies an upper bound on the diffuse gamma-ray flux from high Galactic latitudes(|b|>20 degrees) in the energy range between 100 TeV and 1 PeV. This bound is up to an o… ▽ More Tibet AS-gamma collaboration has recently reported detection of gamma-rays with energies up to Peta-electronvolt from parts of the Galactic plane. We notice that the analysis of gamma-ray flux by the Tibet AS-gamma experiment also implies an upper bound on the diffuse gamma-ray flux from high Galactic latitudes(|b|>20 degrees) in the energy range between 100 TeV and 1 PeV. This bound is up to an order-of-magnitude stronger than previously derived bounds from GRAPES3, KASCADE and CASA-MIA experiments. We discuss the new TibetAS-gamma limit on high Galactic latitude gamma-ray flux in the context of possible mechanisms of multi-messenger (gamma-ray and neutrino) emission from nearby cosmic ray sources, dark matter decays and large scale cosmic ray halo of the Milky Way. △ Less

Submitted 6 September, 2021; v1 submitted 14 July, 2021; originally announced July 2021.

Comments: 4 pages, 3 figures, accepted to A&A

Journal ref: A&A 653, L4 (2021)

Account of the baryonic feedback effect in gamma-ray measurements of intergalactic magnetic fields

Authors: Kyrylo Bondarenko, Alexey Boyarsky, Alexander Korochkin, Andrii Neronov, Dmitri Semikoz, Anastasia Sokolenko

Abstract: Intergalactic magnetic fields in the voids of the large-scale structure can be probed via measurements of secondary gamma-ray emission from gamma-ray interactions with extragalactic background light. Lower bounds on the magnetic field in the voids were derived from the non-detection of this emission. It is not clear a priori what kind of magnetic field is responsible for the suppression of the sec… ▽ More Intergalactic magnetic fields in the voids of the large-scale structure can be probed via measurements of secondary gamma-ray emission from gamma-ray interactions with extragalactic background light. Lower bounds on the magnetic field in the voids were derived from the non-detection of this emission. It is not clear a priori what kind of magnetic field is responsible for the suppression of the secondary gamma-ray flux: a cosmological magnetic field that might be filling the voids or the field spread by galactic winds driven by star formation and active galactic nuclei. We used IllustrisTNG cosmological simulations to study the effect of magnetized galactic wind bubbles on the secondary gamma-ray flux. We show that within the IllustrisTNG model of baryonic feedback, the galactic wind bubbles typically provide energy-independent secondary flux suppression at the level of about 10%. The observed flux suppression effect has to be due to the cosmological magnetic field in the voids. This might not be the case for a special case when the primary gamma-ray source has a hard intrinsic gamma-ray spectrum that peaks in the energy range above 50 TeV. In this case, the observational data may be strongly affected by the magnetized bubble that is blown by the source host galaxy. △ Less

Submitted 3 February, 2022; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: V2: 8 pages, 7 figures, minor changes in the text. Version accepted for publication in A&A

Journal ref: A&A 660, A80 (2022)

Pion decay model of TIBET-AS$γ$ PeV gamma-ray signal

Authors: Sergey Koldobskiy, Andrii Neronov, Dmitri Semikoz

Abstract: Tibet-AS$γ$ collaboration has recently reported a measurement of diffuse $γ$-ray flux from the outer Galactic disk in the energy range reaching PeV. We complement this measurement with the Fermi/LAT measurement of the diffuse flux from the same sky region and study the pion decay model of the combined Fermi/LAT+Tibet-AS$γ$ spectrum. We find that within such a model the average cosmic ray spectrum… ▽ More Tibet-AS$γ$ collaboration has recently reported a measurement of diffuse $γ$-ray flux from the outer Galactic disk in the energy range reaching PeV. We complement this measurement with the Fermi/LAT measurement of the diffuse flux from the same sky region and study the pion decay model of the combined Fermi/LAT+Tibet-AS$γ$ spectrum. We find that within such a model the average cosmic ray spectrum in the outer Galactic disk has the same characteristic features as the local cosmic ray spectrum. In particular, it experiences a hardening at several hundred GV rigidity and a knee feature in the PV rigidity range. The slope of the average cosmic ray spectrum above the break is close to the locally observed slope of the helium spectrum $γ\simeq 2.5$, but is harder than the slope of the local proton spectrum in the same rigidity range. Although the combination of Fermi/LAT and Tibet-AS$γ$ data points to the presence of the knee in the average cosmic ray spectrum, the quality of the data is not yet sufficient for the study of knee shape and cosmic ray composition. △ Less

Submitted 17 August, 2021; v1 submitted 3 May, 2021; originally announced May 2021.

Comments: 5 pages, 2 figures

Journal ref: Phys. Rev. D 104, 043010 (2021)

Possibility of measurement of cosmic ray electron spectrum up to 100 TeV with two-layer water Cherenkov detector array

Authors: A. Neronov, D. Semikoz

Abstract: Measurements of cosmic ray electron+positron spectrum above 10 TeV with ground-based experiments is challenging because of the difficulty of rejection of hadronic extensive air shower background. We study the efficiency of rejection of the hadronic background with water Cherenkov detector array supplemented by muon detection layer. We show that addition of a continuous muon detection layer to the… ▽ More Measurements of cosmic ray electron+positron spectrum above 10 TeV with ground-based experiments is challenging because of the difficulty of rejection of hadronic extensive air shower background. We study the efficiency of rejection of the hadronic background with water Cherenkov detector array supplemented by muon detection layer. We show that addition of a continuous muon detection layer to the experimental setup allows to achieve a ~ 1e-5 rejection factor for hadronic background at 10 TeV and enables measurement of electron spectrum in 10-100 TeV energy range. We show that measurements of electron spectrum in this energy range do not require a high-altitude experiment and can be done with a sea-level detector. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: 5 pages, 5 figures

The POEMMA (Probe of Extreme Multi-Messenger Astrophysics) Observatory

Authors: A. V. Olinto, J. Krizmanic, J. H. Adams, R. Aloisio, L. A. Anchordoqui, A. Anzalone, M. Bagheri, D. Barghini, M. Battisti, D. R. Bergman, M. E. Bertaina, P. F. Bertone, F. Bisconti, M. Bustamante, F. Cafagna, R. Caruso, M. Casolino, K. Černý, M. J. Christl, A. L. Cummings, I. De Mitri, R. Diesing, R. Engel, J. Eser, K. Fang , et al. (51 additional authors not shown)

Abstract: The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to accurately observe ultra-high-energy cosmic rays (UHECRs) and cosmic neutrinos from space with sensitivity over the full celestial sky. POEMMA will observe the extensive air showers (EASs) from UHECRs and UHE neutrinos above 20 EeV via air fluorescence. Additionally, POEMMA will observe the Cherenkov signal from upward-movin… ▽ More The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to accurately observe ultra-high-energy cosmic rays (UHECRs) and cosmic neutrinos from space with sensitivity over the full celestial sky. POEMMA will observe the extensive air showers (EASs) from UHECRs and UHE neutrinos above 20 EeV via air fluorescence. Additionally, POEMMA will observe the Cherenkov signal from upward-moving EASs induced by Earth-interacting tau neutrinos above 20 PeV. The POEMMA spacecraft are designed to quickly re-orientate to follow up transient neutrino sources and obtain unparalleled neutrino flux sensitivity. Developed as a NASA Astrophysics Probe-class mission, POEMMA consists of two identical satellites flying in loose formation in 525 km altitude orbits. Each POEMMA instrument incorporates a wide field-of-view (45$^\circ$) Schmidt telescope with over 6 m$^2$ of collecting area. The hybrid focal surface of each telescope includes a fast (1~$μ$s) near-ultraviolet camera for EAS fluorescence observations and an ultrafast (10~ns) optical camera for Cherenkov EAS observations. In a 5-year mission, POEMMA will provide measurements that open new multi-messenger windows onto the most energetic events in the universe, enabling the study of new astrophysics and particle physics at these otherwise inaccessible energies. △ Less

Submitted 24 May, 2021; v1 submitted 14 December, 2020; originally announced December 2020.

Comments: 66 pages, 53 figures; JCAP accepted version

Journal ref: JCAP 06 (2021) 007

Radio-to-gamma-ray synchrotron and neutrino emission from proton-proton interactions in active galactic nuclei

Authors: Andrii Neronov, Dmitry Semikoz

Abstract: We explore possible physical origin of correlation between radio wave and very-high-energy neutrino emission in active galactic nuclei (AGN), suggested by recently reported evidence for correlation between neutrino arrival directions and positions of brightest radio-loud AGN. We show that such correlation is expected if both synchrotron emitting electrons and neutrinos originate from decays of cha… ▽ More We explore possible physical origin of correlation between radio wave and very-high-energy neutrino emission in active galactic nuclei (AGN), suggested by recently reported evidence for correlation between neutrino arrival directions and positions of brightest radio-loud AGN. We show that such correlation is expected if both synchrotron emitting electrons and neutrinos originate from decays of charged pions produced in proton-proton interactions in parsec-scale relativistic jet propagating through circum-nuclear medium of the AGN. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 6 pages, 1 figure, accepted to JETP Lett

Multiwavelength variability and correlation studies of Mrk 421 during historically low X-ray and $γ$-ray activity in 2015$-$2016

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

Abstract: We report a characterization of the multi-band flux variability and correlations of the nearby (z=0.031) blazar Markarian 421 (Mrk 421) using data from Metsähovi, Swift, Fermi-LAT, MAGIC, FACT and other collaborations and instruments from November 2014 till June 2016. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F… ▽ More We report a characterization of the multi-band flux variability and correlations of the nearby (z=0.031) blazar Markarian 421 (Mrk 421) using data from Metsähovi, Swift, Fermi-LAT, MAGIC, FACT and other collaborations and instruments from November 2014 till June 2016. Mrk 421 did not show any prominent flaring activity, but exhibited periods of historically low activity above 1 TeV (F$_{>1\mathrm{TeV}}<$ 1.7$\times$10$^{-12}$ ph cm$^{-2}$ s$^{-1}$) and in the 2-10 keV (X-ray) band (F$_{2-10 \mathrm{keV}}<$3.6$\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$), during which the Swift-BAT data suggests an additional spectral component beyond the regular synchrotron emission. The highest flux variability occurs in X-rays and very-high-energy (E$>$0.1 TeV) $γ$-rays, which, despite the low activity, show a significant positive correlation with no time lag. The HR$_\mathrm{keV}$ and HR$_\mathrm{TeV}$ show the harder-when-brighter trend observed in many blazars, but the trend flattens at the highest fluxes, which suggests a change in the processes dominating the blazar variability. Enlarging our data set with data from years 2007 to 2014, we measured a positive correlation between the optical and the GeV emission over a range of about 60 days centered at time lag zero, and a positive correlation between the optical/GeV and the radio emission over a range of about 60 days centered at a time lag of $43^{+9}_{-6}$ days.This observation is consistent with the radio-bright zone being located about 0.2 parsec downstream from the optical/GeV emission regions of the jet. The flux distributions are better described with a LogNormal function in most of the energy bands probed, indicating that the variability in Mrk 421 is likely produced by a multiplicative process. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: 45 pages, 30 figures, 6 tables, accepted for publication in MNRAS

The INTEGRAL view of the pulsating hard X-ray sky: from accreting and transitional millisecond pulsars to rotation-powered pulsars and magnetars

Authors: A. Papitto, M. Falanga, W. Hermsen, S. Mereghetti, L. Kuiper, J. Poutanen, E. Bozzo, F. Ambrosino, F. Coti Zelati, V. De Falco, D. de Martino, T. Di Salvo, P. Esposito, C. Ferrigno, M. Forot, D. Götz, C. Gouiffes, R. Iaria, P. Laurent, J. Li, Z. Li, T. Mineo, P. Moran, A. Neronov, A. Paizis , et al. (8 additional authors not shown)

Abstract: In the last 25 years, a new generation of X-ray satellites imparted a significant leap forward in our knowledge of X-ray pulsars. The discovery of accreting and transitional millisecond pulsars proved that disk accretion can spin up a neutron star to a very high rotation speed. The detection of MeV-GeV pulsed emission from a few hundreds of rotation-powered pulsars probed particle acceleration in… ▽ More In the last 25 years, a new generation of X-ray satellites imparted a significant leap forward in our knowledge of X-ray pulsars. The discovery of accreting and transitional millisecond pulsars proved that disk accretion can spin up a neutron star to a very high rotation speed. The detection of MeV-GeV pulsed emission from a few hundreds of rotation-powered pulsars probed particle acceleration in the outer magnetosphere, or even beyond. Also, a population of two dozens of magnetars has emerged. INTEGRAL played a central role to achieve these results by providing instruments with high temporal resolution up to the hard X-ray/soft gamma-ray band and a large field of view imager with good angular resolution to spot hard X-ray transients. In this article, we review the main contributions by INTEGRAL to our understanding of the pulsating hard X-ray sky, such as the discovery and characterization of several accreting and transitional millisecond pulsars, the generation of the first catalog of hard X-ray/soft gamma-ray rotation-powered pulsars, the detection of polarization in the hard X-ray emission from the Crab pulsar, and the discovery of persistent hard X-ray emission from several magnetars. △ Less

Submitted 2 December, 2020; originally announced December 2020.

Comments: Accepted for publication on New Astronomy Reviews as invited contribution

Viscous damping of chiral dynamos in the early universe

Authors: A. Neronov, D. Semikoz

Abstract: Chiral dynamo converting asymmetry between right and left-handed leptons in the early universe into helical magnetic field has been proposed as a possible cosmological magnetogenesis scenario. We show that this mechanism is strongly affected by viscous damping of primordial plasma motions excited by the dynamo. This effect modifies the expected range of strength and correlation length of the chira… ▽ More Chiral dynamo converting asymmetry between right and left-handed leptons in the early universe into helical magnetic field has been proposed as a possible cosmological magnetogenesis scenario. We show that this mechanism is strongly affected by viscous damping of primordial plasma motions excited by the dynamo. This effect modifies the expected range of strength and correlation length of the chiral dynamo field which could have survived till present epoch in the voids of the Large Scale Structure. We show the range of parameters of chiral dynamo field that may have survived in the voids is still consistent with existing lower bounds on intergalactic magnetic field from gamma-ray observations, but only if the right-left lepton asymmetry at the temperature T~80 TeV is very high, close to the maximal possible value. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Comments: 6 pages, 1 figure

NANOGrav signal from MHD turbulence at QCD phase transition in the early universe

Authors: A. Neronov, A. Roper Pol, C. Caprini, D. Semikoz

Abstract: The NANOGrav collaboration has recently reported evidence for the existence of a stochastic gravitational wave background in the 1-100 nHz frequency range. We argue that such background could have been produced by magneto-hydrodynamic (MHD) turbulence at the QCD scale. From the NANOGrav measurement one can infer the magnetic field parameters: comoving field strength close to microGauss and a corre… ▽ More The NANOGrav collaboration has recently reported evidence for the existence of a stochastic gravitational wave background in the 1-100 nHz frequency range. We argue that such background could have been produced by magneto-hydrodynamic (MHD) turbulence at the QCD scale. From the NANOGrav measurement one can infer the magnetic field parameters: comoving field strength close to microGauss and a correlation length close to 10\% of the Hubble radius at the QCD phase transition epoch. We point out that the turbulent decay of a non-helical magnetic field with such parameters leads to a magnetic field at the recombination epoch, which would be sufficiently strong to provide a solution to the Hubble tension problem, as recently proposed. We also show that the MHD turbulence interpretation of the NANOGrav signal can be tested via measurements of the relic magnetic field in the voids of the large scale structure, with gamma-ray telescopes like CTA. △ Less

Submitted 16 March, 2021; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: 6 pages, 2 figures

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

Advances in Understanding High-Mass X-ray Binaries with INTEGRAL and Future Directions

Authors: Peter Kretschmar, Felix Fürst, Lara Sidoli, Enrico Bozzo, Julia Alfonso-Garzón, Arash Bodaghee, Sylvain Chaty, Masha Chernyakova, Carlo Ferrigno, Antonios Manousakis, Ignacio Negueruela, Konstantin Postnov, Adamantia Paizis, Pablo Reig, José Joaquín Rodes-Roca, Sergey Tsygankov, Antony J. Bird, Matthias Bissinger né Kühnel, Pere Blay, Isabel Caballero, Malcolm J. Coe, Albert Domingo, Victor Doroshenko, Lorenzo Ducci, Maurizio Falanga , et al. (26 additional authors not shown)

Abstract: High mass X-ray binaries are among the brightest X-ray sources in the Milky Way, as well as in nearby Galaxies. Thanks to their highly variable emissions and complex phenomenology, they have attracted the interest of the high energy astrophysical community since the dawn of X-ray Astronomy. In more recent years, they have challenged our comprehension of physical processes in many more energy bands… ▽ More High mass X-ray binaries are among the brightest X-ray sources in the Milky Way, as well as in nearby Galaxies. Thanks to their highly variable emissions and complex phenomenology, they have attracted the interest of the high energy astrophysical community since the dawn of X-ray Astronomy. In more recent years, they have challenged our comprehension of physical processes in many more energy bands, ranging from the infrared to very high energies. In this review, we provide a broad but concise summary of the physical processes dominating the emission from high mass X-ray binaries across virtually the whole electromagnetic spectrum. These comprise the interaction of stellar winds with the high gravitational and magnetic fields of compact objects, the behaviour of matter under extreme magnetic and gravity conditions, and the perturbation of the massive star evolutionary processes by presence in a binary system. We highlight the role of the INTEGRAL mission in the discovery of many of the most interesting objects in the high mass X-ray binary class and its contribution in reviving the interest for these sources over the past two decades. We show how the INTEGRAL discoveries have not only contributed to significantly increase the number of high mass X-ray binaries known, thus advancing our understanding of the population as a whole, but also have opened new windows of investigation that stimulated the multi-wavelength approach nowadays common in most astrophysical research fields. We conclude the review by providing an overview of future facilities being planned from the X-ray to the very high energy domain that will hopefully help us in finding an answer to the many questions left open after more than 18 years of INTEGRAL scientific observations. △ Less

Submitted 7 September, 2020; originally announced September 2020.

Comments: 61 pages, 22 figures, 1 table, will be published in New Astronomy Reviews

Snowmass 2021 Letter of Interest: The Probe Of Multi-Messenger Astrophysics (POEMMA)

Authors: A. V. Olinto, F. Sarazin, J. H. Adams, R. Aloisio, L. A. Anchordoqui, M. Bagheri, D. Barghini, M. Battisti, D. R. Bergman, M. E. Bertaina, P. F. Bertone, F. Bisconti, M. Bustamante, M. Casolino, M. J. Christl, A. L. Cummings, I. De Mitri, R. Diesing, R. Engel, J. Eser, K. Fang, G. Fillipatos, F. Fenu, E. Gazda, C. Guepin , et al. (39 additional authors not shown)

Abstract: The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to identify the sources of Ultra-High-Energy Cosmic Rays (UHECRs) and to observe cosmic neutrinos, both with full-sky coverage. Developed as a NASA Astrophysics Probe-class mission, POEMMA consists of two spacecraft flying in a loose formation at 525 km altitude, 28.5 deg inclination orbits. Each spacecraft hosts a Schmidt tele… ▽ More The Probe Of Extreme Multi-Messenger Astrophysics (POEMMA) is designed to identify the sources of Ultra-High-Energy Cosmic Rays (UHECRs) and to observe cosmic neutrinos, both with full-sky coverage. Developed as a NASA Astrophysics Probe-class mission, POEMMA consists of two spacecraft flying in a loose formation at 525 km altitude, 28.5 deg inclination orbits. Each spacecraft hosts a Schmidt telescope with a large collecting area and wide field of view. A novel focal plane is optimized to observe both the UV fluorescence signal from extensive air showers (EASs) and the beamed optical Cherenkov signals from EASs. In POEMMA-stereo fluorescence mode, POEMMA will measure the spectrum, composition, and full-sky distribution of the UHECRs above 20 EeV with high statistics along with remarkable sensitivity to UHE neutrinos. The spacecraft are designed to quickly re-orient to a POEMMA-limb mode to observe neutrino emission from Target-of-Opportunity (ToO) transient astrophysical sources viewed just below the Earth's limb. In this mode, POEMMA will have unique sensitivity to cosmic neutrino tau events above 20 PeV by measuring the upward-moving EASs induced by the decay of the emerging tau leptons following the interactions of neutrino tau inside the Earth. △ Less

Submitted 1 September, 2020; v1 submitted 29 August, 2020; originally announced August 2020.

Comments: 3 pages, no figures, LoI for Snowmass 2021

Sensitivity reach of gamma-ray measurements for strong cosmological magnetic fields

Authors: A. Korochkin, O. Kalashev, A. Neronov, D. Semikoz

Abstract: A primordial magnetic field with the strength in the 1-10 pG range can resolve the tension between different measurements of the Hubble constant and provide an explanation for the excess opacity in the 21 cm line at redshift $15<z<20$, if it is present during the recombination and reionization epochs. This field can also survive in the voids of the large-scale Structure in the present day universe… ▽ More A primordial magnetic field with the strength in the 1-10 pG range can resolve the tension between different measurements of the Hubble constant and provide an explanation for the excess opacity in the 21 cm line at redshift $15<z<20$, if it is present during the recombination and reionization epochs. This field can also survive in the voids of the large-scale Structure in the present day universe. We study the sensitivity reach of the gamma-ray technique for measurement of such relatively strong cosmological magnetic field using deep exposure(s) of the nearest hard spectrum blazar(s) with CTA telescopes. We show that the gamma-ray measurement method can sense the primordial magnetic field with a strength of up to $10^{-11}$~G. Combination of the cosmic microwave background and gamma-ray constraints can thus sense the full range of possible cosmological magnetic fields to confirm or rule out their relevance to the problem of the origin of cosmic magnetic fields, as well as their influence on recombination and reionization epochs. △ Less

Submitted 6 December, 2020; v1 submitted 28 July, 2020; originally announced July 2020.

Comments: 9 pages, 7 figures, accepted to ApJ

Journal ref: 2021, ApJ, 906, 116