Thermodynamic behaviour of a one-dimensional Bose gas at low temperature

Thermodynamic behaviour of a one-dimensional Bose gas at low temperature

Giulia De Rosi INO-CNR BEC Center and Dipartimento di Fisica, Università di Trento, Via Sommarive 14, I-38123 Povo, Italy    Grigori E. Astrakharchik Departament de Física, Universitat Politècnica de Catalunya, 08034 Barcelona, Spain    Sandro Stringari INO-CNR BEC Center and Dipartimento di Fisica, Università di Trento, Via Sommarive 14, I-38123 Povo, Italy
July 6, 2019

We show that the chemical potential of a one-dimensional (1D) interacting Bose gas exhibits a non-monotonic temperature dependence which is peculiar of superfluids. The effect is a direct consequence of the phononic nature of the excitation spectrum at large wavelengths exhibited by 1D Bose gases. For low temperatures , we demonstrate that the coefficient in expansion of the chemical potential is entirely defined by the zero-temperature density dependence of the sound velocity. We calculate that coefficient along the crossover between the Bogoliubov weakly-interacting gas and the Tonks-Girardeau gas of impenetrable bosons. Analytic expansions are provided in the asymptotic regimes. The theoretical predictions along the crossover are confirmed by comparison with the exactly solvable Yang-Yang model in which the finite-temperature equation of state is obtained numerically by solving Bethe-ansatz equations. A 1D ring geometry is equivalent to imposing periodic boundary conditions and arising finite-size effects are studied in details. At we calculated various thermodynamic functions, including the inelastic structure factor, as a function of the number of atoms, pointing out the occurrence of important deviations from the thermodynamic limit.

PACS numbers

I Introduction

It is well known that the thermodynamic behaviour of a superfluid is dominated, at low temperature, by the thermal excitation of phonons Wilks1967 (). This explains, in particular, the peculiar behaviour exhibited by the specific heat as well as by other fundamental thermodynamic functions. A non trivial (and less investigated in the literature) consequence of superfluidity shows up in the non-monotonic behaviour of the chemical potential Papoular2012 (). At low temperature the chemical potential increases with as a consequence of the thermal excitation of phonons. At high temperature, in the ideal gas classical regime, the chemical potential is instead a decreasing function of . This non-monotonic behaviour has been recently measured in a strongly interacting atomic Fermi gas Ku2012 (), where it was shown that the chemical potential exhibits a maximum in the vicinity of the superfluid critical temperature.

It is consequently interesting to explore the low-temperature thermodynamic behaviour of other systems, like one-dimensional (1D) interacting Bose gases, which are known to exhibit a phononic excitation spectrum, despite the fact that they cannot be considered superfluids according to standard definition. By investigating the drag flow caused by a moving external perturbation, Astrakharchik and Pitaevskii Astrakharchik2004 () have in fact shown that 1D Bose gases interacting with contact potential exhibit a traditional superfluid behaviour, characterized by the absence of friction force, only in the weakly interaction regime, where Bogoliubov theory applies and the gas can be locally considered Bose-Einstein condensed, despite the absence of true long range order.

In this work, we investigate the low-temperature expansion of the chemical potential of a 1D Bose gas with contact repulsive interaction for the whole crossover, ranging from the weakly to the strongly interaction limits. A major motivation is given by the possibility of comparing the low-T expansion of the chemical potential with the numerical results now available within the Yang-Yang theory Yang1969 (); Yang1970 (), along the whole interaction strength crossover. Previous comparisons were in fact available only in the case of the Tonks-Girardeau limit Lang2015 (), corresponding to the ideal Fermi gas, where the low-T expansion corresponds to the Sommerfeld expansion. We find that for all intermediate interaction regimes, described at by Lieb-Liniger (LL) theory, the increase of the chemical potential at low temperature follows the law and is actually caused by the phononic nature of the long wavelength elementary excitations, as in usual superfluids Papoular2012 (). The relevant coefficient fixing the law depends on the density derivative of the sound velocity which is calculated using Lieb-Liniger theory. This feature strengthens the analogy with superfluids even in 1D dimension. Importantly, our results can be also generalized to every Luttinger liquid at low temperature whose macroscopic elementary excitations can be described in terms of phonons.

Recently, a ring geometry has been experimentally realized for a microscopic system of atoms Labuhn2016 (). Motivated by the experimental progress, we study in details also the behavior of a gas containing a finite number of atoms in a ring, focusing on the deviations of its thermodynamic behavior from the one in the large limit.

Our system is a uniform gas of bosons interacting with a repulsive contact interaction


where the interaction parameter is related to the 1D coupling constant through , where is the 1D scattering length. The system (1) has been realized experimentally for the whole interaction crossover by suitably tuning the interaction strength Paredes2004 (); Kinoshita2006 (); Cazalilla2011 (), described by the dimensionless parameter


from weak () to strong () interactions Kinoshita2006 (); Haller2009 (); Haller2010 (); Haller2011 (); Guarrera2012 (). The Bogoliubov (BG) perturbative theory can be used in the limit of weak interactions. In the Tonks-Girardeau (TG) limit of strong repulsions the bosons are impenetrable and their wave function can be mapped onto that of an ideal Fermi gas Girardeau1960 ().

The paper is organized as follows.

In Sec. II we derive the low-temperature expansion of the chemical potential, starting from the free energy of an ideal phononic gas. This assumption is fully justified by the low-momenta behavior of the Lieb-Liniger excitation spectrum. The low-temperature expansion exhibits a -dependence on temperature, with the coefficient related to the density derivative of the LL sound velocity at zero temperature. The Bethe-ansatz results for the chemical potential are shown to agree very well with the low-temperature expansion, for the whole BG-TG crossover.

In Sec. III we investigate the BG weakly-interacting gas. By considering the quantum fluctuation contribution in the ground-state energy at , we explore the behavior of the chemical potential and of the sound velocity. While this correction is important at , it does not affect the low-temperature expansion of the chemical potential.

Similarly to Sec. III, we calculate in Sec. IV the first corrections in the interaction parameter to the TG strongly interacting gas. The starting point is the expansion, for large values of , of the ground-state energy of a hard-sphere gas.

In Sec. V we derive the low-temperature expansions of both the adiabatic and the isothermal inverse compressibilities. The coefficients of the laws are studied as a function of the interaction parameter and analytically calculated in the BG and TG limits.

In Sec. VI we consider a ring configuration with a finite number of particles at zero temperature and calculate the finite-size corrections with respect to the thermodynamic limit for the energy, the chemical potential and the sound velocity. Results for the static inelastic structure factor for a finite number of particles are also reported.

In Sec. VII, we draw our final conclusions.

Ii Low-temperature expansion of the chemical potential

It is well known that at the elementary excitations of an interacting 1D Bose gas have a phononic character at small momenta Pitaevskii2016 (); Lieb1963 (), characterized by the linear dispersion relation


At the sound velocity is related to the density dependence of the chemical potential according to the relation


where is the chemical potential and denotes the linear density. The density dependence of the chemical potential at zero temperature can be calculated within the Lieb–Liniger model. The ratio between the sound velocity and the Fermi velocity is known as the Luttinger parameter, , and it plays an important role in defining the long-range properties of one-dimensional systems. Figure 1 shows the dependence of the sound velocity on the interaction parameter for the Lieb-Liniger model, described by Hamiltonian (1). There is a smooth crossover between the mean-field BG value defined as for weak interactions to the Tonks-Girardeau (ideal Fermi gas) value in the limit of strong repulsion.

Figure 1: (Color online) Sound velocity in units of Fermi velocity (solid line) as a function of the interaction parameter , calculated by solving the Lieb-Liniger equations. The Bogoliubov (dotted line, ) and Tonks-Girardeau (dashed line, ) limits, including their first-order corrections (thin solid line and thin dot-dashed line , respectively) are present too, see Secs. III and IV.

For larger momenta the 1D excitation spectrum is characterized by a continuous structure, bounded by two branches of elementary excitations Lieb1963 (); Yang1969 (); Pitaevskii2016 (), which have been the object of recent measurements Meinert2015 (); Fabbri2015 (). For small values of , the Lieb-I particle-like branch corresponds to the Bogoliubov excitation spectrum Lieb1963 (); Kulish1976 (); Pitaevskii2016 (). The Lieb-II hole-like branch is instead associated in the weakly-interacting regime with the dark soliton dispersion predicted by Gross-Pitaevskii theory Pitaevskii2016 (); Ishikawa1980 (); Kulish1976 (). The two branches merge into the phononic spectrum for , Fig. 2.

Figure 2: (Color online) Lieb–Liniger excitation spectrum in the BG regime with (left) and in the deep TG regime with (right). The units are the Fermi energy and the Fermi momentum . The shaded region represents the continuum of the excitations and is delimited by the upper (Lieb I) and the lower (Lieb II) branch of the spectrum. On the left, the Lieb I and II branches are not reported. On the left, the dashed line gives the Bogoliubov dispersion and the dotted line gives the mean field soliton spectrum. In the limit , the Lieb I branch tends to be equal to the Bogoliubov dispersion, while the Lieb II one coincides with the soliton spectrum. The solid line is the Lieb–Liniger phononic spectrum calculated with . On the right, Lieb I and Lieb II branches are reported and they coincide with the particle and hole ideal Fermi gas excitations, respectively. The solid line is the phononic spectrum calculated with the Fermi velocity.

At low temperature () we expect that the thermodynamic behaviour of the system can be calculated in terms of a gas of non interacting phonons. The free energy of this gas is then given by


where is dispersion (3) and we have added the energy calculated at with the Lieb–Liniger theory. Notice that the thermal contribution to is affected by two-body interactions through the dependence of on the interaction parameter . The integral of Eq. (5) yields the following low-temperature expansion for the free energy


which differs from the usual behaviour exhibited by three-dimensional () superfluids Pitaevskii2016 () because of the 1D structure of the integral (5). Starting from result (6) for the free energy, one can calculate the low-temperature expansion of the chemical potential:


where we have introduced the energy scale given by the Fermi energy of a 1D Fermi gas, because it exhibits the same density dependence of the quantum degeneracy temperature of the system. We have also defined the relevant dimensionless parameters of the expansion




which are functions of the interaction parameter and can be calculated at zero temperature using Lieb–Liniger theory. It is worth noticing that the parameter , which is the most relevant because it fixes the leading coefficient of the low- expansion, depends on the density derivative of the sound velocity. The two numerical functions , Eq. (8), and , Eq. (9), have been calculated within LL theory and their values are reported in Figs. 3 and 4 with their BG and TG limits. In particular, the TG limits for and reproduce the low-temperature Sommerfeld expansion of the chemical potential for the 1D ideal Fermi gas, Eq. (21).

Figure 3: (Color online) (solid line) with the leading dependence (dashed line, and dotted line, ) and first order corrections [thin dot-dashed line, and thin solid line, ] for Tonks-Girardeau and Bogoliubov limits, respectively.
Figure 4: (Color online) (solid line) with the Tonks-Girardeau (dashed line, ) and Bogoliubov (dotted line, ) limits.

In Figs. 5-6 we report the temperature dependence of the chemical potential of the system described by Hamiltonian (1) as obtained numerically from the Bethe-ansatz (BA) approach first developed by Yang-Yang Yang1969 (); Yang1970 (); Kheruntsyan2005 (); Lang2015 () for several characteristic values of . The Yang-Yang description has been probed experimentally vanAmerongen2008 (); Vogler2013 () and allows not only to investigate the thermodynamics, but also the Luttinger liquid physics and the quantum criticality of the system Guan2011 (); Wang2013 (); Guan2014 (). The numerical results for the thermodynamics have been derived recently in an analytic fashion by using the polylog functions at finite temperature by Guan and Batchelor Guan2011 (); Guan2014 ().

The crossover from mean-field to Tonks-Girardeau regimes (see Fig. 1) introduces two distinct energy scales. Correspondingly, we rescale the chemical potential in units of the Fermi energy in Fig. 5 and in units of the mean-field zero-temperature chemical potential in Fig. 6. The first choice provides natural units in the TG regime in which strongly repulsive bosons behave similarly to an ideal Fermi gas (IFG) in the limit of . In this regime, the chemical potential as a function of is calculated by inverting the Fermi–Dirac distribution (upper dashed line in Fig. 5):


and, despite the absence of superfluidity, it still exhibits the quadratic low-temperature dependence , which follows from the low-temperature Sommerfeld expansion, Eq. (21).

By reducing the interaction parameter , the system becomes softer and the limit of vanishing interactions, , corresponds to an ideal Bose gas (IBG) with the chemical potential fixed by the relationship (lower dashed line in Fig. 5):


Notice that, because of the absence of Bose-Einstein condensation Mermin1966 (); Hohenberg1967 (), the chemical potential of the 1D ideal Bose gas is always negative and approaches the value as . Remarkably, for all finite interaction strengths the temperature dependence is not monotonic. Moreover, the initial increase is perfectly described by the quadratic low-temperature expansion (7), thereby proving that the model based on a gas of independent phonons well accounts for the thermodynamic behaviour of the 1D interacting Bose gas. This is a non trivial result due to the complex structure of the elementary excitations at larger wave vectors exhibiting a double branch converging into the phonon law (3) only at small momenta. We notice also that the chemical potential for high temperatures, which is a decreasing function of , can be considered as a shift of the ideal Bose chemical potential, Eq. (11), for every value of .

The behavior of the chemical potential in the weakly-interacting regime () is best seen in Fig. 6. For low temperatures the gas behaves like a quasicondensate, exhibiting typical features of superfluids. For , the gas is a thermal degenerate gas, while for the gas behaves classically with . A similar classification of the quantum degeneracy states in 1D trapped configurations was first proposed in Petrov2000 ().

Figure 5: (Color online) Chemical potential as a function of temperature in Fermi units for several values of and at a fixed density . The solid lines represent the Bethe–ansatz (BA) solutions for different values of . The dot-dashed lines are the low–temperature expansions of the chemical potential taking into account only the phononic contribution, Eq. (7). The phononic expansions for are equal to the analytical Sommerfeld expansion of Eq. (21). Both the chemical potentials as a function of for the ideal Fermi (upper dashed line) and ideal Bose (lower dashed line) gas are also reported, Eq. (10) and (11), respectively.

Although there is no phase transition in 1D systems at finite , in the canonical ensemble, there exists a critical point, corresponding to the value of the chemical potential, which separates the vacuum from the filled “Fermi sea” of repulsive bosons at . In particular, a universality class is present in the temperature regime and near the critical point Guan2011 (); Wang2013 (); Guan2014 ().

Figure 6 is similar to Fig. 5, but with the chemical potential expressed in units of the BG chemical potential at zero temperature: and the temperature in units of:


which has been introduced as an appropriate temperature scale for visualizing the behavior of the chemical potential at low temperature. With the new units, the phononic expansion (7) takes the form:

Figure 6: (Color online) Chemical potential as a function of temperature in BG units for several values of . The solid lines represent the Bethe–ansatz (BA) solutions for different values of . The dashed lines are the low–temperature expansions of the chemical potential taking into account only the phononic contribution, Eq. (13). We notice also that the phononic expansion (13) does not hold for very small , like , for which value the low–temperature expansion is not reported.

Figures 5 and 6 point out in a clear way the non-monotonic behavior of the chemical potential as a function of for a fixed value of the density. This is a general feature exhibited by superfluids Papoular2012 () and it is shown here that it characterizes also interacting 1D Bose gas for all finite values of the interaction parameter .

Both figures show also that the phononic expansion describes very well the low-T thermodynamics for all values of , although the region of the applicability of the phononic description depends on . As pointed out in Ref. Lang2017 (), for small values of the interaction parameter , higher-order corrections beyond the linear phononic contribution in the excitation spectrum (3) might be important.

Iii Bogoliubov regime

In the mean-field theory, the chemical potential is linear in density, and the velocity of sound takes the value , see Fig. 1.

The first correction to the mean-field expression for the equation of state comes from the quantum fluctuations Lee1957 (); Lee19572 (); Pitaevskii2016 (). With respect to the 3D case, in 1D this calculation is simpler because it does not require the renormalization of the scattering length due to the absence of ultraviolet divergencies in the calculation of the ground-state energy. Therefore in 1D one can consider all ranges of momenta and one finds Lieb1963 ():




is the Bogoliubov excitation spectrum. By considering the thermodynamic limit of Eq. (14) and by solving the integral in momentum space, one finally finds the first-order correction in the interaction parameter for the ground state energy Lieb1963 ()


The same result can be also found by performing a power series expansion of the Lieb-Liniger equations Kaminaka2011 (); Gudyma20152 (). The correction is negative as it comes from second order perturbation theory and, contrary to the higher dimensions, in 1D there is no renormalization of the coupling constant thus no additional terms have to be added.

Equation (16) allows one to calculate the higher-order corrections for the other thermodynamic quantities at . For the chemical potential, one finds the result


which implies the result


for the expansion of the coefficient , where is the mean-field value. The corresponding result has been plotted in Fig. 3 and well reproduces the exact value of up to values .

From Eq. (4) and Eq. (17), one can calculate also the correction to the sound velocity Lieb1963 ()


which is also reported in Fig. 1, yielding the expression


for the coefficient , Eq. (9) with the Bogoliubov value. Notice that, differently from the case of [see Eq. (18)], the first correction vanishes because of an exact cancellation between the corrections provided by the terms and of Eq. (9). This explains why the Bogoliubov approximation describes correctly the value of for a large interval of values of , up to (see Fig. 4).

Iv Tonks-Girardeau regime

In the TG limit of strong repulsion, , the energetic properties are the same as in an ideal Fermi gas. The thermodynamic quantities do not depend on the coupling constant , but only on the density , encoded in the Fermi energy . This regime can be interpreted as that of a unitary Bose gas with the Bertsch parameter equal to 1 as the chemical potential is equal to the Fermi energy []. Similarly, the sound velocity is equal to the Fermi velocity , see Fig. 1. The low-temperature expansion of the chemical potential in this limit is equal to the first terms of the Sommerfeld expansion (dot-dashed line for in Fig. 5) of the 1D ideal Fermi gas, as already pointed out in Lang2015 ():


which contains the TG limits of and parameters, Figs. 3 and 4.

Leading corrections to the ground-state energy in the TG regime arise from the “excluded volume” and can be obtained from the equation of state of hard spheres (i. e. impenetrable) bosons with diameter  Girardeau1960 ():


In the limit of point-like bosons , Eq. (22) reproduces the ground state energy of the ideal Fermi gas, . Expanding the denominator in Eq. (22) generates a power series with integer coefficients, . It is interesting to notice that for a -interacting potential the momentum-dependent -wave scattering length, , does not affect first and second corrections in but induces a negative correction in front of the third correction. Indeed, the universality of the first and the second corrections becomes evident by comparing low-density expansion of the equation of state for hard spheres, Eq. (22), and contact -potential obtained by solving Bethe equations recursively Astrakharchik2010 (),


The non-universal correction depends on the shape of the potential and for the LL model it has a non-integer coefficient, which qualitatively can be understood by noting that the typical value of the scattering momentum in TG regime is proportional to , which is consistent with terms appearing in expansion (23). The universal terms are the same both in the super Tonks-Girardeau Astrakharchik2005 () () and the strongly repulsive Gudyma20152 () () regimes. From Eq. (23), by introducing the parameter (2) and by considering only the leading term, one finds


From Eq. (24), one easily calculates the correction of the chemical potential at , Eq. (8):


which implies the result Gudyma20152 ()


for , including the first correction to the TG result . Prediction (26) is reported in Fig. 3 for positive values of , its accuracy being good for values of larger than .

From Eq. (4) and Eq. (25), one can calculate also the first correction, at large , to the sound velocity Gudyma20152 (); Valiente2016 ():


which is reported in Fig. 1. For the coefficient , which provides the -correction in the expansion of the chemical potential, we find again an exact cancellation between the correction (provided by the term and entering the expression (9) for , similarly to what happens in the small expansion discussed in the previous Section III in the case of the Bogoliubov gas. We then find that the Tonks-Girardeau expression provides an accurate estimate of for values of larger than (see Fig. 4).

More accurate analytical expressions for the above thermodynamical quantities, which allow to probe the whole range of interaction strength with excellent accuracy, are reported in Ristivojevic2014 (); Jiang2015 (); Prolhac2017 (); Lang2017 ().

V Low-temperature expansion of the inverse compressibility

Here we derive the dependence of the adiabatic and isothermal inverse compressibilities on the interaction parameter in the limit of low temperature.

v.1 Adiabatic inverse compressibility and sound velocity

From the Gibbs–Duhem relation , one finds


where is the entropy density and is the entropy per particle.

At low temperature the entropy per particle of a non–interacting gas of phonons takes the form Pitaevskii2016 ()


which depends on the value (4) of the sound velocity. Use of relation (29) permits to express the dependence of the second contribution to the adiabatic inverse compressibility on the r.h.s. of Eq. (28) on the interaction parameter


in terms of the coefficient , Eq. (9) related to the density derivative of the sound velocity at constant entropy. The first contribution on the r.h.s. of Eq. (28) can be obtained by using Eqs. (7) and (29),


From the above equations one finally finds the low temperature expansion


of the adiabatic inverse compressibility, where


is its value and we have defined the positive quantity


which is reported in Fig. 7 together with its asymptotic limits in the Bogoliubov and Tonks–Girardeau regimes.

Figure 7: (Color online) Value of the dimensionless coefficient of the low–temperature expansion of the adiabatic inverse compressibility (solid line). The BG and TG analytical limits are also shown: (dotted line) and (dashed line).

v.2 Isothermal inverse compressibility

By fixing the temperature in Eq. (28) and by considering the low-temperature expansion of the chemical potential (7), one can also calculate the low-temperature expression for the isothermal inverse compressibility


where we have defined the negative dimensionless coefficient


Notice that the thermal corrections to the isothermal and adiabatic inverse compressibilities have opposite sign, being the coefficient always negative. The absolute value of is reported in Fig. 8 together with the asymptotic limits in the Bogoliubov and Tonks–Girardeau regimes. The negative value of is the consequence of the peculiar temperature dependence of the free energy (6).

Figure 8: (Color online) Absolute value of the dimensionless coefficient of the low–temperature expansion of the isothermal inverse compressibility (solid line). The BG and TG analytical limits are also shown: (dotted line) and (dashed line).

Vi Gas on a ring

The physics in one dimension is unusual in many aspects. The mean–field regime is reached at large densities contrarily to what happens in three dimensions where the weakly–interacting limit corresponds to small densities, according to the limit . For a fixed number of particles the mean–field limit in one dimension, , can be obtained either increasing the linear density , by decreasing the system size , or by increasing the –wave scattering length , i.e. decreasing the coupling constant . Asymptotically, at a certain point, the size of the system will become comparable to the healing length


and finite–size effects will become important. This should be contrasted to the three-dimensional case where the mean–field regime is instead achieved by increasing the system size which consequently becomes larger than the healing length.

Finite–size effects depend on the system geometry. Interestingly, periodic boundary conditions, commonly used as a mathematical tool in the three–dimensional world, in one dimension can be explicitly realized in a ring and have consequently a direct physical interest. This is another peculiarity of the one–dimensional world. In the following we calculate the finite–size dependence of thermodynamic quantities for a gas confined in a ring whose properties are then equivalent to the ones of a linear 1D system satisfying periodic boundary conditions (PBC). If one considers a plane wave and one imposes PBC, one finds that the momentum is quantized according to


where are integers. Moreover, in 1D, all the integrals in momentum space, defined in the thermodynamic limit (, finite), are replaced by a sum over the discretized momenta (38) as:


In the following, we calculate the finite-size corrections in both BG and TG regimes at zero temperature, as well as the static inelastic structure factor for a finite number of particles.

vi.1 Bogoliubov regime at

Let us consider the ground-state energy per particle given by


corresponding to the Bogoliubov regime of small , where is provided by the Bogoliubov spectrum (15). Equation (40) differs from Eq. (14) because it contains the term in the sum. This term has been included in order to avoid self-interaction effects in the leading mean-field term of Eq. (14) which should be replaced by .

By introducing the discretized values of  (38), the energy can be rewritten in the form


where we have introduced the dimensionless variable


depending on the interaction parameter and the function


where the adding of the quantity ensures that the term in the sum is counted just once.

By using the Euler-Maclaurin expansion (see Appendix A), one can calculate the expression for the series (43) for large values of :


In Fig. 9 we report the comparison of the series (43) with its expansion (44). We notice that the two curves agree in an excellent way for . The thermodynamic limit is also reported.

Figure 9: (Color online) Comparison of the numerical series  (43) (solid line) and its analytical expansion (44) (dashed line) holding for . The dot-dashed line represents the thermodynamic value.

For large number of particles, the ground-state energy per particle (41) then takes the form:


and, in the thermodynamic limit, reproduces Eq. (16). The condition is equivalent to requiring that the healing length (37) be smaller than the size of the system.

The ground-state energy contains three contributions: the leading term corresponds to the usual mean field energy, the second contribution arises from the quantum fluctuations and is a one-dimensional analog of the Lee-Huang-Yang correction in 3D, while the last term accounts for finite-size effects and depends explicitly on the interaction parameter .

Finite size corrections can be sizeable, as clearly shown by Fig. 10 where we report the energy per particle as a function of for the thermodynamic limit (16) (dot-dashed line), the Bethe-ansatz (BA) calculation (circle), the Bogoliubov expression (41) (solid line) and the expansion (45) (dashed line). The figure reveals a general good agreement between the BA and the Bogoliubov predictions (41), except for , where Eq. (41), being based on the Bogoliubov approach, is no longer adequate.

Figure 10: (Color online) Comparison of the ground state energy per particle, in BG units, as a function of in the thermodynamic limit of Bogoliubov theory (16) (dot-dashed line), the Bethe-ansatz (BA) calculation (circle), the Bogoliubov expression (41) (solid line) and the expansion (45) (dashed line), for several values of the interaction parameter .

The chemical potential can be obtained by deriving Eq. (40) with respect to , at fixed . One finds


which can be rewritten as , where is provided by Eq. (42) and we have introduced the series


depending on the quantized momenta (38) and such that the zero-momentum term is accounted for once. The Euler-Maclaurin expression, applied to the sum (47), yields


holding in the limit. In Fig. 11 we report the comparison of the series (47) with its expansion (48) holding for . The two curves agree very well for .

Figure 11: (Color online) Comparison of the numerical series  (47) (solid line) and its analytical expansion (48) (dashed line) holding for . The dot-dashed line represents the thermodynamic value.

Using Eq. (48), one can finally write the following expansion for the chemical potential


In Fig. 12 we report the results for the chemical potential as a function of (42) for the thermodynamic limit (17) (dot-dashed line), the Bethe-ansatz calculation (symbols) and the Bogoliubov expression (46) (solid line). The expansion (49) practically coincides with the full series (46). The square symbol corresponds to the “forward” definition of the chemical potential, the star symbol to the “backward” expression , while the circles to the “symmetric” value . While the three definitions of the chemical potential coincide in the thermodynamic limit , they are different in a finite system  nuclear (). In particular, the symmetric definition well agrees with the calculation (46), based on the differential definition , except for the case.

Figure 12: (Color online) Comparison of the chemical potential in BG units as a function of in the thermodynamic limit (17) (dot-dashed line), the Bethe-ansatz (BA) calculation (symbols) and the Bogoliubov expression (46) (solid line), for several values of the interaction parameter . For the BA: (square), (star) and (circle).

From Eq. (46), one can also calculate the sound velocity (4), corresponding to the density derivative of the chemical potential for a fixed value of . The resulting expression,


with the sound velocity defined in the Bogoliubov regime, used in Fig. 1. The above expression can be rewritten as where we have defined the series


after introducing the variable  (42) and the quantized momenta (38). As before, we apply the Euler-Maclaurin formula and we find the expansion


holding in the limit, yielding the asymptotic expansion


for the sound velocity.

vi.2 Tonks–Girardeau regime at

According to Girardeau Girardeau1960 (), the ground-state energy of the gas in the strongly-interacting limit is the same as that of an ideal Fermi gas. The energy for a finite number of particles in a box with periodic boundary conditions is obtained by summing the energy of the single-particle levels in the box,


In the thermodynamic limit, , Eq. (54) results in . The “excluded volume” correction should be present for a finite interaction strength, see the hard–sphere like expression, Eq. (22), and the discussion below it. In order to incorporate the leading finite-size correction close to the Tonks-Girardeau regime we replace with in Eq. (54) resulting in the following expression for the energy per particle


For large values of the interaction parameter one can replace the factor with . In Fig. 13 we report the energy per particle as a function of for the TG regime (54) (solid line), the hard-sphere (HS) like model (55) (dashed and dotted lines) and the Bethe-ansatz solution (symbols) for several values of . We observe a very good agreement between the BA solution and the analytical hard-sphere (55) expression. For the BA results are indistinguishable from the TG limit (54) and they are not reported in the figure. The comparison between Eq. (55) and Eq. (45) reveals that finite-size effects are less important in the TG regime since in the weakly interacting Bogoliubov regime, the correction is amplified by the smallness of .

Figure 13: (Color online) Energy per particle in units of the TG gas energy as a function of . Bethe–ansatz (BA) results (symbols) with different values of the interaction parameter are compared with the TG gas (54) (solid line) and the hard–sphere (HS) model (55) (dashed and dotted lines).

For strong repulsion we obtain the finite-size correction to the chemical potential


and to the sound velocity (4):


It is interesting to note that while the finite-size correction to the energy (55) and the chemical potential (56) scales as with the number of particles, such a correction is instead asymptotically vanishing in the sound velocity (57).

In Fig. 14, we plot the chemical potential with given by Eq. (55) (solid line) as a function of for different values of  muTGring (). In the same figure we plot also the values of and which differ from the symmetric value for small values of nuclear (), similarly to the case of the weakly interacting Bose gas. Differently from the weakly interacting BG gas, the symmetric value however exhibits significant deviations with respect to the differential estimate , for small values of .

Figure 14: (Color online) Chemical potential at in Fermi units as a function of the number of particles in the TG regime (63) for fixed values of (solid line). The dashed lines correspond to the thermodynamic limit of the TG model (64). The symbols correspond to the Bethe-ansatz (BA) calculation: (square), (star) and (circle).

vi.3 Static inelastic structure factor

The ring geometry has a profound effect on the correlation functions. Here we analyze the inelastic static structure factor at zero temperature,


where is the density operator in momentum representation. The static structure factor gives information about two-body correlations and can be measured in experiments by means of Bragg spectroscopy.

In the thermodynamic limit, the static structure factor has a linear behavior at small momenta, , with the slope determined by the sound velocity . The ring geometry introduces both discretization in the allowed momentum and a change in the slope due to the finite-size correction to the sound velocity. The latter effect is rather weak, especially in the Tonks-Girardeau regime, but is important in the context of the finite-size dependence of the Luttinger parameter .

The strongest effect comes from the discretization of the allowed momenta on a ring. For the standing wave values (38), the last term in Eq. (58), corresponding to the square of the so-called elastic form factor, does not contribute. Indeed, one finds , which exactly vanishes for .

Figure 15: (Color online) Static structure factor at in the Tonks-Girardeau limit for different number of particles (solid lines). Values of momenta (38) corresponding to standing waves on a ring are shown with circles. The dashed line represents the phononic law , which coincides with the static structure factor in the thermodynamic limit for in the TG regime.

Figure 15 reports the static structure factor in the Tonks-Girardeau regime. When the probing momentum is equal to a standing wave value (38) in the ring, the value of the static structure factor is exactly the same as in the thermodynamic limit. In this way, discrete points form a linear phononic dependence. As the number of particles is increased, the phononic behavior is better resolved. The absence of the change of the slope means that the finite-size corrections to the sound velocity are negligible in the Tonks-Girardeau regime, confirming the predictions of Eq. (57).

When the probing momentum is different from the allowed values in the ring, the value of depends strongly on the number of particles. Importantly, the small-momentum behavior is no longer linear but rather shows a quadratic dependence on . This qualitative change reflects the change in the structure of the excitation spectrum which becomes discrete. A quadratic dependence on the momentum, , is typical to gapped systems with being the value of the gap. In the discrete case it is not possible to create an excitation with energy smaller than , resulting in a quadratic low-momentum dependence. In the thermodynamic limit and the phononic linear behavior is restored.

Figure 16: (Color online) Static structure factor at and for different number of particles (solid lines). Values of momenta (38) corresponding to standing waves on a ring are shown with circles.

In Fig. 16 we show the static structure factor for , calculated using the diffusion Monte Carlo method. Similarly, the finite-size quadratic behavior at small momenta is replaced by the linear phononic dependence in the thermodynamic limit. Contrarily to the TG case, here the values at depend on the number of particles, although the effect is weak (see, for example, the value at ). In terms of the Luttinger parameter, which in the linear regime corresponds to , this results in its finite-size dependence.

While for the TG regime, the linear dependence extends up to , for weaker interactions the linear regime shrinks (compare Figs. 15-16). Eventually for the linear regime becomes very small and phononic theory cannot provide a good description of the system properties. A similar effect was observed in Figs. 5-6 in the applicability of the phononic theory in the limit of weak interactions.

Vii Conclusion

In this paper we have investigated the low temperature properties of 1D Bose gases along the whole Bogoliubov (BG) — Tonks-Girardeau (TG) crossover. We have shown that, at low temperature, the chemical potential exhibits a typical behavior, which follows from the leading contribution to thermodynamics arising from the thermal excitation of phonons, similarly to what happens in superfluids. The chemical potential is always a decreasing function of at high temperature, thus the increase exhibited by the chemical potential at low temperature is responsible for a typical non-monotonic behavior as a function of . The coefficient of the law has been calculated using the Lieb-Liniger results for the sound velocity and the resulting behavior has been successfully compared with thermodynamic functions obtained from the Yang-Yang theory of 1D interacting Bose gases. We have also presented results for the temperature dependence of the isothermal and adiabatic inverse compressibilities. In particular we have shown that the correction has opposite sign in the two cases.

In the second part of the paper we have focused on the corrections to the thermodynamic functions caused by the finite size of the system. To this purpose, we have considered the useful ring geometry and the mapping with the 1D problem where calculations are carried out using periodic boundary conditions. Explicit results have been obtained in the weakly and strongly interacting regimes where, at zero temperature, the first corrections to the thermodynamic limit, due to finite size effects, can be calculated in analytic form, in excellent agreement with the numerical results provided by the Bethe-ansatz. We have found that finite-size corrections are particularly important in the weakly interacting regime where the healing length can easily become comparable to the size of the system.

Concerning future developments of the analysis carried out in this paper, it is worth mentioning the physical understanding of higher-order corrections (beyond the -law caused by the real excitations of the phononic branch) to the low-temperature thermodynamic behavior. In particular, it is important to understand the temperature corrections arising due to non-symmetric spreading of the phononic branch (different beyond-linear behavior of the lower and upper branches) as well as effects originating from the non-linear behavior of the Bogoliubov spectrum at large momenta. A further perspective of research concerns the finite temperature thermodynamic behavior of 1D Bose gases containing a small number of atoms and confined in a ring of finite size.

Appendix A Euler-Maclaurin expansion for G(y)

In this Appendix, we show the detailed derivation of the expansion holding for  (44) for the series (43).

We use the Euler-Maclaurin expansion which allows to approximate a series as follows Abramowitz2012 ():


where is a continuous function of real numbers in the interval . For , one considers only the first terms in the sum, whose Bernoulli’s numbers are


and are the –derivatives of the function .

By defining the following function


entering the series (43), one estimates the integral


which allows to calculate the thermodynamic limit of the ground-state energy per particle on a ring configuration (41), provided by Eq. (16).

By calculating the first derivative of the function (61), and by using Eq. (59) and Eq. (60), one finally gets the expansion (44) holding for large values of the parameter.

G. De Rosi and S. Stringari would like to acknowledge fruitful and helpful discussions with L. P. Pitaevskii, C. Menotti, S. Giorgini, M. Di Liberto and G. Bertaina. This work has been supported by ERC through the QGBE grant, by the QUIC grant of the Horizon2020 FET program and by Provincia Autonoma di Trento (De Rosi & Stringari). G. De Rosi acknowledges the hospitality of the Computer Simulation in Condensed Matter Research Group (SIMCON) of Universitat Politècnica de Catalunya in Barcelona, where this work was partially done. G. E. Astrakharchik acknowledges partial financial support from the MICINN (Spain) Grant No. FIS2014-56257-C2-1-P. The Barcelona Supercomputing Center (The Spanish National Supercomputing Center – Centro Nacional de Supercomputación) is acknowledged for the provided computational facilities. The authors gratefully acknowledge the Gauss Centre for Supercomputing e.V. ( for funding this project by providing computing time on the GCS Supercomputer SuperMUC at Leibniz Supercomputing Centre (LRZ, The authors would like to acknowledge also G. Lang, X.-W. Guan and the referees of this manuscript for useful comments and suggestions which have allowed some improvements in the revised version of this work.


  • (1) J. Wilks, ”Properties of Liquid and Solid Helium (Monographs on Physics)”, (Clarendon Press, Oxford, 1967).
  • (2) D. J. Papoular, G. Ferrari, L. P. Pitaevskii, and S. Stringari, ”Increasing Quantum Degeneracy by Heating a Superfluid”, Phys. Rev. Lett. 109, 084501 (2012).
  • (3) M. J. H. Ku, A. T. Sommer, L. W. Cheuk, M. W. Zwierlein, ”Revealing the Superfluid Lambda Transition in the Universal Thermodynamics of a Unitary Fermi Gas”, Science 335, 563 (2012).
  • (4) G. E. Astrakharchik and L. P. Pitaevskii, ”Motion of a heavy impurity through a Bose-Einstein condensate”, Phys. Rev. A 70, 013608 (2004).
  • (5) C. N. Yang and C. P. Yang, ”Thermodynamics of a One-Dimensional System of Bosons with Repulsive Delta‐Function Interaction”, J. Math. Phys. 10, 1115 (1969).
  • (6) C. P. Yang, ”One-Dimensional System of Bosons with Repulsive -Function Interactions at a Finite Temperature T”, Phys. Rev. A 2, 154 (1970).
  • (7) G. Lang, F. Hekking, and A. Minguzzi, ”Dynamic structure factor and drag force in a one-dimensional strongly interacting Bose gas at finite temperature”, Phys. Rev. A 91, 063619 (2015).
  • (8) H. Labuhn, D. Barredo, S. Ravets, S. de Léséleuc, T. Macrì, T. Lahaye and A. Browaeys, ”Tunable two-dimensional arrays of single Rydberg atoms for realizing quantum Ising models”, Nature 534, 667 (2016).
  • (9) B. Paredes, A. Widera, V. Murg, O. Mandel, S. Fölling, I. Cirac, G. V. Shlyapnikov, T. W. Hänsch and I. Bloch, ”Tonks-Girardeau gas of ultracold atoms in an optical lattice”, Nature 429, 277 (2004).
  • (10) T. Kinoshita, T. Wenger, D. S. Weiss, ”Observation of a One-Dimensional Tonks-Girardeau Gas”, Science 305, 1125 (2004);
    ”A quantum Newton’s cradle”, Nature (London) 440, 900 (2006).
  • (11) M. A. Cazalilla, R. Citro, T. Giamarchi, E. Orignac, and M. Rigol, ”One dimensional bosons: From condensed matter systems to ultracold gases”, Rev. Mod. Phys. 83, 1405 (2011).
  • (12) E. Haller, M. Gustavsson, M. J. Mark, J. G. Danzl, R. Hart, G. Pupillo, H.-C. Nägerl, ”Realization of an Excited, Strongly Correlated Quantum Gas Phase”, Science 325, 1224 (2009).
  • (13) E. Haller, M. J. Mark, R. Hart, J. G. Danzl, L. Reichsöllner, V. Melezhik, P. Schmelcher, and H.-C. Nägerl, ”Confinement-Induced Resonances in Low-Dimensional Quantum Systems”, Phys. Rev. Lett. 104, 153203 (2010).
  • (14) E. Haller, M. Rabie, M. J. Mark, J. G. Danzl, R. Hart, K. Lauber, G. Pupillo, and H.-C. Nägerl, ”Three-Body Correlation Functions and Recombination Rates for Bosons in Three Dimensions and One Dimension”, Phys. Rev. Lett. 107, 230404 (2011).
  • (15) V. Guarrera, D. Muth, R. Labouvie, A. Vogler, G. Barontini, M. Fleischhauer, and H. Ott, ”Spatiotemporal fermionization of strongly interacting one-dimensional bosons”, Phys. Rev. A 86, 021601(R) (2012).
  • (16) M. Girardeau, ”Relationship between Systems of Impenetrable Bosons and Fermions in One Dimension”, J. Math. Phys. (N. Y.) 1, 516 (1960).
  • (17) E. H. Lieb and W. Liniger, ”Exact Analysis of an Interacting Bose Gas. I. The General Solution and the Ground State”, Phys. Rev. 130, 1605 (1963);
    E. H. Lieb, ”Exact Analysis of an Interacting Bose Gas. II. The Excitation Spectrum”, Phys. Rev. 130, 1616 (1963).
  • (18) L. P. Pitaevskii and S. Stringari, ”Bose–Einstein Condensation and Superfluidity”, (Clarendon Press, Oxford, 2016).
  • (19) F. Meinert, M. Panfil, M. J. Mark, K. Lauber, J.-S. Caux, and H.-C. Nägerl, ”Probing the Excitations of a Lieb-Liniger Gas from Weak to Strong Coupling”, Phys. Rev. Lett. 115, 085301 (2015).
  • (20) N. Fabbri, M. Panfil, D. Clément, L. Fallani, M. Inguscio, C. Fort, and J.-S. Caux, ”Dynamical structure factor of one-dimensional Bose gases: Experimental signatures of beyond-Luttinger-liquid physics”, Phys. Rev. A 91, 043617 (2015).
  • (21) P. P. Kulish, S. V. Manakov, and L. D. Faddeev, ”Comparison of the exact quantum and quasiclassical results for a nonlinear Schrö̈dinger equation”, Theor. Math. Phys. 28, 615 (1976).
  • (22) M. Ishikawa, and H. Takayama, ”Solitons in a One-Dimensional Bose System with the Repulsive Delta-Function Interaction”, J. Phys. Soc. Jpn 49, 1242 (1980).
  • (23) K. V. Kheruntsyan, D. M. Gangardt, P. D. Drummond, and G. V. Shlyapnikov, ”Finite-temperature correlations and density profiles of an inhomogeneous interacting one-dimensional Bose gas”, Phys. Rev. A 71, 053615 (2005).
  • (24) A. H. van Amerongen, J. J. P. van Es, P. Wicke, K. V. Kheruntsyan, and N. J. van Druten, ”Yang-Yang Thermodynamics on an Atom Chip”, Phys. Rev. Lett. 100, 090402 (2008).
  • (25) A. Vogler, R. Labouvie, F. Stubenrauch, G. Barontini, V. Guarrera, and H. Ott, ”Thermodynamics of strongly correlated one-dimensional Bose gases”, Phys. Rev. A 88, 031603(R) (2013).
  • (26) X.-W. Guan and M. T. Batchelor, ”Polylogs, thermodynamics and scaling functions of one-dimensional quantum many-body systems”, J. Phys. A 44, 102001 (2011).
  • (27) M.-S. Wang, J.-H. Huang, C.-H. Lee, X.-G. Yin, X.-W. Guan, and M. T. Batchelor, ”Universal local pair correlations of Lieb-Liniger bosons at quantum criticality”, Phys. Rev. A 87, 043634 (2013).
  • (28) X.-W. Guan, ”Critical phenomena in one dimension from a Bethe ansatz perspective”, Int. J. Mod. Phys. B 28, 1430015 (2014).
  • (29) N. D. Mermin and H. Wagner, ”Absence of Ferromagnetism or Antiferromagnetism in One- or Two-Dimensional Isotropic Heisenberg Models”, Phys. Rev. Lett. 17, 1133 (1966).
  • (30) P. C. Hohenberg, ”Existence of Long-Range Order in One and Two Dimensions”, Phys. Rev. 158, 383 (1967).
  • (31) D. S. Petrov, G. V. Shlyapnikov, and J. T. M. Walraven, ”Regimes of Quantum Degeneracy in Trapped 1D Gases”, Phys. Rev. Lett. 85, 3745 (2000).
  • (32) G. Lang, F. Hekking, and A. Minguzzi, ”Ground-state energy and excitation spectrum of the Lieb-Liniger model: accurate analytical results and conjectures about the exact solution”, arXiv: 1609.08865v3 (2017).
  • (33) T. D. Lee and C. N. Yang, ”Many–Body Problem in Quantum Mechanics and Quantum Statistical Mechanics”, Phys. Rev. 105, 1119 (1957).
  • (34) T. D. Lee, K. Huang, and C. N. Yang, ”Eigenvalues and Eigenfunctions of a Bose System of Hard Spheres and Its Low–Temperature Properties”, Phys. Rev. 106, 1135 (1957).
  • (35) T. Kaminaka, M. Wadati, ”Higher order solutions of Lieb–Liniger integral equation”, Phys. Lett. A 375, 2460 (2011).
  • (36) A. Gudyma, ”Non-equilibrium dynamics of a trapped one-dimensional Bose gas”, Ph.D. Thesis, Université Paris-Saclay, (2015).
  • (37) G. E. Astrakharchik, J. Boronat, I. L. Kurbakov, Y. E. Lozovik, and F. Mazzanti, ”Low-dimensional weakly interacting Bose gases: Nonuniversal equations of state”, Phys. Rev. A 81, 013612 (2010).
  • (38) G. E. Astrakharchik, J. Boronat, J. Casulleras, and S. Giorgini, ”Beyond the Tonks-Girardeau Gas: Strongly Correlated Regime in Quasi-One-Dimensional Bose Gases”, Phys. Rev. Lett. 95, 190407 (2005).
  • (39) M. Valiente, and P. Öhberg, ”Few-Body Route to One-Dimensional Quantum Liquids”, Phys. Rev. A 94, 051606(R) (2016).
  • (40) Z. Ristivojevic, ”Excitation Spectrum of the Lieb-Liniger Model”, Phys. Rev. Lett. 113, 015301 (2014).
  • (41) Y.-Z. Jiang, Y.-Y. Chen and X.-W. Guan, ”Understanding many-body physics in one dimension from the Lieb-Liniger model”, Chin. Phys. B 24, 050311 (2015).
  • (42) S. Prolhac, ”Ground state energy of the -Bose and Fermi gas at weak coupling from double extrapolation”, J. Phys. A: Math. Theor. 50, 144001 (2017).
  • (43) M. Abramowitz and I. A. Stegun, ”Handbook of Mathematical Functions: with Formulas, Graphs, and Mathematical Tables”, (Dover Publications, Mineola, New York, 2012).
  • (44) Important differences between and are known to occur in nuclear physics Bohr2008 (), where they are also employed to identify pairing effects of superfluid nature.
  • (45) A. Bohr and B. R. Mottelson, ”Nuclear Structure: Single Particle Motion”, (World Scientific, 2008), Vol. 1.
  • (46) The explicit expression for the chemical potential is
    which, in the thermodynamic limit, gives
    and for large values of the interaction parameter provides Eq. (56).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description