跳到正文
至真探索 · 核结构

把扔掉的信号捡回来

2026-08-15

探测原子核内部,物理学家最顺手的工具是"听"它放出的伽马射线。可光子钻进探测器时,几乎从不一次把能量交干净。它先撞出一个电子,自己拐个弯飞走,再撞,再拐,一连好几次,直到能量耗尽,或者干脆逃出晶体。这条歪歪扭扭的散射链,在旧式探测器里是最大的噪声来源。Gammasphere 这类老阵列的办法很粗暴:在锗晶体外面围一圈 BGO 闪烁体当盾牌,凡是散射出晶体的光子,一律判成"没测干净",扔掉。伯克利实验室的一篇报道算过一笔账:被扔掉的有用信号约占八成——这正是 Deleplanque 他们后来那篇概念论文要解决的动机([1])。你装了一屋子探测器,真正派上用场的不到两成。

1999 年,劳伦斯伯克利国家实验室(LBNL)的三个人——Michelle Deleplanque、I-Yang Lee、Kai Vetter——提出一个反直觉的方案:别扔了。把每一次碰撞的位置和能量都记下来,再从这些碎片反推出光子原本的路径。这个点子([1])当时像科幻。要记录每次碰撞的精确位置,每块晶体得切成几十个独立读出的小格;还要有算法,能在几微秒里把碎片拼回完整光子。

怎么记录位置?每一块高纯锗晶体被切成 36 段——6 个方位角乘 6 层,像一根切成 36 格的圆铅笔。伽马光子在某一小格附近撞出电子,那一格的电极就跳起一个电流脉冲。脉冲的形状取决于碰撞点到电极的距离和方向。把这段波形跟预先算好的"位置对应波形"表一比,就能反推出碰撞发生在晶体里哪个点,精度到一两毫米。这一步叫脉冲形状分析(pulse shape analysis,PSA),实现时是把观测波形跟预计算基组做拟合,解出每次相互作用的数目、能量、时间和三维位置。

位置有了,接下来是拼图。追踪算法分两派。前向追踪(forward tracking):从第一次碰撞出发,按康普顿散射的运动学规律往外推,看哪条链的能量和角度对得上,用一个打分函数选最像的那条([2])。后向追踪(backtracking):从最后一次碰撞反推回来——最后一次大概率是光电吸收,能量干净([3])。GRETINA 的追踪代码叫 GTAG,由阿贡国家实验室维护,里面塞满了"两个碰撞点太近要不要合并""一条链分叉了怎么办"这类修正。这套办法的完整技术综述,见 Korichi 和 Lauritsen 2019 年的长文([7])。活儿极耗算力:每块晶体每秒要分解两三千个事件,整套装置得配几千核的计算机群。

代价大,回报也大。老阵列扔掉的那八成信号,追踪阵列全留着。GRETINA 是 GRETA 的四分之一规模先导,一开始 7 个四联模块、28 块晶体,覆盖约四分之一的球面,后来扩到 12 个模块。它 2013 年的性能报告([4])给出了关键数字:做完信号分解加追踪,峰总比超过 0.61,同时保住原始效率的 75%。峰总比是"测到的事件里有多少是干净的光峰事件",这个数越高,底下的噪声越少。老阵列 Gammasphere 能干净测到的事件约占 9%,但那靠丢掉八成信号换来的。追踪阵列不丢信号,也能达到相近的峰总比——同一份物理,它能看见弱得多的信号。两套思路的系统对比,见 Lauritsen 等人 2016 年的论文([5])。

真正的革命发生在物理上。先看梨形核。梨形核,就是原子核长得一头大一头小、左右不对称,像一只梨。2013 年,Peter Butler 和 Liam Gaffney 领导的团队在 CERN 的 ISOLDE 装置上,用 Miniball 阵列做库仑激发——让带电粒子擦着核飞过,用电磁场把它激起来——然后量激发时放出的伽马([10])。他们报告:氡-220 的 B(E3)=32±4 韦斯科普夫单位,镭-224 的 B(E3)=42±3。这里 B(E3) 是衡量电八极跃迁强弱的量,韦斯科普夫单位是核物理给跃迁强度定的标准尺;电八极跃迁强,就是"梨形"的直接指纹。镭-224 因此被判成静态梨形,氡-220 被判成八极振动——它也在"梨"和"不梨"之间来回晃。这个"静态还是振动"的区分,至今还是密度泛函理论家们的争论点:不同泛函对哪个核是静态梨形,给不出统一答案。

美国这边,GRETINA 把梨形核的版图推到了中子丰的钡同位素区。2016 年,Brian Bucher 和他的合作者在阿贡国家实验室做实验:用 CARIBU 装置产生钡-144,加速后打在铅-208 靶上做库仑激发,用 GRETINA 测放出的伽马([8])。他们量到基态到第一个 3⁻ 态的 B(E3) = 48(+25/−34) 韦斯科普夫单位——比当时的理论预言普遍偏大。在中子丰的钡同位素区,这是第一个被实验直接钉死的梨形核。第二年他们又对钡-146 做了同样测量,B(E3) = 48(+21/−29),八极形变参数 β₃ 约 0.17——β₃ 是描述"梨"有多歪的参数([9])。更妙的是,钡-144 和钡-146 的内禀电偶极矩差了一个数量级以上,E3 却几乎一样。这说明两边的"梨"差不多大,电偶极的巨大差异来自反射不对称平均场里不同中子轨道怎么占位。教科书里"静态八极形变"这一栏,从此有了干净的实验样本。

再看奇特核。GRETINA 早期最漂亮的成果之一,是飞行中测到镁-32 的第二个 0⁺ 态。镁-32 是"反转岛"上的著名核,中子数 20 的壳层间隙在这里塌掉了。GRETINA 测到那个激发态的寿命,发现它的跃迁强度大到当时的壳模型怎么也算不出来([11])。壳模型,就是把核子按壳层排列的模型。还有硅-34:2017 年 Nature Physics 上的质子敲出实验([25])显示,它的质子 s 轨道近乎空着,相邻的硫-36 里却几乎占满——自旋轨道相互作用对密度和同位旋的依赖,第一次这么直观。2025 年,Alexandra Gade 团队在 FRIB 用 GRETINA 做了第一个在束实验:从铁-64 敲出两个质子,产生铬-62,量到它的形状共存——同一颗核里,两个不同形变的激发带挤在一起([12])。这组数据把壳模型推到了组态混合的极限。

理论那边,最让 GRETA 眼馋的是钙同位素。双重幻核,就是质子和中子数都恰好是幻数、特别稳定的核。钙-60 是不是双重幻核?Gade 团队从铬-62 外推,预言它的第一激发 2⁺ 态约在 1 兆电子伏,跟 N=36/38 的同中子素相当;另一些理论家([13])则对钙-60 的双重幻性提出质疑。两个预言正面冲突。从头算(ab initio,从基本核力出发直接算)还预言硅-48 是双重幻核,第一激发 2⁺ 态约 2.49 兆电子伏,比钙-54 还高约 400 千电子伏([26])——这些预言都还没被实验检验。钙的滴线在哪,吵得更凶。旧的从头算说中子数到 40 就停了,密度泛函普遍说能到 50 上下(各泛函给出的区间在 48 到 56 之间);2025 年的新 N³LO_Texas 手征力([14],手征力是描述核力的有效理论)把从头算的滴线一路推到中子数 51,跟密度泛函和解了,但这个预言同样没被验证。GRETA 的科学案卷里写着,钙-60 的第一次谱学只有它的分辨本领才能实现——这话带点自我宣传的味道,但方向没错。

高自旋物理也一样。超形变带的原子核被拉成橄榄球,长轴差不多是短轴的两倍。但超形变带怎么退回正常形状,一直是个谜:连接跃迁极弱,被统计涨落打碎,旧阵列根本看不清。GRETA 的高效率和高峰总比,正是为这种弱信号准备的。磁转动、反磁转动、手征双重带——磁转动是一群核子像剪刀一样绕主轴转动,手征双重带是同一颗核里出现成对的、互为镜像的转动带——这些靠测磁跃迁和电跃迁的强度比、能级寿命和四极矩,来区分"手征振动还是静态手征"的物理。这类测量要同时看两个光子先后落地、光子往哪个方向飞、以及光子的偏振——这三样恰恰是追踪阵列的强项。

现在说 GRETA 本体。30 个四联模块,120 块 36 段锗晶体,覆盖约八成的球面。设计指标:把散射能量加回去后的全能量效率不低于 30%(1.33 兆电子伏),峰总比不低于 50%,位置分辨好于 2 毫米,能量分辨好于 2.5 千电子伏,每块晶体每秒能扛 5 万次计数。阵列由 LBNL 牵头建造,造价 5800 万美元,2025 年 8 月主体完工,秋天运到密歇根州立大学的 FRIB。2026 年 7 月 28 日,LBNL 宣布它在 FRIB 的再加速束流线上完成首次束流测试,30 个模块全部就位,时间同步和能量分辨都验证过了([20])。8 月 5 日,萨里大学宣布第一个正式科学实验已经开跑([21]):由该校博士后 Jacob Heery 领导,用 JANUS 合作组(FRIB、LBNL、萨里大学)的方案,测量放射性锗同位素的形状。理论一直预言锗是长椭球,近年的线索却暗示可能是扁椭球。这个看似基础的测量,若结果偏向扁椭球,就会挑战现有对质子中子行为的理解。

这套装置还准备搬家。它设计之初就可移动,计划 2028/2029 年从 FRIB 搬到阿贡国家实验室的 ATLAS 装置,接手即将退役的 Gammasphere 的用户群,重点做超形变和超重核。

那"灵敏度高 10-100 倍"的说法是哪来的?得分口径。1999 年那篇概念论文的原话,是 GRETA 会比 Gammasphere、Euroball "强约一千倍"——按分辨本领算的。当代官方口径(FRIB 仪器页、ANL 讲座)是"灵敏度提升约 10-100 倍"。中间还有更细的说法:对高自旋级联这类高多重性事件,提升 100-1000 倍;对低多重性事件,约 10 倍。引用时得说清楚是哪个口径,否则容易吵起来。

欧洲那边,AGATA 是另一种节奏。它跟 GRETA 技术同源:同样 36 段锗晶体,同样脉冲形状分析加追踪([22])。差别在模块:GRETA 用四联模块,AGATA 用三联模块;晶体读出方向、时钟约定也不一样,两套数据格式互不兼容,软件要翻译才能互跑。AGATA 的完整设计是 180 块晶体、60 个三联模块、约 362 千克锗。但它的建造是分期走的:2010-2011 在意大利 LNL 演示,2012-2014 搬到德国 GSI 做相对论束谱学,2014-2021 在法国 GANIL 打了几轮战役、最多时约 60 块晶体(三分之一的球面),2022 年回到 LNL 重组。现在在役的是约 15 个三联模块(45 块晶体),覆盖约四分之一的球面——还没回到 GANIL 时期的规模。规划是先到 2π,再往 3π 和 4π 爬。

AGATA 的困难是钱。2024 年它的指导委员会记录显示,德国教研部(BMBF)对 AGATA 的拨款申请评估推迟,委员会不得不写支持信,GSI 先垫钱买了一台低温恒温器。一期协议(2021-2025)的预算缺口是 12 个探测器胶囊、1 台低温恒温器加配套电子学。这不是小数目。

所以 GRETA 和 AGATA 的关系,是竞争,也是合作。竞争是明摆着的:同一技术概念,同一批物理目标——超重核、奇异形变、远离稳定线的谱学,各自被欧洲的 NuPECC 长期规划和美国核物理长期规划列为旗舰。但合作也是实打实的:两边开系列联合会议,2024 年 11 月第四届在阿贡举行,由法国 CSNSM 的 Amel Korichi 和阿贡的 Torben Lauritsen 组织,议题包括双方装置状态、信号分解、AI/机器学习、模拟基准。GRETA 的信号生成直接用了 AGATA 的探测器库;UCGretina 模拟程序被改造成生成 AGATA 的晶体基函数;两边还真做过"用 AGATA 数据跑 GRETINA 的分解""用 GRETINA 数据跑 AGATA 的脉冲形状分析"。没有一纸正式合作协议,靠的是定期见面和共享软件——同一张图纸上的两座独立工程,互相审计,也互相较劲。

全球版图上没有第三套完整的 4π 追踪阵列。日本 RIKEN 在借力:GT-5 项目打算买一套 GRETINA 四联模块,用在快速放射性束的多普勒修正谱学([23])。加拿大和中国都没有公开的全尺寸追踪阵列计划。

说到中国,画面要分两层。第一层:没有任何中国机构是 GRETA 或 AGATA 的正式成员。用 INSPIRE 数据库检索,24 篇标题含 GRETINA 的论文、101 篇标题含 AGATA 的论文,署中国机构的为零。第二层:中国物理学家并没有缺席这件事的物理。北京大学孟杰团队预言的"多重手征双重带"(MχD,一个核里出现多对手征伙伴带),2013 年在阿贡用 Gammasphere 被证实([16]),孟杰、张双全、赵鹏巍、陈启博、亓斌是共同作者——用的不是追踪阵列,但物理正是 GRETA 未来要做的物理。山东大学(威海)王守宇团队 2024 年在 PRL 报告了溴-74 里手征带和摇摆带共存(实验在南非 iThemba 实验室,8 个 clover 探测器,[17]),也是传统阵列的活儿,同样属于 GRETA 的目标物理。复旦大学王思敏、博士后许志成和 FRIB 的 Witold Nazarewicz 合作,用伽莫夫壳模型解释了钙-36 的 B(E2) 反常([18])——这是理论进入 GRETA 目标物理的现成通道。中国团队还通过美国装置的"用户机制"直接做实验:北京大学李智焕、华辉团队 2021 年领衔在 NSCL 首次观测到四质子非束缚的镁-18([19]),虽然用的不是追踪阵列。

还有一层容易混淆:华裔科学家在美国机构里的参与。Shaofei Zhu 在阿贡,是 GRETINA 项目核心成员,那两篇钡同位素八极论文([8][9])的作者列表里就有他,还有罗切斯特大学的 Ching-Yen Wu。但他们隶属美国机构,不等于中国团队。

中国自己的伽马谱学装置,走的是另一条路线。北京的 CGAC(中国联合伽马阵列,2019 年由原子能院、北大、山大、吉大、清华、深圳大学等共建,[24])有 40 台高纯锗探测器,效率约 3%,是传统阵列的思路。惠州 HIAF 的高能放射性束终端规划了一个 4π 碘化铯闪烁体阵列(CsI(Tl)),走"高流强重离子加相对论在束伽马"的路线——跟 GRETA/AGATA 的锗追踪互补,不是追赶。差距是结构性的:国内没有位置灵敏分段锗的追踪阵列,脉冲形状分析依赖欧洲的方法,高纯锗晶体材料部分靠进口,核电子学 ASIC 刚起步。

开放问题还不少。快束多普勒修正是个极限:原子核以超过三成光速飞行时放出的伽马,能量会被多普勒效应拖宽,修正精度卡在位置分辨上——GRETINA 在快束上的表现,见 Weisshaar 等人的论文([6])。低能段(百千电子伏以下)和极高计数率下的信号堆叠,追踪阵列仍是短板。中子损伤会让锗晶体的位置分辨慢慢变差,4π 规模的长期运维成本没人说得清。AI 能不能让信号分解和追踪更快更稳,第四届联合会议把它列成专题,但还没有定论。还有最根本的物理争论:梨形核里"静态"和"振动"的边界在哪(八极形变的理论全景,见 Cao 等人的 PRC 102, 024311,[15]),钙-60 到底是不是双重幻核,钙的滴线在哪——这些都要等 GRETA 的数据。

说到底,这台机器能做的,是把核物理过去几十年丢掉的那八成信号,重新变成数据。2026 年 8 月,120 块锗晶体已经在 FRIB 开始正式工作;AGATA 还在经费的泥潭里往 2π 爬。两套装置共享算法、互相审计,在同一批核上竞逐。中国人在场——理论、传统阵列、还有在美欧机构的华裔科学家——中国机构却不在两大合作体系里。梨形核的边界、钙-60 的幻性、滴线的位置,这些问题的答案,都要等新数据说话。

*参考文献*

  1. M. A. Deleplanque, I-Y. Lee, K. Vetter, et al., "GRETA: utilizing new concepts in γ-ray detection," *Nucl. Instrum. Methods A* 430, 292 (1999). (LBNL)
  2. G. J. Schmid, M. A. Deleplanque, I-Y. Lee, F. S. Stephens, K. Vetter, et al., "A γ-ray tracking algorithm for the GRETA spectrometer," *Nucl. Instrum. Methods A* 430, 69 (1999). (LBNL/LLNL)
  3. J. van der Marel, B. Cederwall, "γ-ray tracking in germanium: the backtracking method," *Nucl. Instrum. Methods A* 477, 391 (2002). (KTH, Stockholm)
  4. S. Paschalis, et al., "The performance of GRETINA," *Nucl. Instrum. Methods A* 709, 44 (2013). (LBNL)
  5. T. Lauritsen, et al., "Characterization of a gamma-ray tracking array: a comparison of GRETINA and Gammasphere using a ⁶⁰Co source," *Nucl. Instrum. Methods A* 836, 46 (2016). (ANL 等)
  6. D. Weisshaar, et al., "The performance of the γ-ray tracking array GRETINA for γ-ray spectroscopy with fast beams of rare isotopes," *Nucl. Instrum. Methods A* 847, 187 (2017). (NSCL/MSU 等)
  7. A. Korichi, T. Lauritsen, "Tracking γ rays in highly segmented HPGe detectors: A review of AGATA and GRETINA," *Eur. Phys. J. A* 55, 121 (2019). (CSNSM/IN2P3; ANL)
  8. B. Bucher, et al., "Direct evidence of octupole deformation in neutron-rich ¹⁴⁴Ba," *Phys. Rev. Lett.* 116, 112503 (2016); arXiv:1602.01485. (LLNL/ANL/LBNL/U. Rochester/U. Liverpool 等)
  9. B. Bucher, et al., "Direct evidence of octupole deformation in neutron-rich ¹⁴⁶Ba," *Phys. Rev. Lett.* 118, 152504 (2017); arXiv:1703.05268. (含 U. Autónoma de Madrid 理论组)
  10. L. P. Gaffney, P. A. Butler, et al., "Studies of pear-shaped nuclei using accelerated radioactive beams," *Nature* 497, 199 (2013). (U. Liverpool, CERN ISOLDE, KU Leuven 等)
  11. P. Fallon, A. Gade, I-Y. Lee, "GRETINA and Its Early Science," *Annu. Rev. Nucl. Part. Sci.* 66, 321 (2016). (LBNL; NSCL/MSU)
  12. A. Gade, et al., "In-beam spectroscopy reveals competing nuclear shapes in the rare isotope ⁶²Cr," *Nature Physics* 21, 37 (2025), DOI 10.1038/s41567-024-02680-0. (FRIB/MSU 等)
  13. S. Chen, et al., "Level structures of ⁵⁶,⁵⁸Ca cast doubt on a doubly magic ⁶⁰Ca," *Phys. Lett. B* 843, 138025 (2023). (钙-60 双重幻性质疑)
  14. B. Hu, A. Ekström, et al., "The neutron dripline in calcium isotopes from a chiral interaction," arXiv:2512.11723 (2025). (N³LO_Texas 手征相互作用)
  15. Y. Cao, S. E. Agbemava, A. V. Afanasjev, W. Nazarewicz, E. Olsen, "Landscape of pear-shaped even-even nuclei," *Phys. Rev. C* 102, 024311 (2020); arXiv:2004.01319. (MSU/Mississippi State)
  16. A. D. Ayangeakaa, U. Garg, et al., "Evidence for Multiple Chiral Doublet Bands in ¹³³Ce," *Phys. Rev. Lett.* 110, 172504 (2013). (Notre Dame 主导; Argonne Gammasphere; 孟杰/张双全/赵鹏巍等为共同作者)
  17. R. J. Guo, S. Y. Wang, et al., "Evidence for Chiral Wobbler in Nuclei," *Phys. Rev. Lett.* 132, 092501 (2024). (山东大学威海; iThemba LABS)
  18. Z. C. Xu, S. M. Wang, W. Nazarewicz, et al., "Puzzling B(E2; 0⁺→2⁺) strength in the proton dripline nucleus ³⁶Ca," *Phys. Rev. C* 112, L011302 (2025). (复旦大学; FRIB/MSU)
  19. Y. Jin, C. Y. Niu, Z. H. Li, H. Hua, et al., "First Observation of the Four-Proton Unbound Nucleus ¹⁸Mg," *Phys. Rev. Lett.* 127, 262502 (2021). (北京大学; NSCL/MSU)
  20. LBNL News Center, "GRETA Comes Online for a New Era of Nuclear Science," 2026-07-28. newscenter.lbl.gov
  21. University of Surrey, "Surrey leads first nuclear science experiment with world's most advanced gamma-ray detector," 2026-08-05.
  22. AGATA Collaboration, 官方网站 agata.org(成员列表、部署历史、当前配置).
  23. RIKEN RIBF, GT-5 项目提案 (APR057), nishina.riken.jp.
  24. W. Nan, B. Guo, et al., BRIF 伽马谱学综述(含 CGAC),arXiv:2503.17384. (CIAE 等)
  25. M. Mutschler, et al., "A proton density bubble in the doubly magic ³⁴Si nucleus," *Nature Physics* 13, 152 (2017). (含 A. Gade, D. Weisshaar, J. A. Tostevin)
  26. Q. Yuan, et al., "Ab initio study of shell evolution in neutron-rich Si, S, Ar, and Ca isotopes near N = 32 and 34," *Chinese Physics C* 49, 094104 (2025), DOI 10.1088/1674-1137/add5dd. (⁴⁸Si 双重幻核的 VS-IMSRG 预言;中科院近物所等)