Problem of nonlinear conductivity within relaxation time approximation in noncentrosymmetric insulators

技术摘要：非中心对称绝缘体中弛豫时间近似下的非线性电导率

问题陈述
本文探讨了在计算非中心对称绝缘体的非线性电学响应时，广泛使用的弛豫时间近似（RTA）在应用于量子动力学方程时存在的一个关键缺陷。虽然 RTA 是描述非平衡稳态的标准唯象工具，但近期的研究表明，在处理线性响应机制时，它会产生不物理的结果，具体表现为在弱电场下预测绝缘体具有有限的线性电导率。本文研究了类似的致命缺陷是否存在于非线性响应机制（特别是二阶响应）中，并提出了一种修正的理论框架，以准确描述这些现象。

研究方法
作者分析了一个在直流电场 $E$ 下具有有限能隙 $\Delta_k$ 的双带绝缘系统，其中电场通过皮尔斯代换（Peierls substitution）引入。研究分为两个阶段进行：

1. RTA 分析： 作者在标准 RTA 框架下，将密度矩阵 $\rho_k(t)$ 展开至电场 $E$ 的二阶。他们推导了电流的带内（intraband）和带间（interband）贡献的表达式，并将电导率分解为线性（$\sigma^{(1)}$）和非线性（$\sigma^{(2)}$）分量。该分析利用绝热波函数（快照基底）和 Redfield 方程框架来定义跃迁偶极矩矩阵元和贝里联络（Berry connections）。
2. 动力学相位近似 (DPA)： 为了修正识别出的缺陷，作者应用了 DPA 方法，该方法已在其先前的工作（Phys. Rev. B, 109, L180302 (2024)）中建立。这种方法通过 Redfield 方程，从微观角度推导了与热库耦合的量子开放系统的耗散项。其中的关键步骤是在 Redfield 方程积分项中对动力学相位因子进行线性化。针对非线性机制，作者通过引入位移矢量 $R_k$，将跃迁偶极矩矩阵元 $[W_k(t-s)]_{+-}$ 的时间演化展开至电场 $E$ 的二阶，从而完善了该近似。

主要结果

• RTA 在非线性机制中的失效： 分析表明，标准的 RTA 会产生不物理的有限非线性电导率项。具体而言，二阶电导率的带间贡献 $\sigma^{(2)}_{inter}$ 包含与阻尼率 $\tau^{-1}$ 和位移矢量 $R_k$ 成正比的项。即使在不存在费米面的情况下，这些项依然存在，导致在弱电场下对非中心对称绝缘体中非互易输运（nonreciprocal transport）做出了违反直觉的预测。
• 抵消机制： 作者证明，RTA 中的不物理项源于在低场极限下对非对角密度矩阵元素的错误近似。在精确计算（以及 DPA 框架）中，项 $[W_k(t)]_{+-}/\Delta_k(t)$ 得到了修正，这导致了与 $\tau^{-1}$ 成正比的不物理项的抵消。
• DPA 修正： 通过应用 DPA，研究表明不物理的带间线性电导率被完全抵消。将此扩展到非线性机制，改进后的 DPA 近似成功消除了依赖于位移矢量和阻尼率的伪二阶项。因此，该理论框架恢复了线性响应机制中预期的绝缘行为（电导率为零），并提供了一个一致的框架来计算非线性响应，避免了 RTA 固有的伪影问题。

意义与主张
本文声称，识别并修复 RTA 中存在的缺陷，对于准确理解量子材料中的非线性及非平衡现象至关重要。所提出的基于 DPA 的方法提供了一个简单且有效的替代方案，避免了 RTA 在预测绝缘体具有有限电导率方面的致命缺陷。作者将该方法定位为研究非线性及非微扰效应（例如非互易输运和位移电流）的有用工具，特别是在位移矢量起关键作用的缺乏反演对称性的系统中。本研究并非提出新的实验装置或特定的材料应用，而是侧重于建立一个用于计算绝缘体系中非线性电导率的理论一致性框架。