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 A

Condensate fluctuation and thermodynamics of mesoscopic Bose-Einstein condensates: A correlated many-body approach

Satadal Bhattacharyya, Barnali Chakrabarti

Published Feb 23, 2016• 11 citations
Publisher

Abstract

We present a correlated many-body approach to calculate the distribution function and fluctuations for a Bose-Einstein condensate with $N$ interacting atoms in the harmonic confinement. The present formulation uses the recursion relation for the canonical ensemble partition function $(Z)$. $Z$ is calculated from the energy spectrum of the many-body effective potential, which keeps all possible two-body correlations and uses the realistic interatomic interaction. The condensate statistics are in very good agreement with earlier results of an ideal gas for which exact statistical moments for all temperature are known. We also present the numerical results of condensate statistics for real experimental situations. The calculated moments nicely exhibit the mesoscopic effect for a few hundred atoms, whereas the sharp fall in the variance for the large condensate near the critical temperature shows the possibility of phase transition. We also calculate the critical temperature for the mesoscopic regime. Our present calculation mimics the JILA experiment with $^{87}\mathrm{Rb}$ atoms.

DOI
10.1103/physreva.93.023636

Related papers

Mesoscopic Effects in Bose-Einstein Condensate Fluctuations of an Ideal Gas in a Box

Konstantin E. Dorfman

The mesoscopic effects in the quantum trapped gases of the Bose atoms constitute\nthe main subject of the present thesis. These effects are the most difficult for\nthe theoretical analysis in the quantum statistical physics since they can’t be seen by\nneither a standard quantum mechanics of the simple microscopic systems of one or\nvery few atoms nor a standard statistical physics of the macroscopic systems that are\ninfinite in the bulk (thermodynamic) limit.\nMost of the experiments on the cold quantum gases performed in the last decade,\nstarting from the first demonstration of BEC in 1995, involve the mesoscopic systems\nof a finite number of atoms. The mesoscopic effects should manifest themselves most\nclearly and easily near a critical temperature of BEC; however, they could be observed\nalso above and below the critical temperature.\nHere I study the quantum and thermal fluctuations of the Bose-Einstein condensate\n(BEC) in a box with the periodic boundary conditions under a particle-number\nconstraint. The above constraint is the only reason for the BEC and is crucial for the\nmesoscopic effects in the BEC fluctuations, especially in the vicinity of the critical\ntemperature in the Bose gas.\nI employ the particle-number conserving operator formalism of Girardeau and\nArnowitt introduced in 1959 to analyze the canonical ensemble fluctuations. I present\nanalytical formulas and numerical calculations for the central moments of the ground\nstate occupation fluctuations in an ideal Bose gas in a box with a mesoscopic number of particles.\nI present the analysis of the BEC statistics both on a temperature at a fixed\nnumber of particles and on a number of particles at a fixed temperature. Both analyses\nare valid for the purpose of understanding the important mesoscopic effects near the\ncritical temperature. I emphasize the non-Gaussian nature of the fluctuations.\nThe presented formalism can be generalized to the case of a weakly interacting\nBose gas in a box in the framework of the Bogoliubov approximation. The work in\nthis direction is in progress but is not included in the present thesis.

•May 15, 2009•0 citations

Condensate statistics and thermodynamics of weakly interacting Bose gas: Recursion relation approach

Konstantin E. Dorfman, M. Kim, Anatoly A. Svidzinsky

We study condensate statistics and thermodynamics of weakly interacting Bose gas with a fixed total number $N$ of particles in a cubic box. We find the exact recursion relation for the canonical ensemble partition function. Using this relation, we calculate the distribution function of condensate particles for $N=200$. We also calculate the distribution function based on multinomial expansion of the characteristic function. Similar to the ideal gas, both approaches give exact statistical moments for all temperatures in the framework of Bogoliubov model. We compare them with the results of unconstraint canonical ensemble quasiparticle formalism and the hybrid master equation approach. The present recursion relation can be used for any external potential and boundary conditions. We investigate the temperature dependence of the first few statistical moments of condensate fluctuations as well as thermodynamic potentials and heat capacity analytically and numerically in the whole temperature range.

Physical Review A•Mar 14, 2011•8 citations

Fluctuations of a weakly interacting Bose-Einstein condensate

Zbigniew Idziaszek, Łukasz Zawitkowski, Mariusz Gajda +1 more

Fluctuations of the number of condensed atoms in a finite-size, weakly interacting Bose gas confined in a box potential are investigated for temperatures up to the critical region. The canonical partition functions are evaluated using a recursive scheme for smaller systems, and a saddle-point approximation for larger samples that allows to treat realistic size systems containing up to N ∼10 5 particles. We point out the importance of particle-number constraint and interactions between out of condensate atoms for the statistics near the critical region. For sufficiently large systems, the crossover from the anomalous to normal scaling of the fluctuations is observed. The excitations are described in a self-consistent way within the Bogoliubov-Popov approximation, and the interactions between thermal atoms are described by means of the Hartree-Fock method.

arXiv:0708.0092•Europhysics Letters (EPL)•Apr 1, 2009•16 citations

Energy fluctuation of a finite number of interacting bosons: A correlated many-body approach

Satadal Bhattacharyya, M. L. Lekala, Barnali Chakrabarti +1 more

We calculate the energy fluctuation of a truly finite number of interacting bosons and study the role of interaction. Although the ideal Bose gas in thermodynamic limit is an exactly solvable problem and analytic expression of various fluctuation measures exists, the experimental Bose-Einstein condensation (BEC) is a nontrivial many-body problem. We employ a two-body correlated basis function and utilize the realistic van der Waals interaction. We calculate the energy fluctuation $(△{E}^{2})$ of the interacting trapped bosons and plot $\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}$ as a function of $\frac{T}{{T}_{c}}$. In the classical limit $△{E}^{2}$ is related to the specific heat per particle ${c}_{v}$ through the relation $△{E}^{2}={k}_{B}{T}^{2}{c}_{v}$. We have obtained a distinct hump in $\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}$ around the condensation point for three-dimesional harmonically trapped Bose gas when the particle number $N\ensuremath{\simeq}5000$ and above which corresponds to the second-order phase transition. However for finite-size interacting bosons ($N\ensuremath{\simeq}$ a few hundred) the hump is not sharp, and the maximum in $\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}$ can be interpreted as a smooth increase in the scaled fluctuation below ${T}_{c}$ and then a decrease above ${T}_{c}$. To illustrate the justification we also calculate ${c}_{v}$, which exhibits the same feature, which leads to the conjecture that for finite-sized interacting bosons phase transition is ruled out.

Physical Review A•Mar 14, 2016•4 citations

Condensate fluctuations of interacting Bose gases within a microcanonical ensemble

Jianhui Wang, Jizhou He, Yong-li Ma

Based on counting statistics and Bogoliubov theory, we present a recurrence relation for the microcanonical partition function for a weakly interacting Bose gas with a finite number of particles in a cubic box. According to this microcanonical partition function, we calculate numerically the distribution function, condensate fraction, and condensate fluctuations for a finite and isolated Bose-Einstein condensate. For ideal and weakly interacting Bose gases, we compare the condensate fluctuations with those in the canonical ensemble. The present approach yields an accurate account of the condensate fluctuations for temperatures close to the critical region. We emphasize that the interactions between excited atoms turn out to be important for moderate temperatures.

Physical Review E•May 31, 2011•14 citations

Analytical theory of mesoscopic Bose-Einstein condensation in an ideal gas

Vitaly V. Kocharovsky, Vl. V. Kocharovsky

We find the universal structure and scaling of the Bose-Einstein condensation (BEC) statistics and thermodynamics (Gibbs free energy, average energy, heat capacity) for a mesoscopic canonical-ensemble ideal gas in a trap with an arbitrary number of atoms, any volume, and any temperature, including the whole critical region. We identify a universal constraint-cutoff mechanism that makes BEC fluctuations strongly non-Gaussian and is responsible for all unusual critical phenomena of the BEC phase transition in the ideal gas. The main result is an analytical solution to the problem of critical phenomena. It is derived by, first, calculating analytically the universal probability distribution of the noncondensate occupation, or a Landau function, and then using it for the analytical calculation of the universal functions for the particular physical quantities via the exact formulas which express the constraint-cutoff mechanism. We find asymptotics of that analytical solution as well as its simple analytical approximations which describe the universal structure of the critical region in terms of the parabolic cylinder or confluent hypergeometric functions. The obtained results for the order parameter, all higher-order moments of BEC fluctuations, and thermodynamic quantities perfectly match the known asymptotics outside the critical region for both low and high temperature limits. We suggest two- and three-level trap models of BEC and find their exact solutions in terms of the cutoff negative binomial distribution (which tends to the cutoff gamma distribution in the continuous limit) and the confluent hypergeometric distribution, respectively. Also, we present an exactly solvable cutoff Gaussian model of BEC in a degenerate interacting gas. All these exact solutions confirm the universality and constraint-cutoff origin of the strongly non-Gaussian BEC statistics. We introduce a regular refinement scheme for the condensate statistics approximations on the basis of the infrared universality of higher-order cumulants and the method of superposition and show how to model BEC statistics in the actual traps. In particular, we find that the three-level trap model with matching the first four or five cumulants is enough to yield remarkably accurate results for all interesting quantities in the whole critical region. We derive an exact multinomial expansion for the noncondensate occupation probability distribution and find its high-temperature asymptotics (Poisson distribution) and corrections to it. Finally, we demonstrate that the critical exponents and a few known terms of the Taylor expansion of the universal functions, which were calculated previously from fitting the finite-size simulations within the phenomenological renormalization-group theory, can be easily obtained from the presented full analytical solutions for the mesoscopic BEC as certain approximations in the close vicinity of the critical point.

arXiv:1511.01610•Physical Review A•Mar 23, 2010•31 citations
We present a correlated many-body approach to calculate the distribution function and fluctuations for a Bose-Einstein condensate with NNN interacting atoms in the harmonic confinement. The present formulation uses the recursion relation for the canonical ensemble partition function (Z)(Z)(Z). ZZZ is calculated from the energy spectrum of the many-body effective potential, which keeps all possible two-body correlations and uses the realistic interatomic interaction. The condensate statistics are in very good agreement with earlier results of an ideal gas for which exact statistical moments for all temperature are known. We also present the numerical results of condensate statistics for real experimental situations. The calculated moments nicely exhibit the mesoscopic effect for a few hundred atoms, whereas the sharp fall in the variance for the large condensate near the critical temperature shows the possibility of phase transition. We also calculate the critical temperature for the mesoscopic regime. Our present calculation mimics the JILA experiment with 87Rb^{87}\mathrm{Rb}87Rb atoms.
The mesoscopic effects in the quantum trapped gases of the Bose atoms constitute\nthe main subject of the present thesis. These effects are the most difficult for\nthe theoretical analysis in the quantum statistical physics since they can’t be seen by\nneither a standard quantum mechanics of the simple microscopic systems of one or\nvery few atoms nor a standard statistical physics of the macroscopic systems that are\ninfinite in the bulk (thermodynamic) limit.\nMost of the experiments on the cold quantum gases performed in the last decade,\nstarting from the first demonstration of BEC in 1995, involve the mesoscopic systems\nof a finite number of atoms. The mesoscopic effects should manifest themselves most\nclearly and easily near a critical temperature of BEC; however, they could be observed\nalso above and below the critical temperature.\nHere I study the quantum and thermal fluctuations of the Bose-Einstein condensate\n(BEC) in a box with the periodic boundary conditions under a particle-number\nconstraint. The above constraint is the only reason for the BEC and is crucial for the\nmesoscopic effects in the BEC fluctuations, especially in the vicinity of the critical\ntemperature in the Bose gas.\nI employ the particle-number conserving operator formalism of Girardeau and\nArnowitt introduced in 1959 to analyze the canonical ensemble fluctuations. I present\nanalytical formulas and numerical calculations for the central moments of the ground\nstate occupation fluctuations in an ideal Bose gas in a box with a mesoscopic number of particles.\nI present the analysis of the BEC statistics both on a temperature at a fixed\nnumber of particles and on a number of particles at a fixed temperature. Both analyses\nare valid for the purpose of understanding the important mesoscopic effects near the\ncritical temperature. I emphasize the non-Gaussian nature of the fluctuations.\nThe presented formalism can be generalized to the case of a weakly interacting\nBose gas in a box in the framework of the Bogoliubov approximation. The work in\nthis direction is in progress but is not included in the present thesis.
We study condensate statistics and thermodynamics of weakly interacting Bose gas with a fixed total number NNN of particles in a cubic box. We find the exact recursion relation for the canonical ensemble partition function. Using this relation, we calculate the distribution function of condensate particles for N=200N=200N=200. We also calculate the distribution function based on multinomial expansion of the characteristic function. Similar to the ideal gas, both approaches give exact statistical moments for all temperatures in the framework of Bogoliubov model. We compare them with the results of unconstraint canonical ensemble quasiparticle formalism and the hybrid master equation approach. The present recursion relation can be used for any external potential and boundary conditions. We investigate the temperature dependence of the first few statistical moments of condensate fluctuations as well as thermodynamic potentials and heat capacity analytically and numerically in the whole temperature range.
Fluctuations of the number of condensed atoms in a finite-size, weakly interacting Bose gas confined in a box potential are investigated for temperatures up to the critical region. The canonical partition functions are evaluated using a recursive scheme for smaller systems, and a saddle-point approximation for larger samples that allows to treat realistic size systems containing up to N ∼10 5 particles. We point out the importance of particle-number constraint and interactions between out of condensate atoms for the statistics near the critical region. For sufficiently large systems, the crossover from the anomalous to normal scaling of the fluctuations is observed. The excitations are described in a self-consistent way within the Bogoliubov-Popov approximation, and the interactions between thermal atoms are described by means of the Hartree-Fock method.
We calculate the energy fluctuation of a truly finite number of interacting bosons and study the role of interaction. Although the ideal Bose gas in thermodynamic limit is an exactly solvable problem and analytic expression of various fluctuation measures exists, the experimental Bose-Einstein condensation (BEC) is a nontrivial many-body problem. We employ a two-body correlated basis function and utilize the realistic van der Waals interaction. We calculate the energy fluctuation (△E2)(△{E}^{2})(△E2) of the interacting trapped bosons and plot △E2kB2T2\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}kB2​T2△E2​ as a function of TTc\frac{T}{{T}_{c}}Tc​T​. In the classical limit △E2△{E}^{2}△E2 is related to the specific heat per particle cv{c}_{v}cv​ through the relation △E2=kBT2cv△{E}^{2}={k}_{B}{T}^{2}{c}_{v}△E2=kB​T2cv​. We have obtained a distinct hump in △E2kB2T2\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}kB2​T2△E2​ around the condensation point for three-dimesional harmonically trapped Bose gas when the particle number $N\ensuremath{\simeq}5000$ and above which corresponds to the second-order phase transition. However for finite-size interacting bosons ($N\ensuremath{\simeq}$ a few hundred) the hump is not sharp, and the maximum in △E2kB2T2\frac{△{E}^{2}}{{k}_{B}^{2}{T}^{2}}kB2​T2△E2​ can be interpreted as a smooth increase in the scaled fluctuation below Tc{T}_{c}Tc​ and then a decrease above Tc{T}_{c}Tc​. To illustrate the justification we also calculate cv{c}_{v}cv​, which exhibits the same feature, which leads to the conjecture that for finite-sized interacting bosons phase transition is ruled out.
Based on counting statistics and Bogoliubov theory, we present a recurrence relation for the microcanonical partition function for a weakly interacting Bose gas with a finite number of particles in a cubic box. According to this microcanonical partition function, we calculate numerically the distribution function, condensate fraction, and condensate fluctuations for a finite and isolated Bose-Einstein condensate. For ideal and weakly interacting Bose gases, we compare the condensate fluctuations with those in the canonical ensemble. The present approach yields an accurate account of the condensate fluctuations for temperatures close to the critical region. We emphasize that the interactions between excited atoms turn out to be important for moderate temperatures.
We find the universal structure and scaling of the Bose-Einstein condensation (BEC) statistics and thermodynamics (Gibbs free energy, average energy, heat capacity) for a mesoscopic canonical-ensemble ideal gas in a trap with an arbitrary number of atoms, any volume, and any temperature, including the whole critical region. We identify a universal constraint-cutoff mechanism that makes BEC fluctuations strongly non-Gaussian and is responsible for all unusual critical phenomena of the BEC phase transition in the ideal gas. The main result is an analytical solution to the problem of critical phenomena. It is derived by, first, calculating analytically the universal probability distribution of the noncondensate occupation, or a Landau function, and then using it for the analytical calculation of the universal functions for the particular physical quantities via the exact formulas which express the constraint-cutoff mechanism. We find asymptotics of that analytical solution as well as its simple analytical approximations which describe the universal structure of the critical region in terms of the parabolic cylinder or confluent hypergeometric functions. The obtained results for the order parameter, all higher-order moments of BEC fluctuations, and thermodynamic quantities perfectly match the known asymptotics outside the critical region for both low and high temperature limits. We suggest two- and three-level trap models of BEC and find their exact solutions in terms of the cutoff negative binomial distribution (which tends to the cutoff gamma distribution in the continuous limit) and the confluent hypergeometric distribution, respectively. Also, we present an exactly solvable cutoff Gaussian model of BEC in a degenerate interacting gas. All these exact solutions confirm the universality and constraint-cutoff origin of the strongly non-Gaussian BEC statistics. We introduce a regular refinement scheme for the condensate statistics approximations on the basis of the infrared universality of higher-order cumulants and the method of superposition and show how to model BEC statistics in the actual traps. In particular, we find that the three-level trap model with matching the first four or five cumulants is enough to yield remarkably accurate results for all interesting quantities in the whole critical region. We derive an exact multinomial expansion for the noncondensate occupation probability distribution and find its high-temperature asymptotics (Poisson distribution) and corrections to it. Finally, we demonstrate that the critical exponents and a few known terms of the Taylor expansion of the universal functions, which were calculated previously from fitting the finite-size simulations within the phenomenological renormalization-group theory, can be easily obtained from the presented full analytical solutions for the mesoscopic BEC as certain approximations in the close vicinity of the critical point.