Robust and efficient preconditioned Krylov spectral solvers for computing the ground states of fast rotating and strongly interacting Bose-Einstein condensates

JOURNAL OF COMPUTATIONAL PHYSICS

卷:258

Article

Antoine, Xavier^1,2,3  Duboscq, Romain^1,2 

[1] Univ Lorraine, Inst Elie Cartan Lorraine, UMR 7502, F-54506 Vandoeuvre Les Nancy, France

[2] CNRS, Inst Elie Cartan Lorraine, UMR 7502, F-54506 Vandoeuvre Les Nancy, France

[3] INRIA Nancy Grand Est IECL ALICE, Nancy, France

关键词: Gross-Pitaevskii equations;    Bose-Einstein condensates;    Stationary states;    Vortex;    Rotation;    Spectral methods;    Krylov solvers;   

DOI  :  10.1016/j.jcp.2013.10.045

来源: Elsevier

【 摘 要 】

We consider the Backward Euler SPectral (BESP) scheme proposed in [10] for computing the stationary states of Bose-Einstein Condensates (BECs) through the Gross-Pitaevskii equation. We show that the fixed point approach introduced in [10] fails to converge for fast rotating BECs. A simple alternative approach based on Krylov subspace solvers with a Laplace or Thomas-Fermi preconditioner is given. Numerical simulations (obtained with the associated freely available Matlab toolbox GPELab) for complex configurations show that the method is accurate, fast and robust for 2D/3D problems and multi-components BECs. (C) 2013 Elsevier Inc. All rights reserved.

【 授权许可】

Free