Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

Canadian Journal of Physics

Theory of Brueckner and Sawada Applied to a Many-Boson Model

Mohsen Razavy, Edward S. Krebes

Published Feb 1, 1973• 4 citations
Publisher

Abstract

The Bassichis–Foldy model of a simple interacting boson is solved numerically and the results are compared with those obtained by the Bogoliubov approximation and by the Brueckner–Sawada t-matrix formalism. In the normal region, contrary to the widely held view, the Brueckner–Sawada approximation for the energy of the ground state is not reliable for strong, well-behaved, repulsive forces. The Bogoliubov approximation, on the other hand, remains valid for a wide range of values of the coupling constant. In the inverted region, the attractive force causes a population inversion in the levels of the system. For this case a modified Brueckner–Sawada approximation is developed. This method is applied to the calculation of the transition point and the energies of the ground and the first excited states of the system. Here most of the predictions of the modified Brueckner–Sawada approximation are quite accurate. By a simple change in the Bassichis–Foldy model it is shown that even, for two bosons there can be a phase transition. In this model, the derivative of the ground state energy with respect to the coupling constant is discontinuous at the transition point.

DOI
10.1139/p73-037

Related papers

Analysis of a Boson Model Displaying a Phase Transition

William H. Bassichis, Y. Disantnik

The simple boson model proposed by Foldy, which displayed a discontinuity in its energy as a function of interaction strength, is studied using symmetry arguments, the method of linearized equations of motion, and perturbation theory. The particular symmetry responsible for the discontinuity, or phase transition, is pointed out. The method of linearized equations is shown to predict the transition value of the interaction strength and, like the Bogoliubov approximation, to accurately describe the energy of the first excited state. It is further shown that simple perturbation theory is quite accurate except in the transition region. The validity of that feature of the Bogoliubov method in which excited states are described by independent quasiparticles is examined. This description is found to be excellent in the normal region. All comparisons are made with a numerical solution which provided the complete energy spectrum and the associated eigenvectors. It is hoped that these arguments can be applied in considerations of the validity of the various approximations for realistic systems.

Physical Review•Jan 7, 1966•1 citation

Analysis of the Bogoliubov Method Applied to a Simple Boson Model

William H. Bassichis, Leslie L. Foldy

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 ${F}_{T}$, beyond which $\ensuremath{\alpha}$ is constant and large and the usual Bogoliubov approximation is invalid for interactions more attractive than that characterized by ${F}_{T}$. 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.

Physical Review•Feb 17, 1964•41 citations

Effects of the T-matrix center-of-mass approximation in the Brueckner–Sawada theory of liquid helium II

S. K. Bose, D. F. Goble

As a model for liquid helium II, we study the Bose–Einstein gas with a two-body interaction potential of the form ~δ(r − a), where r is the interparticle separation. Excitation spectra for various values of a are calculated using the Brueckner–Sawada approach based on the concept of the T matrix. However, unlike the work of Brueckner and Sawada and many other related works that followed, we take into account the dependence of the T matrix on the center-of-mass momentum of the interacting particles. Excitation spectra calculated with and without this dependence, using two different expressions for the two-particle propagator, indicate the validity of the Brueckner and Sawada center-of-mass approximation for physically interesting values of a. The difference between the two spectra is found to increase as a increases.

Canadian Journal of Physics•Dec 1, 1985•0 citations

Excitation Spectrum of the Bose Liquid

E. Byckling

The assumptions and predictions of the Brueckner-Sawada method (including recent refinements) for the derivation of the excitation spectrum of the high-density boson system are examined using Green's-function techniques and a new method of solving the scattering-matrix equation. In the case of the real interaction potential of ${\mathrm{He}}^{4}$ atoms, the spectrum obtained in our approximation has the correct form but the depletion has a meaningless value. It is pointed out that a probable cause of discrepancies is inconsistent omission of certain self-energy terms.

Physical Review•May 6, 1966•14 citations

Exactly Solvable Many-Boson Model

R. W. Richardson

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.

Journal of Mathematical Physics•Sep 1, 1968•82 citations

Self-consistent-field approximation method for a one-dimensional many-boson system

Oscar Hipólito, R. Lobo

The method of self-consistent-field approximation is applied for the study of the ground-state properties of a one-dimensional system of bosons, of density $\ensuremath{\rho}$, first interacting through finite rigid rods, of length $l$, and again via a repulsive $\ensuremath{\delta}$-function potential. We obtain the energy spectrum and the sound velocity. The results are close to the exact one for $\ensuremath{\rho}l<0.5$. The energy spectrum obtained shows a rotonlike region which is qualitatively similar to the two- and three-dimensional versions. For a Bose gas interacting through a repulsive $\ensuremath{\epsilon}$ potential of strength $2c$, we derive the energy spectrum, the sound velocity, and the ground-state energy. The results are compared with Girardeau and Lieb's exact expressions. In the small-coupling limit our calculations show the same results as obtained by Bogoliubov's perturbation theory. The ground-state energy obtained for small $\ensuremath{\gamma}=\frac{c}{\ensuremath{\rho}}$ is found to agree closely with several different works found in the literature. Unfortunately our theory diverges for the ground-state energy in the limit $\ensuremath{\gamma}\ensuremath{\rightarrow}\ensuremath{\infty}$. The sound velocity in this limit is finite but twice as small as from Girardeau's exact results.

Physical review. B, Condensed matter•Feb 1, 1978•1 citation
The Bassichis–Foldy model of a simple interacting boson is solved numerically and the results are compared with those obtained by the Bogoliubov approximation and by the Brueckner–Sawada t-matrix formalism. In the normal region, contrary to the widely held view, the Brueckner–Sawada approximation for the energy of the ground state is not reliable for strong, well-behaved, repulsive forces. The Bogoliubov approximation, on the other hand, remains valid for a wide range of values of the coupling constant. In the inverted region, the attractive force causes a population inversion in the levels of the system. For this case a modified Brueckner–Sawada approximation is developed. This method is applied to the calculation of the transition point and the energies of the ground and the first excited states of the system. Here most of the predictions of the modified Brueckner–Sawada approximation are quite accurate. By a simple change in the Bassichis–Foldy model it is shown that even, for two bosons there can be a phase transition. In this model, the derivative of the ground state energy with respect to the coupling constant is discontinuous at the transition point.
The simple boson model proposed by Foldy, which displayed a discontinuity in its energy as a function of interaction strength, is studied using symmetry arguments, the method of linearized equations of motion, and perturbation theory. The particular symmetry responsible for the discontinuity, or phase transition, is pointed out. The method of linearized equations is shown to predict the transition value of the interaction strength and, like the Bogoliubov approximation, to accurately describe the energy of the first excited state. It is further shown that simple perturbation theory is quite accurate except in the transition region. The validity of that feature of the Bogoliubov method in which excited states are described by independent quasiparticles is examined. This description is found to be excellent in the normal region. All comparisons are made with a numerical solution which provided the complete energy spectrum and the associated eigenvectors. It is hoped that these arguments can be applied in considerations of the validity of the various approximations for realistic systems.
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 FFF the strength parameter of the interaction). The factor $\ensuremath{\alpha}$ is found to be small and constant over a wide range of FFF and in this range the deviation of WWW from the Bogoliubov value is small. There is, however, a transition point FT{F}_{T}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{F}_{T}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.
As a model for liquid helium II, we study the Bose–Einstein gas with a two-body interaction potential of the form ~δ(r − a), where r is the interparticle separation. Excitation spectra for various values of a are calculated using the Brueckner–Sawada approach based on the concept of the T matrix. However, unlike the work of Brueckner and Sawada and many other related works that followed, we take into account the dependence of the T matrix on the center-of-mass momentum of the interacting particles. Excitation spectra calculated with and without this dependence, using two different expressions for the two-particle propagator, indicate the validity of the Brueckner and Sawada center-of-mass approximation for physically interesting values of a. The difference between the two spectra is found to increase as a increases.
The assumptions and predictions of the Brueckner-Sawada method (including recent refinements) for the derivation of the excitation spectrum of the high-density boson system are examined using Green's-function techniques and a new method of solving the scattering-matrix equation. In the case of the real interaction potential of He4{\mathrm{He}}^{4}He4 atoms, the spectrum obtained in our approximation has the correct form but the depletion has a meaningless value. It is pointed out that a probable cause of discrepancies is inconsistent omission of certain self-energy terms.
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 method of self-consistent-field approximation is applied for the study of the ground-state properties of a one-dimensional system of bosons, of density $\ensuremath{\rho}$, first interacting through finite rigid rods, of length lll, and again via a repulsive $\ensuremath{\delta}$-function potential. We obtain the energy spectrum and the sound velocity. The results are close to the exact one for $\ensuremath{\rho}l<0.5$. The energy spectrum obtained shows a rotonlike region which is qualitatively similar to the two- and three-dimensional versions. For a Bose gas interacting through a repulsive $\ensuremath{\epsilon}$ potential of strength 2c2c2c, we derive the energy spectrum, the sound velocity, and the ground-state energy. The results are compared with Girardeau and Lieb's exact expressions. In the small-coupling limit our calculations show the same results as obtained by Bogoliubov's perturbation theory. The ground-state energy obtained for small $\ensuremath{\gamma}=\frac{c}{\ensuremath{\rho}}$ is found to agree closely with several different works found in the literature. Unfortunately our theory diverges for the ground-state energy in the limit $\ensuremath{\gamma}\ensuremath{\rightarrow}\ensuremath{\infty}$. The sound velocity in this limit is finite but twice as small as from Girardeau's exact results.