Geodesic clustering of zeros of Eisenstein series for congruence groups

We consider a set of generators for the space of Eisenstein series of even weight $k$ for any congruence group $Γ$ and study the set of all of their zeros taken for $Γ(1)$-conjugates of $Γ$ in the standard fundamental domain for $Γ(1)$. We describe (a) an upper bound $κ_Γ+ O(1/k)$ for their imaginary part; (b) a finite configuration of geodesics segments to which all zeros converge in Hausdorff distance as $k \rightarrow \infty$; (c) a finite set containing all algebraic zeros for all weights. The bound in (a) depends on the (non-)vanishing of a new generalization of Ramanujan sums. The proof of (b) originates in a method used to study phase transitions in statistical physics. The proof of (c) relies on the theory of complex multiplication. The results can be made quantitative for specific groups. For $Γ=Γ(N)$ with $4 \nmid N$, $κ_Γ=1$ and the zeros tend to the unit circle, whereas if $4 \mid N$, $κ_Γ=2$ and the limit configuration includes parts of vertical geodesics and circles of radius $2$. In both cases, the only algebraic zeros are at $\mathrm{i}$ and $\exp(2π\mathrm{i}/3)$ for sufficiently large $k$. For $Γ(N)$ with $N$ odd, we use finer estimates to prove a trichotomy for the exact `convergence speed' of the zeros to the unit circle, as well as angular equidistribution of the zeros as $k \rightarrow \infty$.