Search

The self-energy functional theory (SFT) is generalized to describe the real-time dynamics of correlated lattice-fermion models far from thermal equilibrium. This is achieved by starting from a reformulation of the original equilibrium theory in terms of double-time Green's functions on the Keldysh-Matsubara contour. With the help of a generalized Luttinger-Ward functional, we construct a functional Ω[Σ] which is stationary at the physical (nonequilibrium) self-energy Σ and which yields the grand potential of the initial thermal state Ω at the physical point. Nonperturbative approximations can be defined by specifying a reference system that serves to generate trial self-energies. These self-energies are varied by varying the reference system's one-particle parameters on the Keldysh-Matsubara contour. In the case of thermal equilibrium, this approach reduces to the conventional SFT. Contrary to the equilibrium theory, however, "unphysical" variations, i.e., variations that are different on the upper and the lower branches of the Keldysh contour, must be considered to fix the time dependence of the optimal physical parameters via the variational principle. Functional derivatives in the nonequilibrium SFT Euler equation are carried out analytically to derive conditional equations for the variational parameters that are accessible to a numerical evaluation via a time-propagation scheme. Approximations constructed by means of the nonequilibrium SFT are shown to be inherently causal, internally consistent, and to respect macroscopic conservation laws resulting from gauge symmetries of the Hamiltonian. This comprises the nonequilibrium dynamical mean-field theory but also dynamical-impurity and variational-cluster approximations that are specified by reference systems with a finite number of degrees of freedom. In this way, nonperturbative and consistent approximations can be set up, the numerical evaluation of which is accessible to an exact-diagonalization approach. [ABSTRACT FROM AUTHOR]

We extend the nonequilibrium dynamical mean field (DMFT) formalism to inhomogeneous systems by adapting the "real-space" DMFT method to Keldysh Green's functions. Solving the coupled impurity problems using strong-coupling perturbation theory, we apply the formalism to homogeneous and inhomogeneous layered systems with strong local interactions and up to 39 layers. We study the diffusion of doublons and holes created by photoexcitation in a Mott insulating system, the time-dependent build-up of the polarization and the current induced by a linear voltage bias across a multilayer structure, and the photoinduced current in a Mott insulator under bias. [ABSTRACT FROM AUTHOR]

We investigate the possibility to assist the numerically ill-posed calculation of spectral properties of interacting quantum systems in thermal equilibrium by extending the imaginary-time simulation to a finite Schwinger-Keldysh contour. The effect of this extension is tested within the standard maximum entropy approach to analytic continuation. We find that the inclusion of real-time data improves the resolution of structures at high energy, while the imaginary-time data are needed to correctly reproduce low-frequency features such as quasiparticle peaks. As a nonequilibrium application, we consider the calculation of time-dependent spectral functions from retarded Green function data on a finite time interval, and compare the maximum entropy approach to direct Fourier transformation and a method based on Pad6 approximants. [ABSTRACT FROM AUTHOR]

• Dynamical mean-field theory simulations of time-resolved photoemission spectroscopy. Time- and angular-resolved photoemission spectroscopy (trARPES) can directly probe the electronic structure of quantum materials out of equilibrium. This can shed light on the interaction of the electrons with spin, lattice, and orbital degrees of freedom, and help to unravel pathways towards novel out-of-equilibrium phases. Dynamical mean-field theory (DMFT) and its extensions provide a versatile toolbox to interpret such experiments through a theoretical simulation of the underlying microscopic processes. The approach can be applied both to Mott insulators and correlated metals, and it is formulated in terms of non-equilibrium Green's functions, which directly relate to the photoemission spectrum. This article reviews the theoretical description of trARPES within DMFT and related diagrammatic non-equilibrium Green's function techniques. Several applications are discussed, including the photo-induced melting of excitonic order, femtosecond relaxation processes in Mott insulators, and the manipulation of the electronic structure of Mott and charge transfer insulators using photo-doping and strong THz fields. [ABSTRACT FROM AUTHOR]

A fast time propagation method for nonequilibrium Green's functions (NEGF) based on the generalized Kadanoff–Baym Ansatz (GKBA) is applied to a lattice system with a symmetry‐broken equilibrium phase, namely an excitonic insulator (EI). The adiabatic preparation of a correlated symmetry‐broken initial state from a Hartree–Fock wave function within GKBA is assessed by comparing with a solution of the imaginary‐time Dyson equation. It is found that it is possible to reach a symmetry‐broken correlated initial state with nonzero excitonic order parameter by the adiabatic switching (AS) procedure. It is discussed under which circumstances this is possible in practice within reasonably short switching times. [ABSTRACT FROM AUTHOR]

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 [ABSTRACT FROM AUTHOR]

We use a nonequilibrium implementation of the dynamical cluster approximation (DCA) to study the effect of short-range correlations on the dynamics of the two-dimensional Falicov-Kimball model after an interaction quench. As in the case of single-site dynamical mean-field theory, thermalization is absent in DCA simulations, and for quenches across the metal-insulator boundary, nearest-neighbor charge correlations in the nonthermal steady state are found to be larger than in the thermal state with identical energy. We investigate to what extent it is possible to define an effective temperature of the trapped state after a quench. Based on the ratio between the lesser and retarded Green's function, we conclude that a roughly thermal distribution is reached within the energy intervals corresponding to the momentum-patch dependent subbands of the spectral function. The effectively different chemical potentials of these distributions, however, lead to a very hot, or even negative, effective temperature in the energy intervals between these subbands. [ABSTRACT FROM AUTHOR]

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 mapping 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 multilayer 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 Green s function within the MCTDH second quantization representation. [ABSTRACT FROM AUTHOR]