Nonlinear Sciences
See recent articles
Showing new listings for Friday, 1 August 2025
- [1] arXiv:2507.23406 [pdf, html, other]
-
Title: Theory of dispersive shock waves induced by the Raman effect in optical fibersComments: 8 pages, 4 figuresSubjects: Pattern Formation and Solitons (nlin.PS); Optics (physics.optics)
We develop the theory of dispersive shock waves in optical fibers for the case of long-distance propagation of optical pulses, when the small Raman effect stabilizes the profile of the shock. The Whitham modulation equations are derived as the basis for the Gurevich-Pitaevskii approach to the analytical theory of such shocks. We show that the wave variables at both sides of the shock are related by the analogue of the Rankine-Hugoniot condition that follows from the conservation laws of the Whitham equations. Solutions of the Whitham equations yield the profiles of the wave variables that agree very well with the exact numerical solution of the generalized nonlinear Schroedinger equation for propagation of optical pulses.
- [2] arXiv:2507.23469 [pdf, html, other]
-
Title: Exact quasi-periodic solutions to the MKdV equationComments: 21 pages, 1 figureSubjects: Exactly Solvable and Integrable Systems (nlin.SI); Mathematical Physics (math-ph); Algebraic Geometry (math.AG)
In the present paper, a hierarchy of the mKdV equation is integrated by the methods of algebraic geometry. The mKdV hierarchy in question arises on coadjoint orbits in the loop algebra of $\mathfrak{sl}(2)$, and employs a family of hyperelliptic curves as spectral curves. A generic form of the finite-gap solution in any genus is obtained in terms of the $\wp$-functions, which generalize the Weierstrass $\wp$-function. Reality conditions for quasi-periodic wave solutions are completely specified. The obtained solutions are illustrated by plots in small genera.
- [3] arXiv:2507.23525 [pdf, html, other]
-
Title: Wave Turbulence and Cortical DynamicsSubjects: Adaptation and Self-Organizing Systems (nlin.AO)
Cortical activity recorded through EEG and MEG reflects complex dynamics that span multiple temporal and spatial scales. Spectral analyses of these signals consistently reveal power-law behaviour, a hallmark of turbulent systems. In this paper, we derive a kinetic equation for neural field activity based on wave turbulence theory, highlighting how quantities such as energy and pseudo-particle density flow through wave-space (k-space) via direct and inverse cascades. We explore how different forms of nonlinearity, particularly 3-wave and 4-wave interactions, shape spectral features, including harmonic generation, spectral dispersion, and transient dynamics. While the observed power-law decays in empirical data are broadly consistent with turbulent cascades, variations across studies, such as the presence of dual decay rates or harmonic structures, point to a diversity of underlying mechanisms. We argue that although no single model fully explains all spectral observations, key constraints emerge: namely, that cortical dynamics exhibit features consistent with turbulent wave systems involving both single and dual cascades and a mixture of 3- and 4-wave interactions. This turbulence-based framework offers a principled and unifying approach to interpreting large-scale brain activity, including state transitions and seizure dynamics.
New submissions (showing 3 of 3 entries)
- [4] arXiv:2507.23384 (cross-list from physics.bio-ph) [pdf, html, other]
-
Title: Could Living Cells Use Phase Transitions to Process Information?Comments: 8 pages, 5 figuresSubjects: Biological Physics (physics.bio-ph); Disordered Systems and Neural Networks (cond-mat.dis-nn); Soft Condensed Matter (cond-mat.soft); Adaptation and Self-Organizing Systems (nlin.AO); Cell Behavior (q-bio.CB)
To maintain homeostasis, living cells process information with networks of interacting molecules. Traditional models for cellular information processing have focused on networks of chemical reactions between molecules. Here, we describe how networks of physical interactions could contribute to the processing of information inside cells. In particular, we focus on the impact of biomolecular condensation, a structural phase transition found in cells. Biomolecular condensation has recently been implicated in diverse cellular processes. Some of these are essentially computational, including classification and control tasks. We place these findings in the broader context of physical computing, an emerging framework for describing how the native dynamics of nonlinear physical systems can be leveraged to perform complex computations. The synthesis of these ideas raises questions about expressivity (the range of problems that cellular phase transitions might be able to solve) and learning (how these systems could adapt and evolve to solve different problems). This emerging area of research presents diverse opportunities across molecular biophysics, soft matter, and physical computing.
- [5] arXiv:2507.23426 (cross-list from stat.ME) [pdf, html, other]
-
Title: Overcoming error-in-variable problem in data-driven model discovery by orthogonal distance regressionComments: 28 pages, 12 figures, prepared for the Data-driven systems and control: analysis, modelling, optimisation, and stochasticity collection in the journal Mathematics of Control, Signals, and SystemsSubjects: Methodology (stat.ME); Chaotic Dynamics (nlin.CD); Machine Learning (stat.ML)
Despite the recent proliferation of machine learning methods like SINDy that promise automatic discovery of governing equations from time-series data, there remain significant challenges to discovering models from noisy datasets. One reason is that the linear regression underlying these methods assumes that all noise resides in the training target (the regressand), which is the time derivative, whereas the measurement noise is in the states (the regressors). Recent methods like modified-SINDy and DySMHO address this error-in-variable problem by leveraging information from the model's temporal evolution, but they are also imposing the equation as a hard constraint, which effectively assumes no error in the regressand. Without relaxation, this hard constraint prevents assimilation of data longer than Lyapunov time. Instead, the fulfilment of the model equation should be treated as a soft constraint to account for the small yet critical error introduced by numerical truncation. The uncertainties in both the regressor and the regressand invite the use of orthogonal distance regression (ODR). By incorporating ODR with the Bayesian framework for model selection, we introduce a novel method for model discovery, termed ODR-BINDy, and assess its performance against current SINDy variants using the Lorenz63, Rossler, and Van Der Pol systems as case studies. Our findings indicate that ODR-BINDy consistently outperforms all existing methods in recovering the correct model from sparse and noisy datasets. For instance, our ODR-BINDy method reliably recovers the Lorenz63 equation from data with noise contamination levels of up to 30%.
- [6] arXiv:2507.23493 (cross-list from physics.flu-dyn) [pdf, html, other]
-
Title: Fluctuating interfaces in barotropic beta-plane turbulenceSubjects: Fluid Dynamics (physics.flu-dyn); Pattern Formation and Solitons (nlin.PS); Atmospheric and Oceanic Physics (physics.ao-ph)
Zonal jets manifest themselves as bands with sharp interfaces in the vorticity configuration. We develop an algorithm to track these fluctuating vorticity interfaces and systematically investigate their characteristic spatio-temporal behavior. While the interfacial height fluctuations are typically sub-Gaussian, the corresponding $\textit{fluctuation speeds}$ exhibit wider, heavy-tailed distributions reflecting the influence of lateral dispersion induced by the zonal velocity profile along the interfacial contours. The temporal evolution of these fluctuations is further characterized through their power spectrum displaying scale invariance in the frequency domain. The sharp, dense, shock-like features present in the time series of the $\textit{height}$ field suggest a possible lacking of differentiability. We confirm this by calculating the moments of the time-increments of the interfacial height fluctuations. Finally, the fractal nature of these boundaries is investigated systematically through a multifractal approach, revealing the non-trivial, complex statistics of interfaces in such geophysical, turbulent flows.
- [7] arXiv:2507.23498 (cross-list from math.DS) [pdf, html, other]
-
Title: On the lattice property of the Koopman operator spectrumSubjects: Dynamical Systems (math.DS); Functional Analysis (math.FA); Chaotic Dynamics (nlin.CD)
The Koopman operator has become a celebrated tool in modern dynamical systems theory for analyzing and interpreting both models and datasets. The linearity of the Koopman operator means that important characteristics about it, and in turn its associated nonlinear system, are captured by its eigenpairs and more generally its spectrum. Many studies point out that the spectrum of the Koopman operator has a multiplicative lattice structure by which eigenvalues and eigenfunctions can be multiplied to produce new eigenpairs. However, these observations fail to resolve whether the new eigenfunction remains in the domain of the Koopman operator. In this work, we prove that the spectrum of the Koopman operator associated to discrete-time dynamical systems has a multiplicative lattice structure. We further demonstrate that the Koopman operator associated to discrete-time stochastic process does not necessarily have such a structure, demonstrating an important nuance that lies at the heart of Koopman operator theory.
- [8] arXiv:2507.23596 (cross-list from physics.optics) [pdf, html, other]
-
Title: Nonlinear synchronization through vector subharmonic entrainmentSubjects: Optics (physics.optics); Adaptation and Self-Organizing Systems (nlin.AO)
Synchronization is ubiquitous across a wide range of fields. Subharmonic entrainment (SHE) is a nonlinear synchronization phenomenon that results in a locking oscillator at a frequency of an external periodic forcing signal with a fraction of the oscillator frequency. Beyond the fundamentals of nonlinear dynamics, SHE has a range of practical applications from stabilizing ultrafast laser pulses to optimizing control in various engineering and natural systems. However, the vectorial nature of SHE remains elusive. Here, we present the results of a theoretical and experimental study of a vector type of subharmonic entrainment (VSHE) using a passively mode-locked fiber laser as a testbed. We unveil the mechanism of vectorial SHE, in which weak external signals can entrain internal laser dynamics through vectorial coupling. Vectorial SHE presents in the form of synchronization between the subharmonic of mode-locking-driven oscillations and continuous wave (CW) signal through an evolving state of polarization. This CW signal, driven by the internal dynamics of the injected signal, causes VSHE with the frequencies ratios of multiples of ten, resulting in a partially mode-locking regime operation. Our findings offer new control techniques over mode-locking and additional dimension such as polarization states.
- [9] arXiv:2507.23621 (cross-list from cond-mat.stat-mech) [pdf, html, other]
-
Title: How to Incorporate Higher-order Interactions in Analog Ising MachinesComments: 16 pages, 7 figures, including Supplementary MaterialSubjects: Statistical Mechanics (cond-mat.stat-mech); Mathematical Physics (math-ph); Adaptation and Self-Organizing Systems (nlin.AO); Cellular Automata and Lattice Gases (nlin.CG); Applied Physics (physics.app-ph)
Ising machines (IMs) are specialized devices designed to efficiently solve combinatorial optimization problems. Among such problems, Boolean Satisfiability (SAT) is particularly relevant in industrial applications. To solve SAT problems using IMs, it is crucial to incorporate higher-order interactions. However, in analog IMs, interactions of different orders scale unevenly with the continuous spin amplitudes, introducing imbalances that can significantly degrade performance. We present a numerical comparison of methods to mitigate these imbalances, evaluating time-to-solution and success rate on Uniform Random 3-SAT instances from the SATLIB benchmark set. Our results show that the most effective approach employs spin interactions that are proportional to the signs of spins, rather than their continuous amplitudes. This generalizes our previous work, which showed that such interactions best mitigate imbalances induced by external fields in quadratic analog IMs. In this work, its advantage becomes substantially more pronounced, as it naturally mitigates imbalances across all interaction orders. We further demonstrate that smooth approximations of this method make it compatible with analog hardware. Our findings underscore the central role of spin-sign-based interactions in enabling robust and scalable analog IM dynamics.
Cross submissions (showing 6 of 6 entries)
- [10] arXiv:2410.15915 (replaced) [pdf, html, other]
-
Title: Asymptotic behaviors and dynamics of degenerate and mixed solitons for the coupled Hirota system with strong coherent coupling effectsSubjects: Pattern Formation and Solitons (nlin.PS)
In this work, we study the asymptotic behaviors and dynamics of degenerate and mixed solitons for the coupled Hirota system with strong coherent coupling effects in the isotropic nonlinear medium. Using the binary Darboux transformation, we derive the solutions to represent the degenerate solitons with two eigenvalues that are conjugate to each other. We obtain three types of degenerate solitons and provide their asymptotic expressions. Notably, these degenerate solitons exhibit time-dependent velocities, and the relative distance between the two asymptotic solitons increases logarithmically with the higher-order perturbation parameter $|\varepsilon|$ increasing. We also asymptotically reveal four interaction mechanisms between a degenerate soliton and a bell-shaped soliton: (1) elastic interaction with a position shift; (2) inelastic interaction for the degenerate soliton but elastic for the bell-shaped one; (3) elastic interaction for the degenerate soliton but inelastic for the bell-shaped one; and (4) the coherent interaction during a longer interaction region and elastic interaction based on specific parameter conditions. Besides, we analyze a special degenerate vector soliton that exhibits significant coherence effects, and numerically study the relationship between the robustness of such solitons and parameter $\varepsilon$. Our results indicate that $\varepsilon$ significantly affects the coherence of these solitons, and their robustness decreases when $|\varepsilon|$ increases.
- [11] arXiv:2202.12870 (replaced) [pdf, html, other]
-
Title: Correlation functions for open XXZ spin 1/2 quantum chains with unparallel boundary magnetic fieldsComments: 73 pages, several small correctionsJournal-ref: J. Phys. A: Math. Theor. 55 405203 (2022)Subjects: Mathematical Physics (math-ph); Other Condensed Matter (cond-mat.other); High Energy Physics - Theory (hep-th); Exactly Solvable and Integrable Systems (nlin.SI)
In this paper we continue our derivation of the correlation functions of open quantum spin 1/2 chains with unparallel magnetic fields on the edges; this time for the more involved case of the XXZ spin 1/2 chains. We develop our study in the framework of the quantum Separation of Variables (SoV), which gives us both the complete spectrum characterization and simple scalar product formulae for separate states, including transfer matrix eigenstates. Here, we leave the boundary magnetic field in the first site of the chain completely arbitrary, and we fix the boundary field in the last site $N$ of the chain to be a specific value along the $z$-direction. This is a natural first choice for the unparallel boundary magnetic fields. We prove that under these special boundary conditions, on the one side, we have a simple enough complete spectrum description in terms of homogeneous Baxter like $TQ$-equation. On the other side, we prove a simple enough description of the action of a basis of local operators on transfer matrix eigenstates as linear combinations of separate states. Thanks to these results, we achieve our main goal to derive correlation functions for a set of local operators both for the finite and half-infinite chains, with multiple integral formulae in this last case.
- [12] arXiv:2208.10097 (replaced) [pdf, html, other]
-
Title: On correlation functions for the open XXZ chain with non-longitudinal boundary fields : the case with a constraintComments: some small corrections, including corrected signs in the last formula, 29 pagesJournal-ref: SciPost Phys. 16, 099 (2024)Subjects: Mathematical Physics (math-ph); Strongly Correlated Electrons (cond-mat.str-el); High Energy Physics - Theory (hep-th); Exactly Solvable and Integrable Systems (nlin.SI)
This paper is a continuation of [1], in which a set of matrix elements of local operators was computed for the XXZ spin-1/2 open chain with a particular case of unparallel boundary fields. Here, we extend these results to the more general case in which both fields are non-longitudinal and related by one constraint, allowing for a partial description of the spectrum by usual Bethe equations. More precisely, the complete spectrum and eigenstates can be characterized within the Separation of Variables (SoV) framework. One uses here the fact that, under the constraint, a part of this SoV spectrum can be described via solutions of a usual, homogeneous, TQ-equation, with corresponding transfer matrix eigenstates coinciding with generalized Bethe states. We explain how to generically compute the action of a basis of local operators on such kind of states, and this under the most general boundary condition on the last site of the chain. As a result, we can compute the matrix elements of some of these basis elements in any eigenstate described by the homogenous TQ-equation. Assuming, following a conjecture of Nepomechie and Ravanini, that the ground state itself can be described in this framework, we obtain multiple integral representations for these matrix elements in the half-infinite chain limit, generalizing those previously obtained in the case of longitudinal boundary fields and in the case of the special boundary conditions considered in [1].
- [13] arXiv:2406.17552 (replaced) [pdf, html, other]
-
Title: Long-Time and Short-Time Dynamics in a Weighted-Median Opinion Model on NetworksComments: 30 pages, 13 figures, Submitted to SIAM Journal on Applied Dynamical Systems. Version 2 of this manuscript had a mistake in the arxiv title. The manuscript, figure, and appendix of this version is identical to version 2; the only change is the arxiv title which has been changed to align with the title of the manuscriptSubjects: Physics and Society (physics.soc-ph); Social and Information Networks (cs.SI); Dynamical Systems (math.DS); Adaptation and Self-Organizing Systems (nlin.AO)
Social interactions influence people's opinions. In some situations, these interactions eventually yield a consensus opinion; in others, they can lead to opinion fragmentation and the formation of different opinion groups in the form of ``echo chambers''. Consider a social network of individuals with continuous-valued scalar opinions, and suppose that they can change their opinions when they interact with each other. In many models of the opinion dynamics of individuals in a network, it is common for opinion updates to depend on the mean opinion of interacting individuals. As an alternative, which may be more realistic in some situations, we study an opinion model with an opinion-update rule that depends on the weighted median of the opinions of interacting individuals. Through numerical simulations of our median-update opinion model, we investigate how the final opinion distribution depends on network structure. For configuration-model networks, we also derive a mean-field approximation for the asymptotic dynamics of the opinion distribution when there are infinitely many individuals. We numerically investigate its accuracy for short-time opinion dynamics on various networks.
- [14] arXiv:2407.01450 (replaced) [pdf, html, other]
-
Title: Two-Parameter Quantum Groups and $R$-Matrices: Classical TypesJournal-ref: SIGMA 21 (2025), 064, 54 pagesSubjects: Representation Theory (math.RT); High Energy Physics - Theory (hep-th); Mathematical Physics (math-ph); Quantum Algebra (math.QA); Exactly Solvable and Integrable Systems (nlin.SI)
We construct finite $R$-matrices for the first fundamental representation $V$ of two-parameter quantum groups $U_{r,s}(\mathfrak{g})$ for classical $\mathfrak{g}$, both through the decomposition of $V\otimes V$ into irreducibles $U_{r,s}(\mathfrak{g})$-submodules as well as by evaluating the universal $R$-matrix. The latter is crucially based on the construction of dual PBW-type bases of $U^{\pm}_{r,s}(\mathfrak{g})$ consisting of the ordered products of quantum root vectors defined via $(r,s)$-bracketings and combinatorics of standard Lyndon words. We further derive explicit formulas for affine $R$-matrices, both through the Yang-Baxterization technique of [Internat. J. Modern Phys. A 6 (1991), 3735-3779] and as the unique intertwiner between the tensor product of $V(u)$ and $V(v)$, viewed as modules over two-parameter quantum affine algebras $U_{r,s}(\widehat{\mathfrak{g}})$ for classical $\mathfrak{g}$. The latter generalizes the formulas of [Comm. Math. Phys. 102 (1986), 537-547] for one-parametric quantum affine algebras.
- [15] arXiv:2412.09100 (replaced) [pdf, html, other]
-
Title: Generalized Liénard systems and isochronous connectionsComments: v1: Comments are welcome; v2: To appear in IJTP; v3: This version includes the acknowledgments section which was missing in the previous versionJournal-ref: Int. J. Theor. Phys. 64, 212 (2025)Subjects: Quantum Physics (quant-ph); Mathematical Physics (math-ph); Exactly Solvable and Integrable Systems (nlin.SI); Classical Physics (physics.class-ph)
In this paper, we explore some classical and quantum aspects of the nonlinear Liénard equation $\ddot{x} + k x \dot{x} + \omega^2 x + (k^2/9) x^3 = 0$, where $x=x(t)$ is a real variable and $k, \omega \in \mathbb{R}$. We demonstrate that such an equation could be derived from an equation of the Levinson-Smith kind which is of the form $\ddot{z} + J(z) \dot{z}^2 + F(z) \dot{z} + G(z) = 0$, where $z=z(t)$ is a real variable and $\{J(z), F(z), G(z)\}$ are suitable functions to be specified. It can further be mapped to the harmonic oscillator by making use of a nonlocal transformation, establishing its isochronicity. Computations employing the Jacobi last multiplier reveal that the system exhibits a bi-Hamiltonian character, i.e., there are two distinct types of Hamiltonians describing the system. For each of these, we perform a canonical quantization in the momentum representation and explore the possibility of bound states. While one of the Hamiltonians is seen to exhibit an equispaced spectrum with an infinite tower of states, the other one exhibits branching but can be solved exactly in closed form for certain choices of the parameters.
- [16] arXiv:2501.08424 (replaced) [pdf, html, other]
-
Title: Isochronous oscillator with a singular position-dependent mass and its quantizationComments: v1: Comments are welcome; v2: To appear in JMPSubjects: Quantum Physics (quant-ph); Mathematical Physics (math-ph); Exactly Solvable and Integrable Systems (nlin.SI); Classical Physics (physics.class-ph)
In this paper, we present an analysis of the equation $\ddot{x} - (1/2x) \dot{x}^2 + 2 \omega^2 x - 1/8x = 0$, where $\omega > 0$ and $x = x(t)$ is a real-valued variable. We first discuss the appearance of this equation from a position-dependent-mass scenario in which the mass profile goes inversely with $x$, admitting a singularity at $x = 0$. The associated potential is also singular at $x = 0$, splitting the real axis into two halves, i.e., $x > 0$ and $x < 0$. The dynamics is exactly solvable for both the branches and so for definiteness, we stick to the $x > 0$ branch. Performing a canonical quantization in the position representation and upon employing the ordering strategy of the kinetic-energy operator due to von Roos, we show that the problem is isospectral to the isotonic oscillator. Thus, the quantum spectrum consists of an infinite number of equispaced levels. The spacing between the energy levels is found to be insensitive to the specific choices of the ambiguity parameters that are employed for ordering the kinetic-energy operator à la von Roos.
- [17] arXiv:2501.16325 (replaced) [pdf, html, other]
-
Title: Tailored Forecasting from Short Time Series via Meta-learningComments: 23 pages, 12 figuresSubjects: Machine Learning (cs.LG); Chaotic Dynamics (nlin.CD); Computational Physics (physics.comp-ph)
Machine learning models can effectively forecast dynamical systems from time-series data, but they typically require large amounts of past data, making forecasting particularly challenging for systems with limited history. To overcome this, we introduce Meta-learning for Tailored Forecasting using Related Time Series (METAFORS), which generalizes knowledge across systems to enable forecasting in data-limited scenarios. By learning from a library of models trained on longer time series from potentially related systems, METAFORS builds and initializes a model tailored to short time-series data from the system of interest. Using a reservoir computing implementation and testing on simulated chaotic systems, we demonstrate that METAFORS can reliably predict both short-term dynamics and long-term statistics without requiring contextual labels. We see this even when test and related systems exhibit substantially different behaviors, highlighting METAFORS' strengths in data-limited scenarios.
- [18] arXiv:2502.03093 (replaced) [pdf, html, other]
-
Title: Stabilizer Entropy and entanglement complexity in the Sachdev-Ye-Kitaev modelSubjects: Quantum Physics (quant-ph); Strongly Correlated Electrons (cond-mat.str-el); High Energy Physics - Theory (hep-th); Chaotic Dynamics (nlin.CD); Exactly Solvable and Integrable Systems (nlin.SI)
The Sachdev-Ye-Kitaev (SYK) model is of paramount importance for the understanding of both strange metals and a microscopic theory of two-dimensional gravity. We study the interplay between Stabilizer Rényi Entropy (SRE) and entanglement entropy in both the ground state and highly excited states of the SYK4+SYK2 model interpolating the highly chaotic four-body interactions model with the integrable two-body interactions one. The interplay between these quantities is assessed also through universal statistics of the entanglement spectrum and its anti-flatness. We find that SYK4 is indeed characterized by a complex pattern of both entanglement and non-stabilizer resources while SYK2 is non-universal and not complex. We discuss the fragility and robustness of these features depending on the interpolation parameter.