[논문 리뷰] Third virial coefficient of the unitary Bose gas

Yvan Castin, Félix Werner|arXiv (Cornell University)|2012. 12. 21.

Cold Atom Physics and Bose-Einstein Condensates참고 문헌 46인용 수 12

한 줄 요약

이 논문은 조화 퍼텐셜 내에서의 정확한 삼체 문제 해법과 역 잔여법을 사용하여, 공간적으로 균일한 단위 보스 가스의 삼차 비열계수 $a_3$와 삼차 누적량 $b_3$를 해석적으로 결정한다. 삼보스의 질량 각도 $\nu = \pi/6$를 작은 매개변수로 간주함으로써, 저온에서의 에피모프 트리머와 고온에서의 이온화 효과를 모두 고려한 고정밀도 해석적 근사치를 얻는다. 이에 따라 $b_3 = 3\sqrt{3}[S(0^+) + C]$ 형태로 표현되며, 여기서 $C \approx 0.648$는 적분 표현으로부터 유도된 상수로, 온도에 관계없이 $a_3$의 정확한 예측을 가능하게 한다.

ABSTRACT

By unitary Bose gas we mean a system composed of spinless bosons with s-wave interaction of infinite scattering length and almost negligible (real or effective) range. Experiments are currently trying to realize it with cold atoms. From the analytic solution of the three-body problem in a harmonic potential, and using methods previously developed for fermions, we determine the third cumulant (or cluster integral) b_3 and the third virial coefficient a_3 of this gas, in the spatially homogeneous case, as a function of its temperature and the three-body parameter R_t characterizing the Efimov effect. A key point is that, converting series into integrals (by an inverse residue method), and using an unexpected small parameter (the three-boson mass angle nu=pi/6), one can push the full analytical estimate of b_3 and a_3 up to an error that is in practice negligible. ----- Nous entendons par gaz de Bose unitaire un systeme compose de bosons sans spin interagissant dans l'onde s par un potentiel de longueur de diffusion infinie et de portee (reelle ou effective) presque negligeable, systeme pour l'instant abstrait mais dont la tentative de realisation avec des atomes froids est en cours. A partir de la solution analytique connue du probleme a trois corps dans un piege harmonique, et de methodes precedemment developpees pour des fermions, nous determinons le troisieme cumulant b_3, puis le troisieme coefficient du viriel a_3 de ce gaz, dans le cas spatialement homogene, en fonction de sa temperature et du parametre a trois corps R_t caracterisant l'effet Efimov. Un point marquant est qu'en convertissant des series en des integrales (par une methode des residus inverse), puis en utilisant un petit parametre inattendu, l'angle de masse nu=pi/6 des trois bosons, on peut pousser l'estimation completement analytique de b_3 et de a_3 jusqu'a une erreur en pratique negligeable.

연구 동기 및 목표

• 공간적으로 균일한 단위 보스 가스의 삼차 비열계수 $a_3$와 삼차 누적량 $b_3$를 결정함으로써, 실험적·이론적으로 아직 잘 규명되지 않은 상태를 해소하고자 한다.
• 에피모프 물리학을 기술하는 데 필수적인 삼체 매개변수 $R_t$를 단위 보스 가스의 비열 전개에 통합하고자 한다.
• 저온에서의 에피모프 트리머와 고온에서의 이온화 효과를 모두 고려한 $b_3$에 대한 고정밀도 해석적 방법을 개발하고자 한다.
• 정확한 수치 결과와 알려진 극한값과의 비교를 통해 해석적 접근의 정확성을 검증하고, 천분의 일 이내의 정확도를 확보하고자 한다.
• 앞으로의 초냉각 원자 실험, 특히 유한한 $R_t$를 가진 단위 근처에서의 강한 상호작용 보스 입자에 대한 이론적 틀을 제공하고자 한다.

제안 방법

• 조화 퍼텐셜 내 삼체 문제의 정확한 해법을 활용하여, 단위 보스 가스의 클러스터 적분 $b_3$를 유도한다.
• 이산적인 에너지 준위에 대한 합을 적분으로 전환하기 위해 역 잔여법을 적용함으로써, 해석적 계속 및 근사화를 가능하게 한다.
• 삼보스 질량 각도 $\nu = \pi/6$를 작은 매개변수로 간주하여, $b_3 = 3\sqrt{3}[S(0^+) + C]$에서 상수 $C$의 페르뮤티이션 전개를 수행한다.
• 고정밀도 수치 계산과 $\nu = \pi/6$에 대한 3차까지의 페르뮤티이션 전개를 가능하게 하는 $C$의 적분 표현을 유도한다.
• 표준 비열 전개 관계를 통해 $b_3$를 비열계수 $a_3$와 연결함으로써, 실험 측정치와 직접 비교할 수 있도록 한다.
• 정확한 $S(0^+)$와 $C$ 값과의 비교를 통해 근사치의 정확도를 검증하였으며, 오차가 0.1% 이내로 확인되었다.

실험 결과

연구 질문

• RQ1 삼체 매개변수 $R_t$와 온도의 함수로 나타낸 단위 보스 가스의 삼차 누적량 $b_3$의 해석적 형태는 무엇인가?
• RQ2 에피모프 물리학을 지배하는 삼체 매개변수 $R_t$는 어떻게 강한 상호작용 보스 가스의 비열 전개에 일관적으로 통합될 수 있는가?
• RQ3 에피모프 트리머와 그 온도 의존성 인구 수가 존재하는 상황에서도 삼차 비열계수 $a_3$를 고정밀도로 해석적으로 계산할 수 있는가?
• RQ4 저온 영역(에피모프 트리머가 支배하는 영역)과 고온 영역(트리머가 이온화되는 영역)에서 $b_3$의 거동는 어떠한가?
• RQ5 해석적 근사치 $b_3$는 저온에서의 지수형 형태와 고온에서의 로그형 형태와 같은 알려진 극한을 어느 정도 잘 재현할 수 있는가?

주요 결과

• 삼차 누적량은 $b_3 = 3\sqrt{3}[S(0^+) + C]$ 형태로 해석적으로 주어지며, 여기서 $S(0^+)$는 $\lambda/R_t$에 의존하고 $C \approx 0.648$는 유니버설 상수이다.
• 상수 $C$는 고정밀도 수치 계산과 $\nu = \pi/6$에 대한 페르뮤티이션 전개를 가능하게 하는 적분 표현으로부터 도출된다.
• 해석적 근사치 $b_3$는 정확한 값과 천분의 일 이내의 일치를 보이며, 이는 방법의 정확성을 검증한다.
• 저온 영역($k_B T \ll E_t$)에서는 $b_3 \simeq 3\sqrt{3} \, e^{\beta E_t}$로 나타나며, 기저 상태 트리머에서 기대되는 지수적 행동과 일치한다.
• 고온 영역($k_B T \gg E_t$)에서는 $b_3 \simeq 3\sqrt{3} \frac{|s_0|}{2\pi} \ln\left(e^{\gamma + 2\pi C/|s_0|} \beta E_t\right)$로 나타나며, 손실률 분석에서 유도된 로그 의존성과 일치한다.
• 이론은 $k_B T / E_t \approx 112.56$에서 $b_3$의 부호 변화를 예측하지만, 유한한 상호작용 범위를 가진 실제 모델에서는 관측되지 않을 수 있다.