Connection between the GKSL master equation and the Landauer formula

Resumo Técnico: Conexão entre a Equação Mestra GKSL e a Fórmula de Landauer

Definição do Problema
O transporte de elétrons em sistemas mesoscópicos, como junções moleculares e pontos quânticos, é tipicamente analisado utilizando dois frameworks distintos: a abordagem de função de Green (bem adequada para transporte coerente) e a abordagem de equação mestra (comumente usada em regimes de acoplamento fraco envolvendo dissipação e desfasagem). Embora ambos sejam amplamente empregados, a conexão teórica entre eles permanece não totalmente compreendida. Trabalhos anteriores por Jin et al. estabeleceram um vínculo para modelos específicos de sistema-banho sob limites de alta temperatura e alto potencial químico, mostrando convergência para as fórmulas de Meir–Wingreen e Landauer. No entanto, é necessária uma derivação da fórmula de corrente dentro do formalismo da equação mestra de Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) para um sistema não interagente geral, e uma identificação precisa das condições sob as quais ela se reduz à fórmula de Landauer, para esclarecer essa conexão sob condições mais amplas.

Metodologia
Os autores analisam um sistema de elétrons não interagentes conectado a eletrodos esquerdo (L) e direito (R), caracterizados por funções de distribuição de Fermi. A análise procede dentro do framework da equação mestra GKSL sob as seguintes suposições:

1. Aproximação de Markov e Aproximação de Banda Larga (WBLA): O acoplamento entre cada orbital e os eletrodos é suficientemente fraco, permitindo que a distribuição de Fermi seja aproximada como independente da energia dentro da janela de energia relevante.
2. Hamiltoniano: O sistema é descrito por um Hamiltoniano eletrônico geral $\hat{H} = \sum_{ij} \hat{c}^\dagger_i \hat{c}_j$.
3. Dinâmica: A evolução temporal da matriz de densidade $\rho$ é governada pela equação GKSL, incorporando superoperadores de Lindblad que contabilizam a injeção e extração de elétrons dos eletrodos. Os coeficientes $\alpha_{a,i}$ e $\beta_{a,i}$ são derivados das forças de acoplamento $\gamma_{a,i}$ e das distribuições de Fermi $f_{a,i}$.
4. Estado Estacionário: Os autores resolvem para a matriz de densidade de um único elétron em estado estacionário $n$ definindo a derivada temporal como zero. Isso leva a uma equação de Sylvester da forma $AX - XB = Y$.
5. Derivação Matemática: Usando a solução integral geral para a equação de Sylvester envolvendo integração de contorno, os autores expressam a matriz de densidade de estado estacionário $n$ em termos das funções de Green retardada ($G^r$) e avançada ($G^a$). Esta expressão é então substituída na definição de corrente para derivar uma fórmula de corrente expressa via funções de Green.

Principais Contribuições e Resultados
O artigo deriva uma expressão explícita para a corrente elétrica fluindo do eletrodo esquerdo para o sistema dentro do formalismo GKSL. A corrente resultante $I_L$ é dada por:
$$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]$$
Esta formulação permite uma comparação direta com a fórmula padrão de Landauer:
$$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))$$

Os autores identificam duas condições específicas sob as quais a fórmula derivada de GKSL coincide com a fórmula de Landauer:

1. Hamiltoniano Diagonal (Modos Independentes): Se o Hamiltoniano do sistema for diagonal (representando modos multicanais independentes), a corrente GKSL reduz-se a uma forma idêntica à fórmula de Landauer no limite de acoplamento fraco.
2. Distribuições de Fermi Independentes de Energia: Se as funções de distribuição de Fermi $f_L$ e $f_R$ puderem ser tratadas como constantes independentes de energia dentro da faixa de interesse (conforme ilustrado no modelo esquemático), tanto as expressões de Landauer quanto as de GKSL reduzem-se à mesma forma simplificada:
   $$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]$$
   Os autores observam que a condição de distribuições independentes de energia é fisicamente justificada quando todas as energias orbitais de partícula única situam-se estritamente entre os potenciais químicos dos eletrodos ($\mu_R < \epsilon_i < \mu_L$) e a temperatura é suficientemente baixa ($k_B T \ll \mu_L - \mu_R$), de modo que $f_L \approx 1$ e $f_R \approx 0$.

Significância
O artigo fornece uma derivação rigorosa ligando a abordagem da equação mestra GKSL à fórmula de Landauer para sistemas não interagentes. Ao utilizar a solução geral para a equação de Sylvester, os autores demonstram que os dois frameworks não são inerentemente contraditórios, mas tornam-se equivalentes sob restrições físicas específicas. Especificamente, a equivalência mantém-se quando a dependência de energia das distribuições de Fermi pode ser negligenciada dentro da janela de energia relevante. Este trabalho esclarece as condições de contorno teóricas onde a imagem de tunelamento sequencial dissipativo da equação mestra converge com a imagem de transporte coerente do formalismo de Landauer.