Connection between the GKSL master equation and the Landauer formula

Resumen Técnico: Conexión entre la Ecuación Maestra GKSL y la Fórmula de Landauer

Planteamiento del Problema
El transporte de electrones en sistemas mesoscópicos, tales como uniones moleculares y puntos cuánticos, se analiza típicamente utilizando dos marcos distintos: el enfoque de la función de Green (adecuado para el transporte coherente) y el enfoque de la ecuación maestra (comúnmente utilizado en regímenes de acoplamiento débil que involucran disipación y desfase). Aunque ambos son ampliamente empleados, la conexión teórica entre ellos aún no se comprende plenamente. Un trabajo previo realizado por Jin et al. estableció un vínculo para modelos específicos de sistema-baño bajo límites de alta temperatura y alto potencial químico, mostrando convergencia con las fórmulas de Meir–Wingreen y Landauer. Sin embargo, se requiere una derivación de la fórmula de corriente dentro del formalismo de la ecuación maestra de Gorini–Kossakowski–Sudarshan–Lindblad (GKSL) para un sistema no interactuante general, así como una identificación precisa de las condiciones bajo las cuales esta reduce a la fórmula de Landauer, para esclarecer esta conexión bajo condiciones más amplias.

Metodología
Los autores analizan un sistema de electrones no interactuantes conectado a electrodos izquierdo (L) y derecho (R), caracterizados por funciones de distribución de Fermi. El análisis procede dentro del marco de la ecuación maestra GKSL bajo las siguientes suposiciones:

1. Aproximación de Markov y Aproximación de Banda Ancha (WBLA): El acoplamiento entre cada orbital y los electrodos es suficientemente débil, lo que permite que la distribución de Fermi se aproxime como independiente de la energía dentro de la ventana de energía relevante.
2. Hamiltoniano: El sistema se describe mediante un hamiltoniano electrónico general $\hat{H} = \sum_{ij} h_{ij}\hat{c}^\dagger_i \hat{c}_j$.
3. Dinámica: La evolución temporal de la matriz de densidad $\rho$ está gobernada por la ecuación GKSL, incorporando superoperadores de Lindblad que dan cuenta de la inyección y extracción de electrones de los electrodos. Los coeficientes $\alpha_{a,i}$ y $\beta_{a,i}$ se derivan de las fuerzas de acoplamiento $\gamma_{a,i}$ y las distribuciones de Fermi $f_{a,i}$.
4. Estado Estacionario: Los autores resuelven la matriz de densidad de un solo electrón en estado estacionario $n$ estableciendo la derivada temporal en cero. Esto conduce a una ecuación de Sylvester de la forma $AX - XB = Y$.
5. Derivación Matemática: Utilizando la solución integral general para la ecuación de Sylvester mediante integración de contorno, los autores expresan la matriz de densidad en estado estacionario $n$ en términos de las funciones de Green retardada ($G^r$) y avanzada ($G^a$). Esta expresión se sustituye luego en la definición de corriente para derivar una fórmula de corriente expresada a través de funciones de Green.

Contribuciones Clave y Resultados
El artículo deriva una expresión explícita para la corriente eléctrica que fluye desde el electrodo izquierdo hacia el sistema dentro del formalismo GKSL. La corriente resultante $I_L$ viene 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 formulación permite una comparación directa con la fórmula estándar 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))$$

Los autores identifican dos condiciones específicas bajo las cuales la fórmula derivada de GKSL coincide con la fórmula de Landauer:

1. Hamiltoniano Diagonal (Modos Independientes): Si el hamiltoniano del sistema es diagonal (que representa modos multicanal independientes), la corriente GKSL se reduce a una forma idéntica a la fórmula de Landauer en el límite de acoplamiento débil.
2. Distribuciones de Fermi Independientes de la Energía: Si las funciones de distribución de Fermi $f_L$ y $f_R$ pueden tratarse como constantes independientes de la energía dentro del rango de interés (como se ilustra en el esquema del modelo), tanto la expresión de Landauer como la de GKSL se reducen a la misma 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]$$
   Los autores señalan que la condición de distribuciones independientes de la energía está justificada físicamente cuando todas las energías de los orbitales de partícula única se encuentran estrictamente entre los potenciales químicos de los electrodos ($\mu_R < \epsilon_i < \mu_L$) y la temperatura es suficientemente baja ($k_B T \ll \mu_L - \mu_R$), de modo que $f_L \approx 1$ y $f_R \approx 0$.

Significancia
El artículo proporciona una derivación rigurosa que vincula el enfoque de la ecuación maestra GKSL con la fórmula de Landauer para sistemas no interactuantes. Al utilizar la solución general de la ecuación de Sylvester, los autores demuestran que los dos marcos no son inherentemente contradictorios, sino que se vuelven equivalentes bajo restricciones físicas específicas. Específicamente, la equivalencia se mantiene cuando la dependencia de la energía de las distribuciones de Fermi puede despreciarse dentro de la ventana de energía relevante. Este trabajo esclarece las condiciones de contorno teóricas donde la imagen de túnel secuencial y disipativa de la ecuación maestra converge con la imagen de transporte coherente del formalismo de Landauer.