Intersections of Hecke correspondences on modular curves

We compute the arithmetic intersections of Hecke correspondences on the product of integral model of modular curve $\mathcal{X}_0(N)$ and relate it to the derivatives of certain Siegel Eisenstein series when $N$ is odd and squarefree. We prove this by establishing a precise identity between the arithmetic intersection numbers on the Rapoport--Zink space associated to $\mathcal{X}_0(N)^{2}$ and the derivatives of local representation densities of quadratic forms.