In this paper, we report the detection of the very-high-energy (VHE, 100 GeV < E < 100 TeV) and ultra-high-energy (UHE, E > 100 TeV) γ-ray emissions from the direction of the young star-forming region W43, observed by the Large High Altitude Air Shower Observation (LHAASO). The extended γ-ray source was detected with a significance of ∼16 σ by KM2A and ∼17 σ by WCDA, respectively. The angular extension of this γ-ray source is about 0.5 degrees, corresponding to a physical size of about 50 pc. We discuss the origin of the γ-ray emission and possible cosmic ray acceleration in the W43 region using multi-wavelength data. Our findings suggest that W43 is likely another young star cluster capable of accelerating cosmic rays (CRs) to at least several hundred TeV.

We present results over an 11-year Solar cycle of cosmic antiprotons based on 1.1×10^{6} events in the rigidity range from 1.00 to 41.9 GV. The p[over ¯] fluxes exhibit distinct properties. The magnitude of the p[over ¯] flux temporal variation is significantly smaller than those of p, e^{-}, and e^{+}. A hysteresis between the p[over ¯] fluxes and the p fluxes is observed, whereas the p[over ¯] and e^{-} fluxes show a linear correlation. With a model-independent analysis, we found a universal relation between the shape of the rigidity spectrum and the magnitude of flux temporal variation over an 11-year Solar cycle for both positively and negatively charged particles. The simultaneous results on p[over ¯] and p, e^{-}, and e^{+} provide unique information for understanding particle transport in the Solar System as a function of mass, charge, and spectral shape.

We report the properties of precision time structures of cosmic nuclei He, Li, Be, B, C, N, and O fluxes over an 11-year solar cycle from May 2011 to November 2022 in the rigidity range from 1.92 to 60.3 GV. The nuclei fluxes show similar but not identical time variations with amplitudes decreasing with increasing rigidity. In particular, below 3.64 GV the Li, Be, and B fluxes, and below 2.15 GV the C, N, and O fluxes, are significantly less affected by solar modulation than the He flux. We observe that these differences in solar modulation are linearly correlated with the differences in the spectral indices of the cosmic nuclei fluxes. This shows, in a model-independent way, that solar modulation of galactic cosmic nuclei depends on their spectral shape. In addition, solar modulation differences due to nuclei velocity dependence on the mass-to-charge ratio (A/Z) are not observed.

We report the detection of an extended very-high-energy (VHE) γ-ray source coincident with the location of middle-aged (62.4 kyr) pulsar PSR J0248+6021, by using the LHAASO-WCDA data of live 796 d and LHAASO-KM2A data of live 1216 d. A significant excess of γ-ray induced showers is observed both by WCDA in energy bands of 1–25 TeV and KM2A in energy bands of >25 TeV with 7.3σ and 13.5σ, respectively. The best-fit position derived through WCDA data is R.A. = 42.06° ± 0.12° and Dec. = 60.24° ± 0.13° with an extension of 0.69°±0.15° and that of the KM2A data is R.A.= 42.29° ± 0.13° and Dec. = 60.38° ± 0.07° with an extension of 0.37° ±0.07°. No clear extended multiwavelength counterpart of this LHAASO source has been found from the radio band to the GeV band. The most plausible explanation of the VHE γ-ray emission is the inverse Compton process of highly relativistic electrons and positrons injected by the pulsar. These electrons/positrons are hypothesized to be either confined within the pulsar wind nebula or to have already escaped into the interstellar medium, forming a pulsar halo.

The ultra-high-energy (UHE) gamma-ray source 1LHAASO J0007+7303u is positionally associated with the composite SNR CTA1 that is located at high Galactic Latitude b ≈ 10.5°. This provides a rare opportunity to spatially resolve the component of the pulsar wind nebula (PWN) and supernova remnant (SNR) at UHE. This paper conducted a dedicated data analysis of 1LHAASO J0007+7303u using the data collected from December 2019 to July 2023. This source is well detected with significances of 21σ and 17σ at 8–100 TeV and >100 TeV, respectively. The corresponding extensions are determined to be 0.23°±0.03° and 0.17°±0.03°. The emission is proposed to originate from the relativistic electrons accelerated within the PWN of PSR J0007+7303. The energy spectrum is well described by a power-law with an exponential cutoff function dN/dE=(42.4± 4.1)(E 20 TeV)^-2.31± 0.11 exp(-E 110± 25 TeV) TeV−1 cm−2 s−1 in the energy range from 8 to 300 TeV, implying a steady-state parent electron spectrum dN_e/dE_e∝ (E_e 100 TeV)^-3.13± 0.16 exp[(-E_e373± 70 TeV)^2] at energies above ≈ 50 TeV. The cutoff energy of the electron spectrum is roughly equal to the expected current maximum energy of particles accelerated at the PWN terminal shock. Combining the X-ray and gamma-ray emission, the current space-averaged magnetic field can be limited to ≈ 4.5 µG. To satisfy the multi-wavelength spectrum and the γ-ray extensions, the transport of relativistic particles within the PWN is likely dominated by the advection process under the free-expansion phase assumption.

The Water Cherenkov Detector Array (WCDA) is one of the components of Large High Altitude Air Shower Observatory (LHAASO) and can monitor any sources over two-thirds of the sky for up to 7 h per day with >98 per cent duty cycle. In this work, we report the detection of two outbursts of the Fanaroff-Riley I radio galaxy NGC 1275 that were detected by LHAASO-WCDA between 2022 November and 2023 January with statistical significance of 5.2 sigma and 8.3 sigma. The observed spectral energy distribution in the range from 500 GeV to 3 TeV is fitted by a power law with the best-fitting spectral index of alpha = -3.37 +/- 0.52 and -3.35 +/- 0. 29, respectively. The outburst flux above 0.5 TeV was (4.55 +/- 4.21) x10(-11) cm(-2) s(-1) and (3.45 +/- 1.78) x10(-11) cm(-2) s(-1), corresponding to 60 per cent and 45 per cent of Crab Nebula flux, respectively. Variation analysis reveals the variability time-scale of days at the TeV energy band. A simple test by one-zone synchrotron self-Compton model reproduces the data in the gamma-ray band well.

The first source catalog of the Large High Altitude Air Shower Observatory (LHAASO) reported the detection of a very high energy gamma-ray source, 1LHAASO J1219+2915. This Letter presents a further detailed study of the spectral and temporal behavior of this pointlike source. The best-fit position of the TeV source (R.A. = 185.°05 ± 0.°04, decl. = 29.°25 ± 0.°03) is compatible with NGC 4278 within ∼0.°03. Variation analysis shows an indication of variability on a timescale of a few months in the TeV band, which is consistent with low-frequency observations. Based on these observations, we report the detection of TeV γ -ray emissions from this low-luminosity active galactic nucleus. The observation by LHAASO's Water Cherenkov Detector Array during the active period has a significance level of 8.8 σ with a best-fit photon spectral index Γ = 2.56 ± 0.14 and a flux f 1–10 TeV = (7.0 ± 1.1 sta ± 0.35 syst ) × 10 −13 photons cm −2 s −1 , or approximately 5% of the Crab Nebula. The discovery of VHE gamma-ray emission from NGC 4278 indicates that compact, weak radio jets can efficiently accelerate particles and emit TeV photons.

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

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

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-gamma >= 100 TeV) gamma-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 nearby radio galaxy M87 is a very high-energy (VHE) gamma-ray emitter established by observations with ground-based gamma-ray detectors. Here we report the long-term monitoring of M87 from 2021 to 2024 with the Large High Altitude Air Shower Observatory (LHAASO). M87 has been detected by LHAASO with a statistical significance similar to 9 sigma. The observed energy spectrum extends to 20 TeV, with a possible hardening at similar to 20 TeV and then a clear softening at higher energies. Assuming that the intrinsic spectrum is described by a single power law up to 20 TeV, a tight upper bound on the extragalactic background light intensity is obtained. A strong VHE flare lasting 8 days, with a rise time of tau(rise)(r) = 1.05 +/- 0.49 days and decay time of tau(decay)(d) = 2.17 +/- 0.58 days, was found in early 2022. A possible GeV flare is seen also in Fermi Large Area Telescope data during the VHE flare period. The variability time as short as 1 day seen in the LHAASO data suggests a compact emission region with a size of similar to 3 x 10(15) delta cm (delta being the Doppler factor of the emitting region), corresponding to a few Schwarzschild radii of the central supermassive black hole in M87. The continuous monitoring of the source reveals a duty cycle of similar to 1% for VHE flares with a flux above 10(-11) erg cm(-2) s(-1).

The full array of the Large High Altitude Air Shower Observatory（LHAASO） has been in operation since July 2021. For its kilometer-square array（KM2A）, we optimized the selection criteria for very high and ultrahigh energy γ-rays using data collected from August 2021 to August 2022, resulting in an improvement in significance of the detection in the Crab Nebula of approximately 15%, compared with that of previous cuts. With the implementation of these new selection criteria, the angular resolution was also significantly improved by approximately 10% at tens of Te V. Other aspects of the full KM2A array performance, such as the pointing error, were also calibrated using the Crab Nebula. The resulting energy spectrum of the Crab Nebula in the energy range of 10-1000Te V are well fitted by a log-parabola model, which is consistent with the previous results from LHAASO and other experiments.

The measurement of shower muons on an event-by-event basis offers a potent tool for conducting ground-based experiments on gamma rays and cosmic rays due to its sensitivity to primary mass and interaction models. In recent years, underground water Cherenkov detectors as large-area muon counters provide the most powerful way of rejecting cosmic ray background when searching for TeV–PeV gamma rays and cosmic ray electrons, an unprecedented rejection power of 104–105 is achieved. Unburied water Cherenkov detectors are widely used in ground-based gamma astronomy experiments, e.g, Milagro, HAWC, LHAASO-WCDA, etc. However, due to the presence of electromagnetic components, their deployment as event-by-event muon counters has encountered considerable challenges. All the experiments mentioned above reconstruct lateral-distribution-function related parameters to tell a gamma from hadrons. In this work, we first developed a method to utilize the WCDA, to specify muon content in each shower with LHAASO KM2A-WCDA synergy and help LHAASO to gain approximately a 37,650-meter-square effective area as a muon counter.

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

Very recently, the Large High Altitude Air Shower Observatory (LHAASO) reported the observation of the very early TeV afterglow of the brightest-of-all-time GRB 221009A, recording the highest photon statistics in the TeV band ever from a gamma-ray burst. We use this unique observation to place stringent constraints on an energy dependence of the speed of light in vacuum, a manifestation of Lorentz invariance violation (LIV) predicted by some quantum gravity (QG) theories. Our results show that the 95% confidence level lower limits on the QG energy scales are $E_{\mathrm{QG},1}>10$ times of the Planck energy $E_\mathrm{Pl}$ for the linear, and $E_{\mathrm{QG},2}>6\times10^{-8}E_\mathrm{Pl}$ for the quadratic LIV effects, respectively. Our limits on the quadratic LIV case improve previous best bounds by factors of 5--7.

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

The attenuation length of the muon content in extensive air showers provides important information regarding the generation and development of air showers. This information can be used not only to improve the description of such showers but also to test fundamental models of hadronic interactions. Using data from the LHAASO-KM2A experiment, the development of the muon content in high-energy air showers was studied. The attenuation length of muon content in the air showers was measured from experimental data in the energy range from 0.3 to 30 PeV using the constant intensity cut method. By comparing the attenuation length of the muon content with predictions from high-energy hadronic interaction models (QGSJET-II-04, SIBYLL 2.3d, and EPOS-LHC), it is evident that LHAASO results are significantly shorter than those predicted by the first two models (QGSJET-II-04 and SIBYLL 2.3d) but relatively close to those predicted by the third model (EPOS-LHC). Thus, the LHAASO data favor the EPOS-LHC model over the other two models. The three interaction models confirmed an increasing trend in the attenuation length as the cosmic-ray energy increases.

We report the detection of a γ-ray bubble spanning at least 100deg^2 in ultra high energy (UHE) up to a few PeV in the direction of the star-forming region Cygnus X, implying the presence Super PeVatron(s) accelerating protons to at least 10 PeV. A log-parabola form with the photon index Γ(E) = (2.71 ± 0.02) + (0.11 ± 0.02) ×log_10 (E/10 TeV) is found fitting the gamma-ray energy spectrum of the bubble well. UHE sources, `hot spots' correlated with very massive molecular clouds, and a quasi-spherical amorphous γ-ray emitter with a sharp central brightening are observed in the bubble. In the core of ∼ 0.5^∘, spatially associating with a region containing massive OB association (Cygnus OB2) and a microquasar (Cygnus X-3), as well as previously reported multi-TeV sources, an enhanced concentration of UHE γ-rays are observed with 2 photons at energies above 1 PeV. The general feature of the bubble, the morphology and the energy spectrum, are reasonably reproduced by the assumption of a particle accelerator in the core, continuously injecting protons into the ambient medium.

One of the main scientific objectives of the Large High Altitude Air Shower Observatory (LHAASO) is to perform accurate measurements of the energy spectra for different cosmic ray masses, from a few TeV to 100 PeV. As one of the main sub-arrays of LHAASO, the Wide Field-of-View Cherenkov Telescope Array (WFCTA), which consists of 18 telescopes, can help in achieving this goal. The pointing accuracy of each telescope is crucial for reconstructing energy and determining mass-sensitive parameters. In this study, a method for absolute pointing calibration using ultraviolet bright stars was established. The proposed method can achieve a pointing accuracy of approximately 0.02° by using more than five stars. When more stars are used in the telescope's field of view, the accuracy of the pointing calibration can be improved, e.g, a pointing accuracy of 0.01° can be achieved when using more than 15 stars. A high-precision inclinometer with a monitor resolution of 0.003° was installed on the camera to monitor the zenith direction of the telescope at every second. After calibration using bright stars, the absolute pointing accuracy of the inclinometer was 0.02°.

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.