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.
In this work, a direct solution to the few-body Integrodifferential equation using the Lagrange-mesh method is presented.With the Lagrange-mesh method, the few-body Integrodifferential equation is converted into a matrix eigenvalue problem for numerical treatment.The accuracy and stability of the solution is tested by calculating groundstate binding energies for few-boson systems interacting via the Ali-Bodmer, Volkov and Malfliet-Tjon (MTV) nuclear potentials for A-body systems (where A = 4, 5, 6, 7).A rapid convergence of the results is achieved and calculated ground-state energies are in good agreement with those available in the literature, where different approaches were used.The method is promising and could be extended to larger number of particles such as those involved in the Bose-Einstein condensates.
We present a semiclassical study of the spectrum of a few-body system consisting of two short-range interacting bosonic particles in one dimension, a particular case of a general class of integrable many-body systems where the energy spectrum is given by the solution of algebraic transcendental equations. By an exact mapping between δ-potentials and boundary conditions on the few-body wave functions, we are able to extend previous semiclassical results for single-particle systems with mixed boundary conditions to the two-body problem. The semiclassical approach allows us to derive explicit analytical results for the smooth part of the two-body density of states that are in excellent agreement with numerical calculations. It further enables us to include the effect of bound states in the attractive case. Remarkably, for the particular case of two particles in one dimension, the discrete energy levels obtained through a requantization condition of the smooth density of states are essentially in perfect agreement with the exact ones.
Universal low-energy properties are studied for three identical bosons confined in two dimensions. The short-range pairwise interaction in the low-energy limit is described by means of the boundary condition model. The wave function is expanded in a set of eigenfunctions on the hypersphere, and the system of hyperradial equations is used to obtain analytical and numerical results. Within the framework of this method, exact analytical expressions are derived for the eigenpotentials and the coupling terms of hyperradial equations. The derivation of the coupling terms is generally applicable to a variety of three-body problems provided the interaction is described by the boundary condition model. The asymptotic form of the total wave function at a small and a large hyperradius $\ensuremath{\rho}$ is studied, and the universal logarithmic dependence $\ensuremath{\sim}{\mathrm{ln}}^{3}\phantom{\rule{0.2em}{0ex}}\ensuremath{\rho}$ in the vicinity of the triple-collision point is derived. Precise three-body binding energies and the
(2+1)-scattering length are calculated.
We use the hyperspherical adiabatic expansion method to formulate the angular Faddeev-like equation for N identical bosons including effects of two-body correlations. Eigenvalues and eigenfunctions determine the effective potentials in the coupled set of radial equations. The zero-range approximation is implemented on the free solutions as a boundary condition obtained by using the effective-range expansion of the two-body phase shift. This results in a new angular eigenvalue equation expressed as a transcendental equation in the index of the Jacobi functions corresponding to solutions for non-interacting particles. The collapse for attractive potentials is avoided by including scattering length, effective range and shape-parameter terms in the effective-range expansion of the phase shift. Correspondingly the hyperradial dependence in the eigenvalue equation is either a polynomial of first, second, or fourth order. Two-body correlations are included. The arithmetic eigenvalue equation allows derivation of analytic properties such as for example simple scaling relations for the particle number dependence. Excellent agreement with finite-range numerical results is obtained. The couplings between adiabatic channels are discussed. Including only the lowest adiabatic potential is shown to be sufficient for most applications on dilute systems.
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.