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

arXiv:cond-mat/9811344 · Physical Review A

Numerical approach to the ground and excited states of a Bose-Einstein condensed gas confined in a completely anisotropic trap

Barry I. Schneider, David L. Feder

Published Mar 1, 1999• 84 citations
arXivHTMLPDF

Abstract

The ground and excited states of a weakly interacting and dilute Bose-Einstein condensed gas, confined in a completely anisotropic harmonic oscillator potential, are determined at zero temperature within the Bogoliubov approximation. The numerical calculations employ a computationally efficient procedure based on a discrete variable representation (DVR) of the Hamiltonian. The DVR is efficient for problems where the interaction potential may be expressed as a local function of interparticle coordinates. In order to address condensates that are both very large ($\ensuremath{\sim}{10}^{6}$ atoms) and fully anisotropic, the ground state is found using a self-consistent field approach. Experience has demonstrated, however, that standard iterative techniques applied to the solution of the nonlinear partial differential equation for the condensate are nonconvergent. This limitation is overcome using the method of direct inversion in the iterated subspace (DIIS). In addition, the sparse structure of the DVR enables the efficient application of iterative techniques such as the Davidson and/or Lanczos methods, to extract the eigenvalues of physical interest. The results are compared with recent experimental data obtained for Bose-Einstein condensed alkali-metal vapors confined in magnetic traps.

DOI
10.1103/physreva.59.2232

Related papers

Bose condensates in a harmonic trap near the critical temperature

T. H. Bergeman, David L. Feder, Nandor L. Balazs +1 more

The mean-field properties of finite-temperature Bose-Einstein gases confined in spherically symmetric harmonic traps are surveyed numerically. The solutions of the Gross-Pitaevskii (GP) and Hartree-Fock-Bogoliubov (HFB) equations for the condensate and low-lying quasiparticle excitations are calculated self-consistently using the discrete variable representation, while the most high-lying states are obtained with a local-density approximation. Consistency of the theory for temperatures through the Bose condensation point ${T}_{c}$ requires that the thermodynamic chemical potential differ from the eigenvalue of the GP equation; the appropriate modifications lead to results that are continuous as a function of the particle interactions. The HFB equations are made gapless either by invoking the Popov approximation or by renormalizing the particle interactions. The latter approach effectively reduces the strength of the effective scattering length ${a}_{\mathrm{sc}},$ increases the number of condensate atoms at each temperature, and raises the value of ${T}_{c}$ relative to the Popov approximation. The renormalization effect increases approximately with the log of the atom number, and is most pronounced at temperatures near ${T}_{c}.$ Comparisons with the results of quantum Monte Carlo calculations and various local-density approximations are presented, and experimental consequences are discussed.

arXiv:cond-mat/0001079•Physical Review A•May 11, 2000•31 citations

Ground State Properties for Bose-Condensed Gas in Anisotropic Traps

Xu Zhi-Jun, Haibin Li, Xin Xiao-Tian +2 more

Based on the energy functional and variational method, we present a new method to investigate the ground state properties for a weakly interacting Bose-condensed gas in an anisotropic harmonic trap at zero temperature. With this method we are able to find the analytic expression of the ground-state wavefunction and to explore the relevant quantities, such as energy, chemical potential, and the aspect ratio of the velocity distribution. These results agree well with previous ground state numerical solutions of the Gross–Pitaevskii equation given by Dalfovo et al. [Phys. Rev. A 53 (1996) 2477] This new method is simple compared to other methods used to solve numerically the Gross–Pitaevskii equation, and one can obtain analytic and reliable results.

Chinese Physics Letters•May 30, 2006•2 citations

Noncondensate atoms in a trapped Bose gas

Juha Javanainen

We formulate Bogoliubov theory for a weakly interacting Bose-Einstein condensate confined to a spherically symmetric harmonic-oscillator potential. The theory is solved numerically by diagonalizing a Hamiltonian that is second order in boson creation and annihilation operators with the aid of a symplectic transformation. At zero temperature, the number of noncondensate atoms is well approximated by a heuristic expression based on Bogoliubov theory for free atoms.

Physical Review A•Nov 1, 1996•37 citations

Analytical calculations on Landau damping of collective modes in anisotropic Bose-Einstein condensates

Xiaodong Ma, Yong-li Ma, Guoxiang Huang

We develop an analytical approach for calculating the Landau damping of collective modes in the Bose-Einstein condensate trapped by an anisotropic harmonic potential. Based on a variational ground-state wave function obtained by solving the time-independent Gross-Pitaevskii equation beyond Thomas-Fermi approximation, we solve the Bogoliubov--de Gennes equations describing thermal excited quasiparticles and provide divergence-free analytical solutions for the Bogoliubov amplitudes and coupling matrix elements that characterize the interaction between the collective modes and the quasiparticles. With these analytical results we evaluate the Landau damping rates of several collective modes for various anisotropic parameters of the trapping potential in terms of the formulas derived from a time-dependent mean-field theory. In addition, we discuss the dependence of the damping rates on temperature, particle number, trapping frequency, and anisotropic parameter of the system, and compare our theoretical results with experimental and numerical data reported in literature.

Physical Review A•Jan 31, 2007•5 citations

Dark-soliton states of Bose-Einstein condensates in anisotropic traps

David L. Feder, M. S. Pindzola, Lee A. Collins +2 more

Dark soliton states of Bose-Einstein condensates in harmonic traps are studied both analytically and computationally by the direct solution of the Gross-Pitaevskii equation in three dimensions. The ground and self-consistent excited states are found numerically by relaxation in imaginary time. The energy of a stationary soliton in a harmonic trap is shown to be independent of density and geometry for large numbers of atoms. Large-amplitude field modulation at a frequency resonant with the energy of a dark soliton is found to give rise to a state with multiple vortices. The Bogoliubov excitation spectrum of the soliton state contains complex frequencies, which disappear for sufficiently small numbers of atoms or large transverse confinement. The relationship between these complex modes and the snake instability is investigated numerically by propagation in real time.

arXiv:cond-mat/0004449•Physical Review A•Oct 19, 2000•270 citations

Excited states of a dilute Bose-Einstein condensate in a harmonic trap

Alexander L. Fetter, Daniel S. Rokhsar

The low-lying hydrodynamic normal modes of a dilute Bose-Einstein gas in an isotropic harmonic trap determine the corresponding Bogoliubov amplitudes. In the Thomas-Fermi limit, these modes have large low-temperature occupation numbers and they permit an explicit construction of the dynamic structure function $S(\mathbf{q},\ensuremath{\omega})$. The total noncondensate number ${N}^{\ensuremath{'}}(0)$ at zero temperature increases like ${R}^{6}$, where $R$ is the condensate radius measured in units of the oscillator length. The lowest dipole modes are constructed explicitly in the Bogoliubov approximation.

arXiv:cond-mat/9704234•Physical Review A•Feb 1, 1998•62 citations
The ground and excited states of a weakly interacting and dilute Bose-Einstein condensed gas, confined in a completely anisotropic harmonic oscillator potential, are determined at zero temperature within the Bogoliubov approximation. The numerical calculations employ a computationally efficient procedure based on a discrete variable representation (DVR) of the Hamiltonian. The DVR is efficient for problems where the interaction potential may be expressed as a local function of interparticle coordinates. In order to address condensates that are both very large (\ensuremath∼106\ensuremath{\sim}{10}^{6}\ensuremath∼106 atoms) and fully anisotropic, the ground state is found using a self-consistent field approach. Experience has demonstrated, however, that standard iterative techniques applied to the solution of the nonlinear partial differential equation for the condensate are nonconvergent. This limitation is overcome using the method of direct inversion in the iterated subspace (DIIS). In addition, the sparse structure of the DVR enables the efficient application of iterative techniques such as the Davidson and/or Lanczos methods, to extract the eigenvalues of physical interest. The results are compared with recent experimental data obtained for Bose-Einstein condensed alkali-metal vapors confined in magnetic traps.
The mean-field properties of finite-temperature Bose-Einstein gases confined in spherically symmetric harmonic traps are surveyed numerically. The solutions of the Gross-Pitaevskii (GP) and Hartree-Fock-Bogoliubov (HFB) equations for the condensate and low-lying quasiparticle excitations are calculated self-consistently using the discrete variable representation, while the most high-lying states are obtained with a local-density approximation. Consistency of the theory for temperatures through the Bose condensation point Tc{T}_{c}Tc​ requires that the thermodynamic chemical potential differ from the eigenvalue of the GP equation; the appropriate modifications lead to results that are continuous as a function of the particle interactions. The HFB equations are made gapless either by invoking the Popov approximation or by renormalizing the particle interactions. The latter approach effectively reduces the strength of the effective scattering length asc,{a}_{\mathrm{sc}},asc​, increases the number of condensate atoms at each temperature, and raises the value of Tc{T}_{c}Tc​ relative to the Popov approximation. The renormalization effect increases approximately with the log of the atom number, and is most pronounced at temperatures near Tc.{T}_{c}.Tc​. Comparisons with the results of quantum Monte Carlo calculations and various local-density approximations are presented, and experimental consequences are discussed.
Based on the energy functional and variational method, we present a new method to investigate the ground state properties for a weakly interacting Bose-condensed gas in an anisotropic harmonic trap at zero temperature. With this method we are able to find the analytic expression of the ground-state wavefunction and to explore the relevant quantities, such as energy, chemical potential, and the aspect ratio of the velocity distribution. These results agree well with previous ground state numerical solutions of the Gross–Pitaevskii equation given by Dalfovo et al. [Phys. Rev. A 53 (1996) 2477] This new method is simple compared to other methods used to solve numerically the Gross–Pitaevskii equation, and one can obtain analytic and reliable results.
We develop an analytical approach for calculating the Landau damping of collective modes in the Bose-Einstein condensate trapped by an anisotropic harmonic potential. Based on a variational ground-state wave function obtained by solving the time-independent Gross-Pitaevskii equation beyond Thomas-Fermi approximation, we solve the Bogoliubov--de Gennes equations describing thermal excited quasiparticles and provide divergence-free analytical solutions for the Bogoliubov amplitudes and coupling matrix elements that characterize the interaction between the collective modes and the quasiparticles. With these analytical results we evaluate the Landau damping rates of several collective modes for various anisotropic parameters of the trapping potential in terms of the formulas derived from a time-dependent mean-field theory. In addition, we discuss the dependence of the damping rates on temperature, particle number, trapping frequency, and anisotropic parameter of the system, and compare our theoretical results with experimental and numerical data reported in literature.
Dark soliton states of Bose-Einstein condensates in harmonic traps are studied both analytically and computationally by the direct solution of the Gross-Pitaevskii equation in three dimensions. The ground and self-consistent excited states are found numerically by relaxation in imaginary time. The energy of a stationary soliton in a harmonic trap is shown to be independent of density and geometry for large numbers of atoms. Large-amplitude field modulation at a frequency resonant with the energy of a dark soliton is found to give rise to a state with multiple vortices. The Bogoliubov excitation spectrum of the soliton state contains complex frequencies, which disappear for sufficiently small numbers of atoms or large transverse confinement. The relationship between these complex modes and the snake instability is investigated numerically by propagation in real time.
The low-lying hydrodynamic normal modes of a dilute Bose-Einstein gas in an isotropic harmonic trap determine the corresponding Bogoliubov amplitudes. In the Thomas-Fermi limit, these modes have large low-temperature occupation numbers and they permit an explicit construction of the dynamic structure function $S(\mathbf{q},\ensuremath{\omega})$. The total noncondensate number ${N}^{\ensuremath{'}}(0)$ at zero temperature increases like R6{R}^{6}R6, where RRR is the condensate radius measured in units of the oscillator length. The lowest dipole modes are constructed explicitly in the Bogoliubov approximation.