Connection between the GKSL master equation and the Landauer formula

技術要約：GKSLマスター方程式とランダウアー公式の接続

問題提起
分子接合や量子ドットなどのメゾスコピック系における電子輸送は、通常、2つの異なる枠組みを用いて解析される。すなわち、コヒーレントな輸送に適したグリーン関数法と、散逸やデフェージングを伴う弱結合領域で一般的に用いられるマスター方程式法である。これらは共に広く採用されているが、両者の理論的な接続は完全には理解されていない。Jinらによる先行研究では、高温および大きな化学ポテンシャルの極限における特定のシステム・バスモデルにおいて、これらがメイア・ウィングリーン（Meir–Wingreen）およびランダウアー（Landauer）公式へと収束することを示し、両者の関連性を確立した。しかし、一般的な非相互作用系において、ゴリニ・コサック・スダルシャン・リンドブラッド（GKSL）マスター方程式の形式内で電流公式を導出し、どのような条件下でそれがランダウアー公式へと帰着するかを正確に特定することが、より広い条件下でのこの接続を明確にするために必要とされている。

手法
著者らは、フェルミ分布関数によって特徴付けられる左（L）および右（R）電極に接続された、非相互作用電子系を解析している。解析は、以下の仮定の下でGKSLマスター方程式の枠組み内で行われる。

1. マルコフ近似および広帯域近似（WBLA）： 各軌道と電極との結合は十分に弱く、フェルミ分布は関連するエネルギー窓内においてエネルギーに依存しないものと近似できる。
2. ハミルトニアン： システムは、一般的な電子ハミルトニアン $\hat{H} = \sum_{ij} h_{ij}\hat{c}^\dagger_i \hat{c}_j$ によって記述される。
3. ダイナミクス： システムの密度行列 $\rho$ の時間発展は、電極からの電子の注入および抽出を考慮したリンドブラッド・スーパーオペレータを含むGKSL方程式に従う。係数 $\alpha_{a,i}$ および $\beta_{a,i}$ は、結合強度 $\gamma_{a,i}$ とフェルミ分布 $f_{a,i}$ から導出される。
4. 定常状態： 著者らは、時間微分をゼロと置くことで、定常状態の単一電子密度行列 $n$ を解く。これは、$AX - XB = Y$ という形式のシルベスター方程式へと導かれる。
5. 数学的導出： 複素平面上の積分を用いたシルベスター方程式の一般積分解を用いることで、著者らは定常状態の密度行列 $n$ を、遅延（retarded, $G^r$）および先進（advanced, $G^a$）グリーン関数の形式で表現する。この式を電流の定義に代入することで、グリーン関数を用いて表現された電流公式を導出する。

主要な貢献および結果
本論文は、GKSL形式における左電極からシステムへ流れる電流の明示的な表現を導出している。得られる電流 $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から導かれた公式がランダウアー公式と一致する2つの特定の条件を特定した：

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マスター方程式とランダウアー公式を結びつける厳密な導出を提供している。シルベスター方程式の一般解を利用することで、著者らは、これら2つの枠組みが本質的に矛盾するものではなく、特定の物理的制約の下で等価になることを実証している。具体的には、関連するエネルギー窓内においてフェルミ分布のエネルギー依存性を無視できる場合に、この等価性が成立する。本研究は、マスター方程式の散逸的かつ逐次的トンネル現象の描像が、ランダウアー形式のコヒーレント輸送の描像へと収束する理論的な境界条件を明らかにしている。