Papers

Topological signals are dynamical variables supported on higher-order structures such as simplicial or cell complexes. In this work, we investigate the impact of time-delayed interactions on the emergence of global synchronization of topological oscillators, with application to the Stuart-Landau system. We first examine the case where signals are supported on simplexes, or cells, of a given dimension and coupled through the Hodge-Laplace matrix. We derive the Master Stability Function and we prove that for suitable discrete delay values, the stability of the synchronous solution becomes analytically tractable and can be solved by using the Lambert W-function. This analysis yields explicit spectral conditions for Global Topological Synchronization to emerge and shows that admissible delays may either promote or suppress synchronization. We then study delayed interactions between signals supported on simplexes, or cells, of different dimensions coupled by the Dirac operator. Under analogous admissible time delays the variational problem, still related to the Master Stability Function, can again be simplified, however in this case we prove that Global Topological Dirac Synchronization cannot emerge. Numerical simulations on simplicial and cell complexes corroborate and complement these findings.

We examine the spectral dissipation-coherence trade-off inequality conjectured by Oberreiter et al. [Phys. Rev. E 106, 014106 (2022)], stating that the coherent number, defined from the second largest eigenvalue, provides a lower bound on the stationary entropy production per unit oscillation. We disprove this conjecture by explicitly constructing a counterexample, which accompanies finite coherent number with vanishing stationary entropy production. This model also excludes a wide class of thermodynamic bound with the real and imaginary parts of the second largest eigenvalue of the transition rate matrix. Our result implies that the second largest eigenvalue does not necessarily characterize the degree of coherent oscillation.

We study the monic recurrence coefficients of the generalized Charlier polynomials associated with the discrete weight \[ \rho_\mu(k) = \frac{\mu^k}{(k!)^2}, \qquad k\in\mathbb N_0, \qquad \mu>0. \] The coefficients $\beta_n$ and $\gamma_n$ in the recurrence relation \[ P_{n+1}(x) = (x-\beta_n)P_n(x)-\gamma_nP_{n-1}(x) \] are described through complementary moment, determinant, discrete dynamical, and continuous integrable formulations. We give a consistent normalization of the Markov function and its Jacobi continued fraction, and express the recurrence coefficients in terms of Hankel determinants, using a bordered determinant for the diagonal coefficient $\beta_n$. Starting from the original Laguerre–Freud equations, we derive the local recursive system used for computation and relate it to the Toda deformation. We also recall the known Painlevé~V representations associated with the two-parameter generalized Charlier family, in a normalization consistent with the present specialization. Moment/Hankelconstructions provide complementary moment-theoretic realizations of the recurrence coefficients, while the local Laguerre–Freud recursion provides an independent dynamical construction. Their comparison is used for numerical cross-validation. The numerical results are consistent with the leading large-degree behavior \[ \beta_n\sim n, \qquad \gamma_n\sim\mu. \]

We demonstrate the generation of near-visible, broadband Kerr frequency combs using a tapered fibre Fabry-Pérot resonator. By incorporating micron-diameter tapered silica fibre into a pulsed-driven Fabry-Pérot configuration, we exploit geometric dispersion engineering to overcome the strong normal dispersion intrinsic to silica (and nearly all materials) at short wavelengths, achieving net anomalous dispersion in the near-visible regime. The strong mode confinement of the tapered fibre simultaneously enhances the Kerr nonlinearity by over an order of magnitude relative to standard fibre, substantially reducing the threshold drive power for nonlinear structure formation. Through precise control of the cavity fabrication, we experimentally demonstrate and numerically confirm a rich suite of coherent, phase-locked nonlinear structures, including conventional anomalous dispersion cavity solitons, zero dispersion solitons with spectral bandwidths exceeding 16 THz between 750 nm and 800 nm, as well as switching wave combs and broadband Raman solitons, the latter also reported here for the first time in the near-visible regime. Our results establish the tapered fibre Fabry-Pérot resonator as a versatile platform for near-visible and potentially visible coherent frequency comb generation.

In many turbulent flows, the quantities of interest are coarse-grained variables, either because resolving the fine scales is computationally prohibitive or because only coarse observations are available; in both cases, data-driven reduced-order models that evolve only these variables are required. Even when the governing equations are known, their projection onto the coarse variables is not closed, and the central difficulty lies in representing the feedback of the unresolved scales so that the reduced model reproduces the statistical and dynamical properties of the projected high-resolution system. We address this problem with a score-based framework, constructing a reduced-order model of quasi-geostrophic turbulence whose stationary distribution coincides by construction with the distribution learned from data and whose dynamics are calibrated to reproduce the finite-time correlations of the energies in individual wavenumber shells. The resulting model reproduces the stationary statistics of the resolved flow with high accuracy, and its calibration is computationally inexpensive because it requires no forward integration of the model.

The one-dimensional swarmalator model is usually understood through three macroscopic states: synchrony, the phase wave, and asynchrony. We show that the identical model also contains a much larger set of hidden stable equilibria - families of q-twisted and m-clustered states - which random initial conditions almost never reach, because many are embedded in a degenerate neutral manifold of vanishing global basin. They are nonetheless locally robust: under Gaussian perturbations of per-coordinate amplitude {\epsilon} their survival radii shrink only algebraically, as 1/q and 1/m, in sharp contrast to the Gaussian basins and winding ceiling of twisted states on a nearest-neighbour Kuramoto ring, both removed by mean-field coupling. A block-circulant Jacobian fixes the stability in closed form, and far-kick experiments indicate the basins are compact cores rather than tentacled sets. Several patterns reported in extended swarmalator models thus already exist, stable or neutrally stable, in the minimal one.

Reduced descriptions of unresolved fluctuations are commonly based on Gaussian processes, although many realistic forcings have finite correlation times and a renewal structure. For memoryless step (Kubo–Anderson) noise, we show that the reduced dynamics admit two exact and complementary descriptions: a kernel representation governed by the Boolean cumulants of the jump distribution, and a frozen-noise representation adapted to stationary probability densities. For exponentially distributed waiting times, the totally time-ordered $G$-cumulants coincide with the Boolean cumulants of the jump law, and the memory kernel is resummed exactly as the Boolean generator $\eta$ evaluated on a resolvent operator. Second-order closure is exact if and only if the jump law is symmetric Bernoulli; otherwise, the leading closure error is controlled by $(b_4/b_2)\lambda^2$. The Boolean hierarchy, however, acts on the kernel, not on the stationary measure. Stationary densities are approximated by replacing the jump law with its $N$-point Gauss quadrature: each surrogate preserves all multi-time correlations up to order $2N-1$ for any waiting-time law, its kernel is a Padé resummation of the Boolean one, and it is itself an exactly solvable renewal problem. For the linear system with linear multiplicative interaction (LIMI/CAM), the frozen-noise representation reduces the dynamics to a random affine recursion and yields exact support boundaries, singularity exponents and Kesten tail indices, confirmed numerically. Rates, moments and closure errors are governed by the Boolean hierarchy, whereas stationary densities are determined by the geometry of the jump distribution. For memoryless renewal noise, the combinatorial structure controlling finite-correlation reductions is Boolean.

We develop an analytical framework for studying a family of interacting particle systems that we refer to as the multi-particle Oja flow. The common interaction structure underlying this family arises in a variety of applications, including higher-order synchronization models, opinion dynamics, and self-attention dynamics derived from Transformer architectures, motivating a unified mathematical framework for their analysis. We first establish well-posedness for both the finite-particle and the mean-field systems. We also prove convergence of the finite-particle dynamics to the mean-field equation as the number of particles tends to infinity. We then investigate low-dimensional structures underlying the dynamics. In particular, we show that the mean-field solution admits a representation through a time-dependent normalized linear transformation acting on the initial distribution, with the transformation governed by a matrix ODE. Building on this representation, we identify a broad class of parametric families that are invariant under the mean-field dynamics. When the initial distribution belongs to one of these families, the evolution remains within the same family for all time and can therefore be characterized through the dynamics of the corresponding family parameters. As a concrete example, we develop this reduction within the angular central Gaussian family and illustrate how the resulting parameter dynamics provide reduced descriptions of the mean-field evolution for several representative systems.

Criticality in the Ising model is conventionally generated by Hamiltonian dynamics and controlled by temperature. Here we show that critical like behavior can emerge in an Ising system whose dynamics is completely independent of the Ising Hamiltonian and contains no temperature like parameter. We introduce an adaptive cluster dynamics in which spin connectivity is controlled by a local quiet time variable, the time elapsed since a spin was last updated. The Ising Hamiltonian enters only as an observable characterizing the resulting configurations. As the parameter $\alpha$, controlling the connectivity strength, is varied, the system undergoes spontaneous symmetry breaking accompanied by strong collective fluctuations. The cluster size distribution develops an extended power law regime, terminated by a small finite size hump at the largest scales, while the Ising energy exhibits a singular response to the variation of $\alpha$. These results show that memory dependent connectivity alone can generate collective critical behavior, revealing a route to nonequilibrium criticality without energy based dynamics.

We report on the recent results in the theory of discrete integrable equations reduced modulo a prime number $p$ based on the previous results. This article mainly concerns the $q$-discrete analogues of the Painlevé equations. We first give a quick review of the idea of `almost good reduction', which is a property that most integrable systems over a finite field $F_p$ possess when reduced from those defined over the field of $p$-adic numbers $Q_p$. This idea has been introduced in our previous papers and the discrete and q-discrete Painlevé II equations have been proved to have `almost good reduction', which can be interpreted as an arithmetic analogue of the singularity confinement test. In this article we prove that this property holds for several other $q$-discrete Painlevé equations. The theory of Diophantine integrability is reviewed and its relation to the integrability detection by the algebraic entropy is discussed.

Traveling waves of activity in neural tissue can be halted by spatial inhomogeneity in synaptic coupling. We study an integrate-and-fire network in which each neuron fires once and the coupling decays exponentially. For this model the leading-edge firing map reduces exactly to a scalar equation for the wave speed in space, with a slow unstable and a fast stable homogeneous speed, $c_1$ and $c_2$, and a periodic modulation of the coupling enters this equation pointwise. Positive periodic waves terminate in folds. For slowly varying modulation the fold amplitude approaches a plateau set by the trough of the modulation, where the local bottleneck speed is $\sqrt{c_1 c_2}$; the approach to the plateau is set by the local geometry of the trough. For rapidly varying modulation the fold amplitude grows linearly with frequency, with a slope given by a full-amplitude average over the integrated modulation profile; a weak-ripple truncation of that average overestimates the slope about 2.5-fold at the default coupling. A perturbative recursion in the amplitude gives the speed profile explicitly. We distinguish loss of the periodic wave from failure of a particular launch, and the excitatory regime $\epsilon\le 1$ from its sign-changing extension. The results are checked against direct integration of the reduced equation and first-spike firing-map simulations of the network.

The long-run statistics of chaotic dissipative partial differential equations are described by invariant measures, and when the equation carries a parameter these measures form a family whose members can differ in kind, from steady states to sustained spatiotemporal chaos. We ask whether a single generative model can cover such a family. We train a single conditional diffusion model once, on states of the Kuramoto-Sivashinsky equation at 160 values of its hyperviscosity, and evaluate it against held-out data at interleaved parameter values not seen in training, using a classifier two-sample test supported by spectral, geometric and tail statistics. Over most of the family, agreement improves as the dynamics become more chaotic. We identify a structural cause: the weakly chaotic members have invariant measures concentrated on low-dimensional sets, which a model with full support output cannot represent, and the target approaches full dimensionality as chaos increases. Restricting training to the strongly chaotic members reduces the classifier's excess over chance on those members by roughly 40%, and this restricted conditional model matches or exceeds models trained at single parameter values. On those members, generated samples reproduce energy spectra to within a resampling floor, with a small systematic deficit at the largest scales, and match the energy balance of held-out data.

The sudden removal of a potential barrier from a Bose-Einstein condensate (BEC) gives rise to a quantum version of a hydrodynamic dam break and leads to rich wave dynamics that have similarities to other localized defect problems such as domain wall dynamics in spin systems. Denoting $\Delta n$ as the initial difference in density between the upper $n_{1}$ and lower $n_{0}$ reservoirs on either side of the dam, we use the Gross-Pitaevskii equation to study the quasi-one dimensional case in both perturbative ($\Delta n \ll n_{0}$) and non-perturbative ($\Delta n \sim n_{1}$) regimes. In the perturbative regime a pair of outwardly propagating dispersive wavepackets forms which can be viewed as quantum versions of undular tidal bores that have analytic forms at long times in terms of the integrals of Airy functions with a wavelength that grows as $(\hbar^2 t)^{1/3}$. Airy functions are the universal wave functions that dress structurally stable fold caustics where pairs of rays coalesce, and, indeed, we show that the quantum dam break problem has the same ray structure as the naturally occurring caustic phenomenon of a double rainbow, including Alexander's dark band between the two bows where no light is scattered: we identify the intermediate density plateau in the dam break as an analogous `silent band' where only evanescent sound waves can exist. In the opposite regime of a non-perturbative dam break we show that a self-induced sonic horizon occurs when $\Delta n \geq (8/9)n_{1}$. A discussion of possible experimental schemes for amplification of quantum undulations is included as well as an alternative scheme for dispersive wave generation where a constant flow is imprinted on a BEC in box trap.

Skyrmions experience a sideways gyroscopic force when they move. Bound molecules can cancel this force; we call them gyroscopically neutral. We show that the quantum states available on each site of a magnet select the composition needed for cancellation. In two models of $SU(3)$ magnetism, this composition changes from a pair to a trimer. Both models assign exactly the same energy to every configuration of their common classical order parameter, isolating the effect of the local quantum states. Neutral molecules can also unwind when additional local states are allowed. Mean field calculations nevertheless find both compositions bound and locally stable, with the neutral trimer surviving sampled interaction changes. Comparing their motion reveals the changed cancellation rule, and prepared neutral molecules translate over short distances with little deformation.

Epidemic spread and collective opinion formation are tightly coupled processes: individual beliefs about the severity of an outbreak and the value of protective measures shape behavior, while the evolving epidemic itself reshapes public sentiment through risk perception. We propose an agent-based framework that couples opinion dynamics with a compartmental epidemic model, in which each agent carries both a health status and a continuous opinion variable describing adherence to protective measures. Opinion evolution is driven by binary interactions among agents and by an exogenous background field that reflects the current state of the epidemic, while the disease transmission rate depends explicitly on the opinions of interacting agents. Using methods of statistical mechanics, we characterize the resulting macroscopic system for the compartmental mass fractions and mean opinions in both far-from-equilibrium and close-to-equilibrium regimes. The calibration shows that incorporating behavioral opinion dynamics substantially improves the interpretation of epidemic evolution. In particular, the inferred risk perception function evolves consistently with the observed epidemic burden, indicating that behavioral adaptation can explain part of the temporal variation that would otherwise be absorbed into time-dependent transmission parameters. Our results illustrate how integrating opinion heterogeneity in kinetic epidemic models offers quantitative insights into the feedback between behavior and disease spread.

We introduce a new method for constructing Zoll metrics that admit first integrals polynomial in the momenta. The method uses Funk's form of the metric together with an algebraic relation between the first integrals. This reduces the problem to an algebraic system from which all unknown functions can be determined explicitly. In the case of a non-trivial cubic integral, we present an explicit two-parameter family of metrics that covers all such rotationally symmetric Zoll metrics on the $2$-sphere: every such metric is isometric to a member of this family. We compare these metrics with the earlier results of Valent, Duval, and Shevchishin.

We study the Cauchy problem for the Tzitzéica equation, which is an important integrable model arising in affine differential geometry and characterizing proper affine spheres. Recently, Huang, Wang and Zhu (Math. Ann. 395, 18 (2026)) reported the long-time asymptotic result for the Tzitzéica equation for the case of the purely continuous spectrum. Inspired by their work, we conduct in-depth research on the soliton resolution conjecture and asymptotic analysis of the Tzitzéica equation. Compared with the previous work, by developing the $\bar{\partial}$-nonlinear steepest descent method, we reveal that the solution of the Tzitzéica equation exhibits three different asymptotic regions depending on $\xi:=x/t$. For the region $\xi\in(-\infty,-1]\cup[1,\infty)$, where there are no stationary phase points, we rigorously prove that the solution of the Tzitzéica equation decays algebraically to zero with error $O(t^{-1})$. In the region $\xi\in(-1,1)$, the phase function $\theta(z)$ has two stationary phase points. The corresponding asymptotic approximations can be characterized with an $N$-soliton solution as well as an interaction term between soliton solutions and the dispersion term with a residual error of order $O(t^{-3/4})$. In the transition regions $\{(x,t):1-\varepsilon<|\xi|<1\}$, the corresponding asymptotic approximation can be characterized by the $N$-soliton solution together with the interaction terms between the soliton solution and the dispersion term, with a residual error of order $O(t^{-3/4})$. Our results provide a rigorous analysis of the long-time asymptotic behavior of solutions of the Tzitzéica equation in different regions and confirm the soliton resolution conjecture for the Tzitzéica equation.

We study nonlinear holes, bubbles, kink fronts, and quantum-droplet-like states in an elongated Bose mixture governed by an extended Gross-Pitaevskii equation with attractive cubic mean-field interactions and repulsive quartic Lee-Huang-Yang contribution. We classify the constant-amplitude backgrounds: the Bogoliubov-de Gennes spectrum shows that the upper branch is modulationally stable, whereas the lower branch is unstable to long-wavelength perturbations. Grand-potential and pressure determine the vacuum-finite-density coexistence point and the unique background supporting stationary kink and antikink fronts. Symmetric holes on the stable background have a density deficit and integral width that diverge near coexistence as the hole separates into two nonlinear fronts. A multiplicative kink-antikink construction connects bell-shaped states to broad flat-top droplets and approaches the stationary quantum-droplet family as the front separation increases. An antikink-kink sum with equal signs generates bubble-like density depletions without a phase jump, whereas their difference produces dark-hole-like profiles with a $\pi$ phase difference between the asymptotic backgrounds. Bubble depth increases with front separation, distinguishing gray-like and dark-like profiles. Real-time simulations show that deep bubbles persist under weak random perturbations, demonstrating finite-time robustness without establishing spectral stability. In the asymmetric two-component model, localized relative-density modes can induce spinor-driven breakup on a modulationally stable background, while robust holes persist near the kink limit. Collisions of kink-antikink droplet-like states display phase-dependent merging, asymmetric particle transfer, and effective reflection. These results connect nonlinear fronts with bubbles, phase-jump holes, and self-bound droplets in beyond-mean-field Bose mixtures.

Even if it were possible, attempting to model the motion of a fluid in a pipe by tracking each individual molecule would be an extraordinary waste of computation. The emergent laws of fluid mechanics accomplish the same goal on a laptop instead of a datacenter. But a sufficiently complex fluid such as in Earth's atmosphere again requires a datacenter, this time to resolve the many contortions in the fields making up our weather. Might there be some analogous set of emergent laws that would represent the effect of these perturbations without needing their explicit representation on a numerical grid? We argue that the existence of a symmetry principle implies such laws are likely to exist. The atmosphere is scale invariant, implying that the observed structure at one scale is nothing but a stretched version of that at any other scale, at least within the invariant regime. Though some emergent laws built on scale invariance exist, they are built on a challenging ontological foundation, making them hard to develop further. We propose a path forward, first by introducing a new fundamental building block of the atmosphere that we term a “turbulon”, inspired by a lesser-known generalization of turbulence theory that accounts for buoyancy. We then construct the Superposition of Turbulons and Eddies Atmospheric Model (STEAM), and show that simulated atmospheric volumes are broadly plausible when compared to state-of-the-art hydrodynamic simulation output. Some aspects of STEAM clearly need improvement, though other statistics reproduce observations better than hydrodynamic models, and visualizations of simulated clouds are strikingly realistic. We estimate STEAM's computational cost to be up to a million times less than hydrodynamic models, suggesting that the Navier-Stokes equations may not, in fact, be required for simulation and prediction of Earth's weather and climate.

We study pulse-number selection in a stochastic cubic–quintic complex Ginzburg–Landau model of a weakly seeded, normal-dispersion laser in the dissipative-soliton-resonance (DSR) regime. Without noise, a single pulse and pulse pairs persist at the same control parameters, and the total energy of an $N$-pulse state follows a ladder constructed from the single-pulse branch. The final energies of noisy trajectories lie close to the same ladder. Multipulsing therefore does not necessarily lose the single-pulse solution. It can reflect which coexisting state the noisy dynamics reaches. Optical seed injection suppresses energy-dependent multipulsing and, at larger seed power, produces a nonmonotonic DSR energy window. An energy–noise scan shows that energy dominates the multipulse probability, whereas additive noise produces only a modest trend common to all energies, without a noise optimum. The noise-induced formation statistics therefore do not establish canonical stochastic resonance or escape from a pre-existing soliton. Coherent control by a weak monochromatic seed extends predominantly single-pulse operation over a noticeably broader energy window, enhancing dissipative-soliton energy scalability. Within an adiabatic approximation, equal energy sharing among coexisting pulses is stable wherever the single-pulse energy grows less than proportionally with the control energy, as it does over the sampled DSR range. Energy exchange between the pulses then relaxes up to about 200 times more slowly than their total energy.

Influenza A persists by continually changing the antigens that its host population recognizes. How fast it evolves and how many hosts it infects are not independent, because the immunity that selects new variants is produced by the infections the pathogen has already caused. We analyze a model of Andreasen et.\ al.\ that captures this eco-evolutionary feedback, and we extend it to the case where infection may kill the host. In the model infected and recovered hosts are distributions over a one-dimensional antigenic space coupled by cross-immunity. The endemic state is a traveling pulse: the infected distribution keeps a fixed shape and advances through antigenic space at constant speed. In the limit of weak antigenic mutation we construct the pulse asymptotically and find that: (i) the speed of antigenic evolution decreases as prevalence rises, (ii) less standing antigenic variation accompanies either faster or slower antigenic evolution, depending on whether pathogens differ in immune evasion or host lifespan (iii) the total number of infected hosts is largest at an intermediate virulence when the host birth rate is sufficiently low.

The Yukawa-Sachdev-Ye-Kitaev (YSYK) model, a solvable model of non-Fermi liquids with a proposed cavity-QED realization, couples $N$ complex fermions to $M$ bosons of frequency $\omega_0$ through random Yukawa couplings of strength $g$, and its dynamics depends only on $R\equiv\omega_0/g^{2/3}$. For $R\ll1$ the bosons are slow and act as a static random background in which the fermions are free, while for $R\gg1$ they are only virtually excited and generate a random four-fermion interaction of rank $M$, which is generically chaotic for $M\geq2$. Using exact diagonalization at infinite temperature, we follow the emergence of random-matrix statistics between these limits with unfolded spectral statistics and the Krylov complexity of the thermofield-double state. For complex couplings the level-spacing distribution evolves from Poisson to Gaussian-unitary-ensemble (GUE) statistics, the spectral form factor develops a correlation hole and a ramp, and the Krylov complexity evaluated on the unfolded spectrum develops a peak above its late-time plateau whose height reaches the GUE value. All diagnostics locate the onset around $R\approx0.3$, and adding boson modes moves it to smaller $R$. This delineates the parameter regime in which a cavity-QED realization would display random-matrix spectral correlations. On the other hand, real couplings give the same crossover towards Gaussian-orthogonal-ensemble (GOE) statistics, and along an interpolation between the two the peak height follows the mean gap ratio, so the peak identifies the symmetry class as well as the onset of chaos. Without unfolding the peak stays about 20$\%$ below its random-matrix value in both classes, which we trace to the non-semicircular density of states, and at larger $R$ the boson-number bands slow down the spreading of the state in Krylov space.

Modern deep neural networks treat parameters as static floating-point matrices stored in physical memory, incurring Von Neumann memory bottlenecks and representation collapse. We formulate Orbital Error Dynamics (OED), an analytical framework wherein synaptic weights are not stored masses (O(W)), but transient topological resonances (O(1)) derived procedurally from the complex quadratic polynomial map z_{n+1} = z_n^2 + c. We introduce the Bent Sine Wave Hypothesis, demonstrating that non-equilibrium living systems emerge when harmonic waves curl inward through environmental drag toward the cardioid cusp (c = 1/4). We define the Observer Horizon Geometry in parameter space, identifying interior resonance shoulder loci X_upper = (0.25, +0.18) and X_lower = (0.25, -0.18) between the fixed-point basin and the true boundary at c = 0.25 +/- 0.50i. To escape non-convex stagnation without loss zeroing, we introduce a heavy-tailed Biomimetic Perturbed Jump Operator (Omega_tunneling) inspired by mammalian fertilization zinc sparks. We further couple an enteric-cranial Dual-Brain architecture shielded by adaptive CD4+ regulatory immune gating (M_CD4), and project the 4-nucleotide genetic basis (A, T, C, G) across quadrants in C. Multi-seed empirical validation on the Two-Moons manifold (5 seeds, 80/20 train/test split, 32x32 grid, zero test-time updates, zero label leakage) demonstrates that procedural parameterization from a 24-byte coordinate seed achieves 77.67% +/- 5.35% clean test accuracy (within an 8.00-point paired difference of an unconstrained gradient baseline at 85.67% +/- 5.35%, 95% CI: [-1.07%, 17.07%]) and 71.33% +/- 3.80% under distribution shift (N(1.2, 0.4)), alongside conceptual equivalence with an analog optical co-processor.

We explicitly construct the additional integrals of motion, ensuring maximal superintegrability of the discrete-time rational Ruijsenaars-Schneider model. Using them, we investigate the algebraic aspects of superintegrability in both continuous- and discrete-time settings. In particular, we determine the complete structures of the polynomial symmetry algebras associated with both the rational Ruijsenaars-Schneider model and its discretization. We demonstrate that discretization leads to a nontrivial deformation of the continuous symmetry algebra with respect to the discretization parameter, thereby extending recent analogous results from the rational Calogero-Moser system to its relativistic generalization.

Precise conditions are provided for the existence and criticality of Turing bifurcations in a general class of activator-inhibitor reaction-diffusion equations on a one-dimensional infinite domain. The class includes generalised Schnakenberg and Brusselator models, as well as other models with cubic autocatalytic nonlinear terms. Previous numerical work suggests the existence of a bifurcation structure containing localised patterns due to the so-called homoclinic snaking mechanism. This paper provides explicit calculations to justify those results. Two distinct scalings of parameters that lead to tractable normal-form coefficients are considered in the limit that the diffusion ratio $\delta \to 0$. First, under a small-parameter scaling, the Turing bifurcation is shown to be always subcritical. Second, a large-parameter scaling reveals the Turing bifurcation to be supercritical, leading, by continuity, to the existence of a codimension-two degenerate bifurcation. The sign of a 5th-order normal form coefficient is also computed, which is shown to have the correct sign for the local birth of homoclinic snaking. For the case of the Brusselator, two such codimension-two points can be found explicitly, as can the leading-order expression for the Maxwell point, in a parameter wedge about which localised patterns emerge. Numerical results are found to be consistent with the theory.

We show that negative feedback alone generates the diversity of self-organized shapes usually attributed to scale-dependent activation-inhibition. For a broad class of kernels, pairwise competition models only generate hexagonal spot arrays. Higher-order terms eliminate this restriction and promote stripes and gaps. Combining individual-based simulations and nonlinear analysis, we derive the pattern-selection thresholds and construct the full state diagram, including an unusual spots-stripes-spots sequence. Spot patterns are thus not a reliable indicator of proximity to a tipping point.

We consider a hierarchy of models that can be obtained from the original Nyzhnyk system by imposing conditions on parameters, introducing potentials or pseudopotentials, performing limiting processes with respect to a scaling parameter, applying differential substitutions and interpreting parts of independent and/or dependent variables as complex or real. Therefore, among the Nyzhnyk models, one can distinguish between symmetric and asymmetric, dispersive and dispersionless, standard and modified, as well as real, complex, mixed and specific models. One can also distinguish single partial differential equations or systems of such equations, as well as linear or nonlinear Lax representations in the dispersive or dispersionless cases, respectively. For each specified model, we compute the maximal Lie invariance pseudoalgebra and, in the symmetric case, find the point- and contact-symmetry pseudogroups using the megaideal-based version of the algebraic method. It is shown that relations between these pseudoalgebras and between these pseudogroups are induced by relations between the corresponding models. Defining (finite-dimensional) subalgebras are singled out in all these pseudoalgebras in the symmetric and specific cases. Based on the established correspondences between the dispersionless and dispersive models, we completely classify one- and two-dimensional subalgebras of the maximal Lie invariance pseudoalgebra of the (dispersive symmetric potential) Nyzhnyk equation and one-dimensional subalgebras of the maximal Lie invariance pseudoalgebra of its linear Lax representation.

Integer maps are discrete dynamical systems defined on the natural numbers. In this paper, we investigate the dynamics of the shifted Alladi-Erd{\H o}s map, a one-parameter family of integer maps in which each composite number is mapped to the sum of its prime factors, while each prime is mapped to $n+A$, where $A \in \mathbb{N}$ is a fixed shift parameter. By systematically exploring the parameter space for $2 \leq A \leq 10^5$, we uncover a rich bifurcation structure characterized by the emergence, disappearance, and reorganization of attractor cycles as the shift parameter varies. We find that although several attractors may coexist for a given value of $A$, for most values of $A$, almost all natural numbers belong to the basins of only two dominant attractors. Using elementary number theoretic arguments, we explain the observed bifurcation diagram. We further characterize the attractor cycles and quantify the distribution of their basin sizes across the parameter space. Our results reveal an unexpectedly rich landscape of arithmetic dynamics arising from a remarkably simple integer map. %and provide a systematic characterization of how attractors and their basins evolve as shift parameter is varied.

This data article provides the datasets, symbolic derivations, and scripts used to reproduce master diagrams, stationary-phase spectra, windowed first-order coherence functions, and quantum-noise stability maps for strongly chirped dissipative solitons of the cubic-quintic complex Ginzburg-Landau equation in normal and anomalous group-delay dispersion regimes. The repository includes node-regularized normal-dispersion spectra and energies; small-parameter expansions of the branch roots; cavity-map gain-loss update relations; Airy uniformization at the normal-dispersion spectral edge; anomalous-dispersion spectra and coherence calculations; and processed tables for plotting and stability analysis. OriginLab projects are accompanied by open-format .csv/.txt numerical tables to support reuse without proprietary plotting software. Data and code repository: https://doi.org/10.5281/zenodo.22690899.

We obtain a large number of exact hyperbolic, trigonometric, and periodic solutions in terms of Jacobi elliptic functions as well as algebraic solutions with a power law tail of the integrable local Fokas-Lennels equation and integrable nonlocal Fokas-Lennels equation. Further, we consider a one-parameter family of generalized Fokas-Lenells equations and obtain a few of their exact solutions.