  In \cite{BS}, Balasubramanyam and the second named author derived the first
moment of the pair correlation function for Hecke angles lying in small
subintervals of $[0,1]$ upon averaging over large families of Hecke newforms of
weight $k$ with respect to $\Gamma_0(N)$. The goal of this article is to study
higher moments of this pair correlation function. For an integer $r \geq 2$, we
present bounds for its $r$-th power moments. We apply these bounds to record
lower order error terms in the computation of the second and third moments. As
a result, one can obtain the convergence of the second and third moments of
this pair correlation function for suitably small intervals, and under
appropriate growth conditions for the size of the families of Hecke newforms.

Number Theory