돌리 모델 회전파 근사의 정량적 평가

초록

본 논문은 다중 원자 양자 라비 모델인 돌리 모델에 대해 회전파 근사(RWA)의 정확성을 비교적·비교론적으로 비정상적인 경계값을 제공한다. 저자들은 Dicke Hamiltonian과 그 RWA 버전인 Tavis‑Cummings Hamiltonian 사이의 시간 진화 차이를 연산자 노름으로 제한하고, 이 경계가 초기 상태의 총 각운동량 S와 전체 여기수 M, 그리고 결합 강도 λ/ω에 어떻게 의존하는지를 분석한다. 결과는 상태‑의존적이며, 특히 큰 원자 수 N 혹은 열역학적 극한에서도 RWA가 어느 정도 유지되는 조건을 제시한다. 수치 실험은 이론적 경계와 좋은 일치를 보인다.

상세 분석

논문은 먼저 Dicke 모델과 Tavis‑Cummings 모델을 정확히 정의하고, 두 Hamiltonian을 자유 Hamiltonian H₀에 대한 상호작용 그림으로 변환한다. 이 과정에서 H₁(t)와 H₂(t)라는 시간‑의존 연산자를 도입하고, 각각의 유니터리 전이 연산자 U₁(t), U₂(t)를 정의한다. 핵심 정리는 “정리 1”으로, 임의의 초기 상태 ψ∈D₀(유한 광자 수 공간)에 대해 ‖U₁(t)ψ − U₂(t)ψ‖의 상한을 명시적으로 구한다. 상한식은 두 항으로 구성된다. 첫 번째 항은 λ/ω·|sin ωt|·O(√M · S²) 형태이며, 두 번째 항은 λ²/ω·|t|·O(S^{7/2} · M) 형태이다. 여기서 S는 전체 각운동량, M은 전체 여기수(광자 + 스핀 여기)이다. 이 식은 다음과 같은 물리적 의미를 갖는다. (1) 결합 강도 g=λ/ω가 작을수록, 그리고 시스템이 공명(ω=ω₀)일 때 RWA가 더 정확해진다. (2) 광자 수가 많을수록 첫 번째 항이 커져 오차가 증가하지만, 시간에 따라 진동하는 sin ωt 요인 때문에 짧은 시간 구간에서는 억제될 수 있다. (3) 총 각운동량 S가 클수록 두 번째 항이 급격히 커지므로, 다수 원자(N≫1) 시스템에서는 RWA가 제한적일 수 있다. 특히 S∝N/2이므로, N이 커질수록 λ²·N^{7/2}·M/ω·|t| 형태의 오차가 지배적이다. 저자들은 이 결과를 이용해 열역학적 극한(N→∞)에서 RWA가 유지되기 위한 조건을 λ·N^{7/4}≪√ω·|t|^{-1/2} 등으로 제시한다.

수학적 증명은 두 단계로 나뉜다. 첫 번째는 H₁(t)와 H₂(t) 사이의 차이를 직접 계산하고, 그 차이를 시간 적분 형태로 표현한다. 두 번째는 베타 함수와 코시-슈바르츠 부등식을 이용해 연산자 노름을 추정한다. 중요한 점은 D₀가 H₁(t), H₂(t)에 대해 불변 하위공간이라는 사실을 이용해, 무한 차원 Fock 공간 전체가 아니라 유한 차원 부분에만 집중함으로써 계산을 간소화했다는 것이다. 또한, Dicke 상태 |S,m⟩와 Tavis‑Cummings 고유 상태 ψ_{σ,S,M} 사이의 관계를 Bethe 방정식(26)을 통해 명시적으로 연결함으로써, 실제 수치 계산이 가능한 형태로 전환한다.

수치 실험에서는 N=2, 3, 5 등 소규모 원자 수와 다양한 초기 상태(예: |S=S_max, m=−S⟩⊗|n⟩)에 대해 정확한 시간 진화와 Tavis‑Cummings 근사 진화를 비교한다. 그래프는 이론적 상한이 실제 오차를 충분히 포괄함을 보여준다. 특히, 큰 M(광자 수)에서는 첫 번째 항이 지배적이며, 짧은 시간 구간에서 오차가 거의 0에 가까워진다. 반면, 긴 시간 또는 큰 S 경우에는 두 번째 항이 우세해져 오차가 선형적으로 증가한다.

결론적으로, 논문은 Dicke 모델에 대한 RWA의 정확성을 정량적으로 평가하는 최초의 비교론적 경계값을 제공한다. 이 경계는 상태‑의존적이며, 총 각운동량과 전체 여기수가 핵심 파라미터임을 밝힌다. 또한, 다원자 시스템에서 RWA가 언제 유효한지, 그리고 열역학적 극한에서 어떤 스케일링이 필요한지를 명확히 제시함으로써, 실험적 구현(예: 초전도 회로, 양자점 어레이)이나 이론적 모델링에 중요한 지침을 제공한다.