Thermodynamic behaviour of a onedimensional Bose gas at low temperature
Abstract
We show that the chemical potential of a onedimensional (1D) interacting Bose gas exhibits a nonmonotonic 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 zerotemperature density dependence of the sound velocity. We calculate that coefficient along the crossover between the Bogoliubov weaklyinteracting gas and the TonksGirardeau 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 YangYang model in which the finitetemperature equation of state is obtained numerically by solving Betheansatz equations. A 1D ring geometry is equivalent to imposing periodic boundary conditions and arising finitesize 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:
PACS numbersI 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 nonmonotonic 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 nonmonotonic 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 lowtemperature thermodynamic behaviour of other systems, like onedimensional (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 BoseEinstein condensed, despite the absence of true long range order.
In this work, we investigate the lowtemperature 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 lowT expansion of the chemical potential with the numerical results now available within the YangYang theory Yang1969 (); Yang1970 (), along the whole interaction strength crossover. Previous comparisons were in fact available only in the case of the TonksGirardeau limit Lang2015 (), corresponding to the ideal Fermi gas, where the lowT expansion corresponds to the Sommerfeld expansion. We find that for all intermediate interaction regimes, described at by LiebLiniger (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 LiebLiniger 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
(1) 
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
(2) 
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 TonksGirardeau (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 lowtemperature expansion of the chemical potential, starting from the free energy of an ideal phononic gas. This assumption is fully justified by the lowmomenta behavior of the LiebLiniger excitation spectrum. The lowtemperature expansion exhibits a dependence on temperature, with the coefficient related to the density derivative of the LL sound velocity at zero temperature. The Betheansatz results for the chemical potential are shown to agree very well with the lowtemperature expansion, for the whole BGTG crossover.
In Sec. III we investigate the BG weaklyinteracting gas. By considering the quantum fluctuation contribution in the groundstate 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 lowtemperature 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 groundstate energy of a hardsphere gas.
In Sec. V we derive the lowtemperature 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 finitesize 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 Lowtemperature 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
(3) 
At the sound velocity is related to the density dependence of the chemical potential according to the relation
(4) 
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 longrange properties of onedimensional systems. Figure 1 shows the dependence of the sound velocity on the interaction parameter for the LiebLiniger model, described by Hamiltonian (1). There is a smooth crossover between the meanfield BG value defined as for weak interactions to the TonksGirardeau (ideal Fermi gas) value in the limit of strong repulsion.
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 LiebI particlelike branch corresponds to the Bogoliubov excitation spectrum Lieb1963 (); Kulish1976 (); Pitaevskii2016 (). The LiebII holelike branch is instead associated in the weaklyinteracting regime with the dark soliton dispersion predicted by GrossPitaevskii theory Pitaevskii2016 (); Ishikawa1980 (); Kulish1976 (). The two branches merge into the phononic spectrum for , Fig. 2.
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
(5) 
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 twobody interactions through the dependence of on the interaction parameter . The integral of Eq. (5) yields the following lowtemperature expansion for the free energy
(6) 
which differs from the usual behaviour exhibited by threedimensional () superfluids Pitaevskii2016 () because of the 1D structure of the integral (5). Starting from result (6) for the free energy, one can calculate the lowtemperature expansion of the chemical potential:
(7) 
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
(8) 
and
(9) 
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 lowtemperature Sommerfeld expansion of the chemical potential for the 1D ideal Fermi gas, Eq. (21).
In Figs. 56 we report the temperature dependence of the chemical potential of the system described by Hamiltonian (1) as obtained numerically from the Betheansatz (BA) approach first developed by YangYang Yang1969 (); Yang1970 (); Kheruntsyan2005 (); Lang2015 () for several characteristic values of . The YangYang 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 meanfield to TonksGirardeau 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 meanfield zerotemperature 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):
(10) 
and, despite the absence of superfluidity, it still exhibits the quadratic lowtemperature dependence , which follows from the lowtemperature 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):
(11) 
Notice that, because of the absence of BoseEinstein 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 lowtemperature 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 weaklyinteracting 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 ().
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:
(12) 
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:
(13) 
Figures 5 and 6 point out in a clear way the nonmonotonic 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 lowT 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 , higherorder corrections beyond the linear phononic contribution in the excitation spectrum (3) might be important.
Iii Bogoliubov regime
In the meanfield theory, the chemical potential is linear in density, and the velocity of sound takes the value , see Fig. 1.
The first correction to the meanfield 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 groundstate energy. Therefore in 1D one can consider all ranges of momenta and one finds Lieb1963 ():
(14) 
where
(15) 
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 firstorder correction in the interaction parameter for the ground state energy Lieb1963 ()
(16) 
The same result can be also found by performing a power series expansion of the LiebLiniger 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 higherorder corrections for the other thermodynamic quantities at . For the chemical potential, one finds the result
(17) 
which implies the result
(18) 
for the expansion of the coefficient , where is the meanfield 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 ()
(19) 
which is also reported in Fig. 1, yielding the expression
(20) 
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 TonksGirardeau 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 lowtemperature expansion of the chemical potential in this limit is equal to the first terms of the Sommerfeld expansion (dotdashed line for in Fig. 5) of the 1D ideal Fermi gas, as already pointed out in Lang2015 ():
(21) 
which contains the TG limits of and parameters, Figs. 3 and 4.
Leading corrections to the groundstate 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 ():
(22) 
In the limit of pointlike 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 momentumdependent 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 lowdensity expansion of the equation of state for hard spheres, Eq. (22), and contact potential obtained by solving Bethe equations recursively Astrakharchik2010 (),
(23) 
The nonuniversal correction depends on the shape of the potential and for the LL model it has a noninteger 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 TonksGirardeau Astrakharchik2005 () () and the strongly repulsive Gudyma20152 () () regimes. From Eq. (23), by introducing the parameter (2) and by considering only the leading term, one finds
(24) 
From Eq. (24), one easily calculates the correction of the chemical potential at , Eq. (8):
(25) 
which implies the result Gudyma20152 ()
(26) 
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 ():
(27) 
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 TonksGirardeau 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 Lowtemperature 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
(28) 
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 ()
(29) 
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
(30) 
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),
(31) 
From the above equations one finally finds the low temperature expansion
(32) 
of the adiabatic inverse compressibility, where
(33) 
is its value and we have defined the positive quantity
(34) 
which is reported in Fig. 7 together with its asymptotic limits in the Bogoliubov and Tonks–Girardeau regimes.
v.2 Isothermal inverse compressibility
By fixing the temperature in Eq. (28) and by considering the lowtemperature expansion of the chemical potential (7), one can also calculate the lowtemperature expression for the isothermal inverse compressibility
(35) 
where we have defined the negative dimensionless coefficient
(36) 
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).
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
(37) 
and finite–size effects will become important. This should be contrasted to the threedimensional 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
(38) 
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:
(39) 
In the following, we calculate the finitesize 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 groundstate energy per particle given by
(40) 
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 selfinteraction effects in the leading meanfield term of Eq. (14) which should be replaced by .
By introducing the discretized values of (38), the energy can be rewritten in the form
(41) 
where we have introduced the dimensionless variable
(42) 
depending on the interaction parameter and the function
(43) 
where the adding of the quantity ensures that the term in the sum is counted just once.
By using the EulerMaclaurin expansion (see Appendix A), one can calculate the expression for the series (43) for large values of :
(44) 
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.
For large number of particles, the groundstate energy per particle (41) then takes the form:
(45) 
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 groundstate 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 onedimensional analog of the LeeHuangYang correction in 3D, while the last term accounts for finitesize 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) (dotdashed line), the Betheansatz (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.
The chemical potential can be obtained by deriving Eq. (40) with respect to , at fixed . One finds
(46) 
which can be rewritten as , where is provided by Eq. (42) and we have introduced the series
(47) 
depending on the quantized momenta (38) and such that the zeromomentum term is accounted for once. The EulerMaclaurin expression, applied to the sum (47), yields
(48) 
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 .
Using Eq. (48), one can finally write the following expansion for the chemical potential
(49) 
In Fig. 12 we report the results for the chemical potential as a function of (42) for the thermodynamic limit (17) (dotdashed line), the Betheansatz 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.
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,
(50) 
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
(51) 
after introducing the variable (42) and the quantized momenta (38). As before, we apply the EulerMaclaurin formula and we find the expansion
(52) 
holding in the limit, yielding the asymptotic expansion
(53) 
for the sound velocity.
vi.2 Tonks–Girardeau regime at
According to Girardeau Girardeau1960 (), the groundstate energy of the gas in the stronglyinteracting 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 singleparticle levels in the box,
(54) 
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 finitesize correction close to the TonksGirardeau regime we replace with in Eq. (54) resulting in the following expression for the energy per particle
(55) 
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 hardsphere (HS) like model (55) (dashed and dotted lines) and the Betheansatz solution (symbols) for several values of . We observe a very good agreement between the BA solution and the analytical hardsphere (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 finitesize effects are less important in the TG regime since in the weakly interacting Bogoliubov regime, the correction is amplified by the smallness of .
For strong repulsion we obtain the finitesize correction to the chemical potential
(56) 
and to the sound velocity (4):
(57) 
It is interesting to note that while the finitesize 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 .
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,
(58) 
where is the density operator in momentum representation. The static structure factor gives information about twobody 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 finitesize correction to the sound velocity. The latter effect is rather weak, especially in the TonksGirardeau regime, but is important in the context of the finitesize 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 socalled elastic form factor, does not contribute. Indeed, one finds , which exactly vanishes for .
Figure 15 reports the static structure factor in the TonksGirardeau 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 finitesize corrections to the sound velocity are negligible in the TonksGirardeau 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 smallmomentum 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 lowmomentum dependence. In the thermodynamic limit and the phononic linear behavior is restored.
In Fig. 16 we show the static structure factor for , calculated using the diffusion Monte Carlo method. Similarly, the finitesize 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 finitesize dependence.
While for the TG regime, the linear dependence extends up to , for weaker interactions the linear regime shrinks (compare Figs. 1516). 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. 56 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) — TonksGirardeau (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 nonmonotonic behavior as a function of . The coefficient of the law has been calculated using the LiebLiniger results for the sound velocity and the resulting behavior has been successfully compared with thermodynamic functions obtained from the YangYang 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 Betheansatz. We have found that finitesize 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 higherorder corrections (beyond the law caused by the real excitations of the phononic branch) to the lowtemperature thermodynamic behavior. In particular, it is important to understand the temperature corrections arising due to nonsymmetric spreading of the phononic branch (different beyondlinear behavior of the lower and upper branches) as well as effects originating from the nonlinear 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 EulerMaclaurin 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 EulerMaclaurin expansion which allows to approximate a series as follows Abramowitz2012 ():
(59) 
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
(60) 
and are the –derivatives of the function .
Acknowledgements.
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. FIS201456257C21P. 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. (www.gausscentre.eu) for funding this project by providing computing time on the GCS Supercomputer SuperMUC at Leibniz Supercomputing Centre (LRZ, www.lrz.de). 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.References
 (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 BoseEinstein condensate”, Phys. Rev. A 70, 013608 (2004).
 (5) C. N. Yang and C. P. Yang, ”Thermodynamics of a OneDimensional System of Bosons with Repulsive DeltaâFunction Interaction”, J. Math. Phys. 10, 1115 (1969).
 (6) C. P. Yang, ”OneDimensional 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 onedimensional 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 twodimensional 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, ”TonksGirardeau gas of ultracold atoms in an optical lattice”, Nature 429, 277 (2004).

(10)
T. Kinoshita, T. Wenger, D. S. Weiss,
”Observation of a OneDimensional TonksGirardeau 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, ”ConfinementInduced Resonances in LowDimensional 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, ”ThreeBody 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 onedimensional 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 LiebLiniger 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 onedimensional Bose gases: Experimental signatures of beyondLuttingerliquid 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 OneDimensional Bose System with the Repulsive DeltaFunction Interaction”, J. Phys. Soc. Jpn 49, 1242 (1980).
 (23) K. V. Kheruntsyan, D. M. Gangardt, P. D. Drummond, and G. V. Shlyapnikov, ”Finitetemperature correlations and density profiles of an inhomogeneous interacting onedimensional 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, ”YangYang 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 onedimensional Bose gases”, Phys. Rev. A 88, 031603(R) (2013).
 (26) X.W. Guan and M. T. Batchelor, ”Polylogs, thermodynamics and scaling functions of onedimensional quantum manybody 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 LiebLiniger 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 TwoDimensional Isotropic Heisenberg Models”, Phys. Rev. Lett. 17, 1133 (1966).
 (30) P. C. Hohenberg, ”Existence of LongRange 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, ”Groundstate energy and excitation spectrum of the LiebLiniger 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, ”Nonequilibrium dynamics of a trapped onedimensional Bose gas”, Ph.D. Thesis, Université ParisSaclay, (2015).
 (37) G. E. Astrakharchik, J. Boronat, I. L. Kurbakov, Y. E. Lozovik, and F. Mazzanti, ”Lowdimensional 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 TonksGirardeau Gas: Strongly Correlated Regime in QuasiOneDimensional Bose Gases”, Phys. Rev. Lett. 95, 190407 (2005).
 (39) M. Valiente, and P. Öhberg, ”FewBody Route to OneDimensional Quantum Liquids”, Phys. Rev. A 94, 051606(R) (2016).
 (40) Z. Ristivojevic, ”Excitation Spectrum of the LiebLiniger Model”, Phys. Rev. Lett. 113, 015301 (2014).
 (41) Y.Z. Jiang, Y.Y. Chen and X.W. Guan, ”Understanding manybody physics in one dimension from the LiebLiniger 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
(63) (64)