[논문 리뷰] Quantum Friction: Cooling Quantum Systems with Unitary Time Evolution

Aurel Bulgac, Michael McNeil Forbes외 2명|arXiv (Cornell University)|2013. 05. 29.

Advanced Thermodynamics and Statistical Mechanics참고 문헌 4인용 수 5

한 줄 요약

이 논문은 하미르토니안에 시간역행 대칭을 깨는 퍼텐셜을 추가함으로써 양자 다체계를 냉각시키는 보존적이고 기억이 없는 국소 양자 마찰을 소개한다. 이 방법은 $O(N^2ackepsilon N)$ 스케일링을 통해 효율적인 실시간 진동을 가능하게 하여, $N \sim 10^5$에서 $10^6$ 개의 상태를 가진 대규모 페르미온계에서 기존의 $O(N^3)$ 대각화 방법보다 훨씬 빠른 기초 상태 준비를 가능하게 하며, 시간에 따라 변화하는 밀도기반 함수이론(TDDFT) 시뮬레이션의 계산적 병목 현상을 극복한다.

ABSTRACT

We introduce a type of quantum dissipation -- local quantum friction -- by adding to the Hamiltonian a local potential that breaks time-reversal invariance so as to cool the system. Unlike the Kossakowski-Lindblad master equation, local quantum friction directly effects unitary evolution of the wavefunctions rather than the density matrix: it may thus be used to cool fermionic many-body systems with thousands of wavefunctions that must remain orthogonal. In addition to providing an efficient way to simulate quantum dissipation and non-equilibrium dynamics, local quantum friction coupled with adiabatic state preparation significantly speeds up many-body simulations, making the solution of the time-dependent Schrödinger equation significantly simpler than the solution of its stationary counterpart.

연구 동기 및 목표

• 강한 상관성이 있는 페르미온계의 대규모 시간에 따라 변화하는 밀도기반 함수이론(TDDFT) 시뮬레이션을 위한 초기 상태 준비에서 발생하는 계산적 병목 현상을 해결하기 위해.
• 단일 입자 상태 수 $N \sim 10^5$에서 $10^6$에 이르는 시스템에서 자가일관된 기초 상태를 위한 반복적 대각화의 금방 드러나는 $O(N^3)$ 비용을 극복하기 위해.
• 밀도 행렬을 진화시키지 않고 양자 소산을 모델링하기 위해 코사코프스키-린드블라드 마스터 방정식의 대안으로 기억이 없는 보존적 방법을 개발하기 위해.
• 시간에 따라 변화하는 슈뢰딩거 방정식의 해를 구하는 데 있어 정적 해법보다 단순화된 방법을 통해 핵물리학, 초저온 원자, 전자 구조의 비평형 동역학을 효율적으로 시뮬레이션하기 위해.
• 아디아바틱 스위칭과 함께 국소 양자 마찰을 사용하여 3차원 대규모 시스템(예: $16 \times 16 \times 64$ 격자)의 기초 상태를 준비하는 것이 가능함을 보여주기 위해.

제안 방법

• 시간역행 대칭을 직접적으로 깨는 국소 퍼텐셜 $U_t$ 를 도입하여 파동함수의 보존적 진동을 통해 냉각을 유도하기 위해.
• 시간에 따라 변화하는 슈뢰딩거 방정식에 국소 양자 마찰 퍼텐셜을 적용하여, 밀도 행렬 진화 없이도 다체 상태의 정규직교성을 유지하기 위해.
• 국소 양자 마찰과 아디아바틱 상태 준비를 조합하여 간단한 초기 상태에서 목표 기초 상태로 점진적으로 진화하기 위해.
• 수정된 슈뢰딩거 방정식을 사용한 실시간 진동을 통해 TDDFT의 초기 조건을 준비함으로써, 큰 $N \times N$ 하미르토니안의 자가일관된 대각화가 필요 없도록 하기 위해.
• GPU 가속된 연립 편미분방정식(PDE) 해법을 사용하여 대규모 계산 플랫폼(예: UW Hyak 및 NCCS의 Titan)에 이 방법을 구현하기 위해.
• 32,768 × 32,768 단일 입자 하미르토니안과 15,322개의 연립 PDE를 가진 3차원 단위 페르미 가스에서 이 방법을 시험하여, 24개의 GPU에서 약 1시간 내에 수렴을 달성함.

실험 결과

연구 질문

• RQ1 밀도 행렬을 진화시키지 않고도 강한 상관성이 있는 페르미온계에서 양자 소산을 시뮬레이션할 수 있는 보존적이고 기억이 없는 방법을 개발할 수 있는가?
• RQ2 시간역행 대칭을 깨는 국소 퍼텐셜의 포함이 다체 파동함수의 효율적 냉각을 가능하게 하며, 정규직교성을 유지하는가?
• RQ3 반복적 대각화에 비해 국소 양자 마찰이 TDDFT 시뮬레이션의 초기 상태 준비에 있어 계산 비용을 크게 줄일 수 있는가?
• RQ4 국소 양자 마찰을 사용한 실시간 진동이 대규모 시스템의 시간에 따라 변화하는 슈뢰딩거 방정식을 해결하는 데 있어 정적 방법보다 얼마나 뛰어나게 성능을 발휘하는가?
• RQ5 이 방법은 중간핵이나 초저온 페르미 가스와 같은 $N \sim 10^5$에서 $10^6$ 개의 단일 입자 상태를 가진 비평형 동역학을 시뮬레이션하기 위해 스케일업 가능한가?

주요 결과

• 국소 양자 마찰은 밀도 행렬 진화 없이도 시간역행 대칭을 깨는 퍼텐셜을 통해 페르미온 다체계의 보존적 냉각을 가능하게 하며, 메모리 집약적인 마스터 방정식을 피한다.
• 반복적 대각화를 실시간 진동으로 대체함으로써 기초 상태 준비의 계산 비용을 $O(N^3)$에서 $O(N^2\backepsilon N)$으로 감소시킨다.
• UW Hyak 클러스터의 24개 GPU를 사용하여, $16 \times 16 \times 64$ 격자점과 $32,768$ 개의 단일 입자 상태를 가진 3차원 단위 페르미 가스의 기초 상태를 약 1시간 내에 준비하였다.
• 다양한 격자 크기($24^3$, $32^3$, $40^3$)에서 적외선 수렴성이 확인되어 수치적 안정성과 확장 가능성을 입증하였다.
• 아디아바틱 스위칭과 국소 양자 마찰의 조합은 순수 아디아바틱 진동에 비해 스위칭 시간을 극적으로 줄여 수렴 속도를 가속화한다.
• 이 방법은 향후 타이탄(NCCS)과 같은 엑사스케일 시스템(18,680개의 GPU)으로도 확장 가능하며, $100^3$ 격자에서 수백만 개의 연립 PDE를 시뮬레이션할 수 있다.