Observation of universal non-Gaussian statistics of the order parameter across a continuous phase transition

기술 요약: 연속 상전이 과정에서 질서 매개변수의 보편적 비가우시안 통계 관측

문제 제기
2차 상전이는 근본적으로 임계 스케일링과 보편성에 의해 특징지어지며, 여기서 열역학적 특이점은 임계 지수(critical exponents)에 의해 결정된다. 질서 매개변수의 확률 밀도 함수(PDF)가 임계점 근처에서 비가우시안 통계를 나타낸다는 사실은 잘 알려져 있으나, 이는 보편적이며 유니버설리티 클래스(universality class)와 시스템 기하학에 의존하는 특징이다. 하지만 이러한 측면은 실험적으로는 거의 탐구되지 않았다. 거시적 시스템에서 비가우시안 영역은 제어 매개변수의 매우 좁은 구간이며, 시스템 크기가 커짐에 따라 더욱 축소되므로 표준적인 실험적 조절로는 접근이 불가능하다. 또한, 이러한 임계 비가우시안 통계의 특성을 규명하는 것은 본질적으로 어려우며, 기존의 실험 보고는 액정(liquid crystals)과 같은 고전적인 메조스코픽 시스템에 국한되어 왔다. 과제는 양자 다체 계(quantum many-body system)에서 연속 상전이를 가로지르는 질서 매변수의 전체 PDF를 측정함으로써, 표준적인 임계 지수를 넘어 임계 변동의 본질을 조사하는 것이다.

방법론
저자들은 3차원 광격자(optical lattice)에 로딩된 약 4,000개의 메타스테이블 헬륨($^4$He$^*$) 원자로 구성된 초저온 가스를 활용하여 보즈-허바드 모델(Bose-Hubbard model)을 구현하였다. 이 시스템은 열적 변동에 의해 구동되는 초유체-정상 유체(superfluid-to-normal fluid) 상전이를 겪으며, 3D 고전 XY 유니버설리티 클래스에 속한다. 실험은 조화 트랩(harmonic trap) 내에서 수행되어 불균일한 시스템을 생성하며, 이로 인해 서로 다른 상들이 공존하게 된다.

주요 실험 기술은 다음과 같다:

• 단일 원자 해상도 검출: 마이크로 채널 플레이트(MCP) 검출기를 사용하여, 자유 낙하 후 원자들의 전체 운동량 분포를 샷별(shot-by-shot)로 재구성한다. 이를 통해 중심 운동량 모드($k \approx 0$)에 있는 원자의 수 $N_0$로부터 응축 질서 매개변수 진폭 $|\psi| = \sqrt{N_0}$를 측정할 수 있다.
• 전체 계수 통계(Full Counting Statistics): 다양한 상호작용 강도($u = U/J$)에 대해 약 1,200회의 실험을 반복함으로써 완전한 PDF $P(|\psi|)$를 재구성한다.
• 사후 선택(Post-Selection): 총 원자 수 $N$의 샷 간 변동을 완화하기 위해, $N$에 대한 데이터 사후 선택을 수행하여 질서 매개변수의 통계만을 분리한다.
• 이론적 모델링: 측정된 PDF는 란다우 이론(Landau theory)과 유사한 유효 열역학 퍼텐셜 $F(|\psi|) = -\ln P(|\psi|)$를 사용하여 분석된다. 또한, 비가우시안성을 정량화하기 위해 $|\psi|$의 고차 첨은(high-order cumulants)을 계산한다. 수치적 비교는 균질한 3D XY 모델, 고전적 모델, 그리고 저온 양자 로터(low-temperature quantum rotor) 모델과 이루어진다.

주요 기여 및 결과

1. 비가우시안 PDF의 측정: 본 연구는 초유체 상전이를 가로지르는 질서 매개변수 진폭의 전체 확률 분포를 성공적으로 측정하였다. 정렬된(초유체) 상과 무질서(모트 절연체) 상에서는 분포가 가우시안 통계로 잘 근사된다. 그러나 임계점 근근처($u \approx u_c$)에서 PDF는 뚜렷하고 넓은 폭을 가진 비가우시안 형태를 보이며, 이는 질서 매개변수의 강화된 변동을 포착한다.
2. 유효 열역학 퍼텐셜: 저자들은 데이터로부터 유효 퍼텐셜 $F(|\psi|)$를 재구성하고, 이를 저차 다항식 형태 $F(|\psi|) \simeq a_0 + a_2|\psi|^2 + a_4|\psi|^4$로 피팅하였다. 이들은 퍼텐셜의 연속적인 진화를 관찰하였다: 정렬된 상에서는 0이 아닌 $|\psi|$에서 비자명한 극소값을 가지며, 무질서 상에서는 $|\psi|=0$에서 극소값을 가지며, 전이 근처에서는 평탄하고 넓은 분포를 보인다. 계수 $a_2$는 임계값 $u_{c,0}$에서 부호가 변하며 전이를 표시한다.
3. 평균 질서 매개변수의 임계 스케일링: 평균 질서 매개변수 $\langle |\psi| \rangle$는 멱법칙 스케일링 $\langle |\psi| \rangle \propto (u_c - u)^\beta$를 따르며, 이때 임계 지수는 $\beta = 0.35(1)$이다. 이 값은 3D XY 유니버설리티 클래스($\beta_{3DXY} \approx 0.3485$)와 일치하며 평균장 예측($\beta_{MF} = 0.5$)과는 구별되는데, 이는 유한한 크기와 트랩에도 불구하고 시스템이 보편적인 임계 스케일링을 포착했음을 확인시켜 준다.
4. 고차 첨은의 급격한 부호 변화: 시스템이 전이를 가로지를 때 질서 매개변수의 고차 첨은($n \geq 3$인 $\kappa_n$)에서 관찰되는 급격한 부호 변화는 핵심적인 발견이다. 이 첨은들은 정렬된 상과 무질서 상의 깊은 곳(가우시안 통계를 따르는 영역)에서는 소멸하지만, $u_c$ 근처에서는 유의미해진다.
5. 양자 대 고전적 행동: 첨은의 부호 변화는 실험과 균질한 양자 모델의 수치 시뮬레이션 모두에서 관찰되지만, $n \geq 4$ 차수의 첨은의 상세한 거동은 고전적 모델과 다르다. 실험 데이터는 고전적 모델에 의해 재현되지 않으나 저온 양자 모델에 의해 포착된다. 이는 편차가 고전적 모델이 아닌, PDF의 꼬리 부분에 영향을 미치는 양자 변동에 의한 보편적 스케일링 수정(universal corrections to scaling)에서 기인함을 시사한다.
6. 강건성: 첨은의 임계 거동은 두 가지 서로 다른 실험 데이터셋을 비교함으로써 시스템 기하학 및 트랩 구성의 변화에 대해서도 강건함이 입증되었다.

의의
본 논문은 임계 현상과 보편성을 조사하는 데 있어 질서 매개변수 통계의 결정적인 역할을 강조한다. 표준적인 임계 지수 측정을 넘어, 저자들은 전체 PDF와 그 고차 첨은이 임계 변동에 대한 보완적이고 매우 민감한 진단 도구가 될 수 있음을 보여준다. 메조스코픽 불균일 양자 시스템에서 비가우시안 통계와 첨은의 특정 부호 변화를 관측한 것은, 상관 길이가 시스템 크기보다 커질 때 발생하는 중심 극한 정리의 붕괴에 관한 재규격화 그룹(renormalization group)의 예측을 검증한다. 또한, 이러한 통계를 통해 고전적 및 양자적 임계 행동을 구별할 수 있는 능력은 열적 전이에서의 미묘한 양자 효과를 해결하는 데 있어 양자 시뮬레이터의 독특한 역량을 부각시킨다. 이 접근 방식은 임계 변동의 근본적인 본질에 대한 새로운 관점을 제공하며, QCD 상태도 연구를 포함한 복잡한 시스템의 상전이를 식별하기 위한 유망한 관측량을 제시한다.