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

Physical Review

Analysis of a Boson Model Displaying a Phase Transition

William H. Bassichis, Y. Disantnik

Published Jan 7, 1966• 1 citation
Publisher

Abstract

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.

DOI
10.1103/physrev.141.217

Related papers

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

Mohsen Razavy, Edward S. Krebes

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.

Canadian Journal of Physics•Feb 1, 1973•4 citations

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

Chemical Potential and Low-Lying Excitations in the Many-Boson System. II

Ludwig W. Bruch, K. Sawada

The approximate description of a many-boson system by its density fluctuations is considered. The representation of the chemical potential by a scattering matrix is compared with an exact solution in one dimension. It is found to be valid exactly in the low coupling limit and valid to within 10% for intermediate coupling strengths. In the case treated, the ground-state properties are given quite accurately for intermediate coupling strengths. The excitation spectrum is described less accurately. The velocity of sound obtained from the slope of the phonon spectrum is formally different from that obtained using the compressibility. Here, the sound velocities obtained from the two methods are found to agree numerically. Data for actual ${\mathrm{He}}^{4}$ and a simple linear approximation to the density dependence of the scattering matrix are used. The sound velocity obtained is in rough agreement with the experimentally observed velocity in liquid ${\mathrm{He}}^{4}$.

Physical Review•Oct 15, 1963•1 citation

A Soluble Fermi-Gas Model. Validity of Transformations of the Bogoliubov Type

J. des Cloizeaux

A linear soluble Fermi-gas model is used to test the validity of perturbation theory and the criteria of instability obtained by applying transformations of the Bogoliubov type. The model Hamiltonian describing the interaction of the electrons in a simple band is H(ρ)=−Σ(cosk+ρ)ck+ck+ρΣcos q ck+q+cu−q+cuck. For a half-filled band, the energies of the low-lying states are calculated exactly, by using previous results concerning a linear chain of spins with interaction. It is shown that perturbation theory must be valid for −1 < ρ < 1, and, for values of ρ belonging to this range, the slope of the specific heat for T = 0 is calculated by direct application of Landau's theory. On the other hand, anomalous Hartree-Fock states can be built for all values of ρ, by using transformations of the Bogoliubov type, and these states are more stable than the normal Hartree-Fock state. Thus, it appears that the Bogoliubov method may lead to completely erroneous interpretations, when the coupling constants are small.

Journal of Mathematical Physics•Dec 1, 1966•45 citations

Exact versus mean-field description of the Bose-Einstein condensate: a model study

Magdalena A. Załuska–Kotur, Mariusz Gajda, Arkadiusz Orłowski +1 more

We study a system of trapped bosonic particles interacting by model harmonic forces. Our model allows for detailed examination of the notion of an order parameter (a condensate wave function). By decomposing a single particle density matrix into coherent eigenmodes we study an effect of interaction on the condensate. We show that sufficiently strong interactions cause that the condensate disappears even if the whole system is in its lowest energy state. In the second part of our paper we discuss the validity of the Bogoliubov approximation by comparing its predictions with results inferred from the exactl y soluble model. In particular we examine an energy spectrum, occupation, and fluctuations of the condensate. We conclude that Bogoliubov approach gives quite accurate description of the system in the limit of weak interactions.

arXiv:cond-mat/0004469•Apr 27, 2000•4 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
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.
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.
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.
The approximate description of a many-boson system by its density fluctuations is considered. The representation of the chemical potential by a scattering matrix is compared with an exact solution in one dimension. It is found to be valid exactly in the low coupling limit and valid to within 10% for intermediate coupling strengths. In the case treated, the ground-state properties are given quite accurately for intermediate coupling strengths. The excitation spectrum is described less accurately. The velocity of sound obtained from the slope of the phonon spectrum is formally different from that obtained using the compressibility. Here, the sound velocities obtained from the two methods are found to agree numerically. Data for actual He4{\mathrm{He}}^{4}He4 and a simple linear approximation to the density dependence of the scattering matrix are used. The sound velocity obtained is in rough agreement with the experimentally observed velocity in liquid He4{\mathrm{He}}^{4}He4.
A linear soluble Fermi-gas model is used to test the validity of perturbation theory and the criteria of instability obtained by applying transformations of the Bogoliubov type. The model Hamiltonian describing the interaction of the electrons in a simple band is H(ρ)=−Σ(cosk+ρ)ck+ck+ρΣcos q ck+q+cu−q+cuck. For a half-filled band, the energies of the low-lying states are calculated exactly, by using previous results concerning a linear chain of spins with interaction. It is shown that perturbation theory must be valid for −1 < ρ < 1, and, for values of ρ belonging to this range, the slope of the specific heat for T = 0 is calculated by direct application of Landau's theory. On the other hand, anomalous Hartree-Fock states can be built for all values of ρ, by using transformations of the Bogoliubov type, and these states are more stable than the normal Hartree-Fock state. Thus, it appears that the Bogoliubov method may lead to completely erroneous interpretations, when the coupling constants are small.
We study a system of trapped bosonic particles interacting by model harmonic forces. Our model allows for detailed examination of the notion of an order parameter (a condensate wave function). By decomposing a single particle density matrix into coherent eigenmodes we study an effect of interaction on the condensate. We show that sufficiently strong interactions cause that the condensate disappears even if the whole system is in its lowest energy state. In the second part of our paper we discuss the validity of the Bogoliubov approximation by comparing its predictions with results inferred from the exactl y soluble model. In particular we examine an energy spectrum, occupation, and fluctuations of the condensate. We conclude that Bogoliubov approach gives quite accurate description of the system in the limit of weak interactions.
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.