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:1203.4586 · Physical Review A

Occupation numbers of the harmonically trapped few-boson system

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

Published May 14, 2012• 15 citations
arXivHTMLPDF

Abstract

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.''

DOI
10.1103/physreva.85.053614

Related papers

Harmonically trapped four-boson system

D. Blume, Michelle Wynne Sze, John L. Bohn

Four identical spinless bosons with purely attractive two-body short-range interactions and repulsive three-body interactions under external spherically symmetric harmonic confinement are considered. The repulsive three-body potential prevents the formation of deeply bound states with molecular character. The low-energy spectrum with vanishing orbital angular momentum and positive parity for infinitely large two-body $s$-wave scattering length is analyzed in detail. Using the three-body contact, states are classified as universal, quasiuniversal, or strongly nonuniversal. Connections with the zero-range interaction model are discussed. The energy spectrum is mapped out as a function of the two-body $s$-wave scattering length ${a}_{s}, {a}_{s}>0$. In the weakly to medium-strongly interacting regime, one of the states approaches the energy obtained for a hard-core interaction model. This state is identified as the energetically lowest-lying ``BEC state.'' Structural properties are also presented.

arXiv:1802.00129•Physical Review A•Mar 26, 2018•13 citations

Effective renormalized multi-body interactions of harmonically confined ultracold neutral bosons

P R Johnson, D. Blume, Xiangyu Yin +2 more

We calculate the renormalized effective two-, three- and four-body interactions for N neutral ultracold bosons in the ground state of an isotropic harmonic trap, assuming two-body interactions modeled with the combination of a zero-range and energy-dependent pseudopotential. We work to third-order in the scattering length a t (0) defined at zero collision energy, which is necessary to obtain both the leading-order effective four-body interaction and consistently include corrections for realistic two-body interactions. The leading-order, effective three- and four-body interaction energies are and , where ω and σ ( ω ) are the harmonic oscillator frequency and length, respectively, and energies are in units of ℏ ω . The one-standard deviation error ±0.0001 for the third-order coefficient in U 3 ( ω ) is due to numerical uncertainty in estimating a slowly converging sum; the other two coefficients are either analytically or numerically exact. The effective three- and four-body interactions can play an important role in the dynamics of tightly confined and strongly correlated systems. We also performed numerical simulations for a finite-range (FR) boson–boson potential, and it was comparison to the zero-range predictions which revealed that finite-range effects must be taken into account for a realistic third-order treatment. In particular, we show that the energy-dependent pseudopotential accurately captures, through third order, the finite-range physics, and in combination with the multi-body effective interactions gives excellent agreement with the numerical simulations, validating our theoretical analysis and predictions.

arXiv:1201.2962•New Journal of Physics•May 25, 2012•29 citations

Weakly-Interacting Bosons in a Trap within Approximate Second Quantization Approach

Andrij A. Rovenchak

The theory of Bogoliubov is generalized for the case of a weakly-interacting Bose-gas in harmonic trap. A set of nonlinear matrix equations is obtained to make the diagonalization of Hamiltonian possible. Its perturbative solution is used for the calculation of the energy and the condensate fraction of the model system to show the applicability of the method. PACS numbers: 03.75.Hh, 05.30.−d 1.

•Nov 21, 2012•3 citations

Exact diagonalization of the Hamiltonian for trapped interacting bosons in lower dimensions

Tor Haugset, Hårek Haugerud

We consider systems of a small number of interacting bosons confined to harmonic potentials in one and two dimensions. By exact numerical diagonalization of the many-body Hamiltonian we determine the low-lying excitation energies and the ground-state energy and density profile. We discuss the dependence of these quantities on both interaction strength $g$ and particle number $N$. The ground-state properties are compared to the predictions of the Gross-Pitaevskii equation, and the agreement is surprisingly good even for relatively low particle numbers. We also calculate the specific heat based on the obtained energy spectra.

arXiv:cond-mat/9707107•Physical Review A•May 1, 1998•60 citations

Analytic harmonic approach to theN-body problem

J. R. Armstrong, Nikolaj Thomas Zinner, D. V. Fedorov +1 more

We consider an analytic way to make the interacting N -body problem tractable by using harmonic oscillators in place of the relevant two-body interactions. The two-body terms of the N -body Hamiltonian are approximated by considering the energy spectrum and radius of the relevant two-body problem which gives frequency, centre position, and zero point energy of the corresponding harmonic oscillator. Adding external harmonic one-body terms, we proceed to solve the full quantum mechanical N -body problem analytically for arbitrary masses. Energy eigenvalues, eigenmodes, and correlation functions like density matrices can then be computed analytically. As a first application of our formalism, we consider the N -boson problem in two and three dimensions where we fit the two-body interactions to agree with the well-known zero-range model for two particles in a harmonic trap. Subsequently, condensate fractions, spectra, radii, and eigenmodes are discussed as a function of dimension, boson number N , and scattering length obtained in the zero-range model. We find that energies, radii, and condensate fraction increase with scattering length as well as boson number, while radii decrease with increasing boson number. Our formalism is completely general and can also be applied to fermions, Bose–Fermi mixtures, and to more exotic geometries.

arXiv:1011.2453•Journal of Physics B Atomic Molecular and Optical Physics•Feb 18, 2011•48 citations

BOSE-EINSTEIN CONDENSATION IN TRAPS: A QUANTUM MONTE CARLO STUDY

Jonathan L. DuBois

We evaluate the zero temperature properties and Bose-Einstein condensation (BEC) of hard sphere bosons in a harmonic trap using Quantum Monte Carlo methods. The chief purpose of this work is to go beyond the dilute limit, to test the limits of the Gross-Pitaevskii equation and related mean field approximations and to explore the role of interactions in determining the zero temperature properties of the trapped Bose gas. The density is increased by adjusting both the number of trapped bosons, N , and the s-wave scattering length, a, to study systems from the highly dilute region corresponding to early experiments in trapped BEC [2, 3] up to liquid He densities and beyond. Rather than relying on the perturbative Bogoliubov approximation (which assumes a large condensate fraction) to describe the condensate, we use the one-body density matrix (OBDM) formulation of BEC so that the properties of the condensate for systems with arbitrarily large interactions may be studied. In this formulation of BEC, condensate properties are obtained by diagonalizing the OBDM and obtaining the corresponding single particle “natural orbitals” and their occupation numbers for the system. The condensate wave function and condensate fraction are then obtained from the single particle orbital(s) with macroscopic occupation (N0 >> 1). Within this framework, we calculate the effects of interactions and increased density on the ground state energy, the density profile, the momentum distribution, the condensate fraction and condensate “wave-function” and several other properties.

•Jan 1, 2002•1 citation
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⩾2N\ensuremath{\geqslant}2N\ensuremath⩾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.''
Four identical spinless bosons with purely attractive two-body short-range interactions and repulsive three-body interactions under external spherically symmetric harmonic confinement are considered. The repulsive three-body potential prevents the formation of deeply bound states with molecular character. The low-energy spectrum with vanishing orbital angular momentum and positive parity for infinitely large two-body sss-wave scattering length is analyzed in detail. Using the three-body contact, states are classified as universal, quasiuniversal, or strongly nonuniversal. Connections with the zero-range interaction model are discussed. The energy spectrum is mapped out as a function of the two-body sss-wave scattering length as,as>0{a}_{s}, {a}_{s}>0as​,as​>0. In the weakly to medium-strongly interacting regime, one of the states approaches the energy obtained for a hard-core interaction model. This state is identified as the energetically lowest-lying ``BEC state.'' Structural properties are also presented.
We calculate the renormalized effective two-, three- and four-body interactions for N neutral ultracold bosons in the ground state of an isotropic harmonic trap, assuming two-body interactions modeled with the combination of a zero-range and energy-dependent pseudopotential. We work to third-order in the scattering length a t (0) defined at zero collision energy, which is necessary to obtain both the leading-order effective four-body interaction and consistently include corrections for realistic two-body interactions. The leading-order, effective three- and four-body interaction energies are and , where ω and σ ( ω ) are the harmonic oscillator frequency and length, respectively, and energies are in units of ℏ ω . The one-standard deviation error ±0.0001 for the third-order coefficient in U 3 ( ω ) is due to numerical uncertainty in estimating a slowly converging sum; the other two coefficients are either analytically or numerically exact. The effective three- and four-body interactions can play an important role in the dynamics of tightly confined and strongly correlated systems. We also performed numerical simulations for a finite-range (FR) boson–boson potential, and it was comparison to the zero-range predictions which revealed that finite-range effects must be taken into account for a realistic third-order treatment. In particular, we show that the energy-dependent pseudopotential accurately captures, through third order, the finite-range physics, and in combination with the multi-body effective interactions gives excellent agreement with the numerical simulations, validating our theoretical analysis and predictions.
We consider systems of a small number of interacting bosons confined to harmonic potentials in one and two dimensions. By exact numerical diagonalization of the many-body Hamiltonian we determine the low-lying excitation energies and the ground-state energy and density profile. We discuss the dependence of these quantities on both interaction strength ggg and particle number NNN. The ground-state properties are compared to the predictions of the Gross-Pitaevskii equation, and the agreement is surprisingly good even for relatively low particle numbers. We also calculate the specific heat based on the obtained energy spectra.
We consider an analytic way to make the interacting N -body problem tractable by using harmonic oscillators in place of the relevant two-body interactions. The two-body terms of the N -body Hamiltonian are approximated by considering the energy spectrum and radius of the relevant two-body problem which gives frequency, centre position, and zero point energy of the corresponding harmonic oscillator. Adding external harmonic one-body terms, we proceed to solve the full quantum mechanical N -body problem analytically for arbitrary masses. Energy eigenvalues, eigenmodes, and correlation functions like density matrices can then be computed analytically. As a first application of our formalism, we consider the N -boson problem in two and three dimensions where we fit the two-body interactions to agree with the well-known zero-range model for two particles in a harmonic trap. Subsequently, condensate fractions, spectra, radii, and eigenmodes are discussed as a function of dimension, boson number N , and scattering length obtained in the zero-range model. We find that energies, radii, and condensate fraction increase with scattering length as well as boson number, while radii decrease with increasing boson number. Our formalism is completely general and can also be applied to fermions, Bose–Fermi mixtures, and to more exotic geometries.
We evaluate the zero temperature properties and Bose-Einstein condensation (BEC) of hard sphere bosons in a harmonic trap using Quantum Monte Carlo methods. The chief purpose of this work is to go beyond the dilute limit, to test the limits of the Gross-Pitaevskii equation and related mean field approximations and to explore the role of interactions in determining the zero temperature properties of the trapped Bose gas. The density is increased by adjusting both the number of trapped bosons, N , and the s-wave scattering length, a, to study systems from the highly dilute region corresponding to early experiments in trapped BEC [2, 3] up to liquid He densities and beyond. Rather than relying on the perturbative Bogoliubov approximation (which assumes a large condensate fraction) to describe the condensate, we use the one-body density matrix (OBDM) formulation of BEC so that the properties of the condensate for systems with arbitrarily large interactions may be studied. In this formulation of BEC, condensate properties are obtained by diagonalizing the OBDM and obtaining the corresponding single particle “natural orbitals” and their occupation numbers for the system. The condensate wave function and condensate fraction are then obtained from the single particle orbital(s) with macroscopic occupation (N0 >> 1). Within this framework, we calculate the effects of interactions and increased density on the ground state energy, the density profile, the momentum distribution, the condensate fraction and condensate “wave-function” and several other properties.