Large deviations of radial statistics in the two-dimensional one-component plasma

The two-dimensional one-component plasma is an ubiquitous model for several vortex systems. For special values of the coupling constant $\beta q^2$ (where $q$ is the particles charge and $\beta$ the inverse temperature), the model also corresponds to the eigenvalues distribution of normal matrix models. Several features of the system are discussed in the limit of large number $N$ of particles for generic values of the coupling constant. We show that the statistics of a class of radial observables produces a rich phase diagram, and their asymptotic behaviour in terms of large deviation functions is calculated explicitly, including next-to-leading terms up to order 1/N. We demonstrate a split-off phenomenon associated to atypical fluctuations of the edge density profile. We also show explicitly that a failure of the fluid phase assumption of the plasma can break a genuine $1/N$-expansion of the free energy. Our findings are corroborated by numerical comparisons with exact finite-N formulae valid for $\beta q^2=2$.


Introduction and main results
In recent years there has been a considerable interest in the study of systems with logarithmic interactions. The simplest example of models in this class is the two-dimensional one-component plasma (2D-OCP). This system is also known in literature as 'jellium', 2D Dyson's gas or 2D Coulomb gas [4,16,22,25,31,33,44,50]. The 2D-OCP consists of N identical classical pointlike particles, each carrying a charge q (one species of particle) on a twodimensional domain. The Coulomb interaction between any two particles at distance r from each other is −q 2 v(r), where v(r) obeys the Poisson equation. In the planar case v(r) = log(|r|/L), where L is a length scale that fixes the zero of the potential. To ensure charge neutrality, the particles are embedded in a fixed neutralizing background with opposite charge −qN . The canonical distribution of the 2D-OCP at inverse temperature β is P β,N (r 1 , . . . , r N ) = 1 Z N,β e −βH(r1,...,r N ) , In (2), r i = (x i , y i ) ∈ R 2 denotes the position of the i-th particle of the 2D-OCP (i = 1, . . . , N ) and | · | denotes the Euclidean distance. The first term in H (r 1 , . . . , r N ) is the particle-particle contribution to the energy, while the second term is the particle-background contribution (the 2D-OCP experiences the electrostatic potential V generated by the fixed background). The coupling constant βq 2 is often referred to as plasma parameter.
This statistical mechanics fluid model has appeared in several areas of physics and mathematics. Indeed, the logarithmic repulsion in (2) does occur as interaction between vortices and dislocations in real systems such as superconductors [43], superfluids, rotating Bose-Einstein condensates [17,18,42] (we refer to [2] and [46] for detailed reviews). There is also a well-known analogy between the canonical measure (3) of the 2D-OCP and the Laughlin trial wave function [37] (in the symmetric gauge) in the theory of fractional quantum Hall effect [11]. In the large N limit with fixed potential (without the factor N in front of V (r/L)), the 2D-OCP is equivalent to a class of growth models known as 'Laplacian growth' [50,52].
In the following q and L will be set to one for simplicity. For notational convenience, we also identify R 2 C and denote the positions of the particles in the plane by complex numbers z 1 , . . . , z N . With these conventions, (1)-(2) read P β,N (z 1 , . . . , z N ) = 1 Z N,β e −βH(z1,...,z N ) , For β = 2 and 4, (3) turns out to be the eigenvalues joint distribution for normal complex and normal self-dual matrix ensembles respectively [51,35,8,27]. Remarkably, at inverse temperature β = 2, (3) is a determinantal point process. When the background fills uniformly the unit disk D = {|z| ≤ 1}, the plasma particles experience a quadratic confinement V (z) = |z| 2 /2 and, for β = 2, the system is integrable (the partition function is known exactly for finite N ). In this case, (3) describes the eigenvalues distribution of the celebrated complex Ginibre ensemble [29]. The 2D Coulomb gas analogy has been extensively used in random matrix theory to investigate the large dimensional limit of random normal matrices since the seminal works of Dyson [22].
The present paper deals with some macroscopic features of the 2D-OCP. We present explicit results for the disk configuration of the background. Our quantities of interest are the radial moments ∆ (p) where for notational convenience we have introduced polar coordinates z = re iθ (hence r k = |z k |). The different normalizations for p finite and p = ∞ have been chosen in order to have ∆ (p) for all p. The aim of this paper is to derive the asymptotic distribution of this class of radial observables in the limit of large number N of particles at logarithmic scales 1 (i.e. including the large deviation tails). Our findings can be summarized as follows: (i) We show that the probability density function P (∞) where the left (L) and right (R) large deviation functions are computed explicitly (see (35) and (36) below). Note the different exponential rates of suppression for x ≤ 1 (O(N 2 )) and x > 1 (O(N )). In the derivation of our results, it will be clear that: a) The fluctuations for x ≤ 1 are driven by a collective behaviour of the charges; on the contrary, the fluctuations to the right x > 1 are associated to a spontaneous symmetry breaking in the problem (the equilibrium configuration breaks the rotational invariance in the plane). This change of behaviour of the Coulomb gas is at the heart of the change of speed in the large deviations tails; b) Based on the analysis of the next-to-leading order corrections to (6), which we are also able to obtain (see (51)), we show that a genuine 1/N -expansion of the excess free energy for the 2D-OCP fails to exist if the fluid phase assumption is violated (i.e. when the plasma distribution becomes singular). Indeed, our computations show that, if the equilibrium distribution of the Coulomb gas has a singular component, the large-N expansion of log Z N,β includes a non-universal term of order O(N log N ); c) Studying the limit of ψ (∞) L (x) as x → 1 − , we observe a rather unusual feature: the expected matching between the left rate function and the behaviour of the limiting distribution to the extreme left of x = 1, which is customary in this type of problems, does not actually hold in this case.
(ii) For p finite, we will study in detail the large N behaviour of the moment generating function P (p) We compute explicitly the large N limit up to the next-to-leading order term The scaled logarithm of P (p) N (s) is nothing but the excess free energy of the 2D-OCP. This free energy is given by the electrostatic energy term E p (s) at leading order, while the first 1/N -correction is just a configurational entropy term S p (s). The explicit forms of E p (s) and S p (s) are given in (69) and (70). For p = 1 (mean radial displacement) we recover recent results [19]. Amusingly, for 0 < p < 2 the leading order term E p (s) has a weak non-analytic point at s = 0 due to a change in topology of the plasma distribution. Furthermore, the order of the phase transition (i.e. the order of the first non-continuous derivative of E p (s) at s = 0), while being always at least 3, depends on the parameter p (the larger p the weaker the transition and arbitrarily weak singularities at s = 0 are possible). These non-analyticities of the free energy at the ground state of the plasma unveil a non-regular behaviour of the high order cumulants of ∆ (p) N that cannot be detected by a standard Gaussian approximation. The paper is organized as follows. In section 2 we discuss some aspects of the thermodynamic limit of the 2D-OCP and we introduce the main quantities of interest. Section 3 contains a brief review on the role of the 2D-OCP in random matrix models where some finite-N formulae are available for β = 2. These results will be used to substantiate our findings. In section 4 we present our results on the edge density profile ∆ (∞) N . In particular, in subsections 4.1 and 4.2 we demonstrate the split-off mechanism and we show the breakdown of the free energy 1/N -expansion for a constrained 2D-OCP. The large deviation functions of the radial moments ∆ (p) N (p finite) are considered in section 5. We conclude with some final remarks in section 6.

Thermodynamics of the 2D-OCP
In order to describe the large N limit of the 2D-OCP, it is useful to introduce the normalized density of the gas For large N , when the density approaches a macroscopic limit, the partition function Z N,β of the 2D-OCP can be written as a functional integral over densities where A is the so-called action that can be determined as follows. Assuming a fluid phase for the 2D-OCP, one expects that the partition function admits a large-N expansion [49,50]. Under this assumption, the leading order contributions to A are given by the electrostatic energy and the configurational entropy of the Coulomb gas. The energetic contribution comes from where E[ ] is an energy functional and (z) a short-distance cut-off. The energy functional E is the mean-field approximation of H, In the continuous limit, this functional needs to be corrected by introducing a cut-off that regularizes the particles self-interaction. This correction is provided by the second term in (10), where the short-distance cut-off (z) can be chosen as (z) −1 = N (dµ(z)/dz). A second contribution to the action A arises from the Jacobian J of the change of integration variables from the microscopic variables z k 's to µ: k dz k = [Dµ]N !e log J[µ] , where the factor N ! takes into account the symmetry under permutations of the charges. A standard argument [22,50] gives Hence, by combining all the contributions we identify the action as where E[µ] is the mean-field energy functional (11), is the configurational entropy of the gas density, and C β , C β cannot be determined by these macroscopic considerations. Note that the entropic term at order 1/N in the action vanishes at inverse temperature β = 4. The action representation (9) of the partition function is amenable to a saddlepoint approximation; for large N , the value of log Z N,β is given by the action −βN 2 A[µ eq ] evaluated at the equilibrium density µ eq , i.e. the minimizer of the mean-field energy functional E. The leading term of the action is then E[µ eq ]. Leblé and Serfaty [47,Eq. (1.19)] have recently proved that, under regularity assumptions on the equilibrium measure µ eq , the coefficient of the next-to-leading term is universal and equal to C β = −β/4 (the universality of C β was already known for the 2D-OCP constrained on the line [12,23,10]). In section 4 we will see that this is not the case when µ eq has singular components (non-fluid phase). This explicit example shows that the assumption [47, Eq. (2.7)] of µ eq having a density (used to prove the universality of C β ) is pertinent and not only technical. The next correction for large N includes the configurational entropy S[µ eq ]. Note, however, that if the saddle-point density is not in a fluid phase the configurational entropy (14) diverges, as dµ eq (z)/dz is not finite on the set of singular points. The term C β has been recently investigated in [47]; it has the physical interpretation of microscopic free energy (a weighted sum of specific entropy and 'renormalized energy').
The microscopic free energy C β is universal [47, Eq. (1.17)] (independent of V ); however, there are no explicit results on its value for generic β.
From now on, we focus on the harmonic potential case As already mentioned in the introduction, this quadratic potential is generated by a uniform background in the unit disk D = {|z| ≤ 1}. The relevance of this choice in random matrix theory will be discussed later. It is well-known that, in the limit of large number of particles, the 2D-OCP in the harmonic potential (15) concentrates in a bounded region of the plane. Not surprisingly, the limiting distribution of the plasma is the uniform distribution on the unit disk which is exactly the configuration of the 2D-OCP that neutralizes the uniform background of opposite charge. This distribution on the disk is known as circular law in random matrix theory [7,8,30]. The electrostatic energy and the entropy of the gas are then easy to compute, For β = 2, detailed results on the typical fluctuations of radially symmetric observables of this system have been obtained exploiting the aforementioned determinantal structure. In the next section, we summarize some of these results.

The integrable case β = 2 and random matrices
For β = 2, the canonical distribution of the 2D-OCP (3) assumes a determinantal structure and provides a model whose partition function and correlation functions can be computed explicitly [14,35]. The probability measure in the complex plane (corresponding to V (z) = |z| 2 /2), is well-known in the theory of random matrices. Let G = (G ij ) N i,j=1 be a N × N matrix whose entries are independent complex Gaussian with zero mean and unit variance. Ginibre [29] managed to show that the joint distribution of the N eigenvalues z 1 , . . . , z N of N −1/2 G is given by (18) Kostlan in [36] managed to integrate out the angular variables θ j of the eigenvalues z j = r j e iθj , and showed that, up to a random reshuffling, the eigenvalues moduli r j = |z i | are distributed as a collection of independent χ random variables: (r 1 , . . . , r N ) and σ is a random permutation uniformly distributed in S N . For p = ∞, using Kostlan's theorem (19), the following identity holds where γ(k, y) = y 0 dt t k−1 e −t denotes the lower incomplete gamma function. Differentiating (20), one gets a finite-N formula for the probability density of ∆ Using a quantitative version of the central limit theorem, Rider [40] has proved the following remarkable limiting result for the typical fluctuations of ∆ A similar result has been recently established [13] for a general class of radially symmetric external potentials V (z) = V (|z|). Rider's theorem (22) describes the typical fluctuations of the 'top eigenvalue' of complex Ginibre matrices. In this work, we address the question of atypical fluctuations. Few explicit results are available on large deviations for the 2D-OCP, and they are usually valid at the leading order in N , see [5,19,32]. Forrester [26] obtained a finite-N formula for the moment generating function of radially symmetric linear statistics. Specializing his result to ∆ 2 In other words ξ 2 k d = χ 2 2k /2. and from here it is possible to extract the average and variance of ∆ For a rigorous proof of central limit theorems at β = 2 for radially symmetric observables (including ∆ (p) N with average and variance as above) see [41]. Atypical fluctuations are not described, of course, within the Gaussian approximation. Formulae (21), (23) and (24) will be used in the next sections to corroborate numerically our results, as well as to prove some of our claims for β = 2.

Large deviations of the edge density profile
The large N distribution of the plasma µ disk (16) implies that the edge density profile ∆ The above limit can be written as where Z N,β (x) is the partition function of the 2D-OCP with the constraint r k ≤ x for k = 1, . . . , N . As discussed in section 2, (−1/βN 2 ) log Z N,β (x) is dominated for large N by the saddle-point of the action in (9). The saddlepoint is given by the distribution of the plasma that minimizes the electrostatic energy functional with the constraint that all charges are at distance at most x from the origin. Therefore we have where µ x minimizes E[µ] in (27) with the constraints Finding the minimizer of a logarithmic potential under constraints is usually a technically involved task. From electrostatic considerations, we know that the equilibrium distribution of the plasma concentrates on a compact set. In the interior of the support, the distribution has constant density 1/π and the excess of charges accumulates on the boundary of the support 3 . Note that the support of the distribution is itself unknown. In our case, we can take advantage of the rotational symmetry of the problem and eventually we find For x ≥ 1, the saddle-point density is independent of x (the constraint on the 2D-OCP is ineffective). When x < 1, the constraint forces the 2D-OCP in a smaller disk and a fraction of charge accumulates on the boundary of the disk. More precisely, in the interior of the disk of radius x the density is constant and equal to 1/π; the exceeding charge (1 − x 2 ) lies exactly on the boundary of the disk and, by symmetry, is uniformly distributed. Hence, for x < 1 the equilibrium distribution has a singular component (the delta measure at distance x = r in (30)) and we say that the plasma is not in a fluid phase. Evaluating E[µ x ] and using (17) where θ(x) is the Heaviside theta function. Note that (31) is identically zero for x ≥ 1, while for x → 1 − it vanishes as (x − 1) 3 . Therefore the third-order derivative of the large deviation function (31) is discontinuous at x = 1. Similar non-analiticities in the large deviation function of the edge density profile for 2D Coulomb gases in one dimension have been recorded in previous works on Hermitian random matrices [39]. Eq. (31) where P (∞) is the probability that the edge density takes value x. The idea is to compute this limit as energy cost of a configuration with one charge of the Coulomb gas at distance x > 1. The intuition is that, for large N , the most probable configuration of the gas with at least one particle at distance x > 1 is described by the (unperturbed) equilibrium measure µ disk (i.e., the circular law) with one single particle at distance x. Among all possible realizations of the constraint ∆ (∞) N = x (with x > 1), this is the least energetic one 4 . The presence of one particle outside the unit disk corresponds to a spontaneous symmetry breaking in the system: the problem is symmetric under rotations U (1) in the complex plane, however the plasma configuration with the split-off of a single particle is necessarily no longer symmetric. The computation of (32) therefore amounts to estimating the energy cost in pulling one particle away from its equilibrium position inside the unit disk D and relocating it around a generic position xe iθ outside the disk (x > 1). If the plasma distribution is µ disk then the electrostatic energy of a particle at xe iθ is Note that inside the unit disk (x ≤ 1) the electrostatic energy is constant (the gas is at equilibrium). The energy cost to transport a particle outside the unit disk at distance x > 1 is given by the energy difference, and hence We can summarize the statistical fluctuation of the edge density ∆ (∞) N as in (6) with left (L) and right (R) large deviation functions L (x)). Hence, the cumulative distribution and the probability density have the same large-N decay. We notice that the large deviation result (31) does not exhibit a matching with the limiting result by Rider (22) for typical fluctuations. This is due to the fact that the limiting (Gumbel) distribution (22) decays too fast (super-exponentially) towards −∞ to meet the left rate function at an exponential scale. In Fig. 1 we show a comparison of these large deviation functions with the finite-N formula (21) valid for β = 2. In the following subsections, we will corroborate our findings (35)- (36) with an asymptotic expansion of the finite-N formulae (20)-(21) valid for β = 2. Moreover, the following analysis will i) demonstrate the split-off mechanism associated to atypically large values of ∆ (∞) N , and ii) show that the excess free energy of a 2D-OCP constrained on a disk of radius smaller than 1 does not have a genuine 1/N -expansion: the presence of spurious logarithmic terms in the large-N expansion is due to the constrained 2D-OCP being in a non-fluid phase.

Proof of the split-off mechanism for β = 2
In the previous section we have computed the limit (32) using large deviations ideas. A key step in the computation was finding the typical plasma distribution realizing the atypical event ∆ (∞) N = x with x > 1. Guided by intuition, we have guessed that the typical way for the 2D-OCP to satisfy the constraint on the edge density profile is by a split-off mechanism. Evaluating the limit as energy cost for the splitting we found (34). Here we obtain the same result (i.e. Ψ (∞) R (x)) using an asymptotic expansion of a finite-N formula, thus proving the split-off phenomenon.
Let us consider the finite-N formula (21) valid at β = 2: where the Euler-Maclaurin summation formula [1, Eq. 23.1.30] has been used. It is easy to verify that log γ(N t, N for t ≤ 1 < x, and hence we get Using a saddle-point approximation (the exponent achieves its minimum at the boundary t = 1 of the integration domain) we eventually obtain the same result as in (34) lim This proves, at least for β = 2, that the typical configuration of the 2D-OCP with edge density ∆ (∞) N > 1 is given by the unperturbed uniform distribution on the unit disk with one charge located at distance x > 1 from the origin (hence one charge splits from µ disk ).

Large-N expansion for the constrained 2D-OCP
We now consider the asymptotics of the free energy of the 2D-OCP confined in a disk of radius x ≤ 1. This amounts to considering the asymptotics of We have seen that the plasma under the constraint r k ≤ x ≤ 1 is not in a fluid phase (see (30)). At equilibrium, the charges of the 2D-OCP accumulate on the boundary of the disk and the plasma distribution is not absolutely continuous. We want to show here that the excess free energy of the 2D-OCP with hard constraints does not admit a genuine 1/N -expansion. Using the Euler-Maclaurin summation formula we can write up to o(N ) error We have extracted the large-N asymptotics of (42). Up to an error of order o(N ), we have the following asymptotics for large N and x ≤ 1: where J N (x, ξ) in (43) is defined as The expression (47) lends itself to an evaluation in a saddle-point approximation. The largest contribution to the integral comes from the neighborhood of the point that minimizes the exponent g ξ (t) = t − ξ log t. This function achieves its global minimum at t = ξ, which may or may not be in the interior of the integration range [0, x 2 ]. Let us first consider the expansion in the case of 1/N < ξ < x 2 (saddlepoint in the interior of the integration range). Then, we can expand the exponential term around t = ξ and obtain For x 2 < ξ ≤ 1, the function g ξ (t) achieves its minimum at the boundary t = x 2 of the integration range and hence Elementary integrations provide now the asymptotics of (43) using the large-N expansions (48)- (49): Collecting all terms in (50) and (44)- (46) we obtain the following large-N expansion for x ≤ 1 (at β = 2): where Ψ (∞) L (x) is the left rate function (35) and the correction terms are Note the nonzero coefficient of O(log N/N ). This explicit calculation shows that when the 2D Coulomb gas is not in a fluid phase (here for x < 1), the coefficient C β in the large-N asymptotics (13) is not universal and therefore the excess free energy log(Z N,β (x)/Z N,β ) does not admit a genuine 1/Nexpansion. Note that here the hard constraint r k ≤ x has more dramatic consequences on the free energy expansion of the 2D-OCP compared to the one-dimensional case (Hermitian random matrices) [9,15]. For Hermitian matrix models, the large-N eigenvalue density with a constraint on the largest eigenvalue is absolutely continuous [9] and hence the term O(log N/N ) in the large deviation asymptotics is absent. The expansion (51) has been compared with the exact β = 2 formula. See Fig. 2.

Large deviations of radial moments
For finite p, it is more convenient to work in Laplace space and consider the moment generating function P Taking the logarithm, we obtain the generating function of cumulants as excess free energy of the perturbed gas with respect to the unconstrained system where Z (p) N,β (s) is the partition function of the 2D-OCP in a modified external potential Using the large-N representation (9) of the partition function in terms of an action, one evaluates (54) as where A (p) s is the action of the 2D-OCP in the external potential (55) and µ (p) s is the saddle-point density of A (p) s . One readily sees that the saddlepoint density µ (p) s is the equilibrium distribution of the plasma, i.e. µ The asymptotic behaviour of the cumulant generating function is obtained up to the next-to-leading order: where The leading term (59) is the rescaled cumulant generating function of ∆  s (|z|) is sufficiently regular [47]) and hence we have no logarithmic terms in the excess free energy (54)-(58).
To summarize, the statistics of ∆ (p) N is governed at leading order by the electrostatic energy excess E p (s) of the equilibrium distribution in the potential V (p) s (independent of β) with the first correction in 1/N expressed in terms of a configurational entropy term S p (s). Again, the technical challenge is the solution of a constrained optimization problem. In the next section, we solve the variational problem and we show that the statistical behaviour of ∆ (p) N is quite rich and depends on the exponent p > 0.

Solution of the variational problem
Before tackling the optimization problem, notice that the moment generating function P These different ranges correspond to the stability values of s for the 2D-OCP in the effective potential (55). Observe that the equilibrium distribution µ (p) s will inherit the rotational symmetry of the energy functional E (p) s . In fact, the minimization problem can be solved explicitly using the arguments presented in [19]. A suitable modification of those steps provides the solution of the variational problem for s ∈ Ω p as with (p) s (r) given by The equilibrium distribution µ (p) s of the plasma is supported on a radially symmetric domain where r p (s) is the largest positive solution of V s (r p ) = 0 and R p (s) is the unique positive solution of V s (R p )R p = 1. These equations read explicitly Note that as long as s ≥ 0 the effective potential (55) is convex and therefore the plasma distribution is supported on a disk (r p (s) = 0). However, for p < 2, when s ∈ Ω p is negative the convexity of the potential is broken and the inner radius r p (s) becomes positive. Hence, for s < 0 the support Σ p (s) of the plasma distribution is an annulus. This change of topology of the equilibrium distribution has important consequences on the statistical behaviour of the radial moments. The equilibrium measure µ The electrostatic energy excess of µ and eventually one finds for s ∈ Ω p Later we provide more explicit expressions for E p (s) and S p (s) for some values of p. We notice that, for 0 < p < 2: where c p is a constant. Hence, E p (s) is not analytic at s = 0 for 0 < p < 2. Since E p (s) is the leading order in N of the excess free energy (56), this non-analyticities corresponds to phase transitions. Note that the phase transitions occurs at the ground state of the 2D-OCP (at s = 0) and correspond exactly to the disk-to-annulus change in topology of the equilibrium plasma distribution µ (p) s in (62). More precisely the -th derivative of E p (s) is discontinuous when 2 1 − 2 −1 < p ≤ 2 1 − 2 with ≥ 3. For instance, for 0 < p ≤ 2/3 we have a third-order phase transition, for 2/3 < p ≤ 1 a fourth-order phase transition, for 1 < p ≤ 6/5 the transition is of order five, and so on. For p = 2 the function E p (s) is analytic. We stress that, although the underlying mechanism is the same (the change in topology), the order of the singularity in the free energies of the 2D-OCP depends on the particular statistics 5 , i.e. on the exponent p.
Having in mind these weak singularities, one can obtain the first cumulants κ (∆ (p) N at leading order in N as derivatives of E p at s = 0 (for < 4 2−p ). For instance, one finds  (24). The first correction to these results are obtained as derivatives at s = 0 of S p (s).

Explicit formulae for p = 1 and p = 2
For concreteness, we write here the form of E p and S p in terms of elementary functions for some special values of p. First, the expression of E p (s) can be simplified further as For p = 1 (the mean radial displacement) we obtain   Eq. (74) was obtained in [19]. Note that the fourth-order derivative of E 1 (s) is discontinuous at s = 0, in agreement with the previous discussion.
For the moment of inertia ∆ In this case E 2 (s) is real analytic, according to the fact that for p = 2 there is no disk-to-annulus transition in the plasma distribution. In Figs

Concluding remarks
In summary, we have obtained a number of results on radial observables ∆ (p) N (as defined in Eq. (5)) for the two-dimensional one-component plasma in the limit of large number N of particles and for generic values of the coupling constant. For finite p, the moment generating function of the p-th radial moment is obtained with logarithmic accuracy up to the first sub-leading term in the 1/N -expansion. The leading order E p (s) of the excess free energy of the plasma is shown to exhibit a weak non-analytic point at s = 0 for 0 < p < 2: interestingly, the order of this non-analyticity (i.e. the order of the first noncontinuous derivative), while being always larger or equal than 3, is found to depend on p and arbitrarily weak singularities at s = 0 are possible. In the mean-field representation of the problem, these singularities are associated to disk-to-annulus transitions in the equilibrium distribution of the plasma. Since E p (s) is the generating function of cumulants of ∆ (p) N at leading order in N we expect that the non-continuity of the derivatives of E p (s) at s = 0 may be related to a 'non-regular' behaviour of high cumulants as N → ∞. For p = 1 this was numerically observed in [19].
We have considered in detail the edge density profile ∆ (∞) N of the plasma in the plane, i.e. the fluctuations of the farthest particle from the origin. We obtained the rate functions governing the probability of atypical fluctuations to the left and right of the expected position x = 1. The difference in speeds (O(N 2 ) vs. O(N )) between the two large deviation principles is due to different arrangements of the fluid particles realizing an atypical configuration: large fluctuations to the left are realized cooperatively, while to the right by a splitting-off mechanism, which breaks the rotational symmetry of the fluid. This physically intuitive explanation has been also demonstrated for β = 2 by the asymptotic expansion of the finite-N formula (21). To the best of our knowledge, this is one of the few existing proofs of the split-off phenomenon associated to the fluctuations of the maximum of a set of correlated random variables ({r k } in our problem).
It is worth noticing that, at odds with what generally happens in this kind of problems, the large deviation result (31) does not exhibit a smooth matching with the limiting result by Rider (Eq. (22)) for typical fluctuations to the left of the expected value, as the Gumbel distribution decays too fast (super-exponentially) towards −∞ to meet the left rate function at an exponential scale. This problem gives then rise to another of the very few instances where the conventional matching between typical and atypical fluctuations is more delicate (see e.g. [45,48]).
A careful asymptotic analysis, extended to the first two sub-leading orders, also reveals that a genuine 1/N -expansion of the excess free energy requires a fluid-phase assumption; in our problem the equilibrium density µ x in (30) is not absolutely continuous for x < 1, and the large-N expansion of log Pr[∆ (∞) N ≤ x] for β = 2 contains a nontrivial O(N log N ) term (see (51)). This fluid-phase condition should be compared with the so-called one-cut assumption for the 2D-OCP on the line [10]. Note that, strictly speaking, the 1/N -expansion also breaks down in Hermitian matrix models with a hard constraint of type λ max ≤ x (known as 'hard wall'), λ max being the largest eigenvalue. More precisely, one branch in the asymptotics of Pr[λ max ≤ x] contains a nonzero term O(log N ) [9,Eq. (4.35)]. This term is only a 'normalisation constant' independent of x, due to the presence of the hard constraint. Similar spurious O(log N ) contributions of geometric nature are also expected in two-dimensions.
It would be interesting to see if the results reported in this paper can be extended, at least qualitatively, to the 2D-OCP in more generic confining potentials [14,35,53] (other than radially symmetric [3,30,28]) or for the plasma on different planar surfaces (e.g. on a cylinder [11]).