사방형 링에서 비등방성 터널링을 갖는 보손 모델의 적분가능성

초록

저자들은 네 개의 웰이 원형으로 연결된 보손 시스템을 대상으로, 인접 웰 사이의 터널링 강도가 서로 다를 수 있는 일반화된 Bose‑Hubbard Hamiltonian을 제시한다. 양자 역학적 역산법(QISM)과 알제브라적 베트 앤즈 해법을 이용해 정확히 풀 수 있음을 보이고, 여러 개의 의사진공(pseudovacuum) 상태가 필요함을 수치적으로 확인한다. 특히 특정 베트 앤즈 해는 에너지 식이 단순한 닫힌 형태로 표현된다.

상세 분석

이 논문은 네 개의 웰이 원형(4‑well ring)으로 배열된 Bose‑Einstein condensate(BEC) 시스템을 기술하는 새로운 적분가능 모델을 제안한다. 기존의 Bose‑Hubbard 모델에서는 인접 웰 사이의 터널링 계수가 모두 동일하게 가정되는 경우가 많았지만, 여기서는 t₁₂, t₁₄, t₂₃, t₃₄이 각각 독립적인 실수 파라미터 α₁…α₄와 κ에 의해 t_{ij}=−κ α_i α_j 형태로 정의된다. 이 제약은 t₁₂ t₃₄ = t₂₃ t₁₄라는 관계만을 남기면서도 충분히 비등방성(tunneling anisotropy)을 허용한다.

양자 역학적 역산법(QISM)을 적용하기 위해 SU(2) R‑행렬을 시작점으로 삼고, L‑연산자를 L_{i,j}(u)=u I+η N_{i,j} A_{i,j}A_{i,j}†/(α_i²+α_j²)−η I 형태로 정의한다. 여기서 A_{i,j}=α_i a_i+α_j a_j, N_{i,j}=N_i+N_j이며, η는 자유 파라미터이다. 두 개의 L‑연산자를 결합해 T(u)=L_{1,3}(u+ω) L_{2,4}(u−ω) 형태의 모노드리 매트릭스를 구성하고, 그 트레이스 τ(u)=A(u)+D(u) 를 전이 행렬이라 한다. τ(u) 가 서로 다른 스펙트럼 파라미터에 대해 교환 가능함을 보이면, τ(u) 로부터 보존량이 무한히 존재함을 의미한다.

전이 행렬의 계수와 Hamiltonian(1) 사이의 관계를 정리하면
U=κ η⁴/4, μ=κ ω/η, t_{ij}=−κ α_i α_j 로 식별된다. 따라서 H는 τ(u) 의 선형 결합으로 표현될 수 있고, 에너지 스펙트럼은 λ(u) (τ(u)의 고유값) 를 통해 E=−κ λ(u)+… 로 얻어진다.

알제브라적 베트 앤즈 해법을 적용하기 위해 두 개의 상호 교환 연산자 Γ†{1,3}=α₃ a₁†−α₁ a₃†, Γ†{2,4}=α₄ a₂†−α₂ a₄† 를 도입한다. 이 연산자들은 A_{1,3}, A_{2,4} 와 교환하며, 각각 N_{1,3}, N_{2,4} 의 상승 연산자 역할을 한다. 따라서 (Γ†{1,3})^k (Γ†{2,4})^l |0⟩ (k+l≤N) 로 정의되는 다수의 의사진공 |φ_{k,l}⟩ 가 존재한다. 각 의사진공에 대해 B(u) 가 소멸하고 A(u), D(u) 가 고유값을 갖는다는 사실을 이용해, C(u) 를 반복 적용해 베트 앤즈 상태 |ψ_{k,l}⟩ 를 구성한다.

베트 앤즈 방정식은 두 종류로 나뉜다. (i) k+l=N 인 경우, 전이 행렬 고유값 λ_{k,l}(u) 가 단순히 (u+ω+kη)(u−ω+lη)+const 형태이며, 에너지는 E=U(l−k)²+(l−k)μ 로 닫힌 식을 갖는다. 이는 터널링 파라미터 t_{ij} 에 독립적인 특수한 고유상태를 의미한다. (ii) k+l<N 인 경우, λ_{k,l}(u) 에는 추가적인 Bethe roots {v_j} 가 등장하고, η²(v_i+ω+kη)(v_i−ω+lη)=∏_{j≠i}(v_i−v_j±η) 형태의 비선형 방정식을 만족한다.

또한 전이 행렬으로부터는 H와 총 입자수 N 외에 두 개의 추가 보존량 Q_{1,3}= (Γ†{1,3}Γ{1,3})/(α₁²+α₃²), Q_{2,4}= (Γ†{2,4}Γ{2,4})/(α₂²+α₄²) 를 정의할 수 있다. 이들은 서로와 H, N 모두와 교환하므로, 모델은 총 네 개의 독립 보존량을 갖어 완전한 적분가능성을 만족한다.

수치 검증에서는 N=1,2 경우에 대해 전체 힐베르트 공간을 직접 대각화하고, 베트 앤즈 해에서 얻은 스펙트럼과 일치함을 확인한다. 특히 N=1에서는 세 개의 의사진공이 필요하고, N=2에서는 여섯 개의 의사진공이 모두 사용되어야 전체 10개의 에너지 레벨을 복원한다. 이는 “모든 의사진공을 사용해야 완전성(completeness)이 확보된다”는 결론을 뒷받침한다.

전체적으로 이 연구는 비등방성 터널링을 허용하면서도 양자 역학적 역산법을 통해 정확히 풀 수 있는 네 웰 Bose‑Hubbard 모델을 제시하고, 다중 의사진공 구조와 추가 보존량을 통해 기존 모델보다 풍부한 대칭 구조와 해석적 접근성을 제공한다는 점에서 의미가 크다.