In a preceding paper [J. Math. Phys. 3, 475 (1962)], some model Hamiltonians were proposed for quantum-mechanical many-body systems with pair forces. For infinite systems in thermal equilibrium, they led to temperature-domain propagator expansions which were formally summable and expressible by closed equations. These expansions were identical with infinite subclasses of terms from the propagator expansion for the true many-body problem. The two principal models corresponded to ring- and ladder-diagram summations from the true propagator expansion, augmented by infinite classes of self-energy corrections. The model Hamiltonians were called stochastic because they contained parameters whose phases were fixed by random choices. In the present paper, more general models are formulated which yield formally summable propagator expansions for finite systems. The analysis is extended to correlation and Green's functions defined for nonequilibrium ensembles. The nonequilibrium treatment is developed in the Heisenberg representation in such a way that unlinked diagrams do not arise. A basic convergence question associated with the formal closed equations for the model propagators and correlation functions is examined by means of finite-difference integration of the Heisenberg equations of motion. This procedure appears to converge independently of whether the perturbation expansions for the propagators and correlation functions converge. It yields substantial support for the validity of the formal closed model equations.
Identical Particles Second Quantization Independent-Particle Model for Fermions in Finite Systems Two-Particle States and Interactions Noninteracting Bosons and Fermions Propagators in One-Particle Quantum Mechanics Single-Particle Propagator in Many-Particle Systems Perturbation Expansion of the Single-Particle Propagator Dyson Equation and Self-consistent Green's Functions Mean-Field or Hartree-Fock Approximation Beyond the Mean-Field Approximation Interacting Boson Systems Excited States in Finite Systems Excited States in Infinite Systems Excited States in N 2 Systems and In-Medium Scattering Dynamical Treatment of the Self-energy in Infinite Systems Dynamical Treatment of the Self-energy in Finite Systems Bogoliubov Perturbation Expansion for the Bose Gas Boson Perturbation Theory Applied to Physical Systems In-Medium Interaction and Scattering of Dressed Particles Conserving Approximations and Excited States Pairing Phenomena Propagator Framework for Phenomenological Description of Many-Body Systems Finite Temperature Formalism and Applications.
A method for treating nonlinear stochastic systems is described which it is hoped will be useful in both the quantum-mechanical many-body problem and the theory of turbulence. In this method the true problem is replaced by models that lead to closed equations for correlation functions and averaged Green's functions. The model solutions are exact descriptions of possible dynamical systems, and, as a result, they display certain consistency properties. For example, spectral components of Green's functions which must be positive-definite in the true problem automatically are so for the models. The models involve a new stochastic element: Random couplings are introduced among an infinite collection of similar systems, the true problem corresponding to the limit where these couplings vanish. The method is first applied to a linear oscillator with random frequency parameter. The mean impulse-response function of the oscillator is obtained explicitly for two successive models. The results suggest the existence of a sequence of model solutions which converges rapidly to the exact solution of the true problem. Applications then are made to the Schrödinger equation of a particle in a random potential and to Burgers' analog for turbulence dynamics. For both problems, closed model equations are obtained which determine the average Green's function, the amplitude of the mean field, and the covariance of the fluctuating field. The model solutions can be expressed as sums of infinite classes of terms from the formal perturbation expansions of the solutions to the true problems. It is suggested that correspondence to stochastic models may be a useful criterion to help judge the validity of partial summations of perturbation series.
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.
A self consistent field approach is developed for truncating the infinite set of equations obtained in apply ing the equation of motion method to interacting systems. The operator which is assumed to create elementary excita tions in the system is expanded in terms of simpler opera tors and both the coefficients in the expansion and relations independent of the coefficients are found. The method is applied to a system of Bosons at low temperatures in the approximation that the interaction is both tempera ture and density independent. Corrections to the pairing approximation are obtained due both to the inclusion of more terms in the Hamiltonian and the inclusion of collective excitations. It is shown that in the long wavelength limit these corrections can be included by renormalizing both the single particle and pairing energies. The renormalization can be carried out by adding a term to the single particle energy and an extra interaction to the pairing part of the Hamiltonian. This interaction is attractive over part of momentum space and obeys the criterion for a simple Bose condensation. It is shown that no single particle like zero momentum excitation exists. A collective excitation is con sidered for the zero momentum excitation and it is found v vi that its energy must be non-zero if it exists at all. The application to superfluid He11 is discussed.
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.