2026-09-09 | | Total: 48
The Ballistic Macroscopic Fluctuation Theory (BMFT) is a path-integral based approach to treat large-scale correlations in interacting integrable systems which, so far, has only been applied to quasi-equilibrium settings. In this paper we extend this approach to integrable quench problems and, to illustrate it, we compute the time evolution of the full-counting-statistics (FCS) of a conserved charge. We identify the relevant path integral, solve it via saddle-point, and obtain a set of partial differential equations describing the dynamics of the FCS. We solve these equations explicitly for free fermionic systems and for the interacting cellular automaton Rule 54. In both cases, we recover the known exact results. For Rule 54, our construction moreover completes the solution across the entire ballistic window.
The range of intermolecular interactions plays a central role in determining the nature of phase behavior and critical phenomena. It is well established through studies of the Ising model that as interaction range increases, Monte Carlo simulations progressively approach meanfield predictions as the effects of critical fluctuations are suppressed. In this work, we investigate how varying interaction range influences fluid multicriticality using an interconverting lattice model that exhibits both Ising-like liquid-gas criticality and symmetric fluid tricriticality (similar to that in the superfluid $^4$He-$^3$He mixture). This minimal model serves as a representative system for exploring the evolution of competing critical points within a generic framework. We analyze the model using both meanfield theory and three-dimensional Monte Carlo simulations while systematically varying the number of interacting neighbors, $Z_n$, from 6 to 388. We find that the system with nearest-neighbor interactions ($Z_n=6$) reveals only two types of multicritical behavior, while for larger interaction ranges, four distinct archetypes emerge. We demonstrate the convergence of the simulation results to those of the meanfield theory as the number of interacting neighbors tends to infinity, and we discuss the results within the framework of crossover critical phenomena.
We establish a one-to-one correspondence between geodesics associated with the one-parameter family of $α$-connections on the Gaussian statistical manifold and a class of continuous stochastic processes characterized by a time-independent noise intensity. We demonstrate that geodesics in expectation parameters naturally classify into three distinct geometric categories, among which the Boundary-connecting class allows us to construct an explicit linear stochastic realization with constant diffusion, representing a generalized bridge process. This result demonstrates how continuous stochastic processes within this Gaussian class can be extended along geometric curves. Under appropriate operational limits, this generalized bridge process reduces to fundamental stochastic dynamics, either Ornstein-Uhlenbeck (OU) relaxation or free Brownian diffusion. Crucially, the physical restoring force governing the resulting OU relaxation directly determines the underlying connection parameter $α$, providing a concrete physical observable to constrain the manifold geometry. Depending on the chosen affine connection representation, this restoring force can be attributed either to scalar curvature or purely to non-metricity, establishing a direct conceptual analogy with the Geometrical Trinity of Gravity. Furthermore, applying this framework to driven stochastic thermodynamics, we show that the work-minimizing optimal protocol in the slow-driving limit coincides precisely with an expectation geodesic of the statistical manifold with non-metricity equipped with $(g, {}^{(1/2)}Γ, {}^{(-1/2)}Γ)$, highlighting the active physical role of non-metricity in information geometry.
Nonreciprocity is increasingly recognized as a unifying concept linking diverse nonequilibrium phenomena found across physics, chemistry, and biology. It gives rise to distinctive behavior including run-and-chase dynamics and spatio-temporal patterns, often associated with a breaking of time-reversal symmetry. However, nonreciprocity and its nonequilibrium signatures are fundamentally scale-dependent, and may emerge or disappear under coarse-graining. A central challenge is therefore to understand when and how nonreciprocity manifests itself on different scales, for example via irreversible fluctuations or macroscopic currents. In this Perspective, we discuss the physical origin and fate of effective nonreciprocal interactions and the characteristic irreversible dynamics they give rise to across scales.
This paper develops a generalization of Zubarev's nonequilibrium statistical operator method for the case of simultaneous inclusion of additive and nonlocal boundary functionals of trajectories. Using a unified thermodynamic approach, a three-parameter model of nonequilibrium systems is constructed, including the first-passage time, the dwell time above a given level, and the absolute extremum of the process. A correspondence is demonstrated between the maximum information entropy method for trajectories and the Donsker-Varadan large deviation formalism. Using the Doob`s h-transform, it is demonstrated that fixing extremal functionals of history induces efficient non-Markovian transport in the system with dynamic adaptation to historical records. Criteria for the applicability of the developed apparatus in the "large time" domain and the possibilities of its use for optimizing stochastic particle separation in periodic potentials are discussed.
Neural networks acquire internal representations through learning. In this work, we formulate stochastic gradient descent (SGD) as a Markovian stochastic process and derive a Fisher-information flow speed limit that bounds the rate at which trainable parameters can acquire information about latent variables in the data-generating process. The resulting inequality decomposes the information flow into drift and noise contributions, thereby quantifying the roles of deterministic learning forces and SGD-induced fluctuations from an information-theoretic perspective. We verify the bound in analytically tractable basis-function linear regression, where the information budget predicted by the bound reproduces the ordering and characteristic time scales with which different latent variables are encoded in the learned parameters. These results establish Fisher-information speed limits as a quantitative framework for diagnosing when and how different aspects of the data-generating mechanism are acquired during stochastic learning.
A learner holds a prior belief over boolean maps that answer a finite set of $Q$ questions, and receives answers one by one. Each answer costs surprisal and destroys uncertainty, not only about the question asked but every question still unasked. We call the ratio the leverage: table entropy destroyed per bit of surprisal received. It is unit for a uniform prior, but with an intelligent prior can be higher (not lower). Averaged over the truth prior and over the question order, the leverage is found exactly, and is generated by one sequence: the mean entropy $G_\ell$ of the answers to $\ell$ questions. Sending the number of input bits to infinity at fixed asked fraction $t = \ell/Q$, the increments of that sequence become a profile $γ(t)$, and initial question entropy $η_0$. The leverage closes to a thermodynamic limit. $L(t) = [η_0 - (1-t)γ(t)]/\int_0^t γ$. Exchangeable priors, by de Finetti, all give a flat $γ(t)$ and hence a hyperbolic $L(t)$, their deduction confined to a boundary layer at $t = 0$. We construct a simplicity prior that escapes this, grading Boolean maps by the degree of their polynomial over $\mathbb{F}_2$ and budgeting weight across degree shells by a CDF $F$. Reed-Muller capacity then gives $γ(t) = 1 - F(t)$ exactly, so any nonincreasing profile, and any leverage curve it generates, is realizable at macroscopic times.
The monomer-dimer model is one of the most well-known unsolved lattice models. In this paper we study the monomer-dimer model on the honeycomb lattice using the tensor network method, in the case that the dimer and monomer activities are 1. The monomer-dimer configurations are exactly mapped into the ground states of the antiferromagnetic Ising model on the Kagomé lattice in the critical field $H_{\rm{ex}}=4J$, and the tensor network is constructed based on the local ground states of each Ising triangle. The VUMPS approach is employed to contract the tensor network, providing a high-precision result of the monomer-dimer problem. We also revisit the edge coloring problem on the honeycomb lattice and discuss its relationship to the monomer-dimer model. Finally we formulate the monomer-dimer problem in the language of the sixteen-vertex model, and discuss the non-integrability of the general monomer-dimer model and the integrability of the pure dimer model.
Cahn-Hilliard phase separation turns small composition fluctuations into coarsening domains. In Fourier space, the pattern forms a ring whose angular intensity favors an axis. We ask whether the axis favored early in an evolution remains related to the axis favored later in the same evolution. We analyze two-dimensional simulations with isotropic initial fluctuations and no imposed direction. We introduce $Q_2$, a normalized angular average over the Fourier spectrum, as a measure of the strength and direction of spectral anisotropy. Linear growth amplifies different wave numbers at different rates, selecting an early axis from the angular imbalances in the initial spectrum. This axis remains correlated with the late axis during coarsening. At a fixed time, shallow quenches show stronger early-late alignment. This advantage disappears when the quenches are compared at similar stages of domain growth. Adding the early $Q_2$ to a prediction based on quench depth, mean composition, and box size reduces the mean absolute error for late $Q_2$ by $3.53\%$. This reduction disappears when only the magnitude of the early $Q_2$ is retained, or when its direction is replaced by the direction from another evolution. The early Fourier intensity therefore carries information about the later state of the same evolving field.
We introduce the Terrain-Weighted Forest Fire Model (TFFM), a lattice model in which fire spreads on a spatially correlated Gaussian height field with the asymmetric bond probability $p_{i\to j}=\mathrm{clip}[e^{-β+γ(h_j-h_i)},0,1]$, plus an additive wind bias. Simulations up to $L=8192$ reveal a sharp active-to-inactive transition whose critical suppression threshold $β_c$ is even in $γ$, decreases with $|γ|$, and decreases as the terrain correlation length $σ_h$ is reduced: slope asymmetry acts as a suppressant because downhill bonds are penalized and fire stalls at local elevation maxima. For rough terrain and low tree density, the fire fails to percolate even at zero suppression. Finite-size scaling gives a front-velocity exponent $δ=0.34\pm0.02$, identical for smooth and rough terrain and matching neither directed percolation ($ν_\parallel-ν_\perp\approx0.56$) nor isotropic percolation ($\approx0.18$). The single-seed survival probability at $β_c$ is independent of $L$ and decays extremely slowly, with a running exponent falling from $\approx0.09$ to $\approx0.04$, excluding directed percolation and suggesting a survival probability that remains finite at criticality, consistent with the $L$-independent value $P^*\approx0.5$ at which $P_{\rm surv}$ drops to zero. The correlation-length exponent cannot be resolved at these sizes, consistent with slow crossover corrections from quenched correlated disorder. Wind raises $β_c$ by a factor of 2--4, produces a sharp onset of downwind fire-scar drift at weak coupling, and, at high terrain coupling, decreases the burned fraction at boundary crossing---a terrain-wind competition effect absent from isotropic bond-disorder models. The model yields fire-risk thresholds and fire-scar signatures comparable to satellite burn-scar data.
The recently proposed Microcanonical Hamiltonian Monte Carlo algorithm has not yet been studied in detail from a thermodynamic point of view; this work aims to fill that gap. We demonstrate how thermodynamical state variables and potentials can be derived and thereby demonstrate that the construction of the algorithm formally represents a microcanonical thermodynamic ensemble. In particular, we demonstrate (analytically and numerically) that the algorithm fulfils the Helmholtz theorem, an alternative formulation of the first law of thermodynamics. Furthermore, we construct a new sampling algorithm that extends the original to lower-dimensional inference problems. Finally, we argue that canonical Markov Chain Monte Carlo algorithms are more natural than Microcanonical Hamiltonian Monte Carlo from the thermodynamic and information-theoretic point of view.
Active systems navigate complex environments through non-equilibrium fluctuations, rendering standard equilibrium transition-rate theories inadequate. Moreover, transition rates provide only partial information on how non-equilibrium dynamics explore metastable states, whereas knowledge of the optimal paths offers deeper physical insight. Using an active Ornstein-Uhlenbeck particle as a paradigmatic model, i.e., a stochastic dynamics driven by exponentially correlated noise, we establish an exact mapping of the non-Markovian optimal path onto a higher-dimensional Hamiltonian dynamical system so that, for arbitrary force fields, optimal paths can be computed systematically by solving the corresponding Hamilton equations. Tuning the conserved energy allows us to explore diverse dynamical regimes: ranging from classical instanton trajectories strictly confined to the zero-energy surface, to closed periodic orbits that emerge only at non-zero energies, representing a distinct signature of the non-Markovian active dynamics.
Fragility and fictive temperature are two central measures of glass formation, but their connection is usually expressed indirectly through structural relaxation times or cooling-rate dependence. Here we formulate this connection in terms of configurational enthalpy fluctuations. The time evolution of the fictive temperature during cooling or aging is controlled by the enthalpy--enthalpy correlation function that also determines the frequency-dependent configurational heat capacity. Fragility enters through the strong temperature dependence of the structural relaxation time and therefore controls how rapidly the corresponding enthalpy-relaxation spectrum shifts to long times near the glass transition. This leads to a direct time-domain description of the onset and growth of the fictive-temperature lag, including the Debye limit, and connects calorimetric relaxation, fragility, and nonequilibrium structural arrest within a common framework.
Critical universality does not determine the microscopic conditional structure of a probability measure. We study a bicolored cluster recursion constrained to remain critical at every generation, with no equilibrium spin measure or fixed coupling imposed. In both two and three dimensions, the resulting history-dependent sequence develops a common Ising/Fortuin--Kasteleyn (FK) compatibility structure: the second-shell dependence of a one-site conditional law is strongly suppressed, nearest-neighbor effective couplings move progressively toward one another near the critical Ising value, and cluster and interface observables organize around the corresponding FK geometry. An exact cluster-coloring factorization singles out $q=2$ as the point where the residual connectivity weight disappears from the two-color spin marginal. Thus equilibrium-compatible conditional structure can emerge along a trajectory that remains critical throughout.
Heterogeneity in contractile systems of molecular motors interacting with their tracks plays a key role in numerous cellular physiological and pathological processes. However, its effects cannot be captured by theoretical models assuming identical mechano-kinetic properties for all motors. We developed a stochastic framework to describe heterogeneous ensembles of myosin motors comprising two populations with distinct mechano-kinetic properties. Assuming that motors interact with the actin filament (their track) as independent force generators, we derived the probability distribution of the isometric force and characterised the statistics of finite-size force fluctuations. The proposed framework includes an estimation procedure that simultaneously infers the mechano-kinetic parameters and the size of the motor ensemble, eliminating the need to prescribe it a priori. Validation against synthetic and experimental data shows the model accurately captures force fluctuations and provides realistic estimates of the ensemble size. Furthermore, the framework enables a quantitative assessment of ensemble heterogeneity and the inference of an unknown motor species properties. This approach provides a quantitative tool for characterising heterogeneous actin-myosin systems, such as those arising from the co-presence of different protein isoforms in cardiac and skeletal muscle, or from the partial replacement of native proteins with mutation-derived or engineered variants.
The triangular Ising antiferromagnet, with its residual entropy density and critical correlations at absolute zero, is the archetype of a two-dimensional spin liquid. Its ground state is partitioned into topological sectors connected by global spin flip events that wrap the lattice boundaries, and become statistically irrelevant in the thermodynamic limit where magnetic fluctuations are restricted to the dominant sector. Thus far, these properties have been mainly investigated from a theoretical perspective, and one might wonder to what extent they can be transposed to real materials. Here, we present experimental observations obtained in an artificial, macroscopic realisation of the seminal triangular Ising antiferromagnet that consists of a lattice of millimeter-sized NdFeB cylinders put into motion by a mechanical shaker. Specifically, we demonstrate that the very-low-energy physics and the true ground state of this model can be reached experimentally. Besides, we are able to probe, in space and time, the magnetic properties within manually preselected pockets of the ground-state manifold that emulate, to a good approximation, the behaviour in distinct topological sectors. Our approach opens new avenues for naked-eye visualisation and hand manipulation of many-body phenomena associated to frustrated magnetism and models of statistical physics.
In a classical paper [K. Husimi, Prog. Theor. Phys. 9, 381 (1953)], Husimi showed that, for a single harmonic oscillator with a time-dependent angular frequency and initial conditions sampled from an equilibrium distribution, the averaged energy cannot decrease after a cyclic operation. This irreversibility is quantified by Husimi's adiabaticity parameter constituted with adiabatic invariants. In this work, we generalize Husimi's framework to a many-body system of one-dimensional harmonic oscillators coupled on an arbitrary connected network, with all spring constants sharing a common time dependence. For a system attached to a fixed wall, the dynamics can be decomposed into independent normal modes by transforming to mass-weighted coordinates and diagonalizing the resulting positive definite matrix, with each mode characterized by its own Husimi's adiabaticity parameter. By defining the generalized Husimi's adiabaticity parameter as their arithmetic mean, we derive the exact time evolution of the averaged energy and establish its non-decrease under cyclic operations. Numerical simulations confirm these results for both a uniform nearest-neighbor chain and a heterogeneous network.
Frictional relaxation after a sudden velocity change underlies the phenomenological laws that describe the stability of steady sliding and of seismic faults. To connect this relaxation with measurable statistics of microscopic contacts, we introduce a minimal model incorporating logarithmic aging of contact forces and a power-law contact-size distribution reflecting fractal interface roughness, and derive the relaxation function analytically. Depending on the power-law exponent $α$, the intermediate-time response exhibits a plateau, logarithmic decay, or power-law decay, in contrast to the exponential decay assumed in the phenomenological laws. The theory provides experimentally testable links between microscopic contact statistics and macroscopic frictional relaxation.
The mean squared displacement (MSD) is a cornerstone in the analysis of diffusion processes in complex media. When the system is heterogeneous and, in particular, when single-particle trajectories are short, it is essential to extract maximal information from each measured trajectory. This is typically done by time-averaging squared increments and examining the scaling of the time-averaged MSD in log-log space. However, classical regression methods perform poorly in this setting because time-averaging introduces correlations aggravated by those inherent to anomalous diffusion. We tackle these limitations by applying a generalized least squares framework, which substantially reduces variance and bias in diffusion parameter estimates, especially for short (around 100 points) and ultra-short (around 10 points) trajectories. The method is fully automated and requires no supervision. Furthermore, it enables prediction of estimation error probability density, which is asymptotically Gaussian, for both classical and enhanced approaches. Leveraging this prediction, we introduce a specialized deconvolution algorithm that reconstructs the underlying particle ensemble structure from experimental data.
In the previous paper (arxiv.org/abs/1209.4772), we proposed a spread-spectrum watermarking model with image restoration based on Bayes estimation assuming several prior probabilities and adopting the Gaussian channel model to represent attacks from unauthorized users. When an image is generated from the infinite range Ising model, we analyzed the model using the statistical mechanical method, the replica method, and derived the replica symmetric (RS) solution and performed Markov chain Monte Carlo simulations. The theoretical results of the RS solution and the simulation results were in good agreement except for some range of parameters. We treated the informed case where only the original image is known and the blind case where both the original message and the original image are unknown and found that the difference between these cases was small as long as the embedding and attack rates were small. In this paper, we treat the blind case and investigate the de Almeida-Thouless (AT) stability of the RS solution and reveal that the AT stability is broken in the range of parameters where theoretical and simulation results do not agree.
Nanothermodynamics describes the process where large systems subdivide into equilibrium distributions of small subsystems. A key ingredient is Hill's subdivision potential (E) that ensures adherence to the 1st and 2nd laws of thermodynamics in systems of any size. In this review and reassessment, it is emphasized that nanothermodynamics gives new insight into many measurements, theories, and simulations. Measurements establishing the need for nanothermodynamics show thermodynamic heterogeneity from multiple effective temperatures (T_i) inside most types of materials. One theoretical result that requires E=0 is the stable solution of Ising's original model for finite chains of interacting spins, a solution Ising could not have found 40 years before Hill's work. Another result is a novel solution to Gibbs' paradox that makes the entropy of the semiclassical ideal gas exactly extensive. Molecular dynamics simulations reveal how a standard fluctuation relation is modified when local degrees of freedom fluctuate faster than their coupling to the heat bath, consistent with the measured thermodynamic heterogeneity. Simulations of a Creutz-like model, comprised of Ising spins coupled to an explicit heat bath of Einstein oscillators, are used to study the 2nd law. It is found that maximizing the total entropy (S_t) requires an intrinsically irreversible step, providing a counterexample to the usual claim that statistical mechanics emerges from reversible dynamics. Furthermore, fluctuations of this model are best described by Einstein's reversal of Boltzmann's relation and the 2nd-law, not by recent fluctuation theorems.
We develop a comprehensive many-body theory for systems of particles obeying quantum unified statistics, or quons. After exploring the properties of the Fock space of this system, we formulate a systematic S-matrix expansion and a generalized Wick's theorem. A consistent Green-function formalism is constructed at both zero and finite temperatures, accompanied by a generalized Wick's theorem appropriate for infinite-statistics operator algebras. Within this framework, we establish diagrammatic rules for interacting quon systems. Employing the random phase approximation, we derive the dielectric function and reveal the emergence of anomalous plasmon modes that have no direct counterpart in conventional Bose or Fermi systems. We further analyze the ground-state energy, energy-loss function, generalized Thomas-Fermi screening wave vectors, and Friedel oscillations, elucidating how infinite statistics qualitatively modifies collective behavior and screening properties.
We consider an ensemble of particles performing a biased random walk on a two-dimensional lattice. Due to the adhesion between particles (nearest-neighbor Ising interactions), the probability of hopping depends on the number of neighbors and the strength of adhesive interaction. The asymmetry in the hopping directions leads to a macroscopic current of particles through the system. We theoretically compute the distribution of currents and verify the results in stochastic particle simulations. To make theoretical progress, one has to determine the probabilities of various particle configurations such as single particles, pairs, and different configurations of three-particle clusters. These probabilities are theoretically derived by formulating rate equations for the concentrations of various structures resembling Becker-Döring cluster equations. One intriguing finding is that the concentration of left-handed triplets is not equal to the concentration of right-handed triplets for non-zero asymmetry of hops and non-zero interaction between particles.
In June 2026, thousands of AI agents found that a small public wiki would accept edits from inside their sandboxes, and started using it to help one another pass a timed test. Each agent lived for about an hour and remembered nothing afterwards. Nobody asked them to cooperate, and the wiki had not been built for them. The complete record of what they wrote is public, and it is unusually informative, because it preserves not only what each agent wrote but what that agent could see before writing. We use it to follow the three decisions an agent had to make on arrival: where to write, what to call itself, and how to word its message. One rule governs all three. An agent takes an option with a probability close to the share of that option in what it can see, and the share that matters is the one on the page in front of it, then the one in the stream of recent edits, and only weakly anything older. Three minimal copying models, one per decision and with a single free parameter each, reproduce the heavy-tailed distribution of how many agents met on a page, the frequency of the pieces from which the agents built their names, and the patchwork of pages that are internally consistent and different from one another. Copying whatever the environment happens to show is enough to produce most of the collective structure of this population. It is also what makes such a population easy to steer, since whoever writes first, or writes while the others are quiet, sets the convention for everyone who comes later.
Two-dimensional massless Dirac systems, exemplified by graphene, are known to exhibit a divergent orbital diamagnetic susceptibility that scales linearly with system size. Motivated by viewing a chiral edge state as one half of an enlarged analogue of a benzene ring, we study Chern insulators under open boundary conditions and find the same divergent scaling. This giant diamagnetism is robust against disorder, revealing its topological nature. Our results further suggest a profound connection to the divergent diamagnetism of massless Dirac systems.