Papers

We establish rogue wave patterns associated with the fourth Painlevé equation $(\mathrm{P}_{\mathrm{IV}})$ in the multi-component nonlinear Schrödinger and Hirota equations. The generalized Hermite and generalized Okamoto polynomials arise in representations of rational solutions of $\mathrm{P}_{\mathrm{IV}}$, and we show that their roots determine two classes of rogue wave patterns when one of the internal parameters of rogue wave solutions is large. Specifically, the generalized Hermite polynomials arise from rogue wave solutions represented by Schur-polynomial determinants with consecutive indices, whereas the generalized Okamoto polynomials arise from analogous determinants with index jumps of three. Numerical examples for both equations agree with the predictions.

We present a classification of supersymmetric bilinear KdV-type equations admitting unconstrained three-super-soliton solutions. Extending Hirota's classical three-soliton criterion to the supersymmetric setting, we derive the complete bosonic and fermionic compatibility conditions governing the existence of three-super-soliton solutions. We prove that every supersymmetric bilinear KdV-type equation possesses unconstrained one- and two-super-soliton solutions, whereas three-super-soliton solutions exist only when eight integrability conditions are satisfied. These conditions provide a supersymmetric analogue of Hirota's classical integrability criterion and naturally recover the fermionic relations previously introduced by Carstea, revealing their structural origin. As a consequence, we classify the supersymmetric extensions of Hirota bilinear KdV-type equations and show that only a subset of Hietarinta's classical classification remains valid in the unrestricted supersymmetric framework.

We examine a 4-dimensional generalization of predator-prey models describing interactions among four species. We show that these systems exhibit a rich variety of asymptotic behaviors, with solutions that are either integrable, leading to non-sustainable dynamics characterized by unbounded population growth or collapse, or chaotic, where all species coexist over time. Moreover, we review the integrable cases, which are classified into two-parametric families, and examine three-parametric families that behave similarly to the known two-parametric integrable cases. We also analyze chaotic Lotka-Volterra models in terms of their Lyapunov exponents and Poincare' surfaces of section, which turn out to reveal underlying structures that are compatible with what is expected from sustainable dynamics.

Integrable partial differential equations, such as the Korteweg–De Vries (KdV) equation, admit infinite hierarchies of commuting higher symmetries. Their symmetry reductions give rise to integrable finite-dimensional dynamical systems solvable in terms of Abelian functions. This underlies the finite-gap integration method for the KdV equation, introduced by S.P. Novikov and subsequently extended to many other integrable systems. In this paper, I introduce the concept of reduction ideals and a new class of reductions for integrable differential–difference equations, leading to integrable finite-dimensional systems in both commutative and noncommutative settings. This class includes periodic reductions, symmetry reductions, and a broad family of other reductions defined in terms of the associated Lax operator. The reduction constraints define integrable maps that enable solutions of the reduced systems to be extended to solutions of the corresponding differential-difference equations. The construction is illustrated using the Volterra hierarchy.

We construct a class of integrable chiral fermionic quantum field theories with interactions that depend on both space and time. The construction starts from a unitary, regular, difference-form two-body scattering matrix satisfying braiding unitarity and the Yang–Baxter equation. Each freely moving right- or left-moving particle carries a characteristic coordinate that remains constant along its trajectory, and the scattering matrix depends on differences of these coordinates. The many-body amplitudes form a discrete connection over the coordinate-ordering sectors. The Yang–Baxter equation gives local path independence, while periodic winding on a spatial circle produces quantum Knizhnik–Zamolodchikov (qKZ) difference equations expressing global flatness. Once a reference amplitude satisfying the exchange and qKZ constraints is specified, these transport equations determine the exact wavefunction in every fixed-particle-number sector. Independent reparametrizations of the right- and left-moving characteristics give couplings that depend on both space and time. For the chiral $SU(2)$ Gross–Neveu realization, we solve the qKZ equations for this full characteristic family using off-shell Bethe vectors. The resulting amplitudes are written in the characteristic labels, and their pullbacks give exact many-body wavefunctions for the corresponding space-time-dependent couplings. The two-body amplitude is a ratio of Gamma functions that resums the phase shifts from successive scattering processes. Equal-slope affine maps give a spatially homogeneous specialization in which the coupling depends only on time, and the Bethe ansatz with constant coupling is recovered when the qKZ step vanishes.

We develop a matrix pseudodifference operator framework that places theAdler-Pfaff and Pfaff-Toda hierarchies in a common algebraic setting. The bi-infinite Adler-Pfaff hierarchy then becomes a $2\times2$ matrix pseudodifference Lax hierarchy. We introduce the Matrix Pfaff-Toda hierarchy by dressing the matrix Laurent algebra $M_{2N}\bigl(\mathbb{C}[\mathcal{S},\mathcal{S}^{-1}]\bigr)$ using a Pfaff-Toda splitting of the matrix pseudodifference algebra. In the one-component case, its diagonal sector recovers the continuous Pfaff-Toda hierarchy, while the off-diagonal directions supply the missing odd flows, and the powers of a single bare operator reconstruct the full Adler-Pfaff hierarchy. For general $N$, on the regular factorization locus, the multicomponent Pfaff-Toda tau-functions of Savchenko and Zabrodin realize the commuting diagonal sector, with its dressing equations derived directly from the fermionic bilinear identity. Reductions to the block-diagonal and scalar cases recover multicomponent and scalar $2D$ Toda, identifying the even Pfaff hierarchy with an anti-diagonal Toda subhierarchy, and yielding the Krichever-Zabrodin C-Toda hierarchy.

We formulate the linearization of the fourth Painlevé equation as a linearized Hamiltonian system and transform the scalar variational equation to normal form, with the first-derivative term eliminated. From the canonical fundamental solutions of the Jimbo–Miwa Lax pair, we construct quadratic expressions $Q_{IV}$ and $R_{IV}$ that satisfy the identity $\partial_x^2Q_{IV}-U_{IV}Q_{IV}=\partial_\lambda R_{IV}$. Here $U_{IV}$ is the coefficient in the normal-form equation. A rapid-decay cycle with a vanishing boundary term yields an integral solution of the linearized equation. Every nontrivial such cycle has asymptotic branches in at least two distinct Stokes sectors; otherwise the integral vanishes by Cauchy's theorem. A second rapid-decay homology class gives another integral solution. The Wronskian of the two integral solutions depends holomorphically on the monodromy data, so a single nonzero value implies generic linear independence. For one set of complex data, direct quadrature gives numerical evidence of a nonzero Wronskian. We also derive integral formulas for variations of the Stokes multipliers, the connection matrix, and the local monodromy, and verify the contour representation by substitution into the linearized equation.

We consider the application of Sundman transformations to the first four Painlevé equations $P_I$-$P_{IV}$, with $P_V$ and $P_{VI}$ being left for future work. Since the Painlevé equations are nonautonomous, this requires first of all the construction of an equivalent autonomous system, a process which we discuss in detail, and which we characterize as being inverse to the use of Lie symmetries to obtain a reduction of order of an autonomous system. Secondly, we apply a Sundman transformation to this autonomous system, and identify the resulting equation with a corresponding nonautonomous equation. We give an explicit mapping between solutions of this nonautonomous equation and the original Painlevé equation. Interesting special cases of the nonautonomous equation found by Sundman transformation include, for $P_{II}$, an equation with coefficients expressed via Airy functions, and for $P_{IV}$, an equation with coefficients expressed via Weber-Hermite functions. We also discuss an alternative autonomisation process. Finally, a variety of possible directions for future work are identified, along with interesting open problems.

We study the rational quadratic differential equation \[ \frac{\mathrm{d}y}{\mathrm{d}x} = \frac{ay^2+by+cx}{y^2}, \qquad a,b,c\in\C, \] under the non-degeneracy assumptions \[ c\neq 0,\qquad 2ay+b\not\equiv 0. \] Equivalently, after clearing the denominator, we consider the polynomial vector field \[ \dot{x}=y^2,\qquad \dot{y}=ay^2+by+cx. \] We give a complete, directly checkable classification of its Darboux curves. If $a=0$, no non-constant Darboux polynomial exists. If $a\neq 0$, a non-constant Darboux polynomial exists if and only if \[ c=-ab\qquad\text{or}\qquad c=-2ab. \] In these two cases the unique irreducible Darboux polynomials, up to non-zero constant multiples, are respectively \[ y-ax,\qquad y^2-2bx. \] Consequently every non-constant Darboux polynomial is a non-zero constant multiple of a positive integral power of the corresponding irreducible factor. We then place this classification in the framework of Riccati (R-)integrability and rational potentials. Excluding the exceptional branch $c=-2ab$, the analytic integrability theorem identifies $c=-ab$ as the branch admitting a global Riccati-type analytic first integral. For $c=-2ab$, we compute the first four transverse variational groups along the transformed Darboux divisor over the rational function field. Their dimensions are $1,2,3,4$; the fourth is the full group of invertible fourth-order transverse jets. We prove that this exceptional branch admits no Riccati first integral and hence is not R-integrable. Finally, we relate the branch $c=-ab$ to the Quartic Inverse Riccati (QIR) class and discuss an invariant-based classification problem for quartic Abel equations.

We survey a geometric approach to the Navier-Stokes equations developed by the author over the past fifteen years. The central object is a Riemannian space of fourteen dimensions whose metric is Ricci-flat precisely on the solutions of the Navier-Stokes system, so that the behaviour of a viscous incompressible fluid becomes a property of curvature. The space decomposes into a flat six-dimensional part and two dual quadruples of coordinates, Eulerian and Lagrangian, which exhibits the two classical descriptions of a fluid as two halves of a single geometry. We also recall the associated six-dimensional metric, whose integrability condition is the incompressibility of the fluid; the link with the projective geometry of second-order ordinary differential equations in the sense of E. Cartan; and the use of Cartan's invariants to construct a three-dimensional analogue of the Taylor-Green vortex. An application of the same construction to the equations of rotation of a rigid body is indicated.

We study the two-field nonlinear dispersive system \[ u_t=u_{xxx}+6uu_x,\qquad v_t=v_{xxx}+6(uv)_x, \] viewed simultaneously as a coupled wave equation, as the tangent covering of Korteweg–de Vries (KdV), and as a Hamiltonian flow on a tangent Poisson manifold. The main purpose is to make these viewpoints interact at theorem level. First, we prove that the Magri Poisson pencil of KdV admits a complete tangent lift to an explicit compatible pair of matrix Hamiltonian operators. The corresponding recursion operator has triangular tangent form and generates the lifted KdV hierarchy. Second, we identify a five-dimensional point-symmetry algebra together with the infinite hierarchy of tangent generalized symmetries. Third, reduction by the traveling-wave subgroup produces a four-dimensional Hamiltonian system that is the complete tangent lift of the scalar KdV profile dynamics. On the nonsingular elliptic locus, the base profile is written in Weierstrass form and every tangent traveling wave is classified explicitly by derivatives with respect to the energy and integration constants.

We study the polynomial generalized Korteweg–de Vries family \[ u_t+P_m(u)u_x+\kappa u_{xxx}=0,\qquad \kappa\neq0, \] where $P_m$ is a real polynomial of degree $m\ge1$. The purpose is not to generate isolated closed-form waves, but to identify the structural threshold at which pole-type Painlevé behavior and elliptic traveling-wave geometry are simultaneously lost. A Weiss–Tabor–Carnevale dominant-balance calculation gives the universal principal exponent $p=-2/m$. Hence only $m=1$ and $m=2$ can possess a principal Laurent branch with integer pole order; their resonance sets are respectively $\{-1,4,6\}$ and $\{-1,3,4\}$. These two sectors are, after affine and Galilean transformations, the KdV and modified KdV equations, so their compatibility conditions and complete-integrability structures are inherited from the classical hierarchies. In contrast, every degree $m\ge3$ has a fractional leading exponent and therefore fails the strong WTC pole criterion at the first step. The results provide a concise bridge between singularity analysis, Hamiltonian form, algebraic-curve genus and exact nonlinear waves.

We prove that a scalar tau-symmetric Hamiltonian deformation of the Riemann hierarchy in generalized standard form is uniquely determined by the nonzero coefficient of $u_x^2$ and the coefficients of $u_{xx}^g$, $g\ge2$, in its first Hamiltonian density. This determines the original hierarchy up to normal Miura transformations. Combined with the construction of Hodge-class double ramification hierarchies by Buryak–Rossi, this establishes Hodge universality at nonzero dispersion.

We study current fluctuations in a temporally random integrable quantum circuit. Commutativity reduces every drive history exactly to its layer composition, turning annealed fluctuations into a competition between current gain and the large-deviation cost of rare compositions. For each fixed composition, homogeneous thermodynamic Bethe ansatz dressing supplemented by ballistic fluctuation theory yields the conditional current statistics. Their annealed large-deviation contraction predicts a first-order switch between two dominant history classes, terminating on a line of regular cusp endpoints. Finite-time analysis shows how the switch is rounded. Thus temporal randomness can act as an emergent order-parameter-like coordinate in trajectory space.

Interactions of two solitary waves with an up and a down orientation in the modified Korteweg-de Vries equation are shown to produce rogue-like waves. For waves that asymptotically vanish, the maximum ratio between the height of the interaction profile and the height of the tallest incoming wave is 2.41 when the waves have approximately equal speeds, and this ratio decreases to 2 when the speed ratio is 1.5. For waves that approach a non-zero constant at infinity, the same ratio reaches a maximum of 2.65 when the speed ratio of the waves is 6.32.

Quantum reservoir computing (QRC) is a machine learning approach which employs the internal dynamics of a physical system (the reservoir) to encode and process information. In this work, we explore the use of dual-unitary circuits in a brickwork architecture as a platform for QRC, well suited to current noisy intermediate-scale quantum devices. Dual-unitary circuits present both practical and conceptual advantages. Our results indicate that, under appropriate conditions, dual-unitarity can lead to an enhanced regime of operation: we numerically verify that it improves memory effects and nonlinear processing, and shields against finite-shot noise, mitigating exponential concentration. Moreover, dual unitarity offers an intuitive picture of how operator dynamics gives rise to memory and nonlinear processing in circuit-based reservoirs.

The Julia set is defined by the closure of repelling periodic points in chaotic systems. Why does this structure not appear in integrable systems? In this paper, we address this question by demonstrating the existence of the closure of divergences of the periodic equations, which we designate as the "anti-Julia set." We also call the sets before taking the closures the pre-Julia and pre-anti-Julia sets, respectively. We illustrate the transition mechanism by considering a complex map that interpolates between integrable and non-integrable dynamics, by introducing a real deformation parameter $a$. For $0<a\le 1/2$, we show that the Julia set and the anti-Julia set coincide in the complex plane, although the pre-Julia and pre-anti-Julia sets remain completely disjoint. At the integrable limit $a\to 0$, these two dual structures undergo a critical collision and subsequent annihilation, reminiscent of a cosmological Big Crunch. When $1/2<a<1$, on the other hand, the boundary of the pre-anti-Julia set includes the Julia set. We analytically characterize these phenomena, focusing in particular on the asymptotic behavior of the pre-anti-Julia set as it approaches the integrable limit, and provide numerical visualizations that elucidate the underlying mechanisms of this Big Crunch phenomenon.

This paper deals with two kinds of Bäcklund transformations for scalar hyperbolic partial differential equations. We prove that both these types of transformations map solutions of a Darboux integrable equation into solutions of, generally speaking, another but also Darboux integrable equation. The latter fact can be used to roughly check the completeness of a list of Darboux integrable equations. To illustrate this, we apply the above transformations to several equations from a well-known list of Darboux integrable equations and, as a result, obtain a Darboux integrable equation which is absent in this list, but is already known at present. As a generalization of the last equation, we construct a family of Darboux integrable equations that is parametrized by three arbitrary functions, each of which depends on two arguments. This family is probably new.

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.

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.

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.

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.

Symmetry reductions of partial differential equations (PDEs) inherit more geometric structures than other types of reductions: invariant conservation laws, variational structures, and, under suitable conditions, Hamiltonian-type structures of the original model descend to the reduced model through the mechanism of invariant reduction. This paper develops computational algorithms that carry out such reductions explicitly. We use the interpretation of variational $p$-forms as conservation laws of an enlarged system — the (degree-shifted) tangent system — consisting of the original equations together with their linearizations, in which the perturbation variables are treated as anticommuting. This allows us to formulate the algorithms using well-known concepts from the theory of conservation laws, naturally adapted to the graded-commutative setting. The reduction of conservation laws, variational $1$-forms, and presymplectic structures thereby becomes a single algorithmic procedure. We present (i) a homotopy-based reduction algorithm for systems of evolution equations, implemented in Maple; (ii) a descent reduction algorithm applicable to general $\ell$-normal systems for $p>0$; and (iii) a simple reduction algorithm available for point symmetries under suitable conditions, including the description of reductions in terms of systems involving fewer independent variables. Examples include nonlinear evolution equations in one and two spatial dimensions, the Laplace equation, the incompressible Euler equations, and the cotangent system of Pavlov's equation.

We construct isomonodromy transformations for linear $q$-difference systems of the form $Y(qz)=A(z)Y(z)$, where $A(z)=A_0+z A_1$ has diagonal leading coefficient. These transformations shift eigenvalues of the leading coefficient $A_1$ together with roots of $\det A(z)$. They lift compatibility to the right eigenpairs of $A(z)$, yielding a discrete local tau function. The resulting deformation equations preserve the Birkhoff connection matrix, and reduce in their $q\to 1$ limit to the isomonodromic deformation of a meromorphic connection on $\mathbb P^1$ with an irregular singularity of Poincaré rank one at $\infty$.

The exact N-soliton solutions are derived for the coupled negative first order Klein-Gordon (CNKG) equation subject to vanishing boundary conditions by using the inverse scattering method via Gelfand-Levitan-Marchenko equation. Based on the zero-curvature representation, the conservation laws, integrals of motion, and Hamiltonian structure of the aforementioned coupled nonlinear equations are constructed. The Jost functions and their analyticity properties for the Manakov spectral problem are recalled. The integral equations for the eigenfunctions are then used to formulate the Gel'fand-Levitan-Marchenko equations. Solving these equations establishes a direct correspondence between the kernel functions and the potential, yielding the general N-soliton expressions.

Fushchych and Shtelen defined approximate symmetry by expanding the solution of a perturbed differential equation in the small parameter, replacing the original equation by a triangular system for the expansion coefficients. We extend this construction to arbitrary order. A coefficient-map formulation yields the order-$p$ Fushchych–Shtelen (FS) system and establishes, under explicit hypotheses, a correspondence with the Baikov–Gazizov–Ibragimov (BGI) method. The two frameworks are compared for a perturbed cubic wave equation. Up to equivalences, we classify the nonlinearities admitting an FS dilation through second order. The classification consists of a generic power family and logarithmic branches at exceptional exponents. Among perturbations nontrivial at first order, vertical BGI continuations select a proper subfamily of one FS branch, and a second common branch arises when the perturbation first enters at second order; the two methods are related but not interchangeable. Joint reduction by the Lorentz algebra and the inherited FS dilation integrates all resulting FS normal forms; the explicit solutions of the 1989 Fushchych–Shtelen letter are recovered as the first-order members of this reduction. In the part common to both frameworks, the logarithmic corrections arise from the expansion of a power law with a perturbation-dependent exponent. Finally, periodic travelling waves are used to examine the long-scale validity of the FS expansion. The first correction is obtained by quadrature, the second-order secular structure is isolated, and phase renormalization produces bounded order-consistent approximations. For integer-power members of the generic family, the amplitude dependence of the corrected wavenumber agrees with the scaling weights obtained from the symmetry classification.

In this paper, we construct the lattice Kadomtsev-Petviashvili (KP) type eigenfunction equations. A homogeneous nonlocal $\bar{\partial}$ problem is considered, from which we are able to define the eigenfunction of the Lax pair of the lattice KP equation. The eigenfunction together with its expansions at infinity and at a finite analytic point provide formulations of the lattice modified KP equation, the lattice Schwarzian KP equation and the Nijhoff-Quispel-Capel KP (NQC-KP) equation. We also consider an inhomogeneous nonlocal $\bar{\partial}$ problem. It defines the eigenfunction of the Lax pair of the lattice modified KP equation. This eigenfunction generates a direct formulation for the NQC-KP equation, which is different from the previous ones. Explicit solutions of these equations are obtained, from which we can see the difference of the different formulations for same equations.

We present a novel approach to computing the Nijenhuis torsion and Frölicher-Nijenhuis brackets of recursion operators for symmetries, based on their full-fledged forms introduced by Jahnová and Vojčák (2024). In contrast to the conventional approach, which represents recursion operators as maps between shadows of nonlocal symmetries within a given covering, our method allows the Nijenhuis torsion to be computed directly from its defining formula, without requiring any additional mathematical constructions. The same framework can also be used to compute the Frölicher-Nijenhuis bracket of recursion operators and to verify their compatibility directly in their full-fledged forms. The procedure is illustrated by several examples of full-fledged recursion operators for symmetries of differential equations in four independent variables, including a new recursion operator for the four-dimensional universal hierarchy equation that, to the best of our knowledge, has not previously appeared in the literature. A notable feature of this new full-fledged recursion operator is that its shadow component does not appear to admit a reasonable representation in any of the standard conventional forms. Nevertheless, our approach allows us to prove directly that the Nijenhuis torsion of this operator vanishes. For the examples considered, the corresponding Frölicher-Nijenhuis brackets also vanish, confirming the compatibility of the respective pairs of recursion operators. This demonstrates that full-fledged forms provide an effective framework for studying the hereditary and compatibility properties of recursion operators, including highly nonlocal and/or multidimensional cases that are difficult to handle by existing methods.

A class of moving boundary problems of Stefan-type for an extension of the 2+1- dimensional Gardner equation is shown to be amenable to exact solution by application of a Ermakov- Painlevé II symmetry reduction.

A new family of Darboux integrable partial differential equations is constructed. The minimal orders of integrals in both characteristics can simultaneously be arbitrarily high for equations within this family. This suggests that the complete list of Darboux integrable equations may contain substantially more equations than are currently known.