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.

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 of the bright GRB 221009A, which serendipitously occurred within the instrument 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, then rose to a peak about 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 ~0.8°. This is consistent with the core of a structured jet and could explain the high isotropic energy of this GRB.

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.

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.

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 × 10³⁶ 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.

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.

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

The extension of the cosmic-ray spectrum beyond 1 petaelectronvolt (PeV; 10¹⁵ 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 rays¹. 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. ²). Although γ-rays with energies slightly higher than 0.1 PeV have been reported from a few objects in the Galactic plane^3-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.

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.

Measurement of the high-energy gamma-ray emission from the Moon with the Fermi Large Area Telescope

We have measured the gamma-ray emission spectrum of the Moon using the data collected by the Large Area Telescope onboard the Fermi satellite during its first seven years of operation, in the energy range from 30 MeV up to a few GeV. We have also studied the time evolution of the flux, finding a correlation with the solar activity. We have developed a full Monte Carlo simulation describing the interactions of cosmic rays with the lunar surface. The results of the present analysis can be explained in the framework of this model, where the production of gamma rays is due to the interactions of cosmic-ray proton and helium nuclei with the surface of the Moon. Finally, we have used our simulation to derive the cosmic-ray proton and helium spectra near Earth from the Moon gamma-ray data.

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.

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.

Fermi LARGE AREA TELESCOPE OBSERVATIONS OF BLAZAR 3C 279 OCCULTATIONS BY THE SUN

Observations of occultations of bright γ-ray sources by the Sun may reveal predicted pair halos around blazars and/or new physics, such as, e.g., hypothetical light dark matter particles-axions. We use Fermi Gamma-Ray Space Telescope (Fermi) data to analyze four occultations of blazar 3C 279 by the Sun on October 8 each year from 2008 to 2011. A combined analysis of the observations of these occultations allows a point-like source at the position of 3C 279 to be detected with significance of ≈3σ, but does not reveal any significant excess over the flux expected from the quiescent Sun. The likelihood ratio test rules out complete transparency of the Sun to the blazar γ-ray emission at a 3σ confidence level.

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.

Gamma-ray emission concurrent with the nova in the symbiotic binary V407 Cygni

Novae are thermonuclear explosions on a white dwarf surface fueled by mass accreted from a companion star. Current physical models posit that shocked expanding gas from the nova shell can produce x-ray emission, but emission at higher energies has not been widely expected. Here, we report the Fermi Large Area Telescope detection of variable gamma-ray emission (0.1 to 10 billion electron volts) from the recently detected optical nova of the symbiotic star V407 Cygni. We propose that the material of the nova shell interacts with the dense ambient medium of the red giant primary and that particles can be accelerated effectively to produce pi(0) decay gamma-rays from proton-proton interactions. Emission involving inverse Compton scattering of the red giant radiation is also considered and is not ruled out.

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.

Detection of high-energy gamma-ray emission from the globular cluster 47 Tucanae with Fermi

We report the detection of gamma-ray emissions above 200 megaelectron volts at a significance level of 17sigma from the globular cluster 47 Tucanae, using data obtained with the Large Area Telescope onboard the Fermi Gamma-ray Space Telescope. Globular clusters are expected to emit gamma rays because of the large populations of millisecond pulsars that they contain. The spectral shape of 47 Tucanae is consistent with gamma-ray emission from a population of millisecond pulsars. The observed gamma-ray luminosity implies an upper limit of 60 millisecond pulsars present in 47 Tucanae.

Detection of 16 gamma-ray pulsars through blind frequency searches using the Fermi LAT

Pulsars are rapidly rotating, highly magnetized neutron stars emitting radiation across the electromagnetic spectrum. Although there are more than 1800 known radio pulsars, until recently only seven were observed to pulse in gamma rays, and these were all discovered at other wavelengths. The Fermi Large Area Telescope (LAT) makes it possible to pinpoint neutron stars through their gamma-ray pulsations. We report the detection of 16 gamma-ray pulsars in blind frequency searches using the LAT. Most of these pulsars are coincident with previously unidentified gamma-ray sources, and many are associated with supernova remnants. Direct detection of gamma-ray pulsars enables studies of emission mechanisms, population statistics, and the energetics of pulsar wind nebulae and supernova remnants.

Measurement of the cosmic ray e+ +e- spectrum from 20 GeV to 1 TeV with the Fermi Large Area Telescope

Designed as a high-sensitivity gamma-ray observatory, the Fermi Large Area Telescope is also an electron detector with a large acceptance exceeding 2 m;{2} sr at 300 GeV. Building on the gamma-ray analysis, we have developed an efficient electron detection strategy which provides sufficient background rejection for measurement of the steeply falling electron spectrum up to 1 TeV. Our high precision data show that the electron spectrum falls with energy as E-3.0 and does not exhibit prominent spectral features. Interpretations in terms of a conventional diffusive model as well as a potential local extra component are briefly discussed.

The Fermi Gamma-Ray Space Telescope discovers the pulsar in the young galactic supernova remnant CTA 1

Energetic young pulsars and expanding blast waves [supernova remnants (SNRs)] are the most visible remains after massive stars, ending their lives, explode in core-collapse supernovae. The Fermi Gamma-Ray Space Telescope has unveiled a radio quiet pulsar located near the center of the compact synchrotron nebula inside the supernova remnant CTA 1. The pulsar, discovered through its gamma-ray pulsations, has a period of 316.86 milliseconds and a period derivative of 3.614 x 10(-13) seconds per second. Its characteristic age of 10(4) years is comparable to that estimated for the SNR. We speculate that most unidentified Galactic gamma-ray sources associated with star-forming regions and SNRs are such young pulsars.