Moment estimates of the cloud of a planar measure

With a proper function theoretic definition of the {\it cloud} of a positive measure with compact support in the real plane, a computational scheme of transforming the moments of the original measure into the moments of the uniformly distributed mass on the cloud is described. The main limiting operation involves exclusively truncated Christoffel-Darboux kernels, while error bounds depend on the spectral asymptotics of a Hankel kernel belonging to the Hilbert-Schmidt class.


Introduction
This is a continuation of a recent work [26], adding to it some constructive analysis details. The purpose of both articles was/is to compute the moments of the "cloud" of a positive measure with compact support in the plane starting from the moments of the original measure. This moment conversion eliminates the components of the measure which are singular with respect to Lebesgue area, informally regarded as "outliers", and charges the cloud with uniform mass. Two real dimensions are specific in the approximation scheme we propose, twofold: complex variables enter heavily into the picture, as well as a refined spectral analysis for pairs of self-adjoint operators with trace-class commutator.
The present article has two interlaced motivations: to identify a 2D analogue of the classical and still rapidly evolving spectral analysis on the line or the circle, based on moment information, involving orthogonal polynomials and Christofell-Darboux kernels, and second, the current renewed interest of data analysis flavor, consisting in separating outliers from the cloud, in multidimensional point distributions. Our viewpoint is to treat moments as observed data and to extract from them by simple analytic transforms the moments of the cloud. We only mention without expanding in this work the well known, various reconstruction from moments methods.
A landmark contribution to polynomial approximation theory in the complex domain is Thomson's Theorem [31]. It asserts that complex polynomials are either dense in the Lebesgue space L 2 (µ) associated to a positive Borel measure µ with compact support in the plane, or there exist L 2 (µ)-bounded point evaluations for polynomials. More specifically, Thomson's Theorem provides a decomposition of the measure µ into a singular with respect to area part, where polynomials are dense, and measures which carry an open set of point evaluations, see Theorem 2.2 below for the precise statement. The multiplier with the complex variable S µ = M z , acting on the closure of polynomials P 2 (µ), is consequently decomposed into a direct sum of a normal operator and a collections of irreducible subnormal, cyclic operators with a continuum of eigenvalues for their adjoints: S µ = N ⊕ S 1 ⊕ S 2 ⊕ . . . . In very elementary terms, the main idea we exploit in this article is: the self-commutator . . does not "see" the normal part, that is the outliers to the cloud.
The spectral analysis of the analytic Toeplitz operator S µ is rich and provides the numerical procedure we propose to eliminate via traces of commutators its normal direct summand. A second notable result we invoke is due to Berger and Shaw [5], namely the self-commutator [S * µ , S µ ] is trace class. This opens a wide array of tools specific to the notion of principal function, a refined spectral invariant extending the Fredholm index across the spectrum. A trace formula discovered in the 1970-ies by Helton and Howe and an equivalent determinant formula independently obtained by Carey and Pincus provide effective formulas linking the moments of the principal function (in our case the characteristic function of the cloud of the measure µ) and traces of commutators of smooth functions applied to S µ . The third deep result we rely on is due to Carey and Pincus, asserting that the principal function of a subnormal operator is integer valued [8].
This machinery was put in action in our previous article [26] with a proposed computational/transformation scheme: moments of µ −→ moments of the uniform mass on the cloud of µ.
We expand below a quantitative analysis of the algorithm, with error bounds expressed in the entries of the Hessenberg matrix representing S µ with respect to the filtration of P 2 (µ) by finite dimensional subspaces consisting of polynomials. The main transform involves solely limits of some explicit expressions involving truncated Christoffel-Darboux kernels of the original measure µ. And not surprising, the principal error estimate depends on the asymptotics of the singular numbers of Hankel's operator (I − P )M * z P , where P : L 2 (µ) −→ P 2 (µ) denotes the orthogonal projection. The approximation scheme can be adapted to other orthogonal function systems than polynomials, with an appropriate sequence of finite rank projections converging strongly to the identity operator. We remain at an all-inclusive level, applicable to any positive measure with compact support, leaving for future studies the adaptation of this general framework to concrete situations. The examples at the end of the article indicate a few openings in this direction.
To make the present note accessible to a larger group of readers, we (re)expose below the main ingredients, with precise reference to the sources.
Acknowledgements. We are grateful to Bernhard Beckermann for helpful comments on the asymptotics of the Hessenberg matrix associated to a uniform mass distribution in the plane.

Preliminaries
Throughout this note C[z] stands for the algebra of polynomials in one variable with complex coefficients. Let µ be a positive Borel measure with compact support in the complex plane. We assume the support of µ is not finite.
2.1. Orthogonal polynomials. We recall a few facts and conventions referring to complex orthogonal polynomials, as for instance appearing in [29]. The closure of complex polynomials in L 2 (µ) is denoted P 2 (µ), with associated orthogonal projection P . The multiplication M by the complex variable z ∈ C is a bounded linear transform of L 2 (µ) which leaves invariant the subspace P 2 (µ): The linear operator S µ = M | P 2 (µ) : P 2 (µ) −→ P 2 (µ) is called subnormal. It is also a Toeplitz operator: (S µ f )(z) = P (wf (w))(z), f ∈ P 2 (µ). Sometimes we simply write S = S µ . The spectrum of M coincides with the closed support supp(µ) of the measure µ, while the spectrum of S can be larger, containing in addition some connected components of C \ supp(µ).
The constant function 1 is a cyclic vector for S, producing the finite dimensional filtration (Krylov subspaces): We denote by p n (z) the associated complex orthogonal polynomials: The operator S has a distinguished Hessenberg matrix representation with respect to the orthonormal basis (p n ) ∞ n=0 of P 2 (µ): zp j , p k = h kj , j, k ≥ 0, observing the automatic vanishing relations The adjoint operator is represented by the matrix: It is customary to normalize the leading coefficient of p n to be a positive number: The matrix representing S has the form: If the measure µ has support on the real line, then S is a self-adjoint operator, hence the matrix H is symmetric, with real values, with only three non-vanishing diagonals: a well charted Jacobi matrix framework. For every non-negative integer n one can speak of the reproducing kernel (known as the Christoffel-Darboux kernel) of the finite dimensional subspace C n [z]: characterized by the identity: The Christoffel function of order n is These numerical values decrease with n, and the limit is known as the Christoffel function associated to the measure µ, evaluated at the point z.
Incidentally, the leading coefficient of the orthogonal polynomial p n has a similar variational interpretation:

2.2.
Trace and determinant formulas. The analytic subnormal operator S = P M P satisfies the hyponormal commutator inequality [S * , S] ≥ 0. Indeed, for any element f ∈ P 2 (µ): A theorem of Berger and Shaw [5] asserts that the above commutator is trace-class: Tr[S * , S] < ∞. A different proof of Berger-Shaw Theorem is due to Voiculescu [32], in an article putting this very trace bound in the context of general perturbation theory, with a specific link to the modulus of quasi-triangularity of a Hilbert space linear transform. For any polynomials p ∈ C[z, z], one can define an ordered functional calculus (traditionally called hereditary functional calculus) p(S, S * ), by arranging all powers of S * to the left of the powers of S, in every monomial. Then for a pair of polynomials p, q ∈ C[z, z] one proves by degree induction that the commutator [p(S, S * ), q(S, S * )] is trace-class. A remarkable observation due to Helton and Howe [15] Above, dA stands for Lebesgue measure in R 2 .
Definition 2.1. The cloud of a positive Boreal measure µ with compact support in C is the measurable set Σ(µ) appearing in trace formula (2.1).
Note that the set Σ(µ) is only determined up to area null-sets. In practice it is more appropriate to speak about the class [χ Σ(µ) ] of its characteristic function in L 1 (C, µ).
The equivalent formulation as an infinite determinant formula goes back to the origins of Carey and Pincus work: All results above touch the surface of the theory of hyponormal operators, in particular with reference to the principal function of a semi-normal operator. We do not expand here the details, referring to the monograph [18] for complete proofs and historical comments. In the present article we focus on the reconstruction of the cloud Σ(µ) contained in the polynomial convex hull of the support of the measure µ. Note for instance that Tr[S * , S] = 1 π Area Σ(µ), hence the operator S is normal if and only if Area Σ(µ) = 0. We will see shortly that this in turn is equivalent to the density of complex polynomials in L 2 (µ), quite a central topics in approximation theory. For the limited aim of the present work, we confined ourselves to reproduce some well-known computations specific to Hankel operators. The general theory of Hankel operators is superbly exposed in [25]. The following observations appeared in a slightly different form in [7], but they are not uncommon to the general spectral analysis of Toeplitz operators, see for instance [34].
More precisely, the self-commutator of the subnormal operator S = P M P = M P can be factored as: is Hilbert-Schmidt: Tr T * T < ∞. We can express this fact in two different ways: or, there exists a kernel function L ∈ L 2 (C × C, µ ⊗ µ) satisfying and in particular The singular numbers of T will play a central role below.
2.3. Function theory. The fine structure of the subnormal operator S = S µ was elucidated by a rather recent discovery by Jim Thomson. We state in full the main theorem.
Theorem 2.2 (Thomson). Let µ be a positive Borel measure, compactly supported on C. There exists a Borel partition ∆ 0 , ∆ 1 , . . . of the closed support of µ with the following properties: 2) Every operator S µ j , j ≥ 1, is irreducible with spectral picture: The proof appeared in [31], for L p spaces, 1 ≤ p < ∞. A conceptually simpler proof appears in [6]. The central position of Thomson's Theorem was immediately recognized in the monograph [9] (published almost simultaneously with the original article). Corollary 2.3. The cloud Σ(µ) of a positive measure µ with compact support in C is empty if and only if the complex polynomials are dense in L 2 (µ). In case Σ(µ) is non-empty, it contains interior points: The operator S µ 0 = M z ∈ L(L 2 (µ 0 )) is the single normal component of S µ while the summands S µ j collect the non-normal behavior of S µ . For every index j, j ≥ 1, the theorem asserts that S µ j admits a continuum of eigenvalues of multiplicity one, filling the simply connected open set G j : Moreover, it is known that the range of S µ j − λ is closed for such spectral parameters. The corresponding eigenvectors span P 2 (µ j ), to the extent that this functional Hilbert space carries a reproducing kernel. This feature is detected by the local Christoffel function: , f (λ) = 1} > 0, and ultimately by the full Christoffel function: Corollary 2.4. The cloud of a positive measure with compact support µ is non-empty if and only if there exists a point λ ∈ C satisfying Λ µ (λ) > 0.
Szegö's theory of orthogonal polynomials on the unit circle T provides a clear cut picture: if supp(µ) ⊂ T, then Σ(µ) is either empty or equal to the full disk D. The latter case is characterized by Szegö's condition: or equivalently by Christoffel function test: for at least one point λ ∈ D one has Λ µ (λ) > 0, and then the same is true for all z ∈ D. See for instance [1]. Another, related and notable particular case is represented by measures µ subject to the finiteness condition rank[S * , S] < ∞. This situation was fully analyzed by McCarthy and Yang [19,20] with the following conclusion: Thomson's decomposition admits finitely many irreducible summands (that is 0 ≤ j ≤ n) with the normal part represented by a positive measure µ 0 which is singular with respect to harmonic measure on every G j , 1 ≤ j ≤ n, and does not put mass on ∪ n j=1 G j ; in addition every simply connected component G j is a quadrature domain. In their theorem a complete description of the measures µ j , 1 ≤ j ≤ n, is provided. Moreover, under this finite rank condition, the intersection of any two distinct sets G j ∩ G k contains at most one single point, and the union F = ∪ n j=1 G j is polynomially convex. In other terms the complement C \ F is connected.
A bounded open set Ω of the complex plane is called a quadrature domain for analytic functions if there is a distribution of finite support τ in Ω (combination of point masses and their derivatives), such that where L 1 a (Ω, dA) stands for the space of analytic functions in Ω which are Lebesgue integrable. The order of a quadrature domain is the number of nodes, counting multiplicity, in the above cubature formula. A connected quadrature domain has an irreducible real algebraic boundary. The simplest example is of course a disk. Indeed, on the unit disk by Gauss mean value theorem. Any simply connected quadrature domain is a conformal image of the disk by a rational function, for instance a cardiodid or a connected lemniscate of degree four. An informative and accessible survey of quadrature domains, with ramifications to potential theory, inverse geophysical problems, quantum physics and operator theory is [14]. We record these facts under a precise statement.
Proposition 2.5. Let µ be a positive Borel measure with compact support in C. The cloud Σ(µ) appearing in trace formula (2.1) contains every open set G j , j = 0, appearing in Thomson's decomposition of the operator S µ . In case rank[S * µ , S µ ] < ∞ there are only finitely many summands ( 1 ≤ j ≤ n) and the cloud Σ(µ) coincides up to an area null-set with the quadrature domain Ω = The nature of the measure µ is also quite explicit in the case of finite rank selfcommutator, as follows. Let r : D −→ Ω denote a conformal rational map onto a bounded quadrature domain Ω. Let τ be an absolutely continuous measure with respect to harmonic measure ω for Ω, supported on ∂Ω and satisfying Let µ = τ + ν, where ν is a positive, finitely supported measure on Ω. Then complex polynomials are not dense in L 2 (µ) and the operator S µ is a typical irreducible subnormal operator possessing finite rank self-commutator, cf. Theorem 1.12 in [20]. In this case Ω coincides with the cloud Σ(µ).

Finite rank approximation
3.1. Commutator inequalities. Let φ ∈ C ∞ (C). The operator of multiplication by φ on L 2 (µ) is denoted M φ or by the same letter if there is no confusion. Obviously M φ is bounded. Jacobi's identity: ]. In other terms, for an element f ∈ L 2 [µ] we find: From here we deduce at the level of integral kernels: In particular, the commutator [φ, P ] is Hilbert-Schmidt for every smooth function φ. Observe that for φ complex analytic and f ∈ P 2 (µ) one has [φ, P ]f = φP f − P (φf ) = φf − φf = 0.
One step further, we can write explicitly the integral kernel of the commutator appearing in trace formula (2.1). For a polynomial p ∈ C[z, z] either one computes p(S, S * ) via an ordered functional calculus, or one takes P M p P , the result differs by a trace class operator: Hence the trace of [p(S, S * ), S] is not affected by such a rearrangement of terms. Pure algebra yields: Summing up, The above formula extends to all Toeplitz operators T φ = P M φ P with a smooth symbol. In addition, trace formula yields: One can deduce from above a few qualitative properties of the integral kernel L(z, w), by exploiting the distributional sense of the above identity. In the present article we are concerned only with some approximation estimates. The ad-hoc notation z − w is adopted in the following pages, as well as the Besov type semi-norm where F is a closed subset of the complex plane. Identity 3.1 implies, for every smooth function φ: [φ, P ] HS ≤ L 2,µ×µ |φ| supp(µ) .
Let Q denote an orthogonal projection of L 2 (µ) onto a closed subspace of P 2 (µ). For a fixed smooth function φ one finds, along the above computations: The truncated operator [M, P ] * Q has finite rank and it is represented by an integral kernel: More importantly, In conclusion we have proved the following result.
Theorem 3.1. Let µ be a positive Borel measure with compact support in C and let φ ∈ C ∞ (C). For an orthogonal projection Q of L 2 (µ) onto a closed subspace of P 2 (µ) the estimate holds true.
Remark that one can express the distortion factor above in function space norm: The particular case of a finite rank projection Q onto a subspace of polynomials will be of interest in the next sections.
Hankel operator T = [M * , P ] = (I − P )M * P : identifies the singular numbers (κ j ) together with the eigenvectors of its modulus: By convention we include here the null vectors of T , so that (f j ) ∞ j=0 is an orthonormal basis of P 2 (µ). And similarly g j are eigenvectors of T T * , forming an orthonormal basis of P 2 (µ) ⊥ . Courant-Fisher's min-max principle implies the following bound.

This observation is effective in case rank
is an integrable function with respect to µ. These computations lead to a theoretical result, illustrating the key role played by the asymptotics of the singular numbers of the operator T . Then L(z, z) ∈ L 1 (µ) and Regardless to say that there are ample studies and criteria assuring a Hankel operator to be trace-class, see [25].
To this aim, fix an element f ∈ P 2 (µ) and remark that in the weak topology of L 2 (µ × µ).
This observation is consistent with in the even weaker sense of distributions.

Polynomial approximation.
We specialize below the finite rank approximation of the commutator [M * , P ] to the filtration by polynomial subspaces C n [z] labelled by the degree. The notation introduced in the previous sections is unchanged: P n is the orthogonal projection onto C n [z], Christoffel-Darboux kernel of order n is K n (z, w), Hessenberg's matrix associated to the orthogonal polynomials (p j ) is [h jk ]. The rate of convergence in Hilbert-Schmidt norm of [M * , P ]P n to [M * , P ] controls via Theorem 3.1 the finite central convergence in the trace formula. For a fixed degree j ≥ 0 one finds That is, the square norm of the j-th column in the Hessenberg matrix minus the square norm of the j-th row. Therefore We know that as a function of n this is an increasing sequence, converging to Tr[S * , S] = AreaΣ(µ) π . At the level of integral formulas one puts in motion the Christoffel-Darboux kernel. Denoting by L n (z, w) the integral kernel of [M * , P ]P n one finds L n (z, w) = L(z, ζ)K n (ζ, w)dµ(ζ), n ≥ 0, and one step further, approximating the full projection P by P N with N → ∞ we infer a formula in terms only of the moments of the original measure µ. Given a polynomial R(z, z) ∈ C[z, z] it will be convenient to separate the powers of z: Notice thatR is a polynomial in all variables and its restriction to the diagonal z = w coincides with ∂R ∂z (z, z). The main result follows. We denote byσ the polynomial convex hull of a closed set σ ⊂ C.
Theorem 3.6. Let µ be a positive Borel measure with compact support σ in C, with associated Hessenberg matrix (h jk ) ∞ j,k=0 and Christoffel-Darboux kernels K n (z, w), n ≥ 0. The moments of the area measure supported by the cloud of µ can be computed by the formula: For fixed values of n < N the error ǫ N,n = ǫ N,n (µ) in the above limit satisfies: and π ∞ j=0 s j = Area Σ(µ). Proof. By relaxing the two projectors P appearing in the trace formula 1 π Σ(µ) Helton and Howe trace formula (2.1) is coordinate free, and moreover, it is invariant under additive Hilbert Schmidt perturbations S +L of the operator S, provided the selfcommutator [S * +L * , S +L] remains trace-class [32]. In particular, an adapted sequence of finite rank orthogonal projections converging strongly to the identity operator may well improve the estimates in Theorem 3.6. We discuss an example in this direction.
Corollary 3.7. Let µ be a positive measure of compact support and let ν be a finite atomic measure supported by the complement of the polynomial hull of suppµ. There exists a constant ρ > 1 with the property with the second term independent of N > n.
Proof. By finite recurrence and normalization we can assume ν is a Dirac mass at a point a.
The map J : , is isometric and has dense range by Runge approximation theorem. Hence J is a unitary transformation. Denote by π n the orthogonal projection of P 2 (µ + ν) onto C n [z].
The essential Hilbert-Schmidt norm involved in the bound of ǫ N,n (µ + ν) is and an orthonormal basis of the subspace (z − a)C n−1 [z] ⊂ P 2 (µ), simultaneously orthonormal via the unitary map J in P 2 (µ) ⊕ P 2 (ν), produce the estimate Since the point a does not belong to the polynomial convex hull of supp(µ) the variational definition of Christoffel function Λ µ n (a) implies Λ µ n (a) ≤ Cρ −n , for some constants C > 0 and ρ > 1. The identity K µ n (z, a) K µ n (a, a) + K µ n (a, a) 2 2 µ = K µ n (a, a) K µ n (a, a) + K µ n (a, a) 2 = Λ µ n (a) 1 + Λ µ n (a) completes the proof.

Padé type approximation scheme in 2D
The reconstruction of the cloud of a measure µ from its moments can be completed via different paths: geometric tomography, Bergman space methods [13] (if applicable), identification of potential real-algebraic boundary [17], curve fitting along the boundary [2], or by exploiting the same hyponormal operators tools invoked in the previous sections. We reproduce from our preceding article a few details on the latter approximation scheme.
Let g ∈ L 1 comp (C, dA), 0 ≤ g ≤ 1, be a fixed shade function, and think of it as the characteristic function of the shade Σ(µ) of a measure. We consider the moments given. They can be organized in the exponential of a formal generating series: This is of course the power expansion at infinity of the double Cauchy integral appearing in (2.2). Denote in short: We recall a few of the properties of the exponential transform E g : a). The function E g can be extended by continuity to C 2 by assuming the value E g (z, z) = 0 whenever C g(ζ)dA(ζ) |ζ−z| 2 = ∞; b). The function E g (w, z) is analytic in w ∈ C \ supp(g) and antianalytic in z ∈ C \ supp(g); c). The kernel 1 − E g (w, z) is positive semi-definite in C 2 ; d). The behavior at infinity contains as a first term the Cauchy transform of g: The case of a characteristic function g = χ Ω of a bounded domain Ω is particularly relevant for our note. In this case we simply write E Ω instead of E g , and we record the following properties (all proved and well commented in [12]).
The kernel E Ω is also characterized by a Riemann-Hilbert factorization, see [12]. The feature which turns the exponential transform into a suitable shape reconstruction from moments tool is its rationality on quadrature domains, in complete parallelism to the rationality of Cauchy transforms of point masses on the line, or the rationality of exponential of Cauchy transforms of union of intervals.
We reproduce the rational reconstruction procedure. Let d be a fixed integer and let (a kℓ ) d k,ℓ=0 , be a non-negative matrix of potential moments of a "shade function" g(z), 0 ≤ g ≤ 1. Consider the truncated exponential transform is also a non-negative definite matrix, subject to some additional constraints [18]. A necessary and sufficient condition that (a kℓ ) d k,ℓ=0 represent the moments of a quadrature domain of order d is det(b mn ) d m,n=0 = 0, or equivalently the existence of a monic polynomial P (z) of degree d and a rational function of the form such that, at infinity The reader will recognize above a typical 2D Padé approximation scheme. Moreover, for any shade function g, the exponential transform E g coincides with E Ω , where Ω is a quadrature domain if and only if , |z|, |w| ≫ 1.
In this case the zeros of P coincide with the quadrature nodes, while the numerator is the irreducible defining polynomial of the boundary of Ω: c mn z m z n = |P (z)| 2 }.
The above Padé approximation procedure was proposed for the reconstruction of planar shapes in [11], with additional details in [12].

5.1.
Reconstruction of a disk via its exponential transform. As simple and well known the example below might be, it is illustrative for the two dimensional Padé scheme just discussed.
A disk B = {z ∈ C; |z − c| < R} has the exponential transform The truncated exponential transforms is: k,ℓ=0 identify the monic factor P (z) = z − c in the denominator P (z)P (w) of the rational approximant of the full exponential transform. Finally, as in the one dimensional diagonal Padé approximation scheme, one finds: The support of the measure µ contains full circles, hence it is not finite. The complex monomials are orthogonal in L 2 (µ), with norms 1 γ 2 k = z k 2 2,µ = r 2k dρ(r), k ≥ 0.
The orthonormal polynomials are so that the multiplier S = M z acts as a weighted shift: are the moments of a positive measure defined on [0, 1] (specifically dρ( √ r)), and All other entries in the associated Hessenberg matrix H are equal to zero. In short, H is a subnormal weighted shift. The spectrum of H is well understood: it coincides with the closed unit disk. Moreover, the elements of P 2 (µ) are analytic functions in the open disk, subject to a growth condition imposed by the coefficients γ k , cf. [27]. In particular, P 2 (µ) = L 2 (µ) and every λ ∈ D is a bounded point evaluation for P 2 (µ). That is, the cloud of µ is the full disk: Σ(µ) = D.
According to Theorem 3.6, the rate of convergence of the approximation scheme is dictated by the remainder: In general, for Hausdorff moment sequences such as ( 1 γ 2 k ) ∞ k=0 above, the convergence rate of consecutive quotients is known: A multi-fractal gauge, known as the local dimension of a measure, quantifies this asymptotics, cf. [10,23].
The remarkable feature of this class of examples is that all (normalized) rotationally invariant measures share the same cloud.

A uniform mass cloud plus finitely many point masses.
Let Ω ⊂ C be a bounded, connected and simply connected domain with smooth boundary. The uniform area mass ν = χ Ω dA distributed on Ω offers one of the best understood asymptotic analysis of complex orthogonal polynomials, with a century old history, see for instance [30]. In this case P 2 (ν) coincides with the Bergman space L 2 a (Ω), that is the collection of all analytic functions in Ω, square summable with respect to area.The reproducing kernel converges in Ω × Ω to an analytic/anti-analytic positive definite kernel, known as the Bergman kernel of Ω. If φ : Ω −→ D denotes a conformal mapping, then Therefore, the integral kernel representing (I − P )M * P is precisely The asymptotics of the singular numbers of this integral operator (known as a big Hankel operator) were thoroughly studied [?, 16]. In general, one carries to the unit disk all computations, via the inverse conformal mapping φ = ψ −1 . The integral kernel gives rise to a unitarily equivalent integral operator to (I − P )M * P , this time acting on L 2 (D, dA). Within this framework, a theorem due to Nowak [22] asserts that the singular numbers κ j of the Hankel operator (I − P )M * P associated to a domain Ω with smooth boundary satisfy κ j = O( 1 j ). Consequently the best approximation by a sequence of finite rank projections Q n , rank Q n = n, yields the error κ 2 n+1 + κ 2 n+2 + . . ., that is We show that this estimate is sharp for an ellipse, with respect to the polynomial filtration. More precisely, let E be the ellipse whose complement is described by Joukowski's map and parameter c > 0. The associated Hessenberg matrix is three diagonal, due to the special structure of the orthogonal polynomials with respect to area measure on E. Indeed, where ρ = e 2c and U j denotes Chebyshev polynomial of the second kind, see for details [21] pg. 259. We adopt the convention p −1 = 0 as well U −1 = 0. The three term recurrence relation for Chebyshev polynomials In other terms, h jj = 0, In order to evaluate the approximation rate in Theorem (3.6) we remark that j≤n<N <k |h jk | 2 = 0 whenever N > n, hence dictates, up to the stated constants, the rate of convergence. Since Area E = 1 4 (ρ−ρ −1 ) with ρ > 1, we find the error estimate: Conformal and quasi-conformal mapping techniques led recently to sharp estimates of the Hessenberg matrix entries of the subnormal multiplier S ν = M z , acting on P 2 (ν), [4]. This article represents the highest point of several decades of accumulated studies, by many authors. We mention the main setting. Let . . denote the conformal mapping of the exterior of the closed unit disk onto C \ Ω, It is customary to normalize F by the condition c −1 > 0. Theorem 1.2 of [4] asserts that there exists a constant β ≥ 1 so that Hessenberg matrix of S ν is asymptotically close to the Toepliz matrix T with entries c −1 , c 0 , c 1 , . . .. More specifically: In particular, the only non-zero under-diagonal terms satisfy: The value of the constant β depends on the regularity of the boundary of Ω, [4].
Remark that the series ℓ≥1 ℓ|c ℓ | 2 converges via a well known area estimate, in its turn a consequence of Stokes formula: Just for validation: this is nothing else than the trace of the self-commutator of the corresponding Toeplitz matrix T . Recall from Theorem 3.6 that the orthogonal projection P n onto C n [z] satisfies the identity In conclusion, lim n j≤n<k The yet unknown rate of convergence in the latter limit is not expected to be better than O( 1 n ), as the ellipse case shows. One of the first study of estimates of orthogonal polynomials in Bergman space setting is due to Carleman, in the case of real analytic boundaries, see for instance [30]. Without entering into details, we mention that in this scenario there exists ρ > 1, depending on the geometry of ∂Ω (how far Schwarz function of this curve analytically extends inside Ω), such that the decay in formula (3.5) is geometric: where O depends on n.
In general, Corollary 3.7 allows to estimate the error term in the moment approximation formula (3.6) even after adding finitely point masses outside the polynomial convex hull of Ω.

5.4.
Finite rank self-commutator. Even the simplest, finite rank self-commutator scenario raises challenging approximation theory questions. In view of the structural theorem of McCarthy and Yang, we have to focus in this case on a rational conformal map r : D −→ Ω of the disk onto a bounded quadrature domain, and the push forward measure µ = r * (dθ + ν), where dθ is arc length on the unit circle and ν is a finite atomic, positive measure supported by D. Then the operator S µ = P M z P is a cyclic, irreducible subnormal with finite rank self-commutator, and vice-versa [19,20]. Let the Schmidt expansion of (I − P )M * P be: with d finite. Every eigenfunction f j , 0 ≤ j ≤ d, of the self-commutator [S * , S] annihilates a finite codimension ideal of the ring of analytic functions defined on Ω (see [19]). Let a 1 , . . . , a p ∈ Ω denote the support of this ideal. Therefore every f j is a linear combination of the corresponding point evaluation functionals k a 1 , . . . , k ap . Due to the definition of the measure µ, these point evaluation functionals are push forward via r of point evaluations functionals with respect to the measure dθ + ν. And the latter can be explicitly computed recursively (via the so-called Uvarov's transform). We indicate only one step of this transform, corresponding to the measure dθ + δ α with α ∈ D: where C is a constant. More details can be found in [28]. We infer from these formulas that every evaluation functional K dθ+ν (z, w) analytically extends as a function of z, across ∂D, for a fixed value of w ∈ D. The same analytic continuation feature carries to k a 1 , . . . , k ap provided the boundary of Ω is smooth. Recall that the only singular points in the boundary of a quadrature domain are inner cusps [14]. We discuss a generic situation.
Proposition 5.1. Let S be a cyclic subnormal operator with finite-rank self-commutator, so that its spectrum is the closure of a finite union Ω of quadrature domains, plus a finite number of points. If the boundary of Ω is smooth, then there exists ρ > 1, so that the error in Theorem 3.6 satisfies: (I − P )M * (P − P n ) HS = O( 1 ρ n ), n → ∞.
U z n = z n+1 , n ≥ 0, with U * z n+1 = z n and U 1 = 0. One computes immediately the self commutator, in this Hilbert space representation: remarking that every commutator in the last expression annihilates z k , k ≥ 2. Therefore the two eigenfunctions f 0 , f 1 of [S * µ , S µ ] are linear combinations of 1 and z. In view of the preceding proof, the rate of convergence in Theorem 3.6 is a factor of δ n = inf{ z − h((1 − z) 2 ) 2,dθ , h ∈ C n [u]}, clearly equivalent to the best polynomial approximation (in the corresponding Lebesgue space norm) of the function √ w on the boundary of Ω. Remark that the function √ w is analytic and continuous on the closure Ω of the cardiodid, but it is not analytic in a neighborhood of it. The same theorem of Russell and Walsh and a refined Bernstein-Markov Inequality imply that δ n converges to zero at a slower than geometric rate. We do not expand here the technical details related to the cusp singularity in the boundary, see for instance [3]. The exact asymptotic decay of δ n remains unknown to us. It is worth mentioning that a slight perturbation in the definition of the measure µ, for instance µ ǫ = (r ǫ ) * dθ, r ǫ (z) = (z − 1 − ǫ) 2 , ǫ > 0 remains within the class of rank-two self-commutator, this time with a cloud possessing a smooth boundary. 5.5. Non-smooth clouds. Given the full generality of our approach, a vast array of pathologies enters into the picture. We simply make aware the reader of such pitfalls on one of the simplest examples, stressing that the cloud of a measure is only contained in its closed support.
Let Γ be a Jordan curve in the complex plane, with positive area. As for instance constructed by Osgood [24]. Let Ω be the interior of Γ. A celebrated theorem of Carleman, see for instance [30], asserts that complex polynomials are dense in the associated Bergman space. In our notation L 2 a (Ω) = P 2 (χ Ω dA). Denote µ = χ Ω dA and S µ the associated subnormal operator, equal to the multiplication by z on this space. The set of bounded point evaluations for µ is equal to Ω, also equal to the non-essential spectrum of S µ . The spectrum of S µ is equal to Ω. We know that the principal function g of S µ is equal to the characteristic function of Ω, modulo area null-sets. Therefore the cloud Σ(µ) is equal to Ω, and Area(Σ(µ) \ Ω) > 0, while µ(Σ(µ) \ Ω) = 0.