Connection between the GKSL master equation and the Landauer formula

Technische Samenvatting: De verbinding tussen de GKSL-meestervergelijking en de Landauer-formule

Probleemstelling
Elektronentransport in mesoscopische systemen, zoals moleculaire juncties en kwantumstippen, wordt doorgaans geanalyseerd met behulp van twee verschillende kaders: de Green's functie-benadering (geschikt voor coherente transportprocessen) en de meestervergelijking-benadering (veelal gebruikt in regimes met zwakke koppeling waarbij dissipatie en dephasering een rol spelen). Hoewel beide breed worden toegepast, blijft de theoretische verbinding tussen hen niet volledig begrepen. Eerder werk door Jin et al. heeft een link vastgesteld voor specifieke systeem-bad-modellen onder limieten van hoge temperatuur en grote chemische potentiaal, waarbij werd aangetoond dat zij convergeren naar de Meir–Wingreen- en Landauer-formules. Echter, een afleiding van de stroomformule binnen het Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) meestervergelijkingskader voor een algemeen niet-interagerend systeem, en een precieze identificatie van de condities waaronder het reduceert tot de Landauer-formule, is vereist om deze verbinding onder bredere condities te verduidelijken.

Methodologie
De auteurs analyseren een niet-interagerend elektronensysteem dat verbonden is met een linker (L) en rechter (R) elektrode, gekenmerkt door Fermi-distributiefuncties. De analyse verloopt binnen het GKSL-meestervergelijkingskader onder de volgende aannames:

1. Markov-benadering & Wide-Band Limit Approximation (WBLA): De koppeling tussen elk orbitaal en de elektroden is voldoende zwak, waardoor de Fermi-distributie binnen het relevante energievenster als energie-onafhankelijk kan worden benaderd.
2. Hamiltoniaan: Het systeem wordt beschreven door een algemene elektronische Hamiltoniaan $\hat{H} = \sum_{ij} h_{ij}\hat{c}^\dagger_i \hat{c}_j$.
3. Dynamica: De tijdsevolutie van de dichtheidsmatrix $\rho$ wordt beheerst door de GKSL-vergelijking, die Lindblad-superoperatoren bevat die rekening houden met elektroneninjectie en -extractie uit de elektroden. De coëfficiënten $\alpha_{a,i}$ en $\beta_{a,i}$ zijn afgeleid van de koppelingssterkten $\gamma_{a,i}$ en de Fermi-distributies $f_{a,i}$.
4. Stationaire toestand: De auteurs lossen de stationaire enkel-elektronen dichtheidsmatrix $n$ op door de tijdsafgeleide gelijk aan nul te stellen. Dit leidt tot een Sylvester-vergelijking van de vorm $AX - XB = Y$.
5. Wiskundige Afleiding: Met behulp van de algemene integraaloplossing voor de Sylvester-vergelijking via contourintegratie, drukken de auteurs de stationaire dichtheidsmatrix $n$ uit in termen van de geretardeerde ($G^r$) en geavanceerde ($G^a$) Green's functies. Deze expressie wordt vervolgens in de definitie van de stroom ingevuld om een stroomformule via Green's functies af te leiden.

Belangrijkste Bijdragen en Resultaten
Het artikel leidt een expliciete expressie af voor de elektrische stroom die van de linker elektrode naar het systeem vloeit binnen het GKSL-formalisme. De resulterende stroom $I_L$ wordt gegeven door:
$$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]$$
Deze formulering maakt een directe vergelijking mogelijk met de standaard Landauer-formule:
$$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))$$

De auteurs identificeren twee specifieke condities waaronder de via GKSL afgeleide formule samenvalt met de Landauer-formule:

1. Diagonale Hamiltoniaan (Onafhankelijke Moden): Indien de systeemhamiliaan diagonaal is (wat onafhankelijke multichannel-moden vertegenwoordigt), reduceert de GKSL-stroom tot een vorm die identiek is aan de Landauer-formule in het limiet van zwakke koppeling.
2. Energie-onafhankelijke Fermi-distributies: Indien de Fermi-distributiefuncties $f_L$ en $f_R$ als energie-onafhankelijke constanten kunnen worden behandeld binnen het relevante energiebereik (zoals geïllustreerd in het schematische model van het systeem), reduceren zowel de Landauer- als de GKSL-expressies tot dezelfde vereenvoudigde vorm:
   $$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]$$
   De auteurs merken op dat de conditie van energie-onafhankelijke distributies fysiek gerechtvaardigd is wanneer alle enkel-deeltje orbitaalenergieën strikt tussen de chemische potentialen van de elektroden liggen ($\mu_R < \epsilon_i < \mu_L$) en de temperatuur voldoende laag is ($k_B T \ll \mu_L - \mu_R$), zodanig dat $f_L \approx 1$ en $f_R \approx 0$.

Betekenis
Dit artikel biedt een rigoureuze afleiding die de GKSL-meestervergelijking-benadering verbindt met de Landauer-formule voor niet-interagerende systemen. Door gebruik te maken van de algemene oplossing van de Sylvester-vergelijking, tonen de auteurs aan dat de twee kaders niet inherent tegenstrijdig zijn, maar equivalent worden onder specifieke fysieke restricties. Specifiek geldt de equivalentie wanneer de energieafhankelijkheid van de Fermi-distributies verwaarloosd kan worden binnen het relevante energievenster. Dit werk verheldert de theoretische randvoorwaarden waaronder het dissipatieve, sequentiële tunnelingsbeeld van de meestervergelijking convergeert met het coherente transportbeeld van het Landauer-formalisme.