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.
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 present a Gutzwiller von Neumann dynamics (GvND) method for simulating equilibrium and nonequilibrium phenomena in strongly correlated electron systems. Our approach is a real-space formulation of the time-dependent Gutzwiller approximation method. Applying the GvND method to simulate interaction quenches in the Peierls-Hubbard model, we demonstrate the amplification of initial inhomogeneities, which in turn results in the collapse of quench-induced synchronized oscillation. Moreover, we find a dynamical phase transition separating two quasi-stationary regimes with rather distinct spatial distributions of physical quantities after the collapsed oscillation. In particular, in the strong-coupling regime, the system exhibits a dynamic phase separation in the quasi-stationary state. Our results thus underscore the importance of spatial fluctuations in the nonequilibrium dynamics of strongly correlated systems.