© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2111.09166

RTGW2020: A powerful implementation of DFT + Gutzwiller method

Shiyu Peng, Hongming Weng, Xi Dai

Published Nov 17, 2021• 0 citations

Abstract

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 the corresponding "inner loop" where n0 is fixed, the non-interacting wave function and the parameters contained in the Gutzwiller projector are determined by a two-step iterative approach. All derivatives of the implementation process have been analytically derived, which allows us to apply some advanced minimization or root-searching algorithms for both the inner and outer loops leading to the highly efficient convergence. In addition, an atomic diagonalization method taking the point group symmetry into account has been developed for the customized design of the Gutzwiller projector, making it convenient to explore many interesting orders at a lower cost of computation. As benchmarks, several different types of correlated models have been studied utilizing the proposed method, which are in perfect agreement with the previous results by DMFT and multi-orbital slave-boson mean field method. Compared with the linear mixing method, the newton method with analytical derivatives shows much faster convergence for the inner loop. As for the outer loop, the minimization using analytical derivatives also shows much better stability and efficiency compared with that using numerical derivatives.

DOI: 10.48550/arxiv.2111.09166

Related papers

Efficient implementation of the Gutzwiller variational method

Nicola Lanatà, Hugo U. R. Strand, Xi Dai +1 more

We present a self-consistent numerical approach to solve the Gutzwiller variational problem for general multiband 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 density functional theory (DFT). For fixed density the Gutzwiller renormalization matrix is determined as a fixpoint of a proper functional, whose evaluation requires only 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 multiband 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-principles studies (e.g., in combination with DFT) of many complex real materials, where the full intra-atomic interaction is important to obtain correct results.

arXiv:1108.0180•Physical Review B•Jan 31, 2012•87 citations

Gutzwiller variational method for intersite Coulomb interactions: The spinless fermion model in one dimension

Mohammad Sherafati, S. Satpathy, Dix H. Pettey

We study the Gutzwiller method for the spinless fermion model in one dimension, which is one of the simplest models that incorporates the intersite Coulomb interaction. The Gutzwiller solution of this model has been studied in the literature but with differing results. We obtain the Gutzwiller solution of the problem by a careful enumeration of the many-particle configurations and explain the origin of the discrepancy in the literature to be due to the neglect of the correlation between the neighboring bond occupancies. The correct implementation of the Gutzwiller approach yields results different from the slave-boson mean-field theory, unlike for the Hubbard model with on-site interaction, where both methods are known to be equivalent. The slave-boson and the Gutzwiller results are compared to the exact solution, available for the half-filled case, and to the numerical exact diagonalization results for the case of general filling.

arXiv:1302.2108•Physical Review B•Jul 11, 2013•3 citations

Approximation schemes for the study of multi-band Gutzwiller wave functions

Jörg Bünemann, Thorben Linneweber, Florian Gebhard

The minimum of the Gutzwiller energy functional depends on the number of parameters considered in the variational state. For a three-orbital Hubbard model we find that the frequently used diagonal Ansatz is very accurate in high-symmetry situations. For lower symmetry, induced by a crystal-field splitting or the spin-orbit coupling, the discrepancies in energy between the most general and a diagonal Gutzwiller Ansatz can be quite significant. We discuss approximate schemes that may be employed in multi-band cases where a minimization of the general Gutzwiller energy functional is too demanding numerically.

•Jan 1, 2017•1 citation

A rotationally invariant approach based on Gutzwiller wave function for correlated electron systems

Zhuo Ye, Feng Zhang, Yimei Fang +6 more

Abstract We introduce a rotationally invariant approach combined with the Gutzwiller conjugate gradient minimization method to study correlated electron systems. In the approach, the Gutzwiller projector is parametrized based on the number of electrons occupying the onsite orbitals instead of the onsite configurations. The approach efficiently groups the onsite orbitals according to their symmetry and greatly reduces the computational complexity, which yields a speedup of 20 ∼ 50 × in the minimal basis energy calculation of dimers. The computationally efficient approach promotes more accurate calculations beyond the minimal basis that is inapplicable in the original approach. A large-basis energy calculation of F2 demonstrates favorable agreements with standard quantum-chemical calculations Bytautas et al (2007 J. Chem. Phys. 127 164317)

Journal of Physics Condensed Matter•Oct 11, 2022•1 citation

The Out-of-Equilibrium Time-Dependent Gutzwiller Approximation

Michele Fabrizio

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 $Z_2$ symmetry breaking that can be identified as the Mott insulator-to-metal transition.

•Aug 16, 2016•8 citations

Extending the Gutzwiller approximation to intersite interactions

Garry Goldstein, Nicola Lanatà, Gabriel Kotliar

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 ${P}_{R}$ 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.

arXiv:2001.01353•Physical review. B./Physical review. B•Jul 31, 2020•4 citations

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 the corresponding "inner loop" where n0 is fixed, the non-interacting wave function and the parameters contained in the Gutzwiller projector are determined by a two-step iterative approach. All derivatives of the implementation process have been analytically derived, which allows us to apply some advanced minimization or root-searching algorithms for both the inner and outer loops leading to the highly efficient convergence. In addition, an atomic diagonalization method taking the point group symmetry into account has been developed for the customized design of the Gutzwiller projector, making it convenient to explore many interesting orders at a lower cost of computation. As benchmarks, several different types of correlated models have been studied utilizing the proposed method, which are in perfect agreement with the previous results by DMFT and multi-orbital slave-boson mean field method. Compared with the linear mixing method, the newton method with analytical derivatives shows much faster convergence for the inner loop. As for the outer loop, the minimization using analytical derivatives also shows much better stability and efficiency compared with that using numerical derivatives.

We present a self-consistent numerical approach to solve the Gutzwiller variational problem for general multiband 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 density functional theory (DFT). For fixed density the Gutzwiller renormalization matrix is determined as a fixpoint of a proper functional, whose evaluation requires only 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 multiband 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-principles studies (e.g., in combination with DFT) of many complex real materials, where the full intra-atomic interaction is important to obtain correct results.

We study the Gutzwiller method for the spinless fermion model in one dimension, which is one of the simplest models that incorporates the intersite Coulomb interaction. The Gutzwiller solution of this model has been studied in the literature but with differing results. We obtain the Gutzwiller solution of the problem by a careful enumeration of the many-particle configurations and explain the origin of the discrepancy in the literature to be due to the neglect of the correlation between the neighboring bond occupancies. The correct implementation of the Gutzwiller approach yields results different from the slave-boson mean-field theory, unlike for the Hubbard model with on-site interaction, where both methods are known to be equivalent. The slave-boson and the Gutzwiller results are compared to the exact solution, available for the half-filled case, and to the numerical exact diagonalization results for the case of general filling.

The minimum of the Gutzwiller energy functional depends on the number of parameters considered in the variational state. For a three-orbital Hubbard model we find that the frequently used diagonal Ansatz is very accurate in high-symmetry situations. For lower symmetry, induced by a crystal-field splitting or the spin-orbit coupling, the discrepancies in energy between the most general and a diagonal Gutzwiller Ansatz can be quite significant. We discuss approximate schemes that may be employed in multi-band cases where a minimization of the general Gutzwiller energy functional is too demanding numerically.

Abstract We introduce a rotationally invariant approach combined with the Gutzwiller conjugate gradient minimization method to study correlated electron systems. In the approach, the Gutzwiller projector is parametrized based on the number of electrons occupying the onsite orbitals instead of the onsite configurations. The approach efficiently groups the onsite orbitals according to their symmetry and greatly reduces the computational complexity, which yields a speedup of 20 ∼ 50 × in the minimal basis energy calculation of dimers. The computationally efficient approach promotes more accurate calculations beyond the minimal basis that is inapplicable in the original approach. A large-basis energy calculation of F2 demonstrates favorable agreements with standard quantum-chemical calculations Bytautas et al (2007 J. Chem. Phys. 127 164317)

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

Z_2

symmetry breaking that can be identified as the Mott insulator-to-metal transition.

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

{P}_{R}

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.