Knowasiak

Knowledge Social Network

Sign up
ＢＵＬＬＥＴＩＮＳ

Ultrahigh-energy photons up to 1.4 PeV from 12 γ-ray Galactic sources (2021)

Vanic
Vanic
14 February 2022
114
Ultrahigh-energy photons up to 1.4 PeV from 12 γ-ray Galactic sources (2021)

Abstract

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,4,5,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.

This is a preview of subscription content

Access options

Subscribe to Journal

Get full journal access for 1 year

199,00 €

only 3,90 € per issue

Tax calculation will be finalised during checkout.

Buy article

Get time limited or full article access on ReadCube.

$32.00

All prices are NET prices.

Fig. 1: Spectral energy distributions and significance maps.

Data availability

The data supporting the conclusions of this paper are available through the LHAASO web page (http://english.ihep.cas.cn/lhaaso/index.html) in the section ‘Public Data’. All data are in ASCII, and the code used to produce the figures is publicly accessible and listed in the Public Data section.

References

  1. Aloisio, R., Coccia, E. & Vissani, F. (eds) Multiple Messengers and Challenges in Astroparticle Physics (Springer, 2018).

  2. HESS Collaboration. Acceleration of petaelectronvolt protons in the Galactic Centre. Nature 531, 476–479 (2016).

    ADS 

    Google Scholar 

  3. Amenomori, M. et al. First detection of photons with energy beyond 100 TeV from an astrophysical source. Phys. Rev. Lett. 123, 051101 (2019).

    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 

  4. Abeysekara, A. U. et al. Multiple Galactic sources with emission above 56 TeV detected by HAWC. Phys. Rev. Lett. 124, 021102 (2020).

    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 

  5. Abeysekara, A. U. et al. HAWC observations of the acceleration of very-high-energy cosmic rays in the Cygnus Cocoon. Nat. Astron. https://doi.org/10.1038/s41550-021-01318-y (2021).

  6. Amenomori, M. et al. Potential PeVatron supernova remnant G106.3+2.7 seen in the highest-energy gamma rays. Nat. Astron. https://doi.org/10.1038/s41550-020-01294-9 (2021).

  7. Cao, Z. A future project at Tibet: the large high altitude air shower observatory (LHAASO). Chin. Phys C 34, 249–252 (2010).

    ADS 
    CAS 

    Google Scholar 

  8. Aharonian, F. et al. Observation of the Crab Nebula with LHAASO-KM2A – a performance study. Chin. Phys C 45, 025002 (2021).

    ADS 

    Google Scholar 

  9. Aharonian, F. A. Very High Energy Cosmic Gamma Radiation: A Crucial Window on the Extreme Universe (World Scientific, 2004).

  10. Liu, R.-Y. & Yan, H. On the unusually large spatial extent of the TeV nebula HESS J1825–137: implication from the energy-dependent morphology. Mon. Not. R. Astron. Soc. 494, 2618–2627 (2020).

    ADS 
    CAS 

    Google Scholar 

  11. Giacinti, G., et al. Halo fraction in TeV-bright pulsar wind nebulae. Astron. Astrophys. 636, A113 (2020).

    CAS 

    Google Scholar 

  12. Abeysekara, A. U. et al. Extended gamma-ray sources around pulsars constrain the origin of the positron flux at Earth. Science 358, 911–914 (2017).

    ADS 
    CAS 
    PubMed 
    PubMed Central 

    Google Scholar 

  13. Linden, T. et al. Using HAWC to discover invisible pulsars. Phys. Rev. D 96, 103016 (2017).

    ADS 

    Google Scholar 

  14. Sudoh, T., Linden, T. & Beacom, J. F. TeV halos are everywhere: prospects for new discoveries. Phys. Rev. D 100, 043016 (2019).

    ADS 
    CAS 

    Google Scholar 

  15. Gabici, S. & Aharonian, F. A. Searching for Galactic cosmic-ray pevatrons with multi-TeV gamma rays and neutrinos. Astrophys. J. Lett. 665, L131–L134 (2007).

    ADS 
    CAS 

    Google Scholar 

  16. Aharonian, F. A. & Atoyan, A. M. On the emissivity of π0-decay gamma radiation in the vicinity of accelerators of Galactic cosmic rays. Astron. Astrophys. 309, 917–928 (1996).

    ADS 
    CAS 

    Google Scholar 

  17. Cesarsky, C. J. & Montmerle, T. Gamma-rays from active regions in the Galaxy – the possible contribution of stellar winds. Space Sci. Rev. 36, 173–193 (1983).

    ADS 

    Google Scholar 

  18. Bykov, A. M. et al. High-energy particles and radiation in star-forming regions. Space Sci. Rev. 216, 42 (2020).

    ADS 
    CAS 

    Google Scholar 

  19. Aharonian, F., Yang, R. & de Oña Wilhelmi, E. Massive stars as major factories of Galactic cosmic rays. Nat. Astron. 3, 561–567 (2019).

    ADS 

    Google Scholar 

  20. Akaike, H. A new look at the statistical model identification. IEEE Trans. Automat. Contr. 19, 716C723 (1974)

    MathSciNet 

    Google Scholar 

  21. He, H. Design of the LHAASO detectors. Rad. Det. Technol. Meth. 2, 7 (2018).

    Google Scholar 

  22. Fleysher, R., Fleysher, L., Nemethy, P. & Mincer, A. I. Tests of statistical significance and background estimation in gamma-ray air shower experiments. Astrophys. J. 603, 355–362 (2004).

    ADS 

    Google Scholar 

  23. Bartoli, B. et al. TeV gamma-ray survey of the northern sky using the ARGO-YBJ detector. Astrophys. J. 779, 27 (2013).

    ADS 

    Google Scholar 

  24. Cronin, J. W., Gibbs, K. G. & Weekes, T. C. The search for discrete astrophysical sources of energetic gamma radiation. Annu. Rev. Nucl. Part. Sci. 43, 883–925 (1993); erratum 43, 883–937 (1993).

    ADS 
    CAS 

    Google Scholar 

  25. Aharonian, F., Buckley, J., Kifune, T. & Sinnis, G. High energy astrophysics with ground-based gamma ray detectors. Rep. Prog. Phys. 71, 096901 (2008).

    ADS 

    Google Scholar 

  26. Aharonian, F. et al. The Crab nebula and pulsar between 500 GeV and 80 TeV: observations with the HEGRA stereoscopic air Cerenkov telescopes. Astrophys. J. 614, 897–913 (2004).

    ADS 
    CAS 

    Google Scholar 

  27. MAGIC Collaboration. MAGIC very large zenith angle observations of the Crab Nebula up to 100 TeV. Astron. Astrophys. 635, A158 (2020).

    Google Scholar 

  28. Gould, R. J. & Schréder, G. Opacity of the Universe to high-energy photons. Phys. Rev. Lett. 16, 252–254 (1966).

    ADS 
    CAS 

    Google Scholar 

  29. Popescu, C. C. et al. A radiation transfer model for the Milky Way: I. Radiation fields and application to high-energy astrophysics. Mon. Not. R. Astron. Soc. 470, 2539–2558 (2017).

    ADS 
    CAS 

    Google Scholar 

  30. Moskalenko, I. V., Porter, T. A. & Strong, A. W. Attenuation of very high energy gamma rays by the Milky Way interstellar radiation field. Astrophys. J. Lett. 640, L155–L158 (2006).

    ADS 
    CAS 

    Google Scholar 

  31. Abdollahi, S. et al. Fermi Large Area Telescope fourth source catalog. Astrophys. J. Suppl. Ser. 247, 33 (2020).

    ADS 
    CAS 

    Google Scholar 

  32. Dame, T. M., Hartmann, D. & Thaddeus, P. The Milky Way in molecular clouds: a new complete CO survey. Astrophys. J. 547, 792–813 (2001).

    ADS 
    CAS 

    Google Scholar 

  33. Zirakashvili, V. N. & Aharonian, F. Analytical solutions for energy spectra of electrons accelerated by nonrelativistic shock-waves in shell type supernova remnants. Astron. Astrophys. 465, 695–702 (2007).

    ADS 
    CAS 
    MATH 

    Google Scholar 

  34. Khangulyan, D., Aharonian, F. A. & Kelner, S. R. Simple analytical approximations for treatment of inverse Compton scattering of relativistic electrons in the blackbody radiation field. Astrophys. J. 783, 100 (2014).

    ADS 

    Google Scholar 

  35. Kafexhiu, E., Aharonian, F., Taylor, A. M. & Vila, G. S. Parametrization of gamma-ray production cross sections for p p interactions in a broad proton energy range from the kinematic threshold to PeV energies. Phys. Rev. D 90, 123014 (2014).

    ADS 

    Google Scholar 

  36. Manchester, R. N., Hobbs, G. B., Teoh, A. & Hobbs, M. The Australia Telescope National Facility pulsar catalogue. Astron. J. 129, 1993–2006 (2005).

    ADS 

    Google Scholar 

  37. Shan, S. S. et al. Distances of Galactic supernova remnants using red clump stars. Astrophys. J.S 238, 35 (2018).

    ADS 

    Google Scholar 

  38. Acero, F. et al. Constraints on the Galactic population of TeV pulsar wind nebulae using Fermi Large Area Telescope observations. Astrophys. J. 773, 77 (2013).

    ADS 

    Google Scholar 

  39. Ranasinghe, S. & Leahy, D. A. Distances to supernova remnants G20.4 + 0.1, G24.7 − 0.6, and G28.6 − 0.1 and new molecular cloud associations. Mon. Not. R. Astron. Soc. 477, 2243–2250 (2018).

    ADS 
    CAS 

    Google Scholar 

  40. Gotthelf, E. V., Halpern, J. P., Terrier, R. & Mattana, F. Discovery of an energetic 38.5 ms pulsar powering the gamma-ray source IGR J18490-0000/HESS J1849-000. Astrophys. J. Lett. 729, L16 (2011).

    ADS 

    Google Scholar 

  41. Zhang, B. et al. The parallax of W43: a massive star-forming complex near the Galactic bar. Astrophys. J. 781, 89 (2014).

    ADS 

    Google Scholar 

  42. Yang, J., Zhang, J.-L., Cai, Z.-Y., Lu, D.-R. & Tan, Y.-H. Molecular gas distribution around the supernova remnant G40.5 0.5. Chin. J. Astron. Astrophys. 6, 210–216 (2006).

    ADS 
    CAS 

    Google Scholar 

  43. Downes, A. J. B., Pauls, T. & Salter, C. J. G 40.5-0.5: a previously unrecognised supernova remnant in Aquila. Astron. Astrophys. 92, 47–50 (1980).

    ADS 

    Google Scholar 

  44. Gelfand, J. D., Slane, P. O. & Temim, T. The properties of the progenitor supernova, pulsar wind, and neutron star inside PWN G54.1+0.3. Astrophys. J. 807, 30 (2015).

    ADS 

    Google Scholar 

  45. Kirichenko, A. et al. Optical observations of PSR J2021+3651 in the Dragonfly nebula with the GTC. Astrophys. J. 802, 17 (2015).

    ADS 

    Google Scholar 

  46. Azimlu, M. & Fich, M. Study of molecular clouds associated with H II regions. Astron. J. 141, 123 (2011).

    ADS 

    Google Scholar 

  47. Paredes, J. M. et al. Radio continuum and near-infrared study of the MGRO J2019+37 region. Astron. Astrophys. 507, 241–250 (2009).

    ADS 

    Google Scholar 

  48. Bartoli, B. et al. Observation of the TeV gamma-ray source MGRO J1908+06 with ARGO-YBJ. Astrophys. J. 760, 110 (2012).

    ADS 

    Google Scholar 

  49. Chaves, R. C. G. et al. The H.E.S.S. Galactic plane survey. Astron. Astrophys. 612, A1 (2018).

    Google Scholar 

  50. Rygl, K. L. J. et al. Parallaxes and proper motions of interstellar masers toward the Cygnus X star-forming complex. I. Membership of the Cygnus X region. Astron. Astrophys. 539, A79 (2012).

    Google Scholar 

  51. Kothes, R., Uyaniker, B. & Pineault, S. The supernova remnant G106.3+2.7 and its pulsar-wind nebula: relics of triggered star formation in a complex environment. Astrophys. J. 560, 236–243 (2001).

    ADS 

    Google Scholar 

Download references

Acknowledgements

This work is supported in China by National Key R&D programme of China under grants 2018YFA0404201, 2018YFA0404202, 2018YFA0404203 and 2018YFA0404204 and by NSFC (grant numbers 12022502, 11905227, U1931112, 11635011, 11761141001 and U2031105), and in Thailand by RTA6280002 from Thailand Science Research and Innovation. We thank all staff members who work year-round at the LHAASO site at 4,400 m above sea level to maintain the detector and keep the electricity power supply and other components of the experiment operating smoothly. We are grateful to Chengdu Management Committee of Tianfu New Area for constant financial support to research with LHAASO data.

Author information

Affiliations

  1. Key Laboratory of Particle Astrophysics & Experimental Physics Division & Computing Center, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing, China

    Zhen Cao, Y. X. Bai, X. J. Bi, Y. J. Bi, J. F. Chang, X. C. Chang, L. Chen, M. J. Chen, M. L. Chen, S. H. Chen, S. Z. Chen, X. L. Chen, N. Cheng, Y. D. Cheng, H. L. Dai, X. J. Dong, Z. X. Fan, K. Fang, S. H. Feng, B. Gao, L. S. Geng, Q. B. Gou, M. H. Gu, J. G. Guo, Y. Q. Guo, Y. Y. Guo, H. H. He, J. C. He, C. Hou, H. B. Hu, S. C. Hu, Q. L. Huang, F. Ji, X. L. Ji, C. Jin, Cong Li, F. Li, H. B. Li, H. C. Li, J. Li, K. Li, X. R. Li, Y. Z. Li, Zhe Li, C. Liu, J. Liu, J. Y. Liu, W. Liu, H. K. Lv, L. L. Ma, X. H. Ma, M. Y. Qi, X. D. Sheng, J. R. Shi, B. D. Wang, L. Y. Wang, X. J. Wang, Y. D. Wang, Y. J. Wang, Y. P. Wang, Zheng Wang, Y. J. Wei, C. Y. Wu, H. R. Wu, S. Wu, G. Xiao, F. F. Yang, M. J. Yang, Z. G. Yao, L. Q. Yin, X. H. You, Z. Y. You, T. X. Zeng, Z. K. Zeng, M. Zha, X. X. Zhai, S. S. Zhang, X. P. Zhang, Yong Zhang, Yi Zhang, Y. L. Zhang, J. Zhao, B. Zhou, K. J. Zhu & X. Zuo

  2. TIANFU Cosmic Ray Research Center, Chengdu, Sichuan, China

    Zhen Cao, Y. X. Bai, X. J. Bi, Y. J. Bi, J. F. Chang, X. C. Chang, L. Chen, M. J. Chen, M. L. Chen, S. H. Chen, S. Z. Chen, X. L. Chen, N. Cheng, Y. D. Cheng, H. L. Dai, X. J. Dong, Z. X. Fan, K. Fang, S. H. Feng, B. Gao, L. S. Geng, Q. B. Gou, M. H. Gu, J. G. Guo, Y. Q. Guo, Y. Y. Guo, H. H. He, J. C. He, C. Hou, H. B. Hu, S. C. Hu, Q. L. Huang, F. Ji, X. L. Ji, C. Jin, Cong Li, F. Li, H. B. Li, H. C. Li, J. Li, K. Li, X. R. Li, Y. Z. Li, Zhe Li, C. Liu, J. Liu, J. Y. Liu, W. Liu, H. K. Lv, L. L. Ma, X. H. Ma, M. Y. Qi, X. D. Sheng, J. R. Shi, B. D. Wang, L. Y. Wang, X. J. Wang, Y. D. Wang, Y. J. Wang, Y. P. Wang, Zheng Wang, Y. J. Wei, C. Y. Wu, H. R. Wu, S. Wu, G. Xiao, F. F. Yang, M. J. Yang, Z. G. Yao, L. Q. Yin, X. H. You, Z. Y. You, T. X. Zeng, Z. K. Zeng, M. Zha, X. X. Zhai, S. S. Zhang, X. P. Zhang, Yong Zhang, Y. L. Zhang, J. Zhao, B. Zhou, K. J. Zhu & X. Zuo

  3. University of Chinese Academy of Sciences, Beijing, China

    Zhen Cao, X. J. Bi, L. Chen, S. H. Chen, X. L. Chen, J. G. Guo, Y. Y. Guo, H. H. He, J. C. He, H. B. Hu, S. C. Hu, C. Jin, Y. Z. Li, W. W. Tian, Y. P. Wang, Y. J. Wei, G. M. Xiang, Z. Y. You, Z. K. Zeng & K. J. Zhu

  4. Dublin Institute for Advanced Studies, Dublin, Ireland

    F. A. Aharonian

  5. Max-Planck-Institut for Nuclear Physics, Heidelberg, Germany

    F. A. Aharonian

  6. State Key Laboratory of Particle Detection and Electronics, Beijing, China

    Q. An, Zhe Cao, J. F. Chang, M. L. Chen, M. H. Gu, X. L. Ji, K. Jiang, Cheng Li, F. Li, J. Li, X. Li, Z. B. Tang, Zheng Wang, F. F. Yang, T. X. Zeng, L. Zhao & K. J. Zhu

  7. University of Science and Technology of China, Hefei, China

    Q. An, Zhe Cao, K. Jiang, Cheng Li, H. Y. Li, X. Li, B. Liu, Z. B. Tang, J. Xia, R. Z. Yang, R. Zhang & L. Zhao

  8. School of Physical Science and Technology & School of Information Science and Technology, Southwest Jiaotong University, Chengdu, China

    Axikegu, Long Chen, Q. H. Chen, X. L. Guo, Y. He, D. H. Huang, Z. C. Huang, H. Y. Jia, Xin Li, H. Liu, W. J. Long, A. Masood, B. Y. Pang, Q. N. Sun, H. Wang, R. N. Wang, W. X. Wu, S. Q. Xi, J. J. Xia, Y. L. Xin, Y. F. Zhang, B. Zhao, Y. Zheng, X. X. Zhou & F. R. Zhu

  9. College of Physics, Sichuan University, Chengdu, China

    L. X. Bai, J. Chen, S. Hu, Y. Li, Z. X. Liu, Z. H. Wang, C. W. Yang, Y. H. Yao, J. W. Zhang & R. Zhou

  10. School of Astronomy and Space Science, Nanjing University, Nanjing, China

    Y. W. Bao, Y. Chen, Z. G. Dai, R. Y. Liu, X. Y. Wang, B. B. Zhang, H. M. Zhang, X. Zhang & P. Zhou

  11. Center for Astrophysics, Guangzhou University, Guangzhou, China

    D. Bastieri, J. T. Cai, J. H. Fan, S. L. He, Z. Y. Pei, H. G. Wang, H. B. Xiao & L. X. Zhang

  12. School of Physics and Technology, Wuhan University, Wuhan, China

    H. Cai, W. Wang & G. G. Xin

  13. Key Laboratory of Dark Matter and Space Astronomy, Purple Mountain Observatory, Chinese Academy of Sciences, Nanjing, China

    J. Chang, H. L. Dai, Y. Z. Fan, L. Feng, Y. L. Feng, Y. Y. Guo, H. N. He, H. Y. Li, S. M. Liu, D. M. Wei, J. J. Wei, X. F. Wu, J. Xia, Q. Yuan, H. D. Zeng, R. Zhang, Yi Zhang & S. P. Zhao

  14. Hebei Normal University, Shijiazhuang, China

    B. M. Chen, S. W. Cui, B. B. Li, L. Shao, L. Zhang, P. P. Zhang, S. R. Zhang & L. Z. Zhao

  15. Key Laboratory for Research in Galaxies and Cosmology, Shanghai Astronomical Observatory, Chinese Academy of Sciences, Shanghai, China

    Liang Chen, G. M. Xiang, Y. Xing & J. N. Zhou

  16. Key Laboratory of Cosmic Rays (Tibet University), Ministry of Education, Lhasa, Tibet, China

    T. L. Chen,  Danzengluobu, Q. Gao & M. Y. Liu

  17. National Astronomical Observatories, Chinese Academy of Sciences, Beijing, China

    X. H. Cui, W. W. Tian, J. L. Zhang & H. Zhu

  18. School of Physics and Astronomy & School of Physics (Guangzhou), Sun Yat-sen University, Zhuhai, China

    Y. D. Cui, X. B. He, Y. K. Hor, S. J. Lin, J. S. Liu, P. H. T. Tam, W. Wang, J. Y. Yang & L. L. Yang

  19. School of Physics and Astronomy, Yunnan University, Kunming, China

    B. Z. Dai, J. Fang, M. M. Ge, Z. J. Jiang, R. Lu, Z. X. Wang, T. Wen, S. B. Yang, W. Zeng, Li Zhang & P. F. Zhang

  20. Département de Physique Nucléaire et Corpusculaire, Faculté de Sciences, Université de Genève, Geneva, Switzerland

    D. della Volpe, M. Heller & T. Montaruli

  21. Dipartimento di Fisica dell’Università di Napoli “Federico II”, Complesso Universitario di Monte Sant’Angelo, Naples, Italy

    B. D′Ettorre Piazzoli

  22. Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, China

    C. F. Feng, C. D. Gao, W. Gao, W. H. Huang, X. T. Huang, W. L. Li, D. Liu, Y. C. Nan, L. P. Wang, L. Xue, N. Yin, Y. H. Yu, H. Y. Zhang, S. P. Zhao & C. G. Zhu

  23. Department of Engineering Physics, Tsinghua University, Beijing, China

    G. H. Gong, X. J. Hu, Y. N. Liu & Y. M. Ye

  24. School of Physics and Microelectronics, Zhengzhou University, Zhengzhou, China

    Y. A. Han & H. D. Liu

  25. Yunnan Observatories, Chinese Academy of Sciences, Kunming, China

    X. Hou, J. R. Mao, J. C. Wang & D. H. Yan

  26. Institute for Nuclear Research of Russian Academy of Sciences, Moscow, Russia

    D. Kuleshov, K. Levochkin, V. Rulev, O. Shchegolev, Yu. V. Stenkin & V. Stepanov

  27. School of Physics, Peking University, Beijing, China

    Zhuo Li, B. Q. Ma, H. C. Song & R. X. Xu

  28. School of Physical Science and Technology, Guangxi University, Nanning, China

    E. W. Liang, Y. F. Liang, X. N. Sun & X. G. Wang

  29. Tsung-Dao Lee Institute, Shanghai Jiao Tong University, Shanghai, China

    J. L. Liu, J. S. Wang, Zhen Wang, D. L. Xu & H. Zhou

  30. School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai, China

    J. L. Liu, J. S. Wang, Zhen Wang, D. L. Xu & H. Zhou

  31. Department of Physics, Faculty of Science, Mahidol University, Bangkok, Thailand

    W. Mitthumsiri, P. Pattarakijwanich, D. Ruffolo & A. Sáiz

  32. Moscow Institute of Physics and Technology, Moscow, Russia

    O. Shchegolev & Yu. V. Stenkin

  33. National Space Science Center, Chinese Academy of Sciences, Beijing, China

    Z. B. Sun, C. Wang & F. Zheng

Contributions

This work is the result of the contributions and efforts of all members and institutes of the LHAASO Collaboration under the leadership of Zhen Cao, who is the spokesperson of the LHAASO collaboration. Most of the listed 32 institutes participate in the construction of the LHAASO detectors. In particular, H.H.H. leads the design and construction of KM2A, whose scintillator detector array is constructed by teams led by X.D.S. and muon detector array is constructed by teams led by G.X. KM2A is operated by teams led by J.L. and X. Zuo. The calibration of LHAASO detectors is led by H.H.H. with the participation of groups from IHEP, led by H.K.L. and X. Zuo; from Shandong University, led by C.F.F.; from the University of Science and Technology of China, led by C. Li; and from Southwest Jiaotong University, led by H.Y.J. The CR event reconstruction is led by S.Z.C. with the participation of groups from IHEP, Shandong University and Southwest Jiaotong University. The data analysis by groups from IHEP and Shandong University is led by S.Z.C.; from Sun Yat-sen University by P.H.T.T.; and from the University of Science and Technology of China by R.Z.Y. The team at Purple Mt Observatory, led by Yi Zhang, performs an independent cross-check of the analysis. R.Z.Y. and R.Y.L. lead the interpretation of the results, with participation of the groups from Nanjing University, led by Y.C. and X.Y.W.; from Yunan University, led by L. Zhang; from Purple Mt Observatory, led by Q.Y. Zhen Cao, F.A.A. and B.D.P. supervised the data analysis and interpretation and also lead the paper writing. S.M.L. and D.d.V. contribute as co-chairs of the editorial board of LHAASO. Many groups participated in the analysis procedures, including the groups from Mahidol University, led by D.R.; from Shanghai Jiaotong University, led by H. Zhou; from Shanghai Observatory, led by Z.X.W.; from Peking University, led by Zhuo Li; from Guangxi University, led by E.W.L.; from National Observatory, led by W.W.T.; from Wuhan University, led by P. H. T. Tam; and from Yunnan Observatory, led by J.C.W. All other groups—from the University of Chinese Academy of Science and Tianfu Cosmic ray Research Center, led by Zhen Cao; from State Key Laboratory of Particle Detection and Electronics, led by Q.A.; from Tsinghua University, led by Y.N.L.; from the National Space Science Center, led by Z.B.S.; from Guangzhou University, led by J.H.F.; from Hebei Normal University, led by S.W.C.; from Zhengzhou University, led by H.D.L.; from Sichuan University, led by C.W.Y.; from Tibet University, led by T.L.C.; and from the Institute for Nuclear Research and Moscow Institute of Physics and Technology, led by Yu.V.S.—contributed to detector construction and to reviewing the final version before submission.

Corresponding authors

Correspondence to
Zhen Cao, F. A. Aharonian, S. Z. Chen, R. Y. Liu or R. Z. Yang.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature thanks Razmik Mirzoyan and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 Schematic drawing of the LHAASO layout7.

Small red dots indicate the 5,195 scintillator counters, with a spacing of 15 m in the central area of 1 km2 and of 30 m in the skirt area of 0.3 km2 of KM2A. Big blue dots indicate the 1,188 muon detectors distributed in the central area with a spacing of 30 m. The three light-blue rectangles in the centre indicate WCDA, of 78,000 m2 in total. Small black rectangles near WCDA indicate 18 telescopes of WFCTA.

Extended Data Fig. 2 γ-ray energy distributions in SED bins and corresponding bin purity.

a, Distributions of thrown-in energies Etrue in the simulation of events that are in the bins, defined by reconstructed energies Erec above 10 TeV. The fraction on the vertical axis is defined as (dN/dEtrue)/Ni, where Ni is the total number of events in the ith bin, defined by the reconstructed energy Erec. A bin width of Δ(logErec) = 0.2 is selected, according to the energy resolution of 14% above 100 TeV (ref. 8). A power-law spectrum proportional to E−3.09 is assumed here. b, Bin purity, defined as the fraction of events with Etrue in the bin, as a function of Erec in the range from 10 TeV to 2.5 PeV.

Extended Data Fig. 3 Distributions (dots) of events in the reconstructed energy bins from 10 TeV to 10 PeV.

The input (solid grey line) events are generated using the power-law SED determined from the measured SED of the Crab Nebula8. For each case, the input SED has an artificial cutoff (dashed lines) at the Ecut values listed in the key. The distributions demonstrate a clear spillover effect by the events in the bins above Ecut. The effect becomes weaker at higher energies. There is no indication of pollution of the bins above 1 PeV by events at input energies lower than 0.3 PeV.

Extended Data Fig. 4 LHAASO sky map at energies above 100 TeV.

The circles indicate the positions of known very-high-energy γ-ray sources.

Extended Data Fig. 5 Phenomenological fits to the γ−ray observations of LHAASO J1908+0621, and previous observations of potential counterparts.

The inset shows the KM2A significance map, indicating the potential counterparts of the UHE γ-ray source. The colour bar shows the significance ((sqrt{{rm{TS}}})). The green circle indicates the PSF of LHAASO. The Fermi LAT points for LHAASO J1908+0621 analysed in this work, as well as ARGO48, HESS49 and HAWC4 data, are shown together with the LHAASO measurements. The dotted curve shows the leptonic model of radiation, assuming an injection of electron/positron pairs according to the pulsar’s spin-down behaviour, with a breaking index of 2 and an initial rotation period of 0.04 s. A fraction of 6% of the current spin-down power of the pulsar PSR J1907+0602 at a distance of 2.4 kpc is assumed to be converted to e± pairs to support the γ-ray emission. The injection spectrum of electrons is assumed to be (N(E)propto {E}_{e}^{-1.75}exp {-{[{E}_{e}/(800{rm{TeV}})]}^{2}}). The solid curves correspond to the hadronic model of radiation. Two types of energy distributions are assumed for the parent proton population: (i) a single power-law spectrum of parent protons, N(E) ≈ E−1.85exp[−E/(380 TeV)] (thin solid curve); (ii) a broken power-law spectrum with an exponential cutoff of parent protons, with indices 1.2 and 2.7 below and above 25 TeV, respectively, and a cutoff energy of 1.3 PeV (thick solid curve). In the inset sky map, the black diamond shows the position of PSR J1907+0602, the black contours correspond to the location of supernova remnant SNR G40.5-0.5 and the white circle is the position and size of HESS J1908+063. The cyan regions are the dense clumps described in Methods. The average density in the whole γ-ray emission region is estimated to be about 10 cm−3. γ-ray absorption due to photon–photon pair production (see Methods) is taken into account in the theoretical curve.

Extended Data Table 1 Number of on-source events of energy >100 TeV, residual CR background events and corresponding exposure time for the 12 UHE sources
Extended Data Table 2 List of energetic astrophysical objects possibly associated with each LHAASO source

Supplementary information

About this article

Verify currency and authenticity via CrossMark

Cite this article

Cao, Z., Aharonian, F.A., An, Q. et al. Ultrahigh-energy photons up to 1.4 petaelectronvolts from 12 γ-ray Galactic sources.
Nature 594, 33–36 (2021). https://doi.org/10.1038/s41586-021-03498-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41586-021-03498-z

Further reading

  • On the surface brightness radial profile of the extended γ-ray sources

    • Rui-Zhi Yang
    • Bing Liu

    Science China Physics, Mechanics & Astronomy (2022)

  • Regimes of cosmic-ray diffusion in Galactic turbulence

    • P. Reichherzer
    • L. Merten
    • E. G. Zweibel

    SN Applied Sciences (2022)

  • An ultra-high-energy γ-ray telescope at 4,410 m

    • Zhen Cao

    Nature Astronomy (2021)

  • Hunting the strongest accelerators in our Galaxy

    • Petra Huentemeyer

    Nature (2021)

  • LHAASO: A PeVatrons explorer

    • Felix A. Aharonian

    Science China Physics, Mechanics & Astronomy (2021)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Read More

Vanic
WRITTEN BY

Vanic

“Simplicity, patience, compassion.
These three are your greatest treasures.
Simple in actions and thoughts, you return to the source of being.
Patient with both friends and enemies,
you accord with the way things are.
Compassionate toward yourself,
you reconcile all beings in the world.”
― Lao Tzu, Tao Te Ching