Mpemba effect in a two-dimensional bistable potential

TL;DR

This paper introduces an exactly solvable two-dimensional model demonstrating the Mpemba effect, revealing how non-monotonic relaxation behavior depends on initial temperature and potential geometry.

Contribution

It provides an analytical two-dimensional model with explicit solutions showing the Mpemba effect without confining walls, extending previous one-dimensional models.

Findings

The slowest mode coefficient varies non-monotonically with initial temperature.

A crossing condition for the Kullback-Leibler divergence is derived.

The Mpemba effect is demonstrated in specific parameter regimes.

Abstract

We present an exactly solvable model of the Mpemba effect in an overdamped Langevin system confined in a two-dimensional radially symmetric bistable potential. The potential is constructed as a piecewise quadratic-logarithmic function that is continuous and differentiable at the matching radii, enabling an exact mapping of the corresponding Fokker-Planck operator to a Schroedinger-type eigenvalue problem. The relaxation spectrum and eigenmodes are obtained analytically in each region in terms of confluent hypergeometric functions, with eigenvalues determined from matching conditions. Focusing on isotropic equilibrium initial states at inverse temperature $\beta_{\rm ini}$ quenched to a bath at inverse temperature $\beta$, we derive explicit expressions for the mode amplitudes governing long-time relaxation. We demonstrate that the coefficient of the slowest mode exhibits non-monotonic dependence on $\beta_{\rm ini}$ and identify a sufficient crossing condition for the Kullback-Leibler divergence in terms of the two slowest modes, if the global minimum of the potential is located far away from the origin and the second minimum exists near the origin. For corresponding parameters, we demonstrate that the Mpemba effect can be realized. Our results provide a rare example of an analytically tractable two-dimensional model exhibiting anomalous relaxation without any confining walls, extending previous one-dimensional constructions with a hard wall and clarifying the role of radial geometry in nonequilibrium relaxation phenomena.