来源论文: https://arxiv.org/abs/2603.26598v1 生成时间: Apr 01, 2026 15:34

执行摘要

随着纳米技术和微纳加工工艺的飞速发展，热管理已成为限制高性能集成电路、能量转换器件以及量子计算硬件性能的核心瓶颈。传统的辐射传热理论（如普朗克定律）在远场尺度上表现卓越，但在亚波长甚至亚纳米尺度的“近场”区域，由于瞬态场（Evanescent waves）的隧穿效应，热流密度可超越黑体限几个数量级。长期以来，近场辐射传热（NFRHT）主要由基于波动电动力学（FE）的 Rytov 框架描述，但该框架在处理非局域响应、多体耦合以及远离局部热平衡（LTE）的驱动系统时面临严峻挑战。

本综述深入探讨了**非平衡格林函数（NEGF）**作为一种统一量子多体理论框架，如何彻底改变我们对纳米尺度辐射传热的认知。NEGF 不仅能够自然地处理电子、声子与光子的协同输运，还能通过第一性原理计算（DFT+RPA）实现无参数的定量预测。本文将从 NEGF 的哈密顿量构建出发，详细推导 Meir-Wingreen 传热公式，解析其如何解决 FE 框架下的发散难题，并重点阐述其在拓扑热学、扭曲电子学（Twistronics）以及主动冷却调控（如 Floquet 驱动和电流驱动）中的关键作用。

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

1.1 核心科学问题：为何需要 NEGF？

传统的波动电动力学（FE）建立在局部热平衡（LTE）和涨落耗散定理（FDT）之上。其核心假设是物体内部的随机电流波动仅由其局部温度 $T$ 和线性响应性质（介电函数 $\epsilon$）决定。然而，在以下三个场景中，FE 框架失效：

  1. 亚纳米间隙（Sub-nanometer Gaps）： 当间隙接近原子尺度，材料的电子态发生重叠，辐射、电子隧穿和声子导热不再是独立的加和过程。FE 预测的 $1/d^2$ 热流发散在物理上是不合理的。
  2. 强非平衡驱动： 当系统受到外部交变场（Floquet 驱动）或直流电（Current-driven）作用时，系统的准粒子分布不再遵循玻色-爱因斯坦分布，LTE 假设崩溃。
  3. 非局域效应与量子尺寸效应： 在低维材料（如石墨烯、纳米管）中，电荷涨落具有长程相关性，局部介电函数无法描述宏观热流。

1.2 理论基础：哈密顿量与规范选择

NEGF 框架的第一步是构建耦合系统的总哈密顿量：

$$\hat{H}_{tot} = \hat{H}_{mat} + \hat{H}_{em} + \hat{H}_{int}$$

其中 $\hat{H}_{mat}$ 描述电子自由度（通常采用紧束缚模型），$\hat{H}_{em}$ 描述自由电磁场。为了简化计算并确保物理直观，本文采用了时间规范（Temporal Gauge，$\phi=0$）。在这种规范下，标量势为零，矢量势 $\mathbf{A}$ 成为电磁交互的唯一媒介：

$$\mathbf{E} = -\frac{\partial \mathbf{A}}{\partial t}, \quad \mathbf{B} = \nabla \times \mathbf{A}$$

通过 Peierls 替代引入相互作用：

$$H_{mat} + H_{int} = \sum_{\alpha} \sum_{i,j \in \alpha} c_i^\dagger H_{ij} c_j \exp \left( \frac{e_0}{i\hbar} \int_{r_j}^{r_i} \mathbf{A} \cdot d\mathbf{r} \right)$$

这种处理方式完美融合了近场（静电库仑耦合）与远场（横波辐射）交互，将两者置于对等的理论地位。

1.3 技术难点：光子自能与戴森方程

核心挑战在于计算光子自能（Photon Self-energy, $\Pi$），它代表了材料对电磁场的响应。在随机相位近似（RPA）下，$\Pi$ 由电子格林函数的乘积（极化泡图）给出：

$$\Pi_{\mu\nu}(\mathbf{r}\tau; \mathbf{r}'\tau') = -\frac{i}{\hbar} \langle T_c j_\mu(\mathbf{r},\tau) j_\nu(\mathbf{r}',\tau') \rangle_{ir}$$

一旦获得自能，通过求解**戴森方程（Dyson Equation）**即可获得“穿衣”的光子格林函数 $\mathbf{D}$：

$$\mathbf{D} = \mathbf{v} + \mathbf{v} \Pi \mathbf{D}$$

其中 $\mathbf{v}$ 是真空光子格林函数。$\mathbf{D}$ 包含了系统中所有的多次反射、共振模式以及隧穿效应。

1.4 方法细节：Meir-Wingreen 与 Caroli 公式

NEGF 传热理论的最伟大成就之一是将宏观的坡印廷矢量（Poynting vector）与微观格林函数联系起来。物体 $\alpha$ 发射的净热电流 $I_\alpha$ 可表达为：

$$I_\alpha = \int_0^\infty \frac{d\omega}{2\pi} (\hbar\omega) \text{Tr} \left[ \mathbf{D}^< \Pi_\alpha^> - \mathbf{D}^> \Pi_\alpha^< \right]$$

该公式（Meir-Wingreen 公式）具有清晰的物理图像：热流取决于物体的发射（与 $\Pi^>$ 相关）与系统的吸收（与 $\mathbf{D}^<$ 相关）之间的平衡。

在局部热平衡（LTE）限制下，利用涨落耗散定理（FDT），该公式可退化为类似 Landauer 的形式：

$$I_\alpha = \sum_{\beta} \int_0^\infty \frac{d\omega}{2\pi} \hbar\omega \mathcal{T}_{\beta\alpha}(\omega) [N_\beta(\omega) - N_\alpha(\omega)]$$

其中传输系数 $\mathcal{T}$ 由 Caroli 公式给出：

$$\mathcal{T}_{\beta\alpha}(\omega) = \text{Tr} [\mathbf{D}^r \Gamma_\beta \mathbf{D}^a \Gamma_\alpha]$$

这证明了 NEGF 框架不仅是 FE 的量子化推广，更在逻辑上包含了 FE。

2. 关键 Benchmark 体系与计算性能数据

2.1 石墨烯：解决 $1/d$ 发散难题

石墨烯作为典型的 2D 材料，是检验 NEGF 效力的首选体系。早期 FE 计算预测当间隙 $d \to 0$ 时，由于瞬态模式贡献无穷大，热流会趋于发散。然而，NEGF 计算（如图 1a 所示）显示：

  • 饱和效应： 在 $d < 1 \text{ nm}$ 时，热流趋于平缓，最终在接触极限下达到常数。这是由于电子波矢 $k$ 的有限性（量子截断）导致的。
  • 数据对比： 在 $d = 10 \text{ nm}$ 处，NEGF 预测的热流比黑体限高出约 $10^4$ 倍。对于掺杂浓度 $\mu = 0.1 \text{ eV}$ 的石墨烯，热流在 $d^{-2}$ 幂律区间表现稳定，与实验观察高度吻合。

2.2 金属表面：非局域电磁响应的可视化

对于体相金属（如金或铝），FE 依赖的 Drude 模型忽略了表面电荷涨落的长程相关性。NEGF 研究（图 1c）通过计算极化率自能 $\Pi_{zz}(j, j')$ 发现：

  • 相关长度： 表面格点 $j=1$ 与内部格点 $j'$ 的相关性可延伸至约 $100$ 个原子层（约 $20 \text{ nm}$）。
  • 性能影响： 这种非局域性使得极短距离下的热传导系数比局部模型预测的要低，纠正了传统理论对近场热交换的过高估算。

2.3 维度的魔力：1D 与 2D 的对比

NEGF 能够清晰捕捉几何维度对热输运的影响。研究对比了“线-面”（1D-2D，如纳米管与平面）和“面-面”（2D-2D）构型：

  • 缩放定律： 2D-2D 构型热流随距离 $d$ 的衰减遵循 $d^{-2}$；而 1D-2D 构型（图 1d）衰减则遵循极陡的 $d^{-5}$ 规律。这种差异源于 1D 系统中受限电子态产生的独特电荷涨落空间相关性。

2.4 多通道协同效应：亚纳米间隙下的竞争

在 $d < 1 \text{ nm}$ 的金属接头中（图 2a）：

  • 交叉点（Crossover）： 在 $d \approx 0.92 \text{ nm}$ 处，电子隧穿贡献开始超过光子辐射。
  • 非线性增强： 在 $d < 0.6 \text{ nm}$ 时，总热流大于辐射与隧穿贡献的简单算术和。NEGF 揭示了电子隧穿通过改变界面势垒，实际上“增强”了库仑耦合介导的辐射热流，这种非加和性是理解原子级热管理的关键。

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

3.1 总体工作流 (Workflow)

实现无参数的 NEGF 辐射传热模拟通常遵循以下路径：

  1. 第一性原理计算 (DFT)： 使用 VASP 或 Quantum ESPRESSO 计算材料的基础电子能带结构和波函数。
  2. 万尼尔化 (Wannierization)： 利用 Wannier90 将布洛赫波函数转化为局域化的万尼尔基底，提取紧束缚哈密顿量矩阵元 $H_{ij}$。
  3. 自能计算 (Self-energy)： 根据 Adler-Wiser 公式（式 33）计算不可约极化率 $\Pi^0$。对于 2D 体系，这需要对第一布里渊区进行极细密的 $k$ 点采样。
  4. 戴森方程求解： 在频率域内迭代求解光子格林函数。在大体系下，涉及大型稠密矩阵的求逆，通常需要调用 Intel MKL 或 OpenBLAS 库。
  5. 热流积分： 对频率和动量（$q$ 空间）进行积分得到总热流。

3.2 软件包推荐

  • GPAW / BerkeleyGW： 适合计算准粒子自能和 RPA 极化率，是获取 $\Pi$ 的高精度工具。
  • Kwant (定制版)： 虽然 Kwant 主要用于电子输运，但其底层图形算法可借鉴用于构建复杂几何形状的光子 NEGF 框架。
  • JSW-NEGF Code： 由新加坡国立大学 J.-S. Wang 课题组开发的系列代码（常基于 Fortran/Python），是该领域的学术标杆。部分简化版逻辑已在 GitHub 的相关开源学术库中可见。

3.3 复现核心逻辑 (Python 伪代码示例)

import numpy as np

def solve_negf_rht(omega, q, Pi_1, Pi_2, v_bare):
    """
    求解特定频率和动量下的传输系数
    """
    I = np.eye(3)
    # 求解 retarded Green function Dr = (v^-1 - Pi)^-1
    Pi_tot = Pi_1 + Pi_2
    Dr = np.linalg.inv(np.linalg.inv(v_bare) - Pi_tot)
    
    # 计算谱函数 Gamma = i(Pi_r - Pi_a)
    Gamma_1 = 1j * (Pi_1 - Pi_1.conj().T)
    Gamma_2 = 1j * (Pi_2 - Pi_2.conj().T)
    
    # Caroli 公式计算传输系数
    Transmission = np.trace(Dr @ Gamma_1 @ Dr.conj().T @ Gamma_2)
    return Transmission.real

3.4 关键技术陷阱

  • 规范不变性检查： 必须确保数值结果在不同规范（如时间规范与库仑规范）下的等价性，这通常用于验证自能计算的完备性。
  • $k$ 点收敛性： 对于金属体系，极化泡积分对 $k$ 点非常敏感，通常需要 $1000 \times 1000$ 以上的网格才能获得收敛的远场极限。

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

4.1 关键参考文献

  1. Polder & Van Hove (1971): 奠定了 FE 理论基础，近场辐射研究的起点。
  2. Meir & Wingreen (1992): NEGF 输运理论的通用公式，最初用于电子，后被引入光子领域。
  3. J.-S. Wang et al. (2014, 2017): 关键综述与发展，系统地建立了基于时间规范的光子 NEGF 框架。
  4. Zhu et al. (2020): 首次将 DFT+RPA 与 NEGF 结合，实现了第一性原理预测纳米尺度辐射传热。
  5. Venkataram et al. (2020): 揭示了光子与声子通道在亚纳米间隙下的干涉抑制效应。

4.2 局限性评论

尽管 NEGF 展现了强大的普适性，但仍存在以下局限：

  • 计算成本： 格林函数矩阵求逆的复杂度为 $O(N^3)$。对于包含成千上万个原子的复杂纳米结构（如复杂的超材料阵列），计算量巨大，往往需要借助机器学习代理模型（Surrogate models）。
  • 强相关效应缺失： 目前主流研究基于 RPA 近似。在涉及激子、极化激元强耦合或近藤效应（Kondo effect）的系统中，简单的泡图近似可能失效，需要更高级的顶点校正（Vertex corrections）。
  • 多体散射： 虽然 NEGF 框架理论上支持多体交互，但实际计算大多局限于两体交换。如何高效处理多物体间复杂的相互辐射（如粒子阵列中的多级散射）仍是一个开放课题。
  • 声子-光子耦合的深度局域性： 在极性材料中，由于长程偶极交互，格林函数矩阵变得非常稠密，传统的稀疏矩阵加速方法失效。

5. 其他补充：未来愿景与工程应用

5.1 主动驱动：从“被动传热”到“热机控制”

NEGF 框架最激动人心的应用在于主动非平衡调控。图 4 展示了两种范式：

  • Floquet 相位控制（图 4a）： 通过周期性调制两个物体的介电参数，并引入相位差 $\theta$，可以在无温度梯度的情况下产生净热流。这种“等温热传导”为设计全固态热整流器和开关提供了基础。
  • 电流驱动冷却（图 4c）： 利用石墨烯中的“负兰道阻尼（Negative Landau Damping）”效应。当电子漂移速度超过电磁波相速度时，电子动能可转化为光子发射。足够强的电流甚至能驱动热量从冷端流向热端，实现固态辐射制冷。

5.2 拓扑热学与扭曲电子学

  • SSH 模型与边缘态： 在具有拓扑保护边缘态的纳米管系统中，热流展现出非单调的距离依赖性（图 3a）。通过切换系统的拓扑相，可以实现热交换的“开启”与“关闭”，这为热逻辑门（Thermal Logic Gates）的实现奠定了物理逻辑。
  • Twistronics 热信号： 扭曲双层石墨烯（TBG）在“魔角”附近产生的范霍夫奇异点（Van Hove Singularities）会显著增强特定频率的热辐射（图 3d）。这种可调谐的热特征可用于红外隐身或高精度热探测。

5.3 结语：纳米热学的新纪元

NEGF 框架不仅是一个计算工具，更是一种物理思维的转变。它将辐射传热从宏观电动力学的范畴，拉回到了量子多体相互作用的本质。随着量子信息科学与纳米能量转换技术的交汇，NEGF 必将在开发量子热机、纳米级热信息处理器以及高效热电转换设备中发挥不可替代的作用。未来的研究重点将聚焦于如何将 NEGF 与人工智能结合，以探索更广阔的热学材料设计空间。