Investigation of a Spatio-Temporal Fractal Fractional Coupled Hirota System

: This article aims to examine the nonlinear excitations in a coupled Hirota system described by the fractal fractional order derivative. By using the Laplace transform with Adomian decomposition (LADM), the numerical solution for the considered system is derived. It has been shown that the suggested technique offers a systematic and effective method to solve complex nonlinear systems. Employing the Banach contraction theorem, it is confirmed that the LADM leads to a convergent solution. The numerical analysis of the solutions demonstrates the confinement of the carrier wave and the presence of confined wave packets. The dispersion nonlinear parameter reduction equally influences the wave amplitude and spatial width. The localized internal oscillations in the solitary waves decreased the wave collapsing effect at comparatively small dispersion. Furthermore, it is also shown that the amplitude of the solitary wave solution increases by reducing the fractal derivative. It is evident that decreasing the order α modifies the nature of the solitary wave solutions and marginally decreases the amplitude. The numerical and approximation solutions correspond effectively for specific values of time (t). However, when the fractal or fractional derivative is set to one by increasing time, the wave amplitude increases. The absolute error analysis between the obtained series solutions and the accurate solutions are also presented.


Introduction
In recent decades, due to the broad applications of soliton theory in physics, mathematics, and other areas of engineering and applied sciences, the analysis of explicit accurate solutions in the form of solitary wave solutions of evolution equations has played a significant role [1][2][3].In soliton theory, several analytical methods can be applied to calculate the approximate solution to nonlinear partial differential equations [4][5][6].To investigate fractal fractional nonlinear coupled Hirota systems, we consider the Hirota equation [7] ∂ψ(x, t) ∂t where the function ψ(x, t) is a complex-valued function, and δ, γ ∈ R + .The Hirota equation is a modified nonlinear Schrödinger equation (NLSE) that considers higher-order dispersion and time-delay corrections to the cubic nonlinearity [8].It can be observed as a generalization for the NLS equation when the higher dispersion and time-delay changes are considered, as well as complex generalization of the mKdV equation [9].It is also known as soliton equation because this usually has an exact solution where localized moving excitations that preserve their shape in the evolution, similarly to the NLS soliton [10,11].The Hirota equation was initially studied to illustrate the ultra-short pulses experienced from the self-steepening effect and high-order dispersion [8].It has numerous applications, such as localized wave structures, transmission of optical pulses in optical waveguide arrays, rogue waves, solitary waves, and breathers [12,13].It has also been widely studied for rational solitary and rogue wave solutions, interaction properties of complex modified solitons, multi-solitons, dark-bright-rogue wave phenomena, and rogue wave pairs [14][15][16][17].
The Hirota equation has been analytically investigated using different methods containing the Hamiltonian formalism, perturbation theory, Darboux transformation, sine-cosine, and tangent methods [18][19][20][21].The Hirota Equation ( 1) possess more than one parameter solutions in the form solitons, analogous to NLS equation subjected to velocity and amplitude as where A = k 2γ/δ is the amplitude and µ = γ(η 2 − 3k 2 ) represents the wave speed.
The frequency of the phase is given by ω = γ(3η It should be noted that converting the Hirota equation to a coupled Hirota system allows for the study of growth or decay in the system (the real and imaginary parts of the function ψ(x, t)) separately, leading to deeper insights into the solitary wave solutions.Similar to Equation (1), the coupled Hirota system has also been investigated to describe propagation of electromagnetic pulses in optical fibres by employing the inverse scattering transform [24].Several other interesting results have been examined for coupled Hirota system, such as relation amongst dark-bright solitons, Lax pair, Painlevé analyses and rogue waves solutions [25][26][27].Here, we use the concept of fractional fractal derivative of order (α, β) to study the coupled Hirota system where the classical time partial derivatives are replaced with fractional derivatives of the functions of order α (0 < α ≤ 1) in the Caputo sense with power law kernel and fractal derivative β (0 < β ≤ 1).The initial conditions considered herein are The concept of fractional calculus (FC) originated from the challenge that even the concept of a derivative of a classical order could potentially be expanded to stay authentic when the order was not an integer.Following this exceptional theory, the subject of FC caught great attention in mathematics and other areas of scientific research.Since then, FC has been extensively used to examine various physical facts, such as electromagnetic, viscoelastic, and damping theories, artificial intelligence, fluid dynamics, wave propagation theories, chaotic dynamics, heat transfer analysis, appliances, control systems, and various inherent algorithms [28][29][30][31][32][33][34].
In FC, the classical integrals and derivatives are expanded to fractional order as the standard operators are inadequate for understanding many complex systems, because fractional operators produce more realistic and precise results.Numerous fractional derivatives have been derived with distinct kinds of kernels, such as Hilfer, Riemann-Liouville (RL), Caputo, Grünwald, Caputo-Fabrizio, and Atangana-Baleanu derivatives [35,36].Similarly, several novel kinds of fractional operators have proposed to link the ideas of Caputo and other fractional operators [37,38].
In addition to the fractional-order operators, another new theory of differentiation has been suggested, where the derivative is represented as the fractional order with fractal dimension [39].The fundamental concept is called fractal fractional integrals and derivatives [40].If the system under investigation is derivable, the fractal derivative becomes βt β−1 , and the classical operator is regained when the fractal order reaches one [41].In a fractal derivative, the variable is ascended, corresponding to t α .This modified derivative was proposed to model certain real-world physical systems for which conventional physical phenomena are not appropriate and cannot be applied to the non-integral means of fractal dimension [42,43].The theory of fractal fractional integrals and derivatives has been extensively used in engineering applications (fluid motion, aquifer, and crystallization), chaotic processes, and electronic networks [44,45].Here, we will investigate the localized excitations in a coupled Hirota system with time fractional derive in Caputo's sense with fractal dimension.The Laplace transform with Adomian decomposition is a very effective analytical technique, and has been successfully applied to investigate numerous mathematical problems in classical as well as in fractional calculus.We will derive the series solution of the considered system in a systematic manner using LADM.
The article is organized as follows: In Section 2, various fundamental definitions of fractional calculus are included.In Section 3, the general series solution of the coupled Hirota system with order (α, β) is calculated using LADM.The convergence analysis of LADM is examined in Section 4. For confirmation of the obtained findings, a particular example is studied in Section 5 with absolute error analysis.The systematic approximate solutions are compared with numerical results with comprehensive discussions in Section 6.

Preliminaries
Here, we state fundamental definitions associated with fractal fractional calculus and the Laplace and inverse Laplace transforms that will be applied in the next section.Definition 1.Let ψ ∈ H 1 , and α ∈ (0, 1]; we define the Caputo's derivative as , with fractal order β, then the fractal fractional derivative of order α in Caputo sense with power law kernel is given by The above operator can also be represented as where d γ ψ(s) Let ψ(t) be a continuous function on the interval I, then the fractional integral of ψ(t) with fractional order α can be described From the above, the fractal fractional integral can be obtained as [40] FFP t Definition 4. For a function ψ(t), the Laplace transform is symbolized by L and is defined as [46] L Let ψ(t) be a continuous function; the fractal Laplace transform of order α can be written as Definition 6.For the function ψ(x, t), applying Laplace transform to Caputo fractional derivative, we obtain

General Solution of Fractal Fractional Coupled Hirota System
Here, we discuss properties of Laplace transform with order (α, β) and calculate coupled the Hirota system with fractal fractional derivative of order α using the Laplace Adomian decomposition method.In 1980, George Adomian established a novel approach to solving nonlinear functional equations.Since then, this technique has been called the Adomian decomposition method (ADM), and has been applied in numerous studies [47,48].The ADM involves breaking down the equation under consideration into linear and nonlinear parts and generates the series solution, whose terms are determined by a recursive relation using the Adomian polynomials.Numerous fundamental studies on various aspects of the extended version of Adomian's decomposition method have been performed [49][50][51].Similarly, the extended version of the Laplace decomposition method has been used for the solution of different PDEs [52]. Suppose Using the definition presented in [45], we can write Implementing the Laplace transformation to Equation (9) yields This is the Laplace transformation of fractal fractional derivative with power law kernel.
To obtain the general solution of system ( 5) and ( 6), one can apply definition of fractal derivative, which gives Applying the Laplace transformation Using the definition above gives Now, consider ϕ and φ in the series form the non-linear terms are decomposed as The values of A n , B n C n , and D n are known as the Adomian polynomials [53], obtained in the form which gives Applying L −1 to Equations ( 13) and ( 14), we obtain which give , The other terms can be calculated in the similar fashion.

Convergence Analysis
Here, the convergence of the approximate solution obtained through LADM is studied employing the following theorem.

Theorem 1. Let X and Y be Banach spaces, and T
then by Banach contraction theorem, the approximate solutions ( 17) and ( 18) are convergent.
Proof.Let the series ϕ(x, t) = ∑ ∞ n=0 ϕ n obtained by LADM in the form further, we consider that Proof.Mathematical induction can be applied to prove (1) and ( 2).
(1). for n = 1, we have Let the result be true for n − 1, then Similarly, the result can be proven for φ(x, t) = ∑ ∞ n=0 φ n .

Applications
In this section, the numerical examples on the coupled Hirota system is studied with detailed error analysis.

Further, inserting the values gives
In a similar way, other values can be obtained.The complete result can be expressed as Equations ( 22) and ( 23) combined gives the final solution Here, the coupled Hirota system with order (α, β) in the Caputo's sense with power law kernel is calculated via the combination of Adomian decomposition and Laplace transform.To demonstrate the solutions, a specific illustration is considered and obtained the solutions in the form of a series.The convergence of the suggested method also is examined.

Absolute Error Analysis
The error analysis among the exact solutions (2) against solutions ( 22) and ( 23) is provided in the table below.
From Table 1, it follows that the absolute error between accurate and approximately obtained solutions of the Hirota equation reduces when spatial variable increases for particular value of time (t).It is observed that combining in iterations reduces the absolute error.
The numerical solution (24), together with the absolute value of the exact solution (2), are displayed versus the spatial variable x by the dashed and solid curves, as shown in Figure 1a.A pulse-shaped wave profile results from the wave packet's internal oscillations dispersing.The comparison of both curves shows the confirmation of numerical solution.The three-dimensional description to numerical solution (24) is illustrated in Figure 1b against the temporal and spatial variables at δ = 0.25, γ = 1, and η = 1.It reveals that solitary potentials are spatially localized, keeping a stable amplitude for a short time.
The solid curve in left panel of Figure 2 illustrates the real part of Equation (2) versus 10 ≤ x ≤ −10 at time (t) = 0.1.The dashed curved line signifies the real function of the solution (22) for Equation (24) attained by the proposed method (LADM).One can see that the numerical solution accurately meets the obtained systematic approximation.In the right panel of Figure 2, the imaginary value of Equation ( 2) is compared with the numerical result obtained in Equation ( 23) versus 10 ≤ x ≤ −10.It shows the excitations of perturbations, keeping the carrier wave, and as a result establishes a wave packet.The surface plots of analytical solutions ( 22) and ( 23) are depicted in Figure 3.One can see that the internal fluctuations of the waves become constant and sequentially localize the solitary wave.It is important to note that the soliton excitations effect caused by a stability of steepening and the spreading impacts of the nonlinear terms.The nonlinear wave steepening influences the wave breaking for relatively insignificant wave dispersion.Figure 4 illustrates the influence of the fractal variable β by fixing fractional variable α for time (t = 0.1) for obtained solutions ϕ and φ.Similarly, the influence of the variable α with fixed value of β for obtained solutions ϕ and φ is shown in Figure 5.In conclusion, it is found that decreasing the fractal derivative β causes the amplitude to increase.Similarly, reducing the fractional order α modifies the shape of the solitonic solution and reduces its amplitude.Figures 6 and 7 show the behaviour of the solutions ( 22) and ( 23) with various values of α and β for particular value of the spatial variables (x = 1) versus time (t).It is noted that the solitary waves are in quite good agreement for comparatively particular rate of time (t).It is also found that, for α = β = 1, increasing time (t) rapidly increases the wave propagation.Figure 8 constitutes the surface plots of the absolute error of the solutions obtained by LADM and accurate solutions for ϕ and φ for α = β = 1, respectively.Plot for error analysis reported in Table 1.

Conclusions
The coupled nonlinear system was investigated analytically and numerically with fractal fractional derivative with power law kernel.The suggested method (LADM) is applied and calculated the general series solution.The convergence analysis is also examined.For validation, a particular example is considered, and the obtained results are competed with the numerical simulation with physical interpretations.The confinement of the carrier wave and the presence of confined wave packets are observed from the numerical analysis.It is observed that the dispersion nonlinear parameter reduction equivalently influenced the wave amplitude and spatial width.It is revealed that the amplitude of the solitary wave solution increases by decreasing the fractal dimension.It is also observed that decreasing the fractional order α reduces the nature of the solitary wave solutions and slightly decreases the amplitude.The numerical and approximation solutions resemble effectively for specific values of time (t).However, when the fractal or fractional derivative reaches its maximum value by increasing time, the wave amplitude increases.An error analysis is performed between the obtained series solutions versus the accurate solution.It is found that the error between the exact and approximation solutions is minimized for a sufficiently small value of time (t).
As a future work, it will be interesting to study the fractal fractional Hirota equation and coupled Hirota equations with variable coefficients for high-order solitary wave solution, rogue wave solution and dark-bright solutions using the Mellin transform, Chebyshev collocation scheme and Finite difference method reported in [54,55].

Figure 3 .
Figure 3.The wave profiles of analytical solutions (22) and (23) versus spatial and temporal variables.

Figure 4 .
Figure 4.The influence of stable fractional order α with variation in fractal variable β with time (t = 0.1) for obtained solutions ϕ and φ.

Figure 5 .
Figure 5.The influence of stable fractal order β with variation in fractional order α for the obtained solutions ϕ and φ at time (t).

Figure 6 .
Figure 6.The conduct of the functions ϕ(x, t) and φ(x, t) for distinct estimates of β keeping α stable versus time (t).

Figure 8 .
Figure 8.Plot for error analysis reported in Table1.