The Generalized Gutzwiller Method for n=>2 Correlated Orbitals: Itinerant Ferromagnetism in eg-bands

We have generalized the Gutzwiller method to the cases of n=>2 correlated bands and report studies on a degenerate two-band model with Hund's rule type on-site interactions. At half band filling the metal-insulator transitions are usually of first order.

We develop a diagrammatic method for the evaluation of general multi-band Gutzwiller wave functions in finite dimensions. Our approach provides a systematic improvement of the widely used Gutzwiller approximation. As a first application we investigate itinerant ferromagnetism and correlation-induced deformations of the Fermi surface for a two-band Hubbard model on a square lattice.

We use the Gutzwiller variational many-body theory to investigate the stability of orbitally ordered states in a two-band Hubbard-model without spin degrees of freedom. Our results differ significantly from earlier Hartree-Fock calculations for this model. The Hartree-Fock phase diagram displays a large variety of orbital orders. In contrast, in the Gutzwiller approach orbital order only appears for densities in a narrow region around half filling.

In this tutorial presentation, we give a comprehensive introduction into the Gutzwiller variational approach and its merger with the density functional theory. The merits of this method are illustrated by a discussion of results for two-band Hubbard models and realistic multi-band Hamiltonians for iron-pnictide systems.

We introduce Gutzwiller conjugate gradient minimization (GCGM) theory, an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. GCGM uses the Gutzwiller wave function but does not use the commonly adopted Gutzwiller approximation (GA), which is a major source of inaccuracy. Instead, the theory uses an approximation that is based on the occupation probability of the on-site configurations, rather than approximations that decouple the...

We give a comprehensive introduction into an efficient numerical scheme for the minimisation of Gutzwiller energy functionals in studies on multi-band Hubbard models. Our method covers all conceivable cases of Gutzwiller variational wave functions and has been used successfully in previous numerical studies.

We present a self-consistent numerical approach to solve the Gutzwiller variational problem for general multi-band models with arbitrary on-site interaction. The proposed method generalizes and improves the procedure derived by Deng et al., Phys. Rev. B. 79 075114 (2009), overcoming the restriction to density-density interaction without increasing the complexity of the computational algorithm. Our approach drastically reduces the problem of the high-dimensional Gutzwiller min...

The long-standing problem of the effect of correlations on the ferromagnetism of is apparently nearing solution. The ferromagnetism of transition metals compounds, for instance doped manganites, poses a new question: is there some kind of orbital order coexisting with itinerant ferromagnetism? The ideas and techniques introduced by Gutzwiller should be of use again.

In the present paper, we propose an efficient numerical scheme for Gutzwiller method for multi-band Hubbard models with general onsite Coulomb interaction. Following the basic idea of Deng et al. [Phys. Rev. B 79, 075114 (2009)] and extensions by Lanata et al. [Phys. Rev. B 85, 035133 (2012)], the ground state is variationally determined through optimizing the total energy with respect to the variational single particle density matrix (n0), which is called "outer loop". In th...

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 LDA+Gu...