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