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

arXiv:cond-mat/9911210 · Physical Review A

Condensation and interaction range in harmonic boson traps: A variational approach

Jacques Tempere, F. Brosens, L. F. Lemmens, J. T. Devreese

Published Mar 9, 2000• 19 citations
arXivHTMLPDF

Abstract

For a gas of N bosons interacting through a two-body Morse potential a variational bound of the free energy of a confined system is obtained. The calculation method is based on the Feynman-Kac functional projected on the symmetric representation. Within the harmonic approximation a variational estimate of the effect of the interaction range on the existence of many-particle bound states, and on the $N\ensuremath{-}T$ phase diagram is obtained.

DOI
10.1103/physreva.61.043605

Related papers

Bose–Einstein condensate in a double-well potential: Feynman path integral variational approach

Wattana Ratismith, Sigmund Heller, Virulh Sa‐Yakanit +1 more

A Bose–Einstein condensate in a double-well potential is considered by using the Feynman path integral theory combined with a variational approach. The system consists of N interacting bosons confined in the double well which is taken as a harmonic potential with a Gaussian barrier. The calculation is carried out within the first cumulant approximation measured with respect to the harmonic action containing variational parameters. Assuming a separable expression for the trial action corresponds to the usual mean field approximation. Performing the variational calculations, we obtain analytical results for the ground state energy and its condensate wavefunction. Using a projection onto the odd contributions, the first excited state energy and condensate wavefunction are determined, too. We find very good agreement with a full numerical solution of the Gross–Pitaevskii equation over the full range of potential parameters.

Journal of Physics A Mathematical and Theoretical•Apr 3, 2013•1 citation

Variational Monte Carlo analysis of Bose–Einstein condensation in a two-dimensional trap

Zheng Rong-Jie, Jing Jin, Yi Ping Tang

The ground-state properties of a system with a small number of interacting bosons over a wide range of densities are investigated. The system is confined in a two-dimensional isotropic harmonic trap, where the interaction between bosons is treated as a hard-core potential. By using variational Monte Carlo method, we diagonalize the one-body density matrix of the system to obtain the ground-state energy, condensate wavefunction and the condensate fraction. We find that in the dilute limit the depletion of central condensate in the 2D system is larger than in a 3D system for the same interaction strength; however as the density increases, the depletion at the centre of 2D trap will be equal to or even lower than that at the centre of 3D trap, which is in agreement with the anticipated in Thomas–Fermi approximation. In addition, in the 2D system the total condensate depletion is still larger than in a 3D system for the same scattering length.

Chinese Physics•Aug 22, 2006•6 citations

Two-dimensional bosonic droplets in a harmonic trap

Fabian Brauneis, Artem G. Volosniev, H.‐W. Hammer

We investigate a system of bosons in a two-dimensional harmonic trap. In the limit of strong attractive interactions, the bosons make a droplet insensitive to external confinement. For weak interactions, in contrast, the ground state is given by the harmonic trap. In this work, we conduct a variational study of the transition between these two limits. We find that this transition occurs at the critical interaction strength whose value is universal if scaled appropriately with the number of particles. To connect the change in the properties of the system to the classical description of phase transitions, we analyze the static response of the Bose gas related to the isothermal compressibility. Finally, we perform numerically exact calculations for a few particles to demonstrate the effects of finite-range interactions on this transition. We conclude that finite-range effects wash out the point of transition.

arXiv:2410.00570•Physical Review A•Jun 10, 2025•2 citations

Occupation numbers of the harmonically trapped few-boson system

K. M. Daily, Xiangyu Yin, D. Blume

We consider a harmonically trapped dilute $N$-boson system described by a low-energy Hamiltonian with pairwise interactions. We determine the condensate fraction, defined in terms of the largest occupation number, of the weakly interacting $N$-boson system ($N\ensuremath{\geqslant}2$) by employing a perturbative treatment within the framework of second quantization. The one-body density matrix and the corresponding occupation numbers are compared with those obtained by solving the two-body problem with zero-range interactions exactly. Our expressions are also compared with high precision ab initio calculations for Bose gases with $N=2--4$ that interact through finite-range two-body model potentials. Nonuniversal corrections are identified to enter at subleading order, confirming that different low-energy Hamiltonians, constructed to yield the same energy, may yield different occupation numbers. Lastly, we consider the strongly interacting three-boson system under spherically symmetric harmonic confinement and determine its occupation numbers as a function of the three-body ``Efimov parameter.''

arXiv:1203.4586•Physical Review A•May 14, 2012•15 citations

Ground state properties of trapped boson system with finite-range Gaussian repulsion: exact diagonalization study

Mohd. Imran, M. A. H. Ahsan

Abstract We use exact diagonalization to study an interacting system of N spinless bosons with finite-range Gaussian repulsion, confined in a quasi-two-dimensional harmonic trap with and without an introduced rotation. The diagonalization of the Hamiltonian matrix using Davidson algorithm in subspaces of quantized total angular momentum L z is carried out to obtain the N -body lowest eigenenergy and eigenstate. To bring out the effect of quantum (Bose) statistics and consequent phase stiffness (rigidity) of the variationally obtained many-body wavefunction on various physical quantities, our study spans from few-body ( N = 2) to many-body ( N = 16) systems. Further, to examine the finite-range effect of the repulsive Gaussian potential on many-body ground state properties of the Bose-condensate, we obtain the lowest eigenstate, the critical angular velocity of single vortex state and the quantum correlation (measured) in terms of von Neumann entanglement entropy and degree of condensation. It is found that for small values of the range (measured by the parameter σ ) of Gaussian potential, the ground state energy increases for few-boson (2 ⩽ N ⩽ 8) systems but decreases for many-boson ( N > 8) systems. On the other hand for relatively large values of the range of Gaussian potential, the ground state energy exhibits a monotonic decrease, regardless of the number of bosons N . For a given N , there is found an optimal value of the range of Gaussian potential for which the first vortex (with L z = N ) nucleates at a lower value of the rotational angular velocity Ω c1 compared to the zero-range ( δ -function) potential. Further, we observe that the inter-particle interaction and the introduced rotation are competing effects with latter being dominant over the former. With increase in the range of Gaussian potential, the value of von Neumann entropy decreases and the degree of condensation increases implying an enhanced quantum correlation and phase rigidity.

arXiv:1602.03811•Journal of Physics B Atomic Molecular and Optical Physics•Mar 13, 2020•3 citations

Bose-Einstein condensation in a two-dimensional, trapped, interacting gas

Mehmet Bayındır, B. Tanatar

We study the Bose-Einstein condensation phenomenon in a two-dimensional (2D) system of bosons subjected to a harmonic-oscillator-type confining potential. The interaction among the 2D bosons is described by a \ensuremath{\delta} function in configuration space. Solving the Gross-Pitaevskii equation within the two-fluid model we calculate the condensate fraction, ground-state energy, and specific heat of the system. Our results indicate that interacting bosons have similar behavior to those of an ideal system for weak interactions.

arXiv:cond-mat/9807377•Physical Review A•Oct 1, 1998•32 citations
For a gas of N bosons interacting through a two-body Morse potential a variational bound of the free energy of a confined system is obtained. The calculation method is based on the Feynman-Kac functional projected on the symmetric representation. Within the harmonic approximation a variational estimate of the effect of the interaction range on the existence of many-particle bound states, and on the N\ensuremath−TN\ensuremath{-}TN\ensuremath−T phase diagram is obtained.
A Bose–Einstein condensate in a double-well potential is considered by using the Feynman path integral theory combined with a variational approach. The system consists of N interacting bosons confined in the double well which is taken as a harmonic potential with a Gaussian barrier. The calculation is carried out within the first cumulant approximation measured with respect to the harmonic action containing variational parameters. Assuming a separable expression for the trial action corresponds to the usual mean field approximation. Performing the variational calculations, we obtain analytical results for the ground state energy and its condensate wavefunction. Using a projection onto the odd contributions, the first excited state energy and condensate wavefunction are determined, too. We find very good agreement with a full numerical solution of the Gross–Pitaevskii equation over the full range of potential parameters.
The ground-state properties of a system with a small number of interacting bosons over a wide range of densities are investigated. The system is confined in a two-dimensional isotropic harmonic trap, where the interaction between bosons is treated as a hard-core potential. By using variational Monte Carlo method, we diagonalize the one-body density matrix of the system to obtain the ground-state energy, condensate wavefunction and the condensate fraction. We find that in the dilute limit the depletion of central condensate in the 2D system is larger than in a 3D system for the same interaction strength; however as the density increases, the depletion at the centre of 2D trap will be equal to or even lower than that at the centre of 3D trap, which is in agreement with the anticipated in Thomas–Fermi approximation. In addition, in the 2D system the total condensate depletion is still larger than in a 3D system for the same scattering length.
We investigate a system of bosons in a two-dimensional harmonic trap. In the limit of strong attractive interactions, the bosons make a droplet insensitive to external confinement. For weak interactions, in contrast, the ground state is given by the harmonic trap. In this work, we conduct a variational study of the transition between these two limits. We find that this transition occurs at the critical interaction strength whose value is universal if scaled appropriately with the number of particles. To connect the change in the properties of the system to the classical description of phase transitions, we analyze the static response of the Bose gas related to the isothermal compressibility. Finally, we perform numerically exact calculations for a few particles to demonstrate the effects of finite-range interactions on this transition. We conclude that finite-range effects wash out the point of transition.
We consider a harmonically trapped dilute NNN-boson system described by a low-energy Hamiltonian with pairwise interactions. We determine the condensate fraction, defined in terms of the largest occupation number, of the weakly interacting NNN-boson system ($N\ensuremath{\geqslant}2$) by employing a perturbative treatment within the framework of second quantization. The one-body density matrix and the corresponding occupation numbers are compared with those obtained by solving the two-body problem with zero-range interactions exactly. Our expressions are also compared with high precision ab initio calculations for Bose gases with N=2−−4N=2--4N=2−−4 that interact through finite-range two-body model potentials. Nonuniversal corrections are identified to enter at subleading order, confirming that different low-energy Hamiltonians, constructed to yield the same energy, may yield different occupation numbers. Lastly, we consider the strongly interacting three-boson system under spherically symmetric harmonic confinement and determine its occupation numbers as a function of the three-body ``Efimov parameter.''
Abstract We use exact diagonalization to study an interacting system of N spinless bosons with finite-range Gaussian repulsion, confined in a quasi-two-dimensional harmonic trap with and without an introduced rotation. The diagonalization of the Hamiltonian matrix using Davidson algorithm in subspaces of quantized total angular momentum L z is carried out to obtain the N -body lowest eigenenergy and eigenstate. To bring out the effect of quantum (Bose) statistics and consequent phase stiffness (rigidity) of the variationally obtained many-body wavefunction on various physical quantities, our study spans from few-body ( N = 2) to many-body ( N = 16) systems. Further, to examine the finite-range effect of the repulsive Gaussian potential on many-body ground state properties of the Bose-condensate, we obtain the lowest eigenstate, the critical angular velocity of single vortex state and the quantum correlation (measured) in terms of von Neumann entanglement entropy and degree of condensation. It is found that for small values of the range (measured by the parameter σ ) of Gaussian potential, the ground state energy increases for few-boson (2 ⩽ N ⩽ 8) systems but decreases for many-boson ( N > 8) systems. On the other hand for relatively large values of the range of Gaussian potential, the ground state energy exhibits a monotonic decrease, regardless of the number of bosons N . For a given N , there is found an optimal value of the range of Gaussian potential for which the first vortex (with L z = N ) nucleates at a lower value of the rotational angular velocity Ω c1 compared to the zero-range ( δ -function) potential. Further, we observe that the inter-particle interaction and the introduced rotation are competing effects with latter being dominant over the former. With increase in the range of Gaussian potential, the value of von Neumann entropy decreases and the degree of condensation increases implying an enhanced quantum correlation and phase rigidity.
We study the Bose-Einstein condensation phenomenon in a two-dimensional (2D) system of bosons subjected to a harmonic-oscillator-type confining potential. The interaction among the 2D bosons is described by a \ensuremath{\delta} function in configuration space. Solving the Gross-Pitaevskii equation within the two-fluid model we calculate the condensate fraction, ground-state energy, and specific heat of the system. Our results indicate that interacting bosons have similar behavior to those of an ideal system for weak interactions.