Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Communications in Mathematical Sciences

Ground, Symmetric and Central Vortex States in Rotating Bose-Einstein Condensates

Weizhu Bao, Peter A. Markowich, Hanquan Wang

Published Jan 1, 2005• 149 citations
Publisher

Abstract

We study ground, symmetric and central vortex states, as well as their energy and chemical potential diagrams, in rotating Bose-Einstein condensates (BEC) analytically and numerically.We start from the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with an angular momentum rotation term, scale it to obtain a four-parameter model, reduce it to a 2D GPE in the limiting regime of strong anisotropic confinement and present its semiclassical scaling and geometrical optics.We discuss the existence/nonexistence problem for ground states (depending on the angular velocity) and find that symmetric and central vortex states are independent of the angular rotational momentum.We perform numerical experiments computing these states using a continuous normalized gradient flow (CNGF) method with a backward Euler finite difference (BEFD) discretization.Ground, symmetric and central vortex states, as well as their energy configurations, are reported in 2D and 3D for a rotating BEC.Through our numerical study, we find various configurations with several vortices in both 2D and 3D structures, energy asymptotics in some limiting regimes and ratios between energies of different states in a strong replusive interaction regime.Finally we report the critical angular velocity at which the ground state loses symmetry, numerical verification of dimension reduction from 3D to 2D, errors for the Thomas-Fermi approximation, and spourous numerical ground states when the rotation speed is larger than the minimal trapping frequency in the xy plane.

DOI
10.4310/cms.2005.v3.n1.a5

Related papers

DYNAMICS OF THE GROUND STATE AND CENTRAL VORTEX STATES IN BOSE–EINSTEIN CONDENSATION

Weizhu Bao, Yanzhi Zhang

In this paper, we study dynamics of the ground state and central vortex states in Bose–Einstein condensation (BEC) analytically and numerically. We show how to define the energy of the Thomas–Fermi (TF) approximation, prove that the ground state is a global minimizer of the energy functional over the unit sphere and all excited states are saddle points in linear case, derive a second-order ordinary differential equation (ODE) which shows that time-evolution of the condensate width is a periodic function with/without a perturbation by using the variance identity, prove that the angular momentum expectation is conserved in two dimensions (2D) with a radial symmetric trap and 3D with a cylindrical symmetric trap for any initial data, and study numerically stability of central vortex states as well as interaction between a few central vortices with winding numbers ±1 by a fourth-order time-splitting sine-pseudospectral (TSSP) method. The merit of the numerical method is that it is explicit, unconditionally stable, time reversible and time transverse invariant. Moreover, it conserves the position density, performs spectral accuracy for spatial derivatives and fourth-order accuracy for time derivative, and possesses "optimal" spatial/temporal resolution in the semiclassical regime. Finally we find numerically the critical angular frequency for single vortex cycling from the ground state under a far-blue detuned Gaussian laser stirrer in strong repulsive interaction regime and compare our numerical results with those in the literatures.

Mathematical Models and Methods in Applied Sciences•Dec 1, 2005•47 citations

Vortices in a rotating Bose-Einstein condensate: Critical angular velocities and energy diagrams in the Thomas-Fermi regime

Amandine Aftalion, Qiang Du

For a Bose-Einstein condensate placed in a rotating trap and strongly confined along the z axis, we set a framework of study for the Gross-Pitaevskii energy in the Thomas-Fermi regime for an effective two-dimensional (2D) situation in the $x\ensuremath{-}y$ plane. We investigate an asymptotic expansion of the energy, the critical angular velocities of nucleation of vortices with respect to a small parameter $\ensuremath{\varepsilon},$ and the location of vortices. The limit $\ensuremath{\varepsilon}$ going to zero corresponds to the Thomas-Fermi regime. The nondimensionalized energy is similar to the Ginzburg-Landau energy for superconductors in the high-\ensuremath{\kappa} high-field limit and our estimates rely on techniques developed for this latter problem. We also take advantage of this similarity to develop a numerical algorithm for computing the Bose-Einstein vortices. Numerical results and energy diagrams are presented.

arXiv:cond-mat/0103299•Physical Review A•Nov 8, 2001•144 citations

Mathematical Models and Numerical Methods for Bose-Einstein Condensation

Weizhu Bao

The achievement of Bose-Einstein condensation (BEC) in ultracold vapors of alkali atoms has given enormous impulse to the theoretical and experimental study of dilute atomic gases in condensed quantum states inside magnetic traps and optical lattices. This article offers a short survey on mathematical models and theories as well as numerical methods for BEC based on the mean field theory. We start with the Gross-Pitaevskii equation (GPE) in three dimensions (3D) for modeling one-component BEC of the weakly interacting bosons, scale it to obtain a three-parameter model and show how to reduce it to two dimensions (2D) and one dimension (1D) GPEs in certain limiting regimes. Mathematical theories and numerical methods for ground states and dynamics of BEC are provided. Extensions to GPE with an angular momentum rotation term for a rotating BEC, to GPE with long-range anisotropic dipole-dipole interaction for a dipolar BEC and to coupled GPEs for spin-orbit coupled BECs are discussed. Finally, some conclusions are drawn and future research perspectives are discussed.

arXiv:1403.3884•Mar 16, 2014•10 citations

Vortices in Bose-Einstein condensates with dominant dipolar interactions

Marta Abad, M. Guilleumas, R. Mayol +2 more

We present full three-dimensional numerical calculations of single vortex states in rotating dipolar condensates. We consider a Bose-Einstein condensate of $^{52}\text{C}\text{r}$ atoms with dipole-dipole and $s$-wave contact interactions confined in an axially symmetric harmonic trap. We obtain the vortex states by numerically solving the Gross-Pitaevskii equation in the rotating frame with no further approximations. We investigate the properties of a single vortex and calculate the critical angular velocity for different values of the $s$-wave scattering length. We show that, whereas the standard variational approach breaks down in the limit of pure dipolar interactions, exact solutions of the Gross-Pitaevskii equation can be obtained for values of the $s$-wave scattering length down to zero. The energy barrier for the nucleation of a vortex is calculated as a function of the vortex displacement from the rotation axis for different values of the angular velocity of the rotating trap.

arXiv:0906.4259•Physical Review A•Jun 18, 2009•66 citations

Dynamics of Rotating Bose--Einstein Condensates and its Efficient and Accurate Numerical Computation

Weizhu Bao, Qiang Du, Yanzhi Zhang

In this paper, we study the dynamics of rotating Bose--Einstein condensates (BEC) based on the Gross--Pitaevskii equation (GPE) with an angular momentum rotation term and present an efficient and accurate algorithm for numerical simulations. We examine the conservation of the angular momentum expectation and the condensate width and analyze the dynamics of a stationary state with a shift in its center. By formulating the equation in either the two-dimensional polar coordinate system or the three-dimensional cylindrical coordinate system, the angular momentum rotation term becomes a term with constant coefficients. This allows us to develop an efficient time-splitting method which is time reversible, unconditionally stable, efficient, and accurate for the problem. Moreover, it conserves the position density. We also apply the numerical method to study issues such as the stability of central vortex states and the quantized vortex lattice dynamics in rotating BEC.

SIAM Journal on Applied Mathematics•Jan 1, 2006•116 citations

Mathematical theory and numerical methods for Bose-Einstein condensation

Weizhu Bao, Yongyong Cai

In this paper, we mainly review recent results on mathematical theory andnumerical methods for Bose-Einstein condensation (BEC),based on the Gross-Pitaevskii equation (GPE). Starting from the simplest case with one-component BEC of the weakly interacting bosons, we study the reduction of GPE to lower dimensions, the ground states of BEC including the existence and uniqueness as well as nonexistence results, and the dynamics of GPE including dynamical laws, well-posedness of the Cauchy problem as well as the finite time blow-up. To compute the ground state, the gradient flow with discrete normalization (or imaginary time) method is reviewed and various full discretization methods are presented and compared. To simulate the dynamics, both finite difference methods and time splitting spectral methods are reviewed, and their error estimates are briefly outlined. When the GPE has symmetric properties, we show how to simplify the numerical methods. Then we compare two widely used scalings, i.e. physical scaling (commonly used) and semiclassical scaling, for BEC in strong repulsive interaction regime (Thomas-Fermi regime), and discuss semiclassical limits of the GPE. Extensions of these results for one-component BEC are thencarried out for rotating BEC by GPE with an angular momentum rotation, dipolar BEC by GPE with long range dipole-dipole interaction, and two-component BEC by coupled GPEs. Finally, as a perspective, we show briefly the mathematical models forspin-1 BEC, Bogoliubov excitation and BEC at finite temperature.

arXiv:1212.5341•Kinetic and Related Models•Dec 11, 2012•591 citations
We study ground, symmetric and central vortex states, as well as their energy and chemical potential diagrams, in rotating Bose-Einstein condensates (BEC) analytically and numerically.We start from the three-dimensional (3D) Gross-Pitaevskii equation (GPE) with an angular momentum rotation term, scale it to obtain a four-parameter model, reduce it to a 2D GPE in the limiting regime of strong anisotropic confinement and present its semiclassical scaling and geometrical optics.We discuss the existence/nonexistence problem for ground states (depending on the angular velocity) and find that symmetric and central vortex states are independent of the angular rotational momentum.We perform numerical experiments computing these states using a continuous normalized gradient flow (CNGF) method with a backward Euler finite difference (BEFD) discretization.Ground, symmetric and central vortex states, as well as their energy configurations, are reported in 2D and 3D for a rotating BEC.Through our numerical study, we find various configurations with several vortices in both 2D and 3D structures, energy asymptotics in some limiting regimes and ratios between energies of different states in a strong replusive interaction regime.Finally we report the critical angular velocity at which the ground state loses symmetry, numerical verification of dimension reduction from 3D to 2D, errors for the Thomas-Fermi approximation, and spourous numerical ground states when the rotation speed is larger than the minimal trapping frequency in the xy plane.
In this paper, we study dynamics of the ground state and central vortex states in Bose–Einstein condensation (BEC) analytically and numerically. We show how to define the energy of the Thomas–Fermi (TF) approximation, prove that the ground state is a global minimizer of the energy functional over the unit sphere and all excited states are saddle points in linear case, derive a second-order ordinary differential equation (ODE) which shows that time-evolution of the condensate width is a periodic function with/without a perturbation by using the variance identity, prove that the angular momentum expectation is conserved in two dimensions (2D) with a radial symmetric trap and 3D with a cylindrical symmetric trap for any initial data, and study numerically stability of central vortex states as well as interaction between a few central vortices with winding numbers ±1 by a fourth-order time-splitting sine-pseudospectral (TSSP) method. The merit of the numerical method is that it is explicit, unconditionally stable, time reversible and time transverse invariant. Moreover, it conserves the position density, performs spectral accuracy for spatial derivatives and fourth-order accuracy for time derivative, and possesses "optimal" spatial/temporal resolution in the semiclassical regime. Finally we find numerically the critical angular frequency for single vortex cycling from the ground state under a far-blue detuned Gaussian laser stirrer in strong repulsive interaction regime and compare our numerical results with those in the literatures.
For a Bose-Einstein condensate placed in a rotating trap and strongly confined along the z axis, we set a framework of study for the Gross-Pitaevskii energy in the Thomas-Fermi regime for an effective two-dimensional (2D) situation in the $x\ensuremath{-}y$ plane. We investigate an asymptotic expansion of the energy, the critical angular velocities of nucleation of vortices with respect to a small parameter $\ensuremath{\varepsilon},$ and the location of vortices. The limit $\ensuremath{\varepsilon}$ going to zero corresponds to the Thomas-Fermi regime. The nondimensionalized energy is similar to the Ginzburg-Landau energy for superconductors in the high-\ensuremath{\kappa} high-field limit and our estimates rely on techniques developed for this latter problem. We also take advantage of this similarity to develop a numerical algorithm for computing the Bose-Einstein vortices. Numerical results and energy diagrams are presented.
The achievement of Bose-Einstein condensation (BEC) in ultracold vapors of alkali atoms has given enormous impulse to the theoretical and experimental study of dilute atomic gases in condensed quantum states inside magnetic traps and optical lattices. This article offers a short survey on mathematical models and theories as well as numerical methods for BEC based on the mean field theory. We start with the Gross-Pitaevskii equation (GPE) in three dimensions (3D) for modeling one-component BEC of the weakly interacting bosons, scale it to obtain a three-parameter model and show how to reduce it to two dimensions (2D) and one dimension (1D) GPEs in certain limiting regimes. Mathematical theories and numerical methods for ground states and dynamics of BEC are provided. Extensions to GPE with an angular momentum rotation term for a rotating BEC, to GPE with long-range anisotropic dipole-dipole interaction for a dipolar BEC and to coupled GPEs for spin-orbit coupled BECs are discussed. Finally, some conclusions are drawn and future research perspectives are discussed.
We present full three-dimensional numerical calculations of single vortex states in rotating dipolar condensates. We consider a Bose-Einstein condensate of 52Cr^{52}\text{C}\text{r}52Cr atoms with dipole-dipole and sss-wave contact interactions confined in an axially symmetric harmonic trap. We obtain the vortex states by numerically solving the Gross-Pitaevskii equation in the rotating frame with no further approximations. We investigate the properties of a single vortex and calculate the critical angular velocity for different values of the sss-wave scattering length. We show that, whereas the standard variational approach breaks down in the limit of pure dipolar interactions, exact solutions of the Gross-Pitaevskii equation can be obtained for values of the sss-wave scattering length down to zero. The energy barrier for the nucleation of a vortex is calculated as a function of the vortex displacement from the rotation axis for different values of the angular velocity of the rotating trap.
In this paper, we study the dynamics of rotating Bose--Einstein condensates (BEC) based on the Gross--Pitaevskii equation (GPE) with an angular momentum rotation term and present an efficient and accurate algorithm for numerical simulations. We examine the conservation of the angular momentum expectation and the condensate width and analyze the dynamics of a stationary state with a shift in its center. By formulating the equation in either the two-dimensional polar coordinate system or the three-dimensional cylindrical coordinate system, the angular momentum rotation term becomes a term with constant coefficients. This allows us to develop an efficient time-splitting method which is time reversible, unconditionally stable, efficient, and accurate for the problem. Moreover, it conserves the position density. We also apply the numerical method to study issues such as the stability of central vortex states and the quantized vortex lattice dynamics in rotating BEC.
In this paper, we mainly review recent results on mathematical theory andnumerical methods for Bose-Einstein condensation (BEC),based on the Gross-Pitaevskii equation (GPE). Starting from the simplest case with one-component BEC of the weakly interacting bosons, we study the reduction of GPE to lower dimensions, the ground states of BEC including the existence and uniqueness as well as nonexistence results, and the dynamics of GPE including dynamical laws, well-posedness of the Cauchy problem as well as the finite time blow-up. To compute the ground state, the gradient flow with discrete normalization (or imaginary time) method is reviewed and various full discretization methods are presented and compared. To simulate the dynamics, both finite difference methods and time splitting spectral methods are reviewed, and their error estimates are briefly outlined. When the GPE has symmetric properties, we show how to simplify the numerical methods. Then we compare two widely used scalings, i.e. physical scaling (commonly used) and semiclassical scaling, for BEC in strong repulsive interaction regime (Thomas-Fermi regime), and discuss semiclassical limits of the GPE. Extensions of these results for one-component BEC are thencarried out for rotating BEC by GPE with an angular momentum rotation, dipolar BEC by GPE with long range dipole-dipole interaction, and two-component BEC by coupled GPEs. Finally, as a perspective, we show briefly the mathematical models forspin-1 BEC, Bogoliubov excitation and BEC at finite temperature.