宇宙线光谱的扩展超过了1千兆电子伏（PeV；10-15千兆电子伏）表明了所谓的皮瓦特罗宇宙线工厂的存在，这些工厂将粒子加速到PeV能量。我们需要定位和识别这些物体，以找到银河宇宙线的起源。电子和质子质子质子的主要特征是超高能（超过100 TeV）γ辐射。根据对延伸至0.04 PeV（参考文献2）的硬光谱辐射的探测，已经在银河系中心发现了质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子质子。虽然有报道称，在银河系平面3、4、5、6中的一些物体发出了能量略高于0.1PEV的γ射线，但无偏识别和深入探索粒子加速器需要探测能量远高于0.1PEV的γ射线。在这里，我们报告了从12个超高能γ射线源中检测到的超过530个光子，其能量超过100太安电子伏，高达1.4 PeV，其统计显著性大于7个标准偏差。尽管附近有几个潜在的对应物，包括脉冲星风星云、超新星遗迹和恒星形成区，但负责超高能γ射线的粒子尚未被确定和确定（蟹状星云除外），从而为这些极端加速器的起源留下了空白。
我们有一个专门为中国客户服务的网站。请去中国。com订阅这本杂志。
支持本文结论的数据可通过拉萨索网站（http://english.ihep.cas.cn/lhaaso/index.html)在“公共数据”部分。所有数据都是ASCII码，用于生成数字的代码可以公开访问，并在公共数据部分列出。
Aloisio，R.，Coccia，E.&；Vissani，F.（编辑）天体粒子物理学中的多重信使和挑战（Springer，2018）。
赫斯合作。银河系中心的质子加速。《自然》杂志531476-479（2016）。
Amenomori，M.等人首次从天体物理源探测到能量超过100 TeV的光子。菲斯。牧师。莱特。123, 051101 (2019).
Abeysekara，A.U.等人。HAWC检测到的辐射超过56 TeV的多个星系源。菲斯。牧师。莱特。124, 021102 (2020).
Abeysekara，A.U.等人HAWC对天鹅座茧中超高能宇宙线加速的观测。纳特。阿斯顿。https://doi.org/10.1038/s41550-021-01318-y (2021).
阿蒙诺莫里，M.等。潜在的粒子超新星遗迹G106。在最高能量的伽马射线中可以看到3+2.7。纳特。阿斯顿。https://doi.org/10.1038/s41550-020-01294-9 (2021).
曹，Z.西藏未来的一个项目：大型高空空气簇射观测站（拉索）。颏物理C 34249–252（2010）。
Aharonian，F.等人。用拉索-KM2A观测蟹状星云——一项性能研究。颏物理C 45025002（2021年）。
《超高能宇宙伽马辐射：了解极端宇宙的关键窗口》（世界科学出版社，2004年）。
刘若瑜；Yan，H.关于TeV星云HESS J1825–137异常大的空间范围：能量相关形态学的含义。周一。不R.阿斯顿。Soc。494, 2618–2627 (2020).
Giacinti，G.，等。TeV明亮脉冲星风星云中的晕分数。阿斯顿。天体物理学。636，A113（2020年）。
Abeysekara，A.U.等人。脉冲星周围的扩展伽马射线源限制了地球正电子通量的起源。《科学》杂志3581911–914（2017）。
Linden，T.等人利用HAWC发现了看不见的脉冲星。菲斯。牧师。D 96103016（2017年）。
T.苏杜、T.林登和；Beacom，J.F.TeV光环无处不在：新发现的前景。菲斯。牧师。D 100043016（2019年）。
加比奇，S.&；Aharonian，F.A.正在寻找具有多TeV伽马射线和中微子的银河宇宙线粒子。天体物理学。J.莱特。665，L131–L134（2007年）。
阿哈罗尼安，F.A.&；Atoyan，A.M.关于银河宇宙线加速器附近π0衰变伽马辐射的发射率。阿斯顿。天体物理学。309, 917–928 (1996).
塞萨斯基，C.J.&；蒙特默尔，T.来自银河系活跃区域的伽马射线——恒星风的可能贡献。太空科学。牧师。36, 173–193 (1983).
Bykov，A.M.等人，《恒星形成区的高能粒子和辐射》。太空科学。牧师。216, 42 (2020).
阿哈罗尼安·F·杨·R&；德欧尼亚·威廉米（de Oña Wilhelmi，E.）将大质量恒星视为银河宇宙线的主要工厂。纳特。阿斯顿。3, 561–567 (2019).
Akaike，H.统计模型识别的新视角。IEEE Trans。自动装置。控制。19716C723（1974）
弗雷舍，R.，弗雷舍，L.，内梅西，P&；Mincer，A.I.伽马射线空气簇射实验中的统计显著性检验和背景估计。天体物理学。J.603355–362（2004年）。
Bartoli，B.等人。使用ARGO-YBJ探测器对北部天空进行TeV伽马射线测量。天体物理学。J.77927（2013年）。
克罗宁，J.W.，吉布斯，K.G.&；Weekes，T.C.《寻找高能伽马辐射的离散天体物理源》。阿努。牧师。Nucl。部分Sci。43, 883–925 (1993); 勘误表43883-937（1993）。
F.阿哈罗尼安、J.巴克利、T.基富内、；Sinnis，G.地面伽马射线探测器高能天体物理学。众议员Prog。菲斯。71, 096901 (2008).
Aharonian，F.等人。500 GeV和80 TeV之间的蟹状星云和脉冲星：使用HEGRA立体空气切伦科夫望远镜进行观测。天体物理学。J.614897–913（2004年）。
神奇的合作。神奇的大天顶角观测蟹状星云高达100 TeV。阿斯顿。天体物理学。635，A158（2020年）。
古尔德，R.J.&；Schréder，G.宇宙对高能光子的不透明度。菲斯。牧师。莱特。16, 252–254 (1966).
Popescu，C.C.等人。银河系的辐射传输模型：I.辐射场及其在高能天体物理学中的应用。周一。不R.阿斯顿。Soc。470, 2539–2558 (2017).
Moskalenko，I.V.，Porter，T.A.&；银河系星际辐射场对高能伽马射线的强烈衰减。天体物理学。J.莱特。640，L155–L158（2006年）。
Abdollahi，S.等人，《费米大面积望远镜第四源目录》。天体物理学。J.补充。爵士。247, 33 (2020).
达姆，T.M.，哈特曼，D.&；《分子云中的银河：一个新的完整的联合调查》。天体物理学。J.547792–813（2001年）。
Zirakashvili，V.N.&；Aharonian，F.壳型超新星遗迹中非相对论激波加速电子能谱的解析解。阿斯顿。天体物理学。465, 695–702 (2007).
坎格利安，D.，阿哈罗尼安，F.A.&；Kelner，S.R.黑体辐射场中相对论电子逆康普顿散射处理的简单解析近似。天体物理学。J.783100（2014年）。
卡菲休，E.，阿哈罗尼安，F.，泰勒，A.M.&；Vila，G.S.在从运动学阈值到PeV能量的广泛质子能量范围内，p-p相互作用的γ射线产生截面的参数化。菲斯。牧师。D 90123014（2014年）。
R.N.曼彻斯特，G.B.霍布斯，Teoh，A.和；霍布斯，M.澳大利亚望远镜国家设施脉冲星目录。阿斯顿。J.129，1993年至2006年（2005年）。
Shan，S.S.等人。使用红色团块恒星测量银河系超新星遗迹的距离。天体物理学。J.S 238，35（2018年）。
Acero，F.等人利用费米大面积望远镜观测对TeV脉冲星风星云的银河系人口的限制。天体物理学。J.773，77（2013年）。
南部拉纳星河；莱希，华盛顿特区。超新星遗迹的距离。4+0.1，G24。7.− 0.6和G28。6.− 0.1和新的分子云关联。周一。不R.阿斯顿。Soc。477, 2243–2250 (2018).
Gottelf，E.V.，Halpern，J.P.，Terrier，R.&；马塔纳，F.发现一颗高能38.5毫秒脉冲星，为伽马射线源IGR J18490-0000/HESS J1849-000提供动力。天体物理学。J.莱特。729，L16（2011年）。
Zhang，B.等人，《W43的视差：银河棒附近的一个大质量恒星形成复合体》。天体物理学。J.78189（2014年）。
杨，J.，张，J-L.，蔡，Z-Y.，陆，D-R.&；Tan，Y.-H.超新星遗迹G40周围的分子气体分布。5 0.5. 颏J.阿斯顿。天体物理学。6, 210–216 (2006).
唐斯，A.J.B.，保罗，T&；索尔特，C.J.G 40.5-0.5：阿奎拉中一颗先前未被确认的超新星遗迹。阿斯顿。天体物理学。92, 47–50 (1980).
Gelfand，J.D.，Slane，P.O.&；T.特米姆。PWN G54内部的前体超新星、脉冲星风和中子星的性质。1+0.3. 天体物理学。J.807，30（2015年）。
Kirichenko，A.等。用GTC对蜻蜓星云中PSR J2021+3651的光学观测。天体物理学。J.802，17（2015年）。
阿齐姆鲁，M.和；Fich，M.与H II区域相关的分子云研究。阿斯顿。J.141123（2011年）。
Paredes，J.M.等人。MGRO J2019+37区域的射电连续统和近红外研究。阿斯顿。天体物理学。507, 241–250 (2009).
Bartoli，B.等人。用ARGO-YBJ观测TeV伽马射线源MGRO J1908+06。天体物理学。J.760110（2012年）。
Chaves，R.C.G.等人，H.E.S.S.银河平面勘测。阿斯顿。天体物理学。612，A1（2018年）。
Rygl，K.L.J.等人。星际脉泽朝向天鹅座X恒星形成复合体的视差和自行。I.天鹅座X区域的成员资格。阿斯顿。天体物理学。539，A79（2012年）。
科蒂斯，R.，乌亚尼克，B.&；超新星遗迹G106。3+2.7及其脉冲星风星云：复杂环境中触发恒星形成的遗迹。天体物理学。J.560236–243（2001年）。
这项工作在中国得到了国家重点研发中心的支持；D中国的项目，由国家自然科学基金委员会和泰国科学研究与创新基金委员会分别授予2018YFA040201、2018YFA040202、2018YFA040203和2018YFA040204（授予号：12022502、11905227、U1931112、11635011、1176141001和U2031105），以及泰国的RTA6280002。我们感谢常年在海拔4400米的拉索现场工作的所有工作人员，他们维护了探测器，并保持了实验的电力供应和其他部件的顺利运行。我们感谢天府新区成都管委会对拉萨索数据研究的持续资助。
作者
这项工作是拉萨索合作组织所有成员和机构在拉萨索合作组织发言人曹震领导下做出贡献和努力的结果。列出的32个研究所中的大多数都参与了拉索探测器的建设。特别是，H.H.H.领导KM2A的设计和建造，其闪烁体探测器阵列由X.D.S.领导的团队建造，μ子探测器阵列由G.X.领导的团队建造。KM2A由J.L.和X.左领导的团队运营。拉萨索探测器的校准由H.H.H.牵头，IHEP的团队参与，由H.K.L.和X.左牵头；来自山东大学，由C.F.F.领导。；来自中国科技大学，由C. Li领导；来自西南交通大学，由H.Y.J.领导。CR事件重建由S.Z.C.领导，IHEP、山东大学和西南交通大学的团队参与。IHEP和山东大学小组的数据分析由S.Z.C领导。；来自中山大学的P.H.T.T。；来自R.Z.Y.的中国科技大学，由张衣领导的紫色MT天文台的团队对分析进行了独立的交叉检查。由R.Z.Y.和R.Y.L.牵头解释结果，南京大学的小组由Y.C.和X.Y.W.牵头参与。；来自云南大学，由L.Zhang领导；来自紫金山天文台，由曹振英领导，F.A.A.和B.D.P.监督数据分析和解释，并领导论文撰写。S.M.L.和D.D.V.作为《拉索》编委会的联合主席做出贡献。许多小组参与了分析程序，包括由D.R.领导的马希隆大学的小组。；来自上海交通大学，由周先生领导；来自上海天文台，由Z.X.W领导。；来自北京大学，由卓立领导；来自广西大学，由E.W.L.领导。；由W.W.T.领导的国家天文台。；来自武汉大学，由谭宝华领导；来自云南科学院天文台，由中国科学院大学、天府宇宙射线研究中心等所有由甄曺领导的天文台领导；由Q.A.领导的粒子检测与电子学国家重点实验室。；来自清华大学，由Y.N.L.领导。；由Z.B.S.领导的国家空间科学中心。；来自广州大学，由J.H.F.领导。；河北师范大学南华早报。；来自郑州大学，由H.D.L.领导。；来自四川大学，由C.W.Y领导。；来自西藏大学，由T.L.C.领导。；来自核研究所和莫斯科物理与技术研究所，由俞敏洪领导。V.S.——为探测器的建造和提交前的最终版本审查做出了贡献。
红色小点表示5195闪烁体计数器，在1km2的中心区域间隔15m，在KM2A的0.3km2的边缘区域间隔30m。蓝色大圆点表示分布在中心区域的1188个μ子探测器，间距为30米。中心的三个浅蓝色矩形表示WCDA，总共78000平方米。WCDA附近的黑色小矩形显示了WFCTA的18个望远镜。
a、 在模拟容器中的事件时，抛出能量的分布为真，由10 TeV以上的重构能量E rec定义。纵轴上的分数定义为（dn/de true）/ni，其中ni是第i个料仓中的事件总数，由重建的能量E rec定义。根据100 TeV以上14%的能量分辨率，选择了Δ（log E rec）=0.2的料仓宽度（参考文献8）。与E成正比的幂律谱−此处假设为3.09。b、 料仓纯度，定义为料仓中E为真的事件的分数，作为10 TeV至2.5 PeV范围内E rec的函数。
输入（灰色实线）事件是使用根据蟹状星云8的测量SED确定的幂律SED生成的。对于每种情况，输入SED在键中列出的E cut值处都有一个人工截止（虚线）。这些分布表明，E切口上方料仓中的事件具有明显的溢出效应。能量越高，这种效应越弱。没有迹象表明，输入能量低于0.3 PeV时，1 PeV以上的垃圾箱会受到污染。
插图显示了KM2A重要性图，表明UHEγ射线源的潜在对应物。颜色栏显示重要性（\（\sqrt{{\rm{TS}}\）。绿色圆圈表示拉索的PSF。本文分析的拉萨索J1908+0621的费米LAT点，以及ARGO 48、HESS 49和HAWC 4数据，与拉萨索测量结果一起显示。虚线显示了辐射的轻子模型，根据脉冲星的自旋下降行为假设电子/正电子对的注入，断裂指数为2，初始旋转周期为0.04 s。假设脉冲星PSR J1907+0602在2.4 kpc距离处的当前自旋下降功率的6%的一小部分转换为e±对，以支持γ射线发射。电子的注入光谱假定为\（N（E）\propto{E}{E}{E}^{-1.75}\exp \{{{E}{E}/（800{rm{TeV}}}}}}}}}{2}）。固体曲线对应于辐射的强子模型。假设母体质子布居有两种类型的能量分布：（i）母体质子的单一幂律谱，N（E）≈E−1.85exp[− E/（380TeV）]（细实曲线）；（ii）母质子指数截止的破幂律谱，指数分别低于和高于25 TeV的1.2和2.7，截止能量为1.3 PeV（厚固体曲线）。在插入的天空地图中，黑色菱形显示了PSR J1907+0602的位置，黑色轮廓对应于超新星遗迹SNR G40的位置。5-0.5，白色圆圈表示HESS J1908+063的位置和尺寸。青色区域是方法中描述的密集团块。整个γ射线发射区的平均密度估计约为10厘米−3.理论曲线中考虑了光子-光子对产生的γ射线吸收（见方法）。
通过提交评论，您同意遵守我们的条款和社区指南。如果您发现有滥用或不符合我们的条款或指导方针，请将其标记为不适当。