2026-07-28 | | Total: 12
The "butterfly effect", introduced over half a century ago by Edward Lorenz, has shifted from a cornerstone of dynamical systems to a popular metaphor, yet its true physical manifestation in fully developed turbulence spans a spectrum of phenomena from standard chaotic sensitivity to the recently established concept of Eulerian spontaneous stochasticity. This paper presents an attempt at a systematic synthesis that brings these different but interconnected ideas together within the unifying framework of the Finite Size Lyapunov Exponent (FSLE). The FSLE describes the growth rate of perturbations as a function of their scale, enabling a comprehensive characterization of the multiscale physics of turbulent flows. Using the FSLE and the Sabra shell model, extended to include thermal noise, we bridge the classical, small-scale Lyapunov regime with predictability at large scales and its interpretation in terms of Eulerian spontaneous stochasticity. Moreover, using the FSLE and the Kraichnan model, we also illustrate the closely related phenomenon of Lagrangian spontaneous stochasticity. To complete the spectrum of butterfly effects, we also examine the "literal butterfly" scenario of localized, sub-dissipative perturbations. Ultimately, this synthesis clarifies the physical mechanisms that dictate the fundamental boundaries of forecasting in high-Reynolds-number flows.
In this paper, we study an analog of the scalar two-dimensional Fermi-Pasta-Ulam (FPU) lattice. In particular, a variety of dispersive wave structures and localized patterns are numerically identified in the numerical simulations of the FPU lattice, but, to the best of our knowledge, all of these particular wave structures do not admit analytical closed-form expressions. In order to resolve this issue, we perform a dimensional reduction and accordingly derive a modified KdV equation. Based on this reduction, we first take advantage of some of its exact localized solutions to model the associated wave patterns in the FPU lattice. In addition, we explore the two-dimensional generalization of the Riemann problems for the FPU lattice and the corresponding modified KdV reduction, whose evolution dynamics lead to the formation of multiple composite dispersive structures. Moreover, we propose and rigorously derive the KPII limit of the FPU lattice and investigate their associated wedge problems. Finally, all these relevant numerical dynamics are compared to examine the performance of these quasi-continuum long-wave asymptotic limits.
Complex dynamical systems often undergo transitions from periodic to chaotic behaviour as bifurcation parameters vary, making timely detection of these changes essential. Conventional approaches based on the maximal Lyapunov exponent (MLE) generally require either knowledge of the governing equations or sufficiently long, uniformly sampled time series. Their performance degrades when the available data are scanty or contain missing observations, making reliable phase-space reconstruction difficult. We propose a methodology that combines Topological Data Analysis (TDA), specifically 0-D sublevel persistence, with Machine Learning (ML) classifiers to distinguish periodic and chaotic regimes directly from time series. Sublevel persistence extracts topological features by analysing the evolution of minima and maxima, revealing repeating signatures for periodic dynamics and more scattered patterns for chaotic dynamics. These features are used to train logistic regression, support vector machine, and k-nearest neighbour classifiers. Hyperparameters are validated using K-fold cross-validation, yielding average classification accuracies exceeding 90%. The trained classifiers provide binary predictions, identifying periodic and chaotic behaviour in previously unseen data. The proposed methodology is evaluated on the Duffing, Rössler, and Lorenz systems, where the detected transitions closely agree with those identified using the MLE, demonstrating the reliability of the approach. It is further applied to real-world ECG signals to classify normal and abnormal heartbeats, producing encouraging performance across standard statistical metrics. The proposed framework provides an effective alternative for analysing sparse or incomplete time series and is particularly useful in experimental settings where conventional nonlinear time-series methods are limited.
We develop a unified algebraic framework for nonlocal derivative nonlinear Schrödinger (DNLS)-type hierarchies based on loop algebra splittings. Within this framework, nonlocal and reverse space-time reductions are realized through algebraic constraints, leading to the corresponding integrable hierarchies. By constructing simple elements and establishing the associated factorization theory, we derive Darboux transformations (DTs). And apply DTs to construct explicit solutions.
To model high-risk traffic densities a cellular automaton model is constructed, exhibiting Ising model-like properties. The attraction-repulsion forces between vehicles are evaluated as coupling, and the coupling function $p(1-p)=g/8$, where p is the initial distribution of cells with state 1, and g is the Moore neighbor count, which has roots $\cos^2(π/8)$ and $\sin^2(π/8)$ for $g=1$. This is achieved without the use of any trigonometric functions in the code. From the roots, the tangent polynomial $\tan^2(x) + \tan(x)$ emerges. Kramers-Wannier duality is recovered and it is conjectured that the $1/\sqrt{2}$ difference between the roots serves as a fixed point for a projection mechanism from the 2-dimensional Ising model onto a 1-dimensional Ising chain through the Gudermannian function. The sigmoid evaluated at the proposed fixed point is substituted into the derivative of the logistic coupling function, $σ(1/\sqrt{2})(1-σ(1/\sqrt{2}))$, yielding a numerical approximation to the three-dimensional Ising inverse critical coupling. Finally, the results are linked to risk densities in traffic and vehicle types, accounting for the amplification of fatal accidents.
This study presents a data-driven framework for solving nonlinear wave equations, specifically the Gross-Pitaevskii equation (GPE) and the single-component nonlinear Schrodinger equation (NLSE), using Physics-Informed Neural Networks (PINNs). The approach integrates physical constraints directly into the neural network's loss function, enabling efficient training without requiring labelled data. We implement a PINN-based framework for solitons that models a variety of localized wave structures across both equations. Predicted solutions are compared with exact analytical results and show strong agreement with low error. The method effectively captures soliton profiles in both the GPE and NLSE. The accuracy and flexibility of the framework suggest its usefulness for studying nonlinear differential equations relevant to Bose--Einstein condensates and nonlinear optics.
We develop a motif-based framework for spatiotemporal chaos in spatial evolutionary games and use it to map the dynamical phase diagram in the payoff plane. Using Boolean linearization of the imitate-the-best rule, we derive analytical instability thresholds for local motifs including invaders, cooperative pairs, stripe interfaces, and cooperative cores. These thresholds are obtained from payoff balance at contested motif interfaces and recover classical invasion thresholds of spatial evolutionary games, which emerge here as boundaries of the chaotic phase. Combining the Derrida slope with the asymptotic Hamming distance, we obtain a four-region cartography: ordered, transient-chaotic, sustained-chaotic, and subcritical-chaotic dynamics. The phase diagram is organized by density-dependent motif selection: different initial cooperator densities activate different instability mechanisms, yet a small set of motif-instability lines consistently bounds the sustained-chaos region across densities. This cartography reveals a subcritical chaotic phase (Derrida slope $s<1$ but asymptotic Hamming distance $d_\infty>0$), where infinitesimal perturbations decay while finite-amplitude perturbations sustain chaos. The motif-based framework is anchored by an exact benchmark: for homogeneous backgrounds, the Boolean Jacobian yields an exact correspondence between the Derrida slope and spectral radius, linking damage spreading to deterministic instability.
We show that the inclusion of a dimension-six operator in the Higgs potential has a dramatic impact on the stability of oscillons in the $SU(2)$ bosonic sector of the Standard Model, extending their lifetime by orders of magnitude. This happens for the physical value of the ratio between the Higgs and W boson masses, $m_H/m_W=1.556$ and for the dimension-six operator $O_6 = (Φ^\dagger Φ)^3$ whose coupling constant is below the current upper bound.
Locating periodic solutions of chaotic dynamical systems normally requires an initial guess close enough to the target orbit for numerical continuation or gradient-based search to converge. We show that Physics-Informed Neural Networks (PINNs) trained on sparse, noisy observations \emph{without} initial conditions recover periodic orbits of the gravitational three-body problem, including orbit families absent from the training data. The method rests on a second-order ODE formulation, fixed-frequency Fourier features, percentile-based adaptive refinement, and a trainable scaling parameter, each validated on forward problems. Across two 100-seed ensembles, $23$--$25\%$ of runs converge to families not present in the training data. We then ask what determines which family emerges. Two $χ^2$ tests give a consistent answer: changing the training data source significantly shifts the distribution of recovered families ($p < 0.001$, Cramér's $V = 0.339$), whereas switching between the two initialization distributions tested does not ($p = 0.620$, $V = 0.094$). The random seed selects which family a given run recovers; the \emph{distribution} the weights are drawn from does not shift the aggregate frequencies, but the training data does. The evidence is empirical: we do not characterize the loss landscape analytically, and PINNs remain slower than conventional integrators on well-posed initial-value problems. What the experiments establish is that the recovered orbits are verifiable rather than merely plausible: the identified ones refine to genuine periodic solutions, a network trained on Lagrange data recovers the figure-eight choreography (Li--Liao class I.A.1, matched to seven significant digits in $T^*$), and one trained on figure-eight data recovers a Broucke--Hadjidemetriou--Hénon orbit closing to $δ_T < 10^{-9}$.
Photonic reservoir computing is a promising physical machine-learning technique for predicting time-series data. The quantization of the response signal from the reservoir is required for the implementation of photonic reservoir computing, and the number of quantization bits and sampling frequency need to be optimized to achieve high performance and low energy consumption. However, few studies have been reported to investigate the effect of bit quantization and sampling frequency. In this study, we introduce a concept of approximate reservoir computing with a semiconductor laser by quantizing the amplitude of node states in the reservoir and output weights. We evaluate the performance of a chaotic time-series prediction task and energy consumption per sample. We achieve significant reduction of energy consumption by optimizing the number of quantization bits, the sampling frequency, and the injection current of the semiconductor laser, while maintaining the prediction performance.
Photonic reservoir computing has attracted increasing attention as a fast and low-cost approach for time-series prediction. Photonic reservoir computing utilizes the high speed, broad bandwidth, and spatial parallelism of light. However, the effect of the internal connection structure (network topology) on the computing performance has not been investigated for large-scale photonic reservoirs. In this study, we experimentally and numerically demonstrate photonic reservoir computing using a spatial light modulator to systematically evaluate the relationship between the network topology and the performance of reservoir computing. We introduce complex network structures such as small-world and scale-free network topologies of the internal nodes in the reservoir. We perform the memory capacity measurement and the one-step-ahead prediction task of the chaotic time series to compare the performance. We found that the small-world network exhibits the maximum memory capacity and the best prediction performance. Our numerical calculations reveal that the performance of the time-series prediction can be optimized by changing the rewiring probability of the network and the leak rate of the reservoir. We also implement photonic human brain network as a reservoir, which is designed by the connectomes of human brain activities. We found that the network topology strongly affects the performance of reservoir computing, and the small-world network structure outperforms the other configurations.
The $N$-th order rational rogue wave of the nonlinear Schrödinger equation (NLS), which depends on $2 N-2$ real parameters, has been shown to be impossible to generate by the nonlinear superposition formula. We here generate this sequence by a three-term recurrence relation, each step only requiring the computation of three $N-1$-th order determinants of $2 N-2$ variables, while the previous method requires two determinants of order $2 N$ in $2 N$ variables. This allows us to obtain explicitly the seventh wave with its six arbitrary complex parameters. These very compact expressions open the possibility to investigate the possible existence of new patterns in addition to the already observed ones (concentric rings, polygonal configurations, \dots).