REVIEW 3 major objections 4 minor 23 references
Partition function zeros for the Blume-Capel model on a complete graph
T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read For the Blume-Capel model on a complete graph, partition-function zeros show that finite-size systems up to N=1500 still display effective critical exponents far from the thermodynamic-limit mean-field values.
desk verdict A useful finite-size zeros map for the complete-graph Blume-Capel model, but the headline claim rests on power-law fits that need corrections-to-scaling and error bars before I'd trust the exponents. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the exact integral representation of the partition function, Eq. (4), which expresses $Z_N$ as a single integral over a collective variable $x$ for any finite $N$. On top of this, the paper uses partition-function zeros in three complex planes: Fisher zeros (complex temperature), crystal-field zeros (complex $\Delta$), and Lee-Yang zeros (complex magnetic field). These zeros are found by solving the simultaneous equations $\mathrm{Re}\,Z=0$ and $\mathrm{Im}\,Z=0$ (for Lee-Yang zeros, $\mathrm{Re}\,Z=0$ alone). The argument is carried by scaling relations of the form $T_j \sim (j/N)^{1/(2-\alpha)}$, $H_j \sim (j/N)^{\beta\delta/(2-\alpha)}$, and the impact-angle relation (12), which connect the first zero's distance from the real axis to the critical exponents. By fitting the size dependence of the first zero over $N=20$ to $1500$, the paper extracts effective gap exponents and compares them with the mean-field expectations $1/2$, $2/3$, $3/4$ and $5/6$.
What would settle it
A concrete check: compute the first Fisher zero at the tricritical point ($T_t=1/3$, $\Delta_t=2\ln 2/3$) for $N=10^4$ to $10^6$ using the exact integral representation and fit $T_1(N)$ with the asymptotic form $a N^{-2/3}$ plus a logarithmic-correction term. If the effective exponent moves toward $2/3$ once those corrections are included, the paper's claim of non-asymptotic criticality at $N=1500$ is a finite-size fitting artifact; if the exponent remains near $0.4$ even at $N=10^6$, the claim is supported.
Extended reading notes
Core claim
The central claim is that finite-size criticality in this exactly solvable mean-field model is not asymptotic for the sizes studied. The authors analyse the exact integral representation of the partition function and locate its zeros in the complex temperature, crystal-field, and magnetic-field planes. Near the tricritical point at $(T_t,\Delta_t)=(1/3,2\ln 2/3)$, the first Fisher zero scales as $T_1 \sim N^{-0.35}$ for all sizes up to $N=1500$, or $N^{-0.41}$ using only the three largest sizes, compared with the asymptotic mean-field exponent $2/3$. Near the Ising critical point $(T=2/3,\Delta=0)$ the effective exponent is about 0.44, not 0.5. Crystal-field zeros behave well only near the tricritical point and scale with effective exponents close to $2/3$ even when the thermodynamic point is critical ($1/2$). Lee-Yang zeros in the complex magnetic field obey the Lee-Yang circle theorem and their gap exponent approaches $3/4$ at the critical point and $5/6$ at the tricritical point, but along most of the critical line the fitted exponents are skewed toward the tricritical value.
Load-bearing premise
The load-bearing premise is that a pure power law $T_1 \sim N^{-g}$, fitted over the available range of $N$, captures the asymptotic exponent; if the data instead contain logarithmically suppressed corrections, the fitted effective exponents could be artifacts of omitting them.
Editorial extensions
If this is right
- Finite-size scaling analyses that use only a few hundred or thousand sites for this model will report effective exponents that are not the true mean-field values; the gap is especially large at the tricritical point.
- Crystal-field zeros are a reliable probe near the tricritical point but a poor one near the Ising critical point, because the scaling field is perpendicular to the tricritical point rather than to the critical point.
- The point $(T=1/2,\Delta=\ln 2/2)$, which lies on the pseudo-transition line $\Delta=T\ln 2$, shows unusually fast convergence to the critical exponent and may be a useful reference point in finite-size studies.
- Lee-Yang zeros confirm the Lee-Yang circle theorem on the complete graph: all zeros lie on the imaginary magnetic-field axis at the critical line and the tricritical point.
- Effective exponents extracted from the first Fisher zero depend on which system sizes are included (0.35 versus 0.41 at the tricritical point), so reported exponents should carry this dependence.
Reading between the lines
- An implication the authors leave implicit is that the fitted pure power laws may be reconciled with the true mean-field asymptotics by including logarithmic corrections, which are known to appear in mean-field models above the upper critical dimension; if so, the effective exponents would be a crossover phenomenon rather than evidence of a new universality class.
- A testable extension would be to repeat the zero analysis at substantially larger $N$ (for example $10^4$ to $10^5$) and check whether the effective Fisher-zero exponent increases monotonically toward $2/3$; the paper does not report such data.
- The observation that crystal-field zeros are skewed toward tricritical behaviour even along the critical line suggests that fits of thermodynamic quantities near tricritical points in other spin-1 models may similarly be contaminated by tricritical scaling, a caution that could be tested on related models.
- One could also fit the first Fisher zero with explicit correction terms, such as $T_1 = a N^{-2/3}(1 + b N^{-\omega})$, to see whether the effective exponent 0.35 to 0.41 at $N \le 1500$ is a leading-order artefact or a genuine slow crossover.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies finite-size effects in the Blume-Capel model on a complete graph by locating zeros of the exact integral representation of the partition function, Eq. (4), in the complex temperature, crystal-field, and magnetic-field planes for system sizes up to N=1500. The thermodynamic-limit mean-field exponents are recalled, and the finite-size scaling of the first Fisher, crystal-field, and Lee-Yang zeros is compared with those asymptotic exponents. The central claim is that even at the largest accessible sizes the effective critical behaviour read from Fisher zeros is not asymptotic, e.g. T_1 ~ N^{-0.35} or N^{-0.41} at the tricritical point instead of N^{-2/3}. The paper also reports angle estimates, crystal-field-zero scaling along the critical line, and running Lee-Yang gap exponents.
Significance. If the central claim is correct, the paper provides a useful example in which partition-function-zero scaling at accessible finite sizes is governed by effective exponents that differ substantially from the known mean-field asymptotics, which matters for interpreting finite-size data in mean-field and tricritical systems. The work has clear strengths: the starting point is an exact integral representation, the thermodynamic-limit phase diagram is derived analytically, the standard zero-scaling relations are applied, and three different complex-field planes are compared systematically. However, the central quantitative conclusion rests on pure-power-law fits without corrections to scaling, logarithmic terms, or error bars, so the distinction between 'non-asymptotic effective behaviour' and 'large corrections to scaling' is not yet established.
major comments (3)
- [§5.1.3, Eq. (4), Fig. 9] The main claim that criticality is not asymptotic at the tricritical point is based on fitting Im T_1(N) to a pure power law, giving N^{-0.35} or N^{-0.41} versus the expected N^{-2/3}. No corrections to scaling, no logarithmic factors, and no error bars are reported. At a tricritical point the quartic Landau coefficient vanishes, so logarithmic and higher-order corrections are expected; Ref. [8] is precisely about such logarithmic corrections to mean-field scaling. The authors should test fits of the form T_1 = A N^{-g}(1 + B N^{-theta}) or T_1 = A N^{-g}(ln N)^p over the available N range, and report the resulting g-values with uncertainties. If g moves toward 2/3 under such fits, the conclusion reduces to 'corrections are large' rather than 'the asymptotic exponent is not approached.'
- [§5.1.1, Fig. 6] The analogous pure-power-law analysis at the Ising critical point gives g=0.44 or 0.46 versus the expected 0.5, but the fitted range is limited (N=50 to 500, or only the three largest sizes) and no error bars or corrections-to-scaling fits are shown. The statement that the scaling 'approaches the theoretical value, even though slowly' is plausible, but the numerical support is not quantitative. Please provide stability checks under different N-windows and fits with corrections, and report confidence intervals for the effective exponents.
- [§5.3.1, Figs. 18-19] The running Lee-Yang exponent g_h is computed by repeated fits that add one system size at a time, but no uncertainties are quoted, and the conclusion that the exponents 'finally converge' toward 3/4 or 5/6 is drawn from trajectories without error bars. The qualitative trend is visible, but a quantitative statement about convergence requires either error bars or a stated criterion for convergence. This is important because the conclusions section asserts that Lee-Yang scaling 'matched theoretical predictions.'
minor comments (4)
- [Introduction, Section 5] The Introduction promises an 'expanded function approximation' that 'simplifies the calculations without sacrificing accuracy,' but this method is never defined or used anywhere in the paper. The authors should either deliver it or remove the promise from the Introduction.
- [§5.2.1] The claim that crystal-field zeros at the Ising critical point are not well defined because the scaling field is perpendicular to the tricritical point is stated very briefly; a short explanation of why this makes the zeros poorly defined would improve readability.
- [Conclusions] The conclusion states that 'the equivalence between Fisher and crystal field zeros was confirmed analytically,' but the body of the paper only establishes that both obey the same expected scaling relation, Eq. (13) and Eq. (20). The wording is stronger than the presented derivation and should be adjusted.
- [Notation in §5.1.3, Fig. 9] The text writes T_1 ~ N^{0.35} in one place and T_1 ~ (1/N)^{0.35} in another. Since the quantity is the imaginary part of the first zero decreasing with N, the notation should be made uniform to avoid sign ambiguity.
Circularity Check
No significant circularity: finite-size zero scalings are measured from the exact partition function and compared with independently known mean-field exponents; self-citations supply methodology, not the result.
full rationale
The paper's central claim is that the first Fisher-zero coordinate scales as N^{-0.44} at the Ising critical point and N^{-0.35}/N^{-0.41} at the tricritical point, versus the asymptotic mean-field values N^{-1/2} and N^{-2/3}. These asymptotic exponents are not fitted parameters; they are taken from the standard mean-field solution of the Blume-Capel model and from the known scaling relation (13), which is cited to general results [5,8,10,19] rather than derived from the present fits. The finite-size exponents are outputs of power-law fits to zeros obtained by numerically solving the exact integral representation (4), so the measured exponents are not imposed by the asymptotic values being tested. The self-citations in the paper ([8,10,14,15,22]) are used for methodological background, for prior applications of partition-function-zero techniques, and for heuristic expectations about crystal-field zeros; none of them is invoked as a uniqueness theorem or as the source of the numerically fitted exponent values. The discussion of the special point T=1/2, Delta=ln2/2 is a post-hoc observation that this point lies on the pseudo-transition line, but it does not define the scaling exponents through that line. Concerns about neglecting corrections to scaling, logarithmic terms, or error bars are legitimate statistical robustness questions, but they are not circularity: the comparison between a fitted effective exponent and an independently known asymptotic exponent remains an external benchmark. No equation is defined in terms of the quantity it purports to predict, and no fitted parameter is renamed as a prediction. Therefore the derivation chain is self-contained with respect to the central finite-size claim, and no circular step can be exhibited from the text.
Assumptions & free parameters
free parameters (4)
- Effective Fisher-zero exponent near Ising point =
0.44 (all sizes), 0.46 (largest only)
- Effective Fisher-zero exponent near tricritical point =
0.35 (all sizes), 0.41 (largest only)
- Effective crystal-field-zero exponent =
0.62 (T=0.4), 0.69 (T=0.35)
- Effective Lee-Yang gap exponent =
Roughly 0.75 to 0.85 depending on point and fit window
assumptions (4)
- standard math Stratonovich-Hubbard transformation yields exact integral representation of the partition function (Eq. 4).
- domain assumption The scaling relations (12) and (13) connect partition function zero coordinates to critical exponents (known from Lee-Yang and Fisher theory).
- domain assumption Mean-field critical exponents for the Blume-Capel model (Table 3) are the asymptotic values.
- ad hoc to paper The finite-size scaling of the first zero follows a pure power law T1 ~ N^{-g} with no corrections-to-scaling over the fitted range.
Cite this review
Pith. "Pith review of Partition function zeros for the Blume-Capel model on a complete graph." pith.science (2026). https://pith.science/paper/4TTOKQW3
@misc{pith2026250104452,
author = {Pith},
title = {Pith review of: Partition function zeros for the Blume-Capel model on a complete graph},
year = {2026},
howpublished = {\url{https://pith.science/paper/4TTOKQW3}},
note = {Machine review of arXiv:2501.04452}
}
read the original abstract
In this paper we study finite-size effects in the Blume-Capel model through the analysis of the zeros of the partition function. We consider a complete graph and make use of the behaviour of the partition function zeros to elucidate the crossover from effective to asymptotic properties. While in the thermodynamic limit the exact solution yields the asymptotic mean-field behaviour, for finite system sizes an effective critical behaviour is observed. We show that even for large systems, the criticality is not asymptotic. We also present insights into how partition function zeros in different complex fields (temperature, magnetic field, crystal field) give different precision and provide us with different parts of the larger picture. This includes the differences between criticality and tricriticality as seen through the lens of Fisher, Lee-Yang, and Crystal Field zeros.
Figures
Figures from the paper (16 more)
Reference graph
Works this paper leans on
- [8]
-
[1]
M. Blume. Theory of the First-Order Magnetic Phase Change in U O 2. Phys. Rev., 141(2):517– 524, January 1966
work page 1966
-
[2]
H.W. Capel. On the possibility of first-order phase transitions in Ising systems of triplet ions with zero-field splitting. Physica, 32(5):966–988, May 1966
work page 1966
-
[3]
Ian Lawrie and Stephane Sarbach. Theory of Tricritical Points. In Cyril Domb and Joel Louis Lebowitz, editors,Phase transitions and critical phenomena, volume 9. Academic press, London, 1984
work page 1984
-
[4]
Robert B. Pearson. Partition function of the Ising model on the periodic 4 ×4×4 lattice. Phys. Rev. B, 26:6285–6290, Dec 1982. 23
work page 1982
-
[5]
C. Itzykson, R.B. Pearson, and J.B. Zuber. Distribution of zeros in Ising and gauge models. Nuclear Physics B, 220(4):415–433, September 1983
work page 1983
-
[6]
M. L. Glasser, V. Privman, and L. S. Schulman. Complex temperature plane zeros in the mean-field approximation. J Stat Phys, 45(3-4):451–457, November 1986
work page 1986
-
[7]
M. L. Glasser, V. Privman, and L. S. Schulman. Complex-temperature-plane zeros: Scaling theory and multicritical mean-field models. Phys. Rev. B, 35(4):1841–1845, February 1987
work page 1987
Show all 23 references
-
[9]
Krasnytska, B
M. Krasnytska, B. Berche, Yu. Holovatch, and R. Kenna. Violation of Lee-Yang circle theorem for Ising phase transitions on complex networks. EPL, 111(6):60009, September 2015
2015
-
[10]
Partition function zeros for the Ising model on complete graphs and on annealed scale-free networks
M Krasnytska, B Berche, Yu Holovatch, and R Kenna. Partition function zeros for the Ising model on complete graphs and on annealed scale-free networks. J. Phys. A: Math. Theor. , 49(13):135001, April 2016
2016
-
[11]
Ralph Kenna’s scaling relations in critical phenomena
Le ¨ıla Moueddene, Arnaldo Donoso, and Bertrand Berche. Ralph Kenna’s scaling relations in critical phenomena. Entropy, 26(3), 2024
2024
-
[12]
Janke and R
W. Janke and R. Kenna. The Strength of First and Second Order Phase Transitions from Partition Function Zeroes. Journal of Statistical Physics, 102(5/6):1211–1227, 2001
2001
-
[13]
J. J. Ruiz-Lorenzo. Revisiting the Lee-Yang singularities in the four-dimensional Ising model: a tribute to the memory of Ralph Kenna. Condens. Matt. Phys., 27:33604, 2024
2024
-
[14]
Critical and tricritical singularities from small-scale Monte Carlo simulations: the Blume–Capel model in two dimensions
Le ¨ıla Moueddene, Nikolaos G Fytas, Yurij Holovatch, Ralph Kenna, and Bertrand Berche. Critical and tricritical singularities from small-scale Monte Carlo simulations: the Blume–Capel model in two dimensions. J. Stat. Mech., 2024(2):023206, February 2024
2024
-
[15]
Fytas, and Bertrand Berche
Le ¨ıla Moueddene, Nikolaos G. Fytas, and Bertrand Berche. Critical and tricritical behavior of the𝑑= 3 Blume-Capel model: Results from small-scale Monte Carlo simulations. Phys. Rev. E, 110:064144, Dec 2024
2024
-
[16]
T. D. Lee and C. N. Yang. Statistical Theory of Equations of State and Phase Transitions. II. Lattice Gas and Ising Model. Phys. Rev., 87(3):410–419, August 1952. 24
1952
-
[17]
C. N. Yang and T. D. Lee. Statistical Theory of Equations of State and Phase Transitions. I. Theory of Condensation. Phys. Rev., 87(3):404–409, August 1952
1952
-
[18]
W. E. Brittin, editor. Lectures in Theoretical Physics (vol. VIIC) by Fisher M.E. Gordon and Breach, New York, 1968
1968
-
[19]
Berche, R
B. Berche, R. Kenna, and J.-C. Walter. Hyperscaling above the upper critical dimension.Nuclear Physics B, 865(1):115–132, 2012
2012
-
[20]
Butera and M
P. Butera and M. Pernici. The Blume–Capel model for spins S = 1 and 3 / 2 in dimensions d = 2 and 3. Physica A: Statistical Mechanics and its Applications, 507:22–66, October 2018
2018
-
[21]
Rocha-Neto, G
M ´ario J.G. Rocha-Neto, G. Camelo-Neto, E. Nogueira, and S. Coutinho. Thermodynamical behavior of the Blume–Capel model in the vicinity of its tricritical point. Physica A: Statistical Mechanics and its Applications, 629:129145, November 2023
2023
-
[22]
Phase transitions above the upper critical dimension
Bertrand Berche, Tim Ellis, Yurij Holovatch, and Ralph Kenna. Phase transitions above the upper critical dimension. SciPost Phys. Lect. Notes, page 60, 2022
2022
-
[23]
Biskup, C
M. Biskup, C. Borgs, J. T. Chayes, L. J. Kleinwaks, and R. Koteck´ y. General Theory of Lee-Yang Zeros in Models with First-Order Phase Transitions. Phys. Rev. Lett., 84(21):4794–4797, May 2000. 25
2000
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.