Arxiv:Cond-Mat/0405155V3

Universality of the one dimensional Bose gas with delta interaction Luigi Amico(1) and Vladimir Korepin(2) (1) MATIS-INFM & Dipartimento di Metodologie Fisiche e Chimiche (DMFCI), Universitádi Catania, viale A. Doria 6, I-95125 Catania, Italy and (2) C.N. Yang Institute for Theoretical Physics, State University of New York at Stony Brook, NY 11794-3840 We consider several models of interacting bosons in a one dimensional lattice. Some of them are not integrable like the Bose-Hubbard others are integrable. At low density all of these models can be described by the Bose gas with delta interaction. The lattice corrections corresponding to the different models are contrasted. I. INTRODUCTION duces quantum fluctuations that dominate the zero temperature phase diagram. At commensurate number of bosons per site the Bose-Hubbard model predicts a Mott The Bose gas with delta interaction is one of the insulator-superfluid quantum phase transition controlled most famous models for strongly interacting bosons [1, 2] by the ratio between the bosonic tunneling rate t and (the corresponding classical dynamics is described by the Coulomb repulsion. In optical applications the sys- the non-linear Schrödinger equation). In one dimension tem parameters are tuned by changing the laser intensity this model possesses infinitely many, mutually commut- [20]; the criticality of the superfluid-insulator quantum ing integrals of the motion and is therefore called inte- phase transition is substantially altered by the confining grable. The integrability of the Bose gas can be for- harmonic potential [21]. mulated within the Quantum Inverse Scattering Method (QISM) [3] that exploit a crucial property of the two- In this paper we focus on one dimension. We remark body scattering matrix of the system: it satisfy the Yang- that d = 1 is more than an academic case for the physi- Baxter equation. This reduces the many-body dynamics cal systems captured by the models under consideration. to a two body dynamics. In the Appendix we explain In fact bosonic chains are commonly realized both with the QISM & the algebraic Bethe Ansatz [it uses a trans- optical confinement [20] and, using a different technol- fer matrix]. In the last decades many authors studied ogy, with Josephson junctions [17]. Quantum effects are integrable lattice regularization of the continuous Bose strongest in the one-dimensional case. For example, in gas [4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14]. The procedure it the dilute limit of optical lattices quantum many body is not straightforward since a direct quantization of lattice effects cause the breakdown of the Gross-Pitaevskii the- classical integrable theory whose continuous limit is the ory (valid for 3d systems) and drive the system to the classical continuous Bose gas field theory leads to “non- so called “Tonks-Girardeau” regime [22] where the exact local” Hamiltonians (that is the interaction involves far results of Lieb and Liniger are crucial [23]. Exact results neighbors). This inconvenience was circumvented in two, of the Bose-Hubbard model would be desirable to study somewhat complementary ways: keep the scattering ma- physical systems like the above ones, at higher densi- trix of the continuous theory unaltered and change the ties. However the exact solution of the one dimensional transfer matrix and then the Hamiltonian [4, 5, 6, 7]; Bose-Hubbard model is still a challenging problem. For alternatively, change the scattering matrix keeping the infinite range hopping the thermodynamics is fully estab- Hamiltonian formally unaltered [8, 9]. In both cases, lished [25]. For nearest neighbor hopping, the coordinate however, the Hamiltonians though (quasi) local, involve Bethe ansatz fails [24] due to the multi-occupancy char- higher order processes whose physics is hard to read-out. acteristic of the bosonic statistics [26]. Recent numerical studies suggested that the distribution of level spacings Another paradigmatic lattice model for interacting obeys Wigner-Dyson statistics [27], this considered as in- bosons is the Bose-Hubbard model. It was proposed dication of non-integrability [28]. first by Haldane [15] (see also [16]), and widely used arXiv:cond-mat/0405155v3 [cond-mat.stat-mech] 14 Dec 2004 in many different contexts. In mesososcopic physics it In this paper we discuss the integrable continuous describes granular superconductors, Helium films, and field theory of the Bose gas as universal low density Josephson junctions arrays [17]. In quantum optics, the limit of various integrable (Faddeev-Takhtadjan-Tarasov Bose-Hubbard dynamics was proposed to model bosonic [7], Izergin-Korepin [4], and quantum Ablowitz-Ladik atoms in optical lattices [18] and as effective theory for [8, 9, 12, 13, 14]) and non integrable (Bose-Hubbard) lat- higher spin chains [19]. The Bose-Hubbard Hamiltonian tice models. In particular a number of results obtained describes charged bosons hopping in a d-dimensional lat- for the Bose gas can be used for the behavior of the sys- tice and experiencing a Coulomb repulsion. The resulting tems phrased by Bose-Hubbard models. Finally we dis- competition between kinetic and electrostatic energies in- cuss the differences between the existing integrable lat- 2 q tice regularization of the Bose gas and the non-integrable The density, given by D = −q ρ(u)du, depends on the Bose-Hubbard model. chemical potential since q = q(h) [3]. The Bose gas is R The paper is laid out as follow. In the next section we characterized by different phases depending on the sign summarize a number of properties of the integrable field of the coupling constant c. Let us describe the ground theory of the Bose gas. In section III we apply some of state at zero temperature: these properties to the Bose-Hubbard model. The lattice i) For c = 0 we have free bosons. In the ground state all integrable theories are discussed in section IV. In section the particles have zero momentum (Bose-condensation). V we draw our conclusions. In the Appendix we summa- ii) In the repulsive case c > 0 the Pauli principle is rize the main ideas of the QISM. valid (this is the distinctive feature of the one dimensional bosons). The ground state is a Fermi sphere. iii) In the attractive case c < 0 the ground state is a II. THE BOSE GAS WITH DELTA single large bound state of all the particles. All the par- INTERACTION ticles stay close to one another (in configuration space). It is a droplet. One can say that in this case the phases We shall start with continuous field theory. This will separate. be the continuous limit of all bosonic models considered In the repulsive case at zero temperature the asymp- in this paper. The one dimensional Bose gas is described totics (large distances) of correlation functions is [3, 29] by an integrable non-relativistic model of field theory [1, † −1/θ 3].The Hamiltonian is: Ψ (x)Ψ(0) = Ax (4) h i A cos [2πDx] Ψ†(x)Ψ(x)Ψ†(0)Ψ(0) = D2 + + B , (5) = dx (∂ Ψ†)(∂ Ψ)+ cΨ†Ψ†ΨΨ hΨ†Ψ . h i x2 xθ HBG x x − Z (1) where the critical exponent θ depends on the density in Here: [Ψ(x), Ψ†(y)] = δ(x y) and [Ψ†(x), Ψ†(y)] = 0. a remarkably simple way through: The coupling constant is denoted− by c, and h is the chem- θ =2 2(q)=4πD/v . ical potential. It is equivalent to a many body quantum Z F mechanical problem with Hamiltonian: The thermal correlation functions can be obtained by re- vF placing: x sinh (πT x/vF ) in these expression. Here N 2 → πT ∂ vF is the Fermi velocity [for the model under considera- HBG = +2c δ(zj zk) . (2) − ∂z2 − tion, it coincides with the velocity of sound]. So in the j=1 j N≥j>k≥1 X X repulsive case the asymptotic of the thermal correlation The energy is E = N u2 h ; the momenta u’s obey functions is j=1 j − Bethe equations: v −1/θ P Ψ†(x)Ψ(0) = F exp( πT x/v θ) , (6) h iT πT − F N u u + ic 2 iuj L j k 2πT e = − j =1 ...N. (3) † † −2πT x/vF u u ic Ψ (x)Ψ(x)Ψ (0)Ψ(0) T = B2 e 6 j k h i vF kY=j − − 2πT θ Let us consider repulsive case c > 0. In the thermody- + B e−πT θx/vF cos(2πDx) . (7) 3 v namic limit N,L such that N/L = D = const the F ground state energy-density→ ∞ [per unit length] is The coefficients A,B,B2,B3 are related to certain form q q factors [29]. E = ǫ (u)ρ(u)du = (u2 h)ρ(u)du . 0 − Z−q Z−q III. APPLICATION TO THE BOSE-HUBBARD The function ρ(u) describes the distribution of momenta MODELS in the ground state. It is defined by the equation: The Hamiltonian for the one dimensional Bose- q . 2πρ(u)=1+ K(u v)ρ(v)dv = (u) , Hubbard model reads −q − Z Z Ns † † where (u) is called dressed charge. The kernel of the HBH = [U(ni 1) µ]ni t(a ai+1 +a ai) , (8) − − − i i+1 Z 2 2 − integral equation is K(x)=2c/(c +x ). For further pur- i=XNs poses it is convenient to define a function τ(u) through † where the operators ni := ai ai count the number of q † bosons at the site i; operators ai, ai obey the canoni- 2πτ(u)= K(u q) dvK(u v)τ(v) . † − − − cal commutation relations [ai,a ]= δij and and 2Ns +1 Z−q j 3 is the number of sites.

Load more