Weiyang Fang
Recent work gives an all-$β$, all-order stochastic-zeta representation of the correlation functions of the $\Sine_β$ process and determines their leading Vandermonde asymptotics when several variables merge. We compute the first nontrivial correction throughout the regime $mβ>1$. If $a_1,\ldots,a_m$ are distinct real numbers and \[ V(a)=\sum_{1\le i<j\le m}(a_i-a_j)^2, \] then, as $\varepsilon\to0$, \[ ρ^{(m)}_β(\varepsilon a_1,\ldots,\varepsilon a_m) =C^{(m)}_β|\varepsilon|^{β\binom m2}\prod_{i1/2$. The proof combines a finite-$N$ rotational Ward identity, exact Hua--Pickrell trace moments, compact moment bounds for the stochastic-zeta entire function and its derivatives, and a quantitative multivariate expectation--Taylor lemma. As a by-product we evaluate \[ \E_{\HP_{β,mβ/2}}\sum_x x^{-2}=\frac{mβ}{4(mβ-1)(2m+1)}. \] The pole at $mβ=1$ marks the boundary of the present second-moment argument and suggests a transition in the form of the next fusion correction.