The distribution of spacings between quadratic residues

We study the distribution of spacings between squares modulo q, where q is square-free and highly composite, in the limit as the number of prime factors of q goes to infinity. We show that all correlation functions are Poissonian, which among other things, implies that the spacings between nearest neighbors, normalized to have unit mean, have an exponential distribution.