Pith. sign in

REVIEW 2 major objections 3 minor 23 references

Ensemble inequivalence in the Blume-Emery-Griffiths model near a fourth order critical point

T0 review · 2 major / 3 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 microcanonical fourth-order critical point of the BEG model exists, but at different parameters than the canonical one.

desk verdict Microcanonical BEG at K<0 is a solid extension, but the fourth-order point hangs on a branch selection the paper asserts rather than demonstrates. read the letter →

arxiv 1908.07770 v2 pith:BOYVLSEH submitted 2019-08-21 cond-mat.stat-mech

V. V. Prasad , Alessandro Campa , David Mukamel , Stefano Ruffo This is my paper · ORCID
classification cond-mat.stat-mech
keywords ensembleinequivalencemicrocanonicalBlume-Emery-Griffithsmodelfourth-ordercriticalpointlong-rangeinteractionsreentrantphasetransitionstricriticalmean-fieldspin
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 asks whether a fourth-order critical point known to exist in the canonical (fixed-temperature) Blume-Emery-Griffiths model survives when the same model is studied at fixed energy. It establishes that it does survive, but in a different location: the microcanonical fourth-order point sits at higher temperature and a different negative biquadratic coupling than the canonical one. The reason this matters is that a high-order critical point organizes the topology of the whole surrounding phase diagram, so the shift changes which transitions are first order, which are continuous, and where they meet. The microcanonical diagram contains reentrant first-order transitions and temperature discontinuities that the canonical diagram cannot show. A sympathetic reader should take the paper as mapping out this ensemble-dependent multicritical topology in one exactly solvable mean-field model.

What carries the argument

The load-bearing object is the microcanonical entropy per spin $\tilde{s}_+(\epsilon,m)$, obtained by counting configurations with given magnetization $m$ and quadrupole moment $q$, solving the energy relation for $q$ as a function of $m$ and $\epsilon$, and then expanding around the $m=0$ branch: $\tilde{s}_+=s_0+A_m m^2+B_m m^4+C_m m^6+D_m m^8+\cdots$. The vanishing of $A_m$ defines the critical surface, $A_m=B_m=0$ defines the tricritical line, and $A_m=B_m=C_m=0$ with $D_m<0$ defines the fourth-order point; the same coefficient hierarchy in the canonical free energy locates the canonical fourth-order point. This expansion turns the search for the multicritical point into a finite algebraic calculation along the tricritical line.

What would settle it

Compute the exact global maximum of $\tilde{s}_+(\epsilon,m)$ over $m$ at the parameter values of the alleged fourth-order point, especially in a window around $\epsilon_1^*\approx 0.0835$; if the $m=0$ branch is the entropy maximum there, or if the finite-magnetization branch does not outrank it, the claimed location of the microcanonical fourth-order point is wrong, and if the ordering is reversed the claim is supported.

Watch

Extended reading notes

Core claim

The paper's central claim is that the infinite-range Blume-Emery-Griffiths model with negative biquadratic coupling $K<0$ has a fourth-order critical point in the microcanonical ensemble at $(\epsilon_2^*\approx 0.1313,\ \Delta_2^*\approx 0.4369,\ K_2^*\approx -0.0828,\ T^*\approx 0.2924)$, which is distinct from the canonical fourth-order point at $(T^*\approx 0.2402,\ K^*\approx -0.1838,\ \Delta^*\approx 0.399)$. The microcanonical point is found by locating the simultaneous vanishing of the $m^2$, $m^4$ and $m^6$ coefficients in an expansion of the entropy about $m=0$, with the $m^8$ coefficient negative. A second candidate at $\epsilon_1^*\approx 0.0835$ is argued to be preempted by a global entropy maximum away from $m=0$. Around the accepted point the phase diagram has a continuous transition line ending in a critical end point and a reentrant first-order line that enters the ordered phase and separates two ferromagnetic ordered phases; as energy rises at fixed $\Delta$, this produces sequences of first-order, continuous, and again first-order transitions, with temperature discontinuities at the first-order steps.

Load-bearing premise

The identification of the fourth-order point rests on rejecting the lower-energy candidate $\epsilon_1^*\approx 0.0835$ because a nonzero-magnetization entropy maximum is claimed to beat it, but that global entropy comparison is not shown in the paper.

Editorial extensions

If this is right

  • Raising the energy at fixed $\Delta$ and $K$ in the reentrant region takes an ordered state into a disordered one, back into an ordered state through a continuous transition, and then into a second ordered state through another first-order transition.
  • Every microcanonical first-order transition has a temperature jump; at the point where two first-order branches merge the jump disappears.
  • The continuous transition line is the same in both ensembles, but the first-order lines, critical end points, and ordered-phase critical points differ, so ensemble inequivalence is localized precisely where first-order behavior occurs.
  • Since the topology near a high-order critical point persists over a broad parameter range, the existence of two different fourth-order points means the two ensembles disagree throughout a finite neighborhood in $(T,\Delta,K)$ space, not merely at one point.

Reading between the lines

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

  • An immediate testable extension is to compute the canonical phase diagram at the microcanonical fourth-order parameters and verify that no fourth-order singularity appears there; the two ensembles cannot simultaneously host the point.
  • The same entropy-expansion criterion applied to other infinite-range models with high-order multicritical points would predict microcanonical points shifted in temperature relative to canonical ones, and in models with more order parameters the shift could change which ordered phases participate.
  • The reentrant caloric curves imply narrow energy windows with negative specific heat; these should be visible in constant-energy Monte Carlo or molecular-dynamics runs as a decreasing segment in the temperature-versus-energy curve.
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

2 major / 3 minor

Summary. The paper extends the microcanonical analysis of the infinite-range Blume-Emery-Griffiths model to negative biquadratic coupling K, where the canonical ensemble is known to have a fourth-order critical point. The authors derive the microcanonical entropy as a function of energy and magnetization, expand it about m=0, and obtain the critical and tricritical surfaces. Solving Am=Bm=Cm=0 along the tricritical line yields two candidate energies, epsilon_1*≈0.0835 and epsilon_2*≈0.1313; the paper selects epsilon_2* as the microcanonical fourth-order point, with (epsilon_2*, Delta_2*, K_2*)≈(0.1313, 0.4369, -0.0828) and T*≈0.2924. The paper then studies the phase diagram at K=-0.4, reporting reentrant first-order transitions, temperature discontinuities, and critical end points that differ from the canonical ones, and it superimposes the two ensembles' phase diagrams for comparison.

Significance. If the selection of epsilon_2* is correct, the paper provides a concrete example of ensemble inequivalence at a fourth-order critical point, with distinct microcanonical and canonical coordinates and with topological features such as reentrant ordered phases absent in the canonical ensemble. The work is not circular: the multicritical point is solved from the entropy expansion coefficients, no constants are fitted to data, and the printed coefficients give an explicit route to the central result. The quantitative comparison at K=-0.4 is a falsifiable prediction. The main gap is the absence of a displayed global-entropy comparison at the discarded solution epsilon_1*, which is necessary to identify the fourth-order point in the microcanonical ensemble.

major comments (2)
  1. [Sec. III, after Eq. (29) and Fig. 4] The choice of epsilon_2* as the fourth-order point is not established. Figure 4 shows Cm=0 and Dm<0 at both epsilon_1* and epsilon_2*, i.e. local stability of the m=0 branch at both candidates. Microcanonical equilibrium, however, is defined by the global maximum of s+(epsilon,m), as stated earlier in this section. The text asserts that epsilon_1* is 'preempted by a global maximum away from m=0', but no entropy comparison at epsilon_1* is shown. If the m=0 branch happened to be the global maximum at epsilon_1*, the tricritical line would terminate at epsilon_1* and the coordinates in Eq. (29) would be wrong. Please provide an explicit global-maximization check, for example a plot of max_m s+(epsilon,m) and of s+(epsilon,0) along the tricritical line, or a direct numerical comparison at epsilon_1* and epsilon_2*.
  2. [Sec. IV, Figs. 6-9] The quantitative phase diagrams used to illustrate the distinct microcanonical behavior are computed for K=-0.4, which is far from the claimed fourth-order point K_2*≈-0.0828. The reentrant transitions and temperature discontinuities are therefore demonstrated in a different regime than the point identified in Eq. (29). The paper should either include representative calculations at K values close to K_2*, or explicitly justify continuity of the phase-diagram topology from K_2* down to K=-0.4. This does not invalidate the local K=-0.4 results, but it limits the support for the title's claim of behavior 'near a fourth order critical point'.
minor comments (3)
  1. [Sec. II, after Eq. (7)] The sentence saying that the values x and iy which minimize beta * f_tilde correspond to m and q is imprecise; please state explicitly that the saddle-point value of x gives m and the saddle-point value of iy gives q.
  2. [Sec. III, Eq. (27)] Equation (27) uses partial q+/partial epsilon before q+ is reintroduced after Eq. (22); please define q+ once more immediately before this expression to make the notation self-contained.
  3. [Sec. IV, Fig. 7] The points P^{MC}_1, P^{MC}_2 and P^{MC}_3 are defined in the figure caption, but the main text refers to them without a formal definition; consider introducing them explicitly in the text before the first use.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the microcanonical fourth-order point is solved from the entropy expansion rather than fitted or imported.

full rationale

The central microcanonical claim is derived self-contained. The paper expands the exact entropy s̃+(ε,m) about m=0, obtains the coefficients Am, Bm, Cm, Dm, solves Am=Bm=Cm=0 along the tricritical line, and takes Dm<0 as the local stability condition (Sec. III, Eqs. (23)-(28), Fig. 4). The claimed location (ε2*, Δ2*, K2*, T*) is the algebraic solution of these equations, not a fitted parameter, and the canonical fourth-order point (Eq. (17)) is derived independently in Sec. II and used only for comparison. The paper does rely on prior canonical results and on the authors' earlier ensemble-inequivalence framework, but the K<0 microcanonical calculation is presented with its own entropy function, expansion coefficients, and global phase-diagram construction. The text states that ε1* is 'preempted by a global maximum away from m=0' and that 'the only solution which corresponds to a global maximum of the entropy is ε2*,' but no explicit entropy comparison at ε1* is shown. That is an unverified global-maximum step and a potential correctness gap, not a circular reduction: the assertion does not define or fit the fourth-order point. Under the hard rule that circularity must be exhibited as Eq. X = Eq. Y by construction or as a fitted parameter renamed as a prediction, no circular step is present. Therefore the paper receives a non-circular score.

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

The model is fixed by the Hamiltonian; the derivation introduces no fitted parameters. The main auxiliary premises are the q+ branch selection, the validity of the eighth-order entropy expansion, and the unproved rejection of epsilon_1* as a global maximum. The latter is the main source of uncertainty.

assumptions (5)
  • domain assumption The microcanonical equilibrium state maximizes the entropy at fixed energy and is parameterized by m and q.
    Standard microcanonical thermodynamics; the paper uses this to define s+(epsilon)=max_m s~+(epsilon,m).
  • domain assumption For K<0 only the q+ solution of Eq. (22) is physical because q- is negative.
    The paper discards q- on sign grounds; it does not analyze the q<=1 bound or the full energy domain.
  • standard math The entropy expansion Eq. (23) to order m^8, with coefficients Eq. (25), is valid near m=0 and controls stability at the critical point.
    The coefficients are stated without derivation; the identification of the fourth-order point uses signs of C_m and D_m.
  • ad hoc to paper The epsilon_1* about 0.0835 solution of C_m=0 is not the global entropy maximum and can be discarded.
    Asserted in Sec. III with deferred support; the displayed global entropy plots are for K=-0.4, not for epsilon_1*.
  • standard math The canonical free-energy expansion coefficients B_c, C_c, D_c in Sec. II are as stated.
    The canonical fourth-order point is used as the benchmark; these coefficients are presented without full derivation.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Ensemble inequivalence in the Blume-Emery-Griffiths model near a fourth order critical point." pith.science (2026). https://pith.science/paper/BOYVLSEH

@misc{pith2026190807770,
  author       = {Pith},
  title        = {Pith review of: Ensemble inequivalence in the Blume-Emery-Griffiths model near a fourth order critical point},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/BOYVLSEH}},
  note         = {Machine review of arXiv:1908.07770}
}
abstract

The canonical phase diagram of the Blume-Emery-Griffiths (BEG) model with infinite-range interactions is known to exhibit a fourth order critical point at some negative value of the bi-quadratic interaction $K<0$. Here we study the microcanonical phase diagram of this model for $K<0$, extending previous studies which were restricted to positive $K$. A fourth order critical point is found to exist at coupling parameters which are different from those of the canonical ensemble. The microcanonical phase diagram of the model close to the fourth order critical point is studied in detail revealing some distinct features from the canonical counterpart.

Figures

Figures reproduced from arXiv: 1908.07770 by the authors.

Figure 1
Figure 1. FIG. 1. (color online) Schematic plot showing the canonical [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. (color online) The canonical phase diagram in the [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. (color online) The tricritical line obtained from the micro-canonical analysis by solving for [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (5 more)
Figure 4
Figure 4. Figure 4: FIG. 4. (color online) The values of the coefficients [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 6
Figure 6. Figure 6: FIG. 6. (color online) The microcanonical phase diagrams in [PITH_FULL_IMAGE:figures/full_fig_p006_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Panel (A): Schematic phase diagram in the (∆ [PITH_FULL_IMAGE:figures/full_fig_p007_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. The temperature profile (caloric curve) for fixed val [PITH_FULL_IMAGE:figures/full_fig_p007_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. The canonical and microcanonical (∆ [PITH_FULL_IMAGE:figures/full_fig_p008_9.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

23 extracted references · 21 canonical work pages

  1. [1]

    Campa, T

    A. Campa, T. Dauxois, D. Fanelli, and S. Ruffo, Physics of Long-Range Interacting Systems (Oxford University Press, Oxford, 2014)

    work page 2014

  2. [2]

    cbYz1JnEbAXTGe6NJnkM6cWOicE=

    The other solution is preempted by a global max- imum away from m = 0. Thus the fourth order critical point of the microcanonical ensemble takes place at ϵ∗ 2≈ 0.1313, ∆∗ 2≈ 0.4369, K∗ 2≈− 0.0828, (29) 0.06 0.08 0.10 0.12 0.14 ϵ −2 −1 0 1 2 ϵ∗ 1 ϵ∗ 2 Dm Cm FIG. 4. (color online) The values of the coefficients Cm and Dm [see Eq. (23)], plotted as a function ...

    work page

  3. [3]

    Dauxois, S

    T. Dauxois, S. Ruffo, E. Arimondo, and M. Wilkens (Eds.), Dynamics and Thermodynamics of Systems with Long-Range Interactions, Lecture Notes in Physics Vol. 602 (Springer, New York, 2002)

    work page 2002

  4. [4]

    Chavanis, in Dynamics and Thermodynamics of Systems with Long-Range Interactions , Lecture Notes in Physics Vol

    P.-H. Chavanis, in Dynamics and Thermodynamics of Systems with Long-Range Interactions , Lecture Notes in Physics Vol. 602, edited by T. Dauxois, S. Ruffo, E. Ari- mondo, and M. Wilkens (Springer, New York, 2002) p. 208

    work page 2002

  5. [5]

    Padmanabhan, Phys

    T. Padmanabhan, Phys. Rep. 188, 285 (1990)

    1990

  6. [6]

    L. D. Landau and E. M. Lifshitz, Electrodynamics of Continuous Media (Pergamon Press, Oxford, 1960)

    work page 1960

  7. [7]

    D. R. Nicholson, Introduction to Plasma Theory (Krieger Publishing Company, 1992)

    work page 1992

  8. [8]

    Campa, T

    A. Campa, T. Dauxois, and S. Ruffo, Phys. Rep. 480, 57 (2009)

    work page 2009

Show all 23 references
  1. [9]

    Bouchet, S

    F. Bouchet, S. Gupta, and D. Mukamel, Physica A 389, 4389 (2010)

    2010

  2. [10]

    Misawa, Y

    T. Misawa, Y. Yamaji, and M. Imada, J. Phys. Soc. Japan 75, 064705 (2006)

    2006

  3. [11]

    Borgonovi, G

    F. Borgonovi, G. L. Celardo, M. Maianti, and E. Peder- soli, J. Stat. Phys. 116, 1435 (2004)

    2004

  4. [12]

    Mukamel, S

    D. Mukamel, S. Ruffo, and N. Schreiber, Phys. Rev. Lett. 95, 240604 (2005)

    2005

  5. [13]

    Lynden-Bell, Mon

    D. Lynden-Bell, Mon. Not. R. Astron. Soc. 136, 101 (1967)

    1967

  6. [14]

    Chavanis, J

    P.-H. Chavanis, J. Sommeria, and R. Robert, Astrophys. J. 471, 385 (1996)

    1996

  7. [15]

    Latora, A

    V. Latora, A. Rapisarda, and S. Ruffo, Phys. Rev. Lett. 83, 2104 (1999)

    1999

  8. [16]

    Y. Y. Yamaguchi, J Barr´ e, F. Bouchet, T. Dauxois, and S. Ruffo, Physica A 337, 36 (2004)

    2004

  9. [17]

    Blume, V

    M. Blume, V. J. Emery, and R. B. Griffiths, Phys. Rev. A 4, 1071 (1971)

    1971

  10. [18]

    Mukamel and M

    D. Mukamel and M. Blume, Phys. Rev. A 10, 610 (1974)

    1974

  11. [19]

    Krinsky and D

    S. Krinsky and D. Mukamel, Phys. Rev. B11, 399 (1975)

    1975

  12. [20]

    Lajzerowicz and J

    J. Lajzerowicz and J. Sivardi` ere, Phys. Rev. A 11, 2079 (1975)

    1975

  13. [21]

    Hoston and A

    W. Hoston and A. N. Berker, Phys. Rev. Lett. 67, 1027 (1991)

    1991

  14. [22]

    V. V. Hovhannisyan, N. S. Ananikian, A. Campa, and S. Ruffo, Phys. Rev. E 96, 062103 (2017)

    2017

  15. [23]

    Barr´ e, D

    J. Barr´ e, D. Mukamel, and S. Ruffo, Phys. Rev. Lett.87, 030601 (2001)

    2001

Pith tools

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