We study the unitary time evolution of antiferromagnetic order in the Hubbard model after a quench starting from the perfect Neel state. In this setup, which is well suited for experiments with cold atoms, one can distinguish fundamentally different
pathways for melting of long-range order at weak and strong interaction. In the Mott insulating regime, melting of long-range order occurs due to the ultra-fast transfer of energy from charge excitations to the spin background, while local magnetic moments and their exchange coupling persist during the process. The latter can be demonstrated by a local spin-precession experiment. At weak interaction, local moments decay along with the long-range order. The dynamics is governed by residual quasiparticles, which are reflected in oscillations of the off-diagonal components of the momentum distribution. Such oscillations provide an alternative route to study the prethermalization phenomenon and its influence on the dynamics away from the integrable (noninteracting) limit. The Hubbard model is solved within nonequilibrium dynamical mean-field theory, using the density matrix-renormalization group as an impurity solver.
Nonequilibrium dynamical mean-field theory (DMFT) solves correlated lattice models by obtaining their local correlation functions from an effective model consisting of a single impurity in a self-consistently determined bath. The recently developed m
apping of this impurity problem from the Keldysh time contour onto a time-dependent single-impurity Anderson model (SIAM) [C. Gramsch et al., Phys. Rev. B 88, 235106 (2013)] allows one to use wave function-based methods in the context of nonequilibrium DMFT. Within this mapping, long times in the DMFT simulation become accessible by an increasing number of bath orbitals, which requires efficient representations of the time-dependent SIAM wave function. These can be achieved by the multiconfiguration time-dependent Hartree (MCTDH) method and its multi-layer extensions. We find that MCTDH outperforms exact diagonalization for large baths in which the latter approach is still within reach and allows for the calculation of SIAMs beyond the system size accessible by exact diagonalization. Moreover, we illustrate the computation of the self-consistent two-time impurity Greens function within the MCTDH second quantization representation.
The nonequilibrium Dyson (or Kadanoff-Baym) equation, which is an equation of motion with long-range memory kernel for real-time Green functions, underlies many numerical approaches based on the Keldysh formalism. In this paper we map the problem of
solving the Dyson equation in real-time onto a noninteracting auxiliary Hamiltonian with additional bath degrees of freedom. The solution of the auxiliary model does not require the evaluation of a memory kernel and can thus be implemented in a very memory efficient way. The mapping is derived for a self-energy which is local in space and is thus directly applicable within nonequilibrium dynamical mean-field theory (DMFT). We apply the method to study the interaction quench in the Hubbard model for an optical lattice with a narrow confinement, using inhomogeneous DMFT in combination with second-order weak-coupling perturbation theory. We find that, although the quench excites pronounced density oscillations, signatures of the two-stage relaxation similar to the homogeneous system can be observed by looking at the time-dependent occupations of natural orbitals.
We derive an exact mapping from the action of nonequilibrium dynamical mean-field theory (DMFT) to a single-impurity Anderson model (SIAM) with time-dependent parameters, which can be solved numerically by exact diagonalization. The representability
of the nonequilibrium DMFT action by a SIAM is established as a rather general property of nonequilibrium Green functions. We also obtain the nonequilibrium DMFT equations using the cavity method alone. We show how to numerically obtain the SIAM parameters using Cholesky or eigenvector matrix decompositions. As an application, we use a Krylov-based time propagation method to investigate the Hubbard model in which the hopping is switched on, starting from the atomic limit. Possible future developments are discussed.
Extending our previous work, we explore the breathing mode---the [uniform] radial expansion and contraction of a spatially confined system. We study the breathing mode across the transition from the ideal quantum to the classical regime and confirm t
hat it is not independent of the pair interaction strength (coupling parameter). We present the results of time-dependent Hartree-Fock simulations for 2 to 20 fermions with Coulomb interaction and show how the quantum breathing mode depends on the particle number. We validate the accuracy of our results, comparing them to exact Configuration Interaction results for up to 8 particles.
