[论文解读] 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

一句话总结

本文提出了一种用于 Gross-Pitaevskii 基态的保正有限元方法，通过质量集中法保持离散唯一性与保正性。证明了所有非负离散激发态均与基态重合，从而确保梯度流求解器的全局收敛性，并通过严格的先验误差分析建立了最优收敛速率。

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.

研究动机与目标

• 开发一种有限元离散化方法，以保持连续 Gross-Pitaevskii 基态的保正性与符号唯一性。
• 通过证明任意非负离散激发态必与离散基态重合，建立离散唯一性。
• 通过利用离散保正性与唯一性，确保梯度流求解器的全局收敛性。
• 进行严格的先验误差分析，表明所有离散未知量均具有最优收敛速率。

提出的方法

• 该方法在有限元离散化中采用质量集中法，将连续保正性与唯一性性质转移至离散格式。
• 采用 Gross-Pitaevskii 特征值问题的变分形式，作为 Gross-Pitaevskii 能量泛函的约束最小化问题。
• 通过梯度流迭代计算离散基态，其收敛性由所证明的离散唯一性与保正性保证。
• 利用椭圆正则性、节点插值及离散与连续范数之间的等价性，推导先验误差估计。
• 分析依赖于 Riesz 表示定理与嵌入定理，以在网格尺寸上一致地有界离散解的 $ L^\infty $-范数。
• 该方法确保所有非负离散解均为基态，从而可识别迭代求解器的极限。

实验结果

研究问题

• RQ1 有限元离散化能否保持连续 Gross-Pitaevskii 基态的保正性与符号唯一性？
• RQ2 在有限元格式下，每个非负离散激发态是否均与离散基态重合？
• RQ3 在离散保正性与唯一性条件下，能否保证梯度流求解器的全局收敛性？
• RQ4 所提出的非标准有限元方法对所有未知量的最优收敛速率是多少？
• RQ5 在能量有界条件下，如何建立离散解的统一 $ L^\infty $-界？

主要发现

• 离散基态在符号唯一性意义下唯一为正，且所有非负离散激发态均与之重合，从而确保梯度流极限即为基态。
• 该方法在所有范数下均达到最优收敛速率：$ \|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 $。
• 由于离散唯一性与保正性，所提方法确保梯度流求解器全局收敛至离散基态。
• 误差分析证实能量、梯度与状态逼近均达到最优收敛，其界依赖于精确解的正则性。
• 该方法在质量集中下具有鲁棒性，保持了连续问题在离散格式下的关键理论性质。