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

Out-of-equilibrium dynamics of the Bose-Hubbard model in the strong coupling regime

Matthew Fitzpatrick

Published Jul 5, 2019• 1 citation
OpenAlex

Abstract

Experimental advances have made ultracold atoms in optical lattices a favourable setting to study out-of-equilibrium phenomena and attracted considerable attention in recent years.These systems are highly versatile in that experimental parameters can be tuned over a wide range of values in real time.This facilitates the study of quantum quenches, in which parameters in the corresponding Hamiltonian are varied in time faster than the system can respond adiabatically.Such protocols open the door to a rich range of many-body physics and have been studied intensely both theoretically and experimentally.The Bose-Hubbard model (BHM) has been shown to describe interacting ultracold bosons in an optical lattice, allowing the opportunity for experiments to probe the out-of-equilibrium dynamics of the model.The BHM is a particularly convenient context for studying quantum quenches as it displays a quantum phase transition between superfluid and Mott-insulator phases.In this thesis, we develop a strong-coupling approach that allows the study of correlations in space and time in both the superfluid and Mott-insulating phases of the BHM.Specifically, we obtain a two-particle irreducible (2PI) effective action within the contourtime formalism that allows for a description of both equilibrium and out-of-equilibrium phenomena.We derive equations of motion for both the superfluid order parameter and the single-particle many-body Green's functions.First, we assess the accuracy of our formalism by studying the equilibrium solution for the homogeneous BHM and comparing our results to existing strong-coupling methods as well as to exact methods where possible.We then consider homogeneous systems that are initially thermalized in the Mott phase, and which are then subjected to quenches.We solve numerically the equations of motion for this scenario and calculate the single-particle density matrix.We demonstrate a Lieb-Robinson-like maximal propagation velocity for the spreading of single-particle correlations in one, two, and three dimensions.We study the dependence of the maximal propagation velocity on the quench protocol, chemical potential, temperature, and dimensionality.We compare our results to exact methods, existing strong-coupling approaches, and experiments where possible.Lastly, we extend our strong-coupling approach to the disordered BHM and derive equations of motion for the disorder-averaged single-particle Green's function.We discuss how these equations of motion can be used to study the phase stability of many-body localization in the disordered BHM for dimensions higher than one.

Related papers

Contour-time approach to the disordered Bose-Hubbard model in the strong coupling regime

Ali Mokhtari-Jazi, Matthew Fitzpatrick, Malcolm P. Kennett

There has been considerable interest in the disordered Bose Hubbard model (BHM) in recent years, particularly in the context of thermalization and many-body localization. We develop a two-particle irreducible (2PI) strong-coupling approach to the disordered BHM that allows us to treat both equilibrium and out-of-equilibrium situations. We obtain equations of motion for spatio-temporal correlations and explore their equilibrium solutions. We study the equilibrium phase diagram as a function of disorder strength and discuss applications of the formalism to out-of-equilibrium situations. We also note that the disorder strengths where the emergence of non-ergodic dynamics was observed in a recent experiment [Choi et al. Science 352, 1547 (2016)] appear to correspond to the Mott insulator – Bose Glass phase boundary.

arXiv:2304.11260•Nuclear Physics B•Nov 8, 2023•5 citations

Contour-time approach to the Bose–Hubbard model in the strong coupling regime: Studying two-point spatio-temporal correlations at the Hartree–Fock–Bogoliubov level

Matthew Fitzpatrick, Malcolm P. Kennett

We develop a formalism that allows the study of correlations in space and time in both the superfluid and Mott insulating phases of the Bose–Hubbard Model. Specifically, we obtain a two particle irreducible effective action within the contour-time formalism that allows for both equilibrium and out of equilibrium phenomena. We derive equations of motion for both the superfluid order parameter and two-point correlation functions. To assess the accuracy of this formalism, we study the equilibrium solution of the equations of motion and compare our results to existing strong coupling methods as well as exact methods where possible. We discuss applications of this formalism to out of equilibrium situations.

arXiv:1801.01776•Nuclear Physics B•Mar 1, 2018•10 citations

Out-of-Equilibrium Dynamics of the Bose-Hubbard Model

Malcolm P. Kennett

The Bose-Hubbard model is the simplest model of interacting bosons on a lattice. It has recently been the focus of much attention due to the realization of this model with cold atoms in an optical lattice. The ability to tune parameters in the Hamiltonian as a function of time in cold atom systems has opened up the possibility of studying out-of-equilibrium dynamics, including crossing the quantum critical region of the model in a controlled way. In this paper, I give a brief introduction to the Bose Hubbard model, and its experimental realization and then give an account of theoretical and experimental efforts to understand out-of-equilibrium dynamics in this model, focusing on quantum quenches, both instantaneous and of finite duration. I discuss slow dynamics that have been observed theoretically and experimentally for some quenches from the superfluid phase to the Mott insulating phase and the picture of two timescales, one for fast local equilibration and another for slow global equilibration, that appears to characterize this situation. I also discuss the theoretical and experimental observation of the Lieb-Robinson bounds for a variety of quenches and the Kibble-Zurek mechanism in quenches from the Mott insulator to superfluid. I conclude with a discussion of open questions and future directions.

ISRN Condensed Matter Physics•Jun 12, 2013•22 citations

Phase and group velocities for correlation spreading in the Mott phase of the Bose-Hubbard model in dimensions greater than one

Ali Mokhtari-Jazi, Matthew Fitzpatrick, Malcolm P. Kennett

Lieb-Robinson and related bounds set an upper limit on the speed at which information propagates in nonrelativistic quantum systems. Experimentally, light-cone-like spreading has been observed for correlations in the Bose-Hubbard model (BHM) after a quantum quench. Using a two-particle irreducible (2PI) strong-coupling approach to out-of-equilibrium dynamics in the BHM we calculate both the group and phase velocities for the spreading of single-particle correlations in one, two, and three dimensions as a function of interaction strength. Our results are in quantitative agreement with measurements of the speed of spreading of single-particle correlations in both the one- and two-dimensional BHM realized with ultracold atoms. They also are consistent with the claim that the phase velocity rather than the group velocity was observed in recent experiments in two dimensions. We demonstrate that there can be large differences between the phase and group velocities for the spreading of correlations and explore how the anisotropy in the velocity varies across the phase diagram of the BHM. Our results establish the 2PI strong-coupling approach as a powerful tool to study out-of-equilibrium dynamics in the BHM in dimensions greater than one.

arXiv:2008.00015•Physical Review A•Feb 26, 2021•12 citations

Equilibrium and Non-equilibrium Properties of Strongly-Interacting Bosons in Optical Lattices

Wei Xu

This dissertation serves as a summary of my Ph.D. work numerically studying equilibrium and non-equilibrium properties of strongly-interacting one-dimensional (1D) boson systems. This work is motivated by the fact that 1D systems are realizable and highly controllable with ultracold atoms in optical lattice and atom chip experiments. We apply a recent worm algorithmic Monte Carlo approach developed for 1D continuous systems to study their equilibrium properties, both with and without an underlying lattice. We also apply an exact lattice approach based on the Bose-Fermi mapping to check our Monte Carlo results in the Tonks-Girardeau limit, and more importantly, to study far-from-equilibrium expansion dynamics of the systems. We first study the scaling of one-particle correlations of the harmonically trapped Lieb-Liniger gas with changing temperature and interaction strength. Based on the universal behaviors of the density and momentum profiles, we are able to determine the effective parameters needed to fully characterize the system. We also find that the Tonks-Girardeau limit at low temperatures is the ideal regime for the experimental observation of the $1/k^4$ momentum tail. An extra periodic lattice can drive the transition from superfluid to Mott insulator states. Exact and complete phase diagrams for such transitions are available only in the weak interacting and deep lattice limit, in which the system can be described using one-band Bose-Hubbard model. Beyond this limit, we use the worm algorithm in continuous space to map out the phase diagrams at various interaction strengths. We compare our phase diagrams with one-band Bose-Hubbard predictions and identify the regime where the one-band description breaks down. We introduce an inverse confined scattering solution to obtain effective Hubbard parameters, with which the Bose-Hubbard model provides correct results for strong interactions and deep lattices at unit filling. In addition to the equilibrium properties, we also study the expansion dynamics of ultracold atoms in the hard-core limit. Experimentally, this is usually achieved by turning off confining potentials and letting atoms expand in optical lattices. Theoretical studies from initial ground states predicted the occurrence of fermionization of the momentum distribution after long expansion times. In addition, quasicondensation at finite momenta emerges when expanding from Mott insulating domains. Here, we develop a finite-temperature extension of the lattice approach for dynamics. We find the dynamical ferminoization of the momentum distributions at all temperatures. For expansion from initial Mott domains, we observe enhanced correlations reminiscent of dynamical quasicondensation. Surprisingly, we find the systems appear to cool down during the melting of the Mott domains. We use an emergent local Hamiltonian to understand these emergent phenomena.

•Jan 1, 2018•0 citations

Non-equilibrium dynamics in the fully connected Bose-Hubbard model: the emergence of dynamical phase transitions and synchronization

Andrea Pizzi

In the context of out of equilibrium dynamics of isolated interacting quantum many-body systems, we study the mean-field discretized and nonlinear Gross-Pitaevskii equation (GPE) of motion for a fully connected (or all-to-all coupled) Bose-Hubbard model describing a large population of bosons on a lattice of generic dimension V, with potential experimental applications ranging from ultra-cold atomic gases in optical traps to systems of Josephson junctions. Considering a quench on the system Hamiltonian, we describe the time evolution of a system initialized to the Mott insulating phase and taken in the strong superfluid regime, finding a sharp change in the dynamical behaviour when varying the quench intensity across a critical value, that is a dynamical phase transition. Our results generalize the study presented in some previous papers focusing on the particular cases of V = 2, V = 3 and 1D lattice, and place particular emphasis on the relaxation of some dynamical order parameters in the large V limit. On the side of the dynamical phase transition corresponding to the strong superfluid regime, we find such relaxation to consist of a π-synchronization of the mean field bosonic phases, that we find to compete against the eventual site-dependent disorder of the Hamiltonian in what, in analogy with the notorious Kuramoto model for nonlinear coupled oscillators, can be referred to as a synchronization transition. Additionally, in support of the validity of the mean field approximation far from the transition, we report a study beyond mean field based on the Bogolioubov-de Gennes method.

•Jul 26, 2018•0 citations
Experimental advances have made ultracold atoms in optical lattices a favourable setting to study out-of-equilibrium phenomena and attracted considerable attention in recent years.These systems are highly versatile in that experimental parameters can be tuned over a wide range of values in real time.This facilitates the study of quantum quenches, in which parameters in the corresponding Hamiltonian are varied in time faster than the system can respond adiabatically.Such protocols open the door to a rich range of many-body physics and have been studied intensely both theoretically and experimentally.The Bose-Hubbard model (BHM) has been shown to describe interacting ultracold bosons in an optical lattice, allowing the opportunity for experiments to probe the out-of-equilibrium dynamics of the model.The BHM is a particularly convenient context for studying quantum quenches as it displays a quantum phase transition between superfluid and Mott-insulator phases.In this thesis, we develop a strong-coupling approach that allows the study of correlations in space and time in both the superfluid and Mott-insulating phases of the BHM.Specifically, we obtain a two-particle irreducible (2PI) effective action within the contourtime formalism that allows for a description of both equilibrium and out-of-equilibrium phenomena.We derive equations of motion for both the superfluid order parameter and the single-particle many-body Green's functions.First, we assess the accuracy of our formalism by studying the equilibrium solution for the homogeneous BHM and comparing our results to existing strong-coupling methods as well as to exact methods where possible.We then consider homogeneous systems that are initially thermalized in the Mott phase, and which are then subjected to quenches.We solve numerically the equations of motion for this scenario and calculate the single-particle density matrix.We demonstrate a Lieb-Robinson-like maximal propagation velocity for the spreading of single-particle correlations in one, two, and three dimensions.We study the dependence of the maximal propagation velocity on the quench protocol, chemical potential, temperature, and dimensionality.We compare our results to exact methods, existing strong-coupling approaches, and experiments where possible.Lastly, we extend our strong-coupling approach to the disordered BHM and derive equations of motion for the disorder-averaged single-particle Green's function.We discuss how these equations of motion can be used to study the phase stability of many-body localization in the disordered BHM for dimensions higher than one.
There has been considerable interest in the disordered Bose Hubbard model (BHM) in recent years, particularly in the context of thermalization and many-body localization. We develop a two-particle irreducible (2PI) strong-coupling approach to the disordered BHM that allows us to treat both equilibrium and out-of-equilibrium situations. We obtain equations of motion for spatio-temporal correlations and explore their equilibrium solutions. We study the equilibrium phase diagram as a function of disorder strength and discuss applications of the formalism to out-of-equilibrium situations. We also note that the disorder strengths where the emergence of non-ergodic dynamics was observed in a recent experiment [Choi et al. Science 352, 1547 (2016)] appear to correspond to the Mott insulator – Bose Glass phase boundary.
We develop a formalism that allows the study of correlations in space and time in both the superfluid and Mott insulating phases of the Bose–Hubbard Model. Specifically, we obtain a two particle irreducible effective action within the contour-time formalism that allows for both equilibrium and out of equilibrium phenomena. We derive equations of motion for both the superfluid order parameter and two-point correlation functions. To assess the accuracy of this formalism, we study the equilibrium solution of the equations of motion and compare our results to existing strong coupling methods as well as exact methods where possible. We discuss applications of this formalism to out of equilibrium situations.
The Bose-Hubbard model is the simplest model of interacting bosons on a lattice. It has recently been the focus of much attention due to the realization of this model with cold atoms in an optical lattice. The ability to tune parameters in the Hamiltonian as a function of time in cold atom systems has opened up the possibility of studying out-of-equilibrium dynamics, including crossing the quantum critical region of the model in a controlled way. In this paper, I give a brief introduction to the Bose Hubbard model, and its experimental realization and then give an account of theoretical and experimental efforts to understand out-of-equilibrium dynamics in this model, focusing on quantum quenches, both instantaneous and of finite duration. I discuss slow dynamics that have been observed theoretically and experimentally for some quenches from the superfluid phase to the Mott insulating phase and the picture of two timescales, one for fast local equilibration and another for slow global equilibration, that appears to characterize this situation. I also discuss the theoretical and experimental observation of the Lieb-Robinson bounds for a variety of quenches and the Kibble-Zurek mechanism in quenches from the Mott insulator to superfluid. I conclude with a discussion of open questions and future directions.
Lieb-Robinson and related bounds set an upper limit on the speed at which information propagates in nonrelativistic quantum systems. Experimentally, light-cone-like spreading has been observed for correlations in the Bose-Hubbard model (BHM) after a quantum quench. Using a two-particle irreducible (2PI) strong-coupling approach to out-of-equilibrium dynamics in the BHM we calculate both the group and phase velocities for the spreading of single-particle correlations in one, two, and three dimensions as a function of interaction strength. Our results are in quantitative agreement with measurements of the speed of spreading of single-particle correlations in both the one- and two-dimensional BHM realized with ultracold atoms. They also are consistent with the claim that the phase velocity rather than the group velocity was observed in recent experiments in two dimensions. We demonstrate that there can be large differences between the phase and group velocities for the spreading of correlations and explore how the anisotropy in the velocity varies across the phase diagram of the BHM. Our results establish the 2PI strong-coupling approach as a powerful tool to study out-of-equilibrium dynamics in the BHM in dimensions greater than one.
This dissertation serves as a summary of my Ph.D. work numerically studying equilibrium and non-equilibrium properties of strongly-interacting one-dimensional (1D) boson systems. This work is motivated by the fact that 1D systems are realizable and highly controllable with ultracold atoms in optical lattice and atom chip experiments. We apply a recent worm algorithmic Monte Carlo approach developed for 1D continuous systems to study their equilibrium properties, both with and without an underlying lattice. We also apply an exact lattice approach based on the Bose-Fermi mapping to check our Monte Carlo results in the Tonks-Girardeau limit, and more importantly, to study far-from-equilibrium expansion dynamics of the systems. We first study the scaling of one-particle correlations of the harmonically trapped Lieb-Liniger gas with changing temperature and interaction strength. Based on the universal behaviors of the density and momentum profiles, we are able to determine the effective parameters needed to fully characterize the system. We also find that the Tonks-Girardeau limit at low temperatures is the ideal regime for the experimental observation of the 1/k41/k^41/k4 momentum tail. An extra periodic lattice can drive the transition from superfluid to Mott insulator states. Exact and complete phase diagrams for such transitions are available only in the weak interacting and deep lattice limit, in which the system can be described using one-band Bose-Hubbard model. Beyond this limit, we use the worm algorithm in continuous space to map out the phase diagrams at various interaction strengths. We compare our phase diagrams with one-band Bose-Hubbard predictions and identify the regime where the one-band description breaks down. We introduce an inverse confined scattering solution to obtain effective Hubbard parameters, with which the Bose-Hubbard model provides correct results for strong interactions and deep lattices at unit filling. In addition to the equilibrium properties, we also study the expansion dynamics of ultracold atoms in the hard-core limit. Experimentally, this is usually achieved by turning off confining potentials and letting atoms expand in optical lattices. Theoretical studies from initial ground states predicted the occurrence of fermionization of the momentum distribution after long expansion times. In addition, quasicondensation at finite momenta emerges when expanding from Mott insulating domains. Here, we develop a finite-temperature extension of the lattice approach for dynamics. We find the dynamical ferminoization of the momentum distributions at all temperatures. For expansion from initial Mott domains, we observe enhanced correlations reminiscent of dynamical quasicondensation. Surprisingly, we find the systems appear to cool down during the melting of the Mott domains. We use an emergent local Hamiltonian to understand these emergent phenomena.
In the context of out of equilibrium dynamics of isolated interacting quantum many-body systems, we study the mean-field discretized and nonlinear Gross-Pitaevskii equation (GPE) of motion for a fully connected (or all-to-all coupled) Bose-Hubbard model describing a large population of bosons on a lattice of generic dimension V, with potential experimental applications ranging from ultra-cold atomic gases in optical traps to systems of Josephson junctions. Considering a quench on the system Hamiltonian, we describe the time evolution of a system initialized to the Mott insulating phase and taken in the strong superfluid regime, finding a sharp change in the dynamical behaviour when varying the quench intensity across a critical value, that is a dynamical phase transition. Our results generalize the study presented in some previous papers focusing on the particular cases of V = 2, V = 3 and 1D lattice, and place particular emphasis on the relaxation of some dynamical order parameters in the large V limit. On the side of the dynamical phase transition corresponding to the strong superfluid regime, we find such relaxation to consist of a π-synchronization of the mean field bosonic phases, that we find to compete against the eventual site-dependent disorder of the Hamiltonian in what, in analogy with the notorious Kuramoto model for nonlinear coupled oscillators, can be referred to as a synchronization transition. Additionally, in support of the validity of the mean field approximation far from the transition, we report a study beyond mean field based on the Bogolioubov-de Gennes method.