2026-09-07 | | Total: 11
We present a simple algebraic proof of meromorphic non-integrability of a 'mysterious' case of a quadratic Hamiltonian studied by Sokolov and Wolf. The result is obtained by applying the Morales-Ramis method, based on differential Galois theory.
Solutions of the Whitham modulation equations for one-phase periodic waves obeying the sine-Gordon equation are found that describe the evolution of an oscillatory region behind an instability front propagating into the instability region. A simple self-similar solution describes the whole region between two fronts of instability resulting from a localized initial disturbance in the unstable state. Another hodograph solution represents typical waves close to the instability fronts. This theory generalizes the approach used previously for systems obeying the nonlinear Schroedinger equation.
There is a body of research devoted to understanding how complex singularities of solutions of nonlinear partial differential equations (pdes) spontaneously emerge at $t=0^+$ and propagate for $t>0$, and how their behaviour affects the solution on the real axis. Despite the importance of the small-time limit in these studies, there is still a lack of understanding of how complex singularities are born at $t=0^+$, including for initial conditions that are not analytic functions of the spatial variable. In this paper, we use Burgers' equation as a prototype nonlinear pde and study the complex-plane singularities for initial conditions that are piecewise smooth. Using matched asymptotic expansions, we show how infinitely many singularities emerge from points of discontinuity in a pattern that can be described using branches of the Lambert-$W$ function. For various initial conditions, we observe how these singularities rearrange themselves to align with the appropriate exactly-described long-time behaviour, including sigmoid-shaped travelling waves, constant-area (triangular wave) similarity solutions and $N$-wave solutions. In terms of Burgers' equation, our small-time asymptotic analysis of the singularity propagation for piecewise-continuous initial conditions illustrates the types of generic behaviours that arise for inner regions when diffusion dominates advection. More generally, this work is a step towards understanding complex-plane behaviour of solutions of nonlinear partial differential equations with non-analytic initial conditions.
We present a Hessian-based approach to numerically continue periodic orbits in dynamical systems. A loop (periodic orbit candidate) is parametrized as a Fourier series; a loss function is defined based on the deviation of the loop from the physical differential equations. Unlike previous work relying on hand-derived Jacobians, our method automates the process by leveraging automatic differentiation, a common machine learning technique. The continuation direction can be determined by the flat directions of the loss landscapes (directions with zero eigenvalues), making the search of periodic orbits efficient and guided. Our method is integrator-free, precisely initializes oscillations around unstable fixed points, and efficiently detects orbit family intersections and subharmonic bifurcations. As a demonstration, we present full continuations of periodic double pendulum oscillations from fixed points, showing bifurcations along orbit families and categorizing branches of periodic orbits. In particular, we find periodic orbits where both pendulum masses are never simultaneously at rest, which to our knowledge has been missing in the literature.
Chaotic evolution measures how rapidly nearby configurations separate, but the existence of such rapid evolution does not by itself determine what physical information the evolution actually carries. For fixed-coupling JIMWLK evolution with fresh transverse-local, color-isotropic noise, we show that this distinction becomes exact - leading to a no go theorem. At finite rapidity step and finite transverse regulator, the exact linearized update depends on the Wilson-line background only through local adjoint rotations of fresh color-isotropic noise. Once the past is conditioned upon, these rotations leave the complete local noise measure invariant, so the tangent process has the same probability law for any target ensemble evolved with the same regulated fixed-coupling dynamics, including dilute and saturated ensembles. After the exact global color-rotation mode is removed, the reduced tangent process found to have a positive leading Lyapunov exponent which indicate chaotic evolution. Fixed-coupling JIMWLK therefor can be dynamically chaotic while its tangent-only statistics carry no information about the target ensemble. Spatially nonlocal noise however evades the argument because correlations between distinct sites leave relative Wilson lines in the conditional covariance.
We study the operator entanglement of local operators in one-dimensional brickwork circuits whose two-site gate satisfies the braid relation; throughout this work, we call such a gate a Yang--Baxter gate. We establish upper bounds for several structured, overlapping classes of Yang--Baxter gates. We show that the operator Schmidt rank remains uniformly bounded in time for all qubit Yang--Baxter gates and, in arbitrary local dimension, for permutation gates obtained from non-degenerate Yang--Baxter maps. We also show that it grows at most polynomially for involutive dual-unitary Yang--Baxter gates and for arbitrary phase dressings of permutation gates obtained from non-degenerate Yang--Baxter maps. These results imply, respectively, constant and logarithmic upper bounds on the operator entanglement. Conversely, we construct a seven-state involutive Yang--Baxter gate without dual unitarity and a one-site operator whose exact operator Schmidt rank grows exponentially, although the corresponding operator entropies remain undetermined. Entanglement growth in the general Yang--Baxter case remains open. All proofs and selected examples were constructed by ChatGPT 5.6 Sol.
Detecting weak transient signals in noise is a persistent challenge in sensing, communication, and electrophysiology. Here, we demonstrate a physical weak-signal detector based on a semiconductor negative differential resistance (NDR) device operated near a Hopf bifurcation. A coherent input that persists over the response time of the dynamical system can drive a transition between quiescent and oscillatory states, whereas faster stochastic fluctuations are largely suppressed. This nonlinear transformation converts a weak analog threshold crossing into all-or-none voltage spikes and therefore provides asynchronous signal detection without a reference clock. Using a modulated photovoltaic signal, we detect a weak frequency component of 100 Hz as a demonstration, at an input signal-to-noise amplitude ratio as low as 1/500, and reliably recover it in the output spectrum. We then apply the same principle to neuronal multisite extracellular recordings. After standard band-pass filtering, the raw microelectrode signal is transformed by the NDR dynamics, enhancing the distinction between neuronal spikes and background fluctuations. The resulting detected spike times agree closely with those obtained using a traditional spike-detection pipeline. These results establish bifurcation-engineered NDR dynamics as a compact hardware approach to noise-resilient signal discrimination and event-based analog-to-digital conversion that will be useful for neuroprosthetic devices.
Higher-order networks exhibit rich critical phenomena that cannot be captured by traditional pairwise models. Here, we develop an analytical dimension-reduction framework that maps higher-order networked dynamics onto an effective low-dimensional system, allowing accurate prediction of tipping boundaries, bistability regions, and the nature of phase transitions. We demonstrate the power of this framework across a range of dynamical processes, revealing distinct effects of higher-order interactions on transition continuity and hysteresis. Furthermore, we find that system resilience exhibits a profound dependence on the alignment between pairwise and higher-order connectivity, with assortative mixing enhancing tipping toward active states. Our findings establish a general theory for understanding the critical transitions in higher-order networks, offering new insights for anticipating and managing systemic risk in complex systems.
Unstable periodic orbits and steady states are believed to form the backbone of spatiotemporal chaos and turbulence, yet their computation in thermally driven flows remains scarce for transitional regimes. In this work we compute and characterize a steady state and three families of periodic orbits in two-dimensional Rayleigh-Bénard at $\mathrm{Pr}=1$, near the transition to chaos. We find that in its route to chaos, the flow hops between several sets of orbits after becoming quasiperiodic. We use Floquet analysis to study the stability of the orbits obtained and characterize their bifurcations, showing how the appearance of primary and secondary frequencies, as well as phase-locking mechanisms, are all related to the dynamics of the orbits. We study in detail how the flow shadows the orbits found and determine in which regimes each orbit is dynamically relevant or not. Our analysis also reveals two important insights: (1) all symmetries are broken before the onset of chaos, and (2) this onset does not alter the behavior of the heat transport.
We study the stability of collective neutrino oscillations using a distributional approach motivated by the statistical mechanics of Kuramoto synchronization. Treating the ensemble of neutrino flavor polarization vectors in the thermodynamic limit $N\to\infty$, we derive an exact nonlinear Fokker--Planck (continuity) equation for the one-body distribution $F(\vec{\mathbf{S}},ω,t)$ on the flavor sphere. This equation admits a two-parameter family of azimuthally symmetric stationary solutions, whose stability we analyze by linearizing around them. The resulting eigenvalue condition determines the growth or decay rate of small perturbations from \emph{any} initial distribution -- not merely from a state close to full flavor coherence -- thereby going significantly beyond the conventional linear stability analysis of collective modes. In special limits the condition reproduces known synchronization thresholds in the two-beam model, providing a non-trivial check of the framework. We present analytical results for the eigenvalue equation and explore stability phase diagrams for physically relevant frequency distributions.
This paper describes classical soliton scattering in the 2+1d integrable chiral model, related to the Bogomolny monopole equation and self-dual Yang-Mills, as a case study for higher-dimensional integrable scattering. Extended line solitons exhibit nontrivial scattering with lump solitons and with other line solitons, and $N \to N$ interactions can be factorised into a series of $2 \to 2$ interactions. Scattering between parallel line solitons in 3d is related to lump soliton scattering in 2d, allowing for direct comparison with results in 1+1d integrable field theories.