2026-03-24 | | Total: 15
Chimera states are one of the most intriguing phenomena in nonlinear dynamics, characterized by the coexistence of coherent and incoherent behavior in systems of coupled identical oscillators. Despite extensive studies and numerous observations in different settings, the development of reliable and systematic methods to classify chimera states and distinguish them from other dynamical patterns remains a challenging task. Existing approaches are often limited in scope and lack robustness. In this work, we propose a method based on Fourier analysis combined with statistical classification to characterize chimera behavior. The method is applied to a system of topological signals coupled via the Dirac operator, where it successfully captures the rich dynamical regimes exhibited by the model. We demonstrate that the proposed approach is robust with respect to variations in network topology and system parameters. Beyond the specific model considered, the framework provides a general and automated tool for distinguishing different dynamical regimes in complex systems.
Solitons are localized nonlinear wave packets that propagate without spreading because nonlinearity balances dispersion. Their robustness is well understood in effectively one-dimensional systems, but introducing additional spatial dimensions is generally expected to destabilize them or destroy their coherent character. Here we experimentally investigate how deep-water gravity-wave solitons behave when a controlled transverse degree of freedom is introduced through diffraction. Using a large-scale water-wave facility, we generate solitonic wave packets whose transverse structure is imposed across a segmented wavemaker through either a sharp slit or a smooth Gaussian apodization. The resulting two-dimensional wave fields are measured with high spatial resolution. Diffraction reshapes the transverse profile of the wave packet while its longitudinal dynamics retain the characteristic features of a soliton. Nonlinear spectral analysis confirms that the solitonic content is preserved along the direction of propagation, whereas the transverse evolution follows the linear Fresnel laws of diffraction. These observations reveal an unexpected coexistence of nonlinear soliton dynamics and classical wave diffraction.
We study transport and escape in the Stochastic Web Map (SWM), an area-preserving system with phase-space structure controlled by a symmetry parameter $q$ and nonlinearity $K$. By analyzing the survival probability $P_{\text{S}}(n)$ and escape frequency $P_{\text{E}}(\ln n)$, we show that in the chaotic regime escape dynamics is governed by a single time scale $n_{\text{typ}}\propto K^{-2}h^{2}$; here $h$ is the size of the escape horizon. Deviations at large $K$ and small $h$ indicate a breakdown of the quasilinear approximation. Then, upon rescaling the time by $n_{\text{typ}}$, escape statistics becomes universal, independent of $q$. These results demonstrate that escape is controlled by global transport rather than symmetry.
We study the integrability of a two-dimensional Hamiltonian system with a gyroscopic term and a non-homogeneous potential composed of two homogeneous components of different degrees. The model describes the motion of a particle in a plane under the combined influence of a central (Kepler-type) potential, a uniform magnetic field, and a superposition of homogeneous forces. By combining the Levi--Civita regularization with the so-called coupling constant metamorphosis transformation, and employing differential Galois theory, we derive analytical necessary conditions for integrability in the Liouville sense. They put restrictions on the degrees of homogeneity of the potential terms and their values in particular points. The obtained results encompass and generalize several classical galactic and astrophysical models, including the generalized Hill model, the Hénon--Heiles and Armbruster--Guckenheimer--Kim systems, providing a unified framework for studying non-homogeneous Hamiltonians. We demonstrate the effectiveness of the derived integrability obstructions by proving the non-integrability of these models in the presence of a uniform rotational field. The numerical analysis via the Poincaré cross-sections further confirms the analytical results, illustrating the transition from regular to chaotic dynamics as the rotational and non-homogeneous terms are introduced. Moreover, we show that, without the Kepler-type term, a generalized non-homogeneous extension of the exceptional potential remains integrable. The explicit forms of the first integrals are given.
Addressing both natural and societal challenges requires collective cooperation. Studies on collective-risk social dilemmas have shown that individual decisions are influenced by the perceived risk of collective failure. However, existing feedback evolving game models often focus on a single feedback mechanism, such as the coupling between cooperation and risk or between cooperation and cost. In many real-world scenarios, however, the level of cooperation, the cost of cooperating, and the collective risk are dynamically interlinked. Here, we present an evolutionary game model that considers the interplay of these three variables. Our analysis shows that the worst-case scenario, characterized by full defection, maximum risk, and the highest cost of cooperation, remains a stable evolutionary attractor. Nevertheless, cooperation can emerge and persist because the system also supports stable equilibria with non-zero cooperation. The system exhibits multistability, meaning that different initial conditions lead to either sustained cooperation or a tragedy of the commons. These findings highlight that initial levels of cooperation, cost, and risk collectively determine whether a population can avert a tragic outcome.
Recently, a class of spin chains known as ``free fermions in disguise'' (FFD) has been discovered, which possess hidden free-fermion spectra even though they are not solvable via the standard Jordan-Wigner transformation. In this work, we extend this FFD framework to open quantum systems governed by the Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) equation. We establish a general class of exactly solvable open quantum systems within the FFD framework: if the Liouvillian frustration graph is claw-free and has a simplicial clique, the Liouvillian possesses a hidden free-fermion spectrum. In particular, the (even-hole, claw)-free condition automatically guarantees this, enabling exact computation of the Liouvillian gap and an infinite-temperature autocorrelation function. Our results provide the first realization of the FFD mechanism in open quantum systems.
Traditional percolation theory assumes static microscopic rules, limiting its ability to describe real-world complex systems where macroscopic order actively regulates local interactions. Here, we introduce feedback percolation, an unified framework that dynamically couples the microscopic activation probability to the macroscopic size of the giant component. We show that this simple feedback mechanism produces a rich variety of behaviors both analytically and numerically. Depending on the feedback functions, the system exhibits explosive discontinuous jumps, hybrid transitions, limit-cycle oscillations, and routes to chaos, absent in classical percolation. Our findings establish that macroscopic feedback provides a unifying physical mechanism for phenomena ranging from self-regulating oscillations to systemic infrastructure collapse.
We construct the dynamical symmetry of the quantum Calogero model with particle exchange in a confining Coulomb field. This symmetry is governed by the algebra $so(N+1,2)$, deformed by exchange (Dunkl) operators, with its invariant sector generated by the Dunkl angular momentum tensor and the modified Laplace-Runge-Lenz vector. The equidistant analogue of the Hamiltonian, with a linear spectrum, is expressed in terms of the conformal subalgebra $so(1,2)$. In addition, the wave functions of the Calogero-Coulomb Hamiltonian are classified into infinite-dimensional lowest-weight $so(1,2)$ multiplets.
Cooperating first then mimicking the partner's act has been proven to be effective in utilizing reciprocity in social dilemmas. However, the extent to which this, called Tit-for-Tat strategy, should be regarded as equivalent to unconditional cooperators remains controversial. Here, we introduce a biased Tit-for-Tat (T) strategy that cooperates differently toward unconditional cooperators (C) and fellow T players through independent bias parameters. The results show that, even under strong dilemmas in the donation game framework, this three-strategy system can exhibit diverse phase diagrams on the parameter plane. In particular, when T-bias is small and C-bias is large, a ``hidden T phase'' emerges, in which the weakest T strategy dominates. The dominance of the weakened T strategy originates from a counterintuitive mechanism characterizing non-transitive ecological systems: T suppresses its relative fitness to C, rapidly eliminates the cyclic dominance clusters, and subsequently expands slowly to take over the entire population. Analysis in well-mixed populations confirms that this phenomenon arises from structured populations. Our study thus reveals the subtle role of bias regulation in cooperative modes by emphasizing the ``survival of the weakest'' effect in a broader context.
We introduce a framework for the data-driven discovery of stochastic differential equations (SDEs) that unifies, for the first time, the weak-form integration-by-parts approach of Weak SINDy with the stochastic system identification goal of stochastic SINDy. The central novelty is the adoption of spatial Gaussian test functions $K_j(x)=\exp(-|x-x_j|^2/2h^2)$ in place of temporal test functions. Because the kernel weight $K_j(X_{t_n})$ is $\mathcal{F}_{t_n}$-measurable and the Brownian innovation $ξ_n$ is independent of $\mathcal{F}_{t_n}$, every noise term in the projected response has zero conditional mean given the current state -- a property that guarantees unbiasedness in expectation and prevents the structural regression bias that afflicts temporal test functions in the stochastic setting. This design choice converts the SDE identification problem into two sparse linear systems -- one for the drift $b(x)$ and one for the diffusion tensor $a(x)$ -- that share a single design matrix and are solved jointly via $\ell_1$-regularised regression with grouped cross-validation. A two-step bias-correction procedure handles state-dependent diffusion. Validated on the Ornstein--Uhlenbeck process, the double-well Langevin system, and a multiplicative diffusion process, the method recovers all active polynomial generators with coefficient errors below 4\%, stationary-density total-variation distances below 0.01, and autocorrelation functions that faithfully reproduce true relaxation timescales across all three benchmarks.
The Out-of-Time-Order Correlator (OTOC) is a standard algebraic diagnostic of quantum information scrambling, but it offers limited direct geometric intuition. In this note, we propose a Bohmian, trajectory-based framework for constructing a geometric diagnostic of scrambling-related sensitivity using Lagrangian Descriptors (LDs). To avoid the uncertainty-principle obstruction to assigning independent initial position and momentum within a single wave function, we evaluate Bohmian dynamics over a two-dimensional preparation space of localized Gaussian wavepackets labeled by their initial center and momentum kick. For the inverted harmonic oscillator, this construction is analytically tractable: the wavepacket-center dynamics and their dependence on preparation parameters can be written explicitly. In particular, away from the equilibrium origin, the exponential growth of the associated preparation-space stability matrix yields an $\mathcal{O}(e^{ωT})$ bound on the sensitivity of the wavepacket-center LDs, motivating a semiclassical comparison with sensitivity structures associated with OTOC growth. In this sense, the LD provides a geometric indicator of scrambling-related sensitivity. We conclude by discussing how this preparation-space picture suggests a program for future work regarding the distinct microcanonical regimes previously reported for the inverted harmonic oscillator.
In classical systems, chaos is clearly defined via the behavior of trajectories. In quantum systems with a classical analogue one finds that the transition from regular to chaotic dynamics is signified by a change in the spectral statistics. This has been found to remain true for quantum systems with no classical analogue, including many-body systems. Furthermore, quantum chaotic systems explore all the allowed configurations in the Hilbert space, i.e. they are ergodic, while integrable systems, and systems in the many-body localized phase, are restricted to a certain subspace of the available phase space, and hence strongly break ergodicity. In this dissertation, we study the intermediate behavior between ergodicity and localization, i.e. the weak breaking of ergodicity. The model examined is the PXP spin chain model, where spins are allowed to flip only under certain kinetic constraints. We start by reproducing some already established results. First, we explore the eigenstate thermalization hypothesis (ETH) for this model and demonstrate the existence of a small number of states, throughout the PXP spectrum, that violate the ETH. Then we study the level-spacing statistics of the model, a well-known quantum chaos diagnostic, which turns out to be close to semi-Poisson and approach Wigner--Dyson statistics for large system sizes. Moreover, we examine various aspects of the model that have not been studied before. For example, the eigenvector component statistics, another quantum chaos diagnostic, for the PXP model turn out to be non-Gaussian. Finally, we perform a quench, in order to study how the energy spreads throughout the system, and observe ballistic fronts.
Calculating transition probabilities between different states of multistable climate tipping systems is computationally challenging in high-dimensional models. Targeted algorithms, such as the Trajectory-Adaptive Multilevel Splitting (TAMS) method, require an adequate score function to be successful, i.e., to provide an estimate of a transition probability with an acceptable variance when only a relatively small ensemble of model trajectories can be computed. Here, we present a data-driven method to derive a score function based on projecting the model dynamics in a reduced state space. Using a spatially two-dimensional partial differential equation model of the Atlantic Meridional Overturning Circulation, we show that this score function performs better than currently available ones. Using the new score function, transition probabilities can be determined with low variance, even in the case of small noise amplitudes. Besides purely noise-induced transitions, we also consider the scenario of combined stochastic and time-dependent deterministic forcing, presenting a strategy to efficiently simulate AMOC tipping events in global ocean and climate models subject to transient climate change.
Linear Response theory aims to predict how added forcing alters the statistical properties of an unforced system. These kinds of questions have been studied predominantly for autonomous dynamical systems, yet many systems in the physical, natural, and social sciences are inherently nonautonomous, evolving in time under external forcings of various kinds (a canonical example being the climate system). In such settings, one would like to understand how the system's time dependent statistical properties change when additional infinitesimal forcings are applied. This question is of clear practical relevance, but from a rigorous mathematical viewpoint it has been addressed only for a few specific classes of systems/perturbations. Here we provide a rigorous linear response theory for a rather general class of deterministic and random nonautonomous systems satisfying a specific set of assumptions that in some sense extend the standard assumptions used in the autonomous setting. A central ingredient is rapid loss of memory, i.e. sufficiently fast forgetting of initial conditions along the nonautonomous evolution. Our main strategy is to reformulate the sequential dynamics as a fixed-point problem for a global transfer operator acting on an extended sequence space of measures. This yields explicit and readily implementable response formulas for predicting the effect of small perturbations on time-dependent statistical states. We illustrate the theory on two representative classes: sequential compositions of C3 expanding maps and sequential compositions of noisy random maps, where uniform positivity of the noise induces exponential loss of memory.
Under-ice blooms in the Arctic can develop rapidly under conditions where conventional early warning signals based on critical slowing down fail due to strong noise or limited observational records. We analyze noise-induced transitions in a temperature phytoplankton stochastic differential equation model exhibiting bistability between background and bloom states. The committor function defines a stochastic separatrix as its 1/2-isocommittor, and the normal width of the associated transition layer yields a geometric indicator via arc-length averaging. Under systematic variation of noise intensity, this indicator scales linearly with noise strength, while the logarithm of the mean first passage time follows the Freidlin-Wentzell asymptotic law. Eliminating the noise parameter produces an affine scaling between the logarithmic transition time and the inverse square of the geometric indicator. The relation is robust under variations in discretization, neighborhood definition, and diffusion structure, and holds in the weak noise regime where the transition-layer width scales linearly with noise strength. Unlike variance or lag-one autocorrelation, the geometric indicator remains well defined when rapid transitions preclude reliable time-series estimation. These results provide a geometrically interpretable precursor of bloom onset that may support model-based ecological monitoring in high-variability Arctic systems.