Learning objectives / 学习目标
- 从 \(B/A\) 饱和、散射与镜像核系统学说明核力的短程和近似电荷无关性。
- 用自旋投影算符、张量算符和自旋—轨道项区分不同核子通道。
- 建立氘核质心约化问题，计算其束缚尾尺度并解释 \(S\)–\(D\) 混合证据。
- 比较方阱、Yukawa、现实两核子势、手征 EFT 与光学势的用途和边界。
Prerequisites / 先修知识
- 有限深势阱与束缚态。
- 角动量耦合 与 Pauli 矩阵。
- 结合能、分离能与反应阈值。
Scope and conventions / 范围与约定
- “核力”在本课指低能核子自由度之间的有效相互作用；底层 QCD 并未被一个简单位置势精确替代。
- 讨论 \(pp\)、\(nn\)、\(np\) 时，近似电荷无关性必须先剔除主要电磁作用，并允许小的电荷对称破缺。
- 两体相对坐标取 \(\mathbf r=\mathbf r_p-\mathbf r_n\)，约化质量 \(\mu=m_pm_n/(m_p+m_n)\)，质心静止。
- 氘核基态 \(J^\pi=1^+\)、同位旋 \(T=0\)；“主要为 \(^3S_1\)”不等于纯 \(S\) 波。
- 势函数不是直接可观测量；不同正则化、幺正变换或 EFT 截断可给相同相移与束缚能。
Notation and modeling contract / 符号与建模契约
| Symbol | Meaning | SI unit |
|---|---|---|
| \(r\) | nucleon separation | m |
| \(\mu\) | two-body reduced mass | kg |
| \(B_d\) | deuteron binding energy | J |
| \(S\) | dimensionless total two-nucleon spin quantum number; physical spin operator is \(\hbar\mathbf S\) | 1 |
| \(T\) | dimensionless total isospin quantum number (internal label, not spatial angular momentum) | 1 |
| \(S_{12}\) | tensor operator | 1 |
| \(V_C,V_T,V_{LS}\) | central, tensor and spin-orbit potential terms | J |
| \(\delta_\ell(E)\) | scattering phase shift | rad |
1. Short range, saturation, and approximate charge independence / 短程、饱和与近似电荷无关性
Knowledge ID: P10-U02-L01-K01
核子—核子强相互作用的有效吸引主要作用在数 fm 内，更短距离表现强排斥或强短程关联。其有限范围使每个核子只与有限近邻强耦合，从而总结合能主要随 \(A\) 而非 \(A^2\) 增长，这就是饱和的经验含义。饱和并非“势在某距离突然为零”，也不等于核物质不可压缩；它是多体能量与密度的综合性质，包含 Pauli 原理和两体/三体相互作用。
介子交换图像给出长程尾的尺度直觉：质量为 \(m_\pi\) 的最轻交换自由度产生约
\[V(r)\propto\frac{e^{-m_\pi cr/\hbar}}{r}, \qquad R_\pi\sim\frac{\hbar}{m_\pi c}\approx1.4\ \mathrm{fm}.\]
这不意味着整个核势就是单一 Yukawa 标量势；真实一 pion 交换含自旋、张量和同位旋结构，中短程还需更多自由度或接触项。
剔除电磁作用后，低能 \(pp\)、\(nn\)、\(np\) 强相互作用显示近似电荷独立性，可用同位旋形式组织。但“近似”很关键：夸克质量差、电磁修正及电荷对称破缺仍可测。Argonne \(v_{18}\) 就包含主导电荷无关部分、额外电荷依赖/电荷对称破缺项和完整电磁相互作用。
图 1：核子势的定性径向结构；深度和形状依通道与模型而变。
图 2：有限近邻连接给出 \(B\propto A\) 的直觉，而非全连接 \(A(A-1)/2\)。
图 3：去除主要电磁作用后才比较强相互作用通道。
2. Spin and tensor dependence / 自旋与张量依赖
Knowledge ID: P10-U02-L01-K02
两个自旋 \(1/2\) 核子可耦合成 singlet \(S=0\) 或 triplet \(S=1\)。因
\[\boldsymbol\sigma_1\cdot\boldsymbol\sigma_2= 2S(S+1)-3,\]
该算符在 triplet 中本征值为 \(+1\)，在 singlet 中为 \(-3\)；含 \(V_\sigma(r)\boldsymbol\sigma_1\cdot\boldsymbol\sigma_2\) 的势自然区分两通道。Pauli 反对称性还把空间、自旋和同位旋量子数联系起来，不能独立任选。
张量算符写为
\[S_{12}=3(\boldsymbol\sigma_1\cdot\hat{\mathbf r}) (\boldsymbol\sigma_2\cdot\hat{\mathbf r}) -\boldsymbol\sigma_1\cdot\boldsymbol\sigma_2.\]
它依赖自旋相对核子连线的方向，并在给定 \(J\) 下耦合相差 \(\Delta\ell=2\) 的分波。氘核 \(J=1\) 因而可混合 \(^3S_1\) 与 \(^3D_1\)。非零电四极矩是 \(D\) 波/张量关联的重要证据，不过高精度计算还需核子有限尺寸和交换流修正。
现实相互作用还含自旋—轨道 \(\mathbf L\cdot\mathbf S\)、同位旋 \(\boldsymbol\tau_1\cdot\boldsymbol\tau_2\)、动量依赖等算符。把“核力”只画成一条 \(V(r)\) 曲线会丢失这些通道信息。
图 4：左图只画 triplet 的 \(m_S=+1\) 成员；真正的 singlet 是反对称叠加态，两通道的 \(\boldsymbol\sigma_1\cdot\boldsymbol\sigma_2\) 本征值不同。
图 5：张量相互作用不是纯中心力；图仅表示角依赖。
3. Deuteron as the simplest bound nucleus / 氘核作为最简单束缚核
Knowledge ID: P10-U02-L01-K03
氘核由一个质子与一个中子组成，是唯一稳定束缚的两核子系统；双质子与双中子没有稳定束缚基态。其结合能约 \(B_d=2.2246\ \mathrm{MeV}\)，远小于典型核势几十 MeV 的深度，因此是弱束缚系统，波函数长尾延伸到作用区之外。MIT 22.101 的方阱示例强调：单独知道 \(B_d\) 只能约束势深与范围的组合，必须结合散射相移等数据才能进一步确定。
质心分离后相对运动满足
\[\left[-\frac{\hbar^2}{2\mu}\nabla_r^2+V\right]\psi=-B_d\psi.\]
作用区外 \(V\to0\)，\(S\) 波径向尾近似 \(u(r)\propto e^{-\kappa r}\)，
\[\kappa=\frac{\sqrt{2\mu B_d}}{\hbar}.\]
实际基态是 \(^3S_1\) 主成分加小的 \(^3D_1\) 成分，\(J^\pi=1^+\)。氘核非零四极矩、磁矩偏离简单质子加中子之和、以及散射数据共同约束张量力。它没有束缚激发态；把方阱的额外数学能级当成真实氘核态会违背数据。
图 6：氘核基态位于 \(p+n\) 阈值下约 2.2246 MeV。
图 7：弱束缚使氘核波函数显著伸出强相互作用区。
图 8：张量项产生 \(S\)–\(D\) 混合并支持非零四极矩。
4. Effective nuclear potentials / 有效核势
Knowledge ID: P10-U02-L01-K04
有效势按用途分层。有限方阱可训练边界匹配并估计氘核尺度；Yukawa 形式抓住最轻介子决定的长程；Argonne \(v_{18}\) 用 18 类算符项拟合大量 \(pp\)、\(np\) 散射及氘核性质；手征 EFT 按低动量展开并同时生成两体、三体力及误差阶次；核反应光学势则通常为复数、能量依赖的平均场。本课采用时间因子 \(e^{-iEt/\hbar}\)，并写 \(U(r,E)=V(r,E)-iW(r,E)\)、\(W\ge0\)，此时负的虚部表示被消去反应通道造成的通量损失；换用其他符号约定必须同时声明 \(\operatorname{Im}U\) 的符号。
IAEA RIPL 光学模型库明确包含全局和核素特定参数化、球形/振动/转动/色散模型，并提醒远离稳定线时纯现象学外推可能失效。光学势用于散射与传输系数，不应拿来作为束缚氘核的厄米两体 Hamiltonian。类似地，壳层 Woods–Saxon 平均场描述一个核子在核芯中的有效运动，也不是自由空间 \(NN\) 势。
一个常见算符展开示意为
\[V=V_C(r)+V_\sigma(r)\boldsymbol\sigma_1\cdot\boldsymbol\sigma_2 +V_T(r)S_{12}+V_{LS}(r)\mathbf L\cdot\mathbf S+\cdots.\]
势的短程形状会随分辨率和正则化改变；真正要比较的是相移、束缚能、矩和反应截面等可观测量，并在同一 Hamiltonian 中一致处理两体力、三体力与算符。
图 9：方阱与 Yukawa 形式突出不同物理，不应互称精确势。
图 10：不同有效描述通向不同观测任务，不能跨层级无条件混用。
Worked Examples / 完整例题
Worked Example 1 — Pion-range estimate / pion 程估算
系统与参考系： 两个低能核子在质心系，相对坐标沿 \(+r\)。
模型： 只用最轻 pion 的质量估计 Yukawa 指数衰减长度；不计算完整自旋—同位旋势。
已知： 取 \(m_\pi c^2=139.6\ \mathrm{MeV}\)，\(\hbar c=197.32698\ \mathrm{MeV\,fm}\)。
\[R_\pi=\frac{\hbar c}{m_\pi c^2} =\frac{197.32698}{139.6}\ \mathrm{fm}=1.413\ \mathrm{fm}.\]
\[R_\pi=1.413\times10^{-15}\ \mathrm m.\]
检查： 能量单位约掉，留下长度；\(m_\pi\to\infty\) 时作用程趋零。结果只估计长程尾，不能推出核心深度或三体力。
Worked Example 2 — Saturation scaling test / 饱和标度检验
系统与参考系： 比较静止 \(A=16\) 与 \(A=64\) 的假想球形核，质心系；各核原点取在质心，本例只比较标量标度而不依赖坐标方向。
模型 A： 每对核子贡献相同 \(-\epsilon\)；模型 B： 每个核子只有 \(z\) 个有效近邻，每对不重复计数。取 \(\epsilon=1.00\ \mathrm{MeV}\)、\(z=8\)。
全连接模型给
\[\frac{B}{A}=\frac{\epsilon(A-1)}2.\]
故 \(A=16\) 时 \(B/A=7.50\ \mathrm{MeV}\)，\(A=64\) 时 \(31.5\ \mathrm{MeV}\)。近邻模型给
\[\frac{B}{A}\approx\frac{z\epsilon}{2}=4.00\ \mathrm{MeV},\]
与 \(A\) 无关。换成 SI，每核子的三个结果依次为 \(1.20\times10^{-12}\ \mathrm J\)、\(5.05\times10^{-12}\ \mathrm J\) 与 \(6.41\times10^{-13}\ \mathrm J\)。
检查： 全连接模型随 \(A\) 线性发散，与实测 \(B/A\) 数 MeV 且近似饱和冲突；近邻模型只是标度直觉，未包含表面与 Pauli 动能。
Worked Example 3 — Deuteron tail length / 氘核尾长度
系统与参考系： 静止氘核，质心系；\(\mathbf r=\mathbf r_p-\mathbf r_n\)。
模型： 作用区外自由 \(S\) 波指数尾；忽略小 \(D\) 波。
已知： \(B_d=2.2246\ \mathrm{MeV}\)，\(\mu c^2\approx469.46\ \mathrm{MeV}\)。
\[\kappa=\frac{\sqrt{2\mu c^2B_d}}{\hbar c} =\frac{\sqrt{2(469.46)(2.2246)}}{197.32698} =0.2316\ \mathrm{fm^{-1}}.\]
\[\kappa^{-1}=4.318\ \mathrm{fm}=4.318\times10^{-15}\ \mathrm m.\]
\(B_d=3.564\times10^{-13}\ \mathrm J\)。
检查： \(\kappa\) 为逆长度；\(B_d\to0^+\) 时尾长发散，符合弱束缚极限。该长度不是实验 rms 电荷半径，也未含张量混合。
Worked Example 4 — Square-well depth constrained by \(B_d\) / 方阱深度
系统与参考系： 氘核相对运动，质心系，球对称 \(S\) 波。
模型： \(V=-V_0\) for \(r<r_0\)、\(V=0\) outside，取 \(r_0=2.00\ \mathrm{fm}\)；\(B_d=2.2246\ \mathrm{MeV}\)、\(\mu c^2=469.46\ \mathrm{MeV}\)。
坐标： \(r\ge0\)。
边界条件为 \(k\cot(kr_0)=-\kappa\)。由例 3 得 \(\kappa r_0=0.4632\)；最低根数值约 \(kr_0=1.820\)，故 \(k=0.910\ \mathrm{fm^{-1}}\)。势深满足
\[V_0=B_d+\frac{(\hbar c)^2k^2}{2\mu c^2}.\]
\[V_0=2.2246+\frac{(197.32698)^2(0.910)^2}{2(469.46)} \approx36.6\ \mathrm{MeV}=5.86\times10^{-12}\ \mathrm J.\]
检查： \(V_0>B_d\)，且约几十 MeV，符合课堂方阱量级；改变 \(r_0\) 会得到不同 \(V_0\)，说明单一束缚能不能唯一反演势。精确氘核还需散射相移、\(D\) 波和交换流。
Lesson synthesis / 全课综合
短程与饱和解释核结合能的 \(A\) 标度，近似同位旋对称组织 \(pp/nn/np\) 通道，自旋与张量算符解释通道差异和氘核 \(S\)–\(D\) 混合。氘核把束缚能、长尾、磁矩、四极矩与散射连接起来，但仍只约束而不唯一决定势。有效势必须连同分辨率、自由度、截断和拟合数据一起陈述。
Formula-and-decision table / 公式与决策表
| Task | Relation | Use when | Boundary |
|---|---|---|---|
| range estimate | \(R\sim\hbar/(mc)\) | lightest exchanged degree sets tail | not full potential |
| spin channel | \(\boldsymbol\sigma_1\cdot\boldsymbol\sigma_2=2S(S+1)-3\) | two spin-\(1/2\) nucleons | tensor and orbital terms may also act |
| tensor structure | \(S_{12}=3(\boldsymbol\sigma_1\cdot\hat r)(\boldsymbol\sigma_2\cdot\hat r)-\boldsymbol\sigma_1\cdot\boldsymbol\sigma_2\) | anisotropic spin coupling | matrix elements depend on channel |
| deuteron tail | \(\kappa=\sqrt{2\mu B_d}/\hbar\) | outside interaction region | ignores Coulomb and \(D\) component |
| square-well match | \(k\cot(kr_0)=-\kappa\) | pedagogical \(S\)-wave well | not a unique realistic force |
| optical potential | \(U=V-iW\), \(W\ge0\), with \(e^{-iEt/\hbar}\) | elastic/reaction scattering after eliminating channels | negative imaginary part absorbs flux |
Model boundaries / 模型边界
- 近似电荷无关性须在去除主要电磁作用后陈述；AV18 明含小的破缺项。
- 饱和是多体性质，不能由“一条两体势曲线”单独证明。
- 单 pion 交换只控制长程结构，不是整个核势。
- 氘核不是纯 \(^3S_1\)；非零四极矩要求 \(^3D_1\) 混合及精细修正。
- 方阱参数依选定作用程；一个束缚能不能唯一确定势。
- Argonne、手征 EFT、平均场和光学势服务不同尺度；势本身不是可观测量。
- 忽略三核子力可能无法同时描述轻核谱与核物质饱和。
Misconception clinic / 误区诊所
- 核力短程不等于核内势能为零之外有一条硬边界。
- 中子不带净电荷不等于 \(nn\) 与 \(pp\) 强作用自动完全相同。
- “自旋相关”不是经典小球转向，而是量子通道的算符依赖。
- 张量力不是普通中心力；它耦合自旋方向与相对坐标。
- 氘核只含两核子，不等于可以忽略 \(D\) 波、交换流与有限核子尺寸。
- 一个拟合优良的势不是唯一真实势；可观测量等价不要求短程形状相同。
Practice and retrieval / 练习与提取
- [K01] 为什么 \(B/A\) 近似常数反对全对全、等强度的核子吸引？
答案 1
全对全配对数约 \(A^2/2\)，会给 \(B/A\propto A\)；实测近似饱和支持有限范围、Pauli 与近邻相关的共同作用。
- [K01] “电荷无关”比较前为何要去除 Coulomb？
答案 2
Coulomb 只作用于带电质子，会使 \(pp\) 与 \(nn/np\) 观测显著不同；不扣除它就无法测试强作用的同位旋近似。
- [K01] 若交换粒子质量加倍，Yukawa 程怎样变？
答案 3
\(R\sim\hbar/(mc)\)，质量加倍使作用程减半。
- [K02] 求 \(S=0\) 与 \(S=1\) 时 \(\boldsymbol\sigma_1\cdot\boldsymbol\sigma_2\) 的本征值。
答案 4
代入 \(2S(S+1)-3\)：singlet 为 \(-3\)，triplet 为 \(+1\)。
- [K02] 张量项为什么能混合 \(^3S_1\) 和 \(^3D_1\)？
答案 5
\(S_{12}\) 是二阶角张量，在保持 \(J=1,S=1\) 时允许 \(\Delta\ell=2\)，故耦合 \(\ell=0\) 与 \(\ell=2\)。
- [K03] 氘核尾长为何明显大于 pion 程？
答案 6
氘核只有 2.2246 MeV 的弱束缚，\(\kappa^{-1}\propto B_d^{-1/2}\) 约 4.3 fm；作用消失后波函数仍指数延伸。
- [K03] 非零氘核四极矩排除了哪种最简单模型？
答案 7
排除完全纯净、球对称的 \(^3S_1\) 点核模型；需 \(D\) 波/张量关联并在精度上加入其他修正。
- [K04] 为什么光学势可以有虚部？
答案 8
被消去的非弹性与反应通道从弹性通道带走通量，虚部在有效单通道方程中表示这种损失。
- [K04] 方阱拟合 \(B_d\) 后是否已确定真实 \(NN\) 势？
答案 9
没有；不同深度—范围组合可给同一束缚能，还需相移、矩和其他通道数据。
- [K04] 势函数为什么不是直接可观测量？
答案 10
它依自由度、正则化与幺正表示；不同势可产生相同散射矩阵和束缚态可观测量。
Exit ticket / 离堂检验
- 用三个独立观测说明核相互作用不能只是自旋无关中心势。
- 计算给定 \(B_d\) 的 \(\kappa^{-1}\)，并解释为何它不是“核力作用程”。
- 给方阱、AV18、手征 EFT 与光学势各写一个适用问题和一个禁用情境。
References and source notes / 参考资料与来源说明
以下直接来源已于 2026-08-30 实际打开并核读：
- MIT OpenCourseWare 22.101, Lecture 6: The Neutron-Proton System: Bound State of the Deuteron. https://ocw.mit.edu/courses/22-101-applied-nuclear-physics-fall-2006/4e1f3290d7b1b41d293a2a3c0f09e2d6_lec06.pdf
- MIT OpenCourseWare 8.701, Lecture 9.4 Nuclear Force. https://ocw.mit.edu/courses/8-701-introduction-to-nuclear-and-particle-physics-fall-2020/86b99997c9e69c54ff16e5f9521918cb_MIT8_701f20_lec9.4.pdf
- Argonne National Laboratory, Argonne v18 Nucleon-Nucleon Potential. https://www.phy.anl.gov/theory/research/av18/
- Argonne National Laboratory, Nuclear Forces and Nuclear Systems. https://www.phy.anl.gov/theory/research/forces.html
- Argonne National Laboratory, Deuteron Properties. https://www.phy.anl.gov/theory/research/deuteron/
- IAEA Nuclear Data Section, Handbook for Calculations of Nuclear Reaction Data: Reference Input Parameter Library-2, optical-model chapter. https://www-nds.iaea.org/RIPL-3/handbook/ripl2-march2006.pdf
- NIST, 2022 CODATA adjustment — Fundamental Physical Constants: Extensive Listing. https://physics.nist.gov/cuu/pdf/all.pdf
Source-to-content map
| Source | Used for | Knowledge ID / Section | Boundary or convention |
|---|---|---|---|
| MIT 22.101 Lecture 6 | deuteron bound state, square-well matching | K03, K04; Examples 3–4 | elementary central well is pedagogical |
| MIT 8.701 Lecture 9.4 | radial, spin-orbit and charge-independence evidence | K01, K02 | lecture-level summary, not a fitted interaction |
| ANL Argonne v18 | operator-rich \(NN\) potential and phase shifts | K01, K02, K04 | one high-precision representation, not unique |
| ANL Nuclear Forces | charge breaking, \(2N/3N\) and QMC context | K01, K04 | many-body consistency required |
| ANL Deuteron Properties | binding, \(D\) state, moments and wave functions | K02, K03 | calculations depend on selected potential and impulse approximation |
| IAEA RIPL-2 | complex/global/local optical potentials | K04 | scattering effective model; extrapolation limits |
| NIST CODATA 2022 | proton, neutron and deuteron constants | K03; Example 3 | constants carry measurement uncertainty |
复习与拓展
- 前课 Binding Energy and Stability 的饱和趋势在本课得到相互作用解释。
- 下一课 Nuclear Models 把复杂相互作用压缩为液滴、壳层和集体自由度。
- 后续核反应课将使用相移、截面和光学势；必须继续声明实验室系或质心系。