arXiv:1607.03016 · Physical review. B./Physical review. B

Self-consistent Keldysh approach to quenches in the weakly interacting Bose-Hubbard model

Nicola Lo Gullo, Luca Dell’Anna

Published Nov 18, 2016• 8 citations

Abstract

We present a nonequilibrium Green's-functional approach to study the dynamics following a quench in weakly interacting Bose-Hubbard model (BHM). The technique is based on the self-consistent solution of a set of equations which represents a particular case of the most general set of Hedin's equations for the interacting single-particle Green's function. We use the ladder approximation as a skeleton diagram for the two-particle scattering amplitude useful, through the self-energy in the Dyson equation, for finding the interacting single-particle Green's function. This scheme is then implemented numerically by a parallelized code. We exploit this approach to study the correlation propagation after a quench in the interaction parameter, for one and two dimensions. In particular, we show how our approach is able to recover the crossover from the ballistic to the diffusive regime by increasing the boson-boson interaction. Finally we also discuss the role of a thermal initial state on the dynamics both for one- and two-dimensional BHMs, finding that, surprisingly, at high temperature a ballistic evolution is restored.

DOI: 10.1103/physrevb.94.184308

Related papers

Quench Dynamics and Nonequilibrium Phase Diagram of the Bose-Hubbard Model

Corinna Kollath, Andreas Martin Läuchli, Ehud Altman

We investigate the time evolution of correlations in the Bose-Hubbard model following a quench from the superfluid to the Mott insulator. For large values of the final interaction strength the system approaches a distinctly nonequilibrium steady state that bears strong memory of the initial conditions. In contrast, when the final interaction strength is comparable to the hopping, the correlations are rather well approximated by those at thermal equilibrium. The existence of two distinct nonequilibrium regimes is surprising given the nonintegrability of the Bose-Hubbard model. We relate this phenomenon to the role of quasiparticle interactions in the Mott insulator.

arXiv:cond-mat/0607235Physical Review LettersApr 30, 2007569 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, 20191 citation

Light-cone-like spreading of single-particle correlations in the Bose-Hubbard model after a quantum quench in the strong-coupling regime

Matthew Fitzpatrick, Malcolm P. Kennett

We study the spreading of correlations in space and time after a quantum quench in the Bose-Hubbard model. We derive equations of motion for the single-particle Green's function within the contour-time formalism, allowing us to study dynamics in the strong-coupling regime. We discuss the numerical solutions of these equations and calculate the single-particle density matrix for quenches in the Mott phase. We demonstrate light-cone-like spreading of correlations in the Mott phase in one, two, and three dimensions and calculate the propagation velocities in each dimension.

arXiv:1606.04117Physical Review ANov 19, 201820 citations

Quenching the Hubbard Model: Comparison of Nonequilibrium Green's Function Methods

Jan‐Philip Joost, Christopher Makait, Alexander Osterkorn +1 more

ABSTRACT We benchmark nonequilibrium Green's function (NEGF) approaches for interaction quenches in the half‐filled Fermi–Hubbard model in one and two dimensions. We compare fully self‐consistent two‐time Kadanoff–Baym equations (KBE), the generalized Kadanoff–Baym ansatz (GKBA), and the recently developed NEGF‐based quantum fluctuations approach (NEGF‐QF) with second‐order Born (2B) and GW self‐energies. As references, we use DMRG in 1D and time‐dependent perturbation theory at weak coupling. While all methods reproduce the early‐time build‐up of correlations (prethermalization), different deficiencies emerge during the slower relaxation in the long‐time limit: KBE shows overdamped oscillations, GKBA exhibits a negative overshoot of the Fermi‐surface jump, and NEGF‐QF develops a prethermalization plateau with delayed relaxation. In 1D, we find that GKBA + 2B agrees best with DMRG across the explored parameter range, while QF + GW deviates the most, consistent with its reduced self‐consistency. In 2D, all approaches correctly capture the momentum‐dependence of the relaxation dynamics.

Contributions to Plasma PhysicsJan 15, 20261 citation

Lehmann representation of the nonequilibrium self-energy

Christian Gramsch, Michael Potthoff

It is shown that the nonequilibrium self-energy of an interacting lattice-fermion model has a unique Lehmann representation. Based on the construction of a suitable noninteracting effective medium, we provide an explicit and numerically practicable scheme to construct the Lehmann representation for the self-energy, given the Lehmann representation of the single-particle nonequilibrium Green's function. This is of particular importance for an efficient numerical solution of Dyson's equation in the context of approximations where the self-energy is obtained from a reference system with a small Hilbert space. As compared to conventional techniques to solve Dyson's equation on the Keldysh contour, the effective-medium approach allows us to reach a maximum propagation time, which can be several orders of magnitude longer. This is demonstrated explicitly by choosing the nonequilibrium cluster-perturbation theory as a simple approach to study the long-time dynamics of an inhomogeneous initial state after a quantum quench in the Hubbard model on a $10\ifmmode\times\else\texttimes\fi{}10$ square lattice. We demonstrate that the violation of conservation laws is moderate for weak Hubbard interaction and that the cluster approach is able to describe prethermalization physics.

arXiv:1509.05313Physical Review BDec 18, 201514 citations

Space-time correlations in the Bose Hubbard model after a quantum quench

Matthew Fitzpatrick, Malcolm P. Kennett

We study correlations in space and time after a quantum quench in the Bose Hubbard model. We obtain a two particle irreducible effective action within the Schwinger-Keldysh formalism that allows us to study dynamics in both the superfluid and Mott insulating phases. We derive equations of motion for both the superfluid order parameter and two-point correlation functions. We discuss the numerical solutions of these equations and calculate the quasimomentum distribution for quenches in the Mott phase. We demonstrate light-cone like spreading of correlations in the Mott phase in both one and two dimensions and calculate the propagation velocity in each dimension.

Jun 13, 20161 citation