Interaction-induced period shift of the Kohn oscillation in a parity-time-symmetric harmonic trap

技術要約：パリティ時間対称調和トラップにおけるコーン振動の相互作用誘起周期シフト

問題提起
標準的な調和トラップにおいて、ボース＝アインシュタイン凝縮（BEC）の重心（COM）はトラップ周波数で振動する。これはコーン振動として知られる現象である。相互作用が相対座標のみに依存する場合、調和ポテンシャルの定理により、この周波数は相互作用に対して保護されている。先行研究では、非相互作用系に（位置に比例する虚数ポテンシャルである）パリティ時間（PT）対称な利得・損失勾配を導入すると、この周波数が保持され、雲は剛体のように振動しながら原子数が周期的に変調されることが確立されている。

しかし、相互作用の存在はこの保護を打破する。利得および損失のメカニクスは総原子数（$N$）を変調し、これが相互作用エネルギーを変化させる。この変調は雲の幅（$S$）の変化を駆動し、それがさらに、幅の時間微分（$\dot{S}$）および位置・運動量共分散（$C$）を介して重心の動力学へとフィードバックを与える。本研究が取り組む中心的な問いは以下の通りである：この相互作用誘起のフィードバックは振動周期をシフトさせるのか？どの程度シフトさせるのか？そして、どの程度の利得強度までこれらの周期軌道は存続するのか？

手法
著者らは、虚数ポテンシャル項 $i\gamma x$ を含む一次元グロス＝ピタエフスキー（GP）方程式を用いてこの問題を調査している：
$$i\partial_t\psi = \left[ -\frac{1}{2}\partial_x^2 + \frac{1}{2}x^2 + i\gamma x + g|\psi|^2 \right]\psi$$
ここで、$\gamma$ は利得・損失の強さを、$g$ は相互作用の強さを表す。

1. 数値的継続法（Numerical Continuation）： 既知のコーン軌道（$\gamma=0$ の場合）を有限の $\gamma$ へと数値的に継続することで、周期軌道を算出している。この手法は、PT対称性を利用して、特定の位相でトラップ中心を横切る対称な軌道を定義する。また、線形安定性を評価するためにフロケ乗数を計算する。
2. モーメントモデル： 解析的な洞察を得るために、GP方程式から5変数モーメントモデルを導出する。このモデルは、原子数（$N$）、重心位置（$X$）、運動量（$P$）、二乗幅（$S$）、および共分散（$C$）を追跡する。モデルは、階層的なモーメント方程式を閉じるために、固定された密度プロファイル形状（スケーリング仮説）を仮定しており、これには量子圧力および相互作用項が含まれる。
3. 摂動解析： 小さな $\gamma$ に対して解析的な展開を行い、非エルミート・エーレンフェスト定理および双極モードと呼吸モードの結合を利用して、周期の主要項によるシフトを導出する。

主な結果

• 周期シフトのメカニズム： 周期 $T$ は、小さな $\gamma$ において利得強度 $\gamma$ に対して二次的に増加する。このシフトは、利得が $N$ を変調し、それが $S$ を駆動することによって生じる。$S$ が（呼吸モードの共鳴により）位相遅れを伴って $N$ に応答するため、重心の運動方程式における $2\gamma S$ の項が復元力に抗するように働く。
• 解析的予測： トーマス＝フェルミ（TF）領域において、相対的な周期シフトは次のように与えられる：
  $$\frac{T}{2\pi} - 1 \approx 0.79 g^{2/3} \gamma^2$$
  この係数は、雲のサイズ（基底状態の幅 $\sigma_0$ を通じて）のみに依存し、振動振幅には依存しない。
• モデルの精度： 5変数モーメントモデルは、雲のプロファイルが近似的に対称である限り、臨界利得に至るまで、雲の周期と幅を高い精度（数パーセント以内）で再現する。また、モデルは小さな $\gamma$ の極限において二次的なシフトを正確に捉えている。
• 臨界利得と安定性： 各周期軌道の族は、臨界利得 $\gamma_c$ で終了する。この時点で：
  • 軌道はトラップの高エネルギー状態（例：使用されたパラメータにおける $n \approx 475$ および $559$）と共鳴する。
  • 解の族は「折り返し（fold bifurcation）」を起こす。つまり、$\gamma_c$ の直下では、同じパラメータを持つ2つの軌道が存在する。
  • 動的な不安定性（フロケ乗数が単位円を離れる現象）は、多くの場合 $\gamma_c$ の直前またはその時点で発生するが、この不安定性は振動を即座に破壊することなく、雲の形状を歪ませる。
• スケーリング則： 臨界利達は $\gamma_c \propto (A \cdot R_{TF})^{-1/2}$ とスケーリングする（ここで $A$ は振動振幅、$R_{TF}$ はトーマス＝フェルミ半径）。具体的には、計算されたすべての族において $\gamma_c \sqrt{A R_{TF}} \approx 0.280$ である。

意義と主張
本論文は、PT対称トラップにおける「コーン周波数の保護」が、相互作用の存在下では脆弱であることを示すと主張している。利得と損失の相互作用による原子数変調と雲の幅の結合が、剛体的な振動を打破し、測定可能な周期シフトを引き起こす。

著者らは以下の点を強調している：

1. シフトは、利得、損失、および相互作用の相互作用によって駆動される $\gamma$ に関する二次効果である。
2. 5変数モーメントモデルは、臨界点まで、パラメータに依存しないロバストな記述を提供し、周期シフトの物理的メカニズムを検証している。
3. 軌道が高エネルギーのトラップ状態と共鳴する臨界利得が存在することが、これらの周期解の観測可能な範囲を制限している。
4. これらの結果は、利得がコヒーレントであり、かつ雲の大きな変位に対してもトラップが調和的であり続けるという条件下において、光導波路（光の流体）やエンジニアリングされた損失を持つ原子凝縮体の実験的実現に関連している。

本研究は、臨界利得を超えた完全な非線形動力学を解いたり、カオス的挙動を予測したりすることを目的としたものではなく、コーン極限から現れる周期軌道の特性に焦点を当てている。