Calculations of the trimer ground-state energies of three identical bosons interacting via molecular pair potentials are discussed. The accuracies of variational calculations with a Jastrow trial function and of upper and lower bounds in terms of dimer ground-state energies are determined in a harmonic approximation.
Variational calculations of the ground state energy of 4He trimers are reported and are compared with the Hall-Post-Stenschke lower bound on the ground state energy. The variational result for the ground state energy for a range of model pair potentials is in the range −0.05 to −0.2 K; the effect of the triplet nonadditive term is estimated to be less than 1% of the contribution of the pair potentials. A study of the variation of the ground state energy of 3 bosons with pair potential coupling constant is also reported for a Lennard-Jones and a Morse potential model. There is no bound spin 3/2 3He trimer; a lower bound is given for the ground state energy of the spin 1/2 trimer and analogy with the boson results is used to argue that the spin 1/2 trimer is probably not bound. This paper contains computational applications of formal results derived elsewhere.
The 4 He trimer is treated by the Ritz variational method using trial functions whose accuracy has been tested by a similar treatment of the dimer, for which an exact solution exists. These are nodeless functions containing four nonlinear parameters that are varied to give the minimum value for the ground-state energy. Agreement with measurement is obtained for the mean bond lengths of the trimer, but there is some disagreement with other non-variational calculations of the ground-state binding energy.
We calculated the 4He trimer and tetramer ground and excited states with the LM2M2 potential using our Gaussian expansion method for ab initio variational calculations of few-body systems. The method has been extensively used for a variety of three-, four-, and five-body systems in nuclear physics and exotic atomic and molecular physics. The trimer (tetramer) wave function is expanded in terms of symmetric three- (four-) body Gaussian basis functions, ranging from very compact to very diffuse, without assumption of any pair correlation function. The calculated results for the trimer ground and excited states are in excellent agreement with values reported in the literature. The binding energies of the tetramer ground and excited states are obtained as 558.98 and 127.33 mK (0.93 mK below the trimer ground state), respectively. We found that precisely the same shape of the short-range correlation (${r}_{ij}\ensuremath{\lesssim}4$ \AA{}) in the dimer appears in the ground and excited states of the trimer and tetramer. The overlap function between the trimer excited state and the dimer ground state and that between the tetramer excited state and the trimer ground state are almost proportional to the dimer wave function in the asymptotic region (up to $\ensuremath{\sim}$1000 \AA{}). Also, the pair correlation functions of trimer and tetramer excited states are almost proportional to the squared dimer wave function. We then propose a model which predicts the binding energy of the first excited state of 4HeN ($N\ensuremath{\ge}3$) measured from the 4He${}_{N\ensuremath{-}1}$ ground state to be nearly $\frac{N}{2(N\ensuremath{-}1)}{B}_{2}$ where B2 is the dimer binding energy.
A trial wave function describing the ground state of a quantum system of N interacting bosons is written in the Jastrow form, a product of pair functions. With the interaction potential chosen to represent liquid He4, and with the parametrized form of the pair function chosen to include a long-range term which has been found necessary to represent the zero-point motion of the long-wavelength density oscillations, a variational calculation has been performed using a new approximate integral equation for the pair distribution function. This equation, which can also be used for classical fluids, is found to be more accurate for repulsive potentials than the Percus-Yevick equation and comparable to (but much simpler than) the Percus-Yevick 2 equation. The essential results are that including the zero-point motion in the wave function tends to lower the energy, raises the equilibrium density, corrects the behavior of the structure function and the momentum distribution of the particles in the low-wave-number region, and slightly decreases the Bose-Einstein condensate fraction. The value of the lower limit on the wavelength of the density oscillations was determined variationally to be about three interparticle spacings.
The method of correlated basis functions is used in the study of the ground-state properties of a one-dimensional system of bosons interacting through a repulsive $\ensuremath{\delta}$-function potential. Ground-state energies are determined numerically for intermediate values of the coupling constant in three steps: (1) a trial form of the radial distribution function generated by the Bijl-Dingle-Jastrow (BDJ) type of wave function is introduced to obtain approximate ground-state energies, (2) the paired-phonon analysis is applied to improve the BDJ description of the ground state, and (3) leading corrections to the improved approximate energies are evaluted from the second-order perturbation energy generated by the three-phonon vertex. The obtained energy values are found to agree closely with the exact results calculated by Lieb and Liniger.