2026-02-10 | | Total: 24
Motivated by earlier numerical evidence for a percolation-like transition in space-time jamming, we present an analytic description of the transient dynamics of the deterministic traffic model elementary cellular automaton rule 184 (ECA184). By exploiting the deterministic structure of the dynamics, we reformulate the problem in terms of a height function constructed directly from the initial condition, and obtain an equilibrium statistical mechanics-like description over the lattice configurations. This formulation allows macroscopic observables in space-time, such as the total jam delay and jam relaxation time, as well as microscopic jam statistics, to be expressed in terms of geometric properties of the height function. We thereby derive the associated scaling forms and recover the critical exponents previously observed in numerical studies. We discuss the physical implications of this space-time geometric approach.
This paper addresses uphill transport (defined as a regime in which particle flow is opposite to the prescriptions of Fick's diffusion) in drift-diffusion particle transport constrained by volume exclusion. Firstly, we show that the stationary hydrodynamic limit of a multispecies, weakly asymmetric exclusion process (SHDL) naturally predicts precisely characterized uphill regimes in the space of external drivings. Then, with specific reference to systems of oppositely charged particles, we identify well-defined model hypotheses and extensions whereby the SHDL converges to the modified Poisson-Nernst-Planck model, thus bridging the gap between exclusion-based particle models and continuum descriptions commonly used in engineering. The merits and limitations of the models in describing the particle fluxes and predicting uphill transport conditions are investigated in detail with respect to the adopted approximations and simplifications. The results demonstrate the persistence of uphill transport phenomena across modeling scales, clarify the conditions under which they occur, and suggest that uphill transport may play a significant role in nanoscale electrolytes, confined ionic and iontronic devices, and membrane-based technologies.
Numerical simulations of critical lattice systems are fundamentally limited by critical slowing down, as long-range correlations are typically established through slow temporal equilibration. A physically constrained generative framework that replaces temporal relaxation with a spatial projection mechanism for critical systems is proposed. Using the two-dimensional Ising model at criticality as a benchmark, we introduce an energy-constrained kernel that synthesizes large-scale configurations from compact equilibrated seeds by enforcing Hamiltonian-level observables. The generated configurations are projected onto the nearest-neighbor energy manifold, ensuring thermodynamic consistency while retaining universal critical properties. We show that the resulting configurations reproduce scale-invariant spin correlations, Binder cumulants, and isotropic structure factors for lattice sizes exceeding 10,000, without iterative Monte Carlo equilibration. While not a strict renormalization group transformation, and motivated by renormalization ideas, the method provides a practical inverse mapping that retains universal features of criticality and enables efficient GPU-parallel generation of ultra-large critical ensembles.
We derive the exact nonequilibrium steady state of a run-and-tumble particle (RTP) in $d$ dimensions confined in an isotropic harmonic trap $V(\mathbf r)=μr^{2}/2$, with $r=\|\mathbf r\|$. Rotational invariance reduces the problem to the stationary single-coordinate marginal $p_X(x)$, from which the radial distribution $p_R(r)$ and the full joint stationary density follow by explicit integral transforms. We first focus on a generalized trapped RTP in one dimension, where post-tumble velocities are drawn from an arbitrary distribution $W(v)$. Using a Kesten-type recursion, we represent its stationary position in terms of a stick-breaking (or Dirichlet) process, yielding closed-form expressions for its distribution and its moments. Specializing $W(v)$ to the projected velocity law of an isotropic RTP, we reconstruct $p_R(r)$ and the full joint distribution of all the coordinates in $d=1,2,3$. In $d=1$ and $d=2$, the radial law simplifies to a beta distribution, while in $d=3$, we derive closed-form expressions for $p_R(r)$ and the stationary joint distribution $P(x,y,z)$, which differ from a beta distribution. In all cases, we characterize a persistence-controlled shape transition at the turning surface $r=v_0/μ$, where $v_0$ is the self-propulsion speed. We further include thermal noise characterized by a diffusion coefficient $D>0$, showing that the stationary law is a Gaussian convolution of the $D=0$ result, which regularizes turning-point singularities and controls the crossover between persistence- and diffusion-dominated regimes as $D \to 0$ and $D \to \infty$ respectively. All analytical predictions are systematically validated against numerical simulations.
Asymmetric exclusion process (TASEP) along a one-dimensional (1D) open channel sets the paradigm for 1D driven models and nonequilibrium phase transitions in open 1D models. Inspired by the phenomenologies of an open TASEP with Langmuir kinetics (Lk) and with finite resources, we study the stationary densities and phase transitions in a TASEP with Lk connected to a particle reservoir at its both ends. We calculate the stationary density profiles and the phase transitions. The resulting phase diagrams in the plane of the control parameters are significantly different from their counterparts in an open TASEP with Lk. In particular, some of the phases admissible in the open TASEP with Lk model are no longer possible. Intriguingly, our model that is closely related to a TASEP coupled with Lk on a ring with a point defect, admits more phases than the latter. Phenomenological implications of our results are discussed.
We develop a variational method for interacting surface systems with higher-form global symmetries. As a natural extension of the conventional second-quantized Hamiltonian of interacting bosons, we explicitly construct a second-quantized Hamiltonian formulated in terms of a closed surface operator $\hatφ[C_p^{}]$ charged under a $p$-form global symmetry. Applying the variational principle, we derive a functional Schrödinger equation analogous to the Gross-Pitaevskii equation in conventional bosonic systems. In the absence of external forces, the variational equation admits a uniform solution that is uniquely determined by a microscopic interaction potential $U(ψ^*ψ)$ and the chemical potential. This uniform solution describes a uniform gas of bosonic surfaces. Using the obtained energy functional, we show that low-energy fluctuations contain a gapless $p$-form field $A_p^{}$ when the $p$-form global symmetry is $\mathrm{U}(1)$, whereas the $p$-form field becomes massive for discrete symmetries, whose low-energy limit is described by a $\mathrm{BF}$-type topological field theory. As a consequence, the system exhibits abelian topological order with anyonic surface excitations. In the presence of external forces, however, solving the functional equation in full generality remains challenging. We argue, however, that the problem reduces to solving the conventional Gross-Pitaevskii equation when external forces act separately on the center-of-mass and relative motions. In addition, we present analytic solutions for topological defects as analogs of vortex and domain-wall solutions in conventional bosonic systems. Finally, as a concrete microscopic model, we study a $\mathbb{Z}_N^{}$ lattice gauge theory and apply our variational method to this system.
We provide a link between response theory, pullback measures, and optimal fingerprinting method that paves the way for a) predicting the impact of acting forcings on time-dependent systems and b) attributing observed anomalies to acting forcings when the reference state in not time-independent. We first derive formulas for linear response theory for time-dependent Markov chains and diffusions processes. We discuss existence, uniqueness, and differentiability of the pullback measure under general (not necessarily slow or periodic) perturbations of the transition kernels. An explicit Green-Kubo-type formula for the linear response is derived. We analyze in detail the case of periodic reference dynamics, where the unperturbed pullback attractor is periodic but the response is generally not. Our formulas reduce to those of classic linear response if one considers a reference autonomous state. Finally, we show that our results allow for extending the theory of optimal fingerprinting for detection and attribution of climate change (or change in any complex system) for the case of time-dependent background state and for the case where the optimal solution is sought for multiple time slices at the same time. We provide strong numerical support for the findings by applying our theory to a modified version of the Ghil-Sellers energy balance model where we include explicit time dependence in the reference state as a result of natural forcings. We verify the accuracy of response theory in predicting the impact of increases of $CO_2$ in the temperature field even when we discretize the system using Markov state modelling approach. Additionally, we consider a more complex modelling scenario where a localized aerosol forcing is also included in the system and show that the optimal fingerprinting method developed here is able to attribute the climate change signal to the acting forcings.
The Lorentz mirror model provides a clean setting to study macroscopic transport generated solely by quenched environmental randomness. We introduce a hierarchical version that admits an exact recursion for the distribution of left--right crossings, and prove normal transport: the mean conductance scales as (cross-section)/(length) for all length scales if $d\ge3$. A Gaussian approximation, supported by numerics, predicts that, in the marginal case $d=2$, this scaling acquires a logarithmic correction and that the variance-to-mean ratio of conductance converges to the universal value $2/3$ (the ``$2/3$ law'') for all $d\ge2$. We conjecture that both effects persist beyond the hierarchical setting. We finally provide numerical evidence for the $2/3$ law in the original Lorentz mirror model in $d=3$, and interpret it as a universal signature of normal transport induced by random current matching. A YouTube video discussing the background and the main results of the paper is available: https://youtu.be/G1nqKd6MiXo
Entropy is critically examined as a fundamental concept in contemporary science and informatics. Although the typical Shannon entropy provides a proper framework for describing the canonical ensemble, it fails to represent adequately the microcanonical ensemble. This discrepancy manifests additionally in its inability to support a theoretical derivation of the Second Law of thermodynamics.
The establishment of the Wilson-Fisher fixed point (WFP) for $O(n)$ spin models in $d=4-ε$ dimensions stands as a cornerstone of the renormalization group (RG) theory for critical phenomena. However, when long-range (LR) interactions, algebraically decaying as $\propto 1/r^{d+σ}$, are introduced, the fate of the short-range WFP (SR-WFP) has remained a subject of intense debate since the 1970s. We employ two complementary techniques -- the standard field-theoretic RG and a perturbative bootstrap scheme, and perform the $ε$-expansion calculations up to the two-loop level. We show that, as long as $σ<2$, the SR-WFP becomes unstable and a stable LR-WFP emerges, and, in the non-classical regime with $d/2 < σ< 2$, the critical exponents, including the anomalous dimension, are functions of $ε$, $δ=2-σ$ and $n$, which reduce to the exact results in the limiting cases $ε\to 0$, $δ\to 0$ or $n \to \infty$. Our $(4-ε)$-expansion calculations support the scenario that the threshold between the LR- and SR-WFP occurs strictly at $σ_*=2$, well consistent with the recent high-precision numerical study while different from the widely accepted Sak's criterion.
We present implementations of two physically-embedded computation-universal logical operations using a 2-bit logical unit composed of coupled quantum flux parametrons -- Josephson-junction superconducting circuits. To illustrate universality, we investigate NAND gates built from these two distinct elementary operations. On the one hand, Controlled Erasure (CE) is designed using fixed-point analysis and assumes that information must be stored in locally-metastable distributions. On the other, Erasure-Flip (EF) leverages momentum as a computational resource and significantly outperforms the metastable approach, simultaneously achieving higher fidelity and faster computational speed without incurring any additional energetic cost. Notably, the momentum degree of freedom allows the EF to achieve universality by using both nontrivial reversible and irreversible logic simultaneously in different logical subspaces. These results not only provide a practical, high-performance protocol ripe for experimental realization but also underscore the broader potential of momentum-based computing paradigms.
The application of generative modeling to many-body physics offers a promising pathway for analyzing high-dimensional state spaces of spin systems. However, unlike computer vision tasks where visual fidelity suffices, physical systems require the rigorous reproduction of higher-order statistical moments and thermodynamic quantities. While Score-Based Generative Models (SGMs) have emerged as a powerful tool, their standard formulation on Euclidean embedding space is ill-suited for continuous spin systems, where variables inherently reside on a manifold. In this work, we demonstrate that training on the Euclidean space compromises the model's ability to learn the target distribution as it prioritizes to learn the manifold constraints. We address this limitation by proposing the use of Manifold-Aware Score-Based Generative Modeling framework applied to the 64x64 2D XY model (a 4096-dimensional torus). We show that our method estimates the theoretical Boltzmann score with superior precision compared to standard diffusion models. Consequently, we successfully capture the Berezinskii-Kosterlitz Thouless (BKT) phase transition and accurately reproduce second-moment quantities, such as heat capacity without explicit feature engineering. Furthermore, we demonstrate zero-shot generalization to unseen lattice sizes, accurately recovering the physics of variable system scales without retraining. Since this approach bypasses domain-specific feature engineering, it remains intrinsically generalizable to other continuous spin systems.
We investigate the six-state clock universality of the Ising model on the kagome lattice, considering antiferromagnetic nearest-neighbor (NN) and ferromagnetic next-nearest-neighbor (NNN) interactions. Our comprehensive study employs three approaches: the level-spectroscopy method, Monte Carlo simulations, and a machine-learning phase classification technique. In this system, we observe two Berezinskii-Kosterlitz-Thouless (BKT) transitions. We present a phase diagram consisting of three phases: the low-temperature ordered phase with sublattice magnetizations, the intermediate BKT phase, and the high-temperature disordered phase, as a function of the ratio of the NNN interaction to the NN interaction. We verify the six-state clock universality through the machine-learning study, which uses data from the six-state clock model on the kagome lattice for training.
We introduce the pushy random walk, where a walker can push multiple obstacles, thereby penetrating large distances in environments with finite obstacle density. This process gives a more realistic depiction of experimentally observed interactions of active particles in dense media. In one dimension, the walker carves out an obstacle-free cavity whose length grows subdiffusively over time. In two dimensions, increasing obstacle density drives a transition from free diffusion to localized behavior, where the walker is trapped within a cavity whose radius again grows subdiffusively with time.
Functional connectivity estimates are highly sensitive to analysis choices and can be dominated by noise when the number of sampled time points is small relative to network dimensionality. This issue is particularly acute in fMRI, where scan resolution is limited. Because scan duration is constrained by practical factors (e.g., motion and fatigue), many datasets remain statistically underpowered for high-dimensional correlation estimation. We introduce a framework that combines diffusion-based structural coarse-graining with spectral noise filtering to recover statistically reliable functional networks from temporally limited data. The method reduces network dimensionality by grouping regions according to diffusion-defined communication. This produces coarse-grained networks with dimensions compatible with available time points, enabling random matrix filtering of noise-dominated modes. We benchmark three common FC pipelines against our approach. We find that raw-signal correlations are strongly influenced by non-stationary fluctuations that can reduce apparent inter-subject variability under limited sampling conditions. In contrast, our pipeline reveals a broader, multimodal landscape of inter-subject variability. These large-scale organization patterns are largely obscured by standard pipelines. Together, these results provide a practical route to reliable functional networks under realistic sampling constraints. This strategy helps separate noise-driven artifacts from reproducible patterns of human brain variability.
We study the ground-state properties of a family of frustrated spin-1/2 Heisenberg models on two- and three-dimensional decorated lattices composed of connected star-shaped units. Each star is built from edge-sharing triangles with an antiferromagnetic interaction on the shared side and ferromagnetic interactions on the others. At a critical coupling ratio, the ideal star model - defined by equal ferromagnetic interactions - exhibits a macroscopically degenerate ground state, which we map onto a site percolation problem on the Lieb lattice. This mapping enables the calculation of exponential ground-state degeneracy and the corresponding residual entropy for square, triangular, honeycomb, and cubic lattices. Remarkably, the residual entropy remains high for all studied lattices, exceeding 60\% of the maximal value ln(2). Despite a gapless quadratic one-magnon spectrum, the low-temperature thermodynamics is governed by exponentially numerous gapped excitations. For a distorted-star variant of the model, the ground-state manifold is equivalent to that of decoupled ferromagnetic clusters, leading to exponential degeneracy with a lower, yet still substantial, residual entropy. At low temperature the system mimics a paramagnetic crystal of non-interacting spins with high spin value ($s=4$ for a square lattice). The obtained results establish a structural design principle for engineering quantum magnets with a high ground-state degeneracy, suggesting promising candidates for enhanced magnetocaloric cooling and quantum thermal machines.
Homing and navigation are fundamental behaviors in biological systems that enable agents to reliably reach a target under uncertainty. We present a Reinforcement Learning (RL) framework to model adaptive homing in continuous two-dimensional domain. In this framework, the agent's state is given by its angular deviation from home, actions correspond to alignment or stochastic reorientation, and learning is driven by a radial-distance-based cost that penalizes motion away from the target. For a single self-propelled agent moving with constant speed, we find that the mean homing time $\langle T_{\mathrm{home}} \rangle$ exhibits a non-monotonic dependence on the rotational diffusion strength $D_r$, with an optimal noise level $D_r^{*}$, revealing a subtle interplay between exploration and goal-directed correction. Extending to two agents with soft repulsion, one agent consistently reaches home faster than the other, while in multi-agents system, repulsion ensures separation and the fastest agent becomes progressively faster as group size increases. Finally comparing the mean homing time $\langle T_{\mathrm{home}} \rangle$ of an Active Brownian Particle (ABP) and RL agent over an identical range of $D_r$, we find that RL trajectories are shorter, less noisy, and consistently faster. Our results show that cost-driven learning, stochastic reorientation, and inter-agent interactions enable efficient adaptive navigation, linking individual and collective homing. This reinforcement learning framework captures key biological features such as feedback-based route learning, randomness to escape unfavorable orientations, and mutual coordination.
We analyze a variant of the Hopfield model that incorporates an unlearning mechanism based on spin correlations in the high-temperature regime. In the large system limit where extensively many patterns are stored, we employ the replica method under the replica symmetric ansatz to characterize the model analytically. Our analysis provides a systematic and self-consistent framework that yields order-parameter equations and stability conditions at finite temperatures over a wide range of parameter settings. The resulting theory accurately captures the behavior of the signal-to-noise ratio, the memory capacity, and the criteria for selecting optimal hyperparameters, in agreement with the qualitative findings of Nokura (1996 \textit{J. Phys. A: Math. Gen.} \textbf{29} 3871). Moreover, the theoretical predictions show good agreement with numerical simulations, supporting the conclusion that unlearning enhances memory capacity by suppressing spurious memories.
The structure and spatial statistical properties of amorphous ellipsoid assemblies have profound scientific and industrial significance in many systems, from cell assays to granular materials. This paper uses a fundamental theoretical relationship for mixture distributions to explain the observations of an extensive X-ray computed tomography study of granular ellipsoidal packings. We study a size-bi-disperse mixture of two types of ellipsoids of revolutions that have the same aspect ratio of alpha approximately equal to 0.57 and differ in size, by about 10% in linear dimension, and compare these to mono-disperse systems of ellipsoids with the same aspect ratio. Jammed configurations with a range of packing densities are achieved by employing different tapping protocols. We numerically interrogate the final packing configurations by analyses of the local packing fraction distributions calculated from the Voronoi diagrams. Our main finding is that the bi-disperse ellipsoidal packings studied here can be interpreted as a mixture of two uncorrelated mono-disperse packings, insensitive to the compaction protocol. Our results are consolidated by showing that the local packing fraction shows no correlation beyond their first shell of neighbours in the binary mixtures. We propose a model of uncorrelated binary mixture distribution that describes the observed experimental data with high accuracy. This analysis framework will enable future studies to test whether the observed mean-field behaviour is specific to the particular granular system or the specific parameter values studied here or if it is observed more broadly in other bi-disperse non-spherical particle systems.
Although the Transformer architecture has revolutionized artificial intelligence, its underlying mechanisms remain largely heuristic and lack a unified physical theory. In this work, we propose a first-principles framework for information dynamics, treating the attention mechanism as a physical system governed by the principle of least action rather than as an algorithmic optimization. By mapping information states to a Riemannian manifold with the Fisher information metric, we derive the intelligence Lagrangian. We show that the softmax function corresponds to the unique thermodynamic equilibrium state that minimizes the Helmholtz free energy of the information gas. In addition, we identify the query-key interaction as an electrodynamic coupling between an external field and an intrinsic dipole moment. This theory establishes the first law of information thermodynamics, unifying inference (mechanical work) and learning (chemical evolution). It also explains emergent phenomena, such as scaling laws and grokking, as phase transitions characterized by the divergence of specific heat. Finally, we discuss how rotational symmetry breaking in the attention manifold generates massless Goldstone bosons, providing a field-theoretic perspective on rotary positional embeddings (RoPE). Our work connects Statistical Physics and Deep Learning, laying the groundwork for a general theory of physics-based intelligence.
We study the 1d swarmalator model augmented with time delayed coupling. Along with the familiar sync, async, and phase wave states, we find a family of unsteady states where the order parameters are time periodic, sometimes with clean oscillations, sometimes with irregular vacillations. The unsteady states are born in two ways: via a Hopf bifurcation from the phase wave, and a zero eigenvalue bifurcation from the async state. We find both of these boundary curves analytically. A surprising result is that stabilities of the async and sync states are independent of the delay τ; they depend only on the coupling strength.
We investigate operator growth in a Brownian spin Sachdev--Ye--Kitaev (SYK) model with random all-to-all interactions, focusing on the full operator-size distribution. For Hamiltonians containing interactions of order two up to $L$, we derive a closed master equation for the Pauli-string expansion coefficients and recast their dynamics into a generating-function formulation suitable for the large-$N$ limit. This approach allows us to diagonalize the leading-order evolution operator explicitly and obtain exact solutions for arbitrary initial operator distributions, including the effects of decoherence. Going beyond leading order, we develop a systematic $1/N$ expansion that captures higher-order corrections to the operator-size dynamics and the late-time behavior. Our results demonstrate that higher-order effects play a crucial role in operator scrambling and that the full operator-size distribution provides a more refined probe of quantum chaos in Brownian and open quantum systems.
Collective motion in fish schools exemplifies emergent self-organization in active matter systems, yet computational tools for simulating and analyzing these dynamics remain fragmented across research groups. We present dewi-kadita, an open-source Python library implementing the three-dimensional Couzin zone-based model with comprehensive entropy diagnostics tailored for marine collective behavior research. The library introduces seven information-theoretic metrics -- school cohesion entropy, polarization entropy, depth stratification entropy, angular momentum entropy, nearest-neighbor entropy, velocity correlation entropy, and school shape entropy -- that characterize distinct organizational features inaccessible to classical order parameters. These metrics combine into an Oceanic Schooling Index (OSI) providing a single scalar measure of collective disorder. Validation across four canonical configurations (swarm, torus, dynamic parallel, highly parallel) confirms correct reproduction of known phase behaviors: the swarm maintains disorder with polarization $P < 0.1$ and OSI $\approx 0.71$, while the highly parallel state achieves $P = 0.998$ with OSI $= 0.24$ and velocity correlation entropy vanishing to zero. The entropy framework successfully discriminates the torus and dynamic parallel configurations that exhibit comparable order parameter magnitudes through different organizational mechanisms. Numba just-in-time (JIT) compilation accelerates pairwise interaction calculations by $10$--$100\times$, enabling simulations of $150$--$250$ agents over $1000$--$2000$ time steps within five minutes on standard workstation hardware. NetCDF4 output ensures interoperability with oceanographic analysis tools. The library addresses the need for standardized, reproducible infrastructure in collective behavior modeling analogous to established molecular dynamics codes.
The seminal 2009 paper by Bernard, Krauth, and Wilson marked a paradigm shift in Monte Carlo sampling. By abandoning the restrictive condition of detailed balance in favor of the more fundamental principle of global balance, they introduced the Event-Chain Monte Carlo (ECMC) algorithm, which achieves rejection-free, deterministic sampling for hard spheres. This breakthrough demonstrated that persistent, directional dynamics could dramatically accelerate equilibration in dense particle systems. In this commentary, we review this foundational work and elucidate its underlying mechanism using the broader Event-Driven Monte Carlo (EDMC) framework developed in subsequent years. We show how the original hard-sphere concept naturally generalizes to continuous potentials and modern lifted Markov chain formalisms, transforming a surprising specific result into a powerful general class of sampling algorithms.