\documentclass[11pt]{amsart} \usepackage[margin=1.1in]{geometry} \usepackage{amsmath,amssymb,amsthm} \usepackage[colorlinks=true,linkcolor=blue,citecolor=blue,urlcolor=blue]{hyperref} \newtheorem{theorem}{Theorem}[section] \newtheorem{corollary}[theorem]{Corollary} \newtheorem{proposition}[theorem]{Proposition} \newtheorem{lemma}[theorem]{Lemma} \newtheorem{theoremA}{Theorem} \renewcommand{\thetheoremA}{A} \theoremstyle{definition} \newtheorem{definition}[theorem]{Definition} \newtheorem{hypothesis}[theorem]{Hypothesis} \newtheorem{question}[theorem]{Question} \theoremstyle{remark} \newtheorem{remark}[theorem]{Remark} \numberwithin{equation}{section} \newcommand{\C}{\mathbb{C}} \newcommand{\R}{\mathbb{R}} \newcommand{\N}{\mathbb{N}} \newcommand{\Q}{\mathbb{Q}} \newcommand{\Z}{\mathbb{Z}} \newcommand{\D}{\mathbb{D}} \newcommand{\E}{\mathbb{E}} \newcommand{\PP}{\mathbb{P}} \newcommand{\eps}{\varepsilon} \newcommand{\He}{\mathrm{He}} \DeclareMathOperator{\re}{Re} \DeclareMathOperator{\im}{Im} \DeclareMathOperator{\dist}{dist} \DeclareMathOperator{\sech}{sech} \DeclareMathOperator{\arcosh}{arcosh} \title[Explicit onsets, rates and the tail of the collision constant]{Explicit onsets, rates and the tail of the collision constant for heat-evolved Weyl polynomials} \author{B. Baggins} \date{Draft, October 2026} \subjclass[2020]{Primary 26C10, 60G55; Secondary 30C15, 33C45, 35K05, 60G15} \keywords{Hermite polynomials, linear combinations of orthogonal polynomials, nonreal zeros, explicit onset, heat flow, Weyl polynomials, random power series, rate of convergence, collision of zeros, correlation functions of zeros} \begin{document} \begin{abstract} Let $h_k$ be the orthonormal Hermite polynomials for the standard Gaussian weight, and let $a(w)=\sum_{j\ge0}a_jw^j$ be a real power series with radius of convergence greater than $1$, $a(0)\ne0$, no zeros on the unit circle and only simple zeros in the unit disk. For all large $n$ the polynomial $\sum_{j=0}^{n}a_jh_{n-j}$ has exactly as many nonreal zeros as $a$ has in the disk, but the known proof gives the onset by compactness. We make it explicit: if $a$ has no zeros in $\{1-d\le|w|\le1\}$, with $d$ below an explicit threshold, the count is exact as soon as $n$ exceeds a polynomial in $1/d$ and in the inverses of $|a(0)|$, $\min_{|w|=1}|a|$, $\min_{|w|=1-d}|a|$ and $\min_{[-1+d,1-d]}\max(|a|,|a'|)$; no separation of the nonreal zeros of $a$ is needed. By an exact identity this applies to the random Weyl polynomials $G_n$ with independent standard Gaussian coefficients under the backward heat flow $G_{n,t}=\exp(-\frac t{2n}\partial_z^2)G_n$, whose real-rooting time converges in distribution to $T_*=R_*^{-2}$, where $R_*$ is the smallest modulus of a nonreal zero of the Gaussian power series $q(w)=\sum_j\eta_jw^j$. We show that the convergence holds with the rate $n^{-\alpha}$ on $[1+\eps,\infty)$ for every $\alpha<\frac3{50}$, an exponent that does not depend on $\eps$; the dependence on $\eps$ is removed by a super-exponential tail for the number of zeros of $q$ in a disk. Under a coupling of the degrees that shares the top coefficients, the real-rooting time $\Lambda_n$ satisfies $|\Lambda_n-T_*|\le K/n$ for $n\ge n_0$, with $K$ and $n_0$ depending on $q$. We bound $K$ and $n_0$ by polynomials in condition parameters of $q$, with absolute degrees and, when $T_*\ge1+\eps$, coefficients polynomial in $1/\eps$. For the Gaussian power series this gives, on the event $\{T_*\ge1+\eps\}$, a power tail with exponent $1/900$ for the largest scaled displacement $D$ of the final collision of zeros, and without restriction the tail bound $\exp(-c(\log x/\log\log x)^{1/2})$; the inputs include closed forms of the two-point correlation functions of the zeros of $q$. \end{abstract} \maketitle \section{Introduction}\label{sec:intro} \subsection{Setting} This paper makes quantitative three results of the main paper \cite{Main}, to which we refer for background and notation. Let $\xi_0,\xi_1,\dots$ be independent standard real Gaussian variables, let $G_n(z)=\sum_{k=0}^{n}\xi_k(\sqrt n\,z)^k/\sqrt{k!}$ be the real Weyl polynomial of degree $n$, and let $G_{n,t}=\exp(-\frac t{2n}\partial_z^2)G_n$, $t\ge0$, be its evolution under the backward heat flow. The set of times at which $G_{n,t}$ has only real zeros is an interval $[\Lambda_n^{\rm std},\infty)$, and $\Lambda_n^{\rm std}$ is the \emph{real-rooting time}. Let $\eta_0,\eta_1,\dots$ be independent standard real Gaussian variables, let $q(w)=\sum_{j\ge0}\eta_jw^j$ be the \emph{Gaussian power series}, let $R_*$ be the smallest modulus of a nonreal zero of $q$ in $\D=\{|w|<1\}$, and $T_*=R_*^{-2}$. By \cite[Theorem~1.1]{Main}, $\Lambda_n^{\rm std}\to T_*$ in distribution; moreover, for every $t>1$, the number of nonreal zeros of $G_{n,t}$ converges in distribution to the number of nonreal zeros of $q$ in $\{|w|0$ there is $n_0$ such that for $n\ge n_0$ all zeros of $z\mapsto Q_n(\sqrt n\,z)$, where $Q_n=\sum_{j=0}^{n}a_jh_{n-j}$, are simple; exactly one of them lies within $\eps$ of each point $w_i+w_i^{-1}$, and it is real if and only if $w_i$ is real; and the other $n-k$ are real and lie in $(-2-\eps,2+\eps)$. In particular $Q_n$ has exactly as many nonreal zeros as $a$ has in $\D$. By the identity \cite[\textup{(1.6)}]{Main}, $P_{n,t}(z)$ is, up to the factor $t^{n/2}$ and the dilation $z\mapsto z/\sqrt t$, the polynomial $Q_n(\sqrt n\,\cdot)$ for the symbol $a^{(t)}(w)=q(w/\sqrt t)$, which has radius of convergence $\sqrt t>1$ when $t>1$. Both $n_0$ in \cite[Theorem~1.4]{Main} and the constants of the theorem on the final collision \cite[Theorem~9.1]{Main} are obtained by compactness. We bound them explicitly, and we draw two consequences for the Gaussian case: a rate of convergence in \cite[Theorem~1.1]{Main}, and a tail bound for the constant of the final collision. \subsection{An explicit onset} \begin{theorem}[An explicit onset]\label{thm:quant-intro} Let $a(w)=\sum_{j\ge0}a_jw^j$ be a real power series with $|a_j|\le C_0R_0^{-j}$ for all $j\ge0$, where $R_0>1$, such that $a(0)\ne0$ and $a$ has no zeros on the unit circle. Let $01$, for the probability that $G_{n,t}$ has only real zeros and for the law of the number of its nonreal zeros. It rests on Theorem~\ref{thm:quant-intro} and on a bound for the number of zeros of $q$ in a disk (Lemma~\ref{lem:overcrowd}), which makes the exponent independent of $\eps$. The exponent is far from optimal (Remark~\ref{rem:rate}), and we have no rate uniform in $t>1$. \subsection{Explicit constants and the tail of the collision constant} Let $q(w)=\sum_j\eta_jw^j$ now be a deterministic real power series that satisfies the conditions (G1)--(G3) of \cite[Section~9]{Main}: $q(0)\ne0$, the nonreal zeros of $q$ of minimal modulus $R_*$ inside its disk of convergence are one simple conjugate pair $u_*=R_*e^{i\theta_*}$, $\bar u_*$, with $0<\theta_*<\pi$, $q$ has no real zero of modulus $R_*$, and its real zeros of smaller modulus are simple. Let $T_*=R_*^{-2}$, and let $P_{n,t}$ and $\Lambda_n$ be built from the coefficients of $q$ as above. By \cite[Theorem~9.1]{Main} there are $\delta>0$, $K<\infty$ and $n_0$ such that for $n\ge n_0$: $|\Lambda_n-T_*|\le K/n$; all zeros of $P_{n,\Lambda_n}$ are real, and exactly one of them, $z_n$, is multiple, with $|z_n-2R_*^{-1}\cos\theta_*|\le K/n$; for $T_*-\delta\le t<\Lambda_n$ the polynomial $P_{n,t}$ has exactly two nonreal zeros; and $\Lambda_n$ is the only $t\ge T_*-\delta$ at which $P_{n,t}$ has a multiple zero. For the Gaussian power series (G1)--(G3) hold almost surely \cite[Corollary~9.14]{Main}, so that the largest scaled displacement of the collision, \[ D=\sup_{n\ge1}n\bigl(|\Lambda_n-T_*|+|z_n-2R_*^{-1}\cos\theta_*|\bigr), \] where for $nx,\ T_*\ge1+\eps)\le c\,\eps^{-c}x^{-\beta_\eps}\quad\text{with}\quad\beta_\eps\ge\frac{\eps}{35640\log(40/\eps)}, \] and $\PP(D>x,\ T_*\ge1+\eps)\le C_\eps\,x^{-1/900}$, where $\log C_\eps\le c\,\eps^{-c}$. \item[(c)] There are absolute constants $c>0$ and $C$ such that for the Gaussian power series and all $x\ge e^e$ \[ \PP(D>x)\le C\exp\Bigl(-c\Bigl(\frac{\log x}{\log\log x}\Bigr)^{1/2}\Bigr). \] \end{enumerate} \end{theorem} Part~(a) is Theorem~\ref{thm:collision-explicit}, part~(b) is Theorem~\ref{thm:Dtail}, and part~(c) is Corollary~\ref{cor:Dtail}. In (a) the exponents are absolute; the factor exponential in the number of zeros is kept as a parameter because for random $q$ it can be bounded either by Jensen's formula, which gives the exponent $\beta_\eps$ of order $\eps/\log(1/\eps)$ in (b), or through the number of zeros of $q$ in a disk, which gives the exponent $1/900$, independent of $\eps$, but a constant $C_\eps$ that is not polynomial in $1/\eps$. The proof of (b) uses closed forms of the two-point correlation functions of the zeros of $q$ (Proposition~\ref{prop:twopoint}), which give a small-ball bound for the gap $g$ between $R_*$ and the other moduli of zeros (Theorem~\ref{thm:gap}). Part~(c) combines the first bound of (b) with the bound for $\PP(T_*<1+\eps)$ that follows from \cite[Theorem~7.13]{Main}. Whether $\PP(D>x)$ decays like a power of $x$ is open (Question~\ref{q:Dtail}); a small-ball lower bound for $\im u_*$, which we expect but have not proved, would exclude exponents above $1$ (Remark~\ref{rem:Dtail}(3)). \subsection{Ideas of the proofs} \emph{The explicit onset (Section~\ref{subsec:onset}).} The proof of \cite[Theorem~1.4]{Main} deforms $a$, through $a_s(w)=p_{\rm in}(w)\,p_{\rm out}(sw)$, $s\in[0,1]$, into the polynomial $p_{\rm in}$ with the same zeros in $\D$, and excludes real double zeros along the way by a bulk estimate: in the bulk, $(Q_n,Q_n'/\sqrt n)$ is a $2\times2$ matrix, built from $O(\log n)$ steps of the backward three-term recurrence, applied to $(h_n,h_{n-1})$, plus a tail controlled by a Christoffel--Darboux bound. Each step of that estimate has explicit constants once the coefficient bound of $a_s$ and the minimum of $|a_s|$ on the unit circle are controlled uniformly in $s$. Lemma~\ref{lem:outer} controls both: the coefficients of $a_s$ are bounded by a constant times $\Upsilon^k$, and $|a_s|\ge\mathfrak d(a)\,|a|$ on the closed unit disk, where the distortion $\mathfrak d(a)$ of the outer factor $p_{\rm out}$ is at least $\min_{|w|=1}|a|$ times $\Upsilon^{-k}$ divided by a coefficient bound; here $\Upsilon$ depends only on $R_0$, and Jensen's formula bounds the number $k$ of zeros of $a$ in $\D$. In the exterior we compare $Q_n(\sqrt n\,\cdot)/h_n(\sqrt n\,\cdot)$, as a function of $w$ with $z=w+w^{-1}$, with $a$ on the single circle $\{|w|=1-d\}$, and inside it we pass from $a$ to this function by a linear homotopy, along which real double zeros are excluded by $\nu_d$; this is why no separation of the zeros of $a$ enters. \emph{The rate (Section~\ref{subsec:rate}).} For the random symbol $a^{(t)}=q(\cdot/\sqrt t)$, each parameter in the conditions of Theorem~\ref{thm:quant-intro} has a polynomial small-ball bound, uniformly in $t\ge1+\eps$: the coefficient bound, the minimum of $|q|$ on a circle, the near-double real zeros and the zeros near the circle (Lemmas~\ref{lem:coeff}--\ref{lem:annulus}). The factor exponential in the number of zeros of $a^{(t)}$ in $\D$ would cost, through Jensen's formula, an exponent of order $\eps/\log(1/\eps)$; instead we use that the number of zeros of $q$ in a disk has super-exponential tails (Lemma~\ref{lem:overcrowd}). Large times are handled by the tail of $T_*$. \emph{Explicit constants and the tail of $D$ (Section~\ref{sec:explicit}).} Each compactness step of the proof of \cite[Theorem~9.1]{Main} can be made explicit, and Theorem~\ref{thm:quant} replaces \cite[Theorem~1.4]{Main} at the two times $T_*\pm\delta$. This bounds $\delta^{-1}$, $K$ and $n_0$ by powers of $\Pi$, times powers of $1/\eps$. For the Gaussian power series each factor of $\Pi$ has a power tail on $\{T_*\ge1+\eps\}$. For the modulus gap $g$ this is a first-moment bound over pairs of zeros: the two-point correlation functions, which we write in closed form through the Gaussian regression on the zeros, are at most the products of the one-point intensities. The factor exponential in the number of zeros is bounded through the super-exponential tail of that number, as in the proof of Theorem~\ref{thm:rate-intro}. \subsection{Related work} For the literature on zeros of linear combinations of orthogonal polynomials, on the heat flow of random polynomials and on the de Bruijn--Newman constant we refer to the introduction of \cite{Main}. Times $t=t_n$ that tend to $1$ are the subject of \cite{CompB}, and symbols that depend on $n$, varying weights and the Laguerre flow that of \cite{CompC}. The one-point intensities of the zeros of the Gaussian power series go back to Kac \cite{Ka} for the real zeros, and the correlation functions are Pfaffians, by Matsumoto and Shirai \cite{MS} and, for the mixed correlations of real and nonreal zeros, by Forrester \cite{Fo}. We need only the two-point functions, which we write in closed form (Proposition~\ref{prop:twopoint}), and the Rice formulas of Aza\"is and Wschebor \cite{AW}. \subsection{Organisation} Section~\ref{subsec:onset} proves Theorem~\ref{thm:quant-intro} (Theorem~\ref{thm:quant}) and its form with localisation (Theorem~\ref{thm:quant-main}). Section~\ref{subsec:rate} proves Theorem~\ref{thm:rate-intro} (Theorem~\ref{thm:rate}) and bounds the distortion of the outer factor for the random symbol (Proposition~\ref{prop:distortion}). Section~\ref{sec:explicit} proves Theorem~\ref{thm:explicit-intro}: the explicit constants in Section~\ref{subsec:collision-explicit}, the two-point functions and the modulus gap in Section~\ref{subsec:pairs}, small-ball bounds at the collision in Section~\ref{subsec:smallballs}, and the tail of $D$ in Section~\ref{subsec:Dtail}. Section~\ref{sec:remarks} contains remarks and open questions. \section{An explicit onset for the exact count}\label{subsec:onset} The proof of \cite[Theorem~1.4]{Main} in \cite[Section~6]{Main} obtains $n_0$ by compactness. To make it explicit, we control the deformation $a_s(w)=p_{\rm in}(w)\,p_{\rm out}(sw)$ of that proof through the zeros of $a$ in $\D$ (Lemma~\ref{lem:outer}), and we treat the exterior with one circle and one real segment instead of the disks about the points $\zeta_i=w_i+w_i^{-1}$ used there, so that the count needs no separation of the zeros. \begin{definition}[Condition parameters]\label{def:condition} For $C_0<\infty$ and $R_0>1$ let $\mathcal S(C_0,R_0)$ be the set of real power series $a(w)=\sum_{j\ge0}a_jw^j$ with $|a_j|\le C_0R_0^{-j}$ for all $j\ge0$, and put $Q_n[a]=\sum_{j=0}^{n}a_jh_{n-j}$. For $a\in\mathcal S(C_0,R_0)$ and $00$ and $\mu>0$. Let $p_{\rm in}(w)=a(0)\prod_{i=1}^{k}(1-w/w_i)$, $p_{\rm out}=a/p_{\rm in}$ and $G=\alpha/\prod_i|w_i|$. \begin{enumerate} \item[(a)] $|w_i|\ge r_a:=\alpha/(2C_0)$ and $1-|w_i|\ge d_a:=\mu R_0(1-R_0^{-1})^2/C_0$ for every $i$. \item[(b)] $k\log R_1\le\log(\overline M/\alpha)$; hence $\Upsilon^k\le(\overline M/\alpha)^E$. \item[(c)] $p_{\rm out}$ is analytic in $\{|w|1$, $\mu_S>0$, and let $(a_s)_{s\in S}$ be a family in $\mathcal S(C,R)$ such that $|a_s(w)|\ge\mu_S$ for all $s\in S$ and $|w|=1$. Put $\vartheta=R^{-1/2}$ and \begin{gather*} \kappa=\min\Bigl\{1,\ \tfrac12\log R,\ \frac{\mu_SR(1-\vartheta)^2}{6C}\Bigr\},\qquad\Gamma=\frac{8CR^{1/2}}{(1-\vartheta)^2},\\ M_n=\Bigl\lceil\frac{2\log n}{\log R}\Bigr\rceil,\qquad D_H=\frac{C(198/\kappa+4M_n+8)}{(1-\vartheta)^4}. \end{gather*} If $n\ge2$ satisfies \begin{equation}\label{eq:bulk-explicit} \begin{gathered} n\ge\frac{(M_n+2)^2}{\kappa},\qquad n\ge D_H,\qquad n\ge\frac{8(2\Gamma+1)D_H}{\mu_S^2},\\ n^{3/2}\ge\frac{37(\Gamma+1)C}{(1-R^{-2})^{1/2}\mu_S^2}, \end{gathered} \end{equation} then for all $s\in S$ and all real $x$ with $|x|\le2\sqrt n\cosh\kappa$, \[ Q_n[a_s](x)^2+\frac{Q_n[a_s]'(x)^2}n\ \ge\ c^2\bigl(h_n(x)^2+h_{n-1}(x)^2\bigr),\qquad c=\frac{\mu_S^2}{16(\Gamma+1)} . \] \end{lemma} \begin{proof} We follow the proof of \cite[Proposition~4.3]{Main}, with $L=2\cosh\kappa$ and $\kappa_0=\frac12\log R$, so that $e^{\kappa_0}/R=\vartheta$ and $\kappa\le\kappa_0$. \emph{The limit matrix.} By \cite[Lemma~4.2]{Main} with $\kappa_0$, $\det H_{a_s}(z)\ge\mu_S^2$ for $|z|\le2$, and $\det H_{a_s}(z)=a_s(w)a_s(w^{-1})$ for $2<|z|\le L$, $z=w+w^{-1}$, $w$ real, $e^{-\kappa}\le|w|<1$. Since $\max_{|u|\le e^{\kappa_0}}|a_s'(u)|\le C\sum_jjR^{-j}e^{(j-1)\kappa_0}=CR^{-1}(1-\vartheta)^{-2}$, for real $v$ with $e^{-\kappa}\le|v|\le e^\kappa$ and $\varsigma=\operatorname{sign}v$, \[ |a_s(v)-a_s(\varsigma)|\le(e^\kappa-1)\,CR^{-1}(1-\vartheta)^{-2}\le e\kappa\,CR^{-1}(1-\vartheta)^{-2}\le\tfrac e6\,\mu_S<\tfrac12\mu_S . \] As $a_s(\varsigma)$ is real with $|a_s(\varsigma)|\ge\mu_S$, the real numbers $a_s(w)$ and $a_s(w^{-1})$ have the sign of $a_s(\varsigma)$ and modulus at least $\mu_S/2$. So $\det H_{a_s}\ge\mu_S^2/4$ on $[-L,L]$. By \cite[Lemma~3.3]{Main}, $\|G_j(z)\|\le4(j+2)e^{(j+1)\kappa}$, so $\|H_{a_s}(z)\|\le4Ce^\kappa\sum_j(j+2)\vartheta^j\le\Gamma$ for $|z|\le L$. \emph{The head.} The first condition gives $(M_n+2)^2\le\kappa n$ with $\kappa\le1$. So for $j\le M_n$ \cite[Lemma~3.4]{Main}, with $c_1=15/\kappa$, bounds the sum of the absolute values of the entries of the first row of $G_{n,j}-G_j$ by $c_1(j+1)^3e^{j\kappa}/n$ (part (a)) and that of the second row by $(c_1+3)(j+1)^3e^{j\kappa}/n$ (part (b)). Bounding the operator norm by the sum of the absolute values of the entries, $\|G_{n,j}(z)-G_j(z)\|\le(2c_1+3)(j+1)^3e^{j\kappa}/n\le(33/\kappa)(j+1)^3e^{j\kappa}/n$. For $j\ge2$, \cite[Lemma~2.5]{Main} gives $|A_j|\le(j-1)e^{(j-2)\kappa}$, $|B_j|,|A_{j+1}|\le je^{(j-1)\kappa}$ and $|B_{j+1}|\le(j+1)e^{j\kappa}$, so $\|G_j(z)\|\le4(j+1)e^{j\kappa}$. With $e^\kappa/R\le\vartheta$, $\vartheta^{M_n}\le1/n$, $\sum_j(j+1)^3\vartheta^j\le6(1-\vartheta)^{-4}$ and $\sum_{j>M}(j+1)\vartheta^j\le(M+2)\vartheta^{M+1}(1-\vartheta)^{-2}$, \[ \|H_{n,M_n}[a_s](z)-H_{a_s}(z)\|\le\frac{33C}{\kappa n}\sum_{j\ge0}(j+1)^3\vartheta^j+4C\sum_{j>M_n}(j+1)\vartheta^j\le\frac{D_H}n=:\delta_n \] for $|z|\le L$. The second and third conditions give $\delta_n\le1$ and $\delta_n(2\Gamma+1)\le\mu_S^2/8$. As in the proof of \cite[Proposition~4.3]{Main}, $\det H_{n,M_n}[a_s](z)\ge\mu_S^2/8$, $\|H_{n,M_n}[a_s](z)\|\le\Gamma+1$, and hence $s_{\min}(H_{n,M_n}[a_s](z))\ge2c$. \emph{The tail.} The tail estimate in the proof of \cite[Proposition~4.3]{Main} holds with \[ \eps_n=\sqrt2\,C(1-R^{-2})^{-1/2}R^{-M_n}\sqrt{n(1+L/2)}\le2.3\,C(1-R^{-2})^{-1/2}n^{-3/2}, \] since $R^{-M_n}=\vartheta^{2M_n}\le n^{-2}$ and $1+L/2\le1+\cosh1<2.55$. So $\eps_n\le c$ by the last condition, and by \cite[\textup{(4.1)}]{Main} the norm of the vector with entries $Q_n[a_s](x)$ and $Q_n[a_s]'(x)/\sqrt n$ is at least $(2c-\eps_n)\|v_n(x)\|\ge c\|v_n(x)\|$. \end{proof} \begin{lemma}[Exterior error on a disk]\label{lem:ext-disk} Let $a\in\mathcal S(C,R)$ with $R>1$. For $0<|w|<1$ let $\varphi_n(w)=F_n(\mathcal J(w))$, where $\mathcal J(w)=w+w^{-1}$ and $F_n=Q_n[a](\sqrt n\,\cdot)/h_n(\sqrt n\,\cdot)$. Then $\varphi_n$ extends analytically to $\D$, with $\varphi_n(0)=a(0)$ and $\varphi_n(\bar w)=\overline{\varphi_n(w)}$. For $00$ and $\mu>0$, and let $00$ and $\nu_d>0$. Put $C_0'=\overline M\,\Upsilon^k$ and $\mu'=\mathfrak d_0\,\mu$, and let $\kappa$, $\Gamma$, $M_n$ and $D_H$ be the constants of Lemma~\ref{lem:bulk-explicit} for $C=C_0'$, $R=R_1$ and $\mu_S=\mu'$; thus $\kappa=\min\{1,\frac14\log R_0,\mu'R_0^{1/2}(1-\vartheta_0)^2/(6C_0')\}$ with $\vartheta_0=R_0^{-1/4}$. Suppose that {\rm(i)} $d\le\kappa/\sqrt2$, and that $n\ge2$ satisfies \begin{enumerate} \item[(ii)] $n\ge k$, and \eqref{eq:bulk-explicit} holds for these constants; \item[(iii)] $n>2\sqrt2\,C_0'(1-R_0^{-1})^{-1/2}d^{-7/2}/(\mathfrak d_0\mu_d)$; \item[(iv)] $n>64\,C_0(1-R_0^{-2})^{-1/2}d^{-9/2}/\nu_d$. \end{enumerate} Then the number of nonreal zeros of $Q_n[a]$ equals the number of nonreal zeros of $a$ in $\D$, both counted with multiplicity, and all real zeros of $Q_n[a]$ are simple. Every zero $z$ of $Q_n[a](\sqrt n\,\cdot)$ that is not of the form $\mathcal J(w)$ with $|w|<1-d$ is real and lies in $(-2-2d^2,2+2d^2)$. \end{theorem} \emph{Size of the onset.} By Lemma~\ref{lem:outer}(b), $\Upsilon^k\le(\overline M/\alpha)^E$, so all conditions are polynomial in $1/\mu$, $1/\alpha$, $1/\mu_d$, $1/\nu_d$ and $1/d$, with exponents depending on $R_0$ only through $E$. For instance, normalising $C_0=1$ (the conclusion is invariant under $a\mapsto\lambda a$) and taking $d$ of order $\mu^2\Upsilon^{-2k}$, the largest value allowed by (i), conditions (iii) and (iv) read $n>C\Upsilon^{9k}\mu^{-8}\mu_d^{-1}$ and $n>C\Upsilon^{9k}\mu^{-9}\nu_d^{-1}$, with $C$ depending only on $R_0$. Neither the separation of the nonreal zeros of $a$ nor $|a'|$ at its zeros enters. Real zeros enter through $\nu_d$: if $x_10$ let $m(\rho)$ be the minimum of $|a(w)|$ over the $w$ with $|w|\le1$ and $|w-w_i|\ge\rho$ for all $i$. Put $\delta_0'=d_a\min(d_a,\sigma/2)$ and $\beta=1+16C_0^2/\alpha^2$. Let $0<\eps\le\delta_0'/3$, and let \[ \rho=\min\Bigl\{\frac\sigma3,\ \frac{d_a}3,\ \frac{r_a}2,\ \frac\eps\beta\Bigr\},\qquad d=\min\Bigl\{\frac12,\ \frac{d_a}3,\ \frac\kappa{\sqrt2},\ \sqrt{\eps/2}\Bigr\}. \] If $n$ satisfies (ii) of Theorem~\ref{thm:quant} and \[ n>\frac{2\sqrt2\,C_0'}{(1-R_0^{-1})^{1/2}\,d^{7/2}\,\mathfrak d_0\,m(\rho)}, \] then all zeros of $Q_n(\sqrt n\,\cdot)$ are simple, and the conclusions (a) and (b) of \cite[Theorem~1.4]{Main} hold for this $\eps$. \end{theorem} Since $m(\rho)\ge\mu2^{-k}\rho(\sigma/2)^{k-1}$ (Remark~\ref{rem:depth}), the onset is also polynomial in $1/\eps$. \begin{proof} By Lemma~\ref{lem:outer}(a), $|w_i|\le1-d_a<1-d$, and the points of $\{|w|=1-d\}$ have distance at least $d_a-d\ge\rho$ from the zeros. So $\mu_d\ge m(\rho)$, the displayed condition implies (iii), and Steps 1--3 of the proof of Theorem~\ref{thm:quant} apply: $Q_n(\sqrt n\,\cdot)$ has $k$ zeros in $\Omega_{\rm out}$, and its other $n-k$ zeros are real, simple and lie in $(-b,b)\subset(-2-\eps,2+\eps)$, as $2d^2\le\eps$. The disks $B_i=\{|w-w_i|\le\rho\}$ are disjoint, as $\rho\le\sigma/3$; they lie in $\{|w|<1-d\}$, as $|w_i|+\rho\le1-\frac23d_a$; they do not meet $\R$ if $w_i\notin\R$, as $|\im w_i|\ge\sigma/2$, and are symmetric about $\R$ if $w_i\in\R$. Since $|a|\ge m(\rho)>|\varphi_n-a|$ on $\{|w|\le1-d\}\setminus\bigcup_i\operatorname{int}B_i$ (Lemma~\ref{lem:ext-disk}), Rouch\'e's theorem gives exactly one zero of $\varphi_n$ in each $\operatorname{int}B_i$, real if and only if $w_i$ is real, and no other zero in $\{|w|<1-d\}$. On $B_i$, $|w|\ge r_a/2$, so $|\mathcal J'(w)|=|1-w^{-2}|\le\beta$ and $\mathcal J(\operatorname{int}B_i)\subset\{|z-\zeta_i|<\beta\rho\le\eps\}$. Finally, by \cite[\textup{(5.2)}]{Main}, $\dist(\zeta_i,[-2,2])\ge d_a^2$, and for $i\ne i'$ and nonreal $\zeta_i$ \[ |\zeta_i-\zeta_{i'}|=|w_i-w_{i'}|\,\frac{|1-w_iw_{i'}|}{|w_iw_{i'}|}\ge\sigma d_a,\qquad|\im\zeta_i|=|\im w_i|\,\frac{1-|w_i|^2}{|w_i|^2}\ge\frac{\sigma d_a}2 . \] These are all at least $\delta_0'\ge3\eps$, so the disks $\{|z-\zeta_i|<\eps\}$ are disjoint from each other and from $(-2-\eps,2+\eps)$, and from $\R$ when $\zeta_i\notin\R$. This gives (a) and (b) of \cite[Theorem~1.4]{Main}. \end{proof} \begin{remark}[The depth $m(\rho)$]\label{rem:depth} If $\rho\le\sigma/2$, at most one $w_i$ lies within $\sigma/2$ of a point $w$ that enters $m(\rho)$, so $|p_{\rm in}(w)|\ge G\rho(\sigma/2)^{k-1}$, and the proof of Lemma~\ref{lem:outer}(c) gives $m(\rho)\ge\mu2^{-k}\rho(\sigma/2)^{k-1}$. There is no lower bound polynomial in $\mu$, $\alpha$, $\sigma$, $\rho$: for $k\ge2$ and $0\E S+y)\le\exp(-y^2/(2\sum_jc_j^2))$ for $y\ge0$: indeed $\E e^{\lambda|\eta_j|}=2e^{\lambda^2/2}\Phi(\lambda)\le e^{\lambda\E|\eta_j|+\lambda^2/2}$ for $\lambda\ge0$, by the concavity of $\log\Phi$, where $\Phi$ is the standard normal distribution function. Here $\E S_1=\sqrt{2/\pi}\,(1-r_2)^{-2}\le52\eps^{-2}$, $\sum_jj^2r_2^{2j-2}=(1+r_2^2)(1-r_2^2)^{-3}\le2(1-r_2)^{-3}$, $\E S_2=2\sqrt{2/\pi}\,(1-r_2)^{-3}\le818\eps^{-3}$ and $\sum_jj^2(j-1)^2r_2^{2j-4}=4(1+4r_2^2+r_2^4)(1-r_2^2)^{-5}\le24(1-r_2)^{-5}$. Hence, for $0<\delta\le e^{-1}$, \begin{equation}\label{eq:Ytail} \PP\bigl(Y_i>y_i(\delta)\bigr)\le\delta,\qquad y_1(\delta)=\frac{52+46\,\ell_\delta^{1/2}}{\eps^2},\qquad y_2(\delta)=\frac{818+1255\,\ell_\delta^{1/2}}{\eps^3}. \end{equation} \begin{lemma}[Coefficients]\label{lem:coeff} Let $\lambda_\eps=(1+\eps)^{1/4}$ and $C_q=\sup_j|\eta_j|\lambda_\eps^{-j}$. Then $a^{(t)}\in\mathcal S(C_q,\lambda_\eps)$ for every $t\ge1+\eps$, $\PP(C_q>x)\le(7/\eps)e^{-x^2/2}$ for $x\ge1$, and $\PP(|\eta_0|<\delta)\le\delta$ for $\delta>0$. \end{lemma} \begin{proof} For $t\ge1+\eps$, $|\eta_j|t^{-j/2}\le|\eta_j|\lambda_\eps^{-2j}\le C_q\lambda_\eps^{-j}$. For $x\ge1$, $\PP(|\eta_j|>x\lambda_\eps^j)\le e^{-x^2\lambda_\eps^{2j}/2}\le e^{-x^2/2}\lambda_\eps^{-j}$, since $\lambda_\eps^{2j}-1\ge2j\log\lambda_\eps$. Summing over $j$ gives $e^{-x^2/2}/(1-\lambda_\eps^{-1})$, and $1-\lambda_\eps^{-1}\ge(1-2^{-1/4})\eps\ge\eps/7$ by concavity in $\eps$. \end{proof} \begin{lemma}[Small values on a circle]\label{lem:circle-small} For $00$, and \[ \det\operatorname{Cov}\bigl(\re q(u),\im q(u)\bigr)=\frac14\Bigl(\frac1{(1-r^2)^2}-\frac1{|1-u^2|^2}\Bigr)=\frac{r^2\sin^2\theta}{(1-r^2)^2|1-u^2|^2}. \] So $\PP(|q(u)|<2\delta)$ is at most $4\delta/\sqrt\pi$, and also at most $4\pi\delta^2(1-r^2)|1-u^2|/(2\pi r|\sin\theta|)\le4\delta^2/(r|\sin\theta|)$. We use the first bound at the at most $N\delta+4$ points $\theta_k$ with $|\sin\theta_k|<\delta$, and the second at the others, for which $\sum_k|\sin\theta_k|^{-1}\le4/\delta+(2N/\pi)\log(2/\delta)$. With $N\delta\le2\pi y+1$ and $r\le1$, \[ \PP\Bigl(\min_{|u|=r}|q|<\delta\Bigr)\le\delta+\frac{4\delta(2\pi y+5)}{\sqrt\pi}+\frac{16\delta}r+\frac{8\delta(2\pi y+1)}{\pi r}\log\frac2\delta\le\delta\Bigl[1+\frac{(2\pi y+5)(6+3\log\frac2\delta)}r\Bigr], \] and $2\pi y+5\le621\,\ell_\delta^{1/2}\eps^{-2}$ and $6+3\log(2/\delta)\le11.1\,\ell_\delta$ by \eqref{eq:Ytail} and $\ell_\delta\ge1$. \end{proof} \begin{lemma}[Near-double real zeros]\label{lem:nu-small} Let $t\ge1+\eps$ and $01$ and $\beta>0$, with $m_0=4+\max\bigl(e^2L_r^2\Upsilon^{4\beta}/\pi,\ 8e^2A_r^2\Upsilon^{8\beta}\bigr)$, \[ \PP(\mathcal Z_r\ge m)\le2\,\Upsilon^{-\beta m}\quad(m\ge m_0),\qquad\PP\bigl(\Upsilon^{\mathcal Z_r}>s\bigr)\le C(r,\Upsilon,\beta)\,s^{-\beta}\quad(s\ge1), \] with $C(r,\Upsilon,\beta)=\max(2,\Upsilon^{\beta m_0})$. \end{lemma} \begin{proof} Almost surely the zeros of $q$ in $\D\setminus\{0\}$ are simple (\cite[Lemma~7.4]{Main}), $q(0)\ne0$, and $q$ has no zero on $\{|u|=r\}$ (\cite[Lemma~7.5]{Main}(a)). So $\mathcal Z_r=N_\R+2k_r$, where $N_\R$ is the number of real zeros in $[-r,r]$ and $k_r$ the number of zeros in $\D_r^+$; and $\mathcal Z_r\ge m$ implies $N_\R\ge a$ or $k_r\ge b$. We use $\PP(N\ge j)\le\E\binom Nj$ and $j!\ge(j/e)^j$. \emph{Nonreal zeros.} Step~1 of the proof of \cite[Proposition~7.11]{Main} gives $\E\binom{k_r}b\le(\pi r^2C_r)^b(2b)^{b/2}/b!$ with $C_r=2r/(\pi(1-r^2)^3)$, and $\pi r^2C_r=A_r$. So $\PP(k_r\ge b)\le(\sqrt2\,eA_r/\sqrt b)^b$. \emph{Real zeros.} By \cite[Section~6.1, proof of Theorem~2.1]{MS}, the $a$-point correlation function of the real zeros of $q$ is \[ \rho_a(x_1,\dots,x_a)=(2\pi)^{-a/2}(\det\Sigma)^{1/2}\,\E\bigl|q(x_1)\cdots q(x_a)\bigr|,\qquad\Sigma=\Bigl(\frac1{1-x_ix_j}\Bigr)_{i,j\le a}, \] and by Cauchy's determinant formula $\det\Sigma=\prod_i(1-x_i^2)^{-1}\prod_{is)\le\PP(\mathcal Z_r\ge m)\le2\Upsilon^{-\beta m}\le2s^{-\beta}$. For $1\le s<\Upsilon^{m_0}$ the bound $C(r,\Upsilon,\beta)s^{-\beta}\ge1$ is trivial. \end{proof} So $\mathcal Z_r$, whose mean is of order $(1-r)^{-1}$, exceeds $m$ with probability decaying like $m^{-cm}$, and $\Upsilon^{\mathcal Z_r}$ has finite moments of every order. The constants are crude: for $r=0.8$ the first bound is below $1$ only for $m$ of order $3\cdot10^4$, while numerically $\PP(\mathcal Z_{0.8}\ge m)$ falls from $0.32$ at $m=3$ to $5.5\cdot10^{-4}$ at $m=8$. \begin{theorem}[A rate of convergence]\label{thm:rate} Let $q$ be the Gaussian power series and $0<\eps\le1$. For every $\alpha\in(0,\frac3{50})$ there is $C_{\eps,\alpha}<\infty$ such that for all $n\ge2$ \[ \sup_{t\ge1+\eps}\bigl|\PP(\Lambda_n^{\rm std}\le t)-\PP(T_*\le t)\bigr|\le C_{\eps,\alpha}\,n^{-\alpha}. \] For every $\theta\in(0,\frac1{10})$ there is $C_{\eps,\theta}<\infty$ such that for every $t\ge1+\eps$ and $n\ge2$ both $\bigl|\PP(G_{n,t}\text{ has only real zeros})-\PP(T_*\le t)\bigr|$ and the total variation distance between the laws of the number of nonreal zeros of $G_{n,t}$ and of the number of nonreal zeros of $q$ in $\{|w|0$. Let $C$ and $n_{\eps,\theta}$ denote constants depending only on $\eps$ and $\theta$. For $n\ge3$ put \begin{gather*} \delta=n^{-\theta},\qquad x_n=\max\bigl(1,\sqrt{2\theta\log n}\bigr),\qquad Z=\frac{x_n}{\delta(1-\lambda_\eps^{-1/2})},\\ W=Zn^{\theta'},\qquad c_1=\frac{\lambda_\eps^{1/2}(1-\vartheta_0)^2}6,\qquad\bar d_n=\frac{c_1W^{-2}}{\sqrt2}, \end{gather*} and $k_n=\lfloor\theta'\log n/\log\Upsilon_\eps\rfloor$. Let $\mathcal G_n(t)$ be the event that $C_q\le x_n$, $|\eta_0|\ge\delta$, the parameters of Definition~\ref{def:condition} for $a^{(t)}$ satisfy $\mu,\mu_{\bar d_n},\nu_{\bar d_n}\ge\delta$, $a^{(t)}$ has no zero in $\{1-\bar d_n\le|w|\le1\}$, and $\mathcal Z_{r_1}\le k_n$. \emph{On $\mathcal G_n(t)$ the count is exact.} The number $k$ of zeros of $a^{(t)}$ in $\D$ is the number of zeros of $q$ in $\{|u|x_n)+\PP(|\eta_0|<\delta)\le(7/\eps+1)\delta$. The parameters $\mu$ and $\mu_{\bar d_n}$ are minima of $|q|$ on circles with radii in $[T^{-1/2}/2,r_1]$, so Lemma~\ref{lem:circle-small} bounds their contribution by $21000\,\eps^{-2}T^{1/2}\delta\ell_\delta^{3/2}$; Lemma~\ref{lem:nu-small} bounds that of $\nu_{\bar d_n}$ by $(4000\eps^{-3}\ell_\delta^{1/2}+4T)\delta$, and Lemma~\ref{lem:annulus} that of the annulus by $5(1+\eps)^3\eps^{-3}\bar d_n\le C\delta^2$, as $\bar d_n\le CW^{-2}\le CZ^{-2}\le C\delta^2$. Finally, Lemma~\ref{lem:overcrowd} with $r=r_1$, $\Upsilon=\Upsilon_\eps$ and $\beta=\theta/\theta'$ gives $\PP(\mathcal Z_{r_1}\ge m)\le2\Upsilon_\eps^{-\beta m}$ for $m\ge m_0(r_1,\Upsilon_\eps,\beta)$; with $m=k_n+1>\theta'\log n/\log\Upsilon_\eps$ this gives $\PP(\mathcal Z_{r_1}>k_n)\le2n^{-\theta}$ as soon as $k_n+1\ge m_0$, that is, for $\log n\ge m_0\log\Upsilon_\eps/\theta'$. Enlarging $C$ to cover smaller $n$, we get for all $n\ge2$ \begin{equation}\label{eq:rate-T} \sup_{1+\eps\le t\le T}\Bigl(\bigl|\PP(\Lambda_n\le t)-\PP(T_*\le t)\bigr|+\PP\bigl(N_n(t)\ne N_t\bigr)\Bigr)\le C\,T\,n^{-\theta}(\log n)^{3/2}. \end{equation} \emph{Large times.} By the proof of \cite[Lemma~7.9]{Main}, with $r=T^{-1/2}$ and $\arctan\varrho\le\varrho$, $\PP(T_*>T)\le2T^{-3/2}$ for $T\ge2$. For $t\ge T$, by monotonicity, $|\PP(\Lambda_n\le t)-\PP(T_*\le t)|\le|\PP(\Lambda_n\le T)-\PP(T_*\le T)|+2\PP(T_*>T)$. With $T=\max(2,n^{2\theta/5})$ and \eqref{eq:rate-T} the supremum over $t\ge1+\eps$ is at most $Cn^{-3\theta/5}(\log n)^{3/2}$. Given $\alpha<\frac3{50}$, choose $\theta\in(\frac53\alpha,\frac1{10})$; then $n^{-3\theta/5}(\log n)^{3/2}\le Cn^{-\alpha}$. The claims at a fixed $t$ are \eqref{eq:rate-T} with $T=t$, applied with a slightly larger $\theta<\frac1{10}$ to absorb the logarithm, since $P_{n,t}$ has only real zeros exactly when $N_n(t)=0$, and the total variation distance between the laws of $N_n(t)$ and $N_t$ is at most $\PP(N_n(t)\ne N_t)$. \end{proof} \begin{remark}[The exponent]\label{rem:rate} The exponents do not depend on $\eps$: the only place where the argument could lose a factor depending on $\eps$ is the factor $\Upsilon^k$, and its bound $\Upsilon^k\le(\overline M/\alpha)^E$ of Lemma~\ref{lem:outer}(b), with $E\asymp\eps^{-1}\log(1/\eps)$, would give an exponent of order $\eps/\log(1/\eps)$; Lemma~\ref{lem:overcrowd} replaces it. The constants are not explicit and not polynomial in $1/\eps$: the proof needs $\log n\ge m_0\log\Upsilon_\eps/\theta'$, so $\log C_{\eps,\alpha}=\eps^{-O(1/(3/50-\alpha))}$. The limits $\frac1{10}$ and $\frac3{50}$ come from the degree, about $9$ in $1/\mu$, of the explicit onset (Theorem~\ref{thm:quant} and the paragraph after it). Since $\PP(\mu1$. In a Monte Carlo simulation (not reported here; $n\le1024$, $t\in\{1.25,1.5,2,3\}$, $24\,000$ to $48\,000$ samples for each pair $(n,t)$), the probability that exactly one of the events $\{P_{n,t}\text{ has only real zeros}\}$ and $\{T_*\le t\}$ occurs decays with fitted exponents $0.68$ to $0.82$, and $\PP(\Lambda_n\le t)-\PP(T_*\le t)$ is positive at every computed point and decays with fitted exponents $0.4$ to $0.7$. \end{remark} The distortion $\mathfrak d(a^{(t)})$ of the outer factor (Lemma~\ref{lem:outer}(c)) also has a small-ball bound whose exponent does not depend on $\eps$. We need the circle minimum uniformly in the radius. \begin{lemma}[Small values on a circle, uniformly in the radius]\label{lem:circle-uniform} For $t\ge1+\eps$ and $0y_1(x))\le x$, and $\PP(|\eta_0|\delta^{-\eta/4}\}\cup\{\mu<\delta^{1-3\eta/4}\}\cup\{(\Upsilon')^{\mathcal Z_{r_1}}>\delta^{-\eta/2}\} \] we have $\mathfrak d\ge\delta^{1-3\eta/4}\delta^{\eta/4}\delta^{\eta/2}=\delta$. The first event has probability at most $C_\eps\exp(-c_\eps\delta^{-\eta/2})\le C_{\eps,\eta}\delta$ by Lemma~\ref{lem:coeff}, as $\overline M=C_q/(1-\lambda_\eps^{-1/2})$. The third has probability at most $C_{\eps,\eta}\delta$ by Lemma~\ref{lem:overcrowd} with $\Upsilon=\Upsilon'$ and exponent $2/\eta$. For the second, $\mu=\min_{|u|=t^{-1/2}}|q(u)|$: for (a), Lemma~\ref{lem:circle-small} with $r=t^{-1/2}\ge T^{-1/2}$ gives at most $C_\eps T^{1/2}\delta^{1-3\eta/4}\log^{3/2}(1/\delta)$; for (b), Lemma~\ref{lem:circle-uniform} gives at most $C_\eps\delta^{(1-3\eta/4)/2}\log^{3/2}(1/\delta)$. Take $\eta$ small. \end{proof} The exponent in (b) is limited by small radii, where $q(u)\approx\eta_0+\eta_1u$ on $|u|=r$ and the circle minimum is a small-ball problem for $|\eta_0|-r|\eta_1|$; we expect exponents up to $1$ uniformly in $t$, but have not proved it. \begin{remark}[Another deformation]\label{rem:homotopy} Let $a\in\mathcal S(C_0,R_0)$ with $\alpha,\mu>0$, let $p_{\rm in},p_{\rm out}$ be as in Lemma~\ref{lem:outer}, let $\rho_{\rm out}$ be the smallest modulus of a zero of $a$ outside $\overline\D$ (or $R_0$ if it is smaller), and let $a_s^{\rm h}=p_{\rm in}\,p_{\rm out}^s=p_{\rm in}e^{s\log p_{\rm out}}$, with the branch of $\log p_{\rm out}$ that vanishes at $0$; it is analytic in $\{|w|<\rho_{\rm out}\}$ and has real coefficients. On $|w|=1$, $|a_s^{\rm h}|=|p_{\rm in}|^{1-s}|a|^s$, so $\min_{s\in[0,1],|w|=1}|a_s^{\rm h}|=\min(\mu,m_{\rm in})$ with $m_{\rm in}=\min_{|w|=1}|p_{\rm in}|$. Any continuous family from $p_{\rm in}$ to $a$ has circle minimum at most $\min(\mu,m_{\rm in})$ (look at its end points), so this family is optimal in that respect; and $m_{\rm in}\ge\mathfrak d(a)\mu$, since $|p_{\rm out}(0)|=1$. But $m_{\rm in}$ is a product over the zeros of $a$ in $\D$, and its lower bounds cost a factor exponential in $k$; and if a zero $w'$ of $a$ with $|w'|=\rho_{\rm out}0$, that is, $q$ has no multiple zero in $[-r_c,r_c]$; (G3) gives this only in $(-R_*,R_*)$. Let $\delta_1=\min\{\eps/4,\ gT_*^{3/2}/6\}$, let $\mathfrak A$ be the closed annulus $\{(T_*+2\delta_1)^{-1/2}\le|u|\le(T_*-2\delta_1)^{-1/2}\}$, as in \cite[Section~9.3]{Main} with this $\delta_1$, and for $\rho>0$ let $m_{\mathfrak A}(\rho)$ be the minimum of $|q|$ over the points of $\mathfrak A$ at distance at least $\rho$ from $u_*$ and $\bar u_*$. Put \[ c_F=\frac{|\gamma_*|}{8T_*},\qquad\overline B_2=\frac{4000\max(1,C_q)}{\eps^3},\qquad r_1=\min\Bigl\{\frac{c_F}{\overline B_2},\ \delta_1,\ \frac{m_*}2\Bigr\},\qquad\rho_V=\frac{r_1s_*}{6T_*^{3/2}} . \] Finally let $k_-$ be the number of zeros of $q$ in $\{|u|<(T_*-\delta_1)^{-1/2}\}$, $\mathcal U=\Upsilon_c^{k_-}$, and \[ \begin{aligned} \Pi={}&\max(1,|\gamma_*|^{-1})\cdot s_*^{-1}\cdot g^{-1}\cdot\max\bigl(1,m_{\mathfrak A}(\rho_V)^{-1}\bigr)\\ &\cdot\max\bigl(1,(\nu^{\rm glob})^{-1}\bigr)\cdot\max(1,|\eta_0|^{-1})\cdot\max(1,C_q)\cdot T_*\cdot\mathcal U . \end{aligned} \] Each of the nine factors of $\Pi$ is at least $1$; we call them $P_\gamma,P_s,P_g,P_{\mathfrak A},P_\nu,P_\alpha,P_C,T_*,\mathcal U$, in this order. Below, $c$ denotes absolute constants and $c_\eps$ numbers with $c_\eps\le c\,\eps^{-c}$, not the same at each occurrence. \begin{theorem}[Explicit constants]\label{thm:collision-explicit} In this setting the conclusions (a)--(d) of \cite[Theorem~9.1]{Main} hold with $\delta=\delta_1$ and with $K$, $n_0$ such that \[ \delta^{-1}\le6\eps^{-1}\,\Pi,\qquad K\le c_\eps\,\Pi^{5/2},\qquad n_0\le c_\eps\,\Pi^{10}. \] Moreover $\mathcal U\le(\overline M_c/|\eta_0|)^{E_c}$, and $\mathcal U\le\Upsilon_c^{\mathcal Z}$, where $\mathcal Z$ is the number of zeros of $q$ in $\{|u|\le r_c\}$, counted with multiplicity. \end{theorem} The factor $\mathcal U=\Upsilon^{k}$ of Theorem~\ref{thm:quant} is kept as a parameter, because for random $q$ it can be bounded either by Jensen's formula or through the number of zeros in a disk (Theorem~\ref{thm:Dtail}). \begin{lemma}\label{lem:logs} \begin{enumerate} \item[(a)] If $X\ge1$ and $n\ge16X(1+\log X)^2$, then $n\ge X(1+\log n)^2$. In particular this holds for $n\ge64X^2$. \item[(b)] If $P_1,\dots,P_k$ are distinct factors of $\Pi$ and $a_1,\dots,a_k\ge0$, then $\prod_iP_i^{a_i}\le\Pi^{\max_ia_i}$. \end{enumerate} \end{lemma} \begin{proof} (a) The function $n\mapsto n/(1+\log n)^2$ is increasing for $n\ge e$. At $n_1=16X(1+\log X)^2\ge16$, using $\log(1+y)\le y$, we have $1+\log n_1\le1+\log16+3\log X\le4(1+\log X)$, so $n_1/(1+\log n_1)^2\ge X$. Moreover $(1+\log X)^2\le4X$ for $X\ge1$. (b) Each $P_i\ge1$, so $P_i^{a_i}\le P_i^{\max a}$, and a product of distinct factors of $\Pi$ is at most $\Pi$. \end{proof} By Lemma~\ref{lem:logs}(b), a bound by $c_\eps$ times a monomial in the factors of $\Pi$ is a bound by $c_\eps\Pi^e$, where $e$ is the largest exponent of a single factor. \begin{proof}[Proof of Theorem~\ref{thm:collision-explicit}] We go through the proof of \cite[Theorem~9.1]{Main}, which runs from \cite[Section~9.2]{Main} to \cite[Section~9.5]{Main}, and bound each constant. The elementary numerical inequalities used below were also checked on grids by computer. \emph{Step 1: $\delta_1$ and the annulus.} As $g\le1\le T_*$ and $\eps\le1$, $\delta_1\ge\min(\eps/4,g/6)\ge\eps g/6$, so $\delta_1^{-1}\le6\eps^{-1}P_g$. Also $\delta_1\le\eps/40$ and no zero on the unit circle, gives $\mathcal U=\Upsilon_c^{k_-}\le(\overline M_c/|\eta_0|)^{E_c}$. And $(T_*-\delta_1)^{-1/2}0\}$ and $\D_r^+=\D_r\cap\D^+$. The zeros of $q$ have the one-point intensities \[ \rho_r(x)=\frac1{\pi(1-x^2)}\quad(x\in(-1,1)),\qquad\rho_c(z)=\frac{2\im z}{\pi\,|1-z^2|\,(1-|z|^2)^2}=\frac{b_z}{\pi(1-|z|^2)}\quad(z\in\D^+), \] for the real zeros and for the zeros in $\D^+$, with respect to Lebesgue measure \cite[Theorems~2.1 and~2.7]{MS}, \cite[Lemma~7.9]{Main}, where $b_z=2\im z/((1-|z|^2)|1-z^2|)$. Let $\psi(a,b)=|a-b|/|1-a\bar b|$ be the pseudo-hyperbolic distance on $\D$. For $a\in\D$ let $B_a(w)=(w-a)/(1-aw)$ if $a$ is real, and $B_a(w)=(w-a)(w-\bar a)/((1-\bar aw)(1-aw))$ if $a$ is not real. These are Blaschke products with real Taylor coefficients, and $|B_x(y)|=\psi(x,y)$, $|B_z(x)|=\psi(x,z)^2$, $|B_z(w)|=\psi(z,w)\psi(z,\bar w)$ and $|B_z'(z)|=b_z$. Rice formulas are used in the form of \cite[Theorem~6.2]{AW}% , for Gaussian fields $Z:U\to\R^d$ on open sets $U\subset\R^d$ with analytic paths, a nondegenerate law of $Z(t)$ for each $t\in U$, and almost surely no zero $t$ with $\det Z'(t)=0$; the last condition holds below because the zeros of $q$ in $\D\setminus\{0\}$ are simple \cite[Lemma~7.4]{Main}. \begin{lemma}[Regression]\label{lem:regression} Let $a_1,\dots,a_m\in\D$ be distinct points whose set is closed under complex conjugation, and let $\mathfrak B$ be the product of the factors $B_x$ over the real points $x$ and $B_z$ over the pairs $\{z,\bar z\}$. Then the conditional law of $q$ given $q(a_1)=\dots=q(a_m)=0$ (the Gaussian regression) is the law of $\mathfrak B\check q$, where $\check q$ has the law of $q$. In particular, conditionally on these values, $(q'(a_i))_{i\le m}$ has the law of $(\mathfrak B'(a_i)\check q(a_i))_{i\le m}$. The density at $0$ of the real vector formed by the values $q(a_i)$ (real and imaginary parts at the nonreal points, one point of each pair) is the product, over the points taken in any order, of $(2\pi)^{-1/2}(1-x^2)^{1/2}/|\mathfrak B_<(x)|$ for a real point $x$ and $1/(\pi b_z|\mathfrak B_<(z)|^2)$ for a pair $z$, where $\mathfrak B_<$ is the product of the factors of the points taken before. \end{lemma} \begin{proof} For one real point $x$ this is \cite[Lemma~4.2]{MS}: the covariance $1/(1-uv)$ of $q$ minus its regression on $q(x)$ is $B_x(u)B_x(v)/(1-uv)$. For one conjugate pair it is \cite[Lemma~7.2]{MS}, with the same identity for the covariance and the pseudo-covariance. For several points we condition on one point or pair at a time: conditioning a Gaussian vector on a block of coordinates and then on the next block is conditioning on both. Once $q$ has been replaced by $\mathfrak B_A\hat q$, the remaining conditions $q(a)=0$ are equivalent to $\hat q(a)=0$, because $\mathfrak B_A(a)\ne0$; this proves the first claim by induction, also when real points and pairs are mixed. The derivative identity follows by differentiating $\mathfrak B\check q$ at $a_i$, where $\mathfrak B$ vanishes. For the densities, the density at $0$ of a Gaussian vector is the product of the conditional densities at $0$ of its blocks, and by the first claim the conditional law of $q(a)$ given the earlier values is that of $\mathfrak B_<(a)\check q(a)$. Finally $\operatorname{Var}q(x)=(1-x^2)^{-1}$, and $(\det\operatorname{Cov}(\re q(z),\im q(z)))^{1/2}=b_z/2$ (proof of Lemma~\ref{lem:circle-small}). \end{proof} \begin{proposition}[Two-point functions]\label{prop:twopoint} Let $\rho_{rr}$, $\rho_{rc}$ and $\rho_{cc}$ be the two-point correlation functions, with respect to Lebesgue measure, of pairs of real zeros, of a real zero and a zero in $\D^+$, and of pairs of zeros in $\D^+$. For distinct $x,y\in(-1,1)$ and $z,w\in\D^+$, \begin{equation}\label{eq:twopoint} \begin{aligned} \rho_{rr}(x,y)&=\psi\bigl(\psi+\sqrt{1-\psi^2}\arcsin\sqrt{1-\psi^2}\bigr)\,\rho_r(x)\rho_r(y),&&\psi=\psi(x,y),\\ \rho_{rc}(x,z)&=\psi^2(2-\psi^2)\,\rho_r(x)\rho_c(z),&&\psi=\psi(x,z),\\ \rho_{cc}(z,w)&=ab\,(3-a-b)\,\rho_c(z)\rho_c(w),&&a=\psi(z,w)^2,\ b=\psi(z,\bar w)^2 . \end{aligned} \end{equation} In particular $\rho_{rr}\le\rho_r\rho_r$, $\rho_{rc}\le\rho_r\rho_c$ and $\rho_{cc}\le\rho_c\rho_c$. \end{proposition} \begin{proof} \emph{The mixed pair.} The field $Z(x,z)=(q(x),\re q(z),\im q(z))$ on $U=(-1,1)\times\D^+$ has nondegenerate values, since by Lemma~\ref{lem:regression} the density of $Z(x,z)$ at $0$ is finite, and its zeros in $I\times B$ are the pairs formed by a real zero in $I$ and a zero in $B$. Its Jacobian determinant is $q'(x)|q'(z)|^2$, by the Cauchy--Riemann equations. So \cite[Theorem~6.2]{AW} gives $\E[N_\R(I)N_+(B)]=\int_I\int_B\rho_{rc}$ with \[ \rho_{rc}(x,z)=p_{Z(x,z)}(0)\,\E\bigl[|q'(x)|\,|q'(z)|^2\ \big|\ q(x)=0,\ q(z)=0\bigr]. \] By Lemma~\ref{lem:regression} with $\mathfrak B=B_xB_z$: $p_{Z(x,z)}(0)=(2\pi)^{-1/2}(1-x^2)^{1/2}(\pi b_z\psi^2)^{-1}$, since $|B_x(z)|^2=\psi^2$; $q'(x)=B_x'(x)B_z(x)\check q(x)$ with $|B_x'(x)|=(1-x^2)^{-1}$ and $|B_z(x)|=\psi^2$; and $|q'(z)|=\psi b_z|\check q(z)|$. Hence $\rho_{rc}=\psi^2(1-x^2)^{-1/2}b_z(2\pi)^{-1/2}\pi^{-1}M$ with $M=\E[|X|\,|Y|^2]$, $X=\check q(x)$ and $Y=\check q(z)$, viewed in $\R^2$. Write $Y=\mathbf aX+W$ with $W$ independent of $X$; the cross term vanishes because $\E[|X|X]=0$, and $\E|X|^3=2\sqrt{2/\pi}\,\sigma_X^3$. So \[ M=\E|X|\Bigl(\E|Y|^2+\frac{|\E XY|^2}{\E X^2}\Bigr)=\sqrt{\frac2\pi}\,\frac{1+(1-\psi^2)}{(1-x^2)^{1/2}(1-|z|^2)}, \] since $\E X^2=(1-x^2)^{-1}$, $\E|Y|^2=(1-|z|^2)^{-1}$, $\E XY=(1-xz)^{-1}$, $|1-xz|=|1-x\bar z|$ and $(1-x^2)(1-|z|^2)/|1-x\bar z|^2=1-\psi^2$. With $\rho_r(x)\rho_c(z)=b_z/(\pi^2(1-x^2)(1-|z|^2))$ this is the second line of \eqref{eq:twopoint}. \emph{Two nonreal zeros.} Use the field $(\re q(z),\im q(z),\re q(w),\im q(w))$ on $\{(z,w)\in(\D^+)^2:z\ne w\}$ and $\mathfrak B=B_zB_w$, with $|B_z(w)|=\psi(z,w)\psi(z,\bar w)=\sqrt{ab}$. The density at $0$ is $(\pi b_z)^{-1}(\pi b_wab)^{-1}$, and $|q'(z)|^2|q'(w)|^2=b_z^2b_w^2a^2b^2|\check q(z)|^2|\check q(w)|^2$. By Isserlis' formula $\E|X|^2|Y|^2=\E|X|^2\E|Y|^2+|\E XY|^2+|\E X\bar Y|^2$ for complex Gaussian $X,Y$ with real covariance structure, and $(1-|z|^2)(1-|w|^2)/|1-z\bar w|^2=1-a$, $(1-|z|^2)(1-|w|^2)/|1-zw|^2=1-b$; so the ratio of the moment to $\E|X|^2\E|Y|^2$ is $3-a-b$, which gives the third line. \emph{Two real zeros.} Use the field $(q(x),q(y))$ on $\{x\ne y\}$, or \cite[Lemma~6.1]{MS}. The density at $0$ is $(2\pi)^{-1}(1-x^2)^{1/2}(1-y^2)^{1/2}/\psi$, the Jacobian factor is $\psi^2(1-x^2)^{-1}(1-y^2)^{-1}|\check q(x)\check q(y)|$, and $\E|XY|=\frac2\pi\sigma_X\sigma_Y(\sqrt{1-c^2}+c\arcsin c)$ for the correlation $c=\sqrt{1-\psi^2}$ of $X=\check q(x)$, $Y=\check q(y)$. \emph{The bounds.} $t(2-t)\le1$ for all real $t$. On $[0,1]^2$, $\partial_a[ab(3-a-b)]=b(3-2a-b)\ge0$, and symmetrically in $b$, so $ab(3-a-b)\le1$. With $\psi=\cos\varphi$, $\varphi\in[0,\frac\pi2]$, the first factor is $\cos^2\varphi+\varphi\sin\varphi\cos\varphi\le1$, because $\varphi\le\tan\varphi$. \end{proof} The bound $\rho_{rr}\le\rho_r\rho_r$ is \cite[Corollary~2.2]{MS}, and the closed form for $\rho_{rr}$ is equivalent to \cite[(2.6)]{MS}; the closed form for $\rho_{cc}$ is the case $n=2$ of \cite[Theorem~2.7]{MS}, and the bound $\rho_{cc}\le\rho_c\rho_c$ is stated after it. Mixed correlations are not treated in \cite{MS}; Forrester \cite[(30)]{Fo} gives all mixed correlations as a Pfaffian, and $\rho_{rc}$ can also be read off its entries \cite[(49)--(52)]{Fo}. All three closed forms agree with a direct numerical evaluation of the Kac--Rice formula. Only the three bounds are used below. \begin{theorem}[The modulus gap]\label{thm:gap} For $0<\eps\le1$ and $s>0$, \[ \PP\bigl(T_*\ge1+\eps,\ g(1-r_\eps)/2$, the right side exceeds $107\eps^{-3}/8>1$; so let $s\le(1-r_\eps)/2\le0.15$. On the event, $u_*\in\D^+$ is a zero with $|u_*|=R_*\le r_\eps$, and some zero $u'\notin\{u_*,\bar u_*\}$ of $q$ in $\D$ has $\bigl||u'|-R_*\bigr|2/(3\pi)$ there. For $0.71$ for $r\ge0.8$. As $\{R_*<\varrho\}=\{T_*>\varrho^{-2}\}$, \cite[Lemma~7.9]{Main} gives \begin{equation}\label{eq:Rsmall} \PP(R_*<\varrho)\le2\varrho^3\quad(\varrho>0),\qquad\PP(T_*>s)\le2s^{-3/2}\quad(s>1). \end{equation} \cite[Lemma~9.12]{Main} is the qualitative form of Theorem~\ref{thm:gap}. Numerically, $\PP(T_*\ge1.5,\ gy_3(\ell)\bigr)\le e^{-\ell}\quad(\ell\ge0),\qquad y_3(\ell)=\frac{2758+3456\,\ell^{1/2}}{\eps^3}. \end{equation} \begin{lemma}[Small balls for $q'(u_*)$, $\gamma_*$ and $\sin\theta_*$]\label{lem:collision-smallball} For $00$ and $00$ let \[ \mathfrak m(\rho)=\inf\bigl\{|q(u)|:\ |u|\le r_c,\ q\text{ has no zero in }\{|v-u|<\rho\}\bigr\}, \] which is nondecreasing in $\rho$. In the setting of Section~\ref{subsec:collision-explicit}, every point of $\mathfrak A$ has modulus within $g/2$ of $R_*$ (Step 1 of the proof of Theorem~\ref{thm:collision-explicit}), every zero of $q$ other than $u_*,\bar u_*$ has modulus at distance at least $g$ from $R_*$, and $\mathfrak A\subset\overline{\D}_{r_c}$. Hence \begin{equation}\label{eq:depth-annulus} m_{\mathfrak A}(\rho)\ge\mathfrak m\bigl(\min(\rho,g/2)\bigr)\qquad(\rho>0). \end{equation} \begin{lemma}[Depth away from the zeros]\label{lem:depth} For $0<\rho\le\eps/12$ and $0<\delta\le e^{-1}$, with $y=y_3(\log(1/\delta))$, \[ \PP\bigl(\mathfrak m(\rho)<\delta\bigr)\le\delta+\Bigl(\frac{16\delta^2}{\rho^2}+4y\delta\Bigr)\bigl(19+5.4\log(2.04y)\bigr). \] \end{lemma} \begin{proof} Let $\tau=\max(4\delta/\rho,\ 2\sqrt{y\delta})$. If $\tau>1$, the right side exceeds $1$. Otherwise work on $\{Y_2'\le y\}$, whose complement has probability at most $\delta$ by \eqref{eq:Y3tail}. Let $|u|\le r_c$, $|q(u)|<\delta$ and $|q'(u)|\ge\tau$, and put $s=2\delta/|q'(u)|\le\rho/2$. On $|v-u|=s$, $|q(v)-q'(u)(v-u)|\le\delta+ys^2/2\le\frac32\delta<2\delta=|q'(u)(v-u)|$, because $ys^2\le4y\delta^2/\tau^2\le\delta$. By Rouch\'e's theorem $q$ has a zero in $\{|v-u|s)\le c\,\eps^{-c}s^{-1}(1+\log s)$ for $s\ge1$, with an absolute constant $c$. \end{lemma} \begin{proof} Let $\rho_k=2^{-k-1}\eps/12$ and $\delta_k=(2\rho_k)^{3/2}/s\le e^{-1}$ for $k\ge0$. If $V_{\mathfrak m}>s$, then some $\rho\in(\rho_k,2\rho_k]$ has $\mathfrak m(\rho)<\rho^{3/2}/s$, and $\mathfrak m(\rho_k)\le\mathfrak m(\rho)<\delta_k$. By Lemma~\ref{lem:depth}, with $y_k=y_3(\log(1/\delta_k))$ and $\delta_k^2/\rho_k^2=8\rho_k/s^2$, \[ \PP\bigl(\mathfrak m(\rho_k)<\delta_k\bigr)\le\delta_k+\bigl(128\rho_ks^{-2}+4y_k\delta_k\bigr)\bigl(19+5.4\log(2.04y_k)\bigr). \] Here $\delta_k\le2^{-3k/2}/s$, $\sum_k\rho_k\le\eps/12$, and $\log(1/\delta_k)\le\log s+2(k+1)+\frac32\log(12/\eps)$, so $y_k\le c\,\eps^{-3}(1+\log s+k+\log(1/\eps))^{1/2}$. Summing over $k$ gives the claim. \end{proof} \subsection{The tail of \texorpdfstring{$D$}{D}}\label{subsec:Dtail} Let \[ D=\sup_{m\ge1}m\bigl(|\Lambda_m-T_*|+|z_m-2R_*^{-1}\cos\theta_*|\bigr), \] with $z_m=F_m(\eta_0,\dots,\eta_m)$, as in the proof of \cite[Corollary~9.16]{Main}. By \cite[Corollary~9.14]{Main}, $D<\infty$ almost surely. For the first degrees we use two crude bounds. \begin{lemma}[Small degrees]\label{lem:smalldeg} Let $\Sigma_3=\sum_{j\ge3}(j+1)2^{3-j}|\eta_j|$ and $\Lambda^{\sup}=\max\{30,\ (1.42|\eta_2|+45.5\,\Sigma_3)/|\eta_0|\}$. \begin{enumerate} \item[(a)] For every $m\ge2$, $\Lambda_m\le\Lambda^{\sup}$, and $\PP(\Lambda^{\sup}>x)\le365/x$ for $x\ge30$. \item[(b)] For $m\ge2$ and $t\ge0$, every zero $z$ of $P_{m,t}$ satisfies $|z|\le4(1+mt)^{1/2}H_m$, with $H_m=\max\bigl(1,\max_{j\le m}|\eta_j|/|\eta_0|\bigr)$. \item[(c)] $Y'=\sup_{j\ge0}|\eta_j|/(2\log(j+2))^{1/2}$ satisfies $\PP(Y'>y)\le2\cdot2^{-y^2}$ for $y\ge2$. \end{enumerate} \end{lemma} \begin{proof} (a) Let $t\ge30$ and $a_j=\eta_jt^{-j/2}$. In \cite[Lemma~8.2(c)]{Main} with $\kappa=1$ we have $c_\kappa=\coth1+\frac12<1.82$, so for $j\le\sqrt m$ (when $j\le m/2$) $\beta_{m,j}\le2e^{1.82}(j+1)e^j\le12.4(j+1)e^j$, and for $j>\sqrt m$, $\beta_{m,j}\le\sqrt m4e^2$, $e/\sqrt t\le\frac12$, and the left side of \cite[\textup{(8.3)}]{Main} is at most \[ \frac{\sqrt2\,|\eta_2|}t+12.4\,\frac{e^3}{t^{3/2}}\sum_{j\ge3}(j+1)2^{3-j}|\eta_j|\le\frac{1.42|\eta_2|+45.5\,\Sigma_3}t, \] which is less than $|\eta_0|$ for $t>\Lambda^{\sup}$. So $Q_m^{(t)}$, and hence $P_{m,t}$, has only real zeros for every $t>\Lambda^{\sup}$, and $\Lambda_m\le\Lambda^{\sup}$. Conditionally on $(\eta_j)_{j\ge1}$, $\PP(|\eta_0|y(2\log(j+2))^{1/2}\bigr)\le(j+2)^{-y^2}$, and $\sum_{k\ge2}k^{-y^2}\le2\cdot2^{-y^2}$ for $y\ge2$. \end{proof} \begin{theorem}[The tail of $D$]\label{thm:Dtail} Let $q$ be the Gaussian power series, $0<\eps\le1$, and $E_c=\log\Upsilon_c/\log\lambda_{\eps/2}^{1/2}$ as in Section~\ref{subsec:collision-explicit}. There is an absolute constant $c$ such that for all $x\ge1$: \begin{enumerate} \item[(a)] $\PP(D>x,\ T_*\ge1+\eps)\le c\,\eps^{-c}x^{-\beta_\eps}$, where $\beta_\eps=1/(1782E_c)\ge\eps/(35640\log(40/\eps))$; \item[(b)] $\PP(D>x,\ T_*\ge1+\eps)\le C_\eps\,x^{-1/900}$, where $\log C_\eps\le c\,\eps^{-c}$. \end{enumerate} \end{theorem} \begin{proof} Work on $\Omega=\{T_*\ge1+\eps\}$, intersected with the almost sure event on which (G1)--(G3) hold with $\rho_q=1$, $C_q<\infty$, $\nu^{\rm glob}>0$ (as $\eta_0\ne0$ and the zeros of $q$ in $\D\setminus\{0\}$ are simple, \cite[Lemma~7.4]{Main}), and the conclusions of \cite[Corollary~9.14]{Main} hold. Theorem~\ref{thm:collision-explicit} applies. \emph{A deterministic bound.} For $m\ge n_0$ the $m$-th term of $D$ is at most $2K$. For $m=1$ it is $T_*+2T_*^{1/2}$, as $\Lambda_1=0$ and $z_1=0$. For $2\le ms,\ \Omega)\le c\,\eps^{-c}s^{-b}(1+\log s)^2$ for $s\ge1$, with the following $b$: \begin{itemize} \item $P_\gamma$: $b=6/5$, by Lemma~\ref{lem:nearcritical} with $\delta=0$ and $y=y_3(2\log(1/\tau))$, and $\PP(R_*<\varrho)\le2\varrho^3$ (by \eqref{eq:Rsmall}) with $\varrho=\tau^{2/5}$ (or $b=3/2$ by Lemma~\ref{lem:collision-smallball}(b)); \item $P_s$: $b=2$ (Lemma~\ref{lem:collision-smallball}(c)); $P_g$: $b=1$ (Theorem~\ref{thm:gap}); $P_\alpha$: $b=1$; $T_*$: $b=3/2$, by \eqref{eq:Rsmall}; \item $P_C$ and $Y''$: every $b$, by Lemma~\ref{lem:coeff} with $\eps/2$ in place of $\eps$, and Lemma~\ref{lem:smalldeg}(c); $\Lambda^{\sup}$: $b=1$ (Lemma~\ref{lem:smalldeg}(a)); \item $\max(1,V_{\mathfrak m})$: $b=1$ (Lemma~\ref{lem:depth-scalefree}); $P_\nu$: $b=1$, because the proof of Lemma~\ref{lem:nu-small}, with $\eps/2$ in place of $\eps$ and $\delta'=\delta$, gives $\PP(\nu^{\rm glob}<\delta)\le(32000\,\eps^{-3}\log^{1/2}(1/\delta)+4)\delta$; \item $\mathcal U$: for (a), $\mathcal U\le(\overline M_c/|\eta_0|)^{E_c}$ gives $\PP(\mathcal U>s)\le\PP(\overline M_c>s^{1/(2E_c)})+\PP(|\eta_0|s)\le C(r_c,\Upsilon_c,1)\,s^{-1}$, where $\log C(r_c,\Upsilon_c,1)\le c\,\eps^{-15}$ because $\log((1+r_c)/(1-r_c))\le\log(12/\eps)$, $2r_c^3/(1-r_c^2)^3\le54\eps^{-3}$ and $\Upsilon_c\le40/\eps$. \end{itemize} By \eqref{eq:Dpoly}, $D>x$ forces some factor of $\Pi_3$ to exceed $(x/c_\eps)^{1/891}$. Summing the eleven bounds gives (a), the smallest exponent being $1/(2E_c)$ (note $E_c\ge72$) and $E_c\le20\eps^{-1}\log(40/\eps)$, and (b), all exponents being at least $1$ and the logarithms being absorbed by the margin $1/891-1/900$. \end{proof} \begin{corollary}\label{cor:Dtail} There are absolute constants $c>0$ and $C$ such that for $x\ge e^e$ \[ \PP(D>x)\le C\exp\Bigl(-c\Bigl(\frac{\log x}{\log\log x}\Bigr)^{1/2}\Bigr). \] \end{corollary} \begin{proof} By \cite[Theorem~7.13]{Main} there is an absolute $C_0$ with $\PP(T_*<1+\eps)\le C_0e^{-0.8/\eps}$ for $0<\eps\le1$. With Theorem~\ref{thm:Dtail}(a), \[ \PP(D>x)\le C_0e^{-0.8/\eps}+c\,\eps^{-c}\exp\Bigl(-\frac{\eps\log x}{35640\log(40/\eps)}\Bigr). \] Take $\eps=(\log\log x/\log x)^{1/2}$. For large $x$, $\log(40/\eps)\le\log\log x$, so the second term is at most $\exp\bigl(c\log\log x-(\log x/\log\log x)^{1/2}/35640\bigr)$, and the first is $C_0\exp(-0.8(\log x/\log\log x)^{1/2})$. Enlarge $C$ for the remaining $x$. \end{proof} \begin{remark}\label{rem:Dtail} (1) The tail of $D$ on $\{T_*\ge1+\eps\}$ rests on no unproved hypothesis: the small-ball bound for the modulus gap is Theorem~\ref{thm:gap}, and the weighted Rice formula is not needed. (2) In (a) the constant is polynomial in $1/\eps$ and the exponent is of order $\eps/\log(1/\eps)$; it comes from Jensen's bound for $\mathcal U$. In (b) the exponent does not depend on $\eps$, because the number of zeros of $q$ in a disk of radius $r_c$ has super-exponential tails (Lemma~\ref{lem:overcrowd}), but $\log C_\eps$ is of order $\eps^{-O(1)}$. (3) Whether $\PP(D>x)$ has a power bound without the restriction $T_*\ge1+\eps$ is Question~\ref{q:Dtail}. By \cite[Theorem~1.7]{Main} the limit points of $n(\Lambda_n-T_*)$ fill the interval $[\tau_c-W,\tau_c+W]$, so $D\ge W=(2(\im u_*)^2)^{-1}$; a lower bound $\PP(\im u_*<\varrho)\ge c\varrho^2$, which we expect but have not proved, would give $\PP(D>x)\ge c'x^{-1}$, so that no exponent above $1$ is possible. (4) The exponents $1/900$ and $1/(1782E_c)$ come from the exponent $10$ of Theorem~\ref{thm:collision-explicit} and a union bound over eleven factors, and are not optimised. \end{remark} \section{Remarks and questions}\label{sec:remarks} \subsection{Rates}\label{subsec:rates-questions} The exponents in Theorem~\ref{thm:rate}, any $\alpha<\frac3{50}$ and any $\theta<\frac1{10}$ at fixed $t$, do not depend on $\eps$, but the constants are not polynomial in $1/\eps$ (Remark~\ref{rem:rate}). By Remark~\ref{rem:rate} this argument cannot give a fixed-$t$ exponent above $\frac19$; the limit is the degree of the explicit onset of Theorem~\ref{thm:quant}, through the vertex loss of Lemma~\ref{lem:ext-disk} near the edges $\pm2$ and the bulk constant. The distortion of the outer factor has an $\eps$-independent small-ball exponent, uniformly in $t\ge1+\eps$ (Proposition~\ref{prop:distortion}(b)), and the natural alternative deformation is no better (Remark~\ref{rem:homotopy}). \begin{question}\label{q:distortion} (i) For the Gaussian power series, is there a rate $n^{-\alpha}$ in Theorem~\ref{thm:rate} with $\alpha>\frac1{10}$, or with constants polynomial in $1/\eps$? The Monte Carlo exponents of Remark~\ref{rem:rate} are $0.4$ to $0.82$. (ii) Does Proposition~\ref{prop:distortion}(b) hold for every $\beta<1$, uniformly in $t\ge1+\eps$? (iii) Is there a deformation of $a^{(t)}$ into a polynomial with the same zeros in $\D$ along which the minimum of $|a_s|$ on the unit circle and the coefficient bound of $a_s$ are controlled by powers of $\mu$, $|\eta_0|$ and $C_q$ whose exponents do not depend on $\eps$? By Remark~\ref{rem:homotopy}, any such deformation has circle minimum at most $\min_{|w|=1}|p_{\rm in}(w)|$, a product over the zeros of $a^{(t)}$ in $\D$. \end{question} \subsection{The constant in the final collision} \cite[Theorem~9.1]{Main} gives $|\Lambda_n-T_*|=O(1/n)$ with a constant that depends on $q$. By Theorem~\ref{thm:collision-explicit} this constant is bounded, when $T_*\ge1+\eps$, by a polynomial in condition parameters of $q$ with absolute degrees and coefficients polynomial in $1/\eps$, and for Gaussian coefficients $\PP(D>x,\ T_*\ge1+\eps)\le C_\eps x^{-1/900}$ (Theorem~\ref{thm:Dtail}). The constant $C_\eps$ is not polynomial in $1/\eps$, and with a constant polynomial in $1/\eps$ we have only an exponent of order $\eps/\log(1/\eps)$; so for $D$ itself we obtain only $\PP(D>x)\le C\exp(-c(\log x/\log\log x)^{1/2})$ (Corollary~\ref{cor:Dtail}). By Remark~\ref{rem:Dtail}(3), no bound better than $x^{-1}$ can be expected. \begin{question}\label{q:Dtail} Does $\PP(D>x)$ decay like a power of $x$ for the Gaussian power series, without the restriction $T_*\ge1+\eps$? By \cite[Theorem~7.13]{Main} this would follow from Theorem~\ref{thm:Dtail}(b) with $\log C_\eps=o(1/\eps)$, or from Theorem~\ref{thm:Dtail}(a) with an exponent bounded below. What is the true exponent of $\PP(D>x,\ T_*\ge1+\eps)$? \end{question} \section*{Statement on the use of AI} The proofs, computations and text of this paper were produced with AI systems from Anthropic, OpenAI and Google, directed by the author. The proofs were checked by independent AI referee rounds; no human expert has yet read them. \begin{thebibliography}{99} \bibitem[M]{Main} B. Baggins, \emph{The real-rooting time of heat-evolved random polynomials}, preprint (2026). \bibitem[B]{CompB} B. Baggins, \emph{Nonreal zeros of heat-evolved Weyl polynomials near the collapse time}, preprint (2026). \bibitem[C]{CompC} B. Baggins, \emph{Exact counts of nonreal zeros beyond the Weyl weight: varying weights, coefficient profiles and the Laguerre flow}, preprint (2026). \bibitem{An} T. W. Anderson, \emph{The integral of a symmetric unimodal function over a symmetric convex set and some probability inequalities}, Proc. Amer. Math. Soc. \textbf{6} (1955), 170--176. \bibitem{AW} J.-M. Aza\"is, M. Wschebor, \emph{Level Sets and Extrema of Random Processes and Fields}, Wiley, Hoboken, NJ, 2009. \bibitem{Fo} P. J. Forrester, \emph{The limiting Kac random polynomial and truncated random orthogonal matrices}, J. Stat. Mech. Theory Exp. (2010), P12018; doi:10.1088/1742-5468/2010/12/P12018; arXiv:1009.3066. \bibitem{Fuj} M. Fujiwara, \emph{\"Uber die obere Schranke des absoluten Betrages der Wurzeln einer algebraischen Gleichung}, T\^ohoku Math. J. \textbf{10} (1916), 167--171. \bibitem{Ka} M. Kac, \emph{On the average number of real roots of a random algebraic equation}, Bull. Amer. Math. Soc. \textbf{49} (1943), 314--320. \bibitem{LPR} G. Lindgren, K. Podg\'orski, I. Rychlik, \emph{Effective computations of joint excursion times for stationary Gaussian processes}, arXiv:2007.14220 (2020). \bibitem{MS} S. Matsumoto, T. Shirai, \emph{Correlation functions for zeros of a Gaussian power series and Pfaffians}, Electron. J. Probab. \textbf{18} (2013), no.~49, 1--18; doi:10.1214/EJP.v18-2545; arXiv:1212.6108. \end{thebibliography} \end{document}