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.
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.
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 simple model of a boson system is considered for which the exact Schr\"odinger equation reduces to a difference equation which can be solved numerically. Application of the Bogoliubov method and comparison of its results with the numerical results then yield some insight into the validity of the Bogoliubov approximations. The simple model considered has, in the absence of interactions, only a zero energy state and two states of unit energy available to each boson. Initially there is assumed to be no interaction between bosons in the degenerate excited states and all existing interactions are assumed to be repulsive. The ground-state energy is calculated in the Bogoliubov approximation, with and without the depletion effect considered, and in another numerical approximation. These results and the results of the numerical solution of the exact Schr\"odinger equation converge to the Bogoliubov result with depletion ignored when the number of particles in the system approaches infinity. It is surprising to note that, for the interaction strengths considered, the Bogoliubov result is within 3% of the numerical result for as few as 32 particles in the system, and within 10% for only 4 particles. A modification of the above system is considered in which there exists an additional two-body interaction between particles in the degenerate excited states which may be attractive or repulsive. It is shown that the ground-state energy, with this additional interaction present, deviates from the Bogoliubov value linearly with the strength of the added interaction, i.e., $W={W}_{\mathrm{Bog}}+\ensuremath{\alpha}F$ (with F the strength parameter of the interaction). The factor $\ensuremath{\alpha}$ is found to be small and constant over a wide range of F and in this range the deviation of W from the Bogoliubov value is small. There is, however, a transition point FT, beyond which $\ensuremath{\alpha}$ is constant and large and the usual Bogoliubov approximation is invalid for interactions more attractive than that characterized by FT. A Bogoliubov-like approximation is shown to be quite accurate in this region. It is then shown that $F_{T}^{\ensuremath{-}1}$ is a linear function of the number of particles in the system.
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 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 Bethe roots describing the ground-state energy of the integrable one-dimensional (1D) model of interacting bosons with weakly repulsive two-body delta interactions are seen to satisfy the set of Richardson equations appearing in the strong coupling limit of an integrable BCS pairing model. The BCS model describes boson–boson interactions with zero centre of mass momentum of pairs. It follows that the Bethe roots of the weakly interacting boson model are given by the zeros of Laguerre polynomials. The ground-state energy and the lowest excitation are obtained explicitly via the Bethe roots. A direct link has thus been established, in the context of integrable 1D models, between bosons interacting via weakly repulsive two-body delta interactions and strongly interacting BCS boson pairs.
A general system composed of many weakly interacting bosons and/or fermions is considered. The object is to develop a procedure for renormalizing the single-particle creation operators, so as to remove the interactions to successively higher orders of perturbation. The basic idea is that creation operators ${{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$ must satisfy the Hamiltonian commutator equations $[\mathcal{H},{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}]={\ensuremath{\omega}}_{i}{{\ensuremath{\theta}}_{i}}^{\ifmmode\dagger\else\textdagger\fi{}}$, where H is the Hamiltonian and ${\ensuremath{\omega}}_{i}$ are the particle energies. For the many-weakly-interacting-particle system, the zeroth-order boson and fermion creation operators satisfy these equations to zeroth order. It is shown that if the creation operators are renormalized so as to satisfy the Hamiltonian commutator equations to order m, and also to satisfy the appropriate boson commutators and fermion anticommutators to order m, then the problem is solved to order m. In particular, the vectors formed by operating on the ground state with the renormalized creation operators, according to the usual boson and fermion occupation-number representation, are eigenvectors of H and are orthonormal, all to order m. A procedure is otained for finding the (m+1)-order contributions to the particle creation energies, in terms of the m-order operators. Explicit first-order calculations of these general results are provided for a system of bosons and a system of fermions, and these first-order results are shown to include similar results of Rayleigh-Schr\"odinger perturbation theory, the random-phase approximation, and the method of thermodynamic Green's functions. The problem of anharmonic lattice dynamics is studied in detail, and a method of undetermined coefficients is used to renormalize the phonon creation operators to first order. The phonon energies are calculated to second order, and this calculation shows that the interactions between renormalized phonons cannot be removed in second order. Statistical averages of the phonon energies give the energy shifts and lifetimes which have been calculated previously by various propagator techniques. In addition, the renormalized energies are used to calculate the temperature-dependent part of the Helmholtz free energy correct to second order. As a final example, electron-phonon interactions in a normal metal are studied. The electron and phonon creation operators are renormalized to first order, the particle energies are calculated to second order, and statistical averages of the particle energies recover the usual thermodynamic Green's-function results for energy shifts and lifetimes. These examples show the simplicity by which the renormalization procedure obtains a great amount of detailed information about the single-particle nature of a many-particle system.