The many-boson system with repulsive interactions is treated by a variational method based on a variational trial state of an exponential pair-excitation type obtained by a generalization of that of Bogoliubov; the treatment is closely related to an intermediate-coupling approximation with respect to pairs. The nonlinear integral equation which determines the variational ground state is derived, and various properties of this ground state and the associated excited states are examined. The resultant low-lying spectrum lies below that of Bogoliubov by an amount proportional to the total number of particles. The variational principle is shown to produce rigorous energy eigenvalue differences for the pair part of the Hamiltonian. The variational states, however, still exhibit unphysical features characteristic of pair-excitation states: The pair correlation function does not go to zero at zero particle separation and the phonon spectrum exhibits a gap above the ground state. It is suggested that these features can be removed by using states which take into account excitation of momentum-conserving groups of more than two particles.
A many-boson model is formulated and expressions for its exact eigenstates and energies are obtained for both an arbitrary finite and an infinite number of bosons. The Hamiltonian of the model contains interactions between bosons whose momenta have equal magnitudes but opposite directions. The matrix elements of this interaction are taken to be a constant over a range of momenta surrounding k = 0. The ground state of the 2N-particle system is shown to be a product of N pair-creation operators acting on the vacuum state. Each of these pair-creation operators depends upon one of N parameters which are called pair energies. The N pair energies are shown to satisfy a coupled system of nonlinear algebraic equations. The energy of the state is the sum of the pair energies and the occupation probabilities of the single-particle levels are given as simple functions of the pair energies. Similar results are derived for the excited states of the system and for the states of an odd number of particles. These results are valid for both a repulsive and an attractive interaction, since they only depend upon the form of the interaction. The equations are solved algebraically for two model systems. The first of these is one whose single-particle kinetic energy takes on only one value. The equations for this system are solved for an arbitrary interaction strength and it is shown that the pair energies are proportional to the zeros of certain Laguerre polynomials. The second system is one in which the single-particle kinetic energy can take on two values. The equations for this system are solved in the strong repulsive-interaction limit and it is shown that the pair energies are proportional to the zeros of certain Jacobi polynomials. The excitation energies of this second system are shown to be proportional to 1/n and the occupations of the two single-particle levels in the ground state are shown to be proportional to n, where n is the total number of particles. For a repulsive interaction and an arbitrary single-particle spectrum, the algebraic equations for the pair energies are converted into an approximate integral equation for the density of roots which is accurate to order 1/n. This integral equation is solved for a strong interaction which, in the context of this model, means an interaction whose strength is greater than a constant times 1/V⅔ in the limit of a large volume. From this solution, the following results are obtained: (1) the lowest two single-particle levels have occupations of order n; (2) the excitation spectrum is that of a set of noninteracting quasiparticles; and (3) the quasiparticle spectrum has two zeros corresponding to the lowest two single-particle levels. Apart from the presence of two zeros, the quasiparticle spectrum does not differ significantly from that of the noninteracting particles.
An independent-pair ansatz is developed for the many-body wave function of dilute Bose systems. The pair correlation is optimized by minimizing the expectation value of the full Hamiltonian (rather than the truncated Bogoliubov one) providing a rigorous energy upper bound. In contrast with the Jastrow model, hypernetted chain theory provides closed-form exactly solvable equations for the optimized pair correlation. The model involves both condensate and coherent pairing with number conservation and kinetic-energy sum rules satisfied exactly and the compressibility sum rule obeyed at low density. We compute, for bulk boson matter at a given density and zero temperature, (i) the two-body distribution function, (ii) the energy per particle, (iii) the sound velocity, (iv) the chemical potential, (v) the momentum distribution and its condensate fraction, and (vi) the pairing function, which quantifies the off-diagonal long-range order resulting from the structural properties of the two-particle density matrix. The connections with the low-density expansion and Bogoliubov theory are analyzed at different density values, including the density and scattering length regime of interest of trapped-atom Bose-Einstein condensates. Comparison with the available diffusion Monte Carlo results is also made.
This work is concerned with a microscopic model of a nonideal boson gas. The Hamiltonian is treated in a rigorous manner free of all approximations and, in particular, the Bogoliubov approximation of replacing the zero-momentum single-particle creation and destruction operators by a
c number is entirely avoided. Utilizing a method due to Wentzel, the Hamiltonian is replaced by one involving products of two second quantized operators which, in general, can be diagonalized. This replacement is shown to be rigorously valid in the infinite volume limit for all temperatures despite the avoidance of the Bogoliubov approximation. It is found that the assumption of a positive-definite quasiparticle excitation spectrum $\ensuremath{\epsilon}(\mathrm{k})$ becomes untenable for temperatures
T less than a critical value
Tc, if the interparticle potential is described by a certain kernel which possesses positive eigenvalues only, corresponding to a sufficiently repulsive interparticle force. Assuming that $\ensuremath{\epsilon}(\mathrm{k})=0$ only for k=0, it is found that for
T<Tc the system undergoes an Einstein condensation into the k=0 single-particle state, and the integral equations characterizing the system are precisely the same as would be obtained if the Bogoliubov approximation were made, proving within the framework of this model the strict validity of this approximation procedure. The well-known criterion for an Einstein condensation
v(0)>0, where
v(k) is the Fourier transform of the interparticle potential, is shown to be a weaker statement of our eigenvalue criterion. Further, in agreement with Girardeau and Arnowitt, it is found that an energy gap separates $\ensuremath{\epsilon}(0)$ and $\ensuremath{\epsilon}(\mathrm{k})$ for
k>0. An approximation method is developed for solving the integral equations which describe those systems which undergo an Einstein condensation, and it is shown that this method can be justified if $|T\ensuremath{-}{T}_{c}|\ensuremath{\lesssim}\frac{{T}_{c}}{10}$, and in cases of short-range repulsive interparticle forces and systems for which $\ensuremath{\rho}v(0)$ possesses a certain prescribed upper bound, $\ensuremath{\rho}$ being the number of particles per unit volume. The transition temperature is found to be lower than the corresponding quantity for the ideal boson gas of the same density. A detailed discussion of the thermal properties of the system in the vicinity of
Tc is also given. In an Appendix a general type of smeared Einstein condensation is assumed, and it is found that for systems in which this property occurs $\ensuremath{\epsilon}(\mathrm{k})$ is linear for small
k. This ansatz is shown to be tenable only if
v(0)<0. Finally, by assuming an appropriate pseudopotential representation of a hard repulsive core and by including a weak, longranged attractive force between particles, $\ensuremath{\epsilon}(\mathrm{k})$ shows a nonmonotonic behavior of the same qualitative type as observed in liquid
He4.
The model recently proposed in Ref.1 is used to derive linear integro– differential equations whose solutions provide reasonable estimates for the momentum distribution and condensate fraction in interacting many–boson system at zero temperature. An advantage of these equations is that they can be employed in the weak coupling regime and beyond. As an example, analytical treatment of the weak coupling case is given. PACS: 05.30.Jp, 67.90.+z A new way of investigating spatial particle correlations has recently been proposed for the many–boson system in the ground state [1]. It results in integro–differential equations for the pair distribution function. They accurately take into account the short–range boson correlations and, in the weak coupling regime, yield thermodynamics which reasonably agrees
Creation and annihilation operators of quasi-particles are introduced into the description of a finite temperature system of degenerate bosons with the aid of a canonical transformation which generalizes the one due to Bogoliubov. A variational calculation taking into account the interaction between non-zero momentum particles shows that the quasi-particles are phonons in the long wave-length limit. The state of the system, which is approximated by the grand canonical ensemble of the quasi-particles, approaches the Hartree-Fock state when the condensate is depleted completely at the critical temperature.