Cytomegalovirus dynamics model with random behavior

Dong-Mei Li; Bing Chai; Yu-Li Fu; Qi Wang; Dong-Mei Li; Bing Chai; Yu-Li Fu; Qi Wang

doi:10.3934/math.2020410

AIMS Mathematics

2020, Volume 5, Issue 6: 6373-6394. doi: 10.3934/math.2020410

Previous Article Next Article

Research article Special Issues

Cytomegalovirus dynamics model with random behavior

1.
College of science, Harbin University of Science and Technology, 52 xuefu road, Harbin 150080, China
2.
Heilongjiang Hospital, Harbin 150036, China

Received: 13 May 2020 Accepted: 13 July 2020 Published: 10 August 2020
MSC : 60H10, 92B05, 92D30

In view of the particularity of cytomegalovirus infection in infants and considering the uncertainty of infection mode and treatment, a dynamic model of cytomegalovirus with random behavior is established in this paper. The existence and uniqueness of the solution of the model are proved. Sufficient conditions for the existence of asymptotic, ergodic and extinctive solutions are provided. By using numerical simulation, the influence of uncertainty in breast milk handling and treatment on the variation of cytomegalovirus (CMV) are analyzed, which provides theoretical support for the strategy of preventing infant infection and the basis treatment.

Keywords:

Citation: Dong-Mei Li, Bing Chai, Yu-Li Fu, Qi Wang. Cytomegalovirus dynamics model with random behavior[J]. AIMS Mathematics, 2020, 5(6): 6373-6394. doi: 10.3934/math.2020410

Related Papers:

[1]	Naher Mohammed A. Alsafri, Hamad Zogan . Probing the diversity of kink solitons in nonlinear generalised Zakharov-Kuznetsov-Benjamin-Bona-Mahony dynamical model. AIMS Mathematics, 2024, 9(12): 34886-34905. doi: 10.3934/math.20241661
[2]	Muhammad Bilal Riaz, Faiza Naseer, Muhammad Abbas, Magda Abd El-Rahman, Tahir Nazir, Choon Kit Chan . Solitary wave solutions of Sawada-Kotera equation using two efficient analytical methods. AIMS Mathematics, 2023, 8(12): 31268-31292. doi: 10.3934/math.20231601
[3]	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
[4]	Muhammad Bilal, Javed Iqbal, Ikram Ullah, Aditi Sharma, Hasim Khan, Sunil Kumar Sharma . Novel optical soliton solutions for the generalized integrable (2+1)- dimensional nonlinear Schrödinger system with conformable derivative. AIMS Mathematics, 2025, 10(5): 10943-10975. doi: 10.3934/math.2025497
[5]	Naher Mohammed A. Alsafri . Solitonic behaviors in the coupled Drinfeld-Sokolov-Wilson system with fractional dynamics. AIMS Mathematics, 2025, 10(3): 4747-4774. doi: 10.3934/math.2025218
[6]	Naveed Iqbal, Meshari Alesemi . Soliton dynamics in the $(2+1)$ -dimensional Nizhnik-Novikov-Veselov system via the Riccati modified extended simple equation method. AIMS Mathematics, 2025, 10(2): 3306-3333. doi: 10.3934/math.2025154
[7]	Jamshad Ahmad, Zulaikha Mustafa, Mehjabeen Anwar, Marouan Kouki, Nehad Ali Shah . Exploring solitonic wave dynamics in the context of nonlinear conformable Kairat-X equation via unified method. AIMS Mathematics, 2025, 10(5): 10898-10916. doi: 10.3934/math.2025495
[8]	Abeer S. Khalifa, Hamdy M. Ahmed, Niveen M. Badra, Wafaa B. Rabie, Farah M. Al-Askar, Wael W. Mohammed . New soliton wave structure and modulation instability analysis for nonlinear Schrödinger equation with cubic, quintic, septic, and nonic nonlinearities. AIMS Mathematics, 2024, 9(9): 26166-26181. doi: 10.3934/math.20241278
[9]	Saima Noor, Azzh Saad Alshehry, Asfandyar Khan, Imran Khan . Analysis of soliton phenomena in (2+1)-dimensional Nizhnik-Novikov-Veselov model via a modified analytical technique. AIMS Mathematics, 2023, 8(11): 28120-28142. doi: 10.3934/math.20231439
[10]	Amna Mumtaz, Muhammad Shakeel, Abdul Manan, Marouan Kouki, Nehad Ali Shah . Bifurcation and chaos analysis of the Kadomtsev Petviashvili-modified equal width equation using a novel analytical method: describing ocean waves. AIMS Mathematics, 2025, 10(4): 9516-9538. doi: 10.3934/math.2025439

Abstract

1. Introduction

The intriguing quality of nonlinearity in nature has led many scientists to view nonlinear science as the most significant area of research for a fundamental understanding of nature. The study of a variety of classes of NODEs and nonlinear partial differential equations (NPDEs) is crucial to the mathematical modelling of complex processes that vary over time. These models come from a wide range of disciplines, including optical fibres, economics, solid state physics, plasma physics, infectious disease epidemiology, elasticity, neural networks, and population ecology. Due to such applications, the examination of soliton solutions of the before stated phenomena has been an exciting and extremely dynamic topic of research for the past several decades, with the related issue of the construction of closed-form wave solutions to a broad class of nonlinear differential equations, particularlyNPDEs. For such situations, closed-form soliton solutions offer improved internal information. Thus, a great deal of work has been done by mathematicians and physical scientists to obtain closed-form soliton solutions of such NPDEs. A number of effective and powerful methods have been developed, including the first integral method ^[1], Bäcklund transformation method ^[2], the modified simple equation method ^[3], the $(G' /G)$ -expansion method ^[4,5,6], the exp-function method ^[7], the sine-cosine method ^[8], the homogeneous balance method ^[9], the modified Kudryashov method ^[10], the F-expansion method ^[11], the bilinear method ^[12], the tanh-function method ^[13], the projective Riccit equation method ^[14], and the $(G' /G, 1/G)$ -expansion method ^[15] Riccati-Bernoulli Sub-ODE ^[16], Poincaré-Lighthill-Kuo method ^[17], sub-equation method ^[18], Khater method ^[19], Sardar sub-equation method ^[20], EDAM ^{[21,22,23,24,25]}, Hirota bilinear method ^[26], planar dynamic system method ^[27] and Bifurcation analysis method ^[28], among others ^[29,30]. The computation of complex standard eigenvalue problems for laminar-turbulent transition prediction, while the investigates dynamically controlling terahertz wavefronts. In the third reference, structural reliability analysis is performed using an improved wolf pack algorithm, while the fourth reference introduces a class of digital integrators based on trigonometric quadrature rules ^[31,32,33]. Lastly, the fifth reference focuses on opinion formation analysis in group decision-making systems. These references collectively showcase the breadth and depth of research across various domains, highlighting the importance of interdisciplinary approaches in addressing complex scientific and engineering challenges ^[34,35].

One of the most significant, straightforward, and efficient algebraic techniques in the aforementioned analytical methods for obtaining soliton solutions to NPDEs is the m-EDAM. The strategic m-EDAM comprises converting CIE into a NODE by variable transformation, which is then reduced to a set of nonlinear algebraic equations with the assumption of a closed-form solution ^[36,37,38]. Many soliton solutions in terms of exponential, rational, trigonometric, and hyperbolic functions are produced when the resultant system is solved using the Maple tool. A soliton is a single, self-reinforcing wave packet that moves across a medium without altering its shape or velocity. Since soliton solutions to nonlinear fractional partial differential equations (NFPDEs) provide a higher level of detail and coverage than traditional solutions, their study is still important from an academic standpoint ^[39,40,41]. Their resilience and fermion stability make them valuable in many areas of engineering and science. They make it possible for nonlinear systems to maintain long-distance coherence and convey information efficiently ^[42,43].

Numerous researchers have used the proposed m-EDAM to construct a new plethora of soliton solutions for many NPDEs andNFPDEs. For example, Saima et al. in ^[44] constructed soliton solutions for the (2+1)-dimensional Nizhnik-Novikov-Veselov model (NNVM) using a novel variant of m-EDAM called $r$ +mEDAM. Ikram et al. employed the suggested mEDAM to analyse front events and solitons in the fractional Kolmogorov-Petrovskii-Piskunov equation ^[45]. Bilal et al. explored the soliton phenomena in the nonlinear spatiotemporal fractional Higgs system analytically, using an improved version of the ground-breaking analytical technique m-EDAM ^[46]. Lastly, Rashid et al. investigated soliton solutions in fractional quantum mechanical equations of nonlinear systems using the effective m-EDAM ^[47]. Overall, the m-EDAM research revealed that it is a reliable, practical, and easy-to-use tool that offers more general wave solutions for nonlinear models in mathematical physics and engineering.

Motivated by the ongoing research on soliton solutions, this research aims to extract and analyse soliton solutions for (2+1)-dimensional CIE using the m-EDAM. This model is articulated as ^[48]:

$\begin{equation} z_{tx}+(\lambda z^3-z_{xx}-\lambda z)_x+\sigma z_{yy} = 0, \end{equation}$

(1.1)

where $z = z(x, y, t)$ , $\sigma$ is the degradation coefficient, while $\lambda$ is the diffusion coefficient. In a physical environment, the diffusion of a gas in a homogeneous medium is a significant phenomenon, and the CIE offers a helpful model to investigate such phenomena. One famous response to the Duffing equation in the physical sciences is the (2+1)-dimensional CIE ^[49].

Numerous scholars have addressed the CIE by employing various analytical approaches. Using the Hirota bilinear approach, for example, Sualiman et al. produced lump soliton solutions for CIE with variable coefficients ^[50]. To derive the kink and cuspon soliton solutions for CIE, Akbar et al. used the first integral equation in ^[48]. Similarly, a number of new generalised solitonary solutions for the CIE with unique physical structures have been constructed by Sakthivel and Chun ^[49] using the exp-function method. Lastly, Demiray and Bayrakci created soliton solutions for CIE in ^[51] by applying the sine-Gordon expansion approach. However, in this study, we aim to address CIE to construct soliton solutions for it using m-EDAM. Our study reveals that the soliton solutions obtained for CIE are either cuspon and kink soliton solutions, such as dark kink, bright kink, and bright-dark kink. A kink soliton is a kind of solitary wave that is distinguished by an abrupt change or discontinuity in the field or variable it represents. It is also referred to as a topological soliton or domain wall. Kink solitons usually occur in systems that experience a phase transition of a scalar field. Conversely, a cuspon soliton, also known as a cusped soliton, is a different kind of solitary wave that has a cusp, or sharp peak, at the crest. It is a solution to a specific class of nonlinear wave equations in which the propagating solitary wave keeps its shape.

The format of this article is as follows: Section 1 provides an introduction. Section 2 presents the operational methodology of m-EDAM. In Section 3, we build some new families of soliton solutions for CIE, and in Section 4, several illustrations and a graphical description are shown and analysed. Finally, Section 5 provides the conclusion.

2. The operational procedure of m-EDAM

We describe the operating methods of the m-EDAM in this portion of the paper, focusing on addressing the subsequent general NPDE:

$\begin{equation} R(z, z_t, z_{y_1}, z_{y_2},z z_{y_1}, \ldots) = 0, \end{equation}$

(2.1)

where $z = z(t, y_1, y_2, y_3, \ldots, y_j)$ .

The strategy taken to solve Eq (2.1) is given as follows:

(1). Commencing with a variable transformation of the framework $z(t, y_1, y_2, y_3, \dots, y_j) = Z(\vartheta)$ , where $\vartheta$ can be stated in a number of ways, Eq (2.1) proceeds through a transformation, resulting in the post-NODE:

$\begin{equation} Q(Z,Z'Z,Z',\dots) = 0, \end{equation}$

(2.2)

where $Z' = \frac{d Z}{d\vartheta}$ . On some occasions, the NODE becomes vulnerable to the homogeneous balancing principle due to the integration of Eq (2.2).

(2). Next, we suggest a series-based solution to the NODE in Eq (2.1) using the Riccati ODE, which is as follows:

$\begin{equation} Z(\vartheta) = \sum\limits_{l = -N}^{N}\rho_l(\zeta(\vartheta))^l. \end{equation}$

(2.3)

In this context, $\rho_l (l = -N, ..., N)$ denotes the unknown parameters, and $\zeta(\vartheta)$ is the solution to the first-order Riccati ODE stated as:

$\begin{equation} \zeta'(\vartheta) = R(\zeta(\vartheta))^2+Q \zeta(\vartheta)+P, \end{equation}$

(2.4)

where $Q, R$ , and $P$ are invariables.

(3). Applying a homogeneous balancing strategy to Eq (2.2) that involves balancing the greatest nonlinear component and the highest-order derivative will yield the positive integer $N$ utilised in (2.3).

(4). Following that, we place (2.3) into (2.2) or the equation that arises from integrating (2.2), organising terms $\zeta(\vartheta)$ with identical orders. An expression in $\zeta(\vartheta)$ is produced by this process. A set of algebraic equations representing the variables $\rho_l$ (where $l = -N, \ldots, N$ ) and other pertinent parameters is then expressed after the coefficients in this expression are equated to zero.

(5). This system of algebraic equations is then solved using Maple tool.

(6). Then, by calculating and entering the unknown values into Eq (2.3) together with the $\zeta(\vartheta)$ (the solution to (2.4)), analytical soliton solutions for (2.2) are obtained. We can get the resulting families of soliton solutions by using the generic solution of Eq (2.4):

Family 1. For $\Phi < 0$ & $R\neq0$ ,

$\begin{equation*} \zeta_1(\vartheta) = -{\frac {Q}{2R}}+{\frac {\sqrt {-\Phi}\tan \left( \frac{1}{2}\,\sqrt {-\Phi }\vartheta \right) }{2R}}, \end{equation*}$

$\begin{equation*} \zeta_2(\vartheta) = -{\frac {Q}{2R}}-{\frac {\sqrt {-\Phi}\cot \left( \frac{1}{2}\,\sqrt {-\Phi }\vartheta \right) }{2R}}, \end{equation*}$

$\begin{equation*} \zeta_3(\vartheta) = -{\frac {Q}{2R}}+{\frac {\sqrt {-\Phi} \left( \tan \left( \sqrt {-\Phi}\vartheta \right) + \left( \sec \left( \sqrt {-\Phi}\vartheta \right) \right) \right) }{2R}}, \end{equation*}$

$\begin{equation*} \zeta_4(\vartheta) = -{\frac {Q}{2R}}-{\frac {\sqrt {-\Phi} \left( \cot \left( \sqrt {-\Phi }\vartheta \right) + \left( \csc \left( \sqrt {-\Phi}\vartheta \right) \right) \right) }{2R}}, \end{equation*}$

and

$\begin{equation*} \zeta_5(\vartheta) = -{\frac {Q}{2R}}+{\frac {\sqrt {-\Phi} \left( \tan \left( \frac{1}{4}\, \sqrt {-\Phi}\vartheta \right) -\cot \left( \frac{1}{4}\,\sqrt {-\Phi}\vartheta \right) \right) }{4R}}. \end{equation*}$

Family 2. For $\Phi > 0$ & $R\neq0$ ,

$\begin{equation*} \zeta_6(\vartheta) = -{\frac {Q}{2R}}-{\frac {\sqrt {\Phi}\tanh \left( \frac{1}{2}\,\sqrt {\Phi} \vartheta \right) }{2R}}, \end{equation*}$

$\begin{equation*} \zeta_7(\vartheta) = -{\frac {Q}{2R}}-{\frac {\sqrt {\Phi}\coth \left( \frac{1}{2}\,\sqrt {\Phi} \vartheta \right) }{2R}}, \end{equation*}$

$\begin{equation*} \zeta_8(\vartheta) = -{\frac {Q}{2R}}-{\frac {\sqrt {Z} \left( \tanh \left( \sqrt {\Phi}\vartheta \right) +i \left( {\it sech} \left( \sqrt {\Phi}\vartheta \right) \right) \right) }{2R}}, \end{equation*}$

$\begin{equation*} \zeta_9(\vartheta) = -{\frac {Q}{2R}}-{\frac {\sqrt {\Phi} \left( \coth \left( \sqrt {\Phi}\vartheta \right) + \left( {\it csch} \left( \sqrt {\Phi}\vartheta \right) \right) \right) }{2R}}, \end{equation*}$

and

$\begin{equation*} \zeta_{10}(\vartheta) = -{\frac {Q}{2R}}-{\frac {\sqrt {\Phi} \left( \tanh \left( \frac{1}{4}\, \sqrt {\Phi}\vartheta \right) -\coth \left( \frac{1}{4}\,\sqrt {\Phi}\vartheta \right) \right) }{4R}}. \end{equation*}$

Family 3. For $QP > 0$ & $Q = 0$ ,

$\begin{equation*} \zeta_{11}(\vartheta) = \sqrt {{\frac {P}{R}}}\tan \left( \sqrt {RP}\vartheta \right), \end{equation*}$

$\begin{equation*} \zeta_{12}(\vartheta) = -\sqrt {{\frac {P}{R}}}\cot \left( \sqrt {RP}\vartheta \right), \end{equation*}$

$\begin{equation*} \zeta_{13}(\vartheta) = \sqrt {{\frac {P}{R}}} \left( \tan \left( 2\,\sqrt {RP}\vartheta \right) + \left( \sec \left( 2\,\sqrt {RP}\vartheta \right) \right) \right), \end{equation*}$

$\begin{equation*} \zeta_{14}(\vartheta) = -\sqrt {{\frac {P}{R}}} \left( \cot \left( 2\,\sqrt {RP}\vartheta \right) + \left( \csc \left( 2\,\sqrt {RP}\vartheta \right) \right) \right), \end{equation*}$

and

$\begin{equation*} \zeta_{15}(\vartheta) = \frac{1}{2}\sqrt {{\frac {P}{R}}} \left( \tan \left( \frac{1}{2}\,\sqrt {RP}\vartheta \right) -\cot \left( \frac{1}{2}\,\sqrt {RP}\vartheta \right) \right). \end{equation*}$

Family 4. For $RP < 0$ and $Q = 0$ ,

$\begin{equation*} \zeta_{16}(\vartheta) = -\sqrt {-{\frac {P}{R}}}\tanh \left( \sqrt {-RP}\vartheta \right), \end{equation*}$

$\begin{equation*} \zeta_{17}(\vartheta) = -\sqrt {-{\frac {P}{R}}}\coth \left( \sqrt {-RP}\vartheta \right), \end{equation*}$

$\begin{equation*} \zeta_{18}(\vartheta) = -\sqrt {-{\frac {P}{R}}} \left( \tanh \left( 2\,\sqrt {-RP}\vartheta \right)+ \left( i{\it sech} \left( 2\, \sqrt {-RP}\vartheta \right) \right) \right), \end{equation*}$

$\begin{equation*} \zeta_{19}(\vartheta) = -\sqrt {-{\frac {P}{R}}} \left( \coth \left( 2\,\sqrt {-RP}\vartheta \right) + \left( {\it csch} \left( 2\, \sqrt {-RP}\vartheta \right) \right) \right), \end{equation*}$

and

$\begin{equation*} \zeta_{20}(\vartheta) = -\frac{1}{2}\sqrt {-{\frac {P}{R}}} \left( \tanh \left( \frac{1}{2}\,\sqrt {-RP}\vartheta \right) +\coth \left( \frac{1}{2}\,\sqrt {-RP}\vartheta \right) \right). \end{equation*}$

Family 5. For $P = R$ & $Q = 0$ ,

$\begin{equation*} \zeta_{21}(\vartheta) = \tan \left( P\vartheta \right), \end{equation*}$

$\begin{equation*} \zeta_{22}(\vartheta) = -\cot \left( P\vartheta \right), \end{equation*}$

$\begin{equation*} \zeta_{23}(\vartheta) = \tan \left( 2\,P\vartheta \right) + \left( \sec \left( 2\,P\vartheta \right) \right), \end{equation*}$

$\begin{equation*} \zeta_{24}(\vartheta) = -\cot \left( 2\,P\vartheta \right) + \left( \csc \left( 2\,P\vartheta \right) \right), \end{equation*}$

and

$\begin{equation*} \zeta_{25}(\vartheta) = \frac{1}{2}\tan \left( \frac{1}{2}\,P\vartheta \right) -\frac{1}{2}\,\cot \left( \frac{1}{2}\,P\vartheta \right), \end{equation*}$

Family 6. For $R = -P$ & $Q = 0$ ,

$\begin{equation*} \zeta_{26}(\vartheta) = -\tanh \left( P\vartheta \right), \end{equation*}$

$\begin{equation*} \zeta_{27}(\vartheta) = -\coth \left( P\vartheta \right), \end{equation*}$

$\begin{equation*} \zeta_{28}(\vartheta) = -\tanh \left( 2\,P\vartheta \right) + \left( i {\it sech} \left( 2\,P\vartheta \right) \right), \end{equation*}$

$\begin{equation*} \zeta_{29}(\vartheta) = -\coth \left( 2\,P\vartheta \right) + \left( { \it csch} \left( 2\,P\vartheta \right) \right), \end{equation*}$

and

$\begin{equation*} \zeta_{30}(\vartheta) = -\frac{1}{2}\tanh \left( \frac{1}{2}\,P\vartheta \right) -\frac{1}{2}\,\coth \left( \frac{1}{2}\,P\vartheta \right). \end{equation*}$

Family 7. For $\Phi = 0$ ,

$\begin{equation*} \zeta_{31}(\vartheta) = -2\,{\frac {P \left(2+Q\vartheta\, \right) }{{Q}^{2}\vartheta\, }}. \end{equation*}$

Family 8. For $Q = \mu$ , $P = p\mu(p\neq0)$ & $R = 0$ ,

$\begin{equation*} \zeta_{32}(\vartheta) = {e}^{\mu \vartheta}-p. \end{equation*}$

Family 9. For $Q = R = 0$ ,

$\begin{equation*} \zeta_{33}(\vartheta) = P\vartheta\,. \end{equation*}$

Family 10. For $Q = P = 0$ ,

$\begin{equation*} \zeta_{34}(\vartheta) = -{\frac {1}{R\vartheta\,}}. \end{equation*}$

Family 11. For $P = 0$ , $R\neq0$ & $Q\neq0$ ,

$\begin{equation*} \zeta_{35}(\vartheta) = -{\frac {Q}{R \left( \cosh \left( Q\vartheta \right) -\sinh \left( Q\vartheta \right) +1\right) }}, \end{equation*}$

and

$\begin{equation*} \zeta_{36}(\vartheta) = -{\frac {Q \left( \cosh \left( Q\vartheta \right) +\sinh \left( Q\vartheta \right) \right) }{R \left( 1 +\sinh \left( Q \vartheta \right) +\cosh \left( Q\vartheta \right) \right) }}. \end{equation*}$

Family 12. For $Q = \mu$ , $R = p\mu(p\neq0)$ & $P = 0$ ,

$\begin{equation*} \zeta_{37}(\vartheta) = {\frac {{e}^{\mu \vartheta}}{1-p{e}^{\mu \vartheta}}}, \end{equation*}$

where $\Phi = Q^2-4RP$ .

3. The execution of m-EDAM

In this section of the study, soliton solutions for the CIE, as specified in Eq (1.1), are generated using the mEDAM approach. First, apply the following variable transformation as part of the method:

$\begin{equation} z(x,y,t) = Z(\vartheta), \quad \vartheta = y+x-\omega t. \end{equation}$

(3.1)

The following NODE is obtained by applying this transformation to (1.1):

$\begin{equation} \begin{split} -\omega Z''+(\lambda(Z^3-Z)-Z'')'+\sigma Z'' = 0. \end{split} \end{equation}$

(3.2)

The following outcome can be obtained from a single integration on (3.2) with a zero integration constant:

$\begin{equation} (\sigma-\omega) Z'+\lambda(Z^3-Z)-Z'' = 0. \end{equation}$

(3.3)

After creating a homogeneous balancing condition between $Z''$ and $Z^3$ , it is determined that $N = 1$ . Substituting $N = 1$ into Eq (2.3), we obtain the ensuing closed-form solution for (3.3):

$\begin{equation} Z(\vartheta) = \sum\limits_{l = -1}^{1}\rho_l(\zeta(\vartheta))^l. \end{equation}$

(3.4)

We acquire an expression based on the terms $\zeta(\vartheta)$ via including (3.4) in (3.3) and bringing together terms with powers alike $\zeta(\vartheta)$ . By setting the coefficients to zero, this expression can be reduced to a set of algebraic nonlinear equations. After using Maple to solve this system, the following two sets of solutions are obtained:

Case 1.

$\begin{equation} \begin{split} \rho_{{0}} = {\frac {Q}{\sqrt {\Phi}}},\rho_{{1}} = 2\,{\frac {R}{\sqrt { \Phi}}},\rho_{{-1}} = 0,\omega = \sigma,\lambda = \Phi,\sigma = \sigma. \end{split} \end{equation}$

(3.5)

Case 2.

$\begin{equation} \begin{split} \rho_{{0}} = \frac{1}{2}\,{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2},\rho_{{1}} = 0,\rho_{{-1} } = {\frac {P}{\sqrt {\Phi}}},\omega = \sigma-3\,\sqrt {\Phi},\lambda = 2\, \Phi,\sigma = \sigma. \end{split} \end{equation}$

(3.6)

By considering Case 1 and using (3.1), (3.4), and the corresponding general solution to (2.4), we construct the ensuing families of soliton solutions for CIE given in (1.1):

Family 1.1. When $\Phi < 0, \quad R\neq0$ ,

$\begin{equation} \begin{split} z_{1,1}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}+\frac{1}{2}\,{ \frac {\sqrt {-\Phi}\tan \left( \frac{1}{2}\,\sqrt {-\Phi}\vartheta \right) } {R}} \right) , \end{split} \end{equation}$

(3.7)

$\begin{equation} \begin{split} z_{1,2}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{ \frac {\sqrt {-\Phi}\cot \left( \frac{1}{2}\,\sqrt {-\Phi}\vartheta \right) } {R}} \right) , \end{split} \end{equation}$

(3.8)

$\begin{equation} \begin{split} z_{1,3}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}+\frac{1}{2}\,{ \frac {\sqrt {-\Phi} \left( \tan \left( \sqrt {-\Phi}\vartheta \right) +\sec \left( \sqrt {-\Phi}\vartheta \right) \right) }{R}} \right) , \end{split} \end{equation}$

(3.9)

$\begin{equation} \begin{split} z_{1,4}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{ \frac {\sqrt {-\Phi} \left( \cot \left( \sqrt {-\Phi}\vartheta \right) +\csc \left( \sqrt {-\Phi}\vartheta \right) \right) }{R}} \right) , \end{split} \end{equation}$

(3.10)

and

$\begin{equation} \begin{split} z_{1,5}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}+\frac{1}{4}\,{ \frac {\sqrt {-\Phi} \left( \tan \left( \frac{1}{4}\,\sqrt {-\Phi}\vartheta \right) -\cot \left( \frac{1}{4}\,\sqrt {-\Phi}\vartheta \right) \right) }{ R}} \right) . \end{split} \end{equation}$

(3.11)

Family 1.2. When $\Phi > 0, \quad R\neq0$ ,

$\begin{equation} \begin{split} z_{1,6}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{ \frac {\sqrt {\Phi}\tanh \left( \frac{1}{2}\,\sqrt {\Phi}\vartheta \right) }{ R}} \right) , \end{split} \end{equation}$

(3.12)

$\begin{equation} \begin{split} z_{1,7}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{ \frac {\sqrt {\Phi}\coth \left( \frac{1}{2}\,\sqrt {\Phi}\vartheta \right) }{ R}} \right) , \end{split} \end{equation}$

(3.13)

$\begin{equation} \begin{split} z_{1,8}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{ \frac {\sqrt {\Phi} \left( \tanh \left( \sqrt {\Phi}\vartheta \right) +i{\rm sech} \left( \sqrt {\Phi}\vartheta \right) \right) } {R}} \right) , \end{split} \end{equation}$

(3.14)

$\begin{equation} \begin{split} z_{1,9}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{ \frac {\sqrt {\Phi} \left( \coth \left( \sqrt {\Phi}\vartheta \right) +{\rm csch} \left( \sqrt {\Phi}\vartheta \right) \right) }{ R}} \right) , \end{split} \end{equation}$

(3.15)

and

$\begin{equation} \begin{split} z_{1,10}(x,y,t) = & {\frac {Q}{\sqrt {\Phi}}}+{\frac {2R}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{4}\,{ \frac {\sqrt {\Phi} \left( \tanh \left( \frac{1}{4}\,\sqrt {\Phi}\vartheta \right) -\coth \left( \frac{1}{4}\,\sqrt {\Phi}\vartheta \right) \right) }{ R}} \right) . \end{split} \end{equation}$

(3.16)

Family 1.3. When $RP > 0$ and $Q = 0$ ,

$\begin{equation} \begin{split} z_{1,11}(x,y,t) = & -i\tan \left( \sqrt {RP}\vartheta \right) , \end{split} \end{equation}$

(3.17)

$\begin{equation} \begin{split} z_{1,12}(x,y,t) = & i\cot \left( \sqrt {RP}\vartheta \right) , \end{split} \end{equation}$

(3.18)

$\begin{equation} \begin{split} z_{1,13}(x,y,t) = & -i \left( \tan \left( 2\,\sqrt {RP}\vartheta \right) +\sec \left( 2\,\sqrt {RP}\vartheta \right) \right) , \end{split} \end{equation}$

(3.19)

$\begin{equation} \begin{split} z_{1,14}(x,y,t) = & i \left( \cot \left( 2\,\sqrt {RP}\vartheta \right) +\csc \left( 2\,\sqrt {RP}\vartheta \right) \right) , \end{split} \end{equation}$

(3.20)

and

$\begin{equation} \begin{split} z_{1,15}(x,y,t) = & \frac{1}{2i} \left( \tan \left( \frac{1}{2}\,\sqrt {RP}\vartheta \right) -\cot \left( \frac{1}{2}\,\sqrt {RP}\vartheta \right) \right) . \end{split} \end{equation}$

(3.21)

Family 1.4. When $RP < 0$ and $Q = 0$ ,

$\begin{equation} \begin{split} z_{1,16}(x,y,t) = & -\tanh \left( \sqrt {-RP}\vartheta \right) , \end{split} \end{equation}$

(3.22)

$\begin{equation} \begin{split} z_{1,17}(x,y,t) = & -\coth \left( \sqrt {-RP}\vartheta \right) , \end{split} \end{equation}$

(3.23)

$\begin{equation} \begin{split} z_{1,18}(x,y,t) = & - \left( \tanh \left( 2\,\sqrt {-RP} \vartheta \right) +i{\rm sech} \left( 2\,\sqrt {-RP}\vartheta \right) \right) , \end{split} \end{equation}$

(3.24)

$\begin{equation} \begin{split} z_{1,19}(x,y,t) = & -\left( \coth \left( 2\,\sqrt {-RP} \vartheta \right) +{\rm csch} \left( 2\,\sqrt {-RP}\vartheta \right) \right) , \end{split} \end{equation}$

(3.25)

and

$\begin{equation} \begin{split} z_{1,20}(x,y,t) = & -\frac{1}{2} \left( \tanh \left( \frac{1}{2}\,\sqrt {-RP} \vartheta \right) +\coth \left( \frac{1}{2}\,\sqrt {-RP}\vartheta \right) \right) . \end{split} \end{equation}$

(3.26)

Family 1.5. When $R = P$ and $Q = 0$ ,

$\begin{equation} \begin{split} z_{1,21}(x,y,t) = & -i \tan \left( P\vartheta \right) , \end{split} \end{equation}$

(3.27)

$\begin{equation} \begin{split} z_{1,22}(x,y,t) = & i\cot \left( P\vartheta \right) , \end{split} \end{equation}$

(3.28)

$\begin{equation} \begin{split} z_{1,23}(x,y,t) = & -i (\tan \left( 2\,P\vartheta \right) +\sec \left( 2 \,P\vartheta \right) ) , \end{split} \end{equation}$

(3.29)

$\begin{equation} \begin{split} z_{1,24}(x,y,t) = & -i \left( -\cot \left( 2\,P\vartheta \right) -\csc \left( 2 \,P\vartheta \right) \right) , \end{split} \end{equation}$

(3.30)

and

$\begin{equation} \begin{split} z_{1,25}(x,y,t) = & -i \left( \frac{1}{2}\,\tan \left( \frac{1}{2}\,P\vartheta \right) -\frac{1}{2}\, \cot \left( \frac{1}{2}\,P\vartheta \right) \right) . \end{split} \end{equation}$

(3.31)

Family 1.6. When $R = -P$ and $Q = 0$ ,

$\begin{equation} \begin{split} z_{1,26}(x,y,t) = & \tanh \left( P\vartheta \right) , \end{split} \end{equation}$

(3.32)

$\begin{equation} \begin{split} z_{1,27}(x,y,t) = & \coth \left( P\vartheta \right) , \end{split} \end{equation}$

(3.33)

$\begin{equation} \begin{split} z_{1,28}(x,y,t) = & \tanh \left( 2\,P\vartheta \right) + i{\rm sech} \left( 2\,P\vartheta \right) , \end{split} \end{equation}$

(3.34)

$\begin{equation} \begin{split} z_{1,29}(x,y,t) = & \coth \left( 2\,P\vartheta \right) + {\rm csch} \left( 2\,P\vartheta \right) , \end{split} \end{equation}$

(3.35)

and

$\begin{equation} \begin{split} z_{1,30}(x,y,t) = & \frac{1}{2}\,\bigg(\tanh \left( \frac{1}{2}\,P\vartheta \right) +\,\coth \left( \frac{1}{2}\,P\vartheta \right)\bigg) . \end{split} \end{equation}$

(3.36)

Family 1.7. When $P = 0$ , $Q\neq0$ and $R\neq0$ ,

$\begin{equation} \begin{split} z_{1,31}(x,y,t) = & 1-\frac{2}{ \cosh \left( Q\vartheta \right) -\sinh \left( Q \vartheta \right) +1 } , \end{split} \end{equation}$

(3.37)

and

$\begin{equation} \begin{split} z_{1,32}(x,y,t) = & 1-2\,{\frac {\cosh \left( Q\vartheta \right) +\sinh \left( Q \vartheta \right) }{\cosh \left( Q\vartheta \right) +\sinh \left( Q \vartheta \right) +1}} . \end{split} \end{equation}$

(3.38)

Family 1.8. When $Q = \mu$ , $R = p\mu(p\neq0)$ and $P = 0$ ,

$\begin{equation} \begin{split} z_{1,33}(x,y,t) = & 1+2\,{\frac {p \rm e ^{\mu\,\vartheta }}{1-p \rm e ^{\mu\,\vartheta }}}, \end{split} \end{equation}$

(3.39)

where $\vartheta = y+x-\sigma t.$

Now, by considering Case 2 and using (3.1), (3.4), and the corresponding general solution to (2.4), we construct the ensuing families of soliton solutions for CIE given in (1.1):

Family 2.1. When $\Phi < 0, \quad R\neq0$ ,

$\begin{equation} \begin{split} z_{2,1}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}+\frac{1}{2}\,{\frac { \sqrt {-\Phi}\tan \left( \frac{1}{2}\,\sqrt {-\Phi}\vartheta \right) }{R}} \right) ^{-1}+\frac{1}{2}\,{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} , \end{split} \end{equation}$

(3.40)

$\begin{equation} \begin{split} z_{2,2}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{\frac { \sqrt {-\Phi}\cot \left( \frac{1}{2}\,\sqrt {-\Phi}\vartheta \right) }{R}} \right) ^{-1}+\frac{1}{2}\,{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} , \end{split} \end{equation}$

(3.41)

$\begin{equation} \begin{split} z_{2,3}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}+\frac{1}{2}\,{\frac { \sqrt {-\Phi} \left( \tan \left( \sqrt {-\Phi}\vartheta \right) +\sec \left( \sqrt {-\Phi}\vartheta \right) \right) }{R}} \right) ^{-1}+\frac{1}{2}{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} , \end{split} \end{equation}$

(3.42)

$\begin{equation} \begin{split} z_{2,4}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{\frac { \sqrt {-\Phi} \left( \cot \left( \sqrt {-\Phi}\vartheta \right) +\csc \left( \sqrt {-\Phi}\vartheta \right) \right) }{R}} \right) ^{-1}+\frac{1}{2}{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} , \end{split} \end{equation}$

(3.43)

and

$\begin{equation} \begin{split} z_{2,5}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}+\frac{1}{4}\,{\frac { \sqrt {-\Phi} \left( \tan \left( \frac{1}{4}\,\sqrt {-\Phi}\vartheta \right) -\cot \left( \frac{1}{4}\,\sqrt {-\Phi}\vartheta \right) \right) }{R}} \right) ^{-1}+\frac{1}{2}\,{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} . \end{split} \end{equation}$

(3.44)

Family 2.2. When $\Phi > 0, \quad R\neq0$ ,

$\begin{equation} \begin{split} z_{2,6}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{\frac { \sqrt {\Phi}\tanh \left( \frac{1}{2}\,\sqrt {\Phi}\vartheta \right) }{R}} \right) ^{-1}+\frac{1}{2}\,{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} , \end{split} \end{equation}$

(3.45)

$\begin{equation} \begin{split} z_{2,7}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{\frac { \sqrt {\Phi}\coth \left( \frac{1}{2}\,\sqrt {\Phi}\vartheta \right) }{R}} \right) ^{-1}+\frac{1}{2}\,{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} , \end{split} \end{equation}$

(3.46)

$\begin{equation} \begin{split} z_{2,8}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{\frac { \sqrt {\Phi} \left( \tanh \left( \sqrt {\Phi}\vartheta \right) +i{ \rm sech} \left( \sqrt {\Phi}\vartheta \right) \right) }{R}} \right) ^{-1}+\frac{1}{2}\,{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} , \end{split} \end{equation}$

(3.47)

$\begin{equation} \begin{split} z_{2,9}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{2}\,{\frac { \sqrt {\Phi} \left( \coth \left( \sqrt {\Phi}\vartheta \right) +{\rm csch} \left( \sqrt {\Phi}\vartheta \right) \right) }{R}} \right) ^{- 1}+\frac{1}{2}\,{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} , \end{split} \end{equation}$

(3.48)

and

$\begin{equation} \begin{split} z_{2,10}(x,y,t) = & {\frac {P}{\sqrt {\Phi}}} \left( -\frac{1}{2}\,{\frac {Q}{R}}-\frac{1}{4}\,{\frac { \sqrt {\Phi} \left( \tanh \left( \frac{1}{4}\,\sqrt {\Phi}\vartheta \right) - \coth \left( \frac{1}{4}\,\sqrt {\Phi}\vartheta \right) \right) }{R}} \right) ^{-1}+\frac{1}{2}\,{\frac {Q}{\sqrt {\Phi}}}+\frac{1}{2} . \end{split} \end{equation}$

(3.49)

Family 2.3. When $RP > 0$ and $Q = 0$ ,

$\begin{equation} \begin{split} z_{2,11}(x,y,t) = & \frac{1}{2i} \left( \tan \left( \sqrt {RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} , \end{split} \end{equation}$

(3.50)

$\begin{equation} \begin{split} z_{2,12}(x,y,t) = & -\frac{1}{2i} \left( \cot \left( \sqrt {RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} , \end{split} \end{equation}$

(3.51)

$\begin{equation} \begin{split} z_{2,13}(x,y,t) = & \frac{1}{2i} \left( \tan \left( 2\,\sqrt {RP}\vartheta \right) +\sec \left( 2\, \sqrt {RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} , \end{split} \end{equation}$

(3.52)

$\begin{equation} \begin{split} z_{2,14}(x,y,t) = & -\frac{1}{2i} \left( \cot \left( 2\,\sqrt {RP}\vartheta \right) +\csc \left( 2\, \sqrt {RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} , \end{split} \end{equation}$

(3.53)

and

$\begin{equation} \begin{split} z_{2,15}(x,y,t) = & \frac{1}{i} \left( \tan \left( \frac{1}{2}\,\sqrt {RP}\vartheta \right) -\cot \left( \frac{1}{2} \,\sqrt {RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} . \end{split} \end{equation}$

(3.54)

Family 2.4. When $RP < 0$ and $Q = 0$ ,

$\begin{equation} \begin{split} z_{2,16}(x,y,t) = & -\frac{1}{2} \left( \tanh \left( \sqrt {-RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} , \end{split} \end{equation}$

(3.55)

$\begin{equation} \begin{split} z_{2,17}(x,y,t) = & -\frac{1}{2} \left( \coth \left( \sqrt {-RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} , \end{split} \end{equation}$

(3.56)

$\begin{equation} \begin{split} z_{2,18}(x,y,t) = & -\frac{1}{2} \left( \tanh \left( 2\,\sqrt {-RP}\vartheta \right) +i{\rm sech} \left( 2\,\sqrt {-RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} , \end{split} \end{equation}$

(3.57)

$\begin{equation} \begin{split} z_{2,19}(x,y,t) = & -\frac{1}{2} \left( \coth \left( 2\,\sqrt {-RP}\vartheta \right) +{\rm csch} \left( 2\,\sqrt {-RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} , \end{split} \end{equation}$

(3.58)

and

$\begin{equation} \begin{split} z_{2,20}(x,y,t) = & - \left( \tanh \left( 2\,\sqrt {-RP}\vartheta \right) +\coth \left( 2\,\sqrt {-RP}\vartheta \right) \right) ^{-1}+\frac{1}{2} . \end{split} \end{equation}$

(3.59)

Family 2.5. When $R = P$ and $Q = 0$ ,

$\begin{equation} \begin{split} z_{2,21}(x,y,t) = & {\frac {1}{2i\tan \left( P\vartheta \right) }}+\frac{1}{2} , \end{split} \end{equation}$

(3.60)

$\begin{equation} \begin{split} z_{2,22}(x,y,t) = & -{\frac {1}{2i\cot \left( P\vartheta \right) }}+\frac{1}{2} , \end{split} \end{equation}$

(3.61)

$\begin{equation} \begin{split} z_{2,23}(x,y,t) = & {\frac {1}{2i \left( \tan \left( 2\,P\vartheta \right) +\sec \left( 2\,P\vartheta \right) \right) }}+\frac{1}{2} , \end{split} \end{equation}$

(3.62)

$\begin{equation} \begin{split} z_{2,24}(x,y,t) = & {\frac {1}{2i \left( -\cot \left( 2\,P\vartheta \right) -\csc \left( 2\,P\vartheta \right) \right) }}+\frac{1}{2} , \end{split} \end{equation}$

(3.63)

and

$\begin{equation} \begin{split} z_{2,25}(x,y,t) = & {\frac{1}{i \left( \tan \left( \frac{1}{2}\,P \vartheta \right) -\cot \left( \frac{1}{2}\,P\vartheta \right) \right) }}+\frac{1}{2} . \end{split} \end{equation}$

(3.64)

Family 2.6. When $R = -P$ and $Q = 0$ ,

$\begin{equation} \begin{split} z_{2,26}(x,y,t) = & -{\frac {1}{2\tanh \left( P\vartheta \right) }}+\frac{1}{2} , \end{split} \end{equation}$

(3.65)

$\begin{equation} \begin{split} z_{2,27}(x,y,t) = & -{\frac {1}{2\coth \left( P\vartheta \right) }}+\frac{1}{2} , \end{split} \end{equation}$

(3.66)

$\begin{equation} \begin{split} z_{2,28}(x,y,t) = & {\frac {1}{2 \left( -\tanh \left( 2\,P \vartheta \right) -i{\rm sech} \left( 2\,P\vartheta \right) \right) }}+\frac{1}{2} , \end{split} \end{equation}$

(3.67)

$\begin{equation} \begin{split} z_{2,29}(x,y,t) = & {\frac {1}{2 \left( -\coth \left( 2\,P \vartheta \right) -{\rm csch} \left( 2\,P\vartheta \right) \right) }}+\frac{1}{2} , \end{split} \end{equation}$

(3.68)

and

$\begin{equation} \begin{split} z_{2,30}(x,y,t) = & {\frac {1}{ \left( -\tanh \left( \frac{1}{2} \,P\vartheta \right) -\coth \left( \frac{1}{2}\,P\vartheta \right) \right) }}+\frac{1}{2} . \end{split} \end{equation}$

(3.69)

Family 2.7. When $Q = \mu$ , $P = p\mu(p\neq0)$ and $R = 0$ ,

$\begin{equation} \begin{split} z_{2,31}(x,y,t) = & {\frac {p}{ \rm e ^{\mu\,\vartheta }-p}}+1, \end{split} \end{equation}$

(3.70)

where $\vartheta = y+x-(\sigma-3\, \sqrt {\Phi}) t.$

4. Discussion and graphs

In this section, we show schematic illustrations of several wave configurations seen in the system that we are studying. Using m-EDAM, an improved version of EDAM, we discovered and visualised wave patterns in density and 3D graphs. These examples reveal that the soliton solutions generated primarily appear as cuspon or kink solitons. A kink soliton is the perfect example of a lone wave that has a sudden change or discontinuity in the variable or field it represents. It is also known as a domain wall or topological soliton. Kink solitons typically appear in systems that are going through a scalar field phase transition. On the other hand, a cuspon soliton, also called a cusped soliton, is a unique type of single wave with a cusp, or sharp peak, at the apex. It functions as a solution to a particular class of nonlinear wave equations in which the propagating single wave maintains its form.

Particular phenomena inherent to plasma physics account for the predominance of kink and cuspon solitons as principal soliton solutions to the CIE, a fundamental model of ion-acoustic waves to plasma. Particle interactions, magnetic fields, and nonlinear behaviours are all aspects of plasma dynamics that play a role in the creation and spread of single waves. The CIE, which records nonlinear phenomena in plasma, clarifies the creation of isolated waves with characteristic kinks and cusps. The CIE's mathematical foundations and plasma dynamics interact to produce cuspon solitons, which have sharp peaks, and kink solitons, which show sudden changes in fields as a result of phase transitions. As tracers of nonlinear interactions, these solitons provide information about the complex dynamics of ion-acoustic waves. Understanding these solitons helps to better understand plasma physics and provides insight into the wave behaviours of intricate plasma systems.

The m-EDAM has proven to be a reliable, practical, and easy-to-use tool that offers more general wave solutions for nonlinear models in mathematical physics and engineering. Moreover, our results, under some assumptions, agree with the results of other transformation-based methods such as the Riccati equation expansion method, $\left(\frac{G'}{G}\right)$ -expansion method, sub-equation method, F-expansion method, and many more. For instance, all ten solutions by the modified F-expansion method in the study of Seadawy et al. ^[52] can be easily obtained from some of our 37 solutions produced by m-EDAM by assigning values to $P$ , $Q$ , and $R$ presented in our solutions. Our findings thus serve as a generalisation of the outcomes of the various approaches that were previously described. This adds an important new layer to the existing understanding of the CIE and demonstrates the adaptability and broader applicability of our research.

The graphical discussion presents the visualization of various soliton solutions derived from the mathematical models. illustrates the kink soliton solution $z_{1, 6}$ at $t = 0$ for specific parameter values. In , the bright kink soliton solution $z_{1, 9}$ is depicted at $t = 1$ , showcasing its behavior under different parameter settings. displays the bright kink soliton solution $z_{1, 18}$ at $t = 10$ , highlighting its evolution over time. Meanwhile, and focus on the imaginary part of cuspon soliton solutions $z_{1, 21}$ and $z_{1, 24}$ at different time instances, demonstrating their unique characteristics. In , the kink soliton solution $z_{1, 32}$ is visualized at $t = 0$ , revealing its structure under specific parameter conditions. and showcase the dark kink soliton solutions $z_{2, 10}$ and $z_{2, 19}$ at different time points, emphasizing their behavior under varying parameters. presents the bright-dark soliton solution $z_{2, 23}$ at $t = 0$ , while illustrates the kink soliton solution $z_{2, 31}$ at $t = 2$ , providing insights into their dynamics under specific parameter configurations. Overall, these visualizations offer valuable insights into the behavior and characteristics of the soliton solutions under investigation.

Figure 1. The 3D and contour plots of the kink soliton solution

$z_{1, 6}$ stated in (3.12) are plotted for

$P = 1, Q = 3, R = 2, \sigma = 10, t = 0$ .

DownLoad: Full-Size Img PowerPoint

Figure 2. The 3D and contour plots of the bright kink soliton solution

$z_{1, 9}$ stated in (3.15) are plotted for

$P = 8, Q = 10, R = 2, \sigma = 50, t = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 3. The 3D and contour plots of the bright kink soliton solution

$z_{1, 18}$ stated in (3.24) are plotted for

$P = 9, Q = 0, R = -4, \sigma = 20, t = 10$ .

DownLoad: Full-Size Img PowerPoint

Figure 4. The 3D and contour plots of the imaginary part of cuspon soliton solution

$z_{1, 21}$ stated in (3.27) are plotted for

$P = 3, Q = 0, R = 3, \sigma = 30, t = 100$ .

DownLoad: Full-Size Img PowerPoint

Figure 5. The 3D and contour plots of the imaginary part of cuspon soliton solution

$z_{1, 24}$ stated in (3.30) are plotted for

$P = 6, Q = 0, R = 6, \sigma = 40, t = 50$ .

DownLoad: Full-Size Img PowerPoint

Figure 6. The 3D and contour plots of the kink soliton solution

$z_{1, 32}$ stated in (3.38) are plotted for

$P = 0, Q = 10, R = 5, \sigma = 100, t = 0$ .

DownLoad: Full-Size Img PowerPoint

Figure 7. The 3D and contour plots of the dark kink soliton solution

$z_{2, 10}$ stated in (3.49) are plotted for

$P = 2, Q = 10, R = 8, \sigma = 1, t = 5$ .

DownLoad: Full-Size Img PowerPoint

Figure 8. The 3D and contour plots of the kink soliton solution

$z_{2, 19}$ stated in (3.58) are plotted for

$P = 5, Q = 0, R = -5, \sigma = 4, t = 15$ .

DownLoad: Full-Size Img PowerPoint

Figure 9. The 3D and contour plots of the bright-dark soliton solution

$z_{2, 23}$ stated in (3.62) are plotted for

$P = 6, Q = 0, R = 6, \sigma = 100, t = 0$ .

DownLoad: Full-Size Img PowerPoint

Figure 10. The 3D and contour plots of the kink soliton solution

$z_{2, 31}$ stated in (3.70) are plotted for

$P = 10, \mu = 5, p = 2, Q = 5, R = 0, \sigma = 0, t = 2$ .

DownLoad: Full-Size Img PowerPoint

5. Conclusions

In conclusion, our study has thoroughly investigated the propagation of solitons in the CIE, a noteworthy model in plasma sciences, by employing the upgraded m-EDAM. By converting the CIE into a NODE and assuming a closed-form solution, numerous cuspon and kink soliton solutions have been identified, including hyperbolic, trigonometric, exponential, and rational functions. Density charts and 3D models that illustrate the propagating behaviour of specific soliton solutions provide important insights into the behaviour of soliton propagation processes, with direct applicability to domains related to the CIE. The results shown here are remarkable and demonstrate the novelty of using m-EDAM to the CIE, as they have not been covered in academic literature before. These discoveries advance our understanding of temporal evolution processes and nonlinear dynamics, which has significant implications for our comprehension of related physical phenomena. Moreover, the efficacy and uniformity of the employed techniques in this study underscore their broader significance for nonlinear problems across other scientific domains. While the utilisation of m-EDAM has significantly enhanced our understanding of soliton dynamics and their impact on the CIE, it is imperative to acknowledge the constraints of the methodology, particularly with the fair handling of the highest derivative and nonlinear variables. Despite this limitation, the work shows that our understanding of soliton behaviour and nonlinear dynamics will continue to expand and opens up new directions for future research in related fields. That is, the future scope of the current study lies in extending the proposed m-EDAM to study soliton events in nonlinear stochastic models, which will shed further light on how mechanistic and stochastic dynamics interact.

Author contributions

Naveed Iqbal: Conceptualization, Investigation and Writing-review & editing; Muhammad Bilal Riaz: Methodology, Conceptualization, Formal analysis and Administration and Funding acquisition; Meshari Alesemi: Methodology; Taher S. Hassan: Software and Investigation; Ali M. Mahnashi: Validation and Resources; Ahmad Shafee: Formal analysis and Resources. All authors have read and agreed to the published version of the manuscript.

Use of AI tools declaration

The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

This article has been produced with the financial support of the European Union under the REFRESH (Research Excellence For Region Sustainability and High-tech Industries) project number CZ.10.03.01/00/22_003/0000048 via the Operational Programme Just Transition.

The authors are thankful to the Deanship of Graduate Studies and Scientific Research at University of Bisha for supporting this work through the Fast-Track Research Support Program.

Conflict of interest

The authors declare that they have no conflicts of interest.

References

[1]	Q. X. Zeng, J. X. Dong, Y. Meng, et al. Progress in epidemiology of human cytomegalovirus Infection, Shandong medicine, 57 (2017), 1131-1133.
[2]	F. Chiuppesi, T. Kaltcheva, Z. Meng, et al., Identification of a continuous neutralizing epitope within UL128 of human cytomegalovirus, J. Virol., 91 (2017), 1-16.
[3]	S. E. Jackson, G. X. Sedikides, G. M. Mason, et al. Human Cytomegalovirus (HCMV)-Specific CD4⁺ T Cells are polyfunctional and can respond to HCMV-Infected dendritic cells in vitro, J. Virol., 91 (2017), 1-16.
[4]	D. Song, H. Mei, Research Progress of congenital cytomegalovirus infection in newborns, Medical Recapitulate, 23 (2017), 4453-4457.
[5]	D. Zhu, C. Pan, J. Sheng, et al. Human cytomegalovirus reprogrammes haematopoietic progenitor cells into immunosuppressive monocytes to achieve latency, Nat. Microbiol., 3 (2018), 503-513. doi: 10.1038/s41564-018-0131-9
[6]	Y. H. Li, Analysis of the correlation between breastfeeding of HCMV infected mothers and HCMV infection of newborns, MA. Thesis, Qingdao University, 2015.
[7]	D. C. Moylan, S. K. Pati, S. A. Ross, et al. Breast milk HCMV viral load is associated with the establishment of breast milk CMV-pp65-specific CD8 T cells in Human CMV infected mothers, J. Infect. Dis., 216 (2017), 1176-1179. doi: 10.1093/infdis/jix457
[8]	W. F. Wu, Progress in the treatment of CMV infection, Chinese Medical Journal, 38 (2003), 10-12.
[9]	K. Wang, W. Wang, S. Song, Dynamics of an HBV model with diffusion and delay, J. Theor. Biol., 253 (2008), 36-44. doi: 10.1016/j.jtbi.2007.11.007
[10]	G. Alter, D. Heckerman, A. Schneidewind, et al. HIV-1 adaptation to NK-cell-mediated immune pressure, Nature, 476 (2017), 96-100.
[11]	W. O. Kermack, A. G. Mckendrick, A contribution to the mathematical theory of epidemics, B. Math. Biol., 53 (1991), 57-87.
[12]	X. N. Han, The transmission dynamies of SARS, MA. Thesis, PLA Academy of Military Sciences, 2006.
[13]	R. Zhao, H. Wang, X. Wang, et al. Steroid therapy and the risk of osteonecrosis in SARS patients: a dose-response meta-analysis, Osteoporosis International, 28 (2016), 1027-1034.
[14]	L. F. Zhang, Comparison and parameter estimation between deterministic model and stochastic model of infectious disease transmission, MA. Thesis, Southwest Jiaotong University, 2010.
[15]	A. Q. Miao, J. Zhang, T. Zhang, et al. Threshold dynamics of a stochastic SIR model with vertical transmission and vaccination, Comput. Math. Method. M., 2017 (2017), 1-10.
[16]	P. Y. Xia, Dynamic behavior of several random virus models, Ph.D thesis, Northeast Normal University, 2018.
[17]	C. Y. Ji, Asymptotic behavior of stochastic biological model and infectious disease model, Ph.D thesis, Northeast Normal University, 2011.
[18]	Y. Asai, C. Tomás, X. Han, et al. A random model for immune response to virus in fluctuating environments, Springer International Publishing, 2016.
[19]	Y. Wang, J. Liu, Y. Y. Liu, et al. Establishment of mouse brain latent cytomegalovirus activation model, Progress in Modern Biomedicine, 15 (2015), 4414-4418.
[20]	H. Y. Duan, T. Yu, Diagnosis and treatment progress of cytomegalovirus infection, Chinese Journal of Obstetrics and Gynecology, 6 (2010), 68-71.
[21]	M. A. Nowak, C. Bangham, Population dynamics of immune responses to persistent viruses, Science, 272 (1996), 74-79. doi: 10.1126/science.272.5258.74
[22]	M. A. Nowak, S. Bonhoeffer, A. M. Hill, et al. Viral dynamics in hepatitis B virus infection, Proceedings of the National Academy of Sciences, 93 (1996), 4398-4402. doi: 10.1073/pnas.93.9.4398
[23]	X. Q. Niu, W. D. Li, G. F. Zhu, et al. Modeling the transmission dynamics of hepatitis B Virus and data assimilation forecasting, Mathematics in Practice and Theory, 45 (2015), 205-211.
[24]	J. M. Conway, D. Coombs, A stochastic model of latently infected cell reactivation and viral blip generation in treated HIV patients, PLoS Comput. Biol., 7 (2011), 1-24.
[25]	C. Fraser, N. M. Ferguson, R. M. Anderson, et al. The role of antigenic stimulation and cytotoxic T Cell activity in regulating the long-term immunopathogenesis of HIV: mechanisms and clinical implications, Proceedings: Biological Sciences, 268 (2001), 2085-2095. doi: 10.1098/rspb.2001.1777
[26]	C. Fraser, N. M. Ferguson, R. M. Anderson, Quantification of intrinsic residual viral replication in treated HIV-infected patients, Proceedings of the National Academy of Sciences of the United States of America, 98 (2001), 15167-15172. doi: 10.1073/pnas.261283598
[27]	W. Zhang, L. M. Wahl, P. Yu, Viral blips may not need a trigger: how transient viremia can arise in deterministic in-host models, Siam Rev., 56 (2014), 127-155. doi: 10.1137/130937421
[28]	S. Wang, J. Zhang, F. Xu, et al. Dynamics of virus infection models with density‐dependent diffusion, Comput. Math. Appl., 74 (2017), 1-20. doi: 10.1016/j.camwa.2017.05.001
[29]	M. Wei, L. Hu, X. Mao, Neutral stochastic functional differential equations with Lévy jumps under the local Lipschitz condition, Advances in Difference Equations, 2017 (1017), 57.
[30]	R. Khasminskii, Stochastic Stability of Differential Equations, Stochastic stability of differential equations, 1980.
[31]	X. X. Liao, Theory methods and application of sability, 2nd Edition, Huazhong University of Science and Technology Press, Wuhan, 2010.
[32]	N. He, W. D. Wang, A. R. Zhou, et al. Dynamics of stochastic HIV model based on saturation incidence rate, Journal of Southwest University (Natural Science Edition), 40 (2018), 123-125.

This article has been cited by:

1.	Musawa Yahya Almusawa, Hassan Almusawa, Dark and bright soliton phenomena of the generalized time-space fractional equation with gas bubbles, 2024, 9, 2473-6988, 30043, 10.3934/math.20241451
2.	Azzh Saad Alshehry, Safyan Mukhtar, Ali M. Mahnashi, Optical fractals and Hump soliton structures in integrable Kuralay-Ⅱ system, 2024, 9, 2473-6988, 28058, 10.3934/math.20241361

Reader Comments

Your name:*

Email:*
© 2020 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)