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/9511100 · Physical Review Letters

Surface Region of Superfluid Helium as an Inhomogeneous Bose-Condensed Gas

Adam Griffin, S. Stringari

Published Jan 8, 1996• 46 citations
arXivHTMLPDF

Abstract

We present arguments that the low density surface region of self-bounded superfluid ${}^{4}\mathrm{He}$ systems is an inhomogeneous dilute Bose gas, with almost all of the atoms occupying the same single-particle state at $T\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}0$. Numerical evidence for this complete Bose-Einstein condensation was first given by the many-body variational calculations of ${}^{4}\mathrm{He}$ droplets by Lewart, Pandharipande, and Pieper in 1988 [Phys. Rev. B 37, 4950 (1988)]. We show that the low density surface region can be treated rigorously using a generalized Gross-Pitaevskii equation for the Bose order parameter.

DOI
10.1103/physrevlett.76.259

Related papers

Equilibrium and Non-Equilibrium Properties of Dilute Bose Gases

Oldenburg, Jakob Felix

Quantum systems consisting of many particles are both interesting from the physics and the mathematics perspective. One aspect that makes them fascinating from the physics point of view is that they can exhibit quantum effects on macroscopic scales such as Bose-Einstein condensation. This phenomenon, first predicted 100 years ago by Bose [19] and Einstein [44], occurs when a dilute gas of bosons is cooled down to very low temperatures. Below a certain critical temperature a condensate forms, meaning that a macroscopic fraction of the particles is in the same one-particle state, and thus, quantum phenomena become relevant on length scales much larger than the individual atom. The occurrence of Bose-Einstein condensation is also related to superfluidity as suggested in [88]. Bose-Einstein condensation has been achieved experimentally independently by the groups of Cornell and Wieman [3] and Ketterle [40] earning them a Nobel prize in 2001. While Bose and Einstein analyzed theoretically the non-interacting gas, an effective theory describing weakly interacting Bose gases was proposed by Bogoliubov in 1947 with the goal of explaining the occurrence of superfluidity in liquid helium [18]. Using assumptions such as that condensation indeed occurs, his theory yields predictions for the low-energy states of the system. Seen from the mathematics perspective, describing the Bose gas is an interesting challenge due to the strong correlations between particles introduced by the interaction making this case much harder to study than the non-interacting one. In particular, a proof that Bose-Einstein condensation occurs in the thermodynamic limit for continuous systems is still out of reach. However, there has been a lot of work regarding the rigorous analysis of bosonic systems in various relevant regimes. This thesis contains three new contributions establishing rigorous results on bosonic systems in different, physically relevant settings. It is based on the articles [37, 27, 31]. The first result concerns the time-evolution of initially trapped Bose gases in the ultra-dilute Gross-Pitaevskii regime. Under some assumptions on the initial data satisfied in particular by a class of states that are expected to approximate the ground state in a trap, we derive a quasi-free approximation that is close to the true solution in $L^2$-norm. In particular, we show that this approximate wave function is described in terms of a Bose-Einstein condensate in a one-particle state following the time-dependent Gross-Pitaevskii equation together with pair excitations with respect to this condensate. As a corollary, we establish a central limit theorem for bounded one-particle observables. The second result is a simple proof of Bose-Einstein condensation for interacting, ultra-dilute gases. More specifically, we consider scalings interpolating between the Gross-Pitaevskii regime and the thermodynamic limit. While condensation had been shown to hold in these regimes before, the advantage of the method used here is its simplicity. Lastly, the third result establishes an upper bound on the ground state energy of a dilute Bose gas in the thermodynamic limit. It captures rigorously the third order correction to the energy per unit volume for sufficiently low densities as predicted in the physics literature first by Wu and subsequently others. $\textit{Full-text embargoed until: 2026-01-01}$

•Dec 9, 2025•0 citations

Path Integral Monte Carlo Simulations of Helium: From Superfluid Droplets to Quantum Crystals

Erik W. Draeger

Below Tλ = 2.17 K, bulk 4He is a superfluid and has a non-zero Bose-Einstein condensate fraction. This work will focus primarily on how phenomena such as superfluidity, Bose condensation and superfluid vortices are manifested in microscopic, inhomogeneous helium systems. Path Integral Monte Carlo is a powerful method for calculating the equilibrium properties of quantum systems at finite temperature. We have achieved linear scaling of computer time with number of particles through the use of neighbor lists, allowing us to simulate systems of several thousand atoms. We have derived a local superfluid estimator and used it to examine the microscopic superfluid response around a molecule rotating in a helium nanodroplet. We found that the first solvation layer is well-described by a two dimensional superfluid, with thermal excitations occuring at a lower temperature than in bulk helium. The effective moment of inertia of a linear impurity in a helium droplet is calculated, and compared with experimental scattering results. In addition, we calculated the vortex formation energy for both pure droplets and droplets doped with linear impurities, and found that the linear impurities had a negligible impact on the formation energy. A possible spectroscopic signature of vortices in doped helium droplets was suggested. After deriving a local estimator, we calculated the condensate fraction throughout the free helium surface of a semi-infinite slab. These results, along with densitydensity correlation functions, were used to characterize the surface excitations and calculate the extent to which ripplons are present. In addition, the ripplon dispersion relation was calculated using imaginary-time correlation functions, and found to be lll in good agreement with experimental results. Finally, we have calculated the Dey be-Waller factor in solid helium for a range of temperatures and densities, and compared the scaling behavior with the predictions of harmonic theory. The first non-Gaussian contribution to the density distribution was calculated.

•Jan 1, 2001•3 citations

A condition for the Bose-Einstein transition in the superfluid liquid helium

V. N. Minasyan

First, a condition for the Bose-Einstein transition in the superfluid liquid helium is presented due to the creation of a free neutron spinless pairs in a liquid helium and a dilute neutron gas mixture. We proposed a new model of dilute Bose gas, where presence atoms in the condensate is a suppressor for the collective modes as well as a creator for single-particle excitations On other hand, it is shown that the terms, of the interaction between the excitations of the Bose gas and the density modes of the neutron, meditate an attractive interaction via the neutron modes, which in turn leads to a bound state on a spinless neutron pair. The lambda transition is defined by a condition for the Bose-Einstein transition, which transforms reflected neutron pair modes to single neutron modes.

arXiv:0804.1220•Apr 8, 2008•0 citations

Surface tension of Bose-Einstein condensate at finite temperature

Nguyễn Văn Thụ

We examine the influence of nonzero temperature on a bounded surface of a dilute Bose gas confined by a hard wall, utilizing the Gross-Pitaevskii theory and the Bogoliubov-de Gennes equation to describe surface excitations. The theoretical calculations are compared with experimental data for liquid helium II, demonstrating excellent agreement. An empirical relation is found for the temperature-dependence of the surface tension. Furthermore, our findings indicate that the contribution of surface excitations remains below 0.7\% in the most of recent experiments on Bose-Einstein condensate.

arXiv:2502.01104•Feb 3, 2025•0 citations

Bose Condensation at He-4 Interfaces

Erik W. Draeger, David M. Ceperley

Path Integral Monte Carlo was used to calculate the Bose-Einstein condensate fraction at the surface of a helium film at T = 0:77 K, as a function of density. Moving from the center of the slab to the surface, the condensate fraction was found to initially increase with decreasing density to a maximum value of 0.9, before decreasing. Long wavelength density correlations were observed in the static structure factor at the surface of the slab. A surface dispersion relation was calculated from imaginary-time density-density correlations. Similar calculations of the superfluid density throughout He-4 droplets doped with linear impurities (HCN)(sub n) are presented. After deriving a local estimator for the superfluid density distribution, we find a decreased superfluid response in the first solvation layer. This effective normal fluid exhibits temperature dependence similar to that of a two-dimensional helium system.

•Jan 1, 2003•0 citations

Ground state properties and elementary excitation of85Rb Bose-Einstein condensate: analytical method for modified Thomas-Fermi model with large-gas-parameter

Moumita Gupta, Krishna Rai Dastidar

Bose-Einstein condensation in trapped gases of 85 Rb atoms can now be obtained in the large gas-parameter region in presence of Feshbach resonance and the gas density which is characterized by peak gas parameter x pk can attain the value ~ 10 -2 . For such large gas-parameter, the first correction term to the Bogoliubov energy functional is considered to give modified Gross-Pitaevskii (MGP) equation. In this work analytical expressions for the chemical potential, the density of particles and the excitation frequency of the monopole mode have been derived in Thomas-Fermi(TF) approximated approach to the MGP equation (modified TF approach) following method of higher order iteration for Bose-Einstein condensate (BEC) of atoms in spherically symmetric trap. We have also solved the MGP equation by using Runge-Kutta method for integration of nonlinear differential equation considering spherically symmetric trap. The analytical results for chemical potential and density of particles for BEC of 10 3 and 10 4 85 Rb atoms in a spherically symmetric trap are in very good agreement with the numerical results from MGP equation at low as well as at large scattering length i.e, at large gas-parameter region. The frequency of monopole mode calculated using the density of particles in the ground state (found out in higher order iterative method) are also in very good agreement with the existing numerical results. The difference between the analytical and numerical results decreases as the value of scattering length i.e, gas-parameter is increased.

Journal of Physics Conference Series•Sep 1, 2007•4 citations
We present arguments that the low density surface region of self-bounded superfluid 4He{}^{4}\mathrm{He}4He systems is an inhomogeneous dilute Bose gas, with almost all of the atoms occupying the same single-particle state at T=0T\phantom{\rule{0ex}{0ex}}=\phantom{\rule{0ex}{0ex}}0T=0. Numerical evidence for this complete Bose-Einstein condensation was first given by the many-body variational calculations of 4He{}^{4}\mathrm{He}4He droplets by Lewart, Pandharipande, and Pieper in 1988 [Phys. Rev. B 37, 4950 (1988)]. We show that the low density surface region can be treated rigorously using a generalized Gross-Pitaevskii equation for the Bose order parameter.
Quantum systems consisting of many particles are both interesting from the physics and the mathematics perspective. One aspect that makes them fascinating from the physics point of view is that they can exhibit quantum effects on macroscopic scales such as Bose-Einstein condensation. This phenomenon, first predicted 100 years ago by Bose [19] and Einstein [44], occurs when a dilute gas of bosons is cooled down to very low temperatures. Below a certain critical temperature a condensate forms, meaning that a macroscopic fraction of the particles is in the same one-particle state, and thus, quantum phenomena become relevant on length scales much larger than the individual atom. The occurrence of Bose-Einstein condensation is also related to superfluidity as suggested in [88]. Bose-Einstein condensation has been achieved experimentally independently by the groups of Cornell and Wieman [3] and Ketterle [40] earning them a Nobel prize in 2001. While Bose and Einstein analyzed theoretically the non-interacting gas, an effective theory describing weakly interacting Bose gases was proposed by Bogoliubov in 1947 with the goal of explaining the occurrence of superfluidity in liquid helium [18]. Using assumptions such as that condensation indeed occurs, his theory yields predictions for the low-energy states of the system. Seen from the mathematics perspective, describing the Bose gas is an interesting challenge due to the strong correlations between particles introduced by the interaction making this case much harder to study than the non-interacting one. In particular, a proof that Bose-Einstein condensation occurs in the thermodynamic limit for continuous systems is still out of reach. However, there has been a lot of work regarding the rigorous analysis of bosonic systems in various relevant regimes. This thesis contains three new contributions establishing rigorous results on bosonic systems in different, physically relevant settings. It is based on the articles [37, 27, 31]. The first result concerns the time-evolution of initially trapped Bose gases in the ultra-dilute Gross-Pitaevskii regime. Under some assumptions on the initial data satisfied in particular by a class of states that are expected to approximate the ground state in a trap, we derive a quasi-free approximation that is close to the true solution in L2L^2L2-norm. In particular, we show that this approximate wave function is described in terms of a Bose-Einstein condensate in a one-particle state following the time-dependent Gross-Pitaevskii equation together with pair excitations with respect to this condensate. As a corollary, we establish a central limit theorem for bounded one-particle observables. The second result is a simple proof of Bose-Einstein condensation for interacting, ultra-dilute gases. More specifically, we consider scalings interpolating between the Gross-Pitaevskii regime and the thermodynamic limit. While condensation had been shown to hold in these regimes before, the advantage of the method used here is its simplicity. Lastly, the third result establishes an upper bound on the ground state energy of a dilute Bose gas in the thermodynamic limit. It captures rigorously the third order correction to the energy per unit volume for sufficiently low densities as predicted in the physics literature first by Wu and subsequently others. Full-text embargoed until: 2026-01-01\textit{Full-text embargoed until: 2026-01-01}Full-text embargoed until: 2026-01-01
Below Tλ = 2.17 K, bulk 4He is a superfluid and has a non-zero Bose-Einstein condensate fraction. This work will focus primarily on how phenomena such as superfluidity, Bose condensation and superfluid vortices are manifested in microscopic, inhomogeneous helium systems. Path Integral Monte Carlo is a powerful method for calculating the equilibrium properties of quantum systems at finite temperature. We have achieved linear scaling of computer time with number of particles through the use of neighbor lists, allowing us to simulate systems of several thousand atoms. We have derived a local superfluid estimator and used it to examine the microscopic superfluid response around a molecule rotating in a helium nanodroplet. We found that the first solvation layer is well-described by a two dimensional superfluid, with thermal excitations occuring at a lower temperature than in bulk helium. The effective moment of inertia of a linear impurity in a helium droplet is calculated, and compared with experimental scattering results. In addition, we calculated the vortex formation energy for both pure droplets and droplets doped with linear impurities, and found that the linear impurities had a negligible impact on the formation energy. A possible spectroscopic signature of vortices in doped helium droplets was suggested. After deriving a local estimator, we calculated the condensate fraction throughout the free helium surface of a semi-infinite slab. These results, along with densitydensity correlation functions, were used to characterize the surface excitations and calculate the extent to which ripplons are present. In addition, the ripplon dispersion relation was calculated using imaginary-time correlation functions, and found to be lll in good agreement with experimental results. Finally, we have calculated the Dey be-Waller factor in solid helium for a range of temperatures and densities, and compared the scaling behavior with the predictions of harmonic theory. The first non-Gaussian contribution to the density distribution was calculated.
First, a condition for the Bose-Einstein transition in the superfluid liquid helium is presented due to the creation of a free neutron spinless pairs in a liquid helium and a dilute neutron gas mixture. We proposed a new model of dilute Bose gas, where presence atoms in the condensate is a suppressor for the collective modes as well as a creator for single-particle excitations On other hand, it is shown that the terms, of the interaction between the excitations of the Bose gas and the density modes of the neutron, meditate an attractive interaction via the neutron modes, which in turn leads to a bound state on a spinless neutron pair. The lambda transition is defined by a condition for the Bose-Einstein transition, which transforms reflected neutron pair modes to single neutron modes.
We examine the influence of nonzero temperature on a bounded surface of a dilute Bose gas confined by a hard wall, utilizing the Gross-Pitaevskii theory and the Bogoliubov-de Gennes equation to describe surface excitations. The theoretical calculations are compared with experimental data for liquid helium II, demonstrating excellent agreement. An empirical relation is found for the temperature-dependence of the surface tension. Furthermore, our findings indicate that the contribution of surface excitations remains below 0.7\% in the most of recent experiments on Bose-Einstein condensate.
Path Integral Monte Carlo was used to calculate the Bose-Einstein condensate fraction at the surface of a helium film at T = 0:77 K, as a function of density. Moving from the center of the slab to the surface, the condensate fraction was found to initially increase with decreasing density to a maximum value of 0.9, before decreasing. Long wavelength density correlations were observed in the static structure factor at the surface of the slab. A surface dispersion relation was calculated from imaginary-time density-density correlations. Similar calculations of the superfluid density throughout He-4 droplets doped with linear impurities (HCN)(sub n) are presented. After deriving a local estimator for the superfluid density distribution, we find a decreased superfluid response in the first solvation layer. This effective normal fluid exhibits temperature dependence similar to that of a two-dimensional helium system.
Bose-Einstein condensation in trapped gases of 85 Rb atoms can now be obtained in the large gas-parameter region in presence of Feshbach resonance and the gas density which is characterized by peak gas parameter x pk can attain the value ~ 10 -2 . For such large gas-parameter, the first correction term to the Bogoliubov energy functional is considered to give modified Gross-Pitaevskii (MGP) equation. In this work analytical expressions for the chemical potential, the density of particles and the excitation frequency of the monopole mode have been derived in Thomas-Fermi(TF) approximated approach to the MGP equation (modified TF approach) following method of higher order iteration for Bose-Einstein condensate (BEC) of atoms in spherically symmetric trap. We have also solved the MGP equation by using Runge-Kutta method for integration of nonlinear differential equation considering spherically symmetric trap. The analytical results for chemical potential and density of particles for BEC of 10 3 and 10 4 85 Rb atoms in a spherically symmetric trap are in very good agreement with the numerical results from MGP equation at low as well as at large scattering length i.e, at large gas-parameter region. The frequency of monopole mode calculated using the density of particles in the ground state (found out in higher order iterative method) are also in very good agreement with the existing numerical results. The difference between the analytical and numerical results decreases as the value of scattering length i.e, gas-parameter is increased.