来源论文: https://arxiv.org/abs/2602.22789v1 生成时间: Feb 27, 2026 14:06

0. 执行摘要

热电材料能够实现热能与电能的直接转换，在废热回收和固态制冷领域具有巨大的应用潜力。热电转换效率由无量纲优值 $ZT = (GS^2T)/(\kappa_{el} + \kappa_{ph})$ 决定。然而，提升 $ZT$ 面临着相互制约的输运参数之间的微妙平衡难题。近年来，以过渡金属硫族化合物（TMDs）为代表的二维材料因其独特的量子限制效应和尖锐态密度特征，成为热电领域的研究热点。

本研究针对单层二硒化钨（WSe2）纳米带，提出了一种结合单色光辐照（Floquet 工程）与自旋轨道耦合（SOC）效应的协同调控方案。通过构建六轨道紧束缚模型并结合第一性原理计算（DFT/DFPT），研究揭示了辐照引起的跳跃参数重整化如何打破透射谱的对称性，从而将塞贝克系数 $S$ 提高了一个数量级。同时，研究深入探讨了 SOC 对声子动力学的影响，发现其通过增强非谐性散射显著降低了晶格热导率 $\kappa_{ph}$。最终结果显示，在辐照调控下，$WSe_2$ 纳米带的 $ZT$ 值在较宽的温度范围内可远超单位 1，甚至在特定条件下达到 28。这一发现为开发高性能、可逆调谐的二维热电设备提供了重要的理论指导。

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

核心科学问题

在传统的半导体热电体系中，由于电子透射谱在费米能级附近通常表现出近似的对称性，导致电子和空穴对热功率的贡献相互抵消，使得塞贝克系数 $S$ 难以大幅提升。此外，如何在不引入化学掺杂（可能降低迁移率）的前提下，实现输运性质的可逆调谐，是当前凝聚态物理与材料科学的难点。$WSe_2$ 作为 TMDs 的代表，具有强 SOC 和显著的能带分裂，这为调控电子和声子输运提供了基础，但其内在机制仍需系统阐明。

理论基础：Floquet 工程与紧束缚模型

Floquet 工程的核心思想是利用周期性驱动（如光辐照）来调制量子系统的哈密顿量。当系统受到频率为 $\Omega$ 的单色光辐照时，时间周期性哈密顿量 $H(\tau)$ 可以通过 Floquet 理论映射到一个有效静态哈密顿量 $H_{eff}$。在本文研究的高频极限（$\hbar\Omega \gg BW$，其中 $BW$ 为带宽）下，系统主要由零阶 Floquet 扇区主导。光场的影响通过 Peierls 替代引入：

$$\mathbf{k} \to \mathbf{k} + \frac{q}{\hbar}\mathbf{A}(\tau)$$

其中 $\mathbf{A}(\tau)$ 是矢量势。这导致实空间中的跳跃振幅 $t_{nm}$ 被重整化为：

$$\tilde{t}_{nm} = t_{nm} e^{i(p-q)\Theta} J_{p-q}(\Gamma)$$

其中 $J$ 是贝塞尔函数，$\Gamma$ 与光强和跳跃矢量有关。这种重整化是各向异性的，能够改变能带的曲率和对称性。

紧束缚（TB）模型：研究采用了六轨道模型，包括钨（W）的三个 d 轨道（$d_{z^2}, d_{x^2-y^2}, d_{xy}$）和硒（Se）的三个 p 轨道组合。考虑了原子级的自旋轨道耦合项 $\lambda_W$ 和 $\lambda_{Se}$，这对于准确描述 $WSe_2$ 的价带分裂至关重要。

技术难点：多尺度输运模拟的耦合

  1. 电子输运：需要在紧束缚框架下处理大尺寸（30 nm）纳米带的相干输运，计算透射概率 $\mathcal{T}(E)$。利用 Landauer-Büttiker 公式计算电导 $G$、塞贝克系数 $S$ 和电子热导 $\kappa_{el}$。
  2. 声子输运：热导率 $\kappa_{ph}$ 必须通过第一性原理获得。难点在于准确捕捉 SOC 对三阶力常数（IFCs）的影响，这需要昂贵的非谐声子计算。
  3. 协同调控：如何将光致电子能带重塑与 SOC 致声子散射增强效应结合，在热电优值 $ZT$ 公式中寻找最优平衡点。

方法细节

  • DFT 计算：使用 Quantum Espresso 软件包。采用 PBE 泛函和全相对论赝势（显式包含 SOC）。平面波截止动能设为 60 Ry，电荷密度截止为 240 Ry。真空层设为 15 Å 以消除层间相互作用。
  • 声子动力学：通过密度泛函微扰理论（DFPT）获得二阶力常数，并使用有限位移法获得三阶 IFCs。利用 ShengBTE 求解声子玻尔兹曼输运方程（BTE），考虑了 Normal 和 Umklapp 散射过程。
  • 输运模拟：基于 Python 的 Kwant 包实现波函数匹配法。构建了 zigzag 边缘的 $WSe_2$ 纳米带，并将重整化后的 Floquet 紧束缚参数植入其中。

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

Benchmark 体系：单层 Zigzag WSe2 纳米带

研究选择了宽度 $W=15$ nm，长度 $L=10-30$ nm 的纳米带作为基准模型。该体系具有 $D_{3h}$ 点群对称性，在无 SOC 时带隙为 1.63 eV。加入 SOC 后，带隙收缩至 1.33 eV，且在 K 点出现明显的自旋分裂（0.44 eV）。

关键计算数据

  1. 电子态密度（DOS）：

    • SOC 的引入在导带底部诱导出明显的态密度峰。这是由于能带简并度解除以及在窄能级范围内可用电子态密度的增加。
    • 在辐照下，透射概率谱 $\mathcal{T}(E)$ 从连续的平台转变为一系列尖锐的共振峰，且在费米能级附近产生了显著的能隙。这种“能量过滤”效应是提升热电性能的关键。

  2. 塞贝克系数 $S$：

    • 无光照时：$S \approx 60 \, \mu V/K$（在费米能级附近）。
    • 光辐照下（$A_x=1, A_y=2, \theta=\pi/4$）：$S$ 峰值激增至约 $600 \, \mu V/K$。提升了 10 倍。

  3. 电导 $G$ 与电子热导 $\kappa_{el}$：

    • 辐照导致 $G$ 显著下降（约下降 4 倍），因为许多透射通道被关闭。
    • $\kappa_{el}$ 同样被大幅抑制，从 25 nW/K 降至 5 nW/K 以下，这有助于降低总热导 $\kappa$。

  4. 晶格热导率 $\kappa_{ph}$：

    • 在 300 K 时，无 SOC 计算值为 35 W/mK。
    • 含 SOC 计算值为 15 W/mK。SOC 诱导的声子软化（ZA 支在 $\Gamma$ 和 K 点）以及增强的非谐性散射导致热导率下降了约 57%。

最终性能数据：ZT 值

  • 无辐照：$ZT < 0.2$，性能平庸。
  • 辐照调控下：在 $E_F \approx 0.6-0.8$ eV 范围内，$ZT$ 表现出巨大的增强。对于 $L=10$ nm 的纳米带，$ZT$ 最大值可达 28。即使对于较长的纳米带，$ZT$ 依然保持在 10 以上。在空穴输运区（$E_F \approx -0.6$ eV），$ZT$ 也达到了约 8。
  • 温度特性：$ZT_{max}$ 随温度从 100 K 升高至 600 K 呈现单调增加趋势，表明该体系在高温热电应用中具有潜力。

软件栈

  1. Quantum Espresso (QE): 用于结构优化、电荷密度计算和 DFPT 声子计算。 Official Site
  2. ShengBTE: 用于从三阶力常数求解声子 BTE 得到 $\kappa_{ph}$。 GitHub/ShengBTE
  3. Kwant: 用于基于紧束缚模型的量子输运模拟。 Official Site
  4. Wannier90 (可选): 虽然本文使用了文献中的紧束缚参数，但复现时可用此工具从 QE 结果构建自定义 TB 模型。

复现指南

第一阶段：DFT 与声子热导率计算

  1. 执行 vc-relax：获取 $WSe_2$ 的最优晶格常数（本文为 $a=3.32$ Å）。
  2. 执行 scf 与 ph.x：生成二阶力常数。建议使用 $4\times4\times1$ 的 q 点网格。
  3. 执行 third_order.py (ShengBTE 工具)：生成位移超胞。建议使用 $4\times4\times1$ 超胞。计算原子受力获取三阶 IFCs。
  4. 运行 ShengBTE：在 input 文件中开启 SOC 相关修正（如果赝势支持），计算不同温度下的 $\kappa_{ph}$。

第二阶段：Floquet 输运模拟

  1. 在 Python 环境中定义 Kwant 系统：定义 $W$ 和 $Se$ 的位置。根据六轨道模型定义 onsite 和 hopping 矩阵。
  2. 实现 Floquet 算子：
    import numpy as np
def get_floquet_hopping(t_static, d_vec, Ax, Ay, theta):
    # d_vec is the vector between sites m and n
    gamma = np.sqrt((Ax*d_vec[0])**2 + (Ay*d_vec[1])**2 + 2*Ax*Ay*d_vec[0]*d_vec[1]*np.cos(theta))
    return t_static * jv(0, gamma) # jv(0, x) is the 0-th order Bessel function
  3. 计算透射：使用 kwant.smatrix(sys, energy).transmission(1, 0) 获取 $\mathcal{T}(E)$。
  4. 积分求解输运系数：编写小脚本计算 Landauer 积分 $L_n$（公式 12），进而得到 $G, S, \kappa_{el}$。

关键参数复现表（表 I）：

  • $V_{pd\sigma} = 5.803$ eV, $V_{pd\pi} = -1.081$ eV.
  • $\lambda_W = 0.251$ eV, $\lambda_{Se} = 0.439$ eV.

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

关键引用文献

  1. Silva-Guillén et al., Appl. Sci. 6, 284 (2016): 提供了 TMDs 的标准六轨道紧束缚模型框架，是本文电子结构计算的基础。
  2. Groth et al., New J. Phys. 16, 063065 (2014): Kwant 软件包的基石论文，定义了波函数匹配法的标准。
  3. Li et al., Comput. Phys. Commun. 185, 1747 (2014): ShengBTE 算法原文，解决了非谐声子输运的计算问题。
  4. Oka & Aoki, Phys. Rev. B 79, 081406 (2009): Floquet 工程在石墨烯等二维系统中的早期应用，奠定了光致能隙调控的理论。
  5. Tian et al., Phys. Rev. B 85, 184303 (2012): 探讨了 SOC 对重元素体系热导率的影响，是本文声子部分的重要参考。

局限性评论

虽然该研究在理论上展示了惊人的 $ZT$ 值，但在实际量子化学与物理实验中仍存在以下局限性：

  1. 高频极限的物理实现：研究假设 $\hbar\Omega = 20$ eV，这属于极紫外（EUV）波段。在如此高的频率和强场下，材料的损伤阈值、非线性吸收以及多光子激发效应可能非常显著，简单的 Floquet 映射可能失效。
  2. 独立输运假设：公式 13 将电子和声子热导率分开处理。然而在强驱动下，电子-声子耦合（Electron-Phonon Coupling）可能会发生显著变化，导致 Wiedemann-Franz 定律的更复杂违背，独立累加可能导致对 $ZT$ 的过高估计。
  3. 电极接触问题：模型假设了完美的、无反射的半无限大电极。在实际器件中，接触电阻和界面处的声子散射往往是热电效率的主要损失项。
  4. 边缘重构：本研究基于理想的 zigzag 边缘。实际的 $WSe_2$ 纳米带在辐照环境下可能会发生边缘原子重构或氧化，从而改变边缘态并影响能量过滤效应。

5. 其他补充内容

SOC 对声子特性的微观影响机制

本工作的一个亮点是详尽解释了 SOC 如何通过改变电子云分布影响原子间力常数。具体而言：

  • ZA 支软化：$WSe_2$ 的面外弯曲（ZA）模式对长程力常数极为敏感。SOC 引起的电荷重新分布弱化了面外恢复力，导致频率下降。这一现象在含重原子的体系（如 PbTe, SnSe）中已被证实是降低热导率的有效机制。
  • Grüneisen 参数（$\gamma$）的演化：计算显示，含 SOC 系统在低频区域的 $\gamma$ 值显著高于无 SOC 系统（峰值约为 6 vs 4.5）。这意味着格点振动的非简谐程度增加，声子间的碰撞频率提高，直接缩短了声子平均自由程。

热电性能对几何尺寸的依赖性

研究在第 10 页讨论了纳米带宽度 $W$ 对 $ZT_{max}$ 的影响。发现 $ZT_{max}$ 随宽度呈现非单调变化。这是由于几何局限效应（Geometric Confinement）与光致能带重整化之间的竞争。这种非单调性提示我们在设计实际器件时，存在一个最优的“纳米带宽高比”，而非单纯追求更细的纳米线。

未来展望：全光调控热电逻辑电路

基于 Floquet 工程的可逆性，我们可以设想一种“全光开关”热电器件。通过超快脉冲激光调制光强或极化方向，可以在皮秒量级切换器件的热电转换效率，甚至改变热流方向。这为微纳尺度上的热管理和热计算提供了全新的想象空间。

总结：本研究通过严密的跨尺度模拟，论证了外部辐照场作为一种“虚原子”或“虚掺杂”手段，在不改变材料化学组分的情况下，通过重塑量子态分布实现超高性能热电转换的可能性。对于量子化学研究者而言，这不仅是一个材料性能预测工作，更是一次关于非平衡态量子调控的深入探索。