린드블라드 방정식

편집 역사 토론
양자역학 계산물리 물리 마지막 수정: 2026-08-28 04:12:41

1. 개요[편집]

린드블라드 방정식
Lindblad Equation (GKSL)
정식 명칭GKSL 방정식 (고리니-코사코프스키-수다르샨-린드블라드)
발표1976년 (두 팀이 독립적으로)
다루는 대상열린 양자계의 밀도행렬
강제 조건완전양성 · 대각합 보존 · 마르코프 반군
도약연산자 개수최대 d² − 1 개
수치 비용초연산자 d² × d² 또는 궤도 M 벌 × d

닫힌 계는 슈뢰딩거가 풀어 놨다. 문제는 세상에 닫힌 계가 없다는 것이다.

린드블라드 방정식(Lindblad equation, 정식으로는 GKSL 방정식)은 열린 양자계의 밀도행렬이 따를 수 있는 가장 일반적인 마르코프 시간 발전 방정식으로,

dρdt  =  −iℏ[H,ρ]  +  ∑kγk ⁣(LkρLk†−12{Lk†Lk, ρ})\frac{d\rho}{dt} \;=\; -\frac{i}{\hbar}[H,\rho] \;+\; \sum_k \gamma_k\!\left( L_k \rho L_k^\dagger - \tfrac12\{L_k^\dagger L_k,\,\rho\}\right)

의 꼴을 갖는다. HH 는 (램 이동까지 포함한) 유효 해밀토니안, LkL_k 는 도약연산자(jump operator), γk≥0\gamma_k \ge 0 은 그 비율이다.1 1976년 고리니-코사코프스키-수다르샨이 유한 차원에서, 린드블라드가 분리가능 힐베르트 공간의 유계 생성자에 대해 각각 독립적으로 얻었다.2

이 문서의 핵심 주장은 하나다. 이 형태는 누가 고른 모형이 아니라 강제된 결과다. “완전양성이고 대각합을 보존하는 사상의 한 매개변수 반군”이라는 요구만 걸면 생성자가 저 꼴 말고는 없다. 그래서 “감쇠항을 그럴듯하게 하나 붙이자”는 식으로 마스터 방정식을 손으로 짜는 순간, 그것은 십중팔구 GKSL 형태를 벗어나 있고 확률이 음수인 상태를 뱉는다.

기초 개념(밀도행렬·결어긋남·리우빌 초연산자)은 밀도행렬 문서가 이미 정리해 두었으므로, 여기서는 왜 이 형태뿐인가 · 어떻게 유도되며 언제 깨지는가 · 어떻게 수치로 푸는가 세 가지에 집중한다.

2. 형태가 강제된다 — CPTP 반군의 생성자[편집]

시간 발전을 사상 Λt:ρ(0)↦ρ(t)\Lambda_t : \rho(0) \mapsto \rho(t) 로 보자. 물리적으로 요구할 것은 셋이다.

  • 대각합 보존(TP). 확률의 총합이 1로 유지된다.
  • 완전양성(CP). ρ\rho 를 양의 반정부호로 보내는 것만으로는 부족하다. 계가 아무 상관 없는 보조계와 얽혀 있을 수 있으므로 Λt⊗In\Lambda_t \otimes \mathcal{I}_n 이 모든 nn 에 대해 양성이어야 한다. 이 “완전”이 빠지면 전치 사상처럼 혼자서는 멀쩡한데 얽힘을 만나면 음의 고유값을 내는 괴물이 통과한다.
  • 마르코프 반군. Λt+s=ΛtΛs\Lambda_{t+s} = \Lambda_t \Lambda_s, Λ0=I\Lambda_0 = \mathcal{I}. 환경이 과거를 기억하지 않는다는 뜻이다.

이 셋을 만족하는 연속 반군의 생성자 L\mathcal{L}(ρ˙=Lρ\dot\rho = \mathcal{L}\rho)의 일반형을 구하는 것이 GKS와 린드블라드가 한 일이다. 결과는 트레이스 없는 정규직교 연산자 기저 {Fm}m=1d2−1\{F_m\}_{m=1}^{d^2-1} 에 대해

Lρ=−i[H,ρ]+∑m,n=1d2−1amn ⁣(FmρFn†−12{Fn†Fm,ρ})\mathcal{L}\rho = -i[H,\rho] + \sum_{m,n=1}^{d^2-1} a_{mn}\!\left(F_m \rho F_n^\dagger - \tfrac12\{F_n^\dagger F_m, \rho\}\right)

이고, a=(amn)a = (a_{mn}) 가 양의 반정부호 에르미트 행렬(코사코프스키 행렬)이라는 것이 CP의 정확한 내용이다. aa 를 대각화하면 a=u diag(γk) u†a = u\,\mathrm{diag}(\gamma_k)\,u^\dagger, Lk=∑mumkFmL_k = \sum_m u_{mk} F_m 으로 개요의 정준형이 나온다. 즉 “완전양성 = 어떤 행렬이 양의 반정부호” 라는 번역이 이 정리의 알맹이고, 이 구조 때문에 잡음 채널 추정이 반정부호 계획법 문제로 떨어진다.

여기서 실무자가 반드시 알아야 할 함정 하나. 도약연산자는 유일하지 않다.

  • Lk→∑jukjLjL_k \to \sum_j u_{kj} L_j (유니터리 uu) — 생성자가 그대로다.
  • Lk→Lk+ck1L_k \to L_k + c_k \mathbb{1} 과 H→H+12i∑kγk(ck∗Lk−ckLk†)H \to H + \frac{1}{2i}\sum_k \gamma_k (c_k^* L_k - c_k L_k^\dagger) 를 동시에 하면 역시 그대로다.

같은 물리에 서로 다른 {Lk}\{L_k\} 가 무수히 대응한다는 뜻이다. 뒤에 나올 양자 궤적이 “실제로 일어나는 일”이 아닌 이유가 바로 이것이다 — 어떤 도약을 볼지는 환경을 어떻게 측정하느냐(unravelling)에 달렸고, ρ\rho 만이 관측 가능하다.

3. 유도 — 보른-마르코프와 세속근사[편집]

미시적으로는 전체 해밀토니안 Htot=HS+HB+HIH_{\rm tot} = H_S + H_B + H_I 에서 출발해 환경 BB 를 지운다. 상호작용 그림에서 폰 노이만 방정식을 두 번 반복 대입하고 세 단계 근사를 먹인다.

  1. 보른 근사 — 결합이 약하고 환경이 크므로 ρtot(t)≈ρS(t)⊗ρB\rho_{\rm tot}(t) \approx \rho_S(t)\otimes\rho_B. 환경은 계가 뱉은 정보를 되돌려주지 않는다.
  2. 마르코프 근사 — 환경 상관함수의 감쇠시간 τB\tau_B 가 계의 완화시간 τR\tau_R 보다 훨씬 짧다(τB≪τR\tau_B \ll \tau_R). 적분 안의 ρS(s)\rho_S(s) 를 ρS(t)\rho_S(t) 로 바꾸고 적분 하한을 −∞-\infty 로 보낸다. 여기까지가 레드필드 방정식이다.
  3. 세속근사(secular approximation, RWA) — 계 고유진동수의 차 ∣ω−ω′∣|\omega - \omega'| 가 완화율보다 훨씬 클 때, ei(ω−ω′)te^{i(\omega-\omega')t} 로 빠르게 진동하는 교차항을 버린다.

3번을 하지 않은 레드필드 방정식은 완전양성이 보장되지 않는다. 실제로 점유수가 음수로 내려가거나 순수도가 1을 넘는 사고가 재현되며, 이게 “감쇠 넣었더니 확률이 −0.003 나옴” 버그의 표준 원인이다. 세속근사를 거치면 비로소 GKSL 형태가 되고, 비율은 환경 상관함수의 푸리에 변환 γ(ω)\gamma(\omega) 로 주어진다. 열평형 환경이면 KMS 조건 γ(−ω)=e−βℏωγ(ω)\gamma(-\omega) = e^{-\beta\hbar\omega}\gamma(\omega) 가 성립해 정상상태가 깁스 상태로 떨어진다 — 정준 앙상블이 동역학에서 자동으로 나오는 흔치 않은 장면이다.

그리고 이 유도가 깨지는 경우들:

  • 준축퇴. ∣ω−ω′∣|\omega-\omega'| 가 완화율과 비슷하면 세속근사가 근거를 잃는다. 부분 세속근사·통합 리우빌리안 같은 처방이 있다.
  • 구조화된 스펙트럼 밀도. 환경에 뾰족한 모드가 있으면 τB\tau_B 가 길어져 마르코프성이 무너진다. 시간-비국소 나카지마-츠반치히, TCL 전개, 또는 계층운동방정식(HEOM)으로 간다.
  • 강결합·저온. 보른 근사 자체가 죽는다.
  • 국소 vs 대역 린드블라드. 두 부분계가 서로 결합해 있을 때, 각 부분계의 맨 연산자에 소산자를 붙이는 “국소” 방식은 편하지만 평형에서 정상 열류가 흐르는 등 열역학 제2법칙을 위반할 수 있다. 반대로 전체 HH 의 고유연산자로 짜는 “대역” 방식은 열역학은 맞지만 내부 결합이 소산율보다 약할 때 틀린다. 둘 다 근사이고, 어느 쪽이 맞는지는 결합 세기와 소산율의 비교가 정한다.3

4. 2준위계 — 블로흐 방정식과 T1 · T2[편집]

린드블라드가 실무에서 가장 자주 쓰이는 자리는 큐비트다. H=−ℏω02σzH = -\frac{\hbar\omega_0}{2}\sigma_z 에 도약연산자 둘을 얹자.

  • L1=γ1 σ−L_1 = \sqrt{\gamma_1}\,\sigma_- — 들뜬 상태가 바닥으로 떨어진다. 에너지 완화.
  • L2=γϕ σzL_2 = \sqrt{\gamma_\phi}\,\sigma_z — 에너지는 그대로인데 위상만 무작위로 걷는다. 순수 위상소멸.

ρ\rho 를 블로흐 벡터 r=(⟨σx⟩,⟨σy⟩,⟨σz⟩)\mathbf{r} = (\langle\sigma_x\rangle,\langle\sigma_y\rangle,\langle\sigma_z\rangle) 로 옮기면(블로흐 구) 회전좌표계에서

x˙=Δy−xT2,y˙=−Δx−yT2,z˙=−z−z∞T1\dot x = \Delta y - \frac{x}{T_2},\qquad \dot y = -\Delta x - \frac{y}{T_2},\qquad \dot z = -\frac{z - z_\infty}{T_1}

라는 블로흐 방정식이 그대로 나온다. 계수를 맞춰 보면

1T1=γ1,1T2=12T1+2γϕ\frac{1}{T_1} = \gamma_1, \qquad \frac{1}{T_2} = \frac{1}{2T_1} + 2\gamma_\phi

이고, γϕ≥0\gamma_\phi \ge 0 이므로 T2≤2T1T_2 \le 2T_1 이라는 부등식이 즉시 따라온다. 순수 위상소멸이 전혀 없는 이상적인 진폭감쇠에서만 등호가 성립한다. 실험에서 T2>2T1T_2 > 2T_1 이 나왔다면 물리가 아니라 피팅이 틀린 것이다.4

유한 온도 환경이면 도약연산자가 γ(nˉ+1) σ−\sqrt{\gamma(\bar n+1)}\,\sigma_- 와 γnˉ σ+\sqrt{\gamma \bar n}\,\sigma_+ 두 벌이 되고, T1=1/[γ(2nˉ+1)]T_1 = 1/[\gamma(2\bar n+1)], z∞=1/(2nˉ+1)z_\infty = 1/(2\bar n+1) 로 자연스럽게 볼츠만 비율이 재현된다. 온도를 올리면 완화가 빨라지고 정상상태 편극이 줄어든다는, 물리적으로 당연한 결과다.

2준위계의 GKSL 우변을 2×2 복소 밀도행렬 위에서 그대로 조립해 RK4로 적분하고, 블로흐 성분 x=2Re ρ₀₁, y=−2Im ρ₀₁, z=ρ₀₀−ρ₁₁ 은 읽어서만 그린다. 라비 진동하는 벡터가 소산항 때문에 구 표면에서 안으로 말려들어 정상상태 |r|=0.0793(순도 0.5031)으로 수렴하는 동안, 소산만 켠 Ω=0 참조 적분은 해석해와 최대 6.8e−15 만 어긋나고 |Trρ−1| 은 3.1e−15(≈14 ulp) 를 넘지 않으며 ρ 의 두 고유값은 음수로 내려가지 않는다.

블로흐 구 위에서 보면 그림이 명확하다. T1T_1 은 벡터를 북극(또는 열평형 높이) 쪽으로 끌어올리는 세로 방향 힘이고, T2T_2 는 벡터를 축 쪽으로 납작하게 눌러 붙이는 가로 방향 수축이다. 구 표면(순수상태)에서 출발한 점은 내부(혼합상태)로 파고들며, 결국 한 점으로 수렴한다. 그 한 점이 리우빌 초연산자의 영공간이다.

5. 수치 해법 1 — 벡터화와 리우빌 초연산자[편집]

L\mathcal{L} 은 ρ\rho 에 대해 선형이므로, ρ\rho 를 열 방향으로 쌓아 벡터 vec(ρ)∈Cd2\mathrm{vec}(\rho)\in\mathbb{C}^{d^2} 로 펴면 그냥 선형 상미분방정식이 된다. 열 쌓기 규약에서 vec(AρB)=(BT⊗A) vec(ρ)\mathrm{vec}(A\rho B) = (B^{\mathsf T}\otimes A)\,\mathrm{vec}(\rho) 이므로(크로네커 곱)

L^=−i(1⊗H−HT⊗1)+∑kγk[Lˉk⊗Lk−12(1⊗Lk†Lk)−12((Lk†Lk)T⊗1)]\hat{\mathcal{L}} = -i\big(\mathbb{1}\otimes H - H^{\mathsf T}\otimes \mathbb{1}\big) + \sum_k \gamma_k\Big[\bar L_k \otimes L_k - \tfrac12\big(\mathbb{1}\otimes L_k^\dagger L_k\big) - \tfrac12\big((L_k^\dagger L_k)^{\mathsf T}\otimes \mathbb{1}\big)\Big]

이고, 해는 vec(ρ(t))=eL^t vec(ρ(0))\mathrm{vec}(\rho(t)) = e^{\hat{\mathcal{L}}t}\,\mathrm{vec}(\rho(0)) 다. 이 순간 문제가 통째로 표준 수치선형대수로 넘어간다.

하고 싶은 것수치 문제도구
시간 발전크기 d² 선형 ODE크릴로프 지수적분, 경직하면 암시적
정상상태L^\hat{\mathcal{L}} 의 영공간희소 LU, 시프트-역 아놀디
완화 시간리우빌 갭 (고유값 실수부)란초스 알고리즘 계열
선형응답 스펙트럼분해행렬 요소크릴로프 + 되풀이 우변

대가는 잔인하다. nn 큐비트면 d=2nd = 2^n, L^\hat{\mathcal{L}} 은 4n×4n4^n \times 4^n 이다. n=10n=10 이면 조밀 저장에 101210^{12} 개 복소수 — 논외다. 희소행렬 저장과 행렬-벡터 곱만 쓰는 반복법이 사실상 유일한 길이고, 그마저도 n≳12n \gtrsim 12 부터는 버겁다. 다만 정상상태를 정확히 하나 뽑아야 하거나 리우빌 스펙트럼 자체가 궁금하면 이 길밖에 없다. 궤적 방법으로는 스펙트럼이 안 나온다.

특수한 구조가 있으면 훨씬 싸진다. 해밀토니안과 도약연산자가 보손·페르미온 연산자에 대해 이차식이면 리우빌리안이 정확히 대각화된다(프로센의 제3양자화). 가우스 상태는 공분산 행렬 2N×2N2N\times2N 만 굴리면 되므로 지수 폭발이 아예 없다. 광역학의 선형 공동 문제 대부분이 여기 속한다.

6. 수치 해법 2 — 몬테카를로 파동함수 (양자 도약)[편집]

밀도행렬을 안 들고 다니는 우회로가 몬테카를로 파동함수(MCWF, 양자 도약/양자 궤적)다. 달리바르-카스탱-묄머, 둠-촐러-리치가 1992년에, 카마이클이 양자 궤적이라는 이름으로 거의 같은 시기에 내놓았다. 한 궤적은 이렇게 굴린다.

  1. 비에르미트 유효 해밀토니안 Heff=H−i2∑kγkLk†LkH_{\rm eff} = H - \frac{i}{2}\sum_k \gamma_k L_k^\dagger L_k 로 ∣ψ⟩|\psi\rangle 를 결정론적으로 발전시킨다. 노름이 서서히 줄어든다.
  2. 스텝 동안 줄어든 노름 δp=1−∥∣ψ⟩∥2\delta p = 1 - \||\psi\rangle\|^2 가 도약 확률이다. 난수를 뽑아 비교한다.
  3. 도약하면 채널 kk 를 γk⟨Lk†Lk⟩\gamma_k\langle L_k^\dagger L_k\rangle 에 비례해 고르고 ∣ψ⟩→Lk∣ψ⟩|\psi\rangle \to L_k|\psi\rangle 로 갈아치운 뒤 정규화한다. 도약하지 않으면 그냥 정규화한다.
  4. 궤적 MM 개의 ∣ψ⟩⟨ψ∣|\psi\rangle\langle\psi| 를 평균하면 린드블라드 해로 수렴한다.

비용 대비표가 이 방법의 존재 이유를 전부 설명한다.

항목리우빌 직접 적분MCWF
메모리d² (초연산자는 d⁴)d × (동시 궤적 수)
오차결정론적, 적분기 차수통계, M의 제곱근에 반비례
병렬화행렬 연산 수준궤적별로 완전 독립
정상상태·스펙트럼직접 얻음긴 시간 평균으로 우회

즉 정확도 요구가 낮고 차원이 크면 MCWF, 정확도 요구가 높고 차원이 작으면 벡터화다. 통계오차가 1/M1/\sqrt M 로만 줄기 때문에 소수점 넷째 자리를 원하면 궤적을 10410^4 벌 이상 굴려야 하고, 그 지점에서 d2d^2 를 그냥 푸는 게 싸질 수 있다. 손익분기점은 대략 d∼103d \sim 10^3 언저리지만 관측량 종류에 따라 크게 흔들린다. 중요도 표본추출이나 궤적 재가중으로 분산을 줄이는 기법도 있다.

한 가지 개념적 주의. 도약은 물리적 사건이 아니라 계산 장치다. 앞서 본 게이지 자유도 때문에 같은 ρ(t)\rho(t) 를 주는 언래블링이 무수히 많고(도약 대신 연속 확산으로 푸는 양자 상태 확산도 그중 하나), 어떤 것을 볼지는 환경을 어떤 기저로 측정하느냐가 정한다. 평균만이 물리다.

실무 도구는 사실상 QuTiP가 표준이다. mesolve(리우빌 직접), mcsolve(궤적), steadystate(영공간)가 위 세 갈래에 그대로 대응하고, 도약연산자 리스트만 바꿔 끼우는 구조라 모형 실험이 빠르다. 시간의존 계수·플로케 확장·HEOM까지 같은 인터페이스로 붙어 있다.

7. 자주 나는 사고[편집]

  • 소산자를 손으로 짰다가 CP 위반. ρ\rho 에 감쇠항을 “대충” 곱해 넣는 순간 코사코프스키 행렬의 양의 반정부호성이 깨진다. 증상은 음의 점유수 또는 Tr⁡ρ2>1\operatorname{Tr}\rho^2 > 1. 정준형 밖으로 나가지 마라.
  • 대각합이 새는 적분기. 룽게-쿠타법은 Tr⁡ρ=1\operatorname{Tr}\rho = 1 을 자동으로 지켜주지 않는다. 매 스텝 대각합을 감시하고, 새면 시간간격이 아니라 방정식을 의심한다.
  • 경직성. γ\gamma 가 ω0\omega_0 보다 몇 자릿수 작으면 계가 경직계가 된다. 회전좌표계로 옮겨 빠른 위상을 걷어내는 것이 첫 처방이다.
  • 마르코프가 아닌 계에 린드블라드를 쓰기. 비마르코프 계에서 억지로 γ\gamma 를 피팅하면 짧은 시간의 이차 거동(제논 영역)을 절대 못 맞춘다. 감쇠가 지수함수가 아니면 그건 모형이 틀렸다는 신호다.

8. 관련 문서[편집]

9. Footnotes[편집]

  1. 열린 계 이름 짓기의 혼란은 유서가 깊다. GKSL, 린드블라드, 린드블라디언, 마르코프 마스터 방정식, 양자 동역학 반군이 전부 같은 것을 가리킨다. 논문 첫 페이지에서 “Lindblad” 철자를 확인하는 것으로 저자의 국적을 짐작하는 사람도 있다는데, 아마 사실이 아닐 것이다. ↩

  2. 두 논문이 같은 해에 나왔고 가정이 미묘하게 다르다. GKS는 유한 차원 NN 준위계에서 완전양성 반군의 생성자를 분류했고, 린드블라드는 분리가능 힐베르트 공간에서 유계 생성자에 대해 같은 결론을 얻었다. 물리에서 흔히 쓰는 조화진동자의 aa 는 유계가 아니라서 엄밀히는 린드블라드 정리의 적용 범위 밖인데, 다들 아무렇지 않게 쓴다. 수학자가 보면 킹받을 일이지만 물리적으로는 잘 작동한다. ↩

  3. 국소/대역 린드블라드 논쟁은 “편한 근사가 열역학을 위반한다”는 것이 밝혀진 뒤로 양자 열기관 분야에서 꽤 뜨거웠다. 정상 열류가 온도차 없이 흐르는 결과가 나오면 논문이 아니라 영구기관 특허를 내야 하므로, 그 전에 소산자 구성부터 다시 본다. ↩

  4. 그래서 큐비트 논문에는 T2T_2 와 T2∗T_2^* 가 따로 나온다. T2∗T_2^* 는 자유 유도 감쇠에서 그냥 잰 값이라 준정적 주파수 요동(1/f 잡음)까지 섞여 있어 훨씬 짧고, 에코 펄스로 그 성분을 되감아 얻은 것이 T2T_2 다. 논문의 T2T_2 가 유난히 좋아 보이면 펄스 시퀀스를 먼저 확인하는 것이 국룰이다. ↩