来源论文: https://arxiv.org/abs/2605.28761v1 生成时间: May 29, 2026 06:16

探秘二维强关联体系：普通 Van Hove 奇点处 Stoner 铁磁不稳定性的自发规避机制

0. 执行摘要

在凝聚态物理学与量子化学的前沿研究中，二维费米体系的非平庸拓扑与奇异色散结构一直处于核心地位。当体系的费米面（Fermi Surface, FS）与布里渊区（Brillouin Zone, BZ）边界切合于非平庸色散点——即鞍点（Saddle Point，通常称为 Van Hove 奇点 (Van Hove Singularity, VHS)）时，其单粒子态密度（Density of States, DOS）会出现对数发散。在经典的平均场理论（如 Stoner 铁磁不稳定性理论）框架下，这种态密度的发散预示着：对于任意微弱的排斥性哈伯德相互作用 $U > 0$，体系在有限温度 $T_{\text{MF}}$ 下都将不可避免地失稳并迈向铁磁长程序。

然而，这一理论预言与诸多精确的数值模拟（如量子蒙特卡洛模拟）以及实验现象（例如施加单轴压应力的 $\text{Sr}_2\text{RuO}_4$ 在逼近 Van Hove 奇点时并未出现铁磁性，反而展现出超导性）存在尖锐的冲突。为了阐明这一“规避 Stoner 铁磁不稳定性”的深层物理机制，I.S. Tupitsyn 等人在本工作中采用了两种前沿的**图规蒙特卡洛（Diagrammatic Monte Carlo, DiagMC）**方法：

1. 基于自洽一圈重整化和顶点校正的 Bold4+ 方案；
2. 零系统误差、可控截断的组合求和（Combinatorial Summation, CoS） DiagMC 精确计算。

研究表明，Stoner 铁磁不稳定性在极低温度（比平均场预测的 $T_{\text{MF}}$ 低一个数量级）下依然被成功规避。其背后的微观驱动力由双重机制构成：

• 有效相互作用的向下重整化（Effective Coupling Downward Renormalization）：在粒子-粒子（Cooper）通道中的多体散射占据主导，强力压制了粒子-空穴通道中促使铁磁失稳的顶点校正，导致有效相互作用强度 $U^*(T)$ 随温度降低而显著衰减；
• 准粒子残留值（Quasiparticle Residue, $Z$）的抑制与非费米液体（Non-Fermi Liquid, NFL）行为的自发涌现：自能 $\Sigma(\mathbf{Q}_{\text{VH}}, \omega_n)$ 的强频率依赖性导致准粒子权重 $Z \to 0$，静态极化泡（Polarization Bubble） $\Pi(T)$ 的对数发散被极大地削弱或截断，转化为一种非费米液体的幂律行为 $\text{Im}\Sigma \propto \omega_n^{1/2}$。

本工作不仅澄清了数十年来关于二维体系鞍点处 Stoner 铁磁性的理论争议，更为理解强关联过渡金属氧化物（如钌酸盐）和转角二维材料中的非常规超导性铺平了道路。

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

1.1 核心科学问题：平均场的失效与物理图像冲突

传统的 Stoner 铁磁转变是一种连续（二阶）相变，其发生标志为零动量下的静态自旋磁化率 $\chi(T)$ 发生发散。在平均场近似（Mean-Field Approximation, MFA）下，排斥作用为 $U$ 的费米气系统的静态磁化率满足：

其中 $\Pi_0(T) < 0$ 代表自由费米子的静态均匀极化泡（即朗道简并费米子的极化响应函数）。当温度降低至 $T_{\text{MF}}$，满足 Stoner 不稳定性条件：

系统便失稳，自发产生铁磁长程序。当费米面正好切过 ordinary Van Hove 奇点（记为 $\mathbf{Q}_{\text{VH}}$，在 2D 平方晶格中为 $(\pi, 0)$ 或 $(0, \pi)$ 等鞍点），态密度在费米能级 $E_F$ 附近呈对数发散：

（其中 $W$ 为带宽级截断能）。这导致极化泡 $\Pi_0(T)$ 在低温度极限下出现对数发散：

由该公式可推知，无论相互作用 $U > 0$ 多么微弱，只要温度足够低，Stoner 条件必然在某一有限温度 $T_{\text{MF}} > 0$ 被满足。然而，实验与精确计算一致否定了这一理论。核心问题在于：多体修正（即顶点校正与自能效应）究竟如何协同作用，使得系统自发避开了这一看似必然的铁磁失稳？

1.2 理论基础与先前模型（RG 图像及其局限）

为了超越平均场近似，先前的重整化群（RG）研究（如 Ojajärvi 等人）提出，粒子-粒子（pp）通道的散射过程会引起有效排斥相互作用 $U^*(T)$ 的显著向下重整化：

其中 $b \approx 1.88 > 1$ 是由于交叉单圈极化图贡献产生的几何系数。因为 $b > 1$，所以无论 $U$ 有多大，重整化后的耦合强度 $U^*$ 与极化泡的乘积满足：

这预言了 Stoner 不稳定性绝不会发生。然而，这一“RG 图像”存在几个致命的理论局限：

1. 它强行截断并孤立地考虑了特定的交叉图，忽视了其他通道顶点校正的竞争与耦合；
2. 它完全基于费米液体（Fermi Liquid, FL）假设，默认准粒子在 Van Hove 奇点处依然保持良好的相干性，未考虑准粒子权重 $Z$ 降低对极化泡的压制效应；
3. 缺乏无需近似的、数值精确的多体方法来定量验证该理论的正确性。

1.3 技术难点：多通道对数发散的非微扰处理

在 Van Hove 奇点附近，粒子-空穴（ph）通道（主导铁磁和电荷密度波失稳）与粒子-粒子（pp）通道（主导超导配对失稳）发生强烈的相互交织与量子干涉。这在技术上带来了极大挑战：

• 摄动理论的失效：由于态密度的对数发散，高阶费曼图的每一阶都携带对数项 $\ln^n(T)$，导致传统的微扰膨胀在低温度下完全不收敛；
• 常规 QMC 的局限性：非自洽的有限尺寸量子蒙特卡洛（DQC）面临严重的费米子符号问题（Sign Problem），难以在热力学极限下逼近极低温度和鞍点奇点；
• 多通道顶点相互作用的自洽闭合：如何在重整化中平等地对待所有的两体通道，并实现高精度的多圈图自洽计算，是物理学界的长期技术瓶颈。

1.4 本文方法细节

为了克服这些挑战，作者采用了两种互补的高精度图规蒙特卡洛（DiagMC）方案：

1.4.1 各向异性紧束缚哈密顿量

研究基于一个在平方晶格上定义的各向异性费米-哈伯德模型（Fermi-Hubbard Model）：

为了人为提高普通 Van Hove 奇点处的态密度，并使平均场铁磁相变温度 $T_{\text{MF}}$ 处于可控且易于数值观测的范围，作者设计了如下非对称的跳跃振幅（Hopping Amplitudes）：

• 最近邻跳跃：$t_x = t/2, t_y = t$
• 次近邻跳跃：$t'_x = 0.12t, t'_y = -0.24t$
• 单元能量取为 $t=1$。在该色散关系下，自由费米子的能带结构为： $$\epsilon_{\mathbf{k}} = -2 \sum_{\alpha=x,y} \left[ t_{\alpha} \cos k_{\alpha} - t'_{\alpha} \cos 2k_{\alpha} \right]$$ 该色散产生的费米面在 $y$ 方向上被高度“压平”（见图1），在鞍点 $\mathbf{Q}_{\text{VH}} = (\pi, 0)$ 处形成了非常显著的态密度峰值。相互作用强度设为 $U = 2t$。此时平均场 Stoner 临界温度大约为 $T_{\text{MF}} \approx 0.065t$。

1.4.2 Bold4+ 自洽图规方法

Bold4+ 是一种高度精密的自洽图规重整化方法。其核心思想是不仅对单粒子格林函数 $G$ 进行自洽戴森重整化（Dyson Renormalization），还对三个两体通道中的四点顶点函数（Four-Point Vertices）进行全自洽重整化：

• 粒子-空穴通道（Particle-Hole Channel, $ph$）
• 粒子-空穴交换通道（Particle-Hole Exchange Channel, $phx$）
• 粒子-粒子通道（Particle-Particle Channel, $pp$）

它完整包含了自能的所有直至 $U^4$ 阶的费曼图，并且通过几何级数的形式将特定类型的无穷阶图内嵌其中，有效地捕捉了不同物理通道之间的量子干涉与相互竞争。这种方法的优势在于计算效率极高，便于在大范围的温度与化学势参数空间中进行系统扫描。

1.4.3 DiagMC-CoS 精确图规方法

为了彻底摆脱近似，作者引入了基于**组合求和（Combinatorial Summation, CoS）**算法的数值精确图规蒙特卡洛技术。其数学基础是将自能 $\Sigma$ 和静态磁化率 $\chi$ 展开为关于 $U$ 的泰勒级数：

CoS 算法可以高效、并行地在热力学极限下计算第 $n$ 阶所有连通费曼图的代数和 $a_n$（在本项目中截断阶数达到极高的 $N = 8$）。

由于直接级数在低温度下由于对数项的存在而可能发散，作者使用了以下两种数学外推技术对级数进行高精度重构：

1. Dlog Padé 逼近器：对对数导数 $\frac{d}{dU}\ln \chi(U)$ 构造有理分式逼近，能够精确捕捉函数底部的物理奇点；
2. 积分逼近器（Integral Approximants）：利用更广泛的微分方程逼近序列，通过不同逼近器的收敛展开范围，提供对系统外推误差（Systematic Extrapolation Error）的定量、可靠的估计。

2. 关键 Benchmark 体系、数值计算数据与物理发现

2.1 Benchmark 体系构型与费米面行为

为了进行精确的 Benchmark 评估，研究重点考察了以下两种状态下的物理行为（如图1所示）：

1. Van Hove 点共振状态：化学势取为 $\mu = \mu_{\text{VH}} = -0.110t$（对应 $T/t=0.01$ 时相互作用体系的数值），费米面（FS）恰好通过 $\mathbf{Q}_{\text{VH}} = (\pi, 0)$ 点（即图 1 中的 Point A）；
2. 远离 Van Hove 点状态：化学势偏离 $\mu = \mu_{\text{VH}} + 0.2t$，费米面表现为开放结构并与布里渊区边界相交于非鞍点（即图 1 中的 Point B）。

通过计算 $Re G^{-1}(\omega \to 0, \mathbf{k}) = 0$，研究发现即使在 $U=2t$ 的强相互作用下，相互作用诱导的费米面畸变也极其微弱（图 1 中实线与虚线几乎重合）。这表明我们可以安全地忽略由于费米面拓扑结构变化（如 Lifshitz 相变）对 Stoner 不稳定性消失的影响，从而将物理注意力完全集中在纯多体动力学重整化机制上。

2.2 Bold4+ 方案计算数据深度剖析

2.2.1 磁化率与极化泡的温度演化

图 2 展示了 Bold4+ 方案在 $\mu = \mu_{\text{VH}}$ 下计算得到的关键物理量随温度的演化特征：

• 静态均匀磁化率 $\chi(T)$（图 2 左面板）：当温度 $T$ 从平均场转变温度 $T_{\text{MF}} \approx 0.065t$ 降至极低温度 $T = 0.01t$ 时，$\chi(T)$ 呈现单调上升，但显然趋于平缓并最终饱和在约 1.5 左右的有限值，没有任何发散迹象；
• 重整化极化泡 $\Pi(T)$（图 2 中面板）：自由费米子的极化泡 $\Pi_0(T)$ 随着温度降低而呈现出对数发散的特征（黑色虚线）；然而，Bold4+ 中计算出的自洽极化泡 $\Pi(T)$（蓝色虚线点）在极低温度下发生了完全的平流饱和。这一极其关键的物理行为，直接堵死了 Stoner 条件被满足的可能；
• 有效相互作用 $U^*(T)/U$（图 2 右面板）：利用关系式： $$\chi(T) = -\frac{\Pi(T)}{1 + U^*(T)\Pi(T)}$$ 提取出的 $U^*(T)$ 随温度降低显著下降。在 $T \to 0$ 极限下，其行为可极佳地被基于自洽极化泡的代数公式所描述： $$U^*(T)/U = [1 - U^* b \Pi(T)]^{-1} \implies U^*(T) \propto \Pi(T)^{-1/2}$$ 由于 $\Pi(T)$ 趋于饱和，有效相互作用也同样趋于饱和。

2.2.2 准粒子相干性的丧失：自能的非费米液体特征

为了理解极化泡 $\Pi(T)$ 的饱和机制，必须审视体系的自能性质。图 3 绘出了自能虚部在松原频率（Matsubara Frequency）下的演化规律：

• 远离 VHS 点（Point B）：$\text{Im}\Sigma(\mathbf{k}_B, \omega_n)$ 展现出完美的朗道费米液体（Fermi Liquid）特征，即在低频下严格线性依赖于松原频率： $$\text{Im}\Sigma(\omega_n) = -\alpha \omega_n$$
• 处于 VHS 点（Point A）：在极低温度下，这一线性规律被彻底颠覆，自能虚部在 $\omega_n \to 0$ 极限下直接饱和至一个有限的非零常数： $$\text{Im}\Sigma(\mathbf{Q}_{\text{VH}}, \omega_n \to 0) \approx \text{const}$$ 这意味着准粒子的寿命在低能下被强烈压制，准粒子残留值 $Z(\mathbf{Q}_{\text{VH}}) \to 0$。正是由于准粒子相干性的完全丧失，极化泡中两个格林函数的褶积积分被强烈平滑，导致 $\Pi(T)$ 失去了对数发散的动力。

2.3 精确 DiagMC-CoS 方案计算数据深度剖析

Bold4+ 虽然成功解释了 Stoner 不稳定性的规避，但其自能饱和的物理结论带有自洽一圈近似的伪影。接下来，通过数值精确的 DiagMC-CoS 计算（图 4），作者揭示了更为精细、真实的物理图景：

2.3.1 极化泡与有效相互作用的非饱和协同演化

• 静态磁化率 $\chi(T)$（图 4 左面板）：精确计算的 $\chi(T)$（红色点）与 Bold4+ 的定性结果高度吻合，它在极低温度下依然平滑、不发散。在 $T = 0.02t$ 时，$\chi$ 约等于 1.25，比平均场预言的无穷大表现得极为驯服；
• 极化泡 $\Pi(T)$ 的非饱和增长（图 4 中面板）：与 Bold4+ 的饱和预测不同，精确的自洽极化泡 $\Pi(T)$（红色虚线）并未真正饱和，而是继续随着温度降低而增长。但极具启发性的是，其增长斜率相比于自由极化泡 $\Pi_0(T)$ 被大幅度削弱；
• 有效相互作用 $U^*(T)$ 的持续下降（图 4 右面板）：由于 $\Pi(T)$ 还在继续增长，精确提取出的有效相互作用 $U^*(T)/U$ 呈现出单调且未饱和的急剧下降趋势。在最低可达温度下，它的值已从高温下的 1.0 崩塌至 0.52 附近。这意味着多体关联效应将排斥作用削弱了将近一半！

2.3.2 $\omega_n^{1/2}$ 非费米液体幂律的发现

为了揭示为什么精确计算下的 $\Pi(T)$ 会缓慢增长而没有绝对饱和，作者在图 5 中细致分析了 DiagMC-CoS 计算得到的自能虚部频赖关系： 在 $T/t = 0.02$ 且 $\mu = \mu_{\text{VH}}$ 下，$\text{Im}\Sigma(\mathbf{Q}_{\text{VH}}, \omega_n)$ 在小频率范围内不服从费米液体的线性关系，也并非 Bold4+ 预测的常数饱和关系，而是完美契合一个分数幂指数为 $1/2$ 的非费米液体行为：

这一非费米液体行为具有重大理论价值：

1. 由于自能表现为 $\omega_n^{1/2}$，准粒子权重在低能下表现出非简谐的发散： $$Z(\omega_n) \sim [1 - \text{Im}\Sigma/\omega_n]^{-1} \propto \omega_n^{1/2} \to 0$$
2. 当我们将这种非费米液体的格林函数代入极化泡积分中时，由于积分核中存在 $\omega_n^{1/2}$ 的压制，导致原本的对数发散 $\ln(1/T)$ 被软化为一种发散速度极慢的幂律行为。这完满地解释了为什么 $\Pi(T)$ 既不绝对饱和、但又没有导致 Stoner 发散的精细多体机制。

2.4 磁化率的动量空间解析（Fig. 6）与短程铁磁关联

在图 6 中，作者计算了在温度 $T/t = 0.02$ 下，磁化率 $\chi(\mathbf{Q})$ 沿着高对称路径 $(0,0) \to (\pi, 0) \to (\pi, \pi) \to (0,0)$ 的空间分布：

• 短程涨落的主导性：无论在 Bold4+ 还是 CoS 方案中，$\chi(\mathbf{Q})$ 在全动量空间中均表现出非常有限的数值（最高不超过 1.8）。这表明整个体系处于一种极其温和的、仅存在**短程铁磁关联（Short-Range Ferromagnetic Correlations）**的顺磁状态；
• 峰值的非平里位移：尤为惊人的发现是，对于精确的 CoS 计算，静态磁化率的最大值并不是位于代表均匀铁磁性的零动量 $\mathbf{Q}=(0,0)$ 处，而是自发向有限动量方向发生了移动（Peak Shift）。这一非平庸的峰值位移现象，强烈暗示着在低能极限下，体系通过非单调的多体极化响应，自发规避了不稳定的长程铁磁相变，转而表现出非平庸的非均匀电荷/自旋关联特征。

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

由于本文所使用的自洽 Bold4+ 代码和组合求和（CoS）图规蒙特卡洛（DiagMC）算法属于前沿多体计算凝聚态领域的内部专业级代码，本节基于文献 [35-38] 公开的物理逻辑和通用的图规蒙特卡洛技术栈，提供一份面向技术复现的代码架构及工程化设计指南。

3.1 图规蒙特卡洛（DiagMC）的通用数学框架与采样算法

在图规蒙特卡洛中，我们对物理量（如自能 $\Sigma(\mathbf{k}, \tau)$）的费曼图空间进行随机采样。其核心积分表达式为：

其中 $\mathcal{G}_n$ 为 $n$ 阶费曼图的拓扑结构，$\mathcal{D}$ 为该图对应的贡献权重。蒙特卡洛状态空间由三元组组成：$(\mathcal{G}_n, n, \mathcal{X})$，其中 $\mathcal{X}$ 包含了所有的内部环路动量、松原时间等连续变量。为了对该空间进行无偏采样，必须构建满足细致平衡（Detailed Balance）的马尔可夫链更新算子：

1. 阶数增加更新 (Order-Increase Update: $n \to n+1$)： 随机插入一个新的相互作用顶点 $U$，并产生两个新的内部虚线，重新连接传播子。其接受概率满足 Metorpolis-Hastings 准则： $$P_{\text{acc}} = \min\left( 1, \frac{\mathcal{W}(\mathcal{G}_{n+1})}{\mathcal{W}(\mathcal{G}_n)} \frac{p_{\text{prop}}(n \to n+1)}{p_{\text{prop}}(n+1 \to n)} \right)$$
2. 阶数减少更新 (Order-Decrease Update: $n \to n-1$)： 随机删除一个无外线的自能顶点，并重新桥接对应的格林函数线。
3. 动量/时间更新 (Internal Variable Update)： 利用高斯分布或者局域非均匀分布，对某个内部松原时间 $\tau_j$ 或环路动量 $\mathbf{k}_j$ 进行局部微扰，实现快速的相空间混合。

3.2 组合求和（CoS）核心算法的设计模式（Python/C++ 伪代码）

组合求和算法（Combinatorial Summation, CoS）能通过统一生成某阶所有拓扑图的求和来规避对阶乘级（$n!$）个图分别采样的阶乘壁垒。以下给出一个高度简化的 CoS 算法的设计框架：

import numpy as np
from scipy.integrate import nquad

class DiagrammaticCoS:
    def __init__(self, order_limit=8, k_mesh=64, beta=100.0):
        self.N_max = order_limit
        self.k_mesh = k_mesh
        self.beta = beta  # 倒数温度 1/T
        
    def bare_green_function(self, k, tau):
        """
        输入动量 k (2D 向量) 与虚时间 tau，返回自由格林函数 G_0(k, tau)
        含有非对称紧束缚色散关系
        """
        tx, ty = 0.5, 1.0
        tx_prime, ty_prime = 0.12, -0.24
        # 色散关系计算
        epsilon_k = -2.0 * (tx * np.cos(k[0]) - tx_prime * np.cos(2*k[0]) + 
                            ty * np.cos(k[1]) - ty_prime * np.cos(2*k[1]))
        mu_vh = -0.110  # 对应 U=2t 下的 VH 化学势
        xi_k = epsilon_k - mu_vh
        
        # 费米-狄拉克分布与时间演化
        if xi_k > 0:
            return -np.exp(-xi_k * tau) / (1.0 + np.exp(-xi_k * self.beta))
        else:
            return -np.exp(xi_k * (self.beta - tau)) / (1.0 + np.exp(xi_k * self.beta))

    def generate_combinatorial_kernels(self, order):
        """
        核心算法：生成第 n 阶的所有连通图拓扑结构组合算子
        利用分块矩阵行列式或拓扑树，消除重复路径，降低 N! 复杂度
        """
        # 此处在 C++ 实际代码中会通过位掩码（Bitmask）与自同构群（Automorphism Group）生成器实现
        topologies = [] 
        # 示例：返回生成的拓扑函数列表
        return topologies

    def evaluate_coefficient_an(self, n):
        """
        蒙特卡洛积分求解 a_n 的核心循环
        """
        coef_accumulator = 0.0
        # 采样配置：虚时间 tau_1...tau_n， 动量 q_1...q_n
        # 采用多维马尔可夫链 Monte Carlo (MCMC) 进行高维球积分
        return coef_accumulator

3.3 开源工具包与计算生态推荐

目前，在国际强关联多体物理计算社区中，虽然 Tupitsyn 组的 CoS 代码没有完全开源，但有一系列非常成熟的开源框架可以用来复现本文所涉及的物理量计算：

1. TRIQS (Toolbox for Research on Interacting Quantum Systems)

   • 功能：一个非常强大的量子多体计算 C++/Python 框架，拥有极其成熟的自能、格林函数格点映射和松原频率处理库。
   • 链接：https://triqs.github.io
   • 复现路径：可以使用 TRIQS 内置的紧束缚模型模块定义本文的各项异性色散，并调用其强大的松原频率变换工具，以其一圈/二圈自洽循环作为 Bold4+ 的基准起点。

2. ALPS (Algorithms and Libraries for Physics Simulations)

   • 功能：经典的强关联蒙特卡洛计算包，提供高效的并行采样框架和误差自相关时间分析工具。
   • 链接：https://alps.comp-phys.org

3. Simons Collaboration on the Many-Electron Problem 开源 Repo

   • 功能：西蒙斯电子多体协作组在 GitHub 上维护了大量先进的多体计算代码，包括部分简化版的图规蒙特卡洛和极化泡求解器。
   • 链接：https://github.com/simons-collaboration

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

4.1 关键文献引用及导读

1. [1] E. C. Stoner, Proc. R. Soc. Lond. Ser. Math. Phys. Sci. 169, 339 (1939).
   • 重要性：奠定了巡游电子铁磁性（Stoner Ferromagnetism）的经典平均场框架。
2. [30] C. W. Hicks et al., Science 344, 283 (2014) & [31] M. E. Barber et al., Phys. Rev. Lett. 120, 076602 (2018).
   • 重要性：通过对 $\text{Sr}_2\text{RuO}_4$ 施加单轴压应力逼近 Van Hove 奇点的关键实验，证实系统并未出现平均场预言的铁磁相变，而是形成了非常规超导。这是本工作最直接的实验源泉。
3. [34] R. Ojajärvi, A. V. Chubukov, Y.-C. Lee, M. Garst, and J. Schmalian, NJP Quantum Materials 9, 105 (2024).
   • 重要性：提出了粒子-粒子散射通道对有效相互作用进行对数重整化的理论，是本工作所批判和超越的直接理论靶点。
4. [38] E. Kozik, Nature Communications 15, 7916 (2024).
   • 重要性：详细阐述了组合求和（CoS）算法的数学基础，是本工作精确计算的技术基石。

4.2 工作局限性深度评论

尽管本工作是一项具有里程碑意义的研究，但理性的技术作者应当指出其客观存在的几点物理局限性：

4.2.1 阶数截止限制与中等相互作用区间的局限

本工作通过 DiagMC-CoS 方法计算到了 $N=8$ 阶。对于 $U=2t$ 这一弱到中等相互作用区间，利用 Dlog Padé 进行外推是极其可靠的。然而，如果强关联效应进一步增强（如 $U \ge 4t$，这在真实的铜氧化物超导体中很常见），高阶费曼图的级数展开由于非微扰效应将变得极度不稳定。在这类极强关联区间下，有限阶数的重构往往无法保证收敛，这是所有图规蒙特卡洛方法的固有物理软肋。

4.2.2 高阶 Van Hove 奇点（Higher-Order VHS）的普适性尚未验证

本文的研究基于“普通” Van Hove 奇点（即色散关系呈二次鞍点，态密度呈对数发散 $\ln(1/E)$）。在当今热门的转角二维体系（如扭角双层石墨烯、转角过渡金属硫族化合物 TMDs）中，通过微调转角可以实现高阶 Van Hove 奇点（HOVHS）。在 HOVHS 处，能带呈现出更平坦的四次或六次幂指数鞍点，态密度呈现更为剧烈的幂律发散（如 $E^{-1/3}$ 或 $E^{-1/2}$）。在如此剧烈的态密度发散下，有效相互作用的重整化速度能否依然超越这一发散速度、自发规避 Stoner 铁磁性？本工作的方法尚未对此给出明确的物理外推，亟待后续研究跟进。

4.2.3 竞争有序相（超导与电荷密度波）的显式构建缺失

本工作把研究重心完全放在了阐明“为什么铁磁不稳定性会被规避”上，并利用静态磁化率的有限值来论证其正确性。然而，强关联体系在低温度下真正的物理宿命通常是超导态（Superconductivity）或电荷密度波（CDW）。本工作并未在计算体系中显式引入超导配对序参量，亦未能完整绘制出在极低温度下竞争相（如 $p$ 波或 $d$ 波非常规超导）的相图。这使得该工作更像是一个“排除性”的研究，而非一个完整的“建构性”理论。

5. 补充讨论：强关联二维体系中的非费米液体物理学

5.1 从 Landau 费米液体理论的失效谈起

Landau 费米液体理论是整个现代金属物理的基石。其基本假设是：在中微弱相互作用下，多体系统可以通过绝热演化映射为无相互作用的准粒子（Quasiparticles）。这些准粒子继承了自由电子的大部分电荷与自旋量子数，只是质量被重整化，且其寿命 $\tau \propto T^{-2}$，在低温极限下具有近乎无限的稳定性，因而可以被视为良好的单粒子。其自能虚部满足：

这意味着当频率 $\omega_n$ 趋近于零时，自能消失的速度比频率本身还要快，保证了准粒子的相干性：

然而，在本工作的 Van Hove 奇点处，这一基石彻底瓦解。精密的 DiagMC-CoS 数值模拟揭示：

这不仅意味着准粒子权重 $Z$ 在费米面上彻底耗尽（$Z \to 0$），更意味着体系跨入了经典的“非费米液体”或“边缘费米液体”（Marginal Fermi Liquid）范畴。在这一状态下，电荷的非相干散射主导了输运过程。在实验上，这会直接导致一系列非平庸的物理效应：

• 线性电阻行为：由于散射率 $1/\tau \propto T^{\gamma}$（$\gamma < 2$），体系的电阻率不再遵循经典的 $T^2$ 关系，而是呈现出奇妙的 $T^{\gamma}$（甚至线性 $T$）依赖行为；
• 维德曼-弗兰茨定律（Wiedemann-Franz Law）的失效：非相干的电荷输运与热量输运发生脱耦，导致洛伦兹数 $L = \kappa / (\sigma T)$ 偏离经典的 Sommerfeld 极限值。

5.2 粒子-粒子与粒子-空穴通道竞争的物理本质

在强关联物理中，一维系统具有完美的“电荷-自旋分离”和多通道完全干涉特性，可以使用 parquet RG 进行解析描述。而在普通的二维系统（远离 VHS 点）中，粒子-空穴通道与粒子-粒子通道通常是弱耦合的：粒子-粒子通道在零总动量处由于 Cooper instability 产生发散（导致超导），而粒子-空穴通道在有限转移能量处主导电荷/自旋密度波。

然而，Van Hove 奇点的存在，巧妙地将一维的“多通道深度交织”物理注入到了二维体系之中。由于费米面在鞍点处呈高度平坦化，使得粒子-空穴散射和粒子-粒子散射在相同的对数标度上竞争：

• ph 通道试图通过自旋极化来降低动能（Stoner 铁磁）；
• pp 通道则试图利用排斥相互作用下的动量空间符号交替来诱导非常规配对（如 $d$-wave 或 $p$-wave 超导）。

本工作最伟大的物理贡献，在于利用数值上毫无偏置、精确无误的图规蒙特卡洛，毫无争议地表明了：在这个多通道量子干涉的角力场上，粒子-粒子通道的多体散射笑到了最后。它通过强力的顶点相互作用重整化，将排斥相互作用 $U$ 拽离了铁磁失稳的深渊，同时准粒子权重的丧失进一步软化了铁磁响应。这一精妙的多体物理防线，使得系统自发地避开了铁磁性，并将多体能量完整地保留给了更具非凡魅力的基态——非常规超导。这也正是 strained $\text{Sr}_2\text{RuO}_4$ 在高压应变逼近 Van Hove 点时，铁磁不升反降、而超导温度 $T_c$ 出现急剧攀升的极深层次、最底层的微观物理机制。