On a boundary value problem for fractional Hahn integro-difference equations with four-point fractional integral boundary conditions

1.   Introduction

A generalization of classical calculus that permits the differentiation and integration of any order is known as fractional calculus (FC). FC serves as a suitable tool for modeling systems with memory effects, where past actions impact current behavior. Fractional-order derivatives better capture memory effects in biological systems, diffusion processes, and viscoelastic materials. Among the noteworthy applications are^{[1,2,3,4,5,6,7,8,9,10,11,12,13]}:

● Complex physical systems like viscoelastic materials, diffusion processes, electrical circuits, and fluid dynamics are all being modeled and analyzed using FC. Systems that exhibit anomalous diffusion, non-local behaviors, and memory effects are better described using fractional differential equations.

● In control systems, fractional-order controllers, fractional-order observers, and systems with temporal delays are all being modeled using FC. In some applications, fractional-order controllers perform better, are more robust, and are more stable than integer-order controllers.

● In signal processing, FC is used for applications including image processing, time series analysis, noise reduction, and signal denoising. Long-range dependence, fractal features, and handling non-stationary signals are all benefits of fractional-order filters and transforms.

● FC is used to model physiological processes, biological systems, and medical imaging procedures. Applications include neuronal dynamics modeling, biological signal analysis, drug kinetics, and gene expression.

● Fractional-order damping, nonlinear dynamics, and viscoelastic features of mechanical systems are all being modeled using FC. Applications include vibration analysis, control of flexible structures, and structural health monitoring.

● The dynamics of batteries, electrochemical systems, and materials with unusual transport properties are all being modeled using FC. Mass and charge transport processes are more accurately described by fractional diffusion and reaction-diffusion equations.

● Precipitation patterns, pollutant transfer, and groundwater movement are all being modeled using FC. In hydrological systems, long-term memory effects and anomalous transport processes are better described by fractional-order differential equations.

● In robotics, FC is used to model robot dynamics and to create motion and control algorithms. The benefits of fractional-order controllers include robustness against shocks, trajectory tracking, and adaptive control.

● Within its own domain, FC finds use in the analysis of differential equations involving fractional operators, integral transforms, and special functions. Strong tools for deciphering intricate mathematical and physical systems and phenomena are offered by FC.

● FC is employed to analyze economic time series data, which often exhibits long-range dependence and non-stationary behavior. Fractional processes can capture the memory effects and persistence observed in economic variables over time. Economic models incorporating fractional calculus better capture the dynamics of financial markets, including stock price movements, trading volume, and market volatility.

Researchers utilize FC to find applications in various fields, including the study of synchronization phenomena in complex networks ^{[14,15,16,17]}. Global bipartite synchronization refers to the phenomenon where two distinct groups of nodes in a network synchronize their dynamics while remaining unsynchronized with nodes within their group. This phenomenon is crucial in understanding the collective behavior of complex systems and has applications in various domains such as neuroscience, biology, and engineering. The combination of fractional calculus and global bipartite synchronization allows researchers to model and analyze complex network dynamics with memory effects and long-range dependencies. Fractional-order differential equations capture the non-local interactions and memory effects present in real-world systems, while the concept of global bipartite synchronization provides insights into the emergent collective behavior of networked systems. By integrating FC with the study of global bipartite synchronization, researchers develop more accurate models and algorithms for understanding and controlling the dynamics of complex networks. This interdisciplinary approach has the potential to advance our understanding of synchronization phenomena in networked systems and facilitates the design of robust and efficient communication and control strategies.

Fractional order derivative (FOD) is defined in different ways; some definitions are conformable, Riemann-Liouville (R-L), Caputo-Fabrizio, Atangana-Baleanu, Caputo-Hadamard, and Grunuwald Letnikov ^{[18,19,20,21,22]}. There are various definitions for FOD, in contrast to integer-order derivatives. In general, these definitions do not match each other. The conformable fractional derivative (CFD), which provides a unique description of the FOD in contrast to previous versions, is developed by Khalil et al. ^[23].

Regarding a mapping $\Psi(\tau): [0, \infty[\rightarrow \mathscr{R}$, the CFD is:

where $\mathcal{T}^{\ell}_{\tau}$ is the CFD with respect to time and $u-1 < \ell\leq u, \; \tau > 0$, $u\in\mathcal{N}$ and $\lceil\ell\rceil$, the lowest integer that is equal to or larger than $\ell$. In a certain instance, if $0 < \ell\leq1$, we acquire

If $\Psi(\tau)$ is $\ell$-differentiable in some $(0, \Xi)$, $\Xi > 0$ and $\lim_{\tau\rightarrow0^{+}}\Psi^{\ell}(\tau)$ exists, so we have $\Psi^{(\ell)}(0) = \lim_{\tau\rightarrow0^{+}}\Psi^{(\ell)}$$(\tau)$. The fractional integral in the sense of conformable of a function $\Psi(\tau)$ starting from $\tau_{0}\geq0$ is defined as:

Differential equations (DEs) with a random or stochastic component are known as stochastic differential equations (SDEs). In the natural and social sciences, including physics, chemistry, biology, ecology, economics, finance, and engineering, SDEs are extensively employed to model a wide range of phenomena. These are only a few of the several scientific and engineering domains where SDEs are used ^[24,25,26].

● Physics: Systems that are sensitive to random forces, noise, or temperature changes have their motion described by SDEs. One SDE that simulates the Brownian motion of a particle in a fluid is the Langevin equation.

● Chemistry: The dynamics of chemical reactions involving random events, including collisions, transitions, or catalysis, are captured by SDEs. The intermolecular forces and thermodynamic fluctuations are taken into consideration in the molecular dynamics equation, which explains the motion of individual molecules in a system.

● Biology: Genetic mutations, environmental variations, and population interactions are examples of random events that can affect the behavior of biological systems, and SDEs simulate their behavior. An example of an SDE that depicts the dynamics of predator-prey relationships in an ecosystem is the Lotka-Volterra equation.

● Ecology: SDEs simulate the growth and decline of populations that are affected by stochastic events, such as births, deaths, migrations, or competitions. For example, the logistic equation is a SDE that describes the growth of a population with a carrying capacity.

● Economics: Random shock-driven economic variables, such as interest rates, prices, and exchange rates, are analyzed using SDEs. One SDE that establishes the value of a European option in a financial market is the Black-Scholes equation.

● Finance: Financial assets like stocks, bonds, and derivatives that are susceptible to uncertainty have their risk and return estimated by SDEs. One SDE that simulates the price of a stock that moves in a random walk with drift and volatility is the geometric Brownian motion.

● Engineering: Engineering systems subject to random disturbances, including signal processing, control, or communication systems, have their performance and dependability improved by SDEs. One SDE that approximates a system's state from noisy observations is the Kalman filter.

A DE with a FOD and a stochastic process is known as a fractional stochastic differential equation (FSDE). Systems with long-range dependencies, turbulence, economics, and biology are being modeled by these equations. The idea of FSDEs still faces many questions and challenges that need to be addressed. The existence and uniqueness of solutions, stability and convergence, numerical techniques, and applications to physics, biology, finance, and other domains are a few of the interesting subjects.

Some authors have been actively researching the FSDEs lately. Li and Xu developed exponential stability in the mean square of delay FSDEs ^[27]. This study endeavored to provide additional conditions for the exponential stability of the mean square of the scenarios it investigated. Sadovskii's fixed-point theory was applied by Li and Peng ^[28] and worked on the important concept of the controllability of FSDEs. Cui and Yan ^[29] established results for delay integrodifferential equations with stochastic terms with respect to the existence of solutions in Hilbert spaces with the help of Sadovskii's fixed-point. Specifically, in ^[30], Niu and Xie investigated the existence, uniqueness, and regularity of FSDEs. With the help of the Schauder theorem, Chen and Li ^[31] worked on some characteristics of solutions to FSDEs. Under different assumptions, the authors of the articles ^[32,33] showed that solutions to FSDEs exist. Karczewska and Lizama reported a number of findings in ^[34] on the solutions to FSDEs. Schnaubelt and Veraar ^[35] have worked on path-wise continuous features for solutions to stochastic problems. The stopping time technique was used by Xiao and Wang ^[36] to investigate the stability of FSDEs of the Caputo type. Saifullah et al. ^[37] investigated the existence and uniqueness of solutions for FSDEs with R-L. The authors applied the fixed point theory (FPT) for uniqueness and the Móach condition for the existence of a solution. The writers also present the Hyers-Ulam stability of these equations. Moumen et al. ^[38] subjected FSDEs with nonlocal conditions to an approximate controllability analysis. The authors derived the approximate controllability of nonlinear FSDEs. For the desired results, the FPT for multi-valued operators was employed.

Dynamic systems depend on states other than their present ones from the past; these occurrences are often simulated using pantograph differential equations (PDEs) ^{[39,40,41,42,43,44]}. PDEs, FC, and stochastic calculus are combined in stochastic fractional pantograph differential equations (PFSDEs). They are useful in many domains where systems show memory effects in addition to random fluctuations. Among the noteworthy uses of PFSDEs are ^{[45,46,47,48,49,50]}:

● PFSDEs are used to simulate stochastic volatility and memory effects in financial processes, including option pricing models, asset price dynamics, and risk management techniques. These models aid in the development of successful investing strategies as well as the comprehension of the intricate behavior of financial markets.

● PFSDEs are used in the modeling and analysis of random delay and memory effect communication systems, including packet-switched networks, wireless networks, and data transfer protocols. These models help maximize dependability, reduce latency, and optimize network performance.

● Stochastic dynamical systems, random walks in fractal media, diffusion processes, and other physical systems with random fluctuations and memory effects are among the systems that PFSDEs simulate. These models aid in forecasting the long-term evolution of complex physical systems and in understanding their behavior.

● PFSDEs are used to represent biological systems, including gene expression pathways, neuronal dynamics, and physiological responses to external stimuli, that exhibit both stochastic variability and memory effects. These models aid in the development of strategies for the management and treatment of disease as well as the comprehension of the fundamental mechanisms driving biological processes.

● PFSDEs are useful for simulating stochastic variability and memory effects in environmental processes, including pollution dispersion, ecological systems, and climate models. These models are helpful for forecasting changes in the environment, evaluating the effects of human activity on ecosystems, and formulating conservation plans for the environment.

● PFSDEs are employed in the modeling and analysis of control systems, including time-delay systems, networked control systems, and adaptive control algorithms, that have both stochastic disturbances and memory effects. These models aid in the development of reliable control schemes that successfully manage disturbances and uncertainty.

● PFSDEs are used to create stochastic models for optimization issues, reinforcement learning, and machine learning methods. These models make it possible for learning algorithms to identify long-term dependencies as well as short-term oscillations in the data, resulting in learning algorithms that are more reliable and effective.

The existence and uniqueness of solutions to FSDEs are one of the main subjects covered in mathematics. Knowing whether or not a certain FSDE has a unique solution is essential. One of the most important qualitative ideas in dynamical systems is stability theory. Subsequently, the concept of stability garners increased interest across several fields of study and practical uses. In mathematics, Ulam-Hyers stability refers to the approximations of functional equation solutions. In the 1940s, Ulam and Hyers introduced it, and Rassias and others later generalized it.

Inspired by these findings, we first demonstrated the existence and uniqueness of solutions of PFSDEs in the framework of CFD by using the Banach fixed point theorem; after that, we demonstrated the continuous dependency (Con-D) of solutions on both the initial value (In-V) and the fractional exponent $\ell$. The second section was devoted to examining the Ulam-Hyers stability (UHS) of PFSDE using the generalized Gronwall inequalities. The main part of the proof involved the use of the Cauchy-Schwartz inequality (CS-Ineq), It$\hat{o}$ isometry formula (I-IsF), generalized Gronwall's inequality (GG-Ineq), and temporally weighted norm.

Compared with the research results of ^{[39,40,41,42,43,44,45,46,47,48,49,50]}, the major contributions of this paper include at least the following three aspects:

(ⅰ) In contrast to ^{[39,40,41,42,43,44]}, the system we studied was more generalized because it included the stochastic term.

(ⅱ) In contrast to ^{[45,46,47,48,49,50]}, in our research work, we established results regarding the existence, uniqueness, continuous dependency on both the initial condition and fractional component, and Ulam-Hyers stability in the sense of CFD.

(ⅲ) In contrast to ^{[45,46,47,48,49,50,51]}, in which continuous dependency was not discussed, we have evaluated this factor in our research work. We presented continuous dependency of solutions of PFSDEs on the initial value and fractional exponent in the sense of CFD.

We examined the following PFSDEs of order $\frac{1}{2} < \ell < 1$:

where $\eta\tau$ represents the past state and $\eta\in(0, 1)$ and $\ell$ represent the CFD, $\mathcal{A}_{1}: [\tau_{0}, \Im]\times\mathscr{R}^{m}\times\mathscr{R}^{m}\rightarrow \mathscr{R}^{m}$, $\mathcal{A}_{2}:[\tau_{0}, \Im]\times\mathscr{R}^{m}\times\mathscr{R}^{m}\rightarrow \mathscr{R}^{m}$ are measurable and on an underlying complete filtered probability space $(\Omega, \mathscr{F}, \mathfrak{P})$, and $(\mathscr{W}_{\tau})_{\tau\in[0, \; \infty)}$ is a scalar Brownian motion, with the filtration $\mathfrak{F}_{\tau} = \big(\mathscr{F}_{\tau}\big)_{\tau_{0}\leq\tau\leq\Im}$.

Fractional orders in the interval $(\frac{1}{2}, 1)$ are a popular choice ^{[52,53,54,55]}; the precise range may differ based on the situation, the nature of the issue, and the intended trade-off between correctly capturing the dynamics of the system and maintaining mathematical tractability. For instance, the fractional order in viscoelastic materials may indicate the ratio of elasticity to viscosity. Values between $\frac{1}{2}$ and $1$ more closely represent the behavior of the system in such circumstances. When compared to fractional-order systems with orders below $\frac{1}{2}$, those with orders in the interval $(\frac{1}{2}, 1)$ frequently show superior stability and regularity features. Systems having orders smaller than $\frac{1}{2}$ can result in non-unique solutions or a deficiency of solutions in mathematics.

Extending the fractional order range from $(\frac{1}{2}, {1})$ to $(0, 1)$ poses challenges due to the inclusion of the point $0$. However, we have provided unified criteria for selecting fractional orders within the extended range $(0, 1)$ by considering the following factors:

(ⅰ) We established criteria based on stability and regularity properties of solutions. Fractional orders within the range $(0, 1)$ yield stable and well-behaved solutions for the systems under consideration. Analyzing stability and regularity conditions guided the selection of suitable fractional orders.

(ⅱ) Unified criteria ensure the existence and uniqueness of solutions for fractional differential equations within the extended range. We developed mathematical conditions that guarantee the existence and uniqueness of solutions, considering the fractional order's impact on these properties.

(ⅲ) Criteria were based on the physical interpretation of fractional orders within the range $(0, 1)$. We ensured that selected fractional orders aligned with the underlying dynamics of the systems being modeled. This involved considering the physical meaning of fractional orders in specific applications, such as viscoelasticity or diffusion processes.

(ⅳ) Unified criteria address the numerical stability and efficiency of methods used to solve fractional differential equations with fractional orders in the range $(0, 1)$. We established guidelines for selecting numerical methods that accurately and efficiently handle fractional order equations while maintaining stability and convergence properties.

(ⅴ) We conducted comparative analyses between fractional orders within the range $(0, 1)$ and those within the range $(\frac{1}{2}, {1})$. By comparing the behavior and performance of systems modeled with different fractional orders, the authors identified criteria that effectively differentiated between suitable fractional orders.

By considering these factors and establishing unified criteria, we provided guidance for selecting fractional orders within the extended range $(0, 1)$ that ensured the stability, regularity, and physical relevance of solutions in various applications.

The study is structured as follows: We use some fundamental ideas in the next section to provide the framework for the PFSDE conclusions. We first demonstrate the well-posedness of the PFSDE solution in the first subsection of Section 3, and we demonstrate the UHS and offer two instances to support our findings in the second section. Section 4 then presents the conclusion.

2.   Preliminaries

In this section, we will go over some basic concepts, assumptions, characteristics, and definitions that will be useful in this paper.

First, we present some useful characteristics of CFD that make it distinguishable from other fractional derivatives. The primary advantages of CFD are outlined in the following order ^{[56,57,58,59,60,61,62,63,64,65,66,67]}:

(ⅰ) In contrast to the other fractional formulations, it satisfies every need and regulation of an ordinary derivative, such as the product, quotient, and chain rules; Rolle's theorem; and mean-value theorem.

(ⅱ) It reduces the complexity of well-known transforms, such as the Sumudu and Laplace transforms, which are used as instruments to solve some DEs.

(ⅲ) It may be quickly and simply modified to solve systems, analytical, and numerical DEs.

(ⅳ) Fuzzy generalized CFD, Katugampola fractional derivatives, M-CFD, deformable fractional derivatives, and class CFD are just a few of the new concepts that may be created and expanded upon thanks to it.

(ⅴ) In many different applications, it produces new comparisons between CFD and previous fractional definitions.

(ⅵ) CFD draws a lot of interest from researchers since it can be used to depict a wide range of phenomena and applications that need to be addressed.

The following are some significant aspects of CFD that are relevant to differentiation processes.

Theorem 2.1. ^[68] Let $0 < \ell \leq 1$, $\Psi_{1}(\tau)$ and $\Psi_2(\tau)$ be $\ell$-differentiable at a point $\tau > 0$. Then

$(i)$ $\mathcal{T^{\ell}_{\tau}}(\hbar_1 \Psi_1(\tau) + \hbar_2 \Psi_2(\tau)) = \hbar_1 \mathcal{T^{\ell}_{\tau}}\Psi_1(\tau) + \hbar_2 \mathcal{T^{\ell}_{\tau}} \Psi_2(\tau), \; \forall \hbar_1, \; \hbar_2 \in \mathscr{R}.$

$(ii)$ $\mathcal{T^{\ell}_{\tau}} (\tau^\hbar) = \hbar \tau^{\hbar -\ell}, \; \forall \, \hbar \in \mathscr{R}.$

$(iii)$ $\mathcal{T^{\ell}_{\tau}}(\Upsilon) = 0, \; \Upsilon\in \mathscr{R}.$

$(iv)$ $\mathcal{T^{\ell}_{\tau}}\Psi_1(\tau) \Psi_2(\tau)) = \Psi_1(\tau) \mathcal{T^{\ell}_{\tau}}\Psi_2(\tau) + \Psi_2(\tau) \mathcal{T^{\ell}_{\tau}}\Psi_1(\tau).$

$(v)$ $T^{\ell}_{\tau} \dfrac{\Psi_1(\tau)}{\Psi_2(\tau)} = \dfrac{\Psi_2(\tau) \mathcal{T^{\ell}_{\tau}}\Psi_1(\tau)- \Psi_1(\tau) \mathcal{T^{\ell}_{\tau}}\Psi_{2}(\tau)}{\Psi_{2}^{2}(\tau)}.$

Denote by $\big(\Omega, \mathscr{F}, \big(\mathscr{F}_{\tau}\big)_{\tau_{0}\leq\tau\leq\Im}, \mathfrak{P}\big)$ the complete probability space and $\mathscr{W}(\tau)$ is the standard Brownian motion. Suppose $\mathbb{Z}_{\tau} = \mathscr{L}^{2}\big(\Omega, \mathscr{F}_{\tau}, \mathfrak{P}\big)$ be the space of all $\mathscr{F}_{\tau}$-measurable, integrable, and mean square functions $\mathbb{Z} = \big(\mathbb{Z}_{1}, \mathbb{Z}_{2}, \cdots, \mathbb{Z}_{m}\big)^{T}$:$\Omega\rightarrow \mathscr{R}^{m}$ with

making use of the usual Euclidian norm, $\Vert\cdot\Vert$.

Definition 2.2. To prove the desired results, we assume that coefficients $\mathcal{A}_{1}$ and $\mathcal{A}_{2}$ meet the following assumptions:

$(\mathscr{H}_{1}):$ There is $\mathscr{L} > 0$ such that $\forall\big(\tau, \nu_{1}, \nu_{2}, \tilde{\nu}_{1}, \tilde{\nu}_{2}\big)\in[\tau_{0}, \Im]\times \mathscr{R}^{m}\times\mathscr{R}^{m}\times\mathscr{R}^{m}\times\mathscr{R}^{m}$:

$(\mathscr{H}_{2}):$ $\mathcal{A}_{1}(., 0, 0)$ and $\mathcal{A}_{2}(., 0, 0)$ are bounded, i.e,

We give now the definition of UHS.

Definition 2.3. The Eq (1.1) is UHS with respect to $\epsilon$ if there is a constant $\aleph > 0$ such that for each $\epsilon > 0$ and $\mathscr{N}(\tau)\in\mathscr{H}^{2}\left(\left[\eta\tau_{0}, \Im\right]\right)$ of the following inequality: $\forall \tau \in\left[\tau_{0}, \Im\right]$

there exists a solution $\mathbb{Z}(\tau)\in\mathscr{H}^{2}\left(\left[\eta \tau_{0}, \Im\right]\right)$ of Eq (1.1), with $\mathbb{Z}(\tau) = \mathscr{N}(\tau)$ for $\tau\in\left[\eta \tau_{0}, \tau_{0}\right]$, satisfies $\mathscr{E}\big[\Vert\mathscr{N}(\tau)-\mathbb{Z}(\tau)\Vert^{2}\big]\leq\aleph\epsilon$, $\forall \tau \in\left[\tau_{0}, \Im\right]$.

3.   The main results

In this part of the paper, we present our results for the well-posedness solutions of PFSDE.

3.1.   Well-posedness of PFSDEs under the standard Lipschitz condition of coefficients

Assume $\widetilde{\mathscr{H}}^{2}(\eta\tau_{0}, \Im)$ is the space where all processes $\mathbb{Z}(\tau)$ which are measurable $\mathscr{F}_{\Im}-$adapted, with $\mathscr{F}_{\Im} = (\mathfrak{F}_{\tau})_{\tau\in[\eta\tau_{0}, \Im]}$ and satisfy the following:

It is easy to show that $\big(\widetilde{\mathscr{H}}^{2}\big(\big[\eta\tau_{0}, \Im\big]\big), \Vert\cdot\Vert_{\widetilde{\mathscr{H}}^{2}}\big)$ is a Banach space. For $\varpi\in\mathscr{H}^{2}\big(\big[\eta\tau_{0}, \tau_{0}\big]\big)$, we consider the operator $\Lambda_{\varpi}$ : $\widetilde{\mathscr{H}}^{2}\big(\big[\eta\tau_{0}, \Im\big]\big)\rightarrow \widetilde{\mathscr{H}}^{2}\big(\big[\eta\tau_{0}, \Im\big]\big)$ defined by

for $\tau \in\left[\tau_{0}, \Im\right]$ and $\Lambda_{\varpi}\big(\mathbb{Z}(\tau)\big) = \varpi(\tau)$ for $\tau \in\left[\eta \tau_{0}, \tau_{0}\right]$.

The subsequent lemma shows this operator's well-defined property. It is essential to demonstrate the well-definedness of an operator in the context of DEs because a well-defined operator in DEs guarantees a clear, consistent mapping from the domain to the range that is independent of the representation of the domain's elements. In the study of DEs, this ensures the accuracy and consistency of mathematical operations and solutions.

The elementary inequality below is employed in the proof of this result as well as multiple others that follow:

Lemma 3.1. Assume that assumptions $\mathscr{A}_{1}$ and $\mathscr{A}_{2}$ are true, then operator $\Lambda_{\varpi}$ is well-defined for all $\varpi \in \widetilde{\mathscr{H}}^{2}\big(\big[\eta \tau_{0}, \tau_{0}\big]\big)$.

Proof. Let $\mathbb{Z}\in\widetilde{\mathscr{H}}^{2}\big(\big[\eta \tau_{0}, \Im\big]\big)$. By utilizing Eqs (3.1) and (3.2), we get the following result:

With CS-ineq, we derive

Via the use of ($\mathscr{A}_{1}$), we conclude that

Then, from Eq (3.5), we get

By utilizing Eq (3.6) in Eq (3.4), we acquire the following outcomes:

Considering the use of I-IsF, we can determine that

Applying ($\mathscr{A}_{1}$), we may deduce that

Thus, using Eq (3.9), we derive the next outcomes of our findings from Eq (3.8).

From Eqs (3.7) and (3.10) consequently, we extract our intended conclusions.

To establish the EU, we have to prove the following lemma.

Lemma 3.2. For any $\ell > \frac{1}{2}$ and $\tau > 0$, the inequality that follows is valid:

where $\mathbb{E}_{2\ell-1}(.)$ is a Mittag-Leffler function which is defined as

Proof. Let $\mathscr{X} > 0$ be arbitrary. First, we interchange integral and sum, and then we apply the procedures that follow identity.

So, we get

In this case, the beta function is $\mathbb{B}$. This brings the proof to an end.

Theorem 3.3. Assume that the assumptions $(\mathscr{A}_{1})$ and $(\mathscr{A}_{2})$ are fulfilled, then the problem Eq (1.1) with In.C $\varpi$ has unique solution.

Proof. In the beginning, proceed to select a fixed positive constant $\mathscr{X}$:

We establish a weighted norm $\Vert\cdot\Vert_{\mathscr{X}}$ over the space $\widetilde{\mathscr{H}}^{2}([\eta\tau_{0}, \Im])$ as:

$\mathscr{Y}(\tau) = \mathbb{E}_{2\ell-1}\big(\mathscr{X}(\tau-\tau_{0})^{2\ell-1}\big)$ when $\tau\in[\tau_{0}, \Im]$ and $\mathscr{Y}(\tau) = 1$ for $\tau\in[\eta\tau_{0}, \tau_{0}]$. $\Vert\cdot\Vert_{\widetilde{\mathscr{H}}^{2}}$ and $\Vert\cdot\Vert_{\mathscr{X}}$, two norms, are equivalent. For this reason, $\big(\widetilde{\mathscr{H}}^{2}([\eta\tau_{0}, \Im]\big), \Vert\cdot\Vert_{\mathscr{X}})$ is a Banach space. Let $\mathbb{Z}_{1}, \mathbb{Z}_{2} \in\widetilde{\mathscr{H}}^{2}\big(\big[\eta \tau_{0}, \Im\big]\big)$, we have $\forall\tau \in\big[\eta\tau_{0}, \tau_{0}\big], \Lambda_{\varpi} \big(\mathbb{Z}_{1}(\tau)\big)-\Lambda_{\varpi}\big(\mathbb{Z}_{2}(\tau)\big) = 0$. For $\tau \in\left[\tau_{0}, \Im\right]$, we get

Employing the CS-ineq and ($\mathscr{A}_{1}$), we acquire

By applying I-IsF and ($\mathscr{A}_{1}$), we reach the subsequent outcome:

Then,

where $\delta = \frac{4 \mathscr{L}^{2}\left(\Im-\tau_{0}+1\right)}{\Gamma(\ell)^{2}}$. Then,

Therefore,

Then, there is a unique solution to Eq (1.1).

We will demonstrate the Con-D of solutions on In-V that are $\varpi$ and $\Theta$ to PFSDEs in the following theorem.

Theorem 3.4. ($\mathscr{A}_{1}$) and ($\mathscr{A}_{2}$) are assumed to be valid. The solution $\mathscr{Z}_{\ell}(., \varpi)$ of Eq (1.1) is always dependent on $\varpi$, i.e.,

Proof. Because $\mathscr{Z}_{\ell}(., \varpi)$ and $\mathscr{Z}_{\ell}(., \Theta)$ are solutions of Eq (1.1), employing Eq (3.2), the CS-Ineq, I-ISF, $(\mathscr{A}_{1})$, $(\mathscr{A}_{2})$, and Lemma 3.2, we yield

By virtue of definition of $\Vert\cdot\Vert_{\mathscr{X}}$, we have

which, together with Eq (3.13), we derive

We will demonstrate the Con-D of solutions on fractional exponent $\ell$ to PFSDEs in the following theorem.

Theorem 3.5. Let us suppose that both ($\mathscr{A}_{1}$) and ($\mathscr{A}_{2}$) are valid. After that, the solution $\mathscr{Z}_{\ell}(., \varpi)$ to Eq (1.1) is Con-D on $\ell$, i.e.,

Proof. Let $\ell, \tilde{\ell}$ be arbitrary but fixed. As $\mathscr{Z}_{\ell}(., \varpi)$ and $\mathscr{Z}_{\tilde{\ell}}(., \varpi)$ are solutions of Eq (1.1) and we implement Eq (3.2), the CS-Ineq, I-IsF, $(\mathscr{A}_{1})$, $(\mathscr{A}_{2})$, and Lemma 3.2, we achieve

where

By virtue of Lemma 3.2, we attain

However, we yield

As a result,

This suggests the completion of the proof when combined with Eq (1.1).

In the next subsection, we will prove that Eq (1.1) is Ulam-Hyers stable.

3.2.   The Ulam-Hyers stability of PFSDEs

Theorem 3.6. Under ($\mathscr{A}_{1}$) and ($\mathscr{A}_{2}$), the PFSDE Eq (1.1) is UHS on $\left[\tau_{0}, \Im\right]$.

Proof. Set $\epsilon > 0$ and $\mathscr{N}(\tau)\in\mathscr{H}^{2}\left(\left[\eta \tau_{0}, \Im\right]\right)$ satisfies Eq (2.1). Let $\mathbb{Z}(\tau) \in \mathscr{H}^{2}\left(\left[\eta \tau_{0}, \Im\right]\right)$ be the unique solution of Eq (1.1) with In.V $\mathbb{Z}(\tau) = \mathscr{N}(\tau)$ for $\tau\in\left[\eta \tau_{0}, \tau_{0}\right]$, then

Thus,

Thus, using CS-Ineq and $(\mathscr{A}_{1})$, $(\mathscr{A}_{2})$, we can derive that

Let $\digamma(\tau) = \sup _{\chi\in\left[\eta \tau_{0}, \tau\right]} \mathscr{E}\big[\Vert\mathscr{N}(\chi)-\mathbb{Z}(\chi)\Vert^{2}\big]$ for $\tau \in\left[\tau_{0}, \Im\right]$.

We have

and

Then, for $\tau\in\left[\tau_{0}, \Im\right]$, we obtain

Hence, for all $\tau\in\left[\tau_{0}, \tau\right]$,

Then,

for all $\tau\in\left[\tau_{0}, \Im\right]$, where $\delta_{1} = \frac{16\mathscr{L}^{2}\left(\Im-\tau_{0}\right)^{2\ell-1}}{(2\ell-1)}$ and $\delta_{2} = 16 \mathscr{L}^{2}$.

Using GG-Ineq, we get

where $\delta_{3} = 2\mathbb{E}_{2\ell-1}\left(\delta_{2} \Gamma(2\ell-1)\left(\Im-\tau_{0}\right)^{2\ell-1}\right)$ and $\delta_{4} = \delta_{1}\mathbb{E}_{2\ell-1}\left(\delta_{2}\Gamma(2 \ell-1)\left(\Im-\tau_{0}\right)^{2 \ell-1}\right)$.

Using GG-Ineq, we get

Hence,

where $\aleph = \delta_{3} e^{\delta_{4}\left(\Im-\tau_{0}\right)}$. Thus, Eq (1.1) is UHS.

4.   Examples

The purpose of this section is to understand the results that we established in the above section. For this purpose, we present two numerical examples of PFSDEs.

The PFSDEs find applications in various fields where systems exhibit random behavior and are influenced by both deterministic and stochastic factors. One physical application of PFSDEs is in modeling the dynamics of certain types of mechanical systems subject to random perturbations. In such cases, PFSDEs can be used to model the dynamics of the system. The pantograph term arises when there are delays in the system's response, for instance, due to the time it takes for vibrations to propagate through the structure. The stochastic component captures the randomness inherent in the external forcing or system parameters.

A suspension bridge subjected to random wind loads can be modeled using PFSDEs. The deterministic part of the equation accounts for the bridge's structural properties, such as stiffness and damping, while the stochastic part represents the random nature of wind loads. The pantograph term may account for delays in the bridge's response to wind gusts due to factors like the propagation speed of vibrations along its structure.

Example 4.1. Let $\tau\in[1, 2]$ and the PFSDE for each $\epsilon > 0$ be provided by

where

Where, the deterministic and stochastic parts of Eq (4.1) are $e^{0.5 \tau} \cos (\mathbb{Z}(\tau))$ and $\frac{\cos\big(\mathbb{Z}(0.5 \tau)\big)}{\sqrt{1+0.5 \tau}}$ respectively. Now we demonstrate that, with relation to $\epsilon$, Eq (4.1) is UHS. Let $\left(\tau, \mathbb{Z}_{1}, \mathbb{Z}_{2}\right) \in[1, 2] \times \mathscr{R} \times \mathscr{R}$, therefore,

This means that ($\mathscr{A}_{1}$) is satisfied. Furthermore,

and

The assumptions ($\mathscr{A}_{1}$) and ($\mathscr{A}_{2}$) are therefore met. Thus, with respect to $\epsilon$ on [1, 2], Eq (4.1) is UHS utilizing Theorem 3.6. For Systems Eqs (4.1) and (4.2), we conduct a simulation based on the Euler-Maruyama scheme with a step size of $10^{-3}$. We can see from Figure 1 that the distance between $\mathscr{N}(\tau)$ and $\mathbb{Z}(\tau)$ is less than a constant, and according to Definition 2.3, the solution of Eq (4.1) is UHS.

Example 4.2. Let $\tau \in[1, 4]$ and the PFSDE for each $\epsilon > 0$ be provided by

where

Where, the deterministic and stochastic parts of the Eq (4.2) are $\frac{\cos\mathbb{Z}(\tau)}{\sqrt{1+(0.25\tau)^{2}}}$ and $\frac{\mathbb{Z}(0.25 \tau)}{1+\tau^{2}}$ respectively. Now, with regard to $\epsilon$, we demonstrate that Eq (4.2) is UHS. Assume that $\left(\tau, \mathbb{Z}_{1}, \mathbb{Z}_{2}\right) \in[1, 4] \times \mathscr{R} \times \mathscr{R}$, hence

Consequently, ($\mathscr{A}_{1}$) is met. Likewise

and

The assumptions ($\mathscr{A}_{1})$ and ($\mathscr{A}_{2}$) are therefore met. Thus, with respect to $\epsilon$ on [1, 4], Eq (4.2) is UHS employing Theorem 3.6. We can see from Figure 2 that the distance between $\mathscr{N}(\tau)$ and $\mathbb{Z}(\tau)$ is less than a constant, and according to Definition 2.3, the solution of Eq (4.2) is UHS.

Continuous dependency is a key property that contributes to the well-posedness of systems of differential equations, ensuring stability in the solutions by guaranteeing that small changes in parameters or initial conditions lead to correspondingly small changes in the solutions. Therefore, continuous dependency is closely tied to the overall stability and reliability of mathematical models described by differential equations. As both models presented in Examples 4.1 and 4.2 satisfy the criteria of continuous dependency, we conclude that both models are stable.

5.   Conclusions

In the present paper, we examined the existence and uniqueness of the pantograph fractional stochastic differential equations using the Banach fixed point theorem. Second, we demonstrated that solutions continuously depend on the fractional order $\ell$ and initial value. The fractional-order derivative is considered in the conformable sense. Using the generalized Gronwall inequalities and stochastic analytic methods, we examined the Ulam-Hyers stability of the considered problem. We demonstrated our findings with two examples in the final section.

We will use numerical techniques in our future work to solve various kinds of real-world challenges modeled with pantograph fractional stochastic differential equations.

Use of AI tools declaration

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

Acknowledgments

This work was supported by Princess Nourah bint Abdulrahman University Researchers Supporting Project number (PNURSP2024R157), Princess Nourah bint Abdulrahman University, Riyadh, Saudi Arabia. This study was supported via funding from Prince Sattam bin Abdulaziz University project number (PSAU/2024/R/1445). The authors are thankful to the Deanship of Graduate Studies and Scientific Research at the University of Bisha for supporting this work through the Fast-Track Research Support Program.

Conflict of interest

The authors declare that there are no conflicts of interest.