Multi-band Gutzwiller-correlated wave functions reconcile the contrasting concepts of itinerant band electrons versus electrons localized in partially filled atomic shells. The exact evaluation of these variational ground states in the limit of large coordination number allows the identification of quasi-particle band structures, and the calculation of a variational spinwave dispersion. The study of a generic two-band model elucidates the co-operation of the Coulomb repulsion and the Hund's-rule exchange for itinerant ferromagnetism. We present results of calculations for ferromagnetic nickel, using a realistic 18 spin-orbital basis of $4s$, $4p$ and $3d$ valence electrons. The quasi-particle energy bands agree much better with the photo-emission and Fermi surface data than the band structure obtained from spin-density functional theory (SDFT).