Research article

Existence of two homoclinic solutions for a nonperiodic difference equation with a perturbation

  • Received: 25 November 2020 Accepted: 18 February 2021 Published: 25 February 2021
  • MSC : 39A12, 39A23

  • In the present paper, with the combination of a compactness lemma and variational techniques, we establish a sufficient condition to guarantee the existence of two nontrivial homoclinic solutions for a nonperiodic fourth-order difference equation with a perturbation. Our result generalizes and improves some known results.

    Citation: Yuhua Long. Existence of two homoclinic solutions for a nonperiodic difference equation with a perturbation[J]. AIMS Mathematics, 2021, 6(5): 4786-4802. doi: 10.3934/math.2021281

    Related Papers:

    [1] Tianyong Han, Ying Liang, Wenjie Fan . Dynamics and soliton solutions of the perturbed Schrödinger-Hirota equation with cubic-quintic-septic nonlinearity in dispersive media. AIMS Mathematics, 2025, 10(1): 754-776. doi: 10.3934/math.2025035
    [2] Elsayed M. E. Zayed, Mona El-Shater, Khaled A. E. Alurrfi, Ahmed H. Arnous, Nehad Ali Shah, Jae Dong Chung . Dispersive optical soliton solutions with the concatenation model incorporating quintic order dispersion using three distinct schemes. AIMS Mathematics, 2024, 9(4): 8961-8980. doi: 10.3934/math.2024437
    [3] Noha M. Kamel, Hamdy M. Ahmed, Wafaa B. Rabie . Solitons unveilings and modulation instability analysis for sixth-order coupled nonlinear Schrödinger equations in fiber bragg gratings. AIMS Mathematics, 2025, 10(3): 6952-6980. doi: 10.3934/math.2025318
    [4] Yazid Alhojilan, Islam Samir . Investigating stochastic solutions for fourth order dispersive NLSE with quantic nonlinearity. AIMS Mathematics, 2023, 8(7): 15201-15213. doi: 10.3934/math.2023776
    [5] Ninghe Yang . Exact wave patterns and chaotic dynamical behaviors of the extended (3+1)-dimensional NLSE. AIMS Mathematics, 2024, 9(11): 31274-31294. doi: 10.3934/math.20241508
    [6] Islam Samir, Hamdy M. Ahmed, Wafaa Rabie, W. Abbas, Ola Mostafa . Construction optical solitons of generalized nonlinear Schrödinger equation with quintuple power-law nonlinearity using Exp-function, projective Riccati, and new generalized methods. AIMS Mathematics, 2025, 10(2): 3392-3407. doi: 10.3934/math.2025157
    [7] Azzh Saad Alshehry, Safyan Mukhtar, Ali M. Mahnashi . Optical fractals and Hump soliton structures in integrable Kuralay-Ⅱ system. AIMS Mathematics, 2024, 9(10): 28058-28078. doi: 10.3934/math.20241361
    [8] Mohammad Alqudah, Safyan Mukhtar, Albandari W. Alrowaily, Sherif. M. E. Ismaeel, S. A. El-Tantawy, Fazal Ghani . Breather patterns and other soliton dynamics in (2+1)-dimensional conformable Broer-Kaup-Kupershmit system. AIMS Mathematics, 2024, 9(6): 13712-13749. doi: 10.3934/math.2024669
    [9] F. A. Mohammed, Mohammed K. Elboree . Soliton solutions and periodic solutions for two models arises in mathematical physics. AIMS Mathematics, 2022, 7(3): 4439-4458. doi: 10.3934/math.2022247
    [10] Mohammed Aldandani, Abdulhadi A. Altherwi, Mastoor M. Abushaega . Propagation patterns of dromion and other solitons in nonlinear Phi-Four (ϕ4) equation. AIMS Mathematics, 2024, 9(7): 19786-19811. doi: 10.3934/math.2024966
  • In the present paper, with the combination of a compactness lemma and variational techniques, we establish a sufficient condition to guarantee the existence of two nontrivial homoclinic solutions for a nonperiodic fourth-order difference equation with a perturbation. Our result generalizes and improves some known results.



    A soliton is considered a wave packet with a localized form that can maintain its form and velocity through a nonlinear transferring media when the nonlinearity and the dispersion of that soliton are precisely exact balanced [1]. Solitons are important experimentally, in either practical or theoretical studies, due to their robust properties and unique nature, which appear in their dynamic propagation and interactions [2,3]. Numerous scientific studies have been carried out on soliton formation due to their essential dynamical processes [4,5,6,7,8]. These studies display the main concepts of solitons and many of their aspects, especially in optical field, including their models, dynamic features, and nonlinear effects. Many researchers have studied optical solitons to utilize their unique characteristics in technology development. These studies have focused on the contribution of nonlinear solitons in photonics research. In [9,10,11], the phase imprinting method was used with Bose-Einstein condensates (BECs) to generate dark solitons. In [12], Feshbach resonance was used to produce bright matter-wave solitons, converting interactions from repulsion to attraction. Beyond dark and bright solitons, other types of solitons have been studied and generated, such as vortex solitons and surface gap solitons in [13]. The soliton waves are governed through the BECs by the used Gross-Pitaevskii (GP) model, which comes from mean field theory [14]. In [15], an investigation was implemented to obtain bidirectional wave solutions for a newly created two-mode extension of the coupled Korteweg-De Vries (KdV)-Schrödinger equations. Theoretically, aside from the second-order diffraction/dispersion and cubic nonlinearity factors, the nonlinear Schrödinger equation (NLSE) including the parity-time (PT)-symmetric potential mostly describes the dynamics of PT solitons [16]. Many researchers have applied novel and new methods in order to be able to obtain newly created soliton solutions and many other exact solutions, such as the improved modified extended tanh function method [17], the modified extended mapping method [7,18], the extended Fan sub-equation method [19], the F-expansion method [20], and many others [21,22,23].

    In this work, the modified extended direct algebraic method (MEDAM) is applied to a higher-order nonlinearity Schrödinger equation that was used as GP equation in [15,24] but with constant coefficients, as in the following Eq (1.1):

    iΦt+b2Φxx+ϑ3|Φ|2Φ+ϑ5|Φ|4Φ+ϑ7|Φ|6Φ+ϑ9|Φ|8Φ+iνgΦx+VΦ=0, (1.1)

    where the function Φ(x,t) denotes the macroscopic wave, the two variables x relate to space, and t relates to time. b represents the velocity dispersion coefficient constant, and ϑ3, ϑ7, and ϑ9 are the respective real nonlinear constant coefficients. νg represents the group delay coefficient constant, while V denotes the complex PT-symmetric potential constant coefficient. The proposed combination between the used model and the applied MEDAM is discussed for the first time in this study, and it offers a variety of options compared to prior approaches, making it unique. Thus, this study solves the highly dispersive solitons, making them very useful for obtaining different soliton waves with different features. It is a noteworthy addition to the field as it discovers creative traveling wave solutions for the proposed model. The model, which is regarded as a GP equation with constant coefficients, describes the propagation of waves in a higher-order nonlinear BEC system. Compared to earlier research, which mainly concentrated on lower-order nonlinearities, the current work significantly advances the nonlinear Schrödinger equation by including cubic, quintic, septic, and nonic nonlinearities. The authors present new stability regimes and dynamic behaviors not before documented by contrasting their unique solutions with those discovered in the literature. This study differs from previous ones in that it presents novel analysis for managing the complexity of higher-order terms. In addition, we investigated the modulation instability for the suggested model, which provides certain conditions to improve the stability. The discoveries have potential implications in areas where controlling intricate nonlinear interactions is critical, such as optical fiber communications and Bose-Einstein condensates. Nowadays, the trend lies in studying stochastic soliton solutions for stochastic partial differential equations in order to investigate the dynamical behavior of these solitons under the effects of a Wiener process or what is called a white noise effect [25,26,27,28]. As far as we know, no prior studies have tried to solve Eq (1.1) with MEDAM; hence, the solutions obtained in this investigation are novel and have not been documented in any other literature. The following results could undoubtedly help achieve stable localized nonlinear wave transmission in extremely nonlinear media through experimentation. Different types of solutions including dark, bright, and singular solitons; mixed dark-bright solitons; singular periodic wave solutions; and Jacobi elliptic function (JEF) solutions are obtained. 3D illustrations are displayed for clarifying some of the results with their characteristics.

    The derivation of the model's solutions is discussed after a quick overview of the model and its basic mathematical background. The key points of the recommended approach are provided in Section 2. All of the results are displayed in Section 3, which also provides an explanation of the solution's many dynamic wave forms. Section 4 elaborates on the analysis of modulation instability. Graphical 2D, 3D, and contour illustrations of different soliton solutions are offered in Section 5. Some conclusions from the obtained results are given in Section 6.

    This section provides a quick description of the MEDAM's main steps [29,30,31,32].

    Consider the nonlinear partial differential equation (NLPDE) in its most common form, which is displayed below:

    ψ(Φ,Φt,Φx,Φxx,Φxt,Φxxt,...)=0. (2.1)

    Regarding the 2D of space and time coordinates, ψ represents a polynomial function in terms of Φ and some Φ derivatives.

    Procedure (Ⅰ): For obtaining the localized soliton solutions from solving Eq (2.1), we can implement and assume the wave function in Eq (2.2) taking the following form, as in [15]:

    Φ(x,t)=H(ξ) ei(βx+ct+α),   ξ=x+κt,    κ0, (2.2)

    where H acts as the amplitude values of solitons. A soliton's velocity is κ. The soliton wave number is β, the phase constant provides α, and the soliton frequency provides c, where β, α, c, and κ are constants with real values. This work looks to assess them through its processes. The non-linear ordinary differential equation (ODE) is subsequently obtained through placing Eq (2.2) into Eq (2.1):

    P(H,H,H,H)=0. (2.3)

    Procedure (Ⅱ): The following form is supposed as a solution for Eq (2.3):

    H(ξ)=Mn=MγnRn(ξ), (2.4)

    where γn (n=0,1,2,....,M) are parameters with real values that need estimation under the constraint that γM or γM must not be equal to zero simultaneously and R(ξ) follows the prerequisites of the auxiliary equation that is mentioned below (2.5):

    [R(ξ)]2=τ0+τ1R(ξ)+τ2R2(ξ)+τ3R3(ξ)+τ4R4(ξ)+τ6R6(ξ). (2.5)

    Procedure (Ⅲ): Applying the homogeneous balance principle to Eq (2.3) enables determining the values of M that will be used for Eq (2.4).

    Procedure (Ⅳ): By using Eq (2.4) with Eq (2.5) and then by substitution in Eq (2.3), a polynomial of R(ξ) terms is formed. Next, using the Wolfram Mathematica tool, we work out the resulting set of equations that represents a system of nonlinear algebraic equations (NLAEs). By collecting all terms of equal powers and setting their sum equal to zero, all the needed unknown values are obtained. After that, a number of exact solutions for Eq (2.1) may be generated.

    The traveling wave solution of Eq (1.1) is considered to have Eq (2.2) structure.

    With the transition mentioned above, Eq (1.1) will be converted into an ODE, by substituting Eq (2.2) into Eq (1.1). The resulting ODE can be separated into imaginary and real components as follows:

    Real component:

    bH+σH+2ϑ3H3+2ϑ5H5+2ϑ7H7+2ϑ9H9=0, (3.1)

    where σ=2(Vcβνg)bβ2.

    Imaginary component:

    (bβ+νg+κ)H=0, (3.2)

    and from Eq (3.2), we can deduce the wave number of the produced soliton as

    κ=(bβ+νg). (3.3)

    Assume that H=X14, and then Eq (3.1) can be described by

    4bXX3b(X)2+16σX2+32ϑ3X5/2+32ϑ5X3+32ϑ7X7/2+32ϑ9X4=0. (3.4)

    The general solution for Eq (3.4) might be written out using the recommended technique from Section 2 as follows:

    X(ξ)=γ0+γ1R(ξ)+γ1R(ξ), (3.5)

    where γ0,γ1, and γ1 are constants that will be computed according to the restriction that γ1 or γ10 in the same case.

    Plug Eqs (3.5) and (2.5) into Eq (3.4). With the aid of the auxiliary equation, collect coefficients of comparable powers, and equate them all to zero to construct a system of nonlinear algebraic equations. The Wolfram Mathematica program can then be used to solve this system, yielding the results that are shown below:

    Case (1): If τ0=τ1=τ3=τ6=0, the below set of solutions is obtained:

    ϑ3=ϑ5=ϑ7=0,γ1=γ0=0,γ1=145bτ42ϑ9,σ=bτ216.

    By using the acquired set of solutions, Eq (1.1) can be solved, giving its derived analytical solutions as follows:

    (1.1) If τ2>0, τ4<0, and b ϑ9>0, a bright soliton profile is raised as

    Φ1.1(x,t)=(145bτ22ϑ9sech[(x+κt)τ2])14ei(βx+ct+α). (3.6)

    (1.2) If τ2<0, τ4>0, and b ϑ9<0, a singular periodic solution is obtained as

    Φ1.2(x,t)=(145bτ22ϑ9sec[(x+κt)τ2])14ei(βx+ct+α). (3.7)

    (1.3) If τ2<0, τ4>0, and b ϑ9<0, a singular periodic solution is obtained as

    Φ1.3(x,t)=(145bτ22ϑ9csc[(x+κt)τ2])14ei(βx+ct+α). (3.8)

    Case (2): If τ1=τ3=τ6=0 and τ0=τ224τ4, the sets of solutions listed below result:

    (2.1) ϑ3=ϑ5=ϑ7=0,γ1=185bτ22ϑ9τ4,γ0=0,γ1=145bτ42ϑ9,σ=bτ28.

    (2.2) ϑ3=ϑ7=0,ϑ5=32bϑ9τ25,γ1=±τ285b2ϑ9τ4,γ0=185bτ2ϑ9,γ1=0,σ=bτ28.

    (2.3) ϑ3=ϑ7=0,ϑ5=32bϑ9τ25,γ1=0,γ0=185bτ2ϑ9,γ1=±145bτ42ϑ9,σ=bτ28.

    (2.4) ϑ3=ϑ7=0,ϑ5=3bϑ9τ25,γ1=τ285b2ϑ9τ4,γ0=145bτ2ϑ9,γ1=145bτ42ϑ9,σ=bτ22.

    The following structures are the outcome of certain exact solutions to Eq (1.1), as per the set of solutions (2.1):

    (2.1.1) If τ2>0, τ4>0, and b ϑ9<0, the below singular periodic solution is obtained:

    Φ2.1,1(x,t)=(185bϑ9(τ2tan[(x+κt)τ22]+cot[(x+κt)τ22]))14ei(βx+ct+α). (3.9)

    The following are some exact solutions to Eq (1.1) that follow the set of solutions (2.2):

    (2.2.1) If τ2<0, τ4>0, and b ϑ9<0, the below singular soliton solution is reached:

    Φ2.2,1(x,t)=(185bτ2ϑ9(1+coth[(x+κt)τ22]))14ei(βx+ct+α). (3.10)

    The following exact solution structures for Eq (1.1) are obtained based on the set of solutions (2.3) that was mentioned above:

    (2.3.1) If τ2<0, τ4>0, and b ϑ9<0, a dark soliton solution results as

    Φ2.3,1(x,t)=(185bτ2ϑ9(1±tanh[(x+κt)τ22]))14ei(βx+ct+α). (3.11)

    The following exact solution structures for Eq (1.1) are obtained based on the set of solutions (2.4):

    (2.4.1) If τ2<0, τ4>0, and b ϑ9<0, a singular soliton solution is obtained as follows:

    Φ2.4,1(x,t)=(145bτ2ϑ9(1+coth[(x+κt)2τ2]))14ei(βx+ct+α). (3.12)

    Case (3): If τ3=τ4=τ6=0, the following set of solutions is obtained:

    ϑ3=ϑ7=γ0=γ1=0, ϑ5=3τ14bϑ910τ0, γ1=145bτ02ϑ9,σ=bτ216.

    The following are the exact solutions for Eq (1.1) that may be raised by using this set of solutions:

    (3.1) If τ0>0,τ1=0,τ2>0, and b ϑ9<0, a singular soliton solution results as below:

    Φ3.1(x,t)=(145bτ22ϑ9csch[(x+κt)τ2])14ei(βx+ct+α). (3.13)

    (3.2) If τ0>0,τ1=0,τ2<0, and b ϑ9<0, a singular periodic solution is obtained as

    Φ3.2(x,t)=(145bτ22ϑ9csc[(x+κt)τ2])14ei(βx+ct+α). (3.14)

    (3.3) If τ0=τ214τ2,τ1<0,τ2>0, and b ϑ9<0, an exponential solution is obtained as

    Φ3.3(x,t)=(τ225bτ02ϑ9(1τ1+2τ2eτ2(κt+x)))14ei(βx+ct+α). (3.15)

    Case (4): If τ0=τ1=τ2=τ6=0, the following set of solutions results:

    γ1=γ0=ϑ3=ϑ7=σ=0, γ1=145bτ42ϑ9, ϑ5=3τ34bϑ910τ4.

    According to this obtained set of solutions, an exact solution to Eq (1.1) results in the following structure:

    (4.1) If τ3>0, τ4<0, and b ϑ9>0, the following rational wave solution is obtained:

    Φ4.1,1(x,t)=(5bτ42ϑ9(τ34τ4+(x+κt)2τ23))14ei(βx+ct+α). (3.16)

    Case (5): If τ0=τ1=τ6=0, the following set of solutions is obtained:

    γ1=γ0=ϑ3=ϑ7=0, γ1=145bτ42ϑ9, ϑ5=3τ34bϑ910τ4, σ=bτ216.

    By means of the generated set of solutions, Eq (1.1) can be solved, giving the following solutions:

    (5.1) If τ2>0,τ3=2τ2τ4, and b ϑ9<0, the dark soliton and singular solution solutions are obtained in the below structures:

    Φ5.1,1(x,t)=(185bτ22ϑ9(1+tanh[12(x+κt)τ2]))14ei(βx+ct+α), (3.17)

    or

    Φ5.1,2(x,t)=(185bτ22ϑ9(1+coth[12(x+κt)τ2]))14ei(βx+ct+α). (3.18)

    (5.2) If τ2>0,τ32τ2τ4,τ4>0, b ϑ9<0, and (τ32τ2τ4tanh[12(x+κt)τ2])>0, a mixed bright-dark soliton solution is obtained as

    Φ5.2(x,t)=(τ245bτ42ϑ9(sech2[12(x+κt)τ2]τ32τ2τ4tanh[12(x+κt)τ2]))14ei(βx+ct+α). (3.19)

    (5.3) If τ2<0,τ32τ2τ4,τ3<0,τ4>0, and b ϑ9<0, the following singular periodic solution is obtained:

    Φ5.3(x,t)=(τ245bτ42ϑ9(sec2[12(x+κt)τ2]τ32τ2τ4tan[12(x+κt)τ2]))14ei(βx+ct+α). (3.20)

    Case (6): If τ2=τ4=τ6=0, the set of solutions that follow is

    ϑ3=ϑ7=γ1=0, γ1=145bτ02ϑ9,γ0=τ165bτ0ϑ9, ϑ5=τ145bϑ92τ0, τ3=4τ3127τ20,σ=bτ2124τ0.

    Using the revealed set of solutions, Eq (1.1) gives its exact solution that can be expressed as follows:

    (6.1) If τ0<0, τ1<0, τ3>0, and b ϑ9>0, a Weierstrass elliptic doubly periodic solution results as

    Φ6.1(x,t)=(5bϑ9τ0(τ162τ08[(x+κt)τ34;(4τ1τ3,4τ0τ3)]])14ei(βx+ct+α). (3.21)

    Case (7): If τ0=τ1=τ3=0, the set of solutions below results:

    ϑ3=ϑ5=ϑ7=γ1=γ0=τ6=0, γ1=145bτ42ϑ9, σ=bτ216.

    Using the acquired solutions' set, the following solution for Eq (1.1) is arrived at:

    (7.1) If τ2>0 and b ϑ9 τ4<0, the below form of a singular soliton solution is obtained:

    Φ7.1(x,t)=(145bτ22ϑ9τ4csch[(x+κt)τ2])14ei(βx+ct+α). (3.22)

    (7.2) If τ2<0 and b ϑ9 τ4>0, a singular periodic solution is obtained:

    Φ7.2(x,t)=(145bτ22ϑ9τ4sec[(x+κt)τ2])14ei(βx+ct+α). (3.23)

    Case (8): If τ1=τ3=τ6=0 and γ0=0, we obtain the subsequent solutions' sets:

    (8.1) ϑ3=ϑ5=ϑ7=0,γ1=145bτ02ϑ9,γ1=τ285b2ϑ9τ0,τ4=τ224τ0,σ=bτ28.

    (8.2) γ1=τ4=ϑ3=ϑ5=ϑ7=0, γ1=145bτ02ϑ9, σ=bτ216.

    The set of solutions (8.1) indicates that Eq (1.1) has an exact solution, which can be phrased as follows:

    (8.1.1) If τ0=1,τ2=24m2, τ4=1, and b ϑ9<0, a reached singular periodic solution is

    Φ8.1,1(x,t)=(125b2ϑ9csc[2(x+κt)])14ei(βx+ct+α). (3.24)

    The set of solutions (8.2) indicates that Eq (1.1) has certain exact solutions, which are as follows:

    (8.2.1) If τ0=1,τ2=m21,τ4=m2 or τ0=14,τ2=m222,τ4=m44, and b ϑ9<0, two singular periodic solutions can be reached as below:

    Φ8.2,1(x,t)=(145b2ϑ9csc[x+κt])14ei(βx+ct+α), (3.25)

    or

    Φ8.2,2(x,t)=(145b2ϑ9sec[x+κt])14ei(βx+ct+α). (3.26)

    (8.2.2) If τ0=m2,τ2=2m21,τ4=1m2, and b ϑ9>0, a bright soliton solution is reached as

    Φ8.2,3(x,t)=(145b2ϑ9sech[x+κt])14ei(βx+ct+α). (3.27)

    (8.2.3) If τ0=m42m3+m2,τ2=4m,τ4=m2+6m1, b ϑ9<0, and m=322, we obtain the following JEF solution:

    Φ8.2,4(x,t)=(145b(m1)22ϑ9(1+m sn2(x+κt)cn(x+κt)dn(x+κt)))14ei(βx+ct+α). (3.28)

    In this section, we examine the modulation instability (MI) of Eq (1.1), through standard linear stability analysis of Eq (1.1), as shown below [33]:

    Φ=(y(x,t)+R) ei(R x). (4.1)

    In this scenario, R represents the steady-state solution for Eq (1.1). To investigate the development of the disturbance y(x,t), linear stability analysis is utilized. The process involves inserting Eq (3.19) into Eq (1.1) and then linearizing the resulting expression as follows:

    b yxx+2i yt+2i(b R yx+νg) yx+(b R22R νg+2 V+18R4 ϑ9+14R3 ϑ7+10R2 ϑ5+6R ϑ3) (y+y)=0, (4.2)

    where represents the conjugate of the complex function. Let us assume that the solution of Eq (4.2) can be expressed in the manner detailed below:

    y=F1 ei(L xω t)+F2 ei(L xω t), (4.3)

    where ω represents the frequency of perturbation, and L denotes the normal wave number. By substituting Eq (4.3) into Eq (4.2), we derive the dispersion relationship in the subsequent manner:

    bL2[b(L22R2)4V]+8bωLR4ω2+4L[νg(2ωbLR)bLR(R(R(9Rϑ9+7ϑ7)+5ϑ5)+3ϑ3)Lν2g]=0. (4.4)

    Now, when we assess the dispersion relation ω derived from Eq (4.4), the result we obtain is as follows:

    ω=bL R+L νg±12b L2(b L2+2bR24R(νg+9R3 ϑ9+7R2 ϑ7+5R ϑ5+3ϑ3)4V). (4.5)

    Now, we can explore the steady-state stability using the dispersion relation above. If the wave number is real, the steady state transitions to stability against minor disturbances. However, if the wave number is imaginary, the steady-state solution becomes unstable, as the perturbation grows exponentially. Hence, a steady-state solution becomes unstable if b [b L2+2 b R24 R(νg+9 R3 ϑ9+7 R2 ϑ7+5R ϑ5+3 ϑ3)4 V]<0. To obtain the MI gain spectrum G(R), we proceed in the following manner:

    G(R)=2Im[bL R+Lνg±12bL2(bL2+2bR24R(νg+9R3ϑ9+7R2ϑ7+5Rϑ5+3ϑ3)4V)]. (4.6)

    Different sets of solutions were generated by setting specific values for the parameters for Eq (1.1). As a result, this method has led to a number of unique outcomes that have never been reported before. Two- and three-dimension visual graphical aids and contour representations of several specific solutions, are provided to display the recovered solutions and clarify their physical properties.

    The bright soliton solution for Eq (3.6) is shown in Figure 1 with parameters b=0.7, ϑ9=0.5, νg=0.7, β=0.6, c=1.1, α=0.7, τ2=0.8,  1t1,  and 20x20. The centered maximum intensity peak is the predominant feature of the bright soliton. There are several impressive uses for a bright soliton solution in nonlinear fiber optics and optical telecommunications systems. Figure 2 clarifies the singular periodic wave solution for Eq (3.7) with parameters b=0.8, ϑ9=0.6, νg=0.75, β=0.7, c=0.9, α=0.8, τ2=0.8,  5t5,  and 15x15. This type of solution has surprising uses in sophisticated signal processing methods as well as nonlinear optics [34]. Figure 3 displays the dark soliton solution for Eq (3.11) when the parameters are b=0.9, ϑ9=0.4, νg=0.6, β=0.7, c=0.75, α=0.6, τ2=0.9,  0t5,  and 20x20. The centered lowest amplitude is the primary feature of the dark soliton. Because of its stability and resistance to loss, the dark soliton is a great tool for fiber lasers. It was found that a dark soliton is more suitable than a bright soliton for optical communications applications [35]. In addition, the singular soliton solution in Eq (3.13) is shown through Figure 4 according to the applied values of the input parameters as b=0.9, ϑ9=0.5, νg=0.5, β=0.6, c=0.8, α=0.7, τ2=0.9,  0t5,  and 15x15. The singular soliton's centered unbounded amplitude is its primary feature. Such a type of solution is employed to generate useful models by comprehending the severe wave event. Figure 5 presents the spectrum of MI of Eq (4.6) with different values of a and b.

    Figure 1.  Graphical simulation of the bright soliton profile of Eq (3.6).
    Figure 2.  Graphical simulation of the singular periodic profile of Eq (3.7).
    Figure 3.  Graphical simulation of the dark soliton profile of Eq (3.11).
    Figure 4.  Graphical simulation of the singular soliton profile of Eq (3.13).
    Figure 5.  Gain spectrum of MI of Eq (4.6).

    A wide range of highly dispersive soliton solutions of NLPDEs has been discovered in this work. For the suggested model, several exact solutions were found, including periodic, singular periodic, and JEF wave solutions and singular, bright, and dark solitons. This retrieval was made possible using MEDAM techniques, with which we found several solutions that satisfied the governing equation. Additionally, for the physical representation, specific behaviors of the obtained solutions have been visually presented by some graphs. When compared to earlier study findings, the results show the novelty and innovation of the proposed technique, which is accurate, simple to use, and beneficial for optical solitons as well as other mathematical physics applications. All of the outcomes show that a broad spectrum of different solutions for solitary waves that propagate with their shape can be obtained by striking a balance between the effects of dispersion and non-linearity. We also explored the MI analysis of the given presented model.

    Abeer Khalifa: Formal analysis, Software; Hamdy Ahmed: Validation, Methodology; Niveen Badra: Investigation, Writing–review & editing; Wafaa Rabie: Software, Writing–review & editing; Farah M. Al-Askar: Formal analysis, Writing–review & editing; Wael Mohammed: Resources, Writing–review & editing. All authors have read and agreed to the published version of the manuscript.

    The work was supported by Princess Nourah bint Abdulrahman University Researcher Supporting Project number (PNURSP2024R 273), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia.

    The authors declare that they have no conflicts of interest.



    [1] R. P. Agarwal, Difference equations and inequalities: Theory, methods and applications, Marcel Dekker, 1992.
    [2] Y. H. Long, L. Wang, Global dynamics of a delayed two-patch discrete SIR disease model, Commu. Nonlinear Sci., 83 (2020), 105117. doi: 10.1016/j.cnsns.2019.105117
    [3] Y. T. Shi, J. S. Yu, Wolbachia infection enhancing and decaying domains in mosquito population based on discrete models, J. Biol. Dynam., 14 (2020), 679–695. doi: 10.1080/17513758.2020.1805035
    [4] F. M. Atici, G. S. Guseinov, Positive periodic solutions for nonlinear difference equations with periodic coefficients, J. Math. Anal. Appl., 232 (1999), 166–182. doi: 10.1006/jmaa.1998.6257
    [5] Z. M. Guo, J. S. Yu, The existence of periodic and subharmonic solutions of sub-quadratic second order difference equations, J. Lond. Math. Soc., 68 (2003), 419–430. doi: 10.1112/S0024610703004563
    [6] J. Hendersona, R. Luca, Existence of positive solutions for a system of semipositone coupled discrete boundary value problems, J. Differ. Equ. Appl., 25 (2019), 516–541. doi: 10.1080/10236198.2019.1585831
    [7] T. S. He, Y. W. Zhou, Y. T. Xu, C. Y. Chen, Sign-changing solutions for discrete second-order periodic boundary value problems, B. Malays. Math. Sci. So., 38 (2015), 181–195. doi: 10.1007/s40840-014-0012-1
    [8] Y. H. Long, B. L. Zeng, Multiple and sign-changing solutions for discrete Robin boundary value problem with parameter dependence, Open Math., 15 (2017), 1549–1557. doi: 10.1515/math-2017-0129
    [9] Y. H. Long, J. L. Chen, Existence of multiple solutions to second-order discrete Neumann boundary value problems, Appl. Math. Lett., 83 (2018), 7–14. doi: 10.1016/j.aml.2018.03.006
    [10] Y. H. Long, S. H. Wang, Multiple solutions for nonlinear functional difference equations by the invariant sets of descending flow, J. Differ. Equ. Appl., 25 (2019), 1768–1789. doi: 10.1080/10236198.2019.1694014
    [11] Y. H. Long, Existence of multiple and sign-changing solutions for a second-order nonlinear functional difference equation with periodic coefficients, J. Differ. Equ. Appl., 26 (2020), 966–986. doi: 10.1080/10236198.2020.1804557
    [12] M. J. Ma, Z. M. Guo, Homoclinic orbits for second order self-adjoint difference equations, J. Math. Anal. Appl., 232 (2006), 513–521.
    [13] Y. H. Long, Homoclinic orbits for a class of noncoercive discrete hamiltonian systems, J. Appl. Math., 2012 (2012), 1–21.
    [14] Z. Zhou, D. F. Ma, Multiplicity results of breathers for the discrete nonlinear Schrödinger equations with unbounded potentials, Sci. China Math., 58 (2015), 781–790. doi: 10.1007/s11425-014-4883-2
    [15] Y. H. Long, Y. B. Zhang, H. P. Shi, Homoclinic solutions of 2nth-order difference equations containing both advance and retardation, Open Math., 14 (2016), 520–530. doi: 10.1515/math-2016-0046
    [16] L. Erbe, B. G. Jia, Q. Q. Zhang, Homoclinic solutions of discrete nonlinear systems via variational method, J. Appl. Anal. Compt., 9 (2019), 271–294.
    [17] Q. Q. Zhang, Homoclinic orbits for discrete Hamiltonian systems with local super-quadratic conditions, Commun. Pur. Appl. Anal., 18 (2019), 425–434. doi: 10.3934/cpaa.2019021
    [18] J. H. Kuang, Z. M. Guo, Heteroclinic solutions for a class of p-Laplacian difference equations with a parameter, Appl. Math. Lett., 100 (2020), 106034. doi: 10.1016/j.aml.2019.106034
    [19] Y. H. Long, S. H. Wang, J. L. Chen, Multiple solutions of fourth-order difference equations with different boundary conditions, Bound. Value Probl., 2019 (2019), 1–25. doi: 10.1186/s13661-018-1115-7
    [20] S. H. Wang, Y. H. Long, Multiple solutions of fourth-order functional difference equation with periodic boundary conditions, Appl. Math. Lett., 104 (2020), 106292. doi: 10.1016/j.aml.2020.106292
    [21] G. H. Lin, Z. Zhou, Homoclinic solutions of discrete ϕ-Laplacian equations with mixed nonlinearities, Commun. Pur. Appl. Anal., 17 (2018), 1723–1747. doi: 10.3934/cpaa.2018082
    [22] H. Fang, D. P. Zhao, Existence of nontrivial homoclinic orbits for fourth-order difference equations, Appl. Math. Compt., 214 (2009), 163–170. doi: 10.1016/j.amc.2009.03.061
    [23] N. Dimitrov, S. Tersian, Existence of homoclinic solutions for a nonlinear fourth order p-Laplacian difference equation, Discrete Cont. Dyn. B, 25 (2020), 555–567.
    [24] I. Ekeland, Convexity methods in Hamiltonian mechanics, Springer-Verlag Berlin Heidelberg, 1990.
    [25] S. Tersian, J. Chaparova, Periodic and homoclinic solutions of extended Fisher-Kolmogorov equations, J. Math. Anal. Appl., 260 (2001), 490–506. doi: 10.1006/jmaa.2001.7470
    [26] Y. C. Zhou, H. Cao, Y. N. Cao, Difference equations and their applications, Beijing, Science Press, 2014.
  • Reader Comments
  • © 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2349) PDF downloads(124) Cited by(6)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog