English

An Efficient Unconditionally Energy-Stable Numerical Scheme for Bose--Einstein Condensate

Numerical Analysis 2025-11-18 v1 Numerical Analysis

Abstract

A numerical framework is proposed and analyzed for computing the ground state of Bose--Einstein condensates. A gradient flow approach is developed, incorporating both a Lagrange multiplier to enforce the L2L^2 conservation and a free energy dissipation. An explicit approximation is applied to the chemical potential, combined with an exponential time differencing (ETD) operator to the diffusion part, as well a stabilizing operator, to obtain an intermediate numerical profile. Afterward, an L2L^2 normalization is applied at the next numerical stage. A theoretical analysis reveals a free energy dissipation under a maximum norm bound assumption for the numerical solution, and such a maximum norm bound could be recovered by a careful convergence analysis and error estimate. In the authors' knowledge, the proposed method is the first numerical work that preserves the following combined theoretical properties: (1) an explicit computation at each time step, (2) unconditional free energy dissipation, (3) L2L^2 norm conservation at each time step, (4) a theoretical justification of convergence analysis and optimal rate error estimate. Comprehensive numerical experiments validate these theoretical results, demonstrating excellent agreement with established reference solutions.

Keywords

Cite

@article{arxiv.2511.12411,
  title  = {An Efficient Unconditionally Energy-Stable Numerical Scheme for Bose--Einstein Condensate},
  author = {Jing Guo and Cheng Wang and Dong Wang},
  journal= {arXiv preprint arXiv:2511.12411},
  year   = {2025}
}
R2 v1 2026-07-01T07:39:25.805Z