Macroscopic Thermodynamic Framework for the Mpemba Effect
Pith reviewed 2026-05-21 11:13 UTC · model grok-4.3
The pith
A generalized cooling law with an initial-state memory term explains why hotter systems can cool faster.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Starting from linear irreversible thermodynamics, the authors derive the generalized Newton's cooling law dT/dt = −[γ0 + M Q(t)](T − Tr), where the memory term Q(t) is proportional to T(0) − Tr. The coefficient M, generated by the interplay between heat flux and structural evolution, controls anomalous behavior: M > 0 produces the Mpemba effect while M < 0 produces the inverse effect.
What carries the argument
The coefficient M that augments the bare relaxation rate through the initial-state memory term Q(t) in the generalized cooling law.
If this is right
- When M is positive, systems that begin hotter reach the reservoir temperature in less time than systems that begin cooler.
- The sign of M partitions parameter space into regions that exhibit the Mpemba effect, the inverse effect, or ordinary monotonic cooling.
- The same thermodynamic structure supplies testable criteria for the effect in any system whose relaxation involves coupled heat flow and internal structural change.
- Microscopic models of the Mpemba effect can be recovered as special cases once their effective M is computed from the underlying dynamics.
Where Pith is reading between the lines
- The same memory mechanism could govern other relaxation anomalies, such as faster equilibration after larger initial deviations in glassy or biological systems.
- Independent measurement of structural order parameters during cooling would allow direct extraction of M and a check of its thermodynamic origin.
- The framework suggests that engineering the sign of M through material properties could be used to accelerate or suppress thermal equilibration on demand.
Load-bearing premise
The derivation assumes that linear irreversible thermodynamics produces a memory term Q(t) proportional to the initial temperature difference that modifies the bare relaxation rate.
What would settle it
Measure the instantaneous cooling rate for samples prepared at several different initial temperatures and test whether the rate varies with initial temperature difference exactly as predicted by the sign and magnitude of M times Q(t).
Figures
read the original abstract
The counterintuitive Mpemba effect, wherein a hotter system cools faster, critically lacks a general macroscopic theory. Here, starting from linear irreversible thermodynamics, we formulate a generalized Newton's cooling law, $\mathrm{d}T/\mathrm{d}t = -[\gamma_0 + \mathcal{M}Q(t)](T-T_r)$, for a system at temperature $T$ relaxing in a thermal reservoir at $T_r$, where the bare relaxation rate $\gamma_0$ is modified by an initial-state memory term, $Q(t) \propto T(0)-T_r$. Arising from the interplay between heat flux and structural evolution, the coefficient $\mathcal{M}$ governs anomalous relaxation behaviors, where $\mathcal{M} > 0$ ($\mathcal{M} < 0$) induces the (inverse) Mpemba effect. This universal thermodynamic framework maps out phase diagram to provide general criteria for the Mpemba effect in complex systems, offering a macroscopic picture that bridges disparate microscopic approaches.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a macroscopic thermodynamic framework for the Mpemba effect derived from linear irreversible thermodynamics. It introduces a generalized Newton's cooling law dT/dt = −[γ₀ + ℳ Q(t)](T − T_r), where the bare rate γ₀ is modified by an initial-state memory term Q(t) ∝ T(0) − T_r. The coefficient ℳ is claimed to control anomalous relaxation, with ℳ > 0 producing the Mpemba effect and ℳ < 0 the inverse effect, and the framework is said to map a phase diagram for complex systems.
Significance. If the derivation from Onsager relations is completed and the framework is tested against data, the approach could supply a unifying macroscopic description that connects microscopic models of the Mpemba effect through the coupling of heat flux and structural evolution. The introduction of the tunable coefficient ℳ and the resulting phase diagram would then constitute a useful organizing principle.
major comments (2)
- Abstract: the generalized cooling law is stated and attributed to linear irreversible thermodynamics, yet the explicit Onsager matrix, definition of the auxiliary structural variable, and integration that isolates the factor (T(0)−T_r) in Q(t) are not exhibited. Standard LITT flux-force relations produce only dT/dt = −γ₀(T−T_r); the initial-condition dependence therefore rests on an unshown structural postulate whose validity cannot be assessed from the given text.
- Derivation of the memory term (main text): the claim that Q(t) arises from the interplay between heat flux and structural evolution requires the evolution equation for the auxiliary variable and the coupling that makes its contribution depend on the initial rather than the instantaneous temperature difference. Without these steps the sign of ℳ and the phase diagram for Mpemba versus inverse Mpemba remain formal rather than derived.
minor comments (1)
- Notation: the symbol ℳ (script M) is introduced without a clear typographic distinction from other rate coefficients; a brief remark on its physical dimensions would aid readability.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major comment below and indicate the revisions we will make to strengthen the presentation of the derivation.
read point-by-point responses
-
Referee: Abstract: the generalized cooling law is stated and attributed to linear irreversible thermodynamics, yet the explicit Onsager matrix, definition of the auxiliary structural variable, and integration that isolates the factor (T(0)−T_r) in Q(t) are not exhibited. Standard LITT flux-force relations produce only dT/dt = −γ₀(T−T_r); the initial-condition dependence therefore rests on an unshown structural postulate whose validity cannot be assessed from the given text.
Authors: We agree that the abstract summarizes the result without exhibiting the intermediate steps. The main text introduces the auxiliary structural variable and its coupling to heat flux, but does not display the explicit Onsager matrix or the integration yielding Q(t) ∝ T(0)−T_r. We will revise the manuscript by adding a dedicated derivation subsection (or appendix) that presents the Onsager relations, defines the structural variable, and shows the integration step that isolates the initial-temperature-difference factor. This will allow readers to assess the structural postulate directly. revision: yes
-
Referee: Derivation of the memory term (main text): the claim that Q(t) arises from the interplay between heat flux and structural evolution requires the evolution equation for the auxiliary variable and the coupling that makes its contribution depend on the initial rather than the instantaneous temperature difference. Without these steps the sign of ℳ and the phase diagram for Mpemba versus inverse Mpemba remain formal rather than derived.
Authors: We acknowledge that the current text describes the interplay at a conceptual level but does not supply the explicit evolution equation for the auxiliary variable or demonstrate how the coupling produces dependence on the initial (rather than instantaneous) temperature difference. In the revision we will insert the full set of coupled differential equations, show the solution for the structural variable, and derive the resulting form of Q(t) and the sign of ℳ from the thermodynamic fluxes. This will place the phase diagram on a firmer, explicitly derived footing. revision: yes
Circularity Check
Memory term Q(t) defined proportional to initial temperature difference encodes Mpemba classification by construction
specific steps
-
self definitional
[Abstract]
"we formulate a generalized Newton's cooling law, dT/dt = -[γ0 + ℳQ(t)](T-Tr), ... where the bare relaxation rate γ0 is modified by an initial-state memory term, Q(t) ∝ T(0)-Tr. Arising from the interplay between heat flux and structural evolution, the coefficient ℳ governs anomalous relaxation behaviors, where ℳ > 0 (ℳ < 0) induces the (inverse) Mpemba effect."
Q(t) is defined by construction to carry the factor (T(0)-Tr), injecting initial-condition dependence into the effective rate. The sign of the free coefficient ℳ is then stipulated to produce the Mpemba effect (or its inverse). The target phenomenon is thereby recovered by the choice of sign and the inserted proportionality rather than derived from the Onsager matrix or entropy-production structure of linear irreversible thermodynamics.
full rationale
The paper presents a generalized cooling law derived from linear irreversible thermodynamics, but the load-bearing memory term is introduced with explicit proportionality to the initial temperature difference T(0)-Tr. The coefficient M is then assigned a sign that directly determines whether the (inverse) Mpemba effect occurs. This structure makes the anomalous relaxation behavior a direct consequence of the assumed form and parameter choice rather than an independent prediction from the underlying thermodynamic relations. Standard LITT yields constant-rate relaxation; the initial-state dependence requires additional structural postulates whose explicit derivation is not exhibited, rendering the central claim partially circular.
Axiom & Free-Parameter Ledger
free parameters (1)
- M
axioms (1)
- domain assumption Linear irreversible thermodynamics applies and yields a memory term proportional to the initial temperature difference.
invented entities (1)
-
Initial-state memory term Q(t)
no independent evidence
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
generalized Newton’s cooling law … dΔT/dt = −[γ0 + M Q(t)][ΔT − I Q(t)] … M>0 (M<0) induces the (inverse) Mpemba effect
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanalpha_pin_under_high_calibration unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Onsager phenomenological equations … tight-coupling regime (|q|=1) … Δα(t)=k Q(t)
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Forward citations
Cited by 3 Pith papers
-
Nonstabilizerness Mpemba Effects
In U(1)-symmetric random circuits, initial states with lower stabilizer Rényi entropy generate nonstabilizerness faster than those with higher entropy, with the effect also depending on spatial charge structure and ex...
-
Thermodynamic Geometry of Relaxation
A Rayleigh-quotient-based thermodynamic geometry reformulates relaxation dynamics as entropic stiffness versus frictional dissipation, with linear vanishing of the slow relaxation rate near the critical temperature in...
-
Inverse engineering of cooling protocols: from normal behavior to Mpemba effects
Analytical expressions are derived for external temperature protocols that produce any prescribed internal temperature trajectory using Newtonian cooling and microscopic models, including cases with Mpemba effects.
Reference graph
Works this paper leans on
-
[1]
E. B. Mpemba and D. G. Osborne, Cool?, Phys. Educ.4, 172 (1969)
work page 1969
-
[2]
Jeng, The mpemba effect: When can hot water freeze faster than cold?, Am
M. Jeng, The mpemba effect: When can hot water freeze faster than cold?, Am. J. Phys.74, 514 (2006)
work page 2006
-
[3]
J. D. Brownridge, When does hot water freeze faster then cold water? a search for the mpemba effect, Am. J. Phys.79, 78 (2011)
work page 2011
-
[4]
G. S. Kell, The freezing of hot and cold water, Am. J. Phys.37, 564 (1969)
work page 1969
- [5]
-
[6]
Gallear, The bacon–descartes–mpemba phenomenon, Phys
R. Gallear, The bacon–descartes–mpemba phenomenon, Phys. Educ.9, 490 (1974)
work page 1974
-
[7]
D. G. Osborne, Mind on ice, Phys. Educ.14, 414 (1979)
work page 1979
-
[8]
J.Hanneken,Mpembaeffectandcoolingbyradiationtothesky, Phys. Educ.16, 7 (1981)
work page 1981
-
[9]
Auerbach, Supercooling and the mpemba effect: When hot water freezes quicker than cold, Am
D. Auerbach, Supercooling and the mpemba effect: When hot water freezes quicker than cold, Am. J. Phys.63, 882 (1995)
work page 1995
-
[10]
C. A. Knight, The mpemba effect: The freezing times of cold and hot water, Am. J. Phys.64, 524 (1996)
work page 1996
-
[11]
J. I. Katz, When hot water freezes before cold, Am. J. Phys.77, 27 (2009)
work page 2009
-
[12]
M. Vynnycky and S. L. Mitchell, Evaporative cooling and the mpemba effect, Heat Mass Transfer46, 881 (2010)
work page 2010
-
[13]
M. Balazovic and B. Tomasik, Paradox of temperature decreas- ing without unique explanation, Temperature2, 61 (2015)
work page 2015
-
[14]
A.A.Romanovsky,Whichisthecorrectanswertothempemba puzzle?, Temperature2, 63 (2015)
work page 2015
-
[15]
Walker, Hot water freezes faster than cold water
J. Walker, Hot water freezes faster than cold water. why does it do so?, Sci. Am.237, 246 (1977)
work page 1977
-
[16]
Deeson, Cooler-lower down, Phys
E. Deeson, Cooler-lower down, Phys. Educ.6, 42 (1971)
work page 1971
-
[17]
P. K. Maciejewski, Evidence of a convective instability allow- ing warm water to freeze in less time than cold water, J. Heat Transfer118, 65 (1996)
work page 1996
-
[18]
M.VynnyckyandN.Maeno,Axisymmetricnaturalconvection- driven evaporation of hot water and the mpemba effect, Int. J. Heat Mass Transf.55, 7297 (2012)
work page 2012
-
[19]
M. Vynnycky and S. Kimura, Can natural convection alone ex- plain the mpemba effect?, Int. J. Heat Mass Transfer80, 243 (2015)
work page 2015
-
[20]
S.Esposito,R.DeRisi,andL.Somma,Mpembaeffectandphase transitions in the adiabatic cooling of water before freezing, Physica A (Amsterdam)387, 757 (2008). 6
work page 2008
-
[21]
H.C.BurridgeandP.F.Linden,Questioningthempembaeffect: hot water does not cool more quickly than cold, Sci. Rep.6, 37665 (2016)
work page 2016
-
[22]
P.A.Greaney,G.Lani,G.Cicero,andJ.C.Grossman,Mpemba- like behavior in carbon nanotube resonators, Metall. Mater. Trans. A42, 3907 (2011)
work page 2011
-
[23]
Y.-H. Ahn, H. Kang, D.-Y. Koh, and H. Lee, Experimental verifications of mpemba-like behaviors of clathrate hydrates, Korean J. Chem. Eng.33, 1903 (2016)
work page 1903
-
[24]
A. Lasanta, F. V. Reyes, A. Prados, and A. Santos, When the hotter cools more quickly: Mpemba effect in granular fluids, Phys. Rev. Lett.119, 148001 (2017)
work page 2017
-
[25]
Torrente and et al., Large mpemba-like effect in a gas of inelastic rough hard spheres, Phys
A. Torrente and et al., Large mpemba-like effect in a gas of inelastic rough hard spheres, Phys. Rev. E99, 060901 (2019)
work page 2019
-
[26]
M. Baity-Jesi and et al., The mpemba effect in spin glasses is a persistentmemoryeffect,Proc.NatlAcad.Sci.USA116,15350 (2019)
work page 2019
-
[27]
A. Kumar and J. Bechhoefer, Exponentially faster cooling in a colloidal system, Nature584, 64 (2020)
work page 2020
- [28]
-
[29]
N. Giovambattista, P. J. Rossky, and H. E. Stanley, Effect of agingonthephasediagramofasupercooledwatermodel,Phys. Rev. E73, 041604 (2006)
work page 2006
-
[30]
W. Yu, S. L. Simon, and G. A. Voth, Aging and glass transition inapolymermelt: Amoleculardynamicsstudy,J.Chem.Phys. 126, 164505 (2007)
work page 2007
-
[31]
J.-B. Maddalena, M. Tissier, G. Tarjus, and P. Viot, Aging and memory effects in a kinetically constrained spin model, Phys. Rev. Lett.104, 125702 (2010)
work page 2010
- [32]
- [33]
- [34]
-
[35]
Keller and et al., Quenches across the self-organization tran- sition in multimode cavities, New J
T. Keller and et al., Quenches across the self-organization tran- sition in multimode cavities, New J. Phys.20, 025004 (2018)
work page 2018
-
[36]
A.Gijon, A.Lasanta,andE.Hernandez,Pathstowardsequilib- riuminmolecularsystems: thecaseofwater,Phys.Rev.E100, 032103 (2019)
work page 2019
-
[37]
A. Nava and M. Fabrizio, Lindblad dissipative dynamics in the presenceofphasecoexistence,Phys.Rev.B100,125102(2019)
work page 2019
-
[38]
I.Klich,O.Raz,O.Hirschberg,andM.Vucelja,Mpembaindex andanomalousrelaxation,PhysicalReviewX9,021060(2019)
work page 2019
-
[39]
Z.-Y. Yang and J.-X. Hou, Non-markovian mpemba effect in mean-field systems, Physical Review E101, 052106 (2020)
work page 2020
-
[40]
S. Zhang and J.-X. Hou, Theoretical model for the mpemba effectthroughthecanonicalfirst-orderphasetransition,Physical Review E106, 034131 (2022)
work page 2022
-
[41]
T. Van Vu and H. Hayakawa, Thermomajorization mpemba ef- fect, Phys. Rev. Lett.134, 107101 (2025)
work page 2025
-
[42]
G. Teza, J. Bechhoefer, A. Lasanta, O. Raz, and M. Vucelja, Speedups in nonequilibrium thermal relaxation: Mpemba and related effects, Physics Reports1164, 1 (2026), speedups in nonequilibriumthermalrelaxation: Mpembaandrelatedeffects
work page 2026
-
[43]
Prigogine,Étude thermodynamique des phénomènes ir- réversibles(Desoer, Liège, 1947)
I. Prigogine,Étude thermodynamique des phénomènes ir- réversibles(Desoer, Liège, 1947)
work page 1947
-
[44]
S.R.deGrootandP.Mazur,Non-EquilibriumThermodynamics (Dover Publications, New York, 1984)
work page 1984
-
[45]
Van den Broeck, Thermodynamic efficiency at maximum power, Phys
C. Van den Broeck, Thermodynamic efficiency at maximum power, Phys. Rev. Lett.95, 190602 (2005)
work page 2005
-
[46]
S. Sheng and Z. C. Tu, Weighted reciprocal of temperature, weightedthermalflux,andtheirapplicationsinfinite-timether- modynamics, Phys. Rev. E89, 012129 (2014)
work page 2014
-
[47]
Y. Izumida and K. Okuda, Work output and efficiency at maxi- mum power of linear irreversible heat engines operating with a finite-sized heat source, Phys. Rev. Lett.112, 180603 (2014)
work page 2014
-
[48]
H. Yuan, Y. H. Ma, and C. P. Sun, Optimizing thermody- namic cycles with two finite-sized reservoirs, Phys. Rev. E105, L022101 (2022)
work page 2022
-
[49]
M. J. Ondrechen, B. Andresen, M. Mozurkewich, and R. S. Berry, Maximum work from a finite reservoir by sequential carnot cycles, Am. J. Phys.49, 681 (1981)
work page 1981
-
[50]
J. G. Richens, A. M. Alhambra, and L. Masanes, Finite-bath corrections to the second law of thermodynamics, Phys. Rev. E 97, 062132 (2018)
work page 2018
-
[51]
A.M.Timpanaro,J.P.Santos,andG.T.Landi,Landauer’sprin- ciple at zero temperature, Phys. Rev. Lett.124, 240601 (2020)
work page 2020
-
[52]
Y.H.Ma,Effectoffinite-sizeheatsource’sheatcapacityonthe efficiency of heat engine, Entropy22, 1002 (2020)
work page 2020
-
[53]
Y. H. Ma, Simple realization of the polytropic process with a finite-sized reservoir, Am. J. Phys.91, 555 (2023)
work page 2023
-
[54]
I. U. Vakarelski, N. A. Patankar, J. O. Marston, D. Y. C. Chan, and S. T. Thoroddsen, Stabilization of leidenfrost vapour layer bytexturedsuperhydrophobicsurfaces,Nature489,274(2012)
work page 2012
-
[55]
J. G. Leidenfrost,De Aquae Communis Nonnullis Qualitatibus Tractatus(Ovenius, Duisburg, 1756)
-
[56]
J. G. Leidenfrost and C. Wares, On the fixation of water in diverse fire, International Journal of Heat and Mass Transfer9, 1153 (1966)
work page 1966
-
[57]
M. Shirota, M. A. J. van Limbeek, C. Sun, A. Prosperetti, and D.Lohse,Dynamicleidenfrosteffect: Relevanttimeandlength scales, Phys. Rev. Lett.116, 064501 (2016). END MATTER Appendix A: The Generalized Langevin Equation—By defining the uncoupled relaxation rates𝜈𝑇 =𝐿 𝑞𝑞 /( ˜𝐶𝑇 2 r )and 𝜈𝛼 =𝐿 𝛼 𝛼𝑘 𝛼/𝑇r, and substituting the thermodynamic forces into Eq. (4)...
work page 2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.