**1. Introduction**

In the current century, many entropy problems have been expressed by using mathematical models that are nonlinear partial differential equations. New results in the last few years have shown that the relation between non-standard entropies and nonlinear partial differential equations can be applied on new nonlinear wave equations inspired by quantum mechanics. Nonlinear models of the celebrated Klein–Gordon and Dirac equations have been found to admit accurate time dependent soliton-like solutions with the shapes of the so-called q-plane waves. Such q-plane waves are generalizations of the complex exponential plane wave solutions of the linear Klein–Gordon and Dirac equations [1]. Wave progressing of soliton forming and its application in the differential equation has been noticeable in the last few years. The physical phenomena of nonlinear partial differential equations (NLPDEs) may connect to many areas of sciences, for example plasma physics, optical fibers, nonlinear optics, fluid mechanics, chemistry, biology, geochemistry, and engineering sciences. The nonlinear Schrödinger equations describe wave propagation in optical fibers with nonlinear impacts [2–4].

Various numeric and analytic techniques have been used to seek solutions for nonlinear differential equations such as the homotopy perturbation scheme [5], the Adams–Bashforth–Moulton method [6], the shooting technique with fourth-order Runge–Kutta scheme [7–10], the group preserving method [11], the finite forward difference method [12,13], the Adomian decomposition method [14,15], the sine-Gordon expansion method [16–18], the modified auxiliary expansion method [19], the modified exp(−*ϕ* (*ξ*)) expansion function method [20,21], the improved Bernoulli sub-equation method [22,23], the Riccati–Bernoulli sub-ODE method [24], the modified exponential function method [25], the improved tan(*φ* (*ξ*) /2) [26,27], the Darboux transformation method [28,29], the double - *G*- *<sup>G</sup>* , <sup>1</sup> *G*- expansion method [30,31], the - 1 *G*- expansion method [32,33], the decomposition Sumudu-like-integral transform method [34], and the inverse scattering method [35].

In recent years, many researchers have carried out investigations on the governing models in optical fibers. The nonlinear Schrödinger equation, involving cubic and quartic-order dispersion terms, has been investigated to seek the exact optical soliton solutions via the undetermined coefficients method [36], the modified Kudryashov approach [37], the complete discrimination system method [38], the generalized tanh function method [39], the sin-cosine method, as well as the Bernoulli equation approach [40], the semi-inverse variation method [41], a simple equation method [3], and the extended sinh-Gordon expansion method [42].

Now, optical solitons are the exciting research area of nonlinear optics studies, and this research field has led to tremendous advances in their extensive applications. It is identified that the dynamics of nonlinear optical solitons and Madelung fluids are based on the generalized nonlinear Schrödinger dispersive equation and resonant nonlinear Schrödinger dispersive equation. In the research of chirped solitons in Hall current impacts in the field of quantum mechanics, a specific resonant term must be given [43].

Dispersion and nonlinearity are the two key elements for the propagation of solitons over intercontinental ranges. Normally, group velocity dispersion (GVD) leveling with self-phase modulation in a sensitive way allows such solitons to maintain long distance travel. Now, it could occur that GVD is minuscule and therefore completely overlooked, so in this condition, the dispersion impact is rewarded for by third-order (3OD) and fourth-order (4OD) dispersion impacts. This is generally referred to as solitons that are cubic-quartic (CQ). This term was implemented in 2017 for the first time. This model was later extensively researched through different points of view such as the semi-inverse variation principle [41], Lie symmetry [44], conservation rules [45], and the system of undetermined coefficients [37]. Consider the nonlinear Schrödinger and resonant nonlinear Schrödinger equations in the appearance of 3OD and 4OD without both GVD and disturbance. The equations are as follows:

$$
\alpha \| u\_t + i\alpha u\_{xxx} + \beta u\_{xxxx} + cF\left( |u|^2 \right) u = 0,\tag{1}
$$

$$
\delta u\_t + i\alpha u\_{xxx} + \beta u\_{xxxx} + cF\left(|u|^2\right)u + c\_3\left(\frac{|u|\_{xx}}{|u|}\right)u = 0. \tag{2}
$$

In Equations (1) and (2), *u* (*x*, *t*) is the complex valued wave function and *x* (space) and *t* (time) are independent variables. The coefficients *α* and *β* are real constants, while *c*<sup>3</sup> is the Bohm potential that occurs in Madelung fluids. The Bohm potential term of disturbance generates quantum behavior, so that quantum characteristics are closely related to their special characteristics. Therefore, we have the chirped NLSE's disturbance expression giving us the introduction of the theory of hidden variables. Therefore, it will be more crucial to retrieve accurate solutions for the development of quantum mechanics from disturbed chiral (resonant) NLSE [46]. Furthermore, the functional *F* is a real valued algebraic function that represents the source of nonlinearity and *F* - |*u*| 2 *u* : *C* → *C*. In more detail, the function *F* - |*u*| 2 *u* is *p*-times continuously differentiable, so that:

$$F\left(|u|^2\right)u \in \bigcup\_{m,n=1}^{\infty} \mathbb{C}^p\left((-n,n) \times (-m,m) : R^2\right).$$

Suppose that *F* (*u*) = c1*u* + *c*2*u*2, so Equations (1) and (2) can be rewritten as:

$$
\alpha \dot{u} u\_t + \dot{a} u\_{xxx} + \beta u\_{xxxx} + \left(c\_1 |u|^2 + c\_2 |u|^4\right) u = 0,\tag{3}
$$

$$
\delta u\_t + i a u\_{xxx} + \beta u\_{xxxx} + \left(c\_1|u|^2 + c\_2|u|^4\right)u + c\_3\left(\frac{|u|\_{xx}}{|u|}\right)u = 0. \tag{4}
$$

Equation (3) was investigated by making *c*<sup>2</sup> = 0 in [47] via the Kudryashov approach. The conservation laws to obtain the conserved densities for Schrödinger's nonlinear cubic-quarter equation have been analyzed in Kerr and power-law media [45]. The undetermined coefficients method has been employed to construct bright soliton and singular soliton solutions of Equation (1), when nonlinearity has been taken into consideration in the instances of the Kerr law and power law [37]. In this study, we use two methods to investigate soliton solutions of the cubic-quartic nonlinear Schrödinger equation and cubic-quartic resonant nonlinear Schrödinger equation with the parabolic law, namely Equations (3) and (4).
