Multipole and Berezinskii-Kosterlitz-Thouless Transitions in the Two-component Plasma

We study the two-dimensional two-component Coulomb gas in the canonical ensemble and at inverse temperature $\beta>2$. In this regime, the partition function diverges and the interaction needs to be cut off at a length scale $\lambda\in (0,1)$. Particles of opposite charges tend to pair into dipoles of length scale comparable to $\lambda$, which themselves can aggregate into multipoles. Despite the slow decay of dipole--dipole interactions, we construct a convergent cluster expansion around a hierarchical reference model that retains only intra-multipole interactions. This yields a large deviations result for the number of $2p$-poles as well as a sharp free energy expansion as $N\to\infty$ and $\lambda\to0$ with three contributions: (i) the free energy of $N$ independent dipoles, (ii) a perturbative correction, and (iii) the contribution of a non-dilute subsystem. The perturbative term has two equivalent characterizations: (a) a convergent Mayer series obtained by expanding around an i.i.d.\ dipole model; and (b) a variational formula as the minimum of a large-deviation rate function for the empirical counts of $2p$-poles. The Mayer coefficients exhibit transitions at $\beta_p=4-\tfrac{2}{p}$, that accumulate at $\beta=4$, which corresponds to the Berezinskii-Kosterlitz-Thouless transition in the low-dipole-density limit. At $\beta=\beta_p$ the $p$-dipole cluster integrals switch from non-integrable to integrable tails. The non-dilute system corresponds to the contribution of large dipoles: we exhibit a new critical length scale $R_{\beta, \lambda}$ which transitions from $\lambda^{-(\beta-2)/(4-\beta)}$ to $+\infty$ as $\beta$ crosses the critical inverse temperature $\beta=4$, and which can be interpreted as the maximal scale such that the dipoles of that scale form a dilute set.