임의 포텐셜에서 Bose‑Einstein 응축의 고전장 구축 및 비영점 온도에서의 사각 진동

초록

본 논문은 고전장 근사법을 최적화하여 비영점 온도에서 조화 진동 포텐셜에 가두어진 Bose‑Einstein 응축체의 사각(quadrupole) 집단 모드를 시뮬레이션한다. 초기 온도와 외부 교란의 대칭성에 따라 열 구름과 응축부의 진동 특성이 달라지는 것을 재현하고, m=0 모드에서는 온도가 상승함에 따라 응축부가 열 구름의 자연 진동수로 이동하고, m=2 모드에서는 반대 방향으로 이동한다는 실험 결과를 정성적으로 설명한다.

상세 분석

본 연구는 고전장(클래식 필드) 근사를 J. Phys. B 40, R1 (2007)에서 제시된 형태를 기반으로, JILA 실험(Jin et al., Phys. Rev. Lett. 78, 764 (1997))에서 관측된 온도 의존적 집단 모드 변화를 재현한다. 고전장은 매크로스코픽하게 점유된 모드들을 연산자 대신 복소수 진폭으로 대체함으로써, 전체 원자수를 보존하는 비선형 슈뢰딩거 방정식(즉, Gross‑Pitaevskii 형태)으로 기술된다. 여기서 핵심은 ‘컷오프 파라미터 k_max’를 격자 간격에 의해 정의하고, Fourier 변환을 이용해 모드 공간을 제한함으로써 고전장이 실제 물리적 시스템을 충분히 포괄하도록 하는 점이다.

고전장으로부터 얻은 일인자 밀도 행렬을 시간·공간 평균(코그리그 그레이닝)하여, 가장 큰 고유값에 대응하는 고유벡터를 응축 파동함수 ψ₀로 정의한다. 나머지 고유값들은 열 구름 밀도 ρ_T로 해석된다. 이 절차는 실험에서 관측 가능한 컬럼 밀도와 직접 비교 가능하도록 설계되었으며, Penrose‑Onsager 정의와 일치한다.

평형 상태는 초기 상태에 무작위 교란을 가한 뒤, 비선형 방정식을 장시간 진화시켜 에너지와 입자 수가 보존되는 마이크로캐노니컬 집합에 도달하도록 한다. 이후, 고전장에 포함된 모드들의 에너지 분포를 분석하면, 높은 에너지 모드들은 빠르게 진동하고 평균적으로 고전적 에너지 등분배(N_i(ε_i−μ)=k_BT)를 만족한다는 점을 확인한다. 이는 절단 에너지 이하의 모든 점유된 모드가 고전적 열역학에 의해 설명될 수 있음을 의미한다.

동적 응답을 조사하기 위해 외부 교란을 두 가지 대칭(m=0, m=2)으로 가하고, 교란 주파수를 다양하게 스캔한다. 시뮬레이션 결과는 저온에서는 열 구름이 응축부와 동조하여 동일한 진동수를 보이며, 온도가 상승함에 따라 열 구름이 자체 고유 주파수(즉, 비상호작용 가스의 트랩 주파수)에 가까워진다. 특히 m=0 모드에서는 온도 0.65 T_c 부근에서 응축부의 진동수가 급격히 상승해 열 구름 주파수와 일치하는 현상이 재현된다. 반면 m=2 모드에서는 응축부가 열 구름 주파수와 반대 방향으로 이동하며, 임계 온도에 가까워질수록 응축부가 열 구름 주파수에 수렴한다. 이러한 온도 의존적 주파수 이동과 감쇠율 증가는 기존 Zaremba‑Nikuni‑Griffin(ZNG) 이론이나 2차 양자장 이론과 정성적으로 일치하지만, 고전장 근사는 복잡한 두‑성분 모델을 사전에 가정하지 않고도 동일한 현상을 자연스럽게 도출한다는 점에서 의미가 크다.

또한, 고전장 근사는 단일 실현(realization) 수준에서 마이크로캐노니컬 통계에 해당하므로, 실험에서 관측되는 샘플‑투‑샘플 변동성을 직접 시뮬레이션할 수 있다. 이는 기존의 고전적 평균장(Mean‑Field) 접근이 다루기 어려운 비평형 열 구름의 동역학을 포착하는 데 유리하다.

요약하면, 본 논문은 (1) 고전장 근사의 수치 구현 방법을 상세히 제시하고, (2) 코그리그 그레이닝을 통한 응축‑열 구름 분리를 정량화하며, (3) 온도와 교란 대칭성에 따른 집단 모드의 주파수·감쇠 변화를 성공적으로 재현함으로써, 비영점 온도에서의 Bose‑Einstein 응축체 동역학을 이해하는 새로운 시각을 제공한다.