A two-body correlated basis set is used to develop a many-body theory which is valid for any number of bosons in the trap. The formalism incorporates the van der Waals interaction and two-body correlations in an exact way. The theory has successfully been applied to Bose-Einstein condensates---dilute weakly interacting and also dilute but having a large scattering length. Even in the extreme dilute condition, we observe the breakdown of the shape-independent approximation and the interatomic correlation plays an important role in the large particle-number limit. This correlated many-body calculation can handle, within the two-body correlation approximation, the entire range of atom number of experimentally achieved condensates. Next we successfully push the basis function for large scattering lengths where the mean-field results are manifestly bad. The sharp increase in correlation energy clearly shows the beyond-mean-field effect. We also calculate one-particle densities for various scattering lengths and particle numbers. Our many-body calculation exhibits the finite-size effect in the one-body density.
We report the effects of interaction on thermodynamic properties of a repulsive Bose-Einstein condensate confined in a harmonic trap by using the correlated potential harmonics expansion method. This many-body technique permits the use of a realistic interactomic interaction, which gives rise to the effective long-range interaction of the condensate in terms of the
s-wave scattering length. We have computed temperature (
T) dependence of the chemical potential, specific heat, condensate fraction, entropy, pressure, and the average energy per particle of a system containing a large number (
A) of
87Rb atoms in the Joint Institute for Laboratory Astrophysics (JILA) trap. The repulsion among the interacting bosons results in a small but measurable drop of condensate fraction and critical temperature (
Tc), compared to those of a noninteracting condensate. These are in agreement with the experiment. Although all thermodynamic quantities have a strong dependence on
A and to a smaller extent on the interatomic interaction, our numerical calculation appears to show that a thermodynamic quantity per particle follows a universal behavior as a function of
T/Tc. This shows the importance of
Tc for all thermodynamic properties of the condensate. As expected, for
T>Tc, these properties follow those of a trapped noncondensed Bose gas.
We evaluate the zero temperature properties and Bose-Einstein condensation (BEC) of hard sphere bosons in a harmonic trap using Quantum Monte Carlo methods. The chief purpose of this work is to go beyond the dilute limit, to test the limits of the Gross-Pitaevskii equation and related mean field approximations and to explore the role of interactions in determining the zero temperature properties of the trapped Bose gas. The density is increased by adjusting both the number of trapped bosons, N , and the s-wave scattering length, a, to study systems from the highly dilute region corresponding to early experiments in trapped BEC [2, 3] up to liquid He densities and beyond. Rather than relying on the perturbative Bogoliubov approximation (which assumes a large condensate fraction) to describe the condensate, we use the one-body density matrix (OBDM) formulation of BEC so that the properties of the condensate for systems with arbitrarily large interactions may be studied. In this formulation of BEC, condensate properties are obtained by diagonalizing the OBDM and obtaining the corresponding single particle “natural orbitals” and their occupation numbers for the system. The condensate wave function and condensate fraction are then obtained from the single particle orbital(s) with macroscopic occupation (N0 >> 1). Within this framework, we calculate the effects of interactions and increased density on the ground state energy, the density profile, the momentum distribution, the condensate fraction and condensate “wave-function” and several other properties.
The diffusion Monte Carlo method is applied to describe a trapped atomic Bose-Einstein condensate at zero temperature, fully quantum mechanically and nonperturbatively. For low densities, ${n(0)a}^{3}<~2\ifmmode\times\else\texttimes\fi{}{10}^{\ensuremath{-}3}$
[n(0), peak density; a, s-wave scattering length], our calculations confirm that the exact ground-state energy for a sum of two-body interactions depends to a good approximation on only the atomic-physics parameter a, and no other details of the two-body model potential. Corrections to the mean-field Gross-Pitaevskii energy range from being essentially negligible to about 20% for
N=2−−50 particles in the trap with positive s-wave scattering length
a=100−−10000a.u. Our numerical calculations confirm that inclusion of an additional effective potential term in the mean-field equation, which accounts for quantum fluctuations [see, e.g., E. Braaten and A. Nieto, Phys. Rev. B 56, 14 745 (1997)], leads to a greatly improved description of trapped Bose gases.
We calculate the ground-state energy and the collective excitation frequency of trapped bosons at large scattering length interacting via the realistic two-body van der Waals potential. Our many-body method keeps two-body correlations produced by all interacting pairs. When the scattering length is small compared to the trap size and the number of bosons in the trap is of the order of a few thousands, the mean-field results are in good agreement with the many-body results. However for large particle numbers, even when the condensate is sufficiently dilute, the interatomic correlation comes into the picture. When the scattering length is quite large near the Feshbach resonance, the Bose gas becomes highly correlated. The many-body results are close to the Gross-Pitaevskii results for a small number of bosons, however, large deviations are noted in the large particle limit. We also calculate the lowest collective excitation and the interaction energy for large scattering lengths. The monopole excitation frequency exhibits a pronounced dependence on the scattering length. We also observe a universal behavior for the interaction energy at the limit of large scattering length.