We formulate a method to study two-body correlations in a system of N identical bosons interacting via central two-body potentials. We use the adiabatic hyperspherical approach and assume a Faddeev-like decomposition of the wave function. For a fixed hyperradius we derive variationally an optimal integro-differential equation for the hyperangular eigenvalue and wave function. This equation reduces substantially by assuming the interaction range much smaller than the size of the N-body system. At most, one-dimensional integrals then remain. We view a Bose-Einstein condensate pictorially as a structure in the landscape of the potential given as a function of the one-dimensional hyperradial coordinate. The quantum states of the condensate can be located in one of the two potential minima. We derive and discuss properties of the solutions and illustrate with numerical results. The correlations lower the interaction energy substantially. The new multibody Efimov states are solutions independent of details of the two-body potential. We compare with mean-field results and available experimental data.
We investigate systems of identical bosons with the focus on two-body correlations. We use the hyperspherical adiabatic method and a decomposition of the wave function in two-body amplitudes. An analytic parametrization is used for the adiabatic effective radial potential. We discuss the structure of a condensate for arbitrary scattering length. Stability and time scales for various decay processes are estimated. The previously predicted Efimov-like states are found to be very narrow. We discuss the validity conditions and formal connections between the zero- and finite-range mean-field approximations, Faddeev-Yakubovskii formulation, Jastrow ansatz, and the present method. We compare numerical results from present work with mean-field calculations and discuss qualitatively the connection with measurements.
We investigate systems of identical bosons with the focus on two-body correlations and attractive finite-range potentials. We use a hyperspherical adiabatic method and apply a Faddeev type of decomposition of the wave function. We discuss the structure of a condensate as a function of particle number and scattering length. We establish universal scaling relations for the critical effective radial potentials for distances where the average distance between particle pairs is larger than the interaction range. The correlations in the wave function restore the large-distance mean-field behavior with the correct two-body interaction. We discuss various processes limiting the stability of condensates. With correlations we confirm that macroscopic tunneling dominates when the trap length is about half of the particle number times the scattering length.
We study two-body correlations for N identical bosons by use of the hyperspherical adiabatic-expansion method. We use the zero-range interaction and derive a transcendental equation determining the key ingredient of the hyperradial potential. The necessary renormalization is for both repulsive and attractive interactions achieved with an effective range expansion of the two-body phase shifts. Our solutions including correlations provide the properties of Bose-Einstein condensates exemplified by stability conditions as established by mean-field Gross-Pitaevskii calculations. The many-body Efimov states are unavoidable for large scattering lengths.
We study a system of
A identical interacting bosons trapped by an external field by solving ab initio the many-body Schr\"odinger equation. A complete solution by using, for example, the traditional hyperspherical harmonics (HH) basis develops serious practical problems due to the large degeneracy of HH basis. Symmetrization of the wave function, calculation of the matrix elements, etc., become an immensely formidable task as
A increases. Instead of the HH basis, here we use a new basis, called ``potential harmonics'' (PH) basis, which is a subset of HH basis. We assume that the contribution to the orbital and grand orbital [in $3(A\ensuremath{-}1)$-dimensional space of the reduced motion] quantum numbers comes only from the interacting pair. This implies inclusion of two-body correlations only and disregard of all higher-body correlations. Such an assumption is ideally suited for the Bose-Einstein condensate (BEC), which is required, for experimental realization of BEC, to be extremely dilute. Hence three and higher-body collisions are almost totally absent. Unlike the $(3A\ensuremath{-}4)$ hyperspherical variables in HH basis, the PH basis involves only three active variables, corresponding to three quantum numbers---the orbital
l, azimuthal
m, and the grand orbital
2K+l quantum numbers for any arbitrary
A. It drastically reduces the number of coupled equations and calculation of the potential matrix becomes tremendously simplified, as it involves integrals over only three variables for any
A. One can easily incorporate realistic atom-atom interactions in a straightforward manner. We study the ground and excited state properties of the condensate for both attractive and repulsive interactions for various particle number. The ground state properties are compared with those calculated from the Gross-Pitaevskii equation. We notice that our many-body results converge towards the mean field results as the particle number increases.