Hostname: page-component-69cd664f8f-lpdgq Total loading time: 0 Render date: 2025-03-13T05:39:48.133Z Has data issue: false hasContentIssue false

GAUSSIAN HOLOMORPHIC SECTIONS ON NONCOMPACT COMPLEX MANIFOLDS

Published online by Cambridge University Press:  12 March 2025

Alexander Drewitz
Affiliation:
Universität zu Köln, Department Mathematik/Informatik, Weyertal 86–90, 50931 Köln, Germany (adrewitz@uni-koeln.de gmarines@math.uni-koeln.de)
Bingxiao Liu*
Affiliation:
Universität zu Köln, Department Mathematik/Informatik, Weyertal 86–90, 50931 Köln, Germany (adrewitz@uni-koeln.de gmarines@math.uni-koeln.de)
George Marinescu
Affiliation:
Universität zu Köln, Department Mathematik/Informatik, Weyertal 86–90, 50931 Köln, Germany (adrewitz@uni-koeln.de gmarines@math.uni-koeln.de)
Rights & Permissions [Opens in a new window]

Abstract

We provide two constructions of Gaussian random holomorphic sections of a Hermitian holomorphic line bundle $(L,h_{L})$ on a Hermitian complex manifold $(X,\Theta )$, that are particularly interesting in the case where the space of $\mathcal {L}^2$-holomorphic sections $H^{0}_{(2)}(X,L)$ is infinite dimensional. We first provide a general construction of Gaussian random holomorphic sections of L, which, if $H^{0}_{(2)}(X,L)$ is infinite dimensional, are almost never $\mathcal {L}^2$-integrable on X. The second construction combines the abstract Wiener space theory with the Berezin–Toeplitz quantization and yields a Gaussian ensemble of random $\mathcal {L}^2$-holomorphic sections. Furthermore, we study their random zeros in the context of semiclassical limits, including their distributions, large deviation estimates, local fluctuations and hole probabilities.

Type
Research Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2025. Published by Cambridge University Press

1 Introduction

The origins of the study of random polynomials as well as their generalizations in geometry trace back to value distribution theory of analytic functions and quantum chaos, the former being a classic topic. In order to understand the typical distribution of zeros, one forms ensembles of analytic functions by defining the coefficients of their representation in a given basis to be independent random variables with a prescribed distribution and then studies their statistical properties. Importantly, the eigenfunctions of quantum chaotic Hamiltonians are well modeled by random polynomials in terms of the distribution of their zeros, critical points, sup-norms and quantum expectation values (cf. [Reference Bogomolny, Bohigas and Leboeuf11, Reference Nonnenmacher and Voros62, Reference Zelditch79, Reference Shiffman and Zelditch68]).

A natural geometric generalization of polynomials are the holomorphic sections of a holomorphic line bundle. In this paper we study the construction of Gaussian ensembles of $\mathcal {L}^{2}$ -holomorphic sections of a line bundle, focusing on the case when this space is infinite dimensional. We consider a connected complex n-dimensional Hermitian manifold $(X,\Theta )$ without boundary, where $\Theta $ is a smooth Hermitian form compatible with the complex structure of X. We denote by L a holomorphic line bundle over X, and let $h_{L}$ be a smooth Hermitian metric on L. We denote by $H^0(X,L)$ the space of holomorphic sections of L and by $H^0_{(2)}(X,L)$ the space holomorphic sections which are $\mathcal {L}^{2}$ -integrable with respect to the $\mathcal {L}^{2}$ -norm (2.2) induced by $\Theta $ and $h_L$ . If the dimension $d:=\dim _{\mathbb {C}} H^0_{(2)}(X,L)$ is finite, one can construct a Gaussian probability measure on $H^0_{(2)}(X,L)$ induced by the $\mathcal {L}^{2}$ -metric. If, however, $H^0_{(2)}(X,L)$ is infinite dimensional, then this construction is bound to fail.

To tackle this difficulty, we will provide two different approaches for constructing a random holomorphic section from the infinite-dimensional $H^0_{(2)}(X,L)$ . Both of them are natural as extensions of the finite-dimensional case.

The first approach is a direct generalization of the study of random holomorphic functions on $\mathbb {C}^{n}$ to the context of complex geometry. The random holomorphic functions given by power series on $\mathbb {C}$ as well as the distributions of their zeros (or other values) have been studied by Kac [Reference Kac45], Littlewood–Offord [Reference Littlewood and Offord52, Reference Littlewood and Offord53], Offord [Reference Offord63, Reference Offord64, Reference Offord65] and by Edelman–Kostlan [Reference Edelman and Kostlan29, Reference Edelman and Kostlan30], etc. For Gaussian random holomorphic functions, the results have further been extended by Sodin [Reference Sodin71], Sodin–Tsirelson [Reference Sodin and Tsirelson72, Reference Sodin and Tsirelson73, Reference Sodin and Tsirelson74], and then, on $\mathbb {C}^{n}$ , by Zrebiec [Reference Zrebiec82]. In particular, the general Gaussian random holomorphic functions on the domains in $\mathbb {C}$ (also known under the name Gaussian analytic functions, GAFs) have been investigated vastly (cf. [Reference Hough, Krishnapur, Peres and Virág40]) from probabilistic perspectives, serving as examples of point processes on $\mathbb {C}$ .

We construct a Gaussian random section in terms of an orthonormal basis of $H^{0}_{(2)}(X,L)$ such that its distribution is independent of the choice of such basis (cf. Proposition 2.3). More concretely, if $\{S_{j}\}^{d}_{j=1}$ is an orthonormal (Hilbert) basis of $H^{0}_{(2)}(X,L)$ and if $\{\eta _{j}\}_{j=1}^{d}$ denotes a sequence of independent and identically distributed (i.i.d.) standard complex Gaussian variables, we can define a random holomorphic section of L via

(1.1) $$ \begin{align} \psi^{S}_{\eta}:=\sum_{j=1}^{d}\eta_{j}S_{j} \end{align} $$

by using properties of the Bergman kernel associated with $H^{0}_{(2)}(X,L)$ (cf. Proposition 2.1). We will call $\psi ^{S}_{\eta }$ a standard Gaussian random holomorphic section of L. It turns out that when $d=\infty $ , then $\psi ^{S}_{\eta }$ as constructed in Equation (1.1) is almost surely non- $\mathcal {L}^{2}$ -integrable over X (cf. Lemma 2.5). In the case of the Bargmann–Fock space on $\mathbb {C}^{n}$ (cf. Example 2.14), $\psi ^{S}_{\eta }$ is just a Gaussian holomorphic function on $\mathbb {C}^{n}$ as mentioned before. If $d<\infty $ the above construction is equivalent to endowing $H^{0}_{(2)}(X,L)$ with the standard Gaussian probability measure associated to the $\mathcal {L}^{2}$ inner product.

An interesting question which arises is what are the possibilities to randomize $\mathcal {L}^{2}$ -holomorphic sections in a natural manner, or equivalently, how to construct Gaussian probability measures on $H^{0}_{(2)}(X,L)$ in a geometric way. The starting point is the simple observation that for $\mathbf {a}=(a_{j})_{j=1}^{d}\in \ell ^{2}(\mathbb {C})$ , the section

(1.2) $$ \begin{align} \psi^{S}_{\mathbf{a},\eta}:= \sum_{j=1}^{d}\eta_{j}a_{j}S_{j}\, \end{align} $$

is almost surely an $\mathcal {L}^{2}$ -integrable holomorphic section since $\mathbb {P}(\sum _{j}|a_{j}|^{2}|\eta _{j}|^{2}<\infty )=1$ . We introduce the geometric input by taking $(a_{j})_{j=1}^{d}$ to be the spectrum of a Toeplitz operator on $H^{0}_{(2)}(X,L)$ . For any bounded function $f\in \mathscr {C}^\infty (X)$ , the Toeplitz operator with symbol f is the endomorphism of $H^{0}_{(2)}(X,L)$ given by $T_f=PM_f$ , where $M_f$ is the pointwise multiplication with f and P is the orthogonal projection from $\mathcal {L}^{2}(X,L)$ onto $H^{0}_{(2)}(X,L)$ . For our purpose, we consider an injective Hilbert–Schmidt Toeplitz operator $T_{f}$ on $H^{0}_{(2)}(X,L)$ , which defines a measurable norm $\|T_{f}\cdot \|$ on $H^{0}_{(2)}(X,L)$ (cf. Definition 4.1). Taking advantage of the theory of abstract Wiener spaces of Gross [Reference Gross36], we can construct in a unique way a Gaussian probability measure $\mathbb {P}_{f}$ on $H^{0}_{(2)}(X,L)$ associated with $T_{f}$ . A random $\mathcal {L}^{2}$ -holomorphic section following the probability law $\mathbb {P}_{f}$ is exactly given as in Equation (1.2), where each $a_{j}>0$ is an eigenvalue of $T_{f}$ and the orthonormal basis $\{S_{j}\}_{j=1}^{d}$ is such that $T_{f}S_{j}=a_{j}S_{j}$ . For a brief introduction to Gross’ abstract Wiener spaces, we refer to [Reference Janson42, Example 1.25].

On top of the constructions of random holomorphic sections outlined above, we aim to study the distribution of their zeros as $(1,1)$ -currents on X in the framework of semiclassical limits, that is, considering random holomorphic sections of the sequence of high tensor powers $(L^{p},h_{p}):=(L^{\otimes p},h_{L}^{\otimes p})$ , $p\in \mathbb {N}$ , of a given positive Hermitian line bundle $(L,h_{L})$ . As $p\rightarrow \infty $ , the parameter $h:=1/p$ , playing the role of the Planck constant, tends to $0$ . For this purpose, we need to make further assumptions on $(X,\Theta )$ and $(L,h_{L})$ ; see Condition 1.2. For a sequence of random sections $\psi ^{S_{p}}_{\eta }$ constructed as in Equation (1.1) using an orthonormal basis of the Hilbert spaces $H^{0}_{(2)}(X,L^{p})$ , the distributions of the normalized zeros of these sections are expected to converge towards the first Chern form $c_{1}(L,h_{L})$ as $p\rightarrow \infty $ . Thus, zeros tend to concentrate in regions of high curvature.

For compact Kähler manifolds, the asymptotic distribution as $p\rightarrow \infty $ of zeros of sections of $H^0(X,L^p)$ for a positive line bundle L was extensively studied by Shiffman and Zelditch [Reference Shiffman and Zelditch68, Reference Zelditch81]. Pioneering works in the physics literature were Bogomolny, Bohigas and Leboeuf [Reference Bogomolny, Bohigas and Leboeuf11] as well as Nonnenmacher–Voros [Reference Nonnenmacher and Voros62] in the context of quantum chaotic dynamics. A key ingredient in their approach is the asymptotic expansion of the associated Bergman kernel (cf. [Reference Ma and Marinescu56, Reference Tian76, Reference Zelditch80] and the references therein).

Using ideas from complex dynamics, a new method to study the distribution of zeros was introduced by Dinh–Sibony [Reference Dinh and Sibony26], which also provides bounds for the convergence speed. Subsequently, Dinh, Marinescu and Schmidt [Reference Dinh, Marinescu and Schmidt25] extended such results to the noncompact setting, where they assumed that the dimension $d_{p}:=\dim _{\mathbb {C}}H^{0}_{(2)}(X,L^{p})\in \mathbb {N}\cup \{\infty \}$ satisfies $d_{p}=\mathcal {O}(p^{n})$ for $p\gg 0$ (which is the same growth as in the compact case). Along these lines, there are also plenty of generalizations to different geometric (singular Hermitian metrics or singular space X) or probabilistic settings (general probability measures), cf. [Reference Bayraktar, Coman and Marinescu6, Reference Bloom and Levenberg10, Reference Coman and Marinescu15, Reference Coman, Ma and Marinescu13, Reference Coman, Marinescu and Nguyên16, Reference Coman, Marinescu and Nguyên17, Reference Dinh, Ma and Marinescu24, Reference Marinescu and Vu61]. We refer to the survey papers [Reference Bayraktar, Coman, Herrman and Marinescu5, Reference Zelditch81] for more details and references on this topic.

For compact Kähler manifolds, Shiffman, Zelditch and Zrebiec [Reference Shiffman, Zelditch and Zrebiec70] established large deviation estimates for the random zeros of Gaussian holomorphic sections as p grows to infinity. As a consequence, they obtained the exponential decay of the hole probabilities, that is, the probabilities that the Gaussian random holomorphic sections do not vanish on a given domain in X. In our previous paper [Reference Drewitz, Liu and Marinescu27] (see also [Reference Liu and Zielinski54]), we generalized their results to the noncompact setting, especially for the case of Riemann surfaces with cusps, under the assumption $d_{p}=\mathcal {O}(p^{n})$ for $p\gg 0$ . In this paper, we show that these results also extend to the case of Gaussian random holomorphic sections $\psi ^{S_{p}}_{\eta }$ without assuming $d_{p}<\infty $ .

In this semiclassical setting, the use of Toeplitz operators in our construction of the random $\mathcal {L}^{2}$ -holomorphic sections is natural. The family of Toeplitz operators $T_{f,p}=P_pM_f$ on $H^{0}_{(2)}(X,L^{p})$ , $p\in \mathbb {N}$ , is called Berezin–Toeplitz quantization of a given real smooth function $f:X\rightarrow \mathbb {R}$ . The function f can be regarded as a classical observable on the phase space $(X,\Theta )$ and $T_{f,p}$ as the corresponding quantum observable on the quantum space $H^{0}_{(2)}(X,L^{p})$ (cf. [Reference Bordemann, Meinrenken and Schlichenmaier12]). Such operators are central in the study of geometric quantization on Kähler or, more generally, on symplectic manifolds. For more details, we refer to [Reference Bordemann, Meinrenken and Schlichenmaier12, Reference Ma55, Reference Ma and Marinescu57, Reference Ma and Marinescu58, Reference Ma and Marinescu59] and [Reference Ma and Marinescu56, Chapter 7]. In our context, the geometry of $(X,\Theta )$ and $(L,h_L)$ enters the picture through the sequence of Bergman projectors $P_p$ on $H^{0}_{(2)}(X,L^{p})$ and the multiplication with f will have the effect localizing the Gaussian probability measure $\mathbb {P}_{f}$ on the support of f.

Here, we introduce a class of functions f on X such that $T_{f,p}$ is Hilbert–Schmidt for all $p\gg 0$ . Associated to a nonnegative smooth function f in this class, we construct canonically a sequence of probability spaces $(H^{0}_{(2)}(X,L^{p}),\mathbb {P}_{f,p})$ , $p\gg 0$ . We are then interested in investigating the asymptotic behavior of the zeros of random $\mathcal {L}^{2}$ -holomorphic sections as $p\rightarrow \infty $ . Their limit as $(1,1)$ -currents will be given by $c_{1}(L,h_{L})$ on the support of f (where f has zeros of at most second order). One goal of geometric quantization is to recover information about the classical observable f from the action of the quantum observable $T_{f,p}$ in the semiclassical limit $p\to \infty $ . We show here that we can recover the complex Laplacian $\partial \overline \partial \log f^2$ of the logarithm of $f^2$ (see also Ancona–Le Floch [Reference Ancona and Le Floch1] for the compact case in a slightly different setting).

The next level of results concerns the universality of the statistics of zeros of sections of large degree p on small length scales of order the Planck scale $1/\sqrt {p}$ (cf. also [Reference Bleher, Shiffman and Zelditch8, Reference Hannay37, Reference Tao and Vu75]). Note that the Bergman and Toeplitz kernels exhibit their universality at this scale (see Theorem 3.14 and [Reference Ma and Marinescu56, Theorems 4.2.1 & 6.1.1]). By zooming in a small ball $B(x,R/\sqrt {p})$ around an arbitrary point x, one loses track of the special features of the geometrical setting and obtains universal limiting fluctuations of random zeros. We observe two different behaviors depending on whether $f(x)=0$ (with vanishing order $2$ ) or $f(x)\neq 0$ ; see Equation (1.25). A further interesting question is to describe the asymptotic distribution of random zeros outside the support of f.

Our approach to the above results relies on the asymptotic expansion of the on-diagonal Schwartz kernel of the operator $T^{2}_{f,p}=T_{f,p}\circ T_{f,p}$ , as $p\rightarrow \infty $ , whose first terms are computed explicitly in [Reference Ma and Marinescu56, Chapter 7] and in [Reference Ma and Marinescu59], and which also holds on noncompact manifolds cf. [Reference Ma and Marinescu56, Reference Ma and Marinescu57, Reference Ma and Marinescu58, Reference Ma and Marinescu59]. In particular, we can apply them to the case considered in [Reference Ma and Marinescu56, Section 7.5] and the case of bounded geometry discussed in [Reference Ma and Marinescu60] and [Reference Finski33].

A comprehensive study of zeros of Gaussian $\mathcal {L}^{2}$ -holomorphic sections by means of Toeplitz operators associated with a nonsmooth symbol f is given in our follow-up paper [Reference Drewitz, Liu and Marinescu28].

In the next four sections, we provide the setting and formulate our main results.

1.1 Zeros of Gaussian random holomorphic sections

Let us start with a (paracompact) connected complex manifold $(X,J,\Theta )$ , where J denotes the complex structure of X and $\Theta $ is a J-compatible smooth Hermitian form on X. Then $g^{TX}(\cdot ,\cdot )=\Theta (\cdot ,J\cdot )$ is the associated Riemannian metric on $TX$ . In this case, we refer to the triplet $(X,J,\Theta )$ simply as a complex Hermitian manifold. Let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ . Assume that we have $d=\dim _{\mathbb {C}} H^0_{(2)}(X,L)\geq 1$ .

The zero-sets of holomorphic sections are complex analytic sets, and they will be studied through the current of integration they define. We refer to Section 2.3 for basic definitions concerning the currents of integration along zero-divisors of a holomorphic section.

Our first result concerns the expectation of the currents of integration on the zero-divisors of the Gaussian random holomorphic section $\psi ^{S}_{\eta }$ defined in Equation (1.1), as a current on X, that is, of the random $(1,1)$ -current $[\operatorname {\mathrm {Div}}(\psi ^{S}_{\eta })]$ . For any test form $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ , the random variable $\langle [\operatorname {\mathrm {Div}}(\psi ^{S}_{\eta })],\varphi \rangle $ is measurable (cf. [Reference Coman and Marinescu15, proof of Proposition 4.2]). If the random variable $\langle [\operatorname {\mathrm {Div}}(\psi ^{S}_{\eta })],\varphi \rangle $ is integrable for any test form $\varphi $ , then the linear map

$$ \begin{align*}\varphi\mapsto \mathbb{E}\left[\langle [\operatorname{\mathrm{Div}}(\psi^{S}_{\eta})],\varphi\rangle\right]\,,\quad \varphi\in\Omega^{n-1,n-1}_{0}(X),\end{align*} $$

defines a $(1,1)$ -current on X, which we refer to as the expectation of $[\operatorname {\mathrm {Div}}(\psi ^{S}_{\eta })]$ , denoted by $\mathbb {E}[[\operatorname {\mathrm {Div}}(\psi ^{S}_{\eta })]]$ . Next, we define the Fubini–Study current $\gamma (L,h_L)$ on X. For this purpose, let

(1.3) $$ \begin{align} P:\mathcal{L}^2(X,L)\to H_{(2)}^0(X,L) \end{align} $$

be the $\mathcal {L}^{2}$ -orthogonal projection, called the Bergman projection. It has a smooth Schwartz kernel $P(x,y)$ , called the Bergman kernel, cf. Subsection 2.1. The Bergman kernel function $X\ni x\mapsto P(x,x)$ is a nonnegative smooth function on X, and the function $\log P(x,x)$ is locally integrable on X. We set

(1.4) $$ \begin{align} \gamma(L,h_{L})=c_{1}(L,h_{L})+\frac{\sqrt{-1}}{2\pi} \partial\bar{\partial}\log P(x,x), \end{align} $$

where $c_{1}(L,h_{L})$ is the Chern form of $(L,h_{L})$ . See Remark 2.7 for their geometric interpretation.

Theorem 1.1. Assume that $d\geq 1$ . Then the expectation of the random variable $[\operatorname {\mathrm {Div}}(\psi ^{S}_{\eta })]$ exists as a $(1,1)$ -current on X and we have in the sense of currents,

(1.5) $$ \begin{align} \mathbb{E}[[\operatorname{\mathrm{Div}}(\psi^{S}_{\eta})]]=\gamma(L,h_{L})\,. \end{align} $$

In the case $d<\infty $ , Equation (1.5) was already known for line bundles with empty base locus (cf. [Reference Shiffman and Zelditch68, Lemma 3.1], [Reference Ma and Marinescu56, Theorem 5.3.1]) and in several situations when the metric $h_L$ or the base X are singular (see, e.g., [Reference Coman and Marinescu15, Proposition 4.2], [Reference Coman and Marinescu14, Theorem 1.4]). When $d=\infty $ analogs of this result are known in the context of random holomorphic functions on $\mathbb {C}^{m};$ for instance, Edelman and Kostlan [Reference Edelman and Kostlan29, Sections 7 & 8] studied the expectations of complex zeros of random power series (in their paper, they mainly aimed to study the distribution of real zeros). Other interesting examples from complex geometry, where our Theorem 1.1 applies, are given in Subsection 2.5. In particular, Theorem 2.13 gives a more general version of Theorem 1.1 where we allow the Hermitian metric $h_{L}$ to be singular. The equidistribution results for the case of singular Hermitian metric $h_{L}$ with finite-dimensional holomorphic sections were studied in [Reference Coman, Ma and Marinescu13, Reference Coman and Marinescu14, Reference Coman and Marinescu15, Reference Coman, Marinescu and Nguyên16, Reference Coman, Marinescu and Nguyên17, Reference Dinh, Ma and Marinescu24].

1.2 High tensor powers of L: equidistribution and large deviations

We are interested in the semiclassical limit of the zeros of Gaussian holomorphic sections when we replace L by its high tensor powers. For this purpose, we need to make further assumptions on the complex Hermitian manifold $(X,J,\Theta )$ and on $(L,h_{L})$ as follows.

Condition 1.2. The Riemannian metric $g^{TX}$ induced by $\Theta $ is complete and there exist $C,C_{0},\varepsilon>0$ such that on X,

(1.6) $$ \begin{align} \sqrt{-1}R^{L}\geq \varepsilon \Theta,\;\; \sqrt{-1}R^{\det}\geq -C_{0} \Theta,\;\;|\partial \Theta|_{g^{TX}}\leq C\,, \end{align} $$

where $R^{\det }$ be the curvature of the holomorphic connection $\nabla ^{\det }$ on $K_X^*=\det (T^{(1,0)}(X))$ .

If Condition 1.2 holds, the Bergman kernel functions $P_{p}(x,x)$ have an asymptotic expansion in the tensor power p, which is uniform on any given compact subset of X, by [Reference Ma and Marinescu56, Chapter 6]. As a consequence, we have the convergence of currents

(1.7) $$ \begin{align} \frac{1}{p}\gamma(L^{p},h_{p})\rightarrow c_{1}(L,h_{L}) \quad \text{ as }p\rightarrow \infty. \end{align} $$

In the following, we denote by $\psi ^{S_{p}}_{\eta }$ the Gaussian random holomorphic section (as in Equation (1.1)) constructed from an orthonormal basis $S_{p}=\{S^{p}_{j}\}_{j=1}^{d_{p}}$ of $H^{0}_{(2)}(X,L^{p})$ . As is natural, before formulating our concentration estimates, we begin with stating findings for the limit of the expectations $\mathbb {E}[[\operatorname {\mathrm {Div}}(\psi ^{S_{p}}_{\eta })]]$ . While the results are novel in our specific setting and formulated precisely in Theorems 3.1 and 3.7 below, we roughly speaking prove the following results.

Theorem 1.3. Under the above settings, we have, as $p\rightarrow \infty $ ,

$$ \begin{align*}\frac{1}{p}\mathbb{E}[[\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]]\rightarrow c_{1}(L,h_{L}).\end{align*} $$

Moreover, for each $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ , we have that

(1.8) $$ \begin{align} \mathbb{P}\Big(\lim_{p\rightarrow \infty} \frac{1}{p}\langle [\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})],\varphi\rangle= \langle c_{1}(L,h_{L}),\varphi\rangle\Big)=1. \end{align} $$

It is clear that the first point is a consequence of Theorem 1.1 in combination with the limit (1.7). Note that in Theorem 3.1, we actually prove that, on any given compact subset of X, $\frac {1}{p}\mathbb {E}[[\operatorname {\mathrm {Div}}(\psi ^{S_{p}}_{\eta })]]$ admits an asymptotic expansion as $p\rightarrow \infty $ . By our proof of equality (1.8) (see Theorem 3.7), the same conclusion holds for any twice differentiable $(n-1,n-1)$ -form $\varphi $ with compact support, that is, the limit above holds in the sense of currents of order $2$ .

A deeper question on the distribution of zeros is whether sequences of individual sections tend to have equidistributed zeros, that is, if for a random sequence of sections $\psi ^{S_{p}}_{\eta }$ , the normalized currents of integration $p^{-1}[\operatorname {\mathrm {Div}}(\psi ^{S_{p}}_{\eta })]$ converge to $c_{1}(L,h_{L})$ as $p\to \infty $ . We show in Corollary 3.8, that under natural finiteness hypotheses (of the mass of the current $c_{1}(L,h_{L})$ or of the volume of X in the metric given by $c_{1}(L,h_{L})$ ) we can improve the equality (1.8) to

(1.9) $$ \begin{align} \mathbb{P}\Big(\lim_{p\rightarrow \infty} \frac{1}{p} [\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]= c_{1}(L,h_{L})\Big)=1, \end{align} $$

where the limit is taken with respect to the weak convergence of $(1,1)$ -currents on X. In fact, by Remark 3.9 (3), the limit in the equality (1.9) can be taken in the sense of weak convergence of currents of order $0$ , that is, continuous linear functionals on the space of continuous $(n-1,n-1)$ -forms with compact support.

With these distribution results about zeros of random sections at our disposal, a natural next step is to investigate the speed of convergence in terms of large deviation estimates as in [Reference Drewitz, Liu and Marinescu27, Reference Shiffman, Zelditch and Zrebiec70] by allowing the possibility that $d_{p}=\infty $ .

Theorem 1.4. Let $(X,J,\Theta )$ be a connected complex Hermitian manifold, and let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ such that Condition 1.2 holds. For any $\delta>0$ and $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ , there exists a constant $c=c(\delta ,\varphi )>0$ such that for $p\in \mathbb {N}$ , we have

(1.10) $$ \begin{align} \mathbb{P}\Big( \ \Big|\Big\langle\frac{1}{p}[\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})] - c_{1}(L,h_L),\varphi\Big\rangle\Big|>\delta\ \Big)\leq \exp(-c\,p^{n+1}). \end{align} $$

Estimates like (1.10) are called large deviation estimates or concentration of measure estimates [Reference Ledoux50]. Moreover, the proof of the above theorem is independent of equality (1.8); in fact, we can prove the equality (1.8) by using the estimates (1.10) (cf. [Reference Drewitz, Liu and Marinescu27, Section 3.6]).

A natural follow-up question is whether a central limit theorem for the distribution of zeros of $\psi ^{S_{p}}_{\eta }$ as $p\rightarrow \infty $ holds true; we will touch upon this question in Remark 3.17.

Since $c_{1}(L,h_{L})$ is positive, $\frac {1}{n!}c_{1}(L,h_{L})^{n}$ also defines a positive volume element on X, and for $U\subset X$ open we set

(1.11) $$ \begin{align} \operatorname{\mathrm{Vol}}^{L}_{2n}(U)=\int_{U} \frac{1}{n!}{c_{1}(L,h_{L})^{n}}. \end{align} $$

Then for $s_{p}\in H^{0}(X,L^{p})\setminus \{0\}$ we define the $(2n-2)$ -dimensional volume with respect to $c_{1}(L,h_{L})$ of the divisor $\operatorname {\mathrm {Div}}(s_{p})$ (cf. (2.27)) in an open subset $U\subset X$ as

(1.12) $$ \begin{align} \operatorname{\mathrm{Vol}}^{L}_{2n-2}\big(\operatorname{\mathrm{Div}}(s_{p})\cap U\big)=\sum_{V\subset Z(s_{p})}\mathrm{ord}_{V}(s_{p})\int_{V\cap U} \frac{c_{1}(L,h_{L})^{n-1}}{(n-1)!}\,\cdot \end{align} $$

If we use this volume to measure the size of the zeros of $s_{p}$ in U, then Theorem 1.4 leads to the following result.

Theorem 1.5. Let $(X,J,\Theta )$ be a connected complex Hermitian manifold, and let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ . We assume that Condition 1.2 holds. If U is a nonempty relatively compact open subset of X such that $\partial U$ has zero measure in X, then for any $\delta>0$ , there exists a constant $c_{U,\delta }>0$ such that for all p large enough, we have

(1.13) $$ \begin{align} \mathbb{P}\Big( \ \Big|\frac{1}{p}\operatorname{\mathrm{Vol}}^{L}_{2n-2}(\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})\cap U) - n\operatorname{\mathrm{Vol}}^{L}_{2n}(U) \Big|>\delta \ \Big)\leq \exp(-c_{U,\delta}\, p^{n+1})\,. \end{align} $$

In addition, there exists a constant $C_{U}>0$ such that for $p> 0$ ,

(1.14) $$ \begin{align} \mathbb{P}\big(\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})\cap U=\varnothing\big) \leq \exp(-C_{U}p^{n+1})\,. \end{align} $$

The proofs of the above two theorems will be provided in Subsection 3.2. One essential ingredient for these proofs is Proposition 3.11, for which we need a more refined investigation of the local sup-norms of holomorphic sections on X (cf. Subsection 3.3).

The probability in estimate (1.14) is referred to as the hole probability of the random section $\psi ^{S_{p}}_{\eta }$ on the subset U. This estimate then provides us with an upper bound for the hole probabilities for $p> 0$ . In [Reference Shiffman, Zelditch and Zrebiec70, Theorem 1.4] and [Reference Drewitz, Liu and Marinescu27, Proposition 1.7], under additional assumptions on U, a lower bound of the form $\exp (-C^{\prime }_{U}p^{n+1})$ for the hole probabilities was proved. In general, however, such a lower bound remains unclear in the case $d_{p}=\infty $ .

Let us mention here that in the special case of the Bargmann–Fock space, for the standard Gaussian random holomorphic function on $\mathbb {C}^{n}$ (cf. Equation (2.44)), the two-sided bound on the hole probabilities when $U=\mathbb {B}(0,r)$ as $r\rightarrow \infty $ was proved by Sodin–Tsirelson (for $\mathbb {C}$ , [Reference Sodin and Tsirelson73, Theorem 1]) and by Zrebiec (for $\mathbb {C}^{n}$ , [Reference Zrebiec82, Theorem 1.2]). In Subsection 3.4, we will explain how to recover their findings from our general results by specializing to the scaled Bargmann–Fock spaces.

1.3 Random $\mathcal {L}^{2}$ -holomorphic sections and Toeplitz operators

In the setting of Section 1.1, we introduce for a bounded function f on X the associated Toeplitz operator $T_{f}$ defined by $T_{f}: H^{0}_{(2)}(X,L)\ni S\mapsto P(fS)\in H^{0}_{(2)}(X,L)$ , where P is the Bergman projection (1.3) (see Definition 4.4 for further details). If f is smooth and also satisfies

(1.15) $$ \begin{align} \int_{X}|f(x)|P(x,x)\mathrm{dV}(x)<\infty, \end{align} $$

then the operator $T_{f}$ is Hilbert–Schmidt (cf. Proposition 4.7). If in addition f is nonnegative and not identically zero, then $T_{f}$ is injective. For such a function f we get a Hilbert metric $\langle T_{f}\,\cdot ,T_{f}\,\cdot \rangle _{\mathcal {L}^{2}(X,L)}$ on $H^{0}_{(2)}(X,L)$ , which is a measurable norm in the sense of Gross (cf. [Reference Gross36]). Let $\mathcal {B}_{f}(X,L)$ be the completion of $H^{0}_{(2)}(X,L)$ under this measurable norm. The theory of abstract Wiener spaces implies that for f given as above, there exists a unique Gaussian probability measure $\mathcal {P}_{f}$ on $\mathcal {B}_{f}(X,L)$ such that it extends the Gaussian probability measure on any finite-dimensional subspace of $\mathrm {Im}(T_{f})$ associated with the standard $\mathcal {L}^{2}$ -metric. The injective operator $T_{f}$ extends to an isometry of Hilbert spaces

(1.16) $$ \begin{align} \widehat{T}_{f}:\big(\mathcal{B}_{f}(X,L),\|T_{f}\cdot\|\big) \rightarrow \big(H^{0}_{(2)}(X,L),\|\cdot\|_{\mathcal{L}^{2}(X,L)}\big). \end{align} $$

The pushforward $\mathbb {P}_{f}$ of $\mathcal {P}_{f}$ by $\widehat {T}_{f}$ is a Gaussian probability measure on $H^{0}_{(2)}(X,L)$ . We consider the random variable on $(H^{0}_{(2)}(X,L),\mathbb {P}_{f})$ given by $\boldsymbol {s}\mapsto [\operatorname {\mathrm {Div}}(\boldsymbol {s})]$ . Its expectation $\mathbb {E}^{\mathbb {P}_{f}}\left [[\operatorname {\mathrm {Div}}(\boldsymbol {s})]\right ]$ is a $(1,1)$ -current given for $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ by

$$\begin{align*}\big\langle\mathbb{E}^{\mathbb{P}_{f}}\left[[\operatorname{\mathrm{Div}}(\boldsymbol{s})]\right], \varphi\big\rangle= \int_{H^{0}_{(2)}(X,L)}\left(\int_{\operatorname{\mathrm{Div}}(\boldsymbol{s})}\varphi\right) \mathrm{d}\mathbb{P}_f(\boldsymbol{s}) .\end{align*}$$

The on-diagonal restriction $T^{2}_{f}(x,x)$ of the Schwartz kernel of $T^{2}_{f}=T_{f}\circ T_{f}$ is locally integrable on X (cf. Lemma 4.14). Analogously to Equation (1.4), we define the Fubini–Study current associated to f, by

(1.17) $$ \begin{align} \gamma_{f}(L,h_{L})=c_{1}(L,h_{L})+ \frac{\sqrt{-1}}{2\pi}\partial\bar{\partial}\log T^{2}_{f}(x,x)\,. \end{align} $$

It is a closed positive $(1,1)$ -current on X; see Remark 4.16 for their geometric interpretation. In Subsection 4.4, we prove the following result for the expectation of the random zeros of $\mathcal {L}^{2}$ -holomorphic section.

Theorem 1.6. Let $(X,J,\Theta )$ be a connected complex Hermitian manifold, and let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ . The expectation of the current valued random variable $[\operatorname {\mathrm {Div}}(\boldsymbol {s})]$ , where $\boldsymbol {s}$ runs in the Gaussian ensemble $(H^{0}_{(2)}(X,L),\mathbb {P}_{f})$ constructed via Equation (1.16), exists and is given by the Fubini–Study current (1.17), that is,

(1.18) $$ \begin{align} \mathbb{E}^{\mathbb{P}_{f}}\left[[\operatorname{\mathrm{Div}}(\boldsymbol{s})]\right]=\gamma_{f}(L,h_{L}). \end{align} $$

Remark 1.7. The equality (1.18) shows that the expectation of the zeros of $\widehat {T}_f\,\widetilde {\boldsymbol {s}}$ for $\widetilde {\boldsymbol {s}}$ running in the Gaussian ensemble $(\mathcal {B}_{f}(X,L),\mathcal {P}_f)$ , equals the Fubini–Study current (1.17). All the results pertaining to $\mathbb {P}_{f}$ , such as Theorems 1.8 and 1.9, can be viewed as concerning the zeros of $\widehat {T}_f\,\widetilde {\boldsymbol {s}}$ . During the writing of this paper, we became aware of the work of Ancona and Le Floch [Reference Ancona and Le Floch1] on zeros of images $T_{f}\boldsymbol {s}$ of random sections $\boldsymbol {s}$ by a Toeplitz operator $T_f$ , in the case of a compact Kähler manifold X. For compact X, one has $d<\infty $ (we assume that $d>0$ ), and in this case, the random section $\boldsymbol {s}$ in $H^{0}(X,L)$ with respect to the probability measure $\mathbb {P}_{f}$ defined above has the same distribution as the random section $T_{f}\widetilde {\boldsymbol {s}}$ considered by Ancona and Le Floch, where $\widetilde {\boldsymbol {s}}$ is a random section in the standard Gaussian ensemble $(H^{0}(X,L),\mathbb {P}_{\mathrm {st}})$ defined by the inner product $\langle \cdot ,\cdot \rangle _{\mathcal {L}^{2}}$ . In this case, we have $\mathcal {P}_f=\mathbb {P}_{\mathrm {st}}$ , so $\mathbb {P}_f$ is the pushforward of $\mathbb {P}_{\mathrm {st}}$ by $T_f$ , and Equation (1.18) reads

(1.19) $$ \begin{align} \mathbb{E}\left[[\operatorname{\mathrm{Div}}(T_f\,\widetilde{\boldsymbol{s}}\,)]\right]= \int_{H^0(X,L)}[\operatorname{\mathrm{Div}}(T_f\,\widetilde{\boldsymbol{s}}\,)] \mathrm{d}\mathbb{P}_{\mathrm{st}}(\,\widetilde{\boldsymbol{s}}\,)=\gamma_{f}(L,h_{L}); \end{align} $$

that is, the expectation of the zeros of $T_f\,\widetilde {\boldsymbol {s}}$ for $\widetilde {\boldsymbol {s}}$ in the Gaussian ensemble $(H^0(X,L),\mathbb {P}_{\mathrm {st}})$ equals the Fubini–Study current (1.17).

1.4 High tensor powers of L: equidistribution on the support of f

To consider the semiclassical limit in the noncompact setting, we need to impose the same assumptions as in Subsection 1.2. For simplicity, in this subsection we only consider a nontrivial nonnegative smooth function f on X with compact support. Note that our results hold for a general class of nonnegative smooth functions f that are not required to have compact support (cf. Subsections 5.1 and 5.2).

Since f has compact support, condition (1.15) is satisfied for the line bundle $L^{p}$ for each p. This way, we can construct a sequence of probability spaces $(H^{0}_{(2)}(X,L^{p}),\mathbb {P}_{f,p})$ using the corresponding Toeplitz operator $T_{f,p}$ . We denote by $\boldsymbol {s}_{f,p}$ the identity map on the canonical probability space $(H^{0}_{(2)}(X,L^{p}),\mathbb {P}_{f,p})$ .

In Theorems 5.7, 5.9 and 5.11, we prove the general version of the following results.

Theorem 1.8. Let $(X,J,\Theta )$ be a connected complex Hermitian manifold, and let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ such that Condition 1.2 holds. Let U be a (relatively compact) open subset of X such that $f>0$ on $\overline {U}$ . Then there exists $p_{U}\in \mathbb {N}$ such that for all $p\geq p_{U}$ , $\gamma _{f}(L^{p},h_{p})$ is smooth on U, and we have the following expansion of smooth forms on U in any $\mathscr {C}^{\ell }$ -norm for a given $\ell \in \mathbb {N}$ as $p\to \infty $ ,

(1.20) $$ \begin{align} \gamma_{f}(L^{p},h_{p})= pc_{1}(L,h_{L})+\frac{\sqrt{-1}}{2\pi}\partial\overline\partial\log\left( \frac{c_1(L,h_L)^n}{\Theta ^n}f^2\right) + \mathcal{O}\left(p^{-1}\right). \end{align} $$

We have the convergence of smooth forms in any $\mathscr {C}^{\ell }$ -norm (induced by $g^{TX}$ and the associated Levi–Civita connection) on U given by

(1.21) $$ \begin{align} \frac{1}{p}\mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}]=\frac{1}{p}\gamma_{f}(L^{p},h_{p})|_{U}\rightarrow c_{1}(L,h_{L})|_{U}\,,\:\:p\rightarrow \infty. \end{align} $$

Moreover, for any $\varphi \in \Omega ^{n-1,n-1}_{0}(U)$ , we have

(1.22) $$ \begin{align} \mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}\langle[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})], \varphi\rangle=\langle c_{1}(L,h_{L}),\varphi\rangle\right)=1. \end{align} $$

If $\Theta $ is a Kähler form and $\int _{X} c_{1}(L,h_{L})\wedge \Theta ^{n-1}<\infty $ , or if $\int _{X} c_{1}(L,h_{L})^{n}<\infty $ , then

(1.23) $$ \begin{align} \mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}=c_{1}(L,h_{L})|_{U} \right)=1, \end{align} $$

where the limit is taken with respect to the weak convergence of $(1,1)$ -currents on U.

If we allow f to vanish at some points in $\operatorname {\mathrm {supp}} f$ , we observe that the smallest vanishing order of f at a zero is two since $f\geq 0$ . In this situation, we have the following result.

Theorem 1.9. We assume the same conditions for $(X,J,\Theta )$ and $(L,h_{L})$ as in Theorem 1.8 and in addition we assume that the prequantum condition $\Theta =c_{1}(L,h_{L})$ holds. Fix a nontrivial nonnegative smooth function f on X as above. Let U be an open subset of $\operatorname {\mathrm {supp}} f$ such that f only vanishes up to second order in U with nonzero $\Delta f$ at the zeros. Then the convergence (1.21) holds in the sense of $(1,1)$ -currents on U, and the conclusions (1.22), (1.23) of Theorem 1.8 hold.

A general version of the above theorem is provided in Theorem 5.14. One important ingredient in the proofs of the above results is the following identity from Theorem 1.6,

(1.24) $$ \begin{align} \mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{p,f})]]-pc_{1}(L,h_{L})= \frac{\sqrt{-1}}{2\pi}\partial\overline{\partial}\log(T^{2}_{f,p}(x,x)). \end{align} $$

We consider next the fluctuations of zeros on a small geodesic ball $B(x,R/\sqrt {p})$ centered at x, via pairing with a test form $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ . In Subsection 5.3 and under the prequantum condition, our computations (especially by Theorem 5.22 for the case $f(x)=0$ , $\Delta f(x)<0$ ) show that

(1.25) $$ \begin{align} \begin{aligned} \Big\langle\mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]]-pc_{1}(L,h_{L}),\chi_{B\big(x,\frac{R}{\sqrt{p}}\big)}\varphi\Big\rangle =\begin{cases} \mathcal{O}(p^{-n}), \!\!\!\!&\text{if } f(x)>0,\\\mathcal{O}(p^{-n+1}),\!\!\!\!&\text{if } f(x)=0, \Delta f(x)<0, \end{cases} \end{aligned} \end{align} $$

where $\chi _{B(x,R/\sqrt {p})}$ is the indicator function of the set $B(x,R/\sqrt {p})$ , and the first case $f(x)>0$ follows from an analog of [Reference Ancona and Le Floch1, Theorem 1.7] (cf. (1.20)). The coefficients of $p^{-n}$ and of $p^{-n+1}$ in the above estimates can be worked out explicitly. The different powers in Equation (1.25) show that, in the Planck scale, zeros of random sections have higher fluctuations near a zero of f of order $2$ than near points where f does not vanish.

At last, in Subsection 5.4, we consider a not necessarily nonnegative real smooth function f satisfying condition (1.15) for $L^{p}$ , $p\gg 0$ . In this case, $T_{f,p}$ might not be injective and with suitable conditions on the vanishing points of f, we can still extend Theorem 1.9 to this case.

The next four sections of this paper correspond exactly to the above four subsections describing the main results: The first two sections deal with Gaussian random holomorphic sections, and the last two sections deal with random $\mathcal {L}^{2}$ -holomorphic sections using the Toeplitz operators.

2 Gaussian random holomorphic sections

In this section, we define Gaussian ensembles of holomorphic sections of L and study their zeros as a $(1,1)$ -currents on X. While some results proved in this section are not new in the special case of random functions or power series, we were not able to locate these results for holomorphic sections with $d=\infty $ in the literature.

2.1 Holomorphic line bundles and Bergman kernels

Let $(X,J,\Theta )$ be a connected complex Hermitian manifold (without boundary) of complex dimension $n\geq 1$ . To $\Theta $ one can associate a J-invariant smooth Riemannian metric $g^{TX}(\cdot ,\cdot )=\Theta (\cdot ,J\cdot )$ . Let L be a holomorphic line bundle over X, and let $h_{L}$ be a smooth Hermitian metric on L. We denote the corresponding Chern curvature form of L by $R^L$ , and the first Chern form of $(L, h_{L})$ is denoted by

(2.1) $$ \begin{align} c_{1}(L, h_{L})=\frac{\sqrt{-1}}{2\pi}R^{L}. \end{align} $$

Let $\mathscr {C}_0^{\infty }(X,L)$ denote the space of compactly supported smooth sections of L on X. Associated with the metrics $g^{TX}$ and $h_{L}$ , we define the $\mathcal {L}^{2}$ -inner product as follows, for $s_{1}, s_{2}\in \mathscr {C}_0^{\infty }(X,L)$ ,

(2.2) $$ \begin{align} \langle s_1,s_2\rangle_{\mathcal{L}^{2}(X,L)}:=\int_X\langle s_1(x),s_2(x)\rangle_{h_L} \mathrm{dV}(x), \end{align} $$

where $\mathrm {dV}=\frac {1}{n!}\Theta ^n$ is the volume form induced by $\Theta .$ We also let $\mathcal {L}^2(X,L)$ be the separable Hilbert space obtained by completing $\mathscr {C}_0^{\infty }(X,L)$ with respect to the norm $\|\cdot \|_{\mathcal {L}^{2}(X,L)}$ induced by Equation (2.2). Let $H^{0}(X,L)$ denote the vector space of holomorphic sections of L over X. Set

(2.3) $$ \begin{align} H_{(2)}^0(X,L):= \mathcal{L}^2(X,L) \cap H^0(X,L). \end{align} $$

It follows from the Cauchy estimates for holomorphic functions that for every compact set $K\subset X$ there exists $C_K>0$ such that

(2.4) $$ \begin{align} \sup_{x\in K}|s(x)|_{h_{L}} \leq C_K\|s\|_{\mathcal{L}^2(X,L)}\ \ \ \text{for}\ s\in H_{(2)}^0(X,L), \end{align} $$

which in turn implies that $H_{(2)}^0(X,L)$ is a closed subspace of $\mathcal {L}^2(X,L)$ . Moreover, $H_{(2)}^0(X,L)$ is a separable Hilbert space with induced $\mathcal {L}^{2}$ -metric (cf. [Reference Weil78, p. Reference Ma and Marinescu60]).

The evaluation functional $H_{(2)}^0(X,L)\ni S\mapsto S(x)$ is continuous by estimate (2.4), so by Riesz representation theorem for each $x\in X$ there exists $P(x,\cdot )\in \mathcal {L}^2(X,L_x\otimes L^*)$ such that

$$ \begin{align*}s(x)=\int_X P(x,y)s(y)\, \mathrm{dV}(y)\,,\quad \text{for all } s\in H_{(2)}^0(X,L).\end{align*} $$

Set

(2.5) $$ \begin{align} d=\dim H_{(2)}^0(X,L)\in \mathbb{N}\cup\{\infty\}. \end{align} $$

If X is compact, then $d<\infty $ . If $d\geq 1$ , consider an orthonormal basis $\{S_j\}_{j=1}^{d}$ of $H_{(2)}^0(X,L)$ . Then the series $\sum _{j=1}^{d}S_j(x)\otimes (S_j(y))^*$ converges uniformly on every compact together with all its derivatives (cf. [Reference Atiyah2, Proposition 2.4], [Reference Ma and Marinescu56, Remark 1.4.3], [Reference Weil78, p. Reference Offord63]). In particular, $P(x,y)$ is smooth on $X\times X$ . It follows that

(2.6) $$ \begin{align} P(x,y)=\sum_{j=1}^{d}S_j(x)\otimes (S_j(y))^*\,. \end{align} $$

We obtain thus for the Bergman projection (1.3),

$$ \begin{align*}(P s)(x)=\int_X P(x,y)s(y)\,\mathrm{dV}(y),\end{align*} $$

that is, $P(x,y)$ is the integral kernel of the Bergman projection. Recall that the line bundle $L\boxtimes L^*$ on $X\times X$ has fibres $(L\boxtimes L^*)_{(x,y)}:=L_{x}\otimes L^{\ast }_{y}$ for $(x,y)\in X\times X$ . The section $P(\cdot ,\cdot )$ of $L\boxtimes L^*\to X\times X$ is called Bergman kernel.

The canonical identification $L_x\otimes L_x^*=\operatorname {\mathrm {End}}(L_x)=\mathbb {C}$ , $s\otimes s^*\mapsto s^*(s)=|s|^2_{h_L}$ allows to identify $P(x,x)$ as the smooth function

(2.7) $$ \begin{align} P(x,x)=\sum_{j=1}^{d}|S_j(x)|^2_{h_L}, \end{align} $$

called the Bergman kernel function. We deduce that $d=\int _X P(x,x)\, \mathrm {dV}(x)\in \mathbb {N}\cup \{\infty \}.$ Hence, the Bergman kernel function is the dimensional density of $H_{(2)}^0(X,L)$ . If $d=0$ , then the above considerations are trivially true.

2.2 Gaussian random holomorphic sections

The results proved in this subsection are extensions of the well-known results for random power series or random analytic functions on $\mathbb {C}^{n}$ (cf. [Reference Kahane46] or [Reference Edelman and Kostlan29, Section 3]) to the complex geometric setting. We include details of the proofs for the sake of completeness.

Let $\eta =\{\eta _{j}\}_{j\in \mathbb {N}}$ be a sequence of i.i.d. centered real or complex Gaussian random variables, and denote by $\mathbb {P}$ and $\mathbb {E}$ the underlying probability measure and its expectation.

For $d\geq 1$ , let $S=\{S_{j}\}_{j=1}^{d}$ be an orthonormal basis of $H^{0}_{(2)}(X,L)$ . Define

(2.8) $$ \begin{align} \psi^{S}_{\eta}(x)=\sum_{j=1}^{d}\eta_{j}S_{j}(x). \end{align} $$

If $d=0$ , we simply set $\psi ^{S}_{\eta }\equiv 0$ .

Proposition 2.1. The section $\psi ^{S}_{\eta }$ is almost surely a holomorphic section of L on X.

Proof. If d is finite, the claim is clearly true. Hence, it remains to prove it for the case $d=\infty .$ In this case, X is noncompact. Let $\{K_{i}\}_{i\in \mathbb {N}}$ be an increasing sequence of compact subsets of X such that $X=\cup _{i\in \mathbb {N}} K_{i}$ . We can take each $K_{i}$ to be the closure of a relatively compact open subset $U_{i}$ of X. Then to prove this proposition, we only need to show that for each i, $\psi ^{S}_{\eta }$ is almost surely a holomorphic section of L on $U_{i}$ .

Let K be a compact subset of X, and let U be an open relatively compact neighborhood of K. Similarly to estimate (2.4), there exists a constant $C_{U}>0$ such that for $s\in H^{0}_{(2)}(X,L)$ ,

(2.9) $$ \begin{align} \sup_{x\in K} |s(x)|_{h_{L}}\leq C_{U}\|s\|_{\mathcal{L}^{2}(\overline{U},L)}. \end{align} $$

By Equation (2.7), we have

(2.10) $$ \begin{align} \sum_{j=1}^{d}\|S_{j}\|^{2}_{\mathcal{L}^{2}(\overline{U},L)}=\int_{x\in\overline{U}}P(x,x)\mathrm{dV}(x)<+\infty. \end{align} $$

For $j\in \mathbb {N}_{>0}$ , $x\in X$ we consider the $L_{x}$ -valued random variable $X_{j}(x)=\eta _{j}S_{j}(x)$ . Since $\eta _{j}$ is centered, we infer

(2.11) $$ \begin{align} \mathbb{E}[X_{j}(x)]=0\in L_{x}. \end{align} $$

It is then consistent to define the variance as

$$ \begin{align*} \mathrm{Var}(X_{j}(x))=\mathbb{E}[|X_{j}(x)|_{h_{L}}^{2}], \end{align*} $$

and we can compute

(2.12) $$ \begin{align} \mathrm{Var}(X_{j}(x))=\mathbb{E}[|X_{j}(x)|_{h_{L}}^{2}]= \mathrm{Var}(\eta_{j})|S_{j}(x)|^{2}_{h_{L}}= \mathrm{Var}(\eta_{1})|S_{j}(x)|^{2}_{h_{L}}. \end{align} $$

We next prove that for any $k\in \mathbb {N}, N\in \mathbb {N}_{>0}$ and for $r>0$ , we have

(2.13) $$ \begin{align} \mathbb{P}\Big(\sup_{\ell=1,\ldots, N}\sup_{x\in K}\big|\sum_{j=1}^{\ell}X_{k+j}(x)\big|_{h_{L}}>r\Big)<\frac{C^{2}_{U}\mathrm{Var}(\eta_{1})}{r^{2}}\sum_{j=1}^{N}\|S_{k+j}\|^{2}_{\mathcal{L}^{2}(\overline{U},L)}. \end{align} $$

For this purpose, define the stochastic process

(2.14) $$ \begin{align} Y_{\ell}=\big\|\sum_{j=1}^{\ell}X_{k+j}\big\|^2_{\mathcal{L}^{2}(\overline{U},L)}, \quad \ell =1, \ldots, N, \end{align} $$

and observe that by virtue of Equation (2.9), we have

(2.15) $$ \begin{align} \sup_{x\in K}\big|\sum_{j=1}^{\ell}X_{k+j}(x)\big|_{h_{L}}\leq C_{U}Y_{\ell}^{\frac12}. \end{align} $$

As a consequence, we have

(2.16) $$ \begin{align} \mathbb{P}\Big(\sup_{\ell=1,\ldots, N}\sup_{x\in K}\big|\sum_{j=1}^{\ell}X_{k+j}(x)\big|_{h_{L}}>r\Big) \le \mathbb{P}\Big(\sup_{\ell=1,\ldots, N} C_{U}^2 Y_{\ell} >r^2\Big). \end{align} $$

Now, the process $(Y_{\ell })$ , $\ell =1, \ldots , N,$ is a submartingale with respect to the filtration $(\mathcal F_\ell ),$ where $\mathcal F_\ell = \sigma \big (\langle X_{k+i},X_{k+j} \rangle _{\mathcal {L}^{2}(\overline {U},L)},\ i,j =1, \ldots , \ell \big )$ . Therefore, Doob’s submartingale inequality (see, e.g., [Reference Klenke47, Lemma 11.1]) yields

(2.17) $$ \begin{align} \mathbb{P}\Big(\sup_{\ell=1, \ldots, N} Y_{\ell}> \frac{r^2}{C_{U}^2}\Big) \le C_{U}^2 \frac{\mathbb{E}[Y_N]}{r^2}, \end{align} $$

which immediately entails Equation (2.13). Now, letting $N\rightarrow +\infty $ in Equation (2.13), we get

(2.18) $$ \begin{align} \mathbb{P}\Big(\sup_{\ell\in\mathbb{N}_{>0}}\sup_{x\in K}|\sum_{j=1}^{\ell}X_{k+j}(x)|_{h_{L}}>r \Big)\leq \frac{C^{2}_{U}\mathrm{Var}(\eta_{1})}{r^{2}}\sum_{j=1}^{+\infty}\|S_{k+j}\|^{2}_{\mathcal{L}^{2}(\overline{U},L)}. \end{align} $$

Then taking the limit of Equation (2.18) as $k\rightarrow \infty $ , and using Equation (2.10), we infer

(2.19) $$ \begin{align} \mathbb{P}\Big(\limsup_{k\rightarrow+\infty}\sup_{\ell\in\mathbb{N}_{>0}}\sup_{x\in K}|\sum_{j=1}^{\ell}X_{k+j}(x)|_{h_{L}}>r\Big)=0. \end{align} $$

Therefore, a union bound along the sequence of $r = \frac 1N$ immediately supplies us with

(2.20) $$ \begin{align} \mathbb{P}\Big(\limsup_{k\rightarrow+\infty}\sup_{\ell\in\mathbb{N}_{>0}}\sup_{x\in K}|\sum_{j=1}^{\ell}X_{k+j}(x)|_{h_{L}}>0\Big)=0. \end{align} $$

If we take V to be a relatively compact open subset of X, and take $K=\overline {V}$ , then by Equation (2.20), the sum $\sum _{j=1}^{\infty } X_{j}$ is almost surely uniformly convergent on K so that it almost surely defines a holomorphic section on V. This completes the proof of our proposition.

For the purpose of the following definition, we note at this point that a standard complex Gaussian is a random variable having the distribution $\frac {1}{\sqrt {2}} (X + \sqrt {-1}Y)$ , where X and Y are standard real Gaussian variables.

Definition 2.2. The random section $\psi ^{S}_{\eta }$ defined by the equality (2.8) is called a standard Gaussian random holomorphic section of L over X if $\eta =\{\eta _{j}\}_{j\in \mathbb {N}}$ is a sequence of i.i.d. standard complex Gaussian random variables.

Now, we prove that the distribution of a standard Gaussian random holomorphic section $\psi ^{S}_{\eta }$ does not depend on the choice of the orthonormal basis.

Proposition 2.3. Assume that $d\geq 1$ , and assume that $\eta =\{\eta _{j}\}_{j=1}^{d}$ is a sequence of i.i.d. standard complex Gaussian random variables. If $S'=\{S^{\prime }_{j}\}_{j=1}^{d}$ is another choice of orthonormal basis of $H^{0}_{(2)}(X,L)$ , then $\psi ^{S'}_{\eta }$ and $\psi ^{S}_{\eta }$ have the same distribution as random holomorphic sections.

Proof. It is sufficient to find a sequence $\eta '=\{\eta ^{\prime }_{j}\}_{j=1}^{d}$ of i.i.d. standard complex Gaussian random variables such that a.s. $\psi ^{S'}_{\eta }=\psi ^{S}_{\eta '}$ . Let $\ell ^{2}(\mathbb {C})$ denote the Hilbert space of $\ell ^{2}$ -summable complex sequences. If $u=(u_{j})_{j\in \mathbb {N}}\in \ell ^{2}(\mathbb {C})$ , set

(2.21) $$ \begin{align} (\eta,u)_{\ell^{2}}=\sum_{j\in\mathbb{N}} \eta_{j}\bar{u}_{j}. \end{align} $$

By Kolmogorov’s three-series theorem (cf. [Reference Varadhan77]), the series in definition (2.21) is almost surely convergent so that $(\eta ,u)_{\ell ^{2}}$ is a well-defined random variable. By the property of Gaussian random variables, we conclude that $(\eta ,u)_{\ell ^{2}}$ is a centered complex Gaussian random variable with variance $|u|^{2}_{\ell ^{2}}$ . In particular, if $|u|_{\ell ^{2}}=1$ , then $(\eta ,u)_{\ell ^{2}}$ has the same distribution as $\eta _1$ . Moreover, if nonzero $u,v\in \ell ^{2}$ is such that $(u,v)_{\ell ^{2}}=0$ , then $(\eta ,u)_{\ell ^{2}}$ and $(\eta ,v)_{\ell ^{2}}$ are independent. Take $(a_{ij}\in \mathbb {C})_{i,j\in \mathbb {N}}$ such that for each i,

(2.22) $$ \begin{align} S^{\prime}_{i}=\sum_{j\in\mathbb{N}}a_{ij}S_{j}. \end{align} $$

For $j\in \mathbb {N}$ , set $b_{j}=(\bar {a}_{ij})_{i\in \mathbb {N}}$ . Then $b_{j}\in \ell ^{2}(\mathbb {C})$ is with norm $1$ , moreover, if $j\neq j'$ , then $(b_{j},b_{j'})_{\ell ^{2}}=0$ . Now, define

(2.23) $$ \begin{align} \eta^{\prime}_{j}=(\eta,b_j)_{\ell^{2}}. \end{align} $$

Then $\eta '=(\eta ^{\prime }_{j})_{j\in \mathbb {N}}$ is a sequence of i.i.d. centered Gaussian random variables with the same distribution as $\eta $ . By definition, we get that almost surely, $\psi ^{S'}_{\eta }=\psi ^{S}_{\eta '}$ . Therefore, $\psi ^{S'}_{\eta }$ and $\psi ^{S}_{\eta }$ have the same distribution.

Remark 2.4. (a) When $d=\infty $ , note that by taking a sequence of compact subset $\{K_{i}\}_{i\in \mathbb {N}}$ as in the proof of Proposition 2.1, we can define a sequence of seminorms for $H^{0}(X,L)$ , hence a Fréchet distance so that $H^{0}(X,L)$ is a Fréchet space. In Proposition 2.1, we actually prove that $\psi ^{S}_{\eta }$ is a random variable taking values in the Fréchet space $H^{0}(X,L)$ .

(b) When $d\geq 1$ , we have

(2.24) $$ \begin{align} \mathbb{P}(\psi^{S}_{\eta}\equiv 0)=0. \end{align} $$

Indeed, since $\{S_{j}\}_{j=1}^{d}$ is a basis of $H^{0}_{(2)}(X,L)$ , we have

(2.25) $$ \begin{align} \mathbb{P}(\psi^{S}_{\eta}\equiv 0)=\mathbb{P}(\eta_{j}=0, j=1, \cdots )=0. \end{align} $$

(c) In the proof of Proposition 2.1, we do not use the Gaussianity of the $\eta _{j}$ in an essential way. Hence, we can work with any sequence $\eta $ of pairwise uncorrelated centered random variables with uniformly bounded variance. In that case, however, the distribution of the random section $\psi ^{S}_{\eta }$ might depend on the choice of the basis S. Generally, one needs suitable moment conditions on $\eta $ to obtain more results, such as the universality of the zeros of $\psi ^{S}_{\eta }$ , and we refer to [Reference Bayraktar, Coman and Marinescu6, Reference Drewitz, Liu and Marinescu27, Reference Kabluchko and Zaporozhets44] for the related details.

Lemma 2.5. If $d=\infty $ , then with probability one, $\psi ^{S}_{\eta }$ is not $\mathcal {L}^{2}$ -integrable on X.

Proof. The event that $\psi ^{S}_{\eta }$ is $\mathcal {L}^{2}$ -integrable is equivalent to the event $\{\sum _{j=1}^{\infty } |\eta _{j}|^{2}< \infty \}$ . But, for example, by the law of large numbers, we infer

(2.26) $$ \begin{align} \mathbb{P}\Big(\sum_{j=1}^{\infty} |\eta_{j}|^{2}< \infty\Big)=0, \end{align} $$

and the statement of the lemma follows.

2.3 Currents and the Lelong–Poincaré formula

The zero-set of a holomorphic section is a complex analytic set which is in general singular. The analytic tool used to deal with singularities in complex geometry is the theory of currents, introduced by de Rham [Reference de Rham18] (see [Reference Demailly22, Reference Griffiths and Harris35] and especially [Reference Federer31] for complete expositions).

Let X be a complex manifold of dimension n, and let E be a complex vector bundle on X. The space of smooth sections of E is denoted by $\mathscr {C}^\infty (X,E)$ and is endowed with the $\mathscr {C}^\infty $ -topology of uniform convergence of all derivatives on compact sets. The space of smooth sections of E with compact support is denoted by $\mathscr {C}^\infty _0(X,E)$ and is endowed with the topology of inductive limit of spaces of smooth sections with support on a given compact set. In particular, we denote by $\Omega ^{n-1,n-1}_{0}(X)$ the space of of smooth $(n-1,n-1)$ -forms with compact support.

The space of $(1,1)$ -currents on X is the topological dual of the space $\Omega ^{n-1,n-1}_{0}(X)$ (called test forms in this context). In the sequel, we let $\langle T, \varphi \rangle $ be the pairing between a $(1,1)$ -current T and a test form $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ . A $(1,1)$ -current is called of order $k\in \mathbb {N}_0$ if it is continuous in the $\mathscr {C}^k$ -topology; equivalently, it extends as a linear continuous functional to the space of $(n-1,n-1)$ -forms of class $\mathscr {C}^k$ with compact support.

For any analytic hypersurface $V\subset X$ , we define the current of integration $[V]$ on V by

$$ \begin{align*}\varphi\mapsto \int_{V}\varphi:=\int_{V_{\mathrm{reg}}}\varphi, \quad \varphi\in\Omega^{n-1,n-1}_{0}(X),\end{align*} $$

where $V_{\mathrm {reg}}$ is the regular set of V (a complex submanifold of codimension $1$ ). By a theorem of Lelong ([Reference Griffiths and Harris35, p. Reference Fefferman32] [22, III-2.7]) the current of integration on V is a closed positive $(1,1)$ -current. It is clear that $[V]$ is a current of order $0$ on X.

Let L be a holomorphic line bundle on X. For a holomorphic section $s\in H^{0}(X,L)\setminus \{0\}$ , the divisor $\operatorname {\mathrm {Div}}(s)$ of s is defined as the formal sum

(2.27) $$ \begin{align} \operatorname{\mathrm{Div}}(s)=\sum_{V \subset Z(s)} \mathrm{ord}_{V}(s) V, \end{align} $$

where V runs over all the irreducible analytic hypersurfaces contained in $Z(s)$ , and $\mathrm {ord}_{V}(s)$ denotes the vanishing order of s along V. Let $Z(s)$ denote the set of zeros of s, which is a purely $1$ -codimensional analytic subset of X. The current of integration (with multiplicities) on the divisor $\operatorname {\mathrm {Div}}(s)$ is defined by

(2.28) $$ \begin{align} [\operatorname{\mathrm{Div}}(s)]=\sum_{V \subset Z(s)} \mathrm{ord}_{V}(s) [V],\quad \big\langle[\operatorname{\mathrm{Div}}(s)],\varphi\big\rangle= \int_{\operatorname{\mathrm{Div}}(s)}\varphi:= \sum_{V \subset Z(s)} \mathrm{ord}_{V}(s)\int_V\varphi\,. \end{align} $$

Assume that L is endowed with a smooth Hermitian metric $h_L$ . By the Lelong–Poincaré formula [Reference Ma and Marinescu56, Theorem 2.3.3], we have

(2.29) $$ \begin{align} [\operatorname{\mathrm{Div}}(s)]=\frac{\sqrt{-1}}{2\pi} \partial\bar{\partial}\log|s|^{2}_{h_{L}}+ c_{1}(L,h_{L}), \quad\text{for } s\in H^{0}(X,L). \end{align} $$

This important formula is crucial for our purposes. It links the zero-divisor to the curvature and to the logarithm of the pointwise norm of a section, which is analytically easier to tackle and allows the introduction of the Bergman kernel into the picture.

2.4 Expectation of random zeros: proof of Theorem 1.1

In the sequel, we always assume $d=\dim H^{0}_{(2)}(X,L)\geq 1$ . We start with some considerations about the Fubini–Study currents.

Lemma 2.6. Assume that $d\geq 1$ . Then the function $X\ni x\mapsto \log P(x,x) \in \{-\infty \}\cup \mathbb {R}$ is locally $\mathcal {L}^{1}$ -integrable on X with respect to $\mathrm {dV}$ . Thus, $\frac {\sqrt {-1}}{2\pi }\partial \overline {\partial }\log P(x,x)$ defines a $(1,1)$ -current on X.

Proof. Let $e_{L}:U\to L$ be a local holomorphic frame of L. Let $\varphi \in \mathscr {C}^\infty (U)$ be the local weight of $h_L$ with respect to $e_L$ , that is, $|e_L|^2_{h_L}=e^{-2\varphi }$ on U. We consider an orthonormal basis $\{S_j\}_{j=1}^d$ of $H^0_{(2)}(X,L)$ and write

(2.30) $$ \begin{align} S_{j}(x)=f_{j}(x)e_{L}(x), \:\:x\in U,\:\:f_j\in\mathscr{O}(U), \end{align} $$

where $f_{j}$ are nontrivial holomorphic functions on U. Then $P(x,x)=\sum _{j=1}^d|S_j|^2_{h_L}=\sum _{j=1}^d|f_j|^2e^{-2\varphi }$ on U, hence

(2.31) $$ \begin{align} \log P(x,x)=\log\Big(\sum_{j=1}^d|f_j|^2\Big)-2\varphi\,. \end{align} $$

The series $\sum _{j=1}^d|f_j|^2$ converges locally uniformly on U, thus $\log \!\big (\sum _{j=1}^d|f_j|^2\big )$ is a plurisubharmonic function that is not identically $-\infty $ , hence locally integrable.

Lemma 2.6 shows that the Fubini–Study currents (1.4) are well defined. Note that $c_1(L,h_L)|_U=\frac {\sqrt {-1}}{\pi }\partial \overline {\partial }\varphi $ . By applying $\partial \overline {\partial }$ to both sides of equality (2.31) and taking into account formula (1.4), we see that

(2.32) $$ \begin{align} \gamma(L,h_{L})\big|_U=\frac{\sqrt{-1}}{2\pi}\partial\overline{\partial} \log\Big(\sum_{j=1}^d|f_{j}|^{2}\Big)\,, \end{align} $$

thus $\gamma (L,h_{L})$ is a closed positive $(1,1)$ -current. The base locus of $H^{0}_{(2)}(X,L)$ is the proper analytic set

(2.33) $$ \begin{align} \operatorname{\mathrm{Bl}}(X,L):=\big\{x\in X\;|\; s(x)=0 \text{ for all } s\in H^{0}_{(2)}(X,L)\big\}. \end{align} $$

Thus, $\{x\in X:P(x,x)=0\}=\operatorname {\mathrm {Bl}}(X,L)$ . Hence, $\gamma (L,h_{L})$ is a smooth form if $\operatorname {\mathrm {Bl}}(X,L)=\varnothing $ .

Remark 2.7. The Fubini–Study currents have the following geometric interpretation that justifies their name. Note that if X is compact and $\operatorname {\mathrm {Bl}}(X,L)=\varnothing $ , then $\gamma (L,h_{L})$ is the pullback of the Fubini–Study form on the projective space by the Kodaira map defined by $H^0(X,L)$ (see, e.g., [Reference Ma and Marinescu56, (5.1.21)]). If X is noncompact and $H^0_{(2)}(X,L)$ is infinite dimensional, one can still define a Kodaira map to the infinite-dimensional projective space $\mathbb {C}\mathbb {P}^{\infty }$ , cf. [Reference Kobayashi48]. This is a Hilbert manifold $\mathbb {C}\mathbb {P}^{\infty }\!=\ell ^2(\mathbb {C})\setminus \{0\}/\mathbb {C}^*$ modeled on the space $\ell ^2(\mathbb {C})$ . The Fubini–Study metric on $\mathbb {C}\mathbb {P}^{\infty }$ is defined by $\omega _{\mathrm {FS},\infty }= \frac {\sqrt {-1}}{2\pi }\partial \overline \partial \log \|a\|^2$ . Assume that $H^0_{(2)}(X,L)$ has no base locus, and consider an orthonormal basis $\{S_j\}_{j=1}^\infty $ of $H^0_{(2)}(X,L)$ . With the representation (2.30), we define the Kodaira map $\Phi :X\to \mathbb {C}\mathbb {P}^{\infty }$ by $\Phi (x)=[(f_j(x))_j]$ for $x\in U$ . Then $\Phi ^*\omega _{\mathrm {FS},\infty }= \frac {\sqrt {-1}}{2\pi }\partial \overline \partial \log (\sum _j|f_j|^2)$ coincides with formula (2.32).

Now, we are ready to prove Theorem 1.1. Let $\eta =\{\eta _{j}\}_{j=1}^{d}$ is a sequence of i.i.d. standard complex Gaussian random variables. Let $\psi ^{S}_{\eta }$ be the random holomorphic section defined in formula (2.8), and let $[\operatorname {\mathrm {Div}}(\psi ^{S}_{\eta })]$ denote the $(1,1)$ -current given by its zeros (cf. formula (2.27)).

Proof of Theorem 1.1.

We fix a test form $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ and evaluate $\mathbb {E}[\langle [\operatorname {\mathrm {Div}}(\psi ^{S}_{\eta })], \varphi \rangle ]$ . By applying the Lelong–Poincaré formula (2.29) to $\psi ^{S}_{\eta }$ , we get

(2.34) $$ \begin{align} \begin{aligned} \langle [\operatorname{\mathrm{Div}}(\psi^{S}_{\eta})],\varphi\rangle &= \int_{X}\big(\frac{\sqrt{-1}}{2\pi}\partial\bar{\partial} \log|\psi^{S}_{\eta}|^{2}_{h_{L}}+ c_{1}(L,h_{L})\big)\wedge \varphi\\ &=\int_{X}c_{1}(L,h_{L})\wedge\varphi+ \frac{\sqrt{-1}}{2\pi}\int_{X} \log|\psi^{S}_{\eta}|_{h_{L}}^{2} \partial\bar{\partial}\varphi\\ &=\int_{X}c_{1}(L,h_{L})\wedge\varphi+\frac{\sqrt{-1}}{2\pi} \int_{X\backslash \text{Bl}(X,L)}\log|\psi^{S}_{\eta}|_{h_{L}}^{2} \partial\bar{\partial}\varphi. \end{aligned} \end{align} $$

The last identity follows from that $\log |\psi ^{S}_{\eta }|_{h_{L}}^{2}$ is almost surely locally integrable by Proposition 2.1 and Remark 2.4 (b), and that $\text {Bl}(X,L)$ has Lebesgue measure zero in X. For $x\in X\backslash \text {Bl}(X,L)$ , we have $P(x,x)\neq 0$ , and let $e_{L}(x)$ be a unit vector of L at x, define

(2.35) $$ \begin{align} b(x)=\big(P(x,x)^{-1/2}S_{j}(x)/e_{L}(x)\big)_{j\in\mathbb{N}}\in \ell^{2}(\mathbb{C}). \end{align} $$

We have $|b(x)|_{\ell ^{2}}=1$ . Note that

(2.36) $$ \begin{align} P(x,x)^{-1/2}\psi^{S}_{\eta}=(\eta,\overline{b(x)})_{\ell^{2}}\,e_{L}(x). \end{align} $$

Then

(2.37) $$ \begin{align} \mathbb{E}[\log|P(x,x)^{-1/2}\psi^{S}_{\eta}|_{h_{L}}^{2}]=\mathbb{E}[\log\big|(\eta,\overline{b(x)})_{\ell^{2}}\big|^{2}]=\mathbb{E}[\log|\eta_{1}|^{2}]. \end{align} $$

Note that $\mathbb {E}[|\log |\eta _{1}|^{2}|]<\infty $ . By Lemma 2.6, $\log P(x,x)$ is locally integrable on X, then we can apply the Fubini’s theorem to the following integrals so that

(2.38) $$ \begin{align} \begin{aligned} \mathbb{E}\big[\int_{X\backslash \text{Bl}(X,L)} \log|P(x,x)^{-1/2}\psi^{S}_{\eta}|_{h_{L}}^{2} \partial\bar{\partial} \varphi\big] &=\int_{X\backslash \text{Bl}(X,L)} \mathbb{E}\big[\log|P(x,x)^{-1/2}\psi^{S}_{\eta}|_{h_{L}}^{2}\big] \partial\bar{\partial} \varphi\\ &=\mathbb{E}[\log|\eta_{1}|^{2}]\int_{X}\partial\bar{\partial} \varphi=0. \end{aligned} \end{align} $$

Then by formula (2.34), we get

(2.39) $$ \begin{align} \begin{aligned} \mathbb{E}\left[\langle[\operatorname{\mathrm{Div}}(\psi^{S}_{\eta})],\varphi\rangle \right] &=\int_{X}c_{1}(L,h_{L})\wedge\varphi+ \frac{\sqrt{-1}}{2\pi} \int_{X\backslash\text{Bl}(X,L)}\log{P(x,x)}\cdot \partial\bar{\partial}\varphi\\ &=\big\langle c_{1}(L,h_{L})+\frac{\sqrt{-1}}{2\pi}\partial\bar{\partial}\log(P(x,x)), \varphi\big\rangle. \end{aligned} \end{align} $$

This completes the proof.

2.5 Geometric examples

We present in this subsection some interesting examples of Bergman spaces and Fubini–Study currents where our results apply. We start with some simple observations.

(i) If $P(x,x)>0$ (equivalently, x is not in the base locus of $H^0_{(2)}(X,L)$ ), then the $(1,1)$ -form $\sqrt {-1}\partial {\overline \partial }\log P(x,x)$ is smooth in a neighborhood of x, and hence $\gamma (L,h_L)$ , too. In particular, if $\operatorname {\mathrm {Bl}}(X,L)=\varnothing $ , then $\gamma (L,h_L)$ is smooth.

(ii) If $P(x,x)>0$ , let $s_0\in H^0_{(2)}(X,L)$ with $s_0(x)\neq 0$ . Assume that there exist $s_1,\ldots ,s_n\in H^0_{(0)}(X,L)$ such that $d(s_1/s_0)(x),\ldots , d(s_n/s_0)(x)$ are linearly independent (that is, sections of $H^0_{(2)}(X,L)$ give local coordinates at x). Then $\gamma (L,h_L)$ is strictly positive near x.

(iii) Thus, if $\operatorname {\mathrm {Bl}}(X,L)=\varnothing $ and sections of $H^0_{(2)}(X,L)$ give local coordinates at any point in X, then $\gamma (L,h_L)$ defines a Kähler metric on X.

Example 2.8 (Bergman metric).

We consider the case when L is the canonical bundle $K_X$ of X (cf. [Reference Kobayashi48, Reference Weil78]). The space of holomorphic sections of $K_X$ is the space $H^{n,0}(X)$ of holomorphic $(n,0)$ -forms. Such a form can be written in local coordinates $(z_1,\ldots ,z_n)$ as $f(z)dz_1\wedge \ldots \wedge dz_n$ , with f a holomorphic function. We say that a measurable $(n,0)$ -form $\beta $ is an $\mathcal {L}^{2}$ section of $K_X$ if

(2.40) $$ \begin{align} \|\beta\|^2:=2^{-n}(\sqrt{-1})^{n^2} \int_X\beta\wedge\overline\beta<\infty. \end{align} $$

We denote by $H^{n,0}_{(2)}(X)$ the space of $\mathcal {L}^{2}$ holomorphic $(n,0)$ -forms. We have $H^{n,0}_{(2)}(X)=H^{0}_{(2)}(X,K_X)$ , where the right-hand side is defined with respect to an arbitrary metric $\Theta $ on X and the metric on $K_X$ is induced by $\Theta $ .

We assume that $H^{n,0}_{(2)}(X)\neq \{0\}$ and let $\{\beta _j\}_{j=1}^d$ be an orthonormal basis of $H^{n,0}_{(2)}(X)$ . In local coordinates $(U;z_1,\ldots ,z_n)$ write $\beta _j=f_j(z)dz_1\wedge \ldots \wedge dz_n$ . According to (2.32) the Fubini–Study current is given on U by $\gamma (K_X,h_{K_X})|_U=\frac {\sqrt {-1}}{2\pi }\partial {\overline \partial }\log (\sum _j |f_j|^2)$ . If the Fubini–Study current is actually a Kähler metric on X, then it is called the Bergman metric of X. We will denote it by $\omega _B$ . The metric $\omega _B$ is invariant by the group of biholomorphic transformations of X.

If X is an open set in $\mathbb {C}^n$ , the canonical bundle is trivial, so we identify the space $H^{n,0}_{(2)}(X)$ of $\mathcal {L}^{2}$ -holomorphic $(n,0)$ -forms with the space $H^{0}_{(2)}(X)$ holomorphic functions which are $\mathcal {L}^{2}$ with respect to the Lebesgue measure. There is a vast literature on Bergman spaces and kernels on domains in $\mathbb {C}^n$ ; see, for example, [Reference Hedenmalm, Korenblum and Zhu38, Reference Jarnicki and Pflug43] and the references therein.

To give concrete examples let us recall the definition of Stein manifolds, which are interesting due to their rich function-theoretical structure [Reference Grauert and Remmert34]. For a complex manifold X, let $\mathcal {O}(X)$ denote the space of all holomorphic functions on X.

Definition 2.9. A complex manifold X is called Stein if the following two conditions are satisfied: (1) X is homomorphically convex, that is, for every compact subset $K\subset X$ , its holomorphically convex hull $\widehat {K}=\big \{z\in X: |f(z)|\leq \sup _{w\in K}|f(w)|, \forall \, f\in \mathcal {O}(X)\big \}$ is compact. (2) X is holomorphically separable, that is, if $x\neq y$ in X, then there exists $f\in \mathcal {O}(X)$ such that $f(x)\neq f(y)$ .

Let $L\to X$ be a holomorphic line bundle. The cohomology vanishing theorem for coherent analytic sheaves on Stein manifolds (Cartan’s theorem B, cf. [Reference Grauert and Remmert34]) yields the following:

(i) The holomorphic sections $H^0(X,L)$ give local coordinates at each point of X.

(ii) For any closed discrete set $A=\{p_k:k\in \mathbb {N}\}\subset X$ and any family $\{v_k\in L_{p_k}:k\in \mathbb {N}\}$ , there exists $s\in H^0(X,L)$ with $s(p_k)=v_k$ for all $k\in \mathbb {N}$ . In particular, for each $p\in X$ the evaluation map $H^0(X,L)\to L_p$ is surjective and we have $\dim H^0(X,L)=\infty $ .

Example 2.10. Let X be a Stein manifold and $D\Subset X$ be a relatively compact domain. We consider a Hermitian metric on X whose associated $(1,1)$ -form is denoted by $\Theta $ . Let $\mathrm {dV}_\Theta =\Theta ^n/n!$ be the volume form induced by $\Theta $ , where $\dim X=n$ . Let $(L,h_L)$ be a Hermitian holomorphic line bundle. Consider the space $\mathcal {L}^{2}(D,L,h_L,\mathrm {dV}_\Theta )$ of measurable sections S of L over D satisfying $\int _D|S|^2_{h_L}\mathrm {dV}_\Theta <\infty $ , and let $H^0_{(2)}(D,L,h_L,\mathrm {dV}_\Theta )=\mathcal {L}^{2}(D,L,h_L,\mathrm {dV}_\Theta )\cap H^0(X,L)$ . The restriction map $H^0(X,L)\to H^0_{(2)}(D,L,h_L,\mathrm {dV}_\Theta )$ is well defined and injective. We deduce that the space $H^0_{(2)}(D,L,h_L,\mathrm {dV}_\Theta )$ is infinite dimensional, has empty base locus and sections of this space give local coordinates at any point of D. Therefore, $\gamma (L,h_L)$ is smooth on X and if $(L,h_L)$ is semipositive (i.e., $c_1(L,h_L)$ is positive semidefinite), it is a Kähler form.

We deduce from Theorem 1.1 and the discussion from Example 2.8 the following.

Corollary 2.11. For any relatively compact domain $D\Subset X$ in a Stein manifold, the expectation of the zero divisors of the standard Gaussian random holomorphic $(n,0)$ -forms defined from the $\mathcal {L}^{2}$ -holomorphic $(n,0)$ -forms on D is given by the Bergman metric on D. If $D\Subset \mathbb {C}^n$ this is true for standard Gaussian random holomorphic functions defined from the $\mathcal {L}^{2}$ -holomorphic functions on D.

One of the simplest examples is the unit disc $\mathbb {D}\subset \mathbb {C}$ endowed with the Lebesgue measure. Then $P(z,z)=\frac {1}{\pi (1-|z|^2)^2}$ and the Bergman metric

(2.41) $$ \begin{align} \omega_B=\frac{\sqrt{-1}}{\pi}\frac{dz\wedge d\bar{z}}{(1-|z|^{2})^{2}} \end{align} $$

is the hyperbolic metric (up to a constant factor) on the disc. We see on this example that the Bergman metric blows up as $|z|\to 1$ , so the zeros accumulate towards the boundary of $\mathbb {D}$ . This is a more general phenomenon. Assume that $D\Subset \mathbb {C}^n$ is a strictly pseudoconvex domain with smooth boundary. Then the Bergman metric blows up at the boundary in the nontangential directions. More precisely, let $\varrho $ be a defining function of D so that $D=\{\varrho <0\}$ , $d\varrho \neq 0$ on $\partial {D}$ and the Levi form of $\varrho $ is positive definite on the complex tangent space of $\partial {D}$ . For $v\in \mathbb {C}^n$ and $x\in D$ , we denote by $|v|_{D,x}$ the norm of the vector $\sum _{j=1}^nv_j\frac {\partial }{\partial z_j}(x)\in T^{1,0}_x\mathbb {C}^n$ with respect to the Bergman metric. Then for any $x_0\in \partial {D}$ , there exists a neighborhood U of $x_0$ and $C>0$ such that for any $x\in D\cap U$ and $v\in \mathbb {C}^n$ we have $|v|_{D,x}\geq C|\partial \varrho (x)\cdot v|/|\varrho (x)|$ ; see [Reference Diederich23, Reference Fefferman32]. There is a large amount of works about the boundary behavior of the Bergman metric. In addition to the accurate behavior of strictly pseudoconvex domains, there is another feature which indicates that zeros of random $L^2$ holomorphic sections tend to accumulate up to the boundary. This is the completeness of the Bergman metric. It is known for example that the Bergman metric is complete if D is a bounded pseudoconvex domain in $\mathbb {C}^n$ with a $\mathscr {C}^1$ boundary [Reference Ohsawa66] or if D is a bounded hyperconvex domain in $\mathbb {C}^{n}$ (cf. [Reference Błocki and Pflug9, Reference Herbort39]).

Example 2.12 (Singular Hermitian metrics).

As in the case of a smooth Hermitian metric, a singular Hermitian metric on L is defined in terms of an open cover $X=\bigcup U_\alpha $ for which there exist local holomorphic frames $e_\alpha :U_\alpha \to L$ . Consider the transition functions $g_{\alpha \beta }=e_\beta /e_\alpha \in {\mathcal O}^\star _X (U_\alpha \cap U_\beta )$ . If $h_L$ is a singular Hermitian metric on L, then (see [Reference Demailly21], also [Reference Ma and Marinescu56, p. 97]) $|e_\alpha |^2_{h_L}=e^{-2\varphi _\alpha }$ , where the functions $\varphi _\alpha \in L^1_{\mathrm {loc}}(U_\alpha )$ are called the local weights of the metric h. One has $\varphi _\alpha =\varphi _\beta + \log |g_{\alpha \beta }|$ on $U_\alpha \cap U_\beta $ , and the curvature of $h_L$ ,

$$ \begin{align*}c_1(L,h_L)|_{U_\alpha}=\frac{\sqrt{-1}}{\pi}\partial{\overline\partial}\varphi_\alpha,\end{align*} $$

is a well-defined closed (1,1) current on X. The Fubini–Study current $\gamma (L,h_{L})$ is defined in this case by formula (2.32). We say that the metric $h_L$ is (semi)positively curved if $c_1(L,h)$ is a positive current. Equivalently, the weights $\varphi _\alpha $ can be chosen to be plurisubharmonic (psh) functions. In this case, the Bergman kernel function can be defined by formula (2.7) (see [Reference Coman and Marinescu15]), $\log P(x,x)$ is still locally integrable [Reference Coman and Marinescu15, Lemma 3.2] and formula (1.4) holds. Adapting the proof of Theorem 1.1 to this context easily yields:

Theorem 2.13. Let $(X,J,\Theta )$ be a connected complex Hermitian manifold, and let $(L,h_L)$ be a holomorphic line bundle on X with a singular Hermitian metric $h_L$ so that $c_1(L,h_L)$ is a positive current. Assume that $H^0_{(2)}(X,L)\neq \{0\}$ . Then as $(1,1)$ -currents on X we have $\mathbb {E}[[\operatorname {\mathrm {Div}}(\psi ^{S}_{\eta })]]=\gamma (L,h_{L})$ .

Example 2.14 (Bargmann–Fock space: flat Gaussian holomorphic function).

Let L be the trivial line bundle on $\mathbb {C}^{n}$ , but we equip it with the Hermitian metric $h_{L}$ such that $|1|^{2}_{h_{L},z}=e^{-|z|^{2}}$ , $z\in \mathbb {C}^{n}$ . In this case, $R^{L}=\sum _{j=1}^{n}dz_{j}\wedge d\bar {z}_{j}$ . We endow $\mathbb {C}^{n}$ with the flat metric $\Theta =\frac {\sqrt {-1}}{2\pi }\sum _{j=1}^{n}dz_{j}\wedge d\bar {z}_{j}$ , so $\mathrm {dV}_{\Theta }=\pi ^{-n}\Pi _{j=1}^{n}dx_{j}\wedge dy_{j}$ . For a multi-index $\alpha =(\alpha _{1},\ldots ,\alpha _{n})\in \mathbb {N}^{n}$ , we write

(2.42) $$ \begin{align} S_{\alpha}(z)=\frac{z_{1}^{\alpha_{1}}\ldots z_{n}^{\alpha_{n}}}{\sqrt{\alpha_{1}!\ldots\alpha_{n}!}}. \end{align} $$

A straightforward calculation then confirms that $\{S_{\alpha }\}_{\alpha \in \mathbb {N}^{n}}$ forms an orthonormal basis of $H^{0}_{(2)}(\mathbb {C}^{n},L)$ . In this case, we have

(2.43) $$ \begin{align} P(z,z)=\sum_{\alpha\in \mathbb{N}^{n}} \frac{|z_{1}|^{2\alpha_{1}}\ldots |z_{n}|^{2\alpha_{n}}}{\alpha_{1}!\ldots\alpha_{n}!}e^{-|z|^{2}}=1. \end{align} $$

Denoting by $\eta =(\eta _{\alpha })_{\alpha \in \mathbb {N}^{n}}$ a family of i.i.d. standard complex Gaussian random variables, we define the standard Gaussian random holomorphic function on $\mathbb {C}^{n}$ as

(2.44) $$ \begin{align} \psi^{S}_{\eta}=\sum_{\alpha\in\mathbb{N}^{n}} \eta_{\alpha}S_{\alpha}. \end{align} $$

By Theorem 1.1, we have

(2.45) $$ \begin{align} \mathbb{E}[[\operatorname{\mathrm{Div}}(\psi^{S}_{\eta})]]=\gamma(L, h_{L})=\frac{\sqrt{-1}}{2\pi}\sum_{j=1}^{n}dz_{j}\wedge d\bar{z}_{j}. \end{align} $$

3 Equidistribution and large deviation for high tensor powers of line bundles

In the sequel, assume that $\eta =\{\eta _{j}\}_{j\in \mathbb {N}}$ is a sequence of i.i.d. standard complex Gaussian random variables, note that $\mathrm {Var}(\eta _{1})=1$ .

In this section, we consider the setting of Subsection 1.2: let $(X,J,\Theta )$ be a connected complex Hermitian manifold and let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ . We assume Condition 1.2 holds. Let $\dot {R}^{L}\in \mathrm {End}(T^{(1,0)}X)$ such that $x\in X$ , for $u,v\in T_{x}^{(1,0)}X$ ,

(3.1) $$ \begin{align} R^{L}_{x}(u,v)=g^{TX}_{x}(\dot{R}^{L}u,v). \end{align} $$

Set $a_{0}(x):=\det \dot {R}^{L}_{x}$ a smooth function on X, then by condition (1.6), $ a_{0}(x)\geq \varepsilon ^{n}$ .

3.1 Equidistribution of zeros of Gaussian random holomorphic sections

We consider the sequence of Hilbert spaces $H^{0}_{(2)}(X,L^{p})$ , $p\in \mathbb {N}$ large. Set

(3.2) $$ \begin{align} d_{p}=\dim H^{0}_{(2)}(X,L^{p})\in \mathbb{N}\cup\{\infty\}. \end{align} $$

We equip $L^{p}$ with the induced Hermitian metric $h_{p}:=h^{\otimes p}_{L}$ . Let $P_{p}$ denote the orthogonal projection from $\mathcal {L}^{2}(X,L^{p})$ onto $H^{0}_{(2)}(X,L^{p})$ , and let $P_{p}(x,y)$ denote the corresponding Bergman kernel on X with respect to $\mathrm {dV}(x)=\frac {\Theta ^{n}}{n!}$ .

By Ma–Marinescu [Reference Ma and Marinescu56, Theorems 4.1.1 & 6.1.1], we have the following results on the asymptotics of on-diagonal Bergman kernels: There exist coefficients $\boldsymbol {b}_{r}\in \mathscr {C}^{\infty }(X,\mathbb {R})$ , $r\in \mathbb {N}$ , such that for any compact subset $K\subset X$ , any $k,\ell \in \mathbb {N}$ , there exists $C_{k,\ell ,K}>0$ such that for $p\in \mathbb {N}^{\ast }$ ,

(3.3) $$ \begin{align} \left|\frac{1}{p^{n}}P_{p}(x,x)-\sum_{r=0}^{k}\boldsymbol{b}_{r}(x)p^{-r}\right|{}_{\mathscr{C}^{\ell}(K)}\leq C_{k,\ell,K} p^{-k-1}, \end{align} $$

where the $\mathscr {C}^{\ell }$ -norm is induced by $g^{TX}$ . In particular, we have

$$ \begin{align*}\boldsymbol{b}_{0}(x)=\det\left(\frac{\dot{R}^{L}_{x}}{2\pi}\right)=\frac{a_{0}(x)}{(2\pi)^{n}},\end{align*} $$

and $\boldsymbol {b}_{1}$ also has an explicit formula given as in [Reference Ma and Marinescu56, (4.1.9)],

(3.4) $$ \begin{align} \boldsymbol{b}_{1}(x)=\frac{1}{8\pi}\det\left(\frac{\dot{R}^{L}_{x}}{2\pi}\right)\left(\boldsymbol{r}^{X}_{\omega}(x)-2\Delta_{\omega} \big(\log(\det \dot{R}_x^{L})\big) \right), \end{align} $$

where $\boldsymbol {r}^{X}_{\omega }$ , $\Delta _{\omega }$ denote the scalar curvature, and the Bochner Laplacian associated to the Riemannian metric $g^{TX}_{\omega }(\cdot ,\cdot ):=c_{1}(L,h_{L})(\cdot , J\cdot )$ on X. An explicit formula for $\boldsymbol {b}_{2}$ can be found in [Reference Ma and Marinescu59, (0.14) & Remark 0.5].

For $p\in \mathbb {N}_{>0}$ , let $\psi ^{S_{p}}_{\eta }$ be a standard Gaussian random holomorphic section constructed from $H^{0}_{(2)}(X,L^{p}),$ that is, for $\{S^{p}_{j}\}_{j=1}^{d_{p}}$ an orthonormal basis of $H^{0}_{(2)}(X,L^{p})$ with respect to the $\mathcal {L}^{2}$ -metric, and set

(3.5) $$ \begin{align} \psi^{S_{p}}_{\eta}=\sum_{j=1}^{d_{p}}\eta_{j} S^{p}_{j}. \end{align} $$

Theorem 3.1. Let $(X,J,\Theta )$ be a complex Hermitian manifold, and let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ , such that Condition 1.2 holds. Then we have:

  1. (i) There exist functions $\boldsymbol {c}_{r}\in \mathscr {C}^{\infty }(X,\mathbb {R})$ , $r\in \mathbb {N}$ , such that for any compact subset $K\subset X$ and $k,\ell \in \mathbb {N}$ , there exists a $p_{k,\ell ,K}\in \mathbb {N}$ , such that for all $p\geq p_{k,\ell ,K}$ , $\mathbb {E}[[\operatorname {\mathrm {Div}}(\psi ^{S_{p}}_{\eta })]]$ is a smooth real $(1,1)$ -form, and the following expansion holds in the $\mathscr {C}^{\ell }(K)$ -norm induced from $g^{TX}$ and associated Levi–Civita connection,

    (3.6) $$ \begin{align} \frac{1}{p}\mathbb{E}[[\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]]= c_{1}(L,h_{L})+ \frac{\sqrt{-1}}{2\pi}\sum_{r=0}^{k} \frac{\partial\overline{\partial}\boldsymbol{c}_{r}(x)}{p^{r+1}}+\mathcal{O}_{K,\ell,k}(p^{-k-2}). \end{align} $$

    Moreover, each function $\boldsymbol {c}_{r}(x)$ for $r\geq 1$ is a polynomial in terms of the coefficients $\boldsymbol {b}_{s}/\boldsymbol {b}_{0}$ , $s\leq r$ , given in (3.3). We have

    (3.7) $$ \begin{align} \begin{aligned} &\boldsymbol{c}_{0}(x)=\log \boldsymbol{b}_{0}(x),\; \boldsymbol{c}_{1}(x)=\frac{\boldsymbol{b}_{1}(x)}{\boldsymbol{b}_{0}(x)}=\frac{1}{8\pi}\left(\boldsymbol{r}^{X}_{\omega}-2\Delta_{\omega} \big(\log(\det \dot{R}^{L})\big) \right),\\ &\boldsymbol{c}_{2}(x)=\frac{\boldsymbol{b}_{2}(x)}{\boldsymbol{b}_{0}(x)}-\frac{1}{2}\left(\frac{\boldsymbol{b}_{1}(x)}{\boldsymbol{b}_{0}(x)}\right)^2. \end{aligned} \end{align} $$
  2. (ii) As $p\rightarrow +\infty $ , we have the weak convergence

    (3.8) $$ \begin{align} \frac{1}{p}\mathbb{E}[[\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]]\rightarrow c_{1}(L,h_{L}) \end{align} $$
    of $(1,1)$ -currents, that is, for any $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ , as $p\rightarrow +\infty $ ,
    (3.9) $$ \begin{align} \Big \langle\frac{1}{p}\mathbb{E}[[\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]],\varphi\Big \rangle\rightarrow \langle c_{1}(L,h_{L}),\varphi\rangle. \end{align} $$

Proof. It is clear that the second part (ii) of this theorem follows directly from the expansion (3.6). We only need to prove the first part (i). By Theorem 1.1 for $L^{p}$ ,

(3.10) $$ \begin{align} \mathbb{E}[[\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]]=\gamma(L^{p},h_{p}), \end{align} $$

where $\gamma (L^{p}, h_{p})$ is the corresponding Fubini–Study current defined via formula (1.4). To obtain the expansion (3.6), it is enough to study the expansion of $\frac {\sqrt {-1}}{2\pi p}\partial \overline {\partial }\log P_{p}(x,x)$ in p.

For this purpose, we apply the expansion (3.3) to $P_{p}(x,x)$ . Note that $\boldsymbol {b}_{0}$ is uniformly bounded from below by a positive constant. Then for a given compact subset K of X, and $k,\ell \in \mathbb {N}$ , we take the expansion (3.3) but with $\mathscr {C}^{\ell +2}(K)$ -norm. We conclude that, for all sufficiently large p, we have $P_{p}(x,x)>0$ , and

(3.11) $$ \begin{align} \big|\sum_{r=1}^{k}\boldsymbol{b}_{r}(x)p^{-r}+C_{k,\ell+2,K} p^{-k-1}\big|\leq \frac{1}{2}\boldsymbol{b}_{0}(x). \end{align} $$

Then $\frac {\sqrt {-1}}{2\pi p}\partial \overline {\partial }\log P_{p}(x,x)$ is a smooth $(1,1)$ -form on K, and the expansion (3.6) is the consequence of the following identity and the Taylor expansion of $\log (1+x)$ with $|x|<1$ ,

(3.12) $$ \begin{align} \begin{aligned} \log P_{p}(x,x)&=n\log p+\log\big(\sum_{r=0}^{k}\boldsymbol{b}_{r}(x)p^{-r}+R_{k}(p)\big)\\ &=n\log p+\log \boldsymbol{b}_{0}(x)+\log\big(1+\sum_{r=1}^{k}p^{-r}\frac{\boldsymbol{b}_{r}(x)}{\boldsymbol{b}_{0}(x)}+R_{k}(p)\big), \end{aligned} \end{align} $$

where $R_{k}(p)=\mathcal {O}(p^{-k-1})$ denotes the remainder term. In particular, we get formulas (3.7).

The convergence in the limit (3.8) can be improved by imposing further geometric assumptions, for instance the assumption of bounded geometry as in [Reference Ma and Marinescu60].

Definition 3.2 (Bounded geometry).

We say that $(X,J,\Theta )$ , $(L,h_{L})$ have bounded geometry if J, $g^{TX}$ , $R^{L}$ and their derivatives of any order are uniformly bounded on X in the norm induced by $g^{TX}$ , and the injective radius of $(X,g^{TX})$ is strictly positive.

One important class of examples are Galois coverings of compact projective manifolds M endowed with the pull-back of a positive holomorphic line bundle on M.

We recall the following results proved in [Reference Ma and Marinescu60, Theorem 3].

Theorem 3.3 [Reference Ma and Marinescu60].

Let $(X,J,\Theta )$ be a complex Hermitian manifold, and let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ and Condition 1.2 holds. In addition, we assume that they have bounded geometry. Then for any $\ell \in \mathbb {N}$ , we have the following expansion holding in the $\mathscr {C}^{\ell }$ -norm induced from $g^{TX}$ and the associated Levi–Civita connection on X,

(3.13) $$ \begin{align} P_{p}(x,x)=\frac{a_{0}(x)}{(2\pi)^{n}}p^{n}+\mathcal{O}(p^{n-1}) \end{align} $$

Moreover, there exists $p_{0}\in \mathbb {N}$ such that for all $p>p_{0}$ , X is holomorphically convex with respect to the bundle $L^{p}$ and $H^{0}_{(2)}(X,L^{p})$ separates points and gives local coordinates on X.

As a consequence, we get the following results.

Proposition 3.4. Assume the same geometric hypothesis as in Theorem 3.3. Writing $\psi ^{S_{p}}_{\eta }$ for the Gaussian random section constructed from $H^{0}_{(2)}(X,L^{p})$ , then for sufficiently large p, $\mathbb {E}[[\operatorname {\mathrm {Div}}(\psi ^{S_{p}}_{\eta })]]$ is a smooth $(1,1)$ -form on X. Then for any $\ell \in \mathbb {N}$ , we have the following convergence in the $\mathscr {C}^{\ell }$ -norm on X

(3.14) $$ \begin{align} \frac{1}{p}\,\mathbb{E}\left[[\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]\right]\rightarrow c_{1}(L,h_{L}), \:\:\text{as } p\rightarrow +\infty. \end{align} $$

Remark 3.5. Note that under the assumption of bounded geometry and for X noncompact, we have $d_{p}=\infty $ , $p\gg 0$ .

Example 3.6 (Scaled Bargmann–Fock spaces).

We consider the line bundle $(L,h_{L})$ on $\mathbb {C}^{n}$ from Example 2.14, which satisfies the above assumptions. For $p\geq 1$ , an orthonormal basis of $H^{0}_{(2)}(\mathbb {C}^{n},L^{p})$ is given by the family

(3.15) $$ \begin{align} S^{p}_{\alpha}(z)=p^{\frac{n}{2}}S_{\alpha}(\sqrt{p}z), \quad \alpha\in\mathbb{N}^{n}. \end{align} $$

Then the Bergman kernel function is given

(3.16) $$ \begin{align} P_{p}(z,z)\equiv p^{n}. \end{align} $$

Recall the flat Gaussian random holomorphic function $\psi ^{S}_{\eta }$ on $\mathbb {C}$ is defined by the equality (1.5). Then for $p\geq 1$ , we have

(3.17) $$ \begin{align} \psi^{S_{p}}_{\eta}(z)=p^{n/2}\psi^{S}_{\eta}(\sqrt{p}z). \end{align} $$

A direct computation then shows that

(3.18) $$ \begin{align} \frac{1}{p}\mathbb{E}[[\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]]=\mathbb{E}[[\operatorname{\mathrm{Div}}(\psi^{S}_{\eta})]]=\frac{\sqrt{-1}}{2\pi}\sum_{j=1}^{n}dz_{j}\wedge d\bar{z}_{j}. \end{align} $$

Theorem 3.7. Let $(L,h_{L})$ and $(X,\Theta )$ be as in Theorem 3.1. For any given test form $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ , we have

(3.19) $$ \begin{align} \mathbb{P}\Big(\lim_{p\rightarrow+\infty} \frac{1}{p}\big\langle [\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})],\varphi\big\rangle= \langle c_{1}(L,h_{L}),\varphi\rangle\Big)=1. \end{align} $$

Proof. To prove this theorem, we mainly follow the arguments from proof of [Reference Ma and Marinescu56, Theorem 5.3.3], and the possibility of infinite dimension does not lead to complications in this setting. Fix a nontrivial test form $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ . Note that from the proof of Theorem 3.1, we have the convergence

(3.20) $$ \begin{align} \lim_{p\rightarrow\infty}\Big\langle \frac{1}{p}\gamma(L^{p},h_{p}),\varphi\Big\rangle=\langle c_{1}(L,h_{L}),\varphi\rangle. \end{align} $$

Defining the random variable

(3.21) $$ \begin{align} Y_{p}=\frac{1}{p}\Big\langle [\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]-\gamma(L^{p},h_{p}),\varphi\Big\rangle, \end{align} $$

statement (3.19) is equivalent to proving the almost sure convergence

(3.22) $$ \begin{align} Y_{p}\rightarrow 0. \end{align} $$

For any $x\in \operatorname {\mathrm {supp}} \varphi $ , let $e_{L}(x)$ denote a unit vector of $(L_{x},h_{L,x})$ . Set

(3.23) $$ \begin{align} b_{p}(x)=(P_{p}(x,x)^{-1/2}S^{p}_{j}(x)/e_{L}^{\otimes p}(x))\in\ell^{2}(\mathbb{C}). \end{align} $$

Then $\eta \cdot b_{p}(x)$ is a standard complex Gaussian variable. The covariance matrix of the Gaussian vector $(\eta \cdot b_{p}(x),\eta \cdot b_{p}(y))$ depends smoothly on $(x,y)\in \operatorname {\mathrm {supp}}\varphi \times \operatorname {\mathrm {supp}}\varphi $ .

For $v=(v_{1},v_{2})\in \mathbb {C}^{2}$ with $\|v\|=1$ , we consider the integral

(3.24) $$ \begin{align} \rho(v):=\frac{1}{4\pi^{2}}\int_{\mathbb{C}^{2}}e^{-\frac{1}{2}(|z_{1}|^{2}+|z_{2}|^{2})}\big|\log|z_{1}|\cdot\log|v_{1}z_{1}+v_{2}z_{2}|\big|\mathrm{dV}(z). \end{align} $$

The computations in [Reference Ma and Marinescu56, Eqs. (5.3.13) to (5.3.15)] then show that

(3.25) $$ \begin{align} C:=\sup_{v\in\mathbb{C}^{2},\|v\|=1 }\rho(v)<\infty, \end{align} $$

so for $x,y\in \operatorname {\mathrm {supp}}\varphi $ we have

(3.26) $$ \begin{align} \mathbb{E}\Big[\Big|\log\big|P_{p}(x,x)^{-1/2}\sum_{j}\eta_{j}S^{p}_{j}(x) \big|_{h_{p}}\log\big|P_{p}(y,y)^{-1/2} \sum_{j}\eta_{j}S^{p}_{j}(y)\big|_{h_{p}}\Big|\Big]\leq C. \end{align} $$

Note that

(3.27) $$ \begin{align} \mathbb{E}[|Y_{p}|^{2}]=\frac{1}{p^{2}}\mathbb{E}\big[\left|\langle [\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})],\varphi\rangle\right|{}^{2}\big]-\frac{1}{p^{2}}\big|\langle \gamma(L^{p},h_{p}),\varphi\rangle\big|^{2}. \end{align} $$

Then by relations (2.29) and (3.26) and the Fubini–Tonelli theorem we infer that

(3.28) $$ \begin{align} \begin{aligned} \mathbb{E}[|Y_{p}|^{2}] & =\frac{1}{\pi^{2}p^{2}}\int_{X\times X} \big(\partial\overline{\partial}\varphi(x)\big) \big(\overline{\partial\overline{\partial}\varphi(y)}\big)\\ &\quad \mathbb{E}\Big[\log\big|P_{p}(x,x)^{-1/2}\sum_{j}\eta_{j}S^{p}_{j}(x) \big|_{h_{p}}\log\big|P_{p}(y,y)^{-1/2} \sum_{j}\eta_{j}S^{p}_{j}(y)\big|_{h_{p}}\Big]. \end{aligned} \end{align} $$

By estimate (3.26), we conclude $\mathbb {E}[|Y_{p}|^{2}]=\mathcal {O}(\frac {1}{p^{2}})$ . Hence, $\mathbb {E}[\sum _{p\geq 1}|Y_{p}|^{2}]=\sum _{p\geq 1}\mathbb {E}[|Y_{p}|^{2}]<\infty $ , thus $Y_{p}\rightarrow 0$ almost surely.

Corollary 3.8. If $\Theta $ is a Kähler form and $\int _{X} c_{1}(L,h_{L})\wedge \Theta ^{n-1}<\infty $ , or if $\int _{X} c_{1}(L,h_{L})^{n}<\infty $ , then

(3.29) $$ \begin{align} \mathbb{P}\Big(\lim_{p\rightarrow+\infty} \frac{1}{p} [\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]= c_{1}(L,h_{L})\Big)=1, \end{align} $$

where the limit is taken with respect to the weak convergence of $(1,1)$ -currents on X.

Proof. Let’s assume at first that $\Theta $ is a Kähler form with $\int _{X} c_{1}(L,h_{L})\wedge \Theta ^{n-1}<\infty $ . For all $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ , $s_{p}\in H^{0}(X,L^{p})$ , due to the positivity of the current $[\operatorname {\mathrm {Div}}(s_{p})]$ , we have

(3.30) $$ \begin{align} |\langle[\operatorname{\mathrm{Div}}(s_{p})],\varphi\rangle|\leq |\varphi|_{\mathscr{C}^{0}(X)}\langle[\operatorname{\mathrm{Div}}(s_{p})], \Theta^{n-1}\rangle=p|\varphi|_{\mathscr{C}^{0}(X)}\int_{X}c_{1}(L,h_{L})\wedge\Theta^{n-1}. \end{align} $$

The last identity follows from the Poncaré–Lelong formula. We may set $C=\int _{X}c_{1}(L,h_{L})\wedge \Theta ^{n-1}<\infty $ , then

(3.31) $$ \begin{align} \frac{1}{p}|\langle[\operatorname{\mathrm{Div}}(s_{p})],\varphi\rangle|\leq C|\varphi|_{\mathscr{C}^{0}(X)}. \end{align} $$

By considering a countable $\mathscr {C}^{0}$ -dense family of $\varphi $ ’s in $\Omega ^{n-1,n-1}_{0}(X)$ , and applying Theorem 3.7, we get the equality (3.29).

If $\int _{X} c_{1}(L,h_{L})^{n}<\infty $ , we use $c_{1}(L,h_{L})$ instead of $\Theta $ in the above arguments to get the equality (3.29).

Remark 3.9. (1) The additional assumptions in the above corollary are necessary in our approach to the proof of the equality (3.29); however, it is an interesting question whether or not these extra assumptions could be removed.

(2) In the first assumption of Corollary 3.8, we can replace the Kähler condition on $\Theta $ by the condition that $\Theta $ is Gauduchon, which means $\partial \bar {\partial }\Theta ^{n-1}=0$ .

(3) Due to the uniform estimate (3.30) for all holomorphic sections $s_{p}$ with a given p, we can improve Corollary 3.8 as follows: Equation (3.29) holds true when the limit is taken as the weak convergence in the sense of current of order $0$ , that is, for any continuous $(n-1,n-1)$ -form $\varphi $ with compact support on X, we have

(3.32) $$ \begin{align} \lim_{p\rightarrow+\infty}\langle \frac{1}{p} [\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})],\varphi\rangle= \langle c_{1}(L,h_{L}),\varphi\rangle. \end{align} $$

Remark 3.10. For each $p\in \mathbb {N}_{>0}$ , we can take a sequence of i.i.d. standard complex Gaussian random variables $\eta ^{p}=\{\eta ^{p}_{j}\}_{j=1}^{d_{p}}$ and assume that they are mutually independent for different $p.$ We define the flat Gaussian random sections

(3.33) $$ \begin{align} \psi^{S_{p}}_{\eta^{p}}=\sum_{j=1}^{d_{p}}\eta^{p}_{j}S^{p}_{j}, \end{align} $$

where $S_{p}=\{S^{p}_{j}\}_{j=1}^{d_{p}}$ is an orthonormal (Hilbert) basis of $H^{0}_{(2)}(X,L^{p})$ . Then the statements in Theorems 3.1, 3.7, Proposition 3.4 and Corollary 3.8 still hold true for the sequence of random sections $\psi ^{S_{p}}_{\eta ^{p}}, p\geq 1$ .

3.2 Large deviation estimates and hole probability: proofs of Theorems 1.4 and 1.5

In this subsection, we study the large deviation estimates for random zeros in a given domain with respect to the high tensor powers as in [Reference Shiffman, Zelditch and Zrebiec70], [Reference Dinh, Marinescu and Schmidt25] and [Reference Drewitz, Liu and Marinescu27]. In particular, we prove Theorems 1.4 and 1.5. A key intermediate result in the approach to the above theorems is the proposition as follows, whose proof is deferred to the next subsection.

Proposition 3.11. Let $(X,J,\Theta )$ be a connected complex Hermitian manifold, and let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ such that Condition 1.2 holds. Let U be a relatively compact open subset in X. For any $\delta>0$ , there exists $C_{U,\delta }>0$ such that for all $ p\gg 0$ ,

(3.34) $$ \begin{align} \mathbb{P}\Big( \int_{U}\Big|\log{\big|\psi^{S_{p}}_{\eta}(x)\big|_{h_{p}}}\Big| \,\mathrm{dV}(x)\geq \delta p \Big)\leq e^{-C_{U,\delta}p^{n+1}}. \end{align} $$

Proof of Theorem 1.4.

The Poincaré–Lelong formula (2.29) shows that

(3.35) $$ \begin{align} \frac{\sqrt{-1}}{\pi}\partial\overline{\partial} \log{|\psi^{S_{p}}_{\eta}|_{h_{p}}}=[\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})]-pc_{1}(L,h_L) \end{align} $$

as an identity of $(1,1)$ -currents on $X.$ Now, fix $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ , and fix a relatively compact open subset $U\subset X$ such that $\mathrm {supp}\ \varphi \subset U$ . Then

(3.36) $$ \begin{align} \begin{aligned} \Big(\frac{1}{p}[\operatorname{\mathrm{Div}}(s_{p})],\varphi\Big) -\int_{X}c_{1}(L,h_L)\wedge \varphi=\frac{\sqrt{-1}}{p\pi}\int_{X} \log{|\psi^{S_{p}}_{\eta}|_{h_{p}}}\,\partial\overline{\partial}\varphi. \end{aligned} \end{align} $$

Since $\varphi $ has a compact support in U, so has $\partial \overline {\partial }\varphi $ . Set

(3.37) $$ \begin{align} S_{\varphi}=\max_{x\in U}\left|\frac{\sqrt{-1}\partial \overline{\partial}\varphi(x)}{\mathrm{dV}(x)}\right|. \end{align} $$

We can and we may assume that $S_{\varphi }>0$ . Then

(3.38) $$ \begin{align} \begin{aligned} \left|\frac{\sqrt{-1}}{p\pi}\int_{X}\log{|\psi^{S_{p}}_{\eta}|_{h_{p}}}\, \partial\overline{\partial}\varphi\right|\leq \frac{S_{\varphi}}{p\pi}\int_{U} \big|\log{|\psi^{S_{p}}_{\eta}(x)|_{h_{p}}}\big|\,\mathrm{dV}(x). \end{aligned} \end{align} $$

Applying Proposition 3.11 to right-hand side of inequality (3.38) we get estimate (1.10).

Proof of Theorem 1.5.

Estimate (1.14) is a direct consequence of estimate (1.13) by taking $\delta =n\operatorname {\mathrm {Vol}}^{L}_{2n}(U)$ . Hence, it is sufficient to prove estimate (1.13). For this purpose, let $\chi _{U}$ denote the characteristic function of U on X. Let $\delta>0$ be arbitrary, and take $\psi _{1}$ , $\psi _{2}\in \mathcal {C}^{\infty }_{0}(X,\mathbb {R})$ such that $0\leq \psi _{1}\leq \chi _{U}\leq \psi _{2}\leq 1,$ and

(3.39) $$ \begin{align} \int_{X} \psi_{1}\frac{c_{1}(L,h_{L})^{n}}{n!}\geq \operatorname{\mathrm{Vol}}^{L}_{2n}(U)-\delta, \quad \int_{X} \psi_{2}\frac{c_{1}(L,h_{L})^{n}}{n!}\leq \operatorname{\mathrm{Vol}}^{L}_{2n}(U)+\delta\,. \end{align} $$

Note that the existence of such functions is guaranteed by the assumption that $\partial U$ has measure $0$ with respect to $\mathrm {dV}$ , hence also to $\frac {1}{n!}c_{1}(L,h_{L})^{n}$ . For $j\in \{1,2\},$ set $\varphi _{j}=\frac {1}{(n-1)!}{\psi _{j}}c_{1}(L,h_{L})^{n-1}$ . By applying Theorem 1.4 to $\varphi _{j}$ separately, we get exactly estimate (1.13).

3.3 Proof of Proposition 3.11

Let $U\subset X$ be a relatively compact open subset. For $s_{p}\in H^{0}(X, L^{p})$ , we set

(3.40) $$ \begin{align} \mathcal{M}^{U}_{p}(s_{p})=\sup_{x\in U}|s_{p}(x)|_{h_{p}}<+\infty. \end{align} $$

Before proving Proposition 3.11, we need to investigate the probabilities for both, $\mathcal {M}^{U}_{p}(\psi ^{S_{p}}_{\eta })$ taking atypically large and small values, respectively.

Proposition 3.12. For any $\delta>0$ , there exists a constant $C_{U,\delta }>0$ such that for $p\in \mathbb {N}_{>1}$ ,

(3.41) $$ \begin{align} \mathbb{P}\big (\mathcal{M}^{U}_{p}(\psi^{S_{p}}_{\eta})\geq e^{\delta p}\big)\leq e^{-\delta p^{n+1}+C_{U,\delta}p^{n}\log{p}}\,. \end{align} $$

Proof. The basic idea of the proof is that the local sup-norm of a holomorphic function is bounded by its local $\mathcal {L}^{2}$ -norm as in estimate (2.9). We fix $\delta>0$ and let $r>0$ be sufficiently small so that we can choose a finite set of points $\{x_{j}\}_{j=1}^{\ell }\subset U$ such that the geodesic open balls $B^{X}(x_{j},r)$ , $j=1,\ldots , \ell $ form an open covering of $\overline {U}$ . Since r is sufficiently small, then we can assume that each larger ball $B^{X}(x_{j},2r)$ lies in a complex chart (hence viewed as an open subset of $\mathbb {C}^{n}$ ), and that for each j, we can fix a local holomorphic frame $e_{L,j}$ of L on a neighborhood of $B^{X}(x_{j}, 2r)$ with $\sup _{x\in B^{X}(x_{j},2r)} |e_{L,j}(x)|_{h_{L}}=1$ . Set

(3.42) $$ \begin{align} \nu=\min\big\{ \inf_{x\in B^{X}(x_{j},2r)}|e_{L,j}(x)|_{h_{L}} \,:\,j=1 , \ldots, \ell\big\}. \end{align} $$

It is clear that $0<\nu \leq 1$ . By fixing r small enough, we can and do assume that

(3.43) $$ \begin{align} -\log{\nu}\leq \frac{\delta}{4}\,\cdot \end{align} $$

As in estimate (2.9), since U is relatively compact, there exists a constant $C>0$ such that for each $j=1,\ldots , \ell $ , if f is a holomorphic function on a neighborhood of $B^{X}(x_{j},2r)$ , then

(3.44) $$ \begin{align} \sup_{x\in B^{X}(x_{j},r)}|f(x)|\leq C\|f\|_{\mathcal{L}^{2}(B^{X}(x_{j},2r))}, \end{align} $$

where the volume form $\mathrm {dV}(x)$ on X is used in the norm $\|\cdot \|_{\mathcal {L}^{2}(B^{X}(x_{j},2r))}$ . Note that the choices of $x_{j}$ , r, $\ell $ and the constants $\nu $ , C are independent of the tensor power p. Set $\widetilde {U}=\cup _{j} B^{X}(x_{j},2r)\supset U$ . For $p\in \mathbb {N}, s_{p}\in H^{0}(X,L^{p})$ , on each $B^{X}(x_{j},2r)$ , we write

(3.45) $$ \begin{align} s_{p}|_{B^{X}(x_{j},2r)}=f_{j}e_{L,j}^{\otimes p}, \end{align} $$

where $f_{j}$ is a holomorphic function on the chart in $\mathbb {C}^{n}$ corresponding to $B^{X}(x_{j},2r)$ . Then we have

(3.46) $$ \begin{align} \begin{aligned} \mathcal{M}^{U}_{p}(s_{p})=\sup_{x\in U}|s_{p}(x)|_{h_{p}}&\leq \max_{j} \sup_{x\in B^{X}(x_{j},r)}|f_{j}(x)|\\ &\leq C\max_{j} \{\|f_{j}\|_{\mathcal{L}^{2}(B^{X}(x_{j},2r))}\}\\ &\leq \frac{C}{\nu^{p}}\max_{j}\{\|s_{p}\|_{\mathcal{L}^{2}(B^{X}(x_{j},2r),L^{p})}\}\\ &\leq \frac{C}{\nu^{p}}\|s_{p}\|_{\mathcal{L}^{2}(\widetilde{U},L^{p})}. \end{aligned} \end{align} $$

The next step is to estimate the quantity $\mathbb {E}[\|\psi ^{S_{p}}_{\eta }\|^{2p^{n}}_{\mathcal {L}^{2}(\widetilde {U},L^{p})}]$ for $p\geq 2$ . Applying Hölder’s inequality with $\frac {1}{p^{n}}+\frac {p^{n}-1}{p^{n}}=1$ , we get

(3.47) $$ \begin{align} \mathbb{E}\big[\|\psi^{S_{p}}_{\eta}\|^{2p^{n}}_{\mathcal{L}^{2}(\widetilde{U},L^{p})}\big] \leq \operatorname{\mathrm{Vol}}(\widetilde{U})^{p^{n}-1}\mathbb{E}\Big[\int_{\widetilde{U}} |\psi^{S_{p}}_{\eta}(x)|^{2p^{n}}_{h_{p}}(x)\mathrm{dV}\Big]. \end{align} $$

As in estimate (3.44), on a neighborhood of $B^{X}(x_{j}, 2r)$ , write

(3.48) $$ \begin{align} S^{p}_{i}=f^{p}_{i}e_{L,i}^{\otimes p}. \end{align} $$

If $x\in B^{X}(x_{j},2r)$ , set

(3.49) $$ \begin{align} F_{j}(x)=\sum_{i=1}^{d_{p}} \eta_{i} f^{p}_{i}(x). \end{align} $$

Then $F_{j}(x)$ is a complex Gaussian random variable with (total) variance $\sum _{i=1}^{d_{p}}|f^{p}_{i}(x)|^{2}$ . By our assumption on the local frame $e_{L,j}$ , we get

(3.50) $$ \begin{align} \sum_{i=1}^{d_{p}}|f^{p}_{i}(x)|^{2}\leq \frac{1}{\nu^{2p}}P_{p}(x,x). \end{align} $$

Then we have

(3.51) $$ \begin{align} \mathbb{E}\big[|F_{j}(x)|^{2p^{n}}\big]=p^{n}! \Big(\sum_{i=1}^{d_{p}}|f^{p}_{i}(x)|^{2}\Big)^{p^{n}}. \end{align} $$

As a consequence, we get that for $x\in \widetilde {U}$ ,

(3.52) $$ \begin{align} \begin{aligned} \mathbb{E}\big[|\psi^{S_{p}}_{\eta}(x)|_{h_{p}}^{2p^{n}}\big]\leq\frac{p^{n}!}{\nu^{2p^{n+1}}}(P_{p}(x,x))^{p^{n}}. \end{aligned} \end{align} $$

Since we are in the context of $\sigma $ -finite measures and the integrands are nonnegative, Tonelli’s theorem applies so that

(3.53) $$ \begin{align} \mathbb{E}\Big[\int_{\widetilde{U}}|\psi^{S_{p}}_{\eta}(x)|^{2p^{n}}_{h_{p}}\mathrm{dV}(x)\Big]\leq \frac{p^{n}!}{\nu^{2p^{n+1}}}\int_{\widetilde{U}} (P_{p}(x,x))^{p^{n}}\mathrm{dV}(x). \end{align} $$

Moreover, by the on-diagonal estimate for the Bergman kernel on a given compact subset, there exists a constant $C_{\widetilde {U}}>0$ (independent of p) such that for $p\in \mathbb {N}$ , $x\in \widetilde {U}$ ,

(3.54) $$ \begin{align} P_{p}(x,x)\leq C_{\widetilde{U}}p^{n}. \end{align} $$

Combining estimate (3.47) with the above inequalities, we infer that

(3.55) $$ \begin{align} \mathbb{E}\Big[\|\psi^{S_{p}}_{\eta}\|^{2p^{n}}_{\mathcal{L}^{2}(\widetilde{U},L^{p})} \Big]\leq \big(C_{\widetilde{U}}\operatorname{\mathrm{Vol}}(\widetilde{U})\big)^{p^{n}}\frac{p^{n}!}{\nu^{2p^{n+1}}}(p^{n})^{p^{n}}. \end{align} $$

By applying estimate (3.46) to $\psi ^{S_{p}}_{\eta }$ , we get

(3.56) $$ \begin{align} \begin{aligned} \mathbb{E}\big[\mathcal{M}^{U}_{p}(\psi^{S_{p}}_{\eta})^{2p^{n}}\big]&\leq \Big(\frac{C}{\nu^{p}}\Big)^{2p^{n}}\mathbb{E}\big[\|\psi^{S_{p}}_{\eta}\|^{2p^{n}}_{\mathcal{L}^{2}(\widetilde{U},L^{p})}\big] \leq \frac{(\widetilde{C}p^{n})^{2p^{n}}}{\nu^{4p^{n+1}}}, \end{aligned} \end{align} $$

where $C>0$ , $\widetilde {C}>0$ are constants independent of p. Then estimate (3.41) follows from Chebyshev’s inequality and the inequality $\frac {1}{\nu }\leq e^{\frac {\delta }{4}}$ from (3.43).

Remark 3.13. The choice to consider the $p^{n}$ -th moment of $\|\psi ^{S_{p}}_{\eta }\|^{2}$ leads to the exponent $p^{n+1}$ in the exponential of the resulting probability estimate. One can consider arbitrary N-th moments to obtain a more general statement on this probability upper bound.

When X is compact, or if X is noncompact but $d_{p}$ is bounded polynomially in p, then the upper bound $Ce^{-cp^{n+1}}$ can be obtained in a simpler way as in [Reference Drewitz, Liu and Marinescu27, Reference Shiffman, Zelditch and Zrebiec70] (and of course with a much sharper upper bound).

Now, we consider the probabilities of small values of $\mathcal {M}^{U}_{p}(\psi ^{S_{p}}_{\eta })$ , and we will adapt the ideas in [Reference Drewitz, Liu and Marinescu27, Reference Shiffman, Zelditch and Zrebiec70]. At first, we introduce a result on the near-diagonal estimate of Bergman kernel. We fix a point $x\in X$ . Let $\{\boldsymbol {f}_{j}\}_{j=1}^{n}$ be an orthonormal basis of $(T_{x}^{1,0}X, g_{x}^{TX}(\cdot ,\overline {\cdot }))$ such that

(3.57) $$ \begin{align} \dot{R}^{L}_{x}\,\boldsymbol{f}_{j}=\mu_{j}(x)\boldsymbol{f}_{j}, \end{align} $$

where $\mu _{j}(x)$ , $j=1, \ldots , n$ , are the eigenvalues of $\dot {R}^{L}_{x}$ (cf. formula (3.1)). Then by the first inequality in condition (1.6), we have

(3.58) $$ \begin{align} \mu_{j}(x)\geq\varepsilon. \end{align} $$

Set $\boldsymbol {e}_{2j-1}=\frac {1}{\sqrt {2}}(\boldsymbol {f}_{j}+\overline {\boldsymbol {f}}_{j})$ , $\boldsymbol {e}_{2j}=\frac {\sqrt {-1}}{\sqrt {2}}(\boldsymbol {f}_{j}-\overline {\boldsymbol {f}}_{j})$ , $j=1, \ldots , n$ . Then they form an orthonormal basis of the (real) tangent vector space $(T_{x}X, g_{x}^{TX})$ . If $v=\sum _{j=1}^{2n} v_{j}\boldsymbol {e}_{j}\in T_{x}X$ , we can write

(3.59) $$ \begin{align} v=\sum_{j=1}^{n}(v_{2j-1}+\sqrt{-1}v_{2j})\frac{1}{\sqrt{2}} \boldsymbol{f}_{j} + \sum_{j=1}^{n}(v_{2j-1}-\sqrt{-1}v_{2j})\frac{1}{\sqrt{2}} \overline{\boldsymbol{f}}_{j}. \end{align} $$

Set $z=(z_{1},\ldots ,z_{n})$ with $z_{j}=v_{2j-1}+\sqrt {-1}v_{2j}$ , $j=1,\ldots , n$ . We call z the complex coordinate of $v\in T_{x}X$ . Then by formula (3.59),

(3.60) $$ \begin{align} \frac{\partial}{\partial z_{j}} = \frac{1}{\sqrt{2}} \boldsymbol{f}_{j} , \; \frac{\partial}{\partial \overline{z}_{j}} = \frac{1}{\sqrt{2}} \overline{\boldsymbol{f}}_{j} \end{align} $$

so that

(3.61) $$ \begin{align} v=\sum_{j=1}^{m}\Big(z_{j} \frac{\partial}{\partial z_{j}} +\overline{z}_{j} \frac{\partial}{\partial \overline{z}_{j}}\Big). \end{align} $$

Note that $|\frac {\partial }{\partial z_{j}}|^{2}_{g^{TX}}=|\frac {\partial }{\partial \overline {z}_{j}}|^{2}_{g^{TX}}=\frac {1}{2}$ . For $v, v'\in T_{x}X$ , let $z, z'$ denote the corresponding complex coordinates. Define a weighted distance function $\Phi ^{TX}_{x}(v,v')$ as follows,

(3.62) $$ \begin{align} \Phi^{TX}_{x}(v,v')^{2}=\sum_{j=1}^{n}\mu_{j}(x)|z_{j}-z^{\prime}_{j}|^{2}. \end{align} $$

For sufficiently small $\delta _{0}>0$ , we identify the small open ball $B^{X}(x,2\delta _{0})$ in X with the ball $B^{T_{x}X}(0, 2\delta _{0})$ in $T_{x}X$ via the geodesic coordinate. Let $\mathrm {dist}(\cdot ,\cdot )$ denote the Riemannian distance of $(X, g^{TX})$ . There exists $C_{1}>0$ such that for $v,v'\in B^{T_{x}X}(0, 2\delta _{0})$ , we have

(3.63) $$ \begin{align} C_{1}\mathrm{dist}(\exp_{x}(v),\exp_{x}(v'))\geq \Phi^{TX}_{x}(v,v')\geq \frac{1}{C_{1}}\mathrm{dist}(\exp_{x}(v),\exp_{x}(v')). \end{align} $$

In particular,

(3.64) $$ \begin{align} \Phi^{TX}_{x}(0,v)\geq \varepsilon^{1/2} \mathrm{dist}(x,\exp_{x}(v)). \end{align} $$

Moreover, if we consider a compact subset $K\subset X$ , the constants $\delta _{0}$ and $C_{1}$ can be chosen uniformly for all $x\in K$ . For $p\in \mathbb {N}$ , $x,y\in X$ , the normalized Bergman kernel is defined as

(3.65) $$ \begin{align} N_{p}(x,y)=\frac{|P_{p}(x,y)|_{h_{p,x}\otimes h_{p,y}^{\ast}}}{\sqrt{P_{p}(x,x)}\sqrt{P_{p}(y,y)}}. \end{align} $$

The following result was proved in [Reference Drewitz, Liu and Marinescu27, Theorem 5.1], where we use essentially the near-diagonal expansion of Bergman kernel from [Reference Ma and Marinescu56, Theorems 4.2.1 & 6.1.1].

Theorem 3.14. Let $(X,J,\Theta )$ be a connected complex Hermitian manifold, and let $(L,h_{L})$ be a holomorphic line bundle on X with a smooth Hermitian metric $h_{L}$ such that Condition 1.2 holds. Let U be a relatively compact open subset of $X.$ Then the following uniform estimates on the normalized Bergman kernel hold for $x,y\in U$ : For $k\geq 1$ and $b>\sqrt {16k/\varepsilon }$ fixed, we have for $p\gg 0\,$ (such that $b\sqrt {(\log {p})/p}\leq 2\delta _{0}$ ) that

(3.66) $$ \begin{align} N_{p}(x,y)=\begin{cases} &\big(1+o(1)\big)\exp\Big(-\dfrac{p}{4}\,\Phi_{x}(0,v')^{2}\Big), \text{ uniformly for } \mathrm{dist}(x,y)\leq b\sqrt{(\log{p})/p},\\ &\qquad\qquad\qquad\qquad\qquad\qquad\qquad\qquad \;\text{with }y=\exp_{x}(v'), v'\in T_{x}X; \\[4pt] &\mathcal{O}(p^{-k}), \;\; \text{uniformly for } \mathrm{dist}(x,y) \geq b\sqrt{(\log{p})/p}. \end{cases} \end{align} $$

Proposition 3.15. There exist constants $C_{U}>0, C_{U}'>0$ such that for all $\delta>0$ and $p\in \mathbb {N}$ ,

(3.67) $$ \begin{align} \mathbb{P}\big(\mathcal{M}^{U}_{p}(\psi^{S_{p}}_{\eta})\leq e^{-\delta p}\big)\leq e^{-C_{U}\delta p^{n+1}+C_{U}' p^{n}\log{p}}\,. \end{align} $$

Proof. For $x\in X$ , we fix some $\lambda _{x}\in L_{x}$ with $|\lambda _{x}|_{h}=1$ and set

(3.68) $$ \begin{align} \xi_{x}=\frac{\langle \lambda_{x}^{\otimes p}, \psi^{S_{p}}_{\eta}(x)\rangle_{h_{p}}}{\sqrt{P_{p}(x,x)}}\,\cdot \end{align} $$

Then $\xi _{x}$ is a complex Gaussian random variable. Moreover, for any two points $x,y\in X$ , we have

(3.69) $$ \begin{align} \big|\mathbb{E}[\xi_{x}\overline{\xi}_{y}]\big| =N_{p}(x,y). \end{align} $$

Then by the asymptotics (3.66), using arguments similar to those used in [Reference Shiffman, Zelditch and Zrebiec70, Subsection 3.2] or the proof of [Reference Drewitz, Liu and Marinescu27, Theorem 1.13], we can prove a more general version of estimate (3.67) as follows: for a sequence of positive numbers $\{\lambda _{p}\}_{p\in \mathbb {N}}$ with $\lambda _p\rightarrow 0$ , we have

(3.70) $$ \begin{align} \mathbb{P}\big(\mathcal{M}^{U}_{p}(\psi^{S_{p}}_{\eta})\leq \lambda_{p}\big)\leq e^{Cp^{n}\log \lambda_{p} +C'p^{n}\log p},\quad p\gg 0. \end{align} $$

Then, for any $\delta>0$ , choosing $\lambda _{p}=e^{-\delta p}$ in estimate (3.70), we recover (3.67).

Combining Propositions 3.12 and 3.15, we arrive at the following.

Corollary 3.16. For any relatively compact open subset $U\subset X$ , and for $\delta>0$ , there exists a constant $C=C(U,\delta )>0$ such that for $p\gg 1$ ,

(3.71) $$ \begin{align} \mathbb{P}\left(\left|\log{\mathcal{M}^{U}_{p}(\psi^{S_{p}}_{\eta})}\right| \geq \delta p\right) \leq e^{-Cp^{n+1}}. \end{align} $$

Proof of Proposition 3.11.

The proof of Proposition 3.11 follows by combining from the arguments in [Reference Shiffman, Zelditch and Zrebiec70, Subsection 4.1] with Corollary 3.16. Here, we just sketch the proof.

For $t>0$ , set

(3.72) $$ \begin{align} \log^{+}{t}=\max\{\log{t},0\},\; \log^{-}{t}:=\log^{+}(1/t)=\max\{-\log{t},0\}. \end{align} $$

Then

(3.73) $$ \begin{align} |\log{t}|=\log^{+}{t}+\log^{-}{t}. \end{align} $$

Let U be a relatively compact nonempty open subset in $X.$ Then for any nonzero holomorphic section $s_{p}\in H^{0}(X,L^{p})$ , we have that $\big |\log |s_{p}|_{h_{p}}\big |$ is integrable on $\overline {U}$ with respect to $\mathrm {dV}$ . We now start with showing that

(3.74) $$ \begin{align} \mathbb{P}\left( \int_{U}\log^{+}{|\psi^{S_{p}}_{\eta}(x)|_{h_{p}}}\,\mathrm{dV}(x)\geq \frac{\delta}{2} p\right)\leq e^{-C_{U,\delta}p^{n+1}}. \end{align} $$

For this purpose, observe that on U we have

(3.75) $$ \begin{align} \log^{+}{|\psi^{S_{p}}_{\eta}|_{h_{p}}}\leq \big|\log\mathcal{M}^{U}_{p}(\psi^{S_{p}}_{\eta})\big|, \end{align} $$

which then supplies us with

(3.76) $$ \begin{align} \mathbb{P}\left( \int_{U}\log^{+}{|\psi^{S_{p}}_{\eta}(x)|_{h_{p}}}\,\mathrm{dV}(x)\geq \frac{\delta}{2} p\right) \leq \mathbb{P}\left( \left|\log\mathcal{M}^{U}_{p}(\psi^{S_{p}}_{\eta})\right|\geq \frac{\delta}{2 \operatorname{\mathrm{Vol}}(U)} p\right), \end{align} $$

where $\operatorname {\mathrm {Vol}}(U)$ denotes the volume of U with respect to $\mathrm {dV}$ . In combination with Corollary 3.16, this immediately implies estimate (3.74). The next step is to prove that

(3.77) $$ \begin{align} \mathbb{P}\left( \int_{U}\log^{-}{|\psi^{S_{p}}_{\eta}(x)|_{h_{p}}}\,\mathrm{dV}(x)\geq \frac{\delta}{2} p\right)\leq e^{-C_{U,\delta}\,p^{n+1}}. \end{align} $$

Suppose that U contains an annulus $B(2,3):=\{z\in \mathbb {C}^{n}:2<|z|<3\}$ (possibly after rescaling of coordinates) and the line bundle L on $B(4):=\{|z|<4\}$ (still contained in U) has a holomorphic local frame $e_{L}$ . We will mainly work on this annulus $B(2,3)$ instead of a coordinate disc. The reason is that, as in [Reference Shiffman, Zelditch and Zrebiec70, p. 1991 in Subsection 4.1], we will need the two-sided positive bounds for the Poisson kernel of the disc and such bounds will hold on the annulus $B(2,3)$ .

Set $\alpha (x)= \log |e_{L}(x)|^{2}_{h_L}$ . We can then write

(3.78) $$ \begin{align} \psi^{S_{p}}_{\eta}=F_{p}e^{\otimes p}_{L}, \end{align} $$

where $F_{p}$ is a random holomorphic function on $B(4)$ . Then

(3.79) $$ \begin{align} \log|\psi^{S_{p}}_{\eta}|_{h_{p}}=\log|F_{p}|+\frac{p}{2}\alpha. \end{align} $$

In the following estimates, each $K_{i}$ , $i \in \mathbb {N},$ denotes a sufficiently large positive constant. Then by relations (3.73) and (3.76), we have

(3.80) $$ \begin{align} \mathbb{P}\left( \int_{B(2,3)}\log^{+}{|F_{p}|}\,\mathrm{dV}\geq K_{1} p \right)\leq e^{-C_{U,K_{1}}\, p^{n+1}}. \end{align} $$

Using the Poisson kernel and the submean inequality for $\log (|F_{p}|)$ (see [Reference Shiffman, Zelditch and Zrebiec70, Subsection 4.1, pp. 1991]), we can improve estimate (3.80) to get

(3.81) $$ \begin{align} \mathbb{P}\left( \int_{B(2,3)}\big|\log{|F_{p}|}\big|\,\mathrm{dV}\geq K_{2} p \right)\leq e^{-C_{U,K_{2}}\, p^{n+1}}. \end{align} $$

From this point, we proceed as in [Reference Shiffman, Zelditch and Zrebiec70, Section 4.1, p. 1992]. For $\delta \in (0,\frac {1}{2}]$ , we fix a grid in the polar coordinate system of $B(2,3)$ so that, by enlarging a bit the grid cells, we obtain an open covering $\{U_j\}_{j=1}^q$ of $B(2,3)$ consisting of small boxes of diameters $\simeq \delta ^{2n+2}$ . Then for a finite set of points $\{z_{j}\}_{j=1}^{q}$ with $z_j\in U_j$ and for all $s_{p}$ , $p\in \mathbb {N}$ , we have

(3.82) $$ \begin{align} \begin{aligned} &-\int_{B(2,3)} \log|s_{p}|_{h_{p}}\mathrm{dV}\\ &\quad \leq -\sum_{j=1}^{q} \mu_{j} \log|s_{p}|_{h_{p}}(z_{j})+K_{3}\delta\int_{B(2,3)} \big|\log{|f_{p}|}\big|\,\mathrm{dV}+p\delta K_{3} \sup_{z\in B(2,3)}|d\alpha(z)|_{g^{T^{\ast}X}}, \end{aligned} \end{align} $$

where the constant $K_{3}$ does not depend on $\delta $ , but the quantities q and $\mu _{j}>0$ only depend on $\delta $ , and we have $\sum _{j=1}^{q}\mu _{j}\simeq 1$ . Applying the above inequality to $\psi ^{S_{p}}_{\eta }$ and $F_{p}$ , we can use Corollary 3.16 for each $U_j$ : except an event of probability $\leq e^{-c_jp^{n+1}}$ , we can always choose $z_j\in U_j$ such that $\log |\psi ^{S_{p}}_{\eta }|_{h_{p}}(z_{j})> -\delta p$ . Then combining estimate (3.81) with estimate (3.82), we infer that

(3.83) $$ \begin{align} \mathbb{P}\left( -\int_{B(2,3)}\log{|\psi^{S_{p}}_{\eta}|_{h_{p}}}\,\mathrm{dV}\geq K_{4}\delta p\right)\leq e^{-C_{U,\delta}p^{n+1}},\;\forall\;p\gg 0. \end{align} $$

Noting that $\log ^{-}=-\log +\log ^{+}$ and that a finite set of annuli of the form $B(2,3)$ covers U, we can infer estimate (3.77) from estimates (3.76) and (3.83). This completes our proof.

Remark 3.17. Given the information we have for the law of large numbers and large deviations, it is natural to ask about a version of the central limit theorem. Sodin–Tsirelson [Reference Sodin and Tsirelson72, Main Theorem] introduced and proved the asymptotic normality of functionals of the zeros of certain random holomorphic functions on $\mathbb {C}$ or $\mathbb {D}$ . This was extended by Shiffman–Zelditch [Reference Shiffman and Zelditch69, Theorem 1.2] for random holomorphic sections of line bundles on a compact Kähler manifold and for general random polynomials on $\mathbb {C}^{n}$ by Bayraktar [Reference Bayraktar4].

One key ingredient in their approach is the normalized Bergman kernel defined in formula (3.65), viewed as the covariance function of a normalized Gaussian process on $\mathbb {C}$ or X, as constructed in the proof of Proposition 3.15. Using the estimates given in Theorem 3.14 and the seminal result by Sodin and Tsirelson [Reference Sodin and Tsirelson72, Theorem 2.2], one could obtain an extension of [Reference Shiffman and Zelditch69, Theorem 1.2], [Reference Sodin and Tsirelson72, Main Theorem], in the general setting considered here. This direction will be pursued elsewhere.

3.4 Remark on the lower bound for the hole probabilities

To obtain a lower bound of matching order $e^{-cp^{n+1}}$ for the hole probability in estimate (1.14) is generally more complicated. In the case of scaled Bargmann–Fock spaces (cf. Example 3.6), we can give a lower bound and and we sketch its proof below.

Recall that for $p\in \mathbb {N}$ , the family $\{S^{p}_{\alpha }\}_{\alpha \in \mathbb {N}^{n}}$ denotes an orthonormal basis of $H^{0}_{(2)}(\mathbb {C}^{n},L^{p})$ . For $K>0$ , define the index set

(3.84) $$ \begin{align} I(K)=\big \{\alpha=(\alpha_{1},\ldots,\alpha_{n})\in\mathbb{N}^{n}\;:\; 0\leq \alpha_{j}\leq K, j=1,\ldots, n\big\}, \end{align} $$

set $I^{\ast }(K)=I(K)\backslash \{(0,\ldots , 0)\}$ and put

(3.85) $$ \begin{align} q_{p}:=\sharp I(Kp)=(1+\lfloor Kp\rfloor)^{n}=\mathcal{O}(p^{n}). \end{align} $$

For this canonical family of orthonormal basis, we can verify directly the following local concentration condition: For any relatively compact subset $U\subset \mathbb {C}^{n}$ and for any $c>0$ , there exist constants $K=K(U,c)>0, C'=C'(U,c)>0$ such that

(3.86) $$ \begin{align} \sup_{z\in\bar{U}} \sum_{\alpha \notin I(Kp)}|S^{p}_{\alpha}(z)|^{2}_{h_{p}}\leq C' e^{-cp}. \end{align} $$

In fact, to get the above inequality, one can start with the sum over $\alpha _1\in \mathbb {N}$ for the monomial $z_1^{\alpha _1}$ : after applying Stirling’s formula for $\lceil Kp\rceil !$ , $\sum _{\alpha _1> Kp}\frac {|z_1|^{2\alpha _1}}{\alpha _1!}=\mathcal {O}( e^{-c'p})$ , when $z_1$ lies in a given compact subset.

Let $\psi ^{S_{p}}_{\eta }$ be the random holomorphic section (actually, function) on $\mathbb {C}^{n}$ constructed in Example 3.6.

Lemma 3.18. For any relatively compact open subset $U\subset \mathbb {C}^{n}$ , there exists a constant $C^{\prime }_{U}>0$ such that for $p\gg 1$ ,

(3.87) $$ \begin{align} \mathbb{P}\left(\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})\cap U=\varnothing\right)\geq e^{-C^{\prime}_{U}p^{n+1}}. \end{align} $$

Proof. For $U=\varnothing $ the statement is trivial, so assume U nonempty is as in the assumptions. Fix a relatively compact open neighborhood $U'$ of $\overline {U}$ and define the strictly positive quantity

(3.88) $$ \begin{align} \widetilde{M}:=\min_{z\in \overline{U'}} e^{-\frac{|z|^{2}}{2}} \in (0, 1). \end{align} $$

Let the constants K and $C'$ be the ones in estimate (3.86) for the constant $c=-2\log {\widetilde {M}}>0$ and for $U'$ . For $p\in \mathbb {N}$ , write $S^{p}_{0}\equiv p^{n/2}$ for the unit constant section in $H^{0}_{(2)}(\mathbb {C}^{n},L^{p})$ corresponding to $\alpha =(0,\ldots ,0)\in \mathbb {N}^{n}$ . Then

(3.89) $$ \begin{align} \min_{z\in \overline{U'}} |S^{p}_{0}(z)|_{h_{p}}=p^{n/2}\widetilde{M}^{p}. \end{align} $$

Defining the random holomorphic sections

(3.90) $$ \begin{align} \begin{aligned} \psi^{S_{p}}_{\eta,\mathrm{I}}(z):=\sum_{ \alpha \in I^{*}(Kp)}\eta_{\alpha} S^{p}_{\alpha}(z)\,, \quad \psi^{S_{p}}_{\eta,\mathrm{II}}(z):=\sum_{\alpha \notin I(Kp)}\eta_{\alpha} S^{p}_{\alpha}(z), \end{aligned} \end{align} $$

we can decompose

(3.91) $$ \begin{align} \psi^{S_{p}}_{\eta}=\eta_{0} S^{p}_{0} + \psi^{S_{p}}_{\eta,\mathrm{I}} + \psi^{S_{p}}_{\eta,\mathrm{II}}. \end{align} $$

Note that the three random sections on the right-hand side of equality (3.91) are independent from each other. In the rest of the proof, we will regard the above sections as being holomorphic functions on $\mathbb {C}^{n}$ , and let $|\cdot |$ denote the standard modulus on $\mathbb {C}$ (instead of considering the norm $|\cdot |_{h_{p}}$ on line bundle). Applying estimate (2.9) to the function $\psi ^{S_{p}}_{\eta ,\mathrm {II}}$ and using the estimate (3.86), we arrive at the upper bound

(3.92) $$ \begin{align} \begin{aligned} \mathbb{E}\Big[\sup_{z\in U}\Big|\psi^{S_{p}}_{\eta,\mathrm{II}}(z)\Big|^{2}\Big]&\leq C_{U'}\sigma^{2}\int_{U'}\sum_{\alpha \notin I(Kp)}|S^{p}_{\alpha}(z)|^{2}\mathrm{dV}(z)\\ &\leq \widetilde{C}_{U'}\operatorname{\mathrm{Vol}}(U')\sigma^{2} \widetilde{M}^{-2p} e^{-cp}\\ &= \widetilde{C}_{U'}\operatorname{\mathrm{Vol}}(U')\sigma^{2} =:\widetilde{C}', \end{aligned} \end{align} $$

where the last equality follows from our choice $c=-2\log {\widetilde {M}}$ . For any $\lambda>0,$ as a consequence of Chebyshev’s inequality in combination with estimate (3.92), we have

(3.93) $$ \begin{align} \mathbb{P}\Big( \sup_{z\in U}\Big|\psi^{S_{p}}_{\eta,\mathrm{II}}(z)\Big| < \lambda\Big)\geq 1-\frac{\widetilde{C}'}{\lambda^{2}}. \end{align} $$

We define the good event

(3.94) $$ \begin{align} \Omega_{p}=\left\{|\eta_{0}|\geq 1\,;\, |\eta_{\alpha}|\leq \frac{1}{3\sqrt{q_{p}-1}}\widetilde{M}^{p}, \alpha \in I^{\ast}(Kp)\,;\, \sup_{z\in U}\left|\psi^{S_{p}}_{\eta,\mathrm{II}}(z)\right| < \frac{1}{3}p^{n/2} \right\}. \end{align} $$

For all sufficiently large $p\in \mathbb {N}$ , we have

(3.95) $$ \begin{align} \begin{aligned} \mathbb{P}(\Omega_{p})&=\mathbb{P}(|\eta_{0}|\geq 1) \cdot \mathbb{P}\left(\sup_{z\in U}\left|\psi^{S_{p}}_{\eta,\mathrm{II}}(z)\right| < \frac{1}{3}p^{n/2} \right)\\ &\qquad\cdot\mathbb{P}\Bigg(|\eta_{\alpha}|\leq \frac{1}{3\sqrt{q_{p}-1}}\widetilde{M}^{p}, \alpha \in I^{\ast}(Kp)\Bigg)\\ &\geq e^{-1}\Big(1-\frac{9\widetilde{C}'}{ p^{n}}\Big)\cdot \Pi_{\alpha \in I^{\ast}(Kp)}\left(\frac{1}{18(q_{p}-1)}\widetilde{M}^{2p}\right). \end{aligned} \end{align} $$

Then by formula (3.85), there exists $C^{\prime }_{U}>0$ such that for $p\gg 1$ ,

(3.96) $$ \begin{align} \mathbb{P}(\Omega_{p})\geq e^{-C^{\prime}_{U}p^{n+1}}. \end{align} $$

Our lemma then follows once we show the inclusion

(3.97) $$ \begin{align} \Omega_{p}\subset \left\{\operatorname{\mathrm{Div}}(\psi^{S_{p}}_{\eta})\cap U=\varnothing\right\}. \end{align} $$

Indeed, if $|\eta _{\alpha }|\leq \frac {1}{3\sqrt {q_{p}-1}}\widetilde {M}^{p}$ , $\alpha \in I^{\ast }(Kp)$ , then for $z\in U$ ,

(3.98) $$ \begin{align} \begin{aligned} \big|\psi^{S_{p}}_{\eta,\mathrm{I}}(z)\big|^{2}&\leq \Big(\sum_{\alpha\in I^{\ast}(pK)}|\eta_{\alpha}|^{2}\Big)\Big(\sum_{\alpha\in I^{\ast}(pK)}|S^{p}_{\alpha}(z)|^{2}\Big)\\ &\leq \frac{1}{\widetilde{M}^{2p}}\Big(\sum_{\alpha\in I^{\ast}(pK)}|\eta_{\alpha}|^{2}\Big)P_{p}(z,z)\\ &\leq \frac{1}{9}p^{n}. \end{aligned} \end{align} $$

As a consequence, on $\Omega _{p}$ and for $z\in U$ , we get

(3.99) $$ \begin{align} \begin{aligned} \big|\psi^{S_{p}}_{\eta,\mathrm{I}}(z) + \psi^{S_{p}}_{\eta,\mathrm{II}}(z)\big|&\leq \big|\psi^{S_{p}}_{\eta,\mathrm{I}}(z)\big|+ \big|\psi^{S_{p}}_{\eta,\mathrm{II}}(z)\big| \leq \frac{1}{3}p^{n/2}+\frac{1}{3}p^{n/2}\\ &< p^{n/2}\leq \left|\eta_{0}S^{p}_{0}(z)\right|. \end{aligned} \end{align} $$

The above strict inequality implies that the inclusion (3.97) holds. This finishes the proof of the lemma.

We now shortly explain how applying our results to the special case of the Bargmann–Fock space recovers the results by Sodin–Tsirelson (for $\mathbb {C}$ , [Reference Sodin and Tsirelson73, Theorem 1]) and Zrebiec (for $\mathbb {C}^{n}$ , [Reference Zrebiec82, Theorem 1.2]) about the hole probability. They proved that there exist constants $c_{1}\geq c_{2}>0$ such that for $r>0$ large,

(3.100) $$ \begin{align} \exp(-c_{1}r^{2n+2})\leq \mathbb{P}\big (\psi^{S}_{\eta}(z)\neq 0,\; \mathrm{for\;all\;}z\in \mathbb{B}(0,r)\big)\leq \exp(-c_{2}r^{2n+2}), \end{align} $$

where $\mathbb {B}(0,r)=\{z\in \mathbb {C}^{n}\;:\; |z|<r\}$ . Let us now fix $r_{0}>0$ . Then by estimates (1.14) and (3.87), we get

(3.101) $$ \begin{align} \exp(-c \sqrt{p}^{2n+2})\leq \mathbb{P}\big( \psi^{S_{p}}_{\eta}(z)\neq 0,\; \mathrm{for\;all\;}z\in \mathbb{B}(0,r_{0})\big)\leq \exp(-c'\sqrt{p}^{2n+2}). \end{align} $$

By using formula (3.17), the inequality (3.101) is equivalent to

(3.102) $$ \begin{align} \exp(-c \sqrt{p}^{2n+2})\leq \mathbb{P}\big( \psi^{S}_{\eta}(z)\neq 0,\; \mathrm{for\;all\;}z\in \mathbb{B}(0,\sqrt{p}r_{0})\big)\leq \exp(-c'\sqrt{p}^{2n+2}). \end{align} $$

Therefore, we recover the estimates (3.100) by approximating a sufficiently large $r>0$ by $\sqrt {p}r_{0}$ .

Remark 3.19. In the context of a general complete Kähler manifold X, a question related to estimate (3.86) is to find for any relatively compact open subset $U\subset X$ a sequence of orthonormal basis $\{\widetilde {S}^{p}_{j}\}_{j=1}^{d_{p}}$ of $H^{0}_{(2)}(X,L^{p})$ , $p\in \mathbb {N}$ , such that

(3.103) $$ \begin{align} \sup_{x\in\overline{U}} \sum_{j>K' p^{n}}|\widetilde{S}^{p}_{j}(x)|^{2}_{h_{p}}\leq C e^{-cp}, \end{align} $$

where C, $K'$ , c are certain positive constants independent of p, and the sum in the left-hand side is taken to be $0$ if $d_{p}=\dim H^{0}_{(2)}(X,L^{p})\leq K'p^{n}$ . This question is trivial if $d_{p}=\mathcal {O}(p^{n})$ for $p\gg 0$ .

The existence of such a sequence of basis suggests that, on a relatively compact subset, the Bergman projections or Bergman kernels can be approximated by the orthogonal projections or their kernels of a sequence of finite-dimensional subspaces of $H^{0}_{(2)}(X,L^{p})$ . Moreover, one may expect a connection between the number (or dimension of the aforementioned subspace) $K'p^{n}$ and the integration of the dimension density on U,

(3.104) $$ \begin{align} \int_{U}P_{p}(x,x)\mathrm{dV}(x). \end{align} $$

4 Random $\mathcal {L}^{2}$ -holomorphic sections and Toeplitz operators

In this section, we always use the setting in Section 2: $(X,J,\Theta )$ is a connected complex Hermitian manifold (without boundary), and $(L,h_{L})$ is a Hermitian line bundle on X. We do not, however, assume any completeness for $\Theta $ or positivity for $(L,h_{L})$ .

The goal of this section is to introduce a method of ‘canonically randomizing’ the $\mathcal {L}^{2}$ -holomorphic sections of L on X, in particular when $d=\dim H^{0}_{(2)}(X, L)=\infty $ .

As mentioned in the Introduction, this is achieved by the abstract Wiener space construction from probability theory with an approach via Toeplitz operators from the theory of geometric quantization. This induces a Gaussian probability measure on the space of $\mathcal {L}^{2}$ -holomorphic sections.

4.1 Abstract Wiener spaces

To define a Gaussian probability measure on an infinite dimensional Hilbert space, we employ the construction of abstract Wiener spaces introduced by Gross [Reference Gross36]. We also refer to the article of Sheffield [Reference Sheffield67] for further motivation and developments on this topic. For a (complex) vector space $\mathcal {H}$ , a Hermitian norm is a norm on $\mathcal {H}$ which is induced by a Hermitian inner product on it.

Definition 4.1. Let $(\mathcal {H},\|\cdot \|_{\mathcal {H}})$ be a separable Hilbert space of infinite dimension. A Hermitian norm $\|\cdot \|$ is called measurable if for all $\epsilon>0$ , there exists a finite-dimensional subspace $F_{\epsilon }\subset \mathcal {H}$ such that for $F\subset \mathcal {H}$ a subspace of finite dimension with $F\perp F_{\epsilon }$ , one has

(4.1) $$ \begin{align} \mu_{F,\; \|\cdot\|_{\mathcal{H}}} \big(\{x\in F\;:\; \|x\|\geq \epsilon \}\big) < \epsilon, \end{align} $$

where $\mu _{F,\; \|\cdot \|_{\mathcal {H}}}$ denotes the standard Gaussian measure on F with respect to the Hermitian metric associated with $\|\cdot \|_{\mathcal {H}}$ .

Proposition 4.2 (cf. [Reference Gross36],[Reference Kuo49, Chapter I: Theorem 4.3]).

Let $\mathcal {H}$ be a separable Hilbert space with norm $\|\cdot \|_{\mathcal {H}}$ , and $\|\cdot \|$ be a continuous (with respect to $\|\cdot \|_{\mathcal {H}}$ ) Hermitian norm on $\mathcal {H}$ . Then the following two conditions are equivalent:

  1. 1. $\|\cdot \|$ is measurable.

  2. 2. There exists a one-to-one Hilbert–Schmidt operator T on $\mathcal {H}$ such that $\|x\|=\|Tx\|_{\mathcal {H}}$ for $x\in \mathcal {H}$ .

Given a measurable Hermitian norm $\|\cdot \|$ on $\mathcal {H}$ , let $\mathcal {B}$ be the completion of $\mathcal {H}$ with respect to $\|\cdot \|$ . Then $(\mathcal {B},\|\cdot \|)$ is a separable Hilbert space containing $\mathcal {H}$ as a dense subspace. Let $\mathcal {B}^{*}$ be the topological dual space of $\mathcal {B}$ . If $\alpha \in \mathcal {B}^{*}$ , then $\alpha |_{\mathcal {H}}$ is a continuous linear functional on $(\mathcal {H},\|\cdot \|_{\mathcal {H}})$ . If $\alpha $ vanishes identically on $\mathcal {H}$ , then it vanishes on $\mathcal {B}$ . This way, we can regard $\mathcal {B}^{*}$ as a (dense) subspace of $\mathcal {H}^{*}$ , where $\mathcal {H}^{*}$ can be identified with $\mathcal {H}$ via the Hilbert metric associated with $\|\cdot \|_{\mathcal {H}}$ .

By a slight abuse of notation, we denote by $\mathcal {S}$ the Borel $\sigma $ -algebra of $\mathcal {B}$ . Then each $\alpha \in \mathcal {B}^{*}$ is a Borel-measurable function from $\mathcal {B}$ to $\mathbb {C}$ . For $F\subset \mathcal {B}^{*}\subset \mathcal {H}$ an arbitrary finite-dimensional subspace we introduce the notation

(4.2) $$ \begin{align} \begin{aligned} \phi_{F}: \mathcal{B}&\rightarrow F,\; \phi_{F}(b)=\sum_{j=1}^{\dim_{\mathbb{C}} F} (b,v_{j})v_{j}, \end{aligned} \end{align} $$

where $\{v_{j}\}$ is an orthonormal basis of $(F, \|\cdot \|_{\mathcal {H}})$ . Gross [Reference Gross36] proved the following result.

Theorem 4.3. Fix a measurable norm $\|\cdot \|$ on $\mathcal {H}$ as above. There exists a unique probability measure $\mathcal {P}$ on $(\mathcal {B},\mathcal {S})$ such that for $F\subset \mathcal {B}^{*}$ any finite-dimensional subspace,

(4.3) $$ \begin{align} \mathcal{P}(\phi^{-1}_{F}(U))=\mu_{F,\;\|\cdot\|_{\mathcal{H}}}(U), \end{align} $$

for all Borel subset U of F. The triple $(\mathcal {B},\mathcal {S},\mathcal {P})$ is called an abstract Wiener space.

If $\alpha \in \mathcal {B}^{*}$ , then as a function on $\mathcal {B}$ , it is an element of $\mathcal {L}^{2}(\mathcal {B},\mathcal {S},\mathcal {P})$ . We denote this map by

(4.4) $$ \begin{align} \Phi_{0}:\mathcal{B}^{*}\rightarrow \mathcal{L}^{2}(\mathcal{B},\mathcal{S},\mathcal{P}). \end{align} $$

Moreover, for $\alpha \in \mathcal {B}^{*}$ , $\Phi _{0}(\alpha )$ is a Gaussian random variable with zero mean and variance $\|\alpha \|^{2}_{\mathcal {H}}$ . The map $\Phi _{0}$ extends to a continuous linear map

(4.5) $$ \begin{align} \Phi:\mathcal{H}^{*}\simeq \mathcal{H}\rightarrow \mathcal{L}^{2}(\mathcal{B},\mathcal{S},\mathcal{P}), \end{align} $$

where for $y\in \mathcal {H}$ , $\Phi (y)$ is a Gaussian random variable with zero mean and variance $\|y\|^{2}_{\mathcal {H}}$ .

Remark that the above construction is trivial if $\mathcal {H}$ is finite dimensional; indeed, in this case the Hilbert space $\mathcal {B}$ is reduced to $\mathcal {H}$ itself. The probability measure constructed in Theorem 4.3 is the standard Gaussian probability measure on $\mathcal {H}$ with respect to the norm $\|\cdot \|_{\mathcal {H}}$ .

4.2 Toeplitz operators on $H^{0}_{(2)}(X,L)$

Recall that P denotes the orthogonal projection from $\mathcal {L}^{2}(X,L)$ onto $H^{0}_{(2)}(X,L)$ , and $P(x,y)$ , $x,y\in X$ , denotes the corresponding Bergman kernel. W.l.o.g. we may and do always assume that $d=\dim H^{0}_{(2)}(X,L)\geq 1$ in the following.

Definition 4.4. For a bounded function $f\in \mathscr {C}^{\infty }(X,\mathbb {C})$ , set

(4.6) $$ \begin{align} T_{f}: H^{0}_{(2)}(X, L) \rightarrow H^{0}_{(2)}(X,L),\quad T_{f}:=P f, \end{align} $$

where the action of f is the pointwise multiplication by f. The operator $T_{f}$ is called Toeplitz operator associated with f. Equivalently, we consider $T_{f}:\mathcal {L}^{2}(X, L) \rightarrow \mathcal {L}^{2}(X,L)$ , $T_{f}:=P fP$ .

The integral kernel of $T_{f}$ is given by

(4.7) $$ \begin{align} T_{f}(x,x')=\int_{X} P(x,x^{\prime\prime})f(x^{\prime\prime})P(x^{\prime\prime},x') \mathrm{dV}(x^{\prime\prime}). \end{align} $$

Note also that the Hilbert adjoint of $T_{f}$ is $T_{\overline {f}}$ . We introduce a class of bounded smooth functions on X whose associated Toeplitz operators are Hilbert–Schmidt.

Definition 4.5. Let $\mathcal {Q}(X,L\,;\,\mathbb {C})$ be the vector space of bounded smooth complex functions f on X such that

(4.8) $$ \begin{align} \int_{X}|f(x)|P(x,x)\mathrm{dV}(x)<\infty, \end{align} $$

where P is the Bergman kernel of L.

Example 4.6. (1) It is clear that $\mathscr {C}^{\infty }_{c}(X,\mathbb {C})$ is a subspace of $\mathcal {Q}(X,L\,;\,\mathbb {C})$ . In particular, if X is compact, then

(4.9) $$ \begin{align} \mathcal{Q}(X,L\,;\,\mathbb{C})=\mathscr{C}^{\infty}(X,\mathbb{C}). \end{align} $$

(2) Let $\mathscr {C}^{\infty }_{b}(\mathbb {C}^{n},\mathbb {C})$ denote the set of bounded smooth functions on $\mathbb {C}^{n}$ . In the case of the Bargmann–Fock space (see Example 2.14), we have

(4.10) $$ \begin{align} \mathcal{Q}(\mathbb{C}^{n},L\,;\,\mathbb{C})=\mathscr{C}^{\infty}_{b}(\mathbb{C}^{n},\mathbb{C})\cap \mathcal{L}^{1}(\mathbb{C}^{n}, \mathrm{dV}). \end{align} $$

(3) In general, with the assumptions as in Section 3, if we assume further that $(X,J,\Theta )$ , $(L,h_{L})$ have bounded geometry, by [Reference Ma and Marinescu60, Theorem 6], there exist $c>0$ , $C>0$ and $p_{0}\in \mathbb {N}^{*}$ such that for $p\geq p_{0}$ ,

(4.11) $$ \begin{align} cp^{n}\leq \inf_{x\in X} P_{p}(x,x)\leq \sup_{x\in X} P_{p}(x,x)\leq Cp^{n}, \end{align} $$

that is, the Bergman kernel function $P_{p}(x,x)$ is bounded from above and away from zero on X. As a consequence, we get that for $p\geq p_{0}$ ,

(4.12) $$ \begin{align} \mathcal{Q}(X,L^{p}\,;\,\mathbb{C})=\mathscr{C}^{\infty}_{b}(X,\mathbb{C})\cap \mathcal{L}^{1}(X, \mathrm{dV}). \end{align} $$

Proposition 4.7. For $f\in \mathcal {Q}(X,L\,;\,\mathbb {C})$ , the operator $T_{f}$ on $H^{0}_{(2)}(X,L)$ has smooth Schwartz kernel and is Hilbert–Schmidt.

Proof. If $d=\dim H^{0}_{(2)}(X,L)<\infty $ , then the statement is trivial. Hence, we assume $d=\infty $ w.l.o.g. in the sequel. Let $\{S_{j}\}_{j=1}^{\infty }$ be a complete Hilbert basis of $H^{0}_{(2)}(X,L)$ . Note that by [Reference Atiyah2, Proposition (2.4)], for any compact set $K\subset X$ , the series

(4.13) $$ \begin{align} \sum_{j=1}^{\infty}|S_{j}(x)|_{h_{L}}^{2} \end{align} $$

converges uniformly for $x\in K$ . As a consequence, for $K_{1}, K_{2} \subset X$ compact, the series

(4.14) $$ \begin{align} \sum_{j=1}^{\infty}S_{j}(x)\otimes (S_{j}(y))^{\ast} \end{align} $$

converges absolutely and uniformly for $x\in K_{1}$ and $y\in K_{2}$ .

As follows from the properties of holomorphic functions, if we replace $S_{j}(x)$ , $(S_{j}(y))^{\ast }$ by their respective covariant derivatives, then the series (4.13) and (4.14) are still absolutely convergent on any given compact subsets.

By Definition 4.5, for $j\in \mathbb {N}^{\ast }$ , the function $X\ni x\mapsto f(x)|S_{j}(x)|^{2}_{h_{L}}$ is integrable on X with respect to $\mathrm {dV}$ . Furthermore, for $x'\in X$ , $i,j\in \mathbb {N}^{\ast }$ , we have

(4.15) $$ \begin{align} f(x')(S_{i}(x'))^{\ast}S_{j}(x')=f(x')h_{L,x'}\left(S_{j}(x'),S_{i}(x')\right), \end{align} $$

and

(4.16) $$ \begin{align} \int_{X}\left|f(x')\left(S_{i}(x')\right)^{\ast}S_{j}(x')\right| \mathrm{dV}(x')\leq \big\|\sqrt{|f|}S_{i}\big\|_{\mathcal{L}^{2}(X,L)} \cdot\big\|\sqrt{|f|}S_{j}\big\|_{\mathcal{L}^{2}(X,L)}. \end{align} $$

Now, we fix two compact subsets $K_{1}$ , $K_{2}\subset X$ . For $x\in K_{1}$ , $y\in K_{2}$ and $i,j\in \mathbb {N}^{\ast }$ we have

(4.17) $$ \begin{align} \big|S_{i}(x)\otimes (S_{i}(x'))^{\ast}f(x')S_{j}(x')\otimes (S_{j}(y))^{\ast}\big|\leq |S_{i}(x)|_{h_{L}}\cdot \big|f(x')(S_{i}(x'))^{\ast}S_{j}(x')\big|\cdot |S_{j}(y)|_{h_{L}}, \end{align} $$

where the norm in the left-hand side is given by $h_{L, x}\otimes h^{\ast }_{L, y}$ . By inequality (4.16) this entails

(4.18) $$ \begin{align} \begin{aligned} &\int_{X}\big|S_{i}(x)\otimes (S_{i}(x'))^{\ast}f(x')S_{j}(x')\otimes (S_{j}(y))^{\ast}\big|\mathrm{dV}(x')\\&\quad \leq |S_{i}(x)|_{h_{L}}\big\Vert\sqrt{|f|}S_{i}\big\Vert{}_{\mathcal{L}^{2}(X,L)}\cdot\big\Vert \sqrt{|f|}S_{j}\big\Vert{}_{\mathcal{L}^{2}(X,L)} |S_{j}(y)|_{h_L}. \end{aligned} \end{align} $$

Putting things together, we arrive at

(4.19) $$ \begin{align} \begin{aligned} \sum_{i=1}^{\infty}|S_{i}(x)|_{h_{L}}\|\sqrt{|f|}S_{i}\|_{\mathcal{L}^{2}(X,L)} &\leq \Big(\sum_{i=1}^{\infty}|S_{i}(x)|^{2}_{h_{L}}\Big)^{1/2}\Big(\sum_{i=1}^{\infty}\|\sqrt{|f|}S_{i}\|^{2}_{\mathcal{L}^{2}(X,L)} \Big)^{1/2}\\ &=\Big(\sum_{i=1}^{\infty}|S_{i}(x)|^{2}_{h_{L}}\Big)^{1/2}\Big(\int_{X}|f(x')|P(x',x')\mathrm{dV}(x')\Big)^{1/2}\\ &<\infty, \end{aligned} \end{align} $$

and the above estimates still hold if we replace $S_{i}(x)$ by its covariant derivatives at x.

Recalling the Schwartz kernel of $T_{f}$ from (4.7), the above calculations show that $T_{f}(x,y)$ is a smooth section on $X\times X$ . For proving that $T_{f}$ is Hilbert–Schmidt, it only remains to show that

(4.20) $$ \begin{align} \sum_{i,j}\left|\langle T_{f}S_{i},S_{j}\rangle_{\mathcal{L}^{2}(X,L)}\right|{}^{2}<\infty. \end{align} $$

Indeed, by inequality (4.16), we have

(4.21) $$ \begin{align} \left|\langle T_{f}S_{i},S_{j}\rangle_{\mathcal{L}^{2}(X,L)}\right|{}^{2}\leq \left\|\sqrt{|f|}S_{i}\right\|_{\mathcal{L}^{2}(X,L)}^{2}\cdot\left\|\sqrt{|f|}S_{j}\right\|_{\mathcal{L}^{2}(X,L)}^{2}. \end{align} $$

Then

(4.22) $$ \begin{align} \begin{aligned} \sum_{i,j}\left|\langle T_{f}S_{i},S_{j}\rangle_{\mathcal{L}^{2}(X,L)}\right|{}^{2} &\leq\sum_{i,j=1}^{\infty} \left\|\sqrt{|f|}S_{i}\right\|_{\mathcal{L}^{2}(X,L)}^{2}\cdot\left\|\sqrt{|f|}S_{j}\right\|_{\mathcal{L}^{2}(X,L)}^{2}\\ &= \Big(\int_{X} |f(x)|P(x,x)\mathrm{dV}(x)\Big)^{2}<\infty. \end{aligned} \end{align} $$

This completes our proof.

Corollary 4.8. If $f\in \mathscr {C}^{\infty }(X,\mathbb {C})$ is with compact support, then $T_{f}$ is a Hilbert–Schmidt operator on $H^{0}_{(2)}(X, L)$ . Moreover, $T_{f}$ is trace class, and

(4.23) $$ \begin{align} \operatorname{\mathrm{Tr}}[T_{f}]=\int_{X}f(x)P(x,x)\mathrm{dV}(x). \end{align} $$

4.3 Random $\mathcal {L}^{2}$ -holomorphic sections of L

Let $\mathcal {Q}(X,L\,;\,\mathbb {R}_{\geq 0})$ denote the subspace of $\mathcal {Q}(X,L\,;\,\mathbb {C})$ consisting of the functions valued in $\mathbb {R}_{\geq 0}$ . For $f\in \mathcal {Q}(X,L\,;\,\mathbb {R}_{\geq 0})$ , $T_{f}$ is a nonnegative self-adjoint Hilbert–Schmidt (hence compact) operator on $H^{0}_{(2)}(X,L)$ .

Lemma 4.9. For $0\neq f\in \mathcal {Q}(X,L\,;\,\mathbb {R}_{\geq 0})$ , the operator $T_{f}: H^{0}_{(2)}(X,L) \rightarrow H^{0}_{(2)}(X,L)$ is injective.

Proof. Since $f\neq 0$ , there exists an open subset U of X on which f is strictly positive. If $s\in H^{0}_{(2)}(X,L)$ is such that $T_{f}s=0$ , then $0=\langle T_{f}s, s\rangle =\int _{X}f(x)|s(x)|^{2}_{h_L}\mathrm {dV}(x)$ , and hence $s|_{U}=0$ . Since U is open and s is holomorphic on X, we get $s=0$ . This proves the lemma.

Fix $f\in \mathcal {Q}(X,L\,;\,\mathbb {R}_{\geq 0})$ , $f\neq 0$ . If $d<\infty $ , then the above $T_{f}$ is actually an isomorphism on the vector space $H^{0}_{(2)}(X,L)$ . Now, we focus on the case of $d=\infty $ . Since $T_{f}$ is compact and injective, it cannot be surjective. Hence, it does not admit a bounded inverse. Moreover, for any $\lambda \in \mathbb {C}$ , $\lambda \neq 0$ , the operator $T_{f}-\lambda $ is Fredholm with closed range and Fredholm index $0$ . Set $D(T^{-1}_{f})=\mathrm {Range}(T_{f}:H^{0}_{(2)}(X,L) \rightarrow H^{0}_{(2)}(X,L))\subset H^{0}_{(2)}(X,L)$ , which is a dense subspace. The inverse of $T_{f}$ is defined as

(4.24) $$ \begin{align} T_{f}^{-1}:D(T^{-1}_{f})\subset H^{0}_{(2)}(X,L)\rightarrow H^{0}_{(2)}(X,L). \end{align} $$

Let $\sigma (T_{f})\subset \mathbb {R}_{\geq 0}$ denote the spectrum of $T_{f}$ , which is a countable set consisting of two parts: the point spectrum $\sigma _{\mathrm {p}}(T_{f})\subset \mathbb {R}_{> 0}$ (eigenvalues) and the residual spectrum $\sigma _{\mathrm {res}}(T_{f})=\{0\}$ . In this case, the point spectrum of $T_{f}$ is a decreasing sequence of strictly positive eigenvalues of finite multiplicity,

(4.25) $$ \begin{align} \lambda_{1}\geq \lambda_{2}\geq \ldots\geq \lambda_{m}\geq \ldots\rightarrow 0. \end{align} $$

Since any separable (complex) Hilbert space is isometric to the Hilbert space $\ell ^{2}(\mathbb {C})$ by choosing an orthonormal basis, we can choose an orthonormal basis $\{S_{j}\}^{\infty }_{j=1}$ of $H^{0}_{(2)}(X,L)$ with respect to the $\mathcal {L}^{2}$ -metric such that

(4.26) $$ \begin{align} T_{f}S_{j}=\lambda_{j}S_{j}. \end{align} $$

If $S\in H^{0}_{(2)}(X,L)$ , we can write uniquely

(4.27) $$ \begin{align} S=\sum_{j\geq 1}a_{j}S_{j},\quad a_{j}\in\mathbb{C}. \end{align} $$

Then $(a_{j})_{j}\in \ell ^{2}(\mathbb {C})$ , yielding the identification between $H^{0}_{(2)}(X,L)$ and $\ell ^{2}(\mathbb {C})$ .

Since $T_{f}$ is one-to-one and Hilbert–Schmidt, by Proposition 4.2, $\|\cdot \|_{f}:=\|T_{f}\cdot \|$ defines a Hermitian measurable norm on $H^{0}_{(2)}(X,L)$ . We denote by $\mathcal {B}_{f}(X,L)$ the completion of $H^{0}_{(2)}(X,L)$ with respect to $\|\cdot \|_{f}$ and set

(4.28) $$ \begin{align} \ell^{2}_{f}(\mathbb{C})=\Big\{(a_{j}\in\mathbb{C})_{j\geq 1}\;:\;\sum_{j\geq 1}\lambda_{i}^{2}|a_{j}|^{2}<\infty \Big\}. \end{align} $$

It is clearly a separable Hilbert space, and using the basis as in (4.26), we have

(4.29) $$ \begin{align} \mathcal{B}_{f}(X,L)\simeq \ell^{2}_{f}(\mathbb{C}). \end{align} $$

Proposition 4.10. Assume $d=\infty \,$ , $0\neq f\in \mathcal {Q}(X,L\,;\,\mathbb {R}_{\geq 0})$ . Then the operator $T_{f}$ extends uniquely to an isomorphism of Hilbert spaces

(4.30) $$ \begin{align} \widehat{T}_{f}: \big(\mathcal{B}_{f}(X,L),\|\cdot\|_{f} \big)\rightarrow \big (H^{0}_{(2)}(X,L),\|\cdot\|_{\mathcal{L}^{2}(X,L)} \big). \end{align} $$

Given $0\neq f\in \mathcal {Q}(X,L\,;\,\mathbb {R}_{\geq 0})$ , if $d<\infty $ , we set

(4.31) $$ \begin{align} \big(\mathcal{B}_{f}(X,L),\|\cdot\|_{f} \big)= \big(H^{0}_{(2)}(X,L),\|\cdot\|_{f} \big),\quad \text{and} \quad \widehat{T}_{f}:=T_{f}. \end{align} $$

Then we unify our notation for both cases $d<\infty $ and $d=\infty $ .

Definition 4.11. Denote by $\mathcal {P}_{f}$ the probability measure from Theorem 4.3 with the choice $\mathcal {B} = \mathcal {B}_{f}(X,L).$ Let $\mathbb {P}_{f}$ be the Gaussian probability measure on $H^{0}_{(2)}(X,L)$ given by the pushforward of $\mathcal {P}_{f}$ through the isomorphism (4.30). We call $\mathbb {P}_{f}$ the Gaussian probability measure induced by the spectral decomposition of the Toeplitz operator $T_f$ . This way, we randomize the sections in $H^{0}_{(2)}(X,L)$ .

Remark 4.12. When $0<d<\infty $ , then $\mathcal {B}_{f}(X,L)=H^{0}_{(2)}(X,L)$ , and $\mathcal {P}_{f}=\mathbb {P}_{\mathrm {st}}$ is exactly the standard Gaussian probability measure on $H^{0}_{(2)}(X,L)$ with respect to the $\mathcal {L}^{2}$ inner product. So that $\mathbb {P}_{f}$ is the pushforward of $\mathbb {P}_{\mathrm {st}}$ via the isomorphism $T_{f}$ .

Lemma 4.13. Assume $d\geq 1\,$ , $0\neq f\in \mathcal {Q}(X,L\,;\,\mathbb {R}_{\geq 0})$ . For any nonzero $S\in H^{0}_{(2)}(X,L)$ , the random variable on $(H^{0}_{(2)}(X,L), \mathbb {P}_{f})$ defined as $H^{0}_{(2)}(X,L)\ni s\mapsto \langle s, S\rangle _{\mathcal {L}^{2}(X,L)}\in \mathbb {C}$ is a centered complex Gaussian variable with variance $\|T_{f}S\|^{2}_{\mathcal {L}^{2}(X,L)}$ .

Proof. Note that $T_{f}S$ is nonzero in $H^{0}_{(2)}(X,L)$ , the linear form

(4.32) $$ \begin{align} H^{0}_{(2)}(X,L)\ni s'\mapsto \langle s', T_{f}S\rangle_{\mathcal{L}^{2}(X,L)}\in \mathbb{C} \end{align} $$

extends to a bounded linear form on $(\mathcal {B}_{f}(X,L),\|\cdot \|_{f})$ , hence defines an element in $\mathcal {B}_{f}(X,L)^{\ast }$ , denoted by $\Psi _{S}$ . Then by property (4.3), the random variable $\Psi _{S}(s')$ with $s'$ having the law $\mathcal {P}_{f}$ , is a centered complex Gaussian variable with variance $\|T_{f}S\|^{2}_{\mathcal {L}^{2}(X,L)}$ . Put differently, by construction, for $s'\in \mathcal {B}_{f}(X,L)$ ,

(4.33) $$ \begin{align} \Psi_{S}(s')=\langle \widehat{T}_{f}s',S\rangle_{\mathcal{L}^{2}(X,L)}. \end{align} $$

Thus, as a random variable, it is exactly the same as $\langle s, S\rangle _{\mathcal {L}^{2}(X,L)}$ with s having distribution $\mathbb {P}_{f}$ . This completes the proof.

4.4 Zeros of random $\mathcal {L}^{2}$ -holomorphic sections: proof of Theorem 1.6

We assume $d\geq 1$ , and we fix $0\neq f\in \mathcal {Q}(X,L\,;\,\mathbb {R}_{\geq 0})$ . The operator $T^{2}_{f}:= T_{f}\circ T_{f}$ on $H^{0}_{(2)}(X,L)$ is a positive self-adjoint operator of trace class. Let $T^{2}_{f}(x,y)$ denote the Schwartz kernel of $T^{2}_{f}$ .

Lemma 4.14. The function $X\ni x\mapsto \log T^{2}_{f}(x,x)$ is locally integrable on X so that the $(1,1)$ -current $\partial \bar {\partial }\log T^{2}_{f}(x,x)$ is well defined on X.

Proof. Let $\{S_{j}\}_{j=1}^{d}$ be the orthonormal basis of $H^{0}_{(2)}(X,L)$ satisfying the equalities (4.26). Then

(4.34) $$ \begin{align} T^{2}_{f}(x,x)=\sum_{j=1}^{d}\lambda_{j}^{2}|S_{j}(x)|^{2}_{h_{L}}\,, \quad\text{for } x\in X. \end{align} $$

If $d=\infty $ , the above sum is uniformly convergent on any compact subset of X. Similar to the proof of Lemma 2.6, we get that the function $\log T^{2}_{f}(x,x)$ is a quasi-plurisubharmonic function on X, hence locally integrable. This completes our proof.

We can now prove Theorem 1.6 for the zeros of the random $\mathcal {L}^{2}$ -holomorphic sections constructed in last subsection.

Proof of Theorem 1.6.

By (4.34), $T^{2}_{f}(x,x)$ vanishes exactly on $\text {Bl}(X,L)$ . Let $\{S_{j}\}_{j=1}^{d}$ be the orthonormal basis of $H^{0}_{(2)}(X,L)$ satisfying the equalities (4.26). By Lemma 4.13, the complex random variables

(4.35) $$ \begin{align} \eta_{j}:=\frac{1}{\lambda_{j}}\langle \boldsymbol{s}, S_{j}\rangle_{\mathcal{L}^{2}(X,L)},\;j=1,\,2,\,\ldots \end{align} $$

form an i.i.d. sequence of standard centered complex Gaussian variable. As a consequence, we get that for $x\in X$ ,

(4.36) $$ \begin{align} \boldsymbol{s}(x)=\sum_{j}\eta_{j}\lambda_{j}S_{j}(x). \end{align} $$

Recalling the definition of Fubini–Study currents from formula (1.17), we can proceed as in the proof of Theorem 1.1, replacing $P(x,x)$ by $T^{2}_{f}(x,x)$ given in formula (4.34), and we conclude formula (1.18).

Remark 4.15. In the above proof, we observe that the random $\mathcal {L}^{2}$ -holomorphic section $\boldsymbol {s}$ on $(H^{0}_{(2)}(X,L),\mathbb {P}_{f})$ is equivalent to the construction given in formula (4.36), as we explained in the Introduction (cf. Equation (1.2)). More precisely, denote by $S=\{S_{j}\}_{j=1}^{d}$ the orthonormal basis of $H^{0}_{(2)}(X,L)$ as given in formula (4.26), and write $\psi ^{S}_{\eta }$ for the Gaussian random holomorphic section defined by formula (2.8), which can be regarded as a random variable valued in $\mathcal {B}_{f}(X,L)$ . Then the distribution of the identity on $(H^{0}_{(2)}(X,L),\mathbb {P}_{f})$ coincides with the distribution of the random section $\widehat {T}_{f}\psi ^{S}_{\eta }$ .

Remark 4.16. The Fubini–Study currents (1.17) associated to f have the following geometric interpretation, analogous to that of the usual Fubini–Study currents (1.4); see Remark 2.7. Assume that that $H^0_{(2)}(X,L)$ has no base locus, and consider the Toeplitz–Kodaira map $\Phi _f:X\to \mathbb {C}\mathbb {P}^{d-1}$ , given by $\Phi _f(x)=[(T_fS_j(x)/e_L(x))_j]$ , where $e_L$ is a local holomorphic frame around x ( $d=\infty $ is allowed). As in Remark 2.7, formula (4.34) yields $\Phi _f^*\omega _{\mathrm {FS}}=\gamma _{f}(L,h_{L})$ . Thus, the Fubini–Study current $\gamma _{f}(L,h_{L})$ is the pull-back of the Fubini–Study metric by the Toeplitz–Kodaira map $\Phi _f$ .

Remark 4.17. Note that in the above constructions we consider a nonnegative real function f in order to guarantee the injectivity of $T_{f}$ on $H^{0}_{(2)}(X,L)$ . We can also consider another setting. Here, we do not need the injectivity of $T_{f}$ . Let $f\in \mathcal {Q}(X,L;\mathbb {R})$ , which can be negative on a subset of X. Set

(4.37) $$ \begin{align} H^{0}_{(2)}(X,L,f):=(\ker T_{f})^{\perp}=\overline{T_{f}H^{0}_{(2)}(X,L)}\subset H^{0}_{(2)}(X,L). \end{align} $$

It is a Hilbert space, and the sections in $H^{0}_{(2)}(X,L,f)$ are the $\mathcal {L}^{2}$ -holomorphic sections of L detected by f. We consider the (self-adjoint) Hilbert–Schimdt operator

(4.38) $$ \begin{align} T_{f}^{\sharp}:=T_{f}|_{H^{0}_{(2)}(X,L,f)}: H^{0}_{(2)}(X,L,f)\rightarrow H^{0}_{(2)}(X,L,f). \end{align} $$

Then we can proceed as in Subsection 4.3 to construct a Gaussian probability measure $\mathbb {P}_{f}^{\sharp }$ on $H^{0}_{(2)}(X,L,f)$ . Let $\boldsymbol {s}^{\sharp }$ denotes the corresponding random section in $H^{0}_{(2)}(X,L,f)$ given by identity map, then

(4.39) $$ \begin{align} \mathbb{E}^{\mathbb{P}^{\sharp}_{f}}\left[[\operatorname{\mathrm{Div}}(\boldsymbol{s}^{\sharp})]\right]=\gamma_{f}(L,h_{L}), \end{align} $$

where $\gamma _{f}(L, h_{L})$ is given by formula (1.17).

One step further, set $m(f):=\dim \ker T_{f}\in \mathbb {N}\cup \{\infty \}$ , and let $\{S_{j}\}_{j=1}^{m(f)}$ be an orthonormal Hilbert basis of $\ker T_{f}$ , then the Schwartz kernel of the orthogonal projection $P_{\ker T_{f}}$ is given as

(4.40) $$ \begin{align} P_{\ker T_{f}}(x,y)=\sum_{j=1}^{m(f)}S_{j}(x)\otimes(S_{j}(y))^{\ast}. \end{align} $$

Let $\psi _{\ker T_{f}}$ be a standard Gaussian holomorphic section constructed from the Hilbert space $\ker T_{f}$ . We define a Gaussian random holomorphic section as follows,

(4.41) $$ \begin{align} \hat{\boldsymbol{s}}=\psi_{\ker T_{f}}+\boldsymbol{s}^{\sharp}. \end{align} $$

If $m(f)<\infty $ , then $\hat {\boldsymbol {s}}$ is valued in $H^{0}_{(2)}(X,L)$ ; otherwise, it is almost never $\mathcal {L}^{2}$ -integrable on X. Then we have

(4.42) $$ \begin{align} \mathbb{E}\left[[\operatorname{\mathrm{Div}}(\hat{\boldsymbol{s}})]\right]=c_{1}(L,h_{L})+ \frac{\sqrt{-1}}{2\pi}\partial\bar{\partial}\log \big(T^{2}_{f}(x,x) +P_{\ker T_{f}}(x,x)\big). \end{align} $$

Note that since f is bounded on X, then we always have

(4.43) $$ \begin{align} T^{2}_{f}(x,x)\leq T^{2}_{f}(x,x)+P_{\ker T_{f}}(x,x)\leq\max\{\|f\|^{2}_{\infty},1\}P(x,x), \end{align} $$

where $\|f\|_{\infty }$ is the $\mathcal {L}^{\infty }$ -norm of f on X.

We will consider the above different settings in Subsection 5.4 for high tensor powers of a prequantum line bundle on a complete Kähler manifold.

5 Random $\mathcal {L}^{2}$ -holomorphic sections for high tensor powers

By analogy with Section 3, we want to study the asymptotic behavior of the zeros of the random $\mathcal {L}^{2}$ -holomorphic sections for high tensor powers of a given positive line bundle on X. We make the same assumptions for $(X,\Theta )$ and $(L,h_{L})$ as in the beginning of Section 3 (or in Subsection 1.2), in particular, we assume the condition (1.6) and the completeness of $g^{TX}$ .

To construct in a natural way the sequence of random $\mathcal {L}^{2}$ -holomorphic sections of $L^{p}$ , $p\in \mathbb {N}_{>0}$ , we use the Toeplitz operators $\{T_{f,p}\}_{p\in \mathbb {N}_{>0}}$ associated with a suitable positive function f on X. Such operators $\{T_{f,p}\}_{p\in \mathbb {N}_{>0}}$ are already well studied in the context of Berezin–Toeplitz quantization.

5.1 Asymptotics of Toeplitz operators

Recall that $P_{p}$ denotes the orthogonal projection from $\mathcal {L}^{2}(X,L^{p})$ onto $H^{0}_{(2)}(X,L^{p})$ . For a smooth bounded function f on X and $p\in \mathbb {N}_{>0}$ , we set

(5.1) $$ \begin{align} T_{f,p}=P_{p}f P_{p}. \end{align} $$

This defines a bounded linear operator acting on $H^{0}_{(2)}(X,L^{p})$ .

To obtain the asymptotic expansion of the Schwartz kernels of $\{T_{f,p}\}$ , we need further assumptions either on the function f or on the geometry of X and L. We are mainly concerned with the following two cases (note that we always assume $g^{TX}$ to be complete).

  1. (I) The geometric data $(X,J,\Theta )$ and $(L,h_{L})$ are as in Condition 1.2 and the function f be smooth and constant outside a compact subset of X.

  2. (II) The geometric data $(X,J,\Theta )$ and $(L,h_{L})$ are as in Condition 1.2 and have bounded geometry (cf. Definition 3.2), and f is a bounded smooth function on X with (globally) bounded derivatives (with respect to $\nabla ^{TX}$ and $g^{TX}$ ) of any order.

Theorem 5.1 (cf. [Reference Ma and Marinescu56, Chapter 7],[Reference Ma and Marinescu60], [Reference Finski33, Lemmas 3.11, 3.14 & 4.6]).

Assume that either (I) or (II) hods. Then the Toeplitz operator $\{T_{f,p}\}_{p\in \mathbb {N}}$ has the following properties:

(i) For every compact subset $K\subset X$ , every $\epsilon>0$ , and every $\ell ,m\in \mathbb {N}$ , there exists $C_{\ell ,m,\varepsilon }>0$ such that for $p\geq 1$ , $x,x'\in X$ with $d(x,x')>\varepsilon $ , we have

(5.2) $$ \begin{align} |T_{f,p}(x,x')|_{\mathscr{C}^{m}(K\times K)}\leq C_{\ell,m,\epsilon} p^{-\ell}, \end{align} $$

where the $\mathscr {C}^{m}$ -norm is induced by $\nabla ^{TX}$ , and $h_{L}$ , $g^{TX}$ .

(ii) We have the uniform asymptotic expansion as $p\rightarrow \infty $ on any compact subset of X,

(5.3) $$ \begin{align} T_{f,p}(x,x)=\sum_{\ell=0}^{\infty} \boldsymbol{b}_{\ell}(f)(x)p^{n-l} +\mathcal{O}(p^{-\infty}), \end{align} $$

where $\boldsymbol {b}_{\ell }(f)\in \mathscr {C}^{\infty }(X,\mathbb {C})$ , in particular,

(5.4) $$ \begin{align} \boldsymbol{b}_{0}(f)(x)=\boldsymbol{b}_{0}(x)f(x),\, x\in X. \end{align} $$

(iii) The operator norms of $T_{f,p}$ , $p\in \mathbb {N}$ , satisfy

(5.5) $$ \begin{align} \lim_{p\rightarrow \infty} ||T_{f,p}||=||f||_{\infty}. \end{align} $$

(iv) If g is a bounded smooth function on X in the same class as f (Case (I) or (II)), then on any compact subset $K\subset X$ , we have the uniform expansion

(5.6) $$ \begin{align} (T_{f,p}T_{g,p})(x,x)=p^{n}\boldsymbol{b}_{0}(x)f(x)g(x) + \mathcal{O}(p^{n-1}), \end{align} $$

the expansion still holds if we take the derivatives with respect to x of any given order on both sides. In particular, in Case (II), we can refine estimate (5.2) to an exponential decay with respect to $\sqrt {p}$ , and the results (5.2), (5.3) and (5.6) hold uniformly on the whole manifold X.

The above theorem for Case (I) was mainly proved by Ma and Marinescu in [Reference Ma and Marinescu56, Chapter 7], [Reference Ma and Marinescu57]. For Case (II), it can be proved by a variation of the arguments in [Reference Ma and Marinescu56, Chapter 7] by using the exponential estimate for the Bergman kernel obtained in [Reference Ma and Marinescu60]. This is explained by Finski in [Reference Finski33, Sections 3 & 4].

Remark 5.2. The condition that f has bounded derivatives of any order in case (II) is used to guarantee the full expansion (5.3) and can be relaxed for our purposes. If $(X,J,\Theta )$ and $(L,h_{L})$ are as in Condition 1.2 and have bounded geometry, and if we assume that two smooth functions f, g have bounded derivatives up to order $2$ , then the expansion (5.6) still holds uniformly on any given compact subset in the $\mathscr {C}^{2}$ -norm (cf. the proof of [Reference Ma and Marinescu56, Lemma 7.2.4] and [Reference Barron, Ma, Marinescu and Pinsonnault3, Theorem 3.5 and Remark 3.6]).

Remark 5.3. Assume $0\neq f\in \mathscr {C}^{\infty }_{0}(X,\mathbb {R}_{\geq 0})$ to have compact support, hence $T_{f,p}$ is compact for $p\geq 1$ . We use the spectrum of the Toeplitz operators $T_{f,p}$ to define probability measures $\mathbb {P}_{f,p}$ on $H^0_{(2)}(X,L^p)$ (Definitions 4.11 and 5.6). It is thus interesting to know what is the asymptotic distribution of the spectra $\sigma (T_{f,p})$ as $p\to \infty $ . They play a role in several areas of geometric quantization [Reference Berndtsson7, Reference Ioos, Lu, Ma and Marinescu41, Reference Lindholm51]. On $(0, |f|_{\mathscr {C}^0}]$ , the spectral density measure $\mu _{f,p}$ of $T_{f,p}$ is defined as the sum of the Dirac masses at all the eigenvalues (counted with multiplicities) in $\sigma _{\mathrm {p}}(T_{f,p})$ , which is locally finite. A result of [Reference Ioos, Lu, Ma and Marinescu41] shows that as $p\rightarrow +\infty $ , we have the weak convergence of measures on $(0, |f|_{\mathscr {C}^0}]$ ,

(5.7) $$ \begin{align} p^{-n}\mu_{f,p}\rightarrow f_*\left(\frac{1}{n!}c_1(L,h_L)^n\right). \end{align} $$

This extends the results for compact Kähler manifolds or domains in $\mathbb {C}^n$ , such as [Reference Berndtsson7, Reference Lindholm51].

Our results in the sequel will mainly employ the expansion (5.6) with $g=f$ . Note that with further geometric conditions on $(X,\Theta )$ and $(L,h_{L})$ , we have a refined version of formula (5.6). Let $\mathrm {Ric}$ denote the Ricci curvature tensor, and set $\mathrm {Ric}_{\Theta }=:\mathrm {Ric}(J\cdot ,\cdot )$ . Let $\boldsymbol {r}^{X}$ denote the scalar curvature of $(X, g^{TX})$ , and let $\Delta $ be the (positive) Bochner Laplacian associated with $g^{TX}$ acting on the functions. We will use $\langle \cdot ,\cdot \rangle $ to denote the $\mathbb {C}$ -linear extension of the inner product $g^{\Lambda ^{\bullet }T^{*}X}$ . Consider the connection $\nabla ^{T^{*}X}: \mathscr {C}^{\infty }(X, T^{*}X\otimes \mathbb {C})\rightarrow \mathscr {C}^{\infty }(X, T^{\ast }X\otimes T^{\ast }X\otimes \mathbb {C})$ , let $D^{0,1}$ , $D^{1,0}$ denote the its respective $(1,0)$ , $(0,1)$ components.

The following theorem was proved in [Reference Ma and Marinescu59] for a compact Kähler manifold equipped with a prequantum line bundle. It was showed in [Reference Ma and Marinescu56, Reference Ma and Marinescu57, Reference Ma and Marinescu58, Reference Ma and Marinescu59] that the problem can be localized and since the computations are local the result extends to the case of complete (noncompact) Kähler manifolds. In particular, as a consequence of [Reference Ma and Marinescu56, Sections 7.4 & 7.5] (for the Case (I)) and [Reference Ma and Marinescu60] [Reference Finski33, Sections 3 & 4] (for Case (II)), these results hold for both our cases (I), (II).

Theorem 5.4. Assume that $(X,J,\Theta )$ is connected, complete Kähler manifold and that $(L,h_{L})$ is a prequantum holomorphic line bundle on X with smooth Hermitian metric $h_{L}$ , that is, $\Theta =c_{1}(L,h_{L})$ . Let $f,g$ be bounded smooth functions where are constants outside a compact subset (Case (I)), or if in addition $(X,\Theta ), (L,h_{L})$ have the bounded geometry, let $f,g$ be two bounded smooth functions on X such that their derivatives of any order are also bounded on X (Case (II)). Then for $\ell \in \mathbb {N}$ , there exists a smooth function on X, denoted by $\boldsymbol {b}_{\ell }(f,g)$ , which is a polynomial in the derivatives of $f,g$ with coefficients depending only on $\Theta $ and $h_{L}$ , such that on any compact subset $K\subset X$ , we have the uniform expansion as follows ( $N\geq 0$ ),

(5.8) $$ \begin{align} (T_{f,p}T_{g,p})(x,x)=\sum_{\ell=0}^{N} p^{n-\ell}\boldsymbol{b}_{\ell}(f,g)(x)+\mathcal{O}(p^{n-N-1}). \end{align} $$

Furthermore, we have

(5.9) $$ \begin{align} \begin{aligned} \boldsymbol{b}_{0}(f,g)&=fg,\\ \boldsymbol{b}_{1}(f,g)&=\frac{\boldsymbol{r}^{X}}{8\pi}fg-\frac{1}{4\pi}\big((\Delta f)g+f(\Delta g)\big)+\frac{1}{2\pi}\langle \overline{\partial}f,\partial g\rangle,\\ \boldsymbol{b}_{2}(f,g)&=\frac{1}{32\pi^{2}}\Big(f(\Delta^{2}g)+(\Delta^{2}f)g- \boldsymbol{r}^{X}\big(f(\Delta g)+(\Delta f)g\big) \Big)\\ &\quad -\frac{\sqrt{-1}}{8\pi^{2}}\langle \mathrm{Ric}_{\Theta}, f\partial\bar{\partial}g+g\partial\bar{\partial}f\rangle\\ &\quad +\frac{1}{8\pi^{2}}\Big\{\frac{1}{2}\Delta f\cdot\Delta g+\frac{\boldsymbol{r}^{X}}{2}\langle \overline{\partial}f,\partial g\rangle +\langle D^{0,1}\overline{\partial}f,D^{1,0}\partial g\rangle_{g^{T^{\ast}X\otimes T^{\ast}X}}\\ &\quad -\langle \overline{\partial}\Delta f,\partial g\rangle-\langle \overline{\partial}f,\partial \Delta g\rangle\Big\}. \end{aligned} \end{align} $$

Remark 5.5. Actually, by [Reference Ma and Marinescu59, Remark 0.5], a version of Theorem 5.4 holds also, without the prequantum condition $\Theta =c_{1}(L,h_{L})$ . In this case, we have

(5.10) $$ \begin{align} (T_{f,p}T_{g,p})(x,x)=\sum_{\ell=0}^{N} p^{n-\ell}\boldsymbol{b}^{\ast}_{\ell}(f,g)(x)+\mathcal{O}(p^{n-N-1}). \end{align} $$

with explicit formulas for $\boldsymbol {b}^{\ast }_{0}(f,g)$ , $\boldsymbol {b}^{\ast }_{1}(f,g)$ , $\boldsymbol {b}^{\ast }_{2}(f,g)$ (cf. [Reference Ma and Marinescu59, (0.30)]), which generalize the expressions (5.9). We will come back to this point in Remark 5.15.

5.2 Random zeros on the support: proofs of Theorems 1.8 and 1.9

Fix a $p_{0}\in \mathbb {N}_{>0}$ , set

(5.11) $$ \begin{align} \mathcal{Q}_{\geq p_{0}}(X,L;\mathbb{R}_{\geq 0}):=\cap_{p\geq p_{0}} \mathcal{Q}(X,L^{p};\mathbb{R}_{\geq 0}). \end{align} $$

Now, we consider a nontrivial $f\in \mathcal {Q}_{\geq p_{0}}(X,L;\mathbb {R}_{\geq 0})$ . Note that such function always exists, take for instance a nonnegative smooth function with compact support, or in the case of Bargmann–Fock space, a nonnegative Schwartz function on $\mathbb {C}^{n}$ . In the rest of this section, we always consider an integer $p\geq p_{0}$ .

Definition 5.6. Following the construction in Definition 4.11, let $\mathbb {P}_{f,p}$ be the probability measure on $H^{0}_{(2)}(X,L^{p})$ induced by the spectral decomposition of the Toeplitz operator $T_{f,p}$ . We will denote by $\boldsymbol {S}_{f,p}$ the random section in $H^{0}_{(2)}(X,L^{p})$ given by the probability distribution $(H^{0}_{(2)}(X,L^{p}),\mathbb {P}_{f,p})$ .

Let U be an open subset of X, and let $\Omega ^{n-1,n-1}_{0}(U)$ denote the smooth $(n-1,n-1)$ -forms on $\overline {U}$ with compact support in U. For any $(1,1)$ -current $\alpha $ on X, let $\alpha |_{U}$ denote its restriction on U by acting on sections in $\Omega ^{n-1,n-1}_{0}(U)$ .

Theorem 5.7. Assume that $(X,J,\Theta )$ , $(L,h_{L})$ , and f are as in one of the cases (I) or (II), with f a nontrivial function in $\mathcal {Q}_{\geq p_{0}}(X,L;\mathbb {R}_{\geq 0})$ . Let U be an open subset of X such that $f>0$ on U. Then for every compact subset of $K\subset U$ , there exists $p_{K}\in \mathbb {N}$ such that for all $p\geq p_{K}$ , $\gamma _{f}(L^{p},h_{p})$ is smooth on K, and we have the expansion of smooth forms on K in any $\mathscr {C}^{\ell }$ -norm,

(5.12) $$ \begin{align} \gamma_{f}(L^{p}, h_{p})=pc_{1}(L,h_{L})+\frac{\sqrt{-1}}{2\pi}\partial\overline\partial\log(\boldsymbol{b}_{0}f^2) + \mathcal{O}(p^{-1})\,,\quad p\to+\infty. \end{align} $$

Proof. By formula (1.17), we get

(5.13) $$ \begin{align} \gamma_{f}(L^{p},h_{p})=pc_{1}(L,h_{L})+ \frac{\sqrt{-1}}{2\pi}\partial\overline{\partial}\log{T_{f,p}^{2}(x,x)}. \end{align} $$

For any compact subset $K\subset U$ , set $m_{K}:=\max _{x\in K}\boldsymbol{b}_{0}(x)f^2(x)$ , $c_{K}:=\min _{x\in K}\boldsymbol{b}_{0}(x) f^2(x)>0$ . By formula (5.6), for $\ell \in \mathbb {N}$ , the following identity holds uniformly in any $\mathscr {C}^{\ell }(K)$ -norm,

(5.14) $$ \begin{align} T^{2}_{f,p}(x,x)=\boldsymbol{b}_{0}(x)f^{2}(x)p^{n}+\mathcal{O}(p^{n-1}), \:\:x\in K. \end{align} $$

If we are in Case (II), the asymptotics (5.14) hold uniformly over the whole manifold X.

As a consequence, there exists $p_{K}\in \mathbb {N}$ such that for $p\gg p_{K}$ , $x\in K$ , we have

(5.15) $$ \begin{align} 2m_{K}p^{n}\geq T_{f,p}^{2}(x,x)\geq \frac{1}{2}c_{K}p^{n}. \end{align} $$

Then for $p\geq p_{K}$ , $\gamma _{f}(L^{p},h_{p})$ is a smooth form on an open neighborhood of K. Then by the uniform expansion (5.14) and estimate (5.15), we get, in any $\mathscr {C}^{\ell }(K)$ -norm,

(5.16) $$ \begin{align} \frac{\sqrt{-1}}{2\pi }\partial\overline{\partial}\log{T_{f,p}^{2}(x,x)}= \frac{\sqrt{-1}}{2\pi}\partial\overline\partial\log(\boldsymbol{b}_{0}f^2)+\mathcal{O}(p^{-1})\,, \quad p\to+\infty. \end{align} $$

This way, we conclude the expansion (5.12).

Remark 5.8. (1) The expansion (5.12) is a generalization of Tian’s approximation theorem for a given integral Kähler form ([Reference Tian76, Theorem A], [56, Theorem 5.1.4]) for the case of the Toeplitz–Kodaira map (cf. Remark 4.16).

(2) By the arguments in the proof of [Reference Ma and Marinescu56, Lemma 5.1.6] and the conclusions drawn in Theorem 5.7, we infer that for any relatively compact open subset $V\subset X$ such that $f>0$ on $\overline {V}$ , the Toeplitz–Kodaira map restricted to V,

$$ \begin{align*}\Phi_{f,p}|_{V}: V\subset X\rightarrow \mathbb{CP}^{d_{p}-1},\end{align*} $$

is an immersion for p large enough. This remains true even when $d_{p}$ is infinite.

(3) On an open subset where $\boldsymbol{b}_{0}f^2$ is a nonzero positive constant, we have a sharper remainder for the approximations

(5.17) $$ \begin{align} \frac{1}{p}\gamma_{f}(L^{p},h_{p})=c_{1}(L,h_{L})+ \mathcal{O}(p^{-2})\,,\quad p\to+\infty. \end{align} $$

(4) For a compact Kähler manifold X with a prequantum line bundle $(L,h_L)$ , Ancona–Le Floch [Reference Ancona and Le Floch1] proved the approximation theorem in sense of currents under the condition that f vanishes transversally on X:

$$ \begin{align*}\gamma_{f}(L^{p},h_{p})-pc_{1}(L,h_{L})\rightarrow \frac{\sqrt{-1}}{2\pi}\partial\overline\partial\ \log f^2,\; \; p\rightarrow +\infty.\end{align*} $$

Combining Theorem 1.6 with Theorem 5.7, we get the following results.

Theorem 5.9. Assume that $(X,J,\Theta )$ , $(L,h_{L})$ and f are as in one of the cases (I) or (II), with f a nontrivial function in $\mathcal {Q}_{\geq p_{0}}(X,L;\mathbb {R}_{\geq 0})$ . Let U be an open subset of X such that $f>0$ on U. Let V be an open subset of U which is relatively compact in X, then there exists $p_{V}\in \mathbb {N}$ such that for all $p\geq p_{V}$ $\mathbb {E}^{\mathbb {P}_{f,p}}[[\operatorname {\mathrm {Div}}(\boldsymbol {S}_{f,p})]|_{V}]$ is the smooth form $\gamma _{f}(L^{p},h_{p})|_{V}$ on V, and we have the convergence of smooth forms on V in any $\mathscr {C}^{\ell }$ -norms as $p\rightarrow +\infty $ ,

(5.18) $$ \begin{align} \frac{1}{p}\mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{V}] \rightarrow c_{1}(L,h_{L})|_{V}. \end{align} $$

Moreover, we have the expansion in any $\mathscr {C}^{\ell }$ -norms on V for $p\gg 0$ ,

(5.19) $$ \begin{align} \mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{V}]-pc_{1}(L,h_{L})|_{V} =\frac{\sqrt{-1}}{2\pi}\partial\overline\partial\log(\boldsymbol{b}_0f^2)+\mathcal{O}(p^{-1}). \end{align} $$

As a consequence, we have the weak convergence of $(1,1)$ -currents on U as $p\rightarrow +\infty $ ,

(5.20) $$ \begin{align} \begin{aligned} \frac{1}{p}\mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}]&\rightarrow c_{1}(L,h_{L})|_{U},\\ \mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}]- pc_{1}(L,h_{L})|_{U}&\rightarrow \frac{\sqrt{-1}}{2\pi}\partial\overline\partial\log(\boldsymbol{b}_0f^2). \end{aligned} \end{align} $$

Remark 5.10. In line with Remark 5.2, for the case of bounded geometry, to show the weak convergence in terms of currents in the limit (5.20), it is sufficient to require that the function f has bounded derivatives up to order $2$ , since we only need the asymptotics (5.14) to hold in local $\mathscr {C}^{2}$ -norms.

By considering sequences of random sections in the product probability space,

(5.21) $$ \begin{align} (\boldsymbol{S}_{f,p})_{p}\in \prod_{p}\big(H^{0}_{(2)}(X,L^{p}),\mathbb{P}_{f,p}\big), \end{align} $$

we also have the following convergence with probability one.

Theorem 5.11. Let f be as in Theorem 5.9. Let U be an open subset of X such that $f>0$ on U, then for any $\varphi \in \Omega ^{n-1,n-1}_{0}(U)$ , we have

(5.22) $$ \begin{align} \mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}\big\langle[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})], \varphi \big\rangle=\langle c_{1}(L,h_{L}),\varphi\rangle \right)=1. \end{align} $$

If there exists a smooth closed positive $(1,1)$ -form $\alpha $ on X such that $\int _{X}c_{1}(L,h_{L})\wedge \alpha ^{n-1}<\infty $ and that $\alpha $ is strictly positive on $\overline {U}$ , then we have

(5.23) $$ \begin{align} \mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}=c_{1}(L,h_{L})|_{U} \right)=1, \end{align} $$

where the limit is taken with respect to the weak convergence of $(1,1)$ -currents on U.

Proof. Fix a nonzero $\varphi \in \Omega ^{n-1,n-1}_{0}(U)$ . Note that from the proof of Theorem 5.7, we have the convergence

(5.24) $$ \begin{align} \lim_{p\rightarrow\infty}\Big\langle \frac{1}{p}\gamma_{f}(L^{p},h_{p}),\varphi\Big\rangle=\langle c_{1}(L,h_{L}),\varphi\rangle. \end{align} $$

Defining the random variable

(5.25) $$ \begin{align} Y_{f,p}=\frac{1}{p}\Big\langle [\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]-\gamma_{f}(L^{p},h_{p}),\varphi\Big\rangle, \end{align} $$

the statement (5.22) is equivalent to proving that almost surely one has

(5.26) $$ \begin{align} Y_{f,p}\rightarrow 0. \end{align} $$

Note that if we use the construction from the proof of Theorem 1.6, we can write

(5.27) $$ \begin{align} \boldsymbol{S}_{f,p}=\sum_{j=1}^{d_{p}}\eta^{p}_{j}\lambda^{p}_{j}S^{p}_{j}, \end{align} $$

where $\{\eta ^{p}_{j}\}_{j}$ is a sequence of i.i.d. standard complex Gaussian random variables, $\{\lambda ^{p}_{j}\}_{j}$ is the point spectrum of $T_{f,p}$ , and $\{S^{p}_{j}\}_{j}$ is the orthonormal basis of $H^{0}_{(2)}(X,L^{p})$ given by the eigensections of $T_{f,p}$ . Then, as explained in Remark 3.10, we can proceed as in the proof of Theorem 3.7, and we get $\mathbb {E}[|Y_{f,p}|^{2}]=\mathcal {O}(p^{-2})$ , which entails the convergence (5.26), and hence equality (5.22).

Now, we prove (5.23). Note that by our assumptions on $\alpha $ , the analogous arguments in the proof of Corollary 3.8 show that there exists $C_{U,\alpha }>0$ such that for all $\varphi \in \Omega ^{n-1,n-1}_{0}(U)$ ,

(5.28) $$ \begin{align} \frac{1}{p}|\langle[\operatorname{\mathrm{Div}}(s_{p})],\varphi\rangle|\leq C_{U,\alpha}|\varphi|_{\mathscr{C}^{0}(U,\alpha)}, \end{align} $$

where $\mathscr {C}^{0}(U,\alpha )$ denotes the $\mathscr {C}^{0}$ -norm for differential forms on U with respect to the Hermitian metric induced by $\alpha $ . Then the equality (5.23) follows as a consequence of the equality (5.22) by considering a countable $\mathscr {C}^{0}(U,\alpha )$ -dense family of $\varphi $ ’s in $\Omega ^{n-1,n-1}_{0}(U)$ .

It is natural to investigate a relaxation of the assumptions of Theorem 5.11 as follows. For f as above, consider U an open subset of $\operatorname {\mathrm {supp}} f$ . In general, f might vanish at some points in $U,$ and it is a natural and interesting question to understand for which kind of conditions on the zeros of f in U we still can have the equidistribution results for the random zeros on U as above. Since f is nonnegative, if $f(x_0)=0$ , the least possible vanishing order of f at $x_{0}$ is $2$ . In the sequel, we will explain that, if f has only zeros of order $2$ at which $\Delta f$ does not vanish, then the above results still hold. For this purpose, we will employ the results in Theorem 5.4 so that we need the prequantum condition.

Proposition 5.12. Assume that $(X,J,\Theta )$ , $(L,h_{L})$ , f are as in either cases (I) or (II) with f a nontrivial function in $\mathcal {Q}_{\geq p_{0}}(X,L;\mathbb {R}_{\geq 0})$ and the prequantum condition $\Theta =c_{1}(L,h_{L})$ holds. Let U be an open subset of $X.$ If f only vanishes up to order $2$ in U and $\Delta f$ is nonzero at all zeros of f, then for any compact subset K of U, there exists a constant $c_{K}>0$ and $p_{K}\geq p_{0}$ such that for $x\in K$ , $p\geq p_{K}$ ,

(5.29) $$ \begin{align} T^{2}_{f,p}(x,x)\geq c_{K}p^{n-2}. \end{align} $$

Moreover, $\log f^{2}$ is locally integrable on U, and we have the weak convergence of currents on U as $p\rightarrow +\infty $ ,

(5.30) $$ \begin{align} \partial\overline{\partial} \log T^{2}_{f,p}(x,x)\rightarrow \partial\overline{\partial} \log f^{2}. \end{align} $$

Around a point x where $f(x)>0$ , the convergence in (5.30) holds in any local $\mathscr {C}^{\ell }$ -norms.

Proof. With our assumptions, we can apply Theorem 5.4 to $T^{2}_{p,f}$ . Let $x_{0}\in U$ with $f(x_0)=0$ . By our assumption on f, we have

(5.31) $$ \begin{align} -\Delta f(x_{0})\neq 0. \end{align} $$

By taking a geodesic normal coordinate system $(Y=(y_{j})_{j=1}^{2n}\in \mathbb {R}^{2n})$ centered at $x_{0}$ , the Taylor expansion of f near $x_{0}$ reads

(5.32) $$ \begin{align} f(Y)=\sum_{j}c_{j}(x_{0})y_{j}^{2}+\mathcal{O}(|Y|^{3}), \end{align} $$

where $c_{j}(x_{0})\ge 0$ since $f\geq 0$ . Then

(5.33) $$ \begin{align} -\Delta f(x_{0})=\sum_{j}c_{j}(x_{0})>0. \end{align} $$

Now, we compute the terms $\boldsymbol {b}_{\ell }(f,f)$ , $\ell =1,2$ , from (5.9) near $x_{0}$ ,

(5.34) $$ \begin{align} \begin{aligned} &\boldsymbol{b}_{1}(f,f)=\frac{1}{8\pi}(\boldsymbol{r}^{X}f-4\Delta f)f+\frac{1}{2\pi}|\partial f|^{2},\\ &\boldsymbol{b}_{2}(f,f)=\frac{1}{4\pi^{2}}\Big(\sum_{j}c_{j}(x_{0}) \Big)^{2}+\frac{1}{8\pi^{2}}|D^{0,1}\overline{\partial}f(x_{0})|^{2}_{g^{T^{\ast}X\otimes T^{\ast}X}}+\mathcal{O}(|Y|). \end{aligned} \end{align} $$

Setting

(5.35) $$ \begin{align} \begin{aligned} \mu(f,x_{0})=\frac{1}{4\pi^{2}} \Big(\sum_{j}c_{j}(x_{0})\Big)^{2}+\frac{1}{8\pi^{2}}|D^{0,1}\overline{\partial}f(x_{0})|^{2}_{g^{T^{\ast}X\otimes T^{\ast}X}}>0, \end{aligned} \end{align} $$

we can choose a small open neighborhood $V_{x_{0}}$ of $x_{0}$ such that for $x\in V_{x_{0}}$ ,

(5.36) $$ \begin{align} \boldsymbol{r}^{X}_{x}f(x)-4\Delta f(x)\geq 0, \quad \text{and} \quad \boldsymbol{b}_{2}(f,f)(x)\geq \frac{1}{2} \mu(f,x_{0}), \end{align} $$

and so

(5.37) $$ \begin{align} \boldsymbol{b}_{1}(f,f)(x)\geq 0. \end{align} $$

Since $\boldsymbol {b}_{0}(f,f)=f^{2}$ , then from the above computations and estimate (5.15), we get estimate (5.29).

By formula (5.32), on a sufficiently small open neighborhood of $x_{0}$ , we have

(5.38) $$ \begin{align} f(Y)\geq \frac{1}{2}\sum_{j} c_{j}(x_{0})y_{j}^{2}. \end{align} $$

Then it is clear that $\log f^{2}$ is integrable near $x_{0}$ . Then the current $\partial \overline {\partial }\log f^{2}$ is well defined on U. Near a point where f does not vanish, we get the strong convergence in the limit (5.30) by means of asymptotics (5.6) and (5.16).

Next, we focus on a point $x_{0}$ with $f(x_{0})=0$ . Note that

(5.39) $$ \begin{align} p^{-n}T^{2}_{f,p}(x,x)=f^{2}+\boldsymbol{b}_{1}(f,f)p^{-1}+ \boldsymbol{b}_{2}(f,f)p^{-2}+\mathcal{O}(p^{-3}). \end{align} $$

By formulas (5.34), we can take a small open neighborhood $V^{\prime }_{x_{0}}$ of $x_{0}$ such that for $x\in V^{\prime }_{x_{0}}$ , $p\gg 0$ ,

(5.40) $$ \begin{align} \boldsymbol{b}_{1}(f,f)(x)p^{-1}+\boldsymbol{b}_{2}(f,f)(x)p^{-2}+ \mathcal{O}(p^{-3})\geq 0,\quad \text{and} \quad f^{2}(x)\leq p^{-n}T^{2}_{f,p}(x,x)\leq 1. \end{align} $$

Then on $V^{\prime }_{x_{0}}$ , we have

(5.41) $$ \begin{align} \big|\log(p^{-n}T^{2}_{f,p}(x,x))\big|\leq |\log f^{2}(x)|. \end{align} $$

At the same time we have the pointwise convergence of functions as $p\rightarrow \infty $ ,

(5.42) $$ \begin{align} \log(p^{-n}T^{2}_{f,p}(x,x))\rightarrow \log f^{2}(x). \end{align} $$

Since $\log f^{2}$ is integrable near $x_{0}$ , by the dominated convergence theorem, we get the convergence of $(1,1)$ -currents in (5.30) on $V^{\prime }_{x_{0}}$ , hence on U. This completes the proof.

Remark 5.13. In the proof of Proposition 5.12, we see that if f has at least one zero in $K\subset U$ , then the power $(n-2)$ in estimate (5.29) cannot be improved; otherwise, a lower bound of $T_{f,p}^{2}(x,x)$ on K is given in estimate (5.15). When X is compact, this observation indicates that if $f\geq 0$ has only proper zeros of order $2$ and at least one of such vanishing point, then the lowest eigenvalue of $T_{f,p}$ should behave like $\mathcal {O}(\frac {1}{p})$ as p grows. For this kind of results, we refer to the papers [Reference Deleporte19, Reference Deleporte20] of Deleporte. In particular, when X is compact, the lower bound in estimate (5.29) can be deduced from [Reference Deleporte19].

As a direct consequence of Proposition 5.12, we obtain:

Theorem 5.14. Assume that $(X,J,\Theta )$ , $(L,h_{L})$ , f are as in either cases (I) or (II) with f a nontrivial function in $\mathcal {Q}_{\geq p_{0}}(X,L;\mathbb {R}_{\geq 0})$ and the prequantum condition $\Theta =c_{1}(L,h_{L})$ holds. Let U be an open subset of $\operatorname {\mathrm {supp}} f$ such that f only vanishes up to second order in U with nonzero $\Delta f$ at the zeros. Then as $p\rightarrow \infty $ ,

(i) We have the weak convergence of $(1,1)$ -currents on U as $p\rightarrow +\infty $ ,

(5.43) $$ \begin{align} \begin{aligned} \frac{1}{p}\mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}]& \rightarrow c_{1}(L,h_{L})|_{U},\\ \mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}]- pc_{1}(L,h_{L})|_{U}&\rightarrow \frac{\sqrt{-1}}{2\pi}\partial\overline\partial\log f^2\,. \end{aligned} \end{align} $$

For any $\varphi \in \Omega ^{n-1,n-1}_{0}(U)$ , we have

(5.44) $$ \begin{align} \mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}\langle[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})],\varphi\rangle= \langle c_{1}(L,h_{L}),\varphi\rangle\right)=1. \end{align} $$

(ii) If there is a smooth closed positive $(1,1)$ -form $\alpha $ on X with $\int _{X}c_{1}(L,h_{L})\wedge \alpha ^{n-1}<\infty $ such that $\alpha $ is strictly positive on $\overline {U}$ , then we have

(5.45) $$ \begin{align} \mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}=c_{1}(L,h_{L})|_{U} \right)=1, \end{align} $$

where the limit is taken with respect to the weak convergence of $(1,1)$ -currents on U.

Remark 5.15. Following Remark 5.5, we have more general versions of Proposition 5.12 and Theorem 5.14 without the prequantum condition. In fact, let $\boldsymbol {b}^{\ast }_{j}(f,g)$ , $j=0,1,2,\ldots $ , denote the coefficients in the general expansion (5.10), valid without assuming the prequantum condition. We still consider the same function f in Proposition 5.12, and a point $x_{0}\in U$ such that $f(x_{0})=0$ , $-\Delta f (x_{0})>0$ . Then by [Reference Ma and Marinescu59, (0.17), (0.30) and Remark 0.5], we get that, on a very small neighbourhood of $x_{0}$ , we have

(5.46) $$ \begin{align} \begin{aligned} &\boldsymbol{b}^{\ast}_{0}(f,f)(x)=\boldsymbol{b}_{0}(x)f(x)^{2}\geq 0,\\ &\boldsymbol{b}^{\ast}_{1}(f,f)(x)\geq 0,\\ &\boldsymbol{b}^{\ast}_{2}(f,f)(x)=\boldsymbol{b}_{0}(x) \mu(f,x_{0})+\mathcal{O}(|Y|), \end{aligned} \end{align} $$

where $\boldsymbol {b}_{0}(x)$ is the positive function given in (3.3), and Y denote the normal coordinates centered at $x_{0}$ . Note that for $ \boldsymbol {b}^{\ast }_{1}(f,f)(x)\geq 0$ , we use that $-\Delta f(x_{0})$ dominates the terms f, $\partial f/\partial y_{j}$ , $j=1,\ldots , 2n$ , near the point $x_{0}$ . Then we can proceed as in the proof of Proposition 5.12 to get the same conclusions. All the results in Theorem 5.14 still hold except that we need to change the second line in (5.43) to the following limit as $p\to +\infty $ ,

(5.47) $$ \begin{align} \mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}]- pc_{1}(L,h_{L})|_{U}\rightarrow \frac{\sqrt{-1}}{2\pi}\partial\overline\partial\log f^2\,+\frac{\sqrt{-1}}{2\pi}\partial\overline\partial\log\boldsymbol{b}_{0}. \end{align} $$

Example 5.16. We give here simple examples for the situation of Theorem 5.14.

(1) Let $X=\mathbb {D}$ be the unit disc in $\mathbb {C}$ endowed with the hyperbolic metric $\omega :=\omega _{B}$ given by Equation (2.41). This metric induces a metric $h_{K_{\mathbb {D}}}$ on the canonical bundle $K_{\mathbb {D}}$ with $c_1(K_{\mathbb {D}},h_{K_{\mathbb {D}}})=\omega $ , hence $(K_{\mathbb {D}},h_{K_{\mathbb {D}}})$ is a prequantum line bundle on $(\mathbb {D},\omega )$ . For $r\in (0,1)$ consider a cutoff function $\chi \in \mathscr {C}^{\infty }_0(\mathbb {D})$ such that $\chi =1$ on $B_r(0)$ . We define $f\in \mathscr {C}^{\infty }_0(\mathbb {D})$ by $f(z)=\chi (z)|z|^{2}$ . Then we have on $U=B_r(0)$ , with $\delta _0$ the delta distribution at $0$ ,

(5.48) $$ \begin{align} \frac{1}{p}\mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]|_{U}]\rightarrow \omega|_{U},\quad \mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})|_{U}]]- p\omega|_{U}\rightarrow 2\delta_0 \,,\:\:p\rightarrow \infty, \end{align} $$

since by Lelong–Poincaré formula we have $\frac {\sqrt {-1}}{2\pi }\partial \overline \partial \log |z|^2=\delta _0$ .

(2) Let $X=\mathbb {C}\mathbb {P}^1$ endowed with the Fubini–Study metric $\omega _{FS}$ and $(L,h_L)=(\mathscr {O}(1),h_{FS})$ be the hyperplane line bundle endowed with the Fubini–Study metric, so $c_1(L,h_L)=\omega _{FS}$ . We define $f\in \mathscr {C}^{\infty }(\mathbb {C}\mathbb {P}^1)$ in homogeneous coordinates $[z_0,z_1]$ on $\mathbb {C}\mathbb {P}^1$ by $f([z_0,z_1])=|z_0|^{2}/(|z_0|^2+|z_1|^2)$ (note that f is a rescaling of the height function on the sphere $S^2\subset \mathbb {R}^3$ ). Then we have on $\mathbb {C}\mathbb {P}^1$ ,

(5.49) $$ \begin{align} \frac{1}{p}\mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]]\rightarrow \omega_{FS},\quad \mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]]- (p+2)\omega_{FS}\rightarrow 2\delta_{[0,1]} \,,\:\:p\rightarrow \infty, \end{align} $$

since $\frac {\sqrt {-1}}{2\pi }\partial \overline \partial \log f= \delta _{[0,1]}-\omega _{FS}$ , with $\delta _{[0,1]}$ the delta distribution at the point $[0,1]\in \mathbb {C}\mathbb {P}^1$ .

(3) The previous example can be generalized as follows. Let X be a compact Kähler manifold, $(L,h_L)$ be a prequantum holomorphic line bundle on X. Let $D\subset X$ be a smooth complex hypersurface of X, and $F=\mathscr {O}(D)$ be the associated holomorphic line bundle. Let $s_D\in H^0(X,F)$ be a canonical holomorphic section of F, that vanishes at first order on D. It is uniquely determined up to a multiplicative constant. We endow F with a Hermitian metric $h_F$ and consider $f\in \mathscr {C}^\infty (X)$ , $f=|s_D|_{h_F}^2$ . By the Lelong–Poincaré formula (2.29), we have on X as $p\rightarrow \infty $ ,

(5.50) $$ \begin{align} \begin{aligned} \frac{1}{p}\mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]]&\rightarrow c_1(L,h_L),\\ \mathbb{E}^{\mathbb{P}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}_{f,p})]]- pc_1(L,h_L)&\rightarrow 2([D]-c_1(F,h_F))\,, \end{aligned} \end{align} $$

where $[D]$ is the current of integration on D.

These examples illustrate how one can recover the geometric input (the curvature form $c_1(L,h_L)$ ) and the classical observable f in the semiclassical limit $p\to \infty $ . By recovering the zero divisor D of $s_D$ in the limit (5.50), we have determined $s_D$ up to a constant factor.

5.3 Higher fluctuation of random zeros near points of vanishing order two

In this subsection, we work under the assumption that $(X,J,\Theta )$ , $(L,h_{L})$ , f are given in either cases (I) or (II) with f a nontrivial function in $\mathcal {Q}_{\geq p_{0}}(X,L;\mathbb {R}_{\geq 0})$ and with prequantum condition $\Theta =c_{1}(L,h_{L})$ . We aim to examine the random zeros of $\boldsymbol {S}_{f,p}$ near a proper zero of f with vanishing order $2$ , up to Planck scale $1/\sqrt {p}$ . Note that in [Reference Ancona and Le Floch1], for a compact Kähler manifold X and under a different assumption on f, Ancona and Le Floch proved that the random zeros fluctuate more near the zeros of f. We will prove a similar result for our setting. For this purpose, we need to refine the computations in formulas (5.34) in a complex coordinate system centered at $x_{0}$ where f vanishes with order $2$ .

Suppose $f\geq 0$ and that $x_{0}$ is a zero of f with $\Delta f(x_{0})<0$ . Then we can choose a holomorphic coordinate system centered at $x_{0}$ , denoted by $z=(z_{1},\ldots ,z_{n})\in \mathbb {C}^{n}$ such that

(5.51) $$ \begin{align} g^{TX}_{z}=g^{\mathbb{C}^{n}}_{\mathrm{st}}+\mathcal{O}(|z|^{2}), \end{align} $$

where $g^{\mathbb {C}^{n}}_{\mathrm {st}}$ denotes the standard Euclidean metric on $\mathbb {C}^{n}\simeq \mathbb {R}^{2n}$ . Note that we view z as a column vector, and let $(\cdot )^{T}$ denote the transpose of a matrix. In this coordinate system, we can write

(5.52) $$ \begin{align} f(z)=z^{T}A\bar{z}+z^{T}Bz+\bar{z}^{T}\bar{B}\bar{z}+\mathcal{O}(|z|^{3}), \end{align} $$

where the matrix A is Hermitian and semipositive definite, B is symmetric complex matrix and they are determined uniquely by the Hessian of f at $x_{0}$ . Set

(5.53) $$ \begin{align} \hat{f}_{x_{0}}(z)=z^{T}A\bar{z}+z^{T}Bz+\bar{z}^{T}\bar{B}\bar{z}. \end{align} $$

Since $f\geq 0$ , then for any $z\in \mathbb {C}^{n}$ with $\|z\|=1$ ,

(5.54) $$ \begin{align} z^{T}A\bar{z}\geq 2|\Re(z^{T}Bz)|, \end{align} $$

where $\Re (\cdot )$ denotes the real part. In particular, $\hat {f}_{x_{0}}(z)\geq 0$ . Using this complex coordinate system, we compute

(5.55) $$ \begin{align} \begin{aligned} \Delta f(z)&=-4\operatorname{\mathrm{Tr}}[A]+\mathcal{O}(|z|),\\ |\partial f(z)|^{2}&=2|A\bar{z}+2Bz|^{2}+\mathcal{O}(|z|^{3}),\\ |D^{0,1}\overline{\partial}f(z)|^{2}_{g^{T^{\ast}X\otimes T^{\ast}X}}&=16\operatorname{\mathrm{Tr}} [B\bar{B}^{T}]+\mathcal{O}(|z|). \end{aligned} \end{align} $$

Note that $\mu (f,x_{0})$ is defined in (5.35), then we have

(5.56) $$ \begin{align} \mu(f,x_{0})=\frac{1}{\pi^{2}}(\operatorname{\mathrm{Tr}} [A])^{2}+ \frac{2}{\pi^{2}}\operatorname{\mathrm{Tr}} [B\bar{B}^{T}]>0. \end{align} $$

Then we rewrite the computations in formulas (5.34) as follows:

(5.57) $$ \begin{align} \begin{aligned} &\boldsymbol{b}_{0}(f,f)(z)=\hat{f}_{x_{0}}^{2}(z)+\mathcal{O}(|z|^{5}),\\ &\boldsymbol{b}_{1}(f,f)(z)=\frac{2}{\pi}\operatorname{\mathrm{Tr}} [A]\hat{f}_{x_{0}}(z)+\frac{1}{\pi}|A\bar{z}+2Bz|^{2}+\mathcal{O}(|z|^{3}),\\ &\boldsymbol{b}_{2}(f,f)(z)=\mu(f,x_{0})+\mathcal{O}(|z|). \end{aligned} \end{align} $$

Definition 5.17. Associated with the Kähler form $\Theta $ and f near $x_{0}$ , we define a (strictly) positive function on $\mathbb {C}^{n}$ as follows:

(5.58) $$ \begin{align} F_{f,x_{0}}(z)=\hat{f}^{2}_{x_{0}}(z)-\frac{1}{2\pi}(\Delta f)(x_{0}) \hat{f}_{x_{0}}(z)+ \frac{1}{\pi}|A\bar{z}+2Bz|^{2} +\mu(f,x_{0}). \end{align} $$

Note that this function does not depend on the choice of the holomorphic coordinate systems centered at $x_{0}$ satisfying the equality (5.51). Equivalently, we have for $z\in \mathbb {C}^{n}\simeq (T_{x_{0}}X, J_{x_{0}})$ ,

(5.59) $$ \begin{align} F_{f,x_{0}}(z)=\lim_{p\rightarrow\infty} \big\{p^{2}\boldsymbol{b}_{0}(f,f)(z/\sqrt{p})+p\boldsymbol{b}_{1}(f,f)(z/\sqrt{p})+\boldsymbol{b}_{2}(f,f)(z/\sqrt{p})\big\}. \end{align} $$

We also define the following positive quadratic function in $z\in \mathbb {C}^{n}$ ,

(5.60) $$ \begin{align} \widehat{\boldsymbol{b}}_{1}(z)=\lim_{p\rightarrow\infty} p\boldsymbol{b}_{1}(f,f)(z/\sqrt{p})=-\frac{1}{2\pi}(\Delta f)(x_{0}) \hat{f}_{x_{0}}(z)+ \frac{1}{\pi}|A\bar{z}+2Bz|^{2}. \end{align} $$

Proposition 5.18. With the above notation,

(5.61) $$ \begin{align} \beta_{f,x_{0}}:=\partial\overline{\partial} \,\widehat{\boldsymbol{b}}_{1}=\partial\overline{\partial} F_{f,x_{0}}(0)\in \Lambda^{(1,1)}T^{\ast}_{x_{0}}X \end{align} $$

is a positive $(1,1)$ -form on $\mathbb {C}^{n}$ , more precisely,

(5.62) $$ \begin{align} \beta_{f,x_{0}}=(dz)^{T}\wedge K_{f,x_{0}}d\bar{z}, \end{align} $$

where $K_{f,x_{0}}$ is the semipositive definite Hermitian matrix given by

(5.63) $$ \begin{align} K_{f,x_{0}}=\frac{2}{\pi}\operatorname{\mathrm{Tr}} [A]A+\frac{1}{\pi}(A^{2}+4B\bar{B}). \end{align} $$

We have the convergence of $(1,1)$ -forms at $x_{0}$ as $p\rightarrow \infty $ ,

(5.64) $$ \begin{align} \frac{1}{p}\partial\overline{\partial}\log T^{2}_{f,p}(x,x)|_{x=x_{0}}\rightarrow \frac{1}{\mu(f,x_{0})}\beta_{f,x_{0}}=\partial\overline{\partial}\log{ F_{f,x_{0}}}(0). \end{align} $$

Proof. The first part follows directly from the formulae (5.58) and (5.60). We now prove the convergence (5.64). In complex coordinates z centered at $x_{0}$ , we have for $|z|<1$ ,

(5.65) $$ \begin{align} p^{-n}T^{2}_{f,p}(z,z)=p^{-2}F_{f,x_{0}}(\sqrt{p}z)+\mathcal{O}(p^{-3}) +\mathcal{O}(|z|^{5})+p^{-1}\mathcal{O}(|z|^{3})+p^{-2}\mathcal{O}(|z|), \end{align} $$

hence

(5.66) $$ \begin{align} \begin{aligned} &\frac{1}{p}\partial\overline{\partial}\log(p^{-n}T^{2}_{f,p}(z,z))\\ &=\frac{\big(\partial\overline{\partial}F_{f,x_{0}}\big)(\sqrt{p}z)+ \mathcal{O}(p^{-1})+p\mathcal{O}(|z|^{3})+ \mathcal{O}(|z|)}{F_{f,x_{0}}(\sqrt{p}z)+ \mathcal{O}(p^{-1})+p^{2}\mathcal{O}(|z|^{5})+ p\mathcal{O}(|z|^{3})+\mathcal{O}(|z|)}\\ &\quad-\frac{\big(\partial F_{f,x_{0}}\wedge\overline{\partial} F_{f,x_{0}}\big)(\sqrt{p}z)+p^{3}\mathcal{O}(|z|^{7})+ p^{2}\mathcal{O}(|z|^{5})+p\mathcal{O}(|z|^{3})+ \mathcal{O}(|z|)+\mathcal{O}(p^{-1})} {\big\{F_{f,x_{0}}(\sqrt{p}z)+\mathcal{O}(p^{-1})+ p^{2}\mathcal{O}(|z|^{5})+p\mathcal{O}(|z|^{3})+ \mathcal{O}(|z|)\big\}^{2}}. \end{aligned} \end{align} $$

Taking $z=0$ in the equality (5.66) and letting $p\rightarrow \infty $ yields (5.64).

Definition 5.19. For a zero $x_{0}$ of f as above, we define for $R>0$ the linear function

(5.67) $$ \begin{align} \Phi^{R}_{f,x_{0}}:\Lambda^{(n-1,n-1)}_{x_{0}}T^{\ast}X\rightarrow \mathbb{C}, \quad \Phi^{R}_{f,x_{0}}(\alpha):= \frac{\sqrt{-1}}{2\pi}\int_{B^{\mathbb{C}^{n}}(0,R)} \partial\overline{\partial}\log{ F_{f,x_{0}}(z)}\wedge\alpha, \end{align} $$

where $\alpha \in \Lambda ^{(n-1,n-1)}_{x_{0}}$ is viewed as a constant $(n-1,n-1)$ -form on $\mathbb {C}^{n}\simeq (T_{x_{0}}X, J_{x_{0}})$ .

Remark 5.20. The expression of $\Phi ^{R}_{f,x_{0}}(\alpha )$ can be made more concrete by using the formula (5.58), especially if f has a nice shape near $x_{0}$ (for example, if $B=0$ ). We will demonstrate this in Example 5.21, but our expectation is that the calculations in general will be much more complicated, so we will not attempt to do so in this paper.

Example 5.21. Now, we assume f near $x_{0}$ is given by expansion (5.52), where $B=0$ and

(5.68) $$ \begin{align} A=\mathrm{Id}_{n}. \end{align} $$

Then

(5.69) $$ \begin{align} F_{f,x_{0}}(z)=|z|^{4}+\frac{(2n+1)}{\pi}|z|^{2}+\frac{n^{2}}{\pi^{2}}. \end{align} $$

Set $\omega _{0}=\sqrt {-1}\sum _{j}dz_{j}\wedge d\bar {z}_{j}$ . Then we have

(5.70) $$ \begin{align} \begin{aligned} &\sqrt{-1}\partial\overline{\partial}\log{ F_{f,x_{0}}(z)}\wedge\frac{\omega_{0}^{n-1}}{(n-1)!}\\ &=\pi\Big[\frac{(2n-2)\pi^{3}|z|^{6}+(6n^{2}-n-2)\pi^{2}|z|^{4}+(6n^{3}+2n^{2}-3n-1)\pi|z|^{2}+2n^{4}+n^{3}}{\pi^{4}|z|^{8}+(4n+2)\pi^{3}|z|^{6}+(6n^{2}+4n+1)\pi^{2}|z|^{4}+(4n^{3}+2n^{2})\pi|z|^{2}+n^{4}} \Big]\frac{\omega_{0}^{n}}{n!}. \end{aligned} \end{align} $$

In the case $n=1$ , we have

(5.71) $$ \begin{align} \begin{aligned} &\sqrt{-1}\partial\overline{\partial}\log{ F_{f,x_{0}}(z)}=\pi\Big[\frac{3\pi^{2}|z|^{4}+4\pi|z|^{2}+3}{\pi^{4}|z|^{8}+6\pi^{3}|z|^{6}+11\pi^{2}|z|^{4}+6\pi|z|^{2}+1}\Big] \omega_{0}. \end{aligned} \end{align} $$

Theorem 5.22. Assume that $(X,J,\Theta )$ , $(L,h_{L})$ , f are as in either Case (I) or Case (II) with f a nontrivial function in $\mathcal {Q}_{\geq p_{0}}(X,L;\mathbb {R}_{\geq 0})$ and the prequantum condition $\Theta =c_{1}(L,h_{L})$ holds. Let $x_{0}$ be a zero of f with $\Delta f(x_{0})<0$ . Then for any fixed $R>0$ , $\varphi \in \Omega ^{n-1,n-1}_{0}(X)$ , and for all $p> p_{0}$ ,

(5.72) $$ \begin{align} \frac{\sqrt{-1}}{2\pi}\int_{B(x_{0},R/\sqrt{p})}\partial\overline{\partial} \log(T^{2}_{f,p}(x,x))\wedge\varphi= p^{-n+1}\Phi^{R}_{f,x_{0}}\big(\varphi({x_{0}})\big)+ \mathcal{O}(p^{-n+1/2}). \end{align} $$

Proof. For $p\gg 0$ , we identify $B(x_{0},R/\sqrt {p})\simeq B^{\mathbb {C}^{n}}(0,R/\sqrt {p})$ . Then for $z\in B^{\mathbb {C}^{n}}(0,R/\sqrt {p})$ , $l\in \mathbb {N}$ ,

(5.73) $$ \begin{align} p^{l}\mathcal{O}(|z|^{2l+1})=\mathcal{O}(p^{-1/2}). \end{align} $$

Also, note for $z\in B^{\mathbb {C}^{n}}(0,R)$ ,

(5.74) $$ \begin{align} \varphi(z/\sqrt{p})=\varphi(x_{0})+\mathcal{O}(p^{-1/2}). \end{align} $$

Then the asymptotics (5.72) follows from the asymptotics (5.66). This completes the proof.

As explained in Subsection 1.4, formula (5.72) gives different powers of p in the asymptotics (1.25), which exhibit two different behaviors of the fluctuations of the random zeros depending on whether $f(x)=0$ or $f(x)\neq 0$ .

5.4 Case of real functions with negative values

This subsection continues the discussion from Remark 4.17. We examine the equidistribution of random zeros of $\mathcal {L}^{2}$ -holomorphic sections detected by a given real function f, which is not necessarily nonnegative. Now, we consider the case of complete Kähler manifold $(X,\Theta )$ equipped with a prequantum holomorphic line bundle $(L,h_{L})$ . Recall that $\mathcal {Q}(X,L^{p};\mathbb {R})$ is the subspace of $\mathcal {Q}(X,L^{p};\mathbb {C})$ consisting of real-valued functions and that

(5.75) $$ \begin{align} \mathcal{Q}_{\geq p_{0}}(X,L;\mathbb{R}):=\cap_{p\geq p_{0}}\mathcal{Q}(X,L^{p};\mathbb{R}). \end{align} $$

Definition 5.23. Let f be a real smooth function on X, for $x\in X$ , we say f is vanishing properly at x up to order $2$ if one of the following cases holds:

  • $f(x)\neq 0$ , or

  • $f(x)=0$ , $df(x)\neq 0$ , or

  • $f(x)=0,\ df(x)=0,\ \Delta f(x)\neq 0$ with $f\Delta f \leq 0$ on an open neighborhood of x.

For any subset $U\subset X$ , we say f is vanishing properly on U up to order $2$ if it does so for every point in U. Given such a function, we also set

(5.76) $$ \begin{align} \kappa(K):=\max_{x\in K}\mathrm{ord}_{x}(f)\in\{0,1,2\}. \end{align} $$

The following proposition is an extension of Proposition 5.12.

Proposition 5.24. Assume that $(X,J,\Theta )$ , $(L,h_{L})$ , f are given in either cases (I) or (II) with f a nontrivial function in $\mathcal {Q}_{\geq p_{0}}(X,L;\mathbb {R})$ and with prequantum condition $\Theta =c_{1}(L,h_{L})$ . Let U be an open subset of $\operatorname {\mathrm {supp}} f$ such that f vanishes properly on U up to second order. Then for any compact subset K of U, there exists a constant $c_{K}>0$ and $p_{K}\geq p_{0}$ such that for $x\in K$ , $p\geq p_{K}$ ,

(5.77) $$ \begin{align} T^{2}_{f,p}(x,x)\geq c_{K}p^{n-\kappa(K)}. \end{align} $$

Moreover, $\log f^{2}$ is locally integrable on U, and we have the weak convergence of currents on U,

(5.78) $$ \begin{align} \partial\overline{\partial} \log T^{2}_{f,p}(x,x)\rightarrow \partial\overline{\partial} \log f^{2}\,, \quad \text{as } p\rightarrow \infty. \end{align} $$

Around a point x with $f(x)\neq 0$ , the convergence in (5.78) holds in any local $\mathscr {C}^{\ell }$ -norms.

Proof. We start with proving estimate (5.77). For $x_{0}\in U$ , if $f(x_{0})\neq 0$ , then $f^{2}(x_{0})>0$ , so estimate (5.77) holds near $x_{0}$ . If $f(x_{0})=0$ , $df(x_{0})\neq 0$ , then in a sufficiently small neighborhood of $x_{0}$ , there is a constant $c_{x_{0}}>0$ such that have

(5.79) $$ \begin{align} \boldsymbol{b}_{1}(f,f)=\frac{1}{8\pi}(\boldsymbol{r}^{X}f-4\Delta f)f+\frac{1}{2\pi}|\partial f|^{2}\geq c_{x_{0}}|df(x_{0})|^{2}_{g^{T^{\ast}X}_{x_{0}}}>0 \end{align} $$

so that near $x_{0}$ ,

(5.80) $$ \begin{align} T^{2}_{f,p}(x,x)\geq \frac{1}{2}c_{x_{0}}\,p^{n-1}. \end{align} $$

If $\mathrm {ord}_{x_{0}}(f)=2$ , we can adapt the proof of Proposition 5.12. The condition that $\Delta f(x_{0})$ is nonzero with $f\Delta f \leq 0$ near $x_{0}$ implies that on a small neighborhood of $x_{0}$ ,

(5.81) $$ \begin{align} (\boldsymbol{r}^{X}f-4\Delta f)f\geq 0,\; \mu(f,x_{0})>0. \end{align} $$

Then estimate (5.77) still holds near $x_{0}$ . The second part of the proposition follows from arguments analogous to those used to prove Proposition 5.12.

For $f\in \mathcal {Q}_{\geq p_{0}}(X,L;\mathbb {R})$ , the operator $T_{f,p}$ might not be injective so that, in Remark 4.17, we introduce a closed subspace $H^{0}_{(2)}(X,L^{p},f)=(\ker T_{f,p})^{\perp }$ of $H^{0}_{(2)}(X,L^{p})$ and the Gaussian probability measure $\mathbb {P}^{\sharp }_{f,p}$ on it. Consider the sequence of random $\mathcal {L}^{2}$ -holomorphic sections detected by f

(5.82) $$ \begin{align} (\boldsymbol{S}^{\sharp}_{f,p})_{p\geq p_{0}}\in \Pi_{p\geq p_{0}}\big(H^{0}_{(2)}(X,L^{p},f), \; \mathbb{P}^{\sharp}_{f,p}\big), \end{align} $$

and a sequence of random holomorphic sections defined as in the decomposition (4.41),

(5.83) $$ \begin{align} \widehat{\boldsymbol{S}}_{f,p}=\psi_{\ker T_{f,p}}+\boldsymbol{S}^{\sharp}_{f,p},\;\; p\geq p_{0}. \end{align} $$

Recall that the $\mathcal {L}^{2}$ -integrability on X of $\widehat {\boldsymbol {S}}_{f,p}$ depends on the finiteness of $\dim \ker T_{f,p}$ .

From the inequality (4.43) and by Proposition 5.24, we get:

Theorem 5.25. Under the assumptions of Proposition 5.24, let U be an open subset of $\operatorname {\mathrm {supp}} f$ such that f vanishes properly on U up to order $2$ . Then the following assertions hold.

(i) We have the convergence of $(1,1)$ -currents on U as $p\rightarrow \infty $ ,

(5.84) $$ \begin{align} \begin{aligned} &\frac{1}{p}\mathbb{E}^{\mathbb{P}^{\sharp}_{f,p}}[[\operatorname{\mathrm{Div}}(\boldsymbol{S}^{\sharp}_{f,p})]|_{U}] \rightarrow c_{1}(L,h_{L})|_{U},\\ &\frac{1}{p}\mathbb{E}[[\operatorname{\mathrm{Div}}(\widehat{\boldsymbol{S}}_{f,p})]|_{U}] \rightarrow c_{1}(L,h_{L})|_{U}. \end{aligned} \end{align} $$

(ii) For any $\varphi \in \Omega ^{n-1,n-1}_{0}(U)$ , we have

(5.85) $$ \begin{align} \begin{aligned} &\mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}\langle[\operatorname{\mathrm{Div}}(\boldsymbol{S}^{\sharp}_{f,p})],\varphi\rangle= \langle c_{1}(L,h_{L}),\varphi\rangle\right)\\ &=\mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}\langle[\operatorname{\mathrm{Div}}(\widehat{\boldsymbol{S}}_{f,p})],\varphi\rangle= \langle c_{1}(L,h_{L}),\varphi\rangle\right)=1. \end{aligned} \end{align} $$

(iii) If there is a smooth closed positive $(1,1)$ -form $\alpha $ on X with $\int _{X}c_{1}(L,h_{L})\wedge \alpha ^{n-1}<\infty $ such that $\alpha $ is strictly positive on $\overline {U}$ , then we have

(5.86) $$ \begin{align} \begin{aligned} &\mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}[\operatorname{\mathrm{Div}}(\boldsymbol{S}^{\sharp}_{f,p})]|_{U}= c_{1}(L,h_{L})|_{U} \right)=\mathbb{P}\left(\lim_{p\rightarrow \infty}\frac{1}{p}[\operatorname{\mathrm{Div}}(\widehat{\boldsymbol{S}}_{f,p})]|_{U}= c_{1}(L,h_{L})|_{U} \right)=1, \end{aligned} \end{align} $$

where the limit is taken with respect to the weak convergence of $(1,1)$ -currents on U.

Remark 5.26. If X is compact, then $H^{0}_{(2)}(X,L^{p})=H^{0}(X,L^{p})$ , $p\in \mathbb {N}$ , are finite dimensional, and we can take f to be any real smooth function vanishing properly up to order $2$ in the above theorem. If $\kappa (X)\leq 1$ , then the first limit (5.84) is already proved by Ancona and Le Floch [Reference Ancona and Le Floch1]. As mentioned in Subsection 5.3, Ancona and Le Floch studied the fluctuations of the random zeros near a vanishing point of f with order $1$ . Their computations, as well as ours in the proof of Proposition 5.18 for points with vanishing order $2$ , were done by rescaling the expansion (5.8) of $T^2_{f,p}(x,x)$ . Since the problem has local character, their results [Reference Ancona and Le Floch1, Theorem 1.5 and Corollary 1.6] also hold in the noncompact setting for a vanishing point of f of order $1$ .

Acknowledgments

We gratefully acknowledge support of DFG Priority Program 2265 ‘Random Geometric Systems.’ The authors thank Prof. Xiaonan Ma for useful discussions and the referee for carefully reading the paper.

Competing interest

On behalf of all authors, the corresponding author states that there is no competing interest.

Data sharing

Data sharing not applicable to this article as no datasets were generated or analysed during the current study.

Footnotes

G. M. is partially supported by DFG funded projects SFB/TRR 191 (Project-ID 281071066-TRR 191), and the ANR-DFG project QuaSiDy (Project-ID 490843120).

The authors are partially supported by the DFG Priority Program 2265 ‘Random Geometric Systems’ (Project-ID 422743078)

References

Ancona, M. and Le Floch, Y., Berezin–Toeplitz operators, Kodaira maps, and random sections. arXiv:2206.15112 (2022).Google Scholar
Atiyah, M. F., Elliptic operators, discrete groups and von Neumann algebras. in Colloque ‘Analyse et Topologie’ en l’Honneur de Henri Cartan (Orsay, 1974), Astérisque, No. 32-33. (Paris: Soc. Math. France, Paris, 1976), 4372.Google Scholar
Barron, T., Ma, X., Marinescu, G. and Pinsonnault, M., Semi-classical properties of Berezin-Toeplitz operators with ${\mathscr{C}}^k$ -symbol, J. Math. Phys. 55 (2014), 042108, 25.CrossRefGoogle Scholar
Bayraktar, T., Asymptotic normality of linear statistics of zeros of random polynomials. Proc. Amer. Math. Soc. 145 (2017), 29172929.CrossRefGoogle Scholar
Bayraktar, T., Coman, D., Herrman, H. and Marinescu, G., A survey on zeros of random holomorphic sections. Dolomites Res. Notes Approx. 11 (2018), 119.Google Scholar
Bayraktar, T., Coman, D. and Marinescu, G., Universality results for zeros of random holomorphic sections. Trans. Amer. Math. Soc. 373 (2020), 37653791.CrossRefGoogle Scholar
Berndtsson, B., Bergman kernels related to Hermitian line bundles over compact complex manifolds. In Explorations in Complex and Riemannian Geometry, Contemp. Math. vol. 332. (Amer. Math. Soc., Providence, RI, 2003), 117.CrossRefGoogle Scholar
Bleher, P., Shiffman, B. and Zelditch, S., Poincaré–Lelong approach to universality and scaling of correlations between zeros. Comm. Math. Phys. 208 (2000), 771785.CrossRefGoogle Scholar
Błocki, Z. and Pflug, P., Hyperconvexity and Bergman completeness. Nagoya Math. J. 151 (1998), 221225.CrossRefGoogle Scholar
Bloom, T. and Levenberg, N., Random polynomials and pluripotential-theoretic extremal functions. Potential Anal. 42 (2015), 311334.CrossRefGoogle Scholar
Bogomolny, E., Bohigas, O. and Leboeuf, P., Quantum chaotic dynamics and random polynomials. Journal of Statistical Physics 85 (1996), 639679.CrossRefGoogle Scholar
Bordemann, M., Meinrenken, E. and Schlichenmaier, M., Toeplitz quantization of Kähler manifolds and $N\to \infty$ limits. Comm. Math. Phys. 165 (1994), 281296.CrossRefGoogle Scholar
Coman, D., Ma, X. and Marinescu, G., Equidistribution for sequences of line bundles on normal Kähler spaces. Geom. Topol. 21 (2017), 923962.CrossRefGoogle Scholar
Coman, D. and Marinescu, G., Convergence of Fubini-Study currents for orbifold line bundles. Internat. J. Math. 24 (2013), 1350051, 27.CrossRefGoogle Scholar
Coman, D. and Marinescu, G., Equidistribution results for singular metrics on line bundles. Ann. Sci. Éc. Norm. Supér. (4) 48 (2015), 497536.CrossRefGoogle Scholar
Coman, D., Marinescu, G. and Nguyên, V. A., Hölder singular metrics on big line bundles and equidistribution. Int. Math. Res. Not. IMRN, (2016), 50485075.CrossRefGoogle Scholar
Coman, D., Marinescu, G. and Nguyên, V. A.. Approximation and equidistribution results for pseudo-effective line bundles. J. Math. Pures Appl. (9) 115 (2018), 218236.Google Scholar
de Rham, G., Variétés différentiables. Formes, courants, formes harmoniques, Publications de l’Institut de Mathématique de l’Université de Nancago, III. (Hermann, Paris, 1973). Troisième édition revue et augmentée, Actualités Scientifiques et Industrielles [Current Scientific and Industrial Topics], No. 1222b.Google Scholar
Deleporte, A., Low-energy spectrum of Toeplitz operators: The case of wells. J. Spectr. Theory 9 (2019), 79125.CrossRefGoogle Scholar
Deleporte, A., Low-energy spectrum of Toeplitz operators with a miniwell. Comm. Math. Phys. 378 (2020), 15871647.CrossRefGoogle Scholar
Demailly, J. P., Singular Hermitian metrics on positive line bundles. In Complex Algebraic Varieties (Bayreuth, 1990), Lecture Notes in Math., vol. 1507. (Springer, Berlin, 1992) 87104.CrossRefGoogle Scholar
Demailly, J. P., Complex Analytic and Differential Geometry (2012). Available at https://www-fourier.ujf-grenoble.fr/~demailly/documents.html.Google Scholar
Diederich, K., Das Randverhalten der Bergmanschen Kernfunktion und Metrik in streng pseudo-konvexen Gebieten. Math. Ann. 187 (1970), 936.CrossRefGoogle Scholar
Dinh, T. C., Ma, X. and Marinescu, G., Equidistribution and convergence speed for zeros of holomorphic sections of singular Hermitian line bundles. J. Funct. Anal. 271 (2016), 30823110.CrossRefGoogle Scholar
Dinh, T. C., Marinescu, G. and Schmidt, V., Equidistribution of zeros of holomorphic sections in the non-compact setting. J. Stat. Phys. 148 (2012), 113136.CrossRefGoogle Scholar
Dinh, T. C. and Sibony, N., Distribution des valeurs de transformations méromorphes et applications. Comment. Math. Helv 81 (2006), 221258.CrossRefGoogle Scholar
Drewitz, A., Liu, B. and Marinescu, G., Large deviations for zeros of holomorphic sections on punctured Riemann surfaces. Michigan Mathematical Journal (2023), 141. Advance publication.CrossRefGoogle Scholar
Drewitz, A., Liu, B. and Marinescu, G., Toeplitz operators and zeros of square-integrable random holomorphic sections. arXiv:2404.15983 (2024).Google Scholar
Edelman, A. and Kostlan, E.. How many zeros of a random polynomial are real? Bull. Amer. Math. Soc. (N.S.) 32 (1995), 137.CrossRefGoogle Scholar
Edelman, A. and Kostlan, E.. Erratum: ‘How many zeros of a random polynomial are real?’ [Bull. Amer. Math. Soc. (N.S.) 32(1) (1996), 137, 1995; MR1290398 (95m:60082)]. Bull. Amer. Math. Soc. (N.S.) 33, 325.Google Scholar
Federer, H., Geometric Measure Theory, Die Grundlehren der mathematischen Wissenschaften, Band 153. (Springer-Verlag New York, Inc., New York, 1969).Google Scholar
Fefferman, C., The Bergman kernel and biholomorphic mappings of pseudoconvex domains. Invent. Math. 26 (1974), 165.CrossRefGoogle Scholar
Finski, S., Complex embeddings, Toeplitz operators and transitivity of optimal holomorphic extensions. arXiv:2201.04102 (2022).Google Scholar
Grauert, H. and Remmert, R., Theory of Stein Spaces, Classics in Mathematics. (Springer-Verlag, Berlin, 2004).CrossRefGoogle Scholar
Griffiths, P. and Harris, J., Principles of Algebraic Geometry, Wiley Classics Library, (John Wiley & Sons, Inc., New York, 1994).CrossRefGoogle Scholar
Gross, L., Abstract Wiener spaces. In Proc. Fifth Berkeley Sympos. Math. Statist. and Probability (Berkeley, Calif., 1965/66), vol. II: Contributions to Probability Theory, Part 1. (Univ. California Press, Berkeley, CA, 1967), 3142.Google Scholar
Hannay, J. H., Chaotic analytic zero points: exact statistics for those of a random spin state. J. Phys. A 29 (1996), L101L105.CrossRefGoogle Scholar
Hedenmalm, H., Korenblum, B. and Zhu, K., Theory of Bergman Spaces, Graduate Texts in Mathematics, vol. 199. (Springer-Verlag, New York, 2000).CrossRefGoogle Scholar
Herbort, G., The Bergman metric on hyperconvex domains. Math. Z. 232 (1999), 183196.CrossRefGoogle Scholar
Hough, J. B., Krishnapur, M., Peres, Y. and Virág, B., Zeros of Gaussian Analytic Functions and Determinantal Point Processes, University Lecture Series, vol. 51. (American Mathematical Society, Providence, RI, 2009).Google Scholar
Ioos, L., Lu, W., Ma, X. and Marinescu, G., Berezin–Toeplitz quantization of non-compact manifolds. In preparation.Google Scholar
Janson, S., Gaussian Hilbert Spaces, Cambridge Tracts in Mathematics, vol. 29. (Cambridge University Press, Cambridge, 1997).CrossRefGoogle Scholar
Jarnicki, M. and Pflug, P., Invariant Distances and Metrics in Complex Analysis, De Gruyter Expositions in Mathematics, vol. 9. (Walter de Gruyter GmbH & Co. KG, Berlin, 2013).Google Scholar
Kabluchko, Z. and Zaporozhets, D.. Asymptotic distribution of complex zeros of random analytic functions. Ann. Probab. 42 (2014), 13741395.CrossRefGoogle Scholar
Kac, M., On the average number of real roots of a random algebraic equation. Bull. Amer. Math. Soc. 49 (1943), 314320.CrossRefGoogle Scholar
Kahane, J. P., Some Random Series of Functions , 2nd edn., Cambridge Studies in Advanced Mathematics, vol. 5. (Cambridge University Press, Cambridge, 1985).Google Scholar
Klenke, A., Probability Theory, (Springer, Universitext., London, 2014), German edition. Google Scholar
Kobayashi, S., Geometry of bounded domains. Trans. Amer. Math. Soc. 92 (1959), 267290.CrossRefGoogle Scholar
Kuo, H. H., Gaussian Measures in Banach Spaces, Lecture Notes in Mathematics, Vol. 463. (Springer-Verlag, Berlin-New York, 1975).CrossRefGoogle Scholar
Ledoux, M., The Concentration of Measure Phenomenon, Mathematical Surveys and Monographs, vol. 89. (American Mathematical Society, Providence, RI, 2001).Google Scholar
Lindholm, N., Sampling in weighted ${L}^p$ spaces of entire functions in ${\mathbb{C}}^n$ and estimates of the Bergman kernel. J. Funct. Anal. 182 (2001), 390426.CrossRefGoogle Scholar
Littlewood, J. E. and Offord, A. C., On the distribution of the zeros and $\mathrm{a}$ -values of a random integral function. I. J. London Math. Soc. 20 (1945), 130136.CrossRefGoogle Scholar
Littlewood, J. E. and Offord, A. C., On the distribution of zeros and $\mathrm{a}$ -values of a random integral function. II. Ann. of Math. (2) 49 (1948), 885952; errata 50, 990–991 (1949).CrossRefGoogle Scholar
Liu, B. and Zielinski, D., Semipositive line bundles on punctured Riemann surfaces: Bergman kernels and random zeros. arXiv:2407.15106 (2024).CrossRefGoogle Scholar
Ma, X., Geometric quantization on Kähler and symplectic manifolds. In Proceedings of the International Congress of Mathematicians, vol. II. (Hindustan Book Agency, New Delhi, 2010), 785810.Google Scholar
Ma, X. and Marinescu, G., Holomorphic Morse Inequalities and Bergman Kernels, Progress in Mathematics, vol. 254. (Birkhäuser Verlag, Basel, 2007).Google Scholar
Ma, X. and Marinescu, G., Toeplitz operators on symplectic manifolds . J. Geom. Anal. 18 (2008), 565611.CrossRefGoogle Scholar
Ma, X. and Marinescu, G., Berezin–Toeplitz quantization and its kernel expansion. In Geometry and Quantization, Trav. Math., vol. 19. (Univ. Luxemb., Luxembourg, 2011) 125166.Google Scholar
Ma, X. and Marinescu, G., Berezin–Toeplitz quantization on Kähler manifolds. J. Reine Angew. Math. 662 (2012), 156.Google Scholar
Ma, X. and Marinescu, G., Exponential estimate for the asymptotics of Bergman kernels. Math. Ann. 362 (2015), 13271347.CrossRefGoogle Scholar
Marinescu, G. and Vu, D. V., Bergman kernel functions associated to measures supported on totally real submanifolds. J. Reine Angew. Math. 810 (2024), 217251.Google Scholar
Nonnenmacher, S. and Voros, A.. Chaotic eigenfunctions in phase space. J. Statist. Phys. 92 (1998), 431518.CrossRefGoogle Scholar
Offord, A. C., The distribution of the values of an entire function whose coefficients are independent random variables. Proc. London Math. Soc. (3) 14a (1965), 199238.CrossRefGoogle Scholar
Offord, A. C., The distribution of zeros of power series whose coefficients are independent random variables. Indian J. Math. 9 (1967), 175196.Google Scholar
Offord, A. C., The distribution of the values of an entire function whose coefficients are independent random variables. II. Math. Proc. Cambridge Philos. Soc. 118 (1995), 527542.CrossRefGoogle Scholar
Ohsawa, T., A remark on the completeness of the Bergman metric. Proc. Japan Acad. Ser. A Math. Sci. 57 (1981), 238240.CrossRefGoogle Scholar
Sheffield, S., Gaussian free fields for mathematicians. Probab. Theory Related Fields 139 (2007), 521541.CrossRefGoogle Scholar
Shiffman, B. and Zelditch, S., Distribution of zeros of random and quantum chaotic sections of positive line bundles. Comm. Math. Phys. 200 (1999), 661683.CrossRefGoogle Scholar
Shiffman, B. and Zelditch, S., Number variance of random zeros on complex manifolds, II: smooth statistics. Pure Appl. Math. Q. 6 (2010), 11451167.CrossRefGoogle Scholar
Shiffman, B., Zelditch, S. and Zrebiec, S., Overcrowding and hole probabilities for random zeros on complex manifolds. Indiana Univ. Math. J. 57 (2008), 19771997.CrossRefGoogle Scholar
Sodin, M., Zeros of Gaussian analytic functions. Math. Res. Lett. 7 (2000), 371381.CrossRefGoogle Scholar
Sodin, M. and Tsirelson, B., Random complex zeroes. I. Asymptotic normality. Israel J. Math. 144 (2004), 125149.CrossRefGoogle Scholar
Sodin, M. and Tsirelson, B., Random complex zeroes. III. Decay of the hole probability. Israel J. Math. 147 (2005) 371379.CrossRefGoogle Scholar
Sodin, M. and Tsirelson, B., Random complex zeroes. II. Perturbed lattice. Israel J. Math. 152 (2006), 105124.CrossRefGoogle Scholar
Tao, T. and Vu, V., Local universality of zeroes of random polynomials . Int. Math. Res. Not. IMRN (2015), 50535139.CrossRefGoogle Scholar
Tian, G., On a set of polarized Kähler metrics on algebraic manifolds. J. Differential Geom. 32 (1990), 99130.CrossRefGoogle Scholar
Varadhan, S. R. S., Probability Theory, Courant Lecture Notes in Mathematics, vol. 7. New York: New York University, Courant Institute of Mathematical Sciences; (American Mathematical Society, Providence, RI, 2001).Google Scholar
Weil, A., Introduction à l’étude des variétés kählériennes, Publications de l’Institut de Mathématique de l’Université de Nancago, VI. Actualités Sci. Ind. no. 1267. (Hermann, Paris, Paris, 1958).Google Scholar
Zelditch, S., Quantum ergodicity on the sphere. Comm. Math. Phys. 146 (1992), 6171.CrossRefGoogle Scholar
Zelditch, S., Szegö kernels and a theorem of Tian. Internat. Math. Res. Notices (1998) 317331.CrossRefGoogle Scholar
Zelditch, S., From random polynomials to symplectic geometry. In XIIIth International Congress on Mathematical Physics (London, 2000). (Int. Press, Boston, MA, 2001), 367376.Google Scholar
Zrebiec, S., The zeros of flat Gaussian random holomorphic functions on ${\mathbb{C}}^{\mathrm{n}}$ , and hole probability. Michigan Math. J. 55 (2007), 269284.CrossRefGoogle Scholar