나노 규모 시스템의 열평형 도달 메커니즘

초록

본 연구는 스핀‑1/2 입자들로 이루어진 작은 시스템이 동일한 종류의 스핀 환경과 상호작용하면서, 초기 순수 상태에서 시작해 시간에 따라 어떻게 열평형에 도달하는지를 전산 시뮬레이션으로 조사한다. 시스템의 감소밀도행렬 고유값이 시간이 흐를수록 일정한 값으로 수렴하고, 엔트로피도 정착값에 이르는 것이 관찰되었다. 환경 온도가 낮을수록 정착 상태와 정준 분포 사이의 차이가 커지는 특징도 밝혀졌다.

상세 분석

이 논문은 양자 열역학에서 오래된 질문, 즉 “폐쇄된 순수 상태가 어떻게 통계적 혼합 상태를 생성하는가”에 대해 구체적인 수치 실험을 제공한다. 저자들은 전체 시스템을 2‑레벨 스핀(스핀‑1/2)들의 집합으로 모델링하고, 시스템 부분과 환경 부분을 각각 (N_S)와 (N_E)개의 스핀으로 구분한다. 초기 조건은 시스템이 임의의 순수 상태 (|\psi_S(0)\rangle)에 놓이고, 환경은 고전적인 캐노니컬 앙상블을 모방하도록 설계된 순수 상태 (|\phi_E(\beta)\rangle)로 준비된다. 여기서 (\beta=1/k_B T)는 온도 파라미터이며, 환경의 순수 상태는 (\exp(-\beta H_E/2)) 연산자를 무작위 벡터에 적용해 얻는다. 이는 “정준 순수 상태”(canonical typicality) 개념과 일치한다.

시간 전개는 전형적인 시간‑의존 슈뢰딩거 방정식 (i\hbar \partial_t |\Psi(t)\rangle = H |\Psi(t)\rangle)를 직접 수치 적분함으로써 수행된다. Hamiltonian은 시스템 자체 (H_S), 환경 자체 (H_E), 그리고 상호작용 (H_{SE})의 합으로 구성되며, 상호작용은 일반적인 Heisenberg 형태 (\sum_{i\in S, j\in E} J_{ij},\mathbf{\sigma}_i\cdot\mathbf{\sigma}j)를 채택한다. 저자들은 (J{ij})의 크기를 다양하게 조절해 약한 결합부터 강한 결합까지 폭넓은 파라미터 공간을 탐색한다.

핵심 관찰은 감소밀도행렬 (\rho_S(t)=\mathrm{Tr}_E |\Psi(t)\rangle\langle\Psi(t)|)의 고유값 ({\lambda_k(t)})가 시간에 따라 급격히 변동한 뒤, 일정한 정착값 ({\lambda_k^{\infty}})에 수렴한다는 점이다. 이는 엔트로피 (S(t)=-\mathrm{Tr}