电荷交换反应与自旋-同位旋激发

把中子变成质子。不是等它自己衰变——那要等几分钟到几百万年——而是拿一束粒子撞上去，在碰撞的那一瞬间，把中子踢成质子。核物理学家干了这件事干了快五十年。它的正式名字叫电荷交换反应。

听上去像魔术。其实原理不复杂。你拿一个质子去打靶核，质子飞过去的时候，和靶核里的一个中子交换了一个带电的π介子——质子的正电荷传给了中子，中子变成了质子，质子自己变成了中子。结果：靶核的原子序数加了1，出射粒子是一个中子。这就是(p,n)电荷交换反应。

换一个弹丸。(³He,t)反应——氦-3打进去，氚飞出来。两个质子加一个中子变成了一个质子加两个中子。核心还是同一个机制：一个中子和一个质子互换了身份。

换一个方向。(d,²He)反应——氘核打进去，两个质子飞出来。这次是靶核里一个质子变成了中子——靶核的原子序数减1。方向相反，但物理一样。

还有(t,³He)、(⁷Li,⁷Be)、(¹²C,¹²N)、(¹³C,¹³N)。这些都是电荷交换反应的变体。弹丸不同，探针的灵敏度不同——有些对自旋翻转敏感，有些对空间结构敏感——但核心都是同一个东西：把原子核里的一个中子变成质子，或者反过来，然后观察核的响应。

为什么费这么大劲？因为这种"强制变身"能告诉我们原子核最深层的秘密之一：自旋-同位旋响应。

拆开来说。原子核的激发模式不只一种。有些是形状振动——整个核被挤一下，它扁了又圆、圆了又扁。有些是密度振动——核被压缩，密度变高又变低。电荷交换反应激发的，是第三种：同位旋翻转。它的物理图景是：在核里选中一个位置，把那个位置上的中子换成质子。核的其他部分——所有其他核子——对这次"换人"怎么反应？

这就是自旋-同位旋激发。它测量的是核内部一种特殊的"集体刚度"——如果你在一个点上翻转了同位旋，整个核有多少核子会跟着一起翻转？翻转的"波"能传多远？

最能体现这种集体响应的现象，叫伽莫夫-泰勒共振。简称 GT 共振。1975 年，在印第安纳大学回旋加速器上，用(p,n)反应首次在中等质量核中看到了 GT 共振——一个宽大的、集体性的峰，出现在比基态高 8 到 15 MeV 的地方。它的物理很简单：核里的中子和质子同时集体性地翻转自旋，一个中子变成质子，整个核因此升高了大约 10 MeV 的能量。这是一个"共振"——不是某几个核子在动，而是一大群核子同步动。

GT 共振的发现解决了一个四十年的谜：弱相互作用里的"淬火"谜。中子自由衰变时，轴矢量耦合常数 gA ≈ 1.27。但在核内——当你把核视为一个整体时——同样的 gA 表现得好像只有约 1.0。差了约 25%。这个差值的物理原因是什么？有两种答案在 2024-2025 年激烈交锋。

第一种答案：多体关联。2025 年 9 月，一篇基于大规模壳模型计算（arXiv:2509.24542）的论文以锗-76 为例论证：强核力关联——包括形变、跨壳混合、高激发态的混入——把大量 GT 强度推到了极高能区，超出了实验窗口。你测得的那部分——你以为是 100% 的——其实只有 70-80%。不是 gA 变小了，是你没看见剩下的 20-30%。他们的计算显示，手征两体流（来自两个核子同时参与弱相互作用）的贡献只有 5-15%。

第二种答案：QCD 介质效应。Mannque Rho 2024-2025 年在 Modern Physics Letters A 上发表了一系列分析，论证手征有效场论在核介质中漏掉了两个关键效应：标度对称性破缺（trace anomaly / dilaton）和 Brown-Rho 标度——强子质量在核介质中会下降。把这两个效应加回去，gA 的淬火因子精确等于 0.78——无需任何自由参数。Rho 进一步指出，最近的锡-100 超允许 GT 衰变测量如果被证实，可能需要一个更大的、更"基本"的淬火——gA 在介质中降到自由值的一半。

2025 年 3 月，Suzuki 和 Shimizu（东京大学）在 PCM2025 会议上报告了 VS-IMSRG 结果。他们用 VS-IMSRG 的方法从手征核力出发，同时处理壳模型空间截断和两体流效应。在 sd 壳区，从 sd 拓展到 sd⁻²pf² 甚至 sd⁻⁴pf⁴ 的模型空间，逐步回收了约 0.05 的 gA 因子。两体流贡献约 0.07。两者加起来，有效 gA 从约 0.77 升到了约 0.89——距离 1.0 还差一截，但方向是对的。

GT 淬火的答案不是非此即彼。多体关联确实推走了大量强度。两体流确实贡献了一个不可忽略的修正。介质效应确实存在——但它的精确量级仍然受限于我们对核介质中手征微扰论的可靠性评估。这个三方平衡，是 2025 年 GT 淬火研究的核心地带。

GT 共振不是电荷交换反应的唯一产出。更重的伙伴叫自旋-偶极共振——简称 SDR。SDR 比 GT 共振多了一个单位的轨道角动量（ΔL=1 而不是 ΔL=0），所以激发能更高——通常在 20-30 MeV。SDR 对核结构特别敏感。它的总强度分布在三个自旋-宇称道上——0⁻、1⁻、2⁻——这三个道的重心能量差，是张量力的直接探针。

2000 年代以来，RCNP（大阪大学核物理研究中心）的 Grand Raiden 谱仪用 450 MeV 的(³He,t)反应，成为全球高分辨率电荷交换实验的标准。Grand Raiden 的能量分辨约 30 keV——足够在 GT 共振的宽峰里分辨出一个个离散的末态。它的零度测量能力和分散匹配光学系统，至今没有对手。

Grand Raiden 跑了几十年。最近一次重要产出是(³He,t)对锡同位素链的系统扫描。对于¹¹⁶,¹²²Sn(³He,t)反应，Ikeda 求和规则——3(N-Z)——的实测比例只有 29-38%。这意味着在锡区，超过 60% 的 GT 强度"丢了"——移到了更高能区，藏在准自由散射本底下。这也是 GT 淬火的最直接证据之一。

另一个方向也在快速推进——用重离子做电荷交换反应。这一派的物理目标完全不同：测无中微子双贝塔衰变的核矩阵元。

无中微子双贝塔衰变——如果存在——会证明中微子是自己的反粒子。但实验家只能在极深的地下实验室等两个中子同时衰变成两个质子，放出两个电子，没有中微子。等上十年，如果没等到，你只能说：矩阵元小于某个值。要翻译成中微子的有效质量，你需要知道核矩阵元。算不准。

NUMEN 项目——NUclear Matrix Elements for Neutrinoless double beta decay——换了一条路。它不用等两个中子一起衰变。它直接做双电荷交换反应：用重离子（如氧-18）打靶，把靶核里的两个中子同时变成两个质子，测量出射道的氧-18→氖-18 双电荷交换产物的角分布。物理逻辑：双贝塔衰变的末态波函数——初始核和最终核的波函数关系——和双电荷交换反应完全一样。只是弱相互作用换成了强相互作用。

2024 年，NUMEN 在 INFN-LNS（意大利卡塔尼亚）交出了几个关键数据。¹¹⁶Sn(¹⁸O,¹⁸Ne)¹¹⁶Cd 的双电荷交换绝对截面——在 276 MeV 束流能量、极前角——测出来只有 34 nb（纳靶）。这得有多小？一纳靶是十的负三十三次方平方厘米。一百亿亿个氧-18 打过去，只产生几个事件。但正因为截面这么小——且在不同核系统间变化不大——它的系统性测量变得有价值：如果你能在多个 ββ 候选核上测到同样的截面，你就可以对比理论计算和实验数据，提取核矩阵元的相对约束。

NUMEN 也在做单电荷交换和转移反应——同一个¹⁸O+⁴⁸Ti 系统，同时测了(¹⁸O,¹⁸F)、(¹⁸O,¹⁷O)、(¹⁸O,¹⁶O)等所有道，以便从反应机制上分离核结构效应。2024 年 Spatafora、Cappuzzello 等在 PRC 上报道了¹²C(¹⁸O,¹⁸F)¹²B 的道竞争分析——介子交换直接过程和连续转移过程——两者幅度同量级，必须相干叠加才能解释角分布。

NUMEN 的基础设施升级也在 2024-2025 年推进。超导回旋加速器的束流强度要从 10¹⁰ pps 提升到 10¹³ pps——三个数量级。新的反应靶室在 2025 年开始安装。新的前角探测器——基于多层 THGEM 气体倍增器的低气压径迹探测器——完成了原型测试。这些升级的目标很简单：把 34 nb 的截面从"勉强可测"变成"可以系统测"。

回到理论。描述自旋-同位旋响应的标准工具是质子-中子准粒子无规相位近似——简称 pnQRPA。这个名字长得不像人话，但意思不复杂：你先用平均场理论算出核的基态——所有核子填在最好的单粒子轨道上——然后问一个简单的问题：如果我把一个中子变成一个质子（或反过来），这个"粒子-空穴对"会如何演化？是这个对单独振荡，还是会耦合出一大群对的集体运动？

pnQRPA 回答的就是这个问题。它把所有的中子-质子对放在一个矩阵里，对角化这个矩阵，找出"集体模式"——那个把很多很多对同步起来的振动。GT 共振就是 1⁺ 集体模式。SDR 是 1⁻ 集体模式。IAS（同位旋相似态）——把中子变质子但不翻转自旋——是 0⁺ 集体模式。

2024 年 8 月，Ravlić 等人在 Physical Review C 上发表了第一个变形核相对论 pnQRPA。之前的 pnQRPA 都假设核是球形的——对于变形核（橄榄球形的），你只能假装它是圆的。这篇论文打破了球对称限制——在轴向变形下，用相对论点耦合能量密度泛函（DD-PCX）和可分离对力，实现了完整的 pnQRPA 计算。

结果呢？变形导致 GT 强度的剧烈碎裂。在球形极限下，GT 共振是一个干净的峰——能量约 10 MeV，宽度约 2-3 MeV。一旦核变成橄榄球形，这个单一的峰就碎成了一大片——因为变形打破了旋转对称性，不同的 K 投影（角动量在对称轴上的分量）产生了物理上不同的激发模式。你看到的不是"一个 GT 共振"，而是三个——K=0、K=±1——能量分开，强度分散。

同一年，Sagawa 等人发展了减除二次 RPA（SSRPA）——在 pnQRPA 的基础上加了二粒子-二空穴（2p-2h）激发。pnQRPA 只包含 1p-1h——一次变一个人。但在真实核里，一次变了两个人也不是不可能——把两个中子变成两个质子，或者把一个中子变质子、另一个核子跳到更高的轨道。SSRPA 把这种"双倍激发"自洽地加进去。结果是：GT 共振的低能尾巴被显著增强，贝塔衰变半衰期因此被缩短——²²Si、⁶⁸Ni、⁷⁸Ni、¹³²Sn 的半衰期，有张量力比没有张量力差一到两个数量级。

SSRPA 不需要引入自由参数来解释淬火。淬火"自然"出来了——因为 2p-2h 耦合把强度从低能区推到高能区，而你做的实验只看低能区。这不是真正的 gA 变小了，是你的截断。

Robin 和 Litvinova 在 2024 年的另一条线上追同一个问题——他们算的是"电荷交换声子"的耦合。在 pnQRPA 里，集体模式一旦被找到，就被当成独立的声子——彼此不通信。但真实情况是，一个单粒子可以同时和一个 GT 声子耦合，再和一个 SDR 声子耦合。这种耦合——技术上叫粒子-振动耦合——会把单粒子态的能级移来移去，同时改变 GT 共振的形状和位置。Robin 和 Litvinova 的计算显示，加上这种耦合后，不需要调整 gA，就能精确再现⁷⁸Ni 和¹³²Sn 的贝塔衰变半衰期。

这就引出了一个深层问题：GT 淬火到底是一场误会吗？如果我们把所有的多体关联都算进去——跨壳激发、2p-2h、声子耦合、形变效应——而实验窗口永远只能看到总强度的一部分，那所谓"淬火"也许不是 gA 在核介质中变小了，而是我们的理论太小了。

2025 年，两篇论文给这场争论添了新柴。一篇是 Sagawa、Yang、Bai 在 Acta Physica Polonica B 上的 SSRPA 扩展——他们把 SSRPA 的"减除"技术和 Lee-Suzuki 相似变换建立了形式联系。另一篇是 Ejiri（RCNP）在 AIP 会议文集上的 pnQRPA 系统学——用了实验上的 SD 共振能量来标定 pnQRPA 的粒子-空穴耦合强度 gph，得到的 0νββ 核矩阵元随质量数增加而系统下降，反映了自旋-同位旋极化的负相位贡献。

这条线穿透出来，是一条简洁的物理结论：核矩阵元不是核越大就越大的。自旋-同位旋响应的集体性强弱——极化程度——决定了矩阵元的大小。集体性越强——越多核子同步运动——极化越大，矩阵元反而越小。这是一个反直觉但被各种计算反复确认的规律。

现在说中国团队的工作。2025 年 7 月，北航孙保华课题组和张寂潮（现大阪大学博士后）、谷畑勇夫（RIKEN）合作，在 Physical Review X 上发表了一篇被国家自然科学基金委重点报道的论文。他们做了什么事？

电荷拾取反应——电荷交换的一种，弹丸从靶核里"捡"一个质子变成自己的中子——在 GSI 的 900A MeV 束流上，用 24 种 p 壳轻核打碳靶和氢靶。测了它们的电荷交换截面。

结果惊人。从碳-12 到碳-19，截面涨了约 50 倍。按稳定核的经验公式，只应该涨 2.5 倍。五十倍——这是指数增长。原因是：中子越富余，可用的超允许 GT 跃迁通道越多——中子数每增加一个，通道数不是加一个，而是倍增。这个效应在稳定核段几乎看不见——因为中子数太少——但到了滴线附近就爆炸了。

更重要的，孙保华他们提出了一个新的全局经验公式——把质量数和同位旋分开描述——首次在四个数量级的跨度上统一了所有电荷交换截面的测量。GSI 的实验数据、BNL 的老数据、日本的碎片数据——全部落在同一条曲线上。这意味着电荷交换反应有一个普适的标度律——之前没人猜到——而这个标度律的核心就是自旋-同位旋超允许跃迁的态密度。

这个发现不仅对理解中高能重离子反应机制重要。它直接关联到贝塔衰变强度的淬火问题：如果中子富余核的 GT 通道密度这么高——指数增长——那在滴线附近，淬火因子是不是也应该随着同位旋变？当前所有的 r-过程网络计算都假设淬火因子是常数——一个核常数、所有核一个值。如果它不是常数呢？如果在滴线附近，gA 的有效值因为态密度爆炸而变化呢？

没有人敢肯定。但孙保华的数据——50 倍的截面增长——让这个问题的温度升高了。

电荷交换反应的另一个主战场在天体物理。一颗大质量恒星坍缩时，核心温度飙升到几百亿度，电子被质子俘获——反贝塔衰变——中子星的核心因此富中子化。俘获率多大？它由 GT 强度——也就是原子核对"中子变质子"和"质子变中子"两种电荷交换模式的响应——直接影响。铁区核的电子俘获率至今是超新星模拟中最大的核物理不确定性之一。

2025 年，Song、Li 等（中国原子能科学研究院 CIAE）在 Chinese Physics C 上发表了一篇——用核总理论（Gross Theory）系统计算了¹²C、¹⁶O、⁵⁶Fe、²⁰⁸Pb 的中微子-核电荷交换截面，并提出了一个半经验的参数化公式——将中微子有效温度和截面简洁地打包在一起——目的是给天体网络计算提供一个快速数据库。

同年的另一篇，Kasurinen、Pirinen、Suhonen（Jyvaskyla大学）在 PRC 111 上首次把介子交换两体流加进了氙同位素（¹²⁸,¹³⁰,¹³²,¹³⁴,¹³⁶Xe）的超新星中微子截面计算——这是两体流从双贝塔衰变矩阵元计算中的自然延伸。结果显示：含两体流的截面和用唯象淬火因子的截面在大小上相当——但能量分布不同。这个差异可能影响到超新星中微子探测器中重建中微子能量的精度。

中国在电荷交换反应理论方面也在加速。孟杰团队（北大）2024 年在 Physics Letters B 发表了用相对论组态相互作用 DFT 算出的双 GT 矩阵元和 0νββ 矩阵元之间的强线性关联——这个关联是 NUMEN 项目的理论基础之一。2025 年，Chen、Li、Zhang、Niu 把相对论准粒子有限幅度方法（QFAM）拓展到了电荷交换道——这是第一次在变形核中自洽地处理 IAS 和 GT 共振的 K 劈裂——为 Zn 同位素链的电荷交换激发做出了理论预言。

最后，一个物理学上极深的问题，正从电荷交换的数据中慢慢浮出来。

QCD 的基础相互作用是同位旋对称的——上夸克和下夸克之间，强相互作用完全一样。但在核物理里，同位旋对称性破了。破的不是 QCD——是上夸克和下夸克的质量不一样。这个微小的质量差——约 3 MeV——通过复杂的 QCD 动力学放大，产生了核物理里所有电荷交换反应的分支比不对称、镜像核的能级差、以及 IAS 相对于基态的能量位移。

2024 年，Sagawa、Naito、Roca-Maza、Hatsuda 从 QCD 求和规则出发，导出了一个基于第一性原理的电荷对称性破缺（CSB）能量密度泛函。他们证明了 CSB 力有一个独特的径向依赖性——在核表面达到峰值——并且通过与张量力的干涉，显著影响 GT 共振的精细结构和双贝塔衰变的核矩阵元。2025 年 1 月，PRL 134 报道了⁷⁸Zr 的首次激发态测量——一个极端缺中子核——显示的能级结构和基于标准同位旋对称性的预测有显著偏差。

这一切最终指向同一个问题：同位旋对称性在核介质中以什么方式被打破？在一个由强相互作用主导的多体系统中，上夸克和下夸克之间那一点点质量差——占核子质量的 0.3%——是怎样被放大的？电荷交换反应是回答这个问题最锋利的刀子。因为它直接操作同位旋——把一个中子变成质子，测量核的集体响应——所以你测到的不是一两个核子的同位旋破缺，而是整个核的。

一把刀。但刀磨了五十年，还没磨够。

参考文献

C. Robin, E. Litvinova, arXiv:1806.01409 (updated 2024) — Coupling charge-exchange vibrations to nucleons: GT quenching and beta-decay half-lives without gA renormalization (Western Michigan Univ./FRIB/GANIL)
H. Sagawa, C.L. Bai, Z.X. Yang, Acta Phys. Pol. B Proc. Suppl. 18, 2-A32 (2025) — SSRPA for GT and β-decay: 2p-2h coupling + tensor force, Lee-Suzuki connection (RIKEN/Univ. of Aizu)
Suzuki, Shimizu, PCM2025 (Tokyo, March 2025) — GT quenching from VS-IMSRG + extended model spaces: sd⁻²pf² and sd⁻⁴pf⁴, two-body currents contribution ~0.07 (Univ. of Tokyo/CNS)
M. Rho, Mod. Phys. Lett. A 40, 2550048 (2025) — "Quenched gA puzzle in nuclei and nuclear matter and pseudo-conformality in QCD" — Brown-Rho scaling, trace anomaly/dilaton (CEA Saclay)
A. Ravlić, T. Nikšić, Y.F. Niu, P. Ring, N. Paar, Phys. Rev. C 110, 024323 (2024) — Axially deformed pnRQRPA based on point-coupling interactions: first GT, SD, Fermi in deformed nuclei (Univ. of Zagreb/IMP Lanzhou/RIKEN)
J.C. Zhang, B.H. Sun, I. Tanihata et al., Phys. Rev. X 15, 031004 (2025) — Charge pickup reaction cross sections for neutron-rich p-shell isotopes at 900A MeV: 50× exponential growth, new global formula (Beihang Univ./GSI/Osaka Univ.)
NUMEN Collaboration (C. Agodi, F. Cappuzzello et al.), Results in Physics 67, 108563 (2024) — First ¹¹⁶Sn(¹⁸O,¹⁸Ne)¹¹⁶Cd DCE cross section at forward angles: 34 nb (INFN-LNS Catania)
A. Spatafora, F. Cappuzzello et al., Phys. Rev. C 110, 044603 (2024) — ¹²C(¹⁸O,¹⁸F)¹²B single charge exchange at 275 MeV: meson exchange vs sequential transfer competition (INFN-LNS/Univ. Catania)
A. Spatafora, F. Cappuzzello et al., EPJ Web Conf. 304, 01005 (2024) — ¹⁸O+⁴⁸Ti SCE multi-channel analysis (INFN-LNS/MAGNEX)
NUMEN Collaboration (I. Ciraldo, C. Agodi et al.), Nuovo Cimento C 47, 6 (2024) — ²⁰Ne+¹³⁰Te multi-channel campaign: DCE, SCE, transfer for 0νββ NME extraction
Y.K. Wang, P.W. Zhao, J. Meng, Phys. Lett. B 855, 138796 (2024) — Strong linear correlation between 0νββ and double Gamow-Teller NMEs via relativistic CI-DFT (Peking Univ.)
H. Ejiri, AIP Conf. Proc. 3143, 020005 (2025) — Nuclear matrix elements for DBD and spin-dipole nuclear responses: pnQRPA with experimental gph calibration (RCNP Osaka)
H. Sagawa, T. Naito, X. Roca-Maza, T. Hatsuda, Phys. Rev. C 109, L011302 (2024) — CSB energy density functional from QCD sum rules (RIKEN/Univ. of Milan/Univ. of Tokyo)
H. Sagawa, C. Ma, T. Inakura, Prog. Theor. Exp. Phys. 2024, 063D02 (2024) — HF-RPA with tensor interactions: M1, GT, and spin-dipole excitations
N. Song, Z.H. Li et al., Chin. Phys. C 49, 064102 (2025) — Semi-empirical parametrization of ν-nucleus cross sections for supernova networks (CIAE Beijing)
J. Kasurinen, P. Pirinen, J. Suhonen, Phys. Rev. C 111, 014601 (2025) — CC neutrino-nucleus scattering on Xe isotopes with meson-exchange 2BC (Univ. of Jyväskylä)
D. Vale, Y.F. Niu, N. Paar, arXiv:2012.11977 — Spherical pnRQRPA with DD-PCX: GT resonance in ²⁰⁸Pb and T=0 pairing (Univ. of Zagreb/IMP Lanzhou)
A. Tamii, N. Kobayashi, AAPPS Bulletin 34, 7 (2024) — Review of Grand Raiden: dipole responses, PDR, M1, (³He,t) charge exchange at RCNP
K. Yoshida, arXiv:2103.16119 (updated 2025) — Isovector spin susceptibility: isotopic evolution of collectivity in spin-flip and charge-exchange channels (Kyoto Univ.)
D. Carbone et al. (NUMEN), Conf. presentation Sept. 2025 — "Recent results from the NUMEN project": DCE cross sections for 0νββ candidates (INFN-LNS)
C. Eke, I. Ciraldo et al., JINST 20, P02015 (2025) — Multi-layer THGEM prototype for NUMEN focal plane detector
GSI/R3B Collaboration — Charge exchange reaction measurements with relativistic RI beams at R³B (2024-2025)
P. Sarriguren et al., ongoing 2024-2025 — (³He,t) charge-exchange: GT excitations in exotic nuclei and isospin symmetry breaking (CSIC Madrid)
F. Recchia et al., Nature Commun. (2025) — Isospin-symmetric island of inversion in Mo isotopes at N≈Z≈40 (Univ. of Padova/GRETINA/NSCl)
C. J. Lin, Z. H. Li et al., Phys. Rev. C 109, 064610 (2024) — Semimicroscopic description of charge-exchange giant spin-dipole resonances in closed-shell parent nuclei (CIAE Beijing)
S. Beyer, F.M. Nunes, APS DNP Fall Meeting L17.1 (2024) — Bayesian Lane-consistent optical potential from (n,n), (p,p), (p,n) IAS: symmetry energy constraints (FRIB/MSU)
Smith, Hebborn, Nunes, Zegers, Phys. Rev. C 110, 034602 (2024) — Uncertainty quantification in (p,n) IAS reactions: ⁴⁸Ca at 25-160 MeV, 60-100% parametric uncertainties (FRIB/MSU)
D. Frekers, Prog. Part. Nucl. Phys. 57, 217 (2006) — (d,²He) charge-exchange probe: GT strength measurements (historical foundation, Univ. of Münster)