[論文レビュー] Rayleigh functional for nonlinear systems

V. S. Shchesnovich, S. B. Cavalcanti|ArXiv.org|Nov 12, 2004

Advanced Thermodynamics and Statistical Mechanics参考文献 13被引用数 11

ひとこと要約

本稿では、線形系におけるリーマン商の一般化として、非線形系におけるエネルギー汎関数の正規化制約下での数値的最小化を可能にする、レイリー汎関数を導入する。この手法により、汎関数の極値化によって定常解が特定され、第二変種が安定性を決定することが明らかになる。さらに、広く用いられている虚時間緩下法がこのアプローチの特殊な場合であることが示され、ボーズ＝アインシュタイン凝縮や臨界非線形方程式の普遍的解析的恒等式の導出に応用される。

ABSTRACT

We introduce Rayleigh functional for nonlinear systems. It is defined using the energy functional and the normalization properties of the variables of variation. The key property of the Rayleigh quotient for linear systems is preserved in our definition: the extremals of the Rayleigh functional coincide with the stationary solutions of the Euler-Lagrange equation. Moreover, the second variation of the Rayleigh functional defines stability of the solution. This gives rise to a powerful numerical optimization method in the search for the energy minimizers. It is shown that the well-known imaginary time relaxation is a special case of our method. To illustrate the method we find the stationary states of Bose-Einstein condensates in various geometries. Finally, we show that the Rayleigh functional also provides a simple way to derive analytical identities satisfied by the stationary solutions of the critical nonlinear equations. We also show that the functional used in Phys. Rev. A 66, 036612 (2002) is erroneous.

研究の動機と目的

• 非線形偏微分方程式、特に非線形シュレーディンガー（NLS）型方程式の定常解を系統的に数値的に求めるための手法を開発すること。
• レイリー汎関数を用いて正規化制約を自然に組み込む統一的枠組みを提供し、臨時的な正規化手順を回避すること。
• レイリー汎関数と広く用いられている虚時間緩下法との関係を確立し、後者がこのアプローチの特殊な場合であることを示すこと。
• スケール不変性とレイリー汎関数を用いて、臨界非線形方程式の定常解が満たす解析的恒等式を導出すること。
• レイリー汎関数の第二変種分析を用いて、非線形系における局所的最小解と鞍点を区別できるようにすること。

提案手法

• オーダーパラメータのL2ノルムの二乗に対するエネルギー汎関数の比としてレイリー汎関数を定義し、線形系におけるリーマン商の極値性を非線形系へ一般化する。
• レイリー汎関数の変分導関数を用いて、定常解のオイラー＝ラグランジュ方程式を導出し、元のエネルギー最小化問題と整合性を保つ。
• 虚時間における勾配ベースの最小化スキームを実装し、各ステップでエネルギーの負の変分導関数に基づく発展方程式を導出する。
• 2次元および1次元のトラップを含むさまざまな幾何配置におけるボーズ＝アインシュタイン凝縮の基底状態を同定するためにこの手法を適用する。
• エネルギー汎関数のスケール不変性を用いて、臨界NLS方程式における普遍的恒等式（たとえば運動エネルギーと非線形エネルギー項のバランス）を導出する。
• レイリー汎関数の第二変種が、制約付きラグランジアンの第二変種に対応することを示し、解の安定性分類を可能にする。

実験結果

リサーチクエスチョン

• RQ1 正規化制約下でのエネルギー最小化を可能にするために、非線形系におけるリーマン商の一般化を定式化できるか？
• RQ2 レイリー汎関数は、計算物理学で広く用いられている虚時間緩下法とどのように関係しているか？
• RQ3 スケール不変性とレイリー汎関数を用いて、臨界非線形方程式の定常解が満たす解析的恒等式を導出できるか？
• RQ4 レイリー汎関数の第二変種を用いて、非線形系における局所的最小解と鞍点を区別できるか？
• RQ5 レイリー汎関数は、グロス＝ピタエフスキー方程式やNLS方程式などの非線形PDEに対して、数値的および解析的アプローチをどの程度統合できるか？

主な発見

• レイリー汎関数の極値は、エネルギー汎関数から導かれるオイラー＝ラグランジュ方程式の定常解と完全に一致し、非線形系におけるリーマン商の性質が一般化される。
• レイリー汎関数の第二変種は、定常解の安定性を正しく特定でき、局所的最小解と鞍点を区別することが可能になる。
• 虚時間緩下法がレイリー汎関数最小化の特殊な場合であることが形式的に示され、各ステップで正規化が適用される。
• 臨界的な2次元NLS方程式に対して、スケール不変性とレイリー汎関数を用いて、恒等式 $ \int |\nabla\psi|^2 \, d^2x = \frac{1}{2} \int |\psi|^4 \, d^2x $ が普遍的に導出される。
• 1次元臨界的なNLS方程式においても、同じ手法を用いて恒等式 $ -\mu \int |\psi|^2 \, dx = \int \left| \frac{d\psi}{dx} \right|^2 \, dx = \frac{1}{3} \int |\psi|^6 \, dx $ が導出され、普遍性が示される。
• 本手法は、スケール不変性とレイリー汎関数の構造に基づいて、臨界非線形方程式のエネルギーバランス恒等式を体系的かつ解析的に導出する道筋を提供する。