Search

Nonlocal integrable partial differential equations possessing a spatial or temporal reflection have constituted an active research area for the past decade. Recently, more general classes of these nonlocal equations have been proposed, wherein the nonlocality appears as a combination of a shift (by a real or a complex parameter) and a reflection. This new shifting parameter manifests itself in the inverse scattering transform (IST) as an additional phase factor in an analogous way to the classical Fourier transform. In this paper, the IST is analyzed in detail for several examples of such systems. Particularly, time, space, and space-time shifted nonlinear Schr\"odinger (NLS) and space-time shifted modified Korteweg-de Vries (mKdV) equations are studied. Additionally, the semi-discrete IST is developed for the time, space and space-time shifted variants of the Ablowitz-Ladik integrable discretization of the NLS. One soliton solutions are constructed for all continuous and discrete cases., Comment: 33 pages, revised based on reviewer comments

A honeycomb Floquet lattice with helically rotating waveguides and an interface separating two counter-propagating subdomains is analyzed. Two topologically protected localized waves propagate unidirectionally along the interface. Switching can occur when these interface modes reach the edge of the lattice and the light splits into waves traveling in two opposite directions. The incoming mode, traveling along the interface, can be routed entirely or partially along either lattice edge with the switching direction based on a suitable mixing of the interface modes.

In recent years, there has been considerable interest in the study of wave propagation in nonlinear photonic lattices. The interplay between nonlinearity and periodicity has led researchers to manipulate light and discover new and interesting phenomena such as new classes of localized modes, usually referred to as solitons, and novel surface states that propagate robustly. A field where both nonlinearity and periodicity arises naturally is nonlinear optics. But there are other areas where waves propagating on background lattices play an important role, including photonic crystal fibers and Bose-Einstein condensation. In this review article the propagation of wave envelopes in one and two-dimensional periodic lattices associated with additional potential in the nonlinear Schrodinger (NLS) equation, termed lattice NLS equations, are studied. A discrete reduction, known as the tight-binding approximation, is employed to find the linear dispersion relation and the equations governing nonlinear discrete envelopes for two-dimensional simple periodic lattices and two-dimensional non-simple honeycomb lattices. In the limit under which the envelopes vary slowly, continuous envelope equations are derived from the discrete system. The coefficients of the linear evolution system are related to the dispersion relation in both the discrete and continuous cases. For simple lattices, the continuous systems are NLS type equations. In honeycomb lattices, in certain cases, the continuous system is found to be nonlinear Dirac equations. Finally, it is possible to realize so-called topological insulator systems in an optical waveguide setting. The modes supported by these systems are associated with spectral topological invariants and, remarkably, can propagate without backscatter from lattice defects., Comment: Review Article

The propagation of localized solitons in the presence of large-scale waves is a fundamental problem, both physically and mathematically, with applications in fluid dynamics, nonlinear optics and condensed matter physics. Here, the evolution of a soliton as it interacts with a rarefaction wave or a dispersive shock wave, examples of slowly varying and rapidly oscillating dispersive mean fields, for the Korteweg-de Vries equation is studied. Step boundary conditions give rise to either a rarefaction wave (step up) or a dispersive shock wave (step down). When a soliton interacts with one of these mean fields, it can either transmit through (tunnel) or become embedded (trapped) inside, depending on its initial amplitude and position. A comprehensive review of three separate analytical approaches is undertaken to describe these interactions. First, a basic soliton perturbation theory is introduced that is found to capture the solution dynamics for soliton-rarefaction wave interaction in the small dispersion limit. Next, multiphase Whitham modulation theory and its finite-gap description are used to describe soliton-rarefaction wave and soliton-dispersive shock wave interactions. Lastly, a spectral description and an exact solution of the initial value problem is obtained through the Inverse Scattering Transform. For transmitted solitons, far-field asymptotics reveal the soliton phase shift through either type of wave mentioned above. In the trapped case, there is no proper eigenvalue in the spectral description, implying that the evolution does not involve a proper soliton solution. These approaches are consistent, agree with direct numerical simulation, and accurately describe different aspects of solitary wave-mean field interaction

Integrable fractional equations such as the fractional Korteweg-deVries and nonlinear Schr\"odinger equations are key to the intersection of nonlinear dynamics and fractional calculus. In this manuscript, the first discrete/differential difference equation of this type is found, the fractional integrable discrete nonlinear Schr\"odinger equation. This equation is linearized; special soliton solutions are found whose peak velocities exhibit more complicated behavior than other previously obtained fractional integrable equations. This equation is compared with the closely related fractional averaged discrete nonlinear Schr\"odinger equation which has simpler structure than the integrable case. For positive fractional parameter and small amplitude waves, the soliton solutions of the integrable and averaged equations have similar behavior.

A unified method to analyze the dynamics and topological structure associated with a class of Floquet topological insulators is presented. The method is applied to a system that describes the propagation of electromagnetic waves through the bulk of a two-dimensional lattice that is helically-driven in the direction of propagation. Tight-binding approximations are employed to derive reduced dynamical systems. Further asymptotic approximations, valid in the high-frequency driving regime, yield a time-averaged system which governs the leading order behavior of the wave. From this follows an analytic calculation of the Berry connection, curvature and Chern number by analyzing the local behavior of the eigenfunctions near the critical points of the spectrum. Examples include honeycomb, Lieb and kagome lattices. In the nonlinear regime novel equations governing slowly varying wave envelopes are derived. For the honeycomb lattice, numerical simulations show that for relatively small nonlinear effects a striking spiral patterns occurs; as nonlinearity increases localized structures emerge and for somewhat higher nonlinearity the waves collapse., Comment: 30 pages, 11 figures

Nonlinear integrable equations serve as a foundation for nonlinear dynamics, and fractional equations are well known in anomalous diffusion. We connect these two fields by presenting the discovery of a new class of integrable fractional nonlinear evolution equations describing dispersive transport in fractional media. These equations can be constructed from nonlinear integrable equations using a widely generalizable mathematical process utilizing completeness relations, dispersion relations, and inverse scattering transform techniques. As examples, this general method is used to characterize fractional extensions to two physically relevant, pervasive integrable nonlinear equations: the Korteweg-de Vries and nonlinear Schr\"odinger equations. These equations are shown to predict super-dispersive transport of non-dissipative solitons in fractional media.

The inverse scattering transform allows explicit construction of solutions to many physically significant nonlinear wave equations. Notably, this method can be extended to fractional nonlinear evolution equations characterized by anomalous dispersion using completeness of suitable eigenfunctions of the associated linear scattering problem. In anomalous diffusion, the mean squared displacement is proportional to $t^{\alpha}$, $\alpha>0$, while in anomalous dispersion, the speed of localized waves is proportional to $A^{\alpha}$, where $A$ is the amplitude of the wave. Fractional extensions of the modified Korteweg-deVries (mKdV), sine-Gordon (sineG) and sinh-Gordon (sinhG) and associated hierarchies are obtained. Using symmetries present in the linear scattering problem, these equations can be connected with a scalar family of nonlinear evolution equations of which fractional mKdV (fmKdV), fractional sineG (fsineG), and fractional sinhG (fsinhG) are special cases. Completeness of solutions to the scalar problem is obtained and, from this, the nonlinear evolution equation is characterized in terms of a spectral expansion. In particular, fmKdV, fsineG, and fsinhG are explicitly written. One-soliton solutions are derived for fmKdV and fsineG using the inverse scattering transform and these solitons are shown to exhibit anomalous dispersion.

We derive the Whitham modulation equations for the nonlinear Schr\"odinger equation in the plane (2d NLS) with small dispersion. The modulation equations are derived in terms of both physical and Riemann variables; the latter yields equations of hydrodynamic type. The complete 2d NLS Whitham system consists of six dynamical equations in evolutionary form and two constraints. As an application, we determine the linear stability of one-dimensional traveling waves. In both the elliptic and hyperbolic case, the traveling waves are found to be unstable. This result is consistent with all previous investigations of such stability by other methods. These results are supported by direct numerical calculations., Comment: 34 pages, 6 figures

Rogue waves are abnormally large waves which appear unexpectedly and have attracted considerable attention, particularly in recent years. The one space, one time (1+1) nonlinear Schr\"odinger equation is often used to model rogue waves; it is an envelope description of plane waves and admits the so-called Pergerine and Kuznetov-Ma soliton solutions. However, in deep water waves and certain electromagnetic systems where there are two significant transverse dimensions, the 2+1 hyperbolic nonlinear Schrodinger equation is the appropriate wave envelope description. Here we show that these rogue wave solutions suffer from strong transverse instability at long and short frequencies. Moreover, the stability of the Peregrine soliton is found to coincide with that of the background plane wave. These results indicate that, when applicable, transverse dimensions must be taken into account when investigating rogue wave pheneomena., Comment: Supplementary Material available on Publisher's website

The Peierls-Nabarro barrier is a discrete effect that frequently occurs in discrete nonlinear systems. A signature of the barrier is the slowing and eventual stopping of discrete solitary waves. This work examines intense electromagnetic waves propagating through a periodic honeycomb lattice of helically-driven waveguides, which serves as a paradigmatic Floquet topological insulator. Here it is shown that discrete topologically protected edge modes do not suffer from the typical slowdown associated with the Peierls-Nabarro barrier. Instead, as a result of their topological nature, the modes always move forward and redistribute their energy: a narrow (discrete) mode transforms into a wide effectively continuous mode. On the other hand, a discrete edge mode that is not topologically protected does eventually slow down and stop propagating. Topological modes that are initially narrow modes naturally tend to wide envelope states that are described by a generalized nonlinear Schrodinger equation. These results provide insight into the nature of nonlinear topological insulators and their application.

The semi-classical Korteweg-deVries equation for step-like data is considered with a small parameter in front of the highest derivative. Using perturbation analysis Whitham theory is constructed to higher order. This allows the order one phase and the complete leading order solution to be obtained; the results are confirmed by extensive numerical calculations., Comment: 25 pages, 5 figures, a shorter version is accepted for publication in Proc. Roy. Soc. A

Topologically protected waves in the linearly polarized Maxwell's equations with gyrotropic, magneto-optic media were studied a decade ago both computationally and experimentally. This paper develops a robust tight-binding model for this system that makes careful use of Wannier function representations. The model provides very good approximations to the underlying band structure. When solved on a semi-infinite strip, it produces exponentially localized edge modes whose corresponding eigenvalues span the frequency band gaps. A set of coupled differential equations are derived which allows one to find how the electromagnetic field propagates unidirectionally, without backscatter from defects. Furthermore, the discrete model predicts topologically protected edge modes with nontrivial Chern number which are consistent with direct simulation.

Quasi-monochromatic complex reductions of a number of physically important equations are obtained. Starting from the cubic nonlinear Klein-Gordon (NLKG), the Korteweg-deVries (KdV) and water wave equations, it is shown that the leading order asymptotic approximation can be transformed to the well-known integrable AKNS system [6] associated with second order (in space) nonlinear wave equations. This in turn establishes, for the first time, an important physical connection between the recently discovered nonlocal integrable reductions of the AKNS system and physically interesting equations. Reductions include the parity-time, reverse space-time and reverse time nonlocal nonlinear Schrodinger equations., Comment: To appear in J. of Physics A: Mathematical and Theoretical

Dispersive shock waves (DSWs) of the defocusing radial nonlinear Schr\"odinger (rNLS) equation in two spatial dimensions are studied. This equation arises naturally in Bose-Einstein condensates, water waves and nonlinear optics. A unified nonlinear WKB approach, equally applicable to integrable or nonintegrable partial differential equations, is used to find the rNLS Whitham modulation equation system in both physical and hydrodynamic type variables. The description of DSWs obtained via Whitham theory is compared with direct rNLS numerics; the results demonstrate very good quantitative agreement. On the other hand, as expected, comparison with the corresponding DSW solutions of the one-dimensional NLS equation exhibits significant qualitative and quantitative differences., Comment: 58 pages, 8 figures

Topological insulators are studied via tight-binding approximations of longitudinally driven photonic lattices with three lattice sites per unit cell. Two cases are considered in detail: Lieb and Kagome lattices. The lattice is decomposed into three sublattices each of which are allowed move independently of one another. Emphasis is placed on periodic driving induced by laser-etched helical coils along the direction of propagation. The linear Floquet bands are constructed for various inter-sublattice rotation patterns such as: different radii, different frequency, phase offset and quasi one-dimensional motion. Depending on the nature of the band structure, bulk spectral bands with nonzero Chern number are found to support topologically protected edge states which can move unidirectionally. In this case, the modes move scatter-free around defects due to underlying topological protection. Intriguing mode dynamics are found including bi-directional topological modes and bulk-edge leakage i.e. excitation of bulk modes at a defect for edge modes with dispersion frequencies nearby the bulk bands. Finally, certain nonlinear edge modes are also found to propagate unidirectionally and scatter-free around lattice defects., Comment: 17 Figures

The Korteweg-deVries (KdV) equation with step boundary conditions is considered, with an emphasis on soliton dynamics. When one or more initial solitons are of sufficient size they can propagate through the step; in this case the phase shift is calculated via the inverse scattering transform. On the other hand, when the amplitude is too small they become trapped. In the trapped case the transmission coefficient of the associated associated linear Schr\"odinger equation can become large at a point exponentially close to the continuous spectrum. This point is referred to as a {\it pseudo-embedded eigenvalue}. Employing the inverse problem it is shown that the continuous spectrum associated with a branch cut in the neighborhood of the pseudo-embedded eigenvalue plays the role of discrete spectra, which in turn leads to a trapped soliton in the KdV equation., Comment: 8 pages, 2 figures

General soliton solutions to a nonlocal nonlinear Schr\"odinger (NLS) equation with PT-symmetry for both zero and nonzero boundary conditions {are considered} via the combination of Hirota's bilinear method and the Kadomtsev-Petviashvili (KP) hierarchy reduction method. First, general $N$-soliton solutions with zero boundary conditions are constructed. Starting from the tau functions of the two-component KP hierarchy, it is shown that they can be expressed in terms of either Gramian or double Wronskian determinants. On the contrary, from the tau functions of single component KP hierarchy, general soliton solutions to the nonlocal NLS equation with nonzero boundary conditions are obtained. All possible soliton solutions to nonlocal NLS with Parity (PT)-symmetry for both zero and nonzero boundary conditions are found in the present paper., Comment: 26 pages, 3 figures

Whitham modulation theory for certain two-dimensional evolution equations of Kadomtsev-Petviashvili (KP) type is presented. Three specific examples are considered in detail: the KP equation, the two-dimensional Benjamin-Ono (2DBO) equation and a modified KP (m2KP) equation. A unified derivation is also provided. In the case of the m2KP equation, the corresponding Whitham modulation system exhibits features different from the other two. The approach presented here does not require integrability of the original evolution equation. Indeed, while the KP equation is known to be a completely integrable equation, the 2DBO equation and the m2KP equation are not known to be integrable. In each of the cases considered, the Whitham modulation system obtained consists of five first-order quasilinear partial differential equations. The Riemann problem (i.e. the analogue of the Gurevich-Pitaevskii problem) for the one-dimensional reduction of the m2KP equation is studied. For the m2KP equation, the system of modulation equations is used to analyze the linear stability of traveling wave solutions., Comment: 27 pages, 2 figures

Whitham modulation theory for the two dimensional Benjamin-Ono (2DBO) equation is presented. A system of five quasi-linear first-order partial differential equations is derived. The system describes modulations of the traveling wave solutions of the 2DBO equation. These equations are transformed to a singularity-free hyrdodynamic-like system referred to here as the 2DBO-Whitham system. Exact reductions of this system are discussed, the formulation of initial value problems is considered, and the system is used to study the transverse stability of traveling wave solutions of the 2DBO equation., Comment: 8 pages, 1 figure, one reference and final journal information added

A systematic approach for deriving tight-binding approximations in general longitudinally driven lattices is presented. As prototypes, honeycomb and staggered square lattices are considered. Time-reversal symmetry is broken by varying/rotating the waveguides, longitudinally, along the direction of propagation. Different sublattice rotation and structure are allowed. Linear Floquet bands are constructed for intricate sublattice rotation patterns such as counter rotation, phase offset rotation, as well as different lattice sizes and frequencies. An asymptotic analysis of the edge modes, valid in a rapid-spiraling regime, reveals linear and nonlinear envelopes which are governed by linear and nonlinear Schrodinger equations, respectively. Nonlinear states, referred to as topologically protected edge solitons are unidirectional edge modes. Direct numerical simulations for both the linear and nonlinear edge states agree with the asymptotic theory. Topologically protected modes are found; they possess unidirectionality and do not scatter at lattice defect boundaries., Comment: 17 pages, 15 figures

The reverse space-time (RST) Sine-Gordon, Sinh-Gordon and nonlinear Schr\"odinger equations were recently introduced and shown to be integrable infinite-dimensional dynamical systems. The inverse scattering transform (IST) for rapidly decaying data was also constructed. In this paper, IST for these equations with nonzero boundary conditions (NZBCs) at infinity is presented. The NZBC problem is more complicated due to the associated branching structure of the associated linear eigenfunctions. With constant amplitude at infinity, four cases are analyzed; they correspond to two different signs of nonlinearity and two different values of the phase at infinity. Special soliton solutions are discussed and explicit 1-soliton and 2-soliton solutions are found. In terms of IST, the difference between the RST Sine-Gordon/Sinh-Gordon equations and the RST NLS equation is the time dependence of the scattering data. Spatially dependent boundary conditions are also briefly considered., Comment: arXiv admin note: substantial text overlap with arXiv:1612.02726

In 2013 a new nonlocal symmetry reduction of the well-known AKNS scattering problem was found; it was shown to give rise to a new nonlocal $PT$ symmetric and integrable Hamiltonian nonlinear Schr\"{o}dinger (NLS) equation. Subsequently, the inverse scattering transform was constructed for the case of rapidly decaying initial data and a family of spatially localized, time periodic one soliton solution were found. In this paper, the inverse scattering transform for the nonlocal NLS equation with nonzero boundary conditions at infinity is presented in the four cases when the data at infinity have constant amplitudes. The direct and inverse scattering problems are analyzed. Specifically, the direct problem is formulated, the analytic properties of the eigenfunctions and scattering data and their symmetries are obtained. The inverse scattering problem is developed via a left-right Riemann-Hilbert problem in terms of a suitable uniformization variable and the time dependence of the scattering data is obtained. This leads to a method to linearize/solve the Cauchy problem. Pure soliton solutions are discussed and explicit 1-soliton solution and two 2-soliton solutions are provided for three of the four different cases corresponding to two different signs of nonlinearity and two different values of the phase difference between plus and minus infinity. In the one other case there are no solitons.

The genus-1 KP-Whitham system is derived for both variants of the Kadomtsev-Petviashvili (KP) equation (namely, the KPI and KPII equations). The basic properties of the KP-Whitham system, including symmetries, exact reductions, and its possible complete integrability, together with the appropriate generalization of the one-dimensional Riemann problem for the Korteweg-deVries equation are discussed. Finally, the KP-Whitham system is used to study the linear stability properties of the genus-1 solutions of the KPI and KPII equations; it is shown that all genus-1 solutions of KPI are linearly unstable while all genus-1 solutions of KPII {are linearly stable within the context of Whitham theory., Comment: Significantly revised version

A nonlocal nonlinear Schr\"odinger (NLS) equation was recently found by the authors and shown to be an integrable infinite dimensional Hamiltonian equation. Unlike the classical (local) case, here the nonlinearly induced "potential" is $PT$ symmetric thus the nonlocal NLS equation is also $PT$ symmetric. In this paper, new {\it reverse space-time} and {\it reverse time} nonlocal nonlinear integrable equations are introduced. They arise from remarkably simple symmetry reductions of general AKNS scattering problems where the nonlocality appears in both space and time or time alone. They are integrable infinite dimensional Hamiltonian dynamical systems. These include the reverse space-time, and in some cases reverse time, nonlocal nonlinear Schr\"odinger, modified Korteweg-deVries (mKdV), sine-Gordon, $(1+1)$ and $(2+1)$ dimensional three-wave interaction, derivative NLS, "loop soliton", Davey-Stewartson (DS), partially $PT$ symmetric DS and partially reverse space-time DS equations. Linear Lax pairs, an infinite number of conservation laws, inverse scattering transforms are discussed and one soliton solutions are found. Integrable reverse space-time and reverse time nonlocal discrete nonlinear Schr\"odinger type equations are also introduced along with few conserved quantities. Finally, nonlocal Painlev\'e type equations are derived from the reverse space-time and reverse time nonlocal NLS equations., Comment: This paper will be published in Studies in Applied Mathematics in a special volume dedicated to Professor David J. Benny

The generation of rogue waves is investigated via a nonlocal nonlinear Schrodinger (NLS) equation. In this system, modulation instability is suppressed and is usually expected that rogue wave formation would also be limited. On the contrary, a parameter regime is identified where the instability is suppressed but nevertheless the number and amplitude of the rogue events increase, as compared to the standard NLS (which is a limit of the nonlocal system). Furthermore, the nature of these waves is investigated; while no analytical solutions are known to model these events, numerically it is shown that they differ significantly from either the rational (Peregrine) or soliton solution of the limiting NLS equation. As such, these findings may also help in rogue wave realization experimentally in these media.

Rogue waves in birefringent optical fibers are analyzed within the framework of the coupled nonlinear Schr\"odinger (CNLS) system. The generation of rogue waves is frequently associated with modulation instability (MI). It is commonly expected that since the CNLS equations exhibit larger growth rates they should also produce more rogue events than their scalar counterparts. This is found to occur only when both equations are focusing. When at least one component is defocusing, the CNLS equations may still exhibit larger growth rates, compared to the scalar system, but that does not necessarily result in more or larger events. The birefringence angle for which the maximum number of events occurs is also identified and the nature of the rogue wave is described for the different cases.

The appearance of a fundamental long-time asymptotic regime in the two space one time dimensional hyperbolic nonlinear Schr\"odinger (HNLS) equation is discussed. Based on analytical and extensive numerical simulations an approximate self-similar solution is found for a wide range of initial conditions -- essentially for initial lumps of small to moderate energy. Even relatively large initial amplitudes, which imply strong nonlinear effects, eventually lead to local structures resembling those of the self-similar solution, with appropriate small modifications. These modifications are important in order to properly capture the behavior of the phase of the solution. This solution has aspects that suggest it is a universal attractor emanating from wide ranges of initial data., Comment: 36 pages, 26 pages text + 20 figures

Dispersive shock waves (DSWs) in the Kadomtsev-Petviashvili (KP) equation and two dimensional Benjamin-Ono (2DBO) equation are considered using parabolic front initial data. Employing a front tracking type ansatz exactly reduces the study of DSWs in two space one time (2+1) dimensions to finding DSW solutions of (1+1) dimensional equations. With this ansatz, the KP and 2DBO equations can be exactly reduced to cylindrical Korteweg-de Vries (cKdV) and cylindrical Benjamin-Ono (cBO) equations, respectively. Whitham modulation equations which describe DSW evolution in the cKdV and cBO equations are derived in general and Riemann type variables are introduced. DSWs obtained from the numerical solutions of the corresponding Whitham systems and direct numerical simulations of the cKdV and cBO equations are compared with excellent agreement obtained. In turn, DSWs obtained from direct numerical simulations of the KP and 2DBO equations are compared with the cKdV and cBO equations, again with remarkable agreement. It is concluded that the (2+1) DSW behavior along parabolic fronts can be effectively described by the DSW solutions of the reduced (1+1) dimensional equations., Comment: 25 Pages, 16 Figures. The movies showing dispersive shock wave propagation in Kadomtsev-Petviashvili II and Two Dimensional Benjamin-Ono equations are available at https://youtu.be/AExAQHRS_vE and https://youtu.be/aXUNYKFlke0

The propagation of linear and nonlinear edge modes in bounded photonic honeycomb lattices formed by an array of rapidly varying helical waveguides is studied. These edge modes are found to exhibit strong transmission (reflection) around sharp corners when the dispersion relation is topologically nontrivial (trivial), and can also remain stationary. An asymptotic theory is developed that establishes the presence (absence) of edge states on all four sides, including in particular armchair edge states, in the topologically nontrivial (trivial) case. In the presence of topological protection, nonlinear edge solitons can persist over very long distances., Comment: 5 pages, 4 figures. Minor updates on the presentation and interpretation of results. The movies showing transmission and reflection of linear edge modes are available at https://www.youtube.com/watch?v=XhaZZlkMadQ and https://www.youtube.com/watch?v=R8NOw0NvRu0

Pseudospin, an additional degree of freedom inherent in graphene, plays a key role in understanding many fundamental phenomena such as the anomalous quantum Hall effect, electron chirality and Klein paradox. Unlike the electron spin, the pseudospin was traditionally considered as an unmeasurable quantity, immune to Stern-Gerlach-type experiments. Recently, however, it has been suggested that graphene pseudospin is a real angular momentum that might manifest itself as an observable quantity, but so far direct tests of such a momentum remained unfruitful. Here, by selective excitation of two sublattices of an artificial photonic graphene, we demonstrate pseudospin-mediated vortex generation and topological charge flipping in otherwise uniform optical beams with Bloch momentum traversing through the Dirac points. Corroborated by numerical solutions of the linear massless Dirac-Weyl equation, we show that pseudospin can turn into orbital angular momentum completely, thus upholding the belief that pseudospin is not merely for theoretical elegance but rather physically measurable., Comment: 12 pages, 5 figures

The response of a passive mode-locking mechanism, where gain and spectral filtering are saturated with the energy and loss saturated with the power, is examined under the presence of higher order effects. These include third order dispersion, self-steepening and Raman gain. The locking mechanism is maintained even with these terms; mode-locking occurs for both the anomalous and normal regimes. In the anomalous regime, these perturbations are found to affect the speed but not the structure of the (locked) pulses. In fact, these pulses behave like solitons of a classical nonlinear Schrodinger equation and as such a soliton perturbation theory is used to verify the numerical observations. In the normal regime, the effect of the perturbations is small, in line with recent experimental observations. The results in the normal regime are verified mathematically using a WKB type asymptotic theory. Finally, bi-solitons are found to behave as dark solitons on top of a stable background and are significantly affected by these perturbations.

A rogue wave formation mechanism is proposed within the framework of a coupled nonlinear Schrodinger (CNLS) system corresponding to the interaction of two waves propagating in oblique directions in deep water. A rogue condition is introduced that links the angle of interaction with the group velocities of these waves: different angles of interaction can result in a major enhancement of rogue events in both numbers and amplitude. For a range of interacting directions it is found that the CNLS system exhibits significantly more extreme wave amplitude events than its scalar counterpart. Furthermore, the rogue events of the coupled system are found to be well approximated by hyperbolic secant functions; they are vectorial soliton-type solutions of the CNLS system, typically not considered to be integrable. Overall, our results indicate that crossing states provide an important mechanism for the generation of rogue water wave events.

The long-time asymptotic solution of the Korteweg-de Vries equation for general, step-like initial data is analyzed. Each sub-step in well-separated, multi-step data forms its own single dispersive shock wave (DSW); at intermediate times these DSWs interact and develop multiphase dynamics. Using the inverse scattering transform and matched-asymptotic analysis it is shown that the DSWs merge to form a single-phase DSW, which is the `largest' one possible for the boundary data. This is similar to interacting viscous shock waves (VSW) that are modeled with Burgers' equation, where only the single, largest-possible VSW remains after a long time., Comment: 5 pages, 3 figures

A long wave multi-dimensional approximation of shallow water waves is the bi-directional Benney-Luke equation. It yields the well-known Kadomtsev-Petviashvili equation in a quasi one-directional limit. A direct perturbation method is developed; it uses the underlying conservation laws to determine the slow evolution of parameters of two space dimensional, non-decaying web-type solutions to the Benney-Luke equation. New numerical simulations, based on windowing methods which are effective for non-decaying data, are presented. These simulations support the analytical results and elucidate the relationship between the Kadomtsev-Petviashvilli and the Benney-Luke equations and are also used to obtain amplitude information regarding particular web solutions. Additional dissipative perturbations to the Benney-Luke equation are also studied.

Ocean waves are complex and often turbulent. While most ocean wave interactions are essentially linear, sometimes two or more waves interact in a nonlinear way. For example, two or more waves can interact and yield waves that are much taller than the sum of the original wave heights. Most of these nonlinear interactions look like an X or a Y or two connected Ys; at other times, several lines appear on each side of the interaction region. It was thought that such nonlinear interactions are rare events: they are not. Here we report that such nonlinear interactions occur every day, close to low tide, on two flat beaches that are about 2,000 km apart. These interactions are closely related to the analytic, soliton solutions of a widely studied multi-dimensional nonlinear wave equation. On a much larger scale, tsunami waves can merge in similar ways., Comment: 5 pages, 6 figures

Dark soliton formation in mode-locked lasers is investigated by means of a power-energy saturation model which incorporates gain and filtering saturated with energy, and loss saturated with power. It is found that general initial conditions evolve into dark solitons under appropriate requirements also met in the experimental observations. The resulting pulses are well approximated by dark solitons of the unperturbed nonlinear Schr\"{o}dinger equation. Notably, the same framework also describes bright pulses in anomalous and normally dispersive lasers.

A method for approximating dark soliton solutions of the nonlinear Schrodinger equation under the influence of perturbations is presented. The problem is broken into an inner region, where core of the soliton resides, and an outer region, which evolves independently of the soliton. It is shown that a shelf develops around the soliton which propagates with speed determined by the background intensity. Integral relations obtained from the conservation laws of the nonlinear Schrodinger equation are used to approximate the shape of the shelf. The analysis is developed for both constant and slowly evolving backgrounds. A number of problems are investigated including linear and nonlinear damping type perturbations.

Soliton strings in mode-locked lasers are obtained using a variant of the nonlinear Schrodinger equation, appropriately modified to model power (intensity) and energy saturation. This equation goes beyond the well-known master equation often used to model these systems. It admits mode-locking and soliton strings in both the constant dispersion and dispersion-managed systems in the (net) anomalous and normal regimes; the master equation is contained as a limiting case. Analysis of soliton interactions show that soliton strings can form when pulses are a certain distance apart relative to their width. Anti-symmtetric bi-soltion states are also obtained. Initial states mode-lock to these states under evolution. In the anomalous regime individual soliton pulses are well approximated by the solutions of the unperturbed nonlinear Schrodinger equation, while in the normal regime the pulses are much wider and strongly chirped., Comment: 10 pages

The Chazy equation $y''' = 2yy'' - 3y'^2$ is derived from the automorphic properties of Schwarz triangle functions $S(\alpha, \beta, \gamma; z)$. It is shown that solutions $y$ which are analytic in the fundamental domain of these triangle functions, only correspond to certain values of $\alpha, \beta, \gamma$. The solutions are then systematically constructed. These analytic solutions provide all known and one new parametrization of the Eisenstein series $P, Q, R$ introduced by Ramanujan in his modular theories of signature 2, 3, 4 and 6., Comment: to appear in Studies in Applied Mathematics