PeV Energies Research Articles

Search for Joint Multimessenger Signals from Potential Galactic Cosmic-Ray Accelerators with HAWC and IceCube

The origin of high-energy galactic cosmic rays is yet to be understood, but some galactic cosmic-ray accelerators can accelerate cosmic rays up to PeV energies. The high-energy cosmic rays are expected to interact with the surrounding material or radiation, resulting in the production of gamma-rays and neutrinos. To optimize for the detection of such associated production of gamma-rays and neutrinos for a given source morphology and spectrum, a multimessenger analysis that combines gamma-rays and neutrinos is required. In this study, we use the Multi-Mission Maximum Likelihood framework with IceCube Maximum Likelihood Analysis software and HAWC Accelerated Likelihood to search for a correlation between 22 known gamma-ray sources from the third HAWC gamma-ray catalog and 14 yr of IceCube track-like data. No significant neutrino emission from the direction of the HAWC sources was found. We report the best-fit gamma-ray model and 90% CL neutrino flux limit from the 22 sources. From the neutrino flux limit, we conclude that, for five of the sources, the gamma-ray emission observed by HAWC cannot be produced purely from hadronic interactions. We report the limit for the fraction of gamma-rays produced by hadronic interactions for these five sources.

Prospects for a survey of the galactic plane with the Cherenkov Telescope Array

Approximately one hundred sources of very-high-energy (VHE) gamma rays are known in the Milky Way, detected with a combination of targeted observations and surveys. A survey of the entire Galactic Plane in the energy range from a few tens of GeV to a few hundred TeV has been proposed as a Key Science Project for the upcoming Cherenkov Telescope Array Observatory (CTAO). This article presents the status of the studies towards the Galactic Plane Survey (GPS). We build and make publicly available a sky model that combines data from recent observations of known gamma-ray emitters with state-of-the-art physically-driven models of synthetic populations of the three main classes of established Galactic VHE sources (pulsar wind nebulae, young and interacting supernova remnants, and compact binary systems), as well as of interstellar emission from cosmic-ray interactions in the Milky Way. We also perform an optimisation of the observation strategy (pointing pattern and scheduling) based on recent estimations of the instrument performance. We use the improved sky model and observation strategy to simulate GPS data corresponding to a total observation time of 1620 hours spread over ten years. Data are then analysed using the methods and software tools under development for real data. Under our model assumptions and for the realisation considered, we show that the GPS has the potential to increase the number of known Galactic VHE emitters by almost a factor of five. This corresponds to the detection of more than two hundred pulsar wind nebulae and a few tens of supernova remnants at average integral fluxes one order of magnitude lower than in the existing sample above 1 TeV, therefore opening the possibility to perform unprecedented population studies. The GPS also has the potential to provide new VHE detections of binary systems and pulsars, to confirm the existence of a hypothetical population of gamma-ray pulsars with an additional TeV emission component, and to detect bright sources capable of accelerating particles to PeV energies (PeVatrons). Furthermore, the GPS will constitute a pathfinder for deeper follow-up observations of these source classes. Finally, we show that we can extract from GPS data an estimate of the contribution to diffuse emission from unresolved sources, and that there are good prospects of detecting interstellar emission and statistically distinguishing different scenarios.Thus, a survey of the entire Galactic plane carried out from both hemispheres with CTAO will ensure a transformational advance in our knowledge of Galactic VHE source populations and interstellar emission.

The remarkable microquasar S26: A super-Eddington PeVatron

S26 is an extragalactic microquasar with the most powerful jets ever discovered. They have a kinetic luminosity of $L_ j $. This implies that the accretion power to the black hole should be super-Eddington, of the order of $L_ acc j $. However, the observed X-ray flux of this system indicates an apparent very sub-Eddington accretion luminosity of $L_ X We aim to characterize the nature of S26, explain the system emission, and study the feasibility of super-Eddington microquasars as potential PeVatron sources. We first analyze multi-epoch X-ray observations of S26 obtained with XMM-Newton and model the super-Eddington disk and its wind. We then develop a jet model and study the particle acceleration and radiative processes that occur in shocks generated near the base of the jet and in its terminal region. We find that the discrepancy between the jet and the apparent disk luminosities in S26 is caused by the complete absorption of the disk radiation by the wind ejected from the super-Eddington disk. The non-thermal X-rays are produced near the base of the jet, and the thermal X-rays are emitted in the terminal regions. The radio emission observed with the Australia Telescope Compact Array can be explained as synchrotron radiation produced at the reverse shock in the lobes. We also find that S26 can accelerate protons to PeV energies in both the inner jet and the lobes. The ultra-high energy protons accelerated in the lobes of S26 are injected into the interstellar medium with a total power of $ We conclude that S26 is a super-Eddington microquasar with a dense disk-driven wind that obscures the X-ray emission from the inner disk, and that the supercritical nature of the system allows the acceleration of cosmic rays to PeV energies.

Shower formation constraints on cubic Lorentz invariance violation parameters in quantum electrodynamics

We present shower formation constraints on the Lorentz invariance violation (LIV) energy scale for photons with a cubic dispersion relation from recent gamma-ray observations in the 100 TeV–PeV energy range by the Large High Altitude Air Shower Observatory (LHAASO). We assume a Myers–Pospelov effective field theory framework, and calculate the suppression for the Bethe–Heitler process which is mainly responsible for the formation of photon-initiated atmospheric showers. Comparing the high-energy photon events with the suppressed flux predictions, we obtain 95%\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$95\\%$$\\end{document} CL constraints on the LIV energy scale. The obtained shower constraint ELIV≳O1020GeV\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$E_\ extrm{LIV} \\gtrsim \\mathcal {O} \\left( 10^{20} \ ext{ GeV }\\right) $$\\end{document} is significantly weaker than existing birefringence constraints but is independent.

Hadronuclear Interactions in Active Galactic Nuclei Jets as the Origin of the Diffuse High-energy Neutrino Background

The origin of diffuse high-energy neutrinos from TeV to PeV energies detected by the IceCube Observatory remains a mystery. In our previous work, we have shown that hadronuclear (p − p) interactions in active galactic nuclei (AGNs) jets could be important and generate detectable very-high-energy emissions. Here, we further explore these interactions in the AGN jets based on their luminosity function. The diffuse neutrino flux and corresponding γ-ray flux have been calculated and compared with observational data. In our modeling, two beaming patterns are considered separately. To make sure that the corresponding γ-ray flux does not overshoot the diffuse γ-ray background, we find that if the neutrino production region in a jet is opaque to γ-rays, p − p interactions in AGN jets with a small viewing angle (the blazar case) are able to interpret the PeV neutrino background. Similarly, AGN jets with a large viewing angle (the radio galaxy case) may interpret the TeV neutrino background. While, if the neutrino production region is transparent to γ-rays, only blazars have the potential to interpret the diffuse neutrino background around the PeV band. Some caveats are also discussed.

Evidence for particle acceleration approaching PeV energies in the W51 complex

The γ-ray emission from the W51 complex is widely acknowledged to be attributed to the interaction between the cosmic rays (CRs) accelerated by the shock of supernova remnant (SNR) W51C and the dense molecular clouds in the adjacent star-forming region, W51B. However, the maximum acceleration capability of W51C for CRs remains elusive. Based on observations conducted with the Large High Altitude Air Shower Observatory (LHAASO), we report a significant detection of γ rays emanating from the W51 complex, with energies from 2 to 200 TeV. The LHAASO measurements, for the first time, extend the γ-ray emission from the W51 complex beyond 100 TeV and reveal a significant spectrum bending at tens of TeV. By combining the “π0-decay bump” featured data from Fermi-LAT, the broadband γ-ray spectrum of the W51 region can be well-characterized by a simple pp-collision model. The observed spectral bending feature suggests an exponential cutoff at ∼400 TeV or a power-law break at ∼200 TeV in the CR proton spectrum, most likely providing the first evidence of SNRs serving as CR accelerators approaching the PeV regime. Additionally, two young star clusters within W51B could also be theoretically viable to produce the most energetic γ rays observed by LHAASO. Our findings strongly support the presence of extreme CR accelerators within the W51 complex and provide new insights into the origin of Galactic CRs.

Very-high-energy γ-Ray Emission from Young Massive Star Clusters in the Large Magellanic Cloud

The Tarantula Nebula in the Large Magellanic Cloud is known for its high star formation activity. At its center lies the young massive star cluster R136, providing a significant amount of the energy that makes the nebula shine so brightly at many wavelengths. Recently, young massive star clusters have been suggested to also efficiently produce very high-energy cosmic rays, potentially beyond PeV energies. Here, we report the detection of very-high-energy γ-ray emission from the direction of R136 with the High Energy Stereoscopic System, achieved through a multicomponent, likelihood-based modeling of the data. This supports the hypothesis that R136 is indeed a very powerful cosmic-ray accelerator. Moreover, from the same analysis, we provide an updated measurement of the γ-ray emission from 30 Dor C, the only superbubble detected at TeV energies presently. The γ-ray luminosity above 0.5 TeV of both sources is (2–3) × 1035 erg s−1. This exceeds by more than a factor of 2 the luminosity of HESS J1646−458, which is associated with the most massive young star cluster in the Milky Way, Westerlund 1. Furthermore, the γ-ray emission from each source is extended with a significance of >3σ and a Gaussian width of about 30 pc. For 30 Dor C, a connection between the γ-ray emission and the nonthermal X-ray emission appears likely. Different interpretations of the γ-ray signal from R136 are discussed.

A hybrid approach to event reconstruction for atmospheric Cherenkov Telescopes combining machine learning and likelihood fitting

The imaging atmospheric Cherenkov technique provides potentially the highest angular resolution achievable in astronomy at energies above the X-ray waveband. High-resolution measurements provide the key to progress on many of the major questions in high energy astrophysics, including the sites and mechanisms of particle acceleration to PeV energies. The huge potential of the next-generation CTA observatory in this regard can be realised with the help of improved algorithms for the reconstruction of the air-shower direction and energy.Hybrid methods combining maximum-likelihood-fitting techniques with neural networks represent a particularly promising approach and have recently been successfully applied for the reconstruction of astrophysical neutrinos. Here, we present the FreePACT algorithm, a hybrid reconstruction method for IACTs. In this, making use of the neural ratio estimation technique from the field of likelihood-free inference, the analytical likelihood used in traditional image likelihood fitting is replaced by a neural network that approximates the charge probability density function for each pixel in the camera.The performance of this improved algorithm is demonstrated using simulations of the planned CTA southern array. For this setupFreePACT provides significant performance improvements over analytical likelihood techniques, with improvements in angular and energy resolution of 25% or more over a wide energy range and an angular resolution as low as 40′′ at energies above 50TeV for observations at 20° zenith angle. It also yields more accurate estimations of the uncertainties on the reconstructed parameters and significantly speeds up the reconstruction compared to analytical likelihood techniques while showing the same stability with respect to changes in the observation conditions. Therefore, the FreePACT method is a promising upgrade over the current state-of-the-art likelihood event reconstruction techniques.

Studying the role of supernova remnants as cosmic ray PeVatron with LHAASO observations

It is commonly believed that galactic cosmic rays are produced in supernova remnants (SNRs) and accelerated via a diffusive shock acceleration (DSA) mechanism in supernova blast waves driven by expanding SNRs. The latest theoretical advancement of the diffusive shock acceleration hypothesis in SNRs shows that cosmic rays may be accelerated up to the knee energy of the observed cosmic ray spectrum under the amplified magnetic field scenario. There is, however, no empirical evidence to support SNRs as sources of hadrons with energies larger than a few tens of TeV. Very recently, LHAASO observatory reported the very high-energy gamma ray emission between TeV to PeV energy range from two SNRs, Cassiopeia A and IC 443. Above 25 TeV energies, non-detection of gamma-ray flux by LHAASO yields a strong upper limit. In this work, we investigate the implications of the acceleration of cosmic ray protons in the SNR on energetic gamma rays produced in the hadronic interaction of cosmic rays with ambient matter. Our findings imply that when we consider the highest attainable energy of cosmic ray protons in the SNRs to be about 100 TeV, the observed gamma-ray spectra from the two SNRs can be described consistently. Therefore, we conclude that the Cassiopia A and IC 443 SNRs are unlikely to be cosmic ray PeVatrons.

Search for Synchrotron Emission from Secondary Electrons of Proton–Proton Interactions in Galactic PeVatron Candidate HESS J1641–463

HESS J1641−463 is an unidentified gamma-ray source with a hard TeV gamma-ray spectrum, and thus it has been proposed to be a possible candidate for a cosmic-ray (CR) accelerator up to PeV energies (a PeVatron candidate). The source spatially coincides with the radio supernova remnant G338.5+0.1 but has not yet been fully explored in the X-ray band. We analyzed newly taken NuSTAR data, pointing at HESS J1641−463, with 82 ks effective exposure time. There is no apparent X-ray counterpart of HESS J1641−463, while nearby stellar cluster, Mercer 81, and stray-light X-rays are detected. Combined with the archival Chandra data, partially covering the source, we derived an upper limit of ∼6 × 10−13 erg cm−2 s−1 in 2–10 keV (∼3 × 10−13 erg cm−2 s−1 in 10–20 keV). If the gamma-ray emission is originated from the decay of π 0 mesons produced in interactions between CR protons and ambient materials, secondary electrons in the proton–proton interactions can potentially emit synchrotron photons in the X-ray band, which can be tested by our X-ray observations. Although the obtained X-ray upper limits cannot place a constraint on the primary proton spectrum, it will be possible with a future hard X-ray mission.

LHAASO J2108+5157 as a molecular cloud illuminated by a supernova remnant

Context. The search for Galactic PeVatrons – astrophysical accelerators of cosmic rays to PeV energies – has entered a new phase in recent years with the discovery of the first ultra-high-energy (UHE, E > 100 TeV) γ-ray sources by the High Altitude Water Cherenkov (HAWC) observatory and Large High Altitude Air Shower Observatory (LHAASO). Establishing whether the emission is leptonic or hadronic in nature, however, requires multi-wavelength data and modelling studies. Among the currently known UHE sources, LHAASO J2108+5157 is an enigmatic source without clear association to a plausible accelerator, yet spatially coincident with molecular clouds. Aims. We aim to investigate the scenario of a molecular cloud illuminated by cosmic rays accelerated in a nearby supernova remnant (SNR) as an explanation for LHAASO J2108+5157. We aim to constrain the required properties of the SNR as well as which of the clouds identified in the vicinity is the most likely association. Methods. We used a model for cosmic-ray acceleration in SNRs, their transport through the interstellar medium, and subsequent interaction with molecular material to predict the corresponding γ-ray emission. The parameter space of SNR properties was explored to find the most plausible parameter combination that can account for the γ-ray spectrum of LHAASO J2108+5157. Results. In the case that a SNR is illuminating the cloud, we find that it must be young (< 10 kyr) and located within 40–60 pc of the cloud. A SN scenario with a low Sedov time is preferred, with a maximum proton energy of 3 PeV assumed. No SNRs matching these properties are currently known, although an as yet undetected SNR remains feasible. The Galactic CR sea is insufficient to solely account for the observed flux, such that a PeVatron accelerator must be present in the vicinity.

Modeling X-Ray and Gamma-Ray Emission from Redback Pulsar Binaries

We investigated the multiband emission from the pulsar binaries XSS J12270−4859, PSR J2039−5617, and PSR J2339−0533, which exhibit orbital modulation in the X-ray and gamma-ray bands. We constructed the sources’ broadband spectral energy distributions and multiband orbital light curves by supplementing our X-ray measurements with published gamma-ray results, and we modeled the data using intrabinary shock (IBS) scenarios. While the X-ray data were well explained by synchrotron emission from electrons/positrons in the IBS, the gamma-ray data were difficult to explain with the IBS components alone. Therefore, we explored other scenarios that had been suggested for gamma-ray emission from pulsar binaries: (1) inverse-Compton emission in the upstream unshocked wind zone and (2) synchrotron radiation from electrons/positrons interacting with the kilogauss magnetic field of the companion. Scenario (1) requires that the bulk motion of the wind substantially decelerates to ∼1000 km s−1 before reaching the IBS for increased residence time, in which case the formation of a strong shock is untenable, inconsistent with the X-ray phenomenology. Scenario (2) can explain the data if we assume the presence of electrons/positrons with a Lorentz factor ∼ 108 (∼0.1 PeV) that pass through the IBS and tap a substantial portion of the pulsar voltage drop. These findings raise the possibility that the orbitally modulating gamma-ray signals from pulsar binaries can provide insights into the flow structure and energy conversion within pulsar winds and particle acceleration nearing PeV energies in pulsars. These signals may also yield greater understanding of kilogauss magnetic fields potentially hosted by the low-mass stars in these systems.

Does or Did the Supernova Remnant Cassiopeia A Operate as a PeVatron?

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 γ ≥ 100 TeV) γ-rays. In this context, the historical SNR Cassiopeia A (Cas A) is considered one of the most promising targets 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.

The CaloCube calorimeter for high-energy cosmic-ray measurements in space: Response of a large-scale prototype to protons

Direct observation of cosmic rays nuclei is currently limited to energies of the order of hundreds of TeV. In order to extend these observations to higher energies, detectors capable of operating in space with high geometric factor and energy resolution are needed. In particular, highly performing calorimeters based on the CaloCube design can allow to carry out cosmic ray measurements in the PeV energy region. The CaloCube R&D project foresees the installation in space of a homogeneous and isotropic calorimeter composed of cubic scintillator crystals arranged to form a cube of about tons weight, with a high acceptance and capable of detecting particles coming from any direction. A prototype, composed of 5 × 5 × 18 CsI(Tl) crystals, has been tested on high-energy particle beams at CERN SPS accelerator and the results relative to the calorimeter response to protons are reported in this document.

The Maximum Energy of Shock-accelerated Cosmic Rays

Identifying the accelerators of Galactic cosmic ray (CR) protons with energies up to a few PeV (1015 eV) remains a theoretical and observational challenge. Supernova remnants (SNRs) represent strong candidates because they provide sufficient energetics to reproduce the CR flux observed at Earth. However, it remains unclear whether they can accelerate particles to PeV energies, particularly after the very early stages of their evolution. This uncertainty has prompted searches for other source classes and necessitates comprehensive theoretical modeling of the maximum proton energy, , accelerated by an arbitrary shock. While analytic estimates of have been put forward in the literature, they do not fully account for the complex interplay between particle acceleration, magnetic field amplification, and shock evolution. This paper uses a multizone, semianalytic model of particle acceleration based on kinetic simulations to place constraints on for a wide range of astrophysical shocks. In particular, we develop relationships between , shock velocity, size, and ambient medium. We find that SNRs can only accelerate PeV particles under a select set of circumstances, namely, if the shock velocity exceeds ∼104 km s−1 and escaping particles drive magnetic field amplification. However, older and slower SNRs may still produce observational signatures of PeV particles due to populations accelerated when the shock was younger. Our results serve as a reference for modelers seeking to quickly produce a self-consistent estimate of the maximum energy accelerated by an arbitrary astrophysical shock. 1 1 Presented as a thesis to the Department of Astronomy and Astrophysics, The University of Chicago, in partial fulfillment of the requirements for a Ph.D. degree.

EUSO-SPB2: A sub-orbital cosmic ray and neutrino multi-messenger pathfinder observatory

The next generation of ultra-high energy cosmic ray (UHECR) and very-high energy neutrino observatories will address the challenge of the extremely low fluxes of these particles at the highest energies. EUSO-SPB2 (Extreme Universe Space Observatory on a Super Pressure Balloon 2) is designed to prepare space missions to address this challenge. EUSO-SPB2 is equipped with 2 telescopes: the Fluorescence Telescope, which will point downwards and measure fluorescence emission from UHECR air showers with an energy above 2 EeV, and the Cherenkov Telescope (CT), which will point towards the Earth’s limb and measure direct Cherenkov emission from cosmic rays with energies above 1 PeV verifying the technique. Pointed below the limb, the CT will search for Cherenkov emission produced by neutrino-sourced tau-lepton decays above 10 PeV energies and study backgrounds for such events. The EUSO-SPB2 mission will provide pioneering observations and technical milestones on the path towards a space-based multi-messenger observatory.

Cosmic ray mass composition at the knee using azimuthal fluctuations of air shower particles detected at ground by the KASCADE experiment

The presence of hadronic sub-showers causes azimuthal non-uniformity in the particle distributions on the ground in vertical air showers. The LCm parameter, which quantifies the non-uniformity of the signal recorded in detectors located at a given distance on a ring around the shower axis, has been successfully introduced as a gamma/hadron discriminator at PeV energies [23]. In this work, we demonstrate that the LCm parameter can effectively serve as a mass composition discriminator in experiments that employ a compact array of detectors, like KASCADE. We reconstruct the LCm parameter distributions in the energy range lg(E/eV) = [15.0 – 16.0] using measurements from the KASCADE experiment, with intervals of lg(E/eV) = 0.2, which are then fitted with MC templates for five primary nuclei species p, He, C, Si, and Fe considering three hadronic interaction models: QGSjet-II-04, EPOS-LHC and SIBYLL 2.3d. We find that the LCm parameter exhibits minimal dependence on the specific hadronic interaction model considered. The reconstructed fractions of individual species demonstrate a linear decrease in the abundance of protons and He nuclei with increasing energy, while the heavier components become prevalent above the knee as predicted by all three hadronic interaction models. Our findings indicate that the abundance of particle types as a function of energy aligns with different astrophysical models that link the knee to the acceleration and propagation of cosmic rays within the Galaxy. Furthermore, they also demonstrate excellent agreement with three more recent data-driven astrophysical models. These findings suggest that the LCm parameter could be a valuable tool for forthcoming measurements of the LHAASO experiment to enhance our knowledge about the origin and acceleration mechanisms of cosmic rays.

Observation of Gamma Rays up to 320 TeV from the Middle-aged TeV Pulsar Wind Nebula HESS J1849−000

Gamma rays from HESS J1849−000, a middle-aged TeV pulsar wind nebula (PWN), are observed by the Tibet air shower array and the muon detector array. The detection significance of gamma rays reaches 4.0σ and 4.4σ levels above 25 TeV and 100 TeV, respectively, in units of the Gaussian standard deviation σ. The energy spectrum measured between 40 TeV < E < 320 TeV for the first time is described with a simple power-law function of . The gamma-ray energy spectrum from the sub-TeV (E < 1 TeV) to sub-PeV (100 TeV < E < 1 PeV) ranges, including the results of previous studies, can be modeled with the leptonic scenario, i.e., inverse Compton scattering by high-energy electrons accelerated by the PWN of PSR J1849−0001. On the other hand, the gamma-ray energy spectrum can also be modeled with the hadronic scenario in which gamma rays are generated from the decay of neutral pions produced by collisions between accelerated cosmic-ray protons and the ambient molecular cloud found in the gamma-ray-emitting region. The cutoff energy of cosmic-ray protons E p,cut is estimated as , suggesting that protons are accelerated up to the PeV energy range. Our study thus proposes that HESS J1849−000 should be further investigated as a new candidate as a Galactic PeV cosmic-ray accelerator, or “PeVatron.”

The Galactic Population of Pulsar Wind Nebulae and the Contribution of Its Unresolved Component to the Diffuse High-Energy Gamma-ray Emission

In this work, we provide a phenomenological description of the population of galactic TeV pulsar wind nebulae (PWNe) based on suitable assumptions for their space and luminosity distribution. We constrain the general features of this population by assuming that it accounts for the majority of bright sources observed by H.E.S.S. Namely, we determine the maximal luminosity and fading time of PWNe (or, equivalently, the initial period and magnetic field of the pulsar powering the observed emission) by performing a statistical analysis of bright sources in the H.E.S.S. galactic plane survey. This allows us to estimate the total luminosity and flux produced by galactic TeV PWNe. We also evaluate the cumulative emission from PWNe that cannot be resolved by H.E.S.S., showing that this contribution can be as large as ∼40% of the total flux from resolved sources. We argue that also in the GeV domain, a relevant fraction of this population cannot be resolved by Fermi-LAT, providing a non-negligible contribution to the large-scale diffuse emission in the inner galaxy. This additional component could naturally account for a large part of the spectral index variation observed by Fermi-LAT, weakening the evidence of cosmic ray spectral hardening in the inner galaxy. Finally, the same result is obtained for PeV energy, for which the sum of the diffuse component, due to unresolved PWNe, and the truly diffuse emission well saturates the recent Tibet AS-γ data, without the need to introduce a progressive hardening of the cosmic-ray spectrum toward the galactic centre.

Present and future constraints on flavor-dependent long-range interactions of high-energy astrophysical neutrinos

The discovery of new, flavor-dependent neutrino interactions would provide compelling evidence of physics beyond the Standard Model. We focus on interactions generated by the anomaly-free, gauged, abelian lepton-number symmetries, specifically Le – Lμ, Le – Lτ , and Lμ – Lτ , that introduce a new matter potential sourced by electrons and neutrons, potentially impacting neutrino flavor oscillations. We revisit, revamp, and improve the constraints on these interactions that can be placed via the flavor composition of the diffuse flux of high-energy astrophysical neutrinos, with TeV–PeV energies, i.e., the proportion of νe, νμ, and ντ in the flux. Because we consider mediators of these new interactions to be ultra-light, lighter than 10−10 eV, the interaction range is ultra-long, from km to Gpc, allowing vast numbers of electrons and neutrons in celestial bodies and the cosmological matter distribution to contribute to this new potential. We leverage the present-day and future sensitivity of high-energy neutrino telescopes and of oscillation experiments to estimate the constraints that could be placed on the coupling strength of these interactions. We find that, already today, the IceCube neutrino telescope demonstrates potential to constrain flavor-dependent long-range interactions significantly better than existing constraints, motivating further analysis. We also estimate the improvement in the sensitivity due to the next-generation neutrino telescopes such as IceCube-Gen2, Baikal-GVD, KM3NeT, P-ONE, and TAMBO.