Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1107.1354 · Physical Review Letters

Time-Dependent Gutzwiller Theory for Multiband Hubbard Models

Ernst von Oelsen, G. Seibold, Jörg Bünemann

Published Aug 9, 2011• 18 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevlett.107.076402

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

Linear-Response Dynamics from the Time-Dependent Gutzwiller Approximation

Massimo Capone, J. Lorenzana, G. Seibold

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.

•Jan 1, 2014•20 citations

Gutzwiller variational wave function for multiorbital Hubbard models in finite dimensions

Kevin zu Münster, Jörg Bünemann

We develop a diagrammatic method for the evaluation of general multiband 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.

arXiv:1510.07896•Physical review. B./Physical review. B•Jul 26, 2016•13 citations

Linear-response dynamics from the time-dependent Gutzwiller approximation

J Bünemann, M Capone, J. Lorenzana +1 more

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.

arXiv:1303.1665•New Journal of Physics•May 30, 2013•23 citations

Time-dependent Gutzwiller approximation: Recent advances

G. Seibold, Jörg Bünemann, Ernst von Oelsen +1 more

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)

•Sep 4, 2012•0 citations

Time-Dependent Gutzwiller Approximation for the Hubbard Model

G. Seibold, J. Lorenzana

We develop a time-dependent Gutzwiller approximation (GA) for the Hubbard model analogous to the time-dependent Hartree-Fock (HF) method. This new formalism incorporates ground state correlations of the random phase approximation (RPA) type beyond the GA. Static quantities like ground state energy and double occupancy are in excellent agreement with exact results in one dimension up to moderate coupling and in two dimensions for all couplings. We find a substantial improvement over traditional GA and HF+RPA treatments. Dynamical correlation functions can be computed and are also substantially better than HF+RPA ones and obey well behaved sum rules.

arXiv:cond-mat/0006416•Physical Review Letters•Mar 19, 2001•85 citations
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 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.
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.
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)
We develop a time-dependent Gutzwiller approximation (GA) for the Hubbard model analogous to the time-dependent Hartree-Fock (HF) method. This new formalism incorporates ground state correlations of the random phase approximation (RPA) type beyond the GA. Static quantities like ground state energy and double occupancy are in excellent agreement with exact results in one dimension up to moderate coupling and in two dimensions for all couplings. We find a substantial improvement over traditional GA and HF+RPA treatments. Dynamical correlation functions can be computed and are also substantially better than HF+RPA ones and obey well behaved sum rules.