Higher-order terms, convergence, and quantitative validity of the high-activity $z^{-1/2}$ expansion for the hard square lattice gas with nearest- and next-nearest-neighbor exclusion

The $2\times2$ hard square lattice gas (particles excluding both nearest and next-nearest neighbors on $\mathbb{Z}^2$) has a high-activity columnar phase. Ramola and Dhar (2012) constructed a high-activity expansion of its free energy about a perfectly columnar-ordered state and summed the two-rod terms exactly, giving the series through $O(z^{-3/2})$, where $z$ is the particle fugacity. This problem concerns what later literature did with the three limitations they acknowledged. The $O(z^{-2})$ coefficient (three-rod diagrams) has still not been computed; no lower bound on the radius of convergence of the $z^{-1/2}$ series, nor even an analyticity statement for this sliding model's pressure, is known; and low-order truncations are quantitatively useful only for $z\gtrsim150$, far above the critical activity $z_c\approx97.5$. The problem of turning the two-rod calculation into a usable convergent high-activity theory of the columnar phase — coefficients, convergence control, and quantitative validation — remains open.

Background

A lattice gas consists of particles occupying sites of a lattice subject to hard-core exclusion rules; its only parameter is the fugacity $z$ (the weight per particle, $z=e^{\beta\mu}$), so all transitions are entropy-driven. The model here is the 'hard square' gas with exclusion of both nearest and next-nearest neighbors: each particle is a $2\times2$ block on the square lattice, so occupied blocks cannot touch. At maximal density the lattice is tiled by non-touching $2\times2$ blocks on a checkerboard of spacing 2; the close-packing density is $\rho=1/4$ particles per site (one particle per four lattice sites), consistent with the leading term of the series quoted below. Such perfect packings are exponentially numerous: on an $L\times L$ torus there are $4(2^{L/2}-1)$ of them, related by sliding whole rows or columns of blocks by one lattice spacing. This 'sliding' degeneracy is why the standard high-fugacity cluster techniques — Mayer expansions (the usual cluster expansions of statistical mechanics about a reference state) and Pirogov–Sinai theory (a general machinery for proving phase coexistence when the ordered ground states are finite in number) — fail here: a single vacancy can split into two half-vacancies that deconfine along a row or column, making the naive $1/z$ correction extensive in system size. Monte Carlo simulations establish a continuous transition at $z_c\approx97.5$ from a disordered fluid to a columnar phase, in which particles preferentially occupy every other row (or every other column) — four symmetry-related ordered states, whose continuous transition is described in the physics literature by Ashkin–Teller-type universality (the critical behavior of two coupled Ising models). Ramola and Dhar (2012) developed a perturbation expansion about a state with perfect columnar order. The elementary excitations are 'rods': a half-vacancy pair joined by a string of displaced squares, horizontal or vertical, of arbitrary length. Regrouping the double series in the two row-sublattice fugacities $z_A, z_B$ produces a Mayer-like series for a polydisperse gas of interacting vertical rods (the horizontal rods are effectively integrated out, generating a longer-ranged effective interaction between the vertical rods), in which the $n$-th term is of order $z^{-(n+1)/2}$ — fractional powers of $1/z$ are unavoidable for sliding systems. They summed all one- and two-rod terms exactly, obtaining, for example, the density $\rho(z)=\tfrac14-\tfrac18 z^{-1/2}-\tfrac14 z^{-1}-\big(\tfrac92\log\tfrac98+\tfrac{11}{64}\big)z^{-3/2}+O(z^{-2})$, which matches Monte Carlo data at $z=150$ and $z=200$ to three or four digits. Related rigorous work sharpens what is at stake: Jauslin and Lebowitz (2017–2018) proved, via an extension of Pirogov–Sinai theory, that for 'non-sliding' hard-core systems (finitely many perfect coverings, such as hard hexagons or hard crosses) the pressure has a convergent expansion in integer powers of $z^{-1}$ at large $z$; the $2\times2$ hard square model is the paradigmatic sliding counterexample, and one-dimensional analogues show the pressure may be analytic only in $z^{-1/2}$. Most recently, Hadas and Peled (2025) proved — by a novel extension of chessboard estimates to infinite volume plus disagreement percolation — that this very model has four extremal columnar Gibbs measures at sufficiently large fugacity, with anisotropic correlation decay; but their argument produces no analytic expansion and its fugacity threshold is an ineffective constant far above the physical regime. Thus, existence of the ordered phase is now rigorous, yet the analytic structure of the high-activity expansion itself — its coefficients beyond $O(z^{-3/2})$ and its radius of convergence — remains unsettled.

Problem Statement

For the hard square lattice gas with nearest- and next-nearest-neighbor exclusion ($2\times2$ blocks on $\mathbb{Z}^2$) in its high-activity columnar phase, develop the high-activity expansion of Ramola and Dhar beyond the two-rod level into a usable convergent theory. Concretely, the open problems left by that work and unresolved by subsequent literature are: (1) implement the regrouping of the expansion at higher order — an algorithmic enumeration and exact evaluation of all connected defect-cluster (rod) contributions to the regrouped vertical-rod series, beginning with the three-vertical-rod configurations that first contribute at $O(z^{-2})$ — so that the exact coefficients of the free energy, density, and order-parameter series in powers of $z^{-1/2}$ can be pushed beyond the currently known $O(z^{-3/2})$ terms; (2) establish convergence properties of the resulting fractional-power expansion: a lower bound on its radius of convergence in the variable $z^{-1/2}$ (or a determination that the pressure is not analytic in $z^{-1/2}$), together with an effective fugacity above which the series is guaranteed to represent the columnar phase; and (3) quantify, using the higher-order coefficients, how large $z$ must be for truncations of the series to reproduce the thermodynamics of the columnar phase, and how many fractional-order terms are needed, relative to Monte Carlo data and to the critical activity $z_c\approx97.5$. The three components are logically independent, and each states an open problem in its own right. Component (1) requires the exact $O(z^{-2})$ coefficients of the free-energy, density, and order-parameter series from the three-vertical-rod enumeration, with higher orders also in scope. These three components jointly serve one objective: converting the two-rod calculation into a controlled, convergent high-activity theory of the sliding hard-square columnar phase.

References

1. Kabir Ramola and Deepak Dhar, "A high density perturbation expansion for the hard square lattice gas", AIP Conference Proceedings 1447, 135–136 (2012), DOI: 10.1063/1.4709918, https://doi.org/10.1063/1.4709918
2. Kabir Ramola and Deepak Dhar, "High-Activity Perturbation Expansion for the Hard Square Lattice Gas", Physical Review E 86, 031135 (2012), DOI: 10.1103/PhysRevE.86.031135, arXiv:1204.4484, https://arxiv.org/abs/1204.4484
3. Trisha Nath, Joyjit Kundu, and R. Rajesh, "High-Activity Expansion for the Columnar Phase of the Hard Rectangle Gas", Journal of Statistical Physics 160, 1173–1197 (2015), DOI: 10.1007/s10955-015-1285-y, arXiv:1411.7831, https://arxiv.org/abs/1411.7831
4. Trisha Nath and R. Rajesh, "The high density phase of the k-NN hard core lattice gas model", Journal of Statistical Mechanics: Theory and Experiment 2016, 073203 (2016), DOI: 10.1088/1742-5468/2016-07-073203, arXiv:1604.08376, https://arxiv.org/abs/1604.08376
5. Ian Jauslin and Joel L. Lebowitz, "High-fugacity expansion, Lee-Yang zeros and order-disorder transitions in hard-core lattice systems", Communications in Mathematical Physics 364, 655–682 (2018), DOI: 10.1007/s00220-018-3269-7, arXiv:1708.01912, https://arxiv.org/abs/1708.01912
6. Ian Jauslin and Joel L. Lebowitz, "Crystalline Ordering and Large Fugacity Expansion for Hard-Core Lattice Particles", The Journal of Physical Chemistry B 122, 3266–3271 (2018), DOI: 10.1021/acs.jpcb.7b08977, https://doi.org/10.1021/acs.jpcb.7b08977
7. Daniel Hadas and Ron Peled, "Columnar order in random packings of $2\times2$ squares on the square lattice", Journal of the European Mathematical Society (2025), DOI: 10.4171/jems/1680, arXiv:2206.01276, https://arxiv.org/abs/2206.01276
8. N. Vigneshwar, Dipanjan Mandal, Kedar Damle, Deepak Dhar, and R. Rajesh, "Phase diagram of a system of hard cubes on the cubic lattice", Physical Review E 99, 052129 (2019), DOI: 10.1103/PhysRevE.99.052129, arXiv:1902.06408, https://arxiv.org/abs/1902.06408
9. Dipanjan Mandal, Trisha Nath, and R. Rajesh, "Estimating the Critical Parameters of the Hard Square Lattice Gas Model", Journal of Statistical Mechanics: Theory and Experiment 2017, 043201 (2017), DOI: 10.1088/1742-5468/aa650e, arXiv:1702.02332, https://arxiv.org/abs/1702.02332
10. Kabir Ramola, Kedar Damle, and Deepak Dhar, "Columnar Order and Ashkin-Teller Criticality in Mixtures of Hard Squares and Dimers", Physical Review Letters 114, 190601 (2015), DOI: 10.1103/PhysRevLett.114.190601, arXiv:1408.4943, https://arxiv.org/abs/1408.4943