Pith. sign in

REVIEW 1 major objections 4 minor 33 references

Rigorous bound on the aspect ratio for the formation of a nematic phase in hard rod and hard rectangle systems on $\mathbb{Z}^2$

T0 review · 1 major / 4 minor · reviewed 2026-07-14 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Hard rectangles on the square lattice form a nematic phase once their aspect ratio reaches 10^72, the first rigorous quantitative threshold.

desk verdict First explicit (if absurdly large) aspect-ratio threshold for lattice hard-rectangle nematic order; solid quantitative completion of the DG/DGJ two-scale expansion. read the letter →

arxiv 2607.10510 v1 pith:NDQSEKNJ submitted 2026-07-12 math-ph cond-mat.stat-mechmath.MPmath.PR

Qidong He This is my paper · ORCID
classification math-phcond-mat.stat-mechmath.MPmath.PR MSC 82B2082B2682B05
keywords hardrectanglesnematicphaseaspectratioboundclusterexpansiontwo-scalecoarse-graininglatticeliquidcrystalsPirogov-Sinaiorientationalorder
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper proves that hard rectangles of length l and width w on the square lattice, allowed only two orientations, spontaneously break orientational symmetry at intermediate densities when the aspect ratio k = l/w is large enough. The proof converts earlier existential cluster-expansion arguments into an explicit list of inequalities on the parameters and then checks that list for a concrete (non-optimized) choice of constants, obtaining the sufficient bound k ≥ 10^72. The resulting nematic phase is characterized by one-point densities that prefer one orientation over the other while the truncated two-point functions still decay, so translational order is absent. Although the numerical value is far larger than the simulation threshold of 7, the work supplies the first rigorous aspect-ratio estimate for lattice hard-rectangle systems and makes the bookkeeping of the two-scale expansion fully trackable.

What carries the argument

A quantitative two-scale cluster expansion: space is partitioned into tiles of side ≈ l/2, each tile is assigned a spin according to the orientation of rectangles it contains, and the resulting contour model is controlled by Peierls estimates whose constants are tracked all the way to an explicit list of twenty convergence and closure inequalities.

What would settle it

Either recompute the optimization program of Appendix A with tighter free-parameter choices and obtain a feasible aspect ratio strictly smaller than 10^72, or exhibit a concrete (w,l) with k < 10^72 for which one of the twenty inequalities is violated under every admissible parameter tuple.

Watch

Extended reading notes

Core claim

For every aspect ratio k ≥ 10^72 there exists a nonempty open interval of fugacities, written explicitly with the Lambert W function, on which the one-point densities of l imes w hard rectangles on Z^{2} satisfy liminf ho^q ≥ (3/4 - e^ϵ)z for the preferred orientation and limsup ho^{-q} ≤ (1/4 + e^ϵ)z for the opposite orientation, while the truncated two-point function decays exponentially; this establishes orientational symmetry breaking without translational order, uniformly in domain and boundary conditions.

Load-bearing premise

The whole numerical bound rests on a deliberately crude hand-chosen set of free parameters that must make every one of the twenty listed inequalities hold simultaneously; if those particular numbers fail for some admissible rectangle size, the claimed threshold of 10^72 is not justified by the written estimates.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

1 major / 4 minor

Summary. The paper proves that hard l×w rectangles on Z^{2} with two orientations exhibit a nematic phase when the aspect ratio k=l/w is at least 10^{72}. Using a quantitative two-scale cluster expansion (coarse-graining into tiles of side ≈l/2, contour representation of orientational profiles, conversion to a polymer model with multi-body interactions, truncated weights, and inductive Peierls estimates), the authors obtain explicit one- and two-point correlation bounds (Theorem 2) that establish orientational symmetry breaking without translational order on a nonempty intermediate-fugacity interval written with the Lambert W function. The main technical contribution is the complete tracking of all constants and the conversion of the argument into a finite list of closed inequalities (Appendix A.3) whose non-optimized evaluation yields the stated threshold.

Significance. This is the first rigorous, quantitative lower bound on the aspect ratio needed for a nematic phase in lattice hard-rectangle systems. Prior mathematical results (Disertori–Giuliani, Disertori–Giuliani–Jauslin) were purely existential; the present work closes that gap by producing an explicit (albeit enormous) sufficient condition and a fully checkable list of inequalities. The architecture is standard but carefully executed, and the explicit bookkeeping is a genuine advance that makes the method usable for future optimization or computer-assisted improvement. Even though 10^{72} far exceeds the numerical prediction k_min=7, the result supplies a concrete, falsifiable benchmark and a transparent optimization program.

major comments (1)
  1. The central claim of Theorem 2 rests on the simultaneous satisfaction of the full list C1–C20 under the deliberately crude parameter choices of Proposition 25 (§5.1). While the paper correctly labels the evaluation non-optimized and the arithmetic is written so that each step can be verified by hand or machine, a short independent numerical check (or a machine-readable verification script) of the key intermediate bounds (e.g., (207)–(211), (213), (217), (221), (223)) would substantially strengthen confidence that the claimed threshold is indeed justified by the written estimates. Failure of those particular numbers would only force a larger finite threshold, not invalidate the architecture, but the load-bearing numerical step should be made as transparent as possible.
minor comments (4)
  1. The enormous gap between 10^{72} and the numerical value 7 is acknowledged, but a brief remark on which free parameters (the heta_i, the scale factor 1000 in K=zℓ^{2}/1000, K_0=150, etc.) are the most wasteful would help future optimization efforts.
  2. Notation for the many derived rates (κ_i, λ_i, η_i, φ_i) is dense; a short “cheat sheet” table early in §4 or in Appendix A would improve readability.
  3. The open problems in §1.3 are well-posed; a one-sentence pointer to which of them might be approachable with the present polymer technology would be useful.
  4. Typographical consistency: the manuscript occasionally switches between “smoothing squares” and “6L-tiles”; a uniform term would help.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: self-contained quantitative completion of two-scale cluster expansion, with free parameters chosen only to close inequalities.

full rationale

The paper derives a sufficient condition (Theorem 2) for nematic order from the hard-rectangle partition function via an explicit two-scale coarse-graining, contour/polymer representation, and cluster expansion. Prior works [4,5] supply the methodological template (coarse-graining into tiles/smoothing squares, truncated weights, Peierls estimates), but the author tracks every constant, completes implicit steps, and produces a finite list of closed inequalities C1–C20 (Appendix A.3). Proposition 25 then selects crude numerical values (θ_i, K_0=150, τ=K/10, α=K/1000, ν=0, etc.) solely so that those inequalities hold simultaneously; the resulting arithmetic in §5.1 yields the explicit threshold k≳10^72. No parameter is fitted to external data that reappears as a “prediction,” no uniqueness theorem is imported from the author’s own prior work, and no ansatz is smuggled in by self-citation. The bound is acknowledged to be non-optimized and far from the numerical k_min=7; failure of the particular numbers would only enlarge the sufficient k, not collapse the existence claim. The derivation is therefore independent of its inputs by construction.

Assumptions & free parameters 5 free parameters · 5 assumptions · 1 invented entities

The central claim rests on the hard-rectangle model definition, Ueltschi’s abstract polymer cluster expansion, standard lattice geometry, and a long list of hand-chosen numerical parameters that close the Peierls and cluster-expansion inequalities. No new physical entities are postulated; the contours and polymers are standard constructions. The free parameters are the main source of the enormous gap between the rigorous bound and numerics.

free parameters (5)
  • K_0 (threshold for K=zℓ²/1000) = 150
    Set to 150 by hand in Prop. 25 so that K>K_0 forces the Lambert-W interval nonempty and all exponential smallness estimates fire.
  • scale factor 1000 in K=zℓ²/1000 = 1000
    Arbitrary large constant chosen to absorb combinatorial prefactors; directly inflates the final aspect-ratio bound.
  • θ_1=θ_2=θ_3=3/50, θ_4=θ_5=θ_6=1/2 = 3/50 and 1/2
    Free geometric splitting parameters in Lemmas 10 and 12, fixed once for all applications rather than optimized per inequality.
  • τ=K/10, α=K/1000, ϕ_1=z/4, ϕ_2=z², λ_1=K/3000 = as listed in (202)
    Peierls and cluster-expansion auxiliary rates chosen in Prop. 25 to satisfy C8–C20 simultaneously; not derived from a uniqueness principle.
  • ν=0 (with n arbitrary) = 0
    Fugacity-variation parameter set to the boundary of the allowed range; justified only by continuity for small positive ν.
assumptions (5)
  • standard math Ueltschi’s cluster-expansion theorem for abstract polymers with hard-core interactions (Thm. 3 / [29])
    Used throughout §§2.2, 3–4 to expand both the q-oriented partition function and the polymer model; convergence requires the Kotecký–Preiss-type condition (13).
  • domain assumption Hard-core exclusion geometry of l×w rectangles on Z² (Lemma 1 and the overlap conditions (22),(27)–(28))
    Determines the mesoscopic scale ℓ=⌈l/2⌉ and forces same-tile rectangles to share orientation; load-bearing for the spin assignment.
  • domain assumption Two allowed orientations only; particles are axis-parallel hard rectangles of integer sides l≥max{3,2w−1}
    Model definition in §2.1; continuous rotational symmetry is explicitly left open.
  • standard math Inductive recovery of true weights from truncated contour weights under the Peierls bound (Lemma 9 / Presutti [26])
    Breaks the circularity between polymer weights and partition-function ratios; standard in Pirogov–Sinai implementations.
  • ad hoc to paper Simultaneous satisfaction of the full list C1–C20 of convergence and closure conditions under the chosen parameters
    The numerical claim k≥10^72 is exactly the statement that these inequalities hold for the hand-chosen constants of Prop. 25; verified only by direct (crude) estimation in §5.1.
invented entities (1)
  • Two-scale contours and polymers on the 6L smoothing lattice
    purpose: Convert the hard-rectangle model into an effective polymer gas to which cluster expansion and Peierls estimates apply.
    Standard construction adapted from [4,5]; not a new physical object, but the specific 6ℓ smoothing scale and the truncated-weight induction are paper-specific bookkeeping devices.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Rigorous bound on the aspect ratio for the formation of a nematic phase in hard rod and hard rectangle systems on $\mathbb{Z}^2$." pith.science (2026). https://pith.science/paper/NDQSEKNJ

@misc{pith2026260710510,
  author       = {Pith},
  title        = {Pith review of: Rigorous bound on the aspect ratio for the formation of a nematic phase in hard rod and hard rectangle systems on $\mathbbZ^2$},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/NDQSEKNJ}},
  note         = {Machine review of arXiv:2607.10510}
}
abstract

We prove the existence of a nematic phase in a model of $l\times w$ hard rectangles on the square lattice with two allowed orientations and a large aspect ratio $k:=l/w$. The proof is based on a two-scale cluster expansion method developed previously by Disertori--Giuliani for hard rods in 2D and Disertori--Giuliani--Jauslin for hard plates in 3D. Our main contributions lie in explicitly tracking the constants and parameters and completing the arguments left implicit in these works. Hence, the proof produces a sufficient set of quantitative conditions from which estimates for the required aspect ratio can be extracted. A non-optimized evaluation of these conditions yields the bound $k\ge 10^{72}$. Although it vastly overshoots the numerical prediction, $k_{\min}=7$, our result appears to be the first rigorous estimate of the aspect ratio required for the formation of a nematic phase in hard rectangle systems.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

33 extracted references · 2 linked inside Pith

  1. [1]

    A cluster expansion approach to the Heilmann–Lieb liquid crystal model

    Diego Alberici. A cluster expansion approach to the Heilmann–Lieb liquid crystal model. Journal of Statistical Physics, 162(3):761–791, 2016

    2016

  2. [2]

    Jean Bricmont, Koji Kuroda, and Joel L Lebowitz. The structure of Gibbs states and phase co- existence for non-symmetric continuum Widom Rowlinson models.Zeitschrift f¨ ur Wahrschein- lichkeitstheorie und Verwandte Gebiete, 67(2):121–138, 1984

    1984

  3. [3]

    Entropy of fully packed hard rigid rods ond-dimensional hyper- cubic lattices.Physical Review E, 103(4):042130, 2021

    Deepak Dhar and R Rajesh. Entropy of fully packed hard rigid rods ond-dimensional hyper- cubic lattices.Physical Review E, 103(4):042130, 2021

    2021

  4. [4]

    The nematic phase of a system of long hard rods.Communications in Mathematical Physics, 323(1):143–175, 2013

    Margherita Disertori and Alessandro Giuliani. The nematic phase of a system of long hard rods.Communications in Mathematical Physics, 323(1):143–175, 2013

    2013

  5. [5]

    Plate-nematic phase in three dimensions.Communications in Mathematical Physics, 373(1):327–356, 2020

    Margherita Disertori, Alessandro Giuliani, and Ian Jauslin. Plate-nematic phase in three dimensions.Communications in Mathematical Physics, 373(1):327–356, 2020

    2020

  6. [6]

    Onsager’s spherocylinders revisited.Journal of Physical Chemistry, 91(19):4912–4916, 1987

    Daan Frenkel. Onsager’s spherocylinders revisited.Journal of Physical Chemistry, 91(19):4912–4916, 1987

    1987

  7. [7]

    Cambridge University Press, 2017

    Sacha Friedli and Yvan Velenik.Statistical mechanics of lattice systems: a concrete mathe- matical introduction. Cambridge University Press, 2017

    2017

  8. [8]

    On the orientational ordering of long rods on a lattice.Europhysics Letters, 78(2):20003, 2007

    Anandamohan Ghosh and Deepak Dhar. On the orientational ordering of long rods on a lattice.Europhysics Letters, 78(2):20003, 2007

    2007

Show all 33 references
  1. [9]

    Columnar order in random packings of 2×2 squares on the square lattice.Journal of the European Mathematical Society, 2025

    Daniel Hadas and Ron Peled. Columnar order in random packings of 2×2 squares on the square lattice.Journal of the European Mathematical Society, 2025

    2025

  2. [10]

    Extended regime of nematic order in an interacting monomer-dimer model of Heilmann and Lieb.arXiv preprint arXiv:2511.23437, 2025

    Qidong He. Extended regime of nematic order in an interacting monomer-dimer model of Heilmann and Lieb.arXiv preprint arXiv:2511.23437, 2025

    2025

  3. [11]

    Theory of monomer-dimer systems.Communications in Mathematical Physics, 25(3):190–232, 1972

    Ole J Heilmann and Elliott H Lieb. Theory of monomer-dimer systems.Communications in Mathematical Physics, 25(3):190–232, 1972

    1972

  4. [12]

    Lattice models for liquid crystals.Journal of Statistical Physics, 20(6):679–693, 1979

    Ole J Heilmann and Elliott H Lieb. Lattice models for liquid crystals.Journal of Statistical Physics, 20(6):679–693, 1979

    1979

  5. [13]

    Nematic liquid crystal phase in a system of interacting dimers and monomers.Communications in Mathematical Physics, 363(3):955–1002, 2018

    Ian Jauslin and Elliott H Lieb. Nematic liquid crystal phase in a system of interacting dimers and monomers.Communications in Mathematical Physics, 363(3):955–1002, 2018

    2018

  6. [14]

    Phase transitions in a system of hard rectangles on the square lattice.Physical Review E, 89(5):052124, 2014

    Joyjit Kundu and R Rajesh. Phase transitions in a system of hard rectangles on the square lattice.Physical Review E, 89(5):052124, 2014

    2014

  7. [15]

    Asymptotic behavior of the isotropic-nematic and nematic- columnar phase boundaries for the system of hard rectangles on a square lattice.Physical Review E, 91(1):012105, 2015

    Joyjit Kundu and R Rajesh. Asymptotic behavior of the isotropic-nematic and nematic- columnar phase boundaries for the system of hard rectangles on a square lattice.Physical Review E, 91(1):012105, 2015

    2015

  8. [16]

    Phase transitions in systems of hard rectangles with non-integer aspect ratio.The European Physical Journal B, 88(5):133, 2015

    Joyjit Kundu and R Rajesh. Phase transitions in systems of hard rectangles with non-integer aspect ratio.The European Physical Journal B, 88(5):133, 2015. 41

    2015

  9. [17]

    Joyjit Kundu, R Rajesh, Deepak Dhar, and J¨ urgen F Stilck. Nematic-disordered phase transi- tion in systems of long rigid rods on two-dimensional lattices.Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 87(3):032103, 2013

    2013

  10. [18]

    Critical behavior of long straight rigid rods on two-dimensional lattices: Theory and Monte Carlo simulations.The Journal of Chemical Physics, 128(21), 2008

    DA Matoz-Fernandez, DH Linares, and AJ Ramirez-Pastor. Critical behavior of long straight rigid rods on two-dimensional lattices: Theory and Monte Carlo simulations.The Journal of Chemical Physics, 128(21), 2008

    2008

  11. [19]

    High-density hard-core model on triangular and hexagonal lattices.Annales Henri Poincar´ e, 26(9):3321–3381, 2025

    A Mazel, I Stuhl, and Y Suhov. High-density hard-core model on triangular and hexagonal lattices.Annales Henri Poincar´ e, 26(9):3321–3381, 2025

    2025

  12. [20]

    Stability of columnar order in assemblies of hard rectangles or squares.Europhysics Letters, 114(1):10003, 2016

    Trisha Nath, Deepak Dhar, and R Rajesh. Stability of columnar order in assemblies of hard rectangles or squares.Europhysics Letters, 114(1):10003, 2016

    2016

  13. [21]

    High-activity expansion for the columnar phase of the hard rectangle gas.Journal of Statistical Physics, 160(5):1173–1197, 2015

    Trisha Nath, Joyjit Kundu, and R Rajesh. High-activity expansion for the columnar phase of the hard rectangle gas.Journal of Statistical Physics, 160(5):1173–1197, 2015

    2015

  14. [22]

    The effects of shape on the interaction of colloidal particles.Annals of the New York Academy of Sciences, 51(4):627–659, 1949

    Lars Onsager. The effects of shape on the interaction of colloidal particles.Annals of the New York Academy of Sciences, 51(4):627–659, 1949

    1949

  15. [23]

    A condition for long-range order in discrete spin systems with application to the antiferromagnetic Potts model.arXiv preprint arXiv:1712.03699, 2017

    Ron Peled and Yinon Spinka. A condition for long-range order in discrete spin systems with application to the antiferromagnetic Potts model.arXiv preprint arXiv:1712.03699, 2017

    2017 arXiv

  16. [24]

    Long-range order in discrete spin systems.arXiv preprint arXiv:2010.03177, 2020

    Ron Peled and Yinon Spinka. Long-range order in discrete spin systems.arXiv preprint arXiv:2010.03177, 2020

    2010 arXiv

  17. [25]

    Rigidity of proper colorings ofZ d.Inventiones mathematicae, 232(1):79–162, 2023

    Ron Peled and Yinon Spinka. Rigidity of proper colorings ofZ d.Inventiones mathematicae, 232(1):79–162, 2023

    2023

  18. [26]

    Springer, 2009

    Errico Presutti.Scaling limits in statistical mechanics and microstructures in continuum me- chanics. Springer, 2009

    2009

  19. [27]

    Phase transition from nematic to high-density disordered phase in a system of hard rods on a lattice.Physical Review E, 105(3):034103, 2022

    Aagam Shah, Deepak Dhar, and R Rajesh. Phase transition from nematic to high-density disordered phase in a system of hard rods on a lattice.Physical Review E, 105(3):034103, 2022

    2022

  20. [28]

    Boundary-connectivity via graph theory.Proceedings of the American Mathe- matical Society, 141(2):475–480, 2013

    ´Ad´ am Tim´ ar. Boundary-connectivity via graph theory.Proceedings of the American Mathe- matical Society, 141(2):475–480, 2013

    2013

  21. [29]

    Cluster expansions and correlation functions.Moscow Mathematical Journal, 4(2):511–522, 2004

    Daniel Ueltschi. Cluster expansions and correlation functions.Moscow Mathematical Journal, 4(2):511–522, 2004

    2004

  22. [30]

    Generalized Onsager theory of liquid crystals.Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 88(6):062501, 2013

    Xiaobin Xiao and Ping Sheng. Generalized Onsager theory of liquid crystals.Physical Review E—Statistical, Nonlinear, and Soft Matter Physics, 88(6):062501, 2013

    2013

  23. [31]

    An alternate version of Pirogov-Sinai theory.Communications in Mathe- matical Physics, 93(4):559–581, 1984

    Miloˇ s Zahradn´ ık. An alternate version of Pirogov-Sinai theory.Communications in Mathe- matical Physics, 93(4):559–581, 1984. A Constants, parameters, and conditions In this appendix, we collect the constants, parameters, and convergence and closure conditions appearing in ...

    1984

  24. [32]

    10(ii) θ3 (0, 1

    Lem. 10(ii) θ3 (0, 1

  25. [33]

    10(iii) θ4 (0,1) Lem

    Lem. 10(iii) θ4 (0,1) Lem. 12(v) θ5 (0,1) Lem. 12(vi) θ6 (0,1) Lem. 12(vii) α(0,∞) Lem. 24 A.2 Derived quantities Recall the shorthand notationℓ=⌈l/2⌉andϵ=e 1+ν(2w−1)(2l−1)z. The constants and derived parameters are as follows. Expression Location Combinatorial constants cn(ε)...

Pith tools

Reviewed July 14, 2026 · model on record in the stance chip above.