Low degree approximation of random polynomials

TL;DR

This paper demonstrates that random Kostlan polynomials can be approximated by lower-degree polynomials without altering their zero set topology on the sphere, with high probability, and quantifies this relationship.

Contribution

It provides a quantitative link between the degree of approximation and probability, showing low-degree polynomials preserve the topology of Kostlan polynomial zero sets with high probability.

Findings

Zero set topology preserved under low-degree approximation

Probability bounds for complex topological types

Exponential rarity of certain topologies in random polynomials

Abstract

We prove that with "high probability" a random Kostlan polynomial in $n+1$ many variables and of degree $d$ can be approximated by a polynomial of "low degree" without changing the topology of its zero set on the sphere $S^n$. The dependence between the "low degree" of the approximation and the "high probability" is quantitative: for example, with overwhelming probability the zero set of a Kostlan polynomial of degree $d$ is isotopic to the zero set of a polynomial of degree $O(\sqrt{d \log d})$. The proof is based on a probabilistic study of the size of $C^1$-stable neighborhoods of Kostlan polynomials. As a corollary we prove that certain topological types (e.g. curves with deep nests of ovals or hypersurfaces with rich topology) have exponentially small probability of appearing as zero sets of random Kostlan polynomials.

Figures

Equations

$$P(x)=\sum _{|\alpha |=d}{\xi }_{\alpha }\cdot \left({\left(\frac{d!}{{\alpha }_0!\cdots {\alpha }_n!}\right)}^{1/2}{x}_0^{{\alpha }_0}\cdots x_n^{{\alpha }_n}\right),$$

$$\left\{\Vert f-p{\Vert }_{C^1}<\frac{\delta (p)}{2}\right\}⟹(S^n,Z(p))\sim (S^n,Z(f)),$$

$$p{|}_L=\sum _{\ell \le L}{p}_{\ell }.$$

$$\mathbb{P}\left\{{\Vert p-p{|}_L\Vert }_{C^1}<\frac{\delta (p)}{2}\right\}\ge 1-\left(c_5(n)d^{\frac{5n}{2}+2}{L}^{2n-1}{e}^{-\frac{L^2}{d}}{\sigma }^2+\frac{1}{\sigma }\right).$$

$${\Vert \sum _{|\alpha |=d}{\gamma }_{\alpha }{x}_0^{{\alpha }_0}\cdots x_n^{{\alpha }_n}\Vert }_{\mathrm{BW}}={\left(\sum _{|\alpha |=d}{\gamma }_{\alpha }^2\frac{{\alpha }_0!\cdots {\alpha }_n!}{d!}\right)}^{1/2}.$$

$${\mathcal{P}}_{n,d}=\underset{d-\ell \in 2\mathbb{N}}{⨁}\Vert x{\Vert }^{d-\ell }{\mathcal{H}}_{n,\ell }.$$

$${\mathcal{S}}_{n,d}=\{\text{p:Sn→R such that p=P∣Sn with P∈Pn,d}\}={\mathcal{P}}_{n,d}{|}_{S^n}.$$

$${\mathcal{P}}_{n,d}\simeq {\mathcal{S}}_{n,d}\simeq {\mathbb{R}}^N\text{where N=(dn+d)}.$$

$$V_{n,\ell }=\{\text{h:Sn→R such that ΔSnh=−ℓ(ℓ+n−1)h}\}={\mathcal{H}}_{n,\ell }{|}_{S^n}.$$

$$\Vert p{\Vert }_{C^1}=\underset{\theta \in S^n}{\max }|p(\theta )|+\underset{\theta \in S^n}{\max }\Vert {\nabla }_{S^n}p(\theta )\Vert ,$$

$$\Vert p{\Vert }_{L^2}={\left({\int }_{S^n}p(\theta {)}^2\mathrm{d}\theta \right)}^{1/2},$$

$$\Vert p{\Vert }_{H^q}={\left(\sum _{d-\ell \in 2\mathbb{N}}{\ell }^{2q}\Vert p_{\ell }{\Vert }_{L^2}^2\right)}^{1/2}.$$

$${\mathcal{S}}_{n,d}=\underset{d-\ell \in 2\mathbb{N}}{⨁}{V}_{n,\ell }.$$

$$\Vert h_{n,\ell }{\Vert }_{L^2}=w_{n,d}(\ell )\Vert h_{n,\ell }{\Vert }_{\text{BW}},\Vert h_{n,\ell }{\Vert }_{H^q}={\ell }^q{w}_{n,d}(\ell )\Vert h_{n,\ell }{\Vert }_{\text{BW}}\forall h_{n,\ell }\in V_{n,\ell }\subset {\mathcal{S}}_{n,\ell }.$$

$$w_{n,d}(\ell )={\left(\text{vol}(S^n)\Gamma \left(\frac{n+1}{2}\right)\frac{\Gamma \left(\frac{d+\ell }{2}+1\right)}{\Gamma \left(\frac{n+1}{2}+\frac{d+\ell }{2}\right)}\frac{1}{2^d}\left(\genfrac{}{}{0ex}{}{d}{\frac{d-\ell }{2}}\right)\right)}^{1/2}.$$

$$\Vert p{\Vert }_{C^1}\le c_1(n)d^{\frac{1}{2}}\Vert p{\Vert }_{H^q}.$$

$$h_{\ell }(\phi )={\int }_{S^n}{h}_{\ell }(\theta )Z_{\ell }(\phi ,\theta )\mathrm{d}\theta .$$

$$Z_{\ell }({\theta }_1,{\theta }_2)=\sum _{j\in J_{\ell }}{y}_{\ell ,j}({\theta }_1)y_{\ell ,j}({\theta }_2)$$

$$\Vert Z_{\ell }({\theta }_1,\cdot ){\Vert }_{L^2}^2=⟨Z_{\ell }({\theta }_1,\cdot ),Z_{\ell }({\theta }_1,\cdot ){⟩}_{L^2}=Z_{\ell }({\theta }_1,{\theta }_2)=\dim (V_{n,\ell })=O({\ell }^{n-1}),$$

$$|h_{\ell }(\phi )|\le C_1(n){\ell }^{\frac{n-1}{2}}\Vert h_{\ell }{\Vert }_{L^2}\forall \phi \in S^n.$$

$${\partial }_j{h}_{\ell }(\phi )={\int }_{S^n}{h}_{\ell }(\theta ){\partial }_i{Z}_{\ell }(\phi ,\theta )\mathrm{d}\theta ,$$

$$|{\partial }_j{h}_{\ell }(\phi )|$$

$$\le \Vert h_{\ell }{\Vert }_{L^2}{\left({\int }_{S^n}|{\partial }_i{Z}_{\ell }(\phi ,\theta ){|}^2\mathrm{d}\theta \right)}^{1/2}$$

$$\le \Vert h_{\ell }{\Vert }_{L^2}{C}_2(n){\left({\ell }^2{\int }_{S^n}|Z_{\ell }(\phi ,\theta ){|}^2\mathrm{d}\theta \right)}^{1/2}\text{(by cite[cite][@@bibrefSeeley, Theorem 4])}$$

$$\le C_3(n)\Vert h_{\ell }\Vert {\ell }^{\frac{n+1}{2}}.$$

$$\Vert {\nabla }_{S^n}{h}_{\ell }(\phi )\Vert \le C_4(n){\ell }^{\frac{n+1}{2}}\Vert h_{\ell }{\Vert }_{L^2}\forall \phi \in S^n.$$

$$|p(\phi )|+\Vert {\nabla }_{S^n}p(\phi )\Vert$$

$$\le C_5(n)\sum _{d-\ell \in 2\mathbb{N}}{\ell }^{\frac{n+1}{2}}\Vert h_{\ell }{\Vert }_{L^2}$$

$$\le C_6(n){\left(\sum _{d-\ell \in 2\mathbb{N}}{\ell }^{n+1}\Vert h_{\ell }{\Vert }_{L^2}^2\right)}^{1/2}{\left(\frac{d}{2}\right)}^{1/2}$$

$$\le c_1(n)\sqrt{d}\Vert p{\Vert }_{H^q}\text{for q≥2n+1.}$$

$$v_1=\left(d^{\frac{n+1}{2}}\Vert h_d{\Vert }_{L^2},\cdots ,{\ell }^{\frac{n+1}{2}}\Vert h_{\ell }{\Vert }_{L^2},\cdots \right)\text{and}{v}_2=(1,\dots ,1)$$

$$\sum _{d-\ell \in 2\mathbb{N}}{\ell }^{\frac{n+1}{2}}\Vert h_{\ell }{\Vert }_{L^2}=⟨v_1,v_2{⟩}_{{\mathbb{R}}^{[d/2]}}\le \Vert v_1\Vert \Vert v_2\Vert ={\left(\sum _{d-\ell \in 2\mathbb{N}}{\ell }^{n+1}\Vert h_{\ell }{\Vert }_{L^2}^2\right)}^{1/2}{\left(\left[\frac{d}{2}\right]\right)}^{1/2}$$

$$P(x)=\sum _{|\alpha |=d}{\xi }_{\alpha }\cdot \left({\left(\frac{d!}{{\alpha }_0!\cdots {\alpha }_n!}\right)}^{1/2}{x}_0^{{\alpha }_0}\cdots x_n^{{\alpha }_n}\right),$$

Full text

1. Introduction

Over the past few years there has been an intense activity around the field of Random Algebraic Geometry, whose main interest has been studying topological properties of the zero set of random real algebraic equations.

This approach goes back to the classical work of Kac [10], who studied the expected number of real zeroes of a random polynomial in one variable whose coefficients are gaussian random variables, and was later extended and generalized in the 1990s to systems of equations in a sequence of influential papers by A. Edelman, E. Kostlan, M. Shub, S. Smale [3, 24, 4, 12, 26, 25]. More recently, in 2011, P. Sarnak [21] suggested to look at the connected components of a real algebraic curve from the random point of view, proposing a random version of Hilbert’s Sixteenth Problem (to investigate the “number, shape, and position” of the connected components of a real algebraic hypersurface [27]). Since then the area has seen much progress [7, 8, 9, 5, 19, 18, 21, 22, 16, 13, 14, 15], with a focus on the expectation of topological quantities such as the Betti numbers of random algebraic hypersurfaces [9, 8, 5].

In this paper we concentrate on the so called Kostlan model: we sample a random polynomial according to the rule

[TABLE]

with $\{\xi_{\alpha}\}_{|\alpha|=d}$ a family of independent, standard gaussian variables (see Section 3 below for more details). A main feature of this probabilistic model, in the univariate case, is that the expectation of the number of real zeroes of a Kostlan polynomial equals $\sqrt{d}$ [3]. This phenomenon is called “square-root law”: essentially the Kostlan polynomial seems to behave as if its degree is $\sqrt{d}$ rather than $d$. In higher dimensions a similar phenomenon happens to the Betti numbers of its zero set: their expectation is of the order $O(d^{n/2})$, while the deterministic upper bound is $O(d^{n}).$ In this paper we give a further contribution in this direction, by proving the following theorem.

Theorem A** (Low-degree approximation).**

Let $P$ be a random Kostlan polynomial of degree $d$ and $n+1$ many variables and denote by $p=P|_{S^{n}}$ its restriction to the unit sphere $S^{n}\subset\mathbb{R}^{n+1}$ and by $Z(p)\subset S^{n}$ its zero set on the sphere. As $d\to\infty$, with overwhelming probability the pair $(S^{n},Z(p))$ is diffeomorphic to the pair $(S^{n},Z(q))$ where $q$ is the restriction to the sphere of a polynomial of degree $O(\sqrt{d\log d})$.

The idea of the proof of the previous theorem is the following. Thom’s isotopy Lemma implies that, given a function $p:S^{n}\to\mathbb{R}$ whose zero set $Z(p)\subset S^{n}$ is nonsingular, there is a small $C^{1}$ neighborhood (we call it a “stable neighborhood”) such that all functions in this neighborhood have zero sets diffeomorphic to $Z(p)$. However, how large this neighborhood can be depends on $p$ and in Proposition 3 we prove that it contains a $C^{1}$-ball:

[TABLE]

where $\delta(p)$ denotes the distance, in the Bombieri-Weil norm, from $p$ to the set of polynomials with a singular zero set (the “discriminant”, see Section 4). In order to produce a low-degree approximation of $p$, we first write it as $p=\sum_{\ell}p_{\ell}$, where each $p_{\ell}$ denotes the projection of $p$ to the space of spherical harmonics of degree $\ell$, and then take only the part of degree smaller than $L$ of this expansion:

[TABLE]

We will prove that, choosing $L=O(\sqrt{d\log d}),$ with overwhelming probability the difference $p-p|_{L}$ has small enough $C^{1}$-norm to be contained in the above stable neighborhood.

From the technical point of view this last step requires three estimates: we first bound the $C^{1}$-norm of $p-p|_{L}$ with its Sobolev norm (Proposition 1), then the Sobolev norm with the Bombieri-Weil norm of the original polynomial (which is the norm endowing the space of polynomials with the Kosltan gaussian measure, Proposition 2) and finally we estimate the size (i.e. the probability) of a small neighborhood of the discriminant (Proposition 4).

1.1. Consequences

All the previous estimates are quantitative and produce different outcomes for different choices of the degree $L$ to which we truncate the expansion of $p$. The most general bound that we obtain is the following (Theorem 5 below): there exists $c_{5}(n)>0$ such that for every $L,\sigma>1$ we have:

[TABLE]

For example, choosing $L$ to be a fraction of $d=\deg(p)$, the above $\sigma$ can be tuned so that the probability from the statement of Theorem A goes exponentially fast to one as $d\to\infty$.

We use this idea to constraint the typical topology of $(S^{n},Z(p))$ as follows: (i) we identify a “family” of topological types (e.g. hypersurfaces of the sphere $S^{n}$ with more than $\alpha d^{n}$ many components); (ii) we show that we need at least degree $L_{d}$ to realize this topological type (e.g. we need degree at least $c(\alpha)d$ to have $\alpha d^{n}$ many components); (iii) we prove that with “high probability” $p$ can be stably-approximated by a polynomial of degree smaller than $L_{d}$ (which implies its zero set cannot have that topological type). Here are two examples of the application of this strategy.

(Theorem 8 below) The probability that the zero set on $S^{n}$ of a Kostlan polynomial of degree $d$ has total Betti number larger than $\alpha d^{n}$ is bounded by $\gamma_{1}(\alpha)e^{-\gamma_{2}(\alpha)d}$ for some constants $\gamma_{1}(\alpha),\gamma_{2}(\alpha)>0$. This was known for the case $n=1$ (points on $S^{1}$) and for the case $n=2$ (algebraic curves) [6], but only in the case of maximal curves, see Remark 5 below.

(Theorem 9 below) The probability that the zero set on $S^{n}$ of a Kostlan polynomial of degree $d$ contains a nest of depth $\alpha d^{n}$ is bounded by $\gamma_{1}(\alpha)e^{-\gamma_{2}(\alpha)d}$ for some constants $\gamma_{1}(\alpha),\gamma_{2}(\alpha)>0$, see Figure 1.

Remark 1*.*

Since the low-degree approximation from Theorem A is the projection of $p$ to low-degree harmonics, this can be used, in the case $n=1$ and with high probability, to improve the complexity of a certain class of algorithms in real algebraic geometry (e.g. adaptive algorithms for real root isolation), essentially showing that “for most polynomials” the bound on the complexity of these algorithms is better than the absolute deterministic bound. We plan to elaborate on this idea in a forthcoming work.

Acknowledgements

The authors are indebted to Marie-Françoise Roy, who has played a crucial role for the existence of this paper. She is a friend and an intellectual guide.

2. Spaces of polynomials and norms

We denote by $\mathcal{P}_{n,d}=\mathbb{R}[x_{0},\ldots,x_{n}]_{(d)}$ the space of real homogeneous polynomials of degree $d$. We endow $\mathcal{P}_{n,d}$ with the Bombieri-Weil norm, which is defined as follows: writing a homogeneous polynomial in the monomial basis we set:

[TABLE]

For every $\ell=0,\ldots,d$ we will also consider the space $\mathcal{H}_{n,\ell}\subset\mathcal{P}_{n,\ell}$ of homogeneous harmonic polynomials, i.e. polynomials $H$ such that $\Delta_{\mathbb{R}^{n+1}}H=0$. It turns out that the space $\mathcal{P}_{n,d}$ can be decomposed as:

[TABLE]

The decomposition (2.2) has two important properties (see [11]):

• (i)

Given a scalar product which is invariant under the action of $O(n+1)$ on $\mathcal{P}_{n,d}$ by change of variables, the decomposition (2.2) is orthogonal for this scalar product.

• (ii)

The action of $O(n+1)$ on $\mathcal{P}_{n,d}$ preserves each $\mathcal{H}_{n,\ell}$ and the induced representation on the space of harmonic polynomials is irreducible. In particular there exists a unique, up to multiples, scalar product on $\mathcal{H}_{n,\ell}$ which is $O(n+1)$-invariant.

The space $\mathcal{P}_{n,d}$ injects (by taking restrictions of polynomials) into the space $C^{\infty}(S^{n},\mathbb{R})$ of smooth functions on the unit sphere $S^{n}\subset\mathbb{R}^{n+1}$. We denote by

[TABLE]

the image of such injection. In particular the two vector spaces $\mathcal{P}_{n,d}$ and $\mathcal{S}_{n,d}$ are isomorphic:

[TABLE]

We introduce the following convention: given $P\in\mathcal{P}_{n,d}$ we denote by $p=P|_{S^{n}}$ (i.e. we will use capital letters for polynomials in $\mathcal{P}_{n,d}$ and small letters for their restrictions in $\mathcal{S}_{n,d}$). Restricting polynomials in $\mathcal{H}_{n,\ell}$ to the unit sphere we obtain exactly eigenfunctions of the spherical laplacian:

[TABLE]

We will consider various norms on $\mathcal{S}_{n,d}$ (all these norms are in fact defined on $C^{\infty}(S^{n},\mathbb{R})$):

1. (1)

The Bombieri-Weil norm, simply defined for $p=P|_{S^{n}}$ as $\|p\|_{\mathrm{BW}}=\|P\|_{\mathrm{BW}}.$ Note that the same $p:S^{n}\to\mathbb{R}$ can be the restriction of two different $P_{1}\in\mathcal{P}_{n,d_{1}}$ and $P_{2}\in\mathcal{P}_{n,d_{2}}$ (for example: take $P_{2}(x)=\|x\|^{2}P_{1}(x)$), it is therefore important for the computation of the Bombieri-Weil norm to specify the space where $p$ comes from, i.e. its original homogeneous degree. 2. (2)

The $C^{1}$-norm defined for $p\in\mathcal{S}_{n,d}$ as:

[TABLE]

where $\nabla_{S^{n}}p$ denotes the spherical gradient, i.e. the orthogonal projection on the unit sphere of the gradient of $p$. 3. (3)

The $L^{2}$-norm, defined for $p\in\mathcal{S}_{n,d}$ as:

[TABLE]

where “$\mathrm{d}\theta$” denotes integration with respect to the standard volume form of the sphere. In the sequel we will denote by $\{y_{\ell,j}\}_{j\in J_{\ell}}$ a chosen $L^{2}$-orthonormal basis of $V_{n,\ell}$. 4. (4)

The Sobolev $q$-norm, defined for $p=\sum_{\ell}p_{\ell}$ (decomposed as in (2.2)) by:

[TABLE]

(Note that $\|\cdot\|_{H^{0}}=\|\cdot\|_{L^{2}}$.)

The decomposition (2.2) induces a decomposition:

[TABLE]

By property (i) above this decomposition is orthogonal both for the Bombieri-Weil, the $L^{2}$ and the Sobolev scalar products. Moreover, because of property (ii) above, the Bombieri-Weil scalar product, the $L^{2}$ and the Sobolev one are one multiple of the others on $V_{n,\ell}$ (viewed as a subspace of $\mathcal{P}_{n,d}$):

[TABLE]

The rescaling weights are given by (see [5, Example 1]):

[TABLE]

We observe also the following important fact: writing $P=\sum_{\ell}P_{\ell}$ with each $P_{\ell}\in\|x\|^{d-\ell}\mathcal{H}_{n,\ell}$ as in (2.2), when taking restrictions to the unit sphere we have $p=\sum_{\ell}p_{\ell}$ with each $p_{\ell}$ the restriction to $S^{n}$ of a polynomial of degree $\ell$: in other words, the restriction to the unit sphere “does not see” the $\|x\|^{d-\ell}$ factor, which is constant on the unit sphere.

Proposition 1**.**

There exists a constant $c_{1}(n)>0$ such that for every $q\geq\frac{n+1}{2}$ and for every $p\in\mathcal{S}_{n,d}$ we have:

[TABLE]

Proof.

For the proof we use the fact that for every $\ell=0,\ldots,d$ the space $V_{n,\ell}$ with the $L^{2}$-scalar product is a reproducing kernel Hilbert space, i.e. there exists $Z_{\ell}:S^{n}\times S^{n}\to\mathbb{R}$ such that for every $h_{\ell}\in V_{n,\ell}$:

[TABLE]

The function $Z_{\ell}$ (the “zonal harmonic”) is defined as follows: letting $\{y_{\ell,j}\}_{j\in J_{\ell}}$ be an $L^{2}$-orthonormal basis for $V_{n,\ell}$ we set

[TABLE]

(written in this way (2.13) is easily verified). From this it follows that:

[TABLE]

where the last identity follows from [1, Proposition 5.7 (d)] and [1, Proposition 5.8]. We can therefore estimate $h_{\ell}$ using Cauchy-Schwartz in (2.13):

[TABLE]

Similarly, for every orthonormal frame field $\{\partial_{1},\ldots,\partial_{n}\}$ at $\varphi\in S^{n}$, differentiating (2.13) under the integral one obtains:

[TABLE]

and consequently, using again Cauchy-Schwartz:

[TABLE]

From this it follows that:

[TABLE]

Given now $p\in\mathcal{S}_{n,d}$ we write $p=\sum_{\ell}h_{\ell}$ with each $h_{\ell}\in V_{n,\ell}$, as in (2.9). Using (2.16) and (2.22) we can estimate for $\varphi\in S^{n}$:

[TABLE]

In the third inequality we have used Cauchy-Schwartz in $\mathbb{R}^{[d/2]}$ for the vectors:

[TABLE]

so that:

[TABLE]

Taking the supremum over $\varphi\in S^{n}$ in (2.26), this implies the statement.

∎

3. Gaussian measures and random polynomials

The space $\mathcal{P}_{n,d}$ can be turned into a gaussian space by sampling a random polynomial according to the rule:

[TABLE]

with $\{\xi_{\alpha}\}_{|\alpha|=d}$ a family of independent, standard gaussian variables. A random polynomial defined in this way is called a Kostlan polynomial. An alternative way for writing a random Kostlan polynomial is to expand it in the spherical harmonic basis:

[TABLE]

where $\{\xi_{\ell,j}\}_{\ell,j}$ is a family of independent, standard gaussian variables and $\{w_{n,d}(\ell)\}_{d-\ell\in 2\mathbb{N}}$ are given by (2.11).

Remark 2*.*

Observe that both

[TABLE]

are Bombieri-Weil orthonormal bases for $\mathcal{P}_{n,d}.$ More generally, given a basis $\{F_{k}\}_{k=1}^{N}$ for $\mathcal{P}_{n,d}$ which is orthonormal for the Bombieri-Weil scalar product, a random Kostlan polynomial can be defined by:

[TABLE]

where $\{\xi_{k}\}_{k=1}^{N}$ is a family of independent, standard gaussian variables.

Given $L\in\{0,\ldots,d\}$ we consider the projection $\mathcal{S}_{n,d}\to\mathcal{S}_{n,L}$ defined by expanding $p$ in spherical harmonics and taking only the terms of degree at most $L$ of this expansion:

[TABLE]

Proposition 2**.**

There exists a constant $c_{2}(n)>0$ such that for all $t,q\geq 0$ and for every $L\in\{0,\ldots,d\}$ we have:

[TABLE]

Proof.

First observe that, since $\{\left\|p-p|_{L}\right\|_{H^{q}}\leq t\|p\|_{\mathrm{BW}}\}\subset\mathcal{P}_{n,d}$ is a cone, denoting by $S^{N-1}$ the unit sphere in the Bombieri-Weil norm, the required probability equals:

[TABLE]

We will estimate the quantity

[TABLE]

from above using Markov inequality:

[TABLE]

where the expectation is computed sampling a polynomial $p$ uniformly from the unit Bombieri-Weil sphere.

More precisely, expanding $p$ in an $L^{2}$-orthonormal basis $\{y_{\ell,j}\}$ (so that $\{w_{n,d}(\ell)y_{\ell,j}\}$ is a Bombieri-Weil orthonormal basis)

[TABLE]

the condition that $p\in S^{N-1}$ writes $\sum_{\ell,j}\gamma_{\ell,j}^{2}=1$. Consequently, denoting as before “$\mathrm{d}\theta$” the integration with respect to the standard volume form of the sphere, we have:

[TABLE]

We use now the fact that the cardinality of $J_{\ell}$ is $O(\ell^{n-1})$ and that $N\sim\frac{d^{n}}{n!}$, obtaining the estimate:

[TABLE]

Moreover from (2.11) we easily get:

[TABLE]

Substituting (3.17) into (3.16) we get:

[TABLE]

For $y\in\mathbb{R}$ let us denote now by $\{y\}$ the nearest integer to $y$ with the same parity as $d.$ Then we can rewrite:

[TABLE]

We apply now the change of variable $y=x\sqrt{d}$ in the above integral, and obtain:

[TABLE]

In the last line we have used the fact that

[TABLE]

and the convergence is dominated by an integrable function (by De Moivre-Laplace theorem, see also [5, Lemma 6]).

Applying the change of variables $t=x^{2}/2$ we can reduce the last integral to an incomplete Gamma integral:

[TABLE]

For the inequality (3.25) we have used the asymptotic $\Gamma(s,x)\sim x^{s-1}e^{-x}$ for the incomplete Gamma integral:

[TABLE]

Finally, using the estimate (3.26) into (3.11) gives the desired inequality. ∎

Remark 3*.*

The final estimate (3.26) from Proposition 2 takes the following interesting shapes:

If $L=b\sqrt{d}$ with $b>0$, then:

[TABLE]

• -

If $L=\sqrt{bd\log d}$ with $b>0$, then:

[TABLE]

• -

If $L=d^{b}$ with $b\in(\frac{1}{2},1)$, then there exists $c_{1},c_{2}>0$ (depending on $b$) such that:

[TABLE]

• -

If $L=bd$ with $b\in(0,1)$, then there exists $c_{1},c_{2}>0$ (depending on $b$) such that:

[TABLE]

4. Stability

Let us consider the discriminant set $\Sigma_{n,d}\subset\mathcal{S}_{n,d}$ consisting of all those polynomials whose zero set on the sphere is singular:

[TABLE]

Given $p\in\mathcal{S}_{n,d}$ we denote by $\delta(p)$ its distance, in the Bombieri-Weil norm, to $\Sigma_{n,d}:$

[TABLE]

If $Z_{1},Z_{2}\subset S^{n}$ are two smooth hypersurfaces, we will write $(S^{n},Z_{1})\sim(S^{n},Z_{2})$ to denote that the two pairs $(S^{n},Z_{1})$ and $(S^{n},Z_{2})$ are diffeomorphic. Given $f\in C^{1}(S^{n},\mathbb{R})$ we denote by $Z(f)\subset S^{n}$ its zero set. A small perturbation in the $C^{1}$-norm of a function $f\in C^{1}(S^{n},\mathbb{R})$ whose zero set $Z(f)$ is nondegenerate does not change the class of the pair $(S^{n},Z(f))$; next Proposition makes this more quantitative.

Proposition 3**.**

Let $p\in\mathcal{S}_{n,d}\backslash\Sigma_{n,d}$. Given $f\in C^{1}(S^{n},\mathbb{R})$ such that $\|f-p\|_{C^{1}}<\frac{\delta(p)}{2}$, we have:

[TABLE]

Proof.

For $t\in[0,1]$ let us consider now the function $f_{t}=p+t(f-p)$. Since $\|f-p\|_{C^{1}}<\frac{\delta(p)}{2}$, for all $\theta\in S^{n}$ we have:

[TABLE]

Moreover, since $d\geq 1$, from $\|f-p\|_{C^{1}}<\frac{\delta(p)}{2}$ we also deduce $\frac{\|f-p\|_{C^{1}}}{\sqrt{d}}<\frac{\delta(p)}{2},$ which in turn implies for every $t\in[0,1]$ and $\theta\in S^{n}$:

[TABLE]

Recall from [20, Theorem 5.1] the following explicit expression for $\delta(p)$:

[TABLE]

Note that $\left(|p(\theta)|^{2}+\frac{\|\nabla_{S^{n}}p(\theta)\|^{2}}{d}\right)^{1/2}$ equals the distance in $\mathbb{R}^{2}$ between the two vectors $v_{1}(\theta)=(|p(\theta)|,0)$ and $v_{2}(\theta)=\left(0,\frac{\|\nabla_{S^{n}}p(\theta)\|}{\sqrt{d}}\right).$ Observe also that the two vectors $w_{1}(t,\theta)=(|f_{t}(\theta)|,0)$ and $w_{2}(t,\theta)=\left(0,\frac{\|\nabla_{S^{n}}f_{t}(\theta)\|}{\sqrt{d}}\right)$, in virtue of (4.5) and (4.6), satisfy:

[TABLE]

In particular:

[TABLE]

where the strict inequality comes from the fact that $w_{1}$ and $w_{2}$ belong to the interior of the balls.

Taking the minimum over $\theta\in S^{n}$ in the above expression gives:

[TABLE]

In particular the equation $\{f_{t}=0\}$ on $S^{n}$ is regular for all $t\in[0,1]$: whenever $f_{t}(\theta)=0$, then $\nabla_{S^{n}}f_{t}(\theta)$ cannot vanish because of the strict inequality in (4.11). The result follows now from Thom’s Isotopy Lemma.

∎

Next Proposition quantifies how large is the set of stable polynomials in the Bombieri-Weil norm.

Proposition 4**.**

There exists $c_{3}(n),c_{4}(n)>0$ such that for every $s\geq c_{4}(n)d^{2n}$ and for $p\in\mathcal{P}_{n,d}$:

[TABLE]

Proof.

Let $S^{N-1}\subset\mathcal{P}_{n,d}\simeq\mathcal{S}_{n,d}$ be the unit sphere for the Bombieri-Weil norm and consider the algebraic set:

[TABLE]

Observe that there exists a polynomial $Q:\mathbb{C}[z_{0},\ldots,z_{n}]\to\mathbb{C}$ (the discriminant polynomial) which vanishes exactly at polynomials whose zero set in the projective space $\mathbb{C}\mathrm{P}^{n}$ is singular, which has real coefficients and degree $(n+1)(d-1)^{n}$.

Note that if $P\in\Sigma_{n,d}$ then $\{P=0\}\subset\mathbb{C}\mathrm{P}^{n}$ is also singular; it follows that $\overline{\Sigma}$ is contained in $Z(Q)\cap\mathbb{R}[x_{0},\ldots,x_{n}]$ and we can apply [2, Theorem 21.1]. Denoting by $d_{\textrm{sin}}$ the sine distance111Strictly speaking $d_{\textrm{sin}}$ is not a metric on $S^{N-1}$, but rather on $\mathbb{R}\mathrm{P}^{N-1}$. in the sphere, [2, Theorem 21.1] tells that there exists a constant $C_{3}>0$ such that for all $s\geq(2(n+1)(d-1)^{n})N$ we have:

[TABLE]

Taking the cone over the set \{\textrm{p\in S^{N-1}$such that$\frac{1}{d_{\textrm{sin}}(p,\overline{\Sigma})}\geq s}\}, we can rewrite the previous inequality in terms of the Kostlan distribution, obtaining that for all $s\geq(2(n+1)(d-1)^{n})N$:

[TABLE]

Observe now that, since $N={d+n\choose d}$, for some constants $c_{3}(n),c_{4}(n)>0$ we have:

[TABLE]

In particular (4.15) finally implies that for all $s\geq c_{4}(n)d^{2n}$:

[TABLE]

∎

5. Low degree approximation

Theorem 5**.**

There exists $c_{5}(n)>0$ such that for every $L,\sigma>1$ we have:

[TABLE]

Remark 4*.*

Of course the previous statement is interesting if we can choose $L,\sigma>0$ in such a way that $\frac{1}{\sigma}$ goes to zero, but not too fast, and $L$ is significantly smaller than $d$, but not too small, because we still want the exponential term $e^{-\frac{L^{2}}{d}}$ to kill the other factors and make the probability go to one. The choice of $\sigma$ a polynomial in $d$ and $L=O(\sqrt{d\log d})$ is in some sense optimal for our proof, see next Proposition 6.

Proof.

Let $p\in\mathcal{S}_{n,d}$ and $L\in\{0,\ldots,d\}$. We have the following chain of inequalities:

[TABLE]

which hold for every $q\geq\frac{n+1}{2}$, $t>0$ and $s\geq c_{4}(n)d^{2n}$, with probability

[TABLE]

We now make the choices:

[TABLE]

With this choices we have:

[TABLE]

[TABLE]

[TABLE]

where we have set $c_{5}(n)=c_{2}(n)(3c_{1}(n)c_{4}(n))^{2}.$

Because of (5.7) we can apply the estimate in (5.4) which, using (5.8), becomes:

[TABLE]

Using (5.9), the last chain of inequalities holds with probability:

[TABLE]

∎

Proposition 6**.**

For every $a>0$ there exists $b>0$ such that as $d\to\infty$:

[TABLE]

with probability greater than $1-O({d^{-a}}).$

Proof.

Let $\sigma=d^{a}$ and $L=\sqrt{bd\log d}$. Then we have:

[TABLE]

where the last inequality holds for $b>0$ large enough. We apply now the previous Theorem 5 with this choice we have:

[TABLE]

∎

6. Applications to random topology

In this section we show how the previous results can be used to put constraints on the topological type of the pair $(S^{n},Z(p))$ for $p$ a random Kostlan polynomial. The first result is the following.

Theorem 7**.**

Let $p\in\mathcal{S}_{n,d}$ be a random Kostlan polynomial. As $d\to\infty$, with overwhelming probability the pair $(S^{n},Z(p))$ is diffeomorphic to the pair $(S^{n},Z(q))$ where $q$ is a polynomial of degree $O(\sqrt{d\log d})$.

Proof.

This follows from Proposition 6: in fact, by Proposition 3, (5.12) implies that the pairs $(S^{n},Z(p))$ and $(S^{n},Z(p|_{\sqrt{bd\log d}}))$ are diffeomorphic. ∎

6.1. Hypersurfaces with rich topology

For a topological space $X$ we denote by $b(X)$ the sum of its $\mathbb{Z}_{2}$-Betti numbers (sometimes also called the homological complexity of $X$). Recall by [17] that if $P\in\mathbb{R}[x_{0},\ldots,x_{n}]_{d}$, then the zero set of $p=P|_{S^{n}}$ has homological complexity bounded by $b(Z(p))\leq O(d^{n}).$

Theorem 8**.**

For $\alpha>0$ let $M_{\alpha,d}\subset\mathcal{S}_{n,d}$ be the set:

[TABLE]

Then there exist $\gamma_{1}(\alpha),\gamma_{2}(\alpha)>0$ such that:

[TABLE]

Proof.

Observe first that if $q\in\mathbb{R}[x_{0},\ldots,x_{n}]_{L}$ ($q$ is just a polynomial of degree $L$, not necessarily homogeneous), then $b(Z(q))\leq cL^{n}$ for some $c>0$, again by [17]. Hence, if we want $b(Z(q))>\alpha d^{n}$ we must have:

[TABLE]

Arguing as in the proof of Theorem 5, where now we take the projection $\lambda=\lambda_{L}:\mathcal{S}_{n,d}\to\mathcal{S}_{n,L}$ choosing the value $L=\left(\frac{\alpha}{c}\right)^{\frac{1}{n}}d$, we se that for every $t>0$ and $s\geq c_{4}(n)d^{2n}$:

[TABLE]

with probability

[TABLE]

Observe now that:

[TABLE]

for some constants $\gamma_{3}(\alpha),\gamma_{4}(\alpha)>0$, and

[TABLE]

for some constants $\gamma_{5}(\alpha),\gamma_{6}(\alpha)>0$.

Choosing $s$ as in (6.8) and $t$ as in (6.7), for $d>0$ large enough we have $s\geq c_{4}(n)d^{2n}$, and $c_{1}(n)d^{1/2}ts<\frac{1}{2}$; it follows that there exist constants $\gamma_{1}(\alpha),\gamma_{2}(\alpha)>0$ such that

[TABLE]

The condition $b(Z(p))>\alpha d^{n}$ implies that with the choice of $L<\left(\frac{\alpha}{c}\right)^{\frac{1}{n}}d$ we must have $\|p-p|_{L}\|_{C^{1}}\geq\frac{\delta(p)}{2},$ for otherwise the zero set of $p$ would be diffeomorphic to the zero set of $p-p|_{L}$ which, since $\deg(p-p|_{L})<L$, has homological complexity bounded by $b(Z(p-p|_{L}))<cL^{n}<\alpha d^{n}.$ In particular:

[TABLE]

which combined with (6.9) implies the statement.

∎

Remark 5*.*

It is not difficult to derive from Theorem 8 a similar result for random zero projective sets $Z(p)\subset\mathbb{R}\mathrm{P}^{n}$. In this context, the previous result should be compared with [6, Theorem 1], where the authors prove that the Kostlan measure of the set of curves $C\subset\mathbb{R}\mathrm{P}^{2}$ of degree $d$ whose number of components is more than $\frac{(d-1)(d-2)}{2}+1-ad$ is $O(e^{-c_{2}d})$. Theorem 8 is stronger in two senses: it applies to the general case of hypersurfaces in $\mathbb{R}\mathrm{P}^{n}$ and it gives exponential rarefaction for all sets of the form $\{b_{0}(Z(p))\geq\alpha d^{n}\}$ (i.e. not necessarily a linear correction from the maximal bound).

6.2. Depth of a nest

Given $p\in\mathcal{S}_{n,d}\backslash\Sigma_{n,d}$, its zero set $Z(p)\subset S^{n}$ consists of a finite union of connected, smooth and compact hypersurfaces. Fixing a point $y_{\infty}\in S^{n}$ (with $\mathbb{P}=1$ this point does not belong to $Z(p)$), every such component of $Z(p)$ separates the sphere $S^{n}$ into two open sets: a “bounded” one (the open set which does not contain $y_{\infty}$) and an “unbounded” one (the open set which contains $y_{\infty}$). The nesting graph of $Z(p)$ (with respect to $y_{\infty}$) is a graph whose vertices are the components of $Z(p)$ and there is an edge between two components if and only if one is contained in the bounded component of the other. The resulting graph is a forest (a union of trees) and we say that $(S^{n},Z(p))$ has a nest of depth $m$ if this forest contains a tree of depth $m$.

Theorem 9**.**

For $\alpha>0$ let $N_{\alpha d}\subset\mathcal{S}_{n,d}$ be the set:

[TABLE]

Then there exist $c_{1}(\alpha),c_{2}(\alpha)>0$ such that:

[TABLE]

Proof.

The proof is essentially the same as the proof of Theorem 8, after observing that the depth of every nest of the zero set of a polynomial of degree $L$ is smaller than $L$. ∎