Equidistribution of integers represented by standard quadratic form under arithmetic constraints

We study the equidistribution of integers of the form $n= x_1^2 + \cdots + x_d^2$ under the arithmetic constraints given by $(\mathbb{Z}/p\mathbb{Z})^d$. The first step in addressing this problem is to construct modular forms whose Fourier expansion coefficients correspond to the counting problem ov