来源论文: https://arxiv.org/abs/2606.30607v1 生成时间: Jul 04, 2026 05:05

受挫晶格与低密度展开下的巡游铁磁性：Hubbard模型的严格格林函数与相变深度解析

0. 执行摘要

解释过渡金属（如铁、钴、镍）中巡游电子的铁磁性，一直是凝聚态物理和量子化学理论中最具挑战性的经典难题之一。与定域自旋模型（如Ising模型和Heisenberg模型）不同，巡游铁磁性（Itinerant Ferromagnetism）涉及移动电子动能与库仑排斥能之间的微妙竞争。Hubbard模型作为描述强关联电子最简明而深刻的模型，尽管其哈密顿量形式简单，但除了一维体系和极少数特殊极限（如Nagaoka极限、平带铁磁性）外，其基态性质至今缺乏普遍的严格解析解。

本研究探讨了在非定域、受挫晶格体系中，基于严格的“低密度展开”（Low-Density Expansion）理论框架下Hubbard模型的基态铁磁有序化机制。本博文对该论文进行深度学术剖析，核心内容包括：

1. 理论根基：将连续介质中的Lieb-Seiringer-Solovej定理以及简单立方晶格中的Giuliani-Seiringer-Yin严格低密度能量展开公式，拓展至具有强受挫特征的晶格体系（包括次近邻跃迁NNN、对角近邻跃迁DNN和面心立方FCC晶格）。
2. 方法论创新：通过引入解析格林函数（Green’s Functions）和多维布里渊区数值积分技术，精确计算出受挫体系的态密度（DOS）。基于此，推导出了超越Stoner平均场理论的严格基态能量Instability判据，确立了临界态密度 $D_{\text{crit}} = 1/(8\pi a)$ 的解析形式。
3. 核心发现：研究证明，在简单立方晶格中，物理上合理的排斥能 $U$ 无法诱导铁磁性；然而在受挫晶格（如 $t_3 = 0.245t$ 的对角近邻晶格和 $t_2 = -0.25t$ 的FCC晶格）中，铁磁不稳定性在极低的电子密度（如 $\rho_c \approx 0.031$ 和 $\rho_c \approx 0.09$）下即可发生。这一发现不仅证实了“晶格受挫能够显著稳定巡游铁磁性”的物理直觉，而且由于相变发生在低密度区间，极大地提高了该摄动展开理论的物理可信度。

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

1.1 核心科学问题：巡游铁磁性的起源与晶格几何受挫

在金属体系中，传导电子既具有波动性（在空间中自由流动以降低其动能），又具有粒子性（在相遇时产生强烈的库仑排斥能）。巡游铁磁性的核心物理机制源于泡利不相容原理与库仑排斥的协同效应。当体系极化为铁磁态时，所有电子自旋同向，空间波函数必须全反对称，从而迫使电子在空间上自然避开彼此，降低了相互作用能（通常称为交换能作用），但代价是占据了更高的动量状态，导致动能升高。这就是经典的Stoner机制。

然而，Stoner判据：

是在平均场近似（Mean-Field Approximation）下推导出来的，它严重低估了电子间的动态关联效应，导致其预测的铁磁相区在实际材料和更严格的计算中往往不复存在。平均场理论将相互作用强度等同于裸库仑势 $U$，但在真实的低密度极限下，强关联效应会使相互作用被重整化为散射长度（Scattering Length） $a$ 所表征的有效强度。由于有效散射强度受到两体关联的强烈压制，导致在普通简单立方晶格（Simple Cubic）中，低密度下不可能发生巡游铁磁性。

本工作的核心科学问题在于：是否存在某些特定的晶格结构（例如受挫晶格），通过改变单粒子色散关系，在带底附近产生极高甚至发散的态密度（DOS），使得重整化后的有效排斥力在低密度极限下依然能够战胜动能的增加，从而严格诱导基态铁磁长程序？

1.2 理论基础：严格的低密度展开

为了超越平均场近似并给出严格的微观描述，本研究采用了基于多体微扰论的基态能量低密度展开理论。对于连续介质中的三维弱相互作用费米气体，Lieb、Seiringer和Solovej严格证明了在热力学极限下，基态能量密度具有如下严格的渐近形式（简记为LSS定理）：

其中 $\rho_\uparrow$ 和 $\rho_\downarrow$ 分别是自旋向上和自旋向下的电子密度，$a$ 是低能s波两体散射长度。该公式的第一项是无相互作用费米气体的动能项（简记为 $e_0(\rho)$），表现为典型的 $\rho^{5/3}$ 行为；第二项代表有效散射能，其相互作用强度完全由散射长度 $a$ 控制，而非裸库仑势 $U$。

Giuliani（对于上界）以及Seiringer与Yin（对于下界）成功将此定理推广到了格点Hubbard模型。在具有最邻近跃迁的简单立方晶格上，基态能量展开公式具有极其相似的数学结构：

在此处，晶格体系的特征完美地融入了无相互作用项 $E_0(\rho)$ 以及两体有效散射长度 $a$ 中。这正是格点模型相比连续介质的巨大魅力所在：连续介质中的 $e_0(\rho) \sim \rho^{5/3}$ 是机械锁定的，无法进行任何工程设计；而格点上的 $E_0(\rho)$ 则直接依赖于晶格的单粒子色散关系 $\epsilon(\mathbf{k})$。通过人工设计晶格受挫（Frustration），我们可以极大程度地重塑色散关系，从而使动能增加项在低密度下变得极为平缓，进而使有效排斥项占据主导，触发铁磁相变。

1.3 技术难点一：非各向同性晶格色散关系的解析反演与能带平坦化

在连续介质中，等能面是完美的球形，能量和动量具有简单的各向同性二次元关系 $\epsilon(\mathbf{k}) = \frac{\hbar^2 \mathbf{k}^2}{2m}$。这使得将非相互作用能量写为密度的函数变得极其平庸。然而，在受挫三维格点模型中，色散关系 $\epsilon(\mathbf{k})$ 不仅高度非平庸，还具有复杂的空间各向异性（各向异性非等能面）。因此，无法直接得到密度的解析表达式 $\rho(\epsilon_F)$ 和非相互作用能量 $E_0(\rho)$。为了解决这一问题，必须借助格林函数技术，首先计算出全局态密度 $D(E)$。

1.4 技术难点二：散射长度 $a$ 在受挫晶格中的重整化解析定义

晶格系统中的两体散射长度 $a$ 不再是简单的常数，而是紧密依赖于格子上的两体传播子（格林函数）。格点散射长度的精确定义为：

其中，重整化因子 $\gamma$ 是描述在布里渊区（1BZ）中两粒子激发特征的积分项：

其中 $e_{\text{min}}$ 代表能带的底边界（最底能量）。可以明显看出，当裸库仑势 $U \to \infty$ 时，有效散射能 $8\pi a \to 1/\gamma$。这代表了多体理论对局域强关联的极限约束，即由于强关联压制，有效相互作用强度存在一个由能带几何结构决定的坚硬上限，这在经典的平均场理论中是被完全忽略的。

1.5 方法细节：格林函数法、态密度计算与局部失稳判据推导

1.5.1 格林函数与DOS计算

为了实现定量研究，我们首先利用解析延拓技术计算复平面上的单粒子格林函数 $G(z)$：

态密度 $D(E)$ 则是格林函数在实轴上方极微扰动极限下的虚部：

在数值计算中，通常需要选取一个足够小但非零的能量截断 $\eta$（如 $\eta = 10^{-3} \sim 10^{-5}$）来进行三维布里渊区积分。一旦确定了态密度 $D(E)$，无相互作用的总粒子数密度 $\rho$ 和基态动能密度 $E_0$ 便可以通过对费米面下能态的严格积分得到：

1.5.2 基态稳定性分析与临界态密度 $D_{\text{crit}}$ 的推导

为了定量判断体系何时会发生向自旋极化铁磁态的相变，我们引入磁化强度 $P$ 作为变分序参量：

将这一变分关系代入到一阶低密度展开的能量表达式中，我们得到总能量关于自旋极化度 $P$ 的依赖关系 $E_{\text{tot}}(P)$：

现在，我们对序参量 $P$ 在未极化点 $P=0$ 进行二阶泰勒展开，以判断未极化状态（顺磁态）的局部稳定性：

利用关系式 $\frac{d e_F}{d\rho} = \frac{1}{D(e_F)}$，我们可以得到无相互作用能量关于密度的二阶导数：

对总能量函数求关于 $P$ 的二阶导数：

在未极化的顺磁极限点（即自旋对称点 $P=0$），$\rho_\uparrow = \rho_\downarrow = \frac{1}{2}\rho_{\text{tot}}$，费米能级简并为 $e_F(\frac{1}{2}\rho_{\text{tot}})$。二阶导数化简为：

根据变分原理，未极化顺磁基态稳定存在的前提条件是该二阶导数必须严格大于零。若二阶导数小于零，意味着顺磁态处于能量势垒的顶点（局部极大值），体系将自发极化以降低能量。因此，铁磁局部不稳定性发生的严格边界条件为：

这个判据在物理结构上极其优美。它完美重构了Stoner不稳定性，但将Stoner判据中粗糙的裸库仑排斥力 $U$ 严格地替换成了考虑了完整多体动态关联后的重整化散射长度项 $8\pi a$。这是对强关联电子学中巡游铁磁理论的一个极为坚实的理论提升。当且仅当对称状态下的单粒子态密度 $D(e_F)$ 超越临界阈值 $D_{\text{crit}}$ 时，顺磁态失稳，系统展现铁磁有序。若能带底部具有极高或发散的态密度（源于几何受挫带来的能带平坦化），该判据在极低密度极限下便能自发满足！

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

为了检验晶格几何受挫对稳定低密度巡游铁磁性的实际效能，论文精细设计并全面评估了以下四个关键三维晶格体系：

2.1 简单立方晶格（Simple Cubic Lattice）—— 零受挫的基准参照

• 物理特性：最简单、最基础的三维晶格，没有几何受挫。
• 色散关系：其紧束缚近似色散关系为： $$\epsilon(\mathbf{k}) = -2t(\cos k_x + \cos k_y + \cos k_z)$$ 取最邻近跃迁能 $t = 1$。其能带边界为 $[-6, 6]$，在带底 $e_{\text{min}} = -6$ 附近表现出标准的三维自由费米子各向同性行为，色散关系渐近为经典的 quadratic（$k^2$）形式。
• 计算结果：图1展示了其精确计算的DOS，其在带底非常平缓，最大DOS峰值仅在能带中心（由于范霍夫奇点），数值约为 $0.141$。而在低密度限下，带底附近的DOS极小。
• 铁磁稳定性分析：设置排斥能为物理上限 $U = 20$，代入两体散射计算可得，对应的临界态密度 $D_{\text{crit}} = 1/(8\pi a) \approx 0.17$。由于整个低密度填充区间内的真实态密度 $D(E)$ 远远小于这一临界阈值（$0.141 < 0.17$），因此简单立方晶格在整个物理填充范围内（直至 $\rho = 0.8$）均不展现任何铁磁有序（见图2）。若要强行使其产生铁磁性，裸相互作用必须达到非物理强度的 $U > 70$，这超出了绝大多数实际固体材料的物理区限。

2.2 次近邻跃迁立方晶格（Cubic Lattice with NNN Hoppings）—— 适度受挫的探索

• 物理特性：在简单立方晶格的基础上，引入沿着晶胞主轴方向跨越两个晶格常数的次近邻跃迁项 $t_2$（见图3）。这种额外的跃迁路径可以对动能进行调控。
• 色散关系： $$\epsilon_{\text{NNN}}(\mathbf{k}) = -2t(\cos k_x + \cos k_y + \cos k_z) - 2t_2(\cos 2k_x + \cos 2k_y + \cos 2k_z)$$
• 参数区间：$t_2 \in [0, -0.25t]$。当 $t_2 = 0.25t$ 时达到最大受挫。此时，能带底部（$\mathbf{k}=\mathbf{0}$ 附近）的低动量渐近行为由传统的 $k^2$ 陡然转变为 quartic（$k^4$）型扁平化色散关系。
• 计算结果：图4和图5分别定量描绘了 $t_2$ 变温到 $-0.25t$ 时，DOS向极大受挫极限演化的轨迹。在极限受挫 $t_2 = -0.25t$ 时，在带底 $E = -4.5$ 处，DOS不再像简单立方晶格那样趋于0，而是直接发散或保持在极大的有限值（图5中展现了带底的非零跃变）。
• 相变特性：对于 $U = 25$、$t_2 = -0.25$ 的强受挫极限，临界态密度 $D_{\text{crit}} = 0.128$。由于能带中部的范霍夫奇点发散，导致系统在较高填充密度 $\rho_c \approx 0.61$ 处发生一阶铁磁相变（First-order Phase Transition）。如图6所示，随着密度增高，体系的能量曲线 $E(P)$ 从单一极小值点（$P=0$）开始分裂出两个全新的非对称局部极小点（自旋极化态），最终在 $\rho_c \approx 0.61$ 时，这三个极小点能量简并，随后发生一阶跃变跳跃。虽然在此体系中成功诱导了巡游铁磁性，但由于相变临界密度 $\rho_c = 0.61$ 偏大，已经脱离了低密度展开理论最可靠的有效区间（严格低密度一般要求 $\rho \ll 1$）。

2.3 对角近邻跃迁立方晶格（Cubic Lattice with DNN Hoppings）—— 超低密度铁磁性的实现

• 物理特性：引入立方晶格面内对角线方向的最邻近跃迁 $t_3$（见图7）。这种受挫方式在物理上更具破坏性，能够极大地压制电荷的流动，从而压平能带。
• 色散关系： $$\epsilon_{\text{DNN}}(\mathbf{k}) = -2t(\cos k_x + \cos k_y + \cos k_z) - t_3 [\cos(k_x + k_y) + \cos(k_y + k_z) + \cos(k_z + k_x) + \cos(k_x - k_y) + \cos(k_y - k_z) + \cos(k_z - k_x)]$$
• 计算结果与相变：
  • 完全受挫极限（$t_3 = 0.25t$）：带底奇点极其剧烈，色散关系在整条高对称线上完全简并，DOS在带底呈现极具物理争议的发散（图9）。理论分析表明，在零密度极限下系统直接呈现饱和铁磁性，但随后在高密度下相变回顺磁（图11），表现出一种可能非物理的反常行为（这是由于带底非孤立极小值导致低密度展开失效）。
  • 近最大受挫调制（$t_3 = 0.245t$）：为了消除非物理发散，将其稍微偏离最大受挫点。此时，如图8和图10所示，在极低密度下，对称性自发破缺，体系发生二阶巡游铁磁相变。对于温和的排斥能 $U = 5$，在极低的电荷填充度 $\rho_c = 0.031$ 下即可建立稳固的基态铁磁长程序！由于该密度极其微小，高阶剩余项 $\mathcal{O}(\rho^{>2})$ 的贡献微乎其微，因此该体系中铁磁性的存在具有极高的物理和数学可信度，是本论文最具突破性的成果。

2.4 面心立方晶格（Face-Centered Cubic, FCC Lattice）—— 最具物理现实意义的体系

• 物理特性：面心立方（FCC）晶格是自然界中最常见的真实金属结构（如镍 Ni、铂 Pt 等的晶体结构）。它天然具有三维几何受挫特征（由四面体和八面体网格构成）。
• 色散关系（包含次邻近 $t_2$ 跃迁）： $$\epsilon_{\text{FCC}}(\mathbf{k}) = -4t(\cos k_x \cos k_y + \cos k_y \cos k_z + \cos k_z \cos k_x) - t_2(\cos 2k_x + \cos 2k_y + \cos 2k_z)$$
• 参数配置：取跃迁主项 $t = 1$，受挫强度 $t_2 = -0.25$（适度受挫）。此时，带底处的奇点形式为典型的 $D(E) \sim 1/\sqrt{|E - E_0|}$ 形式发散（图13）。
• 相变特性：如图14所示，在排斥能为温和的真实金属范围 $U = 5$ 时，系统在很低的填充密度 $\rho_c = 0.09$ 处触发了二阶铁磁相变。随着密度继续增加，基态自旋极化度不断攀升。这完美地解释了为什么面心立方金属镍（Ni）天然具有强大的巡游铁磁性，而简单立方结构的金属往往只能维持顺磁或反铁磁。这一计算不仅具有高度的学术价值，更与真实凝聚态材料物理有着直接的映射关联。

2.5 数据大汇总与基准对比表

为了直观展示上述所有测试体系的核心定量指标，我们在此重建并扩展了论文中的关键数据表格：

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

为了方便量子化学和关联电子物理研究人员快速复现论文中的核心数据（尤其是态密度 $D(E)$ 的精确数值积分和 $E_{\text{tot}}(P)$ 变分曲线的绘制），以下提供完整的复现指南和基于 Python 科学计算栈（NumPy + SciPy）的严格实现。

3.1 环境配置与依赖软件包

复现本计算不需要依赖复杂的专用量子化学大包，仅需标准的高性能计算科学计算库即可。建议在标准 conda 环境中安装：

conda install numpy scipy matplotlib tqdm

3.2 核心算法实现：以 NNN 立方晶格色散与严格格林函数法 DOS 计算为例

下面给出利用 Python 计算次近邻跃迁（NNN）晶格格林函数，并提取三维态密度 $D(E)$ 的高度优化算法代码。该代码通过引入小虚部常数 $\eta$，直接对三维第一布里渊区（1BZ）进行高能并行网格化求和。

import numpy as np
import matplotlib.pyplot as plt
from tqdm import tqdm

def get_dispersion_nnn(kx, ky, kz, t=1.0, t2=-0.25):
    """
    计算三维 NNN 晶格色散关系 epsilon(k)
    """
    cos_1 = np.cos(kx) + np.cos(ky) + np.cos(kz)
    cos_2 = np.cos(2*kx) + np.cos(2*ky) + np.cos(2*kz)
    return -2 * t * cos_1 - 2 * t2 * cos_2

def calculate_dos_grid(E_range, t=1.0, t2=-0.25, N_k=100, eta=0.05):
    """
    利用三维网格积分法计算态密度 D(E)
    """
    # 在第一布里渊区 [0, pi]^3 生成均匀网格
    # 利用反转对称性，只需计算 1/8 个布里渊区，节省内存并大幅加速
    k_space = np.linspace(0, np.pi, N_k, endpoint=False)
    KX, KY, KZ = np.meshgrid(k_space, k_space, k_space, indexing='ij')
    
    print("正在计算第一布里渊区上的单粒子能谱...")
    epsilon = get_dispersion_nnn(KX, KY, KZ, t, t2)
    epsilon_flat = epsilon.flatten()
    
    dos_results = []
    print("正在进行格林函数虚部数值积分解算...")
    # 遍历能量区间进行复格林函数求和
    for E in tqdm(E_range):
        # G(E + i*eta) = (1 / N) * \sum [1 / (E - epsilon + i*eta)]
        # D(E) = -1/pi * Im[G(E)]
        denominators = E - epsilon_flat + 1j * eta
        green_element = 1.0 / denominators
        g_avg = np.mean(green_element)
        dos = - (1.0 / np.pi) * np.imag(g_avg)
        dos_results.append(dos)
        
    return np.array(dos_results)

# 测试 NNN 极限受挫的态密度
if __name__ == "__main__":
    E_test = np.linspace(-5.0, 8.0, 300)
    # 计算 t2 = -0.25 时的受挫能带 DOS
    dos = calculate_dos_grid(E_test, t=1.0, t2=-0.25, N_k=120, eta=0.03)
    
    plt.figure(figsize=(8, 5))
    plt.plot(E_test, dos, 'r-', lw=2, label=r'$t_2 = -0.25$ (Max Frustrated)')
    plt.xlabel('Energy (E)', fontsize=12)
    plt.ylabel('Density of States D(E)', fontsize=12)
    plt.title('Replicated DOS for 3D NNN Lattice', fontsize=14)
    plt.grid(True, linestyle='--')
    plt.legend()
    plt.savefig('nnn_dos_replicated.png', dpi=300)
    plt.show()

3.3 变分能量 $E_{\text{tot}}(P)$ 的极小化与相变绘图脚本

在获得非零的 $D(E)$ 后，通过下述代码，研究者可以快速通过对自旋自洽磁化强度 $P$ 进行扫描，重构出类似图10和图14中基态自旋极化跃迁的物理曲线：

import scipy.integrate as integrate

def get_fermi_energy_and_E0(rho_target, E_grid, dos_grid):
    """
    给定目标粒子密度，通过对 DOS 积分反求费米能级 e_F 和非相互作用能量 E_0
    """
    # 计算累积积分以确定 N(e_F)
    cum_rho = []
    current_sum = 0.0
    dE = E_grid[1] - E_grid[0]
    for d in dos_grid:
        current_sum += d * dE
        cum_rho.append(current_sum)
    cum_rho = np.array(cum_rho)
    
    # 寻找到最接近目标强度的费米能级标引
    idx = np.abs(cum_rho - rho_target).argmin()
    e_F = E_grid[idx]
    
    # 计算无相互作用动能项 E_0 = \int e * D(e) de
    E_0 = np.sum(E_grid[:idx] * dos_grid[:idx] * dE)
    return e_F, E_0

def scan_energy_versus_magnetization(rho_tot, a_scat, E_grid, dos_grid, P_points=100):
    """
    计算特定自旋极化率 P 在 [-1, 1] 之间的总能量曲线 E_tot(P)
    """
    P_space = np.linspace(0.0, 1.0, P_points)
    E_total_curve = []
    
    for P in P_space:
        rho_up = 0.5 * rho_tot * (1.0 + P)
        rho_down = 0.5 * rho_tot * (1.0 - P)
        
        _, E0_up = get_fermi_energy_and_E0(rho_up, E_grid, dos_grid) if rho_up > 1e-8 else (0, 0)
        _, E0_down = get_fermi_energy_and_E0(rho_down, E_grid, dos_grid) if rho_down > 1e-8 else (0, 0)
        
        # 相互作用能贡献
        E_interaction = 8.0 * np.pi * a_scat * rho_up * rho_down
        E_tot = E0_up + E0_down + E_interaction
        E_total_curve.append(E_tot)
        
    return P_space, np.array(E_total_curve)

3.4 推荐的开源社区高级工具库

如果需要在大规模、多带（Multi-band）或非零温系统下进行更严格的自洽多体计算，本研究推荐以下社区顶级开源工具箱：

1. TRIQS (Toolbox for Research on Interacting Quantum Systems)：开源强关联物理大包。它自带高度优化的多带格林函数积分子、Dyson自洽求解器以及各种先进的量子蒙特卡洛（QMC）杂质求解器。
   • Link: https://triqs.github.io/triqs/
2. Wannier90：如果需要从第一性原理（DFT）计算得到的真实晶格哈密顿量构建低能Hubbard模型跃迁参数 $t_{ij}$，Wannier90 是构建极大限度定域化Wannier轨道的行业标准。
   • Link: http://www.wannier.org/
3. ALPS (Algorithms and Libraries for Physics Simulations)：提供了一系列用于求解强关联晶格模型的量子蒙特卡洛及精确对角化（ED）的标准算法包。
   • Link: https://alps.comp-phys.org/

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

4.1 关键里程碑文献回顾

本研究之所以能建立起如此扎实的物理和数学骨架，主要依赖于以下凝聚态物理及数学物理中的里程碑工作：

1. [2] E. H. Lieb, R. Seiringer, J. P. Solovej, Phys. Rev. A 71, 053605 (2005)
   • 贡献：严格证明了低密度极限下弱相互作用费米气体的能量展开式。这是连续介质中LSS定理的数学基石，为本工作奠定了最坚实的分析学基础。
2. [3] A. Giuliani, J. Math. Phys. 48, 023302 (2007)
   • 贡献：将LSS定理的能量上界严格推广到了三维格点Hubbard模型。这是建立格点低密度展开的里程碑之作。
3. [33] R. Seiringer and J. Yin, J. Stat. Phys. 133, 1139 (2008)
   • 贡献：完成了格点Hubbard模型低密度展开基态能量下界的严格证明，与Giuliani的上界一起，形成了格点二阶能量展开的闭环数学严密性证明。
4. [1] H. Tasaki, “Physics and Mathematics of Quantum Many-Body Systems”, Springer (2020)
   • 贡献：Tasaki教授是强关联物理学严密性证明的宗师。该书详细归纳了平带铁磁性（Flat-band Ferromagnetism）和Nagaoka极限的严格证明，是理解受挫巡游磁性不容错过的圣经。

4.2 局限性深度点评与建设性批判

尽管本工作成功证明了受挫晶格能显著稳定超低密度的巡游铁磁性，提供了一个极其优雅的分析框架，但从量子化学及多体关联物理的严格视角来看，该方案依然存在着以下不可忽视的物理与数学边界局限：

4.2.1 摄动展开截止一阶（高阶关联项的收敛性谜题）

本工作完全忽略了基态能量公式中高阶修正项 $\mathcal{O}(\rho^{>2})$ 的贡献。这是一个极大的物理赌注。在连续介质费米气体中，高阶项包含著名的 Lee-Huang-Yang (LHY) 修正，它倾向于增强对自旋极化的压制。尽管论文辩称，在 $\rho_c = 0.031$ 这样的超低密度下，高阶项在数值上非常微小，但在数学分析上，由于受挫带来的能带扁平化和发散的DOS，往往会导致低密度展开的摄动级数半径急剧收缩，甚至导致级数不收敛。 这一严密性缺失，使得“严格证明”一词在本工作的大部分格点上只能打上引号，仍属于高度严密的物理变分唯象学。

4.2.2 极大受挫极限下带底非孤立奇点导致两体散射理论崩溃

在完全受挫的对角近邻晶格（$t_3 = 0.25t$）和完全受挫的FCC晶格（$t_2 = -0.5t$）中，能带最底部的极小值点不再是空间中的孤立极值点，而是形成了一条连续的三维流形（简并线或简并面）。这直接导致了低能散射理论（s-wave scattering theory）中的散射长度 $a$ 的经典定义发生红外发散而彻底崩溃。 本文无法处理这一发散，只能通过微小偏离最大受挫参数（如 $t_3 = 0.245t$）来人工恢复孤立带底，从而规避这一科学难题。这极大阻碍了该方法被应用于如 Kagome 晶格或 Lieb 晶格等具有严格平带（Flat Band）的物理模型中，这也是本研究最为致命的数学软肋。

4.2.3 空间维度约束（$d=3$ 的机械锁定）

本研究所依赖的格点散射理论公式，其形式在维度上是极其脆弱的。在低维体系（一维 $d=1$ 和二维 $d=2$ 体系）中，低密度展开具有完全不同的数学结构（例如二维体系中的散射能具有 $\ln \rho$ 的对数修正形式，而一维则由于Bethe ansatz解析可解展现出完全不同的无磁性特征）。因此，本套优雅的格林函数局部失稳判据完全不适用于当前低维功能材料化学的热点领域（如单层过渡金属硫族化合物 TMDs，扭角双层石墨烯等）。

5. 补充：关联电子系统的现代设计与多体量子化学的交叉

5.1 受挫 Hubband 模型在分子轨道与金属有机框架（MOFs）设计中的启发

在现代配位化学与多孔材料设计中，如何合成常温分子铁磁体（Molecular Ferromagnets）或磁性金属有机框架（Magnetic MOFs）是一门备受瞩目的尖端前沿。传统的固体物理通常将物理晶格视为既定事实（如Ni的FCC晶格）。然而在量子化学中，研究人员可以通过精心挑选的配体和金属中心，在纳米尺度甚至宏观晶格中“人工组装”出具有特定几何受挫特征的晶格网络。

例如，通过设计三角形、鸽笼形（Kagome）或者具有次近邻轨道重叠的过渡金属配合物分子网络，我们可以在分子层面上重构类似于本文所述的对角近邻（DNN）跃迁路径。本文提供的物理洞察表明：

• 无需寻找极强的局域排斥金属（这通常伴随着重轨道带来的强自旋轨道耦合和复杂的局域退相干）；
• 只需通过精确调控配体场，使分子轨道重叠在低能区产生高度平坦的受挫能带结构。即便配体间的有效排斥能 $U$ 处于中等强度，体系也能在电荷轻微掺杂（低密度 $\rho \sim 0.03$）时，自发跨越临界态密度阻碍，建立高温分子巡游铁磁有序。

这为化学家设计新型高居里温度分子磁体提供了极具操作性的新策略：通过调控跃迁受挫度，寻找极低掺杂下的磁有序跃迁点，避免由于高密度电荷注入引发的电荷密度波（CDW）或超导（SC）等非磁性竞争相。

5.2 巡游铁磁性与非常规超导性的竞争与共存

在强关联物理学中，铁磁不稳定性（Ferromagnetic Instability）通常是不寻常超导配对（如p波三重态自旋超导）的催化剂。在受挫格点Hubbard模型中，由于带底附近巨大的态密度，当磁性受到外场或温度压制而未达到长程序时，强烈的自旋涨落（Spin Fluctuations）极易诱导电子克服库仑阻碍，配对为非常规超导体。本研究对低密度受挫铁磁区限的精确定位，也间接锁定了自旋涨落最剧烈的量子临界点（Quantum Critical Point）。这为探索自旋三重态超导（如在重费米子超导体 $\text{UTe}_2$ 或受挫铁磁体中）提供了极其关键的理论坐标定位，是连接凝聚态物理、先进化学材料合成与强关联多体理论的黄金纽带。