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 (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.

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.

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.

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; 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.

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.

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.

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.

The first pulse of the extremely bright GRB 130427A: a test lab for synchrotron shocks.

Gamma-ray burst (GRB) 130427A is one of the most energetic GRBs ever observed. The initial pulse up to 2.5 seconds is possibly the brightest well-isolated pulse observed to date. A fine time resolution spectral analysis shows power-law decays of the peak energy from the onset of the pulse, consistent with models of internal synchrotron shock pulses. However, a strongly correlated power-law behavior is observed between the luminosity and the spectral peak energy that is inconsistent with curvature effects arising in the relativistic outflow. It is difficult for any of the existing models to account for all of the observed spectral and temporal behaviors simultaneously.

Fermi-LAT observations of the gamma-ray burst GRB 130427A.

The observations of the exceptionally bright gamma-ray burst (GRB) 130427A by the Large Area Telescope aboard the Fermi Gamma-ray Space Telescope provide constraints on the nature of these unique astrophysical sources. GRB 130427A had the largest fluence, highest-energy photon (95 GeV), longest γ-ray duration (20 hours), and one of the largest isotropic energy releases ever observed from a GRB. Temporal and spectral analyses of GRB 130427A challenge the widely accepted model that the nonthermal high-energy emission in the afterglow phase of GRBs is synchrotron emission radiated by electrons accelerated at an external shock.

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.

Detection of the characteristic pion-decay signature in supernova remnants.

Cosmic rays are particles (mostly protons) accelerated to relativistic speeds. Despite wide agreement that supernova remnants (SNRs) are the sources of galactic cosmic rays, unequivocal evidence for the acceleration of protons in these objects is still lacking. When accelerated protons encounter interstellar material, they produce neutral pions, which in turn decay into gamma rays. This offers a compelling way to detect the acceleration sites of protons. The identification of pion-decay gamma rays has been difficult because high-energy electrons also produce gamma rays via bremsstrahlung and inverse Compton scattering. We detected the characteristic pion-decay feature in the gamma-ray spectra of two SNRs, IC 443 and W44, with the Fermi Large Area Telescope. This detection provides direct evidence that cosmic-ray protons are accelerated in SNRs.

Binary millisecond pulsar discovery via gamma-ray pulsations.

Millisecond pulsars, old neutron stars spun up by accreting matter from a companion star, can reach high rotation rates of hundreds of revolutions per second. Until now, all such "recycled" rotation-powered pulsars have been detected by their spin-modulated radio emission. In a computing-intensive blind search of gamma-ray data from the Fermi Large Area Telescope (with partial constraints from optical data), we detected a 2.5-millisecond pulsar, PSR J1311-3430. This unambiguously explains a formerly unidentified gamma-ray source that had been a decade-long enigma, confirming previous conjectures. The pulsar is in a circular orbit with an orbital period of only 93 minutes, the shortest of any spin-powered pulsar binary ever found.