The Mpemba effect in the Descartes protocol: A time-delayed Newton's law of cooling approach

This paper investigates the direct and inverse Mpemba effects within a time-delayed Newton's law of cooling framework by introducing the Descartes protocol, a three-reservoir thermal scheme that analytically characterizes the roles of delay time, waiting time, and temperature to derive exact conditions and optimal parameters for maximizing the effect.

The Gist

The Big Idea: Why Hot Water Might Win the Race to Freeze

You've probably heard the Mpemba effect: the counterintuitive idea that hot water can sometimes freeze faster than cold water. For centuries, people thought this was just a weird trick of nature specific to water. But this paper shows that it's actually a deeper rule of physics that applies to many systems, provided they have a bit of "memory."

The author, Andrés Santos, uses a new experimental setup called the Descartes Protocol to study this. Think of it as a "thermal obstacle course" designed to see which runner (the hot sample or the cold sample) reaches the finish line (equilibrium) first.

---

1. The Rules of the Game: The "Memory" of Heat

Usually, we think of cooling like a car braking: if you hit the brakes (put hot water in a freezer), it slows down immediately based on how fast it's going right now.

But this paper argues that in many real-world systems, heat has memory. It's like a heavy truck with a long suspension. When you hit the brakes, the truck doesn't stop instantly; it keeps lurching forward for a moment because of its past momentum.

In physics terms, this is called a time-delayed cooling law. The temperature of your object at this exact second depends not just on the current room temperature, but on what the temperature was a tiny moment ago (the "delay time," $\tau$).

2. The Three Protocols: Different Ways to Start the Race

The paper compares three different ways to set up the race between two samples, Sample A and Sample B. Imagine they are two runners trying to reach the "Cold Finish Line."

• The Old Way (Two-Reservoir Protocol):

  • Runner A (Hot): Starts hot, then gets dumped into the cold zone.
  • Runner B (Cold): Starts cold, gets heated up to hot, then immediately dumped back into the cold zone.
  • The Problem: This is a bit of a "double dip." Runner B is confused because they were just heated up.

• The Pontus Protocol (The Fisherman's Trick):

  • Both runners start at a "Warm" temperature.
  • Runner A gets heated up to "Hot," then dumped into the cold.
  • Runner B gets cooled down to "Cold," then dumped into the cold.
  • The Issue: It uses three different temperatures, but the setup is a bit complex.

• The Descartes Protocol (The New Star):

  • This is the focus of the paper. It's the simplest and cleanest race.
  • Runner A (The Hot Start): Starts at Hot. Waits for a while. Then gets dumped into the Cold zone.
  • Runner B (The Warm Start): Starts at a Warm temperature (somewhere between Hot and Cold). Waits for a different amount of time. Then gets dumped into the Cold zone.
  • The Twist: Both runners take a single step into the cold. The only difference is when they started waiting and how warm they were to begin with.

3. The Secret Ingredient: The "Wait Time"

The paper discovers that the key to making the Mpemba effect happen isn't just about how hot the water is. It's about timing.

Imagine two people trying to catch a bus (the cold zone).

• If the hot person (A) waits too long, they cool down naturally before the bus arrives, and they lose the advantage.
• If they wait too little, they are still too hot and can't catch up.
• The Sweet Spot: The paper finds that the race is most dramatic when the "waiting time" matches the system's "memory time" (the delay $\tau$).

It's like a drummer hitting a snare drum. If you hit the drum exactly when the echo of the previous hit is fading, you get the loudest, most powerful sound. If you hit it too early or too late, the sound is weak. The author found that the Mpemba effect is "loudest" when the waiting time equals the memory delay.

4. The Results: Who Wins?

The author did the math to find the perfect conditions:

1. The "Goldilocks" Temperature: There is a specific "Warm" starting temperature for Runner B that makes the race most interesting. It's not too hot, not too cold, but just right.
2. The Magnitude: The paper calculates exactly how much faster the hot runner can win.
3. The Surprise Comparison: Even though the Descartes protocol has an extra "knob" to turn (the starting warm temperature), it actually produces a smaller Mpemba effect than the older, messier "Two-Reservoir" protocol.
   • Analogy: It's like having a fancy, high-tech car with more gears (Descartes) vs. a simple, rugged truck (Two-Reservoir). Sometimes, the simpler truck just handles the rough terrain better. The extra control in the Descartes protocol actually limits the maximum speedup you can get.

5. What About Real Life? (Finite-Rate Quenches)

In the real world, you can't instantly drop a cup of water from boiling to freezing. It takes a second to cool down. The paper checks what happens if the "drop" isn't instant.

• The Bad News: If the cooling process is too slow, the "perfect" Mpemba effect disappears because the two runners are no longer experiencing the exact same environment at the same time.
• The Good News: If the cooling is fast enough (but not instant), the effect still happens, just a little less perfectly. It's like running a race in the rain; you still win, but your shoes get a little wet.

Summary: Why This Matters

This paper is important because it strips away the messy details of water (like evaporation or bubbles) and shows that the Mpemba effect is a fundamental property of systems with memory.

• The Takeaway: If you have a system that remembers its past (like a heavy engine, a complex chemical reaction, or even a quantum computer), you can manipulate the timing of your inputs to make a "hot" system cool down faster than a "cold" one.
• The Metaphor: It's not about being "hotter." It's about being strategically positioned in time so that when the environment changes, your internal momentum carries you to the finish line faster than the person who started closer to the goal but had no momentum.

The author calls this the Descartes Protocol because it mimics a recipe from the philosopher René Descartes, who suggested letting boiling water cool to the temperature of spring water before starting an experiment. The paper proves that Descartes was onto something, but the real magic happens when you tune the timing of that cooling perfectly.

Technical Summary

1. Problem Statement

The Mpemba effect is the counterintuitive phenomenon where a system initially at a higher temperature relaxes to equilibrium faster than an identical system starting at a lower temperature. While historically associated with water freezing, it is now recognized as a general non-equilibrium phenomenon in various systems (granular fluids, spin models, quantum systems).

The paper addresses the need for a unified, analytically tractable framework to characterize this effect, specifically focusing on memory effects in thermal relaxation. Previous studies by the author utilized a "two-reservoir protocol" involving double-step quenches. This work introduces a new, simpler scheme called the Descartes protocol to investigate both direct and inverse Mpemba effects within the framework of the time-delayed Newton's law of cooling. The goal is to disentangle the roles of the delay time, waiting time, and initial temperature conditions to determine the precise conditions and magnitude of the effect.

2. Methodology

Theoretical Framework

The study utilizes the time-delayed Newton's law of cooling:
$$\dot{T}(t) = -\lambda [T(t - \tau) - T_b(t)]$$
where:

• $T(t)$ is the system temperature.
• $T_b(t)$ is the time-dependent bath temperature.
• $\tau$ is the delay time (representing memory or finite response time of the environment).
• $\lambda$ is the heat transfer coefficient (set to 1 for time units).

The solution involves a "quasi-exponential" function, $E(t)$, which describes the relaxation dynamics. For long times, this decays exponentially with rates determined by the Lambert $W$ function.

The Descartes Protocol

The authors propose a three-reservoir thermal scheme (illustrated in Fig. 1c):

1. Sample A: Thermalized at a hot temperature ($T_{hot}$) until $t = -t_w$, then quenched to a cold temperature ($T_{cold}$).
2. Sample B: Thermalized at an intermediate "warm" temperature ($T_{warm}$) until $t = 0$, then quenched to the same cold temperature ($T_{cold}$).
3. Key Parameters:
   • $\tau$: Delay time.
   • $t_w$: Waiting time (time difference between the quenches).
   • $\omega$: Normalized warm temperature, defined as $\omega = \frac{T_{warm} - T_{cold}}{T_{hot} - T_{cold}}$.

The condition for the Direct Mpemba Effect is that Sample A (starting hotter) crosses Sample B (starting cooler) at some time $t_\times$ and subsequently remains cooler ($T_A(t) > T_B(t)$ for $t < t_\times$ and $T_A(t) < T_B(t)$ for $t > t_\times$).

The Inverse Mpemba Effect (heating scenario) is analyzed by swapping the roles of hot and cold temperatures.

Analytical Approach

The authors define a temperature difference function $\Delta(t) = \theta_A(t) - \theta_B(t)$. They derive exact conditions for the existence of the effect (bounds on $\omega$) and identify the optimal parameters that maximize the effect's magnitude. They also extend the analysis to finite-rate quenches, where the bath temperature changes exponentially over a characteristic time $\sigma$ rather than instantaneously.

3. Key Contributions

1. Introduction of the Descartes Protocol: A novel experimental design where both samples undergo a single-step quench at different times, utilizing three distinct temperatures ($T_{hot}, T_{warm}, T_{cold}$). This simplifies the parameter space compared to previous double-step protocols while introducing $\omega$ as a tunable control parameter.
2. Exact Conditions for Existence: Derivation of rigorous bounds for the normalized warm temperature $\omega$ required for the Mpemba effect to occur:
   $$e^{-\kappa_0 t_w} < \omega < E(t_w)$$
   where $\kappa_0$ is the dominant decay rate and $E(t_w)$ is the quasi-exponential function evaluated at the waiting time.
3. Optimization of the Effect: Identification of the optimal waiting time $t_w$ that maximizes the effect. The authors prove that the effect is maximized when the waiting time equals the delay time:
   $$t_w^{opt} = \tau$$
4. Quantification and Approximations: Development of accurate, compact analytical approximations for:
   • The optimal $\omega$ ($\tilde{\omega}(t_w)$) that maximizes the effect.
   • The maximal magnitude of the effect ($M_p(t_w)$).
     These approximations avoid the need for solving transcendental equations numerically.
5. Comparison with Two-Reservoir Protocol: A systematic comparison showing that despite having an extra control parameter ($\omega$), the Descartes protocol yields a smaller maximal magnitude of the Mpemba effect compared to the previously studied two-reservoir (double-step) protocol.
6. Finite-Rate Quench Analysis: Demonstration that while finite-rate quenches ($\sigma > 0$) theoretically widen the parameter window for the effect, the requirement for identical bath conditions for both samples after $t=0$ imposes a constraint that strictly forbids a genuine Mpemba effect. However, an approximate effect persists if $\sigma$ is sufficiently small.

4. Key Results

• Phase Space Structure: The existence of the Mpemba effect is confined to a specific region in the 3D parameter space $(\tau, t_w, \omega)$. The width of the allowed $\omega$ interval, $\delta\omega$, vanishes as $t_w \to 0$ or $t_w \to \infty$, peaking at $t_w = \tau$.
• Resonance-like Behavior: The magnitude of the Mpemba effect exhibits a non-monotonic dependence on the waiting time, reaching an absolute maximum when $t_w = \tau$. This suggests a "resonance" between the memory scale of the system and the experimental waiting time.
• Magnitude Comparison: The ratio of the maximal effect in the two-reservoir protocol to that in the Descartes protocol is always greater than 1. The two-reservoir protocol is more efficient at generating a strong Mpemba effect, likely due to the specific trajectory of the double-step quench.
• Inverse Effect: The inverse Mpemba effect (faster heating) follows the same mathematical structure as the cooling case, with the bounds on $\omega$ shifted to $1 - E(t_w) < \omega < 1 - e^{-\kappa_0 t_w}$.
• Finite-Rate Limit: For finite-rate quenches, strict equality of bath conditions ($T_{b,A}(t) = T_{b,B}(t)$ for $t>0$) is impossible unless $\sigma=0$. Consequently, a strict Mpemba effect is precluded, though it remains observable as an approximation for small $\sigma$.

5. Significance

• Unified Framework: The paper provides a unified, analytically solvable model for the Mpemba effect based on non-Markovian (memory) dynamics, bridging the gap between phenomenological observations and rigorous statistical mechanics.
• Experimental Guidance: By identifying $t_w = \tau$ as the optimal condition, the paper offers concrete experimental guidelines for observing the effect in systems with thermal inertia or delayed response (e.g., mesoscopic devices, feedback-controlled systems).
• Protocol Efficiency: The comparison between protocols clarifies that adding more reservoirs or steps does not necessarily maximize the effect; the specific timing and temperature sequencing are critical.
• Memory Effects: The work reinforces the idea that the Mpemba effect is a natural consequence of memory effects in relaxation dynamics, rather than requiring a breakdown of standard cooling laws. It suggests that such effects can be engineered in systems where delays are introduced deliberately.
• Future Directions: The methodology is presented as a tool for analyzing other complex thermal protocols, such as the "Pontus protocol," opening avenues for further research into how protocol structure influences non-equilibrium relaxation pathways.

In summary, this paper rigorously characterizes the Mpemba effect using a time-delayed cooling model, introducing a new experimental protocol that isolates key variables, derives exact existence conditions, and provides practical approximations for maximizing the effect, while highlighting the limitations imposed by finite-rate thermal changes.