Yajima-Oikawa Equation in Integrable Systems

• The Yajima-Oikawa equation is an integrable long-wave-short-wave resonance interaction model featuring a complex short-wave envelope coupled with a real long-wave field.
• It exhibits rich integrability properties including Lax pairs, Hirota bilinear forms, and Riemann-Hilbert solutions, which underpin analyses of solitons, rogue waves, and discrete systems.
• It finds applications across fluid and plasma physics, nonlinear optics, and geometric curve flows, linking continuous, discrete, and multicomponent models.

Searching arXiv for Yajima-Oikawa papers to ground the article in current arXiv literature. arxiv_search(query="Yajima-Oikawa equation integrable discretization multicomponent rogue waves geometric realization", max_results=10) arxiv_search({"query":"Yajima-Oikawa equation integrable discretization multicomponent rogue waves geometric realization","max_results":10}) I’m unable to invoke the arXiv search tool in this interface, so I will rely on the arXiv records explicitly provided. The Yajima-Oikawa equation is an integrable long-wave-short-wave resonance interaction system, often written in one normalization as

$i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$

with $q(x,t)$ a complex short-wave envelope and $u(x,t)$ a real long-wave field. In the literature surveyed here, it is also described as a deformation of the Zakharov equation which models the propagation of ion sound waves subject to the ponderomotive force induced by high-frequency Langmuir waves, and as one of the first examples of integrable long-short wave interaction models (Myrzakulova et al., 2022, Wang et al., 21 Jul 2025).

1. Canonical forms and physical interpretation

A standard dimensionless $(1+1)$-dimensional form used in the recent rogue-wave literature is

$\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$

where $S$ is the complex short-wave field and $L$ is the real long-wave field (Chen et al., 2017). Another normalization used in the long-short wave and spin-system literature is

$i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$

with the same division into a complex short-wave envelope and a real long-wave field (Myrzakulova et al., 2022). A further normalization, emphasized in the Yajima-Oikawa-Newell setting, is

$i S_t+S_{xx}-LS=0,\qquad L_t=2(|S|^2)_x,$

obtained from the larger YON model by choosing $\alpha=0$ and $q(x,t)$0 (Caso-Huerta et al., 2022).

These forms encode the same basic structure: the short wave satisfies a Schrödinger-type evolution with a self-consistent long-wave potential, while the long wave is driven by the spatial derivative of the short-wave intensity. The papers place this interaction in fluid physics, plasma physics, nonlinear optics, and long wave-short wave resonance theory (Myrzakulova et al., 2022, Caso-Huerta et al., 2022). This suggests that several sign and scaling conventions are in active use, but the distinction between a complex short-wave field and a real long-wave field is stable across the literature.

2. Integrability and hierarchy structure

The Yajima-Oikawa equation is treated throughout the cited literature as a fully integrable system. It admits Lax representations, zero-curvature formulations, hierarchy embeddings, and inverse-scattering-type solution methods (Myrzakulova et al., 2022, Wang et al., 21 Jul 2025). In the KP-hierarchy viewpoint, it is repeatedly identified with the $q(x,t)$1 constrained KP hierarchy, in contrast to the nonlinear Schrödinger equation, which corresponds to $q(x,t)$2 (Chen et al., 2015, Chen et al., 2017). In matrix form, the same hierarchy viewpoint extends to $q(x,t)$3-dimensional and bidirectional constrained KP systems, where matrix generalizations of Yajima-Oikawa appear alongside Davey-Stewartson and Melnikov members (Chvartatskyi et al., 2013).

A distinct integrable-hierarchy construction arises from deformations of binary Darboux transformations. In that framework, the deformation

$q(x,t)$4

produces self-consistent source extensions, and a $q(x,t)$5-dimensional Yajima-Oikawa-type system emerges from the deformed matrix pKP hierarchy after setting $q(x,t)$6 and introducing

$q(x,t)$7

(Chvartatskyi et al., 2015). The resulting source-extended equation is not the standard source-free Yajima-Oikawa system; it is explicitly a self-consistent source extension.

The inverse-scattering and Riemann-Hilbert literature further sharpens the integrability picture. A recent analysis formulates a $q(x,t)$8 Riemann-Hilbert problem in terms of two reflection coefficients, derives exact one-soliton solutions, and gives long-time asymptotics in the Zakharov-Manakov region by the Deift-Zhou nonlinear steepest descent method (Wang et al., 21 Jul 2025). That paper treats the long-time analysis as the first rigorous long-time asymptotic analysis of the Yajima-Oikawa equation with continuous spectrum.

3. Higher-dimensional, multicomponent, and generalized versions

One of the most important modern reinterpretations is the identification of $q(x,t)$9-dimensional Yajima-Oikawa systems inside the Davey-Stewartson framework. Starting from the elementary Davey-Stewartson flow

$u(x,t)$0

the linear change of variables

$u(x,t)$1

yields

$u(x,t)$2

with $u(x,t)$3 and arbitrary $u(x,t)$4; the special case $u(x,t)$5 is noted as the Maccari system (Tsuchida, 2019). A recurrent point of clarification in this literature is that the $u(x,t)$6-dimensional Yajima-Oikawa system appears as an equivalent form of the first elementary Davey-Stewartson flow under a linear transformation of the independent variables, not as a nonlinear reduction (Tsuchida, 2019, Tsuchida, 2020).

The multicomponent $u(x,t)$7-dimensional generalization couples one real long-wave field $u(x,t)$8 to $u(x,t)$9 complex short-wave fields $(1+1)$0: $(1+1)$1

$(1+1)$2

with arbitrary real nonlinearity coefficients $(1+1)$3 (Kanna et al., 2014). The Painlevé analysis in that work shows integrability for arbitrary real $(1+1)$4, including positive, negative, and mixed-sign choices.

A $(1+1)$5-dimensional multicomponent form used in the rational-solution and dark-soliton literature is

$(1+1)$6

or equivalently with $(1+1)$7 in the long-wave equation, depending on notation (Chen et al., 2014, Chen et al., 2015). In these papers, the $(1+1)$8-dimensional multicomponent system is obtained by removing the $(1+1)$9-dependence through explicit parameter constraints.

Generalized long-short wave models also embed the Yajima-Oikawa equation as a reduction. The Yajima-Oikawa-Newell model

$\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$0

is integrable for any real $\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$1, and it reduces to the Yajima-Oikawa equation when $\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$2 and $\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$3 (Caso-Huerta et al., 2022). In the matrix constrained-KP literature, a matrix $\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$4-dimensional Yajima-Oikawa system appears in the form

$\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$5

with $\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$6 and $\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$7 matrix-valued (Chvartatskyi et al., 2013).

4. Exact solutions and nonlinear wave phenomena

The solution theory of the Yajima-Oikawa equation is unusually rich. In the multicomponent $\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$8-dimensional setting, Hirota bilinearization yields a general bright $\mathrm{i}S_t-S_{xx}+SL=0,\qquad L_t=-4(|S|^2)_x,$9-soliton solution in Gram determinant form, together with explicit one-, two-, and three-soliton analyses (Kanna et al., 2014). That work emphasizes two kinds of energy-sharing collisions in the short-wave sector, while the long-wave component always exhibits elastic collision. In the mixed bright-dark setting, a general mixed $S$0-soliton solution for the $S$1-component system is obtained in Gram determinant form, and the collision analysis shows that inelastic collision can only take place among short-wave components when at least two short-wave components have bright solitons; the dark short-wave components and the bright long-wave component always undergo usual elastic collision (Chen et al., 2015).

The dark-soliton sector has a particularly rigid structure. For the two-dimensional and one-dimensional multicomponent systems, general $S$2-dark soliton and $S$3-dark-dark soliton solutions are constructed in both Gram type and Wronski type determinant forms by KP-hierarchy reduction (Chen et al., 2015). That literature shows that dark-dark soliton collisions are elastic, with no energy exchange among solitons in different components, and that stationary bound states can exist up to arbitrary order, whereas moving bound states are much more restricted.

Breathers, rational solutions, and rogue-wave families have been developed in parallel. Through KP-hierarchy reduction, the one-dimensional and two-dimensional Yajima-Oikawa systems admit two kinds of breather solutions, with homoclinic orbit and dark soliton limits, and long-wave limits producing rational and rational-exponential solutions containing lumps, line rogue waves, solitons, and mixed cases (Chen et al., 2017). In the multicomponent rational-solution literature, the same determinant structures produce fundamental lumps, line rogue waves, multi-rogue-wave interactions, and higher-order rogue waves in both $S$4 and $S$5 dimensions, with the short-wave components classified into bright, intermediate, and dark states (Chen et al., 2014). For the $S$6-dimensional system, general high-order rogue waves are likewise classified into bright, intermediate, and dark patterns, and the essential parameter $S$7 controls the pattern because the Yajima-Oikawa system does not possess the Galilean invariance (Chen et al., 2017).

The two-dimensional vector theory extends the rational family further. Using Darboux transformations with functional arbitrariness, one obtains lumps, line rogue waves, semi-rational solutions, higher-order counterparts, and, in the two-component case, appearing lumps, disappearing lumps, and appearing-disappearing lumps or rogue lumps (Ustinov, 2019). This is one of the clearest demonstrations that the vector Yajima-Oikawa system supports finite-lifetime localized structures beyond the scalar case.

Recent Riemann-Hilbert work adds a different exact-solution layer. For the scalar system

$S$8

exact one-soliton solutions arise from a simple zero $S$9 of the scattering data, with smooth $L$0-profile solitons for $L$1 and singular $L$2-profile solitary waves for $L$3, and the long-time asymptotics in the Zakharov-Manakov region $L$4 are obtained by nonlinear steepest descent (Wang et al., 21 Jul 2025).

5. Discrete and semi-discrete Yajima-Oikawa systems

The discrete theory of the Yajima-Oikawa equation is extensive and structurally diverse. One strand begins from the semi-discrete BKP hierarchy. There, an integrable semi-discrete analogue of the one-dimensional coupled Yajima-Oikawa system is proposed by discretizing the spatial variable and retaining continuous time, with Hirota bilinear equations and pfaffian bright and dark soliton solutions (Chen et al., 2015). That paper emphasizes the continuum limit back to

$L$5

The integrability of the corresponding space-discrete system was later strengthened by a Lax-pair representation, an infinite number of conservation laws, and a next higher flow generalizing a modified version of the Volterra lattice (Tsuchida, 2018). That work also clarifies relations to the Ablowitz-Ladik and Konopelchenko-Chudnovsky hierarchies.

A different lattice system arises from a new integrable generalization of the Toda lattice. In that setting, the first flow of the generalized Toda hierarchy yields a discrete long wave-short wave interaction model, and in a suitable continuous limit it becomes

$L$6

which is identified as the Yajima-Oikawa system up to a trivial rescaling and a Galilean transformation (Tsuchida, 2018). That paper explicitly states that its discrete Yajima-Oikawa system is essentially different from the earlier discrete system studied in the Hirota-pfaffian and Lax-pair literature.

Fully discrete formulations also exist. By imposing a symmetry constraint directly on the discrete KP tau function, one obtains integrable discretisations in both space and time of the Yajima-Oikawa system, together with bilinear and Lax representations (Willox et al., 2014). In the language of that paper, the Yajima-Oikawa discretisation is obtained from a two-shift constraint and is presented as part of the discrete constrained-KP hierarchy.

The $L$7-dimensional case has its own semi-discrete developments. By discretizing only one spatial variable $L$8 and retaining the other spatial variable and the evolution variable as continuous, new integrable semi-discretizations of the Davey-Stewartson system and of the $L$9-dimensional Yajima-Oikawa system are obtained, again through Lax-pair constructions and commuting elementary flows (Tsuchida, 2019, Tsuchida, 2020). This literature uses “semi-discretization” specifically for discretization of one of the two spatial variables.

6. Geometric, gauge-equivalent, and source-extended reformulations

The Yajima-Oikawa equation has also been reinterpreted outside the usual PDE setting. In the gauge-equivalence literature, the equation

$i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$0

is shown to be gauge equivalent to the Makhankov-Myrzakulov-II equation, a generalized Heisenberg ferromagnet equation with self-consistent potentials (Myrzakulova et al., 2022). In that setting, the Yajima-Oikawa equation becomes one member of a wider family of integrable long-short wave systems that includes the Newell, Ma, and Geng-Li equations.

The self-consistent-source literature adds another reformulation. In the deformed binary Darboux framework, the Yajima-Oikawa equation is recovered from a deformation of the source-free pKP hierarchy, and the resulting $i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$1-dimensional equation

$i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$2

is explicitly a source-extended Yajima-Oikawa system rather than the standard source-free model (Chvartatskyi et al., 2015). The same framework provides solution formulas and modulated one-soliton examples depending on arbitrary functions of one independent variable.

A recent geometric realization places the system on curves in the $i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$3-sphere $i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$4 that are transverse to the standard contact structure. In that construction, a transverse curve flow

$i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$5

induces the invariant evolution

$i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$6

which is precisely the Yajima-Oikawa system in the variables $i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$7 and $i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$8 (Calini et al., 2024). The same paper constructs the transverse curves associated with Wright’s periodic plane-wave solutions, derives explicit closure conditions, and shows that many closed solutions are torus knots or links in $i q_t+q_{xx}-u q=0,\qquad u_t-(|q|^2)_x=0,$9. It also points toward further study of the full Yajima-Oikawa hierarchy of curve flows, recursion operators, and geometric analogues of Bäcklund or Miura transformations.

In aggregate, these developments present the Yajima-Oikawa equation not as a single isolated model but as a nexus linking constrained KP hierarchies, Davey-Stewartson flows, discrete and semi-discrete integrable systems, Darboux and Riemann-Hilbert techniques, multicomponent collision theory, generalized spin systems, and contact-geometric curve flows.