Rotating two-dimensional Bose--Einstein condensates in power-law traps: collapse, Thomas--Fermi scaling, and the inverse-square endpoint

技術要約：べき乗型トラップにおける回転二次元ボース–アインシュタイン凝縮体

問題設定
本研究は、べき乗型トラップポテンシャル $V_s(r) = r^s/s$（指数 $s$ は $-2$から$4$まで）の下での、グロス–ピタエフスキー（GP）方程式によって記述される二次元ボース–アインシュタイン凝縮体（BEC）の基底状態の性質と安定性を調査するものである。本研究では、運動エネルギー、接触相互作用（引力または斥力）、および外部回転の相互作用を取り扱う。中心的な焦点は、引力的相互作用における崩壊閾値付近での凝縮体のサイズのスケーリング挙動、および斥力的相互作用におけるトーマス–フェルミ（TF）極限、特に回転および特異なポテンシャル（具体的には逆二乗ポテンシャル$s = -2$）が存在する場合における挙動である。

手法
著者らは、解析的なスケーリング議論と厳密な数値シミュレーションを組み合わせて用いている：

1. スケーリング解析： スケーリング議論を用いて、凝縮体の幅 $\ell$ およびエネルギーに関するスケーリング則を導出する。ノルム保存型ダイレーション（膨張）の下で、運動エネルギーと接触エネルギーは同一（$\lambda^2$）にスケールし、トラップエネルギーは $\lambda^{-s}$ にスケールすることを確立している。
2. 数値ツール：
   • 2D虚時間伝搬： 引力的および斥力的領域の両方において、回転の有無にかかわらず定常状態を計算するために、Strang分割を用いた正規化勾配流を用いる。
   • 径方向ソルバー： 非回転状態のベンチマークとして、高精度な径方向解（有限体積法、境界付きニュートン反復、リチャードソン補外）を用いる。
   • ボゴリューボフ・スペクトル： 中心状態の安定性閾値（特に双極子周波数）を決定するための励起周波数の計算。
   • 変分境界： ガウス・ブーストされたTownes試行関数によるエネルギーの最小化を行い、基底状態エネルギーの上界を求める。
3. 正則化： 特異なトラップ（$s \le 0$）については、数値的な発散を扱うためにカットオフ $\epsilon$（例：$V_{s,\epsilon} = (r+\epsilon)^s/s$）を用いて正則化を行う。特に逆二乗ポテンシャルの場合における $\epsilon \to 0$ の極限に注意を払う。

主要な貢献と結果

1. スケーリング則と崩壊ダイナミクス

• 引力： 凝縮体は、トラップの形状や回転速度に関わらず（トラップが閉じ込め条件を満たす限り）、Townes結合定数 $g^*$ で崩壊する。崩壊への接近は、$\ell \propto \delta^{1/(s+2)}$（ここで $\delta = 1 - |g|/g^*$）という法則に従う。この法則は、7つの異なる指数（$s \in \{1, 2, 3, 4, -0.5, -1, -1.5\}$）に対して成立する。
• 正則化されたコア： 有限のカットオフ $\epsilon$ を持つ特異トラップの場合、凝縮体の幅が $\epsilon$ と同程度になると、スケーリング則は「コーン則」 $\ell \propto \delta^{1/3}$ へと遷移する。これは、遠方のべき乗則ではなく、正則化されたコアの局所的な線形形状によって決定される。
• 斥力： 斥力的ガスの場合、雲はトーマス–フェルミ半径 $R_{TF} \propto g^{1/(s+2)}$ まで膨張する。TFエネルギー誤差は、密度勾配成分とフロー成分（渦のフローと剛体回転の不一致）に正確に分離されることが示されており、特異なトラップにおいては、フロー成分が勾配成分と同等かそれ以上になる。

2. 回転と対称性の破れ

• 双極子基準： $s > 2$ のトラップにおける引力的な雲について、回転周波数 $\Omega$ が雲の双極子周波数 $\omega_1$ を超えると、中心状態は不安定になる。
  • $s=2$（調和振動子）の場合、$\omega_1$ は一定（Kohnモード）であり、遠心力限界以下では回転は対称性の破れを誘発しない。
  • $s > 2$ の場合、トラップの中心付近は「ソフト」であり、$\omega_1$ は雲が集中するにつれて減少する（$\omega_1 \propto \ell^{(s-2)/2}$）。その結果、弱い引力であっても回転が中心状態を不安定化させることがある。
• 四重極トラップ（$s=4$）： 対称性の破れの閾値は $\delta_c \simeq 0.324 \Omega^6$ であることが判明した。この結果はTownes試行関数から導出され、Bogoliubov計算によって3%以内の精度で確認されている。
• 不安定化後の状態： 一度不安定になると、引力的な雲は軸からずれ、回転系ポテンシャル $V_{eff} = V - \frac{1}{2}\Omega^2 r^2$ の最小環の上に崩壊する。この環に局在する塊の幅は、有効ポテンシャルではなく、元のトラップの曲率によって決定される。
• 斥力的回転： 斥力的なガスにおいて、回転は渦を核生成する。四重極トラップでは、臨界的な減少回転 $\zeta_c = (3/\pi)^{1/6}$ において、TF密度が中心部で空になり、アニュラス（環状構造）を形成する。これは相互作用強度 $g$ に依存しない。

3. 逆二乗の終端（$s = -2$）

• 共変性： $s = -2$ において、正則化されたエネルギーはカットオフ $\epsilon$ と正確に共変である（$E \propto \epsilon^{-2}$）。固有の長さスケールは存在せず、システムは引力的・斥力的な相互作用の両方において中心へ落下する。
• ヘイロー（暈）構造： 斥力的相互作用の場合、凝縮体は高密度のコアと、それを囲む普遍的なヘイローを形成する（密度プロファイル $|\psi|^2 \propto 1/(gr^2)$）。このヘイローの広がり $L_h$ は、結合定数に対して指数関数的に成長する：$L_h \propto \exp(g / [\pi(1+\kappa)])$、ここで $\kappa$ はテールの強度を特徴付ける。
• 普遍性： ヘイローのプロファイルとその指数関数的成長率は普遍的（特定の正則化形式に依存しない）であるが、コアの性質は正則化の詳細に依存する。

4. 崩壊付近の数値的アーティファクト

• 本論文は、崩壊閾値付近における分割ステップ虚時間伝搬の一次の定常バイアスを特定している。このバイアスにより、計算された幅は真の基底状態よりもわずかに大きくなる。
• 決定的なことに、この誤差は標準的な残差テスト（GP方程式の残差を測定するもの）では検出できない。なぜなら、誤差はブリージングモード（呼吸モード）と共にスケールするためである。このバイアスは、閾値付近で $1/\delta$ に比例して増大する。この誤差を排除するには、対称ステップスキームまたはゼロ時間ステップへの補外が必要である。

意義と主張
著者らは、主に以下の3つの新しい結果を主張している：

1. 回転誘起崩壊： 回転が（ソフトなトラップ $s>2$ において）引力的な凝縮体をどのように軸外へ移動させ、環状崩壊へと導くかという正確なメカニズムと閾値（$\delta \propto \Omega^6$）。
2. 逆二乗の終端： カットオフに対するエネルギーの正確な共変性と、斥力的ガスの $s = -2$ 極限における普遍的な指数関数的ヘイローの導出。
3. 数値的バイアス： 虚時間法における、標準的な収束判定では検出できず、かつ幅に大きな影響を与える一次の定常バイアスの特定と特性評価。

本論文は、Townes閾値、一般的なスケーリング則 $\ell \propto \delta^{1/(s+2)}$、あるいは引力的な雲における角運動量の中心質量起源といった既知の結果についての新規性を主張しているのではない。むしろ、これらの既知の法則をべき乗型の全範囲にわたって定量的に検証し、それらを統一的な枠組みへと統合している。本研究は、トラップの幾何学、相互作用の符号、および回転が、特に特異なポテンシャルや臨界崩壊付近において、いかにして2D BECの安定性と構造を共同で決定するかを包括的に記述している。