OSA's Digital Library

Optics Express

Optics Express

  • Editor: Michael Duncan
  • Vol. 12, Iss. 1 — Jan. 12, 2004
  • pp: 42–54
« Show journal navigation

Variational ansatz for the superfluid Mott-insulator transition in optical lattices

J. J. García-Ripoll, J. I. Cirac, P. Zoller, C. Kollath, U. Schollwöck, and J. von Delft  »View Author Affiliations


Optics Express, Vol. 12, Issue 1, pp. 42-54 (2004)
http://dx.doi.org/10.1364/OPEX.12.000042


View Full Text Article

Acrobat PDF (187 KB)





Browse Journals / Lookup Meetings

Browse by Journal and Year


   


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools

Share
Citations

Abstract

We develop a variational wave function for the ground state of a one-dimensional bosonic lattice gas. The variational theory is initally developed for the quantum rotor model and later on extended to the Bose-Hubbard model. This theory is compared with quasi-exact numerical results obtained by Density Matrix Renormalization Group (DMRG) studies and with results from other analytical approximations. Our approach accurately gives local properties for strong and weak interactions, and it also describes the crossover from the superfluid phase to the Mott-insulator phase.

© 2004 Optical Society of America

1. Introduction

During the last years a spectacular development in the storage and manipulation of cold atoms in optical lattices [1

1. M. Greiner, O. Mandel, T. Esslinger, T. W. Hänsch, and I. Bloch, “Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms,” Nature 415, 39 (2002). [CrossRef] [PubMed]

, 2

2. M. Greiner, O. Mandel, T. W. Hänsch, and Immanuel Bloch, “Collapse and revival of the matter wave field of a Bose-Einstein condensate,” Nature 419, 51 (2002). [CrossRef] [PubMed]

] has taken place. Greiner et al. [1

1. M. Greiner, O. Mandel, T. Esslinger, T. W. Hänsch, and I. Bloch, “Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms,” Nature 415, 39 (2002). [CrossRef] [PubMed]

], to name one important example, succeeded in experimentally driving a quantum phase transition between a superfluid and a Mott-insulating phase in bosonic systems. This experimental progress has revived the interest in the Bose-Hubbard model [Eq. (1)] as a generic Hamiltonian for strongly correlated bosons, by which the quantum phase transition can be described [3

3. D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, “Cold Bosonic Atoms in Optical Lattices,” Phys. Rev. Lett. 81, 3108–311 (1998); [CrossRef]

, 4

4. D. Jaksch, V. Venturi, J. I. Cirac, C. J. Williams, and P. Zoller, “Creation of a Molecular Condensate by Dynamically Melting a Mott Insulator,” Phys. Rev. Lett. 89, 040402 (2002) [CrossRef] [PubMed]

]. The Bose-Hubbard Hamiltonian has been used previously in condensed matter physics to study the adsorption of noble gases in nanotubes [5

5. B. Y. Chen, S. D. Mahanti, and M. Yussouff, “Helium atoms in zeolite cages: Novel Mott-Hubbard and Bose-Hubbard systems,” Phys. Rev. Lett. 75, 473–476 (1995). [CrossRef] [PubMed]

], or Cooper pairs in superconducting films with strong charging effects [6

6. Min-Chul Cha, M. P. A. Fischer, S. M. Girvin, M. Wallin, and A. P. Young, “Universal conductivity of two-dimensional films at the superconductor-insulator transition,” Phys. Rev. B 44, 6883–6902 (1991). [CrossRef]

, 7

7. M. P. A. Fischer, G. Grinstein, and S. M. Girvin “Presence of quantum diffusion in two dimensions: Universal resistance at the superconductor-insulator transition,” Phys. Rev. Lett.64, 587–590 (1990);M. P. A. Fischer, “Quantum phase transitions in disordered two-dimensional superconductors,” Phys. Rev. Lett.65, 923–926 (1990). [CrossRef]

]. In this context a lot of work has already been done to characterize the quantum phase transition, the statistics, and the low-energy excitations of the Bose-Hubbbard model [8

8. S. Sachdev, Quantum Phase Transitions, (Cambridge Univ. Press, Cambridge, 2001).

, 9

9. M. P. A. Fischer, P. B. Weichman, G. Grinstein, and D. S. Fisher, “Boson localization and the superfluid-insulator transition,” Phys. Rev. B 40, 546–570 (1989). [CrossRef]

]. However, new interesting questions arise now due to the good tunability of the experiments with optical lattices. In particular, it becomes possible to study time-dependent processes such as driven quantum phase transitions [1

1. M. Greiner, O. Mandel, T. Esslinger, T. W. Hänsch, and I. Bloch, “Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms,” Nature 415, 39 (2002). [CrossRef] [PubMed]

]. A theoretical understanding of such phenomena is challenging, since the characteristics of the superfluid phase—where the atoms tend to delocalize throughout the lattice and large fluctuations in the local density exist —, and the Mott-insulating phase — where the number fluctuations decrease, and a gap in the excitation spectrum opens —, must be covered at the same time. Both regions are separated by a non-analyticity of the spectrum, which implies that a perturbative study [10

10. J. K. Freericks and H. Monien, “Phase diagram of the Bose Hubbard model,” Europhys. Lett. 26, 545–550 (1994). [CrossRef]

, 11

11. J. K. Freericks and H. Monien, “Strong-coupling expansions for the pure and disordered Bose-Hubbard model,” Phys. Rev. B 53, 2691–2700 (1996). [CrossRef]

] works best in strong coupling limit, while a Hartree-Fock-Bogoliubov mean field works best in the superfluid regime. In addition it is possible to develop a mean field theory [12

12. D. van Oosten, P. van der Straten, and H. T. C. Stoof, “Quantum phases in an optical lattice,” Phys. Rev. A 63, 053601 (2001). [CrossRef]

, 13

13. A. M. Rey, K. Burnett, R. Roth, M. Edwards, C. J. Williams, and C. W. Clarck, “Bogoliubov Approach to Super-fluidity of Atoms in an Optical Lattices,” J. Phys. B , 36, 825–841 (2003). [CrossRef]

] based on a Gutzwiller ansatz [3

3. D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, “Cold Bosonic Atoms in Optical Lattices,” Phys. Rev. Lett. 81, 3108–311 (1998); [CrossRef]

, 4

4. D. Jaksch, V. Venturi, J. I. Cirac, C. J. Williams, and P. Zoller, “Creation of a Molecular Condensate by Dynamically Melting a Mott Insulator,” Phys. Rev. Lett. 89, 040402 (2002) [CrossRef] [PubMed]

]: this reproduces the mean field theory in the superfluid limit, as well as the limit of infinite interaction, which raises the hope that the theory also interpolates properly between these limits.

The outline of the paper is as follows: In Section 2 we introduce the quantum rotor model as a possible limit of the Bose-Hubbard Hamiltonian. Next, information about the ground state of the quantum rotor model is obtained variationally as the solution of a Mathieu equation. We can estimate energies, correlation functions and length, and the variance of the density as a function of the only free parameter. A comparison with perturbative estimates demonstrates the accuracy of the method when computing local properties. Since the quantum rotor model is only an approximate description of the optical lattice, in Section 3 we develop a similar variational theory for the Bose-Hubbard Hamiltonian. After bringing the Bose-Hubbard Hamiltonian to an appropriate form, we can estimate the local properties of its ground state. The variational solutions are compared in Sec. 3.3 with the results of DMRG studies of the Bose-Hubbard model. We confirm that the variational method describes very well the local properties of both the Mott-insulator and the superfluid regime, and provides a fairly good interpolation across the phase transition. Finally, in Section 4 we summarize our results.

2. Quantum phase model

2.1. Relation to the Bose-Hubbard model

In this section we show the equivalence of the Bose-Hubbard model

HBH=j=1M[J(aj+1aj+ajaj+1)+U2ajajajaj]U2Mn¯(n¯1),
(1)

and the quantum rotor model for large and integer occupation [6

6. Min-Chul Cha, M. P. A. Fischer, S. M. Girvin, M. Wallin, and A. P. Young, “Universal conductivity of two-dimensional films at the superconductor-insulator transition,” Phys. Rev. B 44, 6883–6902 (1991). [CrossRef]

]. In our notation, M is the number of lattice sites and N=n̄M the number of atoms. Both the Bose-Hubbard model and the quantum rotor model show a phase transition due to the interplay between the kinetic term proportional to J and the interaction term proportional to U. For convenience we have subtracted the ground state energy in the perfect insulator limit U/J→∞.

If we expand a configuration of the lattice using Fock states

ψ=ncnn=ncnn1nM,
(2)

and the number of particles per lattice site is large, nk >1, we may approximate the hopping terms as follows

alajψ=n¯(n¯+1)𝓟Al+Ajn+Δlj
(3)

Here, A±|n〉=|n±1〉 are ladder operators and 𝓟 projects on states with non-negative occupation numbers, nk ≥0. To lowest order the error |Δlj 〉 is

Δlj=ncn(n¯+1)(njn¯)+n(nln)2n¯(n¯+1)n,
(4)

and its norm is bound by

Δljn¯2+(n¯+1)22n¯(n¯+1)σl,
(5)

where σl2 =〈(nl -)2〉 is the variance in the number of particles per lattice site. For the approximation (3) to be valid, the uncertainty in the number of atoms must be small compared to the mean value, ≫σ, and the interaction energy must exceed the neglected terms, Ūn(-1)≫.

HQR=𝓟j[2ρJ(Aj+1+Aj+Aj+Aj+1)+U2(Ajz)2]
(6)

Here ρ=n¯(n¯+1) is approximately the density, Ajz =a†jaj-n̄ is essentially the number operator and we have used that ∑jAjz |ψ〉=0 when we work with states that have a fixed, commensurate number of particles. In the following we will define the energy per lattice site as

ε1MHQR.
(7)

Since the physically interesting states will be concentrated around large occupations, nk =, the usual step now is to drop the projector, 𝓟, and move to the basis of phase states, defined by

nϕ=ein·ϕ(2π)M2,ϕ[π,π]M.
(8)

In doing so, we obtain the identification Aj±e±iϕk and Ajz →-i /∂ϕ j , which produces the usual representation of the quantum rotor model

HQR=j[2ρJcos(ϕjϕj+1)U22ϕj2],
(9)

with the associated state writen as

ψ=(2π)M2ein¯ΣϕkΨ(ϕ)|ϕdMϕ.
(10)

A similar derivation is possible using path integrals [16

16. M. P. A. Fisher and G. Grinstein, “Quantum Critical Phenomena in Charged Superconductors,” Phys. Rev. Lett. 60, 208–211 (1988). [CrossRef] [PubMed]

].

2.2. Variational ansatz

Ψ(ϕ)=j=1Mh(ϕjϕj+1).
(11)
Fig. 1. Instead of working directly with the population of each well, nk , we can use other quantum numbers, wk , defined by the relation nk =wk -wk -1+, and which behave like a set of chemical potentials acting on the barriers that connect neigboring sites.

This representation becomes exact in the Mott-insulating limit, U/J→∞, where hmott (ξ)=1, and in the superfluid limit, U/J→0, where hsf (ξ)=∑n∈ℤδ(ξ-2πn), as can be verified by direct substitution in Eq. (10).

Even though the phase representation is the best one to find a trial wavefunction, it is not the optimal one for performing computations. It is instead more convenient to work with the variables which are conjugate to the phase differences ξj=ϕ j+1-ϕj . These are the new quantum numbers, wk , given by

nk=wkwk1+n¯.
(12)

In terms of these numbers, the ansatz (11) reveals itself as a simple product wavefunction

ψ=h˜(M1)=wh˜w1h˜wM1w1wM1,
(13)

with coefficients given by the Fourier transform

h˜m=h(ξ)eimξdξ.
(14)

As sketched in Fig. 1, the wk play the roles of chemical potentials which are established between different wells: the difference between the potentials on the extremes of a site gives the fluctuations over the mean and commesurate occupation . In this picture

HQR=k=1M1[2ρJ(k++k)+U2(kzk1z)2],
(15)

where ∑±|w〉=|w±1〉 are new ladder operators and ∑ z |w〉=w|w〉.

By minimizing the energy associated with HQR over all states within a given ansatz we can both obtain an upper bound to the energy of the ground state and approximate its wave function. A simple computation with our product ansatz leads to the result [compare (7)]

ε[h˜]4ρJReΣ++U(Σz)2UΣz2,
(16)

where the expected values are computed over a single connection, 〈∑ z 〉=∑ww |hw |2, and the wavefunctions are assumed to be normalized, ∑ w |hw |2=1. Since the stationary states have a well defined parity, (-w)=(-1) Pw , the optimal variational state must satisfy the linear equation

2ρJ(h˜j+1+h˜j1)+Uj2h˜j=εesth˜j,
(17)

which is nothing but the Fourier transform of a Mathieu equation

[U22ξ22ρJcos(ξ)]h(ξ)=εesth(ξ).
(18)
Fig. 2. Estimates for (a) energy energy per lattice site and (b) density fluctuations of the quantum rotor Hamiltonian (6) obtained with the variational method (solid), and perturbative calculations for U≪J (dashed) and UJ (dots).

The estimated ground state energy per site is given by the lowest eigenvalue of either equation.

Using the product ansatz and the same approximations required to derive HQR , we can also compute other properties of the ground state. For instance, the variance of the number of atoms per lattice site

σj2=(ajajn¯)2=2(Σz)2,
(19)

and the correlation functions

aj+1aj=ρjργ1,
(20)
aj+laj=ρk=jj+l1k=ρ(γ1)l,
(21)

which decay exponentially with the distance. This implies that the ansatz (13) only describes properly the decay of the correlations in the Mott-insulating regime, since the correlations in the superfluid regime follow a power law decay. However, as we will see below, local observables (σ, γ1, ε …) are properly estimated even if long–range ones are not.

We have solved Eq. (17) numerically in a truncated space. The results are summarized in Fig. 2, where we also plot reference estimates arising from two other analytical methods. In the limit UJ we compare with a first order perturbative calculation around the solution Ψ(ϕ⃗)=1, which is possible thanks to an energy gap of order 𝒪(U) in the excitation spectrum. In the limit JU we rather use a spin wave or harmonic approximation in which the cosineterms of the Hamiltonian HQR are expanded up to second order in the phase difference between neighboring sites [see Sec. 2

2. M. Greiner, O. Mandel, T. W. Hänsch, and Immanuel Bloch, “Collapse and revival of the matter wave field of a Bose-Einstein condensate,” Nature 419, 51 (2002). [CrossRef] [PubMed]

.3

3. D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, “Cold Bosonic Atoms in Optical Lattices,” Phys. Rev. Lett. 81, 3108–311 (1998); [CrossRef]

]. This approximation is valid in the superfluid regime, where the phase does not vary much between neighboring wells. From the graphical comparison we see that the variational wavefunction provides a fairly accurate description of the ground state of the quantum rotor model in both the superfluid and insulating limits. As a side note, we must remark that this ground state has a divergent fluctuation of the number of particles per site as J→0. This differs from the expected behavior of the ground state of the original Bose-Hubard Hamiltonian, and it reminds us that HQR can only model the atomic ensemble when the variance, σ, is small compared to the mean occupation number, .

2.3. Harmonic approximations to the quantum phase model

In the limit JU it is possible to estimate the ground state of the rotor model (6) analytically. Since we are deep in the superfluid regime, the wavefunction will be concentrated around the line ϕ1=ϕ2=…=ϕM, and we can approximate

HQRj[U22ϕj2ρJ(ϕjϕj+1)2].
(22)

If we remove the periodic boundary conditions on ϕi and change variables, the preceding Hamiltonian may be diagonalized, H=kωk(bkbk+12), with frequencies given by ωk =(8ρJU)1/2 |sin(πk/M)|, where k is an integer in the range -M+1<2k<M and labels the different values of the momentum in the lattice. The ground state energy [Fig. 4(a)] may be estimated as the zero-point energy of the harmonic oscillator. For large M, the sum over k may be replaced with an integral, giving

Eg2Mπ2ρJU.
(23)

σ1π8JρU.
(24)

3. The Bose-Hubbard model

3.1. Coherent states

The phase coherent states |ϕ〉 are defined by

n|ϕ=einϕn!.
(25)

Unlike the phase states defined in Sec. 2, they are not orthogonal to each other, 〈ϕ|θ〉=exp[ei(θ-ϕ) ], but they form a complete basis and an expansion like (10) is still possible. A nice property of the coherent states is that we can rewrite the operators a, a , a a, etc, in terms of the phases very easily. For instance, aϕ=ieiϕϕϕ, and aaϕ=iϕϕ.

Using this representation, we obtain an effective Hamiltonian for the wavefunction Ψ(ϕ⃗), i.e. HBHψ=(2π)M2dMϕein¯Σϕk[HcohtΨ(ϕ)]ϕ, which has the form

Hcoht=Ji,j[2(n¯+1)cos(ϕiϕj)iei(ϕiϕj)ϕj]+U2j(2ϕj2).
(26)

Here Hcoht stands for the tranpose of Hcoh . This operator was already used in Ref. [29

29. J. R. Anglin, P. Drummond, and A. Smerzi, “Exact quantum phase model for mesoscopic Josephson junctions,” Phys. Rev. A 64, 063605 (2001). [CrossRef]

] to study the Bose-Hubbard model with only two sites. On the one hand, it is a non-Hermitian operator1 and we cannot do a simple variational study. On the other hand the Hamiltonian still depends on the phase differences, and it is reasonable to look for approximate eigenstates which have the form (11). This will be done in the following section.

3.2. Variational procedure for non-Hermitian operators

In this section we will find the best variational function which has the product form of Eq. (11). However, as it happened in Section 2, instead of working with phase variables it will be more convenient to develop a representation of the Bose-Hubbard Hamiltonian in terms of connections. This is once more a two-steps process. First we use a similarity transformation suggested by the definition of the coherent states

On=k=1Mnk!n.
(27)

Since OajO -1=A-j and OajO-1=AjzA+j , we find

Hcoh=OHBHO1
=Ji,j(Aiz+n¯)Ai+Aj+U2j(Ajz)2.
(28)

Hcoh=H1+H2,
H1=j[Jn¯Σx+iJjzjy+U(jz)2],
H2=j[J(j1zj+j+1zj)+Ujzj+1z],
(29)

into terms which are local, H 1, and terms which involve pairs of connections, H 2.

For the quantum rotor model we proved that the optimal product wavefunction was an eigenstate of a Hamiltonian which did not couple connections, like H 1. The difference now is that, since the operator Hcoh is not Hermitian, we cannot establish a variational principle and that argument is no longer valid. Nevertheless, we will again propose a variational ansatz which is an eigenstate of the local operator H 1|〉⊗ M =est |〉⊗ N . Using the following equality

ε0=minψ0ψHBHψψ2=minχ0χO2HcohχχO2χ,
(30)

and the product ansatz |χ〉=|〉⊗ M , we arrive to an upper bound for the lowest eigenvalue of the Bose-Hubbard Hamiltonian, expressed in terms of the non-Hermitian one

ε0εest+1Nh˜MO2H2h˜Mh˜MO2h˜Mεest+Δεest.
(31)

1The hermiticity of HBH is maintained due to an implicit projection that takes place when we reconstruct the state HBH|ψ〉 from HtcohΨ(ϕ⃗) (See Ref. [29]).

The way to use this variational principle is as follows. First, for a given J and U we compute the lowest eigenstate of H 1 and this way obtain . After the equivalence (14), finding the ground state of the local Hamiltonian H 1 becomes equivalent to solving a modified Mathieu equation

[U2ξ22J(n¯+1)cos(ξ)2Jsin(ξ)ξ]h=εesth,
(32)

which describes exactly the static properties of a pair of sites with open boundary conditions [29

29. J. R. Anglin, P. Drummond, and A. Smerzi, “Exact quantum phase model for mesoscopic Josephson junctions,” Phys. Rev. A 64, 063605 (2001). [CrossRef]

]. Once we have εest we must still compute the correction Δεest using a rather straightforward expansion which is shown in Sec. 3.4. Surprisingly, Δεest happens to be negative, so that it is actually an improvement over the simple estimate given by εest [See Fig. 5].

From the optimal variational state, |ψ〉=O-1|〉⊗ M , and the estimate for the energy, εvar =εest +Δεest , we may compute other observables. For the density fluctuations and nearest neighbor correlations we use the virial theorem

aj+1aj=Jεvar,
(33)
σ2=Uεvarn¯2,
(34)

whereas for other properties one has to evaluate numerically the matrix products shown in Sec. 3.4. This allows us to prove that for the product states |〉⊗M correlations decay exponentially, opposite to what is expected in the superfluid phase, whose correlations should decay algebraically. Nevertheless, as we will see next, this family of states does estimate accurately the local properties of the optical lattice.

3.3. Comparison to DMRG results

We will now compare the results for the ground state energy, the correlation functions, and the variance of the particle number provided by the two variational ansatz, (18) and (32), and the Gutzwiller ansatz [15

15. A. van Otterlo, K.-H. Wagenblast, R. Baltin, C. Bruder, R. Fazio, and G. Schön, “Quantum phase transitions of interacting bosons and the supersolid phase,” Phys. Rev. B 52, 16176–16186 (1995). [CrossRef]

] with those obtained by DMRG studies of the Bose-Hubbard model. The DMRG, developed 1992 by White [17

17. S. R. White, “Density matrix formulation for quantum renormalization groups,” Phys. Rev. Lett. 69, 2863–2866 (1992). [CrossRef] [PubMed]

, 18

18. S. R. White, “Density-matrix algorithms for quantum renormalization groups,” Phys. Rev. B. 48, 10345–10356 (1993). [CrossRef]

] in the area of condensed matter theory, is a very powerful numerical tool to investigate static and dynamic properties of strongly correlated quasi-one-dimensional spin, fermionic or bosonic quantum systems. The DMRG is an essentially quasi-exact numerical method. The fundamental ideas stem from real space renormalization methods: the system size is grown iteratively while the (exponentially diverging) size of the Hilbert space is kept constant by decimation. Hereby one tries to retain only that subset of states that is essential to describe the physical quantity under consideration. In DMRG these are expectation values with respect to low-lying states (“target states”), and in particular with respect to the ground state wave function.

DMRG builds up the system linearly: at each growth step, suitable density matrices for the target states are derived that yield information on the relevance of Hilbert space states. Building on this information, the states and operators are projected onto Hilbert subspaces of fixed dimension M containing the most relevant states. M is chosen to be small enough to be handled numerically, but large enough to obtain the desired accuracy; numerical results can be extrapolated in M to the exact limit of infinite M in the thermodynamic limit. However, results presented here have converged for the largest M considered and no further extrapolation was necessary.

Fig. 3. (a) The ground state energy per site, ε, (b) nearest neighbor correlation, c 1=〈aj+1aj 〉, and (c) variance of the number of atoms per site, σ 2=〈(nj -)2〉. Plots (b) and (c) use a log-log scale. The results of the DMRG (solid line) are obtained on a system with 128 sites, a maximum occupation number of 9 bosons per site and a reduced space of states of about 200 states. The estimates from the variational theory are plotted using dashed lines. The vertical lines mark the location of the phase transition according to [11]. The mean occupation numbers are denoted with circles (=1), diamonds (=2) and boxes (=3).

Details on the DMRG method may for example be found in [19

19. I. Peschel, X. Wang, M. Kaulke, and K. Hallberg, Density-matrix Renormalization, (Springer-Verlag, Berlin, 1998).

]. In the case of the Bose-Hubbard model the DMRG has been used to study properties of the system [20

20. T. D. Kühner, S. R. White, and H. Monien, “One-dimensional Bose-Hubbard model with nearest-neighbor interaction,” Phys. Rev. B 61, 12474–12489 (2000). [CrossRef]

, 21

21. T. D. Kühner, Diploma work (1997), University of Bonn.

, 22

22. S. Rapsch, U. Schollwöck, and W. Zwerger, “Density matrix renormalization group for disordered bosons in one dimension,” Europhys. Lett. 46, 559–564 (1999). [CrossRef]

]. The results of DMRG agree very well with exact diagonalization results for small systems, with quantum Monte-Carlo simulations e.g. [23

23. G. G. Batrouni, R. T. Scalettar, and G. T. Zimanyi, “Supersolids in the Bose-Hubbard Hamiltonian,” Phys. Rev. Lett. 74, 2527–2530 (1995). [CrossRef] [PubMed]

, 24

24. G. G. Batrouni and R. T. Scalettar, “World-line quantum Monte Carlo algorithm for a one-dimensional Bose model,” Phys. Rev. B 46, 9051 (1992). [CrossRef]

, 25

25. P. Niyaz, R. T. Scalettar, C. Y. Fong, and G. G. Batrouni, “Phase transitions in an interacting boson model with near-neighbor repulsion,” Phys. Rev. B 50, 362–373 (1994). [CrossRef]

, 26

26. N. V. Prokof’ev, B. V. Svistunov, and I. S. Tupitsyn, Phys. Lett. A 238, 253 (1998). [CrossRef]

], and with 13th order perturbation theory [27

27. N. Elstner and H. Monien, “Dynamics and thermodynamics of the Bose-Hubbard model,” Phys. Rev. B 59, 12184–12187 (1999) and ref. therein. [CrossRef]

].

We have used the DMRG to study the properties of the ground state of the Bose-Hubbard model on one-dimensional lattices with 128 sites, and commensurate fillings =1,2, and 3. In Figs. 3(a–c) we show the mean energy per site ε, the nearest neighbor correlations, c 1=〈aj+1aj 〉, and the variance σ of the density, calculated both with the DMRG and with the variational estimates developed above. As expected, there are no indications of the phase transitions in these quantities, neither in the variational results nor in the numerical solutions. Rather, an inflexion of the nearest neighbor correlation points out the location of the superfluid-insulator transition which lies roughly between 3 and 4 (see [20

20. T. D. Kühner, S. R. White, and H. Monien, “One-dimensional Bose-Hubbard model with nearest-neighbor interaction,” Phys. Rev. B 61, 12474–12489 (2000). [CrossRef]

] and ref. therein). The agreement of the two methods is fairly good above the phase transition and below it.

Fig. 4. (a) The ground state energy per site, ε, (b) nearest neighbor correlation, c 1=〈aj+1aj 〉, and (c) variance of the number of atoms per site, σ 2=〈(nj -)2〉. Plot (b) and (c) are in log-log scale. Using filling factor n̄=1, we show results from the variational model for the Bose-Hubbard model using phase coherent states (solid), the quantum rotor model (dashed), the Gutzwiller ansatz for the Bose-Hubbard Hamiltonian (dots) and DMRG (circles). Vertical dash-dot lines mark the location of the phase transition according to [11].

A more detailed comparison is provided in Fig. 4 for the case =1. In this figure we plot together results from the DMRG, the variational ansatz derived above, the quantum rotor model and the well-known Gutzwiller ansatz. The Gutzwiller ansatz [28

28. D. S. Rokhsar and B. G. Kotliar, “Gutzwiller projection for bosons,” Phys. Rev. B 44, 10328–10332 (1991). [CrossRef]

, 3

3. D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, “Cold Bosonic Atoms in Optical Lattices,” Phys. Rev. Lett. 81, 3108–311 (1998); [CrossRef]

, 4

4. D. Jaksch, V. Venturi, J. I. Cirac, C. J. Williams, and P. Zoller, “Creation of a Molecular Condensate by Dynamically Melting a Mott Insulator,” Phys. Rev. Lett. 89, 040402 (2002) [CrossRef] [PubMed]

] is a variational ansatz made of a product of single-site wave functions, |Ψ G 〉=ΠMj=1j〉, where |Φj〈=∑∞m=0 f (j)|mj 〉 and fm(j) are constants. Such a wavefunction cannot be used in the one-dimensional Mott insulator regime, because a perturbative study shows that the corrections of order 𝒪(J/U) are lost and all correlations become zero. However this ansatz gives good results in the superfluid regime, where the long-range order is well described by |Ψ G 〉, and we can use these results and those of the DMRG to assert the accuracy of our variational estimates. As Fig. 4 shows, in the Mott insulator regime, the DMRG results agree perfectly with our variational theory for the Bose-Hubbard model and for the quantum rotor model. Close to the phase transition is the point at which the quantum rotor model no longer describes well the atoms in the optical lattice due to the growth of density fluctuations. At this point we also observe a small disagreement between the DMRG and the coherent states, which is due to the growth of long range correlations and vanishes as we get deeper into the superfluid regime.

3.4. Numerical evaluation of the upper bound

In this section we will show how to compute the corrections to the local energy, Δεest , from Eq. (31). We basically need a method to compute expectation values of the operator O -2 around product states which have the form

ϕ=h˜1h˜k1Ah˜kBh˜k+1h˜k+2h˜M,
(35)

in which at most two contiguous vectors are affected by single-connection operators. For instance, this is the case of the optimal variational state, |ψ〉=O -1|〉⊗ M , where A and B are just the identity, but also of H 2|ψ〉, where A and B are ∑+, ∑- or ∑z. After some manipulations it is possible to write

O2H2ψ=k=1M2utO(HO)k1ZkO(HO)Mk1u,

Zk=J(HzOH+HOHz)U(HzOHz),

with the real matrices and vectors

Hij=h˜i2δij,

(Hα)ij=h˜i(Σαh˜)i2δij,

Oij={[(ij+n¯)!]12,ijn¯0,ij<n¯,

ui=δi0,

i,j,α{+,,z}..

We have used this technique to compute numerically the correction Δεest for different lattice sizes and found small or no differences for more than 30 sites. Intuitively, this is because in the limit of large powers the matrices (HO) k become projectors on the eigenvector with the largest eigenvalue. Using the same type of expansion we may compute other correlations

ak+Δakψ=utO(HO)k1(HO)Δ(HO)Mk+1uutO(HO)Mu.
(36)

For large values of Δ and large lattices, this quotient will decay exponentially as γΔ1 , where γ 1 is the largest eigenvalue of the matrix H_O.

4. Conclusions

In this work we have studied analytically and numerically the properties of the ground state of an ensemble of bosonic atoms in an 1D optical lattice. For the study of the atomic ensemble we have used both the quantum rotor model and the Bose-Hubbard model. Exploiting the fact that in these models there exists only nearest neighbor hopping and local interactions, we have developed a variational wavefunction that may be used to easily estimate local properties, such as the energy per well, the nearest neigbor correlations and the fluctuations of the density. In the case of the quantum rotor model we have verified our results with perturbative calculations around the strongly interacting regime, and with a spin wave approximation around the superfluid regime. In the case of the Bose-Hubbard model we have compared the variational estimates with numerical results obtained using the DMRG technique for a maximum density of three atoms per well. We have concluded that this procedure leads to fairly good estimates of local ground state properties of both Hamiltonians, in both the superfluid and the insulator regime, the largest disagreement being localized around the phase transition. On the other hand, we have also shown that our variational ansatz fails to describe long range properties of the superfluid phase, such as the algebraic decay of the first order correlation function with respect to distance.

Fig. 5. The energy of the product ansatz contains a contribution from each connection, εest , plus the interaction between neighbouring connections, Δεest . In Fig. (a) we show that Δεest (dash) is actually negative, and improves the estimate εest moving it towards the exact value, εDMRG (circles). Everything has been computed for =1. In Fig. (b) we show that the correction Δεest does not change much for large lattices.

Acknowledgments

We thank W. Zwerger for fruitful discussions. J. J. García-Ripoll and J. I. Cirac thank the Deutsche Forschungsgemeinschaft (SFB 631) and the Kompetenznetzwerk Quanteninformationsverarbeitung der Bayerischen Staatsregierung. C. K. and U. S. thank the Hess-Preis of the DFG and the Studienstiftung des deutschen Volkes for financial support. Work at the University of Innsbruck is supported by the Austrian Science Foundation, EU Networks and the Institute for Quantum Information. P.Z. thanks the Max Planck Institut für Quantenoptik for hospitality during his stay in Garching, and thanks the Humboldt Foundation for support during this stay.

References and links

1.

M. Greiner, O. Mandel, T. Esslinger, T. W. Hänsch, and I. Bloch, “Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms,” Nature 415, 39 (2002). [CrossRef] [PubMed]

2.

M. Greiner, O. Mandel, T. W. Hänsch, and Immanuel Bloch, “Collapse and revival of the matter wave field of a Bose-Einstein condensate,” Nature 419, 51 (2002). [CrossRef] [PubMed]

3.

D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, “Cold Bosonic Atoms in Optical Lattices,” Phys. Rev. Lett. 81, 3108–311 (1998); [CrossRef]

4.

D. Jaksch, V. Venturi, J. I. Cirac, C. J. Williams, and P. Zoller, “Creation of a Molecular Condensate by Dynamically Melting a Mott Insulator,” Phys. Rev. Lett. 89, 040402 (2002) [CrossRef] [PubMed]

5.

B. Y. Chen, S. D. Mahanti, and M. Yussouff, “Helium atoms in zeolite cages: Novel Mott-Hubbard and Bose-Hubbard systems,” Phys. Rev. Lett. 75, 473–476 (1995). [CrossRef] [PubMed]

6.

Min-Chul Cha, M. P. A. Fischer, S. M. Girvin, M. Wallin, and A. P. Young, “Universal conductivity of two-dimensional films at the superconductor-insulator transition,” Phys. Rev. B 44, 6883–6902 (1991). [CrossRef]

7.

M. P. A. Fischer, G. Grinstein, and S. M. Girvin “Presence of quantum diffusion in two dimensions: Universal resistance at the superconductor-insulator transition,” Phys. Rev. Lett.64, 587–590 (1990);M. P. A. Fischer, “Quantum phase transitions in disordered two-dimensional superconductors,” Phys. Rev. Lett.65, 923–926 (1990). [CrossRef]

8.

S. Sachdev, Quantum Phase Transitions, (Cambridge Univ. Press, Cambridge, 2001).

9.

M. P. A. Fischer, P. B. Weichman, G. Grinstein, and D. S. Fisher, “Boson localization and the superfluid-insulator transition,” Phys. Rev. B 40, 546–570 (1989). [CrossRef]

10.

J. K. Freericks and H. Monien, “Phase diagram of the Bose Hubbard model,” Europhys. Lett. 26, 545–550 (1994). [CrossRef]

11.

J. K. Freericks and H. Monien, “Strong-coupling expansions for the pure and disordered Bose-Hubbard model,” Phys. Rev. B 53, 2691–2700 (1996). [CrossRef]

12.

D. van Oosten, P. van der Straten, and H. T. C. Stoof, “Quantum phases in an optical lattice,” Phys. Rev. A 63, 053601 (2001). [CrossRef]

13.

A. M. Rey, K. Burnett, R. Roth, M. Edwards, C. J. Williams, and C. W. Clarck, “Bogoliubov Approach to Super-fluidity of Atoms in an Optical Lattices,” J. Phys. B , 36, 825–841 (2003). [CrossRef]

14.

R. Fazio and H. van der Zant, ‘Quantum phase transitions and vortex dynamics in superconducting networks,” Phys. Rep. 355, 235 (2001) and ref. therein. [CrossRef]

15.

A. van Otterlo, K.-H. Wagenblast, R. Baltin, C. Bruder, R. Fazio, and G. Schön, “Quantum phase transitions of interacting bosons and the supersolid phase,” Phys. Rev. B 52, 16176–16186 (1995). [CrossRef]

16.

M. P. A. Fisher and G. Grinstein, “Quantum Critical Phenomena in Charged Superconductors,” Phys. Rev. Lett. 60, 208–211 (1988). [CrossRef] [PubMed]

17.

S. R. White, “Density matrix formulation for quantum renormalization groups,” Phys. Rev. Lett. 69, 2863–2866 (1992). [CrossRef] [PubMed]

18.

S. R. White, “Density-matrix algorithms for quantum renormalization groups,” Phys. Rev. B. 48, 10345–10356 (1993). [CrossRef]

19.

I. Peschel, X. Wang, M. Kaulke, and K. Hallberg, Density-matrix Renormalization, (Springer-Verlag, Berlin, 1998).

20.

T. D. Kühner, S. R. White, and H. Monien, “One-dimensional Bose-Hubbard model with nearest-neighbor interaction,” Phys. Rev. B 61, 12474–12489 (2000). [CrossRef]

21.

T. D. Kühner, Diploma work (1997), University of Bonn.

22.

S. Rapsch, U. Schollwöck, and W. Zwerger, “Density matrix renormalization group for disordered bosons in one dimension,” Europhys. Lett. 46, 559–564 (1999). [CrossRef]

23.

G. G. Batrouni, R. T. Scalettar, and G. T. Zimanyi, “Supersolids in the Bose-Hubbard Hamiltonian,” Phys. Rev. Lett. 74, 2527–2530 (1995). [CrossRef] [PubMed]

24.

G. G. Batrouni and R. T. Scalettar, “World-line quantum Monte Carlo algorithm for a one-dimensional Bose model,” Phys. Rev. B 46, 9051 (1992). [CrossRef]

25.

P. Niyaz, R. T. Scalettar, C. Y. Fong, and G. G. Batrouni, “Phase transitions in an interacting boson model with near-neighbor repulsion,” Phys. Rev. B 50, 362–373 (1994). [CrossRef]

26.

N. V. Prokof’ev, B. V. Svistunov, and I. S. Tupitsyn, Phys. Lett. A 238, 253 (1998). [CrossRef]

27.

N. Elstner and H. Monien, “Dynamics and thermodynamics of the Bose-Hubbard model,” Phys. Rev. B 59, 12184–12187 (1999) and ref. therein. [CrossRef]

28.

D. S. Rokhsar and B. G. Kotliar, “Gutzwiller projection for bosons,” Phys. Rev. B 44, 10328–10332 (1991). [CrossRef]

29.

J. R. Anglin, P. Drummond, and A. Smerzi, “Exact quantum phase model for mesoscopic Josephson junctions,” Phys. Rev. A 64, 063605 (2001). [CrossRef]

OCIS Codes
(000.2690) General : General physics
(000.3860) General : Mathematical methods in physics

ToC Category:
Focus Issue: Cold atomic gases in optical lattices

History
Original Manuscript: November 7, 2003
Revised Manuscript: December 1, 2003
Published: January 12, 2004

Citation
Juan Jose García-Ripoll, J. Cirac, P. Zoller, C. Kollath, U. Schollwöck, and J. von Delft, "Variational ansatz for the superfluid Mott-insulator transition in optical lattices," Opt. Express 12, 42-54 (2004)
http://www.opticsinfobase.org/oe/abstract.cfm?URI=oe-12-1-42


Sort:  Journal  |  Reset  

References

  1. M. Greiner, O. Mandel, T. Esslinger, T. W. Hänsch, I. Bloch, �??Quantum phase transition from a superfluid to a Mott insulator in a gas of ultracold atoms,�?? Nature 415, 39 (2002). [CrossRef] [PubMed]
  2. M. Greiner, O. Mandel, T. W. Hänsch and Immanuel Bloch, �??Collapse and revival of the matter wave field of a Bose�??Einstein condensate,�?? Nature 419, 51 (2002). [CrossRef] [PubMed]
  3. D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner and P. Zoller, �??Cold Bosonic Atoms in Optical Lattices,�?? Phys. Rev. Lett. 81, 3108-311 (1998). [CrossRef]
  4. D. Jaksch, V. Venturi, J. I. Cirac, C. J. Williams and P. Zoller, �??Creation of a Molecular Condensate by Dynamically Melting a Mott Insulator,�?? Phys. Rev. Lett. 89, 040402 (2002). [CrossRef] [PubMed]
  5. B. Y. Chen, S. D. Mahanti and M. Yussouff, �??Helium atoms in zeolite cages: Novel Mott-Hubbard and Bose-Hubbard systems,�?? Phys. Rev. Lett. 75, 473-476 (1995). [CrossRef] [PubMed]
  6. Min-Chul Cha, M. P. A. Fischer, S. M. Girvin, M. Wallin, A. P. Young, �??Universal conductivity of two-dimensional films at the superconductor-insulator transition,�?? Phys. Rev. B 44, 6883-6902 (1991). [CrossRef]
  7. M. P. A. Fischer, G. Grinstein, S. M. Girvin, �??Presence of quantum diffusion in two dimensions: Universal resistance at the superconductor-insulator transition,�?? Phys. Rev. Lett. 64, 587-590 (1990); M. P. A. Fischer, �??Quantum phase transitions in disordered two-dimensional superconductors,�?? Phys. Rev. Lett. 65, 923-926 (1990). [CrossRef]
  8. S. Sachdev, Quantum Phase Transitions, (Cambridge Univ. Press, Cambridge, 2001).
  9. M. P. A. Fischer, P. B. Weichman, G. Grinstein, D. S. Fisher, �??Boson localization and the superfluid-insulator transition,�?? Phys. Rev. B 40, 546-570 (1989). [CrossRef]
  10. J. K. Freericks and H. Monien, �??Phase diagram of the Bose Hubbard model,�?? Europhys. Lett. 26, 545-550 (1994). [CrossRef]
  11. J. K. Freericks and H. Monien, �??Strong-coupling expansions for the pure and disordered Bose-Hubbard model,�?? Phys. Rev. B 53, 2691-2700 (1996). [CrossRef]
  12. D. van Oosten, P. van der Straten and H. T. C. Stoof, �??Quantum phases in an optical lattice,�?? Phys. Rev. A 63, 053601 (2001). [CrossRef]
  13. A. M. Rey, K. Burnett, R. Roth, M. Edwards, C. J. Williams and C. W. Clarck, �??Bogoliubov Approach to Super-fluidity of Atoms in an Optical Lattices,�?? J. Phys. B, 36, 825-841 (2003). [CrossRef]
  14. R. Fazio and H. van der Zant, �??Quantum phase transitions and vortex dynamics in superconducting networks,�?? Phys. Rep. 355, 235 (2001) and ref. therein. [CrossRef]
  15. A. van Otterlo, K.-H. Wagenblast, R. Baltin, C. Bruder, R. Fazio and G. Schön, �??Quantum phase transitions of interacting bosons and the supersolid phase,�?? Phys. Rev. B 52, 16176-16186 (1995). [CrossRef]
  16. M. P. A. Fisher and G. Grinstein, �??Quantum Critical Phenomena in Charged Superconductors,�?? Phys. Rev. Lett. 60, 208-211 (1988). [CrossRef] [PubMed]
  17. S. R. White, �??Density matrix formulation for quantum renormalization groups,�?? Phys. Rev. Lett. 69, 2863-2866 (1992). [CrossRef] [PubMed]
  18. S. R. White, �??Density-matrix algorithms for quantum renormalization groups,�?? Phys. Rev. B. 48, 10345-10356 (1993). [CrossRef]
  19. I. Peschel, X. Wang, M. Kaulke and K. Hallberg, Density-matrix Renormalization, (Springer-Verlag, Berlin, 1998).
  20. T. D. Kühner and S. R. White and H. Monien, �??One-dimensional Bose-Hubbard model with nearest-neighbor interaction,�?? Phys. Rev. B 61, 12474-12489 (2000). [CrossRef]
  21. T. D. Kühner, Diploma work (1997), University of Bonn.
  22. S. Rapsch, U. Schollwöck and W. Zwerger, �??Density matrix renormalization group for disordered bosons in one dimension,�?? Europhys. Lett. 46, 559-564 (1999). [CrossRef]
  23. G. G. Batrouni, R. T. Scalettar and G. T. Zimanyi, �??Supersolids in the Bose-Hubbard Hamiltonian,�?? Phys. Rev. Lett. 74, 2527-2530 (1995). [CrossRef] [PubMed]
  24. G. G. Batrouni and R. T. Scalettar, �??World-line quantum Monte Carlo algorithm for a one-dimensional Bose model,�?? Phys. Rev. B 46, 9051 (1992). [CrossRef]
  25. P. Niyaz, R. T. Scalettar, C. Y. Fong and G. G. Batrouni, �??Phase transitions in an interacting boson model with near-neighbor repulsion,�?? Phys. Rev. B 50, 362-373 (1994). [CrossRef]
  26. N. V. Prokof�??ev, B. V. Svistunov and I. S.Tupitsyn, Phys. Lett. A 238, 253 (1998). [CrossRef]
  27. N. Elstner and H. Monien, �??Dynamics and thermodynamics of the Bose-Hubbard model,�?? Phys. Rev. B 59, 12184-12187 (1999) and ref. therein. [CrossRef]
  28. D. S. Rokhsar and B. G. Kotliar, �??Gutzwiller projection for bosons,�?? Phys. Rev. B 44, 10328-10332 (1991). [CrossRef]
  29. J. R. Anglin, P. Drummond and A. Smerzi, �??Exact quantum phase model for mesoscopic Josephson junctions,�?? Phys. Rev. A 64, 063605 (2001). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.


« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited