Tan's contact of a harmonically trapped one-dimensional Bose gas: strong-coupling expansion and conjectural approach at arbitrary interactions

We study Tan's contact, i.e. the coefficient of the high-momentum tails of the momentum distribution at leading order, for an interacting one-dimensional Bose gas subjected to a harmonic confinement. Using a strong-coupling systematic expansion of the ground-state energy of the homogeneous system stemming from the Bethe-Ansatz solution, together with the local-density approximation, we obtain the strong-coupling expansion for Tan's contact of the harmonically trapped gas. Also, we use a very accurate conjecture for the ground-state energy of the homogeneous system to obtain an approximate expression for Tan's contact for arbitrary interaction strength, thus estimating the accuracy of the strong-coupling expansion. Our results are relevant for ongoing experiments with ultracold atomic gases.


Introduction
Continuous experimental progresses in trapping and cooling atomic gases have led to the realization of one-dimensional (1D) geometries, the observation of the effect of quantum fluctuations [1,2], and the reach of the strongly correlated, Tonks-Girardeau regime [3,4], where many of the properties of bosons are the same as those of free fermions [5]. Several physical properties have been experimentally studied, e.g. the density profiles [3], the two-and three-body correlation functions [6,7], the collective excitation spectrum [8,9], and transport of an impurity in a 1D gas [10] (see e.g. [11] for a comprehensive review).
One of the most common experimental observables is the momentum distribution. This quantity embeds the one-body properties of the gas, in particular it corresponds to the Fourier transform of the one-body density matrix, yielding information on the first-order spatial coherence of the gas. At zero temperature, an interacting 1D Bose gas is predicted to display quasi-off-diagonal long-range order, i.e. an algebraic decay of the one-body density matrix at large distances [12]. High-precision measurements of the momentum distribution are becoming available [13], calling for new theoretical developments.
Recent experimental studies allow to access in particular the high-momentum region of the momentum distribution [14,15]. The high-momentum tails of the momentum distribution n(k) of a gas with contact interactions display a universal n(k) ∝ 1/k 4 decay [16,17], which originates from the zero-range interaction potential: bosons with contact interactions have a cusp in the many-body wavefunction whenever the relative distance of each pair of particles vanishes. This implies a nonanalyticity in the one-body density matrix at short distances, in the form ρ 1 (x) ∼ |x| 3 , which leads to the high-momentum tails.
The weight of the momentum distribution tails, known as Tan's contact , is an important two-body quantity, related to various physical observables, such as the interaction energy, the two-body correlation function at zero distance, the rf-spectroscopy, through Tan's relations [18,19,20,21,22,23]. For a homogeneous 1D Bose gas, Tan's contact at zero temperature can be obtained exactly using the Bethe Ansatz solution for the ground-state energy [24,17,25]. In the experiments, ultracold atomic gases are kept together by an external trap, which is in most cases a harmonic confinement. Tan's contact of an interacting one-dimensional Bose gas under harmonic confinement has been first studied in Ref. [17], where a numerical solution based on the local-density approximation (LDA) of the homogeneous-system result has been provided. Analytically, again in the local-density approximation, a strong-coupling expansion has been derived to zeroth order [17], as well as corrections to first and second order [26]. A comparison with matrix-product state simulations has shown that the local-density approximation works surprisingly well [26], even for a small number of particles.
The study of the equation of state of an interacting 1D Bose gas has received a renewed attention [27,25]. In particular, in Ref. [25] we have derived a strongcoupling expansion for the ground-state energy to an unprecedented accuracy as well as a conjectural expression which is extremely close to the exact numerical solution for a wide range of interaction strengths. Using the above results combined with the local-density approximation, in this work we propose a method to obtain the strongcoupling expansion for Tan's contact of a 1D Bose gas under harmonic confinement to arbitrary order. We also use the conjecture for the ground-state energy to obtain Tan's contact for all interaction regimes, from weak to strong coupling, as well as the lowest-order term in the weak-coupling expansion.
2 Local-density approximation for Tan's contact of a 1D Bose gas We consider N bosonic atoms of mass m at zero temperature, confined by a tight atomic waveguide to a one-dimensional geometry. The atoms interact via a contact potential v(x − x ′ ) = gδ(x − x ′ ), where g is the one-dimensional interaction constant, related to the three-dimensional scattering length of the atoms [28]. In the longitudinal direction the atoms are further confined by an external harmonic confinement V ext (x) = mω 2 0 x 2 /2, describing the optical or magnetic trapping present in the experiments with ultracold atoms. The Hamiltonian of the system reads In the homogeneous system, the interaction strength is measured in terms of the dimensionless parameter γ = mg/(h 2 ρ), where ρ is the average linear density of the homogeneous gas. In the case of harmonic confinement, introducing the harmonicoscillator length a ho = h/(mω 0 ) and the one-dimensional scattering length a 1D = −2h 2 /(gm), the corresponding dimensionless parameter is [29]. According to Tan's sweep relation, Tan's contact in 1D is related to the groundstate energy E of the gas according to [30] In order to determine Tan's contact of the harmonically-confined gas, we employ the density-functional approach developed in [26]. In detail, we define a functional E[ρ] of the density ρ(x) which, in the local-density approximation, reads where ǫ is the ground-state energy density of the homogeneous gas. The ground-state density profile is obtained by minimizing the energy functional, i.e. setting δE/δρ = 0. This yields an implicit equation for the density profile, where R T F = 2µ/(mω 2 0 ) is the Thomas-Fermi radius in a harmonic trap, and the dimensionless average ground-state energy per particle e is such that ǫ(ρ) =h 2 2m ρ 3 e(γ). The chemical potential µ is fixed by imposing the normalization condition N = dx ρ(x). Combining Eqs. (2) and (3), together with the inhomogeneous density profile (4), we obtain Tan's contact within the LDA: This expression readily generalizes the known result for the homogeneous system (see e.g. Ref. [31]): In practice, since e(γ) is not known analytically for the Lieb-Liniger model, one needs to rely on an approximation scheme to obtain Tan's contact in the homogeneous case, and in the trap within the LDA. We proceed by detailing procedures based on various approximations that allow, in principle, to reach excellent accuracy over the whole interaction range.

Scaling relations and methods
In this Section, we combine the various relations obtained above to systematically compute Tan's contact with increasing accuracy. For the energy density ǫ(ρ), we take the solution stemming from the Lieb-Liniger Bethe Ansatz solution, and use various schemes to estimate e(γ), i.e. a strong-coupling expansion where {a k } k≥0 = {π 2 /3, −4π 2 /3, 4π 2 , . . . } is currently known up to order 20 [25], as well as the conjecture where 32 15 , − 96 35 X + 848 315 , . . . } are polynomials explicitly known up to order 6 [25]. Expression (8) is very close to the exact numerical Bethe Ansatz solution for a wide range of interaction strengths. It is also a priori possible to rely on a conjectural expansion in the weakly-interacting regime, where {a ′ k } k≥0 = {1, − 4 3π , 1 6 − 1 π 2 , . . . } is analytically known up to second order at the time being [32]. Higher-order terms were obtained numerically in a very recent work and the exact value of a few of them has been guessed [33].

Strong-coupling expansion for Tan's contact
We first derive the strong-coupling expansion of Tan's contact for a harmonically trapped gas, based on the strong-coupling expansion Eq. (7) for the ground stateenergy of the homogeneous system. Combining Eq. (4), the normalization condition, and the relation g =hω 0 a ho √ N α 0 , the natural rescaled variables are ρ ≡ ρ a ho / √ N , µ ≡ µ/(Nhω 0 ) and x ≡ x/R T F , whereupon we obtain the following set of equations: where A k ≡ (k + 3)a k , and also A straightforward approach to solve the above equations proceeds as follows. One truncates the series in Eq. (10) to order n. Using the values for the coefficients {a k } k≤n from Eq. (7), one can express ρ as the root of a (n + 2) th -degree polynomial as a function of µ and inject into Eq. (11) to find µ, and thus ρ explicitly, consistently to order n by expanding the obtained solution in 1/α 0 . This approach, when carried analytically, suffers from two major drawbacks. First, it is limited to n = 2 since in general, polynomials of order strictly higher than 4 can not be solved by radicals. Second, among the n possible solutions, one has to select the physical one. Using Cardan's method for n = 1 and Ferrari's method for n = 2, putting the roots of the polynomials in a trigonometric form for simplicity, we have used this approach to compute the strong-coupling expansion of Tan's contact to second order in 1/α 0 . This has also been useful to benchmark the more general method described just below. In order to obtain higher-order terms in the strong-coupling expansion, we have developed a more efficient procedure, which allows to systematically compute the expansion of Tan's contact to arbitrary order. The latter relies on the following expansions: µ = +∞ functions, injected into Eqs. (10,11). By analysis, we find that the consistency condition at a given order converts into b j f j (x) = j m=0 b mj (1 − x 2 ) (m+1)/2 , where {b mj } are unkwown coefficients of an upper triangular matrix. Then, synthesis yields: where B is the Euler Beta function. Equations (12) and (13) are the final set of equations. As can be seen, solving the system truncated to order n requires the solution at all lower orders, thus the procedure becomes increasingly lengthy. Moreover, at each step Eq. (12) splits into n+1 independent equations, obtained by equating the coefficients of (1 − x 2 ) (1+m)/2 m=0,...,n in the LHS and RHS. One thus needs to solve a system of n+2 equations to obtain c n and {b mn } m=0,...,n , but fortunately, n of them, giving b mj , m ≥ 1, are fully decoupled.
As a final step, using Eq. (5), we obtain Tan's contact . In natural units imposed by the scaling, i.e. taking C LDA = C LDA a 3 ho /N 5/2 , we have where B k = (k+1)a k+1 . Note that the value is positive because a 1 < 0. At order n, the condition k ′ + k + j ′ = n, where j ′ is the power of α 0 in the integrand, shows that the coefficient of order n is a sum of n+2 n integrals. One of them involves a n+1 , so one needs to know the function e(γ) to order n+1 in 1/γ to obtain the expansion of the Tan's contact to order n in 1/α 0 .

Tan's contact at arbitrary interactions from the conjecture
To check the validity of the strong-coupling method, we have derived another resolution scheme, from Eq. (8). After rescaling and straightforward algebra, one obtains the mixed form Once truncated at order n, it yields a polynomial in 1/(γ + 2) whose roots are found numerically. We proceed as follows. We start from a guessing value of µ and from Eq. (15) we obtain ρ(µ, x). We integrate the density and we use the difference N − ρ dx as control parameter of the accuracy of the initial value for µ. A findroot subroutine exploits this contol parameter to converge to the µ value that ensures the density normalization to N bosons. Once the correct density is computed, we get C LDA from Eq. (5).

Weak-coupling expansion for Tan's contact
At weak iteractions, we derive an expression for Tan's contact using the weak-coupling expansion Eq. (9). Using the same notations as above, we obtain On the other hand, the normalization condition is given by Eq. (11) as before. Here, it is not obvious to what order one should truncate the expressions to obtain a consistent expansion at given order, nor to find the variable in which to expand, as can be seen by evaluating the first orders. Considering only the k = 0 term in the sum, one easily finds ρ(x) = (9/32) 1/3 (1 − x 2 )/α 1/3 0 and µ(α 0 ) = (9/32) 1/3 α 2/3 0 . The expansion to next order is problematic. If one retains terms up to k = 1, corresponding to the Bogoliubov approximation, since the coefficient a ′ 1 is negative, the equation of state at sufficiently large density becomes negative [24]. Then it is not possible to use it to perform the local-density approximation. One may also recall that the local-density approximation breaks down at very weak interactions, where it is not accurate to neglect the quantum tails in the density profile. In this regime, a different scaling parameter becomes relevant [34,35].

Results for Tan's contact
Following the approach presented in Sec. 3.1, the strong-coupling expansion reads: This expression agrees with the zero order one obtained in [17] and with the one derived for a κ-component balanced spinful Fermi gas in [26] to order 2: in fact one may find the bosonic result by taking the limit of very large number of fermionic components κ → ∞ [36]. The weak-coupling expansion obtained from Sec. 3.3 reads: in agreement with [17]. Figure 1 summarizes our results for Tan's contact. Notice that, although the contact is scaled by the overall factor N 5/2 /a 3 ho , it still depends on the number of particles through the factor α 0 /2 = a ho /|a 1D | √ N . First, we notice that the results based on the conjecture are extremely close to the ones obtained from the full solution of the Bethe-Ansatz equation of state in Ref. [17]. Second, by comparing the strong-coupling expansion with the results of the full calculation, we notice that the expansion (17) is valid down to α 0 /2 ≃ 3, and provides an useful analytical expression for Tan's contact in harmonic trap. In order to accurately describe the regime of lower interactions one would need a considerable number of terms in the strong-coupling expansion. The same is true for the series expansion for the ground-state energy of the homogenenous gas [25]. The use of the conjecture (8) is then a valuable alternative with respect to solving numerically the Bethe-Ansatz integral equations, the weak coupling expansion being applicable only for very weak interactions α 0 /2 ≃ 0.05.

Conclusions and outlook
In conclusion, in this work we have determined Tan's contact for a harmonicallytrapped 1D Bose gas. In particular, using an asymptotic expansion of the groundstate energy of the homogeneous system and the local-density approximation, we have developed a general method to obtain a strong-coupling expansion of Tan's contact to arbitrary order, and have provided the coefficients for the expansion up to order four. We have tested it against a calculation of Tan's contact based on a conjecture of the ground-state energy at arbitrary interactions [25] as well as to the full solution of the Bethe-Ansatz equations provided in Ref. [17]. The strong-coupling expansion yields an accurate expression at large interactions, but requires a considerable number of terms to obtain good accuracy at intermediate interactions. In this parameter regime it is then useful to apply the method based on the conjecture for the equation of state. In outlook, the method presented in this work could be used to calculate higher-order terms of the strong coupling expansion with the aim of resumming the terms as was done to obtain the conjectural expression for the homogeneous gas. The local-density approximation could be tested by comparing with ab-initio numerical simulations. It would be interesting to generalize this method to the case of multicomponent 1D gases as well as to finite temperature, beyond the infinitely repulsive Tonks-Girardeau limit of Refs. [37,31].