Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 γ-ray Galactic sources.

The extension of the cosmic-ray spectrum beyond 1 petaelectronvolt (PeV; 1015 electronvolts) indicates the existence of the so-called PeVatrons-cosmic-ray factories that accelerate particles to PeV energies. We need to locate and identify such objects to find the origin of Galactic cosmic rays1. The principal signature of both electron and proton PeVatrons is ultrahigh-energy (exceeding 100 TeV) γ radiation. Evidence of the presence of a proton PeVatron has been found in the Galactic Centre, according to the detection of a hard-spectrum radiation extending to 0.04 PeV (ref. 2). Although γ-rays with energies slightly higher than 0.1 PeV have been reported from a few objects in the Galactic plane3-6, unbiased identification and in-depth exploration of PeVatrons requires detection of γ-rays with energies well above 0.1 PeV. Here we report the detection of more than 530 photons at energies above 100 teraelectronvolts and up to 1.4 PeV from 12 ultrahigh-energy γ-ray sources with a statistical significance greater than seven standard deviations. Despite having several potential counterparts in their proximity, including pulsar wind nebulae, supernova remnants and star-forming regions, the PeVatrons responsible for the ultrahigh-energy γ-rays have not yet been firmly localized and identified (except for the Crab Nebula), leaving open the origin of these extreme accelerators.

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

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±0.05±0.15 PeV. Below the knee, the spectral index is found to be -2.7413±0.0004±0.0050, while above the knee, it is -3.128±0.005±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±0.0003±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.

Measurement of Ultra-High-Energy Diffuse Gamma-Ray Emission of the Galactic Plane from 10 TeV to 1 PeV with LHAASO-KM2A.

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 Letter, 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°10 TeV). The energy spectrum in the inner Galaxy regions can be described by a power-law function with an index of -2.99±0.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 ∼3 than the prediction. A similar spectrum with an index of -2.99±0.07 is found in the outer Galaxy region, and the absolute flux for 10≲E≲60 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.

Exploring Lorentz Invariance Violation from Ultrahigh-Energy γ Rays Observed by LHAASO.

Recently, the LHAASO Collaboration published the detection of 12 ultrahigh-energy γ-ray sources above 100 TeV, with the highest energy photon reaching 1.4 PeV. The first detection of PeV γ rays from astrophysical sources may provide a very sensitive probe of the effect of the Lorentz invariance violation (LIV), which results in decay of high-energy γ rays in the superluminal scenario and hence a sharp cutoff of the energy spectrum. Two highest energy sources are studied in this work. No signature of the existence of the LIV is found in their energy spectra, and the lower limits on the LIV energy scale are derived. Our results show that the first-order LIV energy scale should be higher than about 10^{5} times the Planck scale M_{Pl} and that the second-order LIV scale is >10^{-3}M_{Pl}. Both limits improve by at least one order of magnitude the previous results.

Constraints on Heavy Decaying Dark Matter from 570 Days of LHAASO Observations.

The kilometer square array (KM2A) of the large high altitude air shower observatory (LHAASO) aims at surveying the northern γ-ray sky at energies above 10 TeV with unprecedented sensitivity. γ-ray observations have long been one of the most powerful tools for dark matter searches, as, e.g., high-energy γ rays could be produced by the decays of heavy dark matter particles. In this Letter, we present the first dark matter analysis with LHAASO-KM2A, using the first 340 days of data from 1/2-KM2A and 230 days of data from 3/4-KM2A. Several regions of interest are used to search for a signal and account for the residual cosmic-ray background after γ/hadron separation. We find no excess of dark matter signals, and thus place some of the strongest γ-ray constraints on the lifetime of heavy dark matter particles with mass between 10^{5} and 10^{9} GeV. Our results with LHAASO are robust, and have important implications for dark matter interpretations of the diffuse astrophysical high-energy neutrino emission.

Extended Very-High-Energy Gamma-Ray Emission Surrounding PSR J0622+3749 Observed by LHAASO-KM2A.

We report the discovery of an extended very-high-energy (VHE) gamma-ray source around the location of the middle-aged (207.8 kyr) pulsar PSR J0622+3749 with the Large High-Altitude Air Shower Observatory (LHAASO). The source is detected with a significance of 8.2σ for E>25 TeV assuming a Gaussian template. The best-fit location is (right ascension, declination) =(95.47°±0.11°,37.92°±0.09°), and the extension is 0.40°±0.07°. The energy spectrum can be described by a power-law spectrum with an index of -2.92±0.17_{stat}±0.02_{sys}. No clear extended multiwavelength counterpart of the LHAASO source has been found from the radio to sub-TeV bands. The LHAASO observations are consistent with the scenario that VHE electrons escaped from the pulsar, diffused in the interstellar medium, and scattered the interstellar radiation field. If interpreted as the pulsar halo scenario, the diffusion coefficient, inferred for electrons with median energies of ∼160 TeV, is consistent with those obtained from the extended halos around Geminga and Monogem and much smaller than that derived from cosmic ray secondaries. The LHAASO discovery of this source thus likely enriches the class of so-called pulsar halos and confirms that high-energy particles generally diffuse very slowly in the disturbed medium around pulsars.

Unresolved Gamma-Ray Sky through its Angular Power Spectrum.

The gamma-ray sky has been observed with unprecedented accuracy in the last decade by the Fermi -large area telescope (LAT), allowing us to resolve and understand the high-energy Universe. The nature of the remaining unresolved emission [unresolved gamma-ray background (UGRB)] below the LAT source detection threshold can be uncovered by characterizing the amplitude and angular scale of the UGRB fluctuation field. This Letter presents a measurement of the UGRB autocorrelation angular power spectrum based on eight years of Fermi-LAT Pass 8 data products. The analysis is designed to be robust against contamination from resolved sources and noise systematics. The sensitivity to subthreshold sources is greatly enhanced with respect to previous measurements. We find evidence (with ∼3.7σ significance) that the scenario in which two classes of sources contribute to the UGRB signal is favored over a single class. A double power law with exponential cutoff can explain the anisotropy energy spectrum well, with photon indices of the two populations being 2.55±0.23 and 1.86±0.15.

Search for Spectral Irregularities due to Photon-Axionlike-Particle Oscillations with the Fermi Large Area Telescope.

We report on the search for spectral irregularities induced by oscillations between photons and axionlike-particles (ALPs) in the γ-ray spectrum of NGC 1275, the central galaxy of the Perseus cluster. Using 6 years of Fermi Large Area Telescope data, we find no evidence for ALPs and exclude couplings above 5×10^{-12} GeV^{-1} for ALP masses 0.5≲m_{a}≲5 neV at 95% confidence. The limits are competitive with the sensitivity of planned laboratory experiments, and, together with other bounds, strongly constrain the possibility that ALPs can reduce the γ-ray opacity of the Universe.

Resolving the Extragalactic γ-Ray Background above 50 GeV with the Fermi Large Area Telescope.

The Fermi Large Area Telescope (LAT) Collaboration has recently released a catalog of 360 sources detected above 50 GeV (2FHL). This catalog was obtained using 80 months of data re-processed with Pass 8, the newest event-level analysis, which significantly improves the acceptance and angular resolution of the instrument. Most of the 2FHL sources at high Galactic latitude are blazars. Using detailed Monte Carlo simulations, we measure, for the first time, the source count distribution, dN/dS, of extragalactic γ-ray sources at E>50 GeV and find that it is compatible with a Euclidean distribution down to the lowest measured source flux in the 2FHL (∼8×10^{-12} ph cm^{-2} s^{-1}). We employ a one-point photon fluctuation analysis to constrain the behavior of dN/dS below the source detection threshold. Overall, the source count distribution is constrained over three decades in flux and found compatible with a broken power law with a break flux, S_{b}, in the range [8×10^{-12},1.5×10^{-11}] ph cm^{-2} s^{-1} and power-law indices below and above the break of α_{2}∈[1.60,1.75] and α_{1}=2.49±0.12, respectively. Integration of dN/dS shows that point sources account for at least 86_{-14}^{+16}% of the total extragalactic γ-ray background. The simple form of the derived source count distribution is consistent with a single population (i.e., blazars) dominating the source counts to the minimum flux explored by this analysis. We estimate the density of sources detectable in blind surveys that will be performed in the coming years by the Cherenkov Telescope Array.

Searching for Dark Matter Annihilation from Milky Way Dwarf Spheroidal Galaxies with Six Years of Fermi Large Area Telescope Data.

The dwarf spheroidal satellite galaxies (dSphs) of the Milky Way are some of the most dark matter (DM) dominated objects known. We report on γ-ray observations of Milky Way dSphs based on six years of Fermi Large Area Telescope data processed with the new Pass8 event-level analysis. None of the dSphs are significantly detected in γ rays, and we present upper limits on the DM annihilation cross section from a combined analysis of 15 dSphs. These constraints are among the strongest and most robust to date and lie below the canonical thermal relic cross section for DM of mass ≲100 GeV annihilating via quark and τ-lepton channels.

A limit on the variation of the speed of light arising from quantum gravity effects.

A cornerstone of Einstein's special relativity is Lorentz invariance-the postulate that all observers measure exactly the same speed of light in vacuum, independent of photon-energy. While special relativity assumes that there is no fundamental length-scale associated with such invariance, there is a fundamental scale (the Planck scale, l(Planck) approximately 1.62 x 10(-33) cm or E(Planck) = M(Planck)c(2) approximately 1.22 x 10(19) GeV), at which quantum effects are expected to strongly affect the nature of space-time. There is great interest in the (not yet validated) idea that Lorentz invariance might break near the Planck scale. A key test of such violation of Lorentz invariance is a possible variation of photon speed with energy. Even a tiny variation in photon speed, when accumulated over cosmological light-travel times, may be revealed by observing sharp features in gamma-ray burst (GRB) light-curves. Here we report the detection of emission up to approximately 31 GeV from the distant and short GRB 090510. We find no evidence for the violation of Lorentz invariance, and place a lower limit of 1.2E(Planck) on the scale of a linear energy dependence (or an inverse wavelength dependence), subject to reasonable assumptions about the emission (equivalently we have an upper limit of l(Planck)/1.2 on the length scale of the effect). Our results disfavour quantum-gravity theories in which the quantum nature of space-time on a very small scale linearly alters the speed of light.

Inferred cosmic-ray spectrum from Fermi large area telescope γ-ray observations of Earth's limb.

Recent accurate measurements of cosmic-ray (CR) species by ATIC-2, CREAM, and PAMELA reveal an unexpected hardening in the proton and He spectra above a few hundred GeV, a gradual softening of the spectra just below a few hundred GeV, and a harder spectrum of He compared to that of protons. These newly discovered features may offer a clue to the origin of high-energy CRs. We use the Fermi Large Area Telescope observations of the γ-ray emission from Earth's limb for an indirect measurement of the local spectrum of CR protons in the energy range ∼90 GeV-6 TeV (derived from a photon energy range 15 GeV-1 TeV). Our analysis shows that single power law and broken power law spectra fit the data equally well and yield a proton spectrum with index 2.68±0.04 and 2.61±0.08 above ∼200 GeV, respectively.

A tera-electron volt afterglow from a narrow jet in an extremely bright gamma-ray burst.

Some gamma-ray bursts (GRBs) have a tera-electron volt (TeV) afterglow, but the early onset of this has not been observed. We report observations with the Large High Altitude Air Shower Observatory (LHAASO) of the bright GRB 221009A, which serendipitously occurred within the instrument's field of view. More than 64,000 photons >0.2 TeV were detected within the first 3000 seconds. The TeV flux began several minutes after the GRB trigger and then rose to a peak ~10 seconds later. This was followed by a decay phase, which became more rapid ~650 seconds after the peak. We interpret the emission using a model of a relativistic jet with half-opening angle of ~0.8°. This is consistent with the core of a structured jet and could explain the high isotropic energy of this GRB.

Measurement of separate cosmic-ray electron and positron spectra with the fermi large area telescope.

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Telescope. Because the instrument does not have an onboard magnet, we distinguish the two species by exploiting Earth's shadow, which is offset in opposite directions for opposite charges due to Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using two different methods that produce consistent results. We report the electron-only spectrum, the positron-only spectrum, and the positron fraction between 20 and 200 GeV. We confirm that the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 and 200 GeV are consistent with a fraction that is continuing to rise with energy.

Constraining dark matter models from a combined analysis of Milky Way satellites with the Fermi Large Area Telescope.

Satellite galaxies of the Milky Way are among the most promising targets for dark matter searches in gamma rays. We present a search for dark matter consisting of weakly interacting massive particles, applying a joint likelihood analysis to 10 satellite galaxies with 24 months of data of the Fermi Large Area Telescope. No dark matter signal is detected. Including the uncertainty in the dark matter distribution, robust upper limits are placed on dark matter annihilation cross sections. The 95% confidence level upper limits range from about 10(-26) cm3 s(-1) at 5 GeV to about 5×10(-23) cm3 s(-1) at 1 TeV, depending on the dark matter annihilation final state. For the first time, using gamma rays, we are able to rule out models with the most generic cross section (∼3×10(-26) cm3 s(-1) for a purely s-wave cross section), without assuming additional boost factors.

Peta-electron volt gamma-ray emission from the Crab Nebula.

The Crab Nebula is a bright source of gamma rays powered by the Crab Pulsar's rotational energy through the formation and termination of a relativistic electron-positron wind. We report the detection of gamma rays from this source with energies from 5 × 10-4 to 1.1 peta-electron volts with a spectrum showing gradual steepening over three energy decades. The ultrahigh-energy photons imply the presence of a peta-electron volt electron accelerator (a pevatron) in the nebula, with an acceleration rate exceeding 15% of the theoretical limit. We constrain the pevatron's size between 0.025 and 0.1 parsecs and the magnetic field to ≈110 microgauss. The production rate of peta-electron volt electrons, 2.5 × 1036 ergs per second, constitutes 0.5% of the pulsar spin-down luminosity, although we cannot exclude a contribution of peta-electron volt protons to the production of the highest-energy gamma rays.

Spectrum of the isotropic diffuse gamma-ray emission derived from first-year Fermi Large Area Telescope data.

We report on the first Fermi Large Area Telescope (LAT) measurements of the so-called "extragalactic" diffuse gamma-ray emission (EGB). This component of the diffuse gamma-ray emission is generally considered to have an isotropic or nearly isotropic distribution on the sky with diverse contributions discussed in the literature. The derivation of the EGB is based on detailed modeling of the bright foreground diffuse Galactic gamma-ray emission, the detected LAT sources, and the solar gamma-ray emission. We find the spectrum of the EGB is consistent with a power law with a differential spectral index gamma = 2.41 +/- 0.05 and intensity I(>100 MeV) = (1.03 +/- 0.17) x 10(-5) cm(-2) s(-1) sr(-1), where the error is systematics dominated. Our EGB spectrum is featureless, less intense, and softer than that derived from EGRET data.

Fermi large area telescope search for photon lines from 30 to 200 GeV and dark matter implications.

Dark matter (DM) particle annihilation or decay can produce monochromatic gamma rays readily distinguishable from astrophysical sources. gamma-ray line limits from 30 to 200 GeV obtained from 11 months of Fermi Large Area Space Telescope data from 20-300 GeV are presented using a selection based on requirements for a gamma-ray line analysis, and integrated over most of the sky. We obtain gamma-ray line flux upper limits in the range 0.6-4.5x10{-9} cm{-2} s{-1}, and give corresponding DM annihilation cross-section and decay lifetime limits. Theoretical implications are briefly discussed.

Fermi large area telescope measurements of the diffuse gamma-ray emission at intermediate galactic latitudes.

The diffuse galactic gamma-ray emission is produced by cosmic rays (CRs) interacting with the interstellar gas and radiation field. Measurements by the Energetic Gamma-Ray Experiment Telescope (EGRET) instrument on the Compton Gamma-Ray Observatory indicated excess gamma-ray emission greater, > or approximately equal to 1 GeV relative to diffuse galactic gamma-ray emission models consistent with directly measured CR spectra (the so-called "EGRET GeV excess"). The Large Area Telescope (LAT) instrument on the Fermi Gamma-Ray Space Telescope has measured the diffuse gamma-ray emission with improved sensitivity and resolution compared to EGRET. We report on LAT measurements for energies 100 MeV to 10 GeV and galactic latitudes 10 degrees < or = |b| < or = 20 degrees. The LAT spectrum for this region of the sky is well reproduced by a diffuse galactic gamma-ray emission model that is consistent with local CR spectra and inconsistent with the EGRET GeV excess.

Investigating the Nature of Late-Time High-Energy GRB Emission Through Joint Fermi /Swift Observations.

We use joint observations by the Neil Gehrels Swift X-ray Telescope (XRT) and the Fermi Large Area Telescope (LAT) of gamma-ray burst (GRB) afterglows to investigate the nature of the long-lived high-energy emission observed by Fermi LAT. Joint broadband spectral modeling of XRT and LAT data reveal that LAT non-detections of bright X-ray afterglows are consistent with a cooling break in the inferred electron synchrotron spectrum below the LAT and/or XRT energy ranges. Such a break is sufficient to suppress the high-energy emission so as to be below the LAT detection threshold. By contrast, LAT-detected bursts are best fit by a synchrotron spectrum with a cooling break that lies either between or above the XRT and LAT energy ranges. We speculate that the primary difference between GRBs with LAT afterglow detections and the non-detected population may be in the type of circumstellar environment in which these bursts occur, with late-time LAT detections preferentially selecting GRBs that occur in low wind-like circumburst density profiles. Furthermore, we find no evidence of high-energy emission in the LAT-detected population significantly in excess of the flux expected from the electron synchrotron spectrum fit to the observed X-ray emission. The lack of excess emission at high energies could be due to a shocked external medium in which the energy density in the magnetic field is stronger than or comparable to that of the relativistic electrons behind the shock, precluding the production of a dominant synchrotron self-Compton (SSC) component in the LAT energy range. Alternatively, the peak of the SSC emission could be beyond the 0.1-100 GeV energy range considered for this analysis.