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.

Terrestrial Gamma-ray flashes (TGFs) originating from the Earth's atmosphere, accompanied by thunderstorms and lightning activity, are one of the hot spots in the interdisciplinary of cosmic ray and atmospheric physics. Over the years, satellite experiments have detected thousands of upward TGFa during thunderstorms, while ground-based experiments have observed some downward TGFa. Nowadays, it is widely believed that TGFs accompanying lightning leaders observed by satellite-based and ground-based experiments involve relativistic runaway electron avalanche (RREA) production. Due to triggering the relativistic runaway electron avalanche (RREA) process needing a very large electric field strength and region, it is difficult to study the RREA process through ground-based experiments. In this paper, CORSIKA 7.7410 software package, combined with a vertically uniform electric field model, is adopted to simulate the intensity and energy distribution of RREA electrons in thundercloud with different electric field strengths induced by seed electrons and the secondary electrons in extensive air shower (EAS) from vertical protons with different primary energies. The results show that the number of RREA electrons increases exponentially with the thickness of the thunderclouds increasing, and also increases exponentially with the electric field strength rising. After passing through the atmosphere with an electric field of -3000 V/om and a thickness of 800 m, the number of secondary electrons in RREA process increases by approximately 3x10(4) times. The characteristic length of avalanche (lambda) decreases as the electric field strength increases. When the electric field is -1600 V/com and -3000 V/cm, the lambda is approximately similar to 282 m and similar to 69 m, respectively. The energy spectrum of RREA electrons gradually softens with the increase of layer thickness and strength of electric field, and their average energy increases with the increase of electric field strength, when the thundercloud thickness exceeda 400 m, the mean energy of RREA electrons gradually stabilizes. When secondary particles pass through a thundercloud with an electric field strength of -3000 V/cm and a thickness of 800 m, the mean energy of RREA electrons is approximately 11.7 MeV. Through the Monte Carlo simulations, the RREA process, which is difficult to observe directly in the atmosphere, is successfully simulated. The simulation results provide important information for studying the characteristics of TGF source regions, offer clues for detecting downward TGF in ground-based experiments, and contribute to the research on the triggering mechanism of lightning in the atmosphere. In addition, our simulation results are expected to elucidate the relationship between TGF and lightning activity, promoting interdisciplinary research in the fields of atmospheric physics and cosmic ray physios.

The Large High Altitude Air Shower Observatory (LHAASO) is located at Haizi Mountain (4410 m a. s. l.), Daocheng, Sichuan province, China. Due to its high-altitude location with frequent thunderstorm activities, the LHAASO is suited for studying the effects of near-Earth thunderstorm electric fields on cosmic ray air showers. In this paper, Monte Carlo simulations are performed with CORSIKA and G4KM2A to analyze the flux variations of cosmic ray air showers detected by the kilometer-square array of LHAASO (LHAASO-KM2A) during thunderstorms. The strength, polarity, and layer thickness of the atmospheric electric field (AEF) during thunderstorm are found to be associated with the shower rate variations. The flux of shower events satisfying trigger conditions of the LHAASO-KM2A increases with field intensity, particularly within negative fields, and the enhanced amplitude is more than 5% in -600 V/cm and 12% in AEF of -1000 V/cm, whereas it increases only by 1% and 7% in equivalent positive fields, respectively. While in positive fields ranging from 0 to 400 V/cm, the shower rate decreases with smaller amplitudes. Furthermore, the shower rate increases dramatically with the AEF layer thickness until a certain value, above which the variation trend slows down. The dependence of the trigger rate variation on the primary zenith angle has also been revealed, increasing in lower zenith angle ranges and showing opposite behavior in higher ones. Additionally, we study that the relationship between the trigger rate variations and the primary energies, and find the enhanced amplitude of the shower rate decreases with increasing primary energy. Simultaneously, the shower events with lower primary energy show a significant increase, whereas events with higher primary energy are hardly affected during thunderstorms. These simulation results offer valuable insights into the variation of the trigger rate detected by LHAASO-KM2A during thunderstorms and are beneficial for understanding the acceleration mechanisms of secondary charged particles caused by AEF.

Ultra-high-energy (UHE), exceeding 100 TeV (10^12 electronvolts), γ-rays manifests extreme particle acceleration in astrophysical sources. Recent observations by γ-ray telescopes, particularly by the Large High Altitude Air Shower Observatory (LHAASO), have revealed a few tens of UHE sources, indicating numerous Galactic sources capable of accelerating particles to PeV (10^15 electronvolts) energies. However, discerning the dominant acceleration mechanisms (leptonic versus hadronic), the relative contributions of specific source classes, and the role of particle transport in shaping their observed emission are central goals of modern UHE astrophysics. Here we report the discovery of a giant UHE γ-ray emitter at -17.5° off the Galactic plane - a region where UHE γ-ray sources are rarely found. The emitter exhibits a distinctive asymmetric shape, resembling a giant "Peanut" spanning 0.45° \times 4.6°, indicative of anisotropic particle distribution over a large area. A highly aged millisecond pulsar (MSP) J0218+4232 is the sole candidate accelerator positionally coincident with the Peanut region. Its association with UHE γ-rays extending to 0.7 PeV, if confirmed, would provide the first evidence of a millisecond pulsar powering PeV particles. Such a finding challenges prevailing models, which posit that millisecond pulsars cannot sustain acceleration to PeV energies. The detection reveals fundamental gaps in understanding particle acceleration, cosmic-ray transport, and interstellar magnetic field effects, potentially revealing new PeV accelerator (PeVatron) classes.

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 core-collapse supernova remnant (SNR) Cassiopeia A (Cas A) is one of the brightest galactic radio sources with an angular radius of ∼ 2 . ′ 5 . Although no extension of this source has been detected in the γ -ray band, using more than 1000 days of LHAASO data above ∼0.8 TeV, we find that its spectrum is significantly softer than those obtained with Imaging Air Cherenkov Telescopes (IACTs), and its flux near ∼1 TeV is about 2 times higher. In combination with analyses of more than 16 yr of Fermi-LAT data covering 0.1 GeV–1 TeV, we find that the spectrum above 30 GeV deviates significantly from a single power law and is best described by a smoothly broken power law with a spectral index of 1.90 ± 0.15 stat (3.41 ± 0.19 stat ) below (above) a break energy of 0.63 ± 0.21 stat TeV. Given differences in the angular resolution of LHAASO-WCDA and IACTs, TeV γ -ray emission detected with LHAASO may have a significant contribution from regions surrounding the SNR illuminated by particles accelerated earlier, which, however, are treated as background by IACTs. Detailed modeling can be used to constrain the acceleration processes of TeV particles in the early stage of SNR evolution.

Black holes (BHs), one of the most intriguing objects in the universe, can manifest themselves through electromagnetic radiation initiated by the accretion flow. Some stellar-mass BHs drive relativistic jets when accreting matter from their companion stars, forming microquasars. Non-thermal emission from the radio to teraelectronvolt gamma-ray band has been observed from microquasars, indicating the acceleration of relativistic particles. Here we report detection of four microquasars (SS 433, V4641 Sgr, GRS 1915+105, MAXI J1820+070) of spectra extending to the ultrahigh-energy (UHE; photon energy [Formula: see text] TeV) band, and one microquasar (Cygnus X-1) with a spectrum approaching 100 TeV, using the Large High Altitude Air Shower Observatory. Notably, the total emission associated with SS 433 cannot be interpreted with a single leptonic component. In the UHE band, its emission is in spatial coincidence with a giant atomic cloud, which is consistent with a hadronic origin. An elongated source is discovered from V4641 Sgr with the spectrum continuing up to 800 TeV. The detection of UHE gamma rays demonstrates that accreting BHs and their environments can operate as extremely efficient accelerators of particles up to 1 PeV, suggesting that microquasars are important contributors to Galactic cosmic rays, especially around the 'knee' region.

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.

Large high altitude air shower observatory (LHAASO) is a complex of extensive air shower (EAS) detector arrays, located on the Mt. Haizi (29 degrees 21' N, 100 degrees 08' E) at an altitude of 4410 m a. s. l., Daocheng, Sichuan Province, China. The information about primary cosmic rays can be obtained by using data from secondary particles measured at LHAASO, with photons make up the majority among these secondary particles. During thunderstorms, the atmospheric electric field can affect secondary charged particles (mainly positrons and electrons), thus changing the information of photons on the ground. In this work, Monte Carlo simulations are performed to investigate the effects of near-ground thunderstorm electric fields on cosmic ray secondary photons at LHAASO. A simple model with a vertical and uniform atmospheric electric field in a layer of atmosphere is used in our simulations. During thunderstorms, the number and energy of photons are found to significantly change and strongly depend on the electric field strength. In a field of -1000 V/cm (below the threshold of the relativistic runaway electron avalanche (RREA) process), the number of photons is increased by 23%. Also, the spectrum of photons softens, and the increased number of photons with energy less than 2 MeV exceeds 29%. In an electric field of -1700 V/cm (above the threshold of the RREA process), the number of photons experiences exponential growth, with an increase of 279%. The spectrum of photons becomes softer than that at -1000 V/cm, and the increased number with energy less than 2 MeV is more than 361%. It is consistent with the theory of RREA. For these phenomena of photons at LHAASO, the main factor is that the number of positrons and electrons are increased due to the acceleration of negative electric field on electrons, with increase of 65% in -1000 V/cm and 992% in -1700 V/cm, and the spectrum of positrons and electrons soften. Newborn free positrons/electrons may undergo bremsstrahlung and deposit part of their energy into photons, causing the change of number and energy of photons to follow roughly the same pattern as positrons and electrons. The simulation results can provide the information for understanding the variations of the data detected by LHAASO during thunderstorms and the acceleration mechanisms of secondary charged particles caused by an atmospheric electric field.

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 charged components of cosmic ray secondary particles are deflected by thunderstorm electric fields as they pass through kilometer-scale thunderclouds. As a result, the information on location of secondary particles reaching the observation level will be affected, and the lateral distribution will also be changed. In this paper, the Monte Carlo method is used to simulate the effects of near-earth thunderstorm electric field on the lateral distribution of secondary particles at LHAASO. A vertical and uniform atmospheric electric field model is used in our simulations. The results show that during thunderstorm the lateral distribution of secondary particles widens, and the variation amplitude is not only associated with the strength of electric field, but also dependent upon the primary energy and zenith angle of cosmic rays. In an electric field of - 1000 V·cm−1 (below the threshold of the Relatively Runaway Electron Avalanche, RREA), the variation amplitude of the lateral distribution of secondary particles is about 0.7% for θ= 0∘, and the variation amplitude reaches 4.7% for θ= 50∘. The primary energy of cosmic rays is about 180 GeV, the increasing amplitude is about 0.6%. When the primary energy is about 560 TeV, the variation can be up to 20.1%. In an electric field of - 1700 V·cm−1 (above the threshold of the RREA process), the increasing amplitude of the lateral distribution is greater than that in an electric field of - 1000 V· cm−1. And the variation amplitude is 3.8% for θ= 0∘ and 34% for θ= 50∘, respectively. For the primary energy of about 180 GeV, the increasing amplitude of secondary particles is 9.9%. For the primary energy of about 560 TeV, the variation can be as high as 119%. Our simulation results are helpful to understand the deflection mechanisms of cosmic ray secondary particles generated by the near-earth thunderstorm electric field, as well as the variation of LHAASO data during thunderstorms.

Since the science white paper of the Large High Altitude Air Shower Observatory (LHAASO) published on arXiv in 2019 [e-Print: 1905.02773 (astro-ph.HE)], LHAASO has completed the transition from a project to an operational gamma-ray astronomical observatory LHAASO is a new generation multi-component facility located in Daocheng, Sichuan province of China, at an altitude of 4410 meters. It aims at measuring with unprecedented sensitivity the spectrum, composition, and anisotropy of cosmic rays in the energy range between 10$^{12}$ and 10$^{18}$~eV, and acting simultaneously as a wide aperture (one stereoradiant) continuously operating gamma-ray telescope in the energy range between 10$^{11}$ and $10^{15}$~eV with the designed sensitivity of 1.3\% of the Crab Unit (CU) above 100 TeV. LHAASO's capability of measuring simultaneously different shower components (electrons, muons, and Cherenkov/fluorescence light), will allow it to investigate the origin, acceleration, and propagation of CR through measurement of the energy spectrum, elemental composition, and anisotropy with unprecedented resolution. The remarkable sensitivity of LHAASO will play a key role in CR physics and gamma-ray astronomy for a general and comprehensive exploration of the high energy universe and will allow important studies of fundamental physics (such as indirect dark matter search, Lorentz invariance violation, quantum gravity) and solar and heliospheric physics. The LHAASO Collaboration organized an editorial working group and finished all editorial work of this science book, to summarize the instrumental features and outline the prospects of scientific researches with the LHAASO experiment.

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.