Computing the Ground State Solution of Bose-Einstein Condensates by a Normalized Gradient Flow

In this paper, we present a continuous normalized gradient flow (CNGF) and prove its energy diminishing property, which provides a mathematical justification of the imaginary time method used in the physics literature to compute the ground state solution of Bose--Einstein condensates (BEC). We also investigate the energy diminishing property for the discretization of the CNGF. Two numerical methods are proposed for such discretizations: one is the backward Euler centered finite difference (BEFD) method, the other is an explicit time-splitting sine-spectral (TSSP) method. Energy diminishing for BEFD and TSSP for the linear case and monotonicity for BEFD for both linear and nonlinear cases are proven. Comparison between the two methods and existing methods, e.g., Crank--Nicolson finite difference (CNFD) or forward Euler finite difference (FEFD), shows that BEFD and TSSP are much better in terms of preserving the energy diminishing property of the CNGF. Numerical results in one, two, and three dimensions with magnetic trap confinement potential, as well as a potential of a stirrer corresponding to a far-blue detuned Gaussian laser beam, are reported to demonstrate the effectiveness of BEFD and TSSP methods. Furthermore we observe that the CNGF and its BEFD discretization can also be applied directly to compute the first excited state solution in BEC when the initial data is chosen as an odd function.

[1]  Succi,et al.  Ground state of trapped interacting bose-einstein condensates by an explicit imaginary-time algorithm , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[2]  N. SIAMJ. AN EXPLICIT UNCONDITIONALLY STABLE NUMERICAL METHOD FOR SOLVING DAMPED NONLINEAR SCHRÖDINGER EQUATIONS WITH A FOCUSING NONLINEARITY , 2003 .

[3]  Holland,et al.  Time-dependent solution of the nonlinear Schrödinger equation for Bose-condensed trapped neutral atoms. , 1995, Physical review. A, Atomic, molecular, and optical physics.

[4]  E. M. Lifshitz,et al.  Quantum mechanics: Non-relativistic theory, , 1959 .

[5]  Elliott H. Lieb,et al.  Bosons in a trap: A rigorous derivation of the Gross-Pitaevskii energy functional , 1999, math-ph/9908027.

[6]  Terence Tao,et al.  The nonlinear Schr , 2005 .

[7]  C. J. Pethick,et al.  Vortices in Bose-Einstein-condensed atomic clouds , 1998 .

[8]  Qiang Du,et al.  Ginzburg-Landau vortices: dynamics, pinning, and hysteresis , 1997 .

[9]  R. J. Dodd,et al.  Approximate Solutions of the Nonlinear Schrödinger Equation for Ground and Excited States of Bose-Einstein Condensates , 1996, Journal of research of the National Institute of Standards and Technology.

[10]  P. Markowich,et al.  On time-splitting spectral approximations for the Schrödinger equation in the semiclassical regime , 2002 .

[11]  Gadi Fibich,et al.  Self-Focusing in the Perturbed and Unperturbed Nonlinear Schrödinger Equation in Critical Dimension , 1999, SIAM J. Appl. Math..

[12]  A Gammal,et al.  Improved numerical approach for the time-independent Gross-Pitaevskii nonlinear Schrödinger equation. , 1999, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[13]  S. Adhikari Numerical solution of the two-dimensional Gross-Pitaevskii equation for trapped interacting atoms , 2000, cond-mat/0001361.

[14]  Wolfgang Ketterle,et al.  Bose–Einstein condensation of atomic gases , 2002, Nature.

[15]  T. Hänsch,et al.  Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms , 2002, Nature.

[16]  C. Sulem,et al.  The nonlinear Schrödinger equation : self-focusing and wave collapse , 2004 .

[17]  L. Simon Asymptotics for a class of non-linear evolution equations, with applications to geometric problems , 1983 .

[18]  Sauro Succi,et al.  Particle-inspired scheme for the Gross-Pitaevski equation: An application to Bose-Einstein condensation , 2000 .

[19]  Eric Cornell,et al.  Very Cold Indeed: The Nanokelvin Physics of Bose-Einstein Condensation , 1996, Journal of research of the National Institute of Standards and Technology.

[20]  Edwards,et al.  Numerical solution of the nonlinear Schrödinger equation for small samples of trapped neutral atoms. , 1995, Physical review. A, Atomic, molecular, and optical physics.

[21]  Succi,et al.  Numerical solution of the gross-pitaevskii equation using an explicit finite-difference scheme: An application to trapped bose-einstein condensates , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[22]  E. L. Short,et al.  Quantum Chemistry , 1969, Nature.

[23]  Weizhu Bao,et al.  An Explicit Unconditionally Stable Numerical Method for Solving Damped Nonlinear Schrödinger Equations with a Focusing Nonlinearity , 2003, SIAM Journal on Numerical Analysis.

[24]  Jinchao Xu,et al.  Recent Progress in Computational and Applied PDES , 2002 .

[25]  F. Dalfovo,et al.  Theory of Bose-Einstein condensation in trapped gases , 1998, cond-mat/9806038.

[26]  R. J. Ballagh,et al.  Three-dimensional vortex dynamics in Bose-Einstein condensates , 2000 .

[27]  Qiang Du,et al.  Vortices in a rotating Bose-Einstein condensate: Critical angular velocities and energy diagrams in the Thomas-Fermi regime , 2001 .

[28]  R. J. Ballagh,et al.  Coherent Dynamics of Vortex Formation in Trapped Bose-Einstein Condensates , 1999 .

[29]  B. I. Schnieder,et al.  Numerical approach to the ground and excited states of a Bose-Einstein condensed gas confined in a completely anisotropic trap , 1999 .

[30]  Vladimir V. Konotop,et al.  Solutions of Gross-Pitaevskii equations beyond the hydrodynamic approximation: Application to the vortex problem , 2000 .

[31]  C. Gardiner,et al.  Cold Bosonic Atoms in Optical Lattices , 1998, cond-mat/9805329.

[32]  C. Wieman,et al.  Observation of Bose-Einstein Condensation in a Dilute Atomic Vapor , 1995, Science.

[33]  E. Gross Structure of a quantized vortex in boson systems , 1961 .

[34]  Adhikari Numerical study of the spherically symmetric gross-pitaevskii equation in two space dimensions , 2000, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[35]  Charles S. Adams,et al.  VORTEX FORMATION IN DILUTE INHOMOGENEOUS BOSE-EINSTEIN CONDENSATES , 1998 .

[36]  Daniel S. Rokhsar,et al.  Vortex Stability and Persistent Currents in Trapped Bose Gases , 1997 .

[37]  Shi Jin,et al.  Numerical Study of Time-Splitting Spectral Discretizations of Nonlinear Schrödinger Equations in the Semiclassical Regimes , 2003, SIAM J. Sci. Comput..

[38]  P. Markowich,et al.  Numerical solution of the Gross--Pitaevskii equation for Bose--Einstein condensation , 2003, cond-mat/0303239.

[39]  F. Dalfovo,et al.  Bosons in anisotropic traps: Ground state and vortices. , 1995, Physical review. A, Atomic, molecular, and optical physics.