Assorted Spatial Optical Dynamics of a Generalized Fractional Quadruple Nematic Liquid Crystal System in Non-Local Media

: Nematicons upgrade the recognition of light localization in the reorientation of non-local media. the current research employs a powerful integral scheme using a different procedure, namely, the modified simple equation method (MSEM), to analyze nematicons in liquid crystals from the controlling model. The expanded MSEM is investigated to enlarge the applicability of the standard one. The suggested expansion depends on merging the MSEM and the ansatz method. The new generalized nonlinear n -times quadruple power law is included. With the aid of the symbolic computational package Mathematica, new explicit complex hyperbolic, periodic, and more exact spatial soliton solutions are derived. Moreover, the related existence constraints are obtained. To show the dynamical properties of some of the obtained nematicons, three-dimensional profiles with corresponding contours are depicted with the choice of appropriate values of arbitrary parameters. The fractional impacts in various applicable senses are analyzed to investigate the generality of the considered model.


Introduction
Nonlinear functional problems, including differentials, have wide areas of applications in various fields of science, engineering, medicine, statistics, and artificial intelligence [1][2][3][4]; see also the references included therein.They modeled many real-life phenomena.Therefore, to explain these phenomena and draw conclusions in scientific and practical ways, applied scientists have devoted attention to solving these models.For this purpose, a number of analytic, series-expansion, and numerical schemes were developed, for instance, the improved BPNN method [5], the Adomian decomposition method [6], the variational iteration method [7], the differential transform method [8], the residual power series method [9], the Lie symmetry approach [10], the generalized auxiliary equation method [11], the generalized Riccati equation mapping method [12], the improved Bernoulli sub-equation function technique [13], the G ′ G 2 and F-expansion methods [14], the modified Kudryashov procedure [15], the Sardar sub-ODE method [16], the generalized exponential rational function method [17], and the extended rational sine-cosine, sinh-cosh and sinh-Gordon equation expansion methods [18].
Recently, spatial optical solitons in the applied theory of nonlinear optics, particularly in metasurfaces, metamaterials, liquid crystals, etc., have received considerable attention.In the field of liquid crystals, nematicons, first introduced by Alberucci and Assanto to describe spatial solitons produced by a special type of optical nonlinearity found in nematic liquid crystals, are a class of dielectric media, utilized in processing electronic, thermal, and photo refractive phenomena [19,20].Nematicons present a noteworthy platform for exploring optical processing and switching techniques depending on self and cross localization of light.While the observed effects exhibit inherent slowness, the potential for a faster response can be obtained using other materials or geometries, leading the way towards applied optical communication.
The formal dimensionless soliton dynamics in liquid crystals are given using the coupled system of NPDEs [21]: The dependent variables q(x, t) and p(x, t) symbolize the wave profile and the tilt angle of the liquid crystal molecule, respectively.In Equation ( 1), i = √ −1, and i ∂q ∂t describes the temporal evolution of nematicons, while ∂x 2 stands for the group velocity dispersion (GVD).The nonlinear operator H includes the Kerr, power, parabolic, and dual-power types of nonlinearity.Moreover, α's are assumed to be the parameters.
A novel generalized nematic liquid crystal model, with Equations ( 1) and (2) and the nonlinear nth-quadruple power non-Kerr law, is given by where n is a nonzero integer, and C i 's are assumed to be scalars coefficients.Continuing our previous analysis [31], we analytically processed the liquid crystal governing model, Equations ( 1) and (2), with the quadruple power of nonlinearity in Equation (3), applying the MSEM and its expansion.The outline of this article is as follows: Section 2 includes a brief discussion of the modified simple equation scheme and its expansion for general dimensionless NPDE.Sections 3 and 4 present an analytic treatment of the considered model with formal closed nematicons.In Section 5, the dynamics of some of the obtained nematicons are depicted.To cover more general cases and due to the significance of fractional calculus, the impact of the modified Riemann-Liouville, beta, and modified conformable fractional derivatives have been customized graphically in Section 6.Finally, a concise discussion and the conclusions are provided.

Outline of the Methodology
The main steps of the considered scheme, the MSEM, are demonstrated to process the formal dimensionless nonlinear evolution model given by where F is a polynomial in q(x, t), with its total space x and time t as the partial derivatives.
Step 2. The MSEM is performed by assuming the solution of Equation ( 6) as where ) are constants to be found later.ϕ(η) is an undetermined realvalued function to be calculated later.
Step 3. We compute the positive integer N by implementing the homogeneous balance between terms containing the highest order derivative and nonlinearity in the completely integrated version of Equation ( 6).
Step 4. We replace the assumed ansatz in Equation ( 7) with the value of N in the previous step and its essential derivatives in Equation ( 6).As a result, a polynomial of ϕ −i (η), i = 0, 1, 2, . .., is obtained.We gather the terms of the same power of ϕ −i (η) and make them vanish for each i, and we deduce a mixed algebraic-differential system.By solving this system, the values of ϕ(η) and A i 's are derived.To completely determine the exact solution of Equation ( 4), we place the results into Equation (5).
Consequently, the expansion of the modified simple equation algorithm is based on replacing the unknown function ϕ(η) with a helpful hyperbolic or trigonometric function.Depending on this assumption, the mixed algebraic-differential system should result in some algebraic-functional equations.To clarify, collecting the coefficients of the linearly independent expressions in the obtained equation and equating to zero, a free-functional algebraic system will be obtained.The procedure depends on the assumed ϕ as mentioned in [32].By solving this system for the undetermined parameters, namely the constraints for the solution's existence, and making the backward substitution, the closed form solitons can be be derived.

Mathematical Treatment
Following the procedure in the previous section with application to the coupled system in Equations ( 1) and (2), subject the to non-Kerr operator in Equation (3), we consider the traveling-wave transform η = β(x − ν t) and the solutions of the form q(x, t) = u(η)e i(−κx+ωt+θ 0 ) , ( 8) where ν is the propagation wave speed, and κ, ω, and θ 0 represent the soliton frequency, the wave number of the soliton, and the phase constant, respectively.Substituting Equations ( 8) and ( 9) into Equations ( 1) and (2) and separating the obtained NODEs into real and imaginary parts, we obtain (−2α The imaginary component in Equation ( 12) leads to the constraint condition regarding the soliton speed Thus, Equations ( 10) and ( 11) are deducted to the following completely integrable system of NODEs: We define the following new dependent variables to transform Equations ( 14) and ( 15) into Applying the homogeneous balance principle between the terms containing the highest order derivatives {U U ′′ , V V ′′ } and nonlinear terms V 2 U 2 , U 4 in Equations ( 17) and (18), respectively, results in the following system: This implies that N = M = 1, where M and N represent the corresponding balance integers of U and V, respectively.
Therefore, the formal solutions of Equations ( 17) and ( 18) are In the next section, exact nematicons for our model are derived.

Application
To complete the derivation of the closed-form soliton solutions, we substitute Equation (19) into Equations ( 17) and (18), collect the coefficients of ϕ −i (η), i = 0, 1, 2, 3, 4, with the same power in both obtained sets of differential-algebraic systems, make them vanish, and we obtain The nontrivial solutions with different structures and corresponding existence constraints are listed as follows: Case 1.The unknown function ϕ can be determined directly by solving the linear ordinary differential equations that appear in Equation (22), Equation (23), Equation (27), or Equation (28).With the first choice, we obtain , and where γ 1 , and γ 2 ̸ = 0 are the constants of integration.For the parameters set in Equations (30) and (31), the closed-form spatial soliton solutions of our model will be Case 2. The MSEM provides the possibility of merging two linear differential equations with different orders.In our case, Equation (22) and Equation (28) provide Subsequently, Upon integrating Equation (36), we obtain Consequently, and with Equation ( 35), we obtain That yields where γ 1 , and γ 2 are the constants of integration.The undetermined parameters, as well as the constraints depending on the obtained ϕ, are as follows: provided that α 1 α 2 and α 4 α 5 are negative.The exact nematicons of our problem in Equations ( 1)-(3) will be with Case 3. To expand the use of the MSEM, we assume that ϕ(η) = γ 1 η + γ 0 , with arbitrary constants γ 0 and γ 1 ̸ = 0. We substitute them into the system of Equations ( 20)-( 29) and solve for possible parameters, which gives provided that C 0 = C 1 = C 3 = 0.That is, applying the backward substitution through Equation (19), Equation ( 16), Equation (13), and Equations ( 8) and ( 9), the nematic wave profile and tilt angle of the reduced nth-quadruple system should be Case 4. In this case, the ansatz method is merged with the MSEM by assuming that ϕ(η) = sinh(η).We substitute this into the system of Equations ( 20)-( 29), collect the coefficients of ϕ −i , ϕ −i ϕ ′ , and make them vanish, and the conducted algebraic system is solvable for possible parameters with provided that α 1 α 2 < 0. The nematic wave and tilt angle of our model should be Case 5.As in the previous case, with ϕ(η) = cosh(η), two sets of existence and constraint parameters are found and listed as follows: Set 1.
The common parameters are , provided that α 1 α 2 < 0. The nematicons for the last results of the liquid crystals model in Equations ( 1)-( 3) should be Case 6.By inserting ϕ(η) = sin(η) into the mixed system of Equations ( 20)-( 29) and solving for the included parameters, we obtain Set 1.
The common parameters are , provided that α 1 α 2 < 0. The spatial solitons for the parameters in Set 2 should be Case 7. In the case of ϕ(η) = cos(η), we obtain Set 1.
The common parameters are , provided that α 1 α 2 < 0. The exact spatial soliton solutions for the parameters in Set 2 should be

Graphical Representations
For comparison purposes, the basics of the comparative analysis, showing the dynamical behavior graphically of soliton solutions gives researchers vital indications of the wave motion of the particles under study.Here, several of the spatial solitons derived in the previous section are depicted in 3D.The corresponding contour plots are also listed to further clarify the wave motion.

Fractional Impacts
In the current section, the physical influences of the time-fractional modified Reimann-Liouville, β, and M-truncated derivatives on the solutions to the time-fractional version of the considered nematic liquid crystal model are shown and compared.
The mentioned derivatives are the most recently derived.In addition, the basic properties of the derivative, including the constant, linear, product, quotient, and chain rules, are satisfied.
In what follows, the basic concepts and properties of the used derivatives are listed.
As mentioned before, the linearity, product, quotient, and chain rules are verified in the MRLD case.In addition, we have Definition 2. (Beta derivative (βD)) [46] For µ ∈ (0, 1], the µth-order β operator L is defined as The βD satisfies the basic properties of the integer-order derivative.Moreover, the formula in Equation (68) and provided that f (t) is differentiable, are verified.
The MTD satisfies the main derivative properties.In addition, provided that f (t) is differentiable, is verified.The time-fractional version of the governed liquid crystal model in Equations ( 1)-( 3) is written as i L µ t (q) + α 1 ∂ 2 q ∂x 2 + α 2 p q = 0, (74) We consider the formal solutions in Equations ( 8) and ( 9) with traveling-wave transforms η = β x − ν t µ Γ(µ + 1) , (76) and regarding the MRL, β, and MT fractional derivatives, respectively.The rest of the solution steps are performed similar to that described in Section 3. The obtained nematicons will be as above with η replaced with its value in Equations ( 76)-(78), corresponding to each fractional derivative case.In Figure 9, the 2D profiles of the nematicons in Figure 1 with the effects of various time-fractional operators, the MRLD, βD, and MTD, in comparison to the first derivative (black solid line) are shown.
Consequently, Figures 10-12 present the time-fractional effects on the behaviors of the wave profile q(x, t) and tilt angle p(x, t) that are illustrated in Figure 3, Figure 5, and Figure 7, respectively.65) and (66) at t = 0.5.

Figure 4 .
Figure 4. Contour plots of the (a) wave profile and (b) tilt angle in Figure 3.

Figure 6 .
Figure 6.Contour plots of the (a) wave profile and (b) tilt angle in Figure 5.

Figure 8 .
Figure 8. Contour plots of the (a) wave profile and (b) tilt angle in Figure 7.