  The theory of Maximum Probability Domains (MPDs) is formulated for the
Hubbard model in terms of projection operators and generating functions for
both exact eigenstates as well as Slater determinants. A fast MPD analysis
procedure is proposed, which is subsequently used to analyse numerical results
for the Hubbard model. It is shown that the essential physics behind the
considered Hubbard models can be exposed using MPDs. Furthermore, the MPDs
appear to be in line with what is expected from Valence Bond Theory-based
knowledge.
