来源论文: https://arxiv.org/abs/2606.10976v1 生成时间: Jun 10, 2026 16:47

Hubbard 模型中 Slater-Mott 交叉的低热力学特征：从双重占有数到反铁磁关联长度的深度解析

0. 执行摘要

强关联电子系统中的金属-绝缘体转变以及不同绝缘相之间的过渡，一直是凝聚态物理与理论量子化学研究的核心前沿。在具有嵌套费米面（Nested Fermi Surfaces）的二维和三维半满 Hubbard 模型中，任意微小的排斥相互作用 $U > 0$ 都会在基态诱导出 Néel 型反铁磁（AFM）长程序，从而使系统呈现绝缘行为。然而，这一绝缘态的微观物理机制在不同相互作用强度下有着本质的区别：在弱耦合区域，绝缘性源于费米面嵌套引起的自旋密度波不稳定性，即所谓的 Slater 绝缘体；而在强耦合区域，电荷涨落被极大地抑制，低能物理受局部磁矩之间的超交换作用（Superexchange）支配，即经典的 Mott 绝缘体。由于这两者在 $T = 0$ 时被一个光滑的交叉（Crossover）连接，如何定量、精确且在实验上可行地界定这一 Slater-Mott 交叉边界，一直是一个悬而未决的理论难题。

近期发表的一项杰出工作通过结合精确对角化（Exact Diagonalization, ED）与大尺度辅助场量子蒙特卡洛（Auxiliary-Field Quantum Monte Carlo, AFQMC）模拟，提出并验证了一个简洁、定量受控且在冷原子光晶格实验中高度可测的热力学判据。该研究表明，在极低温度下，双重占有数对温度的导数 $(\partial D / \partial T)_U$ 的符号改变可以完美作为 Slater-Mott 交叉的标志。具体而言：

• 在 Slater 机制下，$(\partial D / \partial T)_U > 0$，反映了电荷涨落起主导作用，低能激发呈电荷主导特征；
• 在 Mott 机制下，$(\partial D / \partial T)_U < 0$，展现出反常的 Pomeranchuk 效应，标志着低能自旋超交换物理的支配地位。

通过对二维半满方晶格 Hubbard 模型进行系统模拟，研究人员成功将该基态交叉边界精准定位在 $U_{\text{cross}}/t = 4.0(2)$。此外，他们独立地利用反铁磁关联长度 $\xi(U)$ 在低温下的极大值行为，得到了完全一致的交叉边界（$U_{\xi}^*/t = 3.94(20)$）。这两者通过麦克斯韦关系（Maxwell Relation）在热力学熵 $s(U)$ 的局部极小值处得到了完美的理论统一。本博客将对该工作的核心科学问题、理论基础、数值方法、基准测试及未来的实验应用进行全方位的深度技术剖析。

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

1.1 核心科学问题：Slater 绝缘体与 Mott 绝缘体的本质区别与交叉

理解相互作用如何改变费米子的物理行为是多体物理学最根本的任务之一。Hubbard 模型的哈密顿量可以写作：

其中 $t$ 代表最近邻跃迁强度，$U$ 代表格点上的库仑排斥能。在半满（Half-filling）状态下，二维方晶格具有完美的费米面嵌套（巢状不稳定性），嵌套矢量为 $\mathbf{Q} = (\pi, \pi)$。这意味着，在 $T=0$ 时，任何微小的相互作用 $U > 0$ 都会导致单粒子激发谱中打开一个能隙，并自发产生 Néel 反铁磁长程序。

尽管对于所有的 $U > 0$，基态都是反铁磁绝缘体，但其背后的物理机制随 $U$ 的增大发生剧变：

1. Slater 机制（弱耦合，$U \ll t$）：绝缘性的成因是平均场水平下的自旋密度波（SDW）。由于费米面嵌套，电子自旋在磁激发下产生非平庸的相干叠加，导致电荷涨落依然显著（双重占有数较多）。其激发谱主要由电荷类型的激发主导，能隙大小 $\Delta_{\text{sp}} \propto t \exp(-2\pi\sqrt{t/U})$。
2. Mott 机制（强耦合，$U \gg t$）：电荷自由度在低能下彻底冻结，每个格点几乎被单个电子占据，双重占有被极大地抑制。系统的物理性质可以用有效海森堡模型（Heisenberg Model）描述，局部磁矩之间的自旋超交换相互作用耦合常数 $J = 4t^2/U$。其激发起于自旋翻转（Spin-flip），属自旋类型激发。

长期以来，实验和理论上界定这一交叉的方法多依赖于谱学性质（如单粒子谱函数的分裂）或系统总能量的二阶导数。这些方法要么对数值近似过于敏感，难以提供干净的定量标准，要么在实验上（如冷原子光晶格）极难直接测量。因此，寻找一个既有严格理论支撑、又能在多体数值计算中精确把控、同时还能直接在实验中观测的热力学判据，成为了该领域的燃眉之急。

1.2 理论基础：双重占有数与热力学熵的麦克斯韦关系

为了构建热力学判据，作者引入了双重占有数（Double Occupancy） $D$：

在热力学体系中，哈密顿量中与相互作用项耦合的正是 $D$。因此，体系的自由能密度 $f = F/N_s$ 关于相互作用强度 $U$ 的偏导数直接对应于双重占有数：

结合热力学基本关系，自由能关于温度 $T$ 的偏导数为熵密度 $s = S/N_s$：

由于自由能是状态函数，其二阶混合偏导数与求导顺序无关，因此我们得到了至关重要的麦克斯韦关系式：

这一公式将热力学熵对相互作用的响应，与双重占有数对温度的响应完美地联系在了一起。正是这个优雅的公式，为定位 Slater-Mott 交叉提供了直接的物理钥匙。

在低温下，由于第三热力学定律，当 $T \to 0$ 时，熵 $s \to 0$。但在有限但极低的温度下，熵密度 $s(U)$ 随 $U$ 的变化呈现出非单调的特征（见图1(a)的插图）。具体表现为：

• 在 $U$ 较小时（Fermi Liquid 及 Bad Metal 区），随着 $U$ 增大，单粒子能隙打开，电荷涨落受到抑制，自旋自由度开始呈现局域化，导致熵呈现增长或高位波动；
• 随着 $U$ 进一步增大进入准反铁磁绝缘区，自旋之间的关联增强。在某一临界相互作用 $U_S$ 处，自旋反铁磁关联达到最强，对应的熵表现出一个局部极小值（Entropy Dip）；
• 超过 $U_S$ 后，由于超交换相互作用 $J = 4t^2/U$ 随 $U$ 增大而减弱，在固定温度 $T$ 下，有效无序度增加（相当于温度相对于 $J$ 升高了，即 $T/J$ 增大），自旋熵反而重新升高。

因此，在 $U = U_S$ 处，$(\partial s / \partial U)_T = 0$。根据麦克斯韦关系，这直接等价于 $(\partial D / \partial T)_U = 0$。这就意味着，双重占有数在温度轴上存在一个极值点 $T_D$，当 $T < T_D$ 时，导数符号发生翻转。

1.3 技术难点与物理机制解析

要在微观上理解为什么 $(\partial D / \partial T)_U$ 的符号在 Slater 和 Mott 区域截然相反，必须诉诸体系的低能激发谱结构。作者通过**精确对角化（ED）**给出了极为直观的微观解释。

在低温下，配分函数可以被截断到仅包含基态（$m=0$）和前几个低能激发态（$m \ge 1$）。双重占有数的温度依赖性可以近似写为：

其中 $E_m$、$D_m$ 和 $g_m$ 分别为第 $m$ 个激发态的能量、双重占有数和简并度，$D_0$ 为基态双重占有数。此式表明，双重占有数的温度行为完全由低能激发态的双重占有数与基态之差 $(D_m - D_0)$ 的正负决定：

• Slater 机制（弱耦合，如 $U/t = 4$）：电荷涨落显著。激发态表现出明显的电荷重新分布特征。随着激发能增高，系统克服能隙，产生更多的电荷涨落，使得低能激发态具有比基态更高的双重占有数，即 $D_m > D_0$。因而，加热会导致双重占有数增加，$(\partial D / \partial T)_U > 0$。
• Mott 机制（强耦合，如 $U/t = 12$）：电荷涨落高度冻结。基态的双重占有数主要来源于虚跃迁过程带来的高阶修正，即 $D_0 \propto (t/U)^2$。当系统受到热激发时，低能激发主要是自旋翻转。自旋翻转会在背景中形成“铁磁键”，根据泡利不相容原理，这会剧烈抑制虚跃迁通道，从而使激发态的双重占有数显着低于基态，即 $D_m < D_0$。因而，加热会抑制虚跃迁，导致双重占有数减小，$(\partial D / \partial T)_U < 0$（这在热力学上对应于著名的 Pomeranchuk 效应）。

这种从“电荷型激发”向“自旋型激发”的转变，完美地对应了 Slater 向 Mott 绝缘体的过渡。

1.4 方法细节：辅助场量子蒙特卡洛（AFQMC）与纠缠通道技术

尽管 ED 提供了清晰的物理解释，但其受限于极小的空间尺度（如 $4 \times 4$），有限尺寸效应极为严重。为了获得热力学极限下的精确物理量，必须采用大尺度数值计算。该工作采用了基于虚时间投影和热力学配分函数集成的**辅助场量子蒙特卡洛（AFQMC）**方法。

为了彻底压制 QMC 模拟中的统计涨落，并精确计算双重占有数和反铁磁关联函数，作者使用了一种自旋通道（HS-$s^z$）与电荷通道（HS-$\hat{n}$）混合实现的哈密顿量离散技术（Hubbard-Stratonovich Transformation）。

• 对于与电荷密度和双重占有数相关的物理量（如 $D$、能量 $E$、熵 $s$），采用 HS-$s^z$ 变换。该变换将相互作用项解耦到局域磁场空间，能极大抑制电荷测量的方差；
• 对于自旋物理量（如自旋-自旋关联函数、反铁磁关联长度 $\xi$），采用电荷通道 HS-$n$ 变换，该变换将相互作用解耦到局域静电势空间，能显著提高自旋关联的测量精度。

此外，针对半满 Hubbard 模型易受到边界条件调制产生有限尺寸能隙的问题，模拟引入了扭曲平均边界条件（Twist-Averaged Boundary Conditions, TABC）。在计算双重占有数及其对温度的导数时，通过在布里渊区内均匀取样大约 $30$ 套扭曲角，极大地消除了有限尺寸效应，使得即使在 $L=20$ 的尺度下，所得结果也已经极度逼近热力学极限。

2. 关键 benchmark 体系，计算所得数据，性能数据

本研究重点对二维半满方晶格 Hubbard 模型进行了系统而精密的基准测试，并进一步拓展到了三维立方晶格体系。以下是该工作中沉淀的核心计算数据与关键物理性能。

2.1 精确对角化（ED）基准数据：$4 \times 4$ 晶格分析

研究者首先对 $L=4$ 的二维半满系统（具有周期性边界条件 PBC）进行了全谱对角化，得到了基态和最低三个激发能级的精确参数。这些基准测试数据如下表所示：

表1：$L=4$ 方晶格 Hubbard 模型在弱耦合 $U/t=4$ 时的低能激发态性质

表2：$L=4$ 方晶格 Hubbard 模型在强耦合 $U/t=12$ 时的低能激发态性质

分析：如上两表所示，对于 $U/t=4$，每一个激发态的 $D_m$ 均大于基态 $D_0$；而对于 $U/t=12$，每一个激发态的 $D_m$ 均小于 $D_0$。这直接证明了前述热力学符号转变机制的正确性。这一微观物理图像的转变点在 $L=4$ 时发生于 $U/t \sim 6 - 7$。然而，由于 $L=4$ 系统尺寸太小，单粒子能隙过大，该转变点严重偏离了热力学极限下的真实交叉边界。

2.2 AFQMC 模拟：二维热力学极限下的双重占有数测量

为了消除尺寸效应并获得精确的边界，作者采用大尺度 AFQMC 模拟了 $L=20$ 的体系并应用了 TABC。通过直接测量关联量 $\langle \hat{H}\hat{D} \rangle - \langle \hat{H} \rangle \langle \hat{D} \rangle$，他们直接得到了 $(\partial D / \partial T)_U$。以下是核心热力学响应数据（如图3所示）：

• $U/t = 3.5$：双重占有数差值 $[D(T) - D(0)]$ 随着温度升高而稳步上升。在 $T \to 0$ 极限下，其斜率明显为正，即 $(\partial D / \partial T)_U > 0$，属于典型的 Slater 绝缘体 行为。
• $U/t = 4.0$：在 $T/t = 0$ 到 $0.06$ 之间，$[D(T) - D(0)]$ 保持在误差范围内几乎不变，斜率几近为零。这标志着该点处于绝缘体交叉的中心界线上。
• $U/t = 4.5$ 及 $4.75$：在低温区（$T/t < 0.06$），$[D(T) - D(0)]$ 随温度升高单调下降，呈现出清晰的负斜率，即 $(\partial D / \partial T)_U < 0$，系统进入 Mott 绝缘体 区域。

表3：不同相互作用下提取出的双重占有数极小值温度 $T_D/t$ 与熵极小值位置 $U_S/t$

通过拟合 $D(T)$ 的局部极小值点 $T_D$ 和熵极小值位置 $U_S$，得到了如下精细的数据点（部分展示于表I）：

通过对 $U_S(T)$ 和 $T_D(U)$ 进行线性外推（见图1中的红虚线），在 $T \to 0$ 时，两者汇聚于同一临界值：

这一极为干净的数据确认了二维半满方晶格 Hubbard 模型在热力学极限下的基态 Slater-Mott 交叉点。

2.3 反铁磁（AFM）关联长度 $\xi(U)$ 极大值的互证

为了使这一热力学判据更具说服力，作者独立地使用反铁磁关联长度 $\xi(U)$ 进行验证。$\xi$ 是通过自旋结构因子 $S_{zz}(\mathbf{q})$ 计算得到的：

其中 $\boldsymbol{\pi} = (\pi, \pi)$，而 $\delta\mathbf{q} = (2\pi/L, 0)$。基于 $L=20$ 和 $L=24$ 的大尺度 PBC-AFQMC 模拟，$\xi(U)$ 在不同温度下的表现如下（图4）：

• 在固定温度 $T/t$ 下，随着相互作用 $U/t$ 的增大，$\xi$ 呈现先上升后下降的单调特征，并在某一处 $U_\xi$ 达到极大值。这一极大值反映了反铁磁长程序的形成与自旋超交换作用减弱之间的竞争。
• 随着温度降低：
  • $T/t = 0.3 \implies U_\xi/t = 5.9(4)$;
  • $T/t = 0.2 \implies U_\xi/t = 5.4(1)$;
  • $T/t = 0.1 \implies U_\xi/t = 4.6(2)$。

对 $U_\xi(T)$ 数据进行线性外推到 $T \to 0$ 极限，得到了独立的外推交叉边界：

该结果与基于热力学双重占有数外推得到的 $U_{\text{cross}}/t = 4.0(2)$ 在误差范围内完美契合！这两种完全不同的物理量（一个是热力学局域量，一个是空间长程关联量）在交叉边界上的高度契合，揭示了该系统中低能物理的自洽性。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

为了在学术研究中复现这一里程碑式的工作，本节将提供详细的代码复现指南、核心算法实现架构，并列出推荐使用的开源多体计算软件包。

3.1 核心算法实现：精确对角化（ED）

由于 ED 能够无系统误差地给出小体系（如 $L=4$）的全部本征态，它是验证微观机制的第一步。在 half-filling 条件下，对于 $4 \times 4$ 晶格（共 16 个格点），自旋向上和向下的电子数均为 8。对应的 Hilbert 空间维度为：

这一维度对于现代台式机而言是完全可以轻松处理的。我们可以利用 Lanczos 算法计算基态以及前几个低能激发态。

Python-SciPy 伪代码复现框架：

import numpy as np
from scipy.sparse.linalg import eigsh
from scipy.sparse import csr_matrix

def build_hubbard_hamiltonian(Lx, Ly, U, t):
    """
    构建二维半满 Hubbard 模型的稀疏矩阵表示。
    这里为了简化说明，采用全基矢构建，实际计算中应充分利用空间对称性(如平移、自旋反转)和粒子数守恒。
    """
    Ns = Lx * Ly
    # 1. 生成所有符合 (N_up, N_down) = (Ns/2, Ns/2) 约束的 Fock 态
    # 2. 构造跃迁项 (Hop-term)
    # 3. 构造库仑排斥项 (Interaction-term)
    # 4. 利用 scipy.sparse.csr_matrix 组装哈密顿量矩阵 H
    pass

# 参数设置
Lx, Ly = 4, 4
t = 1.0
U_values = [4.0, 12.0]

for U in U_values:
    H = build_hubbard_hamiltonian(Lx, Ly, U, t)
    # 求解前 4 个本征值和本征态 (k=4)
    energies, wavefunctions = eigsh(H, k=4, which='SA')
    
    print(f"--- U/t = {U} ---")
    for m in range(4):
        # 计算每个能级对应的双重占有数 D_m
        # D_m = <phi_m| \hat{D} |phi_m>
        psi = wavefunctions[:, m]
        # build_double_occupancy_operator 构建双重占有算符 D_op
        # D_m = np.dot(psi.conj().T, D_op.dot(psi))
        # print(f"m={m}: E_m={energies[m]:.6f}, D_m={D_m:.6f}")

对于 $4 \times 4$ 体系，可直接调用以下成熟的开源精确对角化软件包：

• EDIpack: 用于强关联系统的通用精确对角化和动力学平均场论（DMFT）杂质求解器。
  GitHub Link: https://github.com/QDMRG/EDIpack
• HPhi: 极为强大的通用量子多体哈密顿量精确对角化软件包，由东京大学物性研究所维护，专门针对 Hubbard 模型等进行了极致的并行和优化。
  GitHub Link: https://github.com/issp-center-dev/HPhi

3.2 核心算法实现：辅助场量子蒙特卡洛（AFQMC）

对于 $L \ge 20$ 的大系统，精确对角化完全失效，必须使用量子蒙特卡洛。本文采用的有限温度辅助场量子蒙特卡洛（Finite-T AFQMC）主要基于 Determinant QMC (DQMC) 框架。

复现步骤与技术路线：

1. 时空离散：设定虚时间步长 $\Delta\tau t = 0.02$，投影片数 $M = \beta / \Delta\tau$。采用对称 Trotter 离散将配分函数写为路径积分形式。

2. Hubbard-Stratonovich 变换：

   • 在计算双重占有数等热力学量时，引入 Ising 辅助自旋场 $s_i(\tau) = \pm 1$。公式见原论文公式 (4) 对应的 HS-$s^z$ 通道。
   • 在计算自旋关联和关联长度时，切换至公式 (5) 对应的 HS-$n$ 电荷通道。

3. 扭曲平均边界条件（TABC）的加入： 在单粒子跃迁项中加入相位因子：$t_{ij} \to t_{ij} e^{i \boldsymbol{\theta} \cdot (\mathbf{r}_i - \mathbf{r}_j)}$，其中 $\boldsymbol{\theta} = (\theta_x, \theta_y)$。对多个 $\boldsymbol{\theta}$ 采样（通常在 $[-\pi, \pi)^2$ 区间内进行 $6 \times 6$ 的均匀格点取样），对格林函数进行平均。

4. 无数值微分的热力学量计算： 为了避免对具有随机误差的 $D(T)$ 数据进行数值求导，利用关联测量公式 (11)：

   $$\left(\frac{\partial D}{\partial T}\right)_U = \beta^2 \left( \langle \hat{H}\hat{D} \rangle - \langle \hat{H} \rangle \langle \hat{D} \rangle \right)$$

   该量在 AFQMC 中需要计算四费米子与双重占有算符的关联，对应于格林函数的两点与四点关联测量，需利用 Wick 定理进行收缩。

开源软件包推荐：

• ALF (Algorithms for Lattice Fermions)： 这是目前凝聚态物理界公认最成熟、功能最完备的开源格点费米子量子蒙特卡洛模拟框架（由 Würzburg 大学等国际一流团队开发和维护）。它原生支持自旋和电荷通道的 HS 变换，支持各种复杂的格点几何结构（包括二维方晶格和三维立方晶格），并完美兼容有限温度和基态投影算法，同时具备 Twist-Averaged Boundary Conditions 的实现接口。
  ALF 官方 Git 仓库: https://git.physik.uni-wuerzburg.de/ALF/ALF
  GitHub 镜像: https://github.com/ALF-collaboration/ALF
• Quest (Quantum Electron Simulation Toolbox)： 一个经典的 Fortran 编写的格点 Hubbard 模型 DQMC 代码包，结构简单清晰，非常适合作为复现和魔改的起点。
  GitHub Link: https://github.com/quest-qmc/quest

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键引用文献及物理脉络

该研究的成功并非空中楼阁，它深深根植于以下经典工作之中：

1. Georges et al. (1996) [文献1]：建立了动力学平均场论（DMFT）的完备框架，这是现代强关联物理处理金属-绝缘体转变（Mott Transition）的理论基石。
2. Imada, Fujimori, & Tokura (1998) [文献2]：关于金属-绝缘体转变的百科全书式综述，详细梳理了 Slater 和 Mott 绝缘体在不同物理体系中的竞争。
3. Mazurenko et al. (2017, Nature) [文献24]：首次在冷原子光晶格实验中实现并确认了二维 Fermi-Hubbard 模型中的长程反铁磁 Néel 态，为本文所提热力学标准的实验验证开辟了物理舞台。
4. Shao et al. (2024, Nature) [文献28]：在三维费米-Hubbard 光晶格实验中成功观测到反铁磁相变，这使得本论文中关于三维系统 $U_{\text{cross}}^{\text{3D}}/t = 5.9(3)$ 预测的验证近在咫尺。

4.2 局限性深入评论

尽管这是一项概念清晰、论证严密的高水平研究，但作为面向量子化学和计算材料学研究者的技术视点，我们仍需清醒认识到该方法的以下核心局限性：

1. 严重依赖“无负符号问题”的半满条件： 整个研究的核心计算工具是 AFQMC，而蒙特卡洛能够在极低温度（$T/t \sim 0.05$）下保持极高精度的前提，是系统处于半满（Half-filling）状态。一旦系统引入化学势偏离半满（即掺杂，Doped Hubbard Model），QMC 就会遭遇极其严重的负符号问题（Sign Problem）。这意味着在非半满情况下，无法通过 AFQMC 计算出温度极低时的 $(\partial D / \partial T)_U$。因此，该判据目前极难直接推广到高温超导最核心的掺杂 Hubbard 区域。
2. 外推过程中的系统误差与不确定性： 在确定基态 $U_{\text{cross}}$ 时，研究者采用的是对有限温度下的 $T_D(U)$ 或 $U_S(T)$ 数据点进行线性外推。虽然这一做法在物理直觉上是合理的，但不能完全排除在极低温度极限（如 $T/t < 0.05$）下，系统可能存在非线性的重整化行为或量子临界物理，导致外推值发生微小的系统性偏移。未来的研究需要使用如 Tensor Network（张量网络）等能在 $T=0$ 热力学极限下直接计算的算法进行交叉检验。
3. 实验精度的巨大挑战： 尽管作者乐观地指出最近的冷原子光晶格实验已经实现了低于 0.1% 的双重占有数测量精度，但在真实实验中，光晶格存在不可避免的外置简谐谐振阱（Harmonic Trapping Potential）。这意味着实验样品具有不均匀的密度分布，其中既包含半满区也包含非半满区。外加非均匀势的存在会对双重占有数的局部热力学求导信号造成强烈的平滑化（Smearing Effect），从而极大地增加了在真实实验中直接“可视化”该交叉边界的噪声。

5. 其他必要的补充：Pomeranchuk 效应与自旋-电荷熵竞争的深层图景

作为对本项工作的补充，有必要更深一步探讨 Pomeranchuk 效应 的深层物理本质，以及其在量子化学计算和关联材料设计中的广泛启示。

5.1 费米子的 Pomeranchuk 效应

Pomeranchuk 效应最早由苏联物理学家 Pomeranchuk 于 1950 年提出，用于预测固体 $^3\text{He}$ 在加压加热时会发生反常熔化的物理现象。在普通的液体中，粒子的热运动导致液体熵大于固体熵，因此加热固体通常会导致其熔化。然而，对于 $^3\text{He}$，固相中的原子核自旋是局部化的且高度无序的，贡献了巨大的自旋熵 $k_B \ln 2$。而在液相中，由于费米简并，只有费米面附近的少数原子能贡献自旋涨落，其熵表现为与温度成正比的线性行为（Fermi Liquid $S \propto T$）。因而在某一低温区间，固态 $^3\text{He}$ 的熵反而大于液态。此时，加热系统反而会促进液态向固态转化，或者说在恒压下加热系统，体积反而会收缩，这正对应了负的热膨胀系数。

在强关联 Hubbard 模型的 Mott 绝缘态中，我们也看到了完全相同的物理机制：

• 在 Mott 绝缘区（强耦合），低能自由度完全是自旋式的，形成局部磁矩（类似于固体 $^3\text{He}$）；
• 在 Slater 绝缘区或 Bad Metal 区，电子保留了部分流动电荷特征（类似于液体 $^3\text{He}$）。

因此，在 Mott 区，自旋熵占绝对主导。在温度升高时，局部自旋的无序化释放了熵，根据麦克斯韦关系，系统倾向于将电子更加局域化，导致电荷涨落（双重占有数 $D$）被剧烈抑制。这种“自旋释放熵压制电荷涨落”的过程，是强关联凝聚态物理中最为美妙的景象之一，而 $(\partial D / \partial T)_U < 0$ 正是这一图像最直接、最无争议的判据。

5.2 对量子化学及过渡金属氧化物研究的启示

在现代量子化学中，研究过渡金属氧化物（如高价铜氧化物、镍基超导体催化剂、以及具有重 $\alpha$-d 轨道的配合物）的电子结构时，准确区分体系属于 Slater 型（磁程序诱导能隙）还是 Mott-Heisenberg 型（强定域多体相关，活性空间中电荷涨落完全冻结），对于选择合适的计算方法至关重要：

• 如果体系处于 Slater 机制，密度泛函理论（DFT）辅以适当的杂化泛函（如 B3LYP 或 HSE06）或非限制性 Hartree-Fock（UHF），能够相对准确地描述磁相变和单粒子能隙；
• 如果体系进入 Mott 机制，单决定域的平均场方法会完全失效，即使加入 Hubbard $U$ 修正（DFT+U），也无法正确处理电荷冻结引发的丰富多体关联。此时必须采用多组态自洽场方法（CASSCF/CASPT2）、密度矩阵重整化群（DMRG）或动力学平均场论（DMFT）等强关联工具。

该研究提供的高效热力学边界 $U_{\text{cross}}/t \approx 4.0$ 表明，当活性轨道的有效关联强度超过跃迁积分的 4 倍时，研究者必须坚决放弃传统的单粒子泛函计算，转向多中心多载流子的高阶强关联量子化学算法。这一来自格点 Hubbard 模型的定量标度，无疑为关联分子体系的计算方案选择树立了一座明晰的灯塔。