Expansion of a free Fermi gas released from an isotropic trapping potential

기술 요약: 등방성 트래핑 포텐셜에서 방출된 자유 페르미 가스의 팽창

문제 정의
본 연구는 초기 상태가 등방성 트래핑 포텐셜 $V(\mathbf{x})$의 다체 바닥 상태(many-body ground state)로 준비된 $N$개의 상호작용하지 않는 스핀 없는(또는 스핀 퇴화된) 페르미온 시스템의 비평형 역학을 조사한다. 연구의 초점은 $t=0$에서 포텐셜이 즉각적으로 제거된 이후의 상황에 맞춰져 있다. 저자들은 특히 거대 페르미온 수 $N$의 극한에서 가스의 탄도적 팽창(ballistic expansion)을 분석하여, 후기 시간(late times)에서의 완전한 공간 및 운동량 프로파일을 도출하는 것을 목표로 한다. 주요 목적은 초기 트래핑 포텐셜의 형태(지수 $a$로 특징지어짐)와 공간 차원 $d$가 팽창하는 구름의 점근적 밀도 프로파일과 가장자리 거동에 어떻게 영향을 미치는지 결정하는 것이다.

방법론
저자들은 특정 사례에 대한 정확한 해석적 해법과 일반적인 경우의 대규모 $N$에 대한 준고전적(WKB) 근사를 결性하여 결합하여 사용한다.

1. 정확한 해법 (조화 및 사각 우물 포텐셜):

   • 조화 포텐셜($V \propto x^2$)의 경우, 저자들은 Ermakov-Pinney 방정식에 의해 제어되는 다체 파동함수의 정확한 시간 의존성을 활용한다. 이들은 페르미온들의 재조정된 위치(rescaled positions)가 가우시안 유니터리 앙상블(GUE) 랜덤 행렬의 고유값 분포를 따름을 입증한다.
   • 1차원 사각 우물 포텐셜의 경우, 유한한 $N$에 대한 운동량 분포 $f_N(p)$를 계산하기 위해 정확한 파동함수(2D에서는 베셀 함수, 1D에서는 삼각 함수 포함)를 사용하여 가장자리 거동 및 유한 크기 보정을 분석한다.

2. 준고전적 근사 (일반 포텐셜):

   • 임의의 등방성 멱법칙(power-law) 포텐셜 $V(r) \propto r^a$에 대하여 임의의 차원 $d$에서, 저자들은 전개된 파동함수에 대해 정지 위상 근사(stationary-phase approximation)를 적용한다. 이를 통해 후기 시간에서 페르미온의 위치 $\mathbf{x}$가 초기 운동량 $\mathbf{p}$와 $\mathbf{x} \approx \mathbf{p}t$라는 탄도적 스케일링 관계를 가짐을 확립한다.
   • 운동량 분포 $f_N(p)$는 점유된 단일 입자 고유상태들의 푸리에 변환의 절댓값 제곱의 합으로 유도된다. 거대 $N$ 극한에서, 이 합은 상태 밀도에 대한 토마스-페르미 근사와 준고전적 양자화 규칙을 사용하여 적분으로 변환된다.
   • 고정된 에너지 $E$에서의 미시 정준(microcanonical) 운동량 분포는 고전적 위상 공간 부피를 기반으로 계산되며, 이는 페르미 에너지 $E_F$까지 점유된 에너지 상태들에 대해 적분된다.

주요 기여 및 결과

1. 일반화된 밀도 프로파일:
   본 논문은 임의의 차원 $d$와 임의의 멱법칙 포텐셜 $V(r) \propto r^a$에 대한 명시적인 후기 시간 공간 밀도 프로파일 $\rho_N(r; t)$와 운동량 분포 $f_N(p)$를 유도한다.

   • 운동량 분포: $p < p_F$인 경우, 분포는 다음과 같이 스케일링된다:
     $$f_N(p) \sim \frac{N}{p_F^d} \left( 1 - \frac{p^2}{p_F^2} \right)^{d/a}$$
   • 공간 프로파일: 결과적으로, 구름의 후기 시간 공간 프로파일(for $r < p_F t$)은 다음과 같다:
     $$\rho_N(r; t) \sim \frac{N}{(p_F t)^d} \left( 1 - \frac{r^2}{(p_F t)^2} \right)^{d/a}$$
     여기서 $p_F$는 페르미 운동량이며, $p_F \sim N^{a/((a+2)d)}$로 $N$에 따라 스케일링된다.

2. 보편적 가장자리 지수 (Universal Edge Exponent):
   중심적인 발견은 보편적인 가장자리 지수 $d/a$의 식별이다. 밀도 프로파일은 Wigner semi-circle law(이는 $d=1, a=2$인 특수 사례에 해당함)나 Thomas-Fermi 분포($d=3, a=2$)의 특징인 제곱근 특이점이 아니라, 이동하는 가장자리($r = p_F t$)에서 멱법칙 특이점을 보인다. 이 결과는 모든 멱법칙 포텐셜과 차원에 대해 알려진 법칙들을 일반화한다.

3. 가장자리 거동 및 스케일링 함수:

   • 조화 트랩 ($a=2$): Wigner semi-circle law의 가장자리가 상대적 폭이 $N^{-2/3}$으로 스케일링되는 매끄러운 크로스오버로 확장됨을 보여준다. 가장자리 근처의 프로파일은 Airy 함수의 미분 제곱을 포함하는 스케일링 함수 $A(\xi) = \text{Ai}'(\xi)^2 - \xi \text{Ai}(\xi)^2$에 의해 기술된다. Wigner 구간 외부의 평균 페르미온 수는 작은 유한 상수($\approx 0.061$)이다.
   • 사각 우물 포텐셜 ($a \to \infty$): 무한 우물의 극한에서, 운동량 분포는 페르미 볼 내부에서 균일해진다. 그러나 가장자리에서의 감소는 디감마(digamma) 및 트리감마(trigamma) 함수를 포함하는 다른 스케일링 함수 $F(\xi)$에 의해 기술된다. 조화로운 경우와 달리, 날카로운 포텐셜은 분포에서 긴 대수적 꼬리(algebraic tail) $F(\xi) \sim 1/\xi$를 유발하며, 이로 인해 페르미 구간 외부의 페르미온 수는 $N$에 따라 로그 함수적으로 증가한다 ($N_{out} \sim \ln N$).

4. 고차원 효과:
   $d \ge 2$ 차원에서 궤도 퇴화(orbital degeneracy)의 존재가 다뤄진다. 저자들은 퇴화된 상태의 완전한 껍질(complete shells)에 해당하는 특정 페르미온 수에 대해, 시스템이 유일한 등방성 바닥 상태를 가짐을 보여준다. 2차원 무한 원형 우물의 경우, 적분된 상태 밀도의 진동 성분을 분석하여, 페르미 운동량이 완전한 껍질 수에 제한되더라도 불규칙한 진동을 보임을 입증하며, 이는 고전적 빌리어드(billiard)의 주기 궤도와 연관되어 있다.

의의 및 주장
저자들은 이러한 결과들이 Wigner semi-circle law와 Thomas-Fermi 분포를 일반화하는 포괄적인 분석적 틀을 제공함으로써 독창적이라고 주장한다. 본 연구는 자유 페르미온의 팽창에 대한 상세한 분석적 접근이 상호작용하는 시스템에서는 종종 불가능하다는 점을 강조하며, 초기 구속의 기하학적 구조(지수 $a$를 통해)와 팽창하는 구름의 보편성 클래스 사이의 명확한 연결을 확립한다. 사각 우물 포텐셜에서 페르미 가장자리 외부로의 로그적 입자 증가와 조화 트랩에서의 상수적 입자 수의 차이는 초기 포텐셜의 날카로움이 가장자리 요동의 본질을 어떻게 결정하는지를 보여준다. 이 작업은 귀환 확률(return probabilities)에 대한 이전 연구들의 후속 연구로서, 분석을 전체 밀도 프로파일과 고차원으로 확장한다.