Topology of the nodal set of random equivariant spherical harmonics on $S^3$

TL;DR

This paper investigates the topology of nodal sets of random equivariant spherical harmonics on the 3-sphere, revealing that they typically have a single component and their genus scales cubically with the degree.

Contribution

It provides a probabilistic analysis of the nodal topology of equivariant spherical harmonics, including expected genus calculations and nodal component counts.

Findings

Almost surely a single nodal component

Expected genus proportional to m((N^2 - m^2)/2 + N)

Genus scales as order N^3 when m/N is fixed

Abstract

We show that real and imaginary parts of equivariant spherical harmonics on $S^3$ have almost surely a single nodal component. Moreover, if the degree of the spherical harmonic is $N$ and the equivariance degree is $m$, then the expected genus is proportional to $m \left(\frac{N^2 - m^2}{2} + N\right) $. Hence if $\frac{m}{N}= c $ for fixed $0 < c < 1$, the genus has order $N^3$.

Equations

$${\Pi }_N^m:L^2({\mathbb{S}}^3)\to {\mathcal{H}}_N^m.$$

$${\mathbb{E}}_{N,m}g({\mathcal{Z}}_{u_N^m})$$

$${\mathbb{E}}_{N,m}g({\mathcal{Z}}_{u_N^m})=\frac{1+{\eta }^2}{8\pi }|m|\left(\frac{N^2-m^2}{2}+N\right)-|m|+1,$$

$$\eta =\frac{m/2}{\frac{N^2-m^2}{2}+N},$$

$$\frac{|m|(\#\{f_N^m=0\}-2)}{2}+1.$$

$$\mathbb{E}\#\{f_N^m=0\}=\frac{1+{\eta }^2}{4\pi }\left(\frac{N^2-m^2}{2}+N\right),$$

$$\eta =\frac{m/2}{\frac{N^2-m^2}{2}+N}.$$

$$(z_1,z_2)\mapsto \left(\begin{array}{c}\Re z_1\\ \Im z_1\\ \Re z_2\\ \Im z_2\end{array}\right)\in {\mathbb{S}}^3\subset {\mathbb{R}}^4$$

$$(\alpha ,\phi ,\theta )\mapsto \left(\begin{array}{c}\sin \alpha \cos (\theta +\phi )\\ \sin \alpha \sin (\theta +\phi )\\ \cos \alpha \cos (\theta -\phi )\\ \cos \alpha \sin (\theta -\phi )\end{array}\right)\in {\mathbb{S}}^3\subset {\mathbb{R}}^4,$$

$$e^{i\vartheta }.(z_1,z_2)=(e^{i\vartheta }{z}_1,e^{i\vartheta }{z}_2)$$

$$\pi :(z_1,z_2)\mapsto (2z_1\overline{z_2},|z_1{|}^2-|z_2{|}^2)\in \mathbb{C}\times \mathbb{R}.$$

$$(\alpha ,\phi ,\theta )\mapsto (\alpha ,\phi ,\theta +\vartheta )$$

$$\pi :(\alpha ,\phi ,\theta )\mapsto \left(\begin{array}{c}\sin (2\alpha )\cos (2\phi )\\ \sin (2\alpha )\sin (2\phi )\\ \cos (2\alpha )\end{array}\right)\in S^2\subset {\mathbb{R}}^3.$$

$$\hat{s}(z,\lambda ):={\lambda }^m(s(z)).$$

$${\Delta }_{{\mathbb{S}}^3}={\Delta }_H+\frac{{\partial }^2}{\partial {\theta }^2},$$

$$L^2({\mathbb{S}}^3,dV)=\underset{m\in \mathbb{Z}}{⨁}{\mathcal{H}}^m,$$

$$\widehat{{\nabla }_m^{\ast }{\nabla }_m(f(dz{)}^m)}={\Delta }_H\widehat{(f(dz{)}^m)}.$$

$${\Delta }_{{\mathbb{R}}^4}=4\sum _{j=1}^2\frac{{\partial }^2}{\partial z_j\partial {\bar{z}}_j}.$$

$${\mathcal{H}}_N=\underset{p+q=N}{⨁}{\mathcal{H}}_N^{(p,q)}.$$

$$\Re f_N^m=a_N^m(z),\Im f_N^m=b_N^m(z).$$

$$f_N^m(z)e^{-im\theta }=(a_N^m(z)+i{b}_N^m(z))(\cos m\theta -i\sin m\theta ),$$

$$\left\{\begin{array}[]{l}u_{N}^{m}=a_{N}^{m}\cos m\theta+b_{N}^{m}\sin m\theta,\\ v_{N}^{m}=b_{N}^{m}\cos m\theta-a_{N}^{m}\sin m\theta.\end{array}\right.$$

$${\mathcal{Z}}_{f_N^m}=\{z\in X:f_N^m(z)=0\}.$$

$${\mathcal{Z}}_{{\psi }_N^m}={\pi }^{-1}{\mathcal{Z}}_{f_N^m}.$$

$${\mathcal{N}}_{\Re {\psi }_N^m}\{p\in P_h:\Re {\psi }_N^m(p)=0\},resp.{\mathcal{N}}_{\Im {\psi }_N^m}=\{p\in P_h:\Im {\psi }_N^m(p)=0\}.$$

$${\psi }_N^m=u_N^m+i{v}_N^m.$$

$$\left\{\begin{array}[]{l}\Delta_{G}u_{N}^{m}=-\lambda_{N}^{m}u_{N}^{m},\\ \\ \Delta_{G}v_{N}^{m}=-\lambda_{N}^{m}v_{N}^{m},\\ \\ \frac{\partial}{\partial\theta}u_{j}=mv_{j},\;\;\frac{\partial}{\partial\theta}v_{j}=-mu_{j}.\end{array}\right.$$

$$\chi (X)=\chi (A)+\chi (B)-\chi (A\cap B),$$

$$\chi (X)=n\chi (M).$$

$$Z_{\Re \psi }-{\cup }_{j=1}^k\{(z_j,\theta ):\theta \in [0,2\pi ]\}$$

Full text

Research partially supported by NSF grant DMS-1810747. The first author is partially supported by Sloan Research Fellowship and by NSF grant DMS-1900993.

1. Introduction and statements of the results

In a recent article [JZ18] the authors proved that nodal sets of real or imaginary parts of equivariant (but non-invariant) eigenfunctions of Laplacians $\Delta_{KK}$ of generic ‘Kaluza-Klein’ metrics $g_{KK}$ on unit tangent (or cotangent) bundles $\pi:M\to X$ over Riemann surfaces $(X,g)$ have a single connected component. The generic condition is that [math] be a regular value for the eigenfunctions. The unit sphere ${\mathbb{S}}^{3}\subset{\mathbb{R}}^{4}$ with its standard metric and Laplacian has the standard Hopf fibration $\pi:{\mathbb{S}}^{3}\to{\mathbb{S}}^{2}$ and is an example of a Kaluza-Klein metric. It is a double cover ${\mathbb{S}}^{3}\to SO(3)\simeq U({\mathbb{S}}^{2})$ of the unit tangent bundle of ${\mathbb{S}}^{2}$. Recently, the nodal sets of random wave on $3$-dimensional Euclidean space have been the subject of numerical investigations by A. Barnett, Kyle Konrad, and Matthew Jin [BKJ17], which exhibit a surprising feature: only a small number of nodal components is visible in the computer graphics (Figure 1).

The results of Nazarov-Sodin [NS16] show that in fact there must be $cN^{3}$ distinct nodal components for some (very small) $c>0$, but the other components are evidently too small to be seen in the computer graphics. For this reason, it is conjectured that with probability one, the random spherical harmonic of fixed degree $N$ on ${\mathbb{S}}^{3}$ has one giant component and many much smaller components. In Figure 2, one does see a second small (red) component.

P. Sarnak posed the problem of finding its expected genus, and has proposed that the expected genus is of the order of magnitude of $N^{3}$ [Sar19]. Milnor has proved that the maximal genus of the zero set of a polynomial of degree $N$ is of this order of magnitude, so Sarnak’s proposal is that the nodal sets of random spherical harmonics on ${\mathbb{S}}^{3}$ are rather like Harnack curves (real algebraic curves of degree $N$ and of roughly maximal genus).

The purpose of this note is to link the results of [JZ18] to Sarnak’s proposal. Our main result is that the proposal is true for real/imaginary parts of random equivariant spherical harmonics of degree $N$ on ${\mathbb{S}}^{3}$. That is, we let $\mathcal{H}_{N}^{m}({\mathbb{S}}^{3})$ denote the space of spherical harmonics of degree $N$ which are equivariant of degree $m$ with respect to the $S^{1}$ action defining the Hopf fibration ${\mathbb{S}}^{3}\to{\mathbb{S}}^{2}$. For short, we say that $\psi_{N}^{m}$ is equivariant of degree $(N,m)$. We restrict the natural Gaussian measure on the space $\mathcal{H}_{N}\subset L^{2}({\mathbb{S}}^{3})$ of spherical harmonics of degree $N$ to the subspaces $\mathcal{H}_{N}^{m}({\mathbb{S}}^{3})$. Since $L^{2}({\mathbb{S}}^{3})=\bigoplus_{N=0}^{\infty}V_{N}\otimes V_{N}^{*}$ where $V_{N}$ is the $N$th irreducible representation of ${\mathbb{S}}^{3}$, fixing the weight $m$ of the $S^{1}$ action is the same as fixing one line in $V_{N}^{*}$, so that $\dim\mathcal{H}_{N}^{m}=\dim V_{N}=N+1$. We let $\gamma_{N}^{m}$ be the induced Gaussian measure on $\mathcal{H}_{N}^{m}$. Its covariance function is the Schwartz kernel $\Pi_{N}^{m}(x,y)$ of the orthogonal projection,

[TABLE]

We denote by

[TABLE]

the expected genus of the nodal set of the real part $u_{N}^{m}$ of an equivariant eigenfunction of degree $(N,m)$ with respect to $\gamma_{N}^{m}$.

Theorem 1.1**.**

Let $(\mathcal{H}_{N}^{m},\gamma_{N}^{m})$ be the Gaussian space of equivariant spherical harmonics of degree $(N,m)$. In the following, we assume $m\not=0$.

• (i)

With probability $1$ (w.r.t. $\gamma_{N}^{m}$), the nodal set $\mathcal{Z}_{u_{N}^{m}}$ of the real part $u_{N}^{m}=\Re\psi_{N}^{m}$ (resp. the nodal set $\mathcal{Z}_{v_{N}^{m}}$ of the imaginary part $v_{N}^{m}=\Im\psi_{N}^{m}$) of a random equivariant spherical harmonic $\psi_{N}^{m}\in\mathcal{H}_{N}^{m}$ has a single connected component;

• (ii)

The expected genus of the nodal component is by

[TABLE]

where

[TABLE]

which has modulus less than or equal to $1/2$.

The first statement is almost an application of the main result of [JZ18], where the ‘genericity’ assumption was only used to prove that each real (resp. imaginary) part of an equivariant eigenfunction has [math] as a regular value. Thus, to prove Theorem 1.1 (i) it is only necessary to prove the Bertini-type theorem that [math] is a regular value of $\psi_{N}^{m}\in\mathcal{H}_{N}^{m}$ with probability $1$. This is done in Section 5.4. The second statement has one topological simple part and one probabilistic part. The topological makes use of the identification of equivariant functions in $\mathcal{H}_{N}^{m}$ with sections $f_{N}^{m}e_{L}^{m}$ of the complex line bundle $\mathcal{O}(m)\to{\mathbb{S}}^{2}$; here $e_{L}$ denotes a local frame of $L$ over the affine chart, which we choose to be a holomorphic frame (see Section 2.3).

Lemma 1.2**.**

If [math] is a regular value of $\psi_{N}^{m}$, then the Genus of $\mathcal{Z}_{\Re\psi_{N}^{m}}$ is given by

[TABLE]

where $\psi_{N}^{m}$ is the lift of the section $f_{N}^{m}(z,\bar{z})e_{L}^{m}.$

As mentioned in [JZ18], the key point of the proof is to show that $\pi:{\mathcal{N}}_{u_{m,j}}\to X$ is a kind of ‘helicoid cover’. That is, it is an $m$-fold cover over the complement of the zeros of $f_{N}^{m}$, while the inverse image of a zero is an $S^{1}$ orbit. Thus, it is not a branched cover in the standard sense; rather $\pi$ is locally like the projection of a vertical helicoid onto the horizontal plane.111We thank J. Y. Welschinger for discussions of this local picture and for the Euler characteristic calculation of Lemma 1.2.

The second part is the following Kac-Rice type calculation.

Lemma 1.3**.**

With probability $1$, $f_{N}^{m}$ has isolated non-degenerate zeros, and

[TABLE]

where

[TABLE]

To determine the expected number of zeros of $f_{N}^{m}$ we use the Kac-Rice formula for Gaussian random eigen-sections of a line bundle, which gives an integral formula for the expected number of zeros in terms of the determinant of a matrix formed from values and derivatives of the covariance (two-point) function of the Gaussian random function.

It follows that when $\delta<\frac{m}{N}<1-\delta$ for some small fixed $\delta>0$, the genus of $\mathcal{Z}_{u_{N}^{m_{N}}}$ of order $N^{3}$. The spherical harmonics in $\bigcup_{\begin{subarray}{c}|m|\leq N\\ 2|N-m\end{subarray}}\mathcal{H}_{N}^{m}$ form a set of $N+1$ subspaces of dimension $N+1$ in $\mathcal{H}_{N}$, which intersect the unit sphere $S\mathcal{H}_{N}$ in a union of $N+1$ great spheres of dimension $N$. A random spherical harmonic in $\mathcal{H}_{N}$ is a random linear combination of components in $\mathcal{H}_{N}^{m}$. It is tempting to imagine that one can determine the expected Euler characteristic of the nodal set of a random spherical harmonic in $\mathcal{H}_{N}$ by perturbing the special equivariant spherical harmonics in $\mathcal{H}_{N}^{m}$ (or more precisely, their real and imaginary parts). The genus of the nodal set of $u_{N}^{m}=\Re\psi_{N}^{m}$ is so large that a small perturbation is unlikely to decrease the genus. But this is just speculation.

1.1. Acknowledgements

We thank Peter Sarnak for discussions on the genus of the giant component which prompted this article. We also thank Jean-Yves Welschinger for discussions of the Euler characteristic calculation at the outset. We thank Igor Wigman for helpful discussions.

2. Background and notations

In this section, we first review the results of [JZ18] on eigenfunctions of Kaluza-Klein Laplacians on general $S^{1}$ bundles over Riemann surfaces. We then specialize to the circle bundle ${\mathbb{S}}^{3}\to{\mathbb{S}}^{2}$ and introduce some notation and background concerning eigenfunctions of the standard Laplacian $\Delta_{{\mathbb{S}}^{3}}$on ${\mathbb{S}}^{3}$, on irreducible representations of $SU(2)={\mathbb{S}}^{3}$ and on sections of associated complex line bundles.

2.1. Review of the results of [JZ18]

The main result is [JZ18, Theorem 1.5] pertains to nodal sets of real/imaginary parts of equivariant eigenfunctions of the Kaluza-Klein Laplacian on $S^{1}$ bundles $\pi:M\to X$ over Riemannian surfaces $X$. A Kaluza-Klein metric $G$ is specified by a metric $h$ on $X$ and a connection $\nabla$ on $\pi:M\to X$.

Theorem 2.1**.**

Suppose that the data $(g,h,\nabla)$ of the Kaluza-Klein metric satisfies the generic properties of [JZ18, Theorem 1.4]. Then,

1. (1)

The eigenspace of $\Delta_{G}$ corresponding to $\lambda=\lambda_{m,j}=\lambda_{-m,j}$ is spanned by $\varphi_{m,j}$ and $\varphi_{-m,j}=\overline{\varphi_{m,j}}$. In particular, any real eigenfunction with the eigenvalue $\lambda_{m,j}$ is a constant multiple of $T_{\theta}\left(\Re\varphi_{m,j}\right)$, where $T_{\theta}$ is the $S^{1}$ action on $P$ parameterized by $\theta$. 2. (2)

For $m\neq 0$, the nodal sets of $\Re\varphi_{m,j}$ are connected. 3. (3)

For $m\neq 0$, the number of nodal domains of $\Re\varphi_{m,j}$ is $2$.

The principal ‘generic property’ referred to in the statement is that [math] is a regular value of the (complex-valued) equivariant eigenfunctions. It is also proved that for generic Kaluza-Klein metrics, the eigenspaces are of real dimension $2$ (corresponding to an equivariant eigenfunction and its complex conjugate).

Although the standard metric on ${\mathbb{S}}^{3}$ is Kaluza-Klein with respect to the standard metric on ${\mathbb{S}}^{2}$ and Riemannian connection, the above theorem does not apply to the setting of the present paper. Indeed, the eigenspace of $\Delta_{{\mathbb{S}}^{3}}$ corresponding to spherical harmonics of degree $N$ is of dimension $N^{2}$, quite opposite to the ‘multiplicity two’ property for generic Kaluza-Klein metrics. However, this property was not used to prove (2)-(3), but only to ensure that for generic K-K metrics, all eigenfunctions are equivariant, i.e. the Laplace eigenvalue determines the equivariance degree. This is certainly not true for $\Delta_{{\mathbb{S}}^{3}}$. In Theorem 1.1, however, we restrict to the subspaces of Laplace eigenfunctions with a fixed equivariance degree and prove the result for them, as long as $m\not=0$. We do not allow linear combinations of Laplace eigenfunctions of different equivariance degrees.

The one ‘generic property’ we do need is that [math] is a regular value of the eigenfunctions. This is not necessarily the case for every equivariant eigenfunction, but we prove it is the case almost surely for random equivariant spherical harmonics (Theorem 5.5). Once this is established, the proof of [JZ18, Theorem 1.5] applies with no change to equivariant spherical harmonics, and allows us to conclude that (2)-(3) hold. See Section 6 for a final summary of the argument.

2.2. Coordinates on ${\mathbb{S}}^{3}$ and Hopf fibration

We use two coordinate systems on ${\mathbb{S}}^{3}$:

[TABLE]

where $z_{1},z_{2}\in\mathbb{C}$ and $|z_{1}|^{2}+|z_{2}|^{2}=1$, or

[TABLE]

where $\alpha\in[0,\pi/2]$, $\varphi\in[0,2\pi]$, and $\theta\in[-\pi,\pi]$. In the first coordinate system (2.1), the action of $S^{1}$ is given by

[TABLE]

and the Hopf map $\pi:{\mathbb{S}}^{3}\to{\mathbb{S}}^{2}$ is given by

[TABLE]

In the second coordinate system, (2.3) is equivalent to

[TABLE]

and the Hopf map (2.4) is

[TABLE]

2.3. ${\mathbb{S}}^{3}$ as a Kaluza-Klein 3-fold

The purpose of this section is to introduce the principal notation and terminology, and to explain how ${\mathbb{S}}^{3}$ (equipped with its standard metric) is an example of a Kaluza-Klein metric in the sense of [JZ18]. We also review the relation between equivariant Laplace eigenfunctions on ${\mathbb{S}}^{3}$ and associated eigensections of line bundles over $S^{2}$. In particular, we show that the results of [JZ18] are valid for random spherical harmonics in $\mathcal{H}_{N}^{m}$.

The 3-manifolds $M^{3}$ studied in [JZ18] are $S^{1}$ bundles $M^{3}\to X$ over Riemann surfaces $X$, in particular the unit tangent or cotangent bundles. A Kaluza-Klein metric on $M^{3}$ is a bundle-metric $g$ defined by a connection $\nabla$ on $TM^{3}$ and a Riemannian metric $h$ on $X$. The circular fibers are geodesics of the metric (in particular have a constant length) and the metric on the horizontal spaces is the lift of the metric on $X$. It is evident that the standard metric on ${\mathbb{S}}^{3}$ is Kaluza-Klein with respect to the Hopf fibration $\pi:{\mathbb{S}}^{3}\to{\mathbb{S}}^{2}$. Harmonic analysis on spheres and its relation to the Hopf fibration ${\mathbb{S}}^{3}\to{\mathbb{S}}^{2}$ is elementary and well-known (see e.g. [Fol72]), so we only review two aspects of it: (i) relating eigensections of the Bochner Laplacians on $L^{m}$ to equivariant eigenfunctions of $\Delta_{{\mathbb{S}}^{3}}$; (ii) CR geometry of ${\mathbb{S}}^{3}$.

Associated to the principal $S^{1}$ bundle ${\mathbb{S}}^{3}\to{\mathbb{S}}^{2}$ are the line bundles $L^{m}={\mathbb{S}}^{3}\times_{\chi^{m}}{\mathbb{C}}$ where $\chi(e^{i\theta})=e^{i\theta}$. For $m=1$, $L=\mathcal{O}(1)\to{\mathbb{C}}{\mathbb{P}}^{1}$ in the notation of algebraic geometry [GH78]; it is the spin-bundle, i.e. the square root of the anti-canonical bundle, $K_{{\mathbb{C}}{\mathbb{P}}}^{-1}=(T^{1,0})^{\frac{1}{2}}{\mathbb{C}}{\mathbb{P}}^{1}$. When $m=-1$, $L^{-1}=K_{{\mathbb{C}}{\mathbb{P}}^{1}}\simeq T^{*1,0}{\mathbb{C}}{\mathbb{P}}^{1}$, the canonical bundle. In a standard way, we view ${\mathbb{S}}^{3}\subset L^{*}$ as the unit bundle with respect to the Fubini-Study metric. The connection in this setting is the Chern connection of the Fubini-Study metric; we refer to [GH78, JZ18] for background.

Sections $s$ of $L^{m}\to{\mathbb{S}}^{2}$ naturally lift to $L^{*}$ by

[TABLE]

Thus, the restriction of the lift of $s\in C(X,L^{m})$ ${\mathbb{S}}^{3}$ satisfies $\hat{s}(r_{\theta}x)=e^{im\theta}\hat{s}(x)$. We refer to lifts of sections of $L^{m}$ as ‘equivariant’ functions on ${\mathbb{S}}^{3}$ and denote the space of such functions by $\mathcal{H}^{m}$.

The standard Laplacian $\Delta_{{\mathbb{S}}^{3}}$ is a Kaluza-Klein Laplacian, i.e has the form,

[TABLE]

where $\Delta_{H}$ is the horizontal Laplacian. The fact that the fiber Laplacian is $\frac{\partial^{2}}{\partial\theta^{2}}$ reflects the fact that $S^{1}$ orbits are geodesics isometric to ${\mathbb{R}}/2\pi{\mathbb{Z}}$. It is obvious that $[\Delta_{{\mathbb{S}}^{3}},\frac{\partial^{2}}{\partial\theta^{2}}]=0.$ Since $S^{1}$ acts isometrically on $({\mathbb{S}}^{3},G)$ we may decompose into its weight spaces,

[TABLE]

where $\mathcal{H}^{m}=\{F:{\mathbb{S}}^{3}\to{\mathbb{C}}:F(e^{i\theta}.x)=e^{im\theta}F(x)\}.$ The weight spaces are $\Delta^{H}$-invariant, i.e., $\Delta_{H}:\mathcal{H}^{m}\to\mathcal{H}^{m}$.

Definition 2.2**.**

We define $\mathcal{H}_{N}^{m}$ to be the subspace of degree $N$ spherical harmonics in $\mathcal{H}^{m}$. We call a function belonging to $\mathcal{H}_{N}^{m}$ an equivariant spherical harmonic of degree $(N,m)$.

The lifting map gives a canonical identification $\mathcal{H}_{m}\cong L^{2}(X,L^{m}).$ The Bochner Laplacian $\nabla_{m}^{*}\nabla_{m}$ corresponds to the horizontal Laplacian under this identification, i.e.

[TABLE]

Since ${\mathbb{S}}^{3}$ is a group, $L^{2}({\mathbb{S}}^{3})=\bigoplus_{N=0}^{\infty}V_{N}\otimes V_{N}$ where $V_{N}$ is an irreducible representation of ${\mathbb{S}}^{3}$ of dimension $N+1$. Moreover, $\Delta|_{V_{N}\otimes V_{N}}=N(N+2)=(N+1)^{2}-1$.

2.4. CR structure

Another viewpoint is that ${\mathbb{S}}^{3}=\partial B\subset{\mathbb{C}}^{2}$, i.e. that ${\mathbb{S}}^{3}$ is the boundary of the unit ball, a strictly pseudo-convex domain in ${\mathbb{C}}^{2}$. A defining function for ${\mathbb{S}}^{3}\subset{\mathbb{C}}^{2}$ is the usual Euclidean distance $r=|Z|$ from the origin. Here, $Z\in{\mathbb{C}}^{2}$. Let $\overline{\partial}r$ be the associated $(0,1)$ form. The theory of spherical harmonics on ${\mathbb{S}}^{3}$ has been related to the CR geometry and the representation theory of $SU(2)={\mathbb{S}}^{3}$ in [Fol72].

Spherical harmonics on ${\mathbb{S}}^{3}$ are restrictions of homogeneous harmonic polynomials on ${\mathbb{C}}^{2}$. In complex coordinates $Z=(z_{1},z_{2})$, the Euclidean Laplacian is

[TABLE]

Let $\mathcal{H}_{N}^{(p,q)}$ denote the space of harmonic homogeneous polynomials of degree $N$ on ${\mathbb{C}}^{2}$ which are of degree $p$ in $z_{j}$’s and of degree $q$ in the $\bar{z}_{k}$’s; $N=p+q$. Then $\mathcal{H}^{(p,q)}$ is an irreducible representation of $SU(2)$ and the space of all spherical harmonics of degree $N$ admits the decomposition

[TABLE]

The representation of $U(2)$ on $\mathcal{H}_{N}^{(p,q)}$ is denoted by $\rho(q,-p)$. One has $\rho(q,-p)|_{SU(2)}=\rho(q^{\prime},-p^{\prime})|_{SU(2)}\iff p+q=p^{\prime}+q^{\prime},$ and $\dim\rho(q,-p)=p+q+1$. Hence the decomposition (2.5) is another decomposition of $V_{N}\otimes V_{N}^{*}$ into irreducibles.

The orbits of the Hopf fibration of the action (2.3) define the characteristic directions of the CR manifold and lie in the null space of $\partial r|_{T{\mathbb{S}}^{3}}$. It follows that on polynomials $z^{p}\bar{z}^{q}$ of type $(p,q)$, $S^{1}$ acts by $e^{i(p-q)\theta}$. Thus, the equivariance degree of $\psi_{N}^{p,q}\in\mathcal{H}_{N}^{p,q}$ is $m=p-q$. Since $p+q=N$, $m=2p-N$, i.e. the data $(p,q)$ is equivalent to specifying only $p$ or $q$ or $m$. In particular, $S^{1}$ acts by $e^{iN\theta}$ on the space $\mathcal{H}_{N}^{N,0}$ of holomorphic polynomials $(p,q)=(N,0)$.

Proposition 2.3**.**

$\mathcal{H}_{N}^{(p,q)}=\mathcal{H}_{N}^{2p-N}$. In particular, the dimension of $\mathcal{H}_{N}^{m}$ is $N+1$ if $|m|\leq N$ and $2|N-m$, and [math] otherwise.

2.5. Associated sections of line bundles

It is useful to simultaneously keep in mind the ‘upstairs’ picture of equivariant eigenfunctions of $\Delta_{{\mathbb{S}}^{3}}$ and the ‘downstairs’ picture of sections of complex line bundles, as described in Section 2.2.

We denote the eigensection corresponding to $\psi_{N}^{m}\in\mathcal{H}_{N}^{m}$ as $f_{N}^{m}e_{L}^{m}$ in a local holomorphic frame $e_{L}$ of $L$, so that $f_{N}^{m}$ is a locally defined function on ${\mathbb{C}}{\mathbb{P}}^{1}$, i.e. is a function on the affine chart ${\mathbb{C}}$. Let

[TABLE]

Then,

[TABLE]

so that with $\psi_{N}^{m}=u_{N}^{m}+iv_{N}^{m}$,

[TABLE]

We denote by $\mathcal{Z}_{f_{N}^{m}}$ the zero set of the eigensection $f_{N}^{m}e_{L}^{m}$ on $X$:

[TABLE]

It is easy to see that the zero set $\mathcal{Z}_{\psi_{N}^{m}}$ of $\psi_{N}^{m}$ is the inverse image of $\mathcal{Z}_{f_{N}^{m}}$ under the natural projection $\pi$:

[TABLE]

Usually we study the nodal sets of the real and imaginary parts of the lift, not to be confused with the lifts of the real and imaginary parts of the local expression $f_{N}^{m}$ of the section (since the frame $e_{L}^{m}$ must also be taken into account). In general, it is not obvious whether or not the zero set of $f_{N}^{m}$ is discrete in $X$.

We denote the nodal sets of the real, resp. imaginary parts, of the lift by

[TABLE]

The analysis is the same for real and imaginary parts and we generally work with the imaginary part, following the tradition for quadratic differentials.

Our focus is on the nodal sets of the real or imaginary parts of

[TABLE]

Since $\Delta_{{\mathbb{S}}^{3}}$ is a real operator, the real and imaginary parts (2.7) satisfied the modified eigenvalue system,

[TABLE]

3. Genus of the nodal set: Proof of Lemma 1.2

In this section, we relate the genus of the nodal set of $\Re\psi_{N}^{m}$ to the number of zeros of $f_{N}^{m}$, under the assumption of all zeros of $f_{N}^{m}$ being regular (Lemma 1.2). We first recall a few lemmata:

Lemma 3.1**.**

Let $X$ be a topological space and let $A,B$ be topological subspaces whose interior cover $X$. Then

[TABLE]

where $\chi(\cdot)$ is the Euler characteristic of $\cdot$.

Lemma 3.2**.**

Let $X$ be a $n$-covering of $M$. Then we have

[TABLE]

Now we are ready to prove Lemma 1.2:

Proof.

To simplify the notation, let $\psi=\psi_{N}^{m}$ and $f=f_{N}^{m}$. Let $\{z_{j}\}_{j=1,2,\ldots,k}$ be the complete set of zeros of $f$. We first note that

[TABLE]

is $m$-covering of $k$-punctured sphere. So the Euler characteristic of

[TABLE]

is $m(2-k)$ by Lemma 3.2.

Now we apply Lemma 3.1 with $X=Z_{\Re\psi}$, $A=Z_{\Re\psi}-\cup_{j=1}^{k}\{(z_{j},\theta)~{}:~{}\theta\in[0,2\pi]\}$, and $B$ equal to a sufficiently small open neighborhood of $\bigcup_{j=1}^{k}\{(z_{j},\theta)~{}:~{}\theta\in[0,2\pi]\}$ in $Z_{\Re\psi}$. Then $B$ is homotopic to $\cup_{j=1}^{k}\{(z_{j},\theta)~{}:~{}\theta\in[0,2\pi]\}$, which has Euler characteristic equal to [math], and $A\cap B$ is $m$-covering of a disjoint union of punctured discs, which also has Euler characteristic equal to [math]. This implies that $X=Z_{\Re\psi}$ has Euler characteristic $m(2-k)$, and therefore the conclusion follows. ∎

We note that $c_{1}(L^{m})=m$ (first Chern number, the integral of the first Chern class). By the Hopf theorem $c_{1}(L^{m})$ is the sum over zeros of a smooth section $f_{N}^{m}$ with non-degenerate zeros of the index in ${\mathbb{Z}}_{2}$ of the zero. The index is the degree of the locally defined map $\frac{s(z)}{|s(z)|}$ from a small circle centered at the zero to ${\mathbb{C}}$ in a local trivialization ([BT82, Theorem 11.17]-[BT82, Proposition 12.8].) In particular if $s$ is a holomorphic section, then the indices are all equal to $1$ and $s$ has precisely $m$ zeros (counted with multiplicity). This only occurs in the case $m=N$. Otherwise, the sections in $\mathcal{H}_{N}^{m}$ are smooth eigensections, and they have more than $m$ zeros on average, as is shown in the next section.

4. Gaussian random equivariant spherical harmonics

The space $\mathcal{H}_{N}^{m}$ has thus been identified as the space,

[TABLE]

The basis elements $z^{\alpha}\bar{z}^{\beta}$ are orthogonal on ${\mathbb{S}}^{3}$ but are not of norm $1$. To compute the norms it is advantageous to relate integrals over ${\mathbb{S}}^{3}$ with Gaussian integrals over ${\mathbb{C}}^{2}$, i.e. to use the measure $e^{-|Z|^{2}}dL(Z)$ where $dL$ is Lebesgue measure. The calculations are done in [Fol72, Page 98] and one finds that (using multi-index notation $z^{\alpha}=z_{1}^{\alpha_{1}}z_{2}^{\alpha_{2}}$),

[TABLE]

Henceforth we denote the orthonormalized monomials by $\hat{s}_{N,\alpha,\bar{\beta}}:=\frac{z^{\alpha}\bar{z}^{\beta}}{||z^{\alpha}\bar{z}^{\beta}||}$.

Definition 4.1**.**

The Gaussian random equivariant spherical harmonic $\psi_{N}^{m}\in\mathcal{H}_{N}^{m}=\mathcal{H}_{N}^{p,q},p-q=m$ is defined by the series

[TABLE]

where the coefficients $a_{\alpha,\beta}$ are independent complex normal Gaussians.

When $p=m=N$ these polynomials are known as the $SU(2)$ Gaussian holomorphic polynomials; for general $(p,q)$ we call them $SU(2)^{p,q}$ polynomials. They are sometimes called ‘poly-analytic’ functions. The subspace $\mathcal{H}_{N}^{0}$ consists of invariant eigenfunction pulled back from ${\mathbb{S}}^{2}$. They are real valued and so the equivariant nodal sets are not intersections of two real nodal sets.

4.1. Covariance kernel

The covariance kernel of the Gaussian random equivariant spherical harmonics is the orthogonal projection, $\Pi_{N}^{m}:L^{2}({\mathbb{S}}^{3})\to\mathcal{H}_{N}^{m}$, given by

[TABLE]

where $\Pi_{N}:L^{2}({\mathbb{S}}^{3})\to\mathcal{H}_{N}$ is the covariance kernel of the Gaussian random spherical harmonics.

Explicit calculations to follow are based on the identity,

[TABLE]

where

[TABLE]

is the Chebyshev polynomial of the second kind. (Note: Gegenbauer polynomial for $\dim=3$ is $U_{N}$.)

4.2. Real versus complex Gaussian ensembles

The purpose of this section is to show that that the real parts of complex Gaussian $SU(2)^{p,q}$ are real Gaussian spherical harmonics in the standard sense employed by Nazarov-Sodin, Sarnak etc. ( Taking the real part defines the map,

[TABLE]

Its image is a real subspace we denote by $\mathcal{H}_{N,m}\subset\mathcal{H}_{N}$. If we push forward the complex Gaussian measure $\gamma_{N}^{m}$ on $\mathcal{H}_{N}^{m}$ we get a Gaussian measure on $\mathcal{H}_{N,m}$. Our claim is that this measure coincides with the real Gaussian measure on $\mathcal{H}_{N}$ conditioned on $\mathcal{H}_{N,m}$. We denote by $E_{N}^{m}$ the expectation with respect to $\gamma_{N}^{m}$ and $E_{N,m}$ the conditional expectation of $\gamma_{N}$ conditioned on $\mathcal{H}_{N,m}$./

Lemma 4.2**.**

${\mathbb{E}}_{N}^{p,q}\mathcal{Z}_{\Re\psi_{N}^{m}}={\mathbb{E}}_{N,m}\mathcal{Z}_{\Re\psi_{N}^{m}}$, and ${\mathbb{E}}_{N}^{p,q}g(\mathcal{Z}_{\Re\psi_{N}^{m}})={\mathbb{E}}_{N,m}g(\mathcal{Z}_{\Re\psi_{N}^{m}})$.

Proof.

We assume $p+q=N,p-q=m$. The real part of a complex combination

[TABLE]

equals

[TABLE]

Here, $A_{p,q},B_{p,q}$ are independent $N(0,1)$ random variables. If we condition on $B_{\alpha,\beta}=0$ we get the conditional real Gaussian ensemble.

We consider the measure-valued random variable $\mathcal{Z}_{\Re\psi_{N}^{m}}$ or its Euler characteristic (or any functional of $\Re\psi_{N}^{m})$ as a function $F(A-B)$ where $A,B$ are independent $N(0,1)$ vectors. That is, the Gaussian measure on the coefficients $A,B$ is the product $d\gamma(A)d\gamma(B)$. The conditional Gaussian ensemble is $d\gamma(A)\delta_{0}(B)$. So the Lemma boils down to proving that

[TABLE]

Since $\gamma*\gamma=\gamma$, the left side equals

[TABLE]

as claimed. ∎

5. Kac-Rice formula

In view of Lemma 1.2, the proof of Theorem 1.1 is reduced to the calculation of ${\mathbb{E}}_{N,m}\#\{f_{N}^{m}=0\}$. To this end, we use the Kac-Rice formula. The Kac-Rice formula in the setting of random smooth sections of complex line bundles over Kähler manifolds is proved in [DSZ04]. It makes use of the canonical lift of sections to equivariant functions on the associated $S^{1}$ bundle, and is therefore well adapted to our setting. We briefly review the formula and then move on to the calculation of ${\mathbb{E}}_{N,m}\#\{f_{N}^{m}=0\}$.

5.1. Review of the Kac-Rice formula in our setting

We closely follow the exposition in [DSZ04, Section 5.4].

To state the Kac-Rice result, we need some notation and background. The pullback of the Dirac mass $\delta_{0}$ at $0\in{\mathbb{C}}$ is given by,

[TABLE]

Here, we use that the Jacobian $J_{f}$ of $f:{\mathbb{C}}\to{\mathbb{C}}$ is given by $J_{f}=|df\wedge d\bar{f}|$. The equivarlant lift of $f$ to the circle bundle is $fe^{-m\varphi/2}$ where $\varphi$ is the Kähler potential of the Fubini-Study metric. The pullback of $\delta_{0}$ under this complex-valued function is,

[TABLE]

since, $De^{-m\varphi/2i}f=e^{-m\varphi/2}Df$ at a zero.

Second, we need to recall the joint probability density $D_{N}^{m}(x,\xi;z)$ of the random variables

[TABLE]

and in particular the ‘conditional’ density $D(0,\xi;z)$. The joint probability density is given by,

[TABLE]

where $\Delta=\Delta_{N}^{m}(z)$ is the covariance matrix of $(X_{z},\Xi_{z})$,

[TABLE]

Here, $\nabla^{1}_{z}$, respectively $\nabla^{2}_{z}$, denotes the differential operator on $X\times X$ given by applying $\nabla_{z}$ to the first, respectively second, factor. For notational simplicity, we often drop the super- and sub-scripts $(N,m)$ in what follows.

As discussed in [DSZ04], and as is well-known, $D(0,\xi;z)$ is then given by,

[TABLE]

where

[TABLE]

is the Gaussian density with covariance matrix

[TABLE]

and where

[TABLE]

The formula (5.6) for $D(0,\xi;z)$ simplifies to,

[TABLE]

Proposition 5.1**.**

Let $s=fe$ in a local frame and let $\hat{s}=fe^{-m\varphi/2}$ . Then, ${\mathbb{E}}\mathcal{Z}_{s}$ is the measure on ${\mathbb{C}}{\mathbb{P}}^{1}$ given by

[TABLE]

where $dL$ is Lebesgue measure and where $D(x,\xi;z)$ is the joint probability density of $(f(z),df(z))$, given by (5.10).

Proof.

By definition,

[TABLE]

Here, $\left|d(f(z,\bar{z})e^{-n\varphi/2})\wedge d\bar{(}f(z,\bar{z})e^{-n\varphi/2})\right|$ is a density (the absolute value of a volume form). We then replace the $\delta_{0}(f)$ by the Fourier integral, to get

[TABLE]

In general, if $F$ is a complex Gaussian random field, and $\Phi$ is a (possibly nonlinear) functional,

[TABLE]

Therefore,

[TABLE]

Using that $\int_{{\mathbb{C}}}e^{i\Re x\bar{t}}dL(t)=\delta_{0}(x)$ we get

[TABLE]

It remains to compute $D(x,\xi;z)$, and we outline the calculation in [DSZ04] using the real linear 1-jet map $\mathcal{J}:=J^{1}_{z}$, which is locally written in terms of an orthonormal basis $\{f_{j}\}$ as

[TABLE]

We may regard $\mathcal{J}$ as a map from $a\in{\mathbb{C}}^{N}$ into $(x,\xi)\in{\mathbb{C}}\times{\mathbb{C}}^{2}.$ The joint probability density is the push forward of the measure $e^{-|a|^{2}/2}da$ under $J^{1}_{z}$,

[TABLE]

where $d\dot{a}$ is the surface Lebesgue measure on the subspace $\mathcal{J}^{-1}(x,\xi)$ . This follows from general principles on pushing forward complex Gaussians under complex linear maps $F:{\mathbb{C}}^{d}\to{\mathbb{C}}^{n}$, whereby

[TABLE]

As in [DSZ04], one shows that $\mathcal{J}\mathcal{J}^{*}=\Delta_{N}^{m}$ above. Conditioning on $x=0$ then gives (5.6).

∎

5.2. Symmetries and application to Lemma 1.3

Note that $\Pi_{N}(\Phi(x),\Phi(y))=\Pi_{N}(x,y)$ for any isometry $\Phi\in SO(4)$. However, this is not true for $\Pi_{N}^{m}(x,y)$. To understand the symmetries of $\Pi_{N}^{m}(x,y)$, we first identify $(z_{1},z_{2})\in\mathbb{C}^{2}$ in (2.1) with unit quaternion $p=z_{1}+z_{2}j$. Then the action $r_{\theta}$ (2.3) is equivalent to the left multiplication by $e^{i\theta}$.

We recall that the map from $SU(2)\times SU(2)$ to $SO(4)$ given by mapping a pair of unit quaternions $(p,q)$ to the map $\Phi_{p,q}:x\mapsto\bar{p}xq$ is a surjective homomorphism with the kernel $\{(1,1),(-1,-1)\}$. Observe that among $\Phi_{p,q}$, $\Phi_{1,q}$ for any fixed $q$ commutes with the action $r_{\theta}$, by the associativity of multiplication of quaternions. Hence $\Phi_{1,q}$ leaves $\Pi_{N}^{m}(x,y)$ invariant, because

[TABLE]

We also infer that $\Phi_{1,q}$ is well-defined on the fibers of the Hopf fibration $\pi:{\mathbb{S}}^{3}\to{\mathbb{S}}^{2}$, and therefore induces $SU(2)$-action on ${\mathbb{S}}^{2}$. This induces a surjective homomorphism $SU(2)\to SO(3)$ with the kernel $\{1,-1\}$, and it is well-known that $SO(3)$ acts doubly transitively on ${\mathbb{S}}^{2}$.

Now to prove Lemma 1.3, we need to calculate $\Delta_{N}^{m}$ and $\Lambda$ in (5.7). Equivalently, for a local orthonormal frame $E=\{\partial_{\theta},\mathbf{e}_{1},\mathbf{e}_{2}\}$ near $x\in{\mathbb{S}}^{3}$, we need to calculate

[TABLE]

and

[TABLE]

We then deduce from $\Phi_{1,q}$-invariance of $\Pi_{N}^{m}(x,y)$, and the discussion above that these quantities do not depend on the choice of $x\in{\mathbb{S}}^{3}$ and $E=\{\partial_{\theta},\mathbf{e}_{1},\mathbf{e}_{2}\}$:

Lemma 5.2**.**

$\Delta_{N}^{m}$* and $\Lambda$ in (5.7) are constant matrices.*

5.3. Chebyshev calculations

In this section, we compute $\Delta_{N}^{m}$ and $\Lambda$ in (5.7) explicitly. Firstly, from Lemma 5.2, it is sufficient to compute the matrix at $x=(\alpha,\varphi,\theta)=(\pi/4,0,0)$ in the coordinate system (2.2).

Note that $\Pi_{N}^{m}(x,y)=0$, if $2\nmid N-m$ or if $|m|>N$, by Proposition 2.3. So we assume that $2|N-m$ and $|m|\leq N$ in this section.

Lemma 5.3**.**

Let $U_{N}(x)$ be the Chebyshev polynomial of the second kind. Assume that $m\in\mathbb{Z}$ and $N\in\mathbb{N}$ satisfies $|m|\leq N$.

Then we have

[TABLE]

Proof.

Recall that

[TABLE]

Therefore the first integral is $1$ if $m$ is equal to one of $N,N-2,\ldots,-N$, and [math] otherwise. To compute the second integral, we first differentiate the above equation to get:

[TABLE]

and therefore

[TABLE]

Because the second integral does not depend on the sign of $m$, we assume without loss of generality that $m\geq 0$. Since $\cos(m\theta)\cos\theta=\frac{1}{2}\left(\cos((m+1)\theta)+\cos((m-1)\theta)\right)$, we have

[TABLE]

when $2|N-m$, and [math] otherwise. ∎

Theorem 5.4**.**

We have

[TABLE]

and

[TABLE]

Proof.

Firstly, we have

[TABLE]

by Lemma 5.3. For $\nu=\alpha$ or $\varphi$, we have

[TABLE]

If we write

[TABLE]

and

[TABLE]

then

[TABLE]

which is [math], and

[TABLE]

which simplifies to

[TABLE]

Therefore

[TABLE]

by Lemma 5.3. Likewise,

[TABLE]

and we have

[TABLE]

Therefore

[TABLE]

For the last case,

[TABLE]

Finally, recall that

[TABLE]

where

[TABLE]

and

[TABLE]

where $\{\partial_{\theta},\mathbf{e}_{1},\mathbf{e}_{2}\}$ is a local orthonormal frame. (So $A$, $B$, and $C$ are $1\times 1$, $1\times 2$, and $2\times 2$ complex matrices.) Above computation with (\mathbf{e}_{1},\mathbf{e}_{2})=(\partial_{\theta},\partial_{\varphi})\big{|}_{(\alpha,\varphi)=(\pi/4,0)} imply that

[TABLE]

∎

Proof of Lemma 1.3.

From Proposition 5.1 and Theorem 5.4, we see that

[TABLE]

where

[TABLE]

By change of variables $\Lambda^{-1/2}\xi=\zeta$, we have

[TABLE]

Now let

[TABLE]

and let $v_{1}=\begin{pmatrix}1\\ i\end{pmatrix}$ and $v_{2}=\begin{pmatrix}1\\ -i\end{pmatrix}$. Observe that $\Lambda v_{1}=(\mu+\nu)v_{1}$, and $\Lambda v_{2}=(\mu-\nu)v_{2}$ with $\mu=\frac{1}{N+1}\left(\frac{N^{2}-m^{2}}{2}+N\right)$ and $\nu=\frac{m}{2(N+1)}$. In particular, we have

[TABLE]

where

[TABLE]

Therefore

[TABLE]

and we evaluate the integral as follows:

[TABLE]

In the last equality, we used the fact that $|\eta|\leq\frac{1}{2}$ to evaluate the integral. ∎

5.4. Bertini theorem

We need the following Bertini-type theorem to employ the Kac-Rice formalism and to prove Theorem 1.1:

Proposition 5.5**.**

For all $m\not=0$, [math] is almost surely a regular value of the random equivariant eigenfunction $\psi_{N}^{m}:{\mathbb{S}}^{3}\to{\mathbb{C}}$.

Proof.

We need to show that the derivative $d_{x}\psi_{N}^{m}:T_{x}{\mathbb{S}}^{3}\to{\mathbb{C}}$ is surjective at each point $x$ where $\psi_{N}^{m}(x)=0$ for almost any $\psi_{N}^{m}$. It is sufficient to prove that $\mathcal{H}_{N}^{m}$ has the “1-jet spanning property” that the 1-jet evaluation map,

[TABLE]

is surjective.

The $1$-jet spanning property implies that at each point, $\{d_{x}^{H}\psi_{N}^{m}(x):\psi_{N}^{m}\in\mathcal{H}_{N}^{m}\}$ spans the horizontal tangent space $H_{x}{\mathbb{S}}^{3}$. Since the ensemble is $SU(2)$ invariant, it suffices to prove the spanning property at a single point. Moreover, since $SU(2)$ acts transitively on the unit tangent bundle of ${\mathbb{S}}^{2}$ (or on the horizontal spaces of ${\mathbb{S}}^{3}$), failure to span is equivalent to the existence of $x$ such that $d_{x}^{H}\psi_{N}^{m}(x)=0$ for all $\psi_{N}^{m}$ such that $\psi_{N}^{m}(x)=0$. This is false, since Jacobi polynomials have simple zeros, as can for instance be seen from the Darboux formula.

As explained in [BSZ01, Section 4.1], the $1$-jet spanning property implies that the incidence set $I:=\{(x,\psi_{N}^{m})\in{\mathbb{S}}^{3}\times\mathcal{H}_{N}^{m}:\psi_{N}^{m}(x)=0\}$ is a smooth submanifold, and hence by Sard’s theorem applied to the projection $I\to\mathcal{H}_{N}^{m}$, the zero set

[TABLE]

is a smooth $1$-dimensional submanifold of ${\mathbb{S}}^{3}$ for almost all $\psi_{N}^{m}$. ∎

6. Résumé of the proof Theorem 1.1

Having established all of the ingredients of the proof of Theorem 1.1 outlined in the Introduction, we only review how to assemble the ingredients into a proof. We resume the discussion begun in Section 2.1. To prove (i) of Theorem 1.1 we use the proof of the same statement as Theorem 2.1 ([JZ18, Theorem 1.5]) for general Kaluza-Klein Laplacians on circle bundles. The only point we need to establish to apply the proof is that [math] is almost surely a regular value of the equivariant eigenfunctions, and this is proved in Theorem 5.5. The main new result is therefore (ii) of Theorem 1.1. It is based on a simple formula of Lemma 1.2 relating the genus of the nodal set with the number of zeros of $f_{N}^{m}$; this makes use of the special structure of nodal set of real parts of equivariant eigenfunctions as ‘helicoid covers’ of ${\mathbb{S}}^{2}$. In Section 5, we use the Kac-Rice to calculate the expected number of zeros of the random $f_{N}^{m}$ and prove Lemma 1.3, concluding the proof of Theorem 1.1.