来源论文: https://arxiv.org/abs/2606.07415v1 生成时间: Jun 08, 2026 16:36

量子蒙特卡洛与费米子符号问题的几何破局：李-杨零点在有限温度下的演化及双步外推策略深度解析

0. 执行摘要

费米子符号问题（Fermion Sign Problem, FSP）是凝聚态物理、量子化学和量子统计力学领域中公认的、最具挑战性的“圣杯”问题之一。其本质在于费米子波函数的反对称性，导致在路径积分蒙特卡洛（PIMC）等数值计算中，其统计权重会出现正负交替乃至复数的情况，从而使蒙特卡洛采样的方差随体系尺度和逆温度 $\beta = 1/k_B T$ 呈指数级增长。

近年来，学术界尝试通过引入交换宇称 $\xi$ 作为连续调节参数，建立从无符号问题的玻色系统（$\xi = 1$）或可分辨粒子系统（$\xi = 0$）向有严重符号问题的费米系统（$\xi = -1$）的解析延拓方法。然而，这种延拓方法在低温（Low-T）下的收敛性极其脆弱。北京大学李新征教授团队的系列工作通过将经典统计力学中的李-杨零点（Lee-Yang Zeros）理论引入量子多体系统，为此问题提供了全新的几何视角。在本篇论文中，作者使用了一个可精确求解的非相互作用一维环上粒子模型，系统性地追踪了有限温度下配分函数零点（LY零点）在复 $\xi$ 平面上的运动轨迹。

本篇技术博客将对该论文进行深度解构，探讨以下核心发现：

1. 零点的几何根源：在 $T = 0$ K 时，多费米子体系的 LY 零点分布呈现出普适性的代数规律 $z_k = -1/k$。随着温度升高，这些零点偏离实轴，从而重塑了配分函数的解析结构。
2. 符号问题的静电类比：自由能可以完美类比为复平面的二维静电势，而 LY 零点则是带电线电荷。低温下零点逼近实轴 $\xi = -1$ 处，造成该点附近的“静电势”（自由能）急剧发散，这是任何基于 $\xi$ 直接或间接解析延拓在低温下失效的根本原因。
3. 奇偶粒子数限制：由于配分函数的实数正定约束，偶数粒子和奇数粒子的零点轨迹具有截然不同的宇称对称性保护，导致不同的温度行为。
4. 双步外推策略（Two-Step Strategy）：基于配分函数的代数分解 $Z = (\xi - z_1) Z' + \phi(\beta)$，论文提出了一种巧妙的方法——先在无符号问题的安全区（$\xi \in [0, 1]$）进行高温解析延拓得到高精度的费米子性质，再通过对温标 $\beta$ 的拟合（$T$-fitting）将其延拓至低温区。这为解决实际强关联体系的符号问题提供了切实可行的几何破局方案。

1. 核心科学问题，理论基础，技术难点与方法细节

1.1 核心科学问题：费米子符号问题与解析延拓的失效

在量子统计力学中，具有交换对称性的 $N$ 粒子体系在逆温度 $\beta$ 下的配分函数为：

其中 $\hat{\mathcal{P}}_F$ 和 $\hat{\mathcal{P}}_B$ 分别为费米子和玻色子的投影算符。引入交换宇称 $\xi \in [-1, 1]$ 作为参数，可定义广义配分函数 $Z(N, \beta, \xi)$，它满足如下多项式形式：

这里 $\xi = 1$ 对应玻色子，$\xi = -1$ 对应费米子，$\xi = 0$ 对应可分辨粒子（玻尔兹曼统计）。在数值计算中，$\xi \in [0, 1]$ 的区间是不存在费米子符号问题的。因此，研究者希望通过在 $\xi \in [0, 1]$ 进行精确采样，然后将多项式或有理函数解析延拓至 $\xi = -1$。

然而，这一设想在低温（低 $T$）下会彻底崩溃。主要技术难点在于：随着温度降低，费米子符号因子 $\langle s \rangle_B = Z_F / Z_B$ 呈指数级衰减，这在复平面上表现为配分函数的零点逼近物理实轴。当零点极度靠近 $\xi = -1$ 时，解析延拓的收敛半径被严重压缩，导致多项式拟合或温标外推方法的误差放大数个数量级。要攻克这一难点，首先必须弄清这些零点在复平面上的精确位置、演化轨迹以及它们如何定量影响自由能的解析性。

1.2 理论基础：李-杨零点理论与广义配分函数的多项式性质

李-杨零点理论最初由杨振宁和李政道在1952年提出，用于解释经典格点气体和伊辛模型中的相变。他们证明，当格点气体的逸度（fugacity）在复平面上的零点逼近实轴时，体系发生相变。在本工作中，零点理论被推广用于研究交换对称性空间（$\xi$ 空间）。

对于 $N$ 粒子体系，广义配分函数 $Z(\beta, \xi)$ 可以写为不同共轭类（Conjugate Classes）贡献的代数和。设 $Q$ 代表对称群 $S_N$ 的一个共轭类，则：

其中：

• $Z_Q = \int d\mathbf{R} \langle p \mathbf{R} | e^{-\beta \hat{H}} | \mathbf{R} \rangle$ 是对应于共轭类 $Q$ 的路径积分环形聚合物配分函数。对于非相互作用体系，这个多体积分可以完全分解为单粒子路径积分的乘积。
• $F_Q = \frac{N!}{\prod_m m^{\gamma_m} \gamma_m!}$ 是该共轭类中的元素个数。这里 $1^{\gamma_1} 2^{\gamma_2} \cdots N^{\gamma_N}$ 为循环结构（cycle structure），表示有 $\gamma_m$ 个长度为 $m$ 的循环，满足 $\sum_m m \gamma_m = N$。
• $\sigma(Q) = \sum_{m} (m-1) \gamma_m$ 是将置换还原为恒等置换所需的最小对换（pair permutation）次数，它直接决定了 $\xi$ 的幂次。

由于每个循环对应于一个多粒子环形聚合物，非相互作用体系中长度为 $m$ 的循环在逆温度 $\beta$ 下的贡献可以完全简化为单粒子在有效温度 $T/m$ 下的配分函数 $Z_{11}(m \beta)$：

基于上述代数性质，作者给出了 $N = 1$ 到 $5$ 的广义配分函数的精确代数解析式：

这些多项式的高阶零点即为李-杨（LY）零点。对于一维环上的自由粒子（质量为 $m_e$，环长为 $L$），其单粒子能级为 $E_n = \frac{n^2 h^2}{2 m_e L^2}$，对应的单粒子配分函数为：

结合上式即可在任意温度下精确求解出 LY 零点 $z_k(T)$ 的演化轨迹。

1.3 技术难点与静电学类比

为直观理解零点对自由能的影响，我们分析自由能 $F(T, \xi)$ 的解析表达式：

若将配分函数根据其根（零点）$z_i$ 进行因式分解：

则自由能可写为：

这在物理上具有美妙的静电学类比（Electrostatic Analogy）：

• 复平面 $\xi$ 可以看作一个二维静电场空间。
• 每一个李-杨零点 $z_i$ 扮演着一个带负电的无限长均匀线电荷的角色，其在复平面上 $\xi$ 处产生的二维静电势即对应于 $\ln |\xi - z_i|$ 项。
• 自由能 $F(T, \xi)$ 对应于该静电场在物理实轴 $\xi$ 处的总静电势能。

当零点 $z_i$ 位于复平面内、远离实轴时，其实轴处的势能曲线是极其平滑、解析的（如图4中 $T = 1.4$ 或 $T = 0.8$ 的曲线）。然而，在低温极限下，$z_1$ 极度逼近 $\xi = -1$，根据推导，其渐近行为为：

这意味着一个巨大的带负电荷的“电荷”贴在了 $\xi = -1$ 的皮肤上，导致该处的静电势（即自由能）出现对数发散（Divergence），并在其周围产生极其剧烈的“电场强度”（自由能的极值和跃变）。这直接定量解释了为何在低温下，哪怕在 $\xi \in [0, 1]$ 区间拟合精度达到 $10^{-16}$，解析延拓至 $\xi = -1$ 依然会由于电荷阻隔和势能突变而发生毁灭性的发散，失效区域被 LY 零点彻底锁死。

2. 关键 Benchmark 体系、计算所得数据与性能分析

本研究的核心 Benchmark 体系是无相互作用的一维环上多粒子体系，重点解构了 $N = 4$ 和 $N = 5$ 两个具有代表性的奇偶粒子数体系。

2.1 零点轨迹演化数据与奇偶奇点特征

在 $T = 0$ K 极限下，由公式（2）易知，无论物理参数如何，LY 零点均普适地分布在负实轴上：

• 对于 $N = 4$，三个零点位于： $$ z_1 = -1, \quad z_2 = -1/2, \quad z_3 = -1/3 $$
• 对于 $N = 5$，四个零点位于： $$ z_1 = -1, \quad z_2 = -1/2, \quad z_3 = -1/3, \quad z_4 = -1/4 $$

随着温度 $T$ 的升高（单位为 $\frac{h^2}{2 m_e k_B L^2}$），零点在复平面上开始画出美丽的弧形轨迹：

1. $N = 4$ 体系（偶数粒子）：

   • 在低温下，所有三个零点（红、橙、蓝）均停留在负实轴上。
   • 随着 $T$ 增加，两个原本靠近的零点 $z_2 = -1/2$ 和 $z_3 = -1/3$ 逐渐相向运动。在临界温度 $T \approx 0.278$ 处，它们相撞并形成一个二阶重根（tangency point, $\xi \approx -0.45$）。
   • 当 $T > 0.278$ 时，这对零点脱离实轴，分裂为一对共轭复根，并沿着类似于椭圆的弧线向右上方和右下方运动。
   • 与此同时，最靠近费米子物理点（$\xi = -1$）的零点 $z_1$ 在低温下被推向实轴更左侧（例如 $z_1 < -1$），以满足配分函数在两端的正定性。最后，它在更高温度下独自向右移动，但始终保持在实轴上。

2. $N = 5$ 体系（奇数粒子）：

   • 在低温下，四个零点全部在实轴上。
   • 随着 $T$ 增加，它们两两成对相撞。首先是中间的 $z_2 = -1/2$ 与 $z_3 = -1/3$ 相撞，然后是最边缘的 $z_1 = -1$ 与最内侧的 $z_4 = -1/4$ 分别与其他零点结合。在 $T \approx 1.2$ 以上，所有零点全部脱离实轴，变成两对共轭复根。
   • 在复平面上，零点彻底扫清了负实轴，使得 $\xi = -1$ 处的自由能彻底恢复解析性，此时直接的 $\xi$ 解析延拓方法重新变得极度精确。

2.2 符号因子 $\langle s \rangle_B$ 的定量崩溃行为

费米子符号因子 $\langle s \rangle_B$ 可以严格由所有零点到 $\xi = -1$ 和 $\xi = 1$ 的距离之比表示。由公式（18）可知：

在低温极限下，由于 $z_1 \approx -1 + e^{-\beta E_0}$，上式主导项显然为：

如图5所示，$\ln \langle s \rangle_B$ 与由 $z_1$ 单独贡献的分量 $\ln |\frac{z_1+1}{z_1-1}|$ 具有惊人的重合度（尤其是低温下，两条曲线完全贴合）。这用最无可辩驳的数学事实证明：符号问题的指数级恶化，在几何上完全是由最靠近 $\xi = -1$ 的那一个李-杨零点 $z_1$ 决定的。

2.3 熊恒域（Xiong et al.）等温能线外推法的失效分析

熊恒域等人曾提出一种基于等能线 $\xi_E(T)$ 的隐式外推方法（公式19），其通过拟合高 $T$ 区域的常数能量曲线来推导低温能级。然而本项研究指出，这种方法在本质上依然受到 LY 零点几何分布的强烈制约。

在图6所示的 $\xi-T$ 二维能差谱图中：

• 在 $T$ 较高或能差 $\Delta E$ 较大时，等能线极其平滑，外推工作良好。
• 在低温低能差区域（接近 $T \approx 0$ K, $\xi = -1$ 的左下角区域），由于实轴上密布着红色和灰色的零点轨迹，这些零点像电荷一样使得常能等高线发生极其剧烈的“电场弯曲”与汇聚（bunch up）。
• 如图7所示，作者测试了采用不同阶数（2阶至6阶）的 $\xi_E(T)$ 外推公式对 $N=4$ 系统低温费米能量的预测。在 $T < 0.6$ 的低温区，即使使用极高阶的6阶外推（紫线），其所得能量与真实能量的偏离依然巨大，产生了不可忽视的虚假转折，这正是由于等能线在这里被零点扭曲得无法用多项式拟合。

3. 代码实现细节与复现指南

为了方便量子物理与量子化学计算研究人员快速复现论文中的结果，本节提供一个完整的基于 Python 科学计算库（numpy, scipy, matplotlib）的计算与绘图流程指南。该流程无需商业软件，完全开源。

3.1 算法流程设计

复现的核心在于：

1. 高精度数值计算单粒子一维环配分函数 $Z_{11}(m\beta)$。
2. 根据群论循环结构公式构造 $Z^{(N)}(\beta, \xi)$ 多项式。
3. 在复平面上高精度求解多项式的零点位置。
4. 绘制零点随温度演化的轨迹图（如同图1、图2）。

3.2 核心复现代码（Python）

import numpy as np
from scipy.optimize import root
import matplotlib.pyplot as plt

# 1. 定义物理常数与单粒子配分函数 Z_11
# 设 h^2 / (2 * m * L^2) = 1 为能量单位，k_B = 1
def Z_11(beta, max_n=100):
    """计算自由粒子在一维环上的配分函数"""
    n_array = np.arange(-max_n, max_n + 1)
    energies = n_array**2  # E_n = n^2 (已无量纲化)
    return np.sum(np.exp(-beta * energies))

# 2. 构造多体广义配分函数 Z^(N)(beta, xi)
def get_Z_coefficients(beta, N):
    """根据论文公式(15)返回关于 xi 的多项式系数 [c_0, c_1, ... c_{N-1}]"""
    # 预先计算所需的 Z_11(m * beta)
    z = {m: Z_11(m * beta) for m in range(1, N + 1)}
    
    if N == 1:
        return np.array([z[1]])
    elif N == 2:
        # Z^(2) = 0.5 * (Z_11^2 + xi * Z_11(2beta))
        # 对应 c_0 = 0.5 * Z_11^2, c_1 = 0.5 * Z_11(2beta)
        return 0.5 * np.array([z[1]**2, z[2]])
    elif N == 3:
        # Z^(3) = 1/6 * [ Z_11^3 + 3*xi * Z_11 * Z_11(2beta) + 2*xi^2 * Z_11(3beta) ]
        return (1.0/6.0) * np.array([z[1]**3, 3*z[1]*z[2], 2*z[3]])
    elif N == 4:
        # Z^(4) = 1/24 * [ Z_11^4 + 6*xi*Z_11^2*Z_11(2beta) + xi^2*(8*Z_11*Z_11(3beta) + 3*Z_11^2(2beta)) + 6*xi^3*Z_11(4beta) ]
        c0 = z[1]**4
        c1 = 6 * (z[1]**2) * z[2]
        c2 = 8 * z[1] * z[3] + 3 * (z[2]**2)
        c3 = 6 * z[4]
        return (1.0/24.0) * np.array([c0, c1, c2, c3])
    elif N == 5:
        c0 = z[1]**5
        c1 = 10 * (z[1]**3) * z[2]
        c2 = 20 * (z[1]**2) * z[3] + 15 * z[1] * (z[2]**2)
        c3 = 30 * z[1] * z[4] + 20 * z[2] * z[3]
        c4 = 24 * z[5]
        return (1.0/120.0) * np.array([c0, c1, c2, c3, c4])
    else:
        raise NotImplementedError("仅提供 N <= 5 的解析系数")

# 3. 高精度多项式复零点求解函数
def solve_ly_zeros(beta, N):
    """利用 numpy.roots 求解复多项式零点"""
    coeffs = get_Z_coefficients(beta, N)
    # numpy.roots 需要系数从最高次幂到最低次幂排列，故需要逆序
    roots = np.roots(coeffs[::-1])
    return roots

# 4. 追踪零点随温度的演化 (以 N=4 为例)
N_particles = 4
T_space = np.linspace(0.25, 1.0, 76)  # 温度区间 T = [0.25, 1.0]
beta_space = 1.0 / T_space

all_zeros = []
for beta in beta_space:
    zeros = solve_ly_zeros(beta, N_particles)
    # 对解进行排序，方便画出连续轨迹
    zeros = np.sort_complex(zeros)
    all_zeros.append(zeros)

all_zeros = np.array(all_zeros)

# 5. 绘制 LY 零点轨迹
plt.figure(figsize=(8, 6))
for col in range(N_particles - 1):
    plt.plot(all_zeros[:, col].real, all_zeros[:, col].imag, 'o-', markersize=2, label=f'Root {col+1}')

plt.xlabel('Re $\xi$')
plt.ylabel('Im $\xi$')
plt.title(f'Lee-Yang Zeros Trajectories for N = {N_particles}')
plt.axhline(0, color='black', lw=0.5, ls='--')
plt.grid(True, alpha=0.3)
plt.legend()
plt.show()

3.3 开源相关工具与 Repo 推荐

在更广泛的实际强关联体系研究中，无法通过简单的解析多项式获得配分函数，此时需要利用多体路径积分算法生成数据。以下两个著名的开源量子蒙特卡洛/路径积分分子动力学软件包支持此类计算，并包含丰富的交换置换采样算法：

1. PIMD (Path Integral Molecular Dynamics): 由日本科学家开发的开源分子动力学代码，包含丰富的路径积分算法。 GitHub - pimd-code
2. i-PI: 著名的基于 Python 的分子动力学界面，支持先进的路径积分和玻色子/费米子置换项采样，非常适合用于获取广义交换宇称配分函数数据。 GitHub - ipi-code

4. 关键引用文献与局限性评论

4.1 关键参考文献

本项研究构建在以下几篇极具里程碑意义的学术成果之上：

1. [1] R. C. He, et al. Phys. Rev. E 113, 024115 (2026)：该系列工作的第一篇，系统论证了 $T = 0$ K 下李-杨零点的普适分布规律 $z_k = -1/k$，是本篇论文的理论基石。
2. [7] C. N. Yang and T. D. Lee, Phys. Rev. 87, 404 (1952) 与 [8] T. D. Lee and C. N. Yang, Phys. Rev. 87, 410 (1952)：经典的李-杨相变理论及零点定理，开创了在复平面分析热力学配分函数极点的先河。
3. [22] Y. Xiong and H. Xiong, Phys. Rev. E 107, 055308 (2023)：提出了等能线拟合外推法。本篇工作对其进行了直接的回应和深入的物理解释，指出了其在低温下失效的根本几何原因。
4. [35] B. Hirshberg, V. Rizzi, and M. Parrinello, Proc. Natl. Acad. Sci. USA 116, 21445 (2019)：提出了在路径积分蒙特卡洛中利用递归公式高效率采样多玻色子置换循环的算法，为广义配分函数的数值多项式构建提供了技术支持。

4.2 对本工作局限性与不足的深度学术评论

尽管本工作在一维无相互作用模型上取得了极其优美的解析成果，并且在物理上对符号因子的崩溃给出了颠覆性的几何解释，但若将其作为普适工具推向实际复杂的化学和物理体系，仍存在以下不可忽视的局限性：

1. 强相互作用对零点分布的重塑作用悬而未决： 论文聚焦于非相互作用（Non-interacting）自由粒子体系，这保证了多体积分可以完全因式分解为单粒子贡献（公式12）。然而，量子化学和凝聚态物理的核心在于强电子静电排斥（库仑相互作用）和电子-声子耦合。强关联相互作用会从根本上改变多体配分函数的系数 $c_j(\beta)$。尽管作者在引言中提到将在下一篇论文中讨论相互作用体系，但可以预见，在强关联区，零点在复平面的演化轨迹可能会出现由相互作用诱导的分叉甚至一阶相变，这会极大地增加零点轨迹的不可预测性。

2. 高维连续体系下的零点密集化问题： 一维粒子在有限环上有离散且稀疏的能级，这使得 LY 零点数量有限且轨迹清晰。但在三维大尺度体相系统（如温稠密物质、热致密等离子体）中，系统粒子数 $N \to \infty$，能级分布趋于连续化。此时，LY 零点不再是离散的点，而是会在复平面上连成密集的零点切割线（cuts）。这种连续的零点线将彻底隔离 $\xi \in [0, 1]$ 区域与 $\xi = -1$ 区域，这使得解析延拓的数学难度呈现几何级数增长。

3. 双步外推法对相变/交叉区（Phase Transition/Crossover）的依赖脆弱性： 作者提出的双步外推法的核心假设是：在低温拟合区间内，体系不发生任何物理相变。例如，公式（24）假设 $\ln \phi(\beta)$ 与 $\beta$ 呈高度平滑的线性或低阶多项式关系。然而，如果费米体系在降温过程中发生超导相变（BCS）、磁性相变（如反铁磁有序）或 Wigner 结晶，其自由能本身在该温度点是非解析的。此时，从高 $T$ 外推至低 $T$ 的热力学拟合关系会彻底失效。因此，该方法在探索具有复杂相图的超导、拓扑物态等前沿领域时，依然面临极大挑战。

5. 补充理论延伸与跨学科展望

5.1 交换宇称 $\xi$ 的物理延展：从任意子（Anyons）到分数排除统计（FES）

将交换宇称 $\xi$ 拓展到复平面，不仅是一个聪明的数学技巧，在低维拓扑量子物理中还具有极其深刻的物理实在性：

• 任意子物理（Anyons）：在二维空间中，粒子交换路径的编织群（Braid Group）是无限循环群。当粒子绕彼此旋转一周时，其波函数获得的相位因子可以是任意复数 $e^{i \theta}$。这对应于 $\xi$ 位于复平面单位圆 $|\xi| = 1$ 的物理场景。因此，论文中关于 $\xi$ 在复平面上演化的研究，直接对应于分数统计和任意子体系的热力学性质。
• 哈尔丹分数排除统计（Haldane’s FES）：哈尔丹提出的 exclusion statistics 描述了由于量子态空间被占满导致后续粒子无法进入的统计规律。在复 $\xi$ 平面上的不同点，可以看作代表了不同程度的“排除系数”。LY 零点在复平面上的演化，实际上为我们提供了一个窥探不同分数排除统计物性之间热力学平滑过渡的窗口。

5.2 “双步外推策略”在实际量子化学计算中的应用潜力

在现代量子化学中，对于具有强多配置（Multi-reference）特征的体系（如过渡金属催化剂、过渡态过渡区、多核铁硫簇），传统的密度泛函理论（DFT）往往失效，而精确的辅助场量子蒙特卡洛（AFQMC）又饱受符号问题的困扰。

本工作提出的“双步外推法”提供了一条实用性的折中路线：

1. 第一步（高温无符号区计算）：在较高温度下（此时电子热激发能大于其强关联能，李-杨零点远离实轴，符号问题微弱），利用 PIMC 进行 $\xi \in [0, 1]$ 的无符号蒙特卡洛采样。通过高精度多项式拟合，稳健地外推出高温费米子配分函数 $Z_F(T_\text{high})$ 和内能。
2. 第二步（温标热力学拟合）：结合量子化学在低温下可获得的近似能级谱，构建温标拟合方程（如公式25），将这些高温费米子数据沿温度轴外推。由于避开了在低温下沿 $\xi$ 实轴的外推，这能巧妙绕过零点引发的自由能对数发散，有望为大体系、强多配置电子关联问题的精确求解注入全新的生命力。

总结：本篇论文通过极简的一维粒子模型，用近乎完美的解析推导和深刻的静电学类比，为长久以来困扰学术界的费米子符号问题剥开了繁琐的数学迷雾。它告诉我们，符号问题不是一个单纯的统计涨落问题，其本质上是统计力学复平面上李-杨零点运动的几何必然结果。认识到这一点，并将视角从“对抗符号”转向“绘制零点轨迹”，或将成为未来量子多体计算科学实现飞跃的关键契机。