2026-03-13 | | Total: 18
Living systems are open nonequilibrium systems that continuously exchange energy, matter, and information with their environments, leading to stochastic dynamics with memory and active fluctuations. In this study, we develop a non-Markovian theoretical framework for the entropy dynamics of living systems based on the Keldysh functional formalism and stochastic thermodynamics. The approach naturally incorporates colored environmental noise, memory-dependent dissipation, and many-body interactions, yielding generalized Langevin dynamics and non-Markovian master equations. Within this framework we derive an exact frequency-domain expression for the entropy production rate and show that violations of the fluctuation-dissipation relation provide a direct thermodynamic signature of active biological fluctuations. We further demonstrate that environmental memory enhances low-frequency fluctuations and entropy production, leading to critical slowing down near dynamical instability. These results provide a microscopic physical foundation for the entropy "bathtub" picture of living systems and connect entropy evolution with development, aging, and death in nonequilibrium dynamics.
We introduce a direct Boltzmann inversion method to infer the interaction potential in particle systems using as input particle configurations generated at an arbitrary state point of the system. Unlike iterative Boltzmann inversion, the proposed method does not require performing a new Monte Carlo simulation at each step of the iteration process. It relies instead on enforcing consistency between two independent estimates of the pair correlation function, respectively obtained from interparticle distances and from pairwise forces. As a result, the approach is computationally inexpensive and straightforward to implement. Because it relies on the sole expression of interparticle forces, our method naturally applies to any state point, including when the density is large and alternative methods may fail. Here we present the basic principles of the method and benchmark its performance on a diverse set of test potentials studied using computer simulations. Practical aspects and detailed implementation of the method are also discussed. Owing to its simplicity and generality, the method should be broadly applicable, from the construction of coarse-grained interaction potentials to the inference of effective interactions in non-equilibrium systems.
Accumulation at boundaries represents a widely observed phenomenon in active systems with implications for microbial ecology and engineering applications. To rationalize the underlying physics, we provide analytical predictions for the first-passage properties and spatial distributions of a confined active Brownian particle (ABP). We show that ABPs with absorbing and hard-wall boundary conditions are Siegmund duals, yielding a direct mapping between the propagators of the two problems. We analyze the system across low and high activity regimes -- quantifying persistent motion relative to diffusion -- and show that active motion, together with a favorable initial orientation, typically lowers the mean first-passage time relative to passive diffusion. Notably, the full time-dependent propagator between hard walls approaches a wall-accumulated stationary state given by the derivative of the splitting probability as a consequence of Siegmund duality.
At lower energies, the resonances in scattering experiments are often isolated. In quantum chaotic many-body, disordered or generically stochastic systems, the resonances overlap at larger energies. Eventually, the Ericson regime is reached in which the cross section behaves like a random function. The scattering matrix elements then follow a universal Gaussian distribution. For more than sixty years, the emergence of this robust additional universal behavior on top of the universal system stochasticity awaits a concise analytical treatment. We derive the transition to the Ericson regime in the universal Heidelberg approach and prove the universal Gaussian distribution by a proper asymptotic expansion. We also obtain explicit formulae for the moments of the distributions. We compare with microwave experiments and numerical simulations.
We reassess the concept of transition at minimum work in classical stochastic finite-time thermodynamics, when the system dynamics is modelled by a diffusion process. We show that a well-posed formulation of the optimal control problem corresponding to the minimization of the mean work done on the system during a finite-time transition necessarily requires taking into account speed limits on control protocols. This fact has major qualitative consequences. First, it permits to discriminate between optimal swift engineered equilibration and transitions at minimum work. Second, it shows that in the limit when speed limits are removed, only transitions specified by generalized Schrödinger bridges admit a consistent physical interpretation. To illustrate these points, we focus on the simplest model problem: a levitating particle in a Gaussian moving trap.
Relaxation towards equilibrium is often assumed to be slower when a system starts farther from equilibrium, but this intuition fails in the Mpemba effect. Recent advances in controllable quantum platforms have enabled the exploration of its quantum analogue, the quantum Mpemba effect (QME), yet its microscopic origin remains largely unclear. Here we provide a general framework for understanding the QME in closed quantum many-body chaotic systems by reformulating the equilibration process of local subsystems in terms of Ruelle-Pollicott (RP) resonances. We show that suppressing the initial-state overlap with the dominant RP resonant mode accelerates subsystem equilibration and thereby yields the QME. We further uncover that a novel type of strong QME can occur via complete translation-symmetry breaking of initial states. We substantiate our predictions using the prototypical kicked Ising chain and exotic yet experimentally relevant initial states inspired by number theory. These findings cast the QME in closed many-body systems into a unified framework with open-system analogues and provide experimentally accessible signatures on state-of-the-art quantum platforms.
We show how to systematically construct weak integrability breaking perturbations (WIBs) for classical integrable models on the lattice. These perturbations, which allow quasi-conserved quantities, have mostly been explored in quantum systems, where they are expected to delay the onset of thermalization and diffusive transport to timescales far exceeding those predicted by Fermi's golden rule. However, accessing such long-time dynamics in quantum models is computationally challenging. Classical integrable lattice models offer a complementary setting for probing transport and long-time dynamics under WIBs. In this work, we specialize our general framework to construct several families of WIBs for the Ishimori model, the Toda chain, and the Harmonic Oscillator Chain (HOC). Such constructions can help quantify how WIBs contribute to anomalous transport and serve as a benchmark for thermalization studies in perturbed integrable models. An important example is the Fermi-Pasta-Ulam-Tsingou (FPUT) model: Starting from the HOC, we show that the cubic nonlinearity (the alpha-FPUT interaction) is a genuine WIB perturbation. Using the integrals of motion (IoMs) of the Toda lattice, we explicitly construct corrections to the entire hierarchy of the HOC IoMs, thereby obtaining an infinite tower of quasi-conserved quantities for the alpha-FPUT chain. We further identify the corresponding adiabatic gauge potential (AGP) as a nontrivial trilocal generator in real space, and show that, more generally, any cubic, translationally invariant, momentum-conserving perturbation of the HOC admits such a generator and is therefore a WIB. Together with our transport and AGP-variance studies, our results provide a unified classical framework for weak integrability breaking and for diagnosing anomalous thermalization and transport in nearly integrable Hamiltonian lattice systems.
In this work, we investigate the emergence of topological spin textures in a ferromagnetically coupled bilayer chiral magnet by means of Monte Carlo simulations of a classical spin model including exchange interaction, Dzyaloshinskii-Moriya interaction, magnetic anisotropy, and an external magnetic field. To characterize the topology of the system, we construct a scalar chirality map in the $(K/J,h/J)$ parameter space. Our results reveal several magnetic configurations, including labyrinth structures, skyrmion lattices, ferromagnetic states, and meron-antimeron crystal phases. In particular, we show that the transition from easy-axis to easy-plane anisotropy drives a continuous transformation from skyrmion textures to bimeron-type configurations. The bilayer geometry introduces an additional stabilization mechanism, where interlayer exchange coupling correlates the topological cores in the two layers and increases the energetic cost of defect collapse. These findings provide a systematic topological texture map for bilayer chiral magnets and highlight coupled magnetic layers as a promising platform for stabilizing bimeron-type spin textures in nanoscale spintronic systems.
We study third-order transitions in the two-dimensional Ising and Potts model on regular lattices and Watts--Strogatz small-world networks. Cluster observables are used to track post-critical boundary reorganization and pre-critical cluster breakup. For the Ising model, the critical temperature $T_c$ is calibrated independently from Binder-cumulant crossings and susceptibility peaks, whereas for the Potts model on small-world networks it is identified operationally from the dominant critical peak of $\mathrm d\langle P\rangle/\mathrm dT$. The independent and dependent third-order transitions are identified from the isolated-spin peak and the post-critical structural extremum, respectively. For both lattice and small-world topologies, we find the robust ordering $T_{\mathrm{ind}}<T_c<T_{\mathrm{dep}}$. Increasing the rewiring probability shifts all three characteristic temperatures upward and enhances the visibility of the post-critical transition. The effect is especially clear in the Potts model, where perimeter-based observables are more sensitive to multistate boundary fluctuations. The systematic persistence of the characteristic temperature hierarchy across topologies and finite sizes argues against interpreting these features as incidental finite-size irregularities. Instead, our results support their interpretation as genuine third-order transitions whose structural detectability can be amplified by network topology.
We develop a systematic framework to quantify irreversibility in scalar Model A field theories with a generic non-Hermitian term driving the dynamics. Using the stochastic path-integral formalism, we perform a controlled small-noise expansion, allowing the computation of the entropy production rate (EPR) and violations of the fluctuation-dissipation theorem (FDT). We show that the local EPR is entirely determined by the anti-Hermitian part of the linearised Langevin equation. Around steady states, the non-Hermitian component produces linear corrections to FDT violations and contributes quadratically to the EPR. As an illustration of the applicability of our approach, we analyse a minimal non-Hermitian extension of the Ginzburg-Landau $ψ^4$ theory describing a non-reciprocal Ising model at coarse-grained scales, for which we obtain explicit expressions of the local EPR, showing that it localises at interfaces in non-uniform states. Our results provide a general characterisation of TRS breaking in non-Hermitian scalar field theories.
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.
Reaction-diffusion systems driven far from thermodynamic equilibrium through the injection of energy can support multiple distinct spatial patterns that persist as long-lived dynamical phases. The stability of these metastable phases is not determined by thermodynamics, but by the transition paths connecting them. At finite particle numbers, intrinsic stochasticity induces rare transitions between competing patterns, rendering continuum mean-field descriptions insufficient, while exact stochastic simulations become computationally prohibitive in spatially extended systems. Here, we develop a nonequilibrium instanton framework that enables efficient computation of transition rates between metastable patterns from a single optimal transition path and its fluctuations. Using this theoretical framework, we show that an effective entropy in path space can qualitatively alter stability at finite particle numbers by increasing the exit rates of metastable patterns. By studying models of varying complexity, this work establishes path entropy as a key organizing principle for nonequilibrium pattern formation.
We develop a microscopic theory of thermalisation for a thermometer coupled to a many-body bath beyond standard Markovian and Fermi-golden-rule assumptions. By modeling interaction matrix elements in the non-interacting basis as independent random variables, we derive a diffusion-propagator expression for the reduced dynamics and show that relaxation is controlled by the distribution of interaction-induced level broadenings. The theory predicts a thermalisation timescale set by the inverse typical broadening and yields a non-Markovian generalization of global balance. Exact-diagonalization tests for heavy-tailed Lévy couplings, an all-to-all transverse-field Ising model, and the one-dimensional Imbrie model show good agreement with these predictions.
In quantum theory, equilibrium statistical mechanics is usually formulated through the canonical ensemble, whose privileged status is tied to the Euclidean continuation of time evolution. The microcanonical ensemble, by contrast, is commonly introduced as a separate spectral construction. In this work we show that this asymmetry is representational rather than structural. We formulate the system in an extended Hilbert space in which time is promoted to an auxiliary clock degree of freedom and physical states are selected by a reparametrization-invariant constraint operator $\hat C = \hat P_T + \hat H$. The corresponding projector $δ(\hat C)$ provides a single unified object from which both canonical and microcanonical ensembles emerge as complementary projections in the clock sector. In the clock-time representation, a purely imaginary clock separation yields the Euclidean kernel and the canonical partition function. In the conjugate clock-energy representation, the same projector reduces to the spectral operator $δ(\hat H-E)$ and hence to the microcanonical density of states. The main consequence is structural: canonical and microcanonical statistics need not be introduced as independent constructions, since both are already encoded in the same constrained quantum dynamics.
Synchronization is a hallmark of collective behavior in classical nonlinear systems, yet its realization as a robust many-body phenomenon in coherent quantum systems remains largely unexplored. Here we demonstrate symmetry-protected quantum synchronization and a quantum chimera state in coherent Floquet dynamics on programmable superconducting quantum processors. By implementing stroboscopic evolution of a two-dimensional Heisenberg model on IBM heavy-hex devices, we observe that initially phase-randomized spins spontaneously self-organize into coherent lattice-wide oscillations. On 28 qubits, synchronization persists even for strongly randomized initial states and is stabilized by SU(2) symmetry, as confirmed by explicit symmetry breaking. Scaling up to 156 qubits reveals a qualitatively distinct regime. For weak initial randomness, global synchronization extends across the device. For strong randomness, the system fails to synchronize globally, yet subsets of qubits exhibit robust local phase coherence under homogeneous unitary dynamics. This coexistence of globally desynchronized and locally synchronized regions constitutes a quantum analogue of a classical chimera state. Statevector and matrix-product-state simulations reproduce both the symmetry-protected synchronization and the chimera coexistence, demonstrating that these phenomena arise from the intrinsic Floquet many-body dynamics. Our results establish symmetry-protected synchronization and quantum chimera states as experimentally accessible nonequilibrium dynamical phases in programable many-body quantum systems.
Entropic forces play a fundamental role in nanoscale phenomena, from colloidal self-assembly to biomolecular disaggregation. Here, we develop an exact analytical theory and find general scaling laws for the entropic separation of tether-mediated nanofilament bundles, revealing that a single dimensionless parameter--the ratio of the excluded-volume radius to the tether length--dictates whether filaments are pushed apart or, contrary to the usual expectation, pulled together. This unexpected regime challenges the view that entropic forces invariably promote disaggregation, instead uncovering conditions under which the bundles can remain in attractive metastable states. Brownian dynamics simulations confirm this paradoxical effect, offering predictive insights for applications in biophysics, soft matter physics, and nanotechnology.
One-dimensional self-gravitating systems admit genuine thermodynamical equilibria. For systems with strictly monotonic orbital frequency profile, the Landau and Balescu-Lenard theories predict a relaxation time scaling linearly with the number of particles, $N$, in agreement with simulations. Yet, these theories become ill-posed for degenerate frequency profiles, as is the case in the harmonic potential, where all particles share the exact same mean orbital frequency. Using an exact collision-driven 1D integrator, we investigate numerically the self-consistent relaxation of 1D harmonic self-gravitating systems. We show that harmonic systems relax on a timescale that grows quadratically with $N$. We show that systems that are only partially degenerate display the same quadratic scaling for low $N$, but transition to the linear, non-degenerate behaviour for larger $N$. The larger the fraction of degenerate orbits, the larger the value of $N$ at which this transition of dynamical regime occurs. Finally, we explore the dynamics of fully non-degenerate systems, albeit with finite radial support: we confirm that their relaxation time scales linearly with $N$, though with a substantially larger prefactor than in non-compact systems. Astrophysically, this investigation should offer some new clues on the dynamics of density cores, as in the center of dwarf galaxies.
Integrability conditions on local Hamiltonians for one-dimensional quantum systems to be free and interacting fermions are introduced. The definition of free fermion is the simultaneous satisfaction of the Yang-Baxter equation and Shastry's decorated star-triangle relation by the $R$-matrix, which is more general than the previous `free-fermion algebra' by Maassarani and more special than free fermions as in the context of exactly solvable quantum models or integrable classical two-dimensional vertex models dual to quantum spin chains. Free fermionic $R$-matrices are of the difference form and have a conjugation symmetry. These free Hamiltonians may sometimes be deformed by the conjugation operator to describe an integrable interacting system with non-relativistic $R$-matrices, as are the cases of the Hubbard model and the XY model in a longitudinal field. A further criterion is obtain on precisely when such deformations remain integrable. A practical procedure is proposed to iteratively solve the free fermionic $R$-matrices from local Hamiltonians, which can be used to construct non-relativistic $R$-matrices if the conditions are met.