{'fraction_non_alphanumeric': 0.05125550090603158, 'fraction_numerical': 0.021330572094227285, 'mean_word_length': 3.8075997676926905, 'pattern_counts': {'":': 0, '<': 25, '<?xml version=': 0, '>': 59, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 40, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Local incompressibility estimates for the Laughlin phase

We prove sharp density upper bounds on optimal length-scales for the ground states of classical 2D Coulomb systems and generalizations thereof. Our method is new, based on an auxiliary Thomas–Fermi-like variational model. Moreover, we deduce density upper bounds for the related low-temperature Gibbs states. Our motivation comes from fractional quantum Hall physics, more precisely, the perturbation of the Laughlin state by external potentials or impurities. These give rise to a class of many-body wave-functions that have the form of a product of the Laughlin state and an analytic function of many variables. This class is related via Laughlin’s plasma analogy to Gibbs states of the generalized classical Coulomb systems we consider. Our main result shows that the perturbation of the Laughlin state cannot increase the particle density anywhere, with implications for the response of FQHE systems to external perturbations.


Introduction
The fractional quantum Hall effect (FQHE) [25,29,33,69] is a remarkable feature of the transport properties of 2D electron gases under strong perpendicular magnetic fields and at low temperatures. Soon after its experimental discovery [71], it was recognized in the seminal works of Laughlin [31,32] that the origin of the effect lies in the emergence of a new, strongly correlated, phase of matter. The latter has been argued to host elementary excitations with fractional charge, a fact that was later experimentally confirmed [21,47,63]. Even more fascinating, but still lacking an experimental confirmation, is the possibility that these excitations (quasi-particles) are anyons [6,45,73], i.e. have quantum statistics different from those of bosons and fermions, the only known types of fundamental particles. Due to these exciting prospects, it is an ongoing quest in condensed matter physics to generalize FQHE physics to other, more flexible, contexts than the 2D electron gas [9,11,20,49,72].
One of the distinctive features of strongly correlated FQH states is that they are incompressible liquids, rendering them very robust against perturbations and external fields. Incompressibility, in the form relevant to the FQHE, has in fact two aspects: 1. The proposed strongly correlated FQHE wave functions are approximate ground states for the many-body Hamiltonian of the system at hand. Their energy is separated from the rest of the spectrum by a gap independent of volume and particle number.

2.
Modifications of model FQHE ground states that stay within the (highly degenerate) ground eigenspace of the many-body Hamiltonian cannot increase the local oneparticle density beyond a fixed value.
These two aspects are very non-trivial to check, particularly since we are talking about strongly correlated states of matter for which mean-field descriptions in terms of independent particles are not adequate. The main evidence so far in favor of Properties 1 and 2 has been experimental and numerical and it is an important theoretical challenge to improve on this.
The main paradigms of the FQHE are encoded in Laughlin's wave-function, the simplest FQH state: Here z 1 , . . . , z N ∈ C are the positions of N particles moving in R 2 , identified with the complex plane, and the constant c Lau is a normalization factor (setting the L 2 -norm equal to 1). For fermions, is odd and ≥ 3 (the case = 1 corresponds to noninteracting fermions), while for bosons ≥ 2 is even. The function (1.1) has originally been proposed as a variational ansatz for the ground state of the many-body magnetic Schrödinger Hamiltonian 1 acting on L 2 (R 2N ), the Hilbert space for N 2D particles. Here x ⊥ denotes the vector x ∈ R 2 rotated by π/2 counter-clockwise, so that and thus B 2 x ⊥ is the vector potential of a uniform magnetic field, expressed in symmetric gauge.
To arrive at an ansatz of the form (1.1) one assumes that the energy scales are set, in order of importance, by the strength B of a constant applied magnetic field perpendicular to the plane, the repulsive pair-interaction potential w, the external scalar potential V , representing trapping and/or disorder.
In the sequel we choose units so that the strength of the magnetic field is 2 and the magnetic length therefore 1/ √ 2. Since the magnetic field is the main player, the first reduction is to replace the state space L 2 (R 2N ) N L 2 (R 2 ) by N H where (z is the complex coordinate in the plane) is the lowest Landau level (LLL), ground eigenspace of the magnetic Laplacian −i∇ − B 2 x ⊥ 2 . In other words, one reduces available one-body orbitals to those minimizing the magnetic kinetic energy. Spins are all aligned with the external magnetic field so that one only needs consider symmetric/antisymmetric wave-functions of N H to describe bosons/fermions. Secondly, in order to suppress repulsive interactions, the many-body wave-function of the system should vanish when two particles meet. Combined with the regularity of the LLL orbitals, this forces the inclusion of Jastrow factors (z i − z j ) , ∈ N in the wave-function, leading to (1.1). The integer is chosen to accommodate symmetry (the wave-function must be symmetric/antisymmetric under coordinates exchanges z i ↔ z j to describe bosons/fermions) and to fix the filling factor (number of particles per magnetic flux quantum) equal to 1/ . In the absence of an external potential, there is then no free variational parameter in the ansatz (1.1).
The proposed form (1.1) turns out to be sufficiently robust for serving as the basic ingredient of approximate ground states for large classes of repulsive interaction potentials w and external potentials V .
A possible mathematical formulation of Property 1 is to consider a model pairpotential w whose projection to the N -body LLL has (1.1) as an exact zero-energy ground-state [26,48,70]. The claim would then be that the next energy level is bounded below, independently of N . This problem has remained open so far, see [39,58] for further discussion.
Property 2 is the main subject of the present contribution. To formulate it properly, we note that the arguments leading to choosing (1.1) as a variational ansatz leave the possibility to choose any L 2 -normalized function of the form with F analytic and symmetric under exchange of the z i . This form exhausts the class of functions that minimize the magnetic kinetic energy and at the same time avoid repulsive interactions by vanishing at least as (z i − z j ) as z i and z j come together. In the bosonic case and with = 2 these are exactly the ground states of the contact interaction [58, Section 2.1]. The factor F allows an adaption to an external potential. We shall refer to the class of states of the form (1.4) as fully correlated states. Consider now the one-particle density of a state of the form (1.4). For F = 1, i.e., for the Laughlin state (1.1), it was argued in [31] and proved rigorously in [57,58] that the density takes essentially the constant value 1/(π ) in the disk D(0, √ π N ), and drops quickly to 0 outside of the disk. 2 Property 2 can then roughly be formulated as Some care has to be taken in formulating this rigorously for, in view of existing numerical simulations [18,19], the symbol above cannot actually stand for a microscopic pointwise bound for finite N .
In this paper we improve and generalize results of [59,60] (see also [55, Chapter 3] for further discussion) by proving that (1.6) holds in the sense of local averages over length scales N 1/4 , which are much smaller than the extension of the Laughlin state itself, of order N 1/2 . In view of recent results [7,34], which can be applied to the pure Laughlin state, i.e., F = 1, it is natural to conjecture that (1.6) in fact holds on any length scale N α , α > 0. In order to understand the response of the Laughlin state to external potentials, it is essential, however, to know that the density bound holds for all analytic factors F. The methods of [7,34], in contrast to those we develop here, do not apply in such generality.
The bound (1.6) is an expression of the robustness of the Laughlin state against perturbations by external fields because it forbids a fully-correlated state to accommodate variations of an external potential by concentrating arbitrarily large portions of its mass in energetically favorables places: Any redistribution of the mass must respect the same density bound (1.6) as the Laughlin state.
On the other hand, it has recently been proved in [61] that, if (1.6) is known to hold for all strongly correlated states, then the minimal energy in an external potential within this class of states can asymptotically be attained in states saturating the bound (1.6) and having the simple form where f is a polynomial of a single variable. Physically, this means that it is never favorable to add more correlations to the Laughlin state in order to accommodate an external potential, provided the latter is weak enough so that it does not make the ground state jump across the energy gap assumed in Property 1.
The zeros a j ∈ C of the polynomial f , (1.8) are interpreted as the locations of quasi-holes, each carrying a charge 1/ . If the locations a 1 , . . . , a J are treated quantum-mechanically instead of classically as above, one should expect that the quasi-holes behave as anyons [6,45,73] with statistics parameter −1/ .
The results of the present paper, together with [61], confirm rigorously Laughlin's original intuition [31,32] that the response of the Laughlin state to a disorder potential is to generate uncorrelated quasi-holes, which is also in accord with experimental studies [47].
In the next section we state our main results precisely and discuss them further (they were announced in the short paper [42]). The rest of the paper is devoted to their proofs.

Universal density bounds.
As discussed in the introduction (see also [59,60]), we assume that the magnetic field and repulsive interactions are strong enough, so that the ground state of the system has the form (1.4). We aim at proving that there is a universal upper bound to the local density of all such states, on scales much smaller than the total size of the system. Thus, let be the L 2 -normalized Laughlin function with exponent ∈ N, in units where the magnetic length is 1/ √ 2. As previously explained we consider analytic perturbations thereof: For a function F ∈ L N we denote by 3) the associated one-particle density. Our main result is a universal upper bound on ρ F , holding on mesoscopic length scales: as N → ∞.

Potential energy estimates.
With Theorem 2.1 in hand, we can return to the original physical problem (1.2) in the presence of an external potential V . If the magnetic kinetic energy is frozen by restricting one-body orbitals to the LLL, and if we further assume that the interaction energy is also frozen by restricting to (2.2), the only non-trivial energy term left in the Hamiltonian is set by V . An interesting problem consists in studying the minimization of this energy amongst fully-correlated states: This infimum will not necessarily be attained, but it is easy to see that, for fixed N , the density ρ F is bounded in L ∞ uniformly in F (this follows from the fact that all onebody orbitals are in the LLL [14], see the discussion in [59,Remark 2.2]). Therefore the infimum exists for any V such that its negative part, In view of Theorem 2.1, it is fairly natural to expect a lower bound to (2.6) in terms of the "bathtub energy", defined as [41, Theorem 1.14] We can prove this for a large class of external potentials V , that we now describe. Since the typical extension of the Laughlin state is of order √ N , it is more convenient to state our assumptions in terms of a scaled version of the potential. Assumption 2.2 (The external potential). Let V ∈ C 2 (R 2 ) be a (sequence of) potential(s) that we write in the manner for a (possibly N -dependent) U ∈ C 2 (R 2 ). Assume that there exists a fixed radius R > 0 and an α > 0 such that (1) Local behavior: and independently of N .
A typical example we have in mind is a potential V of the sort where the part U trap represents trapping in a given sample, and satisfies Item 2 of the assumptions, whereas U disorder represents some disorder due to impurities in the sample, is compactly supported in a fixed disk and satisfies Item 1 of the assumption. Item 1 means that we consider a potential V whose variations happen on length scales O(N α ), α > 0, the best scale we may hope for (but cannot quite reach) in our density bounds, Theorem 2.1. Note that this should not be interpreted as saying that the potential energy bounds we are about to state work on the optimal length scale, see Remark 2.4 below. Item 2 is essentially a mild trapping assumption. It is clearly satisfied by any regular enough function U growing at least like c log |x|, c > 2 at infinity.
We shall prove the following corollary, which improves [60, (2.14) The above is obtained by using trial states of the form (1.7). For the special case where V is radially symmetric and monotonously increasing in the radial variable, or is of a Mexican-hat form, such an upper bound was proved earlier in [59,60]. A remarkable consequence is that the Laughlin state stays an approximate ground state in any radial increasing potential, however steep. We refer to [59,60] again and to the companion paper [42] for further discussion of the significance of this.

2.3.
Proof outline: screening properties in the plasma analogy. The backbone of the proof of Theorem 2.1 is Laughlin's plasma analogy, a very fruitful mapping of the many-body density of Lau , and more generally F , to a problem in classical statistical mechanics, originating in [31,32]. As in our previous papers [57][58][59][60] we generalize this to any fully-correlated state (1.4). We rewrite the N -particle probability density of any fully correlated state N in the form of a Boltzmann-Gibbs factor with a classical Hamiltonian To interpret (2.16) as an energy functional for N classical positively charged particles with coordinates z 1 , . . . , z N ∈ C R 2 we recall that Poisson's equation relates the electrostatic potential ϕ generated by a charge distribution ρ in the manner where is the two-dimensional Laplacian. Thus one can see (2.16) as the energy of N positively charged particles subject to the external potential generated by a constant background of opposite charge. Since the potential |z j | 2 can be interpreted as being generated by a uniform charge background of negative density −2/π . pairwise Coulomb repulsion: thus the second term in (2.16) is the 2D Coulomb interaction between particles of identical charge √ 2 . the potential generated by an essentially arbitrary, but positive, (many particle) charge distribution encoded in the (modulus of the) analytic function F. Indeed, for any This is, in fact, the only property of F we use in our method. Equation (2.20) means that − log |F| is superharmonic in every variable, because F is analytic. 4 A plausibility argument for our main results goes as follows: (1) The effective temperature in (2.15) is small compared to the total energy, so that the Boltzmann-Gibbs factor is essentially concentrated around minimizing configurations of the Hamiltonian (2.16).
(2) For F = 1, i.e. in the case of the pure Laughlin state, one should expect that the local density of points in a minimizing configuration is close to 1/( π ). Indeed, this is the condition for the negatively charged fixed constant background (c.f. (2.18)) to efficiently screen the positively charged mobile particles z j , j = 1 . . . N (c.f. (2.19)). (3) If F = 1 we add an additional positive charge distribution to the game because of (2.20). Since the mobile particles z j , j = 1 . . . N are positively charged in this representation, one can expect that such an addition will exert a repelling force on the density distribution of the mobile particles and reduce it, at least on the average. This intuitive picture is, however, not as clear as it may look at first sight. If we imagine the potential term associated to F as being generated by point-like positively charged particles, it is not obvious a priori why several such charges cannot work together and conspire to increase the density locally. Arranged tightly on a circle they could, perhaps, concentrate the density at a high value around the center of the circle. Note also that −2 log |F(z 1 , . . . , z N )| is generated by a mobile charge distribution because z j log |F(z 1 , . . . , z N )| is in general a function of z 1 , . . . , z N and only constant in special cases (essentially when F = f ⊗N is a tensor power, as in (1.7)). Thus, the positions of the additional charges are in general correlated with those of the actual charges. Moreover, as the example F(z 1 , . . . , z N ) = i exp(cz i ) shows, the area occupied by the bulk of the density can be stretched in one direction but compressed in another.
The heuristic argument (1)-(3) above is, in any case, certainly far from a rigorous proof of the bound (2.4). To obtain such a proof, the reduction to considering only ground states of (2.16) (rather than the full Gibbs state) presents, perhaps, the least difficulty. To appreciate this, it is helpful to rescale lengths: we can write (2.21) as a probability measure with Thus, modulo a simple change of length-scales, we can reduce the aim to a statistical mechanics problem in mean-field scaling, the first and second sum in the above being of the same order for large N .
This makes it clear why mostly ground-state configurations are relevant for the proof, although some care is needed to estimate the effect of entropy terms in the classical free energy. In fact, the reduction to classical ground state configurations is responsible for the fact that our main result holds only on length scales N 1/4 . For ground-states instead of Gibbs states, the result holds on scales 1. There is thus certainly room for improvement in this particular part of the proof, see Remark 2.6 below.
With the reduction to ground states of (2.16) out of the way, the main question is now to understand points (2) and (3) of the list above. In our approach, the solution is essentially based on screening properties of ground states of classical Coulomb systems, namely their tendency to achieve local neutrality in order for the total potential generated by the positively and negatively charged particles to compensate one another.
Consider a minimizing configuration Z 0 N = z 0 1 , . . . , z 0 N for the classical Hamiltonian (2.16) (the labeling of the points is irrelevant by symmetry of (2.16) under permutations). To any subset z 0 1 , . . . , z 0 K of the configuration (in fact, any K -tuple of points in R 2 ) we associate an open screening region K such that The above means that the total potential generated by the K positive charges under consideration and a constant background of opposite charge density (π ) −1 contained in K vanishes identically outside K . The background density coming from the first term in (2.16) thus completely screens the K point charges outside of the set K . Note that for (2.23) to hold we must have neutrality, The sets K are constructed in our approach via the minimization of an auxiliary functional bearing some formal resemblance to a 2D version of the usual Thomas-Fermi theory for molecules [40]. Accepting their existence, the rest of the proof has two main steps, going roughly as follows -Given the screening region K associated to a subset z 0 1 , . . . , z 0 K of the minimizing configuration, no other point of a minimizing configuration of particles can lie within K . We refer to this as the exclusion property. It holds for ground states of (2.16), as we shall prove, because the third term herein is superharmonic.
-Any configuration of points satisfying the exclusion property must have its local density bounded above by (π ) −1 . If some region contains a larger density, consider a disk included in it. The screening region associated to the points contained in the disk would have, because of (2.24), to leak outside of the disk. It would then overlap with other points of the configuration, which is impossible because of the exclusion property.
The construction of the screening set, via a variant of Thomas-Fermi theory, is done in Sect. 3. Some care is needed in handling the auxiliary variational problem, because it is of a non-standard type. We next prove the local density bound for ground states of (2.16) in Sect. 4, following the above strategy. The application to low-temperature Gibbs states (2.15) occupies Sect. 5, where we prove Theorem 2.1. Finally we deduce Corollary 2.3 in Sect. 6. Remark 2.6 (Recent related results). The classical Hamiltonian (2.16) and the associated Gibbs state (2.15) in the pure Laughlin case F = 1 have attracted much attention in their own right. Apart from the connection with the FQHE we are primarily concerned about here, one can see them as basic models for trapped plasmas. In this context, a generalization often considered (2D log-gas or 2D β-ensemble) is to replace |z j | 2 → V (z j ) in (2.22), with V : R 2 → R a general trapping potential. Minimizing configurations are then weighted Fekete sets [62], and Gibbs states for special values of the temperature are connected to certain random matrices ensembles [5,22].
For general V the density is not flat: the limiting density of points is proportional to V on its support. One expects the true minimizers and low-temperature Gibbs states to follow this profile on any microscopic length scale N −1/2 (in the scaled variables used in (2.22)). This has recently been proved in [4,53] for ground states (see also [50] for higher dimensional Coulomb and Riesz gases) and in [7,8,34,36] for Gibbs states. See also earlier partial results (not on the optimal scale) in [51,56,65] and [2,3,10,16,17,35,37,64] for related recent literature.
One can see our main result as a generalization of the fine-scale rigidity estimates of [4,7,34,53] to the case of a non-trivial F. Motivated by FQHE physics, this generalization is our chief concern. For general analytic F one should expect only local density upper bounds, but no corresponding lower bound in general. None of the available methods seem to apply here, but the new approach we sketched above might be of interest in the usual context with F = 1 and a general V in (2.22). Note that we obtain the expected density upper bound on essentially optimal length scales for minimizing configurations (see Proposition 4.6 below), but that it remains an open problem to do the same for Gibbs states, i.e. allow any α > 0 instead of only α > 1/4 in Theorem 2.1.

Incompressible 2D Thomas-Fermi Molecules
In this section we discuss a Thomas-Fermi-like variational theory of 2D molecules. We consider the functional where σ is the density distribution of negatively charged "electrons", is the 2D electrostatic potential generated by K point "nuclei" of charge + 1 located at positions x 1 , . . . , x K in the plane and is the electrostatic self-energy of the "electrons". We emphasize that in the application we have in mind, the designations "nuclei"and "electrons" do not refer to real physical entities. The terminology is just chosen for convenience and analogy with standard TF theory for molecules [15,40].
We consider the minimization of E TF over the variational set which means that we are interested in an "incompressible neutral molecule". The term "incompressible" refers to the constraint σ ≤ 1 on the density of electrons and "neutral" to the requirement that the total charge of the electrons equals that of the nuclei.
with a function τ representing the kinetic energy density of the electrons. For τ (σ ) = (const.) · σ 2 this term is the standard local density/semi-classical approximation of the quantum kinetic energy of a 2D electron gas. The minimization problem (3.1)-(3.4), which is the relevant one for our purpose, corresponds formally to taking τ (σ ) = σ p with p → ∞ to enforce the uniform upper bound σ ≤ 1, but we shall not follow this p → ∞ route. Many results we prove below have a natural extension to the case where V nuc is a more general potential, as long as it is generated by a compactly supported charge distribution. One could also consider corresponding 3D models as e.g. in [12], where results parallel to those of this section have been obtained. Note that the 3D setting is in fact simpler to deal with, because of the fall-off at infinity of the 3D Coulomb kernel.
We also mention that there is a kinship with recently studied flocking/swarming models [13,23] where constrained minimization problems similar to the above have been considered and proved to support phase transitions.
The facts we shall need about the variational problem for E TF are collected in the following theorem.
be the total electrostatic potential of the molecule. It is continuous and once continuously differentiable away from the nuclei and tends to zero at infinity. Moreover, for almost all x, and the area of The locations x 1 , . . . , x K of the nuclei lie within TF (x 1 , . . . , x K ) and all connected components of the latter set contain at least one x j .
(v) For a single nucleus, The inclusion in (3.11) says that when a nucleus is added, the density σ TF and the potential TF cannot decrease anywhere. The rest of this section is devoted to the proof of these results.

Existence and uniqueness of a minimizer.
The uniqueness is simpler than the existence and we consider it first: Proof. We only have to show that D(σ, σ ) is strictly convex on M TF . This is a consequence of the positivity of D(μ, μ) for neutral charge distributions μ as in [58,Lemma 3.2]. The general fact is that whenever R 2 μ = 0, see [62, Chapter I, Lemma 1.8]. This implies the claimed strict convexity property: pick two measures μ 1 and μ 2 with μ 1 = μ 2 . Then Equality holds if and only if μ 1 = μ 2 a.e. Uniqueness of the minimizer follows because we minimize a strictly convex functional on a convex set.
Next we consider the boundedness from below and lower semicontinuity of the functional E TF .

Lemma 3.4 (Boundedness from below).
The functional E TF is uniformly bounded from below on M TF . It is moreover lower semi-continuous under strong convergence in L 1 ∩ L p for any p > 1.
Proof. For any σ ∈ M TF we can use the normalization to write The function in parenthesis is clearly uniformly bounded from below outside of a disk D(0, R) that contains all the nuclei. On the disk and non-negative, we can split the integration domain according to whether |x| ≤ R or the other way around, and conclude that both pieces are bounded from below uniformly. The lower semi-continuity follows from the same kind of considerations, using in addition Fatou's lemma and the fact that log | . | ∈ L q loc (R 2 ) for any 1 ≤ q < ∞.
We now prove the existence of a minimizer. Note that if we had been minimizing under the constraint that the mass be equal to N < K , existence of a minimizer would follow from standard arguments discussed e.g. in [62]. At neutrality, when the number of electrons equals that of nuclei, this is slightly more subtle and we use some tricks we learned from [10,27,28].
Proof of Theorem 3.2, Item 1. Uniqueness of the minimizer is already contained in Lemma 3.3. For the existence, consider a minimizing sequence (σ n ) n∈N . In view of the definition of M TF we can extract a (not relabeled) subsequence which converges weakly-as a Radon measure, weakly in any L p (R 2 ) for any 1 < p < ∞, and weaklyin L ∞ (R 2 ) to a candidate minimizer σ TF , i.e.
for any function f ∈ L q (R 2 ), 1 ≤ q < ∞ (the convergence as a Radon measure corresponds to f continuous with compact support). Using Mazur's lemma [41, Theorem 2.13] and the convexity of the functional we could assume without loss that the convergence is strong in any L q (R 2 ), 1 < q < ∞, but we shall not use this. By the lower semicontinuity proved in Lemma 3.4 we obtain Hence we only have to prove that σ TF ∈ M TF . The upper and lower constraints on σ n pass to the limit easily, but it remains to show that This follows from the fact that the sequence (σ n ) n∈N is tight, i.e., no mass escapes to infinity in the limit. To prove this we pick some R > 0 and rewrite the energy as in (3.15), then split the part of the double integral where |x| ≥ R, |y| ≥ R from the rest: The second term is bounded below independently of n and R by the same arguments as in the proof of Lemma 3.4. For the rest we rewrite The function in parenthesis on the second line is bounded below uniformly for |x|, |y| ≥ R large enough, so that, once integrated against 1 |x|≥R 1 |y|≥R σ n (x)σ n (y)dxdy it gives a contribution uniformly bounded from below. For the contribution of the last term we note that All in all we deduce that, the energy of σ n being uniformly bounded above, where C depends neither on n nor on R > 0 provided the latter is large enough. Since log(1 + R 2 ) → +∞ for R → ∞, this implies tightness of the minimizing sequence, namely that for any ε > 0 there exists R ε > 0 such that To study the properties of the minimizer, whose existence and uniqueness we have now established, we next consider the corresponding electrostatic potential.

The Thomas-Fermi potential.
To any electron density σ we associate a potential (3.26) If R 2 σ = K then this is the Coulomb potential generated by a neutral charge distribution. We prove first that it is regular away from the nuclei and decays at infinity. The general case follow by cutting σ at some finite radius R and using the dominated convergence theorem as R → ∞.
For the asymptotic behavior we consider also first the case that σ has support in some disk D(0, R) with R < ∞. For |y| ≤ 1 2 |x| we have log |x − y| = log |x| + log 1 − y x ≤ log |x| + |y/x| 1 − |y/x| ≤ log |x| + 2 |y| |x| (3.28) where the ratio of two vectors in R 2 is interpreted as that of the corresponding complex numbers.
and the result is proved in this case.
Let us now only assume that as in the definition of M TF . This implies that for any ε > 0 there is a R ε such that |y|≥R ε σ (y) log(1 + |y|)dy < ε. (3.31) For every ε > 0 there is a δ > 0 such that log(1 − t) < ε if |t| < δ. Pick some large R and write for |x| ≥ R ≥ R ε (3.32) The first term is bounded by (const.) ε by the choice of δ and because R 2 σ is finite.
(3.34) because of (3.31). In the remaining term (remember that δ · R is large and |x| ≥ δ · R) we split the integral into a disk D(x, 1) of radius 1 around x and the rest. Since the logarithm is square integrable over the unit disk we can use the Cauchy-Schwarz inequality to conclude that The integral over the complement of D(x, 1) is also small because there log |x − y| ≤ 2 log |y|. Summarizing, we have proved that for any ε > 0 there exists R ε > 0 such that for |x| ≥ R ε . Returning to the first equality in (3.28) proves the lemma.
Proof. Equations (3.37)-(3.39) follow from local minimality of the functional (3.1) at σ TF by performing small variations in the standard way, c.f. [24,Lemma 4.1.2] or [13,23]. Here λ is the Lagrange multiplier associated with the mass constraint. That we get only inequalities on the sets {σ TF = 1} and {σ TF = 0} comes from the fact that the variational constraint 0 ≤ σ ≤ 1 allows only one-sided variations on these sets.
The Lagrange multiplier λ we shall refer to as a chemical potential for it is related 5 to variations of the minimal energy as a function of the electronic charge. We have the following: Proof. Assume λ < 0. By Lemma 3.5 TF vanishes at infinity so there is a radius R such that TF (x) > λ for |x| ≥ R. By (3.37)-(3.39) this implies σ TF (x) = 1 a.e. for |x| ≥ R, contradicting the finiteness of σ . Hence λ < 0 is excluded.
Assume λ > 0. As above, there is a radius R such that TF (x) < λ for |x| ≥ R. By Because TF is harmonic on {x : TF (x) < λ}, it cannot have a local minimum on this set so we know that TF (x) ≥ 0 everywhere and hence¯ (R) ≥ 0. By continuity of TF and the assumption λ > 0 there is a radius R 0 such that¯ (R 0 ) > 0, but¯ (R) < λ for R > R 0 . The circle with radius R 0 thus encloses the whole set {σ TF > 0}, up to a set of measure zero. On the other hand, by Newton's Theorem [41, Theorem 9.7],¯ (R 0 ) is proportional to the total charge inside the circle of radius R 0 which is zero and we have a contradiction.
Having excluded both λ < 0 and λ > 0 we conclude that λ = 0. Finally, {x : TF (x) < 0} is empty because TF would be harmonic on this set and zero on its boundary.
With Lemmas 3.6 and 3.7 we have established Item (2) of Theorem 3.2 We next show that σ TF takes only the values 0 or 1 almost everywhere.

Proof of Theorem 3.2, Item (3). We want to show that the set
has zero Lebesgue measure. By the definition (3.6) of TF we have in the sense of distributions and we sketch that argument here.
If a function f is in W 1,1 loc then on the set on which f has any fixed value, its gradient, ∇ f , vanishes except for a set of measure zero, see [1,46,68] or [41,Theorem 6.9]. Away from the nuclei, TF is explicitly seen to be continuous and in L 1 loc , and the same holds for its gradient ∇ TF . On A we thus have ∇ TF = 0. In order to repeat this argument for the Laplacian, we need to show that TF is in W 2,1 loc .
What we do know is that the sum of the second derivatives, TF , is an L 1 loc function (see (3.42)). What we need is that all mixed second derivatives are in L 1 loc . This follows by an exercise involving Fourier transforms, c.f. [23,Lemma 11]. Thus (3.43) does hold a.e. on A, and this contradicts (3.42) unless A has zero Lebesgue measure.
We now turn to the last item of Theorem 3.2, concerning the support of the electron density.

Proof of Theorem 3.2, Item (4).
Step 1. We first prove the estimate (3.8) on the support of σ TF . By Lemma 3.7 we know that λ = 0. There is an R such that |x i | < R for all i. For |x| ≥ R andR > R we introduce the function Then, for R < |x| ≤R, g ≤ 2π. (3.45) Since TF > 0 on this set, which contains no nucleus by assumption, we have σ TF = 1 and TF = 2π on S. Hence (g − TF ) ≤ 0 (3.47) on S, so the function g − TF is superharmonic on S and reaches its minimum at the boundary or at infinity. At infinity the function vanishes. For |x| = R we can make g − TF ≥ 0 by choosingR = R + M R (3.48) where M R is given by (3.9). We conclude that on S we have g − TF ≥ 0 and thus this set must be empty. Then TF ≤ g for |x| ≥ R, so TF = 0 and hence σ TF = (2π) −1 TF = 0 for |x| >R. This yields the desired result, Eq. (3.8).
Step 2. Next we consider the open set TF (x 1 , . . . , x K ) ≡ TF where the potential TF is strictly positive. From Eq. (3.7) we have TF ⊂ {σ TF > 0} while {σ TF > 0} is contained in the closure¯ TF . Since TF is continuous outside the "nuclei", TF is Jordan measurable which means that its boundary has Lebesgue measure zero. Hence TF and {σ TF > 0} are equal up to a null set. Since σ TF takes only the values 1 or 0, it follows that the area of TF is equal to K .
Step 3. Since TF (x) → +∞ when x approaches a location of a nucleus, we have x i ∈ TF for all i. Consider a connected component of TF Clearly, x i / ∈ U for all i = 1, . . . , K . Since TF K ≥ 0 we have K −1 > 0 on U and hence σ K −1 = 1 on U. Thus on U because σ K ≤ 1. Thus K −1 − K is subharmonic on U and U must be empty. We conclude that Step 5. For a single nucleus, uniqueness implies that the minimizer must be rotationally symmetric around x 1 . By Step 3 the support of the density is connected and hence a full disk with σ TF = 1. The Poisson equation then lead to the unique solution (3.53)

Local Density Upper Bounds for Classical Ground States
In this section we discuss the ground state of a generalized classical jellium Hamiltonian as presented in the introduction, and prove the announced local density upper bound on any microscopic scale 1. The argument is split in two steps, corresponding to Sects. 4.1 and 4.2: -First we prove an exclusion rule for ground states configurations: no point can lie inside a screening region defined by any subset of the other points. The latter is given by the open set TF where the TF potential TF corresponding to the chosen subset of the ground state configuration is strictly positive. -Second, we prove that any configuration of points satisfying the above exclusion rule must have a local density bounded above by 1, which is the desired optimal upper bound. In Sect. 4.3 we apply these results to ground states of the Hamiltonian (2.16) entering the plasma analogy, and small perturbations thereof.

Exclusion by screening.
The exclusion rule we shall chiefly rely on in the sequel is as follows: Exclusion by screening). Let X N = {x 1 , . . . , x N } ⊂ R 2 . Take any subset {y 1 , . . . , y K , y K +1 } ⊂ X N of distinct points y i . Define as in (3.10) TF (y 1 , . . . , y K ) to be the open set where the TF potential (3.6), generated by nuclei at locations y 1 , . . . , y K and the corresponding background charge distribution σ TF , is strictly positive. We say that X N satisfies the exclusion rule if y K +1 / ∈ TF (y 1 , . . . , y K ) (4.1) holds for any choice of y 1 , . . . , y K , y K +1 ∈ X N .
We call this "exclusion by screening" because TF defines a region surrounding y 1 , . . . , y K outside of which the background charge density σ TF totally screens the potential generated by y 1 , . . . , y K . Indeed, as we proved in the preceding section, Eq. (3.7), the total potential generated by the background and the point charges together vanishes identically outside TF . Note also that the closure of TF is the essential support of σ TF .
The relevance of this notion to our problem is as follows: Consider a Hamiltonian function on R 2N of the form   H(x 1 , . . . , x N with W symmetric and superharmonic in each variable x i . Equation  Proof. By symmetry of the Hamiltonian we may, without loss, choose Consider fixing all points but x 0 K +1 . The energy to consider is then We claim that if x ∈ TF (x 0 1 , . . . , x 0 K ) ≡ TF then there is anx ∈ ∂ TF such that G(x) < G(x). Thus the minimizing point x 0 K +1 cannot lie in TF . To prove the claim, we add and subtract a term − TF log |x − x |dx to write with a superharmonic function Now, is precisely the TF potential corresponding to "nuclear charges" at x 0 i , . . . x 0 K . Hence, using (3.7), > 0 on TF and zero on the boundary ∂ TF . The first two terms in R are harmonic on TF when taken together. (The Laplacian applied to the first term gives 2π and to the second term −2π on TF .) The other terms are superharmonic on TF . Thus, R takes its minimum on the boundary, so there is ax ∈ ∂ TF with R(x) ≥ R(x) for all x ∈ TF . On the other hand, (x) > 0 = (x) so G(x) > G(x).

Exclusion rule and density bounds.
We now show that any configuration of points satisfying the exclusion rule (4.1) has its local density everywhere bounded above by 1. To this end let R > 0 and define n(R) to be the maximum number of nuclei that a disk D(a, R) of radius R can accommodate while respecting the exclusion rule: It is clear that n(R) is independent of a because a translation of the x i just translates the corresponding exclusion region. Our density bound is as follows:

Theorem 4.3 (Exclusion rule implies a density bound). We have
We first recall from [60] a simpler density bound based on an unpublished theorem of Lieb which is essentially (4.1) for the special case K = 1.  To deduce (4.12), pick a configuration of points {x 1 , . . . , x N } contained in a disk of radius R and smear each point over an open disk of radius 1/(2 √ π). Because of (4.11) these disks do not overlap and must all be contained in a slightly larger disk of radius where the left side is the area covered by the small disks and the right-hand side the area of the large disk. This gives the desired result.
The rough bound Equation (4.12) turns out to be quite useful for the proof of the full Theorem 4.3. The first step is a lemma saying that the lim sup in (4.10) is attained by sequences of configurations with no large vacancies:

Lemma 4.5 (Maximal configurations have no vacancies). Let R k be a sequence of radii
(4.14) Then, for any fixed ε > 0 and any sub-disk D(a, εR k ) ⊂ D(0, R k ), This means that, for a configuration achieving the tightest packing in a given sequence of balls, the density is asymptotically uniform on any length scale comparable to the size of the balls.
Proof. The main idea is that a density lower than max in some sub-disk would have to be compensated for by a density higher than max in another sub-disk. This would contradict the definition of max as the lim sup of the maximal density achievable in any sequence of disks with radii tending to ∞.
Note first that by Lemma 4.4 we have for any R, as in (4.13), Hence max is finite and, in fact, We now tile the disk of radius R k with disjoint smaller disks of radius s := ε R k < 1, labeled by an index i = 1 . . . I . We can achieve this leaving only an area O(ε 2 R 2 k ) untiled. The latter we tile with disjoint still smaller disks of radius t := ε 1/2 R 1/2 k , labeled by j = 1 . . . J . In this way we leave only a domain A h with area |A h | = O(ε R k ) uncovered. Let n (i, s), respectively n ( j, t), be the number of points of the configuration X N k in each of the the smaller and the still smaller disks and let n h be the number of points in A h . By definition of max we have In the second line we used the definition of max as a lim sup and that for each fixed ε > 0, ε R k → ∞. Moreover, we used Lemma 4.4 to bound the density in A h . In the third line we simply used that we tile with disjoint sets. Since by assumption For large R k and fixed ε > 0, the small disk labeled i = 1 must thus satisfy (4.15). For any disk as in the statement, we can construct a tiling such that it is the first small disk and the result follows.
We now describe briefly the main arguments leading to the proof of Theorem 4.3. The area of the "neutralizing region" TF associated with the nuclear charges contained in a disk of radius R is equal to the number of charges in the disk, so the task is to show that this area is asymptotically equal to the area of the disk as R → ∞. On the boundary of the region TF the potential generated by the point charges in the disk and the TF neutralizing density σ TF must vanish because of the exclusion rule. In particular, the potential must vanish at any point of the nuclear charge configuration outside the disk. By the previous lemma, we can assume that the density of such vanishing points is bounded below uniformly, and combining this with a bound on the gradient of the potential we obtain an upper bound on the potential outside the disk. We can now proceed in two different ways (the second one was sketched in [42]): -The bound on the potential outside the disk together with (3.8) implies that we can enclose TF completely in a slightly larger disk whose radius behaves as R (1 + o(1)). Thus the area of TF , which is equal to the number of nuclei in the disk increases at most like the area of the disk plus a small correction. -Using Newton's theorem we obtain a bound on the circular average of the potential which involves the deviation from neutrality for the charge distribution within the disk. A comparison with the previous upper bound again leads to the desired result: Asymptotically the positive charge cannot be larger than the negative charge in the disk.
The details are as follows, the first steps being common to both arguments.
Proof of Theorem 4.3. Consider any maximal density sequence of radii and configurations defined in the previous Lemma 4.5. We use the lemma in two ways. First, taking ε = 1/2, we conclude that the sequence of radii r k := R k /2 and the corresponding number Second, for every ε > 0, the lemma implies that every point in the annulus is at most a distance O(εr k ) from some point in X N k ∩ A as r k → ∞. Let σ TF k be the TF density defined by the point charges {x 1 , . . . , x n k } ⊂ D(0, r k ) and the corresponding TF potential. From Prop. 4.2 we know that TF k vanishes at all points in X N k ∩ A. Moreover, we know that Pick some 1 2 > δ > 0. For |x| ≥ (1 + δ)r k a simple estimate gives a bound on the gradient of the potential: Indeed, if |x| ≥ (1 + δ)r k then |x − x j | ≥ δ · r k for j = 1, . . . , n k and thus Also, splitting the integral into two regions where |x − y| ≥ α = √ n k or vice-versa, because σ TF k = n k and we know from (4.21) that n k ≥ Cr 2 k . Since every point in the annulus A is at most a distance O(εr k ) from a point where TF k vanishes, the gradient estimate (4.25) implies a bound on the potential in the annulus: where in the second inequality we used that n k is in any case smaller than 4r 2 k (1 + o(1)) for large k by Lemma 4.4. With (4.28) at our disposal, we can conclude in two different ways: First proof. By Eq. (3.8) the bound (4.28) implies a bound on the largest radiusr k of a disk D(0,r k ) containing the support of σ TF k : (1)). Choosing δ = ε 1/3 to optimize the above we obtain (1)). (4.30) Since 0 ≤ σ TF k ≤ 1 this implies and thus, since ε > 0 is arbitrary, Thus, picking some fixed δ > δ > 0, because M(r ) ≤ πr 2 is an increasing function. By construction we have, for large k, n k ≥ πr 2 k max (1 + o(1)) and sō Combining with (4.28) we deduce for any ε > 0 and δ > δ > 0 where the o(1) goes to zero when k → ∞. We may thus take the limits k → ∞, ε → 0, δ → 0, δ → 0, in this order to deduce that max ≤ 1 as desired.

Applications.
After scaling, x → z = √ π x, Theorem 4.3 applies to the Hamiltonian for any superharmonic function W of the N variables. To obtain a corresponding result for low-temperature Gibbs states, we shall later use a Feynman-Hellmann argument and thus need to obtain bounds for a perturbed version of (4.33): where ε > 0 is a small enough number.
Proposition 4.6 (Density bound for perturbed plasma ground states). Assume U ∈ C 2 (R 2 ) with U uniformly bounded on R 2 . Let Z 0 N = (z 0 1 , . . . , z 0 N ) be a minimizing configuration for (4.34) and Then, as N → ∞ Proof. For ε = 0 and a disk this is just a combination of Proposition 4.2 and Theorem 4.3, together with the scaling x → √ π z. For the perturbed Hamiltonian, we add and subtract ε 4 and use the fact that U (z) − 1 4 U ∞ |z| 2 is superharmonic. Hence can be absorbed in W and we only need to reproduce the proof by writing (4.34) in the manner To obtain the result for a general domain r , we argue by covering it with disks, as in the "Cheese Theorem" [43,Section 14.4] which quantifies how effectively a region can be approximately covered with non-overlapping disks. We tile r with disks of radius √ r , leaving an open set of area at most O(r ) untiled. In the disks, we argue as previously. In the remaining untiled area we can use the rougher bound on the minimal separation of points, Lemma 4.4 to show that it contains at most O(r ) points. This is of smaller order than the main contribution to (4.36), which comes from the disks, whose union contains O(r 2 ) points, and the result follows.

Local Density Bounds for Classical Gibbs States
Here we prove our main result Theorem 2.1. We proceed in two steps which are roughly as follows: we first obtain the result for fully-correlated states having angular momentum bounded by C N 2 for some C > 0. These have the bulk of their density contained in a region that is not vastly larger than the extension of the Laughlin state. next we eliminate the a priori angular momentum bound by a localization procedure.

Conditional local density bound.
In order to combine neatly with the localization procedure used in Sect. 5.2 below, we must work in a slightly more general setting than discussed in Sect. 2.3. Let μ be a N -body probability density that we can write in the form where L > 0 is a number and log G is subharmonic (see e.g. [41,Chapter 9] for equivalent characterizations): G ≥ 0 and j log G ≥ 0 for any j = 1 . . . N . (5.2) Note that G need not be the modulus of an analytic function here. It could, in fact, be of the form G = |F 1 | + |F 2 | with F 1 and F 2 analytic, and we shall use this later, see Lemma 5.2 below. Marginal probabilities of μ will be denoted as In this subsection we prove the following on the first marginal density: for any α < α where |D| is the area of the disk and o D (1) tends to zero as |D| → ∞.
To see how this relates to our main result, think of the case where μ = | F | 2 is the density of a fully-correlated state, that is where L = +∞ and G = |F| with F analytic. Assuming a total angular momentum of order N 2 Proof. The method is basically the same as in [60,Section 3.1]. We first write μ as a Boltzmann-Gibbs factor: with Z N a normalization constant. Thus μ minimizes the free-energy functional ν log ν (5.8) over all N -particles probability densities ν. We denote by E 0 N = − log Z N the minimum value.
For r, δ > 0 let U r denote the negative characteristic function of a disk with radius r . Let U r,δ be a regularization over a distance δ so that We may choose the regularization in such a way that Consider now a perturbed version of the above classical Hamiltonian H N : along with the associated free-energy functional ν log ν (5.12) whose minimum we denote by E ε N . By definition we have and we will obtain our density estimate from upper and lower bounds to this free energy.

Free-energy upper bound.
We need an upper bound on E 0 N . We use the trial state implies that it cannot increase upon taking an average over disks centered at the z 0 i 's (see e.g. [41,Chapter 9]). For the one-body term we have, integrating in polar coordinates and by the same token, since L ≥ 1 Fixing η we thus deduce Free-energy lower bound. By positivity of the relative entropy (Jensen's inequality) we have, for any pair of probability measures ν, ν t , we deduce that for any probability measure where μ 0 ε is the empirical measure of a minimizing configuration for H ε N . Conclusion. Combining with (5.20) we deduce But, by (5.10) we have Proposition 4.6 applies to μ ε 0 because it is the empirical measure of a minimizing configuration for (5.7) (compare with (4.34)). All in all we thus have D r μ (1) where we have used (5.9) and assumed that r + δ 1. Consider now the case that r = N α with α > 1/4. We can then clearly choose ε, δ depending on r so as to have which implies N −1 r 2 ε −1 log N and 1 εδ −2 (5.32) and thus (5.4) follows.

Localization procedure.
We now want to extend Theorem 5.1 by dropping the angular momentum constraint. We use a rather standard localization procedure, see e.g. [38,54] for its discussion in more general contexts. In order for the localization method to combine efficiently with what we have proved so far, we need the following elementary lemma: Lemma 5.2 (The PL class). The class of functions G satisfying (5.2) is closed under addition. Thus any function G that can be written in the manner G = j α j |F j | with α j positive numbers and F j analytic functions satisfies (5.2).
The rest of the statement follows because any |F| with F analytic satisfies (5.2) as remarked after Eq. (2.20). See also [30,Page 84].
We now proceed to extending Theorem 5.1 to states not necessarily satisfying (5.5). This will conclude the Proof of Theorem 2.1.
Localization procedure. We define a smooth partition of unity on R 2 with L → ∞ a length scale to be optimized over. Consider the classical M-particle state 6 μ M built from μ by localizing exactly M particles and integrating out the other: (5.37) Observe that, using N j=1 χ 2 (z j ) + η 2 (z j ) = 1 and symmetry under particle exchange,  (1) F = N −1 ρ F the one-particle probability density of F . By symmetry under particle exchange again, we have M (z 1 ). (5.40) Finally, since clearly |z − a| 2 χ 2 (z) ≤ L 2+δ χ 2 (z) for arbitrary δ > 0 we have Application. Pick a disk D of center a ∈ R 2 and radius r = N α with α > 1/4 and choose and we expect the right-hand side to be at most of order N −1+2α . To estimate it, we start from (5.40) and first brutally get rid of terms with small M for they have too few particles to contribute to the density at this scale: where we used that and inserted (5.39). We pick γ < 2α to render the first term N γ −1 in the above much smaller than the expected N −1+2α , so that we may focus on the second one.
We thus need to bound the 1-particle probability density μ (1) M , first marginal of the M-particles state μ M . Recalling (5.37) we can write 46) where (c M is a positive constant ensuring normalization of (5.46)) At fixed z M+1 , . . . , z N , the integrand in the above is the squared modulus of an analytic function of the variables z 1 , . . . , z M (they appear only in the second line) and is thus subharmonic. Since the PL class is closed under addition, c.f. Lemma 5.2, we deduce that G M satisfies Assumption (5.2). We can thus apply (5.4) to the normalized version of μ M and deduce for any α < α. Next, we sum the above inequalities from M = N γ to N and insert the result in (5.44) to obtain where we have used that we sum only terms where M 1 to ensure that the o D (1) terms are indeed small. We next insert (5.39) in the first sum and use (5.45) in combination with (5.41) to deal with the second one, obtaining To optimize over γ we take γ = 1 2 + 2α + (1 + δ)β 2 + 2α (5.51) and since we are at liberty to choose α , β arbitrarily close to α and δ arbitrarily small, we can render γ arbitrarily close to 1/2, so that the error terms in (5.50) become O(N −1/2+ε ) for arbitrarily small ε > 0. For α > 1/4 we deduce

Energy in Confining Potentials
We now provide the Proof of Corollary 2.3. This is a variation on the arguments of Sect. 5.1, but it is more convenient to work with scaled variables as in [60,Section 3]. We thus write, with the notation of Assumption 2.2 with the scaled one-particle probability density μ (1) We will deduce the result from the lower bound and estimate (recall the notation (2.21)) Since μ F is the Gibbs state at temperature 1/N associated withH ε N and μ (1) F its first marginal we havẽ 0 N + C(log N + 1) (6.7) using a trial state similar to (5.14) to bound the minimal free energy at temperature 1/N associated withH 0 N , exactly as in [60, Section 3.1]. On the other hand, by Jensen's inequality where we used (2.11) to bound R 2 e −U L . Thus is the empirical measure of a configuration (z 0 1 , . . . , z 0 N ) minimizing the classical Hamil-tonianH ε N (z 1 , . . . , z N )− 1 N N j=1 U L (z j ). We can apply a (re-scaled version of) Proposition 4.6 toμ 0 ε . Thus, for any domain r obtain from a fixed Lipschitz by scaling lengths by a factor r = N α , α > −1/2 where we used (2.9) to estimate U L ∞ . By a Riemann sum approximation on a grid of side length N α we then have whereμ 0 ε is a piecewise constant function satisfying 0 ≤μ 0 ε ≤ 1 π 1 + CεN 1−2α (1 + o(1)) (6.13) and R 2μ 0 ε = 1. (6.14) Since, using (2.9), ∇U L L ∞ ≤ C N 1/2−α we deducẽ for any fixed α > −1/2. Comparing with (6.7) we obtain (1)), where we used that U ≥ U L . Since α > 0 and α can be chosen arbitrarily close to −1/2 in the above, we can choose ε = N −δ with some δ less than, but sufficiently close to, 1 to get The continuity of the bathtub energy as a function of the upper density constraint we have just used is a consequence of the explicit solution of the minimization problem [41, Theorem 1.14]. To conclude the proof of (6.3) we note that, under our assumptions on U , it suffices to take L large enough but fixed to have This follows from the fact that minimizers of the bathtub energy are equal to either 0 or (π ) −1 almost everywhere, and that for potentials increasing at infinity, their support is compact.

Conclusions
By employing a new theorem about multi-particle screening, we have derived sharp local upper bounds on the particle density in ground states of a classical 2D Coulomb gas perturbed by arbitrary many-body potentials which are superharmonic in every variable. We have also derived corresponding bounds for Gibbs states which, via Laughlin's plasma analogy, have direct applications to quantum many-particle states in strong magnetic fields exhibiting the Fractional Quantum Hall Effect (FQHE) and which are fully correlated in order to minimize a repulsive interaction between the particles. Our results hold asymptotically for large particle numbers, on length scales that are small compared to the full system's extension. We expect that there is room for improvement as regards the smallest length scale we can afford. This remains an open problem. The bounds for the classical Coulomb ground states rest on two pillars: First, an "exclusion rule" for ground state configurations which prevents points in such a configuration from entering a screening region defined by the other points in the configuration. Second, a universal density bound that we prove for any configuration of points satisfying the exclusion rule.
For the Gibbs state, the modification of the ground state density caused by the entropy term in the free energy has to be estimated. This is done by perturbing the classical Coulomb Hamiltonian locally and bounding the perturbed free energy, obtaining a density bound via the Feynman-Hellmann principle. This procedure is first carried out under the assumption of an a-priori bound on the angular momentum and then extended to the general case by a localization argument.
Besides the density bounds we also derive lower bounds for the energy of fully correlated FQHE states confined in external potentials. In combination with recently obtained upper bounds [61], this establishes that approximate quantum mechanical ground states in the class of strongly correlated FQHE wave-functions can be obtained by generating uncorrelated quasi-holes on top of the Laughlin state.