Total: 1
We address two central open problems in the theory of anomalous Mpemba-like relaxations: their extension beyond one spatial dimension and their consistent formulation in the thermodynamic limit. Our framework is the antiferromagnetic Ising model on a square lattice under an externally applied magnetic field, which enables us to work in the presence of a phase transition. The rich phase diagram contains two control parameters: temperature and magnetic field. We demonstrate that the standard assumption of relaxation dominated by a single leading exponential is inconsistent for intensive observables exhibiting standard fluctuations. Instead, as the system size increases, a continuous spectrum of time scales emerges. Nevertheless, we make the ansatz that, in the vicinity of the phase transition, the spectral projector onto the slowest time scales can be effectively characterized in terms of an equilibrium thermodynamic quantity: the susceptibility associated with the order parameter of the metastable phase. Combined with the richness of the phase diagram, this ansatz yields qualitative and semi-quantitative predictions for optimal protocols leading to a variety of anomalous relaxation phenomena involving simultaneous variations of temperature and magnetic field. These include direct and inverse Mpemba effects, cooling-heating asymmetries, and faster heating induced by precooling. Careful Monte Carlo simulations validate our theoretical predictions. Furthermore, minimal post-optimization suffices to convert our analytically guided protocols into fully optimal ones that display anomalous relaxations in their most pronounced form.