Universal self-similar evolution of two-dimensional Bose-Einstein condensates in the acoustic regime

問題の所在
本論文は、音響領域において非平衡状態に置かれた二次元（2D）ボース・アインシュタイン凝縮体（BEC）の、非定常かつ自己相似的な進化について調査している。具体的には、初期の広範なスケールを持つ凝縮体が、どのようにしてエネルギーをより小さなスケール（紫外線、UV領域）へと転送する乱流カスケードへと進化するかという問題に取り組んでいる。自己相似的な進化は物理学における一般的な現象（拡散や核爆発など）であるが、非線形波動系におけるスケーリングの性質は、しばしば二つのクラスに分類される。一つは次元解析と保存則によって決定される「第一種」、もう一つは非線形固有値問題によってスケーリング指数が決定される「第二種」である。著者らは、2D BEC、特に音響波が存在する場合の自己相似的なダイナミクスを特徴付け、そのシステムが第一種または第二種の自己相似性を辿るのかどうかを明らかにすることを目的としている。解決すべき特定の課題として、原子BECと比較して、有限の寿命と強い散逸に悩まされるポラリトンBECの挙動が挙げられている。

手法
本研究では、グロス・ピタエフスキー方程式（GPE）と波力学方程式（WKE）の数値シミュレーションを組み合わせた二角的なアプローチを採用している。

1. グロス・ピタエフスキー方程式 (GPE): 著者らは、原子およびポラリトンBECを記述する2D GPEをシミュレートしている。この方程式には、音速 ($c_s$)、補正長 ($\xi$)、外部トラップポテンシャル ($V_{trap}$)、およびポラリトンの損失をモデル化するための線形減衰項 ($\gamma$) が含まれる。

   • シミュレーションには、擬スペクトル法を用いたコード FROST を使用した。
   • 初期条件は、広範なスケールにおける弱いランダムなボゴリューボフ波によって摂動を与えられた基底状態とした。
   • セットアップとして、ユニバーサルな特徴をより良く観察するために、サイズ $75\xi$ の均一な正方形トラップと、トラップのないより大きな周期境界ボックス（サイズ $1024\xi$）の二種類を用いた。

2. 波力学方程式 (WKE): 理論的な洞察を得るために、著者らは弱波乱流理論（WWTT）の枠組みを利用した。彼らは、2Dにおける広範なスケールのボゴリューフ波に対するWKEを導出し、数値的に解いた。

   • WKEは、弱非線形極限におけるエネルギー・スペクトル ($e_k$) の時間進化を記述する。
   • シミュレーションは、4デケードにわたるスペクトル範囲でコード WavKinS を用いて行われた。
   • WKEには、厳密な導出では通常、強制力や散逸がないことを仮定していることを認めつつ、散逸をモデル化するために減衰項がアドホックに追加されている。

3. 解析: 著者らはエネルギー・スペクトルの進化を追跡し、乱流フロントの伝播を定量化するために、特性波数を $k^*(t)$ （エネルギー重み付き平均）として定義した。彼らは、$e_k(t) = E_0 f(t) \Phi(k/k^*(t))$ という形式の自己相似解を求めた。

主な貢献と結果

• 自己相似領域の特定: GPEとWKEの両方のシミュレーションにおいて、フロントが高波数の方向へ伝播し、背後に冪乗則スペクトル $e_k \sim k^{-1}$ を残していく自己相似的な進化が明らかになった。この指数は、2D音響波における直接エネルギーカスケードのコルモゴロフ・ザカロフ（KZ）解に対応している。
• 第二種自己相似性: 「無限容量」を持つシステム（通常、次元解析によって決定される第一種自己相似性を示す）への予想に反して、2D音響系は第二種自己相似性を示す。スケーリング指数は次元的な議論のみから導き出すことはできず、代わりに非線形固有値問題から出現する。
• 普遍定数 $\beta$: 伝播するフロントのダイナミクスは、単一の無次元の普遍定数 $\beta$ によって支配されている。著者らは、フロントの位置 $k^*(t)$ に関する方程式を導出した：
  $$\frac{1}{k^*} \frac{dk^*}{dt} \ln\left(\frac{k^*}{k_0}\right) = \frac{\beta}{\tau} e^{-t/\tau_D}$$
  ここで、$\tau$ は特性時間スケールであり、$\tau_D$ は散逸時間である。
• $\beta$ の数値的決定: $k^*(t)$ の数値的進化を、GPEおよびWKEシミュレーションから理論的予測にフィッティングすることにより、著者らは以下の普遍的な値を決定した：
  $$\beta \approx 8.9$$
• 散逸からの独立性: 驚くべき結果は、自己相似形式 $\Phi(\eta)$ と定数 $\beta$ が、減衰係数 $\gamma$ に依存しないことである。散逸（ポラリトンに関連する）が存在するにもかかわらず、システムは同じ普遍的な自己相似形式を保持している。これは、散逸項の斉次性が、WKEにおける衝突積分と一致していることに起因する。
• スペクトル構造: 自己相似関数 $\Phi(\eta)$ は、小さな $\eta$ において冪乗挙動 $\eta^{-1}$（KZ解と一致）を示し、大きな $\eta$（小さなスケール）において指数関数的減衰 $\exp(-\beta \eta / D)$ を示す。ここで $D$ は拡散係数である。

意義と主張
本論文は、2D音響波のダイナミクスを、退化したタイプの第二種自己相似解の具体的な例として確立したと主張している。その意義は以下の点にある：

1. 普遍的な挙動: システムは、初期条件や散逸の存在に依存しない普遍的な定数 $\beta \approx 8.9$ とペア $(\Phi, \beta)$ によって完全に特徴付けられる。
2. 境界的なシステム: 著者らは、2D音響波が「境界的な無限容量システム」を代表していると示唆している。エネルギー積分は、標準的な無限容量システムにおけるUV発散とは異なり、KZスケーリングにおいては対数的にしか発散しないが、このシステムは第一種と第二種の両方の自己相似性領域の性質を継承している。
3. 理論的検証: 本研究は、弱波乱流理論（WWTT）が、減衰がある場合でもGPEの自己相似的な進化を予測できることを示し、微視的なGPEダイナミクスと巨視的なWKE記述との間の溝を埋めるものである。
4. 実験的関連性: これらの知見は、原子BEC（非粘性の場合）とポラリトンBEC（散逸がある場合）の両方に重要であり、これらにおける非平衡実験、特に小スケール構造の形成と乱流フロントの伝播を解釈するための理論的枠組みを提供する。

著者らは、このシステムにおけるスケーリング指数は、次元的な制約によって課されるのではなく、動的に選択されるものであると結論付けており、これは2D量子流体における非平衡進化の複雑さを浮き彫りにしている。