Off-spectral analysis of Bergman kernels

The asymptotic analysis of Bergman kernels with respect to exponentially varying measures near emergent interfaces has attracted recent attention. Such interfaces typically occur when the associated limiting Bergman density function vanishes on a portion of the plane, the off-spectral region. This type of behaviour is observed when the metric is negatively curved somewhere, or when we study partial Bergman kernels in the context of positively curved metrics. In this work, we cover these two situations in a unified way, for exponentially varying planar measures on the complex plane. We obtain uniform asymptotic expansions of root functions, which are essentially normalized partial Bergman kernels at an off-spectral point, valid in the entire off-spectral component and protruding into the spectrum as well, which allows us to show error function transition behaviour of the original kernel along the interface. In contrast, previous work on asymptotic expansions of Bergman kernels is typically local, and valid only in the bulk region of the spectrum.

1. Introduction 1.1. Bergman kernels and emergent interfaces. This article is a companion to our recent work [13] on the structure of planar orthogonal polynomials. We will make frequent use of methods developed there, and recommend that the reader keep that article available for ease of reference.
The study of Bergman kernel asymptotics has by now a sizeable literature. The majority of the contributions have the flavor of local asymptotics near a given point w 0 , under a positive curvature condition. However, in the study of partial Bergman kernels for the subspace of all functions vanishing to a given order at the point w 0 , the assumption of vanishing has the effect of introducing a negative point mass for the curvature form at w 0 . In addition to the negative curvature which comes from considering partial Bergman kernels defined by vanishing, we allow for the direct effect of patches of negatively curved geometry. Around the set of negative curvature, a forbidden region (or off-spectral set) emerges. This forbidden region is typically larger than the set of actual negative curvature, and may consist of several connectivity components. Recently, the asymptotic behavior of Bergman kernels near the interface at the edge of the forbidden region has attracted considerable attention. In this work, we intend to investigate this in the fairly general setting of exponentially varying weights in the complex plane. The restriction of the Bergman kernel to the diagonal gives us the density of states, which drops steeply at the interface. Indeed, in the forbidden region the density of states vanishes asymptotically, with exponential decay. One of our main results is that the density of states across the interface converges to the error function in a blow-up, provided the interface is smooth.
The key to obtaining the above-mentioned result is in fact our main result. It concerns the expansion of the coherent state k n (z, w 0 ) of depth n at a given off-spectral point w 0 . This is the renormalized reproducing kernel function at the point w 0 for the Bergman space defined by vanishing to order n at w 0 . When n = 0, the coherent state of depth 0 is just the normalized Bergman kernel K(w 0 , w 0 ) − 1 2 K(z, w 0 ). An important feature of the asymptotic expansion is that z and w 0 are allowed to be macroscopically separated. Such truly off-diagonal expansions have, to the best of our knowledge, not appeared elsewhere. It is easy to see that the coherent state of order n = 0, 1, 2, . . ., at w 0 form an orthonormal basis for the Bergman space, which gives an expansion of the density of states, which leads to the error function asymptotics. Similarly, if we want to handle partial Bergman kernels given by vanishing to order n 0 at w 0 , we expand in the basis given by the coherent states of depth n = n 0 , n 0 + 1, . . ..

1.2.
Coherent states and elementary potential theory. We now introduce the objects of study. The Bergman space A 2 mQ is defined as the collection of all entire functions f in with finite weighted L 2 -norm where dA denotes the planar area element normalized so that the unit disk D has unit area, and where Q is a potential with certain growth and regularity properties (see Definition 1.6.1). We denote the reproducing kernel for A 2 mQ by K m , and for a given point w 0 ∈ C, we consider the coherent state (normalized Bergman kernel) (1.2.1) k m,w0 (z) := K m (w 0 , w 0 ) − 1 2 K m (z, w 0 ), which has norm 1 in A 2 mQ . There is a notion of the spectrum S, also called the spectral droplet. This is the closed set defined in terms of the following obstacle problem. Let SH(C) denote the cone of all subharmonic functions on the plane C, and consider the function Q(z) := sup q(z) : q ∈ SH(C), and q ≤ Q on C .
Whenever Q is C 1,1 -smooth and has some modest growth at infinity, it is known thatQ ∈ C 1,1 as well, and it is a matter of definition thatQ ≤ Q pointwise (see, e.g., [10]). Here, C 1,1 denotes the standard smoothness class of differentiable functions with Lipschitz continuous first order partial derivatives. We define the spectrum (or the spectral droplet) as the contact set (1.2.2) S := z ∈ C :Q(z) = Q(z) .
The forbidden region is the complement S c = C \ S. We will need these notions in the context of partial Bergman kernels as well. For a non-negative integer n and a point w 0 ∈ C, we consider the subspace A 2 mQ,n,w0 of A 2 mQ , consisting of those functions that vanish to order at least n at w 0 . It may happen for some n that this space is trivial, for instance when the potential Q has logarithmic growth only, because then the space A 2 mQ consists of polynomials of a bounded degree. We denote its reproducing kernel by K m,n,w0 , and observe that K m,0,w0 = K m . We shall need also the coherent state of depth n at w 0 (or the root function of order n), denoted k m,n,w0 , which is the unique solution to the optimization problem max Re f (n) (w 0 ) : f ∈ A 2 mQ,n,w0 , f mQ ≤ 1 , provided the maximum is positive, in which case the optimizer has norm f mQ = 1.
In the remaining case, the maximum equals 0, and either only f = 0 is possible, or there are several competing optimizers, simply because we may multiply the function by unimodular constants and obtain alternative optimizers. In both remaining instances we declare that k m,n,w0 = 0. When nontrivial, the root function (or coherent state) of order n at w 0 is connected with the reproducing kernel K m,n,w0 : k m,n,w0 (z) = lim ζ→w0 K m,n,w0 (ζ, ζ) −1/2 K m,n,w0 (z, ζ), where the point ζ should approach w 0 not arbitrarily but in a fashion such that the limit exists and has positive n-th derivative at w 0 . The root function k m,n,w0 will play a key role in our analysis, similar to that of the orthogonal polynomials in the context of polynomial Bergman kernels. The root functions k m,n,w0 all have norm equal to 1 in A 2 mQ , except when they are trivial and have norm 0. As a result of the relation (1.2.3), we may alternatively call the root function k m,n,w0 a normalized partial Bergman kernel. The spectral droplet associated to a family of partial Bergman kernels of the above type is defined in Subsection 2.1 in terms of an obstacle problem, and we briefly outline how this is done. For 0 ≤ τ < +∞, let SH τ,w0 (C) denote the convex set SH τ,w0 (C) = q ∈ SH (C) : q(z) ≤ τ log|z − w 0 | + O(1) as z → w 0 , so that for τ = 0 we recover SH(C). We consider the corresponding obstacle problem (1.2.4)Q τ,w0 (z) = sup q(z) : q ∈ SH τ,w0 (C), q ≤ Q on C , and observe that τ →Q τ,w0 is monotonically decreasing pointwise. We define a family of spectral droplets as the coincidence sets (1.2.5) S τ,w0 = z ∈ C : Q(z) =Q τ,w0 (z) .
Due to the monotonicity, the droplets S τ,w0 get smaller as τ increases, starting from S 0,w0 = S for τ = 0. The partial Bergman density ρ m,n,w0 (z) := m −1 K m,n,w0 (z, z) e −2mQ(z) , z ∈ C, may be viewed as the normalized local dimension of the space A 2 mQ,n,w0 (C), and, in addition, it has the interpretation as the intensity of a corresponding (possibly infinite) Coulomb gas. In the case n = 0 we omit the word "partial" and speak of the Bergman density. It is known that in the limit as m, n → +∞ with n = mτ , in the sense of convergence of distributions. In particular, ∆Q ≥ 0 holds a.e. on S. Here, we write ∆ for differential operator ∂∂, which is one quarter of the usual Laplacian. The above convergence reinforces our understanding of the droplets S τ,w0 as spectra, in the sense that the Coulomb gas may be thought to model eigenvalues (at least in the finite-dimensional case). The bulk of the spectral droplet S τ,w0 is the set {z ∈ int(S τ,w0 ) : ∆Q(z) > 0}, where "int" stands for the operation of taking the interior.
1.3. Further background on Bergman kernel expansions. Our motivation for the above setup originates with the theory of random matrices, specifically the random normal matrix ensembles. We should mention that an analogous situation occurs in the study of complex manifolds. The Bergman kernel then appears in the study of spaces of L 2 -integrable global holomorphic sections of L m , where L m is a high tensor power of a holomorphic line bundle L over the manifold, endowed with an hermitian fiber metric h. If {U i } i is a coordinate system on a manifold, then a holomorphic section s to a line bundle L can be written in local coordinates as s = s i e i , where e i are local basis elements for L and s i are holomorphic functions. The pointwise norm of a section s of L m may be then be written as |s| 2 h m = |s i | 2 e −mφi on U i , for some smooth real-valued functions φ i . Along with a volume form on the base manifold, this defines an L 2 space which shares many characteristics with the spaces considered here.
The asymptotic behavior of Bergman kernels has been the subject of intense investigation. However, the understanding has largely been limited to the analysis of the kernel inside the bulk of the spectrum, in which case the kernel enjoys a full local asymptotic expansion. The pioneering work on Bergman kernel asymptotics begins with the efforts by Hörmander [14] and Fefferman [8]. Developing further the microlocal approach of Hörmander, Boutet de Monvel and Sjöstrand [5] obtain a near-diagonal expansion of the Bergman kernel close to the boundary of the given domain. Later, in the context of Kähler geometry, the influential peak section method was introduced by Tian [20]. His results were refined further by Catlin and Zelditch [6,21], while the connection with microlocal analysis was greatly simplified in the more recent work by Berman, Berndtsson, and Sjöstrand [4]. A key element of all these methods is that the kernel is determined by the local geometry around the given point. This property is absent when we consider the kernel near an off-spectral point or near a boundary point of the spectral droplet.
In the recent work [13], we analyze the boundary behavior of polynomial Bergman kernels, for which the corresponding spectral droplet is compact, connected, and has a smooth Jordan curve as boundary. The analysis takes the path via a full asymptotic expansion of the orthogonal polynomials, valid off a sequence of increasing compacts which eventually fill the droplet. By expanding the polynomial kernel in the orthonormal basis provided by the orthogonal polynomials, the error function asymptotics emerges along smooth spectral boundaries.
The appearance of an interface for partial Bergman kernels in higher dimensional settings and in the context of complex manifolds has been observed more than once, notably in the work by Shiffman and Zelditch [19] and by Pokorny and Singer [15]. That the error function governs the transition behavior across the interface was observed later in several contexts. For instance, in [16], Ross and Singer investigate the partial Bergman kernels associated to spaces of holomorphic sections vanishing along a divisor, and obtain error function transition behavior under the assumption that the set-up is invariant under a holomorphic S 1 -action. This result was later extended by Zelditch and Zhou [22], in the context of S 1 -symmetry. More recently, Zelditch and Zhou [23] obtain the same transition for so-called spectral partial Bergman kernels, defined in terms of the Toeplitz quantization of a general smooth Hamiltonian.
Our methods for obtaining asymptotic expansions of coherent states centered at off-spectral points do not easily extend to the higher complex dimensional setting. Illustration of the spectral droplet corresponding to the potential Q(z) = |z| 2 − log(a + |z| 2 ), with a = 0.04. The spectrum is illustrated with a thick line, and appears as the contact set between Q (solid) and the solutionQ to the obstacle function (dashed).
It appears plausible, however, that our analysis would apply when the complex plane C is replaced by a compact Riemann surface with one or possibly several punctures.

1.4.
Off-spectral and off-diagonal asymptotics of coherent states. The main contribution of the present work, put in the planar context, is a non-local asymptotic expansion of coherent states centered at an off-spectral point. For ease of exposition, we begin with a version that requires as few prerequisites as possible for the formulation. We denote by Q(z) an admissible potential, by which we mean the following: (i) Q : C → R is C 2 -smooth, and has sufficient growth at infinity: (ii) Q is real-analytically smooth and strictly subharmonic in a neighborhood of ∂S, where S is the contact set of (1.2.2), (iii) there exists a bounded component Ω of the complement S c = C \ S which is simply connected, and has real-analytically smooth Jordan curve boundary. We consider the case when there exists a non-trivial off-spectral component Ω which is bounded and simply connected, with real-analytic boundary, and pick a "root point" w 0 ∈ Ω. To be precise, by an off-spectral component we mean a connectivity component of the complement S c . This situation occurs, e.g., if the potential is strictly superharmonic in a portion of the plane, as is illustrated in Figure 1.1. In terms of the metric, this means that there is a region where the curvature is negative.
A word on notation. To formulate our first main result, we need the function Q w0 , which is bounded and holomorphic in the off-spectral component Ω and whose real part equals Q along the boundary ∂Ω. To fix the imaginary part, we require that Q w0 (w 0 ) ∈ R. In addition, we need the conformal mapping ϕ w0 which takes Ω onto the unit disk D with ϕ w0 (w 0 ) = 0 and ϕ w0 (w 0 ) > 0. Since the boundary ∂Ω is assumed to be a real-analytically smooth Jordan curve, both the function Q w0 and the conformal mapping ϕ w0 extend analytically across ∂Ω to a fixed larger domain. By possibly considering a smaller fixed larger domain, we may assume that the extended function ϕ w0 is conformal on the larger region. These observations are essential for our first main result, since we want the asymptotics to hold across the interface ∂Ω.
Assuming that Q is an admissible potential, we have the following. Given a positive integer κ and a positive real A, there exist a neighborhood Ω of the closure of Ω and bounded holomorphic functions B j,w0 on Ω for j = 0, . . . , κ, as well as domains Ω m = Ω m,A with Ω ⊂ Ω m ⊂ Ω which meet 1 2 , such that the normalized Bergman kernel at the point w 0 enjoys the asymptotic expansion as m → +∞, where the error term is uniform on Ω m . Here, the main term B 0,w0 is obtained as the unique zero-free holomorphic function on Ω which is smooth up to the boundary with B 0,w0 (w 0 ) > 0, and with prescribed boundary modulus Moreover, if A is big enough, then Ωm |k m,w0 (z)| 2 e −2mQ(z) dA(z) = 1 + O(m −κ ).
As an illustration of this result, we show the Gaussian wave character of the Berezin density k m (w 0 , w 0 ) −1 |k m (z, w 0 )| 2 e −2mQ(z) in Figure 1 Remark 1.4.2. Using an approach based on Laplace's method, the functions B j,w0 may be obtained algorithmically, for j = 1, 2, 3, . . ., see Theorem 3.2.2 below. The details of the algorithm are analogous with the case of the orthogonal polynomials presented in [13].
Commentary to the theorem. The point w 0 may be chosen as any fixed point in Ω, while the point z is allowed to vary anywhere inside the set Ω m , which contains the entire off-spectral component Ω as well as a shrinking neighborhood of the boundary ∂Ω. Hence, we obtain macroscopically off-diagonal asymptotics, where the points z and w 0 are either both off-spectral, or where z is spectral nearboundary and w 0 off-spectral, respectively. To our knowledge, this is the first instance of such asymptotics. Indeed, earlier work covers either diagonal or neardiagonal asymptotics inside the bulk of the spectrum. The analysis of Bergman kernel asymptotics for two macroscopically separated points in the bulk S • appears difficult if we ask for high precision, and the same can be said for the case when w 0 ∈ Ω is off-spectral and z ∈ S • is a bulk point. To address the latter issue, we  1.5. Expansion of partial Bergman kernels in terms of root functions. For a (τ, w 0 )-admissible potential, the partial Bergman kernel K m,n,w0 with the root point w 0 is well-defined and nontrivial. In analogy with Taylor's formula, it enjoys an expansion in terms of the root functions k m,n ,w0 for n ≥ n. Theorem 1.5.1. Under the above assumption of (τ, w 0 )-admissibility of Q, we have that Proof. For n , n ≥ n with n < n , the functions k m,n ,w0 and k m,n ,w0 are orthogonal in A 2 mQ . If one of them is trivial, orthogonality is immediate, while if both are nontrivial, we argue as follows. Let ζ , ζ ∈ C be close to w 0 , and calculate that K m,n ,w0 (ζ , ζ ) − 1 2 K m,n ,w0 (ζ , ζ ) − 1 2 K m,n ,w0 (·, ζ ), K m,n ,w0 (·, ζ ) mQ which tends to 0 as ζ → w 0 . The claimed orthogonality follows. Moreover, since the root functions k m,n ,w0 have unit norm when nontrivial, the expression +∞ n =n k m,n ,w0 (z)k m,n ,w0 (w) equals the reproducing kernel function for the Hilbert space with the norm of A 2 mQ spanned by the vectors k m,n ,w0 with n ≥ n. It remains to check that this is the whole partial Bergman space A 2 mQ,n,w0 . To this end, let f ∈ A 2 mQ,n,w0 be orthogonal to all the the vectors k m,n ,w0 with n ≥ n. By the definition of the space A 2 mQ,n,w0 , this means that where c = 0 is complex. At the same time, the existence of such nontrivial f entails that the corresponding root functions k m,N,w0 is nontrivial as well, and that K m,N,w0 (ζ, ζ) |ζ − w 0 | 2N for ζ near w 0 . On the other hand, the orthogonality between f and k m,N,w0 gives us that where we approach w 0 only in an appropriate direction so that the limit exists. But this contradicts the given asymptotic behavior of f (ζ) near w 0 , since c = 0 tells us that any limit of the right-hand side would be nonzero.
1.6. Off-spectral asymptotics of partial Bergman kernels. Given a point w 0 ∈ C we recall the partial Bergman spaces A 2 mQ,n,w0 , and the associated spectral droplets S τ,w0 (see (1.2.5)), where we keep n = τ m. Before we proceed with the formulation of the second result, let us fix some terminology. Definition 1.6.1. A real-valued potential Q is said to be (τ, w 0 )-admissible if the following conditions hold: (i) Q : C → R is C 2 -smooth and has sufficient growth at infinity: (ii) Q is real-analytically smooth and strictly subharmonic in a neighborhood of ∂S τ,w0 . (iii) The point w 0 is an off-spectral point, i.e., w 0 / ∈ S τ,w0 , and the component Ω τ,w0 of the complement S c τ,w0 containing the point w 0 is bounded and simply connected, with real-analytically smooth Jordan curve boundary. If for an interval I ⊂ [0, +∞[, the potential Q is (τ, w 0 )-admissible for each τ ∈ I and {Ω τ,w0 } τ ∈I is a smooth flow of domains, then Q is said to be (I, w 0 )-admissible.
Generally speaking, off-spectral components may be unbounded. It is for reasons of simplicity that we focus on bounded off-spectral components in the above definition.
A word on notation. We will assume in the sequel that Q is (I, w 0 )-admissible for some non-trivial compact interval I = I 0 . For an illustration of the situation, see Figure 1.3. We let ϕ τ,w0 denote the surjective Riemann mapping which by our smoothness assumption on the boundary ∂Ω τ,w0 extends conformally across ∂Ω τ,w0 . We denote by Q τ,w0 the bounded holomorphic function in Ω τ,w0 whose real part equals Q on ∂Ω τ,w0 and is real-valued at w 0 . It is tacitly assumed to extend holomorphically across the boundary ∂Ω τ,w0 . We now turn to our second main result.
Theorem 1.6.2. Assume that the potential Q is (I 0 , w 0 )-admissible, where the interval I 0 is compact. Given a positive integer κ and a positive real A, there exists a neighborhood Ω τ,w0 of the closure of Ω τ,w0 and bounded holomorphic functions B j,τ,w0 on Ω τ,w0 , as well as domains Ω τ,w0,m = Ω τ,w0,m,A with Ω τ,w0 ⊂ Ω τ,w0,m ⊂ Ω τ,w0 which meet 1 2 , such that the root function of order n at w 0 enjoys the expansion on Ω τ,w0,m as n = τ m → +∞ while τ ∈ I 0 , where the error term is uniform. Here, the main term B 0,τ,w0 is zero-free and smooth up to the boundary on Ω τ,w0 with B 0,τ,w0 (w 0 ) > 0, and with prescribed boundary modulus Moreover, for A big enough, we have as n = τ m → +∞.
Commentary to the theorem. (a) As in Theorem 1.4.1, the point w 0 may be chosen as any fixed point in Ω, while the point z is allowed to vary anywhere inside the set Ω τ,w0,m . Arguing in a fashion analogous to what we did in the commentary following Theorem 1.4.1, we may conclude that k m,n,w0 is small, both pointwise and in the L 2 -sense, away from the set Ω τ,w0,m .
(b) Our Theorems 1.4.1 and 1.6.2 seemingly cover different instances of the root function asymptotics. In fact, in a certain sense, they are equivalent. Indeed, modulo technicalities we may obtain Theorem 1.6.2 from Theorem 1.4.1 by replacing the potential Q(z) withQ τ,w0 (z) = Q(z) − τ log |z − w 0 |. On the other hand, the former theorem is the limit case τ → 0 of the latter.
(c) Theorem 1.6.2 should be compared with Theorem 4.1 in the work of Shiffman and Zelditch [19]. There, an asymptotic expansion of the diagonal restriction K N P (z, z) of the Bergman kernel associated to a family of scaled Newton polytopes N P is obtained, as N → +∞. This expansion is valid deep inside the corresponding forbidden region. In one complex dimension, this is equivalent to studying a certain weighted partial Bergman density for a polynomial Bergman space. Hence there is some overlap with the present work as well as with [13].
(d) The genuinely off-diagonal and off-spectral asymptotic expansion of the coherent states obtained here seem not to have any analogue elsewhere. Although we restrict ourselves to the setting of one complex variable, we compensate for this by making minimal geometric assumptions.

1.7.
Interface transition of the density of states. As a consequence of Theorem 1.6.2, we obtain the transition behavior of the (partial) Bergman densities at emergent interfaces. To explain how this works, we fix a bounded simply connected off-spectral component Ω with real-analytic Jordan curve boundary, associated to either a sequence of partial Bergman kernels K m,n,w0 of the space A 2 mQ,n,w0 with the ratio τ = n m fixed, or, alternatively, a sequence K m of full Bergman kernels. We think of the latter as the parameter choice τ = 0. We assume that the potential Q is real-analytically smooth and strictly subharmonic near ∂Ω. Let z 0 ∈ ∂Ω, and denote by ν ∈ T the inward unit normal to ∂Ω at z 0 . We define the rescaled density where the rescaled variable is defined implicitly by where the convergence is uniform on compact subsets.

1.8.
Comments on the exposition and a guide to the proofs. In Section 2, we explain the proofs of the main results, which are Theorems 1.4.1 and 1.6.2 together with Corollary 1.7.1. Our approach is analogous to that of [13], and we make an effort to explain exactly what needs to be modified for the techniques to apply in the present context. As in [13], the proofs involve the construction of a flow of loops near the spectral boundary, called the orthogonal foliation flow. The flow is constructed in an iterative procedure, which produces both the terms in the asymptotic expansion of the coherent state as well as the successive terms in the expansion of the flow. The flow {γ t } t is constructed to have the following property.
If the normal velocity is denoted by ν, we want the Szegő kernel for the point w 0 with respect to the interior of the curve γ t and the induced weighted arc-length measure e −2mQ νdσ on γ t to be stationary in t. In turn, this allows us to glue together the Szegő kernels to form the Bergman kernel. In order to obtain the coefficient functions of the expansions in a more straightforward fashion, we apply the method developed in [13] which is based on Laplace's method. Another important ingredient, which allows localization to a neighborhood of each off-spectral component, is Hörmander's∂-estimate, suitably modified to the given needs. Since this is the method which allows us to change the geometry drastically, we sometimes refer to it as∂-surgery.
In Section 3, we develop a more general version of the foliation flow lemma, which allows us to introduce a conformal factor in the area form. To avoid unnecessary repetition, both the orthogonal foliation flow and the algorithm for computing the coefficient functions in the asymptotic expansions are explained only in this more general setting. We supply a couple of applications of this extension, including a stability result for the root functions and the orthogonal polynomials under a 1.9. Acknowledgements. We wish to thank Steve Zelditch, Bo Berndtsson, and Robert Berman for their interest in this work. In addition we should thank the anonymous referee for his or her helpful efforts.

2.
Off-spectral expansions of normalized kernels 2.1. A family of obstacle problems and evolution of the spectrum. The spectral droplets (1.2.2) and the partial analogues (1.2.5) were defined earlier. From that point of view, the spectral droplet S is the instance τ = 0 of the partial spectral droplets S τ,w0 . We should like to point out here that the partial spectral droplet S τ,w0 emerges as the full spectrum under a perturbation of the potential Q. To see this, we consider the perturbed potential and observe that the coincidence setS forQ equals the partial spectral droplet S τ,w0 .
The following proposition summarizes some basic properties of the functionQ τ,w0 given by (1.2.4). We refer to [10] for the necessary details.
Proposition 2.1.1. Assume that Q ∈ C 2 (C) is real-valued with the logarithmic growth of condition (i) of Definition 1.6.1. Then for each τ with 0 ≤ τ < τ Q and for each point w 0 ∈ C, the functionQ τ,w0 is a subharmonic function in the plane C which is C 1,1 -smooth off w 0 , and harmonic on C \ (S τ,w0 ∪ {w 0 }). Near the point w 0 we haveQ The evolution of the free boundaries ∂S τ,w0 , which is of fundamental importance for our understanding of the properties of the normalized reproducing kernels, is summarized in the following.
The continuous chain of off-spectral components Ω τ,w0 for τ ∈ I 0 deform according to weighted Laplacian growth with weight 2∆Q, that is, for τ, τ ∈ I 0 with τ < τ , and for any bounded harmonic function h on Ω τ,w0 , we have that Fix a point ζ ∈ ∂Ω τ,w0 , and denote for real ε by ζ ε,w0 the point closest to ζ in the intersection where ν τ (ζ) ∈ T points in the inward normal direction at ζ with respect to Ω τ,w0 . Then we have that and the outer normal n τ −ε,w0 (ζ ε ) satisfies Proof. That the domains deform according to Hele-Shaw flow is a direct consequence of the relation of Ω τ,w0 to the obstacle problem. To see how it follows, assume that h is harmonic on Ω τ,w0 and C 2 -smooth up to the boundary, and apply Green's formula to obtain where the latter integral is understood in the sense of distribution theory. AsQ τ,w0 is a harmonic perturbation of τ times the Green function for Ω τ,w0 , the result follows by writing Ωτ,w 0 \Ω τ ,w 0 h∆QdA as the difference of two integrals of the above form, and by approximation of bounded harmonic functions by harmonic functions C 2 -smooth up to the boundary.
The second part follows along the lines of [13, Lemma 2.3.1].
We turn next to an off-spectral growth bound for weighted holomorphic functions.
Proposition 2.1.3. Assume that Q is admissible and denote by K τ,w0 a closed subset of the interior of S τ,w0 . Then there exist constants c 0 and C 0 such that for any f ∈ A 2 mQ,n,w0 (C) it holds that In case K τ,w0 = ∅, the estimate holds globally.
Proof. This follows immediately by an application of the maximum principle, together with the result of Lemma 2.2.1 in [13], originating from [1].

Some auxilliary functions.
There are a number of functions related to the potential Q that will be useful in the sequel. We denote by Q τ,w0 the bounded holomorphic function on Ω τ,w0 whose real part on the boundary curve ∂Ω τ,w0 equals Q, uniquely determined by the requirement that Im Q τ,w0 (w 0 ) = 0. We also need the functionQ τ,w0 , which denotes the harmonic extension ofQ τ,w0 across the boundary of the off-spectral component Ω τ,w0 . These two functions are connected via Since we work with (τ, w 0 )-admissible potentials Q, the off-spectral component Ω τ,w0 is a bounded simply connected domain with real-analytically smooth Jordan curve boundary. Without loss of generality, we may hence assume that Q τ,w0 ,Q τ,w0 as well as the conformal mapping ϕ τ,w0 extend to a common domain Ω 0 , containing the closureΩ τ,w0 . By possibly shrinking the interval I 0 , we may moreover choose the set Ω 0 to be independent of the parameter τ ∈ I 0 .
2.3. Canonical positioning. An elementary but important observation for the main result of [13] is that we may ignore a compact subset of the interior of the compact spectral droplet S τ associated to polynomial Bergman kernels when we study the asymptotic expansions of the orthogonal polynomials P m,n (with τ = n m ). Indeed, only the behavior in a small neighborhood of the complement S c τ is of interest, and the∂-surgery methods allow us to disregard the rest. The physical intuition behind this is the interpretation of the probability density |P m,n | 2 e −2mQ as the net effect of adding one more particle to the system, and since the positions in the interior of the droplet are already occupied we would expect the net effect to occur near the boundary. The fact that we may restrict our attention to a simply connected proper subset of the Riemann sphereĈ breaks up the rigidity and allows us to apply a conformal mapping to place ourselves in an appropriate model situation.
In the present context, we consider the Riemann mapping ϕ τ,w0 which maps the off-spectral region Ω τ,w0 onto the unit disk D, and has a conformal extension to a neighborhood of Ω τ,w0 . We let the canonical positioning operator for the point w 0 with respect to the off-spectral component Ω τ,w0 be given by and put . An essential property of this operator is that Λ m,n,w0 acts isometrically from L 2 (e −2mRτ,w 0 ) to L 2 (e −2mQ ), wherever it is well-defined.
As for the coherent states, we will analyze them in terms of the canonical positioning operator Λ m,n,w0 . This simplifies the geometry of Ω τ,w0 by mapping it to the unit disk D, and simplifies the weight. Indeed, the function R τ,w0 is flat to order 2 at the unit circle, and consequently the weight e −2mRτ,w 0 behaves like a Gaussian ridge.
We summarize the properties of the operator Λ m,n,w0 in the following proposition. For a potential V and a domain Ω with w 0 ∈ Ω, we denote by A 2 mV,n,w0 (Ω) the space of holomorphic functions on Ω which vanish to order n at w 0 ∈ Ω, endowed with the topology of L 2 (e −2mV , Ω). In case n = 0 we simply denote the space by A 2 mV (Ω).
Proof. The conclusion is immediate by the defining normalizations of the conformal mapping ϕ τ,w0 .
The following definition is an analogue of Definition 3.1.2 in [13]. We denote by Ω 1 a domain containing the closure of the off-spectral component Ω τ,w0 , and let χ 0,τ denote a C ∞ -smooth cut-off function which vanishes off Ω 1 , and equals 1 in a neighborhood of the closure of Ω τ,w0 . (ii) The approximate root functions have norm approximately equal to 1, (iii) The functions F m,n,w0 are approximately real and positive at w 0 , in the sense that the leading coefficient a m,n,w0 = lim z→w0 (z − w 0 ) −n F m,n,w0 (z) satisfies Re a m,n,w0 > 0 and Im a m,n,w0 Re a m,n,w0 = O(m −κ− 1 12 ).
We remark that the exponents in the above error terms are chosen for reasons of convenience, related to the correction scheme of Subsection 2.5.

2.4.
The orthogonal foliation flow. The orthogonal foliation flow {γ m,n,t } t is a smooth flow of closed curves near the unit circle T, originally formulated in [13] in the context of orthogonal polynomials. The defining property is that P m,n should be approximately orthogonal to the lower degree polynomials along the curves Γ m,n,t = φ −1 τ (γ m,n,t ) with respect to the induced measure e −2mQ ν n ds, where ν n denotes the normal velocity of the flow {Γ m,n,t } t and ds denotes normalized arc length measure.
Smoothness classes and polarization of smooth functions. We fix the smoothness class of the weights under consideration, and adapt Definition 4.2.1 in [13] to the the present setting. First, we need the notion of polarization, which applies to real-analytically smooth functions. If R(z) is real-analytic, there exists a function of two complex variables, denoted by R(z, w), which is holomorphic in (z,w) in a neighborhood of the diagonal, with diagonal restriction R(z, z) = R(z). The function R(z, w) is referred to as the polarization of R(z), and it is uniquely determined by its diagonal restriction R(z). If R(z, w) is such a polarization of a function R(z) which is real-analytically smooth near the circle T and quadratically flat there, then R(z) = (1 − |z| 2 ) 2 R 0 (z) and in polarized form while on the annulus A(η −1 , η) := {z ∈ C : η −1 < |z| < η} R is real-analytically smooth and has a polarization R(z, w) which is holomorphic in (z,w) on the 2σ-fattened diagonal annuluŝ A(η, σ) = (z, w) ∈ A(η −1 , η) × A(η −1 , η) : |z − w| ≤ 2σ , and factors as R(z, w) = (1 − zw) 2 R 0 (z, w), where R 0 (z, w) is holomorphic (z,w) on the setÂ(η, σ), and bounded and bounded away from zero there. We say that a subset S ⊂ W(η, σ) is a uniform family, provided that for each R ∈ S, the corresponding R 0 (z, w) is uniformly bounded and bounded away from 0 onÂ(η, σ) while the constant α(R) is uniformly bounded away from 0.
The point with above definition is that it lets us encode uniformity properties of the potentials R τ,w0 with respect to the parameter τ and the point w 0 .
For a polarized function f (z, w), we let f T (z) denote the restriction of f (z, z) for z ∈ T, wherever it is well-defined. We recall from Proposition 4.2.2 of [13] that if f (z, w) is holomorphic in (z,w) on the 2σ-fattened diagonal annulusÂ(η, σ) and if the parameters meet 1 < η ≤ √ 1 + σ 2 + σ, then it follows that the function f T extends holomorphically to the annulus A(η −1 , η). We may need to restrict the numbers η, σ further. Indeed, it turns out that we need that the functions log ∆R, R = √ R (chosen to be positive inside the unit circle and negative outside) as well as log(−z∂ zR ) have polarizations which are holomorphic in (z,w) for (z, w) ∈Â(η, σ) and uniformly bounded there as well. If R belongs to a uniform family of W(η 0 , σ 0 ), then there exist (η 1 , σ 1 ) such that these properties hold for the polarizations with η = η 1 and σ = σ 1 (See Proposition 4.2.3 of [13]), where we moreover require that 1 < η 1 ≤ 1 + σ 2 1 + σ 1 .
Lemma 2.4.2. Let K be a compact subset of each of the domains Ω τ,w0 , where τ ∈ I 0 . Then there exist constants η, σ with η > 1 and σ > 0, such that the collection of weights R τ,w0 with w 0 ⊂ K and τ ∈ I 0 is a uniform family in W(η, σ).
This is completely analogous to the corresponding claim in of [13], which was expressed in the context of an exterior conformal mapping.
The orthogonal foliation flow near the unit circle. The existence of the orthogonal foliation flow around the circle T and the asymptotic expansion of the root functions after canonical positioning are stated in the following lemma (compare with Lemma 4.1.2 in [13]). For the proof, we refer to the sketched proof of Lemma 3.3.2 below, as well as the complete proof of Lemma 4.1.2 in [13], for the case of orthogonal polynomials.
Lemma 2.4.3. Fix an accuracy parameter κ and let R ∈ W(η 0 , σ 0 ). Then, if η 1 is as above, there exist a radius η 2 with 1 < η 2 < η 1 , bounded holomorphic functions f s on D(0, η 1 ) of the form and conformal mappings ψ s,t from D(0, ρ 2 ) into the plane given by s j t lψ j,l such that for s, t small enough, the domains ψ s,t D grow with t, while they remain contained in D(0, η 1 ). Moreover, for ζ ∈ T we have that For small positive s, when t varies in the interval [−β s , β s ] with β s := s 1/2 log 1 s , the flow of loops {ψ s,t (T)} t cover a neighborhood of the circle T of width proportional to β s smoothly. In addition, the first term B 0 is zero-free, positive at the origin, and has modulus |B 0 | = π − 1 4 (∆R) 1 4 on T. The other terms B j are all real-valued at the origin. The implied constant in (2.4.1) is uniformly bounded, provided that R is confined to a uniform family of W(η 0 , σ 0 ).

2.5.∂-corrections and asymptotic expansions of root functions.
In this section, we supply a proof of the main result, Theorem 1.6.2. The proof consists of two parts. First, we construct a family of approximate root function of a given order and accuracy, after which we apply Hörmander-type∂-estimates to correct these approximate kernels to entire functions. The precise result needed for the correction scheme runs as follows.
This is an immediate consequence of Corollary 2.4.2 in [13], and essentially amounts to Hörmander's classical bound for the∂-equation in the given setting.
We turn to the proof of Theorem 1.6.2.
Sketch of proof of Theorem 1.6.2. As the proof is analogous to that of Theorems If we write B j,τ,w0 := (ϕ τ,w0 ) 1 2 B j,τ,w0 • ϕ τ,w0 , it follows that k κ m,n,w0 := m 1 4 Λ m,n,w0 [f κ m,n,w0 ] has the claimed form. It remains to show that k κ m,n,w0 is a family of approximate root functions of order n at w 0 with the stated uniformity property, and to show that it is close to the true normalizing reproducing kernel.
We denote by D m,n,w0 the domain covered by the foliation flow, over the parameter range −δ m ≤ t ≤ δ m , where δ m := m − 1 2 log m. Moreover, we define the domain E τ,w0 as the image of D(0, η ) under ϕ −1 τ,w0 , and let χ 0 = χ 0,τ,w0 denote an an appropriately chosen smooth cut-off function, which takes the value 1 on a neighborhood ofΩ τ,w0 and vanishes off E τ,w0 . If we let χ 1 := χ 0 • ϕ −1 τ,w0 denote the corresponding cut-off extended to vanish off D(0, η ), we may show that which holds for some α 0 > 0 as a consequence of the Gaussian ridge behavior of the function e −2mRτ,w 0 around the unit circle. We now observe that in view of (2.5.2) and (2.5.3), the isometric property of Λ m,n,w0 implies that χ 0 k κ m,n,w0 has mQ,n,w0 be given, and put q = m − 1 4 Λ −1 m,n,w0 [g]. Then, by the isometric property of Λ m,n,w0 and the estimate (2.5.3), it follows that Dm,n,w 0 where we have applied the Cauchy-Schwarz inequality together with the estimate (2.5.3) to obtain the error term. The function f κ m,n,w0 is zero-free up to the boundary in D m,n,w0 provided that m is large enough, as the main term is bounded away from 0 in modulus, and consecutive terms are much smaller. Also, for large enough m, it holds that χ 1 = 1 on D m,n,w0 . We now introduce the function × Re ζ ∂ t ψ m,n,t (ζ)ψ m,n,t (ζ) ds(ζ)dt = 2m where the last step uses the flow equation of Lemma 2.4.3. We now make the crucial observation is that for fixed t, the composition q m,n • ψ m,n,t is holomorphic, so that we may apply the mean value property: Here, in order to obtain the last estimate, we have used (2.5.4) backwards with q m,n replaced by |q m,n |, and the estimate (2.5.2) to obtain Dm,n,w 0 Now if q(0) = 0, that is, if g vanishes to order n + 1 or higher at w 0 , then χ 0 k κ m,n,w0 and g are approximately orthogonal in L 2 mQ . For further details regarding the above computations, we refer to Subsection 4.8 of [13].
The∂-correction scheme. The approximate normalized reproducing kernels are not globally defined, and are consequently not elements of our Bergman spaces of entire functions. However, by applying the Hörmander-type∂-estimate of Proposition 2.5.1, we obtain a solution u = u m,n,w0 to the equation ∂u = k κ m,n,q0∂ χ 0 which exponential decay of the norm of u m,n,w0 in L 2 mQ . By the proposition, it vanishes to order n at the root point w 0 , and has exponentially small norm in L 2 mQ = L 2 (C, e −2mQ ). The function k m,n,w0 := χ 0 k κ m,n,w0 − u m,n,w0 is also an approximate normalized partial Bergman reproducing kernel of the correct accuracy, but this time it at least is an element of the right space, k m,n,w0 ∈ A 2 mQ,n,w0 .
Here, we have the norm estimate (2.5.6) P m,n+1,w0 k m,n,w0 mQ = sup which shows that the correction is very small. By construction, k κ m,n,w0 vanishes precisely to the order n at the root point w 0 . Moreover, the small perturbations u m,n,w0 and P m,n+1,w0 k m,n,w0 vanish at least to order n at w 0 . It follows that k m,n,w0 vanishes precisely to the correct order, that is to say, holds near w 0 for some complex constant C = 0. The constant C is close to being positive real, since the∂-correction u m,n,w0 is small. Indeed, we have C = (1 + O(e −α1m )) C 1 where the constant C 1 > 0 may depend on all the parameters but the parameter α 1 > 0 is a uniform constant. Since the functionk m,n,w0 is automatically orthogonal to A 2 mQ,n+1,w0 , it follows thatk m,n,w0 equals a scalar multiple of the true root function k m,n,w0 : k m,n,w0 = c k m,n,w0 , for some complex constant c = 0. In view of the above, we conclude that c = (1 + O(e −α1m )) c 1 , where c 1 > 0 may depend on all the parameters. Ask m,n,w0 (w 0 ) is approximately real, it follows that c = c γ, where c is real and positive, while γ = 1 + O(m −κ− 1 2 ). It follows from (2.5.6) that we have the estimate and since we obtain that positive constant c 1 has the asymptotics c 1 = 1 + O(m −κ− 1 3 ), which allows to say thatk m,n,w0 and the true root function k m,n,w0 , which differ by a multiplicative constant, are very close. It now follows that so that k m,n,w0 has the desired asymptotic expansion in norm. In view of Proposition 2.1.3, the pointwise expansion is essentially immediate from the L 2 -estimate, at least in the region Ω τ,w0,m where which is where the functionsQ τ,w0 andQ τ,w0 are comparable in the sense that for some fixed positive constant D depending only on Q. The only remaining issue is that the error terms are slightly worse than claimed. However, by replacing κ with an integer larger than κ + 2 + A 2 D and by deriving the expansion with the indicated higher accuracy, we conclude that the desired error terms may be obtained as well.
Turning to the norm control on the set Ω τ,w0,m , we note that by elementary Hilbert space methods, we have that We need to calculate the integral on the right-hand side: where we use that χ 0 = 1 on Ω m,τ,w0 provided that m is big enough. In view of (2.5.7) the first integral on the right-hand side equals 1 + O(m −κ− 1 3 ). For any z ∈ supp(χ 0 ) \ Ω τ,w0,m we have the bound 2m(Q −Q τ,w0 )(z) ≥ A 2 D log m where D is the positive constant encountered previously. Consequently, we have the estimate If A is chosen large enough, it follows that This completes the outline of the proof.
Proof sketch of Theorem 1.4.1. The proof of Theorem 1.4.1 is entirely analogous to the above proof of Theorem 1.6.2, essentially amounting to putting τ = 0 in the latter context. In the setting of Theorem 1.4.1, there exists already a forbidden region around the point w 0 , and hence permits us to consider τ = 0. Indeed, the reason why we required that τ > 0 in the context of Theorem 1.6.2 was to allow for the instance when the off-spectral component Ω τ,w0 shrinks down to the point {w 0 } as τ → 0.

2.6.
Interface asymptotics of the Bergman density. In this section we show how to obtain the error function transition behavior of Bergman densities at interfaces, where the interface may occurs as a result of a region of negative curvature (understood as where ∆Q < 0 holds in terms of the potential Q) or as a consequence of dealing with partial Bergman kernels. Here, we focus on the the partial Bergman kernel analysis. In fact, we may think of the first instance of the full Bergman kernel as a special case and maintain that it is covered by the presented material.
The following Corollary of the main theorem summarizes the asymptotics of normalized off-spectral partial Bergman kernels in a suitable form. The domains Ω τ,w0,m are as in Theorem 1.6.2, for a given positive parameter A chosen suitably large.
Corollary 2.6.1. Under the assumptions of Theorem 1.6.2, we have the asymptotics on the domain Ω τ,w0,m , as n = τ m → +∞ while τ ∈ I 0 , where H Q,τ,w0 is the bounded holomorphic function on Ω τ,w0 whose real part equals 1 4 log(2∆Q) on the boundary, and is real-valued at the root point w 0 .
Proof. In view of the decomposition (2.2.1), this is just the assertion of Theorem 1.6.2 with accuracy κ = 1.
We proceed with a sketch of the error function asymptotics at interfaces, in particular we point out why we may proceed exactly as is done in the proof of Theorem 1.4.1 of [13].
Proof sketch of Corollary 1.7.1. We expand the partial Bergman kernel K m,n,w0 along the diagonal in terms of the root functions k m,n ,w0 , for n ≥ n. We keep τ = n m throughout. In view of Theorem 1.5.1, we have (2.6.1) K m,n,w0 (z m (ξ), z m (ξ)) e −2mQ(zm(ξ)) = +∞ n =n |k m,n ,w0 (z m (ξ))| 2 e −2mQ(zm(ξ)) , where z 0 ∈ ∂Ω τ,w0 and where z m (ξ) gives the rescaled coordinate implicitly by The rescaled Bergman density is then obtained by In view of the assumed (I 0 , w 0 )-admissibility, we may apply the asymptotic expansion in the main result, specifically in the form of Corollary 2.6.1. Since Proposition 2.1.2 tells us how the smooth Jordan curves ∂Ω τ,w0 propagate, a Taylor series expansion of the function Q −Q τ,w0 allows us to write the partial Bergman density approximately as a sum of translated Gaussians where γ 0 = γ z0,w0,Q is a positive constant. As in the proof of Theorem 1.4.1 of [13], we proceed to interpret the above sum (2.6.2) as a Riemann sum for the integral formula for the error function: This proof is complete.
3. The foliation flow for more general area forms 3.1. More general area forms. It will be desirable to obtain some flexibility on the part of the weight e −2mQ in the expansion of Theorem 1.6.2. In particular, in the following subsections we will discuss various situations in which one needs asymptotics for root functions and orthogonal polynomials with respect to measures where V is a positive C 2 -smooth function which is real-analytic in a neighborhood of the fixed smooth spectral interface of interest, which meet the polynomial growth bound for some positive constants C 1 and C 2 and some fixed integer N < +∞. We also require that for some positive constant C 3 , it holds that In particular, this covers working with the spherical area measure dA S (z) := (1 + |z| 2 ) −2 dA(z) in place of planar area measure simply by considering V (z) = (1 + |z| 2 ) −2 . Working with the spherical area measure has the advantage of invariance with respect to rotations and inversion. For a more general conformal factor V , we factor V dA = V S dA S , where V S (z) = (1 + |z| 2 ) 2 V (z), and see that our weighted measure is which has a more invariant appearance. If we write ν(z) = z −1 , the spaces of polynomials of degree at most n with respect to the L 2 -space with measure e −2mQ V S dA S becomes isometrically isomorphic to the L 2 -space of rational functions on the sphere S with a pole of order at most n at the origin, with respect to the L 2 -space with measure e −2mQ•ν V S • ν dA S . This provides an extension of the scale of root functions to zeros of negative order (i.e. poles), and the apparent similarities between orthogonal polynomials and root functions may be viewed in this light. This analogy goes even deeper than that. Assuming that 0 is an off-spectral point for the weighted L 2 -space with measure e −2mQ•ν V S • ν dA S , we may multiply by a suitable power of the conformal mapping from the off-spectral region to the unit disk D, which preserves the origin, to obtain a space of functions holomorphic in a neighborhood of the off-spectral region. Hörmander-type estimates for the∂-equation then permit us to correct the functions so that they are entire, with small cost in norm.
We note that more general area forms appear naturally from working with perturbations of the potential Q. Indeed, if we considerQ = Q − m −1 h for some smooth function h of modest growth, we have that e −2mQ dA = e −2mQ e 2h dA, which corresponds precisely to the conformal factor V = e 2h .
3.2. The asymptotics of root functions and orthogonal polynomials for more general area forms. Our analysis will show that the root function asymptotics of Theorem 1.6.2 holds also in the context of a general area form, with only a slight change in the structure of the coefficients B j,τ,w0 . Let A 2 mQ,V denote the weighted Bergman space of entire functions with respect to the Hilbert space norm The corresponding Bergman kernel is denoted by K m,V . We also need the partial Bergman spaces A 2 mQ,V,n,w0 , consisting of the functions in A 2 mQ,V that vanish at w 0 to order n or higher. These are closed subspaces of A 2 mQ,V which get smaller as n increases: A 2 mQ,V,n+1,w0 ⊂ A 2 mQ,V,n,w0 . The successive difference spaces A 2 mQ,V,n,w0 A 2 mQ,V,n+1,w0 have dimension at most 1. If the dimension equals 1, we single out an element k m,n,w0,V ∈ A 2 mQ,V,n,w0 A 2 mQ,V,n+1,w0 of norm 1, which has positive derivative of order n at w 0 . In the remaining case when the dimension equals 0 we put k m,n,w0,V = 0. As before, we call k m,n,w0,V root functions, and observe that these are the same objects we defined earlier for V = 1 in terms of an extremal problem.
Theorem 3.2.1. Under the assumptions of Theorem 1.6.2 and the above-mentioned assumptions on V , with respect to the interface ∂Ω τ,w0 , we have, using the notation of the same theorem, for fixed accuracy and a given positive real A, the asymptotic expansion of the root function on the domain Ω τ,w0,m which depends on A, where τ = n m , and the implied constant is uniform. Here, the main term B 0,τ,w0 is zero-free and smooth up to the boundary on Ω τ,w0 , positive at w 0 , with prescribed modulus The proof of this theorem is analogous to that of Theorem 1.6.2, given that we have explained how to modify the orthogonal foliation flow with respect to the general area form in Lemma 3.3.2. The lemma is applied with s = m −1 . We omit the necessary details.
We turn next to the computation of the coefficients B j,τ,w0,W in the above expansion. We recall that R τ,w0 is the potential induced by Q in the canonical positioning procedure, and we put analogously For the formulation, we need the orthogonal projection P H 2 0 of L 2 (T) onto the Hardy space H 2 0 of functions f in the Hardy space H 2 that vanish at the origin. If H τ,w0,V denotes the unique bounded holomorphic function on D, whose real part meets with Im H τ,w0,V (0) = 0, the functions B j,τ,w0,V may be obtained algorithmically as B j,τ,w0,V = c j e H τ,w 0 ,V − e H τ,w 0 ,V P H 2 0 eH τ,w 0 ,V F j for some real constants c j = c j,τ,w0,V and real-analytically smooth functions F j = F j,τ,w0,V on the unit circle T. Here, both the constants c j and the functions F j may be computed iteratively in terms of B 0,τ,w0,V , . . . B j−1,τ,w0,V .
One may further derive concrete expressions for the constants c j and the realanalytic functions F j in the above result, in terms of the rather complicated explicit differential operators L k and M k as defined in equation (1.3.4) and Lemma 3.2.1 in [13]. We should mention that the definition of the operator M k contains a parameter l, which is allowed to assume only non-negative values. However, the same definition works also for l < 0, which is necessary for the present application. In terms of the operators M k and L k , we have Here, the index set ‫ת‬ j is defined as ‫ת‬ j = {(i, k, l) ∈ Z 3 +,0 : i, l < j, i + k + l = j}, where we use the notation Z +,0 := {0, 1, 2, . . .}. This theorem is obtained in the same fashion as Theorem 1.3.7 in [13] in the context of orthogonal polynomials, and we do not write down a proof here.
3.3. The flow modified by a conformal factor. We proceed first to modify the book-keeping slightly by formulating an analogue of Definition 2.4.1, which applies to weights after canonical positioning. Definition 3.3.1. Let η and σ be given positive numbers, with η > 1. A pair (R, W ) of non-negative C 2 -smooth weights defined on D(0, η) is said to belong to the class W (η, σ) if R ∈ W(η, σ) and if the weight W meets the following conditions: (i) W is real-analytic and zero-free in the neighborhood A(η −1 , η) of the unit circle T, (ii) The polarization W (z, w) of W extends to a bounded holomorphic function of (z,w) on the 2σ-fattened diagonal annulusÂ(σ, η), which is also bounded away from 0. A collection S of pairs (R, W ) is said to be a uniform family in W (η, σ) if the weights R with (R, W ) ∈ S are confined to a uniform family in W(η, σ), while W (z, w) is uniformly bounded and bounded away from 0 inÂ(η, σ).
Fix a pair (η 0 , σ 0 ). Just as before, we let (η 1 , σ 1 ) denote a possibly more restrictive pair of positive reals with η 1 > 1 such that the relevant polarizations are hermitian-holomorphic and uniformly bounded onÂ(η 1 , σ 1 ). In connection with this definition, we recall the bound which guarantees that that if f (z, w) is holomorphic in (z,w) on the setÂ(σ 1 , η 1 ), then the function f T (z) = f (z,z −1 ) may be continued holomorphically to the annulus A(η −1 1 , η 1 ). We proceed with the main result of this section. such that for s, t small enough it holds that the domains ψ s,t D increase with t, while they remain contained in D(0, η 1 ). Moreover, for ζ ∈ T, we have For small positive s, when t varies in the interval [−β s , β s ] with β s := s 1/2 log 1 s , the flow of loops {ψ s,t (T)} t cover a neighborhood of the circle T of width proportional to β s smoothly. In addition, the main term B 0 is zero-free, positive at the origin, and has modulus |B 0 | = π − 1 4 (∆R) In order to obtain this lemma, we need to modify the algorithm which gives the original result. We proceed to sketch an outline of this modification. The omitted details are available in [13], and we try to guide the reader for easy reading.
We endow the set n with the ordering ≺ induced by the lexicographic ordering, so we agree that (j, l) ≺ (a, b) if j < a or if j = a and l < b.
Proof of Lemma 3.3.2. The conformal mappings ψ s,t are assumed to have the form ψ s,t = ψ 0,t + (j,l)∈ 2κ+1 j≥1 s j t lψ j,l for some bounded holomorphic coefficientsψ j,l and a conformal mapping We make the following initial observation. In the limit case s = 0, the flow equation (3.3.1) of the lemma forces ψ 0,t to be a mapping from D onto the interior of suitably chosen level curves of R, and from this we may obtain the coefficientsψ 0,l . Indeed, if we take logarithms of both sides of the equation and multiply by s we obtain where the J Ψ denotes the associated Jacobian for t in the interval [−β s , β s ], so that J Ψs gets defined on the annulus A(1−β s , 1+β s ), provided that the coefficient functions which define ψ s,t can be found. Assuming some reasonable stability with respect to the variable s as s → 0 + in (3.3.4), we obtain in the limit that In particular, the loop ψ 0,t (T) is a part of the level set where R = 1 2 t 2 . This level set consists of two disjoint simple closed curves, one on either side of T, at least for small enough t and locally near T. For t > 0, we choose the curve outside the unit circle, while for t < 0 we choose the other one. We normalize the mapping ψ 0,t so that it preserves the origin and has positive derivative there. In this fashion, the coefficientsψ 0,l get determined uniquely by the level set condition. We note that the smoothness of the level curves was worked out in some detail in Proposition 4.2.5 in [13]. Moreover, the coefficient functionsψ 0,l are given in terms of Herglotz integrals as in Proposition 4.6.1 of [13], with the obvious modifications.
Our next task is to obtain iteratively the coefficients B j for j = 0, 1, 2, . . . and the higher order corrections to the conformal mapping, given in terms of the coefficientŝ ψ j,l for j = 1, . . . , κ. It turns out to be advantegeous to work instead with h s = log f s as the basic object of study, and make the ansatz for ζ ∈ T, where we recall that R • ψ 0,t = 1 2 t 2 , and put We need to show that (3.3.8) s,t (ζ) = − 1 2 log(4π) + O(|s| κ+ 1 2 + |t| 2κ+1 ), ζ ∈ T when s → 0 while |t| ≤ s 1 2 log 1 s . Since s,t should be smooth in the parameters s and t, by the multivariate Taylor formula, this is equivalent to having the system of equations for ζ ∈ T and (j, l) ∈ 2κ \ (0, 0), fulfilled. If, in addition, we can show that the functions b j andψ j,l remain holomorphic and uniformly bounded in the appropriate domains provided that (R, W ) remains confined to a uniform family of W (η 0 , σ 0 ), the result follows.
How to solve for the unknown coefficient functions. We proceed to solve the system (3.3.9). The approach is to first express the Taylor coefficientsˆ j,l in terms of our unknown coefficient functions b 0 , . . . , b κ andψ j,l for (j, l) ∈ 2κ+1 , and then, as a second step, to insert the system of equations (3.3.9). We determine the unknowns by an iterative procedure. At a given step in the iteration, some of the coefficient functions will be already found. We split the equation (3.3.9) into a term containing precicely one unknown coefficient function and a second term which contains only already determined coefficient functions. Here, we refrain from giving a complete account, which is available in [13] modulo minor modifications to fit the present setup. The higher order Taylor coefficients of the function s,t with respect of s and t are given as follows. When j ≥ 0 and l ≥ 1, we have (3.3.10) 0 =ˆ j,l (ζ) = −2(4∆R(ζ)) 1 2 Re(ζψ j+1,l−1 (ζ)) + T j,l,W (ζ), ζ ∈ T, while for j ≥ 1 and l = 0, we have Here, we inserted the equation (3.3.9) for added convenience. The expressions T j,l,W are real-valued real-analytic functions which are uniformly bounded while (R, W ) remains in a uniform family in W (η 0 , σ 0 ), and may be explicitly written down using the multivariate Faà di Bruno's formula. The crucial point for us is the dependence structure of the functions T j,l,W , which remains the same as in the algorithm for the orthogonal polynomials: (T-i) For l ≥ 1, the function T j,l,W is an expression in terms of the functions b 0 , . . . , b j as wellψ p,q for indices (p, q) ∈ 2κ+1 with (p, q) ≺ (j + 1, l − 1), and also involves R and W ,
The dependence is basically multivariate polynomial dependence, with the need to allow also for partial derivatives of some of the given expressions. This aspect is described in great detail in terms of the polynomial complexity classes introduced in Section 4.8 of [13]. That (T-i)-(T-ii) hold follows by noticing that Propositions 4.2.5, 4.4.1 and 4.6.1 in [13] remain basically unchanged, while Propositions 4.7.1, 4.9.1, 4.10.1 and 4.12.1 in [13] require the obvious modifications related to our replacing the weight e −2mR on the exterior disk D e (0, ρ) by a weight e −2mR W on the disk D(0, ρ). In particular, it is important that the real-analytic function W is strictly positive in a fixed neighborhood of the unit circle T, so that log W is a real-analytic function in the same region. A natural approach to the computations is to write which is essentially the expression which gets expanded in [13], and introduce the modification II s,t = log W • ψ s,t , since then s,t = I s,t + II s,t . If we introduce the notationˆ it follows immediately thatˆ I j,l are as in Proposition 4.12.1 of [13], while the latter may be computed explicitly using the multivariate Faà di Bruno formula [7]. The leading behavior ofˆ j,l comes from the contribution ofˆ I j,l , which is as in [13], while it is easily verified thatˆ II j,l for (j, l) ∈ 2κ is an expression in terms of the coefficient functionsψ p,q with p ≤ j and q ≤ l and of some partial derivatives of log W . But then in particular (p, q) ∈ 2κ with (p, q) ≺ (j, l + 1). It is now immediate that the remainders T j,l,W have the indicated properties.
Sketch of the algorithm. We sketch the solution algorithm below, and indicate where the necessary modifications to the corresponding steps in [13] are required. We need the notation Γ t for the loop with The loop is a perturbation of the unit circle T for t close to 0. However, for t = 0 there are two such loops, one which is inside the circle T, and one which is outside.
For t > 0 we choose the outside loop, whereas for t < 0 we instead choose the inside loop. This way, the domain enclosed by Γ t grows with t.
Step 1. Take as above ψ 0,t to be the conformal mapping ψ 0,t : D → D t , where ψ 0,t (0) = 0 and ψ 0,t (0) > 0, and D t denotes the bounded domain enclosed by the curve Γ t . It follows from the smoothness of the flow of the level curves Γ t (for details see Proposition 4.2.5 in [13]) that we have an expansion which determines the coefficient functionsψ 0,l for all l = 0, 1, 2, . . . (see Proposition 4.6.1 [13]). In particular, we have Step 2. The equation (3.3.9) for (j, l) = (0, 0) together with (3.3.12) gives that which in its turn determines b 0 . Indeed, the only holomorphic function which is real-valued at the origin and meets this equation is where the Herglotz operator H D is given by The definition of the Herglotz operator is extended to the boundary T via nontangential boundary values, and, whenever possible, across the circle T by analytic continuation. As a consequence of the assumptions on the potential Q and the conformal factor V as well as the regularity of the conformal mapping ϕ τ,w0 , the function W is real-analytic and positive in a neighborhood of T, and moreover the pair (R, W ) meets the regularity requirements of Definition 3.3.1. This means that essentially we are in the same setting as explained in [13]. For instance, we may conclude that the function b 0 given by (3.3.13) extends as a bounded holomorphic function on a disk D(0, η 1 ) with radius η 1 > 1.
Step 3. To begin with, we have an integer 1 ≤ j 0 ≤ κ for which we have already successfully determined the coefficient functions b j for 0 ≤ j ≤ j 0 − 1 as well asψ j,l for all (j, l) ∈ 2κ+1 with (j, l) ≺ (j 0 , 0). Note that is known to be so for j 0 = 1, by Steps 1 and 2. In this step, we intend to determine all the coefficient functionsψ j,l with (j, l) ∈ 2κ+1 and (j, l) ≺ (j 0 + 1, 0), and keep track of the equations (3.3.9) with (j, l) ∈ 2κ that get solved along the way. The induction hypothesis includes the assumption that the system of equations (3.3.9) holds for all (j, l) ∈ 2κ with l ≥ 1 and (j, l) ≺ (j 0 − 1, 1) (this is vacuous for j 0 = 1), as well as for (j, l) with 0 < j < j 0 and l = 0. We need only findψ j,l for (j, l) ∈ 2κ+1 with j = j 0 . We do this by induction in the parameter l, starting with l = 0. Assume for the moment that it has been carried out for all l = 0, . . . , l 0 − 1. The coefficientψ j0,l0 which we are looking for appears as the leading term in the equation (3.3.10) corresponding to (j, l) = (j 0 − 1, l 0 + 1), which reads −2(4∆R) 1 2 Re(ζψ j0,l0 ) + T j0−1,l0+1 = 0 on T, where T j0−1,l0+1 is an expression in the already determined coefficient functions, by property (T-i). We solve for the coefficient functionψ j0,l0 in terms of the Herglotz operator:ψ As for the first step l 0 = 0, the above formula applies in that case as well. What is important is that then, the function T j0−1,1,W only depends data known at the beginning of Step 3, in view of (T-i). By the correspondence (3.3.3) and the induction hypothesis, we have made sure that the system of equations (3.3.9) holds for all pairs (j, l) ∈ 2κ with l ≥ 1 and (j, l) ≺ (j 0 , 1) as well as for (j, l) with 0 < j < j 0 and l = 0. This completes Step 3.
Step 4. After having completed Step 3, we find ourselves in the following situation: the coefficient functions b j are known for j < j 0 , whileψ j,l are all known for (j, l) ∈ 2κ+1 with (j, l) ≺ (j 0 + 1, 0). We proceed to determine b j0 using the equation (3.3.11) with index (j, l) = (j 0 , 0). That equation asserts that 2 Re b j0 + T j0,0,W = 0 on T, where T j0,0,W depends on the known data, by property (T-ii). We solve for b j0 using the formula b j0 = − 1 2 H De [T j0,0 ], and we get a function b j0 with b j0 (0) ∈ R. In view of Step 3, this choice makes sure that (3.3.9) holds for all pairs (j, l) with l = 0 and 0 < j ≤ j 0 , as well as for all (j, l) ∈ 2κ with l ≥ 1 and (j, l) ≺ (j 0 , 1). This completes Step 4, and we have extended the set of known data so that we may proceed to Step 3 with j 0 replaced with j 0 + 1. Here, we should insert a word on smoothness. If f is real-analytically smooth along T, then H D [f ] gets to be holomorphic in a larger disk D(0, η) for some η > 1. Hence real-analytic smoothness carries over to the next step in the iterative procedure.
The above algorithm continues until all the unknowns have been determined, up to the point where the whole index set 2κ+1 has been exhausted. In the process, we have in fact solved the system of equations (3.3.9) for all indices (j, l) ∈ 2κ . This means that if we form h s and ψ s,t in terms of the functions b j andψ j,l obtained with the above algorithm, and put f s = exp(h s ), we find by exponentiating the Taylor series expansion of s,t in the parameters s and t that the flow equation (3.3.1) holds. This completes the sketch of the proof of the lemma.

3.4.
Orthogonal polynomial asymptotics for a general area form. The results concerning root functions for more general area forms apply also to the setting of orthogonal polynomials. Although many things are pretty much the same, we make an effort to explain what the precise result is in this context.
We need an appropriate notion of admissibility of the potential Q which applies to the setting of orthogonal polynomials. We recall that the spectral droplet S τ is the contact set S τ = {z ∈ C :Q τ (z) = Q(z)}, which is typically compact, whereQ τ is the function Q τ (z) := sup q(z) : q ∈ SH(C), q ≤ Q on C, q(z) ≤ τ log(|z| + 1) + O(1) .
Definition 3.4.1. We say that the potential Q is τ -admissible if the following conditions are met: (i) Q : C → R is C 2 -smooth, (ii) Q meets the growth bound (iii) The unbounded component Ω τ of the complement of the spectral droplet S τ is simply connected on the Riemann sphereĈ, with real-analytic Jordan curve boundary, (iv) Q is strictly subharmonic and real-analytically smooth in a neighborhood of the boundary ∂Ω τ .
As before, we consider measures e −2mQ V dA, where the conformal factor V is assumed to be nonnegative, positive near the curve ∂Ω τ , and real-analytically smooth in a neighborhood of ∂Ω τ with at most polynomial growth or decay at infinity (3.1.1). We denote by φ τ the surjective conformal mapping φ τ : Ω τ → D e , which preserves the point at infinity and has φ τ (∞) > 0. The function Q τ is defined as the bounded holomorphic function on Ω τ whose real part equals Q on the boundary ∂Ω τ , and whose imaginary part vanishes at infinity.
The orthogonal polynomials P m,n,V have degree n, positive leading coefficient, and unit norm in A 2 mQ,V . They have the additional property that P m,n,V , P m,n ,V mQ,V = 0, n = n .
We will work with τ = n m . Theorem 3.4.2. Suppose Q is τ -admissible for τ ∈ I 0 , where I 0 is a compact interval of the positive half-axis. Suppose in addition that V meets the above regularity requirements. Given a positive integer κ and a positive real A, there exists a neighborhood Ω τ of the closure of Ω τ and bounded holomorphic functions B j,τ,V on Ω τ , as well as domains Ω τ,m = Ω τ,m,κ,A with Ω τ ⊂ Ω τ,m ⊂ Ω τ which meet dist C (Ω c τ,m , Ω τ ) ≥ Am − 1 2 (log m) on Ω τ,m as n = τ m → +∞ while τ ∈ I 0 , where the error term is uniform. Here, the main term B 0,τ,V is zero-free and smooth up to the boundary on Ω τ , positive at infinity, with prescribed modulus In view of Lemma 3.3.2, the construction of approximately orthogonal quasipolynomials may be carried out in the same way as in the case when V = 1.
3.5.∂-surgery in the context of more general area forms. As for the∂surgery technique in the presence of a conformal factor V , we need to mention some modifications that are required. We recall that V meets the conditions (3.1.1) and (3.1.2) and put We note that by (3.1.2), the modified potential Q m is subharmonic and that Q m is strictly subharmonic uniformly in a neighborhood of ∂Ω τ , provided that m is large enough. Here, Ω τ is the unbounded component of the off-spectral set S c τ . Moreover, it follows from the growth and decay bounds (3.1.1) that Q m meets the required growth bound for large enough m. We consider the solutionQ τ,m to the obstacle problemQ τ,m (z) = sup u(z) : u ∈ SH τ (C), u(z) ≤ Q m (z) on C .
The standard regularity theory applies, and implies thatQ τ,m is C 1,1 -smooth. By our smoothness assumptions on S τ , as well as the assumed regularity of the potential Q and the conformal factor V , the coincidence set S τ,m := z ∈ C :Q τ,m (z) = Q(z) moves only very little away from S τ for large m. Indeed, by the main theorem of [18], the free boundary moves with a smooth normal velocity under a smooth perturbation, and it follows that ∂S τ,m is contained in a O(m −1 )-neighborhood of ∂S τ . We note thatQ τ,m is automatically harmonic off S τ,m and that ∆Q τ,m ≥ 0 holds generally.
We now apply Hörmander's classical∂-estimate with the obstacle solutionQ τ,m as potential. Let f ∈ L ∞ (S τ,m ), which vanishes whenever ∆Q τ,m = ∆Q m = 0 on S τ,m . Then there exists a solution u to the problem provided that the right-hand side is finite. Moreover, since u is holomorphic off the compact set S τ,m , this estimate implies a polynomial growth bound u(z) = O(|z| n−1 ) as |z| → +∞.
After these modifications, the∂-surgery may be performed as in the earlier context of the planar area measure dA. For the details, we refer to Section 4.9 in [13] and the comments in Section 2.5 above.