Within a Lagrangian formalism, we derive the time-dependent Gutzwiller approximation for general multi-band Hubbard models. Our approach explicitly incorporates the coupling between time-dependent variational parameters and a time-dependent density matrix from which we obtain dynamical correlation functions in the linear-response regime. Our results are illustrated for the one-band model where we show that the interacting system can be mapped to an effective problem of fermionic quasiparticles coupled to 'doublon' (double occupancy) bosonic fluctuations. The latter have an energy on the scale of the on-site Hubbard repulsion U in the dilute limit but become soft at the Brinkman–Rice transition, which is shown to be related to an emerging conservation law of doublon charge and the associated gauge invariance. Coupling with the boson mode produces a structure in the charge response and we find that a similar structure appears in dynamical mean-field theory.
Within a Lagrangian formalism, we derive the time-dependent Gutzwiller approximation for general multi-band Hubbard models. Our approach explicitly incorporates the coupling between time-dependent variational parameters and a time-dependent density matrix from which we obtain dynamical correlation functions in the linear-response regime. Our results are illustrated for the one-band model where we show that the interacting system can be mapped to an effective problem of fermionic quasiparticles coupled to 'doublon' (double occupancy) bosonic fluctuations. The latter have an energy on the scale of the on-site Hubbard repulsion U in the dilute limit but become soft at the Brinkman–Rice transition, which is shown to be related to an emerging conservation law of doublon charge and the associated gauge invariance. Coupling with the boson mode produces a structure in the charge response and we find that a similar structure appears in dynamical mean-field theory.
Within a Lagrangian formalism, we derive the time-dependent Gutzwiller approximation for general multi-band Hubbard models. Our approach explicitly incorporates the coupling between time-dependent variational parameters and a time-dependent density matrix from which we obtain dynamical correlation functions in the linear-response regime. Our results are illustrated for the one-band model where we show that the interacting system can be mapped to an effective problem of fermionic quasiparticles coupled to 'doublon' (double occupancy) bosonic fluctuations. The latter have an energy on the scale of the on-site Hubbard repulsion U in the dilute limit but become soft at the Brinkman–Rice transition, which is shown to be related to an emerging conservation law of doublon charge and the associated gauge invariance. Coupling with the boson mode produces a structure in the charge response and we find that a similar structure appears in dynamical mean-field theory.
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.
The description of electronic inhomogeneities in transition metal compounds requires a technique which is able to deal with the excitations of an inhomogeneous, strongly correlated system. In this talk I will outline the formalism of the timedependent Gutzwiller approximation [1] for the Hubbard model which allows for the computation of random phase approximation (RPA) like fluctuations on top of the Gutzwiller approximation (GA). The performance of the method is demonstrated by comparing static quantities like ground state energy and double occupancies as well as dynamical correlation functions in the charge, spin and pair channels with HartreeFock+RPA and exact results. I will further show that the method gives a good account of several spectroscopic quantities related to stripes in highTc superconductors. A major part of the talk will address the recent generalization of the method to the case of multiple bands [2,3]. Contrary to the standard RPA, the incorporation of the local multipletstructure within our approach leads to a much larger sensitivity of ferromagnetism on the the ratio between Hund coupling and intraorbital repulsion J/U. Besides, we find that the Brinkman Rice localization transition, which in the paramagnetic GA appears at integer doping is always masked by an antiferromagnetic or an orbitally ordered ground state. In addition, a quantum Lifshitz point is found at larger values of J/U which corresponds to the common boundary of ferro, para and incommensurate magnetic phases. [1] G. Seibold and J. Lorenzana, Phys. Rev. Lett. 86, 2605 (2001) [2] E. von Oelsen, G. Seibold, and J. Bunemann, New J. Phys. 13, 113031 (2011) [3] E. v. Oelsen, G. Seibold, and J. Bunemann, Phys. Rev. Lett. 107, 076402 (2011)
Based on the variational Gutzwiller theory, we present a method for the computation of response functions for multiband Hubbard models with general local Coulomb interactions. The improvement over the conventional random-phase approximation is exemplified for an infinite-dimensional two-band Hubbard model where the incorporation of the local multiplet structure leads to a much larger sensitivity of ferromagnetism on the Hund coupling. Our method can be implemented into local-density approximation and Gutzwiller schemes and will therefore be an important tool for the computation of response functions for strongly correlated materials.
We develop an extension of the Gutzwiller approximation (GA) formalism that includes the effects of Coulomb interactions of arbitrary range (including density density, exchange, pair hopping, and Coulomb-assisted hopping terms). This formalism reduces to the ordinary GA formalism for the multiband Hubbard models in the presence of only local interactions. This is accomplished by combining the 1/z expansion---where z is the coordination number, and only the leading-order terms contribute in the limit of infinite dimensions---with a ${P}_{R}^{\ifmmode\dagger\else\textdagger\fi{}}{P}_{R}\ensuremath{-}I$ expansion, where PR is the Gutzwiller projector on the site R. The method is conveniently formulated in terms of a Gutzwiller Lagrange function. We apply our theory to the extended single-band Hubbard model. Similarly to the usual Brinkman-Rice mechanism, we find a Mott transition. A valence skipping transition is observed, where the occupation of the empty and doubly occupied states for the Gutzwiller wave function is enhanced with respect to the uncorrelated Slater determinant wave function.
We present a simple scheme to evaluate linear response functions including quantum fluctuation corrections on top of the Gutzwiller approximation. The method is derived for a generic multiband lattice Hamiltonian without any assumption about the dynamics of the variational correlation parameters that define the Gutzwiller wavefunction, and which thus behave as genuine dynamical degrees of freedom that add on those of the variational uncorrelated Slater determinant. We apply the method to the standard half-filled single-band Hubbard model. We are able to recover known results, but, as by-product, we also obtain few novel ones. In particular, we show that quantum fluctuations can reproduce almost quantitatively the behaviour of the uniform magnetic susceptibility uncovered by dynamical mean field theory, which, though enhanced by correlations, is found to be smooth across the paramagnetic Mott transition. By contrast, the simple Gutzwiller approximation predicts that susceptibility to diverge at the transition.