Extending the Mpemba effect to the underdamped realm

Technical Summary: Extending the Mpemba Effect to the Underdamped Realm

Problem Statement
The Mpemba effect (ME) is a counterintuitive phenomenon where an initially hotter system equilibrates faster than an otherwise identical, initially colder system. While extensively studied and experimentally verified in classical overdamped systems (where inertia is negligible), the existence and nature of the ME in the underdamped regime (where inertia cannot be neglected) remained unexplored. This is particularly relevant for understanding the ME in water, where molecular dynamics occur on ultrafast timescales where inertial effects are significant, and where the effective potential landscape is often modeled as a double-well structure. The central question addressed is whether the ME persists when the system transitions from the overdamped limit to the underdamped limit, and how the presence of inertia alters the conditions and mechanisms for the effect.

Methodology
The authors investigate the dynamics of a Brownian particle in a one-dimensional spatial potential $V(x)$ coupled to a thermal bath. The system is governed by the Langevin equation, and the evolution of the probability density function (PDF) is described by the Klein-Kramers (KK) equation.

1. Spectral Analysis: The core theoretical approach involves the spectral decomposition of the non-Hermitian KK operator. The relaxation dynamics are governed by the eigenvalues $\lambda_n$ and eigenfunctions $v_n$ (right) and $u_n$ (left) of the operator. The ME is characterized by the non-monotonic dependence of the overlap coefficient $a_2$ (the projection of the initial state onto the slowest relaxation mode) on the initial temperature $T_i$.
2. Perturbative Analysis (Strong Dissipation): In the limit of strong damping ($\gamma \gg m\omega_0$), the authors employ perturbation theory in powers of $1/\epsilon$ (where $\epsilon = \gamma/m\omega_0$) to derive corrections to the eigenfunctions and eigenvalues relative to the overdamped Smoluchowski limit.
3. Energy Diffusion Analysis (Weak Dissipation): In the limit of weak damping ($\gamma \ll m\omega_0$), the fast intra-well oscillations are averaged out, reducing the dynamics to a 1D diffusion process in energy space (Kramers energy-diffusion equation). The authors analyze the nodal and monotonicity properties of the eigenfunctions in this regime.
4. Numerical Simulations: The theoretical predictions are validated using numerical diagonalization of the discretized KK operator and Langevin dynamics simulations. The study utilizes an asymmetric double-well potential (canonical for overdamped ME) and a piecewise-linear double-well potential to explore specific regimes.

Key Contributions and Results

• Persistence in the Strong Dissipation Limit: The authors demonstrate perturbatively that if a robust ME exists in the overdamped limit, it persists in the underdamped regime for sufficiently large but finite damping. The underdamped overlap coefficient $a_2$ is shown to be a perturbative correction to its overdamped counterpart $a_2^{(S)}$. Consequently, the "strong Mpemba temperature" (where $a_2=0$) shifts slightly but remains close to the overdamped value. Numerical simulations confirm that the qualitative features of the ME, including the approach from opposite directions in phase space, survive the inclusion of inertia.

• Absence in Single-Well Potentials (Weak Dissipation): In the ultra-weak-damping limit, the authors prove that the ME cannot occur for smooth, confining single-well potentials with canonical (Boltzmann) initial states. By mapping the energy diffusion equation to a Sturm-Liouville problem, they show that the first excited left eigenfunction $u_2(E)$ is strictly monotonic with respect to energy. This monotonicity ensures that the overlap coefficient $a_2$ is a strictly monotonic function of the initial temperature, precluding the non-monotonicity required for the ME. This contrasts with the overdamped case, where ME can occur in single-well potentials.

• Existence in Double-Well Potentials (Weak Dissipation): The authors show that the ME can arise in the weak-damping limit if the potential is more complex, specifically a double-well potential. In this regime, the energy space becomes a branched structure (a "Y-graph") due to the barrier. The eigenfunction $u_2(E)$ is no longer a single monotonic function of energy because the same energy value can correspond to distinct branches (left and right wells) with different densities of states. This allows the overlap coefficient to exhibit non-monotonic behavior. The authors construct a specific piecewise-linear double-well potential where a strong inverse ME (faster heating) is observed in the weak-damping limit, driven by the distinct scaling of the density of states in the two wells.

• "Mpemba Phases" Diagram: Through numerical decomposition of the KK operator, the authors map out a phase diagram of Mpemba effects as a function of bath temperature ($T_b$) and dissipation rate ($\gamma$). This diagram reveals transitions between different types of effects (weak/strong, direct/inverse) as the damping is varied. Notably, varying the dissipation rate at a fixed temperature can qualitatively alter the nature of the relaxation, such as transitioning from a weak direct effect to a strong direct effect, or even switching between direct and inverse effects.

Significance and Claims
The paper claims to extend the theoretical understanding of the Mpemba effect beyond the overdamped approximation, demonstrating that the effect is robust to the inclusion of inertial dynamics.

• Theoretical Robustness: The study establishes that the ME is not an artifact of overdamped approximations but is a feature that survives into the underdamped regime, provided the potential landscape supports the necessary spectral properties.
• Mechanistic Distinction: It highlights a fundamental difference between the overdamped and weak-damping regimes regarding single-well potentials: while the overdamped regime allows for ME in single wells, the weak-damping regime requires a barrier (multi-well structure) to break the monotonicity of the energy diffusion eigenfunctions.
• Relevance to Water: The authors suggest that their findings are a step toward understanding the ME in water, given that hydrogen-bond rearrangements occur on timescales where inertia is relevant and the effective potential is often modeled as a double-well. However, they modestly note that while these simplified 1D models provide insight, the emergence of an analogous effect in the full many-body system of water is not guaranteed and remains a challenging open problem.

The work does not propose new experimental setups or specific applications for water treatment but rather provides a rigorous theoretical framework and numerical evidence for the existence and conditions of the Mpemba effect in underdamped stochastic systems.