中度关联金属 SrVO3 光谱与输运性质的定量 DFT+DMFT 描述：基于实频多轨道迭代微扰理论的统一框架

来源论文: https://arxiv.org/abs/2607.07378v1 生成时间: Jul 09, 2026 06:01

0. 执行摘要

定量、物质特异性（material-specific）地描述强关联材料的光谱和输运性质，一直是凝聚态物理和量子化学计算领域的核心挑战之一。传统的密度泛函理论（DFT）在处理局域电子关联时存在本质不足，而动力学平均场理论（DMFT）通过将晶格问题映射到自洽嵌入的杂质模型，为处理局域关联提供了强大的工具。然而，由于大多数多轨道 DMFT 计算依赖于连续时间量子蒙特卡洛（CT-QMC）等虚频杂质求解器，面临解析延拓（Analytical Continuation）的不确定性以及低温下的符号问题，这极大限制了理论与输运实验的直接定量对比。

本研究采用了一种高效的、基于实频的多轨道迭代微扰理论（Multi-Orbital Iterative Perturbation Theory, MO-IPT）作为 DMFT 杂质求解器，系统研究了典型中度关联金属 $\text{SrVO}_3$ 的单粒子光谱、直流（DC）电阻率以及光学电导率。研究表明，在合理调节相互作用参数 $(U, J)$ 的前提下，只要保持准粒子权重 $Z$ 恒定（对于 $\text{SrVO}_3$，该实验值为 $Z \approx 0.57$），计算所得的低能光谱特征和输运系数均表现出高度的普适标度性（universal scaling）。计算不仅成功再现了费米液体行为（$\text{Im}\Sigma(\omega) \propto \omega^2$），还精确匹配了超净单晶样品在 $0 - 300\text{ K}$ 范围内的直流电阻率温度曲线，并完美解释了实验中观察到的 $\approx 70\text{ meV}$ 低能间带跃迁（interband transition）特征。该工作不仅证实了 $\text{SrVO}_3$ 中自能的 $\mathbf{k}$ 依赖性较弱，还表明顶点校正（vertex corrections）在其中不占主导地位，为中度关联多轨道体系的定量计算开辟了高效、可靠的新途径。

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

1.1 核心科学问题与物理背景

$\text{SrVO}_3$ 是一种典型的具有钙钛矿结构的 $d^1$ 关联金属，其费米能级附近的电子态完全由孤立的 $\text{V}-3d$ 的 $t_{2g}$ 轨道贡献。由于其结构简单、对称性高（空间群为 $Pm\bar{3}m$）且物性极其稳定，它成为了强关联理论方法事实上的“试金石”（benchmark）。

尽管早期的 DFT+DMFT 研究在再现其质量重正化（$m^*/m \approx 2$）、带宽收缩和哈伯德子带（Hubbard bands）形成方面取得了显著成功，但在定量且统一地描述其光谱与输运性质（尤其是光学电导率和直流电阻率的温度依赖性）方面，仍面临以下关键科学问题：

1. 一致的相干尺度：是否存在一个单一的准粒子重正化因子（或相干尺度） $Z$，能够同时、自洽地定量解释单粒子自能、低温 DC 输运以及高频光学响应？
2. 参数依赖性与普适性：DMFT 计算中常用的相互作用参数（如静电库仑作用 $U$ 和洪特耦合 $J$）在文献中存在较大争议，那么最终的输运性质是否强烈依赖于具体的 $(U, J)$ 取值？
3. 低能光学结构的微观起源：实验中在 $\approx 70\text{ meV}$ 处观察到的红外吸收峰，其微观机制究竟是来自关联诱导的哈伯德物理，还是晶格多轨道杂化导致的能带分裂？

1.2 理论基础：多轨道哈密顿量与 DFT+DMFT 框架

为了描述 $\text{SrVO}_3$ 的低能关联物理，首先通过 DFT 结合最大局域化 Wannier 函数（MLWFs）方法，构建 $t_{2g}$ 轨道的紧束缚哈密顿量，作为无相互作用的基准输入。体系的多轨道 Hubbard 哈密顿量在二次量子化形式下写为：

其中，$\hat{c}_{i\alpha\sigma}^\dagger$（$\hat{c}_{i\alpha\sigma}$）表示在格点 $i$、轨道 $\alpha$、自旋 $\sigma$ 上创建（湮灭）电子的算符，$\hat{n}_{i\alpha\sigma} = \hat{c}_{i\alpha\sigma}^\dagger \hat{c}_{i\alpha\sigma}$ 是粒子数算符。前两项构成单粒子紧束缚哈密顿量 $\hat{H}_0^{\text{tb}}(\mathbf{k})$，表征轨道内能与轨道间跃迁积分。第三项代表局域密度-密度相互作用，其矩阵元 $U_{\alpha\beta}^{\sigma\sigma'}$ 采用对称化的标准形式：

其中 $U$ 是同轨道库仑排斥能，$J$ 是洪特耦合常数。通过单格点 DMFT 近似，自能被假设为纯局域的（$\Sigma(\mathbf{k}, \omega) = \hat{\Sigma}_{\text{imp}}(\omega)$）。自洽循环中的局域格林函数表示为：

1.3 技术难点与 CT-QMC 的局限性

在传统 DFT+DMFT 中，最主流的杂质求解器是连续时间量子蒙特卡洛（CT-QMC）。然而，CT-QMC 具有两大难以逾越的瓶颈：

1. 解析延拓的病态性：CT-QMC 计算在虚频轴（松原频率）进行，计算光谱函数或实频输运性质（如光学电导率 $\sigma_1(\Omega)$）时，必须通过最大熵法（MaxEnt）进行解析延拓。这一步骤对高频和低频细节都极度敏感，容易人为抹平细微能带分裂，给定量计算带来极大系统误差。
2. 计算复杂性与低温限制：在计算低温输运（如 $T \lesssim 100\text{ K}$）时，QMC 面临指数级增长的采样难度，多轨道体系下的符号问题（sign problem）更使其在超低温下几乎不可行。

1.4 解决方案：多轨道迭代微扰理论（MO-IPT）

为了克服上述困难，本工作采用了一种确定性的实频求解器——多轨道迭代微扰理论（MO-IPT）。MO-IPT 基于对局域自能进行二阶微扰展开（基于哈特里极限），其核心自能拟合公式（ansatz）为：

其中，$\Sigma_{\alpha\gamma}^{(2)}(\omega)$ 是二阶成对泡泡（pair-bubble）自能，其表达式包含三重频率卷积积分：

这里 $\rho_\alpha(\varepsilon) = -\frac{1}{\pi} \text{Im}\tilde{G}_\alpha(\varepsilon^+)$ 为哈特里校正后的局域介质谱密度，$\tilde{G}_\alpha$ 为对应的杂质格林函数。未知系数 $A_\alpha$ 和 $B_\alpha$ 通过将拟合公式分别在无穷大频率极限（高频极限）和原子极限下与精确自能渐近线对齐来确定。

此外，对于多轨道体系中复杂的两粒子关联函数 $\langle \hat{n}_m \hat{n}_{m'} \rangle$，本工作引入了对称解耦近似（symmetric decoupling approximation）：

该近似将多自由度的关联解耦，实现了自洽闭合，极大提高了计算效率，且计算量随轨道数呈多项式标度，完美避开了虚频到实频的转换步骤。

1.5 关键算法优化：布里渊区自适应涂宽方案

在计算局域格林函数和输运积分时，布里渊区积分的精度至关重要，特别是对范霍夫奇点（van Hove Singularities, vHs）的刻画。传统的四面体法或固定高斯涂宽法易引入人工展宽。本工作引入了动量依赖的自适应涂宽方案（adaptive smearing scheme），其德尔塔函数逼近公式为：

其中 $H_2(x) = 4x^2 - 2$ 为二阶埃尔米特多项式。其核心在于涂宽因子 $w_{\mathbf{k}}$ 会根据能带在 $\mathbf{k}$ 点的群速度自动调节：

群速度梯度 $\nabla_{\mathbf{k}} \varepsilon_{\mathbf{k}}$ 并非通过数值差分获取，而是通过对紧束缚哈密顿量求解析导数并使用伴随矩阵恒等式严密导出：

该解析方法在保证数值稳定性的同时，完美保留了超精细的谱学分裂。

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

2.1 相互作用参数 $(U, J)$ 空间中的 Iso-$Z$ 轮廓图

在强关联物理中，相互作用参数 $U$ 和 $J$ 的选择往往因具体计算协议而异。本研究的核心发现之一在于：只要准粒子重正化因子 $Z$ 被确定，具体的 $(U, J)$ 组合对低能物理的影响极微。如图 1 所示，作者首先在广泛的区间（$U \in [3.0, 6.0]\text{ eV}$, $J \in [0.2, 0.8]\text{ eV}$）内进行了三维自洽循环，绘制出了等准粒子权重轮廓图（Iso-$Z$ Contours）。

• 库仑排斥项 $U$ 的效应：在固定 $J$ 的情况下，增大 $U$ 会显著增强电子关联，单调减小 $Z$。
• 洪特耦合项 $J$ 的效应：对于 $d^1$ 的 $\text{SrVO}_3$ 体系（处于 1/3 填充，偏离半填充），增大 $J$ 反而会减小有效相干相互作用（通过 $U - 3J$ 通道），从而导致关联效应减弱，使得 $Z$ 增大。因此，等 $Z$ 轮廓在 $U-J$ 平面上表现为正斜率。

通过锁定实验确定的 $\text{SrVO}_3$ 物理重正化尺度 $Z \approx 0.57$（对应质量重正化 $m^*/m = Z^{-1} \approx 1.75$），作者锁定了一条特定的等 $Z$ 轮廓线（对应 $U = 4.5\text{ eV}, J = 0.65\text{ eV}$ 为代表）。

2.2 自能 Benchmark：MO-IPT vs. CT-QMC 及 ARPES 实验

为了验证 MO-IPT 求解器的精确度，作者计算了 $T = 116\text{ K}$ 时自能的虚部，并与精密的 CT-QMC 求解器计算结果以及角分辨光电子能谱（ARPES）提取的实验自能进行了定量对比（图 2）：

• 在低能区（$\omega > -0.2\text{ eV}$），MO-IPT 计算的 $Z|\text{Im}\Sigma(\omega) - \text{Im}\Sigma(0)|$ 曲线与 CT-QMC 计算点几乎完美重合，其斜率与曲率高度一致。
• 理论曲线与由 Aizaki 等人和 Kobayashi 等人测得的低温（$< 20\text{ K}$）ARPES 自能实验数据符合极佳。自能完全遵循 $\text{Im}\Sigma(\omega) \propto \omega^2$ 的经典费米液体规律，证明实频微扰论在远离半填充的多轨道非对称体系中具有不亚于 QMC 的极高精度。

2.3 输运 Benchmark I：温度依赖的直流电阻率 $\rho(T)$

由于多轨道哈伯德哈密顿量只显式包含电子-电子（$e-e$）散射，为了完全模拟实验中直至室温的 DC 电阻率，作者通过现象学 Bloch-Grüneisen (BG) 公式融入了电子-声子（$e-ph$）散射贡献：

其中，$\rho_{e-e}(T) = 1/\sigma_{\text{DC}}(T)$ 纯粹来自 DMFT 自洽电导率。$\rho_{e-ph}(T)$ 表达式为：

拟合优化的输运参数如下：

• 德拜/输运温度：$\Theta_R = 800\text{ K}$
• 电子-声子耦合强度前因子：$\rho_R = 2.5 \times 10^{-4} \ \Omega\cdot\text{cm}$
• 唯象尺度常数：$\sigma_0 = 1.21 \ \Omega^{-1}\text{cm}^{-1}$

计算结果与实验对比分析（图 3）：

• 计算得到的总电阻率 $\rho_{\text{Theory}}(T)$ 在 $0 - 300\text{ K}$ 全区间内与高纯度单晶样品（$SVO-195$，残余电阻比 $RRR = 130$）的实验曲线达成定量匹配。
• 低温区（$T < 25\text{ K}$）：电阻率呈现完美的 $T^2$ 行为（$\rho(T) = \rho_0 + A T^2$），由纯电子-电子局域关联主导。
• 交叉温度（$T \approx 80\text{ K}$）：$\rho_{e-e}(T)$ 与 $\rho_{e-ph}(T)$ 贡献相等，在此温度之上，声子散射迅速超越关联散射，并在高温度区（直至 $300\text{ K}$）表现为准线性温度依赖。

2.4 输运 Benchmark II：高分辨红外光学电导率 $\sigma_1(\Omega)$

利用 Kubo 公式，在无顶点校正的泡泡近似下计算实频光学电导率 $\sigma_1(\Omega)$（图 4）：

• 在 $T \approx 6\text{ K}$ 极低温下，理论曲线重现了双峰结构特征：一是由内能带内跃迁贡献的、尖锐相干的低温德鲁德（Drude）峰；二是在 $\Omega \approx 70\text{ meV}$ 处的陡峭红外间带跃迁峰（interband peak）。
• 随着温度从 $0\text{ K}$ 升高到 $300\text{ K}$（图 5），由于电子关联引起的散射率增加，相干德鲁德峰不断展宽， spectral weight 发生重分布，导致 $70\text{ meV}$ 峰逐渐变宽并被背景吞噬，这与实验中观察到的“相干-不相干”金属转变特征完全一致。

3. 代码实现细节、复现指南与开源工具链

为了使科研同行能够快速复现本项工作中的核心计算流程，以下整理了完整的开源工具链、具体的计算输入配置及实现逻辑：

 ┌─────────────────┐
 │ Quantum ESPRESSO│ ──> PBE-GGA 自洽计算 (提取波函数 & 电荷密度)
 └─────────────────┘
          │
          ▼
 ┌─────────────────┐
 │    Wannier90    │ ──> 投影 V-t2g 轨道 (生成无相互作用紧束缚 H0(k))
 └─────────────────┘
          │
          ▼
 ┌─────────────────┐
 │Custom MO-IPT/DMFT│ ──> 实频自洽迭代自能 Σ(ω) 及局域格林函数 G(ω)
 └─────────────────┘
          │
          ▼
 ┌─────────────────┐
 │  Kubo-Transport │ ──> Peierls 近似计算光学电导率与低温 DC 电阻率
 └─────────────────┘

3.1 步骤 1：DFT 自洽与能带计算 (Quantum ESPRESSO)

• 使用软件：Quantum ESPRESSO (开源官网)
• 主要设置：
  • 泛函：广义梯度近似下的 Perdew-Burke-Ernzerhof (PBE-GGA)。
  • 赝势：采用超软赝势（Ultrasoft Pseudopotentials），截断能设为：波函数截断 60 Ry，电荷密度截断 480 Ry。
  • 晶格常数：采用实验测定的立方钙钛矿常数 $a = 3.84\text{ Å}$，空间群为 $Pm\bar{3}m$。
  • 动量网格：$10 \times 10 \times 10$ 莫恩克霍斯特-派克（Monkhorst-Pack）$k$ 网格。

3.2 步骤 2：最大局域化 Wannier 函数构建 (Wannier90)

• 使用软件：Wannier90 (GitHub Repo)
• 主要步骤：
  • 选取费米能级附近的能带，将其投影至以 $V$ 原子为中心的 $3d$ 轨道之 $t_{2g}$ 空间（$d_{xy}, d_{yz}, d_{zx}$）。
  • 通过最小化高斯传播展宽（spread minimization）完成能带解缠（disentanglement）。生成的紧束缚哈密顿量无相互作用带宽 $W \approx 2.5\text{ eV}$，具有极高的重现精度。

3.3 步骤 3：实频多轨道 DMFT 循环 (MO-IPT Solver)

虽然本项工作使用了定制的 MO-IPT 求解器代码，但研究者可以使用主流强关联计算平台（如 TRIQS 框架或 SolidDMFT）进行类似复现，或将二阶微扰理论（IPT）求解器插件嵌入以下框架：

• TRIQS (Toolbox for Research on Interacting Quantum Systems): GitHub Repo
• TRIQS/DFTTools: GitHub Repo
• MO-IPT 核心自洽循环逻辑（Python/Fortran 混合实现）：
  1. 读取 Wannier90 输出的 wannier90_hr.dat 构建动量空间 $\hat{H}_0^{\text{tb}}(\mathbf{k})$。
  2. 初始化试探自能 $\Sigma_\alpha(\omega) = 0$。
  3. 用式 (3) 计算局域格林函数 $\hat{G}_{\text{loc}}(\omega)$，采用动量自适应涂宽法，设置基础参数 $a = 1.0$。
  4. 提取杂质格林函数 $\tilde{G}_\alpha(\omega) = [G_{\text{loc},\alpha}^{-1}(\omega) + \Sigma_\alpha(\omega)]^{-1}$ 及其谱密度 $\rho_\alpha(\omega)$。
  5. 通过式 (14) 的对称解耦算法计算双粒子局域关联能 $\langle \hat{n}_m \hat{n}_{m'} \rangle$，并解出高频极值参数 $A_\alpha, B_\alpha$。
  6. 用式 (5) 快速计算泡泡图二阶自能 $\Sigma^{(2)}(\omega)$。由于在实频轴计算，该步骤涉及的三重积分可使用快速傅里叶变换（FFT）或自适应高斯求积（Adaptive Gaussian Quadrature）进行极速求和。
  7. 更新并平滑自能：$\Sigma^{\text{new}} = (1-d)\Sigma^{\text{old}} + d \Sigma^{\text{calc}}$，收敛容差设为 $10^{-4}$。

3.4 步骤 4：输运性质积分 (Kubo 形式化)

• 使用软件：woptic 工具包（GitHub Repo）或定制的当前-当前关联函数求和脚本。
• 核心物理公式代码化：
  • Peierls 近似下速度算符 $\hat{v}^x(\mathbf{k})$ 的矩阵元生成：通过对 $\hat{H}_0^{\text{tb}}(\mathbf{k})$ 的解析微商得到（参见式 29）。
  • 执行多重循环计算：

# 核心伪代码：计算光学电导率 sigma_1(Omega)
def compute_optical_conductivity(omega_grid, T, Sigma, H_tb, velocity_x):
    # omega_grid: 实频网格, Sigma: 自洽自能, H_tb: 紧束缚哈密顿量, velocity_x: 速度矩阵
    sigma_1 = np.zeros_like(omega_grid)
    for idx, Omega in enumerate(omega_grid):
        integral = 0.0
        for k in k_points:
            # 计算相互作用谱函数 A(k, w) = -1/pi * Im(G(k, w))
            A_k_w = spectral_function(k, omega_grid, Sigma, H_tb)
            A_k_w_plus_Omega = spectral_function(k, omega_grid + Omega, Sigma, H_tb)
            
            # 费米分布函数差分因子
            df = (fermi(omega_grid, T) - fermi(omega_grid + Omega, T)) / Omega
            
            # 计算速度算符迹项
            trace = np.trace(velocity_x[k] @ A_k_w_plus_Omega @ velocity_x[k] @ A_k_w)
            integral += trace * df
        sigma_1[idx] = sigma_0 * integral
    return sigma_1

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

4.1 关键参考文献及其科学贡献

1. Georges et al., Rev. Mod. Phys. 68, 13 (1996)
   • 贡献：系统奠定了动力学平均场理论（DMFT）的理论基石，详细推导了单格点哈密顿量映射以及 Kubo 输运公式的极限形式。
2. Ahn et al., Phys. Rev. B 106, 085133 (2022)
   • 贡献：通过详尽的光谱学实验，首次揭示了 $\text{SrVO}_3$ 超净薄膜在红外低能区（$\approx 70\text{ meV}$）存在本征间带跃迁，为本工作提供了最关键的实验对比基准。
3. Dasari et al., Eur. Phys. J. B 89, 202 (2016)
   • 贡献：发展并基准化了实频多轨道迭代微扰理论（MO-IPT）杂质求解器，验证了其在偏离半填充非对称多轨道体系中的高度可靠性。
4. Kugler et al., Phys. Rev. B 113, L081105 (2026)
   • 贡献：采用最精密的 CT-QMC 作为求解器，计算了 $\text{SrVO}_3$ 的自能与电阻率，是本项工作最重要的数值 Benchmark 参照物。

4.2 本工作局限性与严厉学术批判

尽管本工作在定量匹配和计算效率上取得了令人瞩目的成绩，但作为前沿学术探讨，其计算方案依然存在以下显著的局限性与改进空间：

1. 单格点 DMFT 的局域性近似限制（非局域自能的缺失）：
   • 批判：单格点 DMFT 假设电子自能在动量空间是绝对均匀的（无 $\mathbf{k}$ 依赖性）。然而，在许多中度关联材料中，非局域动力学涨落（如反铁磁涨落、电荷电荷涨落）会使自能表现出明显的 $\mathbf{k}$ 依赖。这导致计算得到的间带跃迁峰相较于实验明显偏窄（如图 4 所示）。未来必须引入团簇 DMFT（如 CDMFT 或 DCA）或双玻色子（dual boson）等图谱修正方法来缓解该问题。
2. Kubo 公式中顶点校正（Vertex Corrections）的忽略：
   • 批判：作者采用无顶点校正的“泡泡近似”（bubble approximation）来计算光学电导率。虽然在无限维度极限下对于单带 Hubbard 模型顶点校正会严格消失，但在真实的三维多轨道体系中，轨道间的带内和带间非局域顶点项已被证实会对绝对光谱权重及输运散射率产生定量重整。完全忽略顶点校正限制了其在极度强关联体系中的推普能力。
3. Peierls 速度矩阵近似的系统性偏差：
   • 批判：式 (29) 采用的 Peierls 近似完全忽略了能带几何学中重要的、源于 Wannier 基底 $\mathbf{k}$ 依赖性的贝里联络（Berry Connection）贡献。在多轨道钙钛矿结构中，当带间跃迁发生于非平庸的能带交叉或近交叉区（avoided crossings）时，非阿贝尔贝里联络项将显著修正带间吸收强度。未来应采用 woptic 工具中实现的完全规范共变（gauge-covariant）偶极矩阵形式进行修正。
4. 现象学 Bloch-Grüneisen 参数拟合的唯象性：
   • 批判：本工作中的声子物理并未实现首原理意义上的自洽耦合，而是通过经验参数 $\Theta_R$、$\rho_R$ 进行了后处理拟合。这对于一个以“首原理（Ab-initio）定量预测”为卖点的框架而言，略显逊色。未来应当将第一性原理电子-声子耦合自能（DFPT 计算所得的 $\Sigma_{\text{ep}}$）直接反馈至 DMFT 自洽回路中，实现真正无拟合参数的统一输运计算。

5. 补充深度解析：间带跃迁的微观起源与 MIR 峰对比

为了给读者提供更透彻的物理图像，本节对论文附录部分关于 $\approx 70\text{ meV}$ 低能红外峰的微观物理机制以及其与传统中红外（MIR）哈伯德峰的区别进行深度补充解析。

5.1 70 meV 峰的消逝：轨道非对角跃迁的决定性证明

为了探究 $\text{SrVO}_3$ 在 $70\text{ meV}$ 处的红外峰究竟是由何种轨道杂化驱动的，作者做了一个精妙的理论对比实验。在构建的 $3 \times 3$ 的 $t_{2g}$ Wannier 哈密顿量中，其矩阵元素不仅包含轨道内跳跃（对角项，如 $t_{ij}^{\alpha\alpha}$），还包含由于晶格结构畸变和间接杂化导致的轨道间跳跃（非对角项，如 $t_{ij}^{\alpha\beta}, \alpha \neq \beta$）。

作者尝试人为将所有轨道非对角跃迁元素强制设为零（即构建 Diagonal-only Hamiltonian），然后重新计算其非相互作用能带结构与光学电导率（图 11）：

• 能带结构对比：在完整的哈密顿量中，由于非对角跳跃的存在，能带在 $\text{R}-\Gamma$ 和 $\Gamma-\text{M}$ 方向上发生了显著的避免交叉（avoided crossing），产生了明显的带分裂。而在去除了非对角项的“纯对角晶格”中，这些带分裂完全消失，能带发生高度简并。
• 光学电导率对比：如图 11(b) 所示，一旦切断了非对角通道，计算的光学电导率曲线中那颗耀眼的 $70\text{ meV}$ 红外峰瞬间彻底消逝，谱权重几乎跌落为零。
• 物理结论：这无可辩驳地证明了，$\text{SrVO}_3$ 的低能红外跃迁绝非源自于常规的电子局域相互作用或莫特物理，其本质是一种晶格轨道间杂化驱动的、多带分裂间的非对角偶极跃迁。这也解释了为什么传统的单带 Hubbard 模型永远无法再现这一极低能区的光学结构特征。

5.2 70 meV 峰与 1.5 eV 中红外（MIR）峰的本质区别

在过渡金属氧化物（如 $\text{SrVO}_3$、$\text{CaVO}_3$）的光谱研究中，文献中经常提及一个位于 $\approx 1.5\text{ eV}$ 处的宽阔“中红外峰”（Mid-Infrared Peak）。读者极易将本工作探讨的 $70\text{ meV}$ 峰与之混淆。以下是两者物理本质的深度剖析：

这一对比再次强调了本工作的核心科学贡献：通过引入基于实频的高效多轨道 DMFT 计算，本研究第一次在全频段范围内（从极低能的间带杂化分裂能标 $\approx 70\text{ meV}$，到高能的莫特-哈伯德关联能标 $\approx 1.5\text{ eV} - 5.0\text{ eV}$）对关联金属进行了无缝、定量且统一的物理刻画，为强关联材料光谱输运的理性设计（如高性能透明导电氧化物 TCO 的开发）奠定了坚实的理论范式。