2026-07-17 | | Total: 19
We introduce a thermodynamic theory of voting and show that it provides a good description of distribution of party votes in EU elections. The theory traces parallels between system energies of coupled nonlinear oscillators and party vote fractions. Such a classical system evolution is characterized by the conservation of total energy and probability norm that leads to the Rayleigh-Jeans (RJ) thermalization and condensation at low energy states. A similar thermalization also describes the wealth inequality in society. This feature belongs to the phenomena of constraint driven condensation known in statistical mechanics. We show that the RJ theory well depicts the Lorenz and Pareto curves obtained from the EU vote results. The theory also recovers the dispersion of votes between candidates of first round presidential elections in France.
Starting from a \emph{probabilistic numerics} approach to the Fermion sign problem in path integral Monte Carlo, we recast the arithmetic calculation of a Fermionic observable as a statistical inference problem. We develop approaches that learn the behaviour of Fermion exchange cycles binned by the conjugacy class of the permutation group (which we term `permutation family'). This extends the work of DuBois, Brown and Alder\cite{dubois2017overcoming} to inhomogeneous and more complex systems. Monte Carlo samples are used to train models for both the probability of a permutation family and the energy of this set of exchange permutations. The overall Fermionic energy is then directly inferred from these models, without using a direct ratio estimator on the Monte Carlo samples. By imposing physical understanding as inductive priors, we produce accurate and useful fits that remain robust even in regimes with severe sign problems. We generalise the linear (ideal-gas style) models of DuBois et al. with Bayesian priors that enforce the intuitive models of Feynman\cite{Feynman1953A} at their asymptotic limits. These linear models serve as the baseline for a Long Short-Term Memory (LSTM) neural network, which is tasked with learning only the residual many-body \emph{correlations} on top of the physical model. We develop active important sampling methods driven by these models, which direct the Monte Carlo chains toward undersampled permutation regions, to efficiently reduce the variance in the observable. We apply this framework to small experiments on benchmark systems: the spin-polarised uniform electron gas, and electrons in a 2D harmonic confining potential. In both cases we demonstrate that this inference-based framework can extract stable energies in regimes where direct Monte Carlo sampling fails due to the sign problem.
We derive a closed-form analytical expression for the distribution of eccentricities (DoE) in random regular graphs (RRGs) that consist of $N$ nodes of degree $c$. The DoE is given by the tail distribution $P(E > \ell) \simeq 1 - \exp \left[ - \exp \left( - \frac{ e^{b \ell} - μ}β \right) \right]$, where the distance $\ell$ takes integer values, $b = \ln (c-1)$ is the shape parameter, $β= \frac{c-2}{c} N$ is the scale parameter and $μ= \frac{c-2}{c} N \ln N$ is the location parameter. By providing the full distribution rather than a single characteristic length scale, we present a detailed view of the large-scale structure. In spite of the fact that the degrees of all the nodes are the same, their eccentricities exhibit non-trivial variations. We derive a closed-form expression for the mean eccentricity, which is given by $\langle E \rangle \simeq \frac{\ln N}{\ln (c-1)} + \frac{\ln \ln N}{\ln (c-1)} - \frac{ \ln c - \ln (c-2) }{ \ln (c-1) } + \frac{1}{2}$. We calculate the mode of the DoE, which exhibits a staircase profile as a function of the network size. Interestingly, the mode is given by $E_{\rm mode} ={\rm Round} \left( \langle E \rangle \right)$, where ${\rm Round}( x )$ is the nearest integer to $x$. We also calculate the variance ${\rm Var}(E)$ and show that it exhibits oscillations as a function of the network size $N$. The results presented in this paper may serve as benchmarks for algorithmic approaches to eccentricity calculations in large sparse networks. The eccentricities are important in practical applications such as broadcasting and global dissemination, where the network performance is determined by the longest delay times.
Nonlinear thermodynamic computers based on Langevin dynamics exploit thermal fluctuations as a physical substrate for computation. Recent work has shown that quartic-confined fluctuating degrees of freedom can act as thermodynamic neurons capable of nonlinear function approximation at finite observation times. Here we extend this paradigm from mean-only readout to moment-resolved readout. Instead of representing each driven reservoir solely by its first moment, we construct a response vector from the elementwise raw polynomial moments \(\mathbb{E}[\bm{x}]\), \(\mathbb{E}[\bm{x}^{\odot 2}]\), and \(\mathbb{E}[\bm{x}^{\odot 4}]\). These observables combine displacement and central-shape contributions and are naturally aligned with the linear, quadratic, and quartic terms of the local driven dynamics. We further introduce a heterogeneous multi-reservoir architecture in which three reservoirs with distinct initialization and training histories form a joint \(2304\)-dimensional response representation. Under the fixed MNIST \(60000/10000\) reproduction protocol, feature-level fusion achieves the best observed accuracy of \(9695/10000=96.95\%\), compared with \(9682/10000=96.82\%\) for the strongest single-reservoir model and \(9684/10000=96.84\%\) for equal-weight logit averaging. An exact paired McNemar test does not establish a statistically significant improvement over the strongest single reservoir, but the ablation and wrong-set overlap results provide suggestive evidence of complementary classification errors. These results motivate higher-order polynomial-moment readout and reservoir heterogeneity as candidate design principles for finite-time Langevin computing.
Universality is usually associated with asymptotic behavior becoming independent of details of the initial state. Using a unified renormalization-group (RG) framework for self-similar dynamics, we show that retaining relevant length scales leads naturally to distinct classes of memory-retaining fixed points. A modified density-dependent diffusion model reveals a general scaling structure in which initial-condition information remains asymptotically relevant, with the Barenblatt equation emerging as a special case. These results provide a new perspective on universality and anomalous scaling.
We study a system of two ferromagnetic one-dimensional Ising chains coupled to a thermal bath, which are driven out of equilibrium by being moved past one another at a constant speed. We show that even at modest speeds, magnetic friction between the two chains significantly increases the ability of the system to order. In particular, at inverse temperature $β$, Ising coupling $J$, and sliding speed $v$, the dynamics retains memory of its initial magnetization for a time that increases from $\exp(O(βJ))$ at $v = 0$ to $\exp(O((βJ)^2v\ln v))$ when $v>v_c$, where $v_c$ is a small constant. Magnetic friction thus provides a simple mechanism for parametrically slowing down thermalization in a one-dimensional magnet.
Thermal rectification (TR) typically relies on structural asymmetry or material heterogeneity. Here, we show that gravity alone can induce and modulate TR in gaseous systems. Using a minimal model of a single gas particle confined in a two-dimensional channel between heat baths at different temperatures, we analytically demonstrate that gravitational fields generate TR, enabling perfect unidirectional heat conduction across broad gravitational parameter ranges. This effect exhibits an intrinsic trade-off between rectification efficiency and heat power. Extending beyond the single-particle limit, numerical simulations confirm that gravitationally induced TR persists in interacting many-particle systems. Notably, in interacting gas mixtures, the rectification direction can be reversed. Gravity-mediated control of heat currents thus provides new fundamental insights into thermal transport and suggests design principles for gaseous thermal diodes.
Deep thermalisation characterises the emergence of universal quantum state ensembles on subsystems due to projective measurements on their complement. We study the notion of locality, or lack thereof, in this phenomenon by considering a subsystem partitioned into two disjoint subregions which remain causally disconnected at all times under unitary dynamics. We show that the onset of deep thermalisation in this geometry is fundamentally bounded by measurement-induced entanglement teleportation between the subregions. While measurements on the environment generate entanglement across the disconnected partitions -- suggesting an apparent non-locality -- we demonstrate that generic locally interacting systems exhibit an emergent locality. Specifically, the timescales for both deep thermalisation and entanglement teleportation scale logarithmically with the distance separating the subregions. Exceptions to this include special circuits where the randomness of the measurement outcomes is perfectly transmitted to the ensemble of states of the subsystem, conditioned on the outcomes; in such cases the timescale for deep thermalisation is finite leading to genuine non-locality.
We propose an experimentally motivated model in which the Mpemba effect can be controlled through the system's inertia. The model describes a polymer undergoing a denaturation transition whose force-extension curve contains a plateau that slows relaxation. When the system is initially prepared at a higher temperature or under a weaker stretching force, the bead accumulates greater kinetic energy, allowing it to cross the plateau more rapidly and thereby producing the Mpemba effect. Increasing the bead mass broadens the range of initial temperatures over which this mechanism operates. A similar mechanism also generates the inverse Mpemba effect.
Stochastic binary networks are widely used to describe collective dynamics in complex systems and to perform neuromorphic computation, yet realistic networks often contain both asymmetric interactions and finite signal propagation times that fall outside conventional theories. Here we study stochastic binary networks with asymmetric and time-delayed interactions motivated by experimental observations in coupled superparamagnetic tunnel junctions. We find that time delay fundamentally reshapes the dynamics induced by anti-symmetric couplings, producing strong oscillatory temporal correlations consistent with experiment. At the same time, sufficiently long delays drive the steady-state probabilities toward equal state occupations even in strongly coupled systems. These apparently featureless probability distributions coexist with pronounced temporal correlations, distinguishing them from equilibrium high-temperature behavior. We further show analytically that delay-induced uniform distributions emerge in a broad class of stochastic networks, while symmetry-breaking bias fields restore interaction-dependent steady states with qualitatively modified behavior. Simulations of networks with five coupled spins demonstrate that these effects persist beyond minimal systems with only two spins. Our results establish a unified framework for stochastic binary networks in the intermediate regime between symmetric instantaneous interactions and asymmetric or time-delayed interactions, and suggest that asymmetry and delay can be exploited as functional resources in neuromorphic hardware and complex network dynamics.
Quantum many-body scars provide an exception to generic thermalising dynamics, but their relation to integrability remains unclear. Here we apply the periodic-orbit framework to the integrable XYZ/XXZ spin-1/2 chain, developing an energy-resolved approach that enables us to locate and track periodic orbits across the spectrum. We identify families of scarred orbits and resolve their supporting eigenstates in terms of Bethe-ansatz quantum numbers, revealing a common string core dressed by zero-momentum magnon pairs. This allows us to reconstruct both the scarred eigenstates and the associated towers directly within the Bethe ansatz, and explain the equidistant tower spacing analytically from the decoupling of zero-momentum magnon pairs in the Bethe equations, providing a microscopic realisation of scar phenomenology in an integrable setting. We further show that this structure persists upon breaking integrability with a transverse field, with the periodic orbits continuing to govern the dynamics. Our results establish a direct connection between the periodic-orbit picture of quantum scarring and the algebraic structure of integrable models, showing that key features of scar dynamics can be understood analytically within the Bethe ansatz framework.
The stochastic process model of rough surface contact, widely known as Persson's theory of contact, serves as a representative multi-scale model that has been extensively applied across various fields of tribology. In this chapter, we briefly introduce the background of the development of Persson's theory of contact. We thoroughly discuss Persson's theory for purely normal elastic contact, with a special focus on solving the probability density of the contact pressure and the interfacial gap using partial differential equations. Subsequent applications of these fundamental results in addressing more complex interfacial properties in other fields of tribology are also examined. Finally, several recommendations regarding future studies of Persson's theory are proposed. This review article is expected to assist researchers in quickly familiarizing themselves with the current state of the art of Persson's theory and to attract more attention from tribologists and solid mechanicians, thereby contributing to the development and application of Persson's theory of contact.
In two dimensions, the geometric constraints due to Pauli blocking and conservation laws lead to the even-odd effect exhibited by the electron-electron scattering lengths: electron-electron collisions are more efficient at relaxing the even angular harmonics of the distribution function than the odd ones. Inspired by a recent experiment on the magnetotransport in the Corbino disk geometry, we numerically analyze the electron flows in this geometry across all the regimes. We predict a clear signature of the even-odd effect - enhancement of the resistance sensitivity $\partial R/\partial(B^2)$ at small magnetic fields $B\rightarrow 0$. This enhancement is most prominent at the crossover from the ballistic to the tomographic regime, and gradually disappears when the temperature is further increased. Our estimates suggest that in the temperature range of the experiment, the effect should be small. This implies that the attribution of the anomalous scaling of the kinematic viscosity, that was observed in the experiment, to the even-odd effect might need more careful consideration. As a side note, we show how the method of characteristics can be extended to treat the long-lived odd harmonics, which allows one to recast the linearized Boltzmann equation as a system of integral ones.
We numerically analyze a discretized version of the quantum contact process, a non-equilibrium model having a competition between a coherent infection process and a dissipative recovery process. Previously, the continuous-time version has been discussed to have a phase transition in a novel quantum universality class, while experimental investigations of a related process have pointed towards the transition being in the classical directed percolation class. We demonstrate that these differences stem from details of the discretization of the process and construct a systematic way to continuously interpolate between these results. Finally, we provide evidence that the modification of the universality class is due to destructive interference effects, showing a clear quantum-mechanical origin.
Slender structures are highly flexible, spanning several orders of magnitude in length scale. Their deformation depends on the slenderness of their cross sections, highlighting that the elasticity and geometry of structures are intrinsically coupled. The deformation of the cross-section becomes significant, particularly when tubes and pipes are subjected to bending, known as the Brazier instability. Although the bending performance of slender structures is quantified experimentally using a canonical three-point bending test, their numerical counterparts remain under-explored because complex contact mechanics must be implemented in simulations. In this study, we develop a computational framework to simulate experimental three-point bending tests using a hybrid material point method (hybrid-MPM) approach, which integrates Lagrangian finite element and Eulerian finite difference frameworks. We adapt our framework to elastic tubes and tape springs as canonical examples that exhibit characteristic bending deformation in which the cross-sectional and lengthwise bending are coupled. The predictions of numerical simulations are validated against desktop experiments and classical theory. The excellent agreement between the simulation and the experiments implies that the hybrid-MPM framework provides a robust computational framework for predicting the large deformation of structures involving complex contact, such as soft robots and deployable structures.
Community structures are common in real networks, and extracting them provides valuable insight in applications ranging from drug discovery to market segmentation. Overlapping community detection (OCD) is the task of clustering networked data in which nodes may belong to multiple clusters. Existing OCD algorithms often struggle to achieve a suitable balance between detection quality and scalability. We, therefore, propose Highway, a scalable OCD algorithm that exploits the sparse backbone of the input network to perform efficient community inference. We used 728 Lancichinetti-Fortunato-Radicchi benchmark networks to compare Highway and its ablated version against 10 existing OCD algorithms. Our results, based on five performance measures, demonstrate a competitive performance for Highway. It ranks first in overlapping normalized mutual information with a 6.9% improvement over the strongest baseline. It also ranks second in all the other four performance measures. These comparative results suggest that Highway coupled with its backbone procedure offers a suitable accuracy-efficiency trade-off. The Highway algorithm is open-source and available as part of the CDlib library.
Understanding how complex systems transition between order and chaos is a central challenge of nonequilibrium physics. While weak perturbations of classical integrable systems give rise to a mixed phase space of coexisting regular and chaotic trajectories, analogous behavior in interacting quantum many-body systems has remained elusive. Here we develop and experimentally implement a hybrid quantum-classical feedback protocol that autonomously discovers and stabilizes long-lived regular trajectories in a superconducting quantum processor. Each iteration combines short-time quantum evolution with classical optimization that projects the dynamics back onto a low-entanglement variational manifold, effectively distilling coherence from chaotic evolution. The stabilized trajectories reveal a quantum many-body mixed phase space emerging from nonlinear variational dynamics, without a direct analogue in classical or few-body quantum systems. Our results establish a versatile framework for algorithmic discovery and control of coherent dynamics previously inaccessible to experiment.
Magic and entanglement are independent quantum resources, yet their exact relation in many-body dynamics has remained elusive. We uncover two structural principles. First, at any stabilizer state, the curvature of the second stabilizer Rényi entropy under an arbitrary Hermitian generator equals the quantum Fisher information up to a fixed normalization, creating a bidirectional bridge between computational and metrological resources. Second, for commuting Ising evolution on any forest graph, a Clifford pruning circuit yields the full stabilizer-Rényi family at arbitrary size and in any spatial embedding, thereby furnishing a graph-theoretic construction of families with finite magic density and vanishing entanglement density in the thermodynamic limit. We solve two paradigmatic one-dimensional realizations central to quantum simulation -- an Ising quench and a kicked Floquet chain -- exactly for arbitrary system size and directly in the thermodynamic limit, revealing finite magic density with vanishing entanglement density, distinct magic and entanglement revival periods, and Clifford points with zero magic but finite bipartite entanglement. The same tangent geometry fixes initial growth, perturbative revival lifting, and stability of thermodynamic magic minima. Large-scale Pauli-basis matrix-product-state calculations verify all predictions, and the tangent bridge yields a concrete protocol for detecting magic through established quantum-Fisher-information measurements.
We consider inhomogeneous random graphs in which vertices are assigned i.i.d.\ random weights, pairs of distinct vertices are connected by an edge independently with a probability that is a bi-variate function of the weights of the vertices, and single vertices are connected to themselves by a self-loop independently with a probability that is a uni-variate function of the weight of the vertex. We apply a renormalisation transformation in which vertices are aggregated into groups of equal size according to a greedy algorithm, namely, distinct groups of aggregated vertices are connected by an aggregated edge if and only if there is at least one edge connecting two constituent vertices across the groups, while a group of aggregated vertices is connected to itself by an aggregated self-loop if and only if there is at least one self-loop at an internal vertex or one edge connecting a pair of distinct internal vertices. We analyse what happens when the renormalisation transformation is iterated. In particular, we show that, starting from appropriately scaled connection functions, the iterated renormalised graphs converge to a two-parameter family of random graphs, acting as an attractor in a universality class. We consider a light-tailed regime, for which the scaling limit is a homogeneous Erdős--Rényi random graph, and a heavy-tailed regime, for which the scaling limit is an inhomogeneous random graph with stable infinite-mean random weights and an exponential disconnection function. Different scalings are needed for the two regimes. Which of the two regimes prevails depends on the choice of the connection functions and the choice of the law of the random weights.