nlin updates on arXiv.org
nlin updates on the arXiv.org e-print archive.
Robustness of mKdV Breathers: A Numerical Perspective
oai:arXiv.org:2609.27080v1
arXiv:2609.27080v1 Announce Type: new
Abstract: We systematically investigate breather solutions in the modified Korteweg-de Vries (mKdV) equation via numerical simulations. We show that the breather solutions are stable under a variety of perturbations, including amplitude changes, modifications of internal parameters, and perturbations of the nonlinearity of the equation. Our results are not only consistent with the known analytical studies on the stability of breathers in integrable systems, extending them further, but also provide numerical evidence that such breather-type structures can persist and remain stable over the simulated time scales in some non-integrable settings with symmetric potentials.
Multiprecision computation of bright and dark solitons in the discrete nonlinear Schr\"odinger equation
oai:arXiv.org:2609.27663v1
arXiv:2609.27663v1 Announce Type: new
Abstract: We study the spectral stability of bright and dark solitons in the discrete nonlinear Schr\"odinger (DNLS) equation using multiprecision arithmetic. The eigenvalues governing stability are exponentially small in the lattice spacing and cannot be resolved with standard double precision. To address this, we develop a computational framework combining multiprecision arithmetic, an exact Jacobian for the stationary problem, and a squared-operator formulation for spectral analysis. This enables accurate resolution of exponentially small eigenvalues and direct comparison with exponential-asymptotic predictions. Our results show that onsite bright solitons are spectrally stable, whereas intersite bright solitons and both onsite and intersite dark solitons are unstable. Bright solitons require only a few eigenvalues and allow efficient large-scale computations, while dark solitons demand higher precision due to their proximity to the continuous spectrum. Simulations up to \(N=65{,}250\) grid points (31.7 GB RAM) highlight the necessity of multiprecision arithmetic for capturing beyond-all-orders spectral effects.
Noise-Induced Predictability Redistribution Across Forecast Horizons of Extreme Events in Chaotic Dynamics
oai:arXiv.org:2609.27877v1
arXiv:2609.27877v1 Announce Type: new
Abstract: Extreme events (EEs) in chaotic dynamics are rare broad excursions whose forecastability can be altered by dynamical noise. We investigate how noise changes EE occurrence and prediction skill across forecast horizons in a third-order autonomous chaotic flow. A single clean-data threshold is frozen for all realizations, broad events are defined by one maximum per excursion, and a future window W=15 is predicted from a 15-time-unit history using HistGradientBoosting with chronological data separation. As the forecast gap G between the observed history and the future event window increases, the clean Matthews correlation coefficient (MCC) decreases from 0.641 at G=0 to 0.165 at G=15. Noise dependence is evaluated with ten paired realizations at eight amplitudes. The mean short-horizon score increases from 0.456 in clean data to 0.546 at sigma=0.007; the paired gain is 0.0895 (95% CI 0.0494-0.1295; Holm-adjusted p=0.0234). Noise strongly increases EE occurrence while event amplitude and width remain comparatively stable. Equalizing positive training counts across noise levels substantially attenuates the short-horizon gain, whereas strong noise reduces intermediate-horizon skill. We term this horizon-dependent, nonuniform change in forecast skill noise-induced predictability redistribution (NIPR).
Heterogeneity-enhanced stochastic resonance improves liquid-state computing in delayed spiking neural networks
oai:arXiv.org:2609.27896v1
arXiv:2609.27896v1 Announce Type: new
Abstract: We investigate the joint effects of stochastic forcing and quenched structural heterogeneity on stochastic resonance (SR) and liquid-state computation in a small-world network of excitable FitzHugh--Nagumo neurons. Heterogeneity is introduced separately through coupling strengths and time delays drawn from Gaussian, bimodal, and shifted-exponential distributions. Weak periodic forcing resolves the stochastic-resonance landscape, whereas weak aperiodic driving probes noise-assisted signal encoding and Liquid State Machine (LSM) forecasting. Heterogeneity is not generically beneficial, but reorganizes the resonance landscape and computational performance in a distribution-dependent manner. Gaussian coupling disorder broadens the strong-response region, whereas bimodal disorder produces a more pronounced enhancement and shifts resonance toward weaker noise. Shifted-exponential coupling disorder behaves differently: increasing its scale broadens and shifts the coupling distribution toward less responsive regions. Time-delay heterogeneity can enhance or suppress SR depending on how the delay distribution samples the structured delay-response landscape. Under aperiodic forcing, stronger input--output coherence accompanies lower LSM prediction error, showing that SR can improve LSM performance and that suitable coupling heterogeneity can further enhance this computational benefit; Gaussian and bimodal coupling disorder shift the optimum toward weaker noise and reduce the noise-optimized root-mean-square prediction error, with the largest reduction obtained for the bimodal case. These results identify stochastic forcing and quenched heterogeneity as coupled control parameters and show that the computational enhancement of LSMs depends on disorder structure rather than magnitude alone.
General conditions for Turing and wave instabilities in reaction-diffusion systems
oai:arXiv.org:2609.28036v1
arXiv:2609.28036v1 Announce Type: new
Abstract: Necessary and sufficient conditions are provided for a diffusion-driven instability of a stable equilibrium of a reaction-diffusion system with $n$ components and a diagonal diffusion matrix. These can be either Turing or wave instabilities. Known necessary and sufficient conditions are reproduced for there to exist diffusion rates that cause a Turing bifurcation of a stable homogeneous state in the absence of diffusion. The method of proof here though, which is based on a study of dispersion relations in the contrasting limits in which the wavenumber tends to zero and to $\infty$, gives a constructive method for choosing diffusion constants. The results are illustrated on a model for the dispersion of malaria, a 3-component FitzHugh-Nagumo-like model proposed to study excitable wavetrains, and for two different coupled Brusselator systems with 4 components
Amplitude equations for wave bifurcations in reaction-diffusion systems
oai:arXiv.org:2609.28037v1
arXiv:2609.28037v1 Announce Type: new
Abstract: A wave bifurcation is the counterpart to a Turing instability in reaction-diffusion systems, but where the critical wavenumber corresponds to a pure imaginary pair rather than a zero temporal eigenvalue. Such bifurcations require at least three components and give rise to patterns that are periodic in both space and time. Depending on boundary conditions, these patterns can comprise either rotating or standing waves. Restricting to systems in one spatial dimension, complete formulae are derived for the evaluation of the coefficients of the weakly nonlinear normal form of the bifurcation up to order five, including those that determine the criticality of both rotating and standing waves. The formulae apply to arbitrary $n$-component systems ($n\geq 3$) and their evaluation is implemented in software which is made available as supplementary material. The theory is illustrated on two different versions of three-component reaction-diffusion models of excitable media that were previously shown to feature super- and subcritical wave instabilities and on a five-component model of two-layer chemical reaction. In each case, two-parameter bifurcation diagrams are produced to illustrate the connection between complex dispersion relations and different types of Hopf, Turing, and wave bifurcations, including the existence of several codimension-two bifurcations.
Computation of Turing bifurcation normal form for n-component reaction-diffusion systems
oai:arXiv.org:2609.28039v1
arXiv:2609.28039v1 Announce Type: new
Abstract: General expressions are derived for the amplitude equation valid at a Turing bifurcation of a system of reaction-diffusion equations in one spatial dimension, with an arbitrary number of components. The normal form is computed up to fifth order, which enables the detection and analysis of codimension-two points where the criticality of the bifurcation changes. The expressions are implemented within a Python package, in which the user needs to specify only expressions for the reaction kinetics and the values of diffusion constants. The code is augmented with a Mathematica routine to compute curves of Turing bifurcations in a parameter plane and automatically detect codimension-two points. The software is illustrated with examples that show the versatility of the method including a case with cross-diffusion, a higher-order scalar equation, and a four-component system.
PR-Smoother: Simulator-Preserving Non-Gaussian Smoothing for Data Assimilation
oai:arXiv.org:2609.26890v1
arXiv:2609.26890v1 Announce Type: cross
Abstract: Many physical data assimilation (DA) workflows require smoothing methods that represent non-Gaussian posteriors over physical state variables, scale to high-dimensional simulators, train from observation windows alone, and remain compatible with calibration of the prescribed simulator. We introduce PR-Smoother, a simulator-preserving amortized smoother designed for this prescribed-simulator DA regime. Its key design principle is to keep the prescribed simulator explicit in both the evidence lower bound and the variational family: rather than learning replacement dynamics or a learned trajectory prior, PR-Smoother learns only future-conditioned corrections around the prescribed rollout. This yields an explicit non-Gaussian smoothing distribution over physical trajectories and supports joint state, parameter, and sensor-bias learning from observations alone. The variational family contains the exact smoother in deterministic and linear-Gaussian limits. Empirically, PR-Smoother captures multimodal posteriors in 4-dimensional Lorenz-96, remains accurate under ambiguous nonlinear observations and process noise in 40-dimensional Lorenz-96, and scales to joint state-parameter-bias inference in 16,384-dimensional Kolmogorov flow.
A catalog of exact steady states in Kolmogorov flow
oai:arXiv.org:2609.27062v1
arXiv:2609.27062v1 Announce Type: cross
Abstract: Two-dimensional Kolmogorov flow is a common proving ground for the dynamical systems picture of turbulence. Hundreds of unstable periodic orbits have been found for this system across many Reynolds numbers and domain sizes, although the most well explored is Re$=40$ on a square $2\pi\times 2\pi$ domain with a forcing wavenumber of $k=4.$ Here I present a catalog of 8952 exact steady states of the flow in this canonical parameter regime. These solutions provide a fixed constellation informing the dynamics of both turbulence and transients.
Effects of the Internal Transport Barrier in the plasma confinement with Resonant Magnetic Perturbations
oai:arXiv.org:2609.27164v1
arXiv:2609.27164v1 Announce Type: cross
Abstract: In the pursuit of steady-state fusion energy, Advanced Tokamak regimes rely on Internal Transport Barriers (ITBs) with reversed magnetic shear. While Resonant Magnetic Perturbations (RMPs) are commonly applied to control Edge Localized Modes (ELMs) by intentionally making the plasma edge chaotic, the RMPs can also affect the reversed-shear cores. To address this compatibility, we propose a new, adjustable analytical model for plasma current that explicitly incorporates the core, ITB, and edge pedestal components. By applying this current formulation to a Hamiltonian map, we compare how easily magnetic chaos reaches the core in ITB-driven reversed-shear (Non-Twist) plasmas. We find that the localized ITB current creates "twin" inner resonances and pushes outer magnetic island chains closer to the edge. These results highlight a critical operational trade-off: while the Non-Twist topology makes ELM control easier at lower coil currents, it severely shrinks the safety margin against a complete, RMP-induced loss of global confinement.
Shift-register synchronization of dual-wavelength pulse trains in resonators with quadratic nonlinearity
oai:arXiv.org:2609.27655v1
arXiv:2609.27655v1 Announce Type: cross
Abstract: Self-organization of nonlinear waves underpins phenomena ranging from hydrodynamic pattern formation to optical frequency-comb generation. We introduce and numerically demonstrate a shift-register synchronization mechanism for pulse trains in quadratic nonlinear resonators. The synchronized state persists despite substantial group-velocity mismatch through continuous parametric energy exchange, resulting in deterministic discrete pulse shifts after each cavity round trip. Nonlinear parametric interaction modifies the effective walk-off, allowing the pulse trains to satisfy the shift-register synchronization conditions. The proposed mechanism reveals a previously unexplored synchronization regime in quadratic resonators and provides a route toward multicolour frequency combs, optical memories, and ultrafast photonic information processing.
Quantized six-vertex model at roots of unity: Frobenius property and free-parafermion spectra
oai:arXiv.org:2609.27740v1
arXiv:2609.27740v1 Announce Type: cross
Abstract: We study the quantized six-vertex (Q6V) model on arbitrary admissible diagrams at $q=\varepsilon$, where $\varepsilon$ is a primitive $N$th root of unity with $N$ odd. For the commuting layer transfer matrices $\mathbb{T}(y)$ with mixed boundary conditions, we establish the quantum Frobenius property \[ \mathbb{T}(y)\mathbb{T}(\varepsilon y)\cdots \mathbb{T}(\varepsilon^{N-1}y) = \mathscr{P}_N(y^N)\mathbb{I}, \] where $\mathscr{P}_N$ is a scalar Laurent polynomial obtained directly from the monomial expansion of $\mathbb{T}(y)$. Our proof exploits the underlying three-dimensional integrability, in particular the tetrahedron equations and their consequences. The quantum Frobenius relation essentially determines the joint eigenvalues of the associated Hamiltonians up to multiplicities and shows that they take a free-parafermion form, thereby extending such structures from previously known one-dimensional examples to a broad class of genuinely two-dimensional quantized six-vertex models. The construction also admits independent local gauge choices and recovers the long-range-interacting $\tau_{\; 2}$ model, which includes the original free-parafermion model as a special case. We also apply the same framework to the relativistic quantum Toda chain at odd roots of unity and obtain the free-parafermion spectra for the relevant Q6V transfer matrices with fixed boundary states.
Bifurcation and Multistability in a Collective-Risk Public Goods Game with Quorum-Activated Protection
oai:arXiv.org:2609.28132v1
arXiv:2609.28132v1 Announce Type: cross
Abstract: The danger of collective risk in a public goods game could provide an effective escape route from the tragedy of the commons. Existing studies usually embed risk mitigation in the contribution itself, whereas many real systems rely on a separate costly and non-excludable protective activity that requires a participation quorum. Whether such protection can evolve and persist remains unclear. Here, we formulate an $N$-player public-goods game with ordinary cooperators, defectors, and protective cooperators. Defectors endogenously increase a saturating probability of collective failure, whereas protective cooperators incur an additional cost and reduce this risk for all group members once their number reaches the quorum. Replicator analysis and numerical bifurcation calculations yield three main findings. First, protective cooperation gains an advantage when the focal participant completes the quorum and the resulting expected loss reduction exceeds its additional cost. This advantage varies non-monotonically with group composition. Second, when a single protective cooperator cannot activate protection, protection effectiveness does not enter its first-order invasion fitness at full defection; greater effectiveness cannot overcome the rare-invasion barrier when full defection is locally asymptotically stable. Third, a boundary saddle--node bifurcation creates a saddle and a locally asymptotically stable coexistence equilibrium of defectors and protective cooperators. When full cooperation and full defection are also locally asymptotically stable, three attractors coexist, and the numerically observed evolutionary outcome depends on the initial population composition. These results distinguish the effectiveness of activated protection from its ability to become established from rarity.
Optimal limits on weak integrability breaking and protected thermal memory near qutrit exchange
oai:arXiv.org:2609.28237v1
arXiv:2609.28237v1 Announce Type: cross
Abstract: Although integrability does not universally require a continuous one-site symmetry, we rigorously prove that every jointly analytic, regular Yang-Baxter deformation of the qutrit exchange interaction necessarily retains a nontrivial, analytically varying one-site charge. Breaking this local symmetry imposes a fundamental physical constraint on approximate conservation, governed by the optimal uniform bound $\delta^3 \le C\varepsilon$ that explicitly relates the minimal one-site symmetry defect $\delta$ to the local current-conservation residual $\varepsilon$. While breaking all one-site charges strictly forbids an exact integrable completion, an optimally compensated nearest-neighbor interaction saturates this cubic limit and anomalously extends the guaranteed infinite-temperature energy-current correlation window to order $|\lambda|^{-3}$ in the perturbation strength $\lambda$. Furthermore, we reveal a fundamental resonance obstruction for intrinsic conversion perturbations that strictly prevents any exact first-order repair of a broken one-site charge on any finite ring. Nevertheless, we demonstrate that the complete eight-dimensional charge memory matrix remains thermodynamically protected and approaches the identity for timescales $t=o(|\lambda|^{-3/2})$, a robust feature of the full infinite-temperature dynamics when the thermodynamic limit is taken before weak coupling.
Transfer Dynamics and Spectral Cascades in Graph-Coupled Kuramoto Networks
oai:arXiv.org:2609.28432v1
arXiv:2609.28432v1 Announce Type: cross
Abstract: This work introduces a spectral transfer formalism for nonlinear synchronisation dynamics to resolve how internal dynamical activity redistributes across network scales, moving beyond traditional global observables like the Kuramoto order parameter. By projecting Kuramoto phase dynamics onto the eigenbasis of the graph Laplacian, we construct a time-dependent transfer matrix and derive novel network metrics, including the spectral flux, which quantifies directional redistribution between structural dynamical scales rather than physical transport across nodes. Numerical simulations on modular and hierarchical networks reveal highly structured, intermittent spectral transfer episodes characterised by forward and inverse spectral cascades, directional reversals, and complex interaction geometries. Crucially, we demonstrate a clear separation between macroscopic coherence and microscopic spectral interactions, showing that highly organised, fluctuating internal redistribution dynamics persist even beneath nearly stationary and featureless global synchronisation states. This framework extends conventional synchronisation analysis and establishes a multiscale approach for interpreting network dynamics as an evolving process of spectral reorganisation.
Large-space and Large-time Asymptotics for the Focusing Nonlinear Schr\"{o}dinger Soliton Gas
oai:arXiv.org:2605.21091v2
arXiv:2605.21091v2 Announce Type: replace
Abstract: We investigate the large-space and large-time asymptotic behavior of a soliton gas for the focusing nonlinear Schr\"odinger equation. The soliton gas is constructed as the continuum limit of pure $N$-soliton solutions as $N\to\infty$, with the discrete spectrum confined to two segments $\Sigma_1$ and $\Sigma_2$. In particular, our framework does not require the discrete spectrum to be confined to the imaginary axis. By combining the nonlinear steepest descent method with an appropriate $g$-function mechanism, we show that, as $x\to-\infty$, the soliton gas is asymptotically described by a finite-gap elliptic solution with constant coefficients. In the large-time regime $t\to+\infty$, we assume that the endpoint $F$ lies on the trajectory of $H(\xi)$ with $\xi=\frac{x}{2t}\in(-E_1-\sqrt{2}E_2,-E_1)$, namely, $F=H(\hat{\xi})$, $\hat{\xi}\in (-E_1-\sqrt{2}E_2,-E_1)$. Under this assumption, we prove that the solution exhibits distinct asymptotic behaviors in different regions of the variable $\xi=\frac{x}{2t}$. More precisely, there exist an exponentially decaying region $\xi\in(-E_1,+\infty)$, a modulated elliptic-wave region $\xi\in(\hat{\xi},-E_1)$, and an unmodulated elliptic-wave region $\xi\in(-\infty,\hat{\xi})$.
On singular solitons of the KP equation and the Go-diagrams
oai:arXiv.org:2606.04663v2
arXiv:2606.04663v2 Announce Type: replace
Abstract: It has been proven that real and regular soliton solutions of the KP equation are classified in terms of the totally nonnegative Grassmannian. It is well known that vertex operators can be used to construct soliton solutions. In this paper, we consider several regular soliton solutions and study their combinations through products of vertex operators. In general, the resulting solutions become singular. Totally nonnegative elements are parametrized by the Le-diagrams introduced by Postnikov. We show that the resulting singular solutions can be parametrized by Go-diagrams, which extend Le-diagrams and arise in the Deodhar decomposition of the Grassmannian.
Volterra-Bogoyavlensky lattices and solutions of $A_{2n}^{(1)}$ invariant Painlev\'e equations
oai:arXiv.org:2606.13413v4
arXiv:2606.13413v4 Announce Type: replace
Abstract: The objective of this work is to develop a framework that exploits the lattice structure of the $k$-th Volterra--Bogoyavlensky equations ($k\in\mathbb N$, $k>1$) to generate rational solutions of higher symmetric Painlev\'e equations.
For $k=2$, we show that the Volterra lattice, equipped with suitable initial conditions, exactly models the one- and two-dimensional orbits generated by half-translation operators of the $A_2^{(1)}$ symmetric Painlev\'e IV equations. This correspondence yields explicit closed-form expressions for all solution components in terms of generalized Okamoto polynomials and leads to new algebraic recurrence relations among these polynomials.
We present two generalizations of the above Volterra lattice. One is derived from a fractional translation of the $A_{4}^{(1)}$ symmetric Painlev\'e equations. It generalizes Volterra lattice structure in the multi-component setup of the affine $A_{4}^{(1)}$ group and it is shown to generate solutions of the $A_{4}^{(1)}$ symmetric Painlev\'e equations from the seed solution invariant under dihedral group $D_{5}$. The other is the $k=3$ Bogoyavlensky lattice structure. It satisfies recurrence relations that naturally extend recurrence relations of the Volterra lattice.
These results shed light on connection between Volterra--Bogoyavlensky lattices, dihedral symmetries, and rational solutions of higher Painlev\'e systems.
Multi-kinks and composite oscillons in a commensurable and non degenerate double sine-Gordon model
oai:arXiv.org:2407.09682v3
arXiv:2407.09682v3 Announce Type: replace-cross
Abstract: In this paper, we introduce a commensurable and non-degenerate double sine-Gordon model, in which a partial breaking of vacuum degeneracy provides a mechanism for the emergence of static multi-kinks. These multi-kinks $K_n$ are stable field configurations with internal structure, consisting of $n$ localized energy packets with well-defined separations. The properties of the multi-kinks are thoroughly analyzed, including the novel phenomenology that arises during their collisions. In particular, we observe the emergence of long-lived composite oscillons that reflect the original structure of the multi-kinks. The sub-kink's positions and their vibration modes provide collective coordinates that are used to construct a phenomenological model, which offers a good qualitative explanation of the observed oscillon properties. Linear radiation effects are consistently incorporated into the model, revealing their important role in the synchronization of the oscillon's vibrational components. In contrast, oscillon decay and the sudden annihilations observed in some collisions are driven by a novel staccato-like mechanism, whereby repeated crossings of the mass threshold by the dominant frequency generate a sequence of radiation bursts.
Benchmarking Proton Tunneling Splittings with a Wavefunction-Based Double-Well Model: Application to the Formic Acid Dimer
oai:arXiv.org:2511.01903v2
arXiv:2511.01903v2 Announce Type: replace-cross
Abstract: Proton tunneling across hydrogen bonds is a fundamental quantum effect with implications for spectroscopy, catalysis, and biomolecular stability. While state-of-the-art instanton and path-integral methods provide accurate multidimensional tunneling splittings, simplified one-dimensional models remain valuable as conceptual and benchmarking tools. Here we develop a wavefunction-based framework for tunneling splittings using a Cornell-type double-well potential and apply it as a benchmark for hydrogen-bond tunneling. Analytical WKB estimates and numerical finite-difference solutions are compared across a range of barrier parameters, showing consistent agreement. As a test case, we map the formic acid dimer (FAD) barrier onto a quartic double-well model parameterized to reproduce the reported barrier height of $V_b \\approx 2848~\\text{cm}^{-1}$. The resulting tunneling splitting of about $0.037~\\text{cm}^{-1}$ matches the reduced-dimensional calculations of Qu and Bowman. The close agreement between numerical and semiclassical results highlights the pedagogical and diagnostic value of one-dimensional models, while comparison with molecular benchmarks clarifies their limitations relative to full multidimensional quantum treatments.
Quantum inversion of the vacuum hierarchy and multi-kink unbinding in the non-degenerate double sine-Gordon model
oai:arXiv.org:2609.24861v2
arXiv:2609.24861v2 Announce Type: replace-cross
Abstract: We study the quantum corrections to the vacuum energies and to the binding of the multi-kinks of the non-degenerate double sine-Gordon model. With normal ordering defined at the meson mass of the primary vacuum, the one-loop calculation predicts an inversion of the classical vacuum hierarchy, in which quantum fluctuations favour the central secondary vacuum above a coupling set by the classical splitting of the vacua and their meson mass ratio. Matrix-product-state calculations for $n=4$ confirm the inversion, and extrapolation of the bulk energy densities in the lattice spacing supports its persistence in the continuum. The crossing lies below the one-loop prediction, by $22$ percent at $\alpha=0.3$, and is reproduced closely by the Gaussian effective potential. Finite-chain calculations at $\alpha=0.3$ show that the central region of the multi-kink of topological charge $Q=4$ expands as the difference of the vacuum energies falls, as a collective-coordinate description predicts, and that beyond the crossing it grows with the size of the system. The composite therefore approaches the loss of its binding with a correction to its total mass still under a tenth of the classical value.