Open access peer-reviewed chapter

Water Waves and Light: Two Unlikely Partners

Written By

Georgios N. Koutsokostas, Theodoros P. Horikis, Dimitrios J. Frantzeskakis, Nalan Antar and İlkay Bakırtaş

Submitted: 19 March 2020 Reviewed: 08 December 2020 Published: 07 January 2021

DOI: 10.5772/intechopen.95431

From the Edited Volume

Nonlinear Optics - From Solitons to Similaritons

Edited by İlkay Bakırtaş and Nalan Antar

Chapter metrics overview

550 Chapter Downloads

View Full Metrics


We study a generic model governing optical beam propagation in media featuring a nonlocal nonlinear response, namely a two-dimensional defocusing nonlocal nonlinear Schrödinger (NLS) model. Using a framework of multiscale expansions, the NLS model is reduced first to a bidirectional model, namely a Boussinesq or a Benney-Luke-type equation, and then to the unidirectional Kadomtsev-Petviashvili (KP) equation – both in Cartesian and cylindrical geometry. All the above models arise in the description of shallow water waves, and their solutions are used for the construction of relevant soliton solutions of the nonlocal NLS. Thus, the connection between water wave and nonlinear optics models suggests that patterns of water may indeed exist in light. We show that the NLS model supports intricate patterns that emerge from interactions between soliton stripes, as well as lump and ring solitons, similarly to the situation occurring in shallow water.


  • Kadomtsev-Petviashvili equation
  • nonlocal nonlinear Schrödinger equation
  • line solitons
  • soliton interactions
  • lump solitons
  • patterns of light

1. Introduction

For many decades, solitons, namely robust localized waveforms that propagate undistorted in nonlinear dispersive media, have been a topic of particular interest in physics [1, 2] and applied mathematics [3, 4]. These waveforms have a variety of fascinating properties that should be mentioned: due to their particle-like nature, solitons collide elastically, thus preserving their shape after the collision process. In some cases, as e.g. in shallow water, line soliton collisions give rise to the emergence of various wave patterns, including X-, H-, or Y-shaped waves, as well as other, even more complicated, waveforms – see, e.g., Ref. [5] for a set of remarkable examples and observations. Such patterns can be described by multidimensional line soliton solutions of the Kadomtsev-Petviashvili II (KPII) equation, a variant of the KP equation [6] with negative dispersion – as is the case of water waves with small surface tension. In particular, the KP equation can be expressed as:


where σ=±1: for σ=+1 (negative dispersion/small surface tension) Eq. (1) is a KPII, while for σ=1 (positive dispersion/large surface tension) Eq. (1) is a KPI. We should mention that both versions of the KP belong to the – rather limited – class of completely integrable equations in 2+1-dimensions [3].

Important to the collision-induced emergence of different types of wave patterns in shallow water is the fact that quasi (line) solitons of the KPII are stable. Furthermore, of paramount importance is the effect of soliton resonances [7, 8, 9, 10], whereby two (or more) colliding solitons resonate under certain conditions, thus creating novel stable structures. It is, thus, not surprising that the collision dynamics of many robust solitons, can give rise to a wealth of complex wave patterns [11, 12, 13, 14, 15, 16] (see also [1]).

The above discussion motivates the following question: can solitons, similar to those governed by the KPII equation and feature rich collision dynamics (as in shallow water waves), be also predicted – and, possibly, observed – in nonlinear optics? In this context, key model is the nonlinear Schrödinger (NLS) equation, which is known to support soliton solutions (see, e.g., Ref. [17] and references therein). As was shown some time ago [18], the 2+1-dimensional defocusing NLS can be asymptotically reduced to the KPI equation, which allows one to approximate dark soliton stripes of the NLS as line solitons of KPI; then, the fact that line KPI solitons are unstable [3], can be used to investigate the transverse instability of dark solitons [18, 19, 20]. On the other hand, as was recently shown [21], the presence of nonlocality introduces an effective surface tension in the system. Then, in the case where the surface tension is small – corresponding to the case of strong nonlocality – it can be shown that the nonlocal NLS model can be asymptotically reduced to a KPII equation. This result allows for the derivation of approximate stable line soliton solutions of the nonlocal NLS, which behave similarly to the stable line solitons of KPII, i.e., they undergo elastic collisions and form patterns similar to those observed in shallow water.

The scope of this work is to review the above ideas and present a variety of soliton solutions and their dynamics in the framework of a defocusing nonlocal NLS model. These solutions stem from the soliton solutions of KP models, which are derived from the nonlocal NLS via a two-step multiscale expansion method. We show that, similarly to the water wave case, the nonlocal NLS can also support complicated wave patterns, arising from “complex interactions” between two or more antidark soliton stripes that are supported in the strong nonlocality regime. In addition, we present other soliton solutions, namely dark lumps (pertinent to the KPI equation), as well as ring dark and antidark solitons (pertinent to the cylindrical KP (cKP), alias Johnson’s Equation [22, 23]). In particular, a summary of our main results, as well as the organization of this work, are as follows.

First, in Section 2, we introduce the nonlocal NLS model, present its continuous-wave (cw) solution, and study its stability. Then, we perform a two-step multiscale analysis. First, at an intermediate stage, we derive a Boussinesq – or a Benney-Luke (BL) [24] – type equation, which is a long-wave, multi-dimensional and bi-directional, shallow water wave model. Then, from the Boussinesq/BL-type equation, we derive its far field, pertinent to the long time behavior, which is the KP equation, both in Cartesian and in cylindrical geometry. The reduction to the KP model allows us to construct approximate soliton solutions of the nonlocal NLS, which are presented in Section 3. In the strong nonlocality regime, we derive antidark soliton stripes, satisfying an effective KPII equation, which are shown to form patterns similar to those observed in shallow water; these include X-, Y, and H-shaped waveforms, as well as other complicated patterns arising from resonant interactions between more than two solitons. We also present dark lump solitons, as well as ring dark and antidark soliton solutions of the NLS, which are numerically found to propagate undistorted in the framework of the NLS model. Finally, in Section 4, we briefly discuss the conclusions of this work.


2. Nonlocal NLS and multiscale analysis

2.1 Introduction of the model and linear regime

The model under consideration is a two-dimensional (2D) defocusing NLS, with a nonlocal nonlinearity, which is expressed in the following dimensionless form:


where subscripts denote partial derivatives, uC, nR,


is the Laplacian in Cartesian or polar coordinates respectively, while the parameter dR+ measures the degree of nonlocality (see below). The above system finds a number of physical applications. For instance, in the context of optics, this model describes the evolution of the complex electric field envelope u, via a nonlinear wave equation in the paraxial regime, coupled with a diffusion-type equation for the nonlinear correction to the refractive index n depending on the intensity I=u2 (in this case, t represents the propagation coordinate) [25, 26]. In this context, the system (2)–(3) has been used to model experiments on liquid solutions exhibiting thermal nonlinearities [27, 28]. In addition, it has also been used in studies of plasmas (in this case, n represents the relative electron temperature perturbation) [29, 30], as well as in nematic liquid crystals (in this case, n is the optically induced angle perturbation) [31, 32].

As mentioned above, parameter d represents the degree of nonlocality: for d=0, one recovers the local limit, whereby the system (2)–(3) reduces to the defocusing 2D NLS equation with a local cubic (Kerr-type) nonlinearity:


while for d0, Eqs. (2) and (3) feature a spatially nonlocal nonlinearity. Notice that in our perturbative analytical approach below, d will be treated as a free parameter, of order O1, thus allowing the model to acquire an arbitrary degree of nonlocality.

Our analytical approach relies on the analysis of the hydrodynamic form of the model. This can be derived upon using the Madelung transformation:


where u0 is a complex constant, while the unknown real functions ρ and ϕ denote the amplitude and phase of the field u, respectively. Substituting into Eqs. (2) and (3), we obtain the following hydrodynamical system of equations:


where xy or r1rθ is the gradient operator in Cartesian or in polar coordinates respectively.

The simplest, non trivial, solution of Eqs. (2) and (3) is of the form:


which corresponds to the continuous wave (cw) solution u=u0exp(iu02t) and the constant function n=u02. This solution will serve as a “pedestal”, on top of which we will seek nonlinear excitations, namely solitons. It is thus necessary to investigate whether the solution (7) is modulationally stable. This can be done upon introducing to Eqs. (7)(9) the following perturbation ansatz:


where 0<ε1 is a formal small parameter. This way, we obtain, at Oε, the following linear system:


The above system can straightforwardly be decoupled as follows: solving Eq. (12) for n˜, as well as Eq. (13) for Δϕ˜, and substituting back in Eq. (13), the following linear equation for ρ˜ is obtained:


Since the above solution is linear we can examine its stability in first order by examining the dispersion relation of plane wave solutions, i.e., expikrωt, of wavenumber k=k and frequency ω (here, e.g., in Cartesian geometry, r=xy and k=kxky). Then, it can readily be found that these plane waves are characterized by the following dispersion relation:


From the above equation we can obtain the following information. First, it is observed that ωRkR, which indicates that the steady-state solution is modulationally stable. This result is important since, below, we will seek for soliton solutions on top of the stable cw background (7). Second, in the long-wavelength limit, such that dk21, Eq. (16) can be reduced to the following Bogoliubov-type dispersion relation:




is the squared “speed of sound”, namely the velocity of linear waves propagating on top of the cw solution. In addition, parameter α is given by:


and plays the role of an effective surface tension for our original system. Indeed, as discussed in Ref. [21], this can be inferred by the fact that Eq. (17) is reminiscent of the dispersion relation for shallow water waves [4]:


Here, c02=gh is the velocity (g is the acceleration of gravity and h the depth of water at rest), while T̂=T/ρgh2, where ρ is the density and T the surface tension. Comparing Eqs. (17) and (20), one can identify the correspondence: 3T̂4du02, which indeed suggests that a surface tension analogue in our problem is du02, depending on the nonlocality parameter d and the cw intensity u02.

Obviously, regarding the magnitude of the effective surface tension, one may identify the following two regimes:

  • α<04du02>1: weak surface tension – strong nonlocality (for fixed u02),

  • α>04du02<1: strong surface tension – weak nonlocality (for fixed u02).

Below it will be shown that these two regimes correspond, respectively, to a KP-II and a KP-I equation, with the former giving rise to complex interactions of line solitons (that we will study in detail).

2.2 The nonlinear regime – asymptotic analysis

2.2.1 The intermediate stage – Bousinesq/Benney-Luke equation

We now proceed by analyzing the fully nonlinear regime. We start by seeking solutions of Eqs. (7)(9) in the form of the following asymptotic expansions:


where 0<ε1 is a formal small parameter and the unknown functions ρj, Φ and nj (j=1,2,) depend on the following stretched variables:


for the Cartesian and polar geometry, respectively (the angular coordinate θ in the polar geometry remains unchanged). By substituting the expansions (21)–(23) into Eqs. (7)(9), and equating terms of the same order in ε, we obtain the following results. First, Eq. (13), at Oε3/2 and Oε5/2, yields:


Furthermore, Eq. (14) at Oε and Oε2 leads to the equations:


which connect the amplitudes ρ1,2 and n1,2 with the phase Φ; note that, here,


in Cartesian and polar coordinates, respectively. In addition, Eq. (14) yields:


To this end, combining Eq. (29) with Eqs. (26) and (27), yields the following biderectional nonlinear dispersive wave equation for Φ:


where C2 is the squared speed of sound given by Eq. (18) and the parameter α is given by Eq. (19). Here it is worth observing the following. First, at the leading-order in ε, Eq. (30) is the standard second-order wave equation, while at order Oε, the linear part of Eq. (30) corresponds to the dispersion relation (16) for the small-amplitude linear waves of Eqs. (2) and (3) propagating on top of the steady state with u=u0 and n=u02. The full Eq. (30), incorporates fourth-order dispersion and quadratic nonlinear terms, resembling the Boussinesq and Benney-Luke [24] equations. These models, are used to describe bidirectional shallow water waves, in the framework of small-amplitude and long-wave approximations [4].

2.2.2 Long-time behavior – KP equation

Using a multiscale expansion method, similar to the one employed in the water wave problem [4], we now derive the KP equation, which is obtained under the additional assumptions of quasi-two-dimensionality and unidirectional propagation. In particular, we introduce the asymptotic expansion:


where the unknown functions Φ (=0,1,) depend on the variables:


for the Cartesian and polar geometry, respectively. Substituting the expansion (31) into Eq. (30), at the leading-order in ε, we obtain the equations:


in Cartesian and polar coordinates respectively. The above equations imply that Φ0 can be expressed as a superposition of two waves. In the Cartesian geometry, these waves are a right-going one, Φ0R, which depends on χ, and a left-going one, Φ0L, depending on χ˜. Similarly, in the polar case, Φ0 can be expressed as superposition of a radially expanding wave (depending on ϱ), and a radially contracting one (depending on ϱ˜). Thus, the solutions of Eqs. (34) and (35) read:


In addition, at order Oε, we obtain the following equations:

4C2Φ1χχ˜= CΦ0χχRΦ0χ˜LΦ0χRΦ0χ˜χ˜L+χ2CΦ0TR+α4Φ0χχχR3C2Φ0χR2C2Φ0YYR+χ˜2CΦ0TL+α4Φ0χ˜χ˜χ˜L+3C2Φ0χ˜˜L2C2Φ0YYL,E38

for the Cartesian geometry, and

4C2Φ1ϱϱ˜= CΦ0ϱϱRΦ0ϱ˜LΦ0ϱRΦ0ϱ˜ϱ˜L+ρ2CΦ0TR+α4Φ0ϱϱϱR3C2Φ0ϱR2CTΦ0ϱ˜R1T2Φ0ΘΘR+ϱ˜2CΦ0TL+α4Φ0ϱ˜ϱ˜ϱ˜L+3C2Φ0ϱ˜L2CTΦ0ϱ˜L1T2Φ0ΘΘR,E39

for the polar geometry. Here, it is important to observe the following. Once Eq. (38) is integrated in χ or χ˜ [and, similarly, Eq. (39) in ϱ or ϱ˜], secular terms arise: the secular therms are those in the square brackets in the right-hand side, because are functions of χ or χ˜ (ϱ or ϱ˜ in polar) alone, not both. Removal of these secular terms leads to two uncoupled nonlinear evolution equations for Φ0R and Φ0L in both geometries. Furthermore, using the equation ΦT=n1 (i.e., the linear part of Eq. (29), together with Eqs. (26) and (27), it is found that the amplitude ρ1 can also be decomposed to a left- and a right-going wave; this means that


with the fields ρ1R and ρ1L satisfying a pair of two uncoupled KP equations. In the Cartesian geometry, the KP equations are of the form:


while in the polar geometry, the relevant equations take the form of the so-called cylindrical KP (cKP), or Johnson’s equations (see Refs. [22, 23]):


2.2.3 Versions of the KP equations

It is now convenient to further simplify the equations above, namely the KP and the cKP equations, in order to express them in their standard form [4]. We consider right-going waves in the Cartesian geometry, and radially expanding ones in the polar geometry, and introduce the following transformations,


which cast the respective KP and cKP equations into the following form:


Here, the parameter σ2 is defined as:


It is now clear that, in the Cartesian case, Eq. (45) includes both versions of the KP equation, KPI and KPII [33]. Indeed, for σ2=1α<0, i.e., for u02>1/4d or d>1/4u02, Eq. (45) is a KPII equation; this corresponds to a small effective surface tension. On the other hand, for σ2=1α>0, i.e., u02<1/4d or d<1/4u02, Eq. (45) is a KPI equation; this corresponds to the case of large effective surface tension (see end of Section 2.1 and discussion below).

In other words, for a fixed degree of nonlocality d, a larger (smaller) background amplitude u0, as defined by the sign of α, corresponds to KPII (KPI); similarly, for a fixed background amplitude u0, a strong (weak) nonlocality, as defined by the above regimes of d, corresponds to KPII (KPI). Notice that in the local limit of d=0, the asymptotic analysis leads only to the KPI model, in which line solitons are unstable; this fact was used to better understand self-focusing and transverse instability of plane dark solitons of the defocusing NLS Equation [19, 20]. Notably, the same parameter α can also been shown to distinguish solutions in 1D [34] and radially symmetric [35] systems; a similar case, pertinent to the polar geometry, will be considered below.

Once again, it is important to highlight the fact that the existence of these regimes resembles the situation occurring in shallow water. In this context, weak surface tension corresponds to σ2=1 in Eq. (45) (i.e., KPII), while strong surface tension is pertinent to σ2=1 (i.e., KPI). Thus, there exists an immediate connection between the original problem with the one of shallow water waves: relatively large (small) background amplitude or nonlocality corresponds to weak (strong) surface tension, leading to KPII (KPI).


3. Approximate soliton solutions

We now employ the results of the above analysis and construct the solution of our original problem. This can be done as follows. Once a solution of Eq. (45) [or (46)] is known, an approximate solution of Eqs. (2) and (3) has the form:


It is important to note that Eq. (48) describes two different types of solitons, namely dark and antidark ones. Indeed, for α<0, i.e., for solutions satisfying the KPII equation (corresponding to the strong nonlocality regime) the solution (48) has the form of a hump on top of the cw background and is, thus, an antidark soliton. On the other hand, for α>0, i.e., for solutions satisfying the KPI equation (corresponding to the weak nonlocality regime) the solution (48) has the form of a dip on top of the cw background and is, thus, a dark soliton. Note that in the local limit, with d=0, it turns out that α=1 and the solutions have always the form of dark solitons. In other words, antidark solitons are only supported by a sufficiently strong nonlocality.

Below, we will use known solutions of the KP and cKP equations and present corresponding solutions of the nonlocal NLS. In addition, we will perform direct numerical simulations to examine the evolution and – in some cases – the interaction dynamics of the approximate soliton solutions (48).

3.1 Antidark stripe solitons and their interactions

We start with the Cartesian geometry and consider, in particular, the simplest soliton solutions of the KP equation, the so-called line solitons. Specifically, the one-line soliton solution of Eq. (45), traveling at an angle to the Y-axis, reads [33]:


where κ, λ and δ are free parameters. We are particularly interested in the case where the corresponding solution of the NLS takes the form of an antidark soliton, with U obeying the KPII equation: since, in this case, line solitons are stable, we may expect that their interactions will give rise to patterns resembling those observed in shallow water [5]. An example of such an antidark soliton is ahown in Figure 1, for t=0. For this example, the following parameter values were used: u0=1 and d=1/3 (leading to α=1/3<0), as well as κ=λ=1, δ=0, and ε=0.2. These values will also be used below.

Figure 1.

A typical single antidark soliton solution of the nonlocal NLS (uxy0), with κ=λ=1, δ=0; other parameter values are: u0=C=1, d=1/3, and ε=0.2.

In what follows, we choose line solitons with specific parameters, so that their angle during interaction will determine the type of pattern. In the simulations, we evolve the initial configuration up to t=600, so that solitons have enough time to interact and generate interesting interaction patterns. In all cases, the numerical integration of Eqs. (2) and (3) is performed via a high-accuracy pseudo-spectral method in space and fourth-order in time. To adjust to the periodic boundary conditions, the scheme requires we place the solutions on top of a wide super-Gaussian background that decays at infinity and makes the initial condition periodic. As it will be seen below, interacting antidark NLS solitons follow closely the KP dynamics, giving rise to patterns that are usually observed in water.

We start with the interaction of two solitons, and express the two-line soliton solution of KPII as follows:


with ZiκiX+λiYκi2+3λi2T+δi.

We first focus on interacting line solitons that form Y- and X-shaped waveforms, and proceed with more complex ones later in the text. In Figure 2, we show a typical Y-shaped pattern, whereby the resonant two-soliton interaction gives rise to the emergence of one soliton, with its maximum height being four times that of the incoming solitons (this resonant Y-shaped solution was originally found by Miles [7, 8]). This feature is clearly observed in the figure. Parameter values in this case are: κ1=κ2/2=1/2, λ1=3λ2=3/4 (and for all cases, also pertaining to X-type interactions below, δ1=δ2=0).

Figure 2.

A typical Y-type interaction, uxy0, whereby out of two incoming colliding solitons, a single line soliton emerges. Parameter values are: κ1=κ2/2=1/2, λ1=3λ2=3/4.

Next, we proceed with the X-type interactions, which can be discriminated according to the resulting “stem”: an interaction with a short stem (Figure 3-left), an interaction with a long stem where the stem height is higher than the incoming line solitons (Figure 3-middle and right respectively).

Figure 3.

X-type interactions, uxy0, with: κ1=κ2=1/2, λ1=λ2=2/3 (left), κ1=κ2=1/2, λ1=1/40.01, λ2=3/4 (middle), κ1=κ2=1/2, λ1=1/2, λ2=1/21010 (right). Observe that the length of the stem is increased from left to right (see also text).

In the same category falls an X-type interaction with a long stem, with the stem height being lower than the tallest incoming line soliton, which we refer to here as an H-interaction. A typical such interaction is given in Figure 4, and parameter values: κ1=κ2/2=1/2, λ1=1/2107, and λ2=0.

Figure 4.

An H-type interaction, uxy0, with κ1=κ2/2=1/2, λ1=1/2107, λ2=0.

An immediate generalization of the above is a three-wave interaction, where now




so that the resulting interaction is as in Figure 5, for t=200. Here, κ1=1, κ2=2, κ3=3, λ1=λ2/2=λ3/5=1/3 and, as before, δ1=δ2=δ3=0.

Figure 5.

A 2–3 wave interaction, uxy0, with κ1=1, κ2=2, κ3=3, λ1=λ2/2=λ3/5=1/3 and δ1=δ2=δ3=0.

More intricate structures may still form whose mathematical description is more complex. Such novel weblike structures for the KP equation, termed N-in M-out, have been found using Wronskian methods [11, 13, 15] and the so-called τ-function. We illustrate the basic construction of these solitons but do not go into further detail as this is beyond the scope of this work. The most convenient way to write the line soliton solutions of the KPII system is through the Wronskian, or τ-function, which replaces the F function above. Now, the solutions of the KPII equation are written in a similar manner as above, as


where this new τ-function represents the Wronskian determinant


Here superscripts denote differentiation with respect to X, and the set of fn functions constitute the set of linearly independent solutions of the system:


In particular, for line solitons these are defined as:


where θm=kmX+km2Y+km3T+θ0m with distinct real parameters km ‘s with the property: k1<k2<<kM. The parameters θ0m are real constants. Importantly, the coefficients anm define an N×M matrix of rank N, due to the linear independence of the functions fn, such that anmA. Below we provide two examples, and their evolution for a 2–2 and 3–3 interaction (see Figures 6 and 7, respectively).

Figure 6.

A typical 2–2 wave evolution for t=0 (left), t=200 (middle) and t=400 (right). All figures refer to the variable uxyt.

Figure 7.

A typical 3–3 wave evolution, for t=0 (left), t=400 (middle), and t=600 (right). All figures refer to the variable uxyt.

We start with a 2–2 interaction, by defining the matrix A, such that


and the relative real parameters k1=1, k2=0, k3=1 and k4=2 corresponding to the relative θmm=14. Then, the evolution of this interaction is depicted in Figure 6.

Similarly, for a 3–3 wave interaction one has:


and k1=1, k2=1/2, k3=0, k4=1/2, k5=1 and k6=3/2 to produce the following interaction pattern of Figure 7.

It is important to note here that the matrices A for the 2–2 and 3–3 solitons, are not arbitrarily chosen. The N×N maximal minors of A must be non-negative in order for the KP solution q to be regular for all X, Y and T. In the literature, such matrices are called totally non-negative matrices.

3.2 Dark lump solitons

So far, we have focused on structures obeying the effective KPII equation. These structures, namely the antidark stripe solitons, are quasi one-dimensional (1D) states, which are stable. The respective states that satisfy the effective KPI equation, namely dark stripe solitons, are expected to be transversely unstable – see, e.g., Ref. [36] for a stability analysis of KP line solitons.

Nevertheless, the KPI equation, that corresponds to the weakly nonlocal regime, supports stable, purely 2D, soliton states known as lumps. These solitons are weakly localized, namely they decay algebraically as x,y. The respective solitons of the nonlocal NLS, which are of the dark type, i.e., they have the form of dark lumps, can be effectively described as solutions of Eq. (45) for the field U:


where q and k are free real parameters. Notice that solutions of this type have not yet been observed in water due to the fact that the surface tension is small.

Similarly to the case of antidark stripe solitons, once Eq. (60) is substituted into Eq. (48) gives rise to the approximate dark lump soliton of the nonlocal NLS. Using the relevant analytical expression, it is straightforward to prepare the corresponding initial condition by setting t=0, and numerically integrate Eqs. (2) and (3) to examine the evolution of dark lumps. In Figure 8, shown is the result of such a simulation for the dark lump; for this simulation, we have used the parameter values u0=1, d=1/5 (so that α=1/5>0, corresponding to dark soliton states), as well as k=0 and q=1. It can readily be observed that the dark lump solitons propagate undistorted up to t=20, which is the “numerical horizon” for this particular simulation. This result indicates that, similarly to antidark stripe solitons, dark lumps can also be supported by the nonlocal NLS system.

Figure 8.

A typical dark lump solution, at t=0, with k=0, q=1 (top panel) and its evolution for t=0, 10 and 20. The initial condition shown in the top panel is given by Eq. (60).

3.3 Ring dark and antidark solitons

We now turn our attention to the cKP Eq. (46), and focus on radially symmetric solutions that do not depend on Θ. In such a case, Eq. (46) reduces to a quasi-1D equation, namely the cylindrical Korteweg-de Vries (cKdV), which is also a completely integrable model. The cKdV is a generic equation that can been used to describe cylindrical solitons in shallow water, plasmas, etc. [1]. As was shown in Ref. [37], the cKdV equation admits an exact soliton solution which, in the framework of Eq. (46), can be expressed as follows:


where η and ϱ0 are free parameters of the soliton. It is clear that Eq. (61) has the form of a sech2-pulse on top of a rational background. A clearer picture regarding the structure of the cKdV equation, can be obtained by means of an asymptotic asymptotic analysis [38, 39]. Indeed, as shown in these works, to leading-order, in the regime Tϱ, the primary wave UϱT, decays to zero at both upstream and downstream infinity. The primary wave has a form similar to that of Eq. (50), but with the following change of variables: Xϱ and δϱ0. However, there is a very important difference: κ now becomes a slowly-varying function of T, due to the presence of the term U/2T in the cKdV. Following the analysis of [38, 39], and using the original coordinates, it can be concluded that:


where κ02 is a constant setting the solitary wave amplitude at t=t0. We can know express an approximate solution of Eqs. (2) and (3), for the polar chase, and for the primary solitary wave. This is of the form of Eq. (48), with the soliton amplitude and velocity varying as t2/3, and the width varying as t1/3, as follows from Eqs. (50) and (62).

Notice that this approximate solution is a ring-shaped solitary wave, on top of cw background, which can be either dark (for α>0) or anti-dark (for α<0). Notice, also that ring dark solitons were predicted to occur in optical media exhibiting either Kerr [40] or non-Kerr [41] nonlinearities, and were later observed in experiments [42]. On the other hand, ring anti-dark solitons were only predicted to occur in non-Kerr – e.g., saturable media [41, 43]. This picture is complemented by this analysis, according to which a relatively strong [i.e., d>1/4u02] nonlocal nonlinearity can also support ring anti-dark solitary waves.

As in the Cartesian case, we have performed direct numerical simulations to examine the evolution of ring-shaped solitons in the nonlocal NLS model. First, in Figure 9, shown are 3D plots depicting the profiles of the ring dark solitons and ring antidark solitons. Here, we use d=1/5 for the dark soliton (α>0) and d=1/3 for the antidark (α<0), while the rest of the parameters are set equal to unity except ε=0.2 as before. Furthermore, both solitons have an initial radius of r0=10.

Figure 9.

(Color online) typical ring dark (left) and anti-dark (right) solitons.

In addition, in Figure 10, we show the evolution of the ring solitons depicted in Figure 9. Once again, it is observed, that both types of ring solitons propagate undistorted up to end of the simulation (t=10), which indicates that these states are also supported by the nonlocal NLS.

Figure 10.

(Color online) the evolution of the dark (top) and anti-dark (bottom) ring solitons of Figure 9.


4. Conclusions and discussion

In conclusion, we have studied a 2D defocusing NLS model with a spatially nonlocal nonlinearity. The considered model is relevant to optical beam propagation in a variety of physical settings, including thermal media, plasmas, and nematic liquid crystals. We have analytically treated this nonlocal NLS by means of a two-stage multiscale analysis. This led to a Boussinesq/Benney-Luke-type equation, and a KP equation (in both Cartesian and polar geometries), which are models arising in the context of shallow water waves. This analysis revealed that the nonlocal NLS is characterized by an effective surface tension, which is related to the degree of nonlocality: small (large) effective surface tension corresponds to the strong (weak) nonlocality regime, whereby the asymptotic reduction leads to a KPII (KPI) equation.

Apart from the type of KP (KPII or KPI), the effective surface tension also controls the type of soliton: antidark (dark) solitons on top (off) the cw pedestal can only be formed in the strong (weak) nonlocality regime. We have thus found a number of different soliton states that can be supported by the nonlocal NLS. These include antidark stripe solitons, dark lump solitons, as well as ring dark and antidark solitons (in the polar case). We also used numerical simulations to examine the existence, evolution, as well as interaction dynamics. For the latter, we focused on the case of antidark solitons, which were shown to form a plethora of patterns coming out of their collisions. These include Y-, X-, H-shaped waveforms, as well as more complicated patterns involving more than two solitons; all these, have been observed in shallow water [5], a fact that highlights the deeper connection between water waves and light! In addition, we examined, by means of direct simulations, the existence and propagation of other soliton states, namely dark lump solitons, as well as ring dark and antidark solitons. We found that all these states do exist in the framework of the nonlocal NLS, and can propagate undistorted up to times set by the simulations.

At this point, it is also relevant to comment on the possibility of observing the predicted patterns in real experiments. In that regard, first we consider the case of antidark solitons. The observation of these states took place some time ago, as it was first reported in Ref. [44]; this fact indicates that the necessary experimental setting is already available. As concerns the experimental set up needed for the observation of the Y-, X- and H-waves and other related patterns, it may be similar to that of Ref. [44]. In particular, one may employ at first a cw laser beam, which is split into two parts via a beam-splitter. One branch goes through a cavity system to form a pulse (as happens in typical pulsed lasers); this pulse branch undergoes phase-engineering, i.e., passes through a phase mask so that the characteristic phase jump of the antidark soliton is inscribed. Then, the cw and the phase-engineered pulse are incoherently coupled inside the nonlocal medium, e.g., a nematic liquid crystal. This process forms one antidark soliton, as in Ref. [44]. To observe the Y-, X- or H-patterns predicted above, two such antidark solitons have to be combined inside the crystal. The angle between the two incident beams, which should be appropriately chosen so that a specific pattern be formed, can be controlled by a rotating mirror in one of the branches of the beam-splitter.

As concerns the possibility of the experimental observation of the other soliton states that were predicted by our analysis, we note the following. First, regarding the weakly localized dark lump solitons, they can be experimentally observed as follows. In a real experiment, it is straightforward to create a dark soliton by using, e.g., a proper phase mask (see, e.g., the review [45] and references therein); this can also be done even in the case of ring dark solitons [46]. Such structures are prone to the transverse modulational instability in 2D, and are known to decay into vortices [45]. However, as was predicted in Ref. [19], sufficiently weak gray solitons can decay to dark lumps satisfying an effective KPI equation. This result suggests that, to observe dark lumps, one only need to embed a shallow dark soliton (like the ones we considered herein) into a 2D space, e.g., in a planar waveguide; then, the transverse instability will set in, and dark lump solitons will emerge. Finally, regarding the ring solitons, as mentioned above, they have already been observed in a variety of experiments (see, e.g., Ref. [46] and the review [45]) and, thus they can readily be created, e.g., in thermal media or in nematic liquid crystals. Thus, a variety of patterns that can be observed in water may also emerge in optics.

The obtained results relied on the formal reduction of the nonlocal NLS model to the KP (cKP) equation. It would be interesting to investigate if other types of nonlocality, arising, e.g., in the context of chromium Bose-Einstein condensates [47], or in defocusing colloidal media [48] can support relevant soliton states. Such studies are in progress and pertinent results will be presented elsewhere.


  1. 1. Infeld E, Rowlands G. Nonlinear Waves, Solitons, and Chaos. Cambridge University Press; Cambridge, 1990
  2. 2. Dauxois T, Peyrard M. Physics of Solitons. Cambridge University Press; Cambridge, 2006
  3. 3. Ablowitz MJ, Clarkson PA. Solitons, nonlinear evolution equations and inverse scattering. Cambridge University Press; Cambridge, 1991
  4. 4. Ablowitz MJ. Nonlinear dispersive waves: Asymptotic analysis and solitons. Cambridge University Press; Cambridge, 2011
  5. 5. Ablowitz MJ, Baldwin DE. Nonlinear shallow ocean-wave soliton interactions on flat beaches. Phys Rev E. 2012;86:036305
  6. 6. Kadomtsev BB, Petviashvili VI. On the stability of solitary waves in weakly dispersing media. Sov Phys Dokl. 1970;15:539–541
  7. 7. Miles JW. Obliquely interacting solitary waves. J Fluid Mech. 1977;79:157–169
  8. 8. Miles JW. Resonantly interacting solitary waves. J Fluid Mech. 1977;79:171–179
  9. 9. Newell AC, Redekopp LG. Breakdown of Zakharov-Shabat Theory and Soliton Creation. Phys Rev Lett. 1977;38:377–380
  10. 10. Ohkuma K, Wadati M. The Kadomtsev-Petviashvili Equation: the Trace Method and the Soliton Resonances. J Phys Soc Jpn. 1983;52:749–760
  11. 11. Biondini G, Kodama Y. On a family of solutions of the Kadomtsev-Petviashvili equation which also satisfy the Toda lattice hierarchy. J Phys A: Math Gen. 2003;36:10519
  12. 12. Biondini G, Chakravarty S. Soliton solutions of the Kadomtsev-Petviashvili II equation. J Math Phys. 2006;47:033514
  13. 13. Chakravarty S, Kodama Y. Classification of the line-soliton solutions of KPII. J Phys A: Math Theor. 2008;41:275209
  14. 14. Chakravarty S, Kodama Y. Soliton Solutions of the KP Equation and Application to Shallow Water Waves. Stud Appl Math. 2009;123:83–151
  15. 15. Chakravarty S, Lewkow T, Maruno KI. On the construction of the KP line-solitons and their interactions. Applicable Analysis. 2010;89:529–545
  16. 16. McDowell T, Osborne M, Chakravarty S, Kodama Y. On a class of initial value problems and solitons for the KP equation: A numerical study. Wave Motion. 2017;72:201–227
  17. 17. Kivshar YS, Agrawal GP. Optical Solitons: From Fibers to Photonic Crystals. Academic Press; New York, 2003
  18. 18. Kuznetsov EA, Turitsyn SK. Instability and collapse of solitons in media with a defocusing nonlinearity. JETP. 1988;67:1583–1588
  19. 19. Pelinovsky DE, Stepanyants YA, Kivshar YS. Self-focusing of plane dark solitons in nonlinear defocusing media. Phys Rev E. 1995;51:5016–5026
  20. 20. Kivshar YS, Pelinovsky DE. Self-focusing and transverse instabilities of solitary waves. Phys Rep. 2000;331:117–195
  21. 21. Horikis TP, Frantzeskakis DJ. Light Meets Water in Nonlocal Media: Surface Tension Analogue in Optics. Phys Rev Lett. 2017;118:243903
  22. 22. Johnson RS. Water waves and Korteweg-de Vries equations. J Fluid Mech. 1980;97:701–719
  23. 23. Johnson RS. A modern introduction to the mathematical theory of water waves. Cambridge University Press; 1997
  24. 24. Benney DJ, Luke JC. On the interactions of permanent waves of finite amplitude. J Math and Phys. 1964;43:309–313
  25. 25. Rotschild C, Cohen O, Manela O, Segev M, Carmon T. Solitons in nonlinear media with an infinite range of nonlocality: first observation of coherent elliptic solitons and of vortex-ring solitons. Phys Rev Lett. 2005;95:213904
  26. 26. Krolikowski W, Bang O, Nikolov NI, Neshev D, Wyller J, Rasmussen JJ, et al. Modulational instability, solitons and beam propagation in spatially nonlocal nonlinear media. J Opt B: Quantum Semiclass Opt. 2004;6:S288–S294
  27. 27. Ghofraniha N, Conti C, Ruocco G, Trillo S. Shocks in nonlocal media. Phys Rev Lett. 2007;99:043903
  28. 28. Conti C, Fratalocchi A, Peccianti M, Ruocco G, Trillo S. Observation of a gradient catastrophe generating solitons. Phys Rev Lett. 2009;102:083902
  29. 29. Litvak AG, Mironov VA, Fraiman GM, Yunakovskii AD. Thermal self-effect of wave beams in a plasma with a nonlocal nonlinearity. Sov J Plasma Phys. 1975;1:60–71
  30. 30. Yakimenko AI, Zaliznyak YA, Kivshar YS. Stable vortex solitons in nonlocal self-focusing nonlinear media. Phys Rev E. 2005;71:065603(R)
  31. 31. Conti C, Peccianti M, Assanto G. Route to nonlocality and observation of accessible solitons. Phys Rev Lett. 2003;91:073901
  32. 32. Assanto G. Nematicons: Spatial Optical Solitons in Nematic Liquid Crystals. New Jersey: Wiley-Blackwell; 2012
  33. 33. Ablowitz MJ, Clarkson PA. Solitons, nonlinear evolution equations and inverse scattering. Cambridge University Press; Cambridge, 1991
  34. 34. Horikis TP. Small-amplitude defocusing nematicons. J Phys A: Math Theor. 2015;48:02FT01
  35. 35. Horikis TP, Frantzeskakis DJ. Ring dark and antidark solitons in nonlocal media. Opt Lett. 2016;41:583–586
  36. 36. Ablowitz MJ, Segur H. Solitons and the inverse scattering transform. SIAM; 1981
  37. 37. Hirota R. Exact solutions to the equation describing “cylindrical solitons”. Phys Lett A. 1979;71:393–394
  38. 38. Johnson RS. A note on an asymptotic solution of the cylindrical Korteweg-de Vries equation. Wave Motion. 1999;30:1–16
  39. 39. Ko K, Kuehl HH. Cylindrical and spherical Korteweg-deVries solitary waves. Phys Fluids. 1979;22:1343–1348
  40. 40. Kivshar YS, Yang X. Ring dark solitons. Phys Rev E. 1994;50:R40–R43
  41. 41. Frantzeskakis DJ, Malomed BA. Multiscale expansions for a generalized cylindrical nonlinear Schrödinger equation. Phys Lett A. 1999;264:179–185
  42. 42. Dreischuh A, Neshev D, Paulus GG, Grasbon F, Walther H. Ring dark solitary waves: Experiment versus theory. Phys Rev E. 2002;66:066611
  43. 43. Nistazakis HE, Frantzeskakis DJ, Malomed BA, Kevrekidis PG. Head-on collisions of ring dark solitons. Phys Lett A. 2001;285:157–164
  44. 44. Coskun TH, Christodoulides DN, Kim YR, Chen Z, Soljacic M, M S. Bright Spatial Solitons on a Partially Incoherent Background. Phys Rev Lett. 2000;84:2374–2377
  45. 45. Kivshar YS, Luther-Davies B. Dark optical solitons: physics and applications. Phys Rep. 1998;298:81–197
  46. 46. Dreischuh A, Neshev D, Paulus GG, Grasbon F, Walther H. Ring dark solitary waves: Experiment versus theory. Phys Rev E. 2002;66:066611
  47. 47. Lahaye T, Menotti C, Santos L, Lewenstein M, Pfau T. The physics of dipolar bosonic quantum gases. Rep Prog Phys. 2009;72:126401
  48. 48. An X, Marchant TR, Smyth NF. Optical dispersive shock waves in defocusing colloidal media. Physica D. 2017;342:45–56

Written By

Georgios N. Koutsokostas, Theodoros P. Horikis, Dimitrios J. Frantzeskakis, Nalan Antar and İlkay Bakırtaş

Submitted: 19 March 2020 Reviewed: 08 December 2020 Published: 07 January 2021