개방 양자계

양자역학 계산물리 물리 마지막 수정: 2026-08-29 04:13:26

1. 개요[편집]

개방 양자계 Open Quantum System
다른 이름	열린 양자계, 산일 양자계
상태 기술	밀도행렬 (상태벡터로는 불가)
시간 발전	비유니터리 — CPTP 사상
일반형	크라우스(연산자합) 표현
크라우스 연산자 개수	최대 d² 개
마르코프 극한	린드블라드 방정식

학부 양자역학은 우주에 계가 하나뿐이라고 가정하고 시작한다. 실험실에는 계 하나와 그 계를 끊임없이 건드리는 나머지 우주가 있다.

개방 양자계(open quantum system)는 자신보다 훨씬 큰 환경과 에너지·정보를 주고받는 양자계를 말하며, 그 시간 발전은 슈뢰딩거 방정식이 주는 유니터리 변환이 아니다. 큐비트 한 개, 공동 속 광자 모드, 색소 분자 하나 — 실제로 우리가 측정하는 대상은 예외 없이 여기에 속한다.

닫힌 계의 유니터리 진화 $|\psi(t)\rangle = e^{-iHt/\hbar}|\psi(0)\rangle$ 는 우아하지만 두 가지를 원리적으로 예측하지 못한다. 첫째, 들뜬 원자가 스스로 바닥으로 떨어지는 일이 없다(유니터리는 가역이다). 둘째, 순수상태는 영원히 순수상태로 남는다( $\operatorname{Tr}\rho^2$ 가 보존된다). 그런데 실험에서 원자는 붕괴하고 큐비트는 죽는다. 모형이 틀린 게 아니라 계의 경계를 잘못 그은 것이다.

이 문서는 그 경계를 다시 긋는 방법 — 계와 환경을 함께 놓고 환경을 지웠을 때 남는 구조 — 를 다룬다. 마르코프 극한에서 나오는 미분방정식의 유도와 수치해법은 린드블라드 방정식에, 비대각 성분이 죽는 물리적 의미는 결어긋남에 몰아 두었다.

2. 전체는 유니터리인데 부분은 아니다[편집]

출발점은 정직하다. 계 $S$ 와 환경 $E$ 를 합친 전체는 닫힌 계이므로 유니터리로 발전한다.

\rho_{SE}(t) = U(t)\,\big(\rho_S(0)\otimes\rho_E\big)\,U^\dagger(t), \qquad U = e^{-iH_{\rm tot}t/\hbar}

우리가 관측할 수 있는 것은 $S$ 뿐이므로 환경을 부분대각합으로 지운다.

\rho_S(t) = \operatorname{Tr}_E\!\left[U\,(\rho_S(0)\otimes\rho_E)\,U^\dagger\right] \;\equiv\; \mathcal{E}_t\big[\rho_S(0)\big]

$\mathcal{E}_t$ 를 양자 채널 또는 동역학 사상이라 부른다. 여기가 결정적이다 — $U$ 가 유니터리여도 부분대각합을 씌운 $\mathcal{E}_t$ 는 일반적으로 비유니터리다. 정보가 $S$ 에서 $E$ 로 새는 것이 수식으로는 저 대각합 기호 하나에 들어 있다.¹

그래서 개방 양자계에서는 상태벡터를 아예 포기해야 한다. 전체가 순수상태여도 $S$ 와 $E$ 가 얽히는 순간 $\rho_S$ 는 혼합상태가 되고, 혼합상태는 $|\psi\rangle$ 하나로 적히지 않는다(양자 얽힘). “계를 열었다”와 “밀도행렬을 써야 한다”는 같은 말이다.

3. 크라우스 표현 — 채널의 일반형[편집]

위 정의는 환경을 명시적으로 들고 다녀야 해서 계산에 못 쓴다. 환경 기저 $\{|e_k\rangle\}$ 를 잡고 대각합을 풀어 쓰면 환경이 사라진다. $\rho_E = |e_0\rangle\langle e_0|$ 인 경우(환경이 처음에 순수상태) $K_k \equiv \langle e_k|U|e_0\rangle$ 로 두면

\mathcal{E}(\rho) = \sum_k K_k\,\rho\,K_k^\dagger, \qquad \sum_k K_k^\dagger K_k = \mathbb{1}

이 된다. 이것이 크라우스 표현(Kraus representation) 또는 연산자합 표현이다. $K_k$ 는 계의 힐베르트 공간 위 연산자일 뿐이라 환경 차원이 아무리 커도 상관없고, 개수는 $d^2$ 개를 넘지 않는다( $d$ = 계의 차원).

거꾸로도 성립한다는 것이 이 이야기의 알맹이다. 물리적으로 허용되는 사상이 만족해야 할 조건은 둘뿐이다.

대각합 보존(TP): $\operatorname{Tr}\mathcal{E}(\rho)=1$ . 확률이 새지 않는다.
완전양성(CP): 임의의 보조계에 대해 $\mathcal{E}\otimes\mathcal{I}_n$ 이 양성. 계가 다른 것과 얽혀 있어도 음의 확률이 나오지 않는다.

이 둘(CPTP)을 만족하는 사상은 정확히 크라우스 형태를 갖는 것들이고, 또한 정확히 “보조계를 붙여 유니터리를 돌린 뒤 지운 것”들이다. 앞 명제가 최의 정리, 뒤 명제가 스타인스프링 확장정리다. 즉 채널·크라우스·환경 딸림(dilation)은 같은 물건의 세 얼굴이다. 실무적 함의가 크다 — 어떤 잡음 모형을 크라우스 꼴로 적을 수 있다면 그것을 실현하는 환경이 반드시 존재하고, 반대로 크라우스 꼴로 못 적으면 그 모형은 비물리적이다.

수치적으로 CP를 판정하는 표준 도구는 최 행렬(Choi matrix)이다. 최대얽힘상태 $|\Omega\rangle=\frac{1}{\sqrt d}\sum_i|ii\rangle$ 에 채널을 반쪽만 걸어

J(\mathcal{E}) = (\mathcal{I}\otimes\mathcal{E})\big[\,|\Omega\rangle\langle\Omega|\,\big]

를 만들면, $\mathcal{E}$ 가 CP인 것과 $J(\mathcal{E})\succeq 0$ 인 것이 동치다(최-야미오우코프스키 동형). 채널이 $d^2\times d^2$ 행렬 하나로 바뀌므로 “채널이 물리적인가”라는 질문이 고유값 문제로, 채널 추정·최적화가 반정부호 계획법으로 떨어진다. $J$ 의 랭크가 곧 필요한 크라우스 연산자 개수이며, 랭크 1이면 그 채널은 유니터리다.

4. 대표 채널 네 개[편집]

큐비트 채널 몇 개만 손에 익혀 두면 개방 양자계 논문의 절반이 읽힌다. 블로흐 구에서 공을 어떻게 찌그러뜨리는지로 외우는 것이 빠르다.

채널	크라우스 연산자	블로흐 벡터에 하는 일	물리적 정체
진폭감쇠	$\mathrm{diag}(1,\sqrt{1-\gamma})$ , $\sqrt{\gamma}\,\lvert0\rangle\langle1\rvert$	가로 $\sqrt{1-\gamma}$ 배, 세로 $z\to\gamma+(1-\gamma)z$	자발방출, $T_1$
위상감쇠	$\mathrm{diag}(1,\sqrt{1-\lambda})$ , $\sqrt{\lambda}\,\lvert1\rangle\langle1\rvert$	가로만 $\sqrt{1-\lambda}$ 배, 세로 불변	순수 위상소멸
탈분극	$\sqrt{1-3p/4}\,\mathbb{1}$ , $\sqrt{p/4}\,\sigma_{x,y,z}$	공 전체를 $(1-p)$ 배 등방 수축	”아무 일이나 일어남”
비트반전	$\sqrt{1-p}\,\mathbb{1}$ , $\sqrt{p}\,\sigma_x$	$x$ 축 둘레로 납작하게	고전 비트 오류의 양자판

진폭감쇠 줄에 이 분야에서 제일 자주 쓰이는 부등식이 숨어 있다. 가로 방향은 $\sqrt{1-\gamma}$ , 세로 방향은 $(1-\gamma)$ 로 줄어드니 가로가 세로의 제곱근 속도로 죽는다. 시간 상수로 옮기면 정확히 $T_2 \le 2T_1$ 이다. 자세한 것은 결어긋남 참고.

탈분극 채널은 양자 오류 정정 문헌의 기본 잡음 모형인데, 이유는 물리적 정확성이 아니라 파울리 연산자만으로 적히기 때문이다. 파울리 채널은 안정자 형식 안에서 고전적으로 효율적으로 추적된다(고테스만-닐 정리). 진짜 잡음을 파울리 채널로 억지로 비트는 파울리 트위얼링이 실험에서 표준 전처리가 된 것도 같은 이유다 — 물리를 잘 맞춰서가 아니라 시뮬레이션이 가능해져서다.

5. 마르코프 근사와 그 청구서[편집]

크라우스 표현은 한 시각의 사상을 완벽히 기술하지만 미분방정식이 아니다. $\{\mathcal{E}_t\}$ 를 미분방정식으로 적으려면 조건이 더 필요하다.

\mathcal{E}_{t+s} = \mathcal{E}_t\,\mathcal{E}_s

즉 반군 성질이다. 물리적으로는 “환경이 방금 받은 정보를 계에 되돌려주지 않는다”는 뜻이고, 미시적으로는 환경 상관함수의 감쇠시간 $\tau_B$ 가 계의 완화시간 $\tau_R$ 보다 훨씬 짧다는 요구( $\tau_B \ll \tau_R$ )로 나타난다. 이 조건 아래 생성자의 형태가 강제로 결정되는 것이 GKSL 정리이고, 그 유도·형태·수치해법 전부는 린드블라드 방정식 문서에 있다.

청구서를 정확히 적어 두자. 마르코프 근사를 쓰면 아래를 포기한다.

짧은 시간의 이차 거동. 진짜 붕괴는 $t\to0$ 에서 $1-(t/\tau)^2$ 로 시작한다(그래서 양자 제논 효과가 존재한다). 마르코프 해는 처음부터 지수함수라 이 영역을 원리적으로 못 맞춘다.
정보의 역류. 구조화된 환경(뾰족한 스펙트럼 밀도, 광결정 밴드갭, 저온 포논)에서는 계가 뱉은 결맞음이 되돌아온다. 마르코프 해는 단조 감쇠만 낼 수 있다.
초기 계-환경 상관. 유도는 $\rho_{SE}(0)=\rho_S\otimes\rho_E$ 를 깔고 시작한다. 강결합계에서 이미 상관이 있으면 사상이 CP가 아닐 수도 있고, 그러면 크라우스 표현 자체가 성립하지 않는다.

비마르코프성을 정량화하는 척도도 표준화되어 있다. 두 상태의 대각합 거리 $D(\rho_1,\rho_2)$ 는 CPTP 사상 아래 절대 늘지 않으므로, $\frac{d}{dt}D>0$ 인 구간의 넓이를 합치면 “구별가능성이 되돌아온 총량”이 된다(BLP 척도). 다른 계열은 사상의 분해 $\mathcal{E}_{t,0}=\mathcal{E}_{t,s}\mathcal{E}_{s,0}$ 에서 중간 사상 $\mathcal{E}_{t,s}$ 가 CP인지를 묻는다(CP-분해가능성, RHP 척도). 둘은 동치가 아니며 어떤 계는 한쪽 기준으로만 비마르코프다.²

6. 비마르코프 계를 수치로 푸는 법[편집]

마르코프가 깨지면 선택지는 셋이다.

1. 시간-비국소 방정식을 그대로 푼다. 나카지마-츠반치히 방정식은 정확하지만 기억핵 $\mathcal{K}(t-s)$ 에 대한 적분미분방정식이라 모든 과거를 들고 있어야 한다. 시간-국소 전개(TCL)는 미분방정식 꼴을 유지하는 대신 계수가 시간에 의존하고, 차수를 자르면 어느 구간부터 조용히 틀린다.

2. 환경 일부를 계로 편입한다. 뾰족한 환경 모드 몇 개를 계 쪽으로 끌어와 “확장된 계 + 마르코프 잔여 환경”으로 다시 나누는 처방이다(유사모드, 반응좌표 사상). 계 차원이 커지는 대신 다시 린드블라드를 쓸 수 있어서 실무에서 가성비가 좋다.

3. 환경을 진짜로 시뮬레이션한다. 계층운동방정식(HEOM)은 조화 환경 + 특정 스펙트럼 밀도라는 조건 아래 무한 계층을 잘라 수치적으로 정확한 답을 준다. 사슬 사상(TEDOPA)은 환경을 1차원 사슬로 바꿔 텐서 네트워크로 굴린다. 정확한 대신 비싸고, 계층 절단 차수나 결합차원이 곧 수렴 파라미터다.

메모리 벽은 어느 길로 가든 같다. 큐비트 $n$ 개면 $\rho$ 가 $4^n$ 개 복소수이고 초연산자는 $16^n$ 이다. 개방 계 $n$ 큐비트가 닫힌 계 $2n$ 큐비트와 같은 값을 부른다. 그래서 밀도행렬을 아예 들지 않는 몬테카를로 파동함수나, $\rho$ 를 텐서 네트워크로 압축하는 행렬곱 밀도연산자가 큰 계의 실질적 출구다.

7. 실무에서 나는 사고[편집]

감쇠항을 손으로 붙인다. $\dot\rho$ 우변에 $-\Gamma\rho$ 같은 항을 그냥 얹는 순간 CP가 깨지고, 증상은 며칠 뒤 음의 점유수나 $\operatorname{Tr}\rho^2>1$ 로 나타난다. 크라우스 또는 GKSL 정준형 밖으로 나가지 마라.
채널을 곱해 놓고 마르코프라고 부른다. 같은 크라우스 집합을 $N$ 번 반복 적용하는 것은 반군이지만, 실험에서 잰 $\mathcal{E}_{t_1},\mathcal{E}_{t_2}$ 가 서로 합성되는지는 별도로 확인해야 한다. 확인 없이 지수 피팅부터 하는 것이 국룰인데, 그래서 잔차가 이상한 논문이 나온다.
환경 온도를 잊는다. 유한 온도면 올라가는 도약( $\sigma_+$ )도 있어야 하고, 없으면 정상상태가 절대영도로 수렴한다. 상온 실험을 $T=0$ 모형으로 맞추려다 결합상수를 비물리적으로 키우는 사고가 흔하다.
환경 차원을 유한하게 잘라 놓고 “정확히 풀었다”고 한다. 조화 진동자 모드를 $n_{\max}$ 에서 자르면 강한 구동에서 반드시 새고, 증상은 에너지가 슬금슬금 사라지는 것이다. 절단 차원을 두 배로 올렸을 때 결과가 안 움직이는지 보는 것이 유일한 방어다.³

8. 관련 문서[편집]

정보가 “샌다”는 표현은 은유가 아니다. 전체가 순수상태로 남는 한 폰 노이만 엔트로피의 총합은 0이고, $S(\rho_S)$ 가 늘어난 만큼이 정확히 $S$ 와 $E$ 사이 얽힘이다. 즉 계의 엔트로피 증가는 사라진 정보가 아니라 주소가 환경으로 옮겨간 정보다. 물론 환경에서 되찾아 올 방법이 없다는 점에서 실용적으로는 사라진 것과 같다. 세금과 비슷하다. ↩
“이 계는 비마르코프다”라는 문장은 척도를 명시하지 않으면 정보량이 거의 없다. BLP로 재면 마르코프인데 RHP로 재면 비마르코프인 모형을 만드는 것은 연습문제 수준이라, 리뷰어가 제일 먼저 묻는 것도 어느 척도냐다. ↩
이건 개방 계 전용 사고가 아니라 절단이 들어가는 모든 수치해석의 공통 죄목이다. 다만 여기서는 증상이 “발산”이 아니라 “그럴듯한데 틀린 정상상태”로 나와서 훨씬 늦게 들킨다. 발산은 친절한 버그다. ↩