2025-04-08 | | Total: 27
Stabilizing thermodynamically unstable phases in many-body systems -- such as suppressing pathological neuronal synchronization in Parkinson's disease or maintaining magnetic order across broad temperature ranges -- remains a persistent challenge. In traditional approaches, such phases are stabilized through intervening in the dynamics of all system constituents or introducing additional interactions. Here, we offer a hitherto-unexplored alternative -- subsystem resetting, whereby intervention in the dynamics of only a part of the system, and that too only occasionally in time, is implemented through resetting its state to a reset configuration. Just playing with a few parameters, e.g., the nature of the reset configuration and the size of the reset subsystem, one achieves a remarkable and robust control over the phase diagram of the bare dynamics. We demonstrate that these universal effects span a wide variety of scenarios, including equilibrium and non-equilibrium, mean-field and non-mean-field dynamics, with and without quenched disorder. Despite the challenges posed by memory effects, we obtain explicit analytical predictions, validated by simulations.
Stochastic storage models based on essentially non-Gaussian noise are considered. The stochastic description of physical systems based on stochastic storage models is associated with generalized Poisson (or shot) noise, in which the jump values can be quite large. Stochastic storage models have a direct physical meaning: some elements enter the system and leave it. Storage processes fit into the general scheme of dynamic systems subject to the additive influence of a random process. The main relationships of storage models are described, and the possibilities of applying the mathematical provisions of stochastic storage processes to various physical problems are indicated. A number of examples of applying the stochastic storage model are considered.
We consider the finite-time quench dynamics in the quantum transverse field Ising model which exhibits a second order phase transition from a paramagnetic to a ferromagnetic phase, as the transverse magnetic field is decreased. These dynamics have been thoroughly investigated in previous studies when the critical point is crossed during the quench; here, we quench the system from deep in the paramagnetic phase to just above the critical field so that the system remains in the gapped phase throughout the quench duration. On linearly quenching the infinitely large system, we find that the behavior of mean longitudinal defect density and mean transverse magnetization at the end of the quench falls into three distinct scaling regimes as the quench time is increased. For sufficiently small quench times, these observables remain roughly constant, but for larger quench times, a crossover occurs from the Kibble-Zurek scaling law to the quadratic quench rate law when the Kibble-Zurek time is of the order of relaxation time at the final quench field. These results are shown analytically using power series and uniform asymptotic expansions of the exact solution of the model, and also compared with an adiabatic perturbation theory in the third regime. We find that the above mentioned scaling regimes hold for quenches within the ferromagnetic phase also, and provide a general scaling argument for crossover from the Kibble-Zurek regime to an adiabatic regime for slow quenches within a gapped phase.
The observation of the Fermi-Pasta-Ulam-Tsingou (FPUT) paradox, namely the lack of equipartition in the evolution of a normal mode in a nonlinear chain on unexpectedly long times, is arguably the most famous numerical experiment in the history of physics. Since seventy years after its publication, most studies in FPUT chains still focus on long wavelength initial states similar to the original paper. It is shown here that all characteristic features of the FPUT paradox are rendered even more striking if short(er) wavelength modes are evolved instead. Since not every normal mode leads to equipartition, we also provide a simple technique to predict which modes, and in what order, are excited starting from an initial mode (root) in α-FPUT chains. The excitation sequences associated with a root are then shown to spread energy at different speeds, leading to prethermalization regimes that become longer as a function of mode excitation number. This effect is visible in observables such as mode energies and spectral entropies and, surprisingly, also in the time evolution of invariant quantities such as Lyapunov times and Kolmogorov-Sinai entropies. Our findings generalize the original FPUT experiment, provide an original look at the paradox's source, and enrich the vast literature dedicated to studying equipartition in classical many-body systems.
We investigate the response of a system of hard spheres to two classes of perturbations over a range of densities spanning the fluid, crystalline, and glassy regimes within a molecular dynamics framework. Firstly, we consider the relaxation of a "thermal inhomogeneity," in which a central region of particles is given a higher temperature than its surroundings and is then allowed to evolve under Newtonian dynamics. In this case, the hot central "core" of particles expands and collides with the cold surrounding material, creating a transient radially-expanding "compression wave," which is rapidly dissipated by particle-particle collisions and interaction with periodic images at the boundary, leading to a rapid relaxation to equilibrium. Secondly, we consider a rapid compression of the spheres into a disordered glassy state at high densities. Such rapidly compressed systems exhibit very slow structural relaxation times, many orders of magnitude longer than thermalization times for simple temperature inhomogeneities. We find that thermal relaxation of the velocity distribution is determined simply by the total collision rate, whereas structural relaxation requires coordinated collective motion, which is strongly suppressed at high density, although some particle rearrangement nevertheless occurs. We further find that collisions propagate significantly faster through glassy systems than through crystalline systems at the same density, which leads to very rapid relaxation of velocity perturbations, although structural relaxation remains very slow. These results extend the validity of previous observations that glassy systems exhibit a hybrid character, sharing features with both equilibrium and non-equilibrium systems. Finally, we introduce the hard sphere causal graph, a network-based characterization of the dynamical history of a hard sphere system, which encapsulates several useful...
We study stationary fluctuations of conserved slow modes in a two-lane model of hardcore particles which are expected to show universal behaviour. Specifically, we focus on the properties of fluctuations at a special umbilic point where the characteristic velocities coincide. At large space and time scales, fluctuations are described by a system of stochastic Burgers equations studied recently in [13]. Our data suggest coupling-dependent scaling functions and, even more surprisingly, coupling-dependent dynamical scaling exponents, distinct from KPZ scaling exponent typical for surface growth processes.
Supervised machine learning (ML) methods are emerging as valid alternatives to standard mathematical methods for identifying knots in long, collapsed polymers. Here, we introduce a hybrid supervised/unsupervised ML approach for knot classification based on a variational autoencoder enhanced with a knot type classifier (VAEC). The neat organization of knots in its latent representation suggests that the VAEC, only based on an arbitrary labeling of three-dimensional configurations, has grasped complex topological concepts such as chirality, unknotting number, braid index, and the grouping in families such as achiral, torus, and twist knots. The understanding of topological concepts is confirmed by the ability of the VAEC to distinguish the chirality of knots 9_{42} and 10_{71} not used for its training and with a notoriously undetected chirality to standard tools. The well-organized latent space is also key for generating configurations with the decoder that reliably preserves the topology of the input ones. Our findings demonstrate the ability of a hybrid supervised-generative ML algorithm to capture different topological features of entangled filaments and to exploit this knowledge to faithfully reconstruct or produce new knotted configurations without simulations.
Studied are coupled KPZ equations with three control parameters X,Y,T. These equations are used in the context of stretched polymers in a random medium, for the spacetime spin-spin correlator of the isotropic quantum Heisenberg chain, and for exciton-polariton condensates. In an earlier article we investigated merely the diagonal X=Y, T=1. Then the stationary measure is delta-correlated Gaussian and the dynamical exponent equals z = \tfrac{3}{2}. We observed that the scaling functions of the dynamic correlator change smoothly when varying X. In this contribution, the analysis is extended to the whole X-Y-T plane. Solutions are stable only if XY \geq 0. Based on numerical simulations, the static correlator still has rapid decay. We argue that the parameter space is foliated into distinct universality classes. They are labeled by X and consist of half-planes parallel to the Y-T plane containing the point (X,X,1).
In this paper we use singular perturbation theory to solve the 2D narrow capture problem for a set of partially accessible targets \calU_k, k=1,\ldots,N, in a bounded domain \Omega\subset \R^2. In contrast to previous models of narrow capture, we assume that when a searcher finds a target by attaching to the partially adsorbing surface \partial \calU_k it does not have immediate access to the resources within the target interior. Instead, the searcher remains attached to the surface for a random waiting time \tau, after which it either gains access to the resources within ({\em surface absorption}) or detaches and continues its search process ({\em surface desorption)}. We also consider two distinct desorption scenarios -- either the particle continues its search from the point of desorption or rapidly returns to its initial search position. We formulate the narrow capture problem in terms of a set of renewal equations that relate the probability density and target flux densities for absorption to the corresponding quantities for irreversible adsorption. The renewal equations, which effectively sew together successive rounds of adsorption and desorption prior to the final absorption event, provide a general probabilistic framework for incorporating non-Markovian models of desorption/absorption and different search scenarios following desorption. We solve the general renewal equations in two stages. First, we calculate the Laplace transformed target fluxes for irreversible adsorption by solving a Robbin boundary value problem (BVP) in the small-target limit using matched asymptotic analysis. We then use the inner solution of the BVP to solve the corresponding Laplace transformed renewal equations for non-Markovian desorption/absorption, which leads to explicit Neumann series expansions of the corresponding target fluxes.
All digital devices have components that implement Boolean functions, mapping that component's input to its output. However, any fixed Boolean function can be implemented by an infinite number of circuits, all of which vary in their resource costs. This has given rise to the field of circuit complexity theory, which studies the minimal resource cost to implement a given Boolean function with any circuit. Traditionally, circuit complexity theory has focused on the resource costs of a circuit's size (its number of gates) and depth (the longest path length from the circuit's input to its output). In this paper, we extend circuit complexity theory by investigating the minimal thermodynamic cost of a circuit's operation. We do this by using the mismatch cost of a given circuit that is run multiple times in a row to calculate a lower bound on the entropy production incurred in each such run of the circuit. Specifically, we derive conditions for mismatch cost to be proportional to the size of a circuit, and conditions for them to diverge. We also use our results to compare the thermodynamic costs of different circuit families implementing the same family of Boolean functions. In addition, we analyze how heterogeneity in the underlying physical processes implementing the gates in a circuit influences the minimal thermodynamic cost of the overall circuit. These and other results of ours lay the foundation for extending circuit complexity theory to include mismatch cost as a resource cost.
Spin systems exposed to the influence of random magnetic fields are paradigmatic examples for studying the effect of quenched disorder on condensed-matter systems. In this context, previous studies have almost exclusively focused on systems with Ising or continuous symmetries, while the Potts symmetry, albeit being of fundamental importance also for the description of realistic physical systems, has received very little attention. In the present study, we use a recently developed quasi-exact method for determining ground states in the random-field Potts model to study the problem with four states. Extending the protocol applied for the three-state model, we use extensive finite-size scaling analyses of the magnetization, Binder parameter, energy cumulant, specific heat, and the connected as well as disconnected susceptibilities to study the magnetic ordering transition of the model. In contrast to the system in the absence of disorder, we find compelling evidence for a continuous transition, and we precisely determine the critical point as well as the critical exponents, which are found to differ from the exponents of the three-state system as well as from those of the random-field Ising model.
We extend the chemical-potential-based free-energy lattice Boltzmann (LB) model of Li et al. [Phys. Rev. E 103, 013304 (2021)] by integrating generalized equilibria, originally formulated for the color-gradient LB model using sixth-order Hermite polynomials [Saito et al., Phys Rev E 108, 065305 (2023)], into a thermodynamically consistent framework. Our model is formulated on a three-dimensional D3Q27 lattice with a central-moment collision scheme, simplifying implementation and improving Galilean invariance. Numerical tests, including flat-interface equilibrium, stationary and moving droplets in free space, and wetting on solid surfaces, confirm the model's capability to accurately simulate multiphase phenomena while maintaining strict thermodynamic consistency.
We investigate non-local conserved charges in perturbed two-dimensional conformal field theories from the point of view of the 3d SymTFT of the unperturbed theory. In the SymTFT we state a simple commutation condition which results in a pair of compatible bulk and defect perturbations, such that the perturbed line defects are conserved in the perturbed CFT. In other words, the perturbed defects are rigidly translation invariant, and such defects form a monoidal category which extends the topological symmetries. As examples we study the A-type Virasoro minimal models M(p,q). Our formalism provides one-parameter families of commuting non-local conserved charges for perturbations by a primary bulk field with Kac label (1,2), (1,3), or (1,5), which are the standard integrable perturbations of minimal models. We find solutions to the commutation condition also for other bulk perturbations, such as (1,7), and we contrast this with the existence of local conserved charges. There has been recent interest in the possibility that in certain cases perturbations by fields such as (1,7) can be integrable, and our construction provides a new way in which integrability can be found without the need for local conserved charges.
We present efficient and practical protocols to measure the second Rényi entropy (RE), whose exponential is known as the purity. We achieve this by establishing a direct connection to a Loschmidt echo (LE) type measurement sequence, applicable to quantum many-body systems. Notably, our approach does not rely on random-noise averaging, a feature that can be extended to protocols to measure out-of-time-order correlation functions (OTOCs), as we demonstrate. By way of example, we show that our protocols can be practically implemented in superconducting qubit-based platforms, as well as in cavity-QED trapped ultra-cold gases.
Martensites subjected to quasistatic deformation are known to exhibit power law distributed acoustic emission in a broad range of scales. However, the origin of the observed scaling behavior and the mechanism of self organization towards criticality remains obscure. Here we argue that the power law structure of the fluctuations spectrum can be interpreted as an effect of inertia. The general insight is that inertial dynamics can become a crucial player when the underlying mechanical system is only marginally stable. We first illustrate the possibility of inertia-induced criticality using an elementary example of mass points connected by bi-stable springs. We then explore the effects of inertia in the fully realistic two and three dimensional continuum models of specific elastic phase transitions in crystals.
Links and knots are exotic topological structures that have garnered significant interest across multiple branches of natural sciences. Coherent links and knots, such as those constructed by phase or polarization singularities of coherent light, have been observed in various three-dimensional optical settings. However, incoherent links and knots - knotted or connected lines of coherence singularities - arise from a fundamentally different concept. They are hidden in the statistic properties of a randomly fluctuating field, making their presence often elusive or undetectable. Here, we theoretically construct and experimentally demonstrate such topological entities of incoherent light. By leveraging a state-of-the-art incoherent modal-decomposition scheme, we unveil incoherent topological structures from fluctuating light speckles, including Hopf links and Trefoil knots of coherence singularities that are robust against coherence and intensity fluctuations. Our work is applicable to diverse wave systems where incoherence or practical coherence is prevalent, and may pave the way for design and implementation of statistically-shaped topological structures for various applications such as high-dimensional optical information encoding and optical communications.
Stabilizing cooperation among self-interested individuals presents a fundamental challenge in evolutionary theory and social science. While classical models predict the dominance of defection in social dilemmas, empirical and theoretical studies have identified various mechanisms that promote cooperation, including kin selection, reciprocity, and spatial structure. In this work, we investigate the role of localized imitation in the evolutionary dynamics of cooperation within an optional Public Goods Game (PGG). We introduce a model where individuals belong to distinct groups and adapt their strategies based solely on comparisons within their own group. We identify different dynamical regimes, including stable fixed points, limit cycles, and Rock-Scissors-Paper-type oscillations. Our analysis, grounded in a replicator-type framework, reveals that such group-level imitation can stabilize cooperative behavior, provided that groups are not initially polarized around a single strategy. In other words, restricting imitation to group-level interactions mitigates the destabilizing effects of global competition, providing a potential explanation for the resilience of cooperation in structured populations.
Topological gapless phases of matter have been a recent interest among theoretical and experimental condensed matter physicists. Fermionic chains with extended nearest neighbor couplings have been observed to show unique topological transition at the multicritical points between distinct gapless phases. In this work, we show that such topological gapless phases and the transition between them can be simulated in a quantum walk. We consider a three-step discrete-time quantum walk and identify various critical or gapless phases and multicriticalities from the topological phase diagram along with their distinguished energy dispersions. We reconstruct the scaling theory based on the curvature function to study transition between gapless phases in the quantum walk. We show the interesting features observed in fermionic chains, such as diverging, sign flipping and swapping properties of curvature function, can be simulated in the quantum walk. Moreover, the renormalization group flow and Wannier state correlation functions also identify transition at the multicritical points between gapless phases. We observe the scaling law and overlapping of critical and fixed point properties at the multicritical points of the fermionic chains can also be observed in the quantum walk. Furthermore, we categorize the topological transitions at various multicritical points using the group velocity of the energy eigenstates. Finally, the topological characters of various gapless phases are captured using winding number which allows one to distinguish various gapless phases and also show the transitions at the multicritical points.
We study the fluctuations of generic currents in multi-terminal, multi-channel quantum transport settings. In the quantum regime, these fluctuations and the resulting precision differ strongly depending on whether the device is of fermionic or bosonic nature. Using scattering theory, we show that the precision is bounded by constraints set by the entropy production and by the activity in the spirit of thermodynamic or kinetic uncertainty relations, valid for fermionic and bosonic quantum systems and even in the absence of time-reversal symmetry. Furthermore, we derive a combined thermodynamic kinetic uncertainty relation, which is tight over a wide range of parameters and can hence predict the reachable precision of a device. Since these constraints can be expressed in terms of observables accessible in transport measurements, such as currents and bandwidth, we foresee that the tight thermodynamic kinetic uncertainty-like bounds are also useful as an inference tool: they can be exploited to estimate entropy production from transport observables, such as the charge current and its noise, which are more easily accessible in experiment.
A system of exclusive fermions occurs when two fermions of opposite spin are prohibited from occupying the same quantum level. We derive the distribution of exclusive fermions via the employment of the grand canonical ensemble. Salient features of its statistical properties, compared to the free electron gases, include: larger Fermi energy, higher degeneracy pressure, but the same Pauli paramagnetism and Landau diamagnetism. In particular, higher degeneracy pressure leads to an inflation of the Chandrasekhar limit to 1.6 times when applied to white dwarf stars and neutron stars.
The importance of Fisher information is increasing in nonequilibrium thermodynamics, as it has played a fundamental role in trade-off relations such as thermodynamic uncertainty relations and speed limits. In this work, we investigate temporal Fisher information, which measures the temporal information content encoded in probability distributions, for both classical and quantum systems. We establish that temporal Fisher information is bounded from above by physical costs, such as entropy production in classical Langevin and Markov processes, and the variance of interaction Hamiltonians in open quantum systems. Conversely, temporal Fisher information is bounded from below by statistical distances (e.g., the Bhattacharyya arccos distance), leading to classical and quantum speed limits that constrain the minimal time required for state transformations. Our work provides a unified perspective of speed limits from the point of view of temporal Fisher information in both classical and quantum dynamics.
In this article, we investigate the impact of self-alignment and anti-self-alignment on collective phenomena in dense active matter. These mechanisms correspond to effective torques that align or anti-align a particles orientation with its velocity, as observed in active granular systems. In the context of motility-induced phase separation (MIPS) - a non-equilibrium coexistence between a dense clustered phase and a dilute homogeneous phase - both self- and anti-self-alignment are found to suppress clustering. Specifically, increasing self-alignment strength first leads to flocking within the dense cluster, and eventually to the emergence of a homogeneous flocking phase. In contrast, anti-self-alignment induces a freezing phenomenon, progressively reducing particle speed until MIPS is suppressed and a homogeneous phase is recovered. These results are supported by scaling arguments and are amenable to experimental verification in high-density active granular systems exhibiting self- or anti-self-alignment.
Spin-glass models of associative memories are a cornerstone between statistical physics and theoretical neuroscience. In these networks, stochastic spin-like units interact through a synaptic matrix shaped by local Hebbian learning. In absence of self-interactions (i.e., autapses), the free energy reveals catastrophic forgetting of all stored patterns when their number exceeds a critical memory load. Here, we bridge the gap with biology by considering networks of deterministic, graded units coupled via the same Amari-Hopfield synaptic matrix, while retaining autapses. Contrary to the assumption that self-couplings play a negligible role, we demonstrate that they qualitatively reshape the energy landscape, confining the recurrent dynamics to the subspace hosting the stored patterns. This allows for the derivation of an exact overlap-dependent Lyapunov function, valid even for networks with finite size. Moreover, self-interactions generate an auxiliary internal field aligned with the target memory pattern, widening the repertoire of accessible attractor states. Consequently, pure recall states act as robust associative memories for any memory load, beyond the critical threshold for catastrophic forgetting observed in spin-glass models -- all without requiring nonlocal learning prescriptions or significant reshaping of the Hebbian synaptic matrix.
Dynamic processes on networks are fundamental to understand modern day phenomena such as information diffusion and opinion polarization on the internet or epidemic spreading. However, they are notoriously difficult to study broadly as small changes in initial conditions, the process or the network can lead to very different evolution trajectories. Here we apply the information-theoretic framework of maximum caliber to study the statistics of such systems analytically in a general way. We verify the dynamics deduced from maximum caliber by using simulations of different processes on different networks, introduce an approximation of the dynamics which significantly simplifies the problem and show that the approximation can be used to recover well-established models of population dynamics that are typically not thought of as taking place on a network. This provides a theoretical tool that allows to study dynamic processes on networks broadly without loosing sight of known results.
We investigate the spin-1/2 Heisenberg model with competing ferromagnetic and antiferromagnetic interactions on diamond-decorated lattices. Tuning the exchange interactions to the boundary of the ferromagnetic phase, we analyze the models with two types of diamond units: distorted and ideal diamonds. In the distorted diamond model, flat bands in the magnon spectra indicate the localized states confined to small regions (`trapping cells') of the lattice. Remarkably, these trapping cells can host up to five and seven localized states for square and cubic lattices, respectively, leading to the macroscopic ground state degeneracy and high value of residual entropy. The problem of calculating ground state degeneracy reduces to that of non-interacting spins, whose spin value equal to half the number of localized magnons in the trapping cell. In contrast, ideal diamond models feature ground states composed of randomly distributed isolated diamond diagonal singlets immersed in a ferromagnetic background. Counting the ground state degeneracies here maps onto the percolation problem in 2D and 3D lattices. Our analysis shows that ideal diamond models possess even greater ground state degeneracy than their distorted counterparts. These findings suggest that synthesizing diamond-decorated-type compounds holds great promise for low-temperature cooling applications.