REVIEW 2 major objections 5 minor 31 references
Mixed phases in feedback Ising models
T0 review · 2 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read When spin-spin coupling grows with magnetization, a mean-field Ising model develops stable mixed phases at zero temperature and a four-way classification of the transitions among them.
desk verdict A clean, checkable mean-field study of magnetization-dependent coupling; the zero-temperature sliding-dynamics gap is real but patchable. 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 central object is the zero-temperature Glauber dynamics $dm/dt = -m + \mathrm{sgn}(h - g(m))$, where $g(m) \equiv -m f(m) - \tfrac{1}{2} m^2 f'(m)$, together with the equilibrium curve $C_0 = \{(h,m): h = g(m)\}$. The stability of a mixed-phase point $m_0$ on $C_0$ is decided by the slope $g'(m_0)$: a positive slope makes the point super-stable, with linear decay rates $-m_0 \pm 1$, while a negative slope makes it super-unstable, with linear growth. For linear feedback, $g(m)$ is a parabola whose vertex at $m = -1/(3\gamma)$ creates the stable and unstable branches when $\gamma > 1/3$. The trapping condition of Eq. (18), bounding the field derivative between $(-m_0-1)g'(m_0)$ and $(-m_0+1)g'(m_0)$, is what guarantees that perturbations decay linearly and that the system stays on a stable mixed branch.
What would settle it
Solve or simulate the finite-temperature Glauber equation $dm/dt = -m + \tanh(\beta(h - g(m)))$ for the linear FIM with $\gamma > 1/3$ and any small positive $\beta$; if the stable mixed branch $m_1$ disappears or perturbations decay exponentially rather than linearly, then the paper's super-stability claim is confined to exactly $T = 0$.
Extended reading notes
Core claim
In the mean-field feedback Ising model with Hamiltonian $\hat{H}(m) = -hm - \tfrac{1}{2} f(m) m^2$ and linear coupling $f(m) = 1 + \gamma m$, the zero-temperature equilibrium condition is $h = g(m)$ with $g(m) = -m(\tfrac{3}{2}\gamma m + 1)$. The central claim is that for $\gamma > 1/3$, this equilibrium curve has a vertex at $m = -1/(3\gamma)$, giving a stable mixed-phase branch $m_1$ below the vertex and an unstable mixed-phase branch $m_2$ above it, even though the coupling remains ferromagnetic over most of the magnetization range. Stable mixed phases are always super-stable, meaning perturbations decay linearly and reach zero in finite time; the trapping condition of Eq. (18) determines when a time-dependent field keeps the system on such a branch. The paper also shows that stable mixed phases can be true ground states only when $\gamma > 1$, and that the linear FIM exhibits all four types of phase transitions at its branch endpoints, including a three-stage sequence $m_- \to m_1 \to m_+$.
Load-bearing premise
The entire stability classification is carried out at exactly zero temperature, where the thermal transition function is replaced by a sign function, and it assumes the trajectory has a well-defined sliding motion on the discontinuous equilibrium curve; if smoothing the discontinuity at finite temperature changes which mixed phases are stable, the central claim holds only in the $\beta \to \infty$ limit.
Editorial extensions
If this is right
- A fully ferromagnetic coupling can support stable mixed phases at zero temperature, which the classical Curie-Weiss model cannot (there, stable mixed phases require antiferromagnetic coupling).
- Stable mixed phases are super-stable, so under a slowly varying field the system follows the stable branch until the field derivative violates the trapping condition in Eq. (18), giving a concrete criterion for when a driven system is thrown off equilibrium.
- The linear FIM exhibits all four transition types and the three-stage sequence $m_- \to m_1 \to m_+$, providing a minimal model for stage-wise transformations such as those described by Comte's law of three stages.
- Because Eq. (20) reconstructs a feedback function $f(m)$ from any single-valued bifurcation diagram, the framework offers a data-driven route to modeling multistable systems from observed phase diagrams alone.
- Maxwell constructions show that stable mixed phases become true ground states in the linear FIM only when $\gamma > 1$, at which point the coupling is partly antiferromagnetic.
- A finite-temperature simulation of the Glauber equation for the linear FIM with $\gamma > 1/3$ would settle whether the stable mixed branch and its linear decay persist away from $\beta = \infty$.
Reading between the lines
- The paper's stability criterion is stated in terms of the sign function, so the entire super-stability result is framed at exactly zero temperature; a natural extension would be to compute finite-$\beta$ corrections and check whether the fold and the linear decay rate survive thermal smoothing.
- Because the same slope condition $g'(m)>0$ is derived for general feedback functions, any empirically reconstructed bifurcation diagram with a stable intermediate branch can be used to predict the feedback strength needed to stabilize that branch, which may be useful in ecological or socio-economic tipping-point applications.
- The three-stage interpretation suggests a testable minimal model: if a community's coupling increases linearly with the fraction of positive states, the intermediate stage loses stability once that fraction reaches roughly $1/3$, a prediction that could in principle be compared with historical or observational transition data.
- The equivalence between the linear FIM and a model with two- and three-spin interactions means the stable mixed phases could be realized in spin systems with engineered multispin couplings, providing a concrete experimental or numerical target.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces a mean-field Ising model with a magnetization-dependent coupling f(m), concentrates on the linear feedback case f(m)=1+γm, and studies the zero-temperature limit in which the Glauber dynamics reduces to dm/dt=−m+sgn(h−g(m)). It identifies a curve C0 of mixed-phase (MP) equilibria, classifies their stability by the sign of g′(m), and claims that stable MPs are super-stable, meaning that perturbations decay linearly in time and vanish in finite time. It also derives a trapping condition for time-varying magnetic fields, classifies four types of phase transitions, computes Maxwell points, and concludes that stable MPs can be ground states only for γ>1. The paper is analytic, with numerical simulations used to illustrate the predicted trajectories.
Significance. If the main claims hold, the feedback Ising model is a useful exactly solvable extension of the Curie-Weiss model: it shows that stable mixed phases can exist at zero temperature with partly ferromagnetic coupling, something the classical mean-field Ising model does not allow. The derivations are explicit and checkable: the piecewise-linear flow in Eq. (13), the trapping condition in Eq. (18), and the Maxwell-point formulas in the End Matter are all transparent and internally consistent. The model has no fitted parameters; γ is a modeling parameter, and the dynamical-simulation parameters are fully specified in the End Matter. The main caveat is that the central stability statements are made on the discontinuity curve C0 of the sgn function, where the governing equation is not a well-defined ODE without an additional convention; this is a missing well-posedness statement rather than a demonstrated instability.
major comments (2)
- [Equilibria and their stability, Eq. (9)] Eq. (9) is not a well-defined ordinary differential equation on the curve C0, where h=g(m) and sgn(0) is undefined. With the common convention sgn(0)=0, only m=0 on C0 is an equilibrium of Eq. (9), so the entire stable-MP branch would not exist. The paper treats the whole curve C0 as a set of equilibria and derives the trapping condition Eq. (18); this is valid only under a Filippov sliding interpretation (or an equivalent convex-selection convention for sgn(0)). Please state this convention explicitly, justify it as the β→∞ limit of the smooth Glauber dynamics, and provide the Filippov sliding vector field from which Eq. (18) follows. This point is load-bearing because the existence and super-stability of mixed phases are the central claims of the paper.
- [Equilibria and their stability, finite-β regularization] The claimed linear-in-time, finite-time decay of perturbations to a stable MP is a singular-limit artifact that does not survive at any finite β. For finite β, near a point m0 with g′(m0)>0 and h=g(m0), the linearization of Eq. (4) has derivative F′(m0)=−1−β(1−m0^2)g′(m0)<0, so perturbations decay exponentially with a rate proportional to β, and the exact trapping condition Eq. (18) becomes only approximate. The stable branch itself does survive regularization, so the central existence claim is unaffected, but the paper should either state explicitly that super-stability is a property of the exact β→∞ sliding solution, quantify the finite-β crossover, or soften the abstract and main-text wording so that the finite-time linear decay is not presented as a finite-temperature phenomenon.
minor comments (5)
- [Abstract / Finite-dimensional FIMs] The abstract states that the paper discusses basic properties of finite-dimensional FIMs, but the body contains only two brief, speculative sentences about spatial patterns and no actual finite-dimensional analysis; please either add a short section on finite-dimensional FIMs or temper the abstract's claim.
- [Phase transitions] The terms 'second-order' and 'third-order' are used for the inexact type-1 and type-4 transitions without a formal definition of the order in terms of derivatives of a thermodynamic quantity; a brief definition or an asymptotic derivation would make the classification verifiable.
- [Inverse problem, Eq. (20)] The reconstruction formula f(m)=−(2/m^2)∫_0^m g(u)du is stated without naming the regularity conditions on g that guarantee a nonsingular f; please specify the mild assumptions mentioned in the text.
- [References] Reference [29] contains a formatting artifact in the author name ('Ho/suppress lyst'); please correct it.
- [Applications paragraph] The discussion of Comte's law of three stages is clearly labeled a possible application, but it is disconnected from the preceding mathematics; consider condensing it or explicitly mapping the model variables to the social stages in the text.
Circularity Check
Central model predictions are self-contained; the only circularity is a naming-level equivalence between 'stable MP' and 'super-stable'.
-
self definitional
[Abstract; Section 'Equilibria and their stability', Eq. (18)]
"Moreover, stable MPs are always super-stable, meaning that perturbations decay linearly in time. ... However, m0 is 'super-stable' when g′(m0) > 0 and h′(t0) satisfies a trapping condition h′(t0) ∈ ((−m0 − 1)g′(m0), (−m0 + 1)g′(m0)), (18) since perturbation decays linearly in t with rates..."
The paper defines 'super-stable' by exactly the condition g′(m0)>0 plus the trapping condition (18). For a static field h′(t0)=0, Eq. (18) is satisfied precisely when g′(m0)>0, which is the same condition that makes a C0 branch attracting (stable) in this zero-temperature Filippov picture. Thus the advertised result 'stable MPs are always super-stable' is true by construction: the term 'super-stable' was introduced to mean 'stable with linear-in-time decay', and the decay rates in Eq. (19) are the defining content of that term. This is a semantic/naming tautology rather than an independently derived theorem. It does not affect the substantive predictions, which are obtained by direct calculation from g(m).
full rationale
The paper's substantive derivation chain is self-contained. The linear FIM bifurcation diagram follows by substituting f(m)=1+γm into Eq. (6), giving g(m)=−m((3/2)γm+1); the vertex at m=−1/(3γ), h=1/(6γ) yields a stable MP branch for γ>1/3 by the sign of g′(m). No parameter is fitted to data and no external benchmark is used, so nothing is 'fitted input called prediction'. The phase-transition classification, trapping condition, and Maxwell points are derived in closed form from Eq. (9) and the Hamiltonian, not from any fitted quantity. Self-citations [4,5,12,25] are background or textbook references and are not load-bearing for the model's claims; there is no imported uniqueness theorem, no ansatz smuggled in via citation, and no renaming of an external empirical pattern as unification. The only identified circular element is the abstract's statement that stable MPs are always super-stable, which is effectively a restatement of the definition of super-stable; this is a naming-level issue, not a logical loop in the predictions. The separate well-posedness gap at the discontinuity h=g(m) in Eq. (9) (the paper does not state the Filippov convention that makes C0 a set of equilibria) is a rigor/correctness concern, not circularity, and is therefore not scored as a circular step.
Assumptions & free parameters
free parameters (1)
- γ (feedback strength) =
not fitted (model parameter)
assumptions (4)
- domain assumption Mean-field (fully connected) limit with N→∞
- domain assumption Zero-temperature limit β→∞ replacing tanh(βx) with sgn(x)
- domain assumption Glauber dynamics as the microscopic evolution rule
- ad hoc to paper Linear feedback ansatz f(m)=1+γm
Cite this review
Pith. "Pith review of Mixed phases in feedback Ising models." pith.science (2026). https://pith.science/paper/VJPNM34Y
@misc{pith2026250609025,
author = {Pith},
title = {Pith review of: Mixed phases in feedback Ising models},
year = {2026},
howpublished = {\url{https://pith.science/paper/VJPNM34Y}},
note = {Machine review of arXiv:2506.09025}
}
read the original abstract
We study mean-field Ising models in which the coupling depends on the magnetization via a feedback function. We identify mixed phases (MPs) and show that they can be stable at zero temperature for sufficiently strong feedback. Moreover, stable MPs are always super-stable, meaning that perturbations decay linearly in time. Feedback Ising models (FIMs) provide a useful framework for phase transformations between aligned phases via stable and unstable intermediate phases in multistable systems. We also analyze the dynamical behavior of FIMs driven by a varying magnetic field and discuss basic properties of finite-dimensional FIMs.
Figures
Reference graph
Works this paper leans on
-
[1]
J. P. Sethna, Statistical Mechanics: Entropy, Order Pa- rameters, and Complexity (Oxford University Press, New York, 2006)
work page 2006
-
[2]
Time-dependent statistics of the Ising model,
R. J. Glauber, “Time-dependent statistics of the Ising model,” J. Math. Phys. 4, 294–307 (1963)
work page 1963
-
[3]
Theory of phase-ordering kinetics,
A. J. Bray, “Theory of phase-ordering kinetics,” Adv. Phys. 43, 357–459 (1994)
work page 1994
-
[4]
Freezing in Ising ferromagnets,
V. Spirin, P. L. Krapivsky, and S. Redner, “Freezing in Ising ferromagnets,” Phys. Rev. E 65, 016119 (2001)
work page 2001
-
[5]
Zero- temperature relaxation of three-dimensional Ising ferro- magnets,
J. Olejarz, P. L. Krapivsky, and S. Redner, “Zero- temperature relaxation of three-dimensional Ising ferro- magnets,” Phys. Rev. E 83, 051104 (2011)
work page 2011
-
[6]
Superdiffusion-like behavior in zero-temperature coars- ening of the d = 3 Ising model,
D. Gessert, H. Christiansen, and W. Janke, “Superdiffusion-like behavior in zero-temperature coars- ening of the d = 3 Ising model,” Sci. Rep. 13, 13270 (2023)
work page 2023
-
[7]
Yasser Roudi, Joanna Tyrcha, and John Hertz, “Ising model for neural data: Model quality and approximate methods for extracting functional connectivity,” Phys. Rev. E 79, 051915 (2009)
work page 2009
-
[8]
Major depres- sion as a complex dynamic system,
Ang´ elique O. J. Cramer, Claudia D. van Borkulo, Erik J. Giltay, Han L. J. van der Maas, Kenneth S. Kendler, Marten Scheffer, and Denny Borsboom, “Major depres- sion as a complex dynamic system,” PLOS ONE 11 (2016), 10.1371/journal.pone.0167490
Show all 31 references
-
[9]
Stochastic and mesoscopic models for tropi- cal convection,
Andrew J. Majda and Boualem Khouider, “Stochastic and mesoscopic models for tropi- cal convection,” Proceedings of the National Academy of Sciences 99, 1123–1128 (2002), https://www.pnas.org/doi/pdf/10.1073/pnas.032663199
2002 doi
-
[10]
Stability of money: phase transitions in an ising economy,
Stefan Bornholdt and Friedrich Wagner, “Stability of money: phase transitions in an ising economy,” Phys- ica A: Statistical Mechanics and its Applications 316, 453–468 (2002)
2002
-
[11]
Neural networks and physical systems with emergent collective computational abilities,
J. J. Hopfield, “Neural networks and physical systems with emergent collective computational abilities,” Proc. Natl. Acad. Sci. USA 79, 2554–2558 (1982)
1982
-
[12]
Ising model for melt ponds on arctic sea ice,
Y.-P. Ma, I. Sudakov, C. Strong, and K. M. Golden, “Ising model for melt ponds on arctic sea ice,” New J. Phys. 21, 023004 (2019)
2019
-
[13]
social temperature
and collective decision-making [14] display analogous loops that either reinforce or oppose the prevailing order. Our original goal was to enrich the Curie-Weiss model with feedback. The equilibrium and dynamical charac- teristics of the feedback Ising model (FIM) are non-triv...
2025 arXiv
-
[14]
Spatial patterns of tree yield explained by endogenous forces through a correspondence between the ising model and ecology,
A. E. Noble, T. S. Rosenstock, P. H. Brown, J. Machta, and A. Hastings, “Spatial patterns of tree yield explained by endogenous forces through a correspondence between the ising model and ecology,” Proceedings of the National Academy of Sciences of the United States of America...
2018
-
[15]
Ising-based model of opinion formation in a complex network of interpersonal interactions,
A. Grabowski and R. A. Kosi´ nski, “Ising-based model of opinion formation in a complex network of interpersonal interactions,” Physica A: Statistical Mechanics and its Applications 361, 651–664 (2006)
2006
-
[16]
Marc Mezard and Andrea Montanari, Information, physics, and computation (Oxford University Press, 6 2009)
2009
-
[17]
Gapless spin-fluid ground state in a random quantum heisenberg magnet,
Subir Sachdev and Jinwu Ye, “Gapless spin-fluid ground state in a random quantum heisenberg magnet,” Physical review letters 70, 3339 (1993)
1993
-
[18]
The soft mode in the sachdev-ye-kitaev model and its gravity dual,
Alexei Kitaev and S Josephine Suh, “The soft mode in the sachdev-ye-kitaev model and its gravity dual,” Jour- nal of High Energy Physics 2018, 1–68 (2018)
2018
-
[19]
Energy gaps in quantum first-order mean-field–like transitions: The problems that quantum annealing cannot solve,
Thomas J¨ org, Florent Krzakala, Jorge Kurchan, An- thony C Maggs, and Justine Pujos, “Energy gaps in quantum first-order mean-field–like transitions: The problems that quantum annealing cannot solve,” Euro- physics Letters 89, 40004 (2010)
2010
-
[20]
Reverse annealing for the fully connected p-spin model,
Masaki Ohkuwa, Hidetoshi Nishimori, and Daniel A Li- dar, “Reverse annealing for the fully connected p-spin model,” Physical Review A 98, 022314 (2018)
2018
-
[21]
Dynamics of reverse annealing for the fully connected p-spin model,
Yu Yamashiro, Masaki Ohkuwa, Hidetoshi Nishimori, and Daniel A Lidar, “Dynamics of reverse annealing for the fully connected p-spin model,” Physical Review A 100, 052321 (2019)
2019
-
[22]
Fluctuations of the magnetization in the p- spin curie–weiss model,
Somabha Mukherjee, Jaesung Son, and Bhaswar B Bhat- tacharya, “Fluctuations of the magnetization in the p- spin curie–weiss model,” Communications in Mathemat- ical Physics 387, 681–728 (2021)
2021
-
[23]
Expectation bubbles in a spin model of markets,
S. Bornholdt, “Expectation bubbles in a spin model of markets,” Int. J. Mod. Phys. C 12, 667–674 (2001)
2001
-
[24]
Opinion for- mation model for markets with a social temperature and fear,
Sebastian M Krause and Stefan Bornholdt, “Opinion for- mation model for markets with a social temperature and fear,” Physical Review E—Statistical, Nonlinear, and Soft Matter Physics 86, 056106 (2012)
2012
-
[25]
Feedback-induced self-oscillations in large interacting systems subjected to phase transitions,
D. De Martino, “Feedback-induced self-oscillations in large interacting systems subjected to phase transitions,” J. Phys. A: Math. Theor. 52, 015001 (2019)
2019
-
[26]
Krapivsky, Sidney Redner, and Eli Ben-Naim, A Kinetic View of Statistical Physics (Cambridge Uni- versity Press, 2010)
Pavel L. Krapivsky, Sidney Redner, and Eli Ben-Naim, A Kinetic View of Statistical Physics (Cambridge Uni- versity Press, 2010)
2010
-
[27]
Auguste Comte, The positive philosophy of Auguste Comte (Blanchard, 1858)
-
[28]
Control of multistability,
Alexander N Pisarchik and Ulrike Feudel, “Control of multistability,” Physics Reports 540, 167–218 (2014)
2014
-
[29]
Spatial localization in dissipative systems,
E Knobloch, “Spatial localization in dissipative systems,” Annual Review of Condensed Matter Physics 6, 325–359 (2015)
2015
-
[30]
The ising model celebrates a cen- tury of interdisciplinary contributions,
Michael W. Macy, Boleslaw K. Szymanski, and Janusz A. Ho/suppress lyst, “The ising model celebrates a cen- tury of interdisciplinary contributions,” npj Complex 1 (2024), 10.1038/s44260-024-00012-0
2024 doi
-
[31]
Physics of the cryosphere,
Alison F. Banwell, Justin C. Burton, Claudia Cenedese, et al. , “Physics of the cryosphere,” Nature Reviews Physics 5, 446–449 (2023). End Matter Algorithm and parameters for dynamics. —In Fig. 3, the governing equation (4) is solved using the MATLAB solver ode15s for stiff ord...
2023
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Sign in to comment.