[論文レビュー] Positivity preserving finite element method for the Gross-Pitaevskii ground state: discrete uniqueness and global convergence

Moritz Hauck, Yizhou Liang+1 more|arXiv (Cornell University)|May 27, 2024

Nonlinear Photonic Systems被引用数 4

ひとこと要約

本稿では、質量積分法を用いて離散的一意性および正性を保つようにしたグロス＝ピタエフスキー地面状態の正性を保つ有限要素法を提案する。非負の離散的励起状態がすべて地面状態に一致することを証明し、勾配流れソルバーのグローバル収束を保証するとともに、厳密な事前誤差解析により最適収束率を確立する。

ABSTRACT

We propose a positivity preserving finite element discretization for the nonlinear Gross-Pitaevskii eigenvalue problem. The method employs mass lumping techniques, which allow to transfer the uniqueness up to sign and positivity properties of the continuous ground state to the discrete setting. We further prove that every non-negative discrete excited state up to sign coincides with the discrete ground state. This allows one to identify the limit of fully discretized gradient flows, which are typically used to compute the discrete ground state, and thereby establish their global convergence. Furthermore, we perform a rigorous a priori error analysis of the proposed non-standard finite element discretization, showing optimal orders of convergence for all unknowns. Numerical experiments illustrate the theoretical results of this paper.

研究の動機と目的

• 連続的なグロス＝ピタエフスキー地面状態の正性および符号を除く一意性を保つ有限要素離散化を開発すること。
• 非負の離散的励起状態が離散的地面状態に一致することを証明することで、離散的一意性を確立すること。
• 離散的正性および一意性を活用して、勾長流れソルバーのグローバル収束を保証すること。
• すべての離散未知数について最適収束率を示す厳密な事前誤差解析を実施すること。

提案手法

• 連続的な正性および一意性の性質を離散的設定に移すために、有限要素離散化において質量積分法を採用する。
• グロス＝ピタエフスキー固有値問題の変分定式化を、グロス＝ピタエフスキーエネルギー汎関数の制約付き最小化として用いる。
• 離散的地面状態は勾配流れ反復法により計算され、証明された離散的一意性および正性により収束が保証される。
• 楕円的正則性、ノード内挿法、離散的および連続的ノルム間の等価性を用いて、事前誤差推定値を導出する。
• リーマン表現および埋め込み定理に依拠して、メッシュサイズに依存しない離散解の $ L^∞ $-ノルムを有界化する。
• すべての非負の離散解が地面状態に一致することを保証し、反復的ソルバーの極限の特定を可能にする。

実験結果

リサーチクエスチョン

• RQ1 有限要素離散化は、連続的なグロス＝ピタエフスキー地面状態の正性および符号を除く一意性を保てるか？
• RQ2 有限要素設定におけるすべての非負の離散的励起状態は、離散的地面状態に一致するか？
• RQ3 離散的正性および一意性のもとで、勾長流れソルバーのグローバル収束を保証できるか？
• RQ4 提案された非標準的有限要素法のすべての未知数について、最適収束率は何か？
• RQ5 エネルギーが有界である条件下で、離散解の均一な $ L^∞ $-有界性をどのように確立できるか？

主な発見

• 離散的地面状態は符号を除いて一意に正であり、すべての非負の離散的励起状態がこれに一致するため、勾長流れの極限が地面状態に収束することが保証される。
• すべてのノルムで最適収束率を達成する：$ \|u - u_h\|_{H^1} \lesssim h $、$ \|u - u_h\|_{L^2} \lesssim h^2 $、および $ \|u - u_h\|_{L^\infty} \lesssim h^{2-d/2} $（$ d \in \{1,2,3\} $）。
• エネルギーが有界である条件下で、離散解の均一な $ L^\infty $-有界性が確立され、$ \|u_h\|_{L^\infty} \lesssim 1 $ および $ \|\hat{u}_h\|_{L^\infty} \lesssim 1 $ が成り立つ。
• 離散的一意性および正性を活用することで、勾長流れソルバーが離散的地面状態にグローバルに収束することが保証される。
• 誤差解析により、エネルギー、勾配、状態近似について最適収束が確認され、正確な解の正則性に依存する有界性が得られる。
• 質量積分法に対しても本手法はロバストであり、連続問題の主要な理論的性質が離散的設定でも保たれる。