Connection between the GKSL master equation and the Landauer formula

技术摘要：GKSL 主方程与兰道尔公式之间的联系

问题陈述
介观系统（如分子结和量子点）中的电子输运通常使用两种截然不同的框架进行分析：格林函数方法（适用于相干输运）和主方程方法（常用于涉及耗散和去相干的弱耦合机制）。虽然两者都被广泛应用，但它们之间的理论联系尚未得到充分理解。Jin 等人的前人工作在高温和大化学势极限下，为特定的系统-库（system-bath）模型建立了联系，表明其收敛于 Meir–Wingreen 和兰道尔（Landauer）公式。然而，仍需在广义非相互作用系统的 Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) 主方程形式内，推导出电流公式，并精确识别其还原为兰道尔公式的条件，以阐明这种联系在更广泛条件下的情况。

方法论
作者分析了一个连接到左（L）和右（R）电极的非相互作用电子系统，电极由费米分布函数表征。分析在满足以下假设的 GKSL 主方程框架下进行：

1. 马尔可夫近似与宽带近似 (WBLA)： 每个轨道与电极之间的耦合足够弱，使得费米分布在相关能量窗口内可以被近似为与能量无关。
2. 哈密顿量： 系统由一个通用的电子哈密顿量 $\hat{H} = \sum_{ij} h_{ij}\hat{c}^\dagger_i \hat{c}_j$ 描述。
3. 动力学： 系统密度矩阵 $\rho$ 的时间演化由 GKSL 方程控制，该方程包含了描述从电极注入和提取电子的 Lindblad 超算符。系数 $\alpha_{a,i}$ 和 $\beta_{a,i}$ 是由耦合强度 $\gamma_{a,i}$ 和费米分布 $f_{a,i}$ 推导而来的。
4. 稳态： 作者通过令时间导数为零来求解稳态单电子密度矩阵 $n$。这导致了一个形式为 $AX - XB = Y$ 的 Sylvester 方程。
5. 数学推导： 利用涉及围道积分的 Sylvester 方程的一般积分解，作者利用有理度（retarded, $G^r$）和反理度（advanced, $G^a$）格林函数来表达稳态密度矩阵 $n$。随后将该表达式代入电流定义中，以推导出通过格林函数表示的电流公式。

主要贡献与结果
本文推导了在 GKSL 形式下从左电极流向系统的电流 $I_L$ 的显式表达式。所得电流 $I_L$ 为：
$$I_L = \frac{q}{h} \int_{-\infty}^{\infty} d\epsilon \, \text{Tr} \left[ f_L \gamma_L G^a(\epsilon) \gamma_R G^r(\epsilon) - f_R \gamma_R G^a(\epsilon) \gamma_L G^r(\epsilon) \right]$$
该公式允许与标准的兰道尔公式进行直接比较：
$$I = \frac{q}{h} \int_{-\infty}^{\infty} d\epsilon \, \text{Tr} \left[ \gamma_L G^r(\epsilon) \gamma_R G^a(\epsilon) \right] (f_L(\epsilon) - f_R(\epsilon))$$

作者确定了两个使 GKSL 推导出的公式与兰道尔公式一致的具体条件：

1. 对角哈密顿量（独立模式）： 如果系统哈密顿量是对角的（代表独立的多通道模式），则 GKSL 电流在弱耦合极限下简化为与兰道尔公式相同的形式。
2. 能量无关的费米分布： 如果费米分布函数 $f_L$ 和 $f_R$ 在感兴趣的能量范围内可以被视为能量无关的常数（如模型示意图所示），则兰道尔和 GKSL 表达式都简化为相同的简化形式：
   $$I = \frac{q}{h} (f_L - f_R) \int_{-\infty}^{\infty} d\epsilon \, \text{Tr} \left[ \gamma_L G^r(\epsilon) \gamma_R G^a(\epsilon) \right]$$
   作者指出，当所有单粒子轨道能量严格位于电极的化学势之间（$\mu_R < \epsilon_i < \mu_L$）且温度足够低（$k_B T \ll \mu_L - \mu_R$）时，能量无关分布的条件在物理上是合理的，此时 $f_L \approx 1$ 且 $f_R \approx 0$。

意义
本文为连接非相互作用系统的 GKSL 主方程方法与兰道尔公式提供了严谨的推导。通过利用 Sylvester 方程的一般解，作者证明了这两种框架并非本质上相互矛盾，而是在特定的物理约束下趋于等价。具体而言，当费米分布的能量依赖性在相关能量窗口内可以忽略不计时，这种等价性成立。这项工作阐明了主方程的耗散、顺序隧穿图景与兰道尔形式的相干输运图景在何处发生收敛的理论边界条件。