Asymptotic long-range order for the XY-model on random geometric graphs

We study the classical $XY$-model on random geometric graphs $\mathcal{G}_{n, \varepsilon}$, which are obtained by sampling $n \in \mathbb{N}$ independent points in a finite domain $Ω\subset \mathbb{R}^d$, $d \geq 2$, and connecting two points by and edge if their distance is of order $\varepsilon >