Voronoi limit measures for iterates of constant-coefficient differential operators on rational functions with simple poles

Abstract: Bøgvad and Hägg proved that for a rational function with simple poles, the zeros of successive derivatives accumulate on the Voronoi diagram of the pole set, and the normalized zero-counting measures converge to a canonical probability measure supported on this diagram. We extend this result from pure derivatives to iterates of an arbitrary monic constant-coefficient differential operator. Let $h(z)=A(z)/B(z)$ be a reduced rational function, where $B$ is monic of degree $b\ge2$ with distinct zeros $S={z_1,\dots,z_b}$, and let $P(D)=\sum_{j=0}^m c_jD<sup>j$ be a monic constant-coefficient differential operator of order $m\ge1$. After clearing denominators, we can write $P(D)^{n(h)=\widetilde} A_n/B<sup>{mn+1}$ and study the zeros of the numerator polynomials $\widetilde A_n$. If $r:=\min{j:c_j\neq0}$, then (after passing to the proper part of $h$ when $r&gt;0$) the associated zero-counting measures converge vaguely to $\frac{m(b-1)}{bm-r}\,μ<em>S,$ where $μ_S$ is the Bøgvad--Hägg probability measure supported on the Voronoi diagram $V_S$. In particular, the limit is a probability measure exactly when $P(D)=D<sup>m$; otherwise a proportion $\frac{m-r}{bm-r}$ of zeros escapes to infinity (in the sense of vague convergence). When $r&lt;m$, the unshifted logarithmic potentials diverge, but an explicit factorial renormalization yields $L¹</em>{\mathrm{loc}}(\mathbb C)$ convergence to a subharmonic limit with Riesz measure $\frac{m(b-1)}{bm-r}\,μ_S$. Apart from this scalar factor, the limiting measure is determined solely by the pole configuration; the coefficients of $P(D)$ affect only an additive constant in the limiting potential.