[論文レビュー] On soliton (in-)stability in multi-dimensional cubic-quintic nonlinear Schrödinger equations

Rémi Carles, Christian Klein+1 more|arXiv (Cornell University)|Dec 21, 2020

Advanced Mathematical Physics Problems参考文献 33被引用数 5

ひとこと要約

本稿は、立方項が集束的で、五次項が拡散的である2次元および3次元の立方-五次非線形シュレーディンガー方程式における基底状態の軌道的安定性を調査する。時間分割スペクトル法を用いた数値シミュレーションにより、作用最小化基底状態は、その質量が臨界閾値を上回る場合に軌道的安定であることを示している。一方、臨界未満の質量を有する摂動は完全に分散的解へと進化する。主な発見は、過剰な質量を有する不安定な基底状態が、放射を介して安定で、周波数が低い孤立波へと進化することであり、これは異なる非線形モデル間で普遍的な不安定性シナリオを確認するものである。

ABSTRACT

We consider the nonlinear Schrödinger equation with a focusing cubic term and a defocusing quintic nonlinearity in dimensions two and three. The core of this article is the notion of stability of solitary waves. We recall the two standard notions of orbital stability in the context of nonlinear Schrödinger equations, and show that they must be considered as independent from each other. We investigate numerically the notion of orbital stability of ground states in the radially symmetric case, confirming existing conjectures or leading to new ones.

研究の動機と目的

• 多次元立方-五次NLSにおける基底状態の安定性分類に関する長年の曖昧さを解消し、エネルギー最小化解と作用最小化解を区別すること。
• 特に、臨界閾値に対する質量の役割に注目して、回転対称設定における最小作用基底状態の軌道的安定性を数値的に調査すること。
• 特に、摂動を受ける孤立波の長期的ダイナミクスに関する既存の予想を検証または精緻化すること、特に安定基底状態への漸近的収束について。

提案手法

• 本研究では、次元 d=2 および d=3 における立方-五次非線形シュレーディンガー方程式を数値的に解くために、時間分割スペクトル法を用いる。
• 基底状態は、周波数 ω ∈ (0, 3/16) のさまざまな値に対して、定常方程式 $-\frac{1}{2}\Delta\phi + \omega\phi - |\phi|^2\phi + |\phi|^4\phi = 0$ の解として計算される。
• 軌道的安定性は、既知の作用最小化基底状態 $Q_\omega$ を、小さな局所的 $H^1$-正則摂動で摂動し、時間経過に伴う $L^\infty$-ノルムの変化を追跡することで評価される。
• 解の最終状態は、事前に計算された基底状態のライブラリ $Q_{\underline{\omega}}$ と比較され、安定孤立波の漸近的周波数 $\underline{\omega}$ が同定される。
• 初期データの質量は、不安定な $Q_\omega$ と最終的な安定状態 $Q_{\underline{\omega}}$ の質量と比較され、放射損失が定量的に評価される。

実験結果

リサーチクエスチョン

• RQ1 多次元立方-五次NLSにおける作用最小化基底状態は、どのような条件下で軌道的安定であるか？
• RQ2 初期データの質量が、摂動を受けた孤立波が局在的であるか、分散的であるかを決定づける要因となるか？
• RQ3 過剰な質量を有する不安定な基底状態は、周波数が低く、質量が小さい安定基底状態に進化できるか？
• RQ4 放射は、この系における摂動を受ける孤立波の長期的ダイナミクスにおいて、どのような役割を果たすか？
• RQ5 漸近的安定孤立波の最終周波数 $\underline{\omega}}$ に対する普遍的な選択機構は存在するか？

主な発見

• 初期データの質量が不安定基底状態 $Q_\omega$ の質量を上回る場合、解は $\underline{\omega} < \omega$ を満たす安定孤立波 $e^{i\underline{\omega}t}Q_{\underline{\omega}}$ へと進化し、放射による質量損失が示唆される。
• 初期質量が $M(Q_\omega)$ より小さい場合、解は完全に分散し、純粋な分散的挙動が確認される。
• 数値的シミュレーションでは、解の最終 $L^\infty$-ノルムと安定基底状態 $Q_{\underline{\omega}}$ の $L^\infty$-ノルムとの間に良好な一致が得られ、ω = 0.01 の場合 $\underline{\omega} \approx 0.047$ となる。
• 最終的な安定状態の質量は常に初期質量より低く、ω = 0.01 の場合 $M(Q_{\underline{\omega}}) \approx 77.05$ であるのに対し、$M(Q_{\omega=0.01}) \approx 79.44$ であるため、顕著な放射が確認される。
• 不安定性シナリオは、対称的および非対称なガウス型摂動の両方に対して、$\omega < \omega_{\text{crit}}$ のさまざまな値で成立する。
• 最終的な漸近的状態は、常に初期の不安定状態よりも周波数が低く、質量が小さい安定基底状態であるため、予想5.1（普遍的な動的緩和）を支持する。