2025-07-15 | | Total: 27
Plasmas and other systems with long-range interactions are commonly found in non-equilibrium steady states that are outside traditional Boltzmann-Gibbs statistics, but can be described using generalized statistical mechanics frameworks such as superstatistics, where steady states are treated as superpositions of canonical ensembles under a temperature distribution. In this work we solve the problem of inferring the possible steady states of a composite system AB where subsystem A is described by superstatistics and EAB=EA+EB. Our result establishes a closure property of superstatistics, namely that A is described by superstatistics if and only if AB and B are also superstatistical with the same temperature distribution. Some consequences of this result are discussed, such as the impossibility of local thermal equilibrium (LTE) for additive subsystems in non-canonical steady states.
We analyse the transitions between established phases of the El Niño Southern Oscillation (ENSO) by surveying the daily data of the Southern Oscillation Index from an entropic viewpoint using the framework of stochastic Statistical Physics. We evaluate the variation of entropy produced due to each recorded path of that index during each transition as well as taking only into consideration the beginning and the end of the change between phases and verified both integral fluctuation relations. The statistical results show that these entropy variations have not been extreme entropic events; only the transition between the strong 1999−2000 La Niña to the moderate 2002−2003 El Niño is at the edge of being so. With that, the present work opens a long and winding avenue of research over the application of stochastic Statistical Physics to Climate Dynamics.
Many have speculated whether nonequilibrium systems obey principles of maximum or minimum entropy production. In this work, we use stochastic thermodynamics to derive the condition for the minimum entropy production state (MEPS) for continuous-time Markov chains (CTMCs), even far from equilibrium. We show that real nonequilibrium steady states (NESS) generally violate both the MINEP and MAXEP principles. However, through numerical sampling of large interconnected CTMCs, we find that as system size increases, the steady-state entropy production tends to converge toward the minimum. This suggests that large nonequilibrium systems may self-organize to make efficient use of thermodynamic resources, offering a nuanced perspective on the longstanding debate between MAXEP and MINEP.
Continuous spin models with long-range interactions of the form r−σ, where r is the distance between two spins and σ controls the decay of the interaction, exhibit enhanced order that competes with thermal disturbances, leading to a richer variety of phases and types of phase transitions. In-depth research in this area not only aids in comprehending the complex behaviors in theoretical models but also provides valuable insights into the diverse phase transitions observed in real materials. Here, we identify that the true long-range ordered phase encompasses distinct scaling regimes, which we term Enhance Long-Range Ordered (EnLRO) and Reduce Long-Range Ordered (ReLRO) regimes. In the one-dimensional XY model, the crossover from EnLRO to ReLRO regimes occurs around σ≈1.575, while in two dimensions, the crossover happens near σ≈3.2. Applying finite-size scaling analysis, we extract the critical exponents that characterize the order-to-disorder phase transitions in the EnLRO and ReLRO regimes, constructing comprehensive phase diagrams. The analysis is further extended to the 1D and 2D long-range Heisenberg models, where we find the EnLRO-ReLRO crossover at σ≈1.575 and σ≈3.22, respectively. The similar crossover points suggest that the distinction between EnLRO and ReLRO regimes is a generic feature in continuous spin models with long-range interactions.
Random graph models have been instrumental in characterizing complex networks, but chemical reaction networks (CRNs) are better represented as hypergraphs. Traditional models of random CRNs often reduce CRNs to bipartite graphs, representing species and reactions as distinct nodes, or simpler derived graphs, which can obscure the relationship between the statistical properties of these representations and the physical characteristics of the CRN. We introduce a straightforward model for generating random CRNs that preserves their hypergraph structure as well as atomic composition, enabling the direct study of chemically relevant features. Notably, our approach distinguishes two notions of connectivity that are equivalent in graphs but differ fundamentally in hypergraphs. These notions exhibit percolation-like phase transitions, which we analyze in detail. The first type of connectivity has relevance to steady-state synthesis and transduction, determining the effective reactions an open CRN can perform at steady state. The second type is suitable to identify which species can be produced from a given initial set of species in a closed CRN. Our findings highlight the importance of hypergraph-based modeling for uncovering the complex behaviors of CRNs.
The possibility that spatial geometry may emerge from the entanglement structure of a quantum many-body system is a subject of fundamental interest. Here, we propose and numerically test a candidate distance metric in 1D, d_E, defined purely from quantum mutual information (I) via the relation d_E = K_0 / sqrt(I). Using large-scale density-matrix renormalization group (DMRG) simulations, we compute this emergent distance for the ground state of the 1D spin-1/2 XXZ chain, a canonical model system. Our simulations show that in the quantum critical phase at anisotropy Δ=1.0, the mutual information exhibits a power-law decay consistent with the emergence of a valid metric space. In stark contrast, within the gapped, antiferromagnetic phase (Δ=2.0), where mutual information decays exponentially, the emergent distance grows exponentially, a behavior inconsistent with the triangle inequality. These results provide numerical evidence that this information-theoretic definition can yield a well-behaved geometry in critical systems, offering a quantitative tool for probing quantum phases and motivating further analytical investigation into the foundations of emergent space.
The analytical model of a glass-forming system is formulated within the formalism analogous to gauge theory constructions in quantum field theory. This work explores the scope of the proposed approach and investigates the equilibrium behavior of the model under the mean-field approximation. The analysis reveals three possible equilibrium scenarios, only one of which exhibits strong glass-forming ability. Additionally, an upper limit for the initial temperature is identified, beyond which quenching into the glassy state becomes impossible.
From bird flocks and fish schools to migrating cell sheets, collective motion is a ubiquitous biological phenomenon that inspires quantitative modeling through self-propelled particle (SPP) frameworks. Conventional SPP models prescribe either distance-based (metric) or rank-based (topological) interactions; however, empirical studies indicate that real groups may blend both types of interaction. Motivated by this graded perception, we introduce a new three-dimensional SPP model in which metric and topological alignments act simultaneously and are weighted by a single tunable mixing parameter called the interaction parameter. Large-scale simulations spanning a wide ranges of interaction parameters and densities revealed rich dynamics. Even when the global order parameter is low, cluster-level analysis with HDBSCAN shows that particles self-organize into several spatially distinct but internally well-aligned sub-flocks, exposing a hidden layer of order. Most importantly, an intermediate balance in which metric and topological cues contribute almost equally maximizes the global order parameter and markedly improves robustness to density variations. Numerical experiments supported by linear stability analysis demonstrate that activating metric and topological interactions in concert bridges the traditional modeling dichotomy and furnishes a more adaptive and resilient framework for collective motion. Therefore, the proposed model therefore provides a versatile platform for exploring mixed-interaction effects in biological and engineered multi-agent systems.
Finding a target in a complex environment is a fundamental challenge in nature, from chemical reactions to sperm reaching an egg. An effective strategy to reduce the time needed to reach a target is to deploy many searchers, increasing the likelihood that at least one will succeed by using the statistics of rare events. When the underlying stochastic process involves broadly distributed step sizes, rare long jumps dominate the dynamics, making the use of multiple searchers particularly powerful. We investigate the statistics of extreme events for the mean first passage time in a system of N independent walkers moving with jumps distributed according to a power law, where target-reaching is governed by single, large fluctuations. We show that the mean first passage time of the fastest walker scales as ⟨τN⟩∼1/N, representing a dramatic speed-up compared to classical Brownian search strategies. From this, we derive a scaling law relating the number of walkers required to reach a target within a given time to the size X of the search region. As an application, we model biological fertilization, predicting how the optimal number of spermatozoa scales with uterus size across species. Our predictions match empirical data, suggesting that evolution may have exploited rare-event dynamics and broadly distributed motion to optimize reproductive success. This theory applies broadly to any population of searchers operating within a region of size X, providing a universal framework for efficient search in disordered or high-dimensional environments.
We consider the finite-timestep Brownian dynamics of a single particle confined in one dimension, with a nonuniform temperature profile. In such an open-boundary scenario, one cannot observe any net probability current in the nonequilibrium steady state (NESS). On the other hand, the nonequilibrium nature of this system is exhibited through the asymmetry in forward and backward transition probabilities, as is reported in this work through the stochastic simulation analysis and theoretical arguments. The irreversibility becomes prominent nearby the temperature interface. We propose that the observed irreversibility can be accounted for via a virtual-gyration scenario, while the collapse of virtual gyrations upon the one-dimensional coordinate leads to the absence of probability current.
The spectrum of Markov generators encodes physical information beyond simple decay and oscillation, which reflects irreversibility and governs the structure of correlation functions. In this work, we prove an ellipse theorem that provides a universal thermodynamic geometric constraint on the spectrum of Markov rate matrices. The theorem states that all eigenvalues lie within a specific ellipse in the complex plane. In particular, the imaginary parts of the spectrum, which indicate oscillatory modes, are bounded by the maximum thermodynamic force associated with individual transitions. This spectral bound further constrains the possible values of correlation functions of two arbitrary observables. Finally, we compare our result with a previously proposed conjecture, which remains an open problem and warrants further investigation.
Ferromagnetic exponential random graph models (ERGMs) are nonlinear exponential tilts of Erdős-Rényi models, under which the presence of certain subgraphs such as triangles may be emphasized. These models are mixtures of metastable wells which each behave macroscopically like new Erdős-Rényi models themselves, exhibiting the same laws of large numbers for the overall edge count as well as all subgraph counts. However, the microscopic fluctuations of these quantities remained elusive for some time. Building on a recent breakthrough by Fang, Liu, Shao and Zhao [FLSZ24] driven by Stein's method, we prove quantitative central limit theorems (CLTs) for these quantities and more in metastable wells under ferromagnetic ERGMs. One main novelty of our results is that they apply also in the supercritical (low temperature) regime of parameters, which has previously been relatively unexplored. To accomplish this, we develop a novel probabilistic technique based on the careful analysis of the evolution of relevant quantities under the ERGM Glauber dynamics. Our technique allows us to deliver the main input to the method developed by [FLSZ24], which is the fact that the fluctuations of subgraph counts are driven by those of the overall edge count. This was first shown for the triangle count by Sambale and Sinulis [SS20] in the Dobrushin (very high temperature) regime via functional-analytic methods. We feel our technique clarifies the underlying mechanisms at play, and it also supplies improved bounds on the Wasserstein and Kolmogorov distances between the observables at hand and the limiting Gaussians, as compared to the results of [FLSZ24] in the subcritical (high temperature) regime beyond the Dobrushin regime. Moreover, our technique is flexible enough to also yield quantitative CLTs for vertex degrees and local subgraph counts, which have not appeared before in any parameter regime.
The Hubbard model at U→∞ has recently been shown to have resonant valence bond (RVB) ground states on the corner-sharing sawtooth and pyrochlore lattices in the dilute doping limit of a single vacancy. In an effort to further generalize those results, I study how the ground state is modified when not all corners are shared between two tetrahedra as in the quasi-1D lattices of a pyrochlore stripe, and how to approach the problem in the case of finite doping. Using a non-Abelian version of the flux inequality, the tetrahedron chain is shown to have degenerate RVB-like ground states. The Bethe ansatz (BA) is adapted to solve the sawtooth chain with spinless or spin-polarized fermions and multiple holons, which is the first example of applying BA to a quasi-1D lattice.
Time delay is ubiquitous in many experimental and real-world situations. It is often unclear whether time delay plays a significant role in observed phenomena, and if it does, how long the time lag really is. This would be invaluable knowledge when analyzing and modeling such systems. Hitherto, no universal method is available by which the time delay can be inferred. To address this problem, we propose and demonstrate two different methods to infer time delay in overdamped Langevin systems with delayed feedback. In the first part, we focus on the power spectral density based on the positional data and use a characteristic signature of the time delay to infer the delay time. In limiting cases, we establish a direct relation of the observations made for nonlinear time-delayed feedback forces to analytical results obtained for the linear system. In other situations despite the absence of this direct relation, the characteristic signature remains and can be exploited by a semiautomatic method to infer the delay time. Furthermore, it may not always desirable or possible to observe a system for a long time to infer dependencies and parameters. Thus, in the second part, we propose for the first time a probing method combined with a neural network to infer the delay time, which requires only short observation time series. These proposed methods for inferring time delays in stochastic systems may prove to be valuable tools for gaining deeper insight into the role of delay across a wide range of applications -- from the behavior of individual colloidal particles under feedback control to emergent collective phenomena such as flocking and swarming.
The nonequilibrium dynamics of a quantum impurity immersed with a finite velocity in a one-dimensional system of weakly interacting bosons is studied. We uncover and characterize different regimes of relaxation dynamics. We find that the final impurity velocity remains constant in a large interval of sufficiently big and realistic initial velocities. The underlying physical mechanism is the emission of the dispersive density shock waves that carry away the excess of the initial impurity momentum, while locally the system remains in the same stationary state. In contrast, a heavy impurity with the same coupling constant relaxes differently and the regime of constant final velocity disappears. Furthermore, a fast heavy impurity exhibits damped velocity oscillations in time before reaching a stationary state. This process is accompanied by the oscillations of the local depletion of the boson density around the impurity, until their positions coincide and they continue the motion together. Decreasing the impurity-boson coupling or increasing the strength of repulsion between bosons, the oscillations get amplified. In the case of a heavy impurity with the mass bigger than the critical one, the ground state energy as a function of momentum exhibits cusps and metastable branches. We show that they manifest themselves by a soliton emission, a considerable slowing down of the relaxation, and a change of the impurity direction of motion with respect to the initial one.
Neural networks storing multiple discrete attractors are canonical models of biological memory. Previously, the dynamical stability of such networks could only be guaranteed under highly restrictive conditions. Here, we derive a theory of the local stability of discrete fixed points in a broad class of networks with graded neural activities and in the presence of noise. By directly analyzing the bulk and outliers of the Jacobian spectrum, we show that all fixed points are stable below a critical load that is distinct from the classical \textit{critical capacity} and depends on the statistics of neural activities in the fixed points as well as the single-neuron activation function. Our analysis highlights the computational benefits of threshold-linear activation and sparse-like patterns.
Understanding how a quantum many-body state is maintained stably as a nonequilibrium steady state is of fundamental and practical importance for exploration and exploitation of open quantum systems. We establish a general equivalent condition for an open quantum many-body system governed by the Gorini-Kossakowski-Sudarshan-Lindblad dynamics under local drive and/or dissipation to have a quantum independent and identically distributed (i.i.d.) steady state. We present a sufficient condition for a system to have a quantum i.i.d. steady state by identifying a set of operators that commute with arbitrary quantum i.i.d. states. In particular, a set of quantum i.i.d. states is found to be an invariant subset of time evolution superoperators for systems that satisfy the sufficient condition. These findings not only identify a class of models with exactly solvable steady states but also lead to a no-go theorem that precludes quantum entanglement and spatial correlations in a broad class of quantum many-body steady states in a dissipative environment.
We derive a thermodynamically consistent quantum master equation that satisfies locality for quadratic systems coupled to independent and identical baths at each site. We show that the quasi-local Redfield equation coincides exactly with the Davies equation, which satisfies the detailed-balance condition, due to cancellation of quantum coherence generated by each bath. This derivation does not rely on the secular approximation, which fails in systems with vanishing energy-level spacings. We discuss generalizations of our result to slowly driven quadratic systems and generic quantum many-body systems. Our result paves the way to a thermodynamically consistent description of quantum many-body systems.
The maximum energy of the EPR model on a weighted graph is known to be upper-bounded by the sum of the total weight and the value of maximum-weight fractional matching~(MWFM). Recently, Apte, Parekh and Sud~(APS) conjecture that the bound could be strengthened by replacing MWFM with maximum weight matching~(MWM). Here we test this conjecture on a special class of regular graphs that Henning and Yeo constructed many years ago. On this class of regular graphs, MWMs achieve tight lower bounds. As for the maximum energy of the EPR model, we have recently devised a new algorithm called Fractional Entanglement Distribution~(FED) based on quasi-homogeneous fractional matchings, which could achieve rather high accuracy. Applying the FED algorithm to the EPR model on Henning-Yeo graphs, we could thus obtain energy as high as possible and matching value as low as possible, and then make high-precision tests of the APS conjecture. Nevertheless, our numerical results do not show any evidence that the APS conjecture could be violated.
In recent years, graph-based machine learning techniques, such as reinforcement learning and graph neural networks, have garnered significant attention. While some recent studies have started to explore the relationship between the graph structure of neural networks and their predictive performance, they often limit themselves to a narrow range of model networks, particularly lacking mesoscale structures such as communities. Our work advances this area by conducting a more comprehensive investigation, incorporating realistic network structures characterized by heterogeneous degree distributions and community structures, which are typical characteristics of many real networks. These community structures offer a nuanced perspective on network architecture. Our analysis employs model networks such as random and scale-free networks, alongside a comparison with a biological neural network and its subsets for more detailed analysis. We examine the impact of these structural attributes on the performance of image classification tasks. Our findings reveal that structural properties do affect performance to some extent. Specifically, networks featuring coherent, densely interconnected communities demonstrate enhanced learning capabilities. The comparison with the biological neural network emphasizes the relevance of our findings to real-world structures, suggesting an intriguing connection worth further exploration. This study contributes meaningfully to network science and machine learning, providing insights that could inspire the design of more biologically informed neural networks.
Non-Bloch band theory serves as a cornerstone for understanding intriguing non-Hermitian phenomena, such as the skin effect and extreme spectral sensitivity to boundary conditions. Yet this theory hinges on translational symmetry and thus breaks down in disordered systems. Here, we develop a real-space Lyapunov formulation of band theory that governs the spectra and eigenstates of disordered non-Hermitian systems. This framework yields universal non-Hermitian Thouless relations linking spectral density and localization to Lyapunov exponents under different boundary conditions. We further identify an exact topological criterion: skin modes and Anderson-localized modes correspond to nonzero and zero winding numbers, respectively, revealing the topological nature of the skin-Anderson transition. This transition is dictated by an essential Lyapunov exponent and gives rise to novel unidirectional critical states. Our formulation provides a unified and exact description of spectra and localization in generic one-dimensional non-Hermitian systems without translational symmetry, offering new insights into the interplay among non-Hermiticity, disorder, and topology.
Paper fragments in free fall constitute a simple yet paradigmatic mechanical system exhibiting remarkably complex motions. Despite a long history of investigation, this system has defied comprehensive first-principles modeling, motivating the development of phenomenological and experimental approaches to classify the free-fall dynamics of small paper fragments. Here we apply the Bandt-Pompe symbolization method to extract high-dimensional features corresponding to ordinal-pattern transitions (so-called ordinal networks) from observed area time series of video-recorded falling papers shaped as circles, squares, hexagons, and crosses. We then represent each trajectory as a node in a weighted similarity network, with edges encoding pairwise dynamical similarity, and identify motion classes via community detection. Our method automatically clusters trajectories into tumbling and chaotic falls in excellent agreement with expert visual classification. Notably, it outperforms previous approaches based on classical physical features derived from complete three-dimensional trajectories -- especially for cross-shaped papers -- and requires no prior specification of the number of motion classes. We further find that trajectories diverging from expert classifications occupy more central positions in the similarity network, suggesting more complex and ambiguous dynamic behavior.
We investigate defects in the two-dimensional transverse-field Ising ferromagnet on periodic L×L lattices after quantum annealing from high to vanishing field. With exact numerical solutions for L≤6, we observe the expected critical Kibble-Zurek (KZ) time scale ∝Lz+1/ν (with z=1 and 1/ν≈1.59) at the quantum phase transition. We also observe KZ scaling of the ground-state fidelity at the end of the process. The excitations evolve by coarsening dynamics of confined defects, with a time scale ∝L2, and interface fluctuations of system-spanning defects, with life time ∝L3. We build on analogies with classical simulated annealing, where we characterize system-spanning defects in detail and find differences in the dynamic scales of domain walls with winding numbers W=(1,0)/(0,1) (horizontal/vertical) and W=(1,1) (diagonal). They decay on time scales ∝L3 (which applies also to system-spanning domains in systems with open boundaries) and ∝L3.4, respectively, when imposed in the ordered phase. As a consequence of L3.4 exceeding the classical KZ scale Lz+1/ν=L3.17 the probability of W=(1,1) domains in SA scales with the KZ exponent even in the final T=0 state. In QA, also the W=(1,0)/(0,1) domains are controlled by the KZ time scale L2.59. The L3 scale can nevertheless be detected in the excited states, using a method that we develop that should also be applicable in QA experiments.
We obtain a form of the local thermodynamic equilibrium density operator which is invariant under pseudo-gauge transformations of the stress-energy and the spin tensors. This operator is an excellent candidate to describe the dynamics of a system which is assumed to achieve local equilibrium from a pseudo-gauge invariant quantum state, a situation which is believed to occur, for instance, in nuclear collisions at very high energy. As a consequence of pseudo-gauge invariance, the ambiguity affecting the predictions of mean values of observables from a local equilibrium state can be removed.
In these notes, we present a rigorous and self-contained introduction to the fundamental concepts and methods of quantum many-body theory. The text is designed to provide a solid theoretical foundation for the study of interacting quantum systems, combining clarity with mathematical precision. Core topics are developed systematically, with detailed derivations and comprehensive proofs that aim to make the material accessible to graduate students and beginning PhD students. Special attention is given to formal consistency and pedagogical structure, so as to guide the reader through both the conceptual and technical aspects of the subject. This work is intended as a reliable starting point for further exploration and research in modern quantum many-body physics.