Free energy expansion of determinantal Coulomb gases in the quadratic fields with a point charge

This paper derives the explicit free energy expansion up to the constant term for determinantal Coulomb gases in quadratic fields with a point charge, identifying the constant term with the Liouville action and utilizing a deformation framework combined with the foliation flow method to extend isotropic results to anisotropic settings.

The Gist

The Big Picture: A Dance of Particles

Imagine a crowded dance floor (the complex plane) filled with thousands of tiny, energetic dancers (particles). These dancers have a very specific rule: they really dislike being too close to one another. They push each other away, like magnets with the same pole facing each other. This is what physicists call a Coulomb gas.

However, the dance floor isn't empty. There is a "music" playing (an external potential) that tries to pull the dancers toward the center or shape them into a specific formation. The paper studies what happens when you have a huge number of these dancers ($N$) and you want to predict the total "energy" or "effort" of the whole system as the crowd gets infinitely large.

The Special Ingredients

The authors are looking at a very specific type of dance floor with two unique features:

1. The Elliptic Shape (The Anisotropy): Usually, the music pulls dancers equally in all directions, forming a perfect circle. But in this paper, the music is "stretched." It pulls harder in one direction than the other, turning the circle into an ellipse. The parameter $\tau$ controls how stretched this ellipse is.
2. The Point Charge (The VIP): There is a special "VIP" standing at a specific spot ($a$) on the floor. This VIP has a strong gravitational pull (a logarithmic singularity) that attracts the dancers. The strength of this pull is controlled by $c$.

The Three Ways the Crowd Can Arrange

Depending on how strong the VIP is ($c$), how far away they stand ($a$), and how stretched the floor is ($\tau$), the crowd forms three different shapes (called "droplets"):

• Regime I (The Donut): The crowd forms a ring with a hole in the middle. The VIP is inside the hole, and the dancers surround them but don't touch the center.
• Regime II (The Solid Blob): The crowd forms a solid, filled-in shape (like a squashed circle). The VIP is either outside the crowd or the hole has been filled in.
• Regime III (The Two Islands): The crowd splits into two separate, disconnected islands. (The authors note this paper focuses on the first two shapes, not the split islands).

The Main Goal: Counting the Energy

The authors want to calculate the Free Energy of this system. Think of Free Energy as the "total cost" of organizing this massive dance.

They are looking for a formula that predicts this cost as the number of dancers ($N$) goes to infinity. They know the cost is made up of several layers:

• The Big Layer ($N^2$): The main cost, which grows very fast.
• The Middle Layer ($N \log N$): A secondary cost.
• The Small Layer ($N$): A smaller correction.
• The Tiny Layer ($\log N$): Even smaller.
• The Constant Layer ($O(1)$): The final, tiny adjustment that doesn't change with the number of dancers.

The Breakthrough: While previous researchers could calculate the big layers, this paper successfully calculates the Constant Layer (the final tiny adjustment) for this specific, stretched, VIP-influenced scenario.

The Secret Sauce: How They Did It

To find this final number, the authors used a clever trick called Deformation.

Imagine you have a complex, knotted rope (the current system with the VIP and the stretch). It's hard to untangle and measure directly. Instead, the authors slowly "morphed" the rope:

1. They slowly moved the VIP to a different spot.
2. They slowly un-stretched the floor until it was a perfect circle again.

By tracking how the "cost" changed during these slow movements, they could work backward to find the exact cost of the original, complicated shape.

The Mathematical Tools:

• Orthogonal Polynomials: They used a special set of mathematical "rulers" (polynomials) that are perfectly balanced against the crowd's arrangement. By looking at the first few numbers (coefficients) of these rulers, they could deduce the total energy.
• Liouville Action: This is a fancy geometric term they use to describe the "shape cost." They found that the final constant term in their energy formula is directly linked to this geometric shape cost. It's like saying the final price tag of the dance depends on the curvature of the dance floor's edge.

Why This Matters (According to the Paper)

• Connecting Geometry and Physics: The paper shows that the tiny, constant part of the energy isn't just a random number; it's deeply connected to the geometry of the shape the particles form.
• A New Map: They created a new method to solve these problems that doesn't rely on the old, heavy tools (like Riemann-Hilbert problems) used in simpler cases. Instead, they used a "foliation flow" method, which is like tracing the flow of water over a landscape to understand its shape.
• Random Matrices: The results also help predict the behavior of "characteristic polynomials" in elliptic random matrices (a type of complex number grid used in physics and engineering).

What They Didn't Do

The paper explicitly states that they did not solve the case where the crowd splits into two separate islands (Regime III). They also did not apply these results to clinical uses or specific engineering devices; the work remains purely theoretical, focused on understanding the mathematical behavior of these particle systems.

In a nutshell: The authors figured out the exact "price tag" for a massive, stretched-out crowd of repelling particles with a VIP guest, by slowly morphing the system into a simpler shape and using advanced geometry to track the changes.

Technical Summary

Technical Summary: Free Energy Expansion of Determinantal Coulomb Gases in Quadratic Fields with a Point Charge

Problem Statement
This paper investigates the large-$N$ asymptotic expansion of the free energy (log-partition function) for a determinantal Coulomb gas in the complex plane. The system is governed by an external potential $Q(z)$ that combines an anisotropic quadratic field with a logarithmic point charge insertion:
$$Q(z) = \frac{1}{1-\tau^2}(|z|^2 - \tau \text{Re } z^2) - 2c \log |z-a|,$$
where $\tau \in [0, 1)$ represents the anisotropy (non-Hermiticity), $c \ge 0$ is the charge strength, and $a \ge 0$ is the location of the point charge. The study focuses on regimes where the associated equilibrium droplet (the support of the macroscopic density) is either simply connected (Regime II) or doubly connected (Regime I). The primary objective is to derive the free energy expansion up to and including the constant term ($O(1)$), explicitly computing all coefficients, and identifying the constant term with the Liouville action associated with the droplet geometry.

Methodology
The authors employ a deformation framework that relates variations of the free energy to the refined asymptotics of planar orthogonal polynomials. The methodology diverges from previous approaches in the isotropic case ($\tau=0$) by avoiding the Riemann–Hilbert method in favor of the "foliation flow" theory developed by Hedenmalm and Wennman.

The proof strategy proceeds in three main stages:

1. Deformation of the Free Energy: The authors establish differential formulas (Proposition 2.3) linking the derivatives of the free energy with respect to the parameters $a$ and $\tau$ to the leading and subleading coefficients ($A_{n,N}$ and $B_{n,N}$) of the monic planar orthogonal polynomials $P_{n,N}(z)$. Unlike previous works that relied on isomonodromic tau-functions or complex Riemann–Hilbert analysis, this paper provides an elementary proof based on the orthogonality of polynomials and one-point observables.
2. Asymptotic Analysis of Orthogonal Polynomials: Utilizing the general theory of planar orthogonal polynomials for algebraic Hele-Shaw potentials, the authors derive explicit asymptotic expansions for the coefficients $A_{n,N}$ and $B_{n,N}$.
   • In the simply connected regime, the analysis relies on the first subleading correction term ($F_1$) in the Hedenmalm–Wennman expansion, which is explicitly computed in terms of conformal geometric functionals (Liouville action and curvature).
   • In the doubly connected regime, a "hole-filling" reduction (Proposition 6.2) is used. This observation allows the authors to relate the orthogonal polynomials of the multiply connected droplet to those of a simply connected reference domain (an ellipse), leveraging the fact that the outer boundary remains fixed. This reduces the problem to the asymptotics of Hermite polynomials.
3. Integration and Reference Models: The free energy is reconstructed by integrating the derived variations along specific paths in the parameter space $(a, \tau)$. The integration constants are determined using exactly solvable reference models: the induced Ginibre ensemble (for $a=\tau=0$) and the elliptic Ginibre ensemble (for $a \to \infty$).

Key Contributions and Results

• Explicit Free Energy Expansion: The paper derives the full expansion of $\log Z_N(Q)$ up to the constant term $C_5$ for both simply and doubly connected droplets. The expansion takes the form:
  $$\log Z_N(Q) \sim C_1 N^2 + C_2 N \log N + C_3 N + C_4 \log N + C_5 + O(N^{-1}).$$
  All coefficients are computed explicitly in terms of the parameters $a, c, \tau$.
• Identification of the Constant Term: A central result is the identification of the constant term $C_5$ with the Liouville action $L[K_Q]$ of the droplet, alongside topological terms involving the Euler characteristic $\chi$ and the spectral determinant. Specifically, the term $C_5$ is shown to be:
  $$C_5 = \chi \zeta'(-1) + \frac{\log(2\pi)}{2} + \frac{1}{24}L[K_Q] - \frac{1}{24\pi} \int_{\partial K_Q} \log(\Delta Q) \kappa \, ds.$$
• Generalization of the Isotropic Case: The results extend the findings of Byun et al. (specifically [39] in the paper's bibliography) from the isotropic case ($\tau=0$) to the anisotropic elliptic case. The paper demonstrates that the deformation approach is robust enough to handle the loss of radial symmetry and the breakdown of duality relations that exist in the isotropic setting.
• Connection to Characteristic Polynomials: The results provide the large-$N$ asymptotics for the moments of characteristic polynomials of the elliptic Ginibre ensemble, given by the ratio of partition functions $Z_N(Q)/Z_N(Q_e)$.

Significance and Claims
The authors claim that this work establishes a general framework connecting free energy expansions, refined asymptotics of planar orthogonal polynomials, and conformally invariant geometric functionals. Key aspects of their contribution include:

• Methodological Shift: The paper moves away from the Riemann–Hilbert approach, which was dominant in isotropic studies, to a deformation-based method utilizing the foliation flow theory. This allows for a systematic tracking of the underlying conformal geometry.
• Elementary Proof of Deformation Formulas: The derivation of the deformation formulas (Proposition 2.3) is presented as an elementary proof based on orthogonality, contrasting with the technical Riemann–Hilbert constructions used in prior literature.
• Geometric Interpretation: The work explicitly links the constant term of the free energy to the Liouville action, offering a transparent geometric description that avoids the delicate regularisation required for spectral determinants on unbounded domains.
• Limitations and Scope: The authors modestly note that their results are currently restricted to simply and doubly connected droplets. They acknowledge that the multi-component regime (Regime III) remains open, as the current orthogonal polynomial techniques do not apply, and a different approach (potentially involving Riemann–Hilbert analysis) would be required. Furthermore, while the method is suggested to be extensible to general algebraic Hele-Shaw potentials, the full implementation for a broad class of potentials requires further development of deformation paths and finer asymptotic expansions.

The paper concludes that the developed strategy successfully bridges the gap between random matrix theory and conformal geometry, providing a concrete realization of the conjectured free energy expansion in a non-trivial, anisotropic setting.