We formulate a general system of kinetic equations for a non-stationary two-dimensional gas of elastically interacting line solitons and apply it to the description of a soliton gas governed by the Kadomtsev-Petviashvili II (KPII) equation. We then verify the predictions of the kinetic theory in two analytically tractable problems: the oblique interaction of a KPII line soliton with a one-dimensional soliton condensate of the Korteweg-de Vries equation, and the interaction of a trial KPII soliton with a monochromatic KPII soliton gas. In both cases, we compare the analytical results with direct numerical simulations obtained by constructing two-dimensional soliton gases via exact KPII N-soliton solutions for large N, using appropriately chosen random distributions of soliton parameters. The comparison demonstrates excellent agreement, thereby providing strong validation of the proposed kinetic theory of 2D non-equilibrium soliton gases.

We study the formation and interaction of dispersive shock waves (DSWs) in the two-dimensional Toda lattice subject to wedge-type initial conditions, and show that their interaction gives rise to a discrete analog of Mach reflection for dispersive shock waves in discrete systems. The initial jump across each leg of the wedge acts locally as a Riemann problem for the one-dimensional Toda lattice, producing two oblique DSWs whose leading-edge soliton amplitude is determined explicitly by the one-dimensional Whitham modulation theory. The two-dimensional nature of the problem manifests when these oblique DSWs meet along the symmetry axis. We show that, for compressive wedges (i.e., when the initial conditions are such that two oblique DSWs that are generated propagate toward each other), a critical slope q_cr separates two regimes: in the subcritical regime (qq_cr) the interaction is ordinary and produces a localized peak whose amplitude is determined analytically. We also show qualitatively that a similar dichotomy between two regimes exists for expansive wedges (i.e., when the initial conditions are such that the two oblique DSWs propagate away from each other). We confirm all analytical predictions by comparing them with the results of direct numerical simulations. Finally, we show that the continuum limit of the result is consistent with the analogous theory for the Kadomtsev-Petviashvili equation, providing an independent validation of the analytical framework.

The semiclassical limit of the derivative nonlinear Schrodinger equation with periodic initial conditions is studied analytically and numerically. The spectrum of the associated scattering problem for a certain class of initial conditions, referred to as periodic single-lobe potentials, is numerically computed, and it is shown that the spectrum becomes confined to the real and imaginary axes or the spectral parameter in the semiclassical limit. A formal Wentzel-Kramers-Brillouin expansion is computed for the scattering eigenfunctions, which allows one to obtain asymptotic expressions for the number, location and size of the spectral bands and gaps. The results of these calculations suggest that, in the semiclassical limit, all excitations in the spectrum become effective solitons. Finally, the analytical predictions are compared with direct numerical simulations as well as with numerical calculations of the Lax spectrum, and the results are shown to be in excellent agreement.

We derive and analyze, analytically and numerically, two first-order continuum models to approximate the nonlinear dynamics of granular crystal lattices, focusing specifically on solitary waves, periodic waves, and dispersive shock waves. The dispersive shock waves predicted by the two continuum models are studied using modulation theory, DSW fitting techniques, and direct numerical simulations. The PDE-based predictions show good agreement with the DSWs generated by the discrete model simulation of the granular lattice itself, even in cases where no precompression is present and the lattice is purely nonlinear. Such an effective description could prove useful for future, more analytically amenable approximations of the original lattice system.

A (2+1)-dimensional hyperbolic system of four quasi-linear partial differential equations is derived that describes the modulations of lump solutions of the Kadomtsev-Petviashvili I (KPI) equation in the presence of a mean field. The system is then shown to satisfy the necessary conditions for integrability of hydrodynamic chains. Moreover, a suitable reduction of the resulting modulation system is applied to study the interactions between lumps and a rarefaction wave for the mean field. Precise conditions are derived that describe how the lump parameters change as a result of the interaction, and which in particular determine whether the lump is transmitted through or trapped inside the rarefaction wave. The theoretical predictions are compared to direct numerical simulations of the KPI equation, showing excellent agreement.

A Lax pair (L,P) is sometimes thought of as a structural certificate, in that the spatial operator L carries the spectral data of an integrable system, and its isospectral evolution under ∂_t L = [L,P] encodes the nonlinear dynamics. Yet, experience shows that the correspondence between equations and Lax pairs is much more nuanced than this picture suggests. Equations can admit Lax pairs that fail to encode the expected integrable structure. This paper probes that anomalous corner of the Lax pair landscape through five case studies (the Euler top, the free Schrödinger equation, the inviscid Burgers equation, the shallow water system, and the Korteweg–de Vries equation), each illustrating a different way the link to integrability can be distorted. The approach combines analytical calculations with the Sparse Identification of Lax Operators (SILO) framework, which proved useful throughout, in some cases confirming the textbook pair and in others surfacing alternatives worth understanding on their own terms. The recurring lesson across the five cases is that compatibility underdetermines the Lax representation, so that anomalous pairs are regular features of the landscape rather than pathologies. Notably, we show that a spectrally degenerate Korteweg–de Vries Lax pair, classified as fake by standard criteria, still generates the full conservation hierarchy through its operator algebra, which shows that a blunt dichotomy between true and fake Lax pairs can be too reductive.

In the present work we analyze traveling and dispersive shock waves of a two-dimensional Fermi-Pasta-Ulam-Tsingou lattice. In the first part of the paper, using variational techniques we prove the existence of both periodic and solitary traveling waves for convex potentials. In the case of unimodal profiles we are able to remove the assumption of convexity. The variational formulation also provides a natural algorithm for the numerical computation of traveling waves, which we use to explore both solitary and periodic traveling waves. The numerical computations are compared with analytical approximations based on the derivation of the KdV equation for quasi-one-dimensional propagation. In the second part of the paper, we focus on dispersive shock waves (DSWs), which are expanding modulated waves that connect states of different amplitude. In particular, we focus on line DSWs, which are constant along one direction and propagate in the direction orthogonal to which it is constant. Such solutions form when subject to quasi-one-dimensional jump initial data. We find that while the shape of the DSW depends on the direction of travel, properties such as the speed and amplitude do not. The systematic numerical study of the line DSWs is then compared to those predicted by the KdV equation along the line of propagation. Key characteristics of the DSWs, such as the speeds of the trailing and leading edges, are investigated for various jump heights, yielding good agreement between simulation and KdV approximation in the limit of vanishing jump height. Finally, we apply the DSW fitting method to study the trailing and leading edge characteristics of the DSW, finding even better agreement to the numerics when compared to the KdV prediction. The KdV prediction and DSW fitting predictions agree in the limit of small jump height.

We study an initial-boundary-value problem (IBVP) for a system of coupled Maxwell-Bloch equations (CMBE) that model two colors or polarizations of light resonantly interacting with a degenerate, two-level, active optical medium with an excited state and a pair of degenerate ground states. We assume that the electromagnetic field approaches non-vanishing plane waves in the far past and future. This type of interaction has been found to underlie nonlinear optical phenomena including electromagnetically induced transparency, slow light, stopped light, and quantum memory. Under the assumptions of unidirectional, lossless propagation of slowly-modulated plane waves, the resulting CMBE become completely integrable in the sense of possessing a Lax Pair. In this paper, we formulate an inverse scattering transform (IST) corresponding to these CMBE and their Lax pair, allowing for the spectral line of the atomic transitions in the active medium to have a finite width. The scattering problem for this Lax pair is the same as for the Manakov system. The main advancement in this IST for CMBE is calculating the nontrivial spatial propagation of the spectral data and determining the state of the optical medium in the distant future from that in the distant past, which is needed for the complete formulation of the IBVP. The Riemann-Hilbert problem is used to extract the spatio-temporal dependence of the solution from the evolving spectral data. We further derive and analyze several types of solitons and determine their velocity and stability, as well as find dark states of the medium which fail to interact with a given pulse.

We present an analytical model of integrable turbulence in the focusing nonlinear Schrödinger (fNLS) equation, generated by a one-parameter family of finite-band elliptic potentials in the semiclassical limit. We show that the spectrum of these potentials exhibits a thermodynamic band/gap scaling compatible with that of soliton and breather gases depending on the value of the elliptic parameter m of the potential. We then demonstrate that, upon augmenting the potential by a small random noise (which is inevitably present in real physical systems), the solution of the fNLS equation evolves into a fully randomized, spatially homogeneous breather gas, a phenomenon we call breather gas fission. We show that the statistical properties of the breather gas at large times are determined by the spectral density of states generated by the unperturbed initial potential. We analytically compute the kurtosis of the breather gas as a function of the elliptic parameter m, and we show that it is greater than 2 for all nonzero m, implying non-Gaussian statistics. Finally, we verify the theoretical predictions by comparison with direct numerical simulations of the fNLS equation. These results establish a link between semiclassical limits of integrable systems and the statistical characterization of their soliton and breather gases.

Modulational instability (MI) is a fundamental phenomenon in the study of nonlinear dynamics, spanning diverse areas such as shallow water waves, optics, and ultracold atomic gases. In particular, the nonlinear stage of MI has recently been a topic of intense exploration and has been shown to manifest, in many cases, in the generation of dispersive shock waves (DSWs). In this Letter, we experimentally probe the MI dynamics in an immiscible two-component ultracold atomic gas with exclusively repulsive interactions, catalyzed by a hard-wall-like boundary produced by a repulsive optical barrier. We analytically describe the expansion rate of the DSWs in this system, generalized to arbitrary intercomponent interaction strengths and species ratios. We observe excellent agreement among the analytical results, an effective 1D numerical model, full 3D numerical simulations, and experimental data. Additionally, we extend this scenario to the interaction between two counterpropagating DSWs, which leads to the production of Peregrine soliton structures. These results further demonstrate the versatility of atomic platforms toward the controlled realization of DSWs and rogue waves.

The oblique collisions and dynamical interference patterns of two-dimensional dispersive shock waves are studied numerically and analytically via the temporal dynamics induced by wedge-shaped initial conditions for the Kadomtsev-Petviashvili II equation. Various asymptotic wave patterns are identified, classified, and characterized in terms of the incidence angle and the amplitude of the initial step, which can give rise to either subcritical or supercritical configurations, including the generalization to dispersive shock waves of the Mach reflection and expansion of viscous shocks and line solitons. An eightfold amplification of the amplitude of an obliquely incident flow upon a wall at the critical angle is demonstrated. Applications of the results include bore interactions in geophysical fluid dynamics.

We report the first experimental observation of breather gases (BGs) in optics, realized in a recirculating fiber loop enabling virtually lossless propagation over 1200 km. Initiated by a slowly modulated optical background perturbed by noise, the BGs form through a nonlinear fission process and demonstrate spatiotemporal dynamics that align closely with numerical simulations of the focusing one-dimensional nonlinear Schrödinger equation. The minimal dissipation in our setup enables a statistical characterization of the BGs and confirms the theoretically predicted doubling of kurtosis during the evolution of the BGs. These results open new avenues for experimental studies of integrable turbulence involving solitons on finite background.

The inverse spectral theory for a self-adjoint one-dimensional Dirac operator associated periodic potentials is formulated via a Riemann-Hilbert problem approach. The resulting formalism is also used to solve the initial value problem for the nonlinear Schrodinger (NLS) equation. A uniqueness theorem for the solutions of the Riemann-Hilbert problem is established, which provides a new method for obtaining the potential from the spectral data. Two additional, scalar Riemann-Hilbert problems are also formulated that provide conditions for the periodicity in space and time of the solution generated by arbitrary sets of spectral data. The formalism applies for both finite-genus and infinite-genus potentials. The formalism also shows that only a single set of Dirichlet eigenvalues is needed in order to uniquely reconstruct the potential of the Dirac operator and the corresponding solution of the defocusing NLS equation, in contrast with the representation of the solution of the NLS equation via the finite-genus formalism, in which two different sets of Dirichlet eigenvalues are used.

In this paper we focus on a discrete physical model describing granular crystals, whose equations of motion can be described by a system of differential difference equations (DDEs). After revisiting earlier continuum approximations, we propose a regularized continuum model variant to approximate the discrete granular crystal model through a suitable partial differential equation (PDE). We then compute, both analytically and numerically, its traveling wave and periodic traveling wave solutions, in addition to its conservation laws. Next, using the periodic solutions, we describe quantitatively various features of the dispersive shock wave (DSW) by applying Whitham modulation theory and the DSW fitting method. Finally, we perform several sets of systematic numerical simulations to compare the corresponding DSW results with the theoretical predictions and illustrate that the continuum model provides a good approximation of the underlying discrete one.

We study two-dimensional stationary soliton gas in the framework of the time-independent reduction of the Kadomtsev-Petviashvili (KPII) equation, which coincides with the integrable two-way "good" Boussinesq equation in the xy plane. This (2 + 0)D reduction enables the construction of the spatial analog of the kinetic equation for the stationary gas of KP solitons by invoking recent results on (1 + 1)D bidirectional soliton gases and generalized hydrodynamics of the Boussinesq equation. We then use the spectral kinetic theory to analytically describe two basic types of 2D soliton gas interactions: (i) refraction of a line soliton by a stationary soliton gas, and (ii) oblique interference of two soliton gases. We verify the analytical predictions by numerically implementing the corresponding KPII soliton gases via exact N-soliton solutions with N-large and appropriately chosen random distributions for the soliton parameters. We also explicitly evaluate the long-distance correlations for the two-component interference configurations. The results can be applied to a variety of physical systems, from shallow water waves to Bose-Einstein condensates.

The Maxwell-Bloch system of equations with inhomogeneous broadening is studied, and the local and global well-posedness of the corresponding initial-boundary value problem is established by taking advantage of the integrability of the system and making use of the corresponding inverse scattering transform (IST). A key ingredient in the analysis is the L2{L}^{2}-Sobolev bijectivity of the direct and IST established by Xin Zhou for the focusing Zakharov-Shabat problem.

In the present work we revisit the shock wave dynamics in a granular chain with precompression. By approximating the model by an alpha-Fermi-Pasta-Ulam-Tsingou chain, we leverage the connection of the latter in the strain variable formulation to two separate integrable models, one continuum, namely the KdV equation, and one discrete, namely the Toda lattice. We bring to bear the Whitham modulation theory analysis of such integrable systems and the analytical approximation of their dispersive shock waves in order to provide, through the lens of the reductive connection to the granular crystal, an approximation to the shock wave of the granular problem. A detailed numerical comparison of the original granular chain and its approximate integrable-system-based dispersive shocks proves very favorable in a wide parametric range. The gradual deviations between (approximate) theory and numerical computation, as amplitude parameters of the solution increase are quantified and discussed.

The inverse scattering transform is developed to solve the Maxwell-Bloch system of equations that describes two-level systems with inhomogeneous broadening, in the case of optical pulses that do not vanish at infinity in the future. The direct problem, which is formulated in terms of a suitably-defined uniformization variable, combines features of the formalism with decaying as well as non-decaying fields. The inverse problem is formulated in terms of a 2× 2 matrix Riemann-Hilbert problem. A novel aspect of the problem is that no reflectionless solutions can exist, and solitons are always accompanied by radiation. At the same time, it is also shown that, when the medium is initially in the ground state, the radiative components of the solutions decay upon propagation into the medium, giving rise to an asymptotically reflectionless states. Like what happens when the optical pulse decays rapidly in the distant past and the distant future, a medium that is initially excited decays to the stable ground state as t→∞ and for sufficiently large propagation distances. Finally, the asymptotic state of the medium and certain features of the optical pulse inside the medium are considered, and the emergence of a transition region upon propagation in the medium is briefly discussed.

We experimentally realize the Peregrine soliton in a highly particle-imbalanced two-component repulsive Bose-Einstein condensate in the immiscible regime. The effective focusing dynamics and resulting modulational instability of the minority component provide the opportunity to dynamically create a Peregrine soliton with the aid of an attractive potential well that seeds the initial dynamics. The Peregrine soliton formation is highly reproducible, and our experiments allow us to separately monitor the minority and majority components, and to compare with the single component dynamics in the absence or presence of the well with varying depths. We showcase the centrality of each of the ingredients leveraged herein. Numerical corroborations and a theoretical basis for our findings are provided through three-dimensional simulations emulating the experimental setting and via a one-dimensional analysis further exploring its evolution dynamics.

Two-dimensional reductions of the KP-Whitham system, namely the overdetermined Whitham modulation system for five dependent variables that describe the periodic solutions of the Kadomtsev-Petviashvili equation, are studied and characterized. Three different reductions are considered corresponding to modulations that are independent of x , independent of y , and of t (i.e., stationary), respectively. Each of these reductions still describes dynamic, two-dimensional spatial configurations since the modulated cnoidal wave generically has a nonzero speed and a nonzero slope in the x y plane. In all three of these reductions, the properties of the resulting systems of equations are studied. It is shown that the resulting reduced system is not integrable unless one enforces the compatibility of the system with all conservation of waves equations (or considers a reduction to the harmonic or soliton limit). In all cases, compatibility with conservation of waves yields a reduction in the number of dependent variables to two, three and four, respectively. As a byproduct of the stationary case, the Whitham modulation system for the Boussinesq equation is also explicitly obtained.