  Lattice measurements of spatial distributions of the light quark bilinear
densities in static mesons allow to test directly and in detail the wave
functions of quark models. These distributions are gauge invariant quantities
directly related to the spatial distribution of wave functions. We make a
detailed comparison of the recent lattice QCD results with our own quark
models, formulated previously for quite different purposes. We find a striking
agreement not only between our two quark models, but also with the lattice QCD
data for the ground state in an important range of distances up to about 4/GeV.
Moreover the agreement extends to the L=1 states [j^P=(1/2)^+]. An explanation
of several particular features completely at odds with the non-relativistic
approximation is provided. A rather direct, somewhat unexpected and of course
approximate relation between wave functions of certain quark models and QCD has
been established.
