Communications in Mathematical Physics Large Deviations for Gibbs Measures with Singular Hamiltonians and Emergence of Kähler – Einstein Metrics

In the present paper and the companion paper (Berman, Kähler–Einstein metrics, canonical random point processes and birational geometry. arXiv:1307.3634, 2015) a probabilistic (statistical-mechanical) approach to the construction of canonical metrics on complex algebraic varieties X is introduced by sampling “temperature deformed” determinantal point processes. The main new ingredient is a large deviation principle for Gibbs measures with singular Hamiltonians, which is proved in the present paper. As an application we show that the unique Kähler–Einstein metric with negative Ricci curvature on a canonically polarized algebraic manifold X emerges in the many particle limit of the canonical point processes on X. In the companion paper (Berman in 2015) the extension to algebraic varieties X with positive Kodaira dimension is given and a conjectural picture relating negative temperature states to the existence problem for Kähler–Einstein metrics with positive Ricci curvature is developed.


Introduction
In the present paper and the companion paper [9] a probabilistic approach to the construction of canonical metrics on a complex algebraic varieties X is introduced, by sampling random point processes defined in terms of algebro-geometric data, canonically attached to X. The processes are "positive temperature deformations" of determinantal (fermionic) point processes and the main new ingredient is a large deviation principle for Gibbs measures with singular Hamiltonians which is proved in the present paper. As an application we show that the unique Kähler-Einstein metric with negative Ricci curvature on a canonically polarized algebraic manifold X emerges in the many particle limit of the canonical point processes on X. More generally, in the presence of a stress-energy tensor on X it is shown that the unique Kähler metric solving Einstein's equation on X with negative cosmological constant (in Euclidean signature) emerges in the many particle limit.
The generalization to the construction of canonical metrics and measures on a general algebraic variety X of positive Kodaira dimension are given in the companion paper [9], by exploiting the global pluripotential theory and variational calculus in [12,21,14,7]. This leads to a new probabilistic link between algebraic geometry on one hand (in particular the Minimal Model Program) and Kähler-Einstein geometry on the other. A conjectural picture is also developed describing the relation between the existence of negative temperature states and the existence problem for Kähler-Einstein metrics with positive Ricci curvature. In particular, relations to algebro-geometric stability properties, as in the Yau-Tian-Donaldson conjecture are described in [9]. See also [8,42] for connections to optimal transport in the real setting (corresponding to the case when X is toric and abelian variety, respectively) and [6] for connections to physics. A classical problem is to establish conditions for the existence of a macroscopic limit of the empirical measures δ N in the many particle limit N → ∞. More precisely, the problem is to show that the random measures δ N admit a deterministic limit µ β ∈ M 1 (X) in the sense that the law (1.2) Γ N := (δ N ) * µ (N ) β of δ N , defining a probability measure on M 1 (X), converges, as N → ∞, weakly to a Dirac mass concentrated on some µ β in M 1 (X). Equivalently, the marginals (µ (N ) weakly as probability measures on X j as N → ∞, which in the terminology of Kac and Snitzmann [53] means that the sequence µ (N ) β is chaotic. A stronger exponential notion of convergence of δ N , with an explicit speed and rate functional, is offered by the theory of large deviations, by demanding that the laws Γ N satisfy a Large Deviation Principle (LDP) with speed r N and a rate functional F, symbolically expressed as Γ N (µ) ∼ e −rN F (µ) , N → ∞ and assuming that F admits a unique minimizer µ β in M 1 (X). Loosely speaking this means that the probability of finding a cloud of N points x 1 , ..., x N on X such that the corresponding measure 1 N i δ xi approximates a volume form µ is exponentially small unless µ is the minimizer µ β of F β .
Our main general result establish such a LDP for a class of singular Hamiltonians: converges, as N → ∞, to a Gateaux differentiable functional F (u) on C 0 (X) • H (N ) is uniformly quasi-superharmonic, i.e. ∆ x1 H (N ) (x 1 , x 2 , ...x N ) ≤ C on X N Then, for any fixed β > 0, the measures (δ N ) * (e −βH (N ) dV ⊗N ) on M 1 (X) satisfy, as N → ∞, a large deviation principle (LDP) with speed βN and good rate functional where the functional E(µ) is the Legendre-Fenchel transform of −F (−·) and D dV (µ) is the entropy of µ relative to dV. In particular, the empirical measures δ N of the corresponding random point processes on X converge in law to the deterministic measure given by the unique minimizer µ β of F β . Moreover, if the equation (1.4) dF |u = e βu dV X e βu dV on C 0 (X) admits a solution u β , then the corresponding differential µ β := dF |u β is the minimizer of F β .
It follows from the previous theorem that the LDP indeed also holds for the corresponding Gibbs measures with the rate functional F β + C β , where C β is the following constant: It should be stressed that even the convergence of the first marginals of µ (N ) β , implied by the previous theorem, appears to be a new result.
As explained in Section 4.1 the asymptotics in the first assumption of the theorem may be replaced by the weaker assumption that there exists a functional E(µ) on M 1 (X) such that H (N ) (x 1 , ..., x N )/N → E(µ) in the sense of Gamma convergence. Moreover, Theorem 1.1 can be viewed as a generalization of the Gärtner-Ellis theorem in the setting of Gibbs measures (see Section4.2). Let us also point out that that the restriction that X be compact can be removed if suitable growth-assumptions of H (N ) "at infinity" are made. But since our main application concerns the case of compact complex manifolds, we have, for simplicity taken X to be compact. It may be illuminating to point out that in thermodynamical terms the content of Theorem 1.1 can be heuristically expressed as follows. Imagine that we know the macroscopic ground state (i.e. the state of zero energy E) of a system of a large number N of particles in thermal equilibrium at zero temperature (i.e. at β = ∞). If we can rule out any first order phase transitions at zero-temperature (which essentially means that the macroscopic equilibrium states is unique), then increasing the temperature (i.e decreasing β) leads to a new macroscopic equilibrium state, minimizing the corresponding free energy functional E − S/β, where S is the physical entropy (i.e. S = −D with our sign conventions). In fact, in the complex geometric setting to which we next turn. the zero-temperature limit β → ∞ is reminiscent of a (second order) gas-liquid phase transition [10].
1.2. Application to Kähler-Einstein geometry. Let nowX be an n−dimensional complex algebraic projective variety of positive Kodaira dimension. This means that the plurigenera N k of X are increasing: where H 0 (X, kK X ) denotes, as usual, the complex vector space of all pluricanonical (holomorphic) n−forms of X at level k, i.e. H 0 (X, kK X ) is the space of all global holomorphic sections of the k tensor power of the canonical line bundle K X := Λ n (T * X) of X (using additive notation of tensor powers). In terms of local holomorphic coordinates z 1 , ..., z n on X this simply means that the elements s (k) of H 0 (X, kK X ) may be represented by local holomorphic functions s (k) on X, such that |s (k) | 2/k transforms as a density on X and thus defines a measure on X. To any such algebraic variety X we can associate the following canonical sequence of probability measures µ (N k ) on X N k : totally antisymmetric (and thus defined up to a multiplicative complex number) and Z N k is the normalizing constant. The probability measure µ (N k ) thus defined is symmetric, i.e. invariant under the natural action of the permutation group S N k , independent of the choice of generator det S (k) and hence defines a canonical random point process on X with N k points. As shown in the companion paper [9], it follows from Theorem 1.1, combined with the asymptotics in [12] that the corresponding empirical measures δ N k converge in law, as k → ∞, towards a deterministic measure µ can on X, which is thus canonically attached to X. In fact, using the pluripotential theory and variational calculus in [14,7] the limiting measure µ can is shown to coincide with the canonical measure of Song-Tian [54] and Tsuji [55] previously defined in terms of Kähler-Einstein geometry or equivalently as solutions to certain complex Monge-Ampère equations. In the present paper we will show how to apply Theorem 1.1 in the special case when K X is positive (i.e. ample) to deduce the following Theorem 1.2. Let X be a compact complex manifold such with positive canonical line bundle K X . Then the empirical measures δ N k of the corresponding canonical random point processes on X converge in law, as N k → ∞, towards the normalized volume form dV KE of the unique Kähler-Einstein metric ω KE on X. More precisely, the law of δ N k satisfies a large deviation principle with speed N k whose rate functional may be identified with Mabuchi's K-energy functional on the space of Kähler metrics in c 1 (K X ).
By the celebrated Aubin-Yau theorem [3,60] (consisting of Kähler forms, for k sufficiently large) converges, as k → ∞, to the Kähler-Einstein metric ω KE in the weak topology of currents on X.
Theorem 1.2 fits into a more general setting of "temperature deformed" determinantal point processes attached to a polarized manifold (X, L), i.e. a compact complex manifolds X endowed with a positive line bundle L (Theorem 5.7). More precisely, in the general setting the point processes are attached to the data ( · , dV, β k ) consisting of a Hermitian metric · on a L, a volume form dV on X and a sequence of positive numbers β k → β ∈]0, ∞]. Then the corresponding probability measures on X N k are defined by where det S (k) is a generator of the top exterior power Λ N k H 0 (X, kL). Concretely, the corresponding LDP is equivalent to the following asymptotics for the L 2β k /k −norm of the generator det S (k) of the determinant line of H 0 (X, kL) which is orthonormal with respect to the L 2 −product determined by ( · , dV ) : (by Lemma 4.7). In this general setting the limiting deterministic measure µ β minimizing F β is the volume form of the unique Kähler metric ω β in the first Chern class of L solving the twisted Kähler-Einstein equation where the twisting form η is explicitly determined by ( · , dV, β). The point is that when L = K X any given volume form dV naturally defines a metric · dV on L and the probability measures on X N k attached to ( · dV , dV, 1) are precisely the canonical ones defined by formula 1.6. Moreover, in this special case η vanishes and the equation 1.10 thus reduces to the the usual Kähler-Einstein equation 1.7. The more general twisted version of the equation has previously appeared in various situations in Kähler geometry [38,54,55]. From the physics point of view the twisting form η corresponds to the (trace-reversed) stress-energy tensor in Einstein's equations on X (with Euclidean signature). The Hamiltonians corresponding to the probability measures 1.9 are strongly non-linear unless X is a Riemann surface, i.e. unless n = 1. In fact, in the simplest latter case, i.e. when X is the Riemann sphere, where W is the Green function of the corresponding Laplace operator and then the corresponding functional E(µ) is the Dirichlet energy (Remark 5.11).
In general, the connection to the Kähler-Einstein geometry of (X, L) will be shown to arise from the fact that the equation 1.4 is intimately related to the complex Monge-Ampère equation where ω 0 is the normalized curvature two form of the given metric · on L. More precisely, the two equations coincide for smooth functions u such that ω 0 + i∂∂u is a Kähler form (i.e. smooth and positive). In this complex geometric setting the strong non-linearity of the Hamiltonians H (N ) when n ≥ 2 is reflected in the non-linearity of the complex Monge-Ampère operator appearing in the left hand side of equation 1.12 (coinciding with the Laplacian when n = 1). Furthermore, the singularity of H (N ) (which is present for any dimension n) is a reflection of the fact that solutions to the (generalized) Calabi-Yau equation are, in general, singular when µ is a probability measure on X (as is clear already for the Laplace equation appearing when n = 1). Finally, let us point out that the extension to general complex algebraic manifolds X with positive Kodaira dimension, established in the companion paper [9], relies on an extension of Theorem 5.7 to line bundles L, which are big (but not necessarily positive); see Section 5.4.

1.3.
Comparison with previous results. First a comment on relations to the physics literature: in the case n = 1 (i.e. in two real dimensions) the quasi-linear Laplace type equation 1.12 arises as the macroscopic equilibrium equation in a range of statistical mechanical models of mean field type: it is called the Joyce-Montgomery equation in Onsager's vortex model for 2D turbulence, the Poisson-Boltzmann equation in the Debye-Hückel theory of plasmas and electrolytes and the Lane-Emden equation in stellar physics (see [36]). But the Monge-Ampère equation (n > 1) does not seem to have a appeared in any statistical mechanical model before. On the other hand, in the case when β k := k the density of the corresponding probability measure has a natural quantum mechanical interpretation: it is the squared amplitude of the Slater determinant representing a maximally filled many particle state of N free fermions on X, subject to an exterior magnetic field (the corresponding single particle wave functions are elements of H 0 (X, kL) and represent the corresponding lowest Landau levels). The case when β k = 1 ν k, for a given positive integer ν, also appears in the fractional Quantum Hall Effect, where the corresponding probability density is the squared amplitude of the Laughlin state (see the review [46] and references therein).
1.3.1. Large deviations. The LDP in Theorem 1.1 in the case when H (N ) is uniformly equicontinuous is essentially well-known in the setting of mean field models [35,8] (it then also applies to the case of negative β, by replacing H (N ) with −H (N ) ). But the key feature of the previous theorem is that it applies to a large class of singular Hamiltonians and in particular H (N ) is allowed to be strongly repulsive in the sense that it blows up, as two points merge (and hence the Gibbs measure may be ill-defined when β is negative). It seems that the only previous class where a convergence result as in Theorem 4.6 has been established for singular Hamiltonians is in the "linear" case when H (N ) is a sum of pair interactions with a mean field scaling: where the pair interaction W is allowed to be singular along the diagonal, as long it is lower semi-continuous and in L 1 loc (this is indeed a mean field interaction in the sense that each particle x i is exposed to the average of the pair interactions W (x i , x j ) for the N − 1 remaining particles). Then the asymptotics of the partitions functions 1.5 can be obtained using the method of Messer-Spohn [49], which is based on the Gibbs variational principle and which crucially relies on the the existence of the mean energyĒ(µ) corresponding to H (N ) (see [43,27] for the case of a logarithmic singularity which is motivated by Onsager's vortex model for 2D turbulence [50,36]). A similar argument applies in the case of "finite order", i.e. when H (N ) is a sum of j−point interactions for a uniformly bounded j (then E(µ) depends polynomially on µ). However, the main point of the previous theorem is to avoid the latter assumption which is not satisfied in the application to Kähler-Einstein geometry (apart from the classical lowest dimensional setting of Riemann surfaces). In particular, the present proof bypasses the problem of the existence of the limiting mean energies. Instead the main idea of the proof is to exploit the Riemannian orbifold geometry of the space of configurations of N points on X, viewed as the singular quotients X N /S N , where S N is the symmetric group acting on X N by permuting the factors. The key result is a submean inequality for positive quasi-subharmonic functions on X N /S N with a distortion coefficient which is sub-exponential in the dimension (Theorem 2.1), which is closely related to an inequality of Li-Schoen [48].
There is also another approach to large deviation principles for mean field Hamiltonians of the form 1.14 originating in the literature on random matrices and Coulomb gases [18,19,29,51], which as explained in [51], is closely related to the notion of Gamma convergence (see also [57,58] for applications to univariat random polynomials). This approach seems to be limited to the case when β N ≫ log N and in particular β = ∞ so that the entropy contributions can be neglected. 1 . See also [34] for a general LDP for Hamiltonians of the form 1.14 using weak convergence methods.
Let us also point out that the role of (det S (k) )(x 1 , x 2 , ...x N k ) appearing in formula 1.9 is played by the classical Vandermonde determinant in the random matrix literature (for Example 5.6. In fact, there is a non-compact analogue of Theorem 5.7 in Euclidean C n which specializes to the setting of random matrix theory and the 2D log gas when n = 1 and β = ∞ and to the 2D vortex model (for n = 1 and β < ∞) and which can be proved by supplementing the proof of Theorem 5.7 with a tightness estimate, as in the non-compact setting considered for β = ∞ in [5] (see also [20] for the case β = ∞). Details will appear elsewhere.
1.3.2. Kähler geometry. A statistical mechanics approach has previously been applied to conformal geometry [44], as opposed to the present complex-geometric setting. The role of the "determinantal" Hamiltonian 1.11 is in the conformal setting played by a mean field Hamiltonian of the form 1.14 with a logarithmic pair interaction and the role of the fully non-linear complex Monge-Ampère operator is played by a linear conformally invariant operator, which is zero-order perturbation of a power of the Laplacian (the Paneitz operator). Accordingly previous results in [43,27] concerning such Hamiltonians can be applied in the conformal setting (compare the discussion above) in the conformal setting, while the present setting seems to require new methods.
The present probabilistic should be viewed in the light of the pervasive philosophy in Kähler geometry, going back to Yau [61], of approximating metrics on a complex algebraic manifold with algebraically defined Bergman metrics, which may me identified with elements of the symmetric space GL(N, C)/U (N ). For example, the quasi-explicit Kähler metrics ω k in formula 1.8, approximating the Kähler-Einstein metric ω KE on a canonically polarized manifold X, are analogs of Donaldson's balanced metrics in GL(N, C)/U (N ) [33]. One advantage of the present approach is that, as shown in the companion paper [9], the approximation also applies when, for example, X is of general type, where the role of ω KE is played by the the canonical Kähler-Einstein current on X (which is singular along a subvariety of X) [21,14]. In another direction it would be interesting to see if the present approach can be implemented to construct numerical simulations of Kähler-Einstein metrics, using Monte Carlo type methods, complementing the different numerical approaches in [33,32] (see [4] for relations between Monte Carlo simulations and similar polynomial determinantal point processes).
Even if the connection between canonical random point processes on a complex algebraic manifold X does not seem to have been studied before, there are some connections to previous work on random polynomials/holomorphic sections in a given back-ground geometry [52]; in particular in the one-dimensional setting where an LDP was obtained in [57,58]. Another probabilistic approach to the space of Kähler metrics has been introduced in a a series of papers by Ferrari, Klevtsov and Zelditch [37], motivated by Quantum Field Theory. The approach aims at approximating random random Kähler metrics with random Bergman metrics. Accordingly, the role of the N −particle space X N /S N is in [37] played by the symmetric space GL(N.C)/U (N ). In conclusion, it would be very interesting to understand the precise connections between [37] and the present setting, as well as the connection to Donaldson's balanced metrics [33].
Acknowledgment. It is a pleasure to thank Sebastien Boucksom, David Witt-Nyström, Vincent Guedj and Ahmed Zeriahi for the stimulating collaborations [12,13,14], which paved the way for the present work. I am also grateful to Bo Berndtsson for infinitely many fruitful discussions on complex analysis and Kähler geometry over the years. Thanks, in particular, to Sebastien Boucksom for illuminating discussions on Lemma 3.6. The present paper, together with the companion paper [9], supersedes the first arXiv version of the paper [9]. This work was supported by grants from the ERC and the KAW foundation.
Organization. In section 2 we prove the submean inequality in large dimensions, which plays a key role in the subsequent section 3 where the general LDP in Theorem 1.1 is proved. In Section 4 we make a digression on relations to previous methods and notions used in the literature on large deviations. The applications to Kähler-Einstein geometry are given and Section 5. For the convenience of readers lacking background in Kähler geometry we start the section by giving a reasonably self-contained account of the Kähler geometry setup (including some rudiments of pluripotential theory). The article is concluded with an outlook in Section 6 on some open problems and an appendix where the dimension dependence on the constant in the Cheng-Yau gradient estimate is obtained, by tracing through the usual proof.

Submean inequalities in large dimension
2.1. Setup. Let (X, g) be a n−dimensional Riemannian manifold and assume that Ric g ≥ −κ 2 (n − 1)g for some positive constant κ (sometimes referred to as the normalized lower bound on the Ricci curvature). Let G a finite group acting by isometries on X and denote by M := X/G the corresponding quotient equipped with the distance function induced by the metric g, i.e.
where d X is the Riemannian distance function on (X, g). Even though the quotient M is not a manifold in general (since G will in general have fixed points) it still comes with a smooth structure in the following sense. Denote by p the natural projection map from X N to M. Using the projection p we can identify a function f on M with G−invariant function p * f on X and accordingly we say that f is smooth if p * f is. Similarly, there is a natural notion of Laplacian ∆on the quotient M : the Laplacian ∆u of a locally integrable function u on M is the signed Radon measure defined byˆM for any smooth function f on M. More generally, by localization, this setup naturally extends to the setting of Riemannian orbifolds (see [22]), but the present setting of global quotients will be adequate for our purposes.

Statement of the submean inequality.
Theorem 2.1. Let (X, g) be a Riemannian manifold of dimension n such that Ric g ≥ −κ 2 (n − 1)g and G a finite group acting by isometries on X. Denote by M := X/G the corresponding quotient equipped with the distance function induced by the metric g and let v be a non-negative function on M such that ∆ g v ≥ −λ 2 v for some non-negative constant λ. Then, for any δ > 0 and ǫ ∈]0, 1] there exist constants A and C such that where C only depends on an upper bound on κ and A only depends on δ and ǫ (assuming that the balls above are contained in a compact subset of M ).
Note that by the G−invariance we may as well replace the functional v and the balls on M with their pull-back to X.

2.3.
Proof of the submean inequality in Theorem 2.1. We will follow closely the elegant proof of Li-Schoen [48] of a similar submean inequality. But there are two new features here that we have to deal with: • We have to make explicit the dependence on the dimension n of all constants and make sure that the final contribution is sub-exponential in n • We have to adapt the results to the singular setting of a Riemannian quotient Before turning to the proof we point out that it is well-known that submean inequalities with a multiplicative constant C(n) do hold in the more general singular setting of Alexandrov spaces (with a strict lower bound −κ on the sectional curvature). But it seems that the current proofs (see for example [40]), which combine local Poincaré and Sobolev inequalities with the Moser iteration technique, do not give the subexponential dependence on C(n) that we need.
We recall that the two main ingredients in the proof of the result of Li-Schoen referred to above is the gradient estimate of Cheng-Yau [30] and a Poincaré-Dirichlet inequality on balls. Let us start with the gradient estimate that we will need: for some absolute constant C (in particular, independent of n, κ and a).
Proof. In the smooth case this is the celebrated Cheng-Yau gradient estimate [30]. The result is usually stated without an explicit estimate of the multiplicative constant C n in terms of n, but tracing through the proof in [30] gives C n ≤ Cn (see the appendix in the present paper and also [2] for a probabilistic proof providing an explicit constant). We claim that the same estimate holds in the present setting using a lifting argument. To see this recall that the usual proof of the gradient estimate proceeds as follows (see the appendix). Set φ(x) := |∇ log u| (= |∇u| /u) and F (x) := φ(x)(ρ x0 − a) 2 . Then F attains its maximum in a point x 1 in the interior of B a (x 0 ) (otherwise |∇u| vanishes identically and then we are trivially done). Hence, F (x) ≤ F (x 1 ) on some neighborhood U of x 1 . Now, in case F (or equivalently ρ x0 ) is smooth on U we get ∆F ≤ 0 and ∇F = 0 at x 1 . Calculating ∆F and using Bochner formula and Laplacian comparison then gives which is the desired estimate. In the case when ρ x0 is not smooth on U, i.e. x 1 is contained in the cut locus of x 0 one first replaces ρ x0 with a smooth approximation ρ (ǫ) x0 of ρ x0 (which is a local barrier for ρ x0 ) and then lets ǫ → 0 to get the same conclusion as before. In the singular case M = X/G we proceed as follows. First we identify F with a G−invariant function on the inverse image of B R (x 0 ) in X (and x 0 and x 1 with a choice of lifts in the corresponding G−orbits) and set F : . By definitionρ x0 ≥ ρ x0 on X and, after possibly changing the lift of the point x 1 we may assume thatρ x0 = ρ x0 at x = x 1 and henceρ x0 < a (after perhaps shrinking U ). In particular,F ≤ F on U andF = F at x 1 and henceF also has a local maximum at x 1 . But then the previous argument in the smooth case gives that 2.1 holds with ρ x0 replaced bỹ ρ x0 . But since the two functions agree at x 1 this concludes the proof in the general case.
Proof. Set v := log h and fix x ∈ B ǫδ (x 0 ). Integrating along a minimizing geodesic connecting x 0 and x and using the gradient estimate in the previous proposition gives In particular, for any two points which implies the proposition after renaming the constant C.
The second key ingredient in the proof of Theorem 2.1 is the following Poincaré-Dirichlet inequality: Proposition. Let f be a smooth function on B δ (x 0 ) vanishing on the boundary. ThenˆB where the constant C only depends on an upper bound on κ.
Proof. We follow the proof in [48] with one crucial modification (compare the remark below). To fix ideas we first consider the case of a Riemannian manifold. Fix a point p in the boundary of the ball B 1 (x 0 ) and denote by r 1 (x) the distance between x ∈ M and p. From the standard comparison estimate for the Laplacian we get (in the weak sense and point-wise away from the cut locus of p). In particular, for any positive number a we deduce the following inequality on B δ (x 0 ) (using that g(∇r 1 , ∇r 1 ) = 1) a.e.) Hence, setting a := n( 1 Multiplying by |f | and integrating once by parts (and using that ∇r Estimating e −ar1 ≤ e −a(1−δ) in the left hand side above and rearranging giveŝ |f |dV, (using that g(∇r 1 , ∇r 1 ) ≤ 1 in the sense of upper gradients). This shows that the L 1 −version of the Poincaré inequality in question holds with the constant ( 1 , which for δ sufficiently small is bounded from above by e n(4+2κ)δ . The general Riemannian L 2 −Poincare inequality now follows from replacing |f | with |f | 2 and using Hölder's inequality. Finally, in the case of the a Riemannian quotient M we can proceed exactly as above using that the Laplacian comparison estimate in formula 2.2 is still valid. Indeed, the pull-back p * r 1 of r 1 to X is an infimum of functions for which the corresponding estimate holds (by the usual Laplacian comparison estimate and the assumption that G acts by isometries). But then the estimate also holds for the function p * r 1 , by basic properties of Laplacians. More generally, the required Laplacian comparison estimate was shown in [22] for general Riemannian orbifolds.
Remark 2.4. The only difference from the argument used in [48] is that we have taken the point p to be of distance 1 from x 0 rather than distance 2δ, as used in [48]. For δ small this change has the effect of improving the exponential factor from e n(1+δκ) to e n(δ+δκ) , which is crucial as we need a constant in the Poincare inequality which has subexponential growth in n as δ → 0.
Next, by the triangle inequalitŷ But h is the solution to a Dirichlet problem and as such minimizes the Dirichlet norm´B δ(x 0 ) |∇h| 2 over all subharmonic functions with the same boundary values as h. Accordingly,ˆB (as is seen by multiplying with a suitable smooth function χ supported on B 2δ such that χ = 1 on B δ ). All in all this concludes the proof of Theorem 2.1 in the case λ = 0.
Finally, to handle the general case (i.e. λ = 0) we set N := M ×] − 1, 1[ equipped with the standard product metric and apply the previous case to the function ve λt to get But restricting the sup in the left hand side to B ǫδ (x 0 ) × {0} and using that which concludes the proof of the general case (after a suitable rescaling).

Proof of the large deviation principle for Gibbs measures
Given a compact topological space X we will denote by C 0 (X) the space of all continuous functions u on X, equipped with the sup-norm and by M(X) the space of all signed (Borel) measures on X. The subset of M(X) consisting of all probability measures will be denoted by M 1 (X). We endow M(X) with the weak topology, i.e. µ j is said to converge to µ weakly in M(X) if µ j , u j → µ, u for any continuous function u on X, i.e. for any u ∈ C 0 (X), where u, µ denotes the standard integration pairing between C 0 (X) and M(X) (equivalently, the weak topology is precisely the weak*-topology when M(X) is identified with the topological dual of C 0 (X)). A functional F on C 0 (X) will be said to be Gateaux differentiable if it is differentiable along affine lines and for any u in C 0 (X) there exists an element dF |u in M(X), called the differential of F at u, such that for any A random point process with N particles is by definition a probability measure µ (N ) on the N −particle space X N which is symmetric, i.e. invariant under permutations of the factors of X N . The empirical measure of a given random point process is the following random measure ). By definition the law of δ N is the push-forward of µ (N ) to M 1 (X) under the map δ N , which thus defines a probability measure on M 1 (X). Now fix a back-ground measure µ 0 on X and let H (N ) be a given N −particle Hamiltonian, i.e. a symmetric function on X N , which we will assume is lower semi-continuous (and in particular bounded from below, since X is assumed compact). Also fixing a positive number β the corresponding Gibbs measure (at inverse temperature β) is the symmetric probability measure on X N defined as is called the (N −particle) partition function. In our setting we will take µ 0 to be the volume form dV of a fixed Riemannian metric. Given a continuous function u on X we will also write where u has been identified with the following function on the product X N : (ii) A sequence Γ k of measures on P satisfies a large deviation principle with speed r k and rate function I if lim sup

Preliminaries on Large Deviation Principles and
for any open subset G of P.
Remark 3.2. The LDP is said to be weak if the upper bound is only assumed to hold when F is compact. Anyway, we will be mainly interested in the case when P is compact and hence the notion of a weak LDP and an LDP then coincide (and moreover any rate functional is automatically good).
We will be mainly interested in the case when Γ k is a probability measure (which implies that I ≥ 0 with infimum equal to 0). Then it will be convenient to use the following alternative formulation of a LDP (see Theorems 4.1.11 and 4.1.18 in [31]): Proposition 3.3. P be a metric space and denote by B ǫ (µ) the ball of radius ǫ centered at µ ∈ P. Then a sequence Γ N of probability measures on P satisfies a weak LDP with speed r N and a rate functional I iff We note the following simple lemma which allows one to extend the previous proposition to the non-normalized measures (δ N ) * e −βH (N ) µ ⊗N

Proof. SetΓ
By assumption C N,β is uniformly bounded and we denote by C β a given limit point of the sequence obtained by replacing N with a subsequence N j .
showing that C β is independent of the subsequence. Hence, the whole sequence converges towards C β , which proves one direction in the Lemma. The converse is proved in a similar way.
We will also use the following classical result of Sanov, which is the standard example of a LDP for point processes [31] (the result follow, for example, from the Gärtner-Ellis theorem; see Section 4.2). We recall that the relative entropy D µ0 (also called the Kullback-Leibler divergence or the information divergence in probability and information theory) is the functional on M 1 (X) defined by when µ has a density µ µ0 with respect to µ 0 and otherwise D µ0 (µ) := ∞. When µ 0 is a probability measure,D µ0 (µ) ≥ 0 and D µ0 (µ) = 0 iff µ = µ 0 (by Jensen's inequality). in terms of the canonical pairing between V and V * . In the present setting we will take V = C 0 (X) and V * = M(X), the space of all signed Borel measures on a compact topological space X. We will use the following variant of the Brøndsted-Rockafellar property A * [26]: Lemma 3.6. Let f be function on C 0 (X) which is Gateaux differentiable. Then, for any µ ∈ M(X) such that f * (µ) < ∞ there exists a sequence of u j ∈ C 0 (X) such that Proof. First recall that a convex function g on a topological vector space E is said to be subdifferentiable at x ∈ E if g(x) < ∞ and g admits a subgradient x * at x, i.e. an element x * in the topological dual E * such that for any y ∈ E g(y) ≥ g(x) + (y − x), x * The set of all such subgradients is denoted by (∂g)(x). Now assume that g = f * for a convex function f on a Banach space V. Then g is a lower semi-continuous function convex function on the topological vector space E := V * equipped with its weak topology. According to [26,Thm 2] any element µ ∈ E * such that f * (µ) < ∞ has the property that there exists a sequence µ j → µ in V * such that f * is subdifferentiable at µ j with a subgradient in E. Equivalently, this means that there exists u j ∈ E such that µ j ∈ (∂f )(u j ) (as follows from the definition the Legendre-Fenchel transform). Finally, setting E := C 0 (X) and observing that if f is Gateaux differentiable at u ∈ V, then (∂f )(u) = {df |u } (as is seen by restricting f to any affine line) thus concludes the proof.
Remark 3.7. By convexity, if µ = df |u for some u ∈ V := C 0 (X), then f * (µ) = u, df |u −f (u), which is essentially the classical definition of the Legendre transform of f at µ. Accordingly, the previous lemma may be reformulated as the statement that the Legendre-Fenchel transform is the greatest lower semi-continuous extension to all of V of the Legendre transform, originally defined on (df )(V ) ⊂ V * .
3.3. The proof of Theorem 5.7. We start with the following simple Lemma 3.8. Assume that H (N ) satisfies the quasi-superharmonicity assumption in the second point of Theorem 1.1. Then, for any sequence of positive numbers Proof. The inequality ≤ is trivial and to prove the reversed inequality we fix a sequence of x (N ) ∈ X N realizing the infimum appearing the right hand side above. Then replacing the integral of X N with an integral over the L ∞ −ball for a fixed number ǫ and a fixed metric g with distance function d g , and using the classical submean inequality in each variable with a fixed multiplicative constant C giveŝ To handle the case when β N = β + o(1) for a finite β we will need to use the subexponential dependence on the dimensions of the multiplicative constant appearing in Theorem 2.1. To this end we first recall that, since X is assumed compact, the weak topology on M 1 (X) is metrized by the Wasserstein 2-metric d induced by a given Riemannian metric g on X, where d(µ, ν) 2 := inf Γ∈Γ(µ,ν)ˆd where Γ(µ, ν) is the space of all couplings between µ and ν, i.e. all probability measures Γ on X × X such that the push forward of Γ to the first and second factor is equal to µ and ν respectively. Proposition 3.9. For any given ǫ > 0 there exists a positive constant C ǫ such that the following submean inequality holds on X N , for any N : where B r (x (N ) ) denotes the inverse image in X N , under the map δ N , of the Wasserstein ball of radius r centered at δ N (x (N ) ) Proof. First observe that the pull-back of d on M 1 (X) to the quotient space X (N ) := X N /S N under the map δ N defined by the empirical measure (formula 3.1) coincides with 1/N 1/2 times the quotient distance function on X (N ) , induced by the product Riemannian metric on X N : Indeed, this is well-known and follows from the Birkhoff-Von Neumann theorem which gives that for any symmetric function c(x, y) on X × X we have that if where S N is the symmetric group on N letters. Now consider the metric space (X (N ) , d (N ) ) which is the quotient space defined with respect to the finite group S N acting isometrically on the Riemannian manifold (X N , g N ), where g N denotes 1/N times the product Riemannian metric. By assumption H (N ) is S N −invariant and ∆ gN H (N ) ≤ C on X N (using the obvious scaling of the Laplacian). Moreover, since X is compact there exists a non-negative number k such that Ric g ≥ −kg on X and hence rescaling gives Ric g N ≥ −kN g N on (X N , g N ).
But the dimension of X N is equal to nN and hence setting κ 2 := k/n + 1 shows that, for N large, the assumptions in Theorem 2.1 are satisfied for u := e −βH (N ) and (X, g) replaced by (X N , g N ). Applying the latter theorem with δ = ǫ and using the pull-back property in formula 3.7 then shows that the submean property 3.6 indeed holds.
We will also rely on the following simple but very useful lemma (which was used in the similar context of Fekete points in [13]). Lemma 3.10. Fix u * ∈ C 0 (X) and assume that x (N ) * ∈ X N is a minimizer of the function (H (N ) + u * )/N on X N . If the corresponding large N − limit F (u) exists for all u ∈ C 0 (X) and F is Gateaux differentiable at u * , then δ N (x Note that f is a concave function in t (since it is defined as an inf of affine functions) and f N (t) is affine in t. But then it follows from the differentiability of f at t = 0 which thus concludes the proof of the lemma.
The upper bound in the LDP. By Lemma 3.4 it will be enough to establish the LDP for the non-normalized measures Γ N in formula 3.3. To prove the upper bound of the integrals appearing in the equivalent formulation of the LDP in Prop 3.3 we fix a function u ∈ C 0 (X) and rewrite Then, trivially, for any fixed ǫ > 0, Hence, replacing the sup over B ǫ (µ) with the sup over all of X N k and applying Sanov's theorem relative to the tilted volume form µ u gives lim ǫ→0 lim sup using that D e βu dV (µ) = −β´uµ + D dV (µ). Finally, taking the infimum over all u ∈ C 0 (X) shows that the lim sup in the previous formula is bounded from above by −F (µ), Remark 3.11. In the argument above dV can be replaced by any finite measure µ 0 on X.
The lower bound in the LDP. As usually the proof of the lower bound in the LDP is the hardest. We first assume that µ = dF |u for some u ∈ C 0 (X). Denote by x (N ) ∈ X N a sequence of minimizers of H (N ) + u. By Lemma 3.10 we have that δ(x (N ) ) → µ weakly and hence for N sufficiently large where δ(ǫ) is a modulus of continuity for u on X tending to zero with ǫ (by the compactness of X). Next, applying the submean inequality 3.6 gives Since δ(x (N ) → µ we may, for N sufficiently large, assume that B ǫ 2 /2 (µ) ⊂ B ǫ 2 (δ(x (N ) )) and hence letting N → ∞ and using Sanov's theorem (i.e. Prop 3.5) for ǫ fixed and the assumed convergence of (H (N ) + u)(x (N ) )/N gives Since µ is a candidate for the inf in the right hand side above the inf in question may be estimated from above by D dV (µ) and hence letting ǫ → 0 concludes the proof under the assumption that µ := dF |u for some u ∈ C 0 (X). To prove the general case we invoke Lemma 3.6 to write µ as a weak limit of µ j := dF |uj for u j ∈ C 0 (X). We may then replace u in the previous argument with u j for a fixed j and replace µ with µ j in the previous argument to get, for j ≥ j δ , lim inf N →∞ But for j sufficiently large µ j is in the ball B ǫ 2 /2 (ν) and hence the inf above is bounded from above by D dV (µ) giving Letting first ǫ → 0 and then j → ∞ gives Finally, by Lemma 3.6 we may assume that E(µ j ) → E(µ) and that concludes the proof.
The equation for the minimizer µ β . Finally, the equation 1.4 follows immediately from the following general convex analytical result: Lemma 3.12. Let X be a compact topological space and f and g be Gateaux differentiable convex functionals on C 0 (X) such that the differentials dg and df takes values in M 1 (X). Then • The following identity holds: • if the sup in the right hand side above is attained at some u 0 in C 0 (X) (i.e. if −f (−u)− g(u) admits a critical point u 0 ), then, setting F (u) := −f (−u), the measure µ 0 := dF |u0 minimizes the functional f * + g * on M 1 (X).
Proof. First observe that f and g are Lipschitz continuous on the Banach space C 0 (X). Indeed, setting u t : and similarly for g. The first point in the lemma is then obtained as a special case of the Fenchel-Rockafeller duality theorem which only requires that f and g be convex on a Banach space V and that f and g be finite at some point u where f is moreover assumed continuous [25, Thm 1.12]. To prove the second point we let u 0 be a critical point of F (u) − g(u) on C 0 (X), i.e.
(3.10) dF |u0 = dg |u0 , which, by convexity, means that u 0 realizes the sup in the right hand side of formula 3.9. We rewrite, which concludes the proof, using the first point.

Relations to Γ−convergence, the Gärtner-Ellis theorem and mean energy
Before turning to the applications of Theorem 1.1 in the complex geometric setting we explore some relations to previous results and methods in the literature.

4.1.
Relations to Gamma convergence. We recall that a sequence of functions E N on a topological space P is said to Γ−converge to a function E on P if (such a sequence µ N is called a recovery sequence); see [24]. It then follows that E is lower semi-continuous on P. In the present setting we take, as before, P = M(X) Proof. Using the characterization of a LDP in Proposition 3.3, the upper bound in the LDP follows almost immediately from the liminf property of the Gammaconvergence together with Sanov's theorem. To prove the lower bound fix µ ∈ M 1 (X) and take a recovery sequence µ N corresponding to a sequence x (N ) ∈ X N . Then, using the same notation for the balls as in the proof of Theorem 1.1, we have, for ǫ > 0 fixed and N largê using the submean inequality in Theorem 2.1 in the last inequality. Letting first N → ∞ and then ǫ → 0 then concludes the proof, using Sanov's theorem again.
It should be stressed that, in general, the functional E(µ) in the previous theorem will not be convex and hence the subset C β ⊂ M 1 (X) consisting of the minima of F β will, in general, consist of more than one element. By general principles the LDP then implies that any limit point Γ ∞ ∈ M 1 (M 1 (X)) of the laws Γ N is concentrated on C β (in the terminology of statistical mechanics Γ ∞ is thus a mixed state defined as a superposition of the pure states δ µ where µ ∈ C β ).
Remark 4.2. The proof of the previous theorem in the case β = ∞ is much simpler as it is does not require the sub-exponential dependence on the dimension in the submean inequality in Theorem 2.1. Indeed, the rough exponential bound used in in the proof of Lemma 3.8 is enough. Moreover, all that is used in the proof for β < ∞ is that ∆ x1 (e −βN H (N ) ) ≥ −λ β e −βN H (N ) for a constant λ β independent on N (but the assumption that ∆ x1 H (N ) ≤ C is a convenient way of ensuring that the previous inequality holds for any β).

Example 4.3.
In the case when X = R n equipped with the Euclidean distance it is known that the mean field Hamiltonian with pair interaction of the form W (x, y) = w(|x − y|) (formula 1.14) Γ−convergences towards E(µ) :=´X 2 W µ ⊗ µ, if w is lower semi-continuous and increasing close to 0 (see [ [18,16,29] for similar results). The proof exploits the explicit nature of E(µ) and a similar argument applies on a compact manifold when W is continuous away from the diagonal with a singularity of the local form w(|x − y|) close to the diagonal (compare [57,58]).
In contrast to the previous example, for the "determinantal" Hamiltonian 5.16 appearing in the complex geometric setting there is no explicit candidate for a limit E(µ). Instead the Gamma convergence is a consequence of the following dual criterion.

4.1.1.
A criterion for Gamma convergence using duality. Next we separate out the convex analysis used in the proof of Theorem1.1 to get the following criterion for Γ−convergence: where f is a Gateaux differentiable convex function on C 0 (X). Then E N converges to E := f * in the sense of Γ−convergence on the space M 1 (X), equipped with the weak topology. (1) and hence taking the sup over all µ ∈ M 1 (X) gives Finally, letting first N → ∞ and then taking the sup over all u ∈ C 0 (X) concludes the proof of the lower bound for E N (µ N ).
To prove the existence of a recovery sequence we first assume that µ = df |uµ for some u µ ∈ C 0 (X). Then, by since, by assumption, f N (u) = f (u) + o(1). Now, by the weak compactness of M 1 (X) the sup defining f N is attained at some µ N ∈ M 1 (X) and hence Next, by a minor generalization of Lemma 3.10 µ N → µ(:= df |uµ ) and hence f * (µ) = 0 + E N (µ N ) + o(1), as desired. Finally, the proof of the existence of recovery sequence for any µ such that E(µ) < ∞ is concluded by a simple diagonal argument based on Lemma 3.6 applied to E := f * .
Thanks to the previous proposition the first assumption in Theorem 1.1 thus implies (also using Lemma 3.8) that E N → E in the sense of Γ−convergence on M(X). Accordingly we recover Theorem 1.1 from Theorem 4.1.
Remark 4.5. In general, if E N gamma converges to a function E on M 1 (X), then it follows (almost directly) that E * N → E * point-wise on C 0 (X). Hence, the point of the previous proposition is that it gives a converse statement under the assumption that E * is Gateaux differentiable. By basic convex duality it thus follows from the previous proposition that E N converges to a strictly convex functional E on M 1 (X) iff E * N → E * point-wise on C 0 (X), with E * Gateaux differentiable.

Relations to the Gärtner-Ellis theorem. First observe that
where Γ N is the measure (δ N ) * (e −βN H (N ) dV ⊗N ) on M 1 (X) and Γ N denotes its Laplace transform on C 0 (X). The Gärtner-Ellis theorem may, applied to the sequence of measures Γ N on M 1 (X), viewed as a subset of the locally convex Hausdorff topological vector spaceM(X), may then be formulated as follows (see (see [31,Cor 4.6.14, p. 148] and references therein): Theorem 4.6. Let H (N ) be a sequence of Hamiltonians on X N and β N a sequence of positive numbers such that β N → β ∈]0, ∞]. Assume that, for any u ∈ C 0 (X), as N → ∞, where F is a Gateaux differentiable function. Then the measures Γ N := (δ N ) * (e −βH (N ) dV ⊗N ) on M 1 (X) satisfy, as N → ∞, an LDP with speed β N N and good rate functional Compared with the Gärtner-Ellis theorem the main point of Theorem 1.1 is thus that, under the quasi-subharmonicity assumption in the second point of the theorem, the assumption that the convergence of the partition functions in formula 4.5 holds for β = ∞ is enough to ensure that one gets an LDP for any β < ∞. Moreover, as a consequence, the convergence then also hold for any β < ∞ with −F β (·) defined as the Legendre-Fenchel transform of the rate functional F β appearing in Theorem 1.1. In fact, the latter convergence is equivalent to the LDP in question, as made precise by the following Lemma 4.7. Let H (N ) be a sequence of Hamiltonians on X N and β N a sequence of positive numbers such that β N → β ∈]0, ∞[. Assume that, for any given volume form dV, the corresponding partition functions Z N.βN satisfy with E(µ) convex. Then the measures (δ N ) * (e −βH (N ) dV ⊗N ) on M 1 (X) satisfy, as N → ∞, an LDP with speed β N N and good rate functional F β . Moreover, if the asymptotics above also holds for β = ∞ with E(µ) strictly convex, then the LDP holds for β = ∞, as well.
Proof. Fixing a volume form dV and applying the asymptotics in the lemma to the volume forms e −βu dV for any u ∈ C 0 (X) reveals that the asymptotics 4.5 hold with f β given by the Legendre-Fenchel transform of E + D dV /β. Now, if E is convex, then E + D dV /β is strictly convex (since D dV is) and hence it follows from basic convex duality that f β is Gateaux differentiable. In fact, the differential µ u := df β|u is the unique minimizer attaining the sup defining f β (u), viewed as the Legendre-Fenchel transform of E + D dV /β. Equivalently, µ u is the unique minimizer of the strictly convex functional µ → E(µ) + u, µ + D dV (µ)/β.
Remark 4.8. Let β N be sequence tending to ∞. By convex duality the Gärtner-Ellis theorem may in the present setting, be formulated as follows (also using Varadhan's lemma [31] in the converse): let E N be a sequence of functions on M 1 (X). Then in the sense of a LDP, with E(µ) strictly convex iff β N N times the log of the Laplace transform of e −βN N EN (δ N ) * (dV ⊗N ) converges to the Gateaux differentiable function E * on C 0 (X).

4.3.
Relations to the existence of the mean energy. Given a sequence of Hamiltonians H (N ) on X N we set If the limit as N → ∞ exists then we will call it the mean energy of µ, denoted bȳ E(µ).
Example 4.9. If H (N ) is the mean field Hamiltonian associated to the pair interaction potential W (formula 1.14) then, trivially, E(µ) =Ē N (µ) for any µ such that W ∈ L 1 (µ).
It follows immediately from the definition that if the limit of E * N (:= f N ) appearing in formula 4.4 exists thenĒ (µ) ≥ f * (µ). (but, in general this is a strict inequality, for example ifĒ(µ) is not convex). In particular, under the assumptions in Theorem 1.1 we haveĒ(µ) ≥ E(µ), where E(µ) appears as the rate functional in Theorem 1.1 for β = ∞ (using Lemma 3.8). Motivated by the complex geometric applications discussed in Section 6 this leads one to consider the following Problem 4.10. Show that the assumptions on H (N ) in Theorem 1.1 imply that the corresponding mean energyĒ(µ) exists when µ is a volume for (perhaps under additional appropriate assumptions on H (N ) ).
As illustrated by the following lemma this problem turns out to be related to the asymptotics of the Gibbs measures with β negative: Lemma 4.11. Assume that there exists some negative β 0 such that for any β ≥ β 0 the corresponding Gibbs measures are well-defined, i.e. Z N.βN < ∞ for N sufficiently large. Moreover, assume that there exists a functional E(µ) such that for any volume form dV. Then the mean energyĒ(µ) exists for any volume form µ andĒ(µ) = E(µ).
Proof. First observe that, by Jensen's inequality, the number f N (β) := − 1 N log Z N.β appearing in the right hand side above for N is concave in β (and, by assumption, finite). Moreover, ∂f N (β)/∂β =Ē N (dV ) at β = 0. Further more, the finite function f (β) defined by the right hand side in formula 4.6 is also concave, as it is an infimum of a family of linear functions and for β = 0 the infimum attained is precisely at µ = dV. Hence, by basic convex analysis, the derivative of f at β = 0 exists and is given by E(dV ). Finally, the proof is concluded by using that if f N is a sequence of convex functions converging point-wise to convex function f such that f N and f are differentiable at 0 then the corresponding derivatives at 0 also converge.
It should, however, be stressed that, if H (N ) is too singular then the partition function Z N.β is equal to ∞, for any β < 0 (even if H (N ) is quasi-superharmonic as in the assumptions of Theorem 1.1). Indeed, for the mean field Hamiltonian corresponding to a pair interaction W this happens as soon as W has a repulsive power-law singularity, i.e. W (x, y) ∼ |x − y| α with α < 0 close to the diagonal. On the other hand, in the case of a logarithmic singularity Z N.βN is indeed finite for β 0 < 0 and sufficiently close to 0 (see [23] for the corresponding LDP in the setting of the 2D vortex model).
Using the Gibbs variational principle some converses to Lemma 4.11 can be established [11], where the existence of the mean energy is assumed (and some additional assumptions), by extending the approach of Messer-Spohn [49]. However is should be stressed that the main point of our proof of Theorem 1.1 is that it does note rely on the existence of the mean energyĒ(µ), which, as pointed above, is an open problem in the present setting.

Applications to Kähler-Einstein geometry
In this section we will apply Theorem 1.1 to complex manifolds X equipped with a line bundle L, assuming that L is positive. The extension to big line bundles (and varieties of positive dimension) is given in the companion paper [9], using the full power of the pluripotential theory developed in [21,14,7] (see the discussion in Section 5.4).

5.1.
Kähler geometry setup. Let X be an n−dimensional compact complex manifold and denote by J the corresponding complex structure viewed as an endomorphism of the real tangent bundle satisfying J 2 = −I.

5.1.1.
Kähler forms/metrics. On a complex manifold (X, J) anti-symmetric two forms ω and symmetric two tensors g on T X ⊗ T X, which are J−invariant, may be identified by setting g := ω(·, J·) Such a real two form ω is said to be Kähler if dω = 0 and the corresponding symmetric tensor g is positive definite (i.e. defines a Riemannian metric) 2 . Conversely, a Riemannian metric g is said to be Kähler if it arises in this way (in Riemannian terms this means that parallel transport with respect to g preserves J). By the local ∂∂− lemma a two form ω is closed, i.e. dω = 0 if and only if ω may be locally expressed as ω = i 2π ∂∂φ, in terms of a local smooth function φ (called a local potential for ω). In real notation this means that (and hence ω is Kähler iff φ is strictly plurisubharmonic). The normalization above ensures that dd c log |z| 2 is a probability measure on C. We will denote by [ω] ∈ H 2 (X, R) the de Rham cohomology represented by ω. If ω 0 is a fixed Kähler form then, according to the global ∂∂− lemma, any other Kähler metric in [ω 0 ] may be globally expressed as where ϕ is determined by ω 0 up to an additive constant. We set The association ϕ → ω ϕ thus allows one to identify H(X, ω)/R with the space of all Kähler forms in [ω 0 ].

5.1.2.
Metrics on line bundles and curvature. Let L be a holomorphic line bundle on X and · a Hermitian metric on L. The normalized curvature two form of · may be (locally) written as in terms of a given local trivialization holomorphic section s of L. The corresponding cohomology class [ω] is independent of the metric · on L and coincides with the first Chern class c 1 (L) in H 2 (X, R) ∩ H 2 (X, Z) (conversely, any such cohomology class is the first Chern class of line bundle L). A line bundle L is said to be positive if it admits a metric with positive curvature, i.e. such that the curvature form ω is Kähler. Fixing a reference metric · on L with curvature form ω 0 any other metric on L may be expressed as · e −u/2 , for u ∈ C ∞ (X) and its curvature is positive iff u ∈ H(X, ω). When L is the canonical line bundle K X , i.e. the top exterior power of the holomorphic cotangent bundle of X : any volume form dV on X induces a smooth metric · dV on K X , by locally setting dz dV := dz/dV, where dz := dz 1 ∧ · · · ∧ dz n in terms of local holomorphic coordinates. When dV is the volume form of a given Kähler metric ω on X, i.e. dV = ω n /n!, then its curvature form may be identified with minus the Ricci curvature of ω, i.e.
where c n dz ∧ dz is a short hand for the local Euclidean volume form i 2 dz 1 ∧ dz 1 ∧ · · · ∧ i 2 dz n ∧ dz n . By a slight abuse of notation we will also write Ric (dV ) for the right hand side in formula 5.1.

Twisted Kähler-Einstein metrics.
A Kähler metric ω β is said to be a twisted Kähler-Einstein metric if it satisfies the twisted Kähler-Einstein equation where the form η is called the twisting form. Since ω β is Kähler the form η is necessarily closed and J−invariant. The corresponding equation at the level of cohomology classes is Fixing, once and for all, a volume form dV on X gives the following one-to-one corresponds between twisting forms η and Kähler forms ω 0 solving the cohomological equation above: The following lemma then follows directly from the expression 5.1 for the Ricci curvature of a Kähler metric: Lemma 5.1. Let X be a compact complex manifold endowed with a J−invariant and closed form η. Then a Kähler form ω β solves the corresponding twisted Kähler-Einstein equation 5.2 iff ω β := ω 0 + dd c ϕ β for a unique ϕ β ∈ H(X, ω) solving the PDE The celebrated Aubin-Yau theorem may now be formulated as follows: Example 5.3. A complex manifold X admits a Kähler-Einstein metric with negative Ricci curvature iff K X is positive (and the metric is unique). Indeed, if K X is positive then, by the very definition of positivity, we can take ω 0 := −RicdV for some volume form on X, ensuring that η = 0 above, with β = 1 (and the converse is trivial).
Remark 5.4. When n ≥ 2 the equation is 5.2 precisely the trace-reversed formulation of Einstein's equations on X (with Euclidean signature): −β is the cosmological constant and η is the trace-reversed stress-energy tensor. Here we are only concerned with the solutions which are Kähler metrics.

5.1.4.
The projection operator P ω0 to the space P SH(X, ω 0 ). Next, we recall the definition of the operator P introduced in [12] (which turns out to be related to the limit as β → ∞ of the equations 5.4). Given u ∈ C 0 (X) we set {ϕ(x) : ϕ ≤ u} which defines an operator P : C 0 (X) → P SH(X, ω 0 ) from C 0 (X) to the space P SH(X, ω 0 ) of all ω 0 −psh functions on X, i.e. all upper semi-continuous functions ϕ in L 1 (X) such that ω ϕ ≥ 0 in the sense of currents. In fact, the operator P preserves C 0 (X) and hence defines a projection operator from C 0 (X) onto P SH(X, ω 0 ) ∩ C 0 (X). More generally, if u ∈ C ∞ (X), the current dd c (P u) has coefficients in L ∞ loc , i.e. which is hence a positive number independent of ϕ ∈ C ∞ (X). Up to a trivial scaling we may and will assume that V = 1. When n = 1 the operator M A may be identified with the Laplacian, but when n ≥ 2 it is fully non-linear. The oneform on C ∞ (X) defined by M A is closed and hence admits a primitive E, i.e. a functional on C ∞ (X) whose differential is given by The functional E is only determined up to an additive constant which may be fixed by the normalization condition E(0) = 0. Using pluripotential theory [21,14] the operator M A can be extended from H ω0 to all of P SH(X, ω 0 ) giving a positive measures satisfyinĝ Similarly the functional E also extends from H(X, ω 0 ) to an increasing lower-semi continuous functional on P SH(X, ω 0 ). We then set (5.9) F (u) := (E • P )(u), which by [12] defines a Gateaux differentiable functional on C 0 (X). More precisely, This setup leads to a direct variational approach for solving complex Monge-Ampère equations, including the Aubin-Yau equation 5.4, in the more general setting of big cohomology classes and singular volume forms dV [14] (compare Section 5.4). However, in the present setting where L is positive the pluripotential theory can be dispensed with by observing that M A(ϕ) is a well-defined probability measure as long as ω ϕ is in L ∞ loc (using that M A(ϕ) is point-wise defined almost everywhere on X ). Then F (u) may be defined by first taking u to be in C ∞ and using the regularity result 5.6 for P u. One then defines F on C 0 (X) as the unique continuous extension of F from C ∞ (X), using that F (u) is Lipschitz continuous on C ∞ (X) with respect to the C 0 −norm (as follows form general principles; see the beginning of the proof of Lemma 3.12).

5.2.
The "temperature deformed" determinantal point processes on X. Let (X, L) be a polarized manifold, i.e. an n−dimensional complex compact manifold X endowed with a positive holomorphic line bundle L. We will denote by H 0 (X, kL) the space of all global holomorphic sections with values in the k th tensor power of L (using additive notation for tensor powers). By the Hilbert-Samuel theorem To the data ( · , dV, β k ) consisting of a Hermitian metric · on L, a volume form dV on X and a sequence of positive number β k we can associate the following sequence of symmetric probability measures on X N k : where det S (k) is a generator of the top exterior power Λ N k H 0 (X, kL), viewed as a one-dimensional subspace of H 0 (X N k , (kL) ⊠N k ) under the usual isomorphism between H 0 (X N k , (kL) ⊠N k ) and the n fold tensor product of H 0 (X, L). The number Z N k ,β is the normalizing constant By homogeneity the probability measure µ (N k ,β) is independent of the choice of generator det S (k) and thus only depends on the data ( · , dV, β k ). We will refer to to the corresponding random point processes on X, as the temperature deformed determinantal point processes on X attached to ( · , dV, β k ) (the special case β k = k defines a bona fide determinantal point process, as recalled below).
Remark 5.5. Since the transformation ( · , dV, β k ) → ( · e −u/2 , e uβ k dV, β k ), for u ∈ C 0 (X), leaves the probability measure 5.11 invariant, the processes above only depend on the data ( · , dV, β k ) through the corresponding two form η, defined by formula 5.3. Moreover, to any twisting form η such that the cohomology class ([η] + c 1 (K X ))/β k defines a positive class in H 2 (X, R) ∩ H 2 (X, Z), i.e. is the first Chern class of a line bundle L, arises from a choice of data ( · , dV, β k ) (compare Section 5.1.3).
It will be convenient to take det S (k) to be the generator determined by a basis s 1 , ..., s N k in H 0 (X, kL) which is orthonormal with respect to the L 2 −product determined by ( · , dV ) for any fixed volume form dV on X : We then take (det S (k) )(x 1 , x 2 , ..., x N ) := taking X to be a non-singular algebraic variety of P m and L as the restriction to X of O(1) gives a polarized where the elements in H 0 (X, kL) are, for k sufficiently large, the restrictions to X of homogeneous polynomials of degree k on P m (in fact, by the Kodaira embedding theorem any polarized manifold (X, L) may, after replacing L with a sufficiently high tensor power, be concretely realized as (X, O(1) |X )). In the case of X = P 1 (=the Riemann sphere) with · denoting the Fubini-Study metric on O(1) whose curvature form ω 0 is the invariant measure on P 1 one can take the base {s i } to consist of monomials and factorize where N = k + 1 and X has been identified with the unit-sphere in Euclidean R 3 and where In the physics literature the corresponding ensemble appears as a Coulomb gas of N unit-charge particles (i.e a one component plasma) confined to the sphere in a neutralizing uniform background ω (see for example [28]). More generally, on any Riemann surface of genus g the bosonization formula [1] gives (5.14) det where G is the Green function of the Laplacian induced by the metric ω 0 and where the second term r appearing above vanishes for genus g = 0, while for g > 0 it may be expressed in terms of the Riemann eta function on the Jacobian torus of the Riemann surface X (giving a contribution which is lower order than the first term; see [58] and references therein). However, when n > 1 it should be stressed that there is no tractable formula for det S (k) (x 1 , ...x N ), even to the leading order.
When β k = k the probability measure µ (N k ,β k ) in formula 5.11defines a determinantal point process i.e. its density can be written as where K (k) (x, y) denotes the kernel of the orthogonal projection onto the space H 0 (X, kL) viewed as a subspace of the space C ∞ (X, kL) of all smooth sections equipped with the L 2 −norm determined by ( · , dV ) [41,5].
The following result generalizes the LDP in [5] for determinantal point processes (or more generally for the case β = ∞) to the general case where β k → β ∈]0, ∞] : Theorem 5.7. Let (X, L) be a polarized manifold and assume given the data ( · , dV, β k ) consisting of a Hermitian metric · on L, a volume form dV on X and a sequence of positive number β k → β ∈]0, ∞]. Then the law of the empirical measures δ N k of the corresponding deformed determinantal point processes with N k particles satisfies a LDP with speed β k N k and rate functional where E ω0 (µ) is the pluricomplex energy of µ with respect to the curvature form ω 0 of · and In particular, δ N k converges in law to the deterministic measure given by the unique minimizer µ β of F β . Moreover, when β < ∞ the measure µ β is the normalized volume form ω β of the twisted Kähler-Einstein metric corresponding to the twisting form η := βω 0 + RicdV.
In fact, the Kähler form ω β may be recovered directly from the limiting volume form µ β by differentiation twice (as follow from the very definition of the twisted Kähler-Einstein equation 5.2): Using basic compactness properties of the space P SH(X, ω 0 ) one then arrives at the following corollary (see [9] for the proof): Corollary 5.8. Given data as in the previous theorem with β ∈]0, ∞[, the following sequence of Kähler forms on X converges to the unique solution ω β of the the twisted Kähler-Einstein metric corresponding to the twisting form η := βω 0 + RicdV.
Remark 5.9. The previous corollary yields a quasi-explicit way of approximating the solution ω β to the twisted KE equation in question (or equivalently the solution ϕ β of the corresponding complex Monge-Ampère equation 5.4), by performing integrals over the spaces X N k −1 of increasing dimension. The procedure becomes explicit as soon as one has constructed bases in the spaces H 0 (X, kL), for k sufficiently large.

5.2.1.
The canonical random point processes on X. We start by recalling the basic fact that, by the very definition of the canonical line bundle K X , any holomorphic section s k of the k th tensor power of K X (i.e. s k ∈ H 0 (X, kK X ) induces a measure on X, symbolically denoted by (s k ∧s k ) 1/k . Concretely, given an open set U ⊂ X with holomorphic coordinates (z 1 , ..., z n ) and writing s k|U = f k dz ⊗k for a holomorphic function f k on U, where dz := dz 1 ∧ · · · ∧ dz n trivializes K X over U, (s k ∧s k ) 1/k |U = |f k | 2/k i n 2 dz ∧ dz, which is independent of U and thus defines a global measure on X (using any holomorphic atlas on X). We also recall that any volume form dV on X induces a metric · dV on the canonical line bundle K X with the property that, if s k ∈ H 0 (X, kK X ) then (s k ∧s k ) 1/k |U may be expressed as dV dV, as follows immediately from the definitions. Now, fixing a volume form dV on X we can apply the relation 5.15 to X N equipped with the induced volume form dV ⊗N and the corresponding metric on L and deduce that the canonical probability measure µ (N k ) on X N k defined by formula 1.6 coincides with the probability measured in formula 5.11 corresponding to the data ( · dV , dV, 1) Hence, Theorem 1.2 is indeed a special case of Theorem 5.7 (also using that η = 0 for this particular data).

5.3.
Proof of Theorem 5.7. To apply Theorem 1.1 in the present setting first note that the Hamiltonian is given by where det S (k) is defined by formula 5.13. The validity of the first assumption in Theorem 1.1 is then a consequence of the following result from [12], where β N k = k : Theorem 5.10. [12]. Let L → X be a positive line bundle equipped with a smooth Hermitian metric · on L with curvature form ω 0 and dV a volume form on X. Then where F is the Gateaux differentiable functional defined by formula 5.9 To verify the second assumption in Theorem 1.1, concerning quasi-superharmonicity, we first observe that we may as well assume that dV is the volume form dV g of the metric g defined by the Kähler form ω 0 . Indeed, dV = e −uβ dV g for some smooth function u and hence changing dV corresponds to changing the metric · to · e −u/2 . Next, we recall that, in general, log s 2 is kω−psh for any holomorphic section s of kL → X (where ω is the curvature form of · ). Hence, we get, for some positive constant λ. Applying the latter inequality to det(s (k) (·, x 2 , ..., x N ) for x 2 , ..., x N thus shows that Theorem 1.1 can be applied to get the LDP in Theorem 5.7. Next, we will show that the unique minimizer µ β of the rate functional F β appearing in Theorem 1.1 coincides with the normalized volume form ω β of the corresponding twisted Kähler-Einstein metric, by applying the general Lemma 3.12. It should however be stressed that while the infimum in the left hand side of formula3.9 is always attained at some µ 0 ∈ M 1 (X) (by weak compactness and lower-semi continuity) this is not so for the right hand side, in general. But in the present setting the sup is attained, when L is assumed to be positive, thanks to the Aubin-Yau theorem. Indeed, first setting g(u) = β −1 logˆe βu dV, for a given β ∈]0, ∞[ gives g * (µ) = β −1 D dV (µ) if µ ∈ M 1 (X) and g * (µ) = ∞ otherwise, as is well-known [31] (and follows from Jensen's inequality applied to the log). Moreover, by the dominated convergence theorem Letting F be the functional on C 0 (X) defined by formula 5.9 the critical point equation 3.10 thus becomes when u is smooth, say. Up to replacing u by u + C we may as well assume that the denominator above is equal to 1. In particular, when u ∈ H(X, ω) the equation when µ ∈ M 1 (X), which coincides with the pluricomplex energy of µ, with respect to ω 0 in [14] (using the notation in [7]). More concretely, a direct calculation reveals that when µ is a volume form where ϕ µ ∈ H(X, ω 0 ) is the solution to the Calabi-Yau equation 1.13, which in Aubin's notation [3] means that E(µ) = c n (I − J)(ϕ µ ) (using [21] the formula above holds for any µ such that E(µ) < ∞, by letting ∧ denote the non-pluripolar products [21]). Thus E(µ) is a generalization of the classical Dirichlet energy on a Riemann surface. The relation F β (ω n ) = κ(ω), where κ denotes the twisted version of Mabuchi's K-energy then follows from the Chen-Tian formula for the K-energy (see [7] and [9] for a direct proof using convex analysis). Moreover, the restriction to H(X, ω 0 ) of the dual functional f (−u) + g(u) appearing in Lemma 3.12 coincides with the Ding functional in Kähler geometry [7]. An alternative proof of the fact that ω n β minimizes F β on M 1 (X) can then be given by using that ω β is a critical point of κ and hence, by convexity, minimizes κ on H(X, ω 0 ). Accordingly, the Calabi-Yau isomorphism ω → ω n shows that ω n β minimizes the restriction of F β to the subspace of all volume forms in M 1 (X). However, showing that the infimum of F β over all of M 1 (X) coincides with the infimum over the subspace of volume forms requires the following non-trivial fact: any µ such that E(µ) < ∞ can be written as a weak limit of volume forms µ j such that E(µ j ) → E(µ) and D dV (µ j ) → D dV (µ) (see [17] where more general results are obtained). 5.4. The generalization to big line bundles and varieties of positive Kodaira dimension. Let us briefly give some indications about the extension of Theorem 5.7 to line bundles L which are merely assumed big, established in the companion paper [9]. In analytic terms L is big iff c 1 (L) contains a positive current on X which is strictly positive in the sense that it is bounded from below by a Kähler form. However, in general, there is a proper open subset Ω ⊂ X such that all positive currents in c 1 (L) are equal to −∞ on the complement X −Ω (which can be taken to be a complex subvariety of X). Fixing a reference smooth Hermitian metric · on L with curvature form ω 0 in c 1 (L) the space of positive currents in c 1 (L) gets identified, as before, with the space P SH(X, ω 0 ) of all ω 0 −psh functions, modulo constants (however, in general all elements in P SH(X, ω 0 ) will be singular along the subvariety X − Ω). Moreover, the non-pluripolar Monge-Ampère operator can be defined on P SH(X, ω 0 ), by restricting to Ω [21]. Then the functional F can be defined essentially as before and Theorem5.7 still holds (again using [12] to verify the first assumption in Theorem 1.1) Invoking, the general Theorem 1.1 thus establishes an LDP with a rate functional F β , admitting a unique minimizer µ β as before. However, one new difficulty is to show that µ β can be written as M A(ϕ β ) for the solution to the equation 5.4 with minimal singularities, whose existence is provided by the general results in [21,14]. The problem is that Lemma 3.12 cannot be applied as it is not clear that ϕ β is of the form P u for some u in C 0 (X) (even if u can be taken to be in L ∞ (X)). But using the variational calculus in [14,7] shows that µ β is of the desired form.
In particular, when K X is big, i.e. X is a variety of general type, the corresponding positive current ω β is the canonical Kähler-Einstein current in X [21,14]. In the general case of a variety of positive Kodaira dimension κ ≤ n (where κ = n iff K X is big) one can use the Ithaka fibration X → Y to represent K X as the pull-back of a big line bundle L on the κ−dimensional manifold Y. Using the Fujino-Mori canonical bundle formula this reduces the proof of the convergence on X to the application of a generalization of Theorem 5.7 concerning big line bundles on Y endowed with a singular volume form dV. As shown in [9] this realizes the corresponding canonical limiting current ω β as the pull-back to X of a (singular) Kähler form on Y solving a twisted Kähler-Einstein equation of the form 5.2, where η is a current on Y determined by the geometry of X (the canonical current ω β first appeared in a different geometric context in [54,55]).
where the constant C k ensures that´X ϕ (k) dV = 0.
The conjectural formula above can be seen as a generalization to the non-linear complex Monge-Ampère operator of the classical Green's formula for the solution of the Poisson equation for the Laplacian on a Riemann surface. Indeed, when X is a Riemann surface the limit ϕ above is precisely given by the Green formula in question (as follows from the bosonization formula 5.14). It turns out that the validity of the conjecture above would follow from the existence of the corresponding mean energyĒ(µ), for any volume form µ (see Problem 4.10). This is shown precisely as in the setting of the real Monge-Ampère operator considered in [8,42] where the analog of the previous conjecture was established using permanents as a replacements of the determinants appearing in the present setting. In particular, when X is a Calabi-Yau manifold, i.e. K X is trivial, the conjecture would imply a quasi-explicit formula for the unique Ricci flat Kähler metric ω ∈ c 1 (L), i.e. solving the Kähler-Einstein equation with vanishing cosmological constant, Λ = 0.
6.2. β < 0. By Lemma 4.11 the existence of the mean energy (and thus the resolution of the conjecture above) would follow if one could establish the asymptotics in formula 4.6 of the corresponding partition functions Z N k ,β/k (assumed finite) for all β > β 0 , for some negative number β 0 . It can be shown that Z N,βN is indeed finite for for some negative β 0 , sufficiently close to zero. In fact, both sides of formula 4.6 are finite when β > β 0 (where the critical negative β 0 depends on (X, L)). This motives the following Conjecture 6.2. Let (X, L) be a polarized manifold and assume given the data ( · , dV ) consisting of a Hermitian metric · on L, a volume form dV on X. For a given negative number β 0 the following is equivalent: • For any β > β 0 the partition functions Z N k ,β are finite for k sufficiently large • For any β > β 0 the functional βF β admits a minimizer on M 1 (X) • For any β > β 0 the measures (δ N ) * e −βH (N k ) dV ⊗N k on M 1 (X) satisfy a LDP with speedN and rate functional βF β (µ) = βE ω0 (µ) + D dV (µ) where E ω0 (µ) is the pluricomplex energy of µ with respect to the curvature form ω 0 of · .
In particular, if the conjectural LDP above holds then the functional βF β is lower semi-continuous and the large N −limit of the laws of δ N k for the corresponding random point processes is concentrated on the (non-empty) set of minimizers of βF β . By [7] any such minimizer is the volume form of a Kähler metric ω β solving the twisted Kähler-Einstein equation 5.2 corresponding to the data (ω 0 , dV, β) and βF β may be identified with the corresponding twisted K-energy functional. Moreover, if the LDP holds then it follows that Z N k ,β ≤ C N β , when β > β 0 . The conjecture should be contrasted with the fact that, in general, βF β is unbounded from below if β is sufficiently negative and even when βF β is bounded from below there exist, in general, twisted Kähler-Einstein metrics whose volume forms do not minimizer of βF β .
In the case when L is the dual −K X of the canonical line bundle, i.e. X is a Fano manifold (which equivalently means that η can be taken to be zero) the equivalence between the two points in the conjecture above can be seen as a probabilistic analog of the Yau-Tian-Donaldson conjecture saying that a Fano manifold X admits a Kähler-Einstein metric with positive Ricci curvature ((i.e. Λ > 0) iff X is Kstable in the algebro-geometric sense; see the companion paper [9] for more detailed explanations of these relations.