2026-02-10 | | Total: 15
Chaos in both dissipative systems and conservative systems is investigated on the approach of renormalization group. It is found that the chaos is regarded as the critical phenomenon of equilibrium statistics in phase space. The two Feigenbaum constants in the period-doubling bifurcation systems correspond to two independent critical exponents, which are universal and can be adopted to distinguish the classes of chaos. For the conservative systems, due to the critical nature of the chaos, the isolated systems with different parameters are correlated in the phase space, and therefore the isolated system is no longer isolated in the phase space. The information of conservative systems is irreversibly lost over time, which leads to the increase entropy in an isolated system, and the contradiction between the second law of thermodynamics and the reversibility of isolated systems can be resolved.
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 discuss the classical elliptic Toda chain introduced by Krichever and the elliptic Ruijsenaars-Toda chain introduced by Adler, Shabat and Suris. It is shown that these models can be obtained as particular cases of the elliptic Ruijsenaars chain. We explain how the classical $r$-matrix structures are derived for these chains. Also, as a by-product, we prove that the elliptic Ruijsenaars-Toda chain is gauge equivalent to discrete Landau-Lifshitz model of XYZ type. The elliptic Toda chain is also gauge equivalent to XYZ chain with special values of the Casimir functions at each site.
Parrondo's paradox (PP) is a fundamental principle in nonlinear science where the alternation of individually losing strategies leads to a winning outcome. In this topical review, we provide the first systematic panorama of the synergy between PP and chaos. We observe a bidirectional connection between the two areas. The first direction is the translation of PP into the interplay between Order and Chaos through either Chaos + Chaos $\to$ Order (CCO) or Order + Order $\to$ Chaos (OOC). In this vein, many quantifiers, such as Lyapunov Exponents, $λ$, and entropic measures, are used. Second, we note that chaos can be used to engineer switching protocols that can lead to nontrivial effects in diverse PP cases. Our review clarifies the universality of PP and highlights its robust theoretical and practical applications across several areas of science and technology. Finally, we delineate key open questions, emphasizing the unresolved theoretical limits, the role of high-dimensional maps and continuous flows, and the critical need for more experimental verification of the dynamic PP in chaotic systems. For completeness, we also provide a full Python code that allows the reader to observe the many facets of the PP.
We investigate elastic, inelastic, and coalescent collisions between two-dimensional flat-top solitons supported by the cubic-quintic nonlinear Schrödinger equation. Numerical simulations reveal distinct collision regimes ranging from nearly elastic scattering to strongly inelastic interactions leading to long-lived merged states. We demonstrate that the transition between these regimes is primarily controlled by the relative phase of the solitons at the collision point, with out-of-phase collisions suppressing overlap and in-phase collisions promoting strong interaction. Kinetic-energy diagnostics are introduced to quantitatively characterize collision outcomes and to identify phase- and separation-dependent windows of elasticity. To interpret the observed dynamics, we extract effective phase-dependent interaction potentials from collision trajectories, providing a mechanical picture of attraction and repulsion between flat-top solitons. The stability of merged states formed after strongly inelastic collisions is explained by their lower energetic cost, arising from interfacial energetics, where a balance between internal pressure and edge tension plays a central role. A variational analysis based on direct energy minimization supports this picture by revealing robust energetic minima associated with stationary two-dimensional flat-top solitons.
In this paper we demonstrate that the phase space arclength of a trajectory, quantified by the time-averaged Lagrangian descriptor, is a robust and self-contained chaos indicator. By invoking Birkhoff's Ergodic Partition Theorem, we show that this scalar function distinguishes dynamical regimes via its power spectral density: for regular motion it converges to a delta function, whereas for chaotic trajectories the spectrum exhibits an inverse power-law $(1/ω)$ driven by the phenomenon of dynamical stickiness. With this approach, we avoid the computation and simulation of the variational equations and the usage of neighboring orbits, making it the simplest geometrical chaos indicator derivable from Lagrangian descriptors. Its computational efficiency enables the study of high-dimensional systems and allows the generation of large datasets of classified initial conditions, ideal for training Machine Learning models. We validate these findings using the Hénon-Heiles and the Fermi-Pasta-Ulam systems. By linking the geometrical properties of phase space to spectral analysis, this work provides the mathematical justification to establish Lagrangian descriptors as a rigorous, self-sufficient framework for the global analysis of chaos and regularity in Hamiltonian systems.
We study the dynamics of discrete breathers -- spatially localized and time-periodic solutions -- inside the bandgap of a nonlinear honeycomb lattice where the dispersion landscape approaches a so-called semi-Dirac point in which the bands cross linearly in one direction and quadratically in the orthogonal direction. By studying breather dynamics in two opposing asymptotic regimes, near the continuum and anti-continuum limits, we capture the features of hybrid coherent structures on the lattice that are highly discrete at the breather's central peak and have tails well approximated by exact separable solutions to an effective long-wave PDE theory at spatial infinity. We find that breathers are dynamically stable over a wide range of parameters and locate an instability transition. Finally, we analyze the Floquet stability of spatially extended nonlinear plane waves bifurcating from the zero solution at the edges of the gap and how they shape breather profiles inside the gap.
Localized growth in driven materials is often governed by intermittent failure, yet how a material's history biases failure sites remains poorly understood. Using pause-restart experiments on chemical precipitate membranes, we quantify the probability of age-dependent breaching. We show that the kinetics follow a stretched-exponential aging law with parameters that obey one-parameter scaling. As the system approaches a critical concentration, the stretching exponent $β$ tends to zero, signaling a crossover to scale-free, power-law behavior. A stochastic cellular automaton based on this aging rule reproduces the emergent filaments and their concentration-dependent thickening. Our findings identify aging-controlled failure with long-lived but decaying memory as a general route to pattern formation in far-from-equilibrium systems.
For three-dimensional piecewise-smooth systems of ordinary differential equations, this paper characterises the stability of points that belong to a switching surface and are equilibria of exactly one of the two neighbouring pieces of the system. Stability is challenging to characterise when nearby orbits repeatedly switch between regular motion on one side of the switching surface, and sliding motion on the switching surface, as defined via Filippov's convention. We prove that in this case stability is governed by the behaviour of a global reinjection mechanism of a four-parameter family of piecewise-linear hybrid systems, and perform a detailed numerical study of this family.
Flow matching (FM) is increasingly used for time-series generation, but it is not well understood whether it learns a general dynamical structure or simply performs an effective "trajectory replay". We study this question by deriving the velocity field targeted by the empirical FM objective on sequential data, in the limit of perfect function approximation. For the Gaussian conditional paths commonly used in practice, we show that the implied sampler is an ODE whose dynamics constitutes a nonparametric, memory-augmented continuous-time dynamical system. The optimal field admits a closed-form expression as a similarity-weighted mixture of instantaneous velocities induced by past transitions, making the dataset dependence explicit and interpretable. This perspective positions neural FM models trained by stochastic optimization as parametric surrogates of an ideal nonparametric solution. Using the structure of the optimal field, we study sampling and approximation schemes that improve the efficiency and numerical robustness of ODE-based generation. On nonlinear dynamical system benchmarks, the resulting closed-form sampler yields strong probabilistic forecasts directly from historical transitions, without training.
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.
A systematic understanding of integrability breaking in translationally invariant spin chains with genuine three-site interactions remains lacking. In this work, we introduce and classify minimal nonintegrable spin-$1/2$ Hamiltonians, defined as models that saturate injectivity while admitting no nontrivial local conserved charges beyond the Hamiltonian. We first rigorously establish the nonintegrability of the deformed Fredkin spin chain with periodic boundary conditions by mapping it to a nearest-neighbor composite-spin representation and excluding all admissible $3$-local conserved charges. Guided by its structure, we then construct five classes of spin-$1/2$ models with genuine three-site interactions. One class is integrable, while the remaining four contain exactly two interaction terms and constitute the minimal nonintegrable three-site models. Our results delineate a sharp boundary between integrability and nonintegrability beyond the nearest-neighbor paradigm.
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.
Kolmogorov-Arnold Networks (KANs) shift neural computation from linear layers to learnable nonlinear edge functions, but implementing these nonlinearities efficiently in hardware remains an open challenge. Here we introduce a physical analog KAN architecture in which edge functions are realized in materia using reconfigurable nonlinear-processing units (RNPUs): multi-terminal nanoscale silicon devices whose input-output characteristics are tuned via control voltages. By combining multiple RNPUs into an edge processor and assembling these blocks into a reconfigurable analog KAN (aKAN) architecture with integrated mixed-signal interfacing, we establish a realistic system-level hardware implementation that enables compact KAN-style regression and classification with programmable nonlinear transformations. Using experimentally calibrated RNPU models and hardware measurements, we demonstrate accurate function approximation across increasing task complexity while requiring fewer or comparable trainable parameters than multilayer perceptrons (MLPs). System-level estimates indicate an energy per inference of $\sim$250 pJ and an end-to-end inference latency of $\sim$600 ns for a representative workload, corresponding to a $\sim$10$^{2}$-10$^{3}\times$ reduction in energy accompanied by a $\sim$10$\times$ reduction in area compared to a digital fixed-point MLP at similar approximation error. These results establish RNPUs as scalable, hardware-native nonlinear computing primitives and identify analog KAN architectures as a realistic silicon-based pathway toward energy-, latency-, and footprint-efficient analog neural-network hardware, particularly for edge inference.
On the contrary to the common intuition that a steep expulsive potential makes quantum states widely delocalized, we demonstrate that one- and two-dimensional (1D and 2D) Schrödinger equations, which include expulsive potentials that are \emph{steeper than the quadratic} (anti-harmonic-oscillator) ones, give rise to \emph{normalizable} (effectively localized) eigenstates. These states constitute full continuous spectra in the 1D and 2D cases alike. In 1D, these are spatially even and odd eigenstates. The 2D states may carry any value of the vorticity (alias magnetic quantum number). Asymptotic approximations for wave functions of the 1D and 2D eigenstates, valid far from the center, are derived analytically, demonstrating excellent agreement with numerically found counterparts. Special exact solutions for vortex states are obtained in the 2D case. These findings suggest an extension of the concept of bound states in the continuum, in quantum mechanics and paraxial photonics. Gross-Pitaevskii equations are considered as the nonlinear extension of the 1D and 2D settings. In 1D, the cubic nonlinearity slightly deforms the eigenstates, maintaining their stability. On the other hand, the quintic self-focusing term, which occurs in the photonic version of the 1D model, initiates the dynamical collapse of states whose norm exceeds a critical value.