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.
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.
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.
Motivated by recent experiments, we study the relaxation dynamics and thermalization in the one-dimensional Bose-Hubbard model induced by a global interaction quench. Specifically, we start from an initial state that has exactly one boson per site and is the ground state of a system with infinitely strong repulsive interactions at unit filling. Using exact diagonalization and the density-matrix renormalization-group method, we compute the time dependence of such observables as the multiple occupancy and the momentum distribution function. Typically, the relaxation to stationary values occurs over just a few tunneling times. The stationary values are identical to the so-called diagonal ensemble on the system sizes accessible to our numerical methods and we further observe that the microcanonical ensemble describes the time averages of many observables reasonably well for small and intermediate interaction strength. The expectation values of observables in the canonical ensemble agree quantitatively with the time averages obtained from the quench at small interaction strengths, and qualitatively provide a good description even in parameter regimes where the microcanonical ensemble is not applicable due to finite-size effects. We discuss our numerical results in the framework of the eigenstate thermalization hypothesis. Moreover, we also observe that the diagonal and the canonical ensembles are practically identical for our initial conditions already on the level of their respective energy distributions for small interaction strengths. Finally, we discuss implications of our results for the interpretation of a recent sudden expansion experiment [Phys. Rev. Lett. 110, 205301 (2013)], in which the same interaction quench was realized.
We study the nonequilibrium dynamics in the fermionic Hubbard model after a sudden change of the interaction strength. To this scope, we introduce a time-dependent variational approach in the spirit of the Gutzwiller ansatz. At the saddle-point approximation, we find at half filling a sharp transition between two different regimes of small and large coherent oscillations, separated by a critical line of quenches where the system is found to relax. Any finite doping washes out the transition, leaving aside just a sharp crossover. In order to investigate the role of quantum fluctuations, we map the model onto an auxiliary quantum Ising model in a transverse field coupled to free fermionic quasiparticles. Remarkably, the Gutzwiller approximation turns out to correspond to the mean-field decoupling of this model in the limit of infinite coordination lattices. The advantage is that we can go beyond mean field and include Gaussian fluctuations around the non-equilibrium mean-field dynamics. Unlike at equilibrium, we find that quantum fluctuations become massless and eventually unstable before the mean-field dynamical critical line, which suggests they could even alter qualitatively the mean-field scenario.
We review the recently proposed extension of the Gutzwiller approximation, M. Schiro' and M. Fabrizio, Phys. Rev. Lett. 105, 076401 (2010), designed to describe the out-of-equilibrium time-evolution of a Gutzwiller-type variational wave function for correlated electrons. The method, which is strictly variational in the limit of infinite lattice-coordination, is quite general and flexible, and it is applicable to generic non-equilibrium conditions, even far beyond the linear response regime. As an application, we discuss the quench dynamics of a single-band Hubbard model at half-filling, where the method predicts a dynamical phase transition above a critical quench that resembles the sharp crossover observed by time-dependent dynamical mean field theory. We next show that one can actually define in some cases a multi-configurational wave function combination of a whole set of mutually orthogonal Gutzwiller wave functions. The Hamiltonian projected in that subspace can be exactly evaluated and is equivalent to a model of auxiliary spins coupled to non-interacting electrons, closely related to the slave-spin theories for correlated electron models. The Gutzwiller approximation turns out to be nothing but the mean-field approximation applied to that spin-fermion model, which displays, for any number of bands and integer fillings, a spontaneous Z2 symmetry breaking that can be identified as the Mott insulator-to-metal transition.
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.
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.