  For fixed $\alpha \in [0,1]$, consider the set $S_{\alpha,N}$ of dilated
squares $\alpha, 4\alpha, 9\alpha, \dots, N^2\alpha \, $ modulo $1$. Rudnick
and Sarnak conjectured that for Lebesgue almost all such $\alpha$ the
gap-distribution of $S_{\alpha,N}$ is consistent with the Poisson model (in the
limit as $N$ tends to infinity). In this paper we prove a new estimate for the
triple correlations associated to this problem, establishing an asymptotic
expression for the third moment of the number of elements of $S_{\alpha,N}$ in
a random interval of length $L/N$, provided that $L > N^{1/4+\varepsilon}$. The
threshold of $1/4$ is substantially smaller than the threshold of $1/2$ (which
is the threshold that would be given by a na\"{i}ve discrepancy estimate).
  Unlike the theory of pair correlations, rather little is known about triple
correlations of the dilations $(\alpha a_n \, \text{mod } 1)_{n=1}^{\infty} $
for a non-lacunary sequence $(a_n)_{n=1}^{\infty} $ of increasing integers.
This is partially due to the fact that second moment of the triple correlation
function is difficult to control, and thus standard techniques involving
variance bounds are not applicable. We circumvent this impasse by using an
argument inspired by works of Rudnick--Sarnak--Zaharescu and Heath-Brown, which
connects the triple correlation function to some modular counting problems.
  In an appendix we comment on the relationship between discrepancy and
correlation functions, answering a question of Steinerberger.
