Strong Mpemba Effect Through a Reentrant Phase Transition
Pith reviewed 2026-05-07 05:34 UTC · model grok-4.3
The pith
Quenches to the paramagnetic phase in the antiferromagnetic Ising model produce a strong Mpemba effect via a staggered slow mode.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
We show that the strong direct and inverse Mpemba effects arise when quenches terminate in the paramagnetic phase of the antiferromagnetic Ising model in a magnetic field. These anomalous relaxation phenomena originate from the selective excitation of the slowest relaxation mode, which in the paramagnetic phase is purely staggered. Consequently, quenches starting from the paramagnetic phase have zero overlap with the slow mode and exhibit a strong (inverse) Mpemba effect. Quenches from the antiferromagnetic phase excite the staggered mode and display a slow-relaxation tail. By varying the lattice coordination number we show that the strong Mpemba effect disappears in the absence of reentranc
What carries the argument
The staggered relaxation mode that becomes slowest in the paramagnetic phase after the reentrant transition, with its overlaps computed under the pair approximation.
If this is right
- Quenches terminating in the paramagnetic phase have zero overlap with the slowest mode and therefore relax without the slow tail.
- Quenches that begin in the antiferromagnetic phase excite the staggered mode and therefore exhibit a pronounced slow-relaxation tail.
- The strong Mpemba effect vanishes when the coordination number is changed so that reentrance disappears.
- The pair approximation is sufficient to establish the direct link between the reentrant transition and the anomalous relaxation behavior.
Where Pith is reading between the lines
- The same selective excitation of a staggered slow mode could appear in other lattice models or real materials that possess reentrant transitions, offering a route to phase-controlled relaxation speeds.
- Experiments on antiferromagnetic samples in applied fields could directly test whether starting temperature determines the presence of the slow tail after a quench.
- Mapping equilibrium phase diagrams may become a practical tool for predicting and designing fast-relaxation protocols in systems far from equilibrium.
Load-bearing premise
The pair approximation accurately captures the relaxation modes and their overlaps for the antiferromagnetic Ising model across the reentrant transition.
What would settle it
Numerical simulations of the full dynamics of the antiferromagnetic Ising model that find nonzero overlap between paramagnetic initial conditions and the slowest mode, or that retain the strong Mpemba effect even when reentrance is removed, would falsify the claim.
Figures
read the original abstract
We investigate temperature quenches across the reentrant phase transition of the antiferromagnetic Ising model in a magnetic field and show that the strong direct and inverse Mpemba effects arise when quenches terminate in the paramagnetic phase. These anomalous relaxation phenomena originate from the selective excitation of the slowest relaxation mode, which in the paramagnetic phase is purely staggered. Consequently, quenches starting from the paramagnetic phase have zero overlap with the slow mode and exhibit a strong (inverse) Mpemba effect. Quenches from the antiferromagnetic phase excite the staggered mode and display a slow-relaxation tail. By varying the lattice coordination number we show that the strong Mpemba effect disappears in the absence of reentrance. Our results provide the first demonstration of the strong (inverse) Mpemba effect in the antiferromagnetic Ising model based on the pair-approximation, and establish a link between anomalous relaxation and the equilibrium phase behavior.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates temperature quenches across the reentrant phase transition in the antiferromagnetic Ising model in a magnetic field within the pair approximation. It claims that strong direct and inverse Mpemba effects occur precisely when quenches terminate in the paramagnetic phase, because the slowest relaxation mode there is purely staggered (zero overlap with uniform magnetization). Quenches from the antiferromagnetic phase excite this mode and show slow tails, while the effect vanishes when reentrance is suppressed by reducing coordination number. The results are derived from linearizing the closed dynamical equations for one- and two-point functions and analyzing eigenvalue overlaps.
Significance. If the pair-approximation spectrum is faithful, the work supplies a concrete mechanism linking the strong Mpemba effect to equilibrium phase structure via selective mode excitation, and it is the first such demonstration for the antiferromagnetic Ising model. The coordination-number control provides a falsifiable link between reentrance and anomalous relaxation that could be tested in other mean-field or lattice models.
major comments (3)
- [§3.2] §3.2 (pair-approximation closure and mode analysis): The central claim that the slowest eigenvalue in the paramagnetic phase has strictly zero overlap with the uniform magnetization (hence exact strong inverse Mpemba for paramagnetic initial conditions) rests on the pair closure. Because the closure truncates the BBGKY hierarchy at the two-spin level, it can artificially decouple staggered and uniform sectors; the reentrant transition itself is also a mean-field artifact of the same closure. No comparison to the exact master equation on small clusters (e.g., 2×2 or 4-site lattices) or to Monte Carlo trajectories is reported, leaving the zero-overlap result unverified.
- [§4.3] §4.3 (coordination-number scan): The disappearance of the strong Mpemba effect for coordination numbers that eliminate reentrance is presented as evidence that reentrance is necessary. However, the same scan simultaneously changes the accuracy of the pair approximation (higher z improves mean-field character), so the observed correlation could be an artifact of the closure rather than a generic feature of the model. A direct comparison at fixed z between pair-approximation and Monte Carlo relaxation spectra would be required to separate these effects.
- [§5] §5 (quantitative Mpemba measures): The reported crossing times and relaxation-rate ratios that quantify the “strong” Mpemba effect are obtained from the linearized pair dynamics without error bars or sensitivity analysis to the closure parameters. Because the pair approximation is uncontrolled for dynamics, it is unclear whether the reported effect sizes survive restoration of higher-order correlations.
minor comments (3)
- [Figure 2] Figure 2: the eigenvector components plotted for the slowest mode should include an explicit statement that the uniform-magnetization projection is numerically zero to machine precision, together with the tolerance used.
- [§2 and §3] Notation: the symbols for the staggered and uniform magnetizations are introduced inconsistently between the equilibrium phase diagram (§2) and the dynamical equations (§3); a single consistent notation table would help.
- [Introduction] The abstract states “first demonstration … based on the pair-approximation”; a brief sentence in the introduction citing prior Mpemba studies on Ising models (even if they used different approximations) would clarify the novelty claim.
Simulated Author's Rebuttal
We thank the referee for the thorough review and valuable comments. We will revise the manuscript to clarify the limitations of the pair approximation and to include additional analysis addressing the concerns raised. Our responses to each major comment are provided below.
read point-by-point responses
-
Referee: §3.2 (pair-approximation closure and mode analysis): The central claim that the slowest eigenvalue in the paramagnetic phase has strictly zero overlap with the uniform magnetization (hence exact strong inverse Mpemba for paramagnetic initial conditions) rests on the pair closure. Because the closure truncates the BBGKY hierarchy at the two-spin level, it can artificially decouple staggered and uniform sectors; the reentrant transition itself is also a mean-field artifact of the same closure. No comparison to the exact master equation on small clusters (e.g., 2×2 or 4-site lattices) or to Monte Carlo trajectories is reported, leaving the zero-overlap result unverified.
Authors: We agree that the pair approximation truncates the BBGKY hierarchy and that the reentrant transition is a mean-field artifact. Within this approximation, the linearized dynamics in the paramagnetic phase decouple the uniform and staggered sectors due to symmetry, resulting in zero overlap of the slowest (staggered) mode with the uniform magnetization. We will add to the revised manuscript a comparison of the pair-approximation spectrum with the exact master equation on a 4-site cluster to verify the decoupling at least qualitatively for small systems. A full Monte Carlo study on larger lattices is beyond the present scope. revision: partial
-
Referee: §4.3 (coordination-number scan): The disappearance of the strong Mpemba effect for coordination numbers that eliminate reentrance is presented as evidence that reentrance is necessary. However, the same scan simultaneously changes the accuracy of the pair approximation (higher z improves mean-field character), so the observed correlation could be an artifact of the closure rather than a generic feature of the model. A direct comparison at fixed z between pair-approximation and Monte Carlo relaxation spectra would be required to separate these effects.
Authors: We acknowledge that the coordination-number scan affects both reentrance and the accuracy of the pair approximation. The results show that within the pair approximation the strong Mpemba effect vanishes when reentrance is eliminated. In the revision we will add explicit discussion of this limitation and note that a fixed-z comparison with Monte Carlo would be needed to confirm the effect is not an artifact of the closure. The mechanism we identify remains valid within the approximation employed. revision: partial
-
Referee: §5 (quantitative Mpemba measures): The reported crossing times and relaxation-rate ratios that quantify the “strong” Mpemba effect are obtained from the linearized pair dynamics without error bars or sensitivity analysis to the closure parameters. Because the pair approximation is uncontrolled for dynamics, it is unclear whether the reported effect sizes survive restoration of higher-order correlations.
Authors: Since the measures are computed from deterministic linearized equations, no error bars apply. The pair closure has no free parameters. We will include a sensitivity analysis with respect to quench parameters and initial conditions in the revised §5. We will also add a statement noting that quantitative values may change with higher-order correlations, while the qualitative link to the reentrant transition persists. revision: yes
Circularity Check
Derivation follows directly from pair-approximated dynamical equations; no circularity
full rationale
The paper closes the master equation at the pair level, linearizes the resulting ODEs for the two-point correlations, and solves the eigenvalue problem for the relaxation spectrum. The claim that the slowest mode in the paramagnetic phase is purely staggered (zero overlap with uniform magnetization) is obtained by explicit diagonalization of that matrix; the strong Mpemba effect is then a direct algebraic consequence of the overlap integral being zero for paramagnetic initial conditions. No parameter is fitted to Mpemba data and then re-used as a prediction, no self-citation supplies a uniqueness theorem that forces the result, and the reentrant transition itself is an independent output of the same closure. The logical chain is therefore self-contained within the stated approximation and does not reduce to a tautology.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption The pair approximation provides a quantitatively reliable description of both statics and dynamics for the antiferromagnetic Ising model in a field.
- domain assumption The slowest relaxation mode in the paramagnetic phase is purely staggered.
Reference graph
Works this paper leans on
-
[1]
S. R. de Groot and P. Mazur,Non-equilibrium Thermody- namics, 2nd ed. (North-Holland, Amsterdam, 1962)
1962
-
[2]
Onsager, Phys
L. Onsager, Phys. Rev.37, 405 (1931)
1931
-
[3]
Onsager, Phys
L. Onsager, Phys. Rev.38, 2265 (1931)
1931
-
[4]
R. Kubo, M. Yokota, and S. Nakajima, J. Phys. Soc. Jpn. 12, 1203 (1957)
1957
-
[5]
E. M. Bertin, J.-P. Bouchaud, J.-M. Drouffe, and C. Godr` eche, J. Phys. A: Math. Gen.36, 10701 (2003)
2003
-
[6]
Bouchbinder and J
E. Bouchbinder and J. S. Langer, Soft Matter6, 3065 (2010)
2010
-
[7]
Santos, Phys
A. Santos, Phys. Rev. E109, 044149 (2024)
2024
-
[8]
Lahini, O
Y. Lahini, O. Gottesman, A. Amir, and S. M. Rubinstein, Phys. Rev. Lett.118, 085501 (2017)
2017
-
[9]
Van Vu and H
T. Van Vu and H. Hayakawa, Phys. Rev. Lett.134, 107101 (2025)
2025
-
[10]
E. B. Mpemba and D. G. Osborne, Phys. Educ.4, 172 (1969)
1969
-
[11]
G. Teza, J. Bechhoefer, A. Lasanta, O. Raz, and M. Vucelja, Phys. Rep.1164, 1 (2026)
2026
-
[12]
Lu and O
Z. Lu and O. Raz, Proc. Natl. Acad. Sci. U.S.A.114, 5083 (2017)
2017
-
[13]
Klich, O
I. Klich, O. Raz, O. Hirschberg, and M. Vucelja, Phys. Rev. X9, 021060 (2019)
2019
-
[14]
A. Kumar, R. Ch´ etrite, and J. Bechhoefer, Proc. Natl. Acad. Sci.119(2022), 10.1073/pnas.2118484119
-
[15]
G. Teza, R. Yaacoby, and O. Raz, Phys. Rev. Lett.131, 017101 (2023)
2023
-
[16]
Deg¨ unther and U
J. Deg¨ unther and U. Seifert, Europhys. Lett.139, 41002 (2022)
2022
-
[18]
Van Vu and Y
T. Van Vu and Y. Hasegawa, Phys. Rev. Res.3, 043160 (2021)
2021
-
[19]
S. K. Manikandan, Phys. Rev. Res.3, 043108 (2021)
2021
-
[20]
Meibohm, D
J. Meibohm, D. Forastiere, T. Adeleke-Larodo, and K. Proesmans, Phys. Rev. E104, L032105 (2021)
2021
-
[21]
Dieball, G
C. Dieball, G. Wellecke, and A. Godec, Phys. Rev. Res. 5, L042030 (2023)
2023
-
[22]
Ib´ a˜ nez, C
M. Ib´ a˜ nez, C. Dieball, A. Lasanta, A. Godec, and R. A. Rica, Nat. Phys.20, 135 (2024)
2024
-
[23]
Thermal relaxation asymmetry persists under inertial effects,
C. Dieball and A. Godec, “Thermal relaxation asymmetry persists under inertial effects,” (2026), arXiv:2603.18721 [cond-mat.stat-mech]
-
[24]
Baity-Jesi, E
M. Baity-Jesi, E. Calore, A. Cruz, L. A. Fernandez, J. M. Gil-Narvi´ on, A. Gordillo-Guerrero, D. I˜ niguez, A. Lasanta, A. Maiorano, E. Marinari, V. Martin-Mayor, J. Moreno-Gordo, A. M. Sudupe, D. Navarro, G. Parisi, S. Perez-Gaviro, F. Ricci-Tersenghi, J. J. Ruiz-Lorenzo, S. F. Schifano, B. Seoane, A. Taranc´ on, R. Tripiccione, and D. Yllanes, Proc. Na...
2019
-
[25]
Zhang, G
J. Zhang, G. Xia, C.-W. Wu, T. Chen, Q. Zhang, Y. Xie, W.-B. Su, W. Wu, C.-W. Qiu, P.-X. Chen, W. Li, H. Jing, and Y.-L. Zhou, Nat. Commun.16, 301 (2025)
2025
-
[26]
L. K. Joshi, J. Franke, A. Rath, F. Ares, S. Murciano, F. Kranzl, R. Blatt, P. Zoller, B. Vermersch, P. Calabrese, C. F. Roos, and M. K. Joshi, Phys. Rev. Lett.133, 010402 (2024)
2024
-
[27]
Nava and M
A. Nava and M. Fabrizio, Phys. Rev. B100, 125102 (2019). 8
2019
-
[28]
F. Ares, P. Calabrese, and S. Murciano, Nat. Rev. Phys. 7, 451 (2025)
2025
-
[29]
Summer, M
A. Summer, M. Moroder, L. P. Bettmann, X. Turkeshi, I. Marvian, and J. Goold, Phys. Rev. X16, 011065 (2026)
2026
-
[30]
Aharony Shapira, Y
S. Aharony Shapira, Y. Shapira, J. Markov, G. Teza, N. Akerman, O. Raz, and R. Ozeri, Phys. Rev. Lett.133, 010403 (2024)
2024
-
[31]
Biswas, V
A. Biswas, V. V. Prasad, O. Raz, and R. Rajesh, Phys. Rev. E102, 012906 (2020)
2020
-
[32]
Jeng, Am
M. Jeng, Am. J. Phys.74, 514 (2006)
2006
-
[33]
Lasanta, F
A. Lasanta, F. Vega Reyes, A. Prados, and A. Santos, Phys. Rev. Lett.119, 148001 (2017)
2017
-
[34]
Kumar and J
A. Kumar and J. Bechhoefer, Nature584, 64 (2020)
2020
-
[35]
Zhang and J.-X
S. Zhang and J.-X. Hou, Phys. Rev. E106, 034131 (2022)
2022
-
[36]
L. Li and J.-X. Hou, Entropy28(2026), 10.3390/e28010100
-
[37]
Yang and J.-X
Z.-Y. Yang and J.-X. Hou, Phys. Rev. E101, 052106 (2020)
2020
-
[38]
Chatterjee, S
S. Chatterjee, S. Ghosh, N. Vadakkayil, T. Paul, S. K. Singha, and S. K. Das, Phys. Rev. E110, L012103 (2024)
2024
-
[39]
Ghosh, P
S. Ghosh, P. Pathak, S. Chatterjee, and S. K. Das, Com- mun. Phys.8, 359 (2025)
2025
-
[40]
J. A. Sauer and H. N. V. Temperley, Proc. R. Soc. Lond. Ser. A. Math. Phys. Sci.176, 203 (1940)
1940
-
[41]
C. G. B. Garrett, J. Chem. Phys.19, 1154 (1951)
1951
-
[42]
J. M. Ziman, Proc. Phys. Soc.64, 1108 (1951)
1951
-
[43]
Kasteleijn, Physica22, 387 (1956)
P. Kasteleijn, Physica22, 387 (1956)
1956
-
[44]
Vives, T
E. Vives, T. Cast´ an, and A. Planes, Am. J. Phys.65, 907 (1997)
1997
-
[45]
M¨ uller-Hartmann and J
E. M¨ uller-Hartmann and J. Zittartz, Z. Phys. B27, 261 (1977)
1977
-
[46]
Binder and D
K. Binder and D. P. Landau, Phys. Rev. B21, 1941 (1980)
1941
-
[47]
M. E. Fisher, Proc. R. Soc. Lond. A. Math. Phys. Sci. 254, 66 (1960)
1960
-
[48]
M. E. Fisher, Proc. R. Soc. Lond. A. Math. Phys. Sci. 256, 502 (1960)
1960
-
[49]
Holtzman and O
R. Holtzman and O. Raz, Commun. Phys.5, 280 (2022)
2022
-
[50]
Gal and O
A. Gal and O. Raz, Phys. Rev. Lett.124, 060602 (2020)
2020
-
[51]
Aruga Katori, T
H. Aruga Katori, T. Goto, S. Ebii, and A. Ito, J. Magn. Magn. Mater.104-107, 1639 (1992)
1992
-
[52]
L. I. Glazman and A. S. Ioselevich, Z. Phys.80, 133 (1990)
1990
-
[53]
Fujii, T
H. Fujii, T. Okamoto, T. Shigeoka, and N. Iwata, Solid State Commun.53, 715 (1985)
1985
-
[54]
M. J. P. Gingras and E. S. Sørensen, Phys. Rev. B57, 10264 (1998)
1998
-
[55]
C. K. Thomas and H. G. Katzgraber, Phys. Rev. E84, 040101 (2011)
2011
-
[56]
Morita, J
T. Morita, J. Phys. A.19, 1701 (1986)
1986
-
[57]
Kitatani, S
H. Kitatani, S. Miyashita, and M. Suzuki, Phys. Lett. A.108, 45 (1985)
1985
-
[58]
Knafo, T
W. Knafo, T. Thebault, S. Raymond, P. Manuel, D. D. Khalyavin, F. Orlandi, E. Ressouche, K. Beauvois, G. Lapertot, K. Kaneko, D. Aoki, D. Braithwaite, and G. Knebel, Phys. Rev. X15, 021075 (2025)
2025
-
[59]
Unraveling anomalous relaxation effects in the thermodynamic limit,
E. Pomares, V. Mart´ ın-Mayor, A. Lasanta, and G. ´Alvarez, “Unraveling anomalous relaxation effects in the thermodynamic limit,” (2026), arXiv:2603.11326 [cond-mat.stat-mech]
-
[60]
Burley, Physica27, 768 (1961)
D. Burley, Physica27, 768 (1961)
1961
-
[61]
Shirley, Phys
T. Shirley, Phys. Lett. A42, 183 (1972)
1972
-
[62]
Hui, Phys
K. Hui, Phys. Rev. B38, 802 (1988)
1988
-
[63]
J. L. Monroe, Physica A206, 218 (1994)
1994
-
[64]
M. A. Neto and J. R. de Sousa, Phys. Rev. B70, 224436 (2004)
2004
-
[65]
M. A. Neto, R. A. dos Anjos, and J. R. de Sousa, Phys. Rev. B73, 214439 (2006)
2006
-
[66]
Katsura and M
S. Katsura and M. Takizawa, Prog. Theor. Phys.51, 82 (1974)
1974
-
[67]
D. A. Huckaby and D. A. Huckaby, Rom. J. Phys.50, 453 (2005)
2005
-
[68]
K. Blom, N. Ziethen, D. Zwicker, and A. Godec, Phys. Rev. Res.5, 013135 (2023)
2023
-
[69]
K. Blom,Pair-Correlation Effects in Many-Body Sys- tems: Towards a Complete Theoretical Description of Pair-Correlations in the Static and Kinetic Description of Many-Body Systems(Springer Nature, Cham, Switzer- land, 2023)
2023
-
[71]
K. Blom, U. Thiele, and A. Godec, Phys. Rev. E111, 024207 (2025)
2025
-
[72]
Saito and R
Y. Saito and R. Kubo, J. Stat. Phys.15, 233 (1976)
1976
-
[73]
See Supplemental Material at
-
[74]
Falasco and M
G. Falasco and M. Esposito, Rev. Mod. Phys.97, 015002 (2025)
2025
-
[75]
Blom and A
K. Blom and A. Godec, Phys. Rev. X11, 031067 (2021)
2021
-
[76]
Quenches terminating in the antiferromagnetic phase are not considered, because if the dynamics starts withs(0) = 0 it remains trapped at this value and therefore cannot relax toward the true equilibrium state withs eq ̸= 0
-
[77]
Lapolla and A
A. Lapolla and A. Godec, Phys. Rev. Lett.125, 110602 (2020)
2020
-
[78]
The free energy difference is not a metric in the strict mathematical sense, as it does not satisfy the triangle in- equality [22]
-
[79]
K. Blom, U. Thiele, and A. Godec, SciPost Phys.20, 005 (2026)
2026
-
[80]
Penrose, J
O. Penrose, J. Stat. Phys.63, 975–986 (1991)
1991
-
[81]
Ising Model and Kinetic Mean Field The- ories,
F. Ducastelle, “Ising Model and Kinetic Mean Field The- ories,” inTheory and Applications of the Cluster Varia- tion and Path Probability Methods, edited by J. L. Mor´ an- L´ opez and J. M. Sanchez (Springer US, Boston, MA, 1996) pp. 21–35. 1 Supplementary Material for: Anomalous Relaxation Through A Reentrant Phase Transition Kristian Blom1,2, Doron Beny...
1996
-
[82]
R. J. Glauber, J. Math. Phys.4, 294 (1963)
1963
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.