Canonical and grand canonical partition functions of Dyson gases as tau-functions of integrable hierarchies and their fermionic realization

The partition function for a canonical ensemble of 2D Coulomb charges in a background potential (the Dyson gas) is realized as a vacuum expectation value of a group-like element constructed in terms of free fermionic operators. This representation provides an explicit identification of the partition function with a tau-function of the 2D Toda lattice hierarchy. Its dispersionless (quasiclassical) limit yields the tau-function for analytic curves encoding the integrable structure of the inverse potential problem and parametric conformal maps. A similar fermionic realization of partition functions for grand canonical ensembles of 2D Coulomb charges in the presence of an ideal conductor is also suggested. Their representation as Fredholm determinants is given and their relation to integrable hierarchies, growth problems and conformal maps is discussed.

💡 Research Summary

The paper establishes a deep connection between the statistical mechanics of two‑dimensional Coulomb gases (Dyson gases) and integrable hierarchies, using a fermionic operator formalism. Starting from the canonical ensemble of N identical charges confined to the complex plane and interacting via a logarithmic potential, the authors introduce free fermionic creation and annihilation operators ψ_n and ψ_n† (n∈ℤ) together with the current operators J_k = Σ_n :ψ_n ψ†{n+k}:. By expanding the external background potential V(z) in a Laurent series V(z)=∑{k>0}(t_k z^k + s_k \bar z^k), they construct a group‑like element

 G(V)=exp(∑{k>0} t_k J_k) · exp(∑{k>0} s_k J_{‑k})

and show that its vacuum expectation value reproduces the canonical partition function:

 ⟨0| G(V) |0⟩ = Z_N(V).

This identity identifies Z_N(V) as a τ‑function of the two‑dimensional Toda lattice hierarchy. Consequently, the Hirota bilinear equations governing the Toda hierarchy become equivalent to the recursion relations satisfied by the Dyson gas partition function, providing an algebraic integrable‑systems perspective on the Coulomb gas.

The authors then examine the dispersionless (quasiclassical) limit, where N→∞ while the effective Planck constant ℏ→0. In this regime the logarithm of the τ‑function converges to a harmonic function Φ(z, \bar z) whose level set defines an analytic curve Γ in the complex plane. The coefficients {t_k, s_k} become the moments (or “times”) of Γ, establishing a one‑to‑one correspondence between the Toda flows and deformations of the curve. This reproduces the well‑known inverse potential problem: given V(z), the equilibrium charge density is encoded in the shape of Γ, and the curve itself can be generated by a parametric conformal map from the unit disk. Thus the integrable structure of the Dyson gas is precisely the integrable structure underlying the theory of analytic curves and conformal mappings.

In the grand‑canonical setting, the system is coupled to an ideal conductor that creates image charges. The authors derive a kernel K(z,w) that captures both particle–particle and particle–image interactions. The grand‑canonical partition function can be written as a Fredholm determinant

 Z_{GC}=Det(1+K),

which, again, admits a fermionic representation as a vacuum expectation value of a suitably modified group‑like element. This shows that the grand‑canonical partition function is also a τ‑function of the Toda hierarchy. The Fredholm determinant form links the problem to growth processes such as Laplacian growth, the Laplace equation with moving boundaries, and stochastic Loewner evolution, where the evolution of the domain is governed by the same hierarchy.

Overall, the paper provides a unified framework that brings together random matrix theory, Coulomb gas statistics, integrable hierarchies, and complex analysis. By exploiting the free‑fermion realization, the authors obtain explicit operator expressions, clarify the role of dispersionless limits, and reveal how Fredholm determinants encode both canonical and grand‑canonical ensembles. The results open avenues for extending the approach to multi‑component gases, non‑isotropic potentials, and quantum Hall systems, where similar integrable structures are expected to emerge.