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

Physical Review A

Quantum dynamics of Bose-Hubbard Hamiltonians beyond the Hartree-Fock-Bogoliubov approximation: The Bogoliubov back-reaction approximation

Igor Tikhonenkov, J. R. Anglin, Amichay Vardi

Published Jan 10, 2007• 36 citations
Publisher

Abstract

We formulate a method for studying the quantum field dynamics of ultracold Bose gases confined within optical lattice potentials, within the lowest Bloch-band Bose-Hubbard model. Our formalism extends the two-sites results of Vardi and Anglin [Phys. Rev. Lett. 86, 568 (2001)] to the general case of $M$ lattice sites. The methodology is based on mapping the Bose-Hubbard Hamiltonian to an $\mathrm{SU}(M)$ pseudospin problem and truncating the resulting hierarchy of dynamical equations for correlation functions, up to pair correlations between $\mathrm{SU}(M)$ generators. Agreement with few-site exact many-particle calculations is consistently better than the corresponding Hartree-Fock-Bogoliubov approximation. Moreover, our approximation compares favorably with a more elaborate two-particle irreducible effective action formalism, at a fraction of the analytic and numerical effort.

DOI
10.1103/physreva.75.013613

Related papers

Nonequilibrium dynamics of condensates in a lattice with the two-particle-irreducible effective action in the 1∕N expansion

Kristan Temme, Thomas Gasenzer

The dynamical evolution of a Bose-Einstein condensate trapped in a one-dimensional lattice potential is investigated theoretically in the framework of the Bose-Hubbard model. The emphasis is set on the far-from-equilibrium evolution in a case where the gas is strongly interacting. This is realized by an appropriate choice of the parameters in the Hamiltonian, and by starting with an initial state, where one lattice well contains a Bose-Einstein condensate while all other wells are empty. Oscillations of the condensate as well as noncondensate fractions of the gas between the different sites of the lattice are found to be damped as a consequence of the collisional interactions between the atoms. Functional integral techniques involving self-consistently determined mean fields as well as two-point correlation functions are used to derive the two-particle-irreducible (2PI) effective action. The action is expanded in inverse powers of the number of field components $\mathcal{N}$, and the dynamic equations are derived from it to next-to-leading order in this expansion. This approach reaches considerably beyond the Hartree-Fock-Bogoliubov mean-field theory, and its results are compared to the exact quantum dynamics obtained by Rey et al. [Phys. Rev. A 69, 033610 (2004)] for small atom numbers.

arXiv:cond-mat/0607116•Physical Review A•Nov 2, 2006•16 citations

Nonequilibrium dynamics of optical-lattice-loaded Bose-Einstein-condensate atoms: Beyond the Hartree-Fock-Bogoliubov approximation

Ana María Rey, B. L. Hu, Esteban Calzetta +2 more

In this work a two-particle irreducible (2PI) closed-time-path (CTP) effective action is used to describe the nonequilibrium dynamics of a Bose-Einstein condensate selectively loaded into every third site of a one-dimensional optical lattice. The motivation of this work is the recent experimental realization of this system. Patterned loading methods may be useful for quantum computing with trapped atoms. This system also serves to illustrate many basic issues in nonequilibrium quantum-field theory pertaining to the dynamics of quantum correlations and fluctuations which goes beyond the capability of a mean-field theory. By numerically evolving in time the initial-state configuration using the Bose-Hubbard Hamiltonian an exact quantum solution is available for this system in the case of few atoms and wells. One can also use it to test various approximate methods. Under the 2PI CTP scheme with this initial configuration, three different approximations are considered: (a) the Hartree-Fock-Bogoliubov (HFB) approximation, (b) the next-to-leading-order $1∕\mathcal{N}$ expansion of the 2PI effective action up to second order in the interaction strength, and (c) a second-order perturbative expansion in the interaction strength. We present detailed comparisons between these approximations and determine their range of validity by contrasting them with the exact many-body solution for a moderate number of atoms and wells. As a general feature we observe that because the second-order 2PI approximations include multiparticle scattering in a systematic way, they are able to capture damping effects exhibited in the exact solution, which a mean-field collisionless approach fails to produce. While the second-order approximations show a clear improvement over the HFB approximation, our numerical results show that they fail at late times, when interaction effects are significant.

arXiv:cond-mat/0308305•Physical Review A•Mar 22, 2004•100 citations

Quantum Dynamics of Strongly Correlated Ultracold Bose Gases in Optical Lattices

Markus Hild

Ultracold bosonic gases in optical lattices are strongly correlated quantum systems similar to solids. The strong correlation between the electrons in a solid on the one hand, and the bosonic atoms in optical lattice on the other, exhibit various quantum phenomena like insulation, conductivity, localization of electrons and atoms, respectively. Controlled by the intensity of the lattice laser, the ultracold bosonic gas can be transferred from a regime with superfluid character for shallow lattices into a regime of strong correlations, the Mott insulator. As an additional external parameter besides the lattice depth, one can generate spatial inhomogeneities by superimposing an additional standing wave (so-called two-color superlattices), which gives rise to localization effects or the formation of a Bose-glass phase. In the present work, numerical simulations are employed in order to investigate characteristic signatures of the quantum phases in the low-energy excitation spectrum of one-dimensional systems. We simulate temporal small amplitude modulations of the optical lattice in analogy to experiments, and evaluate the response of the system from the time-evolved initial state. The lattice systems are described in the framework of the Bose-Hubbard model. For the evaluation of the time-evolved state, we employ several numerical methods. We analyze systems of small size (6 particles on 6 sites) using an exact time-evolution by integration of the time-dependent Schrödinger equation. The formulation of an importance truncation scheme enables us to retain only the relevant components of the model space in the strongly correlated regime and, thus, allows for the investigation of systems with 10 particles on 10 sites using exact time-evolution. Based on this method, we present results of the Mott-insulating regime as well as for the Bose-glass phase. Furthermore, we employ particle-hole methods, which allow for the treatment of systems with experimental lattice sizes and particle numbers. Starting from the equation of motion method we adapt the Tamm-Dancoff approximation as well as the random-phase approximation for the occupation number representation of the Bose-Hubbard model. We present results of simulations of up to 50 particles on 50 sites and discuss the impact of the lattice depth on the low-energy excitations (U-resonance). Moreover, the impact of a two-color superlattice and the variation of its amplitude is investigated.

•Jan 1, 2010•0 citations

Dynamics of Quantum Statistical Correlations in Ultracold Bose Gases

Cédric Bodet

The dynamical evolution of an ultracold Bose gas distributed across the sites of an optical lattice is investigated theoretically in the framework of the Bose-Hubbard model. First, the focus is set on the evolution of squeezing correlations in the two mode system. It is shown that the eigenstates of the Hamiltonian do not exploit the full region of possible squeezing allowed by Heisenbergs uncertainty relation for number and phase fl uctuations. The development of nonclassical correlations and relative number squeezing is studied at the transition from the Josephson to the Fock regime. Comparing the full quantum evolution with classical statistical simulations allows us to identify quantum aspects of the squeezing formation. In the quantum regime, the measurement of squeezing allows us to distinguish even and odd total particle number states. Then, a far from equilibrium quantum field theory method, the so-called two-particle-irreducible effective action approach, is presented for the description of the dynamics in larger lattices. The resulting dynamics is compared to the classical statistical time evolution. The validity of the quantum field evolution is probed for various initial conditions in the classical regime.

•Jan 1, 2011•0 citations

Finite-temperature theory of superfluid bosons in optical lattices

D. Baillie, P. B. Blakie

A practical finite-temperature theory is developed for the superfluid regime of a weakly interacting Bose gas in an optical lattice with additional harmonic confinement. We derive an extended Bose-Hubbard model that is valid for shallow lattices and when excited bands are occupied. Using the Hartree-Fock-Bogoliubov-Popov mean-field approach, and applying local-density and coarse-grained envelope approximations, we arrive at a theory that can be numerically implemented accurately and efficiently. We present results for a three-dimensional system, characterizing the importance of the features of the extended Bose-Hubbard model and compare against other theoretical results and show an improved agreement with experimental data.

arXiv:0906.4606•Physical Review A•Sep 21, 2009•24 citations

Bose-Hubbard ground state: Extended Bogoliubov and variational methods compared with time-evolving block decimation

Ippei Danshita, Pascal Naidon

We determine the ground-state properties of a gas of interacting bosonic atoms in a one-dimensional optical lattice. The system is modeled by the Bose-Hubbard Hamiltonian. We show how to apply the time-evolving block decimation method to systems with periodic boundary conditions and employ it as a reference to find the ground state of the Bose-Hubbard model. Results are compared with recently proposed approximate methods, such as Hartree-Fock-Bogoliubov (HFB) theories generalized for strong interactions and the variational Bijl-Dingle-Jastrow method. We find that all HFB methods do not bring any improvement to the Bogoliubov theory and therefore provide correct results only in the weakly interacting limit, where the system is deeply in the superfluid regime. On the other hand, the variational Bijl-Dingle-Jastrow method is applicable for much stronger interactions but is essentially limited to the superfluid regime as it reproduces the superfluid--Mott-insulator transition only qualitatively.

arXiv:0809.1486•Physical Review A•Apr 2, 2009•24 citations
We formulate a method for studying the quantum field dynamics of ultracold Bose gases confined within optical lattice potentials, within the lowest Bloch-band Bose-Hubbard model. Our formalism extends the two-sites results of Vardi and Anglin [Phys. Rev. Lett. 86, 568 (2001)] to the general case of MMM lattice sites. The methodology is based on mapping the Bose-Hubbard Hamiltonian to an SU(M)\mathrm{SU}(M)SU(M) pseudospin problem and truncating the resulting hierarchy of dynamical equations for correlation functions, up to pair correlations between SU(M)\mathrm{SU}(M)SU(M) generators. Agreement with few-site exact many-particle calculations is consistently better than the corresponding Hartree-Fock-Bogoliubov approximation. Moreover, our approximation compares favorably with a more elaborate two-particle irreducible effective action formalism, at a fraction of the analytic and numerical effort.
The dynamical evolution of a Bose-Einstein condensate trapped in a one-dimensional lattice potential is investigated theoretically in the framework of the Bose-Hubbard model. The emphasis is set on the far-from-equilibrium evolution in a case where the gas is strongly interacting. This is realized by an appropriate choice of the parameters in the Hamiltonian, and by starting with an initial state, where one lattice well contains a Bose-Einstein condensate while all other wells are empty. Oscillations of the condensate as well as noncondensate fractions of the gas between the different sites of the lattice are found to be damped as a consequence of the collisional interactions between the atoms. Functional integral techniques involving self-consistently determined mean fields as well as two-point correlation functions are used to derive the two-particle-irreducible (2PI) effective action. The action is expanded in inverse powers of the number of field components N\mathcal{N}N, and the dynamic equations are derived from it to next-to-leading order in this expansion. This approach reaches considerably beyond the Hartree-Fock-Bogoliubov mean-field theory, and its results are compared to the exact quantum dynamics obtained by Rey et al. [Phys. Rev. A 69, 033610 (2004)] for small atom numbers.
In this work a two-particle irreducible (2PI) closed-time-path (CTP) effective action is used to describe the nonequilibrium dynamics of a Bose-Einstein condensate selectively loaded into every third site of a one-dimensional optical lattice. The motivation of this work is the recent experimental realization of this system. Patterned loading methods may be useful for quantum computing with trapped atoms. This system also serves to illustrate many basic issues in nonequilibrium quantum-field theory pertaining to the dynamics of quantum correlations and fluctuations which goes beyond the capability of a mean-field theory. By numerically evolving in time the initial-state configuration using the Bose-Hubbard Hamiltonian an exact quantum solution is available for this system in the case of few atoms and wells. One can also use it to test various approximate methods. Under the 2PI CTP scheme with this initial configuration, three different approximations are considered: (a) the Hartree-Fock-Bogoliubov (HFB) approximation, (b) the next-to-leading-order 1∕N1∕\mathcal{N}1∕N expansion of the 2PI effective action up to second order in the interaction strength, and (c) a second-order perturbative expansion in the interaction strength. We present detailed comparisons between these approximations and determine their range of validity by contrasting them with the exact many-body solution for a moderate number of atoms and wells. As a general feature we observe that because the second-order 2PI approximations include multiparticle scattering in a systematic way, they are able to capture damping effects exhibited in the exact solution, which a mean-field collisionless approach fails to produce. While the second-order approximations show a clear improvement over the HFB approximation, our numerical results show that they fail at late times, when interaction effects are significant.
Ultracold bosonic gases in optical lattices are strongly correlated quantum systems similar to solids. The strong correlation between the electrons in a solid on the one hand, and the bosonic atoms in optical lattice on the other, exhibit various quantum phenomena like insulation, conductivity, localization of electrons and atoms, respectively. Controlled by the intensity of the lattice laser, the ultracold bosonic gas can be transferred from a regime with superfluid character for shallow lattices into a regime of strong correlations, the Mott insulator. As an additional external parameter besides the lattice depth, one can generate spatial inhomogeneities by superimposing an additional standing wave (so-called two-color superlattices), which gives rise to localization effects or the formation of a Bose-glass phase. In the present work, numerical simulations are employed in order to investigate characteristic signatures of the quantum phases in the low-energy excitation spectrum of one-dimensional systems. We simulate temporal small amplitude modulations of the optical lattice in analogy to experiments, and evaluate the response of the system from the time-evolved initial state. The lattice systems are described in the framework of the Bose-Hubbard model. For the evaluation of the time-evolved state, we employ several numerical methods. We analyze systems of small size (6 particles on 6 sites) using an exact time-evolution by integration of the time-dependent Schrödinger equation. The formulation of an importance truncation scheme enables us to retain only the relevant components of the model space in the strongly correlated regime and, thus, allows for the investigation of systems with 10 particles on 10 sites using exact time-evolution. Based on this method, we present results of the Mott-insulating regime as well as for the Bose-glass phase. Furthermore, we employ particle-hole methods, which allow for the treatment of systems with experimental lattice sizes and particle numbers. Starting from the equation of motion method we adapt the Tamm-Dancoff approximation as well as the random-phase approximation for the occupation number representation of the Bose-Hubbard model. We present results of simulations of up to 50 particles on 50 sites and discuss the impact of the lattice depth on the low-energy excitations (U-resonance). Moreover, the impact of a two-color superlattice and the variation of its amplitude is investigated.
The dynamical evolution of an ultracold Bose gas distributed across the sites of an optical lattice is investigated theoretically in the framework of the Bose-Hubbard model. First, the focus is set on the evolution of squeezing correlations in the two mode system. It is shown that the eigenstates of the Hamiltonian do not exploit the full region of possible squeezing allowed by Heisenbergs uncertainty relation for number and phase fl uctuations. The development of nonclassical correlations and relative number squeezing is studied at the transition from the Josephson to the Fock regime. Comparing the full quantum evolution with classical statistical simulations allows us to identify quantum aspects of the squeezing formation. In the quantum regime, the measurement of squeezing allows us to distinguish even and odd total particle number states. Then, a far from equilibrium quantum field theory method, the so-called two-particle-irreducible effective action approach, is presented for the description of the dynamics in larger lattices. The resulting dynamics is compared to the classical statistical time evolution. The validity of the quantum field evolution is probed for various initial conditions in the classical regime.
A practical finite-temperature theory is developed for the superfluid regime of a weakly interacting Bose gas in an optical lattice with additional harmonic confinement. We derive an extended Bose-Hubbard model that is valid for shallow lattices and when excited bands are occupied. Using the Hartree-Fock-Bogoliubov-Popov mean-field approach, and applying local-density and coarse-grained envelope approximations, we arrive at a theory that can be numerically implemented accurately and efficiently. We present results for a three-dimensional system, characterizing the importance of the features of the extended Bose-Hubbard model and compare against other theoretical results and show an improved agreement with experimental data.
We determine the ground-state properties of a gas of interacting bosonic atoms in a one-dimensional optical lattice. The system is modeled by the Bose-Hubbard Hamiltonian. We show how to apply the time-evolving block decimation method to systems with periodic boundary conditions and employ it as a reference to find the ground state of the Bose-Hubbard model. Results are compared with recently proposed approximate methods, such as Hartree-Fock-Bogoliubov (HFB) theories generalized for strong interactions and the variational Bijl-Dingle-Jastrow method. We find that all HFB methods do not bring any improvement to the Bogoliubov theory and therefore provide correct results only in the weakly interacting limit, where the system is deeply in the superfluid regime. On the other hand, the variational Bijl-Dingle-Jastrow method is applicable for much stronger interactions but is essentially limited to the superfluid regime as it reproduces the superfluid--Mott-insulator transition only qualitatively.