Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1111.5085 · Physical Review B

Projection operator approach to the Bose-Hubbard model

Anirban Dutta, C. Trefzger, K. Sengupta

Published Aug 24, 2012• 35 citations
arXivHTMLPDF

Abstract

We develop a projection operator formalism for studying both the zero-temperature equilibrium phase diagram and the nonequilibrium dynamics of the Bose-Hubbard model. Our work, which constitutes an extension of that of Trefzger and Sengupta [Phys. Rev. Lett. 106, 095702 (2011)], shows that the method provides an accurate description of the equilibrium zero-temperature phase diagram of the Bose-Hubbard model for several lattices in two and three dimensions. We show that the accuracy of this method increases with the coordination number ${z}_{0}$ of the lattice and reaches to within $0.5%$ of quantum Monte Carlo data for lattices with ${z}_{0}=6$. We compute the excitation spectra of the bosons using this method in the Mott and the superfluid phases and compare our results with mean-field theory. We also show that the same method may be used to analyze the nonequilibrium dynamics of the model both in the Mott phase and near the superfluid-insulator quantum critical point where the hopping amplitude $J$ and the on-site interaction $U$ satisfy ${z}_{0}J/U\ensuremath{\ll}1$. In particular, we study the nonequilibrium dynamics of the model both subsequent to a sudden quench of the hopping amplitude $J$ and during a ramp from ${J}_{i}$ to ${J}_{f}$ characterized by a ramp time $\ensuremath{\tau}$ and exponent $\ensuremath{\alpha}$: $J(t)={J}_{i}+({J}_{f}\ensuremath{-}{J}_{i}){(t/\ensuremath{\tau})}^{\ensuremath{\alpha}}$. We compute the wave function overlap $F$, the residual energy $Q$, the superfluid order parameter $\ensuremath{\Delta}(t)$, the equal-time order parameter correlation function $C(t)$, and the defect formation probability $P$ for the above-mentioned protocols and provide a comparison of our results to their mean-field counterparts. We find that $Q$, $F$, and $P$ do not exhibit the expected universal scaling. We explain this absence of universality and show that our results for linear ramps compare well with the recent experimental observations.

DOI
10.1103/physrevb.86.085140

Related papers

Nonequilibrium Dynamics of the Bose-Hubbard Model: A Projection-Operator Approach

C. Trefzger, K. Sengupta

We study the phase diagram and nonequilibrium dynamics involving ramp of the hopping amplitude J(t)=Jt/τ with ramp time τ of the Bose-Hubbard model at zero temperature using a projection-operator formalism which allows us to incorporate the effects of quantum fluctuations beyond mean-field approximations in the strong-coupling regime. Our formalism yields a phase diagram which provides a near exact match with quantum Monte Carlo results in three dimensions. We also compute the residual energy Q, the superfluid order parameter Δ(t), the equal-time order parameter correlation function C(t), and the wave function overlap F which yields the defect formation probability P during nonequilibrium dynamics of the model. We find that Q, F, and P do not exhibit the expected universal scaling. We explain this absence of universality and show that our results compare well with recent experiments.

arXiv:1008.1285•Physical Review Letters•Feb 28, 2011•65 citations

A projection operator approach for ultracold bosons coupled to a cavity

Sarbarish Sen, Bimalendu Deb, K. Sengupta

We employ a projection operator technique to study the equilibrium phases and the non-equilibrium quench and ramp dynamics of ultracold bosons described by a two-dimensional (2D) Bose-Hubbard model coupled to a high-finesse cavity with competing local short-range and cavity-mediated long-range interactions. Our analysis relies on a systematic elimination of high-energy virtual hopping processes; this enables us to capture the effects of short-range quantum fluctuations on the phase diagram of the model while treating the cavity-induced long-range interaction within mean-field theory. We find that the phase boundaries between the superfluid (SF), supersolid (SS), charge-density-wave (CDW) and Mott insulating (MI) phases are significantly modified by the presence of these fluctuations. We use this approach to study the dynamics of bosons following a quench/ramp from the CDW to the SS and SF phases; our analysis indicates the presence of oscillatory dynamics characterized by a single frequency due to a quench which takes the system near the critical point. In contrast, for a quench deep inside the SF phase, the dynamics involves multiple frequencies. An explicit comparison with Gutzwiller mean-field theory shows the importance of quantum fluctuations in shaping such post-quench dynamics; this feature may be tested experimentally as we discuss. We also study the ramp dynamics of these bosons and identify a bound on the ramp rate above which the present method is expected to produce accurate results.

arXiv:2609.34908•Sep 28, 2026•0 citations

Dynamical mean-field theory for the Bose-Hubbard model

Wen-Jun Hu, Ning-Hua Tong

The dynamical mean-field theory (DMFT), which is successful in the study of strongly correlated fermions, was recently extended to boson systems [K. Byczuk and D. Vollhardt, Phys. Rev. B 77, 235106 (2008)]. In this paper, we employ the bosonic DMFT to study the Bose-Hubbard model which describes on-site interacting bosons in a lattice. Using exact diagonalization as the impurity solver, we get the DMFT solutions for the Green's function, the occupation density, as well as the condensate fraction on a Bethe lattice. Various phases are identified: the Mott insulator, the Bose-Einstein condensed (BEC) phase, and the normal phase. At finite temperatures, we obtain the crossover between the Mott-type regime and the normal phase, as well as the BEC-to-normal phase transition. Phase diagrams on the $\ensuremath{\mu}/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane and on the $T/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane are produced ($\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}$ is the scaled hopping amplitude). We compare our results with the previous ones, and discuss the implication of these results to experiments.

arXiv:0907.2928•Physical Review B•Dec 11, 2009•49 citations

Quantum dynamics of disordered bosons in an optical lattice

Chien-Hung Lin, Rajdeep Sensarma, K. Sengupta +1 more

We study the equilibrium and nonequilibrium properties of strongly interacting bosons on a lattice in the presence of a random bounded disorder potential. Using a Gutzwiller projected variational technique, we study the equilibrium phase diagram of the disordered Bose-Hubbard model and obtain the Mott insulator, Bose glass, and superfluid phases. We also study the nonequilibrium response of the system under a periodic temporal drive where, starting from the superfluid phase, the hopping parameter is ramped down linearly in time, and back to its initial value. We study the density of excitations created, the change in the superfluid order parameter, and the energy pumped into the system in this process as a function of the inverse ramp rate $\ensuremath{\tau}$. For the clean case the density of excitations goes to a constant, while the order parameter and energy relax as $1/\ensuremath{\tau}$ and $1/{\ensuremath{\tau}}^{2}$ respectively. With disorder, the excitation density decays exponentially with $\ensuremath{\tau}$, with the decay rate increasing with the disorder, to an asymptotic value independent of the disorder. The energy and change in order parameter also decrease as $\ensuremath{\tau}$ is increased.

arXiv:1208.4852•Physical Review B•Dec 26, 2012•12 citations

Extended Bose Hubbard model of interacting bosonic atoms in optical lattices: From superfluidity to density waves

G. Mazzarella, S. M. Giampaolo, Fabrizio Illuminati

For systems of interacting, ultracold spin-zero neutral bosonic atoms, harmonically trapped and subject to an optical lattice potential, we derive an Extended Bose Hubbard (EBH) model by developing a systematic expansion for the Hamiltonian of the system in powers of the lattice parameters and of a scale parameter, the lattice attenuation factor. We identify the dominant terms that need to be retained in realistic experimental conditions, up to nearest-neighbor interactions and nearest-neighbor hoppings conditioned by the on-site occupation numbers. In the mean field approximation, we determine the free energy of the system and study the phase diagram both at zero and at finite temperature. At variance with the standard on site Bose Hubbard model, the zero-temperature phase diagram of the EBH model possesses a dual structure in the Mott insulating regime. Namely, for specific ranges of the lattice parameters, a density wave phase characterizes the system at integer fillings, with domains of alternating mean occupation numbers that are the atomic counterparts of the domains of staggered magnetizations in an antiferromagnetic phase. We show as well that in the EBH model, a zero-temperature quantum phase transition to pair superfluidity is, in principle, possible, but completely suppressed at the lowest order in the lattice attenuation factor. Finally, we determine the possible occurrence of the different phases as a function of the experimentally controllable lattice parameters.

arXiv:cond-mat/0509161•Physical Review A•Jan 26, 2006•76 citations

Exact solution of the Bose-Hubbard model on the Bethe lattice

Guilhem Semerjian, Marco Tarzia, Francesco Zamponi

The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.

arXiv:0904.3075•Physical Review B•Jul 31, 2009•42 citations
We develop a projection operator formalism for studying both the zero-temperature equilibrium phase diagram and the nonequilibrium dynamics of the Bose-Hubbard model. Our work, which constitutes an extension of that of Trefzger and Sengupta [Phys. Rev. Lett. 106, 095702 (2011)], shows that the method provides an accurate description of the equilibrium zero-temperature phase diagram of the Bose-Hubbard model for several lattices in two and three dimensions. We show that the accuracy of this method increases with the coordination number z0{z}_{0}z0​ of the lattice and reaches to within 0.50.5%0.5 of quantum Monte Carlo data for lattices with z0=6{z}_{0}=6z0​=6. We compute the excitation spectra of the bosons using this method in the Mott and the superfluid phases and compare our results with mean-field theory. We also show that the same method may be used to analyze the nonequilibrium dynamics of the model both in the Mott phase and near the superfluid-insulator quantum critical point where the hopping amplitude JJJ and the on-site interaction UUU satisfy z0J/U\ensuremath≪1{z}_{0}J/U\ensuremath{\ll}1z0​J/U\ensuremath≪1. In particular, we study the nonequilibrium dynamics of the model both subsequent to a sudden quench of the hopping amplitude JJJ and during a ramp from Ji{J}_{i}Ji​ to Jf{J}_{f}Jf​ characterized by a ramp time \ensuremathτ\ensuremath{\tau}\ensuremathτ and exponent \ensuremathα\ensuremath{\alpha}\ensuremathα: J(t)=Ji+(Jf\ensuremath−Ji)(t/\ensuremathτ)\ensuremathαJ(t)={J}_{i}+({J}_{f}\ensuremath{-}{J}_{i}){(t/\ensuremath{\tau})}^{\ensuremath{\alpha}}J(t)=Ji​+(Jf​\ensuremath−Ji​)(t/\ensuremathτ)\ensuremathα. We compute the wave function overlap FFF, the residual energy QQQ, the superfluid order parameter \ensuremathΔ(t)\ensuremath{\Delta}(t)\ensuremathΔ(t), the equal-time order parameter correlation function C(t)C(t)C(t), and the defect formation probability PPP for the above-mentioned protocols and provide a comparison of our results to their mean-field counterparts. We find that QQQ, FFF, and PPP do not exhibit the expected universal scaling. We explain this absence of universality and show that our results for linear ramps compare well with the recent experimental observations.
We study the phase diagram and nonequilibrium dynamics involving ramp of the hopping amplitude J(t)=Jt/τ with ramp time τ of the Bose-Hubbard model at zero temperature using a projection-operator formalism which allows us to incorporate the effects of quantum fluctuations beyond mean-field approximations in the strong-coupling regime. Our formalism yields a phase diagram which provides a near exact match with quantum Monte Carlo results in three dimensions. We also compute the residual energy Q, the superfluid order parameter Δ(t), the equal-time order parameter correlation function C(t), and the wave function overlap F which yields the defect formation probability P during nonequilibrium dynamics of the model. We find that Q, F, and P do not exhibit the expected universal scaling. We explain this absence of universality and show that our results compare well with recent experiments.
We employ a projection operator technique to study the equilibrium phases and the non-equilibrium quench and ramp dynamics of ultracold bosons described by a two-dimensional (2D) Bose-Hubbard model coupled to a high-finesse cavity with competing local short-range and cavity-mediated long-range interactions. Our analysis relies on a systematic elimination of high-energy virtual hopping processes; this enables us to capture the effects of short-range quantum fluctuations on the phase diagram of the model while treating the cavity-induced long-range interaction within mean-field theory. We find that the phase boundaries between the superfluid (SF), supersolid (SS), charge-density-wave (CDW) and Mott insulating (MI) phases are significantly modified by the presence of these fluctuations. We use this approach to study the dynamics of bosons following a quench/ramp from the CDW to the SS and SF phases; our analysis indicates the presence of oscillatory dynamics characterized by a single frequency due to a quench which takes the system near the critical point. In contrast, for a quench deep inside the SF phase, the dynamics involves multiple frequencies. An explicit comparison with Gutzwiller mean-field theory shows the importance of quantum fluctuations in shaping such post-quench dynamics; this feature may be tested experimentally as we discuss. We also study the ramp dynamics of these bosons and identify a bound on the ramp rate above which the present method is expected to produce accurate results.
The dynamical mean-field theory (DMFT), which is successful in the study of strongly correlated fermions, was recently extended to boson systems [K. Byczuk and D. Vollhardt, Phys. Rev. B 77, 235106 (2008)]. In this paper, we employ the bosonic DMFT to study the Bose-Hubbard model which describes on-site interacting bosons in a lattice. Using exact diagonalization as the impurity solver, we get the DMFT solutions for the Green's function, the occupation density, as well as the condensate fraction on a Bethe lattice. Various phases are identified: the Mott insulator, the Bose-Einstein condensed (BEC) phase, and the normal phase. At finite temperatures, we obtain the crossover between the Mott-type regime and the normal phase, as well as the BEC-to-normal phase transition. Phase diagrams on the $\ensuremath{\mu}/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane and on the $T/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane are produced ($\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}$ is the scaled hopping amplitude). We compare our results with the previous ones, and discuss the implication of these results to experiments.
We study the equilibrium and nonequilibrium properties of strongly interacting bosons on a lattice in the presence of a random bounded disorder potential. Using a Gutzwiller projected variational technique, we study the equilibrium phase diagram of the disordered Bose-Hubbard model and obtain the Mott insulator, Bose glass, and superfluid phases. We also study the nonequilibrium response of the system under a periodic temporal drive where, starting from the superfluid phase, the hopping parameter is ramped down linearly in time, and back to its initial value. We study the density of excitations created, the change in the superfluid order parameter, and the energy pumped into the system in this process as a function of the inverse ramp rate $\ensuremath{\tau}$. For the clean case the density of excitations goes to a constant, while the order parameter and energy relax as $1/\ensuremath{\tau}$ and $1/{\ensuremath{\tau}}^{2}$ respectively. With disorder, the excitation density decays exponentially with $\ensuremath{\tau}$, with the decay rate increasing with the disorder, to an asymptotic value independent of the disorder. The energy and change in order parameter also decrease as $\ensuremath{\tau}$ is increased.
For systems of interacting, ultracold spin-zero neutral bosonic atoms, harmonically trapped and subject to an optical lattice potential, we derive an Extended Bose Hubbard (EBH) model by developing a systematic expansion for the Hamiltonian of the system in powers of the lattice parameters and of a scale parameter, the lattice attenuation factor. We identify the dominant terms that need to be retained in realistic experimental conditions, up to nearest-neighbor interactions and nearest-neighbor hoppings conditioned by the on-site occupation numbers. In the mean field approximation, we determine the free energy of the system and study the phase diagram both at zero and at finite temperature. At variance with the standard on site Bose Hubbard model, the zero-temperature phase diagram of the EBH model possesses a dual structure in the Mott insulating regime. Namely, for specific ranges of the lattice parameters, a density wave phase characterizes the system at integer fillings, with domains of alternating mean occupation numbers that are the atomic counterparts of the domains of staggered magnetizations in an antiferromagnetic phase. We show as well that in the EBH model, a zero-temperature quantum phase transition to pair superfluidity is, in principle, possible, but completely suppressed at the lowest order in the lattice attenuation factor. Finally, we determine the possible occurrence of the different phases as a function of the experimentally controllable lattice parameters.
The exact solution of a quantum Bethe lattice model in the thermodynamic limit amounts to solve a functional self-consistent equation. In this paper we obtain this equation for the Bose-Hubbard model on the Bethe lattice, under two equivalent forms. The first one, based on a coherent-state path integral, leads in the large connectivity limit to the mean-field treatment of Fisher et al. [Phys. Rev. B 40, 546 (1989)] at the leading order, and to the bosonic dynamical mean field theory as a first correction, as recently derived by Byczuk and Vollhardt [Phys. Rev. B 77, 235106 (2008)]. We obtain an alternative form of the equation using the occupation number representation, which can be easily solved with an arbitrary numerical precision, for any finite connectivity. We thus compute the transition line between the superfluid and Mott insulator phases of the model, along with thermodynamic observables and the space and imaginary-time dependence of correlation functions. The finite connectivity of the Bethe lattice induces a richer physical content with respect to its infinitely connected counterpart: a notion of distance between sites of the lattice is preserved, and the bosons are still weakly mobile in the Mott insulator phase. The Bethe lattice construction can be viewed as an approximation to the finite-dimensional version of the model. We show indeed a quantitatively reasonable agreement between our predictions and the results of Quantum Monte Carlo simulations in two and three dimensions.