Dynamical thermalization, Rayleigh-Jeans condensate, vortexes and wave collapse in quantum chaos fibers and fluid of light

技術要約：量子カオス・ファイバーおよび光の流体における動的熱化、レイリー・ジーンズ凝縮、渦、および波の崩壊

問題提起
本論文は、カオス的なD型ビリオド（ビリヤード）幾何学内において、非線形シュレディンガー方程式（NSE）によって支配される非線形場の時間発展を調査している。このモデルは、カオス的な断面を持つマルチモード光ファイバーにおける縦方向の光伝搬の記述と、カー非線形媒体（原子蒸気）における「光の流体」のダイナミクスのモデリングという、二重の物理的目的を担っている。中心となる問題は、外部の熱浴が存在しないハミルトン系における動的熱化のメカニズムを扱うことである。具体的には、弱非線形性におけるコルモゴロフ・アーノルド・モーザー（KAM）理論に支配される準可積分ダイナミクスのレジームと、特定の「カオス境界」を超えた際に発生するカオス的熱化のレジームを区別することを目的としている。主要な目標は、結果として得られる熱分布が量子ボース＝アインシュタイン（BE）統計に従うのか、あるいは古典的なレイリー・ジーンズ（RJ）統計に従うのかを判断し、これらの条件下での基底状態への凝縮形成を分析することである。

手法
本研究では、解析的な導出と広範な数値シミュレーションの組み合わせを用いている。

• モデル: システムは、D型ビリオド（円から直線部分を切り取った形状）におけるディリクレ境界条件を持つNSEで定義される。非線形パラメータ $\beta$ は相互作用の強さを制御し、$\beta > 0$ は斥力性（defocusing）、$\beta < 0$ は引力性（focusing）の非線形性を表す。
• 数値積分: 2つの手法が用いられている。短時間ダイナミクス（$t \le 40$）に対しては、ノルムとエネルギーを高精度に保存するトロッター分解法が使用される。長時間発展（$t$ が最大1000まで）については、長期にわたる安定性と精度を確保するため、最大 $1.36 \times 10^5$ 個のメッシュノードを用いた有限要素法によるFREEFEMパッケージを採用している。
• 初期条件: シミュレーションは通常、パリティ対称性を破り、カオス的な混合を誘発するために、線形固有モードの重ね合わせ（例：$(\psi_m + \psi_{m+1})/\sqrt{2}$）から開始される。
• 解析的枠組み: 著者らは、マイクロカノニカルアンサンブル統計と鞍点近似を用いて、等間隔なエネルギー準位を持つ簡略化されたモデル（RJSモデル）におけるレイリー・ジーンズ凝縮率の解析的な推定値を導出し、モンテカルロシミュレーションと比較している。

主な貢献と結果

1. 動的熱化とカオス境界:
   著者らは、非線形強度が特定のカオス境界を超えると、システムが動的熱化を起こすことを示している。この境界を下回る場合、ダイナミクスは準可積分（KAMレジーム）に留まる。境界を上回ると、システムは以下の古典的なレイリー・ジーンズ熱分布によって記述される定常状態へと進化する：
   $$\rho_m = \frac{T}{E_m - \mu}$$
   ここで、$T$ は温度、$\mu$ は化学ポテンシャルである。これは、熱化が外部の散逸ではなく、純粋にハミルトン的なカオスによる結果であることを裏付けている。

2. レイリー・ジーンズ（RJ）凝縮:
   主要な知見は、全確率（または粒子数）の約80〜90%が基底状態（$m=1$）に蓄積する、RJ凝縮の出現である。

• メカニズム: この凝縮は、モード数 $N$ が大きいシステムにおいて、低い実効温度（$T \ll E_m$）で発生する。非凝縮確率の割合 $W_{ex}$ は $(\ln N)/N$ に比例してスケールし、$N \to \infty$ のとき消失する。
• フロリッヒ凝縮との区別: 本論文は、このRJ凝縮をフロリッヒ凝縮とは明確に区別している。フロリッヒモデルは高温（$T \gg E_m$）と外部からのポンピングおよび散逸を必要とする。対照的に、RJ凝縮は、外部からのポンピングなしに、保存的なハミルトン系において低温度で発生する。
• ボース＝アインシュタイン（BE）凝縮との区別: 数値結果は、BE分布とRJ分布の劇的な違いを示している。BE統計であれば、検討されたパラメータにおいて基底状態の占有率は $\approx 0.38$ と予測されるが、RJ統計は $\approx 0.98$ を予測する。シミュレーションはRJの予測と一致しており、古典的な場の熱化の性質を裏付けている。

3. エントロピーの進化:
   研究では、古典的なボルツマン・エントロピー（$S_B$）と量子フォン・ノイマン・エントロピー（$S_q$）の両方を追跡している。

• $S_B(t)$ は、ボルツマンのH定理と一致して単調に増加する。
• $S_q(t)$ は非単調な挙動を示す。まずエネルギーが高次モードへ拡散するにつれて増加し、極大に達した後、システムがRJ凝縮状態に落ち着くにつれて大幅に減少する。この減少は、確率が基底状態に高度に集中することで、実効的な占有モード数が減少することに起因する。

4. 波の崩壊（引力性非線形性）:
   引力性のケース（$\beta < 0$）において、著者らは波の崩壊を調査している。無限空間における崩壊が有限の時間で起こるのに対し、D型ビリオドの有限サイズとカオスの存在は、崩壊と熱化の間の相互作用をもたらす。本研究は、開いた空間のシナリオとは異なる挙動として、十分に高い正のエネルギーにおいても波の崩壊が起こり得ることを示している。

5. 渦のダイナミクス（斥力性非線形性）:
   強い非線形性を伴う斥力性のレジーム（$\beta > 0$）において、システムは長寿命の渦を特徴とする超流動レジームを示す。これにより、D型ファイバーが、カオス的な「光の流体」における渦のダイナミクスを研究するためのプラットフォームとなることが確立された。

6. 実験的関連性:
   論文では、マルチモード光ファイバーや原子蒸気システムにおいてこれらの効果を実現するためのパラメータについて論じている。従来の光ファイバーにおける「セルフクリーニング」に関する実験は、KAM理論が成立しない縮退したモードスペクトルで行われていた可能性があるが、D型ファイバーはKAMレジームが有効な一般的なケースを提供しており、カオス境界とそれに続くRJ熱化の観測を可能にしている。

意義
本論文は、一般的な非線形振動子系における動的熱化に関する基礎的な理解を提供すると主張している。D型ビリオドを用いることで、スペクトルの縮退からカオスの影響を分離し、熱化が量子ボース＝アインシュタイン分布ではなく、古典的なレイリー・ジーンズ分布をもたらすことを実証している。RJ凝縮を、フロリッヒ凝縮とは異なる低温度の保存現象として特定したことは、古典的な場の統計力学を明確にするものである。さらに、本研究は量子カオス理論、非線形光学、および流体力学の間の架け橋となり、マルチモードファイバーや原子蒸気における熱化、凝縮、および渦ダイナミクスの実験的研究のための理論的基礎を提供するものである。