Rotman and Eisenberg (2009) showed that the Mayer cluster coefficients of a hard-core lattice gas can be reorganized into a tridiagonal symmetric R matrix whose spectrum determines the Yang-Lee zeros and the critical singularities, and derived exact formulas for the singular exponents sigma (nonphysical singularity) and sigma-prime (= 1-alpha of the ordering transition) from the asymptotic approach of the R elements to constants, validated on the exactly solved hard-hexagon model. The paper itself identified the method's open limits: the strong-correction regime where the formulas give imaginary exponents because the spectral edge shifts (realized by the triangular-lattice N2 model, a four-state-Potts criticality), the contradiction between its square-lattice 4-NN prediction sigma-prime = 0.28(6) and Monte Carlo, the unproven reality of R for purely repulsive models, and, as a numerical caveat, the exponential error-sensitivity of the recursive construction. Later Monte Carlo work (Nath and Rajesh 2014) established that the 4-NN model has an Ising-class first transition plus a second two-color Ashkin-Teller transition, sharpening the explicitly open question of why the cluster analysis fails there. This record poses the surviving problem: complete the asymptotic theory - characterize the spectrum-edge-shift regime, diagnose or correct the multiple-transition failure, and settle the reality of R - so that the validity range of the only series method able to extract critical behavior of hard-core gases from low-activity expansions alone is established.

Background

Hard-core lattice gases - lattice models whose only interaction is excluded volume (each particle blocks a fixed set of neighboring sites), so there is no energy scale and any ordering transition is entropy-driven - are the standard discrete laboratories for freezing. Their equation of state in terms of the activity zz is the Mayer expansion ρ(z)=∑n≥1n bnzn,\rho(z)=\sum_{n\ge1} n\,b_n z^n, where bnb_n are the Mayer cluster integrals (cluster coefficients), computable exactly by transfer-matrix enumeration of finite clusters.

Rotman and Eisenberg (2009), extending Eisenberg and Baram (2006, 2007), reorganized this series into a spectral problem: one constructs a tridiagonal real symmetric matrix RR (a Jacobi matrix) whose powers reproduce the cluster coefficients, (Rn)11=(−1)n(n+1)bn+1(R^n)_{11}=(-1)^n(n+1)b_{n+1}, by an explicit element-by-element recursion. The density then becomes a sum over the eigenvalues λ\lambda of RR, ρ(z)=∑λψ1(λ)2z−1+λ,\rho(z)=\sum_\lambda \frac{\psi_1(\lambda)^2}{z^{-1}+\lambda}, so the reciprocals of the eigenvalues are the Yang-Lee zeros (the zeros of the grand partition function in the complex activity plane, which control the singularities of the series). In every purely repulsive model examined, the elements of RR converge rapidly to constants - BB on the diagonal, AA off the diagonal - with computable corrections, e.g. Rn,n=B+b/n2+b′cos⁡(qn)/nR_{n,n}=B+b/n^2+b'\cos(qn)/n and Rn,n+1=A+a/(n+1/2)2+a′cos⁡(q(n+1/2))/(n+1/2)R_{n,n+1}=A+a/(n+1/2)^2+a'\cos(q(n+1/2))/(n+1/2). The critical activities are the values of zz at which −z−1-z^{-1} meets the spectrum edges B±2AB\pm 2A: the physical ordering point zt=(2A−B)−1z_t=(2A-B)^{-1} and a nonphysical singularity −z0=−(2A+B)−1-z_0=-(2A+B)^{-1} on the negative real axis (the reason the standard series-extrapolation techniques (ratio analysis and Padé approximants, which extrapolate a truncated power series to its nearest singularity) fail near the transition). A renormalization-group-style decimation of the matrix then yields exact expressions for the singularity exponents of ρ(z)\rho(z), defined by ρ(z)−ρc∼(1−z/zt)σ′\rho(z)-\rho_c\sim(1-z/z_t)^{\sigma'} and its counterpart at −z0-z_0: σ′=121−4(2a−b)A−[2a′cos⁡(q/2)−b′]2[1−cos⁡(q)]A2,σ=121−4(2a+b)A−[2a′cos⁡(q/2)+b′]2[1−cos⁡(q)]A2.\sigma'=\frac{1}{2}\sqrt{1-\frac{4(2a-b)}{A}-\frac{[2a'\cos(q/2)-b']^2}{[1-\cos(q)]A^2}},\qquad \sigma=\frac{1}{2}\sqrt{1-\frac{4(2a+b)}{A}-\frac{[2a'\cos(q/2)+b']^2}{[1-\cos(q)]A^2}}. Here σ′=1−α\sigma'=1-\alpha, where α\alpha is the specific-heat exponent of the corresponding universality class, and σ\sigma describes the repulsive (Lee-Yang-type) nonphysical singularity. The method is benchmarked on Baxter's exactly solved hard-hexagon model (triangular lattice, nearest-neighbor exclusion), reproducing the exact σ′=2/3\sigma'=2/3 and σ=1/6\sigma=1/6 to within one per cent (the paper's fits give 0.6662 and 0.1655; Baxter 1980; Joyce 1988), and was then applied to the kk-NN family, in which a particle excludes all sites up to the kk-th nearest neighbors: the triangular-lattice N2 model (k=2k=2; established by several numerical studies to have a single continuous transition in the four-state Potts universality class, with σ′=1/3\sigma'=1/3 at μc=1.75682(2)\mu_c=1.75682(2), i.e. zt=eμc=5.794z_t=e^{\mu_c}=5.794; Zhang and Deng 2008, Jaleel et al. 2022) and the square-lattice N4 model (k=4k=4; at the time believed from Monte Carlo to have a single Ising-class transition with σ′=1\sigma'=1; Fernandes et al. 2007).

The source paper identified three unresolved limits of this analysis, together with one numerical caveat. (i) Strong-correction regime: when the fitted amplitudes are large enough that the square-root arguments above turn negative - formally imaginary σ\sigma or σ′\sigma' - the spectrum edge shifts away from 2A±B2A\pm B and the exponent formulas are invalid; the triangular N2 model realizes this (4(2a−b)/A≃2.7>14(2a-b)/A\simeq2.7>1, predicted zt=9.35z_t=9.35 versus the measured 5.7945.794), and no systematic treatment of this regime exists. (ii) The method's N4 prediction σ′=0.28(6)\sigma'=0.28(6) excluded the Ising class and contradicted Monte Carlo; the paper called for detailed numerical studies to settle the discrepancy. (iii) The entire construction presumes that RR is real for purely repulsive models - verified in every computed example but with no proof in general. (iv) Numerical caveat: the recursive construction of RR from the bnb_n is exponentially sensitive to coefficient errors; the source's Appendix A accordingly prescribes that the cluster integrals used be exact or at least known to high accuracy, and that the construction be carried out in high-precision arithmetic.

Problem Statement

Complete the asymptotic theory of the R-matrix (Mayer cluster-coefficient) method for hard-core gases with purely repulsive interactions, by resolving the three failure modes and gaps identified by Rotman and Eisenberg (2009), which together block the single objective of extracting the location, order, and critical exponents of entropy-driven ordering transitions from low-activity series alone.

(a) Strong-correction (spectrum-edge-shift) regime. Characterize analytically the spectral edge and the density singularity when the correction amplitudes are strong enough that the published exponent formulas would give imaginary σ\sigma or σ′\sigma': determine where the edge of the R-matrix spectrum moves, what replaces the exponent formulas in this regime, and a criterion delimiting their validity. The triangular-lattice N2 model, with its 17 published exact cluster coefficients, is the canonical test: a correct theory must be consistent with that model's established four-state-Potts criticality (zt=5.794z_t=5.794, σ′=1/3\sigma'=1/3), where the current analysis predicts zt=9.35z_t=9.35 and breaks down.

(b) Failure under multiple transitions. Diagnose and correct the method's failure for the square-lattice 4-NN model: the cluster analysis gives σ′=0.28(6)\sigma'=0.28(6), excluding the Ising class, while Monte Carlo now firmly establishes a first continuous Ising-class transition at μc=4.705\mu_c=4.705 followed by a second continuous transition in the two-color Ashkin-Teller class (a one-parameter universality family interpolating between two decoupled Ising models and the four-state Potts model). Determine which assumption fails - the single-dominant-singularity hypothesis, the fitted asymptotic regime, or the exponent formulas themselves - and provide the modified analysis, or a proof that the finite cluster series cannot decide the question. This is explicitly posed as open by Nath and Rajesh (2014).

(c) Reality of R. Prove, or exhibit a counterexample to, the assumption that the recursively constructed R matrix is real for every purely repulsive model.

Deliverables (a) and (b) are the physics core - the two concrete regimes in which the asymptotic exponent theory is known to fail - and (c) secures the mathematical foundation the whole construction presumes. A solution is complete when each is resolved within the formalism of the source paper (hard-core gases with purely repulsive interactions in the R-matrix representation, which applies equally to continuum and lattice models in all dimensions), and when the corrected theory is consistent with the hard-hexagon validation benchmark, where the exact exponents σ′=2/3\sigma'=2/3 and σ=1/6\sigma=1/6 are reproduced by the method's fits to within about one per cent (0.6662 and 0.1655).

References

  1. Z. Rotman and E. Eisenberg, 'Critical exponents from cluster coefficients', Phys. Rev. E 80, 031126 (2009). DOI: 10.1103/PhysRevE.80.031126; arXiv:0907.4037. https://doi.org/10.1103/PhysRevE.80.031126
  2. E. Eisenberg and A. Baram, 'Estimation of critical exponents from cluster coefficients and the application of this estimation to hard spheres', Proc. Natl. Acad. Sci. USA 104, 5755-5758 (2007). DOI: 10.1073/pnas.0700778104. https://doi.org/10.1073/pnas.0700778104
  3. E. Eisenberg and A. Baram, 'Analysis of the ordering transition of hard disks through the Mayer cluster expansion', Phys. Rev. E 73, 025104(R) (2006). DOI: 10.1103/PhysRevE.73.025104. https://doi.org/10.1103/PhysRevE.73.025104
  4. R. J. Baxter, 'Hard hexagons: exact solution', J. Phys. A: Math. Gen. 13, L61-L70 (1980). DOI: 10.1088/0305-4470/13/3/007. https://doi.org/10.1088/0305-4470/13/3/007
  5. G. S. Joyce, 'On the hard-hexagon model and the theory of modular functions', Philos. Trans. R. Soc. Lond. A 325, 643-702 (1988). DOI: 10.1098/rsta.1988.0077. https://doi.org/10.1098/rsta.1988.0077
  6. H. C. M. Fernandes, J. J. Arenzon, and Y. Levin, 'Monte Carlo simulations of two-dimensional hard core lattice gases', J. Chem. Phys. 126, 114508 (2007). DOI: 10.1063/1.2539141. https://doi.org/10.1063/1.2539141
  7. W. Zhang and Y. Deng, 'Monte Carlo study of the triangular lattice gas with first- and second-neighbor exclusions', Phys. Rev. E 78, 031103 (2008). DOI: 10.1103/PhysRevE.78.031103. https://doi.org/10.1103/PhysRevE.78.031103
  8. T. Nath and R. Rajesh, 'Multiple phase transitions in extended hard-core lattice gas models in two dimensions', Phys. Rev. E 90, 012120 (2014). DOI: 10.1103/PhysRevE.90.012120; arXiv:1404.6902. https://doi.org/10.1103/PhysRevE.90.012120
  9. A. A. A. Jaleel, D. Mandal, J. E. Thomas, and R. Rajesh, 'Freezing phase transition in hard-core lattice gases on the triangular lattice with exclusion up to seventh next-nearest neighbor', Phys. Rev. E 106, 044136 (2022). DOI: 10.1103/PhysRevE.106.044136; arXiv:2206.04985. https://doi.org/10.1103/PhysRevE.106.044136
  10. Z. Rotman and E. Eisenberg, 'Ideal glass transition in a simple two-dimensional lattice model', Phys. Rev. E 80, 060104(R) (2009). DOI: 10.1103/PhysRevE.80.060104; arXiv:0908.2595. https://doi.org/10.1103/PhysRevE.80.060104
  11. Z. Rotman and E. Eisenberg, 'Direct measurements of the dynamical correlation length indicate its divergence at an athermal glass transition', Phys. Rev. Lett. 105, 225503 (2010). DOI: 10.1103/PhysRevLett.105.225503. https://doi.org/10.1103/PhysRevLett.105.225503
  12. S. Darjani, J. Koplik, S. Banerjee, and V. Pauchard, 'Liquid-hexatic-solid phase transition of a hard-core lattice gas with third neighbor exclusion', J. Chem. Phys. 151, 104702 (2019). DOI: 10.1063/1.5123231; arXiv:1908.05555. https://doi.org/10.1063/1.5123231
  13. S. Darjani, J. Koplik, V. Pauchard, and S. Banerjee, 'Glassy dynamics and equilibrium state on the honeycomb lattice: Role of surface diffusion and desorption on surface crowding', Phys. Rev. E 103, 022801 (2021). DOI: 10.1103/PhysRevE.103.022801. https://doi.org/10.1103/PhysRevE.103.022801