Efficient implementation of the Gutzwiller variational method

TL;DR

This paper introduces an efficient, self-consistent numerical method for solving the Gutzwiller variational problem in multi-band models, improving computational stability and symmetry handling, suitable for first-principles material studies.

Contribution

It generalizes previous approaches to include arbitrary on-site interactions without increasing computational complexity, simplifying the high-dimensional minimization process.

Findings

Method is highly efficient and stable

Able to incorporate full rotationally invariant Hund's interactions

Suitable for first-principles calculations in complex materials

Abstract

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 minimization by mapping it to a minimization only in the variational density matrix, in the spirit of the Levy and Lieb formulation of DFT. For fixed density the Gutzwiller renormalization matrix is determined as a fixpoint of a proper functional, whose evaluation only requires ground-state calculations of matrices defined in the Gutzwiller variational space. Furthermore, the proposed method is able to account for the symmetries of the variational function in a controlled way, reducing the number of variational parameters. After a detailed description of the method we present calculations for multi-band Hubbard models with full (rotationally invariant) Hund's rule on-site interaction. Our analysis shows that the numerical algorithm is very efficient, stable and easy to implement. For these reasons this method is particularly suitable for first principle studies -- e.g., in combination with DFT -- of many complex real materials, where the full intra-atomic interaction is important to obtain correct results.