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

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

Andrea Pizzi

Published Jul 26, 2018• 0 citations
OpenAlex

Abstract

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.

Related papers

Quench-induced dynamical phase transitions and π -synchronization in the Bose-Hubbard model

Andrea Pizzi, Fabrizio Dolcini, Karyn Le Hur

We investigate the nonequilibrium behavior of a fully connected (or all-to-all coupled) Bose-Hubbard model after a Mott to superfluid quench, in the limit of large boson densities and for an arbitrary number $V$ of lattice sites, with potential relevance in experiments ranging from cold atoms to superconducting qubits. By means of the truncated Wigner approximation, we predict that crossing a critical quench strength the system undergoes a dynamical phase transition between two regimes that are characterized at long times either by an inhomogeneous population of the lattice (i.e., macroscopical self-trapping) or by the tendency of the mean-field bosonic variables to split into two groups with phase difference $\ensuremath{\pi}$, that we refer to as $\ensuremath{\pi}$-synchronization. We show the latter process to be intimately connected to the presence, only for $V\ensuremath{\ge}4$, of a manifold of infinitely many fixed points of the dynamical equations. Finally, we show that no fine-tuning of the model parameters is needed for the emergence of such $\ensuremath{\pi}$-synchronization, that is in fact found to vanish smoothly in presence of an increasing site-dependent disorder, in what we call a synchronization crossover.

arXiv:1810.12414•Physical review. B./Physical review. B•Mar 4, 2019•11 citations

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

Matthew Fitzpatrick

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.

•Jul 5, 2019•1 citation

Mean-Field Dynamics of the Bose–Hubbard Model in High Dimension

Shahnaz Farhat, Denis Périce, Sören Petrat

Abstract This paper is a first step toward rigorously establishing the validity of mean-field approximations for the dynamics of quantum systems in a large dimensional limit. Our result is a derivation of the large d mean-field limit of the Bose–Hubbard model on a square lattice. The Bose–Hubbard model effectively describes bosons on a lattice with on-site interactions and nearest-neighbour hopping, serving as a foundational framework for understanding strong particle interactions and the superfluid to Mott insulator phase transition. We prove a trace norm estimate between the one-lattice-site reduced density of the Schrödinger dynamics and the mean-field dynamics in the limit of large dimension. Here, the mean-field approximation is in the hopping amplitude and not in the interaction, leading to a very rich and non-trivial mean-field equation. This mean-field equation does not only describe the condensate, as is the case when the mean-field description comes from a large particle number limit averaging out the interaction, but it allows for a phase transition since it contains the full non-trivial interaction. Our work is a rigorous justification of a simple case of the highly successful dynamical mean-field theory (DMFT) for bosons, which somewhat surprisingly yields many qualitatively correct results in three dimensions.

arXiv:2501.05304•Annales Henri Poincaré•Oct 10, 2026•0 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

Phases and non-equilibrium dynamics of ultracold bosons in an optical quasicrystal

Reeve, Leanne

This thesis concerns the physics of ultracold bosons in quasicrystalline lattices, studied by means of quantum simulation using ultracold 39K atoms loaded to an eight-fold rotationally symmetric optical quasicrystal. Quasicrystalline and quasiperiodic systems, thanks to their long-range order and lack of periodicity, exhibit a wealth of interesting properties such as a dense and self-similar momentum space, and a localization transition to Anderson insulator or Bose glass states (in non-interacting or interacting cases respectively) reminiscent of that in disordered systems, whilst lacking the rare regions that typically exist in such systems. Here we will cover experiments into this system’s behaviour in both equilibrium and non-equilibrium settings. First, we focus on mapping out its phase diagram in both weakly and strongly interacting regimes. We observe consistency with numerical predictions for the superfluid-Bose glass phase transition using coherence measurements, as well as finding a region of roughly zero compressibility around where the numerics predicts the existence of a Mott insulator ground state. We also study which regimes of the phase diagram can be ramped through adiabatically, finding that entropy is created upon traversing the Bose glass phase irrespective of ramp duration, which we attribute to its localized nature. Next we look into the dynamics of the system after a quench of the lattice depth in the weakly interacting regime. Upon quenching from the superfluid regime to the Bose glass regime, we observe a quantum walk to higher order diffraction peaks in momentum space whose maximum order is unbounded by kinetic energy thanks to the dense, self-similar nature of the quasicrystal’s momentum space. Quenching in the opposite direction we observe instead a walk between sites in real space for short times by looking at the position (as opposed to momentum) space autocorrelation function of the cloud, finding an expansion rate in agreement with an intuitive estimate based on the Lieb-Robinson bounds of a periodic system. We also look into the timescales of these dynamics in the context of what we might expect from the underlying Hamiltonian, before finally investigating the loss of coherence exhibited by the system after a slower quench of the lattice depth from an initially superfluid cloud. We observe that coherence is lost sooner after the quench for deeper final lattices, with the loss becoming almost instantaneous for quenches that cross the localization transition across the whole cloud. After this we explore the system’s transport dynamics by means of a quench to an excited state at the same time as an additional confining potential is switched off. We observe that unlike in the case of a periodic square lattice, where ballistic expansion that scales with the tunnelling strength occurs in the non-interacting case before becoming slower upon the introduction of interactions thanks to diffusion, in the quasicrystal the transport rate no longer scales with tunnelling strength. Instead we observe a drop relative to the tunnelling rate by 2 − 3 orders of magnitude as the lattice depth is increased from 1Er to 5Er with little dependence on the interaction strength. We also note that the finite rate measured for the deeper lattices is potentially an overestimate since the finite loss rate of atoms makes mass transport indistinguishable from density-dependent losses at that point. Finally, we discuss use of a quantum gas magnification scheme with which to magnify the real space density distribution of our atoms such that it can be measured with single-site resolution. Some preliminary results are presented regarding how the density distribution changes with lattice depth in the weakly interacting regime, in which we observe that as lattice depth increases a higher proportion of atoms occupy deeper sites close to the centre of the quasicrystal’s configuration space, as we might expect.

•Jun 23, 2025•0 citations

Bose-Hubbard model in the canonical ensemble: a beyond mean-field approach

Tista Banerjee

Ultracold atoms in optical lattices are versatile testbeds to study and manipulate equilibrium and out-of-equilibrium aspects of quantum many-body systems whose behavior can be described by Hubbard-type Hamiltonians. In this paper, we consider an ansatz wave-function which respects total particle-number conservation for such systems and goes beyond mean-field theory; this wave-function has the same complexity in the number of parameters as the mean-field Gutzwiller ansatz, and captures quantum correlations and entanglement via projection onto an effective low-energy manifold. This ansatz can be exploited to study quantum phases observed in a large class of systems realizable in such experimental platforms and is useful to study quantum dynamics. We show that the relaxation dynamics of various out-of-equilibrium initial states under sudden quench of Hamiltonian parameters can be studied with this ansatz wavefunction within the framework of time-dependent variational principle. We present a quantitative comparison with small-scale exact diagonalization results in the 1D Bose-Hubbard model with and without external trapping potentials.

arXiv:2508.01692•Aug 3, 2025•0 citations
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.
We investigate the nonequilibrium behavior of a fully connected (or all-to-all coupled) Bose-Hubbard model after a Mott to superfluid quench, in the limit of large boson densities and for an arbitrary number VVV of lattice sites, with potential relevance in experiments ranging from cold atoms to superconducting qubits. By means of the truncated Wigner approximation, we predict that crossing a critical quench strength the system undergoes a dynamical phase transition between two regimes that are characterized at long times either by an inhomogeneous population of the lattice (i.e., macroscopical self-trapping) or by the tendency of the mean-field bosonic variables to split into two groups with phase difference $\ensuremath{\pi}$, that we refer to as $\ensuremath{\pi}$-synchronization. We show the latter process to be intimately connected to the presence, only for $V\ensuremath{\ge}4$, of a manifold of infinitely many fixed points of the dynamical equations. Finally, we show that no fine-tuning of the model parameters is needed for the emergence of such $\ensuremath{\pi}$-synchronization, that is in fact found to vanish smoothly in presence of an increasing site-dependent disorder, in what we call a synchronization crossover.
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.
Abstract This paper is a first step toward rigorously establishing the validity of mean-field approximations for the dynamics of quantum systems in a large dimensional limit. Our result is a derivation of the large d mean-field limit of the Bose–Hubbard model on a square lattice. The Bose–Hubbard model effectively describes bosons on a lattice with on-site interactions and nearest-neighbour hopping, serving as a foundational framework for understanding strong particle interactions and the superfluid to Mott insulator phase transition. We prove a trace norm estimate between the one-lattice-site reduced density of the Schrödinger dynamics and the mean-field dynamics in the limit of large dimension. Here, the mean-field approximation is in the hopping amplitude and not in the interaction, leading to a very rich and non-trivial mean-field equation. This mean-field equation does not only describe the condensate, as is the case when the mean-field description comes from a large particle number limit averaging out the interaction, but it allows for a phase transition since it contains the full non-trivial interaction. Our work is a rigorous justification of a simple case of the highly successful dynamical mean-field theory (DMFT) for bosons, which somewhat surprisingly yields many qualitatively correct results in three dimensions.
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.
This thesis concerns the physics of ultracold bosons in quasicrystalline lattices, studied by means of quantum simulation using ultracold 39K atoms loaded to an eight-fold rotationally symmetric optical quasicrystal. Quasicrystalline and quasiperiodic systems, thanks to their long-range order and lack of periodicity, exhibit a wealth of interesting properties such as a dense and self-similar momentum space, and a localization transition to Anderson insulator or Bose glass states (in non-interacting or interacting cases respectively) reminiscent of that in disordered systems, whilst lacking the rare regions that typically exist in such systems. Here we will cover experiments into this system’s behaviour in both equilibrium and non-equilibrium settings. First, we focus on mapping out its phase diagram in both weakly and strongly interacting regimes. We observe consistency with numerical predictions for the superfluid-Bose glass phase transition using coherence measurements, as well as finding a region of roughly zero compressibility around where the numerics predicts the existence of a Mott insulator ground state. We also study which regimes of the phase diagram can be ramped through adiabatically, finding that entropy is created upon traversing the Bose glass phase irrespective of ramp duration, which we attribute to its localized nature. Next we look into the dynamics of the system after a quench of the lattice depth in the weakly interacting regime. Upon quenching from the superfluid regime to the Bose glass regime, we observe a quantum walk to higher order diffraction peaks in momentum space whose maximum order is unbounded by kinetic energy thanks to the dense, self-similar nature of the quasicrystal’s momentum space. Quenching in the opposite direction we observe instead a walk between sites in real space for short times by looking at the position (as opposed to momentum) space autocorrelation function of the cloud, finding an expansion rate in agreement with an intuitive estimate based on the Lieb-Robinson bounds of a periodic system. We also look into the timescales of these dynamics in the context of what we might expect from the underlying Hamiltonian, before finally investigating the loss of coherence exhibited by the system after a slower quench of the lattice depth from an initially superfluid cloud. We observe that coherence is lost sooner after the quench for deeper final lattices, with the loss becoming almost instantaneous for quenches that cross the localization transition across the whole cloud. After this we explore the system’s transport dynamics by means of a quench to an excited state at the same time as an additional confining potential is switched off. We observe that unlike in the case of a periodic square lattice, where ballistic expansion that scales with the tunnelling strength occurs in the non-interacting case before becoming slower upon the introduction of interactions thanks to diffusion, in the quasicrystal the transport rate no longer scales with tunnelling strength. Instead we observe a drop relative to the tunnelling rate by 2 − 3 orders of magnitude as the lattice depth is increased from 1Er to 5Er with little dependence on the interaction strength. We also note that the finite rate measured for the deeper lattices is potentially an overestimate since the finite loss rate of atoms makes mass transport indistinguishable from density-dependent losses at that point. Finally, we discuss use of a quantum gas magnification scheme with which to magnify the real space density distribution of our atoms such that it can be measured with single-site resolution. Some preliminary results are presented regarding how the density distribution changes with lattice depth in the weakly interacting regime, in which we observe that as lattice depth increases a higher proportion of atoms occupy deeper sites close to the centre of the quasicrystal’s configuration space, as we might expect.
Ultracold atoms in optical lattices are versatile testbeds to study and manipulate equilibrium and out-of-equilibrium aspects of quantum many-body systems whose behavior can be described by Hubbard-type Hamiltonians. In this paper, we consider an ansatz wave-function which respects total particle-number conservation for such systems and goes beyond mean-field theory; this wave-function has the same complexity in the number of parameters as the mean-field Gutzwiller ansatz, and captures quantum correlations and entanglement via projection onto an effective low-energy manifold. This ansatz can be exploited to study quantum phases observed in a large class of systems realizable in such experimental platforms and is useful to study quantum dynamics. We show that the relaxation dynamics of various out-of-equilibrium initial states under sudden quench of Hamiltonian parameters can be studied with this ansatz wavefunction within the framework of time-dependent variational principle. We present a quantitative comparison with small-scale exact diagonalization results in the 1D Bose-Hubbard model with and without external trapping potentials.