Properties of Gutzwiller wave functions for multiband models

We analyze the Mott transition in multiband Hubbard models with the inclusion of multiplet exchange splittings as it arises in infinite dimensions by using the generalized Gutzwiller wave function introduced by Bunemann, Weber, and Gebhard [Phys. Rev. B 57, 6896 (1998)]. We also present an extension of that variational wave function to account for broken-symmetry solutions, which still allows an exact analytical treatment. Our analysis reveals some drawbacks of the variational wave function, which, in our opinion, imply that Gutzwiller-type wave functions do not properly characterize quasiparticles close to a Mott transition even in the limit of infinite dimensions.