来源论文: https://arxiv.org/abs/2603.18998v1 生成时间: Mar 20, 2026 12:00

掺杂哈伯德模型热熵计算的统一框架：基于辅助场量子蒙特卡罗（AFQMC）的深度解析

0. 执行摘要

在强关联电子系统的研究中，热熵（Thermal Entropy）是刻画量子临界性、相变边界以及理解诸如莫特绝缘体（Mott Insulator）等复杂物态的核心热力学量。然而，在实验（如冷原子光晶格）中直接测量温度极其困难，熵反而成为更基本的物理准则。在理论计算领域，尤其是在偏离半满（Half-filling）的掺杂哈伯德模型中，传统的熵计算方法往往面临数值不稳定性或极高的计算成本。

近期，Song等人在其研究中提出了一套针对掺杂费米-哈伯德模型的高效熵计算框架。该框架基于巨正则系综（Grand-canonical Ensemble），利用辅助场量子蒙特卡罗（AFQMC）方法，通过温度、相互作用强度和化学势的三维参数空间路径积分，实现了对热熵的精确评估。本文将详尽解析该工作的理论基石、算法实现、Benchmark性能以及其在理解强关联系统物理图像中的关键作用。

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

1.1 核心科学问题

哈伯德模型（Hubbard Model）是描述强关联电子系统的标准模型，其哈密顿量包含动能项（hopping）和局域库仑排斥项（U）。尽管半满情形下的热熵计算已相对成熟，但在掺杂（Doped）情形下，由于引入了额外的自由度——化学势（μ）或费米填充数（n），系统展现出极为丰富的物理现象，如伪能隙（Pseudogap）、费米液体 crossover 以及可能的量子临界点（QCP）。如何在高维参数空间（T, U, μ）中高效、无偏地计算热熵，并利用熵的非单调行为揭示相图背后的物理，是本研究解决的核心科学问题。

1.2 理论基础：巨正则系综下的热力学势

研究的理论出发点是巨正则势密度 $\Omega$。根据定义：

$$\Omega = e - Ts + \mu n$$

其中 $e$ 是总能量密度，$s$ 是热熵密度，$n$ 是费米填充。根据热力学基本方程，其全微分形式为：

$$d\Omega = -s dT + h_I dU + n d\mu$$

这里 $h_I$ 对应于相互作用项的期望值（与双占据数 $D$ 相关）。通过对 $\Omega$ 进行数学上的全微分展开并与物理定义对比，作者推导出了热熵密度的全微分关系式：

$$ds = \frac{de_k + U dh_I + \mu dn}{T}$$

这一公式是整个计算框架的基石，它允许我们在参数空间中通过不同的积分路径来求解 $s$。

1.3 技术难点

  1. 高填充灵敏度：在巨正则系综中，填充数 $n$ 对化学势 $\mu$ 和温度 $T$ 极度敏感，尤其在低温区，细微的参数变动会导致填充数的显著偏离。
  2. 数值稳定性：传统的温度积分方案（Integration over T）在 $T \to 0$ 极限下容易发散，且需要高密度的采样点以保证积分精度。
  3. 负符号问题（Sign Problem）：在 Repulsive Hubbard 模型偏离半满时，蒙特卡罗采样会出现严重的符号问题，限制了可达到的低温范围。

1.4 四种计算方案详解

作者提出了四种互补的路径积分方案，涵盖了几乎所有的模拟场景：

  • 方案(a)：固定 $U$ 和 $n$ 随 $T$ 变化 这是最常用的路径。通过引入无穷大温度极限下的已知熵 $s_n(\infty)$（由填充数 $n$ 确定的原子极限熵），将积分分为高低温两段处理。利用分部积分法，作者将熵表达为关于内能 $e_{in}$ 的积分。此方案的优势在于利用了高温极限的解析解，减少了对 $T \to \infty$ 区域的模拟需求。

  • 方案(b)：固定 $U$ 和 $\mu$ 随 $T$ 变化 在某些计算方法中，固定 $\mu$ 比固定 $n$ 更容易实现。此时，化学势项 $\mu dn$ 会贡献到微分中。作者推导出类似的积分公式，通过总能量 $e_{tot}$（包含 $\mu n$ 项）来求解熵。

  • 方案(c)：固定 $U$ 和 $T$ 随 $\mu$ 变化 这是表征掺杂依赖性的关键路径。通过从半满（$\mu=0$）或空穴极限（$\mu \to \infty$）出发，对费米填充 $n$ 随化学势的变化进行积分。这在 AFQMC 中非常实用，因为填充数 $n(\mu)$ 通常具有很高的计算精度。

  • 方案(d)：固定 $T$ 和 $n$ 随 $U$ 变化 该方案通过 Hellmann-Feynman 定理，利用双占据数 $D$ 对 $U$ 的积分来求解熵的变化。这对于研究金属-绝缘体转变（Mott crossover）过程中的熵演变至关重要。

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

2.1 3D 非相互作用体系验证

首先，作者在 $U=0$ 的 3D 简立方晶格（$L=4$）上验证了方案(a)和(b)。由于非相互作用体系具有解析的费米-狄拉克分布，可以得到熵的精确解。计算结果显示，数值积分（蓝方块/红圆圈）与解析曲线（黑色虚线）在整个温度区间完美符合，验证了公式的数学正确性。

2.2 3D 相互作用体系的跨方案一致性

在 $U/t = 6$ 的强关联区域，作者选取了一个参考点 $(U/t=6, T/t=0.3, n=0.875)$，分别从温度路径、相互作用路径和化学势路径汇合于此。图7显示，四种不同路径得到的熵值在误差范围内高度一致。这种“闭合环路验证”是物理模拟中最严苛的 Benchmark，证明了该框架在强关联区域的稳健性。

2.3 2D 体系的大规模 Benchmark 与 DiagMC 对比

为了证明方法的实际应用价值，作者对 2D 平方晶格（$L=12, 16$）进行了大规模 AFQMC 模拟，并与图蒙特卡罗（DiagMC）的数据进行了对比（图9）。

  • 数据一致性：在 $U/t=2$ 和 $U/t=4$ 时，AFQMC 得到的热熵随化学势的变化曲线与 DiagMC 的结果几乎重合。
  • 精度优势：特别是在 $U/t=6$ 的情况下，DiagMC 由于级数收敛问题难以给出稳定的结果，而 AFQMC 凭借该框架能够给出高精度、低噪声的熵曲线，揭示了熵在 $\mu/t \approx 1.4$ 附近的峰值结构。

2.4 物理发现：熵峰与量子临界点（QCP）

计算数据表明，在特定掺杂浓度下，熵密度 $s(\mu)$ 展现出一个明显的峰值。通过麦克斯韦关系 $(\partial s / \partial n)_{U,T} = -(\partial \mu / \partial T)_{U,n}$ 的分析，作者指出这个熵峰对应于化学势随温度变化的驻点。这与动态平均场理论（DMFT）预测的掺杂诱导量子临界点密切相关，反映了系统从伪能隙状态向费米液体状态的演变。

3. 代码实现细节，复现指南与工具链

3.1 核心算法：有限温度 AFQMC

本工作的计算核心是有限温度辅助场量子蒙特卡罗。其基本流程如下：

  1. Trotter-Suzuki 分解：将配分函数 $Z = Tr(e^{-\beta H})$ 在虚时方向离散化为 $M$ 个切片，$\Delta\tau = \beta/M$。本研究采用了对称分解，误差阶数为 $O(\Delta\tau^2)$。
  2. Hubbard-Stratonovich (HS) 变换：将四费米子相互作用项通过辅助场（Auxiliary Fields）转换为单体项。作者采用了 spin-$s^z$ 通道的变换方式。
  3. 辅助场采样：使用带有延迟更新（Delayed Update）技术的马尔可夫链蒙特卡罗采样辅助场。
  4. 观测量测量：测量动能密度 $e_k$、双占据数 $D$ 和费米填充 $n$。

3.2 化学势 $\mu$ 的自动调节（Reproduction Guide）

为了在巨正则系综中实现固定填充 $n$ 的模拟，复现者需实现以下自洽调节算法：

  • 第一步：区间寻优。通过二分法（Bisection）在粗糙的 $\mu$ 区间内找到目标填充 $n_{target}$ 附近的值。
  • 第二步：泰勒展开修正。利用电荷压缩率 $\kappa = -\partial n / \partial \mu$（可通过密度-密度关联函数直接计算），通过公式 $\mu_{target} = \mu_0 + (n(\mu_0) - n_{target}) / \kappa(\mu_0)$ 进行快速收敛。通常只需 2-3 次迭代即可达到 $10^{-4}$ 级别的精度。

虽然论文未直接给出官方 Repo，但该算法可基于以下成熟的开源框架实现：

  • QUEST (Quantum Electron Simulation Tool)：经典的有限温度 AFQMC 实现，支持哈伯德模型。 Link
  • ALPS/QUEST：集成在 ALPS 项目中的版本，具有良好的并行性。
  • SmoQyQMC.jl：基于 Julia 语言的现代 QMC 软件包，非常适合复现此类算法路径积分。 Link

3.4 复现关键参数提示

  • 时步误差：对于双占据数 $D$ 的测量，Trotter 误差显著。建议采用 $\Delta\tau t = 0.015$ 甚至更小，或进行 $\Delta\tau \to 0$ 的外推。
  • 重加权（Reweighting）：在远离半满时，需开启符号重加权以处理平均符号 $\langle S \rangle < 1$ 的情况。

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

4.1 关键引用文献

  1. [13, 14] He et al. (2025)：作者团队之前在半满情形下的工作，奠定了相互作用积分法的基础。
  2. [53] Lenihan et al. (2021)：利用 DiagMC 研究了掺杂 2D 哈伯德模型，是本文主要的 Benchmark 对象。
  3. [50, 51] Mikelsons et al. (2009)：早期利用动力学团簇近似（DCA）探讨量子临界点的工作。
  4. [24, 54] Paiva et al. (2010)：确立了温度积分法在冷原子系统熵估算中的地位。

4.2 工作局限性评论

尽管本工作提供了一个极其强大且通用的框架，但在实际应用中仍存在以下局限：

  1. 符号问题的本质限制：虽然框架本身是通用的，但在极低温（$T/t < 0.2$）和显著掺杂区，AFQMC 的平均符号会迅速衰减至零。这意味着在最感兴趣的超导相区（Superconducting phase），该方法仍面临巨大挑战。
  2. Trotter 误差的敏感性：正如附录 A 所讨论的，熵的积分对双占据数的精度要求极高。对于复杂的模型（如包含 $t'$ 或各向异性项），Trotter 误差的控制将消耗大量的计算资源。
  3. 单带模型局限：本框架目前仅在单带哈伯德模型上验证，对于更接近真实材料的多带模型（Multiorbital models），全微分项将变得极其复杂，参数空间的路径选择难度也会增加。

5. 其他必要补充：物理图像与未来展望

5.1 麦克斯韦关系的深远意义

本工作不仅给出了计算公式，还系统推导并验证了多条麦克斯韦关系。例如：

$$\left( \frac{\partial s}{\partial U} \right)_{n,T} = -\left( \frac{\partial D}{\partial T} \right)_{n,U}$$

这条等式将微观的量子关联量（双占据数随温度的变化）与宏观的热力学量（熵随相互作用的变化）直接联系起来。这意味着通过观察双占据数随温度的升降（即 Pomeranchuk 效应），我们就可以预判系统热熵的演化趋势。

5.2 对冷原子实验的指导作用

在光晶格冷原子实验中，总原子数和绝热过程（即熵保持不变）是实验控制的基准。该框架允许实验学家根据测得的局域观测量（如密度分布、双占据率），反推出系统的绝热曲线（Isentropes），从而精确判定实验系统所处的温度范围。这对于寻找高温超导的冷原子模拟证据至关重要。

5.3 框架的扩展性

作者在结论中提到，该框架可以轻松扩展到包含次近邻跳迁（$t'$）或各向异性耦合（$t_z$）的广义哈伯德模型。通过在巨正则势 $\Omega$ 中添加相应的自然变量（如 $\partial \Omega / \partial t'$），可以构建类似的积分路径。这为未来研究三维哈伯德模型的相图提供了强有力的理论武器。

5.4 总结

Song等人的这项工作通过将热力学势的严谨理论与现代数值算法（AFQMC）相结合，终结了掺杂哈伯德模型熵计算中路径零散、精度不足的局面。其提出的四种方案不仅是计算工具，更是理解关联电子系统复杂热力学行为的物理指南。