Open access peer-reviewed chapter - ONLINE FIRST

Bright, Dark, and Kink Solitary Waves in a Cubic-Quintic-Septic-Nonical Medium

By Mati Youssoufa, Ousmanou Dafounansou and Alidou Mohamadou

Submitted: February 3rd 2020Reviewed: May 13th 2020Published: June 19th 2020

DOI: 10.5772/intechopen.92819

Downloaded: 33

Abstract

In this chapter, evolution of light beams in a cubic-quintic-septic-nonical medium is investigated. As the model equation, an extended form of the well-known nonlinear Schrödinger (NLS) equation is taken into account. By the use of a special ansatz, exact analytical solutions describing bright/dark and kink solitons are constructed. The existence of the wave solutions is discussed in a parameter regime. Moreover, the stability properties of the obtained solutions are investigated, and by employing Stuart and DiPrima’s stability analysis method, an analytical expression for the modulational stability is found.

Keywords

  • higher-order nonlinear Schrödinger equation
  • spatial solitons
  • stability analysis method
  • modulational instability
  • optical fibers

1. Introduction

The study of spatial solitons in the field of fiber-optical communication has attracted considerable interest in recent years. In a uniform nonlinear fiber, soliton can propagate over relatively long distance without any considerable attenuation. The formation of optical solitons in optical fibers results from an exact balancing between the diffraction and/or group velocity dispersion (GVD) and the self-phase modulation (SPM). The theorical prediction of a train of soliton pulses from a continuous-wave (CW) light in optical fibers was first suggested by Hasegawa and Tappert [1, 2] and first experimentally demonstrated by Mollenauer et al. [3] in single-mode fibers in the case of negative GVD, in liquid CS2by Barthelemy et al. in 1985 [4]. In nonlinear optic, optical solitons are localized electromagnetic waves that transmit in nonlinear Kerr or non-Kerr media with dispersion or (and) diffraction without any change in shapes. In nonlinear media, the dynamics of spatial optical solitons is governed by the well-known nonlinear Schrödinger (NLS) equation. Depending on the signs of GVD, the NLS equation admits two distinct types of soliton, namely, bright and dark solitons. The bright soliton exists in the regime of anomalous GVD, and the dark soliton arises in the regime of normal GVD. The physics governing the soliton differs depending on whether one considers a bright or a dark soliton and accordingly features distinct applications [5, 6, 7, 8]. The unique property of optical solitons, either bright or dark, is their particle-like behavior in interaction [9].

In addition to fundamental bright and dark solitons, various other forms and shapes of solitary waves can appear in nonlinear media. Kink solitons, for example, are an important class of solitons which may propagate in nonlinear media exhibiting higher-order effects such as third-order dispersion, self-steepening, higher-order nonlinearity, and intrapulse stimulated Raman scattering. In the setting of nonlinear optics, a kink soliton represents a shock front that propagates undistorted inside the dispersive nonlinear medium [10]. This type of solitons has been studied extensively, both analytically and numerically [11, 12, 13]. These spatial soliton solutions can maintain their overall shapes but allow their widths and amplitudes and the pulse center to change according to the management of the system’s parameters, such as the dispersion, nonlinearity, gain, and so on [14].

The cubic nonlinear Schrödinger equation (CNLSE) has been widely used to model the propagation of light pulse in material’s systems involving third-order susceptibility χ3, though, for moderate pulse intensity, the higher-order nonlinearities are related to higher-order nonlinear susceptibilities (nonlinear responses) of a material. For example, the cubic-quintic-nonlinear Schrödinger equation (CQNLSE) models materials with fifth-order susceptibility χ5. This kind of nonlinearity (cubic-quintic CQ) is named as parabolic law nonlinearity and existing in nonlinear media such as the p-toluene sulfonate (PTS) crystals. The parabolic law can closely describe the nonlinear interaction between the high-frequency Langmuir waves and the ion acoustic waves by ponderomotive forces [15, 16], in a region of reduced plasma density, and the nonlinear interaction between Langmuir waves and electrons. In addition, CQ was experimentally proposed as an empirical description of special semiconductor (e.g., AlGaAs, CdS, etc.) waveguides and semiconductor-doped glasses, particularly for the CdSxSe1x-doped glass, which exhibit a significant fifth-order susceptibilities χ5as experimentally reported earlier [17, 18]. Moreover, using high laser intensity, the saturation of nonlinearity has been established experimentally in many materials such as nonlinear organic polymers, semiconductor-doped glasses, and so on, which have the property that their absorption coefficient decreases [19]. More generally, a self-defocusing χ5usually accounts for the saturation of χ3.

In recent years, many influential works have devoted to construct exact analytical solutions of CQNLSE, such as the pioneering work of Serkin et al. [20]. In particular, Dai et al. [21, 22, 23, 24, 25] obtained exact self-similar solutions (similaritons), their nonlinear tunneling effects of the generalized CQNLSE, and their higher-dimensional forms with spatially inhomogeneous group velocity dispersion, cubic-quintic nonlinearity, and amplification or attenuation.

Since the measurement of third-, fifth-, and seventh-order nonlinearities of silver nanoplatelet colloids using a femtosecond laser [26], an extension of nonlinear Schrödinger equation including the cubic-quintic-septic nonlinearity was used to model the propagation of spatial solitons. In [27], for example, the authors performed numerical calculations based on higher-order nonlinearity parameters including seventh-order susceptibility χ7(a chalcogenide glass is an example). This seeds several motivations to discover new features of solitons with combined effects of higher-order nonlinear parameters. In this regard, Houria et al. [28] constructed dark spatial solitary waves in a cubic-quintic-septic-nonlinear medium, with a profile in a functional form given in terms of “sech23”. They have also investigated chirped solitary pulses for a derivative nonical-NLS equation on a CW background [29]. It is obvious to notice that the contributions of the higher-order nonlinearities can give way to generate stable solitons in homogeneous isotropic media and influence many aspects of filamentation in gases and condensed matters [30, 31, 32, 33].

Recently, the study of modulational instability (MI) in non-Kerr media has receiving particular attention. MI is a fundamental and ubiquitous process that appears in most nonlinear systems [6, 9, 34, 35, 36, 37]. This instability is referred to as modulation instability because it leads to a spontaneous temporal modulation of the CW beam and transforms it into a pulse train. During this process, small perturbations upon a uniform intensity beam grow exponentially due to the interplay between nonlinearity and dispersion or diffraction. As a result, under specific conditions, a CW light often breaks up into trains of ultrashort solitons like pulses [9]. To date, there has not been any report of MI in the cubic-quintic-septic-nonical-nonlinear Schrödinger equation (CQSNNLSE).

Our study will be focused on the analysis of solitary wave’s solutions of systems described by the higher-order NLSE named CQSNNLSE. We will discuss the model with higher-order nonlinearities and explore the dynamics of bright, dark, and kink soliton solutions. Finally, the linear stability analysis of the MI is formulated, and the analytical expression of the gain of MI is obtained. Moreover, the typical outcomes of the nonlinear development of the MI are reported.

2. Model equation

The dynamics of (1 + 1)-dimensional (one spatial and one temporal variables) spatial optical solitons is the well-known nonlinear Schrödinger equation. If we consider the higher-order effects, an extended model is required, and the propagation of optical pulses through the highly nonlinear waveguides can be described by the CQSNNLSE:

Ez=iα1Ett+iα2|E|2E+iα3|E|4E+iα4|E|6E+iα5|E|8E,E1

where Eztis the slowly varying envelope of the electric field, the subscripts zand tare the spatial and temporal partial derivatives in the frame moving with the pulsed solutions, α1is the parameter of diffraction or dispersion, and α2, α3, α4, and α5are the cubic, quintic, septic, and nonical nonlinear terms, respectively. This model is relevant to some applications in which higher-order nonlinearities are important.

For example, Eq. (1) with α1=12, α4=1, and α5=0was used to study numerically the stability conditions of one-dimensional spatial solitons [38]. Recently, Eq. (1) with α5=0was analyzed for systems that are valid for several types of septic nonlinear materials [28]. Here, we consider arbitrary parameters αjj=12345for the sake of a general analysis that is valid for several types of nonical media.

To obtain the exact analytic optical solitary-wave solutions of Eq. (1), we can employ the following transformation:

Ezt=θt+βzeikzωt=θζeikzωt.E2

Here, θζis a real function and βis a real constant to be determined.

Upon substituting Eq. (2) into Eq. (1) and separating the real and imaginary parts, one obtains

β=2α1ω,E3
θζζ=k+α1ω2α1θα2α1θ3α3α1θ5α4α1θ7α5α1θ9,E4

Eq. (4) represents the evolution of an anharmonic oscillator with an effective potential energy V[28] defined by

Vθ=k+α1ω22α1θ2+α24α1θ4+α36α1θ6+α48α1θ8+α510α1θ10.E5

Integrating Eq. (4) yields

θζ2=a1θ2a2θ4a3θ6a4θ8a5θ10+2ξ,E6

where

a1=k+α1ω22α1,an=αnnα1n=2345,E7

and ξis the constant of integration, which can represent the energy of the anharmonic oscillator [39].

In order to get the exact soliton solutions, we first rewrite Eq. (6) in a simplified form by using transformation:

θζ=u12ζ.E8

By substituting Eq. (8) into Eq. (6), we obtain a new auxiliary equation possessing a sixth-degree nonlinear term:

14uζ2=a1u2a2u3a3u4a4u5a5u6+2ξ.E9

To solve Eq. (9), we will employ three types of localized solutions named bright, dark, and kink solitons. In the following, we solve Eq. (9) by using appropriate ansatz and obtain alternative types of solitary-wave solutions on a CW background and investigate parameter domains in which these optical spatial solitary waves exist.

3. Exact solitary-wave solutions

In this section, we find bright-, dark-, and kink-solitary-wave localized solutions of Eq. (9), by using a special ansatz:

3.1 Bright solitary-wave solutions

The bright solitary solutions of Eq. (9) have the form:

ubζ=Ab1+Nbcoshαbζ,E10

where Ab, Nb, and αbare real constants which represent wave parameters (Aband αbrelated to the amplitude and pulse width of the bright wave profiles, respectively) to be determined by the physical coefficients of the model.

Substituting the ansatz Eq. (10) into Eq. (9), we obtain the unknown parameters Ab, Nb, αb, and energy ξ:

Ab=2a1a3αb=a14Nb=a32+4a1a5a32ξ=a1a22a3,E11

with parametric conditions

a2a3+4a1a4=0,a1>0,a3>0,a32+4a1a5>0.E12

Thus, the exact bright solitary-wave solutions on a CW background of Eq. (1) are of the form:

Ebzt=2a1a31+a32+4a1a5a32cosha14t+2α1ωz14×eikzωt.E13

3.2 Dark solitary-wave solutions

The dark solitary solutions of Eq. (9) take the form [40]:

udζ=Adsinhαdζ1+Ndsinh2αdζ.E14

Here Ndis a real constant supposed to be positive. Real parameters Adand αdare related to the amplitude and pulse width of the dark wave profiles, respectively.

By substituting the ansatz Eq. (14) into Eq. (9), we get the unknown parameters Ad, αd, and energy ξ:

Ad=2a1Nda3αd=a12ξ=a12a3,E15

subject to the parametric conditions

2a2a3+2a1a4a3=0,a1>0,a3>0,a32+4a1a5>0.E16

The exact dark solitary-wave solutions on a CW background of Eq. (1) are of the form:

Edzt=2a1Nda3sinh2a12t+βz1+Ndsinh2a12t+2α1ωz14×eikzωt.E17

3.3 Kink solitary-wave solutions

The kink solitary solutions of Eq. (9) are in the following form:

ukζ=Ak1+tanhαkζ,E18

where Akand αkare real parameters related to the amplitude and pulse width of the kink wave profiles, respectively.

Substituting Eq. (18) into Eq. (9), we get

Ak=a34a5αk=a324a5ξ=a224a4,E19

under the parametric conditions

2a2a5a3a4=0,a32+4a1a5=0,a3>0,a5<0.E20

Thus, the exact bright solitary-wave solutions on a CW background of Eq. (1) are of the form:

Ekzt=a34a51+tanha324a5t+2α1ωz14×eikzωt.E21

The previous three exact solitary-wave solutions (13), (17), and (21) exist for the governing nonical-NLS model due to a balance among diffraction (or dispersion) and competing cubic-quintic-septic-nonical nonlinearities. For better insight, we plot in Figure 1 the intensity profile on top of the related first two exact solution solitons named bright and dark, corresponding to the CQNLS models (with α4=0,α5=0) that is available in the current literature. As we can see from Eq. (21), the kink solitons exist only if a50, consequently α50; thus, we cannot plot the corresponding CQNLS kink solution.

Figure 1.

Intensity |Ej(z, t)|2 distribution of the (a) bright and (b) dark solitons given by Eqs. (14) and (17), respectively, with the parameter values corresponding to CQNLS models as α1=0.5, α2=0, α3=1, α4=0, α5=0, k=1, and ω=1.

4. Modulational instability of the CW background

One of the essential aspects of solitary waves is their stability on propagation, in particular their ability to propagate in a perturbed environment over an appreciable distance [41]. Unlike the conventional pulses of different forms, the solitons are relatively stable, even in an environment subjected to external perturbations.

The previous three exact solitary-wave solutions given by the expressions (13), (17), and (21) are sitting on a CW background, which may be subject to MI. If this phenomenon occurs, then the CW background will be quickly destroyed, which will inevitably cause the destruction of the soliton. It is therefore of paramount importance to verify whether the condition of the existence of the soliton can be compatible with the condition of the stability of the CW background. Since MI properties can be used to understand the different excitation patterns on a CW in nonlinear systems, in this section, we perform the standard linear stability analysis [9, 34] on a generic CW:

E0zt=P0eiϕnl,E22

in the system modeled by Eq. (1), where ϕnl=P0α2+α3P0+α4P02+α5P03zis the nonlinear phase shift induced by self-phase modulation and non-Kerr quintic-septic-nonical nonlinear terms, P0being the initial power inside a medium exhibiting optical nonlinearities up to the ninth order. A perturbed nonlinear background plane-wave field for the CQSNNLSE (Eq. (1)) can be written as

Ezt=P0+azteiϕnl,E23

where aztis a small perturbation field which is given by collecting the Fourier modes as

azt=a+eiKzΩt+aeiKzΩt,E24

a+and aare much less than the background amplitude P0, and Ωrepresents the perturbed frequency. Here, the complex field |azt|P0. Thus, if the perturbed field grows exponentially, the steady state (CW) becomes unstable. Inserting the expression of a perturbed nonlinear background Eq. (23) into Eq. (1), with respect to Eq. (24), we obtain after linearization the following dispersion relation:

K=|α1||Ω|Ω2sgnα1Ωc2,E25

where Ωc=1α12α2P0+4α3P02+6α4P03+8α5P04and sgnα1=±1depending on the sign of α1sgnα1=+1forα1>0andsgnα1=1forα1<0. The dispersion relation (25) shows that the steady-state stability depends critically on whether the light experiences normal or anomalous GVD inside the fiber. In the case of normal GVD (α1<0), the wave number Kis real for all Ω, and the steady state is stable against small perturbations. By contrast, in the case of anomalous GVD (α1>0), Kbecomes imaginary for Ω<Ωc, and the perturbation aztgrows exponentially with z. As a result, the CW solution E0is inherently unstable for α1>0. This instability is referred to as modulation instability because it leads to a spontaneous temporal modulation of the CW beam and transforms it into a pulse train. Similar instabilities occur in many other nonlinear systems and are often called self-pulsing instabilities [9, 42, 43, 44, 45]. Then, by setting sgnα1=1, one can obtain the MI gain G=2ImK, where the factor 2 convert Gto power gain. The gain exists only if for |Ω|<Ωcand is given by

GΩ=2|α1Ω|Ωc2Ω2.E26

The gain attains its peak values when the modulated frequency reaches its optimum value, i.e., its optimum modulation frequency (OMF). The OMF corresponding to the gain spectrum (26) is given by

Ωop=±Ωc2,E27

and the peak value given by

Gop=GΩop=|α1Ωc2|.E28

In Figure 2, we have shown the variation of OMF, computed from Eq. (27) as a function of the GVD parameter (α1). The parameter values we have used are given as [34]

Figure 2.

Variation of optimum modulation frequency Ωop as a function of second-order dispersion α1.

P0=15W,α2=2736W1/km,α3=2.63W2/km,α4=9.12×104W3/km,α5=0.5W4/km.E29

We can observe that the OMF increases (respectively decreases) with the increasing α1<0(respectively with the increasing α1>0).

Figure 3 shows the variation of MI gain as a function of the nonic nonlinearity α5. The MI gain increases with the decreasing nonic nonlinearity. In Figure 4, as the input power increases, the maximum gain also increases.

Figure 3.

Variation of the MI gain G as a function of the nonic nonlinearity α5, with the same parameter values as in Figure 2.

Figure 4.

Variation of the MI gain Gkm−1 as a function of frequency ΩHz, at a four-power level P0 for an optical fiber. The other parameters are α5=0.5ps2/km,α2=2736W−1/km,α3=2.63W−2/km,α4=−9.12×10−4W−3/km,α5=0.5W−4/km.

The MI gain spectrum in Figure 5 is a constitutive of two symmetrical sidebands which stand symmetrically along the line Ω=0. The maximum gain is nil at the zero perturbation frequency Ω=0; thus, there is no instability at the zero perturbation frequency.

Figure 5.

Variation of the MI gain Gkm−1 as a function of frequency Ω and the GVD α1. The other parameters are P0=15W,α5=−0.5ps2/km,α2=2736W−1/km,α3=2.63W−2/km,α4=−9.12×10−4W−3/km,α5=5W−4/km.

5. Conclusion

In this chapter, we have investigated the higher-order nonlinear Schrödinger equation involving nonlinearity up to the ninth order. We have constructed exact solutions of this equation by means of a special ansatz. We showed the existence of a family of solitonic solutions: bright, dark, and kink solitons. The conditions on the physical parameters for the existence of this propagating envelope have also been reported. These conditions show a subtle balance among the diffraction or dispersion, Kerr nonlinearity, and quintic-septic-nonical non-Kerr nonlinearities, which has a profound implication to control the wave dynamics. Moreover, by employing Stuart and DiPrima’s stability analysis method, an analytical expression for the MI gain has been obtained. The outcomes of the instability development depend on the nonlinearity and dispersion (or diffraction) parameters. Results may find straightforward applications in nonlinear optics, particularly in fiber-optical communication.

Conflict of interest

The authors declare no conflict of interest.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Mati Youssoufa, Ousmanou Dafounansou and Alidou Mohamadou (June 19th 2020). Bright, Dark, and Kink Solitary Waves in a Cubic-Quintic-Septic-Nonical Medium [Online First], IntechOpen, DOI: 10.5772/intechopen.92819. Available from:

chapter statistics

33total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us