pith. sign in

arxiv: 2603.19887 · v2 · pith:G6QXSYR4new · submitted 2026-03-20 · ❄️ cond-mat.stat-mech · physics.class-ph· quant-ph

Macroscopic Thermodynamic Framework for the Mpemba Effect

Pith reviewed 2026-05-21 11:13 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech physics.class-phquant-ph
keywords Mpemba effectgeneralized cooling lawirreversible thermodynamicsmemory effectsanomalous relaxationphase diagramnon-equilibrium thermodynamics
0
0 comments X

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.

The paper develops a macroscopic thermodynamic description of the Mpemba effect using linear irreversible thermodynamics. It starts from the standard picture of relaxation to a reservoir and adds a memory term Q(t) proportional to the initial temperature difference that modifies the cooling rate. A single coefficient M arising from the coupling of heat flow to structural changes then determines whether relaxation is faster or slower for hotter initial states. The resulting framework supplies a phase diagram and general criteria for when the effect appears in complex systems, linking many earlier microscopic accounts under one thermodynamic umbrella.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2603.19887 by Yu-Han Ma, Yun-Qian Lin, Z. C. Tu.

Figure 1
Figure 1. Figure 1: FIG. 1. (a) Schematic of Example I. A system is surrounded by [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. The occurrence of the ME. (a) The exact cooling time [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

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)
  1. 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.
  2. 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)
  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

2 responses · 0 unresolved

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
  1. 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

  2. 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

1 steps flagged

Memory term Q(t) defined proportional to initial temperature difference encodes Mpemba classification by construction

specific steps
  1. 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

1 free parameters · 1 axioms · 1 invented entities

The framework rests on the domain assumption of linear irreversible thermodynamics and introduces the memory coefficient M together with the term Q(t) whose functional form is postulated to capture initial-state dependence.

free parameters (1)
  • M
    Coefficient that sets the magnitude and sign of the memory correction to the relaxation rate; its value determines whether the Mpemba effect appears.
axioms (1)
  • domain assumption Linear irreversible thermodynamics applies and yields a memory term proportional to the initial temperature difference.
    Invoked as the starting point for deriving the generalized Newton's law.
invented entities (1)
  • Initial-state memory term Q(t) no independent evidence
    purpose: To modify the cooling rate according to the system's starting temperature.
    Postulated to encode history dependence arising from heat flux and structural evolution.

pith-pipeline@v0.9.0 · 5705 in / 1320 out tokens · 56304 ms · 2026-05-21T11:13:57.820204+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Nonstabilizerness Mpemba Effects

    quant-ph 2026-05 unverdicted novelty 7.0

    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...

  2. Thermodynamic Geometry of Relaxation

    cond-mat.stat-mech 2026-04 unverdicted novelty 6.0

    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...

  3. Inverse engineering of cooling protocols: from normal behavior to Mpemba effects

    cond-mat.stat-mech 2026-04 conditional novelty 6.0

    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

57 extracted references · 57 canonical work pages · cited by 3 Pith papers

  1. [1]

    E. B. Mpemba and D. G. Osborne, Cool?, Phys. Educ.4, 172 (1969)

  2. [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)

  3. [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)

  4. [4]

    G. S. Kell, The freezing of hot and cold water, Am. J. Phys.37, 564 (1969)

  5. [5]

    Firth, Cooler?, Phys

    I. Firth, Cooler?, Phys. Educ.6, 32 (1971)

  6. [6]

    Gallear, The bacon–descartes–mpemba phenomenon, Phys

    R. Gallear, The bacon–descartes–mpemba phenomenon, Phys. Educ.9, 490 (1974)

  7. [7]

    D. G. Osborne, Mind on ice, Phys. Educ.14, 414 (1979)

  8. [8]

    Educ.16, 7 (1981)

    J.Hanneken,Mpembaeffectandcoolingbyradiationtothesky, Phys. Educ.16, 7 (1981)

  9. [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)

  10. [10]

    C. A. Knight, The mpemba effect: The freezing times of cold and hot water, Am. J. Phys.64, 524 (1996)

  11. [11]

    J. I. Katz, When hot water freezes before cold, Am. J. Phys.77, 27 (2009)

  12. [12]

    Vynnycky and S

    M. Vynnycky and S. L. Mitchell, Evaporative cooling and the mpemba effect, Heat Mass Transfer46, 881 (2010)

  13. [13]

    Balazovic and B

    M. Balazovic and B. Tomasik, Paradox of temperature decreas- ing without unique explanation, Temperature2, 61 (2015)

  14. [14]

    A.A.Romanovsky,Whichisthecorrectanswertothempemba puzzle?, Temperature2, 63 (2015)

  15. [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)

  16. [16]

    Deeson, Cooler-lower down, Phys

    E. Deeson, Cooler-lower down, Phys. Educ.6, 42 (1971)

  17. [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)

  18. [18]

    M.VynnyckyandN.Maeno,Axisymmetricnaturalconvection- driven evaporation of hot water and the mpemba effect, Int. J. Heat Mass Transf.55, 7297 (2012)

  19. [19]

    Vynnycky and S

    M. Vynnycky and S. Kimura, Can natural convection alone ex- plain the mpemba effect?, Int. J. Heat Mass Transfer80, 243 (2015)

  20. [20]

    S.Esposito,R.DeRisi,andL.Somma,Mpembaeffectandphase transitions in the adiabatic cooling of water before freezing, Physica A (Amsterdam)387, 757 (2008). 6

  21. [21]

    Rep.6, 37665 (2016)

    H.C.BurridgeandP.F.Linden,Questioningthempembaeffect: hot water does not cool more quickly than cold, Sci. Rep.6, 37665 (2016)

  22. [22]

    P.A.Greaney,G.Lani,G.Cicero,andJ.C.Grossman,Mpemba- like behavior in carbon nanotube resonators, Metall. Mater. Trans. A42, 3907 (2011)

  23. [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)

  24. [24]

    Lasanta, F

    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)

  25. [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)

  26. [26]

    Baity-Jesi and et al., The mpemba effect in spin glasses is a persistentmemoryeffect,Proc.NatlAcad.Sci.USA116,15350 (2019)

    M. Baity-Jesi and et al., The mpemba effect in spin glasses is a persistentmemoryeffect,Proc.NatlAcad.Sci.USA116,15350 (2019)

  27. [27]

    Kumar and J

    A. Kumar and J. Bechhoefer, Exponentially faster cooling in a colloidal system, Nature584, 64 (2020)

  28. [28]

    Blatt, P

    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,Observingthequantummpembaeffectinquantum simulations, Phys. Rev. Lett.133, 010402 (2024)

  29. [29]

    Giovambattista, P

    N. Giovambattista, P. J. Rossky, and H. E. Stanley, Effect of agingonthephasediagramofasupercooledwatermodel,Phys. Rev. E73, 041604 (2006)

  30. [30]

    W. Yu, S. L. Simon, and G. A. Voth, Aging and glass transition inapolymermelt: Amoleculardynamicsstudy,J.Chem.Phys. 126, 164505 (2007)

  31. [31]

    Maddalena, M

    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)

  32. [32]

    Zhang, Y

    X. Zhang, Y. Huang, Z. Ma, Y. Zhou, J. Zhou, W. Zheng, Q. Jiang, and C. Q. Sun, Hydrogen-bond memory and water- skin supersolidity resolving the mpemba paradox, Phys. Chem. Chem. Phys.16, 22995 (2014)

  33. [33]

    Jin and W

    J. Jin and W. A. Goddard, Mechanisms underlying the mpemba effect in water from molecular dynamics simulations, J. Phys. Chem. C119, 2622 (2015)

  34. [34]

    Lu and O

    Z. Lu and O. Raz, Nonequilibrium thermodynamics of the markovian mpemba effect and its inverse, Proc. Natl Acad. Sci. USA114, 5083 (2017)

  35. [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)

  36. [36]

    A.Gijon, A.Lasanta,andE.Hernandez,Pathstowardsequilib- riuminmolecularsystems: thecaseofwater,Phys.Rev.E100, 032103 (2019)

  37. [37]

    Nava and M

    A. Nava and M. Fabrizio, Lindblad dissipative dynamics in the presenceofphasecoexistence,Phys.Rev.B100,125102(2019)

  38. [38]

    I.Klich,O.Raz,O.Hirschberg,andM.Vucelja,Mpembaindex andanomalousrelaxation,PhysicalReviewX9,021060(2019)

  39. [39]

    Yang and J.-X

    Z.-Y. Yang and J.-X. Hou, Non-markovian mpemba effect in mean-field systems, Physical Review E101, 052106 (2020)

  40. [40]

    Zhang and J.-X

    S. Zhang and J.-X. Hou, Theoretical model for the mpemba effectthroughthecanonicalfirst-orderphasetransition,Physical Review E106, 034131 (2022)

  41. [41]

    Van Vu and H

    T. Van Vu and H. Hayakawa, Thermomajorization mpemba ef- fect, Phys. Rev. Lett.134, 107101 (2025)

  42. [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

  43. [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)

  44. [44]

    S.R.deGrootandP.Mazur,Non-EquilibriumThermodynamics (Dover Publications, New York, 1984)

  45. [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)

  46. [46]

    Sheng and Z

    S. Sheng and Z. C. Tu, Weighted reciprocal of temperature, weightedthermalflux,andtheirapplicationsinfinite-timether- modynamics, Phys. Rev. E89, 012129 (2014)

  47. [47]

    Izumida and K

    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)

  48. [48]

    H. Yuan, Y. H. Ma, and C. P. Sun, Optimizing thermody- namic cycles with two finite-sized reservoirs, Phys. Rev. E105, L022101 (2022)

  49. [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)

  50. [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)

  51. [51]

    A.M.Timpanaro,J.P.Santos,andG.T.Landi,Landauer’sprin- ciple at zero temperature, Phys. Rev. Lett.124, 240601 (2020)

  52. [52]

    Y.H.Ma,Effectoffinite-sizeheatsource’sheatcapacityonthe efficiency of heat engine, Entropy22, 1002 (2020)

  53. [53]

    Y. H. Ma, Simple realization of the polytropic process with a finite-sized reservoir, Am. J. Phys.91, 555 (2023)

  54. [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)

  55. [55]

    J. G. Leidenfrost,De Aquae Communis Nonnullis Qualitatibus Tractatus(Ovenius, Duisburg, 1756)

  56. [56]

    J. G. Leidenfrost and C. Wares, On the fixation of water in diverse fire, International Journal of Heat and Mass Transfer9, 1153 (1966)

  57. [57]

    Shirota, M

    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)...