Pith. sign in

REVIEW 3 major objections 4 minor 17 references

New results on the domain of analyticity of the free energy for the Ising model

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

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

Pith's one-line read The Ising pressure is analytic for β≤0.322 in two dimensions, the paper claims, using a cluster expansion built on cycle decompositions.

desk verdict The high-temperature theorem fails on two independent grounds—the polymer gas (2.32) is not the high-temperature expansion, and Lemma 3.2 contradicts the stated interaction—but the strong-field and low-temperature sections have salvageable content. read the letter →

arxiv 2608.08396 v1 pith:3RG6XUNZ submitted 2026-08-09 math-ph math.COmath.MP

Nguyen Tong Xuan , Nguyen Dang Thien Thu This is my paper · ORCID
classification math-phmath.COmath.MP MSC 82B2082B05
keywords IsingmodelfreeenergyanalyticityclusterexpansionFernandez-ProcaccicriterionVeblen'stheoremgeneratingfunctionsPeierlscontoursself-avoidingpolygons
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

The paper claims that for the Ising model on Z^d with an external field, the free-energy function is analytic on strictly larger parameter regions than previously known: a strong-field region, a high-temperature region at zero field, and a low-temperature region. The mechanism is a cluster expansion for polymer and contour representations, with convergence controlled by the Fernandez–Procacci criterion. The headline claim is the high-temperature region: in dimension two the pressure p^∅(β) is analytic for β≤0.322, well beyond the classical published bound and the value 0.151 previously reported for d=2. The authors obtain the improvement by rewriting the high-temperature expansion using Veblen's theorem, which decomposes even subgraphs into edge-disjoint cycles, and by counting polygons through known self-avoiding-polygon bounds.

What carries the argument

The load-bearing object is the polymer and cluster expansion for the Ising partition function, converted into a convergence problem by the Fernandez–Procacci criterion: the series converges if a Gruber–Kunz condition of the form sup_x ∑_{S∋x} w(S)$e^{{a|S|}}$ ≤ e^a−1 holds. Around this, the paper uses three counting devices: a generating function p(X)=X(1+p(X))^{2d} for connected subsets in the strong-field regime; Veblen's theorem in the high-temperature regime, which turns the even-subgraph expansion into a gas of primitive cycles (polygons) with activity (tanh β)^{|p|}; and a floor-stack generating function h($X^{{1/d}}$Y)≤g(X,Y)+(d−1)\tilde g(X,Y) for primitive contours at low temperature. The sharpened analyticity regions are obtained by optimizing the auxiliary parameter a in each Gruber–Kunz condition.

What would settle it

Look for a triangle in the lattice graph defined by ||i−j||_∞=1: in d=2 the three sites (0,0), (1,0), and (0,1) are pairwise adjacent, so the graph contains an odd cycle. Since Lemma 3.2 asserts every cycle—and hence every primitive polygon—has an even number of edges, this single configuration disproves the lemma and removes the justification for summing only even powers in the high-temperature convergence condition.

Watch

Extended reading notes

Core claim

On the paper's own terms, the discovery is an improved analyticity theorem for the Ising pressure in all three regimes. In a strong external field, analyticity holds whenever 2h≥φ_st(d) with φ_st(d)=(4d+1)log(1+1/(4d))+log(4d), which is asymptotically much smaller than the standard textbook threshold. At zero field and high temperature, the pressure p^∅(β) is analytic for β≤β_N(d)=$tanh^{{-1}}$[$φ_1^{{high}}$(\bar a_{high})], computed numerically as 0.322 in d=2; at low temperature, analyticity holds for β above a threshold φ_low(m,κ) whose d=2 value is 0.8226. The proof combines the Fernandez–Procacci convergence criterion with Veblen's theorem (every finite connected graph all of whose vertices have even degree is an edge-disjoint union of simple cycles) and with generating-function bounds on the number of connected subsets, contours, and polygons.

Load-bearing premise

The argument's load-bearing step is the assertion that every primitive cycle in the lattice has even length, which is needed to reduce the convergence sum to even powers but fails for the sup-norm adjacency rule already in d=2.

Editorial extensions

If this is right

  • In dimension two, the high-temperature analyticity threshold becomes β≤0.322 instead of β≤0.151, so the technique roughly doubles the domain.
  • The ratio of the new strong-field threshold to the standard one decays exponentially in d, so the gain grows with dimension in the magnetic-field regime.
  • The low-temperature threshold in d=2, β≥0.8226, is smaller than the previously known 0.94, so analyticity extends closer to the region where order appears.
  • All three regimes are treated through one framework, so the same Fernandez–Procacci plus generating-function strategy can be applied to other lattice spin systems with polymer representations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the high-temperature argument survives a correction for odd cycles, then counts of self-avoiding polygons with even and odd length would both enter the convergence sum; because the number of odd polygons in the sup-norm graph grows similarly (2d−1)^{k−1}, the final threshold would likely be smaller than 0.322 but still above the classical bound.
  • The generating-function bounds appear to transfer to any lattice with bounded coordination number; one could test the strong-field and low-temperature thresholds on the triangular lattice, where the bipartite assumption already fails.
  • The d=2 high-temperature value 0.322 remains below the known critical inverse temperature of the square-lattice Ising model, so the method does not locate the phase transition; a natural next step is to check whether the same cluster expansion can be pushed closer to the critical point with anisotropic weights.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper studies the domain of analyticity of the free energy of the Ising model in three regimes: strong external field, high temperature, and low temperature. The authors use the Fernandez--Procacci cluster-expansion criterion, generating-function bounds for the number of polymers/contours, and a graph-theoretic high-temperature expansion based on Veblen's theorem. The main advertised results are an improved strong-field threshold 2h ≥ φ_st(d), a high-temperature analyticity domain β ≤ β_N(d) that is claimed to exceed the classical Simon bound and, in d=2, to give β ≤ 0.322 versus Procacci's 0.151, and a low-temperature domain β ≥ min φ_low(m,κ) that is claimed to improve on Balister--Bollobás and Procacci. The proof strategy is clearly organized, and the Fernandez--Procacci criterion is stated in the appendix in a standard form. However, the central high-temperature identity is invalid, and the strong-field counting does not match the model defined by the stated interaction.

Significance. If the advertised results were correct, they would constitute a substantial quantitative improvement over existing analyticity bounds for the Ising model, and the combination of generating functions with the Fernandez--Procacci criterion would be a useful contribution. The paper also makes a genuine effort to compare constants and to present the cluster-expansion framework self-contained. Unfortunately, the high-temperature theorem, which is explicitly described as one of the main contributions, rests on a false equality between the even-subgraph expansion and a polymer gas over primitive polygons. The strong-field theorem is also not established for the model as defined, because the counting argument uses the wrong graph degree. These are load-bearing errors, not presentation issues, and they invalidate the headline quantitative claims.

major comments (3)
  1. [§2.3, Eq. (2.32); §3.2, proof of Theorem 2.4] The passage from the even-subgraph expansion (2.29) to the primitive-polygon gas (2.32) is invalid. A finite even subgraph need not be a disjoint union of pairwise compatible primitive polygons: on Z^2 with nearest-neighbor edges, the union of two unit squares sharing exactly one vertex is an even subgraph with eight edges, but under Definition 2.2 it decomposes into two primitive square polygons that share a site. Since d(p1,p2)=0, the compatibility factor ζ(p1,p2) in (2.32) vanishes and the tuple is forbidden, whereas the original sum (2.29) contains this graph with weight (tanhβ)^8. Thus (2.32) is not equal to the high-temperature partition function. This failure is independent of the bipartiteness issue and directly invalidates Theorem 2.5 and the claimed value β_N(2)=0.322.
  2. [Lemma 3.2 and Eq. (3.37), §3.2] Lemma 3.2 is false for the interaction defined in (2.2). With ||i−j||∞=1, the graph contains triangles for every d≥2 (for example, the vertices (0,0), (1,0), and (0,1) in d=2 are pairwise at L∞ distance 1), so not every cycle has even length. The proof's statement that 'Z^d is hypercube graph' is also not correct for this adjacency. Consequently, the even-power restriction in the convergence condition (3.37) is unjustified, and the analyticity region in Theorem 2.5 is not established even if the representation problem identified in the previous comment were resolved.
  3. [§3.1, Eq. (3.21); Theorem 2.2] The strong-field counting uses the recurrence p_{L+1}(X)=X(1+p_L(X))^{2d}, which is the correct branching bound for the standard nearest-neighbor lattice with degree 2d. Under the model defined by (2.2), however, adjacency is d(i,j)=||i−j||∞≤1, so each site has 3^d−1 neighbors. The recurrence therefore does not bound the number of connected sets for the model actually under study, and the threshold φ_st(d) in Theorem 2.2 is not justified. In addition, the comparison with Friedli--Velenik is not a comparison with the same model, since Friedli--Velenik treat the standard nearest-neighbor Ising model.
minor comments (4)
  1. [Eqs. (2.34) and (2.38)] The pressure formula in (2.34) should contain log coshβ, not β, and the thermodynamic-limit formula in (2.38) should contain d log coshβ rather than d coshβ; the displayed expressions are missing the logarithm.
  2. [Throughout] The notation for the lattice edges is inconsistent: (2.2) and (2.11) use ||i−j||∞=1, while (2.27) uses |x−y|=1. Since the two choices define different models, the authors should state one convention and use it consistently.
  3. [§2.4, Eq. (2.52)] The partition function Ξ^LT_Λ is called the 'large-field polymer partition function' in the sentence before (2.52); this should be 'low-temperature polymer partition function'.
  4. [Eq. (2.43)] The numerical value 1.0783 is asserted as the result of a maximization over a, but the maximizing argument is not displayed; the authors should provide the calculation or a reference for this value.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the derivation relies on external convergence criteria and combinatorial bounds, with only a minor non-load-bearing self-citation.

full rationale

No circular step reduces a target result to an input. The claimed analyticity domains are sufficient conditions obtained by applying the Fernández–Procacci convergence criterion (Proposition A.1, from ref. [5]) to polymer and contour gases, then bounding polygon/contour counts by generating-function arguments from refs. [1], [8], and [10], and finally optimizing the free auxiliary parameter a>0. The comparisons with Simon and Balister–Bollobás are posterior consistency checks, not inputs to the derivation. The only self-citation, ref. [12], is cited alongside refs. [2,5] as a general reference for cluster-expansion theory and is not used as a load-bearing theorem; no uniqueness claim from the authors' prior work is invoked. The score of 2 reflects only the presence of that minor, non-load-bearing self-citation, not actual circularity. The paper's serious weaknesses are correctness issues rather than circularity: the equality between the high-temperature even-subgraph expansion (2.29) and the polygon polymer gas (2.32) is not established because edge-disjoint Veblen cycles may share vertices, which the compatibility factor (2.13) would forbid, and Lemma 3.2's bipartite assumption conflicts with the L-infinity interaction (2.2). These concerns do not make the derivation equivalent to its own assumptions, so they do not raise the circularity score.

Assumptions & free parameters 3 free parameters · 6 assumptions · 0 invented entities

The central claims depend on the Fernandez-Procacci criterion and combinatorial counting bounds. No physical constants are fitted to experimental data; the auxiliary parameters a, m, kappa are optimization variables in the proof. The notable load-bearing assumptions are the bipartiteness of the lattice, which is false under (2.2), and the uniqueness of primitive-polygon decompositions, which is assumed without proof.

free parameters (3)
  • a (strong-field auxiliary) = a_st = log(1+1/(4d))
    Auxiliary parameter in the Fernandez-Procacci test function; optimized to minimize phi_st(a). Not a physical constant and not fitted to data.
  • a_high (high-temperature auxiliary) = unique maximizer a_bar_high of phi_1^high(a), no closed form given
    Scaling parameter in the Gruber-Kunz condition for the polygon gas; chosen to maximize the analyticity domain.
  • (m, kappa) (low-temperature auxiliary pair) = not explicitly optimized in general; for d=2 the resulting minimum is 0.822614
    Parameters in the floor-stack contour bound; constrained to set L and minimized to yield the low-temperature threshold.
assumptions (6)
  • standard math Fernandez-Procacci criterion for convergence of cluster expansions (Proposition A.1)
    Quoted from [5] and used as the core convergence tool in all three regimes.
  • standard math Veblen's theorem: a finite graph is an edge-disjoint union of cycles iff every vertex has even degree (Lemma 2.3)
    Used to decompose even subgraphs in the high-temperature expansion.
  • domain assumption Generating function p(X)=X(1+p(X))^{2d} bounds the number of connected subsets containing the origin
    Standard tree-function bound for lattice animals, derived via p_{L+1}=X(1+p_L)^{2d} in Section 3.1, using coordination number 2d.
  • ad hoc to paper The lattice graph with the paper's adjacency is bipartite, so all cycles have even length
    Invoked in Lemma 3.2 to restrict polygon sums to even lengths; contradicts the L-infinity adjacency in (2.2), which yields triangles in d>=2.
  • ad hoc to paper Each even subgraph has a unique decomposition into edge-disjoint primitive polygons
    Unstated, used implicitly in passing from (2.29) to (2.32); not proven and likely false, leading to overcounting.
  • standard math Isoperimetric bound |∂S| >= d V_d(1)^{1/d} |S|^{(d-1)/d} for connected subsets S
    Used in Lemma 3.1 to drop the beta-dependence in the strong-field convergence condition.

how reviews work

0 comments
Cite this review

Pith. "Pith review of New results on the domain of analyticity of the free energy for the Ising model." pith.science (2026). https://pith.science/paper/3RG6XUNZ

@misc{pith2026260808396,
  author       = {Pith},
  title        = {Pith review of: New results on the domain of analyticity of the free energy for the Ising model},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3RG6XUNZ}},
  note         = {Machine review of arXiv:2608.08396}
}
read the original abstract

We investigate the analyticity of the free energy of the Ising model in the presence of a non-zero external magnetic field, at high temperature, and at low temperature. Using the Fernandez--Procacci convergence criterion for cluster expansions, together with generating-function techniques and graph-theoretical methods, we derive improved convergence conditions in all three regimes. In particular, the generating-function approach yields sharper estimates for polymers and contours in the strong-field and low-temperature regimes, while a new high-temperature expansion based on Veblen's theorem provides a substantially larger analyticity region than the classical results in the literature.

Figures

Figures reproduced from arXiv: 2608.08396 by the authors.

Figure 1
Figure 1. A configuration of the Ising model. Each connected component of the shaded area delimits one of the polymers S1, . . . , S6. partition function, let us introduce the definitions of compatible and incompatible objects as follow. Definition 2.1. Let us define S, S′ to be compatible, and denote S ∼ S ′ , if d(S, S′ ) ≥ 2. Otherwise S and S ′ are incompatible and we denote S ≁ S ′ . Denote (2.13) ζ(S, S′ ) =  1 if S ∼ … view at source ↗
Figure 2
Figure 2. A comparison with the result of Friedli and Velenik (The green and red lines represent the Friedli–Velenik bound and our result, respectively.) [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. The rate between our bound and the Friedli–Velenik bound 2.3. The Ising model at high-temperature without magnetic field (h = 0). In the case of high temperatures and a vanishing external field, for simplicity in computation, we consider the Ising model with free boundary conditions. Its partition function is defined as: (2.27) Z ∅ Λ (β, 0) = X σΛ exp ( β X x,y∈Λ |x−y|=1 σxσy ) . We use the identity (2.28) exp{βσxσy… view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: A comparison with Simon’s result (green line presented for Simon’s result and red line presented for our result) [PITH_FULL_IMAGE:figures/full_fig_p010_4.png]
Figure 5
Figure 5. Figure 5: The rate between our bound and Simon’s bound As noted in [14], Procacci established an analyticity domain for p ∅(β, 0) in the case d = 2, which is valid for β ≤ 0.151. However, Procacci’s bound is considerably weaker than our result, which extends the limit to β ≤ βN …
Figure 6
Figure 6. Figure 6: A configuration of the two-dimensional Ising model in a finite box Λ with + boundary condition. At low temperature, the lines separating regions of + and − spins are expected to be short and sparse, leading to a positive magnetization in Λ. Let Vγ be the set containing…
Figure 7
Figure 7. Figure 7: Comparison with the result of Balister and Bollob´as (the green and red lines represented our result, Balister–Bollob´as result, respectively) In particular, for d = 2, Procacci [14] obtain an analyticity domain for (2.67) β ≥ 0.94. In this case, we obtain a better dom…
Figure 8
Figure 8. Figure 8: The rate between our bound and Balister and Bollob´as bound (i) Let us start with the following direct consequence of from Proposition A.1 that |Γ|S (wβ,h) converges if for each S ∈ P, (3.1) 1 +X n≥1 X (S1,...,Sn)∈Pn S≁Si , Si∼Sj , 1≤i,j≤n Yn i=1 wβ,h(Si)ea(Si) ≤ e a(S…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

17 extracted references · 17 canonical work pages

  1. [1]

    Communications in Mathematical Physics, 273 (2), 305-315 (2007)

    Balister, P., Bollob´ as, B.:Counting regions with bounded surface area. Communications in Mathematical Physics, 273 (2), 305-315 (2007)

    work page 2007

  2. [2]

    598–617 (2010)

    Bissacot R., Fern´ andez R., Procacci A.:On the convergence of cluster expansions for polymer gases.Journal of Statistical Physics 139, pp. 598–617 (2010)

    work page 2010

  3. [3]

    Brydges D. C. :A short cluster in cluster expansions.In Critical Phenomena, Random Systems, Gauge Theories, Osterwalder, K. and Stora, R. (eds.), Elsevier, 129–83 (1984)

    work page 1984

  4. [4]

    Dobrushin R. L.:Estimates of semiinvariants for the Ising model at low temperatures.Topics in Statistics and Theoretical Physics, American Mathematical Society Translations (2) 177, 59–81 (1996)

    work page 1996

  5. [5]

    New bounds from an old approach

    Fern´ andez R., Procacci A.:Cluster expansion for abstract polymer models. New bounds from an old approach. Communications in Mathematical Physics 274(1), 123-140 (2007)

    work page 2007

  6. [6]

    Cambridge University Press (2017)

    Friedli S., Velenik Y.:Statistical Mechanics of Lattice Systems: A Concrete Mathematical Introduction. Cambridge University Press (2017)

    work page 2017

  7. [7]

    Gruber C., Kunz H.:General properties of polymer systems.Communications in Mathematical Physics 22,133–61 (1971)

    work page 1971

  8. [8]

    (eds) Polygons, Polyominoes and Polycubes

    Guttmann, A.J.:History and Introduction to Polygon Models and Polyominoes.In: Guttmann, A.J. (eds) Polygons, Polyominoes and Polycubes. Lecture Notes in Physics, vol 775. Springer, Dordrecht (2009)

    work page 2009

Show all 17 references
  1. [9]

    D.:Complex zero-free regions at large —q— for multivariate Tutte poly- nomials (alias Potts-model partition functions) with general complex edge weights

    Jackson B., Procacci A., Sokal A. D.:Complex zero-free regions at large —q— for multivariate Tutte poly- nomials (alias Potts-model partition functions) with general complex edge weights. Journal of Combinatorial Theory Series B 103, 21–45 (2013)

    2013

  2. [10]

    L., Mazel A

    Lebowitz J. L., Mazel A. E.:Improved Peierls Argument for High-Dimensional Ising Models. Journal of Statistical Physics 90, Issue 3–4, 1051–1059 (1998)

    1998

  3. [11]

    Ndreca S., Procacci A., Scoppola B.:Improved bounds on coloring of graphs.European Journal Combinatorics 33, 592-609 (2012)

    2012

  4. [12]

    X.; Fern´ andez R.:Convergence of Cluster and Virial expansions for Repulsive Classical Gases

    Nguyen T. X.; Fern´ andez R.:Convergence of Cluster and Virial expansions for Repulsive Classical Gases. Journal of Statistical Physics (2020), 179, 448–484

    2020

  5. [13]

    Penrose O.:Convergence of fugacity expansions for classical systems.In Statistical mechanics: foundations and applications. A. Bak (ed.), Benjamin, New York (1967)

    1967

  6. [14]

    Procacci A.:Cluster expansion methods in rigorous statistical mechanics.arXiv:2308.06380v1

    arXiv

  7. [15]

    Veblen O.:An application of modular equations in analysis situs. Ann. Math. 14 (1912), 86–94

    1912

  8. [16]

    Simon B.:The statistical mechanics of lattice gases. Vol. I. Princeton Series in Physics. Princeton University Press, Princeton, NJ, (1993)

    1993

  9. [17]

    E., Mayer M

    Mayer J. E., Mayer M. G.:Statistical Mechanics. John Wiley and Sons Inc, New York (1940). Quy Nhon University, 170 An Duong Vuong street, Quy Nhon Nam ward, Gia Lai province, Vietnam Email address:nguyentongxuan@qnu.edu.vn Quy Nhon University, 170 An Duong Vuong street, Quy Nh...

    1940

Pith tools

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