The wave functions of boson and fermion gases are known even when the particles have harmonic interactions. Here we generalize these results by solving exactly the N–body Schrödinger equation for potentials V that can be any function of the sum of the squares of the distances of the particles from one another in three dimensions. For the harmonic case that function is linear in r2. Explicit N–body solutions are given when U(r) = –2Mh–2V(r) = 𝛇r–1 – 𝛇2r–2. Here M is the sum of the masses and r2 = 1/2 M–2ΣmImJ(xI – xJ)2. For general U(r) the solution is given in terms of the one– or two–body problem with potential U(r) in three dimensions. The degeneracies of the levels are derived for distinguishable particles, for bosons of spin zero and for spin 1/2 fermions. The latter involve significant combinatorial analysis which may have application to the shell model of atomic nuclei. For large N the fermionic ground state gives the binding energy of a degenerate white dwarf star treated as a giant atom with an N–body wave function. The N–body forces involved in these extraordinary N–body problems are not the usual sums of two–body interactions, but nor are forces between quarks or molecules. Bose–Einstein condensation of particles in three dimensions interacting via these strange potentials can be treated by this method.
The wave functions of boson and fermion gases are known even when the particles have harmonic interactions. Here we generalize these results by solving exactly the N–body Schrödinger equation for potentials V that can be any function of the sum of the squares of the distances of the particles from one another in three dimensions. For the harmonic case that function is linear in r2. Explicit N–body solutions are given when U(r) = –2Mh–2V(r) = 𝛇r–1 – 𝛇2r–2. Here M is the sum of the masses and r2 = 1/2 M–2ΣmImJ(xI – xJ)2. For general U(r) the solution is given in terms of the one– or two–body problem with potential U(r) in three dimensions. The degeneracies of the levels are derived for distinguishable particles, for bosons of spin zero and for spin 1/2 fermions. The latter involve significant combinatorial analysis which may have application to the shell model of atomic nuclei. For large N the fermionic ground state gives the binding energy of a degenerate white dwarf star treated as a giant atom with an N–body wave function. The N–body forces involved in these extraordinary N–body problems are not the usual sums of two–body interactions, but nor are forces between quarks or molecules. Bose–Einstein condensation of particles in three dimensions interacting via these strange potentials can be treated by this method.
We consider an analytic way to make the interacting N -body problem tractable by using harmonic oscillators in place of the relevant two-body interactions. The two-body terms of the N -body Hamiltonian are approximated by considering the energy spectrum and radius of the relevant two-body problem which gives frequency, centre position, and zero point energy of the corresponding harmonic oscillator. Adding external harmonic one-body terms, we proceed to solve the full quantum mechanical N -body problem analytically for arbitrary masses. Energy eigenvalues, eigenmodes, and correlation functions like density matrices can then be computed analytically. As a first application of our formalism, we consider the N -boson problem in two and three dimensions where we fit the two-body interactions to agree with the well-known zero-range model for two particles in a harmonic trap. Subsequently, condensate fractions, spectra, radii, and eigenmodes are discussed as a function of dimension, boson number N , and scattering length obtained in the zero-range model. We find that energies, radii, and condensate fraction increase with scattering length as well as boson number, while radii decrease with increasing boson number. Our formalism is completely general and can also be applied to fermions, Bose–Fermi mixtures, and to more exotic geometries.
We extend a method previously presented for the three-body problem with attractive interparticle potentials to the case of repulsive interparticle potentials, with periodic boundary conditions on each particle. As before, we decompose the wave function into three parts, and from the Schrödinger equation write equations for these parts. We observe that, in one-dimension and for δ-function potentials,these equations can be easily solved numerically, and we present these solutions for potential“strengths” ranging from zero to infinity. We then discuss a more general method for solving these equations, which method involves expansions in a certain set of two-body functions. This general method is not necessarily limited to either one dimension or δ-function potentials, but, as a check on it, we do apply it to that case and get good agreement with the previous numerical results. As a further exploration of the method, we apply it to square well potentials in one dimension. In an Appendix we discuss the set of two-body functions that we use.
The quantum-mechanical problems of N 1-dimensional equal particles of mass m interacting pairwise via quadratic (``harmonical'') and/or inversely quadratic (``centrifugal'') potentials is solved. In the first case, characterized by the pair potential ¼mω2(xi − xj)2 + g(xi − xj)−2, g > −ℏ2/(4m), the complete energy spectrum (in the center-of-mass frame) is given by the formula E=ℏω(12N)12[12(N−1)+12N(N−1)(a+12)+ ∑ l=2Nlnl],with a = ½(1 + 4mgℏ−2)½. The N − 1 quantum numbers nl are nonnegative integers; each set {nl; l = 2, 3, ⋯, N} characterizes uniquely one eigenstate. This energy spectrum can also be written in the form Es = ℏω(½N)½ [½(N − 1) + ½N(N − 1)(a + ½) + s], s = 0, 2, 3, 4, ⋯, the multiplicity of the sth level being then given by the number of different sets of N − 1 nonnegative integers nl that are consistent with the condition s=∑l=2Nlnl. These equations are valid independently of the statistics that the particles satisfy, if g ≠ 0; for g = 0, the equations remain valid with a = ½ for Fermi statistics, a = −½ for Bose statistics. The eigenfunctions corresponding to these energy levels are not obtained explicitly, but they are rather fully characterized. A more general model is similarly solved, in which the N particles are divided in families, with the same quadratic interaction acting between all pairs, but with the inversely quadratic interaction acting only between particles belonging to the same family, with a strength that may be different for different families. The second model, characterized by the pair potential g(xi − xj)−2, g > −ℏ2/(4m), contains only scattering states. It is proved that an initial scattering configuration, characterized (in the phase space sector defined by the inequalities xi ≥ xi.1, i = 1, 2, ⋯, N = 1, to which attention may be restricted without loss of generality) by (initial) momenta pi, i = 1, 2, ⋯, N, goes over into a final configuration characterized uniquely by the (final) momenta pi′, with pi′=pN+1−i. This remarkably simple outcome is a peculiarity of the case with equal particles (i.e., equal masses and equal strengths of all pair potentials).
The resources required to solve the general interacting quantum N-body problem scale exponentially with N, making the solution of this problem very difficult when N is large. In a previous series of papers we develop an approach for a fully interacting wave function for a confined system of identical bosons with a general two-body interaction. This method tames the N scaling by developing a perturbation series that is order-by-order invariant under a point group isomorphic with SN. Group theory and graphical techniques are then used to solve for the wave function exactly and analytically at each order, yielding a solution for the general N-body problem which scales as N0 at any given order. Recently this formalism has been used to obtain the first-order fully interacting wave function for a system of harmonically confined bosons interacting harmonically. In this paper, we report the application of this N-body wave function to a system of N fully interacting bosons in three dimensions. We derive an expression for the density profile for a confined system of harmonically interacting bosons. Choosing this simple interaction is not necessary or even advantageous for our method, however this choice allows a direct comparison of our exact results through first order with exact results obtained in an independent solution. Our density profile for the wave function through first order in three dimensions is indistinguishable from the first-order exact result obtained independently and shows strong convergence to the exact result to all orders.
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.
Exact solutions to many-body interacting systems of both bosonic and fermionic particles confined to harmonic potential in an arbitrary number of dimensions are given. Energy levels and their degeneracies for trapped identical particles interacting via harmonic forces are calculated. This specific form of the interaction allows for analytical solutions. The mutual interaction, attractive or repulsive, modifies significantly the properties of the considered system. For a large number of particles the interaction essentially results in a frequency shift. Statistical properties (e.g., microcanonical and grand canonical partition functions) as well as some illustrative, physically relevant examples are discussed. Our results give an unusual opportunity for further studies of interacting systems in the framework of the exactly soluble model.