Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2606.06362 · Journal of Statistical Mechanics Theory and Experiment

A closed system setting for quantum thermalisation in free fermions

Purvaash Panduranghan-Udhayashankar, Filiberto Ares, Pasquale Calabrese

cond-mat.stat-mechquant-phSubmitted Jun 4, 2026• Announced Jun 5, 2026• 0 citations• Revised Oct 15, 2026 (v2)
v1 · Jun 4, 2026v2 · Oct 15, 2026
arXivHTMLPDF

Abstract

We study thermalisation and the possible occurrence of the Mpemba effect in a closed quantum setting that mimics the interaction of a system with thermal reservoirs coupled only at its boundaries. Specifically, we consider a tripartite geometry in which a finite chain, initially prepared at a finite temperature, is suddenly connected on both sides to two semi-infinite chains of the same nature held at a different temperature. These outer chains act as thermal baths, while the full system evolves unitarily under a homogeneous Hamiltonian. This setup provides a simple quantum realisation of a temperature quench and closely resembles the original scenario in which the classical Mpemba effect was first observed. We focus on two paradigmatic free-fermion models, the XX chain and the transverse-field Ising chain, which respectively preserve and break the global $U(1)$ particle-number symmetry. As a probe of relaxation, we consider the Frobenius distance between the time-evolved reduced density matrix of the central subsystem and its stationary state, which is the thermal state at the bath temperature. Exploiting the free-fermionic structure of both models, the dynamics remains Gaussian and the Frobenius distance can be expressed exactly in terms of two-point correlation functions. Combining this representation with generalised hydrodynamics, we derive analytical predictions for the Frobenius distance in the hydrodynamic limit, providing a complete characterisation of the thermalisation process. Using these results, we investigate the possible occurrence of the Mpemba effect. We find that, despite the genuine non-equilibrium dynamics displayed by the system, no Mpemba effect arises in this setting. Our analysis identifies a broad class of boundary-driven thermalisation protocols in which relaxation is fully characterised analytically and exhibits no anomalous acceleration of equilibration.

Journal reference
J. Stat. Mech. (2026) 103104
DOI
10.1088/1742-5468/aea559
Comments
39 pages, 9 figures

Related papers

Quantum Mpemba effect in free-fermionic mixed states

Filiberto Ares, Vittorio Vitale, Sara Murciano

Recently, a novel probe to study symmetry breaking, known as entanglement asymmetry, has emerged and has been utilized to explore how symmetry is dynamically restored following quantum quenches. Interestingly, it has been shown that, in certain scenarios, greater initial symmetry breaking leads to faster restoration, akin to a quantum Mpemba effect. This study focuses on investigating the effect of mixed initial states and nonunitary dynamics on symmetry restoration. The mixedness of a state can arise from different sources. We consider dephasing or dissipative processes affecting initial pure states or unitary dynamics of initially thermal states. In the former case, the stationary state after the quench is independent of the initial configuration, resembling the phenomenology of the classical Mpemba effect. Investigating the XY spin chain model, through a combination of analytical calculations and numerical simulations, we identify the conditions for the occurrence of the quantum Mpemba effect. It turns out that this phenomenon still occurs in the presence of dissipation or at finite temperature, even though it will be eventually suppressed as the state becomes more mixed.

arXiv:2405.08913•Physical review. B./Physical review. B•Mar 19, 2025•43 citations

Mpemba effect and super-accelerated thermalization in the damped quantum harmonic oscillator

Stefano Longhi

The behavior of systems far from equilibrium is often complex and unpredictable, challenging and sometimes overturning the physical intuition derived from equilibrium scenarios. One striking example of this is the Mpemba effect, which implies that non-equilibrium states can sometimes relax more rapidly when they are further from equilibrium. Despite a rich historical background, the precise conditions and mechanisms behind this phenomenon remain unclear. Recently, there has been growing interest in investigating accelerated relaxation and Mpemba-like effects within quantum systems. In this work, we explore a quantum manifestation of the Mpemba effect in a simple and paradigmatic model of open quantum systems: the damped quantum harmonic oscillator, which describes the relaxation of a bosonic mode in contact with a thermal bath at finite temperature T . By means of an exact analytical analysis of the relaxation dynamics based on the method of moments in both population and coherence subspaces, we demonstrate that any initial distribution of populations with the first r moments exactly matching those of the equilibrium distribution shows a super-accelerated relaxation to equilibrium at a rate linearly increasing with r , leading to a pronounced Mpemba effect. In particular, one can find a broad class of far-from-equilibrium distributions that relax to equilibrium faster than any other initial thermal state with a temperature T ′ arbitrarily close to T . The super-accelerated relaxation effect is shown to persist even for a broad class of initial states with non-vanishing coherences, and a general criterion for the observation of super-accelerated thermalization is presented.

arXiv:2411.09589•Quantum•Mar 26, 2025•20 citations

Thermalisation and temporal relaxation in closed quantum systems

Sam Genway

This thesis approaches questions concerning the thermalisation of subsystems of closed quantum systems, prepared in pure states of definite energy but far from equilibrium, under exact unitary evolution. Taking motivation from experiments in the field of ultracold atoms, an extensive study of relaxation to a thermal state in the Hubbard model is presented. The study of small local subsystems in Hubbard-model lattice clusters has led to some interesting findings. Explored are the effects of interactions between fermions, the initial-state energy and the energy uncertainty in the initial state and their effects on relaxation dynamics and thermalisation. The most significant finding is that while subsystem thermalisation is seen for a large range of subsystem-bath coupling strengths, the temporal form of the relaxation varies markedly from exponential decay for weak couplings with a crossover to Gaussian behaviour with increased coupling strength. This is found to hold more generally for random couplings between the subsystem and bath and for bosons as well as fermions, thus demonstrating generality. As well as being demonstrated numerically, this behaviour is derived for a generic class of bi-partite quantum systems which may be described with the use of random matrices. A Brownian motion model is employed to show the exponential to Gaussian crossover when the subsystem-bath coupling matrix takes a banded form. This result agrees well with numerical Hubbard-model results, and yields identical results at short times to those from straight-forward perturbative methods. It is demonstrated that the non-Markovian Gaussian behaviour should also be observable in the limit of macroscopic baths.

•Jun 1, 2011•1 citation

Quantum quench from a thermal tensor state: Boundary effects and generalized Gibbs ensemble

Mario Collura, Dragi Karevski

We consider a quantum quench in a noninteracting fermionic one-dimensional field theory. The system of size $L$ is initially prepared into two halves $\mathcal{L}$ $([\ensuremath{-}L/2,0])$ and $\mathcal{R}$ $([0,L/2])$, each of them thermalized at two different temperatures ${T}_{\mathcal{L}}$ and ${T}_{\mathcal{R}}$, respectively. At a given time, the two halves are joined together by a local coupling and the whole system is left to evolve unitarily. For an infinitely extended system $(L\ensuremath{\rightarrow}\ensuremath{\infty})$, we show that the time evolution of the particle and energy densities is well described via a hydrodynamic approach which allows us to evaluate the correspondent stationary currents. We show, in such a case, that the two-point correlation functions are deduced, at large times, from a simple nonequilibrium steady state. Otherwise, whenever the boundary conditions are retained (in a properly defined thermodynamic limit), any current is suppressed at large times, and the stationary state is described by a generalized Gibbs ensemble, which is diagonal and depends only on the post-quench mode occupation.

arXiv:1402.1944•Physical Review B•Jun 30, 2014•94 citations

Scaling and front dynamics in Ising quantum chains

Thierry Platini, Dragi Karevski

We study the relaxation behaviour of the quantum Ising chain, focusing our attention onto the non-equilibrium dynamics of the transverse magnetization. The initial states, from which the magnetization relaxes, are product states such as those generated by setting in contact several systems, each initially equilibrated at a given temperature. Due to the free fermionic structure of the chain, the dynamics of the transverse magnetization is easily expressed in a compact form. In the completely factorized initial state, corresponding to a situation where all the spins are thermalized independently, we obtain in the scaling limit the Green function associated to the transverse magnetization. The relaxation behaviour is also considered in the system-bath case, where part of the chain called the system is thermalized at a temperature Ts and the remaining part is at a temperature Tb. The magnetization profiles show a scaling behaviour. Moreover, in the extreme case Tb=∞ and Ts=0, it is shown that the magnetization relaxes in quantized steps in the strong transverse field region.

•Jan 1, 2005•34 citations

Role of electron–electron interaction in the Mpemba effect in quantum dots

Juliane Graf, Janine Splettstoesser, Juliette Monsel

The Mpemba effect has initially been noticed in macroscopic systems-namely that hot water can freeze faster than cold water-but recently its extension to open quantum systems has attracted significant attention. This phenomenon can be explained in the context of nonequilibrium thermodynamics of Markovian systems, relying on the amplitudes of different decay modes of the system dynamics. Here, we study the Mpemba effect in a single-level quantum dot coupled to a thermal bath, highlighting the role of the sign and magnitude of the electron-electron interaction in the occurrence of the Mpemba effect. We gain physical insights into the decay modes from a dissipative symmetry of this system called fermionic duality. Based on this analysis of the relaxation to equilibrium of the dot, we derive criteria for the occurrence of the Mpemba effect using two thermodynamically relevant measures of the distance to equilibrium, the nonequilibrium free energy and the dot energy. We furthermore compare this effect to a possible exponential speedup of the relaxation. Finally, we propose experimentally relevant schemes for the state preparation and explore different ways of observing the Mpemba effect in quantum dots in experiments.

arXiv:2412.18456•Journal of Physics Condensed Matter•Mar 21, 2025•9 citations
We study thermalisation and the possible occurrence of the Mpemba effect in a closed quantum setting that mimics the interaction of a system with thermal reservoirs coupled only at its boundaries. Specifically, we consider a tripartite geometry in which a finite chain, initially prepared at a finite temperature, is suddenly connected on both sides to two semi-infinite chains of the same nature held at a different temperature. These outer chains act as thermal baths, while the full system evolves unitarily under a homogeneous Hamiltonian. This setup provides a simple quantum realisation of a temperature quench and closely resembles the original scenario in which the classical Mpemba effect was first observed. We focus on two paradigmatic free-fermion models, the XX chain and the transverse-field Ising chain, which respectively preserve and break the global U(1)U(1)U(1) particle-number symmetry. As a probe of relaxation, we consider the Frobenius distance between the time-evolved reduced density matrix of the central subsystem and its stationary state, which is the thermal state at the bath temperature. Exploiting the free-fermionic structure of both models, the dynamics remains Gaussian and the Frobenius distance can be expressed exactly in terms of two-point correlation functions. Combining this representation with generalised hydrodynamics, we derive analytical predictions for the Frobenius distance in the hydrodynamic limit, providing a complete characterisation of the thermalisation process. Using these results, we investigate the possible occurrence of the Mpemba effect. We find that, despite the genuine non-equilibrium dynamics displayed by the system, no Mpemba effect arises in this setting. Our analysis identifies a broad class of boundary-driven thermalisation protocols in which relaxation is fully characterised analytically and exhibits no anomalous acceleration of equilibration.
Recently, a novel probe to study symmetry breaking, known as entanglement asymmetry, has emerged and has been utilized to explore how symmetry is dynamically restored following quantum quenches. Interestingly, it has been shown that, in certain scenarios, greater initial symmetry breaking leads to faster restoration, akin to a quantum Mpemba effect. This study focuses on investigating the effect of mixed initial states and nonunitary dynamics on symmetry restoration. The mixedness of a state can arise from different sources. We consider dephasing or dissipative processes affecting initial pure states or unitary dynamics of initially thermal states. In the former case, the stationary state after the quench is independent of the initial configuration, resembling the phenomenology of the classical Mpemba effect. Investigating the XY spin chain model, through a combination of analytical calculations and numerical simulations, we identify the conditions for the occurrence of the quantum Mpemba effect. It turns out that this phenomenon still occurs in the presence of dissipation or at finite temperature, even though it will be eventually suppressed as the state becomes more mixed.
The behavior of systems far from equilibrium is often complex and unpredictable, challenging and sometimes overturning the physical intuition derived from equilibrium scenarios. One striking example of this is the Mpemba effect, which implies that non-equilibrium states can sometimes relax more rapidly when they are further from equilibrium. Despite a rich historical background, the precise conditions and mechanisms behind this phenomenon remain unclear. Recently, there has been growing interest in investigating accelerated relaxation and Mpemba-like effects within quantum systems. In this work, we explore a quantum manifestation of the Mpemba effect in a simple and paradigmatic model of open quantum systems: the damped quantum harmonic oscillator, which describes the relaxation of a bosonic mode in contact with a thermal bath at finite temperature T . By means of an exact analytical analysis of the relaxation dynamics based on the method of moments in both population and coherence subspaces, we demonstrate that any initial distribution of populations with the first r moments exactly matching those of the equilibrium distribution shows a super-accelerated relaxation to equilibrium at a rate linearly increasing with r , leading to a pronounced Mpemba effect. In particular, one can find a broad class of far-from-equilibrium distributions that relax to equilibrium faster than any other initial thermal state with a temperature T ′ arbitrarily close to T . The super-accelerated relaxation effect is shown to persist even for a broad class of initial states with non-vanishing coherences, and a general criterion for the observation of super-accelerated thermalization is presented.
This thesis approaches questions concerning the thermalisation of subsystems of closed quantum systems, prepared in pure states of definite energy but far from equilibrium, under exact unitary evolution. Taking motivation from experiments in the field of ultracold atoms, an extensive study of relaxation to a thermal state in the Hubbard model is presented. The study of small local subsystems in Hubbard-model lattice clusters has led to some interesting findings. Explored are the effects of interactions between fermions, the initial-state energy and the energy uncertainty in the initial state and their effects on relaxation dynamics and thermalisation. The most significant finding is that while subsystem thermalisation is seen for a large range of subsystem-bath coupling strengths, the temporal form of the relaxation varies markedly from exponential decay for weak couplings with a crossover to Gaussian behaviour with increased coupling strength. This is found to hold more generally for random couplings between the subsystem and bath and for bosons as well as fermions, thus demonstrating generality. As well as being demonstrated numerically, this behaviour is derived for a generic class of bi-partite quantum systems which may be described with the use of random matrices. A Brownian motion model is employed to show the exponential to Gaussian crossover when the subsystem-bath coupling matrix takes a banded form. This result agrees well with numerical Hubbard-model results, and yields identical results at short times to those from straight-forward perturbative methods. It is demonstrated that the non-Markovian Gaussian behaviour should also be observable in the limit of macroscopic baths.
We consider a quantum quench in a noninteracting fermionic one-dimensional field theory. The system of size LLL is initially prepared into two halves L\mathcal{L}L $([\ensuremath{-}L/2,0])$ and R\mathcal{R}R ([0,L/2])([0,L/2])([0,L/2]), each of them thermalized at two different temperatures TL{T}_{\mathcal{L}}TL​ and TR{T}_{\mathcal{R}}TR​, respectively. At a given time, the two halves are joined together by a local coupling and the whole system is left to evolve unitarily. For an infinitely extended system $(L\ensuremath{\rightarrow}\ensuremath{\infty})$, we show that the time evolution of the particle and energy densities is well described via a hydrodynamic approach which allows us to evaluate the correspondent stationary currents. We show, in such a case, that the two-point correlation functions are deduced, at large times, from a simple nonequilibrium steady state. Otherwise, whenever the boundary conditions are retained (in a properly defined thermodynamic limit), any current is suppressed at large times, and the stationary state is described by a generalized Gibbs ensemble, which is diagonal and depends only on the post-quench mode occupation.
We study the relaxation behaviour of the quantum Ising chain, focusing our attention onto the non-equilibrium dynamics of the transverse magnetization. The initial states, from which the magnetization relaxes, are product states such as those generated by setting in contact several systems, each initially equilibrated at a given temperature. Due to the free fermionic structure of the chain, the dynamics of the transverse magnetization is easily expressed in a compact form. In the completely factorized initial state, corresponding to a situation where all the spins are thermalized independently, we obtain in the scaling limit the Green function associated to the transverse magnetization. The relaxation behaviour is also considered in the system-bath case, where part of the chain called the system is thermalized at a temperature Ts and the remaining part is at a temperature Tb. The magnetization profiles show a scaling behaviour. Moreover, in the extreme case Tb=∞ and Ts=0, it is shown that the magnetization relaxes in quantized steps in the strong transverse field region.
The Mpemba effect has initially been noticed in macroscopic systems-namely that hot water can freeze faster than cold water-but recently its extension to open quantum systems has attracted significant attention. This phenomenon can be explained in the context of nonequilibrium thermodynamics of Markovian systems, relying on the amplitudes of different decay modes of the system dynamics. Here, we study the Mpemba effect in a single-level quantum dot coupled to a thermal bath, highlighting the role of the sign and magnitude of the electron-electron interaction in the occurrence of the Mpemba effect. We gain physical insights into the decay modes from a dissipative symmetry of this system called fermionic duality. Based on this analysis of the relaxation to equilibrium of the dot, we derive criteria for the occurrence of the Mpemba effect using two thermodynamically relevant measures of the distance to equilibrium, the nonequilibrium free energy and the dot energy. We furthermore compare this effect to a possible exponential speedup of the relaxation. Finally, we propose experimentally relevant schemes for the state preparation and explore different ways of observing the Mpemba effect in quantum dots in experiments.