Connection between the GKSL master equation and the Landauer formula

기술 요약: GKSL 마스터 방정식과 란다우어 공식 간의 연결성

문제 제기
분자 접합이나 양자점과 같은 메조스코픽 시스템에서의 전자 수송은 일반적으로 두 가지 구별된 프레임워크를 사용하여 분석된다: 그린 함수 접근법(결맞음 수송에 적합)과 마스터 방정식 접근법(소산 및 결맞음 해제와 관련된 약한 결합 영역에서 흔히 사용됨). 두 방법 모두 널리 사용되지만, 이들 사이의 이론적 연결은 여전히 완전히 이해되지 않은 상태이다. Jin 등의 선행 연구는 고온 및 큰 화학 퍼텐셜 한계 하의 특정 시스템-배스 모델에 대해 연결성을 확립하여, Meir–Winggin 및 란다우어 공식으로의 수렴을 보여준 바 있다. 그러나 일반적인 비상호작용 시스템에 대하여 GKSL(Gorini–Kossakowski–Sudarshan–Lindblad) 마스터 방정식 형식 내에서 전류 공식을 유도하고, 이 식이 란다우어 공식으로 환원되는 조건을 정확히 식별하는 과정이 필요하며, 이는 더 넓은 조건 하에서 이 연결을 명확히 하기 위함이다.

방법론
저자들은 페르미 분포 함수로 특징지어지는 좌측(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. 정상 상태: 저자들은 시간 미분 항을 0으로 설정함으로써 정상 상태 단일 전자 밀도 행렬 $n$을 구한다. 이는 $AX - XB = Y$ 형태의 실베스터 방정식(Sylvester equation)으로 이어진다.
5. 수학적 유도: 경로 적분을 포함하는 실베스터 방정식의 일반적인 적분 해를 사용하여, 저자들은 후퇴(retarded, $G^r$) 및 전진(advanced, $G^a$) 그린 함수를 통해 정상 상태 밀도 행렬 $n$을 표현한다. 이 표현식은 전류 정의에 대입되어 그린 함수를 통해 표현된 전류 공식을 도출한다.

주요 기여 및 결과
본 논문은 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로부터 유도된 공식이 란다우어 공식과 일치하게 되는 두 가지 특정 조건을 식별하였다:

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 마스터 방정식을 란다우어 공식과 연결하는 엄밀한 유도를 제공한다. 저자들은 실베스터 방정식의 일반 해를 활용함으로써, 두 프레임워크가 본질적으로 모순되는 것이 아니라 특정 물리적 제약 조건 하에서 동등해짐을 입증하였다. 구체적으로, 이 등가성은 관련 에너지 창 내에서 페르미 분포의 에너지 의존성을 무시할 수 있을 때 성립한다. 이 연구는 소산적이고 순차적인 터널링 묘사인 마스터 방정식이 결맞는 수송 묘사인 란다우어 형식론과 수렴하는 이론적 경계 조건을 명확히 한다.