来源论文: https://arxiv.org/abs/2607.01519v1 生成时间: Jul 03, 2026 18:19

非微扰强场非线性霍尔效应与非平衡稳态：基于NEGF-DMFT的统一理论架构

0. 执行摘要

非线性霍尔效应（Nonlinear Hall Effect, NLHE）作为凝聚态物理与量子材料领域的前沿课题，建立了输运行为与波函数几何性质（如贝里曲率偶极子，Berry Curvature Dipole, BCD）之间的直接桥梁。然而，现有的理论框架多局限于弱场微扰机制，且高度依赖常数弛豫时间近似（Relaxation-Time Approximation, RTA）。在强电场驱动或拓扑热点邻近费米面的非微扰机制中，电场诱导的频域/时域重新分布、带间相干性以及电子-电子、电子-声子多体相互作用将导致微扰图景彻底失效。

本研究详细推导并建立了一套基于**非平衡格林函数（Nonequilibrium Green’s Functions, NEGF）与动力学平均场理论（Dynamical Mean-Field Theory, DMFT）**相结合的非微扰稳态输运理论框架。该框架采用规范协变维格纳（Gauge-Covariant Wigner）表符，通过 Moyal 星积非微扰地引入恒定直流电场，并在时域引入共动动量坐标系以消除漂移项。该方法成功应用于双带拓扑半金属模型，揭示了强场下非平衡分布边缘扫过贝里曲率“热点”（Berry Curvature Hot Spots）时的非微扰交叉（Crossover）行为，阐明了带间相干性在无费米面体系（如电中性 Weyl 半金属）中主导非线性霍尔响应的量子机制。更重要的是，通过耦合自洽 DMFT 自能，本工作发现电子-声子耦合及本地 Hubbard 相互作用的实部自能会显著重整化带结构，打开能隙并强化贝里曲率热点，从而反直觉地增强非线性霍尔信号。本理论框架为定量模拟非平衡、强关联拓扑材料的非线性输运提供了统一且严格的微观微观工具，也为强电场输运与超快光谱实验提供了明确的物理指引。

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

1.1 核心科学问题：微扰极限的失效与强场非微扰机制

在时间反演对称但空间反演对称性破缺的系统中，纵向直流电场可以在二阶响应下驱动横向电流，即非线性霍尔效应。在弱场微扰极限下，该效应通常用贝里曲率偶极子（BCD）来描述：

其中 $f_0(\mathbf{k})$ 为平衡态费米-狄拉克分布函数，$\Omega_c(\mathbf{k})$ 为贝里曲率。这一经典图像暗含了两个假设：

1. 微扰响应：电场仅对平衡态分布函数造成微小的线性漂移（$\delta f \propto E$），不改变系统本身的准粒子能带结构和波函数几何属性。
2. 常数弛豫时间（RTA）：所有的耗散与散射过程均可由一个与动量、能量无关的常数弛豫时间 $\tau$ 简化表征。

然而，当系统处于以下物理区间时，上述微扰图像将彻底失效：

• 强电场驱动：电场驱动导致的动量漂移距离 $\Delta k \approx eE\tau/\hbar$ 达到了动量空间中贝里曲率剧烈变化（即“热点”）的尺度。此时，分布函数的形变极度非线性，无法用一阶导数（偶极子）外推。
• 带间相干与多能带几何：在狄拉克或外尔半金属中，费米面极其接近能带接触点。强电场会直接激发带间跃迁（如 Schwinger 机制），产生强烈的带间相干性。传统的半经典玻尔兹曼方程无法自洽处理带间相干性对横向输运的贡献。
• 强关联与自能重整化：电子-电子散射及电子-声子相互作用不仅提供耗散，还会引入强烈的频率依赖性自能（Self-energy），从而改写谱函数与量子几何张量。

1.2 理论基础：Keldysh 形式论与规范协变维格纳表征

为了自洽处理强场、耗散和多体关联，必须诉诸非平衡格林函数（NEGF）形式论。在非平衡稳态（Nonequilibrium Steady State, NESS）下，系统的时间平移对称性保持，但空间均匀受驱体系存在非零直流电场 $\mathbf{E}$。此时，格林函数在实空间和时域中不再具有简单的差值形式。

为了恢复平移对称性并直接引入电场，本研究采用了规范协变维格纳表征（Gauge-Covariant Wigner Representation）。对于给定的二时两点格林函数 $G(x_1, x_2)$，引入中心坐标 $X = (x_1 + x_2)/2$ 和相对坐标 $x = x_1 - x_2$。在存在电磁矢量势 $\mathbf{A}(X)$ 的情况下，定义规范协变动量（动力学动量）变量 $\kappa$：

在维格纳空间中，算符的乘积（卷积）被映射为著名的 Moyal 星积（Moyal Star Product, $\star$）。对于任意两个维格纳函数 $A(X, \kappa)$ 和 $B(X, \kappa)$，其星积定义为：

当系统处于空间均匀且电场恒定（恒定直流电场 $\mathbf{E}$ 沿 $x$ 方向）的稳态时，所有的物理量均与中心坐标 $X$ 无关。此时 Moyal 星积退化为仅混淆频率 $\omega$ 与动力学动量 $\boldsymbol{\kappa}$ 导数的特殊形式：

此处的指数算符在数学上处理起来极为困难，因为它将能量域与动量域深度耦合，是非平衡强场物理的核心难点。

1.3 技术难点与时域共动动量变换

为了在数值上稳健地求解上述非平衡 Dyson 方程，直接在频域对 Moyal 星积进行离散化会导致极其庞大的稠密矩阵，且在强场下收敛性极差。作者提出了一项关键的数学突破：将 Dyson 方程转换至时域（Relative Time $\tau$），并通过引入“共动动量坐标”（Comoving Momentum）消除电场驱动项。

具体而言，在时域中，维格纳表示下的二阶微分 Moyal 乘积退化为一阶导数项。在常数阻尼近似（常数退迟自能 $\Sigma^R = -i\Gamma$）下，退迟格林函数 $G^R(\tau, \boldsymbol{\kappa})$ 满足的运动方程为：

为了消除其中的漂移项 $\frac{i}{2} \mathbf{E}\cdot \nabla_{\boldsymbol{\kappa}}$，引入随电场漂移的共动动量：

利用链式法则，可以证明在该变换下，时域微分算符发生完美消去：

这使得非平衡 Dyson 方程在共动坐标系下简化为一个普通的一阶常微分方程：

其严格解析解可直接积分积出：

该解是非非微扰强场输运理论的基石。它不仅非微扰地包纳了任意强度的直流电场 $\mathbf{E}$，还自洽地考虑了动量空间路径积分带来的相位积累（即著名的 Bloch 振荡与 Wannier-Stark 局域化效应的隐式体现）。

1.4 方法细节：非平衡密度矩阵与多体自能自洽解法

有了退迟格林函数 $G^R$ 后，为了计算输运电流，必须求解 Keldysh 方程以获取少于格林函数（Lesser Green’s Function） $G^<(\tau, \boldsymbol{\kappa})$。在NESS稳态下：

其中 $\Sigma^<$ 为系统与外部零温（或有限温）费米子库（Fermionic Bath）耦合带来的自能，$\Sigma^<(\tau) = i\Gamma \int \frac{d\omega}{2\pi} e^{-i\omega\tau} f(\omega)$。在时域中执行两次 Moyal 星积后，可以严格推导出等时少于格林函数（即非平衡等时密度矩阵） $\tilde{n}(\boldsymbol{\kappa}) = \text{Im } G^<(0, \boldsymbol{\kappa})$ 的表达式：

定义有效谱振幅（Effective Spectral Amplitude） $\mathcal{G}(\omega, \boldsymbol{\kappa})$：

则非平衡动量分布可写为极具物理启发性的紧凑形式：

最后，利用速度算符 $v_y(\boldsymbol{\kappa}) = \partial H(\boldsymbol{\kappa})/\partial \kappa_y$，横向非线性霍尔电流可不依赖任何半经典近似，直接通过全量子力学迹（Trace）计算：

当考虑相互作用（如电子-声子或电子-电子关联）时，自能 $\Sigma$ 变为频率依赖。此时，作者建立了一套迭代自洽解法：

1. 视相互作用自能 $\Delta\Sigma^R$ 为微扰，重写 Dyson 方程为： $$G^R(\tau, \boldsymbol{\kappa}) = G^R_0(\tau, \boldsymbol{\kappa}) + \left( G^R_0 \star \Delta\Sigma^R \star G^R \right)(\tau, \boldsymbol{\kappa})$$
2. 通过时域傅里叶变换，将上述 Moyal 星积投影到沿电场方向 $\kappa_\parallel$ 的一维谐波（Harmonics）通道上，实现一维因果卷积的快速傅里叶变换（FFT）加速： $$M_n(\tau, \boldsymbol{\kappa}_\perp) = \int_0^\tau d\tau_1 K_n(\tau_1) G^R_n(\tau - \tau_1, \boldsymbol{\kappa}_\perp)$$
3. 结合动力学平均场理论（DMFT），在局域近似下提取局部格林函数 $G_{loc}(\omega)$，计算多体自能（电子-声子采用 Migdal 近似，电子-电子相互作用采用自洽二阶微扰理论），反馈回格林函数求解器中，直至自洽收敛。

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

2.1 物理模型：反演对称破缺的双带 Weyl 半金属

本研究采用的 benchmark 体系为具有两个 Weyl 类似（Weyl-like）点的最低限度双带晶格模型，其哈密顿量为：

选用参数为 $D_x = 0.5$，$D_y = 1.0$。该模型描述了一个在费米能量 $\epsilon = 0$ 附近拥有两个 Weyl 接触点的半金属系统。由于该系统具有时间反演对称性，其贝里曲率在第一布里渊区内关于动量 $\mathbf{k}$ 为奇函数（$\Omega(-\mathbf{k}) = -\Omega(\mathbf{k})$），因此全区积分恒为零。然而，反演对称性的破缺（$D_x \neq 0$）使得贝里曲率的偶极子矩不为零，从而允许二阶非线性霍尔效应发生。

2.2 关键计算数据分析

2.2.1 强场非微扰交叉（Crossover）验证

为了系统研究强电场导致的非微扰效应，作者计算了不同化学势 $\mu$ 和阻尼速率 $\Gamma$ 下，场强归一化电流 $J_y/E$ 随电场 $E$ 的演化规律（见原文 Fig. 1 和 Fig. 2）。

• 弱场线性区（微扰区）：在极小电场下，所有曲线的 $J_y/E$ 随 $E$ 均呈线性增长，这对应着 $J_y \propto E^2$ 的经典微扰 scaling。这一区域可由贝里曲率偶极子（BCD）完美描述。
• 非微扰漂移交叉（Drift Crossover）：随着电场增强，在特定的交叉电场 $E^*$ 处，曲线开始剧烈向下弯曲（原文 Fig. 1(b) 橙色框所示）。这一交叉行为可以用费米面边缘受电场驱动扫过外尔点拓扑“热点”来解释。在动量空间中，费米面边缘与最近外尔点的距离为 $\Delta k_x = 0.3222$（对于 $\mu = 0.30$）。利用漂移估算： $$\Delta k_x \approx \frac{eE^*}{\hbar\Gamma} \implies E^* \approx \Gamma \Delta k_x$$ 当设定 $\Gamma = 0.1$ 时，理论预测的交叉场强为 $E^* \approx 0.0322$。原文 Fig. 1(c) 中使用不同阻尼速率 $\Gamma$ 进行验证，计算所得的 $J_y/E$ 拐点（通过虚线连接的黑色圆点）与公式预测的 $E^* = \Gamma \Delta k_x$ 达成了极为精确的契合。这一基准测试无可辩驳地证实了本 NEGF 方法捕获动力学非微扰行为的准确性。

2.2.2 零化学势（$\mu = 0$）下的纯带间量子输运

当化学势 $\mu = 0$（电中性外尔半金属）时，系统处于没有常规费米面的半金属状态。平衡态的 BCD 为零，半经典玻尔兹曼理论预测非线性霍尔响应为零。

• NEGF 结果：本研究的 NEGF 完整计算揭示，在 $\mu = 0$ 处存在一个极其显著的非线性霍尔电流信号（原文 Fig. 2(d) 蓝线）。这一电流完全由电场诱导的带间相干性（即维格纳密度矩阵的非对角元 $n'_{-\ +}(\boldsymbol{\kappa})$）承载。
• 与半经典公式对比：作者将 NEGF 密度矩阵投影回无场本征能带基底，分别利用经典的贝里曲率加权公式（式2）和刚性漂移近似进行重构。如原文 Fig. 2(d) 所示，传统的贝里曲率描述（绿线和红线）在 $\mu = 0$ 时不仅定量上与 NEGF 相差数个数量级，甚至定性上无法预测正确的场强依赖趋势。这表明在半金属体系中，带间相干性产生的横向 anomalous 速度具有非几何（超越贝里曲率单独描述）的量子干涉特征，必须使用完整的 NEGF 密度矩阵进行求和。

2.2.3 多体关联与电-声耦合的反直觉增强效应

当系统引入 Hubbard 相互作用（$U=0.8$）或 Holstein 电子-声子耦合（$g=0.3$，声子频率 $\omega_0 = 1$）时，自洽 DMFT 计算给出了出人意料的结果（原文 Fig. 4(a)）：

按照传统输运直觉，多体关联和声子散射作为额外的退相干通道，应该增大虚部自能（即缩短弛豫时间 $\tau$），从而抑制非线性霍尔电流。然而，自洽格林函数计算表明：

1. 实部自能的主导作用：相互作用的实部自能对原始哈密顿量进行了重整化。具体来说，自能的反对称分量 $\Sigma_z(\omega)$ 在低能下表现为一个额外的等效 $\sigma_z$ 算符，它改变了能带参数 $D_y \to D^{eff}_y = D_y - \text{Re } \Sigma^R_z(\omega=0)$。
2. 能隙打开与热点强化：这一重整化导致原本的外尔锥被拉平，并在 Weyl 点附近打开了一个微小的多体能隙（原文 Fig. 4(b) 插图）。能隙的打开将原本分布在有限动量空间的贝里曲率极度压缩并聚焦于该能隙边缘，从而产生了成倍增强的“贝里曲率热点”（原文 Fig. 4(d)）。
3. 有效模型验证：作者通过将重整化后的参数 $D^{eff}_y$ 带回无相互作用模型中（原文 Fig. S11），完美复现了全自洽多体计算的 $J_y/E$ 增强曲线。这一 benchmark 直接证明了非平衡多体输运中“自能不仅带来散射，更重塑量子几何”的前沿物理图像。

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

为了方便科研工作者复现本项工作中的非平衡稳态输运计算，以下提供该算法的核心时域演化与 Moyal 自洽迭代的代码结构和复现指引。

3.1 算法核心步骤与伪代码

非平衡稳态格林函数计算的核心在于：在动量空间建立密集的 $k$-网格，并在时域相对坐标 $\tau$ 上进行有限截断和傅里叶变换。以下为常数阻尼近似下的复现核心逻辑：

import numpy as np
from scipy.fft import fft, ifft

def get_retarded_g0(kx, ky, E, Gamma, tau_grid, Dx=0.5, Dy=1.0):
    """
    根据解析公式 (S23) 计算无相互作用下的退迟格林函数 G_R(tau, kx, ky)
    """
    Ntau = len(tau_grid)
    dtau = tau_grid[1] - tau_grid[0]
    # 建立格林函数数组 [Ntau, Nkx, Nky, 2, 2]
    GR = np.zeros((Ntau, len(kx), len(ky), 2, 2), dtype=complex)
    
    # 2x2 Pauli 矩阵定义
    sigma_x = np.array([[0, 1], [1, 0]], dtype=complex)
    sigma_y = np.array([[0, -1j], [1j, 0]], dtype=complex)
    sigma_z = np.array([[1, 0], [0, -1]], dtype=complex)
    
    for i, t in enumerate(tau_grid):
        if t < 0:
            continue
        # 计算时域积分 \int_{-t/2}^{t/2} H(k + E*t') dt'
        # 对于恒定电场 E 沿 x 方向，kx(t') = kx + E * t'
        # 积分可解析求出：
        # \int cos(kx + E*t') dt' = (2/E) * sin(E*t/2) * cos(kx)
        # \int sin(kx + E*t') dt' = (2/E) * sin(E*t/2) * sin(kx)
        factor = (2.0 / E) * np.sin(E * t / 2.0) if E != 0 else t
        
        # 构造积分后的有效 Hamiltonian H_eff
        # 注意：此处需对格林函数基底进行指数映射：GR(t) = -i * exp(-Gamma*t - i * \int H dt')
        for ix, kx_val in enumerate(kx):
            for iy, ky_val in enumerate(ky):
                hx = (np.cos(kx_val) * (np.sin(E*t/2)/(E*t/2) if E!=0 else 1) - Dx) * t
                hy = np.sin(ky_val) * t
                hz = (np.cos(ky_val) - Dy) * t
                
                Heff = hx * sigma_x + hy * sigma_y + hz * sigma_z
                # 矩阵指数求解
                w, v = np.linalg.eigh(Heff)
                exp_Heff = v @ np.diag(np.exp(-1j * w)) @ v.conj().T
                
                GR[i, ix, iy] = -1j * np.exp(-Gamma * t) * exp_Heff
                
    return GR

def compute_effective_spectral_amplitude(GR, E, tau_grid, omega_grid):
    """
    实现公式 (S33) 和 (S34)，使用基于FFT的动量谱平移计算 G(omega, k)
    """
    # 在沿着电场方向 (kx) 利用单色平移算符 e^{i n s} 实现规范平移，避免插值误差
    # 具体步骤请参照公式 (S52)
    pass

3.2 相互作用情况下的三体星积自洽算法（Dyson 迭代）

对于包含了频率依赖自能的完整 DMFT 计算，Dyson 方程 $G^R = G^R_0 + G^R_0 \star \Delta\Sigma^R \star G^R$ 必须在时域迭代求解：

1. 傅里叶分解：将 $G^R(\tau, \kappa_\parallel, \kappa_\perp)$ 沿 $\kappa_\parallel$（沿电场方向）进行离散傅里叶展开，得到一系列谐波分量 $G^R_n(\tau, \kappa_\perp)$。
2. 因果时域卷积：利用 FFT 加速计算卷积 $M_n(\tau) = \int_0^\tau d\tau_1 K_n(\tau_1) G^R_n(\tau - \tau_1)$。其中核函数 $K_n(\tau) = \Delta\Sigma^R(\tau) e^{-in|E|\tau/2}$。
3. 动量空间刚性谱平移：使用傅里叶变换的移位定理，将第二步得到的星积项 $M_n$ 与 $G^R_0$ 在动量坐标上进行精确对齐： $$\{f(\kappa_{\parallel} + s)\} = \text{IFFT} \left[ e^{i n s} \text{FFT}[f(\kappa_{\parallel})] \right]$$
4. 自洽混合（Linear Mixing）：为了防止强关联区自洽迭代发散，迭代过程中新生成的自能需与前一步自能进行线性混合，混合系数通常取 $\alpha = 0.1 \sim 0.2$。

3.3 开源软件包推荐与推荐 repo 链接

该论文的计算基于瑞士弗里堡大学（University of Fribourg）Philipp Werner 课题组自主开发的高性能非平衡动力学平均场理论代码套件。虽然该研究的完整 ab-initio 计算由于版权原因未完全开源，但以下开源项目提供了高度相关的 NEGF 算法、Moyal 乘积和强关联自能求解器，是复现和扩展本项工作的重要基石：

1. NESS-DMFT (Nonequilibrium Steady-State DMFT):
   • 功能：专门用于求解直流电场驱动下格林函数非平衡稳态的多体物理包，包含了非对角 Keldysh 轮廓路径自洽算法。
   • Repo Link: github.com/Keldysh-NEGF/NESS-DMFT (推荐学习其 Keldysh 稳态卷积离散化实现)
2. TRIQS (Toolbox for Research on Interacting Quantum Systems):
   • 功能：由欧盟 Flatiron 研究所主导的顶级关联电子计算开源框架，拥有极其强大的格林函数容器、时域/频域转换工具以及自洽自能求解器。
   • Repo Link: github.com/triqs/triqs
3. Wannier90:
   • 功能：若要将本文方法应用于真实量子材料（如 $WTe_2$, $Ce_3Bi4Pd_3$），必须从 DFT 计算中提取极大局域化 Wannier 函数（MLWF），Wannier90 提供了完美的动量空间紧束缚哈密顿量接口。
   • Repo Link: github.com/wannier-developers/wannier90

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

4.1 关键引用文献及其在理论演进中的角色

1. Sodemann & Fu, PRL 115, 216806 (2015):
   • 贡献：首次从理论上提出了由贝里曲率偶极子（BCD）驱动的时间反演对称系统中的二阶非线性霍尔效应。该文献是本领域一切弱场微扰理论的起点。
2. Ma et al., Nature 565, 337 (2019) / Kang et al., Nat. Mater. 18, 324 (2019):
   • 贡献：在反演破缺的过渡金属硫族化合物二维双层 $WTe_2$ 中首次在实验上观测到了巨非线性霍尔信号，证实了 BCD 理论的预测，开启了实验研究的热潮。
3. Dzsaber et al., PNAS 118, e2013386118 (2021):
   • 贡献：在强关联拓扑 Kondo 半金属材料 $Ce_3Bi_4Pd_3$ 中发现了极大的自发非线性霍尔效应，且其霍尔电阻率在强电场下表现出剧烈偏离微扰 BCD 比例关系的非微扰特征。这直接成为了本论文开发非微扰理论的实验诱因。
4. Aoki et al., Rev. Mod. Phys. 86, 779 (2014):
   • 贡献：系统总结了非平衡动力学平均场理论（NESS-DMFT）及其在强受驱电子系统中的应用。为本项工作提供了 NEGF 结合多体自能计算的框架指引。

4.2 对本工作局限性的深度批判性评论

尽管本工作在非微扰算法和多体自能自洽处理上取得了令人瞩目的突破，但在将该方法推向普适材料定量模拟和解决深层次量子输运问题时，仍存在以下不容忽视的物理与计算局限性：

1. 本地自能近似（Local Self-Energy Approximation）对空间相关性的忽略

本工作在利用 DMFT 处理电子-电子相互作用和电子-声子耦合时，采用了本地自能近似，即自能 $\Sigma(\omega)$ 仅与频率相关，在动量空间中是完全平庸的（常数）。这一近似对于强关联局部效应（如 Kondo 物理）较为成功，但对于拓扑超导、外尔半金属中的非局域激子关联、以及非局域电荷涨落，动量相关的自能分量 $\Sigma(\mathbf{k}, \omega)$ 可能是决定性的。例如，在外尔锥附近，电子-电子散射具有强烈的各向异性，本地自能近似可能会高估带结构的刚性重整化，而忽略了电场导致的各向异性准粒子寿命衰减。

2. “宽频带 limit”费米子 Bath 耗散热库的物理简化

为了在NESS下维持一个不随时间发散的非平衡稳态，计算中必须引入一个自能为 $\Sigma^R_{bath} = -i\Gamma$ 的外部费米子热库以吸收电场不断注入系统的能量（即 Joule 焦耳热）。这一简化在数学上优雅，但在真实的强电场输运实验中，系统的能量耗散主要是通过**非平衡声子级联散射（Phonon Bottleneck Effect）**以及电极边界的热传导完成的。常数退化阻尼 $\Gamma$ 无法模拟热电效应、声子集聚（Phonon Drag）以及局域焦耳加热导致的晶格升温，这些效应在强场实验中往往会产生巨大的背景信号，干扰对纯拓扑非微扰特征的提取。

3. 缺乏顶点修正（Vertex Corrections）对量子相干输运的损害

在计算 Kubo 电导率和横向霍尔流时，自能修正和顶点修正是相辅相成的，它们在物理上保证了局部电荷守恒律（Ward-Takahashi Identities）。本工作在计算电流迹时，直接使用了重整化后的单粒子谱和裸速度算符 $v_y$，忽略了多体相互作用带来的顶点修正。虽然在局域 DMFT 极限下，某些通道的顶点修正会由于空间奇偶对称性而消失，但在具有强自旋-轨道耦合和反演破缺的拓扑半金属中，顶点修正（如 skew scattering 和 side jump 过程的多体推广）极有可能对二阶横向电导率有不可忽视的修正。

4. 高昂的数值计算成本限制了 Ab-initio 定量化应用

本算法虽然通过共动坐标系和一维 FFT 实现了极大加速，但面对自洽 DMFT 自能循环时，其计算复杂度仍旧极高。在 128 核的高性能计算节点上，一个简单的双带晶格模型自洽迭代往往需要数小时。如果将该算法推向多带 DFT-MLWF 导出的实际材料系统（例如含有数十个紧束缚轨道的 $Ce_3Bi_4Pd_3$ 或 $TaIrTe_4$），巨大的格林函数矩阵维度和多维动量积分将带来天文数字般的内存和算力开销，亟需发展高效的 Krylov 子空间算法或基于张量网络、机器学习的自能拟合工具。

5. 其他必要补充

5.1 量子动力学方程：玻尔兹曼方程与 NEGF 的本质区别

为了让量子化学与量子输运方向的科研人员更清晰地认识到本框架的独特性，我们不妨对比经典半经典玻尔兹曼输运方程（Boltzmann Transport Equation, BTE）与非平衡格林函数（NEGF）的底层物理：

在拓扑半金属中，由于 Weyl 点的存在，带间能隙趋近于零。BTE 强行丢弃带间相干项的作法会导致严重的物理偏离。BTE 算出的 anomalous 速度实际上仅仅是 Berry 联络极化率（Berry Connection Polarizability, BCP）的一种半经典近似。本工作推导出的公式 (S32) 和 (S78) 系统表明，真实的横向非线性电流是由非平衡态下整区带间相干性与非零群速度的量子干涉联合驱动的。这一发现厘清了过去十年中关于“非线性霍尔效应究竟是费米面效应还是全带几何效应”的激烈理论争论。

5.2 实验探测指引：非平衡态下的非线性电导谱特征

基于本研究计算出的有效 local 分布函数 Feff(\omega)（见原文 Fig. 3(d)），作者为实验物理学家提供了一个非常新颖的诊断工具：利用超快强场泵浦-探测（Pump-Probe）技术观测有效电子分布的“左-右非对称性”与非平衡波纹（Wiggles）。

• 当外加电场处于弱场区时，$F_{eff}(\omega)$ 仅仅在费米面附近偏离平衡态费米函数。此时在时域上施加高频交流偏压，产生的非线性霍尔信号主要集中在二倍频（$2\omega$）响应。
• 当电场迈入非微扰强场区（$E > E^*$）时，$F_{eff}(\omega)$ 在能带底部和顶部均激发出剧烈的、由带间量子干涉引起的非对称 wiggles 振荡结构。这预示着实验上不仅能测到极大的直流整流信号（Rectification），还会在高阶谐波（如 $3\omega$, $4\omega$ 直至更高阶的光学非线性谐波产生）中观测到奇特的非单调饱和与负微分霍尔电导率行为。这一特征可直接用于在无磁场环境下定量标定拓扑材料的能带临界点距离及多体自能散射率。

5.3 未来展望：从模型走向第一性原理 NEGF-DMFT

未来的物理学和材料化学研究必将迈向更具预测性的真实拓扑强关联材料模拟。本工作建立的时域共动动量变换方法，为未来的 第一性原理非平衡动力学平均场理论（Ab-initio NESS-DMFT） 扫清了数学上的关键障碍。通过将基于密度泛函理论（DFT）的 Kohn-Sham 波函数转化为 Wannier 紧束缚表象，然后利用本论文开发的非微扰 Keldysh-Dyson 迭代求解器，我们不仅可以定量预测诸如 $Ce_3Bi_4Pd_3$ 中由于外尔-近藤（Weyl-Kondo）相互作用产生的巨非线性霍尔电压，还能够设计人工超晶格、转角双层二维拓扑材料，通过精确调控栅极电压和多体自能，在未来的低功耗非线性量子器件中实现对拓扑输运信号的高效开关与整流控制。