2026-08-18 | | Total: 15
We report new quasilinear and fully-nonlinear 5th-order Schwarzian evolution equations. These are symmetry-integrable evolution equations in 1+1 dimensions, i.e. equations that admit Lie-Bäcklund symmetries, whereby it is required that the equations are kept invariant under the Möbius transformation for their dependent variable.
In this paper, we study the emergence and dynamics of pure-quartic optical dispersive shock waves governed by the nonlinear Schrödinger (NLS) equation with self-defocusing Kerr nonlinearity and fourth-order dispersion. The corresponding dispersionless hydrodynamic system reveals a distinct mechanism for wave breaking, driven by the nonlinear self-steepening of the hydrodynamic velocity. We illustrate this mechanism through optical dam-break configurations described by Riemann problems and show that, in contrast to the classical quadratic NLS equation, wave breaking occurs prior to wave splitting even in the small amplitude regime. Numerical simulations of the pure-quartic NLS (PQNLS) equation and its hydrodynamic approximation confirm these qualitatively distinct dynamical regimes.
We construct novel higher-order lump solutions of the modified Kadomtsev-Petviashvili-I (mKP-I) equation by applying the generalized long-wave limit method with spectral perturbations. By tuning the phase parameters in the soliton solution, we obtain second- and third-order lump solutions, which, to the best of our knowledge, are reported here for the first time for the mKP-I equation. A detailed asymptotic analysis reveals that these degenerate lumps exhibit anomalous scattering analogous to that in the KP-I equation. However, in contrast to the KP-I case, no complete energy equipartition occurs after the collision, and the two lumps remain distinguishable in the subleading terms of their asymptotic amplitudes. This distinction highlights a difference in the interaction dynamics between the modified and standard KP-I equations.
Largest Lyapunov exponents (LLEs) quantify exponential sensitivity, but data-driven estimates are often obtained from finite-precision trajectories. We show that increasing the length of a single reduced-precision chaotic record can eventually degrade a forecast-error LLE estimate. Using the logistic map at r=4, an ESP32 single-precision trajectory is reproduced bit-for-bit by NumPy float32. Across 10,000 random float32 initial conditions, every trajectory reaches an exact recurrence before iteration 7612. For one long float32 record, the estimated LLE changes from 0.6853 at N=15,000 to 0.1827 at N=20,000 and approximately zero at N=30,000 as exact train-test histories saturate. At N=100,000, the long float32 record gives 0.0016, whereas independently restarted length-100 trajectories give 0.6917; matched float64 controls remain near ln(2)=0.6931. The collapse is reproduced for 28 representative initial conditions, and its onset is strongly correlated with the recurrence scale set by transient length and digital period (Pearson r=0.982). Thus, finite-state recurrence can turn additional samples into duplicate futures rather than new dynamical information, while independent restarts substantially delay this saturation.
In this article we consider the Drinfeld-Sokolov hierarchy of type $A$ and its similarity reduction. Then the similarity reduction is expressed as $2$ types of Hamiltonian systems, one is with a Poisson bracket and another a Poisson bracket of canonical coordinates. We also investigate a connection between the similarity reduction and the isomonodromy deformation system.
We investigate stochastic resonance in a hierarchy of delayed excitable systems, from a single FitzHugh--Nagumo neuron to single-layer and duplex networks with autaptic, pairwise, triadic higher-order, and interlayer interactions. The response to weak periodic forcing is characterized by the spectral amplitude at the driving frequency, restricting the stochastic-resonance analysis to deterministically subthreshold regimes. We show that the order and organization of interactions generate distinct resonance regimes. Among the autapse-free single-layer architectures, pairwise coupling produces the largest attainable response, whereas weak triadic coupling minimizes the noise amplitude required for resonance. Delayed autaptic feedback primarily shifts the optimal noise level and, within the configurations investigated, yields a genuine enhancement only when pairwise and triadic interactions coexist, revealing a non-additive coupling effect. In duplex networks, interlayer coupling acts predominantly as a response-equalization mechanism: it enhances the weaker layer, generally at the expense of the stronger one, with the redistribution controlled mainly by the uncoupled resonance-capacity gap. Increasing the interlayer delay generally suppresses this collective response. These results identify interaction order, resonance mismatch, and delay as key control parameters for noise-assisted signal processing in multilayer excitable systems.
We investigate integrability of Euler-Lagrange equations associated with 3D second-order Lagrangians of the form \begin{equation*} \int f(u_{xy},u_{xt},u_{yt})\ \text{d}x\text{d}y\text{d}t. \end{equation*} It is demonstrated that there are exactly four different types of such Lagrangian densities $f$: the first one is given by the formula $f=\sqrt{u_{xy}u_{xt}u_{yt}}$, the second and the third are more complicated (although still representable in elementary functions), whereas the most generic fourth one is expressible in terms of the Lobachevsky function, revealing unexpected links to spherical/hyperbolic trigonometry and Schläfly-type formulas. Dispersionless Lax pairs and integrable dispersive deformations of the corresponding Euler-Lagrange equations are also constructed.Remarkably, dispersive deformation of the first Lagrangian density coincides with the Lagrangian of the classical Darboux system arising in the theory of triply-orthogonal coordinate systems in $\mathbb{R}^3$. Dispersive deformations of the three other cases provide Lagrangian formulation of semi-discrete and fully discrete versions of the Darboux system, with one, two and three discrete variables, respectively.
A scuba diver holding constant depth balances on an unstable equilibrium: the gas carried in the suit and buoyancy compensator compresses with depth, so the buoyant force falls as the diver sinks and rises as the diver ascends. We represent the diver as a proportional-derivative controller that regulates this compressible-buoyancy saddle after a finite reaction delay, and we derive the governing delay differential equation from the vertical force balance and the isothermal gas law, reducing it to a damping ratio, two control gains, and a dimensionless delay. The characteristic spectrum, obtained by pseudospectral collocation of the semigroup generator and checked against a direct Newton solution of the characteristic equation, locates the Hopf boundary that separates stable hovering from sustained porpoising; for the baseline diver the critical reaction delay is 3.36 s and the onset period is 28.8 s. The bifurcation is supercritical, and because the saturating force is the quadratic hydrodynamic drag, the limit-cycle amplitude grows in proportion to the delay excess rather than as its square root. The safe-operating envelope shows that runaway ascent is triggered by saturation of the compensator, not by loss of linear stability, so a stable and an unstable diver can share the same escape threshold. As onset is approached, the lag-one autocorrelation and variance rise while the fitted recovery rate falls and matches the spectral abscissa, giving an eigenvalue-exact early warning of the transition.
New types of soliton and periodic waves are identified for a nonlinear dissipative medium where the pulse propagation is governed by the cubic complex Ginzburg-Landau equation. We find that the dynamical equation for the pulse amplitude supports two distinct types of kink and antikink solitons with different functional forms. It is found that the obtained kink and antikink soliton waveforms occur under the same fixed inverse velocity. The results also indicate that the periodic waves can propagate with variety of wave forms such as sn, cn, dn and their rational forms as well. It is also shown that in the long-wave limit, the derived periodic waves degenerate into different bright and dark soliton pulses. The stability analysis based on the theory of dispersive waves in nonlinear optics is developed. It is shown that some elliptic and soliton solutions are quasi-stable in the context of passive mode locking lasers described by complex Ginzburg-Landau equation.
Abrupt transitions in complex systems are often preceded by early warning signals. However, most indicators rely on the notion of critical slowing down and do not generally extend to rate-induced tipping where transitions can occur without local loss of stability. This is problematic in stochastic, nonautonomous systems where internal variability and time-varying variables interact to shape tipping onset. We use Koopman operator theory to develop a unified early warning framework for both bifurcation and rate-induced tipping in stochastic systems. Our approach builds on residual Koopman mode decomposition that measures discrepancies between dynamics and their finite-dimensional approximation, and extends it to the control setting by augmenting the observable space with time-varying control variables. In idealized examples, the resulting indicators recover expected signatures near bifurcation points and improve detection in rate-induced regimes where classical indicators fail. We further show that learned embeddings through deep learning outperform prescribed dictionaries, especially in a high-dimensional setting. Applied to simulations of the Atlantic Meridional Overturning Circulation, our Koopman-based indicators distinguish tipping from non-tipping trajectories and reveal interpretable spectral signatures prior to critical transition.
Neural autoregressive models have rapidly emerged as powerful emulators of high-dimensional chaotic systems, yet their long-term instability and error growth remain poorly understood, leading to ad-hoc solutions. Here, we develop an eigenanalysis framework that reveals the dynamical origin of this error growth. By analyzing the Jacobian of the learned one-step update map with respect to the state, we show how inference-time error growth, and thus model stability, is governed by its spectral radius. Direct-step architectures (models that predict the next state from the previous one) generically admit unstable eigenvalues with magnitudes exceeding one, explaining the rapid divergence of these widely used models. In contrast, integration-constrained models (where the time derivative is estimated and integrated with a higher-order integrator) collapse their eigenspectrum onto the unit circle, yielding neutral stability and a universal linear error-scaling law. The largest eigenvalue of this Jacobian provides an architecture-agnostic, a priori diagnostic of short-term skill, long-term stability, and spectral bias, without requiring an expensive rollout. Leveraging this theory, we introduce a stability-promoting loss that explicitly regularizes Jacobian-driven error amplification, improving both forecast accuracy and dynamical robustness. Demonstrated across $29$ models spanning two architectures, several explicit and implicit integrators, and multiple loss functions on the Kuramoto-Sivashinsky system, our results establish a theoretical foundation for the design and evaluation of neural emulators of chaotic multi-scale dynamics. More broadly, our framework is a step toward the kind of a priori stability analysis that numerical analysis provides for discretizations of differential equations and that scientific machine learning currently lacks.
Adjoint methods provide a powerful route for gradient-based optimization, but their physical implementation is obstructed in generic nonlinear systems because the adjoint dynamics requires backward-time evolution, Jacobian transposition, and terminal-value constraints. Here we show that nonlinear parity-time ($\mathcal{PT}$)-symmetric systems overcome this obstruction. Using a class of nonlinear non-Hermitian resonator networks, we establish symmetry relations that map the formal adjoint dynamics onto experimentally accessible forward-time evolutions supplemented by controlled injections. This construction enables exact in-situ evaluation of adjoint gradients without requiring explicit backward propagation or matrix transposition. We demonstrate the approach in nonlinear $\mathcal{PT}$-symmetric resonator chains, where the resulting optimization protocol autonomously discovers parameter configurations that realize prescribed spatio-temporal functionalities, including uniform energy redistribution and targeted wave transport at predefined time windows. Our results identify $\mathcal{PT}$ symmetry as a resource for implementing computational sensitivities within physical systems and establish a route toward self-optimizing nonlinear machines.
We compute the dynamic critical exponent $z$ of the Yang--Lee edge singularity in relaxational dynamics using perturbative renormalization-group methods to the three-loop order. The calculation combines diagram reduction, analytic two-loop evaluation via parametric integration with hyperlogarithms, and numerical evaluation of three-loop integrals using the Sector Decomposition method. The perturbative series obtained is resummed using Padé and Padé-Borel-Leroy methods to produce estimates of $z$ in various spatial dimensions. The resulting values are consistent with previous perturbative and functional renormalization-group calculations.
Spatial early warning signals (EWSs) seek evidence of an approaching tipping point from a single snapshot of many interacting elements. Existing theory largely assumes spatial homogeneity, whereas networks introduce systematic differences among nodes that may obscure fluctuation-based warning signals. We develop a mathematical framework for spatial EWSs in stochastic dynamical systems on networks. We find that the expected spatial variance, a popular spatial EWS, decomposes exactly into a structural contribution from heterogeneity in the equilibrium state and a fluctuation contribution determined by the stationary covariance. Near a simple steady-state bifurcation, the potentially divergent covariance concentrates along the critical eigendirection: the left eigenvector determines how strongly noise excites the critical fluctuation, while the right eigenvector determines its spatial pattern. Consequently, the spatial variance has a divergent fluctuation contribution when the limiting critical eigendirection is noise-excited and spatially nonuniform after centering. In contrast, the spatial coefficient of variation generally saturates, while skewness, kurtosis, and Moran's $I$ approach network-dependent limits without a universal warning direction. We also derive results for homogeneous networks, node-wise baseline subtraction as preprocessing, and Hopf bifurcations, for which the limiting distributions are qualitatively different. These results clarify when spatial EWSs provide reliable warnings and why their performance depends on network structure, noise, and preprocessing.
In the late 1980s, Baulieu and Grossman demonstrated that the supersymmetric formulation of Langevin stochastic differential equations (SDEs), proposed earlier by Parisi and Sourlas, belongs to the family of Witten-type topological field theories (TFTs). From a certain angle, this finding may appear puzzling: TFTs have no local degrees of freedom, whereas SDEs do exhibit local fluctuations. In this paper, we address this apparent contradiction in the context of the supersymmetric theory of stochastic dynamics (STS), a generalization of the Parisi-Sourlas-Baulieu-Grossman approach to SDEs of arbitrary form. We further discuss how, within STS, dynamical chaos can be identified with the spontaneous breakdown of the corresponding topological supersymmetry, while 1/f noise can be understood as a consequence of the Goldstone theorem. The butterfly effect, in turn, calls for an effective field-theoretic description that may itself possess a hidden topological structure, as we speculate on the basis of the Ginzburg-Landau approach and insights from AdS/CFT duality.