In this paper we wish to study hadronic wave functions using an instan
ton model for the QCD vacuum. The wave functions are defined in terms
of gauge-invariant Bethe-Salpeter amplitudes which we have determined
numerically using a Monte Carlo simulation of the instanton ensemble.
We find that the pion and the proton, as well as the rho meson and the
DELTA, have very similar wave functions but observe a sizable splitti
ng between mesons or baryons with different spin. We compare our resul
ts with data obtained in lattice gauge simulations.