Statistical solution and piecewise Liouville theorem for the impulsive discrete Zakharov equations

1.   Introduction

The derivatives of arbitrary order, where the integer-order differentiation and n-fold integration are unified and generalized, the theory of integrals is known as Fractional Calculus. Abel in 1823, described the first application of derivative of order. In 1965, the fractional derivative (FD) was introduced by Lebinitz as a generalization of the Integral order derivative. Later it was reconsidered by Euler, Abel, Riemann Liouville, Grunwald and Letnikov. In several field of research, the topic of fractional calculus plays a major role in the real world problems. Among the diverse fields of science, the fractional calculus has the great application, say in physics, thermodynamics, viscoelasticity, biology, control theory, electrochemistry ^[1,2,3] and acts a controller model for population growth, practically etc. In recent years, fractional differentiation has been drawing increasing attention in the study of social and physical behaviors where scaling power law of fractional order appears universal as an empirical description of such complex phenomena^[4]. Over the past few decades, numerous analyses, real-world problems, and numerical methods were resolved by fractional derivatives and integrals^[5]. Fractional operators used to illustrate better the reality of real-world phenomena with the hereditary property^[6]. The determination of curve's shape such that the time of descent of non-friction point mass sliding down through the curve under the action of gravity is not independent of the starting point is dealt with by the tautochrome problem in which it is the relation of application of fractional calculus with the solution of integral equations. Fractional derivative are amazing tool for illustrating the memory and hereditary properties of diverse materials and processes. On finding the approximation of the solution of the system, numerical analysis have played a major role^[7].

The Caputo fractional differential operator is introduced by the Italian mathematician Caputo in 1967^[8]. Some problems of visco-elasticity are formulated and solved by M.Caputo ^[9] with his own definitions of fractional differentiation. The relationship between the Caputo Fractional Derivative(CFD) with Riemann-Liouville(RL) fractional derivative, Atangana-Baleanu(AB) fractional derivative has been very strong which describes the generalized Mittag-Leffler(ML) functions among their kernels, using certain mathematical model to obtain the results in betterment^[10].

Neutral differential equations occur when max$\{n_1, n_2, ......, n_k\} = n$. The past and present values of the function is dependent by the neutral differential equations, which is similar to retarded differential equations, but the neutral differential equations also depends on derivatives with delays. Neutral type differential equations ^[11,12,13] acts as a model for elastic network arise in high speed computers. That is, for the use of interconnection of switching circuit. Neutral differential equation occur in various branches of applied mathematics, as a result, seeking major heed in recent decades. The development of neutral differential system have been done by many researchers, mentioning the diverse fixed point approaches, mild solutions, and nonlocal conditions^[14,15]. The delay differential equations has the major application in the behaviour of real populations. The systems with impulses are utilized for studying the dynamics of processes subject to abrupt changes at discrete moments ^[16,17].

A set with no empty closed subset of $\mathscr{R}$ is called Time scale. The time scale has been introduced to federate and enhance the theory of differential equation, and many other defined difference systems. The differential equations on time scales for the existence and uniqueness of IVP has been stated by Hilger with some applications. The union of disjoint closed real intervals on time scales acts as an excellent framework for the study of population. In the last few years, differential equations in time scale is majorly developed (see for example^{[18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]}).

In ^[33], by applying fixed point theorems the authors discussed the existence, uniqueness and stability for the non linear fractional differential equations with non linear integral boundary condition on time scales. In our work, we discuss the existence and uniqueness solution to the neutral functional sequential integro differential equations with Caputo fractional derivative on time scale's $\mathbb{T}$ Cauchy problem,

Here, ${}^{c}{\Delta^{\rho}}$, ${}^{c}{\Delta^{\varrho}}$ are CFD. The given functions are $0 < \rho, \rho < 1,$ $\phi :\mathscr{J}\times C([-\epsilon, 0]_\mathbb{T}, \mathscr{R})\rightarrow \mathscr{R}$ and $\varsigma\in C([-\epsilon, 0]_\mathbb{T}, \mathscr{R})$. The function ${\mathit{p}}$ in $[-\epsilon, T]_\mathbb{T}$ and for ${\nu}\in\mathscr{J}$ and $\epsilon > 0$, the element of $C_\epsilon$ : = $C([-\epsilon, 0]_\mathbb{T}, \mathscr{R})$ and we denote by ${\mathit{p}}_\nu$ as,

2.   Preliminaries

Definition 2.1. ^[8] The CFD of order $\beta$ is defined by, for $\beta > 0, t > 0, \beta, {\mathit{x}}, t\in\mathscr{R}$. The fractional operator is,

Definition 2.2. The mapping $\Sigma$: $\mathbb{T}\rightarrow \mathbb{T}$, defined as $\Sigma(t) = inf\{s\in \mathbb{T}: s > t\}$ having inf $\phi = sup \mathbb{T}$. The forward jump operator is defined as a Time scale with an arbitrary non empty closed subset of $\mathscr{R}$ and is denoted by $\mathbb{T}$. (i.e., $\Sigma$(N) = N if $\mathbb{T}$ has a maximum N).

The mapping $\Omega$: $\mathbb{T}\rightarrow \mathbb{T}$, defined by $\Omega(t) = sup\{s\in \mathbb{T}: s < t\}$ with sup $\phi = inf \mathbb{T}$. (i.e., $\Omega$(N) = N if $\mathbb{T}$ has a maximum N) is called the backward jump operator.

Here, the symbol $\phi$ denotes the empty set.

Remark 2.3. In definition 2.2, if $\mathbb{T}$ contains $\bar{n}$ a maximum, then inf $\phi = sup \mathbb{T}$ (i.e., $\Sigma(\bar{n}) = \bar{n}$) and if $\mathbb{T}$ contains $\underline{n}$ a minimum, then sup $\phi = inf \mathbb{T}$ (i.e $\Omega(\underline{n}) = n$), $\phi$ is denoted as the empty set.

Definition 2.4.^[30] A continuous function ${\mathit{u}}$ : $\mathbb{T}\rightarrow \mathscr{R}$ at right dense points is called rd-continuous and in left-dense points left sided limit exists.

Definition 2.5.^[33] A function U: [${\mathit{x}}, {\mathit{y}}$]$\rightarrow \mathscr{R}$ is known as $\Delta$ anti-derivative of function ${\mathit{u}}$: [${\mathit{x}}, {\mathit{y}}$]$\rightarrow \mathscr{R}$, where U is continuous on [${\mathit{x}}, {\mathit{y}}$], $\Delta$ differentiable on (${\mathit{x}}, {\mathit{y}}$), and U$^{\Delta}$(t) = ${\mathit{u}}$(t) $\forall$ t$\in$(${\mathit{x}}, {\mathit{y}}$), where [${\mathit{x}}, {\mathit{y}}$] is the closed bounded interval in $\mathbb{T}$.

The $\Delta$-integral of ${\mathit{u}}$ from ${\mathit{x}}$ to ${\mathit{y}}$ is,

Definition 2.6.^[34] The fractional integral of order $\beta$ of ${\mathit{g}}$ is,

where $\Gamma$ is the Gamma function, $\mathbb{T}$ be a time scale with interval [${\mathit{x}}, {\mathit{y}}$] of $\mathbb{T}$, and ${\mathit{g}}$ be an integrable function on [${\mathit{x}}, {\mathit{y}}$].

Definition 2.7.^[35] Suppose that $\mathbb{T}$ is a time scale. The CFD of order $\beta$ of ${\mathit{g}}$ is defined by

where m = [$\beta$]+1 and [$\beta$] denotes the integer part of $\beta$.

Theorem 2.8.^[36] Let $\beta$, $\beta^{'}$ and ${\mathit{g}}$ be an integrable function on [${\mathit{x}}, {\mathit{y}}$], then,

Lemma 2.9. (Nonlinear alternative for Leray-Schauder type)^[37,38,39]. Consider $U$ a mapping from $\bar{F}$ to $C$, where $C$ is a closed, convex subset of $E$, which is a Banach space and $F$, an open subset of $C$ with 0$\in$F. Then either of the conditions hold:

(ⅰ) $\bar{F}$ has a fixed point,

(ⅱ) ${\mathit{f}} = \Lambda F({\mathit{f}})$, ${\mathit{f}}\in\partial F$ and $\Lambda \in (0, 1)$.

Lemma 2.10. (Krasnoselskii fixed point theorem)^[37]. Let $W$ be closed, bounded, convex and non-empty subset of a Banach space $S$. Let M, N be the operators such that (a) $M{\mathit{u}} +N{\mathit{v}} \in W$ whenever ${\mathit{u}}, {\mathit{v}} \in W$; (b) M is continuous and compact; (c) N is a contraction mapping. Then there exists $z\in W$ such that $z = Mz + Nz$.

Lemma 2.11.^[20] Let $\mathbb{T}$ be a time scale and ${\mathit{g}}$ be an increasing continuous function on the interval $[{\mathit{x}}, {\mathit{y}}]$with the time-scale. Let $G$ be the extension of ${\mathit{g}}$ to the real interval $[{\mathit{x}}, {\mathit{y}}]$, then

then,

3.   Main sequels

We need the following assumptions

(A1) There exists $\Lambda > 0$ such that, $|\psi(\nu, {\mathit{z}})-\psi(\nu, {\mathit{\bar{z}}})| \leq \Lambda||{\mathit{z}}-{\mathit{\bar{z}}}||_{c}$ for $\nu \in \mathscr{J}$ and every ${\mathit{z}}, {\mathit{\bar{z}}} \in C_{\epsilon}$.

(A2) There exists a non '- ve\textquoteright constant $\kappa$ such that, $|\phi(\nu, {\mathit{z}})-\phi(\nu, {\mathit{\bar{z}}})| \leq \kappa||{\mathit{z}}-{\mathit{\bar{z}}}||_{c}$ for $\nu \in \mathscr{J}$ and every ${\mathit{z}}, {\mathit{\bar{z}}} \in C_{\epsilon}$.

(A3) The function $\psi, \phi$ are continuous.

(A4) There exists a continuous non-decreasing function $\psi:[0, \infty]_{\mathbb{T}}\rightarrow (0, \infty)$ and a function $\varUpsilon\in C(\mathscr{J}, \mathscr{J}^{+})$ such that, $|\psi(\nu, {\mathit{z}})| \leq \varUpsilon(\nu)\psi({\mathit{||z||}}_{c})$ for each $(\nu, {\mathit{z}})\in \mathscr{J}\times C_{\epsilon}$.

(A5) There exists a constant ${\bf{L}} > 0$ such that,

where, ${\bf{H}} = \Gamma(\varrho+1){||\varsigma||}_{c} +[|\phi|+\kappa{||\varsigma||}_{c}+2\Phi_{o}]T^{\varrho}$

(A6) There exists a continuous function $\psi, \phi:\mathscr{J}\times C_{\epsilon}\rightarrow \mathscr{R}$, For each

(A7) For each ${\mathit{p}}_{s}, {\mathit{q}}_{s}\in\mathscr{J}$ such that

Remark 3.1. By (A2) for each $(\nu, {\mathit{z}}) \in \mathscr{J}$,

where $\Phi_{o} = {sup_{s\in[0, T]}}_{\mathbb{T}}|\Phi(s, 0)|.$

Remark 3.2. A function ${\mathit{p}} \in$ ${\rm{\tilde E}}$ is said to be a solution of the problem

If ${\mathit{p}}$ satisfies the equation ${}^{c}{\Delta^{\rho}}\Big[{}^{c}{\Delta^{\varrho}}{\mathit{p}}({\nu})-{\Phi\Big({\nu}, {\mathit{p}}_{\nu}, \int_{0}^{t}k_{1}(t, s, {\mathit{p}}_{s})ds\Big)}\Big] = {\psi\Big({\nu}, {\mathit{p}}_{\nu}, \int_{0}^{t}k_{2}(t, s, {\mathit{p}}_{s})ds\Big)}$ on $\mathscr{J}$, the condition ${\mathit{p}}({\nu}) = {\varsigma({\nu})}$ on $[-{\epsilon, 0}]_\mathbb{T}$ and ${}^{c}{\Delta^{\varrho}}{\mathit{p}}(0) = {\phi}$.

Theorem 3.3. The function ${\mathit{p}}\in$ ${\rm{\tilde E}}$ is the solution of the problem.

If

Proof. Using,

we get,

where $\beta \in \mathscr{R}$.

$D^{\varrho}{\mathit{p}}(0) = \phi$, $\beta = \phi - \Phi(0, \varsigma(0))$, is given as,

Thus,

we find $\beta = 0$ and hence the proof.

Theorem 3.4. Assume that (A1) and (A2) holds if,

then there exists a unique element for Initial value problem (1.1) in ${\rm{\tilde E}}$. Proof:. Choose $E: C([-\epsilon, T]_{\mathbb{T}}, \mathscr{R})\rightarrow C([-\epsilon, T]_{\mathbb{T}}, \mathscr{R})$ by,

Let ${\mathit{p}}, {\mathit{q}} \in C([0, T]_{\mathbb{T}}, \mathscr{R})$. Then by (A1) and (A2) we get,

By Lemma 2.11 we get,

Thus,

The operator is contraction $E$. Hence, in problem (1.1) By using Banach contraction principle $E$ is a unique solution on $[-\epsilon, T]_{\mathbb{T}}$ and it has a unique fixed point.

Theorem 3.5. Assume the hypothesis (A2)–(A5) hold. If $\dfrac{\kappa T^{\varrho}}{\Gamma(\varrho+1)} < 1$, then the IVP (1.1) must contain atleast one solution on ${\rm{\tilde E}}$.

Proof. Let us prove the operator $E:C(-\epsilon, T)_{\mathbb{T}, \mathscr{R}}\rightarrow C(-\epsilon, T)_{\mathbb{T}, \mathscr{R}}$ is continuous and completely continuous.

Claim 1: Let $U_{n}\rightarrow y$ in $C(-\epsilon, T)_{\mathbb{T}, \mathscr{R}}$ then $E$ is continuous, where $\{U_{n}\}$ is the sequence. Then,

By Lemma 2.11 we get,

Since $\psi$ and $\Phi$ are continuous functions,

as $n\rightarrow \infty$.

Claim 2: E maps bounded sets into bounded sets in $C({[-\epsilon, T]}_{\mathbb{T}}, \mathscr{R})$. It is necessary to prove for that $k > 0$, and '+ ve\textquoteright constant $\tilde{\Lambda}$ such that, for every ${\mathit{p}} \in {B}_{k} = {{{\mathit{p}} \in C({[\epsilon, T]}_{\mathbb{T}}, \mathscr{R})}:{||{\mathit{p}}||}_{\infty}\leq k}$. we have, ${||E({\mathit{p}})||}_{\infty}\leq \tilde{\Lambda}$.

Using (A4) and (A5), for every $\nu\in\mathscr{J}$, we have,

By Lemma 2.11 we get,

Thus,

Claim 3: E is mapped from bounded to equicontinuous sets of $C({[\epsilon, T]}_{\mathbb{T}}, \mathscr{R})$. $\nu_{1}, \nu_{2} \in \mathscr{J}, \nu_{1} < \nu_{2}, B_{k}$ is the bounded set of $C({[\epsilon, T]}_{\mathbb{T}}, \mathscr{R})$.Let ${\mathit{p}} \in B_{k}$.

By Lemma 2.11 we get,

Thus,

As $\nu_{1}\rightarrow \nu_{2}$, the R.H.S of the above inequality $\rightarrow$ 0. The equicontinuity for the case $\nu_{1}\leq \nu_{2}\leq 0$ & $\nu_{1}\leq 0 \leq \nu_{2}$ is obvious. By using Arzela Ascoli theorem, $C({[-\epsilon, T]}_{\mathbb{T}}, \mathscr{R})\rightarrow C({[-\epsilon, T]}_{\mathbb{T}}, \mathscr{R})$ is continuous and completely continuous.

Claim 4: A set $U \subseteq C({[-\epsilon, T]}_{\mathbb{T}}, \mathscr{R})$ with ${\mathit{p}} \neq {\mathit{m}}\, E({\mathit{p}})$ and for ${\mathit{m}} \in (0, 1)$and ${\mathit{p}} \in \partial {\mathit{p}}$.

Let $U \in C({[-\epsilon, T]}_{\mathbb{T}}, \mathscr{R})$ & ${\mathit{p}} = {\mathit{m}}E({\mathit{p}})$ for some $0 < {\mathit{m}} < 1$. Then, for each $\nu \in \mathscr{J}$, we have,

Considering, for every $\nu \in \mathscr{J}$, we obtain,

Thus,

which can be expressed as

There exists ${\bf{L}}$ and such that ${{||{\mathit{p}}||}_{-\epsilon, T}}_{\mathbb{T}} \neq L$, Setting for ${\mathit{p}}\in \{{\mathit{p}} \in C([-\epsilon, T]_{\mathbb{T}}, \mathscr{R}):{||{\mathit{p}}||_{[-\epsilon, T]}}_{\mathbb{T}}\leq {\bf{L}}\}$. From ${\mathit{p}}$, there is no ${\mathit{p}} \in \partial {\mathit{p}}$ such that ${\mathit{p}} = {\mathit{m}}E {\mathit{p}}$ there exists ${\mathit{m}}\in(0, 1)$. E has a fixed point $u \in \bar{U}$ by using the nonlinear alternative for Leary-Schauder type as the solution of the problem (1.1).

Theorem 3.6. Consider that (A2)–(A3) and $\dfrac{\kappa{T}^{\varrho}}{\Gamma(\varrho+1)} < 1$ holds, and (A6) $|\psi(\nu, {\mathit{x}})|\leq \mathcal{X}_{1}(\nu), |\Phi(\nu, {\mathit{x}})|\leq \mathcal{X}_{2}(\nu)$ for all $({\mathit{p}}, {\mathit{x}}) \in \mathscr{J}\times\mathscr{R},$ where $\mathcal{X}_{1}, \mathcal{X}_{2} \in (\mathscr{J}, \mathscr{R}^{+})$. On $\tilde{E}$ the problem (1.1) has atleast one solution defined.

Proof. The operator $\mathcal{Q}_{1}$ and $\mathcal{Q}_{2}:$

Put

and

and define $D_{\omega} = {{\mathit{p}} \in C([\epsilon, T]_{\mathbb{T}}, \mathscr{R}) = ||{\mathit{p}}||_{\infty}\leq \omega}$

Claim 1: Any sort of ${\mathit{p}}, {\mathit{q}} \in D_{\omega} : \mathcal{Q}_{1}{\mathit{p}} +\mathcal{Q}_{2}{\mathit{q}} \in D_{\omega}$: Any sort of ${\mathit{p}}, {\mathit{q}} \in D_{\omega}$. From (2.5), (2.6) and Lemma 2.11, we have,

This shows that $\mathcal{Q}_{1} y+ \mathcal{Q}_{2} z = D_{\omega}$.

Claim 2: On $D_\omega$, $\mathcal{Q}_{1}$ is a contraction mapping:

Let ${\mathit{p}}, {\mathit{q}}\in D_\omega$. By (2.5) and Lemma 2.11 we have,

Thus,

Claim 3: $\mathcal{Q}_{2}$ and $\psi$ is continuous so $\mathcal{Q}_{2}$ is continuous. $\mathcal{Q}_{2}$ is uniformly bounded on $D_{\omega}$. i.e.,

Claim 4: $\mathcal{Q}_{2}$ is equicontinuous. Defining

For ${\nu}_{1}, {\nu}_{2} \in {[0, T]}_{\mathbb{T}}, {\nu}_{1}, {\nu}_{2}$ by (2.6) and Lemma 2.11 we have,

As $\nu_{1}\rightarrow \nu_{2}$ the R.H.S of the above inequality $\rightarrow$ 0. Therefore $\mathcal{Q}_{2}$ is equicontinuous, where on $D_{\omega}$ $\mathcal{Q}_{2}$ is relatively compact. Therefore, by Arzela-Ascoli theorem $\mathcal{Q}_{2}$ is compact on $D_{\omega}$.

The problem (1.1) has atleast one solution ${[-\epsilon, T]}_{\mathbb{T}}$.

4.   Example

Let us consider the fractional functional integro-differential equation on a Time scale T.

Let

For ${\mathit{p}}, {\mathit{q}} \in [0, \infty]$ and $\nu \in \mathscr{J}$, we have

and

And put $k_{1} = k_{2},$

Hence the conditions (A1) and (A2) hold with $\Lambda = \kappa$ = $\frac{1}{500}$. Thus the problem has a unique solution on $[-\in, 1]_\mathbb{T}$.

5.   Conclusions

In this work, we obtain the existence and uniqueness solution to the integro differential equations for the Caputo fractional derivative on Time scale. The solution of the neutral fractional differential equations along the finite delay condition is derived by using the fixed point theory.

In future we look forward more on circuit analysis, in particular by using Sequential fractional order Neutral functional Integro differential equations on time scales with Caputo fractional operator over Banach spaces.

Acknowledgments

The authors extend their appreciation to the Deputyship for Research & Innovation, Ministry of Education in Saudi Arabia for funding this research work through the project number (IF2/PSAU/2022/01/21626).

Conflict of interest

The authors declare that they have no conflicts of interest.