⚛️ 核理论

Pfaffian and minor formulas in number-conserved pairing theory

本文提出了用于计算奇偶系统数守恒配对理论中范数和密度矩阵的闭合形式 Pfaffian 及行列式公式，并将其应用于 $^{70}$ Ge 和 $^{69}$ Ge 的角动量投影 BCS 计算，以解释其低能集体带。

原作者： G. J. Fu

发布于 2026-10-21

📖 1 分钟阅读🧠 深度阅读

原作者： G. J. Fu

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

在原子核的心脏地带，质子和中子并非仅仅静止地堆积在一起，而是进行着一场复杂的、协调的量子力学之舞。为了理解这些微小粒子的行为，科学家们经常依赖于将它们视为成对运动的理论，就像舞者手牵手一样。这种配对是一种基本的相互作用力，有助于将原子核凝聚在一起并决定其形状。然而，该领域面临的一个主要挑战是，最常用的描述这些配对的数学工具往往违反了一项自然的关键规则：它们允许粒子总数发生波动。虽然这种近似法在处理大型系统时表现良好，但在研究像原子核这样微小的有限物体时却会失效，因为在这些物体中，精确计数每一个粒子对于准确性至关重要。长期以来，研究人员一直需要一种既能描述这些配对系统，又能严格保持粒子数量固定的方法，而这项任务需要穿越复杂的计算迷宫，以预测这些原子核如何振动、旋转和改变形状。

上海同济大学的物理学家傅刚（G. J. Fu）开展的一项新研究，为解决这一特定问题提供了一套强大的新数学工具。研究人员开发了精确的闭合形式公式，使科学家能够在严格守恒粒子数的情况下计算原子核的性质。这些公式充当了桥梁，连接了原子核的不同理论状态，并允许研究人员将它们混合在一起，以找到对现实最准确的描述。与以往难以处理破裂配对或奇数粒子情况的方法不同，这种新方法可以无缝衔接奇数和偶数个粒子。它能够处理配对完整以及配对破裂的情况，提供了一种统一的方法来计算“模”（即发现原子核处于某种特定状态的概率）以及描述原子核对各种力做出反应的矩阵元。通过使用这些公式，科学家现在可以应用先进的技术（例如投影原子核的总自旋），在不丢失粒子守恒这一关键细节的前提下，揭示这些微观世界的真实结构。

为了测试这些新公式的力量，该研究将其应用于锗元素的两种特定同位素：锗-70和锗-69。锗-70是一个偶数核，意味着它拥有相等数量的质子和中子，这些粒子都可以配对。相比之下，锗-69有一个奇数的中子，留下了一个未配对的粒子。研究人员使用了一种称为角动量投影数守恒BCS理论的方法，该方法将新公式与一种强制要求原子核具有特定、明确自旋的技术相结合。目标是观察这种方法是否能准确重现已知能量级和原子核的形状，特别是它们的集体能带，即那些像家族一样协同运动的激发态序列。

对于锗-70，计算表明该原子核并不只有一种简单的形状。相反，最低能量状态（即基态带和伽马带）主要源于一种三轴形，即原子核在三个不同方向上被拉伸，呈现出一种不是完美圆形的、略微被挤压的球体形状。相比之下，高能贝塔带则被发现源于一种扁平的轴对称形状。研究表明，基态通过混合这些不同的形状获得了显著的稳定性，而以往的方法可能会忽略这一过程。结果与实验数据吻合得惊人地好，准确预测了能级以及它们之间跃迁的强度。例如，该模型成功描述了原子核在旋转时如何释放能量，区分了不同类型的转动能带，并解释了为什么某些跃迁较强而另一些较弱。

研究还将注意力转向了带有未配对中子的奇数质量核——锗-69。这种情况更为复杂，因为未配对的粒子会以复杂的方式与配对的核心发生相互作用。研究人员构建了一个模型，其中未配对的中子与锗-68核心的两种不同底层形状相耦合。计算表明，锗-69的低能态并非由单一形状描述，而是由这两个核心的混合体构成。随着原子核自旋的增加，其状态的特征会发生变化，从由一种核心形状主导转变为由另一种核心形状主导。这种混合对于理解正宇称态至关重要，因为这些状态需要来自两种底层结构的贡献才能得到准确描述。该模型成功重现了观测到的能级和跃迁强度，证实了新公式可以在不损失精度的情况下处理未配对粒子带来的复杂性。

这项工作的意义不仅限于这两个同位素。通过提供一个适用于奇数和偶数粒子数、以及完整和破裂配对的严谨数学框架，该研究为更精确地模拟广泛的核现象打开了大门。它允许研究人员以一种此前难以实现的细节水平，去探索量子相变和核形状的演化。这些公式是生成坐标法的多功能工具，该技术用于混合不同的核构型以寻找真实的基态。本质上，这项研究精炼了物理学家观察原子核的“透镜”，确保即使在探索物质最微观尺度下最复杂和最具集体性的行为时，也不会违反粒子守恒这一基本规则。研究结果表明，不同形状之间的相互作用以及未配对粒子的存在，是解锁这些量子系统完整行为的关键，为理解宇宙的构建模块提供了更清晰的路径。

技术摘要：配对理论中 Pfaffian 与行列式小项（Minor）公式的研究

问题陈述
配对理论，特别是 BCS（Bardeen-Cooper-Schrieffer）和 HFB（Hartree-Fock-Bogoliubov）方法，是模拟量子多体系统的标准工具。然而，这些方法本质上破坏了 U(1) 规范对称性，导致粒子数不守恒。虽然在处理大系统时，这种违反在正则系综内通常可以忽略不计，但对于原子核等有限尺寸系统而言，准确描述内在性质和量子相变需要严格的粒子数守恒。

现有的恢复粒子数的方法，例如通过对规范角进行数值积分的粒子数投影法，计算量巨大。另一种选择是直接守恒粒子数的方法，例如守恒粒子数对凝聚态，这提供了一个更直接的框架。然而，计算这些态的范数（norm）和矩阵元——特别是在处理破对配置（奇数粒子数）或在生成坐标方法（GCM）中涉及不同正则基组混合时——仍然非常复杂。此外，当底层波函数缺乏良好的角动量时，通常需要进行角动量投影，这要求能够高效地评估重叠（overlap）和哈密顿量矩阵。

方法论
本研究推导了具有精确粒子数之多体算符的范数与矩阵元的闭合解析公式。该形式化方法同时处理了偶数和奇数粒子数系统，并兼容了破对配置（配对数 $\nu > 0$ ）。

该方法的核心包括：

广义破对态： 该形式化方法将标准的 $N$ 对凝聚态 $|c\rangle = (\hat{P}^\dagger)^N |0\rangle$ 扩展为破对态 $|\alpha\rangle = (\hat{P}^\dagger)^N \hat{c}^\dagger_{\alpha_1} \cdots \hat{c}^\dagger_{\alpha_\nu} |0\rangle$ ，其中 $\nu$ 代表未配对粒子的数量。
代数推导： 利用广义 Kronecker delta 以及正则基组间酉变换的性质，将重叠 $\langle \alpha | \hat{R} | \beta' \rangle$ （其中 $\hat{R}$ 是旋转算符）表示为对配对指标组合的求和。
Pfaffian 与行列式小项（Minor）表述： 推导出的表达式以两种等价形式呈现：
- 基于行列式小项（Minor-based）： 范数被表示为由变换矩阵 $X = U^* W'^\top$ 构成的子矩阵的行列式（minors）之和。
- 基于 Pfaffian（Pfaffian-based）： 通过将库珀对产生算符用矩阵 $P$ 表示，重叠被表示为涉及 $P$ 和 $Q$ （旋转后的对矩阵）的子矩阵的 Pfaffian，并乘以单粒子变换矩阵的行列式。
算符矩阵元： 该形式化方法被推广到 $m$ 体算符，提供了不同破对态及不同正则基组之间矩阵元的紧凑表达式。这些公式适用于 GCM 中的配置混合。

主要贡献

闭合形式公式： 本文提供了对偶数和奇数粒子数以及具有破对态有效的范数与矩阵元的显式闭合形式表达式。
统一框架： 推导出的公式既适用于对凝聚态，也适用于含有未配对粒子的配置，从而促进了在同一理论框架下对奇质量原子核的处理。
计算效率： 通过将重叠和矩阵元表示为 Pfaffian 和行列式小项，这项工作为对称性恢复技术（特别是角动量投影）提供了一条计算可行路径，而无需依赖于对规范角的数值积分。

结果与应用
作者将推导的形式化方法应用于角动量投影的守恒粒子数 BCS（PNBCS）理论，利用 JUN45 有效相互作用在 $1p1/2, 1p3/2, 0f5/2, 0g9/2$ 壳模型空间内研究了 $^{70}\text{Ge}$ 和 $^{69}\text{Ge}$ 两个核素。

$^{70}\text{Ge}$ （偶-偶核）：
- 计算混合了两个内在 Hartree-Fock (HF) 配置：一个三轴形解 $(\beta, \gamma) = (0.25, 36^\circ)$ 和一个轴对称扁平解 $(\beta, \gamma) = (0.13, 60^\circ)$ 。
- 结果重现了低能集体带（基态带、 $\gamma$ 带和 $\beta$ 带）。研究发现，基态带和 $\gamma$ 带主要源于三轴形配置，而 $\beta$ 带则与扁平配置相关。
- 计算的 $B(E2)$ 跃迁强度和激发能与实验数据以及全配置相互作用（FCI）壳模型计算结果吻合良好。
- 研究强调了角动量投影和配置混合带来的能量增益，解释了为什么尽管在投影前扁平 HF 解能量更低，但基态 $\gamma$ 带的 $0^+$ 态能量却低于 $\beta$ 带的 $0^+$ 态。
$^{69}\text{Ge}$ （奇质量核）：
- 多体基组通过将单个未配对中子耦合到 $^{68}\text{Ge}$ 的两个对凝聚态（代表形状共存）来构建。
- PNBCS 计算成功描述了六个低能集体带（ $5/2^-, 3/2^-, 1/2^-, 11/2^+$ 以及两个 $9/2^+$ 带）。
- 对核心重叠（ $W_i$ ）的分析揭示了波函数结构的自旋依赖演化。对于负宇称带，其特征表现为：在低自旋时由其中一个内在核心主导，而在高自旋时转向由另一个核心主导。
- 正宇称带头显示出两个核心配置之间显著的混合，表明单一的内在结构不足以描述这些态。

意义
论文声称，所提出的形式化方法为在守恒粒子数配对理论中进行对称性恢复计算提供了一个严谨且高效的工具。通过提供涉及破对和不同正则基组的范数与矩阵元的解析公式，该工作促进了生成坐标方法和角动量投影在处理包括奇质量系统在内的计算中的应用。对 $^{70}\text{Ge}$ 和 $^{69}\text{Ge}$ 的成功应用证明了该方法在描述中质量原子核的集体结构和形状共存方面的能力，为标准对称性破缺方法或计算量巨大的 FCI 方法提供了可行的替代方案。

类似论文