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.
We give a comprehensive introduction into a diagrammatic method that allows for the evaluation of Gutzwiller wave functions in finite spatial dimensions. We discuss in detail some numerical schemes that turned out to be useful in the real-space evaluation of the diagrams. The method is applied to the problem of d-wave superconductivity in a two-dimensional single-band Hubbard model. Here, we discuss in particular the role of long-range contributions in our diagrammatic expansion. We further reconsider our previous analysis on the kinetic energy gain in the superconducting state.
In this work we analyze the variational problem emerging from the Gutzwiller approach to strongly correlated systems. This problem comprises the two main steps: evaluation and minimization of the ground state energy $W$ for the postulated Gutzwiller Wave Function (GWF). We discuss the available methods for evaluating $W$, in particular the recently proposed diagrammatic expansion method. We compare the two existing approaches to minimize $W$: the standard approach based on the effective single-particle Hamiltonian (EH) and the so-called Statistically-consistent Gutzwiller Approximation (SGA). On the example of the superconducting phase analysis we show that these approaches lead to the same minimum as it should be. However, the calculations within the SGA method are easier to perform and the two approaches allow for a simple cross-check of the obtained results. Finally, we show two ways of solving the equations resulting from the variational procedure, as well as how to incorporate the condition for a fixed number of particles.
Partially-projected Gutzwiller variational wavefunctions are used to describe the ground state of disordered interacting systems of fermions. We compare several different variational ground states with the exact ground state for disordered one-dimensional chains, with the goal of determining a minimal set of variational parameters required to accurately describe the spatially-inhomogeneous charge densities and spin correlations. We find that, for weak and intermediate disorder, it is sufficient to include spatial variations of the charge densities in the product state alone, provided that screening of the disorder potential is accounted for. For strong disorder, this prescription is insufficient and it is necessary to include spatially inhomogeneous variational parameters as well.
We use the Gutzwiller variational theory to calculate the ground-state phase diagram and quasi-particle bands of LaOFeAs. The Fe3d--As4p Wannier-orbital basis obtained from density-functional theory defines the band part of our eight-band Hubbard model. The full atomic interaction between the electrons in the iron orbitals is parameterized by the Hubbard interaction U and an average Hunds-rule interaction J. We reproduce the experimentally observed small ordered magnetic moment over a large region of (U,J) parameter space. The magnetically ordered phase is a stripe spin-density wave of quasi-particles.
We develop the Gutzwiller approximation method to obtain the renormalized Hamiltonian of the SU(4) $t$-$J$ model, with the corresponding renormalization factors. Subsequently, a mean-field theory is employed on the renormalized Hamiltonian of the model on the honeycomb lattice under the scenario of a cooperative condensation of carriers moving in the resonating valence bond state of flavors. In particular, we find the extended $s$-wave superconductivity is much more favorable than the $d+id$ superconductivity in the doping range close to quarter filling. The pairing states of the SU(4) case reveal the property that the spin-singlet pairing and the spin-triplet pairing can exist simultaneously. Our results might provide new insights into the twisted bilayer graphene system.