跳到正文
返回核反应
至真探索 · 核反应

少体物理与Faddeev方程

2026-08-25

1960年，列宁格勒的Ludvig Faddeev写下一组耦合积分方程，把量子三体问题从一团乱麻变成数学上可解的东西。两个粒子怎么散射，Lippmann-Schwinger方程早就解决；可一加到三个粒子，方程里会冒出"非连接图"——某些项只让其中两个粒子相互作用、第三个晾在一边不参与，算出来的东西发散。Faddeev的招是把三体波函数拆成三个分量，每个分量对应一对粒子加一个旁观者，三个分量耦合在一起却各自满足正确的边界条件。非连接图没了，方程能算。这套方程从此成为核少体物理的标准工具，名字就叫Faddeev方程。Faddeev本人在圣彼得堡国立大学和Steklov数学所圣彼得堡分所工作了一辈子，2017年去世，除了三体方程，还和Victor Popov一起给规范场路径积分定了Faddeev-Popov鬼场的规矩。

六十多年过去，这套方程还在被刷新。2024年这一年，少体圈同时冒出三件值得记的事：相对论Faddeev出了3D无t矩阵的新写法，氘核打核的三体散射换上微观光学势，氚和氦-3的结合能被用来给三核力的电荷依赖成分加新约束。把这三件放一起看，能看清Faddeev方程这条线走到哪了。

先说清楚方程长什么样。Faddeev把三体波函数Ψ拆成Ψ=Ψ₁+Ψ₂+Ψ₃，每个Ψᵢ满足一个积分方程，里头出现的不是裸的两体相互作用V，而是两体散射的t矩阵——也就是两体散射问题先解出来、再塞进三体方程。这是经典Faddeev的形式：t矩阵依赖型。1967年Alt、Grassberger和Sandhas（AGS）把方程重写成转移算符的形式，更光滑，成为散射计算的事实标准。Coulomb力是个老大难：它是长程的，Faddeev原始框架只处理短程。Merkuriev在1971到1980年间想出办法，把每对粒子的Coulomb切成短程和长程两部分，长程部分单独处理，恢复各通道的渐近解耦——这就是Faddeev-Merkuriev方程。四体呢？Yakubovsky 1967年把Faddeev的思路推到任意多体，给出Faddeev-Yakubovsky方程，四体是两套耦合分量（3+1链和2+2链）。

相对论版是后话。三体要做Poincaré不变，就要处理boosted两体t矩阵——把静止系的两体散射"抬升"到运动系，伴随Wigner转动。Witala、Golak、Skibiński（克拉科夫Jagiellonian大学）、Glöckle（波鸿鲁尔大学）、Kamada（九州工业大学）、Polyzou（爱荷华大学）在2011年把相对论三核Faddeev方程扩到包含三核力（arXiv:1101.4053），这套一直是标杆。麻烦在于boosted t矩阵算起来又慢又重。2024年1月，PTEP刊出Mohammadzadeh、Radin和Hadizadeh的一篇（PTEP 2024, 013D03，arXiv:2312.05132），给出一个"t矩阵无关"的相对论3D Faddeev方程：直接把boosted两体相互作用塞进方程，跳过两体t矩阵这一整步。思路继承自Golak等2013年的非相对论3D形式（Few-Body Syst. 54, 2427）。他们用Malfliet-Tjon和Yamaguchi两种自旋无关势验证，无t矩阵版精确复现了传统t矩阵依赖版的3体质量本征值，还省一大笔算力。关键是3D——不做偏波分解（partial wave decomposition），直接在矢量动量空间解，避开角动量代数在高自旋时的指数爆炸。同一团队还把这套无t矩阵思路搬到四体Yakubovsky上（arXiv:2311.08713），思路一脉相承。补一句：这篇PTEP论文由NSF资助（2000029），Hadizadeh是少体3D Faddeev/Yakubovsky的长期玩家；三位作者的精确机构在公开摘要里没单列，从arXiv提交记录和资助号看是北美线的研究组。

第二件是氘核打核的三体散射。2024年Few-Body Systems第65卷第68篇（DOI:10.1007/s00601-024-01937-x，arXiv:2502.17351），Deltuva、Jurčiukonis、Likandrovas和Torres Fernandez四人写的。他们解Faddeev型三体转移算符方程，算氘核在¹²C、¹⁶O、²⁴Mg上的弹性散射和破裂，给出微分截面和分析本领，和已有实验数据吻合得不错。新意在于：把光学势从唯象的换成微观全局光学势——也就是从核子-核散射从头算出来的势，而不是靠拟合参数凑的。结果跟唯象势一比，某些运动学区有系统性偏差，说明微观势还有地方要改。Deltuva是Faddeev-Merkuriev+Coulomb这条路的核心人物，早年从里斯本大学核物理中心起家，现在在维尔纽斯大学理论物理与天文研究所，四作者全在这个所。这条路的本事是把长程Coulomb塞进动量空间Faddeev方程：用screening加renormalization的招数（Deltuva 2005起），把Coulomb软化成短程近似、算完再重整。这套方法不只是三核，也用到氘核打⁴He到⁵⁸Ni、质子打弱束缚的¹¹Be/¹⁵C/¹⁷O，相当于拿Faddeev方程给传统的CDCC（连续离散耦合道）方法当标尺，量一量近似到底准不准。

第三件是氚和氦-3。这俩是三核力的试金石。故事老：拿高精度两核力势——AV18、CD-Bonn、Nijmegen——解Faddeev方程算³H和³He结合能，结果比实验低0.5到1 MeV，欠束缚。补一个2π交换型三核力就回来了。三核力的源流要追到1957年Fujita和Miyazawa，后来Tucson-Melbourne、Urbana IX都是这条路。真正系统化是手征有效场论（χEFT）：3N力在次次领头阶（NNLO）第一次出现，带两个待定参数c_D和c_E。2002年Epelbaum、Nogga、Glöckle、Kamada、Meißner、Witała那篇经典（PRC 66, 064001，arXiv:nucl-th/0208023）把c_D、c_E用氚结合能和nd双重态散射长度定标，第一次完整分析了含3NF的nd散射。这六个人的机构是手征核力那个黄金组合：Epelbaum在Jülich加波鸿、Nogga那时在亚利桑那、Glöckle在波鸿、Kamada在九州工大、Meißner在Jülich加格拉茨、Witała在克拉科夫。往后N3LO的长程3NF也塞进了氚的Faddeev计算（Skibiński等九人，arXiv:1107.5163）。2024年12月，Gapchenko、Goleva、Egorov在Phys. Atom. Nucl. 87, 682（DOI:10.1134/S1063778824700698，俄文Yadernaya Fizika的英译）上给了一个新算法：直接积分法解齐次Faddeev方程求³H/³He结合能，同样不做偏波分解，两体t矩阵由local和nonlocal势生成。他们估了Coulomb排斥和3N力对结合能的影响，得出一个有意思的结论：³H和³He缺失束缚之间那点143 keV的skew，必须靠电荷依赖的NN相互作用来解释——给3NF的电荷依赖成分下了新约束。三位作者的精确机构在公开摘要里没标，期刊属性是俄文核物理，团队属俄罗斯。同年11月，Bovermann、Epelbaum、Krebs和Lee（arXiv:2411.19613）在格点上做了非微扰N3LO手征三核子模拟，c_D由氚β半衰期定、c_E由氚结合能定，预言的氦-3（helion）能量落在实验误差内，氚和氦-3的电荷半径比实验低几个百分点——指向两体电荷密度交换流算符还缺。Epelbaum和Krebs在波鸿，Lee是格点核物理的能手。

四体这边，Yakubovsky的方程在2001年由Nogga、Kamada、Glöckle（NPA 689, 357，arXiv:nucl-th/0010005）用真实NN加各种3N力系统解过3N和4N bound态。结论是：所有3N力模型都能消掉氚和α粒子的欠束缚，理论值和实验吻合得相当好，给四核力留下的余地很小。3D Yakubovsky含3BF的形式是Hadizadeh等人2007年做的（arXiv:0704.2056），把方程写成矢量Jacobi动量的耦合三维积分，六变量，省掉偏波代数。到halo核，Ahmadi Pouya和Rajabi（伊朗Shahrood理工大学物理系，arXiv:1904.12451）把8核子Yakubovsky形式套到⁸He上，用α+n+n+n+n的五体子团模型，约化掉8和7的束缚核子贡献后退化成⁶He的α+n+n三体——halo核这种松散结构，正是少体精确解能直接咬住的对象。

把这些串起来看，Faddeev方程这条线六个轴在同时转。一是从短程到长程：Merkuriev和Deltuva把Coulomb塞进来，质子-氘核散射从此能和氘核-核散射同精度算。二是从非相对论到相对论：Witala/Polyzou那套标杆，再到2024年Mohammadzadeh-Hadizadeh的t矩阵无关版，省算力是实在收益。三是从偏波分解到3D：Golak 2013开了头，Hadizadeh接到四体，Gapchenko 2024接到氚氦结合能——3D不是花哨，是因为高自旋和多体时偏波通道数会爆炸，3D直接拿矢量解更利落。四是从t矩阵依赖到t矩阵无关：传统Faddeev要先解两体t矩阵再塞进去，新写法直接用相互作用，省一步；但要注意，t矩阵无关目前只对bound state可靠，散射还得靠t矩阵。五是从3体到4体到halo子团：Yakubovsky 1967给框架，Nogga-Kamada-Glöckle用真实力算，Ahmadi Pouya-Rajabi套到⁸He。六是3N力的层级演进：Fujita-Miyazawa的2π、Tucson-Melbourne和Urbana IX的唯象、χEFT的c_D/c_E、N3LO长程项、再到2024年格点非微扰N3LO——每一代都用氚和氦-3的结合能当定标锚点，³H/³He一直是少体第一性原理的基准线。

2024这一年特别能说明问题：三条原本各走各的路——相对论形式化、三体散射的微观势、³H/³He对3NF的约束——在同一年各自交了答卷，背后是同一个共识：少体精确解是把核力理论钉死的唯一办法。两核力再精，三体里总会差那么0.5-1 MeV，这缺口正是3N力的签名；3N力的参数又只能靠³H/³He/⁴He这些小系统来定，而不是靠重核拟合。Faddeev方程的价值就在这：它给"从夸克到核"那条第一性原理链提供了一个能精确验算的结点。Faddeev本人2017年走了，他1960年那组方程还在替全人类算着最干净的那几个原子核。

参考文献：

  1. L.D. Faddeev, "Scattering theory for a three-particle system," Sov. Phys. JETP 12, 1014 (1961)；及专著 Mathematical Aspects of the Three-Body Problem in Quantum Scattering Theory (Steklov Inst., 1965).
  2. E.O. Alt, P. Grassberger, W. Sandhas, "Reduction of the three-particle scattering problem to multi-channel two-particle Lippmann-Schwinger equations," Nucl. Phys. B2, 167 (1967)（AGS方程）.
  3. S.P. Merkuriev, "Coulomb scattering theory for three charged particles," Ann. Phys. (N.Y.) 98, 282 (1976)；及早期 1971/1974 系列.
  4. O.A. Yakubovsky, "On the integral equations in the theory of N particles," Sov. J. Nucl. Phys. 5, 937 (1967)（Faddeev-Yakubovsky 四体方程）.
  5. M. Mohammadzadeh, M. Radin, M.R. Hadizadeh, "Relativistic Faddeev 3D Equations for Three-Body Bound States Without Two-Body t-Matrices," Prog. Theor. Exp. Phys. 2024, 013D03 (2024). DOI: 10.1093/ptep/ptad153；arXiv:2312.05132（NSF资助2000029；作者精确机构见原文）.
  6. A. Deltuva, D. Jurčiukonis, D. Likandrovas, J. Torres Fernandez（维尔纽斯大学理论物理与天文研究所）, "Three-Body Calculation Of Deuteron–Nucleus Scattering Using Microscopic Global Optical Potential," Few-Body Syst. 65, 68 (2024). DOI: 10.1007/s00601-024-01937-x；arXiv:2502.17351.
  7. A. Gapchenko, O. Goleva, M. Egorov, "Binding Energies of ³H and ³He Nuclei in Three-Body Faddeev Equations with Direct Integration," Phys. Atom. Nucl. 87, 682–696 (2024). DOI: 10.1134/S1063778824700698（俄Yadernaya Fizika英译；作者属俄罗斯团队，精确机构见原文）.
  8. E. Epelbaum, A. Nogga, W. Glöckle, H. Kamada, U.-G. Meißner, H. Witała, "Three-Nucleon Forces from Chiral Effective Field Theory," Phys. Rev. C 66, 064001 (2002). arXiv:nucl-th/0208023（机构：Jülich/波鸿/亚利桑那/九州工大/格拉茨/克拉科夫Jagiellonian）.
  9. R. Skibiński, J. Golak, K. Topolnicki, H. Witała, E. Epelbaum, W. Glöckle, H. Krebs, A. Nogga, H. Kamada, "The triton with long-range chiral N3LO three-nucleon forces," arXiv:1107.5163（Jagiellonian Kraków + Bochum + Jülich + Arizona + Kyushu IT）.
  10. H. Witała, J. Golak, R. Skibiński, W. Glöckle, H. Kamada, W.N. Polyzou, "Three-nucleon force in relativistic three-nucleon Faddeev calculations," arXiv:1101.4053（Jagiellonian Kraków + Bochum + Kyushu IT + Iowa）.
  11. A. Nogga, H. Kamada, W. Glöckle（Ruhr U. Bochum）, "Solution of the Faddeev-Yakubovsky equations using realistic NN and 3N interaction," Nucl. Phys. A 689, 357 (2001). arXiv:nucl-th/0010005.
  12. J. Golak, K. Topolnicki, R. Skibiński, W. Glöckle, H. Kamada, A. Nogga, "A Three-Dimensional Treatment of the Three-Nucleon Bound State," Few-Body Syst. 54, 2427 (2013). DOI: 10.1007/s00601-012-0472-5；arXiv:1206.3192（3D无偏波分解先声）.
  13. M.R. Hadizadeh et al., "Bound State Calculations of the Three-Dimensional Yakubovsky Equations with the inclusion of Three-Body Forces," arXiv:0704.2056 (2007)；及四体t矩阵无关版 arXiv:2311.08713 (2023).
  14. E. Ahmadi Pouya, A.A. Rajabi（Shahrood University of Technology, Semnan, Iran，物理系）, "Implementation of the 8-Nucleon Yakubovsky Formalism for Halo Nucleus ⁸He," Am. J. Mod. Phys. 8, 53 (2019)；arXiv:1904.12451.
  15. L. Bovermann, E. Epelbaum, H. Krebs, D. Lee, "Non-perturbative three-nucleon simulation using chiral lattice EFT," arXiv:2411.19613 (2024)（Epelbaum/Krebs: Bochum；Lee: 格点核物理）.
  16. A. Deltuva, A.C. Fonseca, R. Lazauskas, "Faddeev equation approach for three-cluster nuclear reactions," arXiv:1201.4979 (2012)（里斯本大学核物理中心 + IPHC Strasbourg；含Coulomb的screening-renormalization方法综述）.