2D vs 3D tracking in bacterial motility analysis

Jacqueline Acres; Jay Nadeau; Jacqueline Acres; Jay Nadeau

doi:10.3934/biophy.2021030

AIMS Biophysics

2021, Volume 8, Issue 4: 385-399. doi: 10.3934/biophy.2021030

Previous Article Next Article

Research article Special Issues

2D vs 3D tracking in bacterial motility analysis

Jacqueline Acres ,
Jay Nadeau ^,

Department of Physics, Portland State University, 1719 SW 10th Ave., Portland, OR 97201, USA

Received: 02 October 2021 Accepted: 24 November 2021 Published: 09 December 2021

Digital holographic microscopy provides the ability to observe throughout a large volume without refocusing. This capability enables simultaneous observations of large numbers of microorganisms swimming in an essentially unconstrained fashion. However, computational tools for tracking large 4D datasets remain lacking. In this paper, we examine the errors introduced by tracking bacterial motion as 2D projections vs. 3D volumes under different circumstances: bacteria free in liquid media and bacteria near a glass surface. We find that while XYZ speeds are generally equal to or larger than XY speeds, they are still within empirical uncertainties. Additionally, when studying dynamic surface behavior, the Z coordinate cannot be neglected.

Keywords:

Citation: Jacqueline Acres, Jay Nadeau. 2D vs 3D tracking in bacterial motility analysis[J]. AIMS Biophysics, 2021, 8(4): 385-399. doi: 10.3934/biophy.2021030

Related Papers:

[1]	Qun Dai, Shidong Liu . Stability of the mixed Caputo fractional integro-differential equation by means of weighted space method. AIMS Mathematics, 2022, 7(2): 2498-2511. doi: 10.3934/math.2022140
[2]	Thabet Abdeljawad, Sabri T. M. Thabet, Imed Kedim, Miguel Vivas-Cortez . On a new structure of multi-term Hilfer fractional impulsive neutral Levin-Nohel integrodifferential system with variable time delay. AIMS Mathematics, 2024, 9(3): 7372-7395. doi: 10.3934/math.2024357
[3]	Choukri Derbazi, Zidane Baitiche, Mohammed S. Abdo, Thabet Abdeljawad . Qualitative analysis of fractional relaxation equation and coupled system with Ψ-Caputo fractional derivative in Banach spaces. AIMS Mathematics, 2021, 6(3): 2486-2509. doi: 10.3934/math.2021151
[4]	J. Vanterler da C. Sousa, E. Capelas de Oliveira, F. G. Rodrigues . Ulam-Hyers stabilities of fractional functional differential equations. AIMS Mathematics, 2020, 5(2): 1346-1358. doi: 10.3934/math.2020092
[5]	Kaihong Zhao, Shuang Ma . Ulam-Hyers-Rassias stability for a class of nonlinear implicit Hadamard fractional integral boundary value problem with impulses. AIMS Mathematics, 2022, 7(2): 3169-3185. doi: 10.3934/math.2022175
[6]	Jiqiang Zhang, Siraj Ul Haq, Akbar Zada, Ioan-Lucian Popa . Stieltjes integral boundary value problem involving a nonlinear multi-term Caputo-type sequential fractional integro-differential equation. AIMS Mathematics, 2023, 8(12): 28413-28434. doi: 10.3934/math.20231454
[7]	Hasanen A. Hammad, Hassen Aydi, Hüseyin Işık, Manuel De la Sen . Existence and stability results for a coupled system of impulsive fractional differential equations with Hadamard fractional derivatives. AIMS Mathematics, 2023, 8(3): 6913-6941. doi: 10.3934/math.2023350
[8]	Dongming Nie, Usman Riaz, Sumbel Begum, Akbar Zada . A coupled system of $p$ -Laplacian implicit fractional differential equations depending on boundary conditions of integral type. AIMS Mathematics, 2023, 8(7): 16417-16445. doi: 10.3934/math.2023839
[9]	Thanin Sitthiwirattham, Rozi Gul, Kamal Shah, Ibrahim Mahariq, Jarunee Soontharanon, Khursheed J. Ansari . Study of implicit-impulsive differential equations involving Caputo-Fabrizio fractional derivative. AIMS Mathematics, 2022, 7(3): 4017-4037. doi: 10.3934/math.2022222
[10]	Subramanian Muthaiah, Dumitru Baleanu, Nandha Gopal Thangaraj . Existence and Hyers-Ulam type stability results for nonlinear coupled system of Caputo-Hadamard type fractional differential equations. AIMS Mathematics, 2021, 6(1): 168-194. doi: 10.3934/math.2021012

Abstract

1. Introduction

Fractional differential equations have played an important role and have presented valuable tools in the modeling of many phenomena in various fields of science and engineering ^{[6,7,8,9,10,11,12,13,14,15,16]}. There has been a significant development in fractional differential equations in recent decades ^{[2,3,4,5,26,23,33,37]}. On the other hand, many authors studied the stability of functional equations and established some types of Ulam stability ^{[1,17,18,19,20,21,22,24,27,28,29,30,31,32,33,34,35,36,37]} and references there in. Moreover, many authors discussed local and global attractivity ^{[8,9,10,11,34]}.

Benchohra et al. ^[13] established some types of Ulam-Hyers stability for an implicit fractional-order differential equation.

A. Baliki et al. ^[11] have given sufficient conditions for existence and attractivity of mild solutions for second order semi-linear functional evolution equation in Banach spaces using Schauder's fixed point theorem.

Benchohra et al. ^[15] studied the existence of mild solutions for a class of impulsive semilinear fractional differential equations with infinite delay and non-instantaneous impulses in Banach spaces. This results are obtained using the technique of measures of noncompactness.

Motivated by these works, in this paper, we investigate the following initial value problem for an implicit fractional-order differential equation

$\begin{align} \left\{ \begin{array}{c} ^C\mathfrak{D}^\alpha \left[ \mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\right] = \mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\mathfrak{\beta} \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))) \; \; \; \mathfrak{t}\in {J}, \; 1 \lt \alpha \leq 2 , \; \alpha \geq \beta, \\\\ \bigg(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(t))\bigg)\bigg|_{\mathfrak{t} = 0} = 0\; \; \; \; \mbox{and}\; \; \; \; \; \frac{d}{d\mathfrak{t}}\big[\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big]_{\mathfrak{t} = 0} = 0\end{array} \right. \end{align}$

(1.1)

where $\; ^C\mathfrak{D}^\alpha\;$ is the Caputo fractional derivative, $\; \mathfrak{h}:J \times R \longrightarrow R, \; \mathfrak{g} _{1}: J \times R \times R\longrightarrow R \;$ and $\mathfrak{g} _{2}:J \times R \longrightarrow R$ are given functions satisfy some conditions and $\; J = [0, T]$ .

we give sufficient conditions for the existence of solutions for a class of initial value problem for an neutral differential equation involving Caputo fractional derivatives. Also, we establish some types of Ulam-Hyers stability for this class of implicit fractional-order differential equation and some applications and particular cases are presented.

Finally, existence of at least one mild solution for this class of implicit fractional-order differential equation on an infinite interval $\; J = [0, +\infty), \;$ by applying Schauder fixed point theorem and proving the attractivity of these mild solutions.

By a solution of the Eq (1.1) we mean that a function $x \in C^2(J, R)$ such that

(i) the function $\mathfrak{t} \to \left[\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\right]\in C^2(J, R)$ and

(ii) $\mathfrak{x}$ satisfies the equation in (1.1).

2. Preliminaries

Definition 1. ^[23] The Riemann-Liouville fractional integral of the function $\mathfrak{f} \in L^1([\mathfrak{a}, \mathfrak{b}])$ of order $\alpha \in R_+$ is defined by

$\mathfrak{I}_{\mathfrak{a}}^{\alpha}\; \mathfrak{f}(\mathfrak{t}) = \int_{\mathfrak{a}}^\mathfrak{t}\; \frac{(\mathfrak{t}-\mathfrak{s})^{\alpha-1}}{\Gamma(\alpha)} \; \mathfrak{f}(\mathfrak{s}) \; d\mathfrak{s}.$

and when $\mathfrak{a} = 0,$ we have $\mathfrak{I}^{\alpha}\; \mathfrak{f}(\mathfrak{t}) = \mathfrak{I}_{0}^{\alpha}\; \mathfrak{f}(\mathfrak{t}).$

Definition 2. ^[23] For a function $\; \mathfrak{f}:[\mathfrak{a}, \mathfrak{b}]\rightarrow R\;$ the Caputo fractional-order derivative of $\; \mathfrak{f}, \;$ is defined by

$^C\mathfrak{D}^\alpha \mathfrak{h}(\mathfrak{t}) = \frac{1}{\Gamma(\mathfrak{n}-\alpha)}\int_\mathfrak{a}^\mathfrak{t} \frac{\mathfrak{h}^{(\mathfrak{n})}(\mathfrak{s})}{(\mathfrak{t}-\mathfrak{s})^{\mathfrak{n}-\alpha-1}}\; d\mathfrak{s},$

where where $\; \mathfrak{n} = [\alpha]+ 1\;$ and $\; [\alpha]\;$ denotes the integer part of the real number $\; \alpha.$

Lemma 1. ^[23]. Let $\; \alpha \geq 0\;$ and $\; \mathfrak{n} = [\alpha]+ 1.\;$ Then

$\mathfrak{I}^\alpha(^C\mathfrak{D}^\alpha \mathfrak{f}(\mathfrak{t})) = \mathfrak{f}(\mathfrak{t})- \sum\limits_{k = 0}^{\mathfrak{n}-1}\frac{\mathfrak{f}^k(0)}{k!}\mathfrak{t}^k$

Lemma 2. Let $\; \mathfrak{f} \in L^1([\mathfrak{a}, \mathfrak{b}])\;$ and $\; \alpha \in (0, 1], \;$ then

(i) $^C\mathfrak{D}^{\alpha}\mathfrak{I}^\alpha \mathfrak{f}(\mathfrak{t}) = \mathfrak{f}(\mathfrak{t}).\; \;$

(ii) The operator $\; \mathfrak{I}^\alpha\;$ maps $\; L^1([\mathfrak{a}, \mathfrak{b}])\;$ into itself continuously.

(iii) For $\; \gamma, \; \beta > 0, \;$ then

$\mathfrak{I}_\mathfrak{a}^\beta \; \mathfrak{I}_\mathfrak{a}^\gamma \; \mathfrak{f}(\mathfrak{t})\; = \; \mathfrak{I}_\mathfrak{a}^\gamma \mathfrak{I}_\mathfrak{a}^\beta \; \mathfrak{f}(\mathfrak{t})\; = \; \mathfrak{I}_\mathfrak{a}^{\gamma+\beta}\; \mathfrak{f}(\mathfrak{t}),$

For further properties of fractional operators (see ^[23,25,26]).

3. Main results

Consider the initial value problem for the implicit fractional-order differential Eq (1.1) under the following assumptions:

(i) $\mathfrak{h}:J \times R \longrightarrow R$ is a continuous function and there exists a positive constant $K_1$ such that:

$\mid \mathfrak{h}(\mathfrak{t}, \mathfrak{x})- \mathfrak{h}(\mathfrak{t}, \mathfrak{y})\mid \leqslant K_1 \mid \mathfrak{x}-\mathfrak{y}\mid \; \; \mbox{ for each }\; \; \mathfrak{t}\in J \; \mbox{ and } \; \mathfrak{x}, \mathfrak{y}\in R.$

(ii) $\mathfrak{g} _{1}: J \times R \times R\longrightarrow R$ is a continuous function and there exist two positive constants $K, H$ such that:

$\mid \mathfrak{g} _{1}(\mathfrak{t}, \mathfrak{x}, \mathfrak{y}) - \mid \mathfrak{g} _{1}(\mathfrak{t}, \widetilde{\mathfrak{x}}, \widetilde{\mathfrak{y}}) \mid \leqslant K \mid \mathfrak{x}- \widetilde{\mathfrak{x}} \mid + H \mid \mathfrak{y}- \widetilde{\mathfrak{y}}\mid \; \mbox{ for each }\; \mathfrak{t}\in J \; \mbox{ and }\; \mathfrak{x}, \widetilde{\mathfrak{x}}, \mathfrak{y}, \widetilde{\mathfrak{y}} \in R$

(iii) $\mathfrak{g }_{2}:J \times R \longrightarrow R$ is a continuous function and there exists a positive constant $K_2$ such that:

$\mid \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x})- \mathfrak{g}_2(\mathfrak{t}, \mathfrak{y})\mid \leqslant K_2 \mid \mathfrak{x}-\mathfrak{y}\mid \; \mbox{ for each } \; \mathfrak{t}\in J \; \mbox{ and} \; \mathfrak{x}, \mathfrak{y}\in R.$

Lemma 3. Let assumptions (i)-(iii) be satisfied. If a function $\mathfrak{x} \in C^2(J, R)$ is a solution of initial value problem for implicit fractional-order differential equation $(1.1)$ , then it is a solution of the following nonlinear fractional integral equation

$\begin{equation} \mathfrak{x}(\mathfrak{t}) = \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s} \end{equation}$

(3.1)

Proof. Assume first that $\mathfrak{x}$ is a solution of the initial value problem (1.1). From definition of Caputo derivative, we have

$\mathfrak{I}^{2-\alpha}\mathfrak{D}^2\left(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\right) = \mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))).$

Operating by $\; \mathfrak{I}^{\alpha-1}\;$ on both sides and using Lemma 2, we get

$\mathfrak{I}^{1}\mathfrak{D}^2\left(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\right) = \mathfrak{I}^{\alpha-1}\mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))).$

Then

$\frac{d}{d\mathfrak{t}}\big(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big)- \frac{d}{d\mathfrak{t}}\big(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big)\bigg|_{\mathfrak{t} = 0} = \mathfrak{I}^{\alpha-1}\mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))).$

Using initial conditions, we have

$\frac{d}{d\mathfrak{t}}\big(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big) = \mathfrak{I}^{\alpha-1}\mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))).$

Integrating both sides of (1.1), we obtain

$\big(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big)-\big(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big)\bigg|_{\mathfrak{t} = 0} = \mathfrak{I}^{\alpha}\mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))).$

Then

$\mathfrak{x}(\mathfrak{t}) = \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}$

Conversely, assume that $\mathfrak{x}$ satisfies the nonlinear integral Eq (3.1). Then operating by $\; ^C\mathfrak{D}^\alpha \;$ on both sides of Eq (3.1) and using Lemma 2, we obtain

$\begin{eqnarray*} ^C\mathfrak{D}^\alpha \big(\mathfrak{x}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big) & = &^C\mathfrak{D}^\alpha \mathfrak{I}^\alpha \mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t})))\\ & = &\mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))). \end{eqnarray*}$

Putting $\mathfrak{t} = 0$ in (3.1) and since $\; \mathfrak{g}_1\;$ is a continuous function, then we obtain

$\big(\mathfrak{x}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big)\bigg|_{\mathfrak{t} = 0} = I^\alpha \mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t})))\bigg|_{\mathfrak{t} = 0} = 0.$

Also,

Then we have

$\frac{d}{d\mathfrak{t}}\big(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big)\bigg|_{\mathfrak{t} = 0} = \mathfrak{I}^{\alpha-1}\mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}(\mathfrak{t})))\bigg|_{\mathfrak{t} = 0} = 0.$

Hence the equivalence between the initial value problem (1.1) and the integral Eq (3.1) is proved. Then the proof is completed.

Definition 3. The Eq (1.1) is Ulam-Hyers stable if there exists a real number $\mathfrak{c}_{\mathfrak{f}} > 0$ such that for each $\epsilon > 0$ and for each solution $\mathfrak{z}\in C^{2}(J, R)$ of the inequality

$\mid ^C\mathfrak{D}^\alpha \left[ \mathfrak{z}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{z}(\mathfrak{t}))\right] - \mathfrak{g}_1( \mathfrak{t}, \mathfrak{z}(\mathfrak{t}), \mathfrak{I}^\beta\mathfrak{ g}_2(\mathfrak{t}, \mathfrak{z}(\mathfrak{t}))) \mid \leqslant \epsilon, \; \mathfrak{t}\in J,$

there exists a solution $\mathfrak{y}\in C^{2}(J, R)$ of Eq (1.1) with

$\mid \mathfrak{z}(\mathfrak{t})-\mathfrak{y}(\mathfrak{t})\vert \leqslant \mathfrak{c}_{\mathfrak{f}}\epsilon, \; \mathfrak{t}\in J .$

Definition 4. The Eq (1.1) is generalized Ulam-Hyers stable if there exists $\psi_{\mathfrak{f}}\in C(R_{+}, R_{+}), \; \psi_{\mathfrak{f}}(0) = 0,$ such that for each solution $\mathfrak{z}\in C^{2}(J, R)$ of the inequality

$\mid ^C\mathfrak{D}^\alpha \left[ \mathfrak{z}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{z}(\mathfrak{t}))\right] - \mathfrak{g}_1( \mathfrak{t}, \mathfrak{z}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{z}(\mathfrak{t}))) \mid \leqslant \epsilon , \; \mathfrak{t}\in J,$

there exists a solution $\mathfrak{y}\in C^{2}(J, R)\;$ of Eq (1.1) with

$\mid \mathfrak{z}(\mathfrak{t})-\mathfrak{y}(\mathfrak{t})\vert \leqslant \psi_{\mathfrak{f}} (\epsilon), \; \mathfrak{t}\in J .$

Definition 5. The Eq (1.1) is Ulam-Hyers-Rassias stable with respect to $\varphi\in C(J, R_{+})$ if there exists a real number $\mathfrak{c}_{\mathfrak{f}} > 0$ such that for each $\epsilon > 0$ and for each solution $\mathfrak{z}\in C^{2}(J, R)$ of the inequality

$\mid ^C\mathfrak{D}^\alpha \left[ \mathfrak{z}(\mathfrak{t})-\mathfrak{ h}(\mathfrak{t}, \mathfrak{z}(\mathfrak{t}))\right] - \mathfrak{g}_1( \mathfrak{t}, \mathfrak{z}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{z}(\mathfrak{t}))) \mid \leqslant \epsilon \varphi(\mathfrak{t}) , \; \mathfrak{t}\in J,$

there exists a solution $\mathfrak{y}\in C^{2}(J, R)$ of Eq (1.1) with

$\mid \mathfrak{z}(\mathfrak{t})-\mathfrak{y}(\mathfrak{t})\vert \leqslant \mathfrak{c}_\mathfrak{f}\epsilon\varphi(\mathfrak{t}), \; \mathfrak{t}\in J.$

Definition 6. The Eq (1.1) is generalized Ulam-Hyers-Rassias stable with respect to $\varphi\in C(J, R_{+})$ if there exists a real number $\mathfrak{c}_{\mathfrak{f}, \varphi} > 0$ such that for each solution $\mathfrak{z}\in C^{2}(J, R)$ of the inequality

$\mid ^C\mathfrak{D}^\alpha \left[ \mathfrak{z}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{z}(\mathfrak{t}))\right] - \mathfrak{g}_1( \mathfrak{t}, \mathfrak{z}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{z}(\mathfrak{t}))) \mid \leqslant \varphi(\mathfrak{t}) , \; \mathfrak{t}\in J,$

there exists a solution $\mathfrak{y}\in C^{2}(J, R)$ of Eq (1.1) with

$\mid \mathfrak{z}(\mathfrak{t})-\mathfrak{y}(\mathfrak{t})\vert \leqslant \mathfrak{c}_{\mathfrak{f}, \varphi} \varphi (\mathfrak{t}) , \; \mathfrak{t}\in J.$

Now, our aim is to investigate the existence of unique solution for (1.1). This existence result will be based on the contraction mapping principle.

Theorem 1. Let assumptions (i)-(iii) be satisfied. If $\; K_1+ \dfrac{K\; T^\alpha}{\Gamma(\alpha+1)}+\dfrac{K_2\; H\; T^{\alpha+\beta}}{\Gamma(\beta+1)\Gamma(\alpha+1)} < 1, \;$ then there exists a unique solution for the nonlinear neutral differential equation of fractional order.

Proof. Define the operator $\; N\;$ by:

$N\mathfrak{x}(\mathfrak{t}) = \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}, \; \mathfrak{t}\in J.$

In view of assumptions (i)-(iii), then $\; N: C^2(J, R) \rightarrow C^2(J, R)\;$ is continuous operator.

Now, let $\mathfrak{x}$ and $, \widetilde{\mathfrak{x}} \in C^2(J, R), \;$ be two solutions of (1.1)then

$\begin{eqnarray*} \mid N\mathfrak{x}(\mathfrak{t})-N\widetilde{\mathfrak{x}}(\mathfrak{t})\mid& = & \bigg| \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\\ &+&\dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}(\theta))d\theta\bigg)d\mathfrak{s}\\ \; \; \; \; \; \; &-&\mathfrak{h}(\mathfrak{t}, \widetilde{\mathfrak{x}}(\mathfrak{t}))\\ &-&\dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_{1}\bigg(\mathfrak{s}, \widetilde{\mathfrak{x}}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \widetilde{\mathfrak{x}}(\theta))d\theta\bigg)d\mathfrak{s}\bigg|\\ &\leqslant& K_1|\mathfrak{x}(\mathfrak{t})-\widetilde{\mathfrak{x}}(\mathfrak{t})|\\ &+& \dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg| \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}(\theta))d\theta\bigg)\\ \; \; \; \; \; \; &-& g_{1}\bigg(\mathfrak{s}, \widetilde{\mathfrak{x}}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \widetilde{\mathfrak{x}}(\theta))d\theta\bigg)\bigg| d\mathfrak{s}\\ \; \; \; \; \; \; &\leqslant&K_1|\mathfrak{x}(\mathfrak{t})-\widetilde{\mathfrak{x}}(\mathfrak{t})|\\ &+& \dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1} K\mid \mathfrak{x}(\mathfrak{s})-\widetilde{\mathfrak{x}}(\mathfrak{s})\vert d\mathfrak{s} \\ &+&H\dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1} \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mid \mathfrak{g}_{2}(\theta, \mathfrak{x}(\theta))- \mathfrak{g}_{2}(\theta, \widetilde{\mathfrak{x}}(\theta))\vert d\theta d\mathfrak{s}\\ \; \; \; \; \; \; &\leqslant&K_1|\mathfrak{x}(\mathfrak{t})-\widetilde{\mathfrak{x}}(\mathfrak{t})|+ \dfrac{K}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1} \mid \mathfrak{x}(\mathfrak{s})-\widetilde{\mathfrak{x}}(\mathfrak{s})\vert d\mathfrak{s} \\ &+& \dfrac{H}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\dfrac{K_2}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mid \mathfrak{x}(\theta)- \widetilde{\mathfrak{x}}(\theta)\vert d\theta d\mathfrak{s}. \end{eqnarray*}$

Then

$\begin{eqnarray*} || N\mathfrak{x}(\mathfrak{t})-N\widetilde{\mathfrak{x}}(\mathfrak{t})||&\leqslant&K_1||\mathfrak{x}-\widetilde{\mathfrak{x}}||+ \dfrac{K|| \mathfrak{x}-\widetilde{\mathfrak{x}}|\vert}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1} d\mathfrak{s} \\ &+& || \mathfrak{x}- \widetilde{\mathfrak{x}}\vert| \dfrac{H}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\dfrac{K_2}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1} d\theta d\mathfrak{s}\\ &\leqslant&K_1||\mathfrak{x}-\widetilde{\mathfrak{x}}||+ \dfrac{K|| \mathfrak{x}-\widetilde{\mathfrak{x}}|\vert\; T^\alpha}{\Gamma(\alpha+1)} \\ &+& || \mathfrak{x}- \widetilde{\mathfrak{x}}\vert|\dfrac{K_2\; T^\beta}{\Gamma(\beta+1)} \dfrac{H}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1} d\mathfrak{s}\\ &\leqslant&K_1||\mathfrak{x}-\widetilde{\mathfrak{x}}||+ \dfrac{K|| \mathfrak{x}-\widetilde{\mathfrak{x}}|\vert\; T^\alpha}{\Gamma(\alpha+1)} \\ &+& || \mathfrak{x}- \widetilde{\mathfrak{x}}\vert|\dfrac{K_2\; T^\beta}{\Gamma(\beta+1)} \dfrac{H\; T^\alpha}{\Gamma(\alpha+1)}\\ &\leq& \left[K_1+ \dfrac{K\; T^\alpha}{\Gamma(\alpha+1)}+\dfrac{K_2\; H\; T^{\alpha+\beta}}{\Gamma(\beta+1)\Gamma(\alpha+1)} \right]|| \mathfrak{x}-\widetilde{\mathfrak{x}}|\vert \end{eqnarray*}$

Since $\; K_1+ \dfrac{K\; T^\alpha}{\Gamma(\alpha+1)}+\dfrac{K_2\; H\; T^{\alpha+\beta}}{\Gamma(\beta+1)\Gamma(\alpha+1)} < 1. \;$ It follows that $N$ has a unique fixed point which is a solution of the initial value problem (1.1) in $\; C^2(J, R).\;$

4. Stability of solutions of the IVP (1.1)

4.1. Ulam-Hyers stability

Theorem 2. Let assumptions of Theorem 1 be satisfied. Then the fractional order differential Eq (1.1) is Ulam-Hyers stable.

Proof. Let $\; \mathfrak{y}\in C^2(J, R)\;$ be a solution of the inequality

$\begin{equation} \mid ^C\mathfrak{D}^\alpha \left[ \mathfrak{y}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))\right] - \mathfrak{g}_1( \mathfrak{t}, \mathfrak{y}(\mathfrak{t}), \mathfrak{I}^\beta\mathfrak{ g}_2(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))) \mid \leqslant \epsilon, \; \epsilon \gt 0, \; \mathfrak{t}\in J. \end{equation}$

(4.1)

Let $\; \mathfrak{x}\in C^2(J, R)\;$ be the unique solution of the initial value problem for implicit fractional-order differential Eq (1.1). By using Lemma 3, The Cauchy problem (1.1) is equivalent to

Operating by $\; \mathfrak{I}^{\alpha-1}\;$ on both sides of (4.1) and then integrating, we get

$\begin{eqnarray*} \bigg| \mathfrak{y}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))&-&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))d\theta\right)d\mathfrak{s}\bigg| \\ &\leqslant& \frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\epsilon \; d\mathfrak{s}, \\ &\leq & \frac{\epsilon\; T^\alpha}{\Gamma(\alpha+1)}. \end{eqnarray*}$

Also, we have

$\begin{eqnarray*} |\mathfrak{\mathfrak{y}}(\mathfrak{\mathfrak{t}})-\mathfrak{\mathfrak{x}}(\mathfrak{t})|& = & \bigg|\mathfrak{y}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))-\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{ x}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ & = &\bigg|\mathfrak{y}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))-\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\\ &+& \mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))d\theta\right)d\mathfrak{s}\\ &-& \mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))-\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ &\leq& \bigg|\mathfrak{y}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))-\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{ y}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ &+& | \mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))-\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))| \\ &+&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1} \bigg|\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))d\theta\right)\\ &-&\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)\bigg|d\mathfrak{s}\\ &\leq& \frac{\epsilon\; T^\alpha}{\Gamma(\alpha+1)}+K_1|\mathfrak{y}(\mathfrak{t})-\mathfrak{x}(\mathfrak{t})|\\ &+&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1} \left[K|\mathfrak{y}(\mathfrak{t})-\mathfrak{x}(\mathfrak{t})|+\frac{H}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}|\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))-\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))|d\theta\right]d\mathfrak{s} \end{eqnarray*}$

$\begin{eqnarray*} ||\mathfrak{y}-\mathfrak{x}||&\leq& \frac{\epsilon \; T^\alpha}{\Gamma(\alpha+1)}+K_1||\mathfrak{y}-\mathfrak{x}||\\ &+&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1} \left[K||\mathfrak{y}-\mathfrak{x}||+\frac{H\; K_2||\mathfrak{x}-\mathfrak{y}||T^\beta}{\Gamma(\beta+1)}\right]d\mathfrak{s}\\ &\leq& \frac{\epsilon\; T^\alpha}{\Gamma(\alpha+1)}+K_1||\mathfrak{y}-\mathfrak{x}||\\ &+&\frac{K\; T^\alpha||\mathfrak{y}-\mathfrak{x}||}{\Gamma(\alpha+1)} +\frac{H\; K_2T^{\beta+\alpha} ||\mathfrak{x}-\mathfrak{y}||}{\Gamma(\beta+1)\Gamma(\alpha+1)}. \end{eqnarray*}$

Then

$||\mathfrak{y}-\mathfrak{x}||\leq \frac{\epsilon \; T^\alpha}{\Gamma(\alpha+1)}\left[1- \left( K_1+\frac{K\; T^\alpha}{\Gamma(\alpha+1)} +\frac{H\; K_2\; T^{\beta+\alpha} }{\Gamma(\beta+1)\Gamma(\alpha+1)}\right)\right]^{-1} = \mathfrak{c}\; \epsilon,$

thus the intial value problem (1.1) is Ulam-Heyers stable, and hence the proof is completed. By putting $\psi(\varepsilon) = \mathfrak{c}\; \varepsilon, \; \psi(0) = 0 \;$ yields that the Eq (1.1) is generalized Ulam-Heyers stable.

4.2. Ulam-Hyers-Rassias stability

Theorem 3. Let assumptions of Theorem 1 be satisfied, there exists an increasing function $\varphi\in C(J, R)$ and there exists $\lambda_{\varphi} > 0$ such that for any $\mathfrak{t}\in J,$ we have

$\mathfrak{I}^{\alpha} \varphi (\mathfrak{t})\leqslant \lambda_{\varphi} \varphi(\mathfrak{t}),$

then the Eq (1.1) is Ulam-Heyers-Rassias stable.

Proof. Let $\; \mathfrak{y}\in C^2(J, R)\;$ be a solution of the inequality

$\begin{equation} \mid ^C\mathfrak{D}^\alpha \left[ \mathfrak{y}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))\right] - \mathfrak{g}_1( \mathfrak{t}, \mathfrak{y}(\mathfrak{t}), \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))) \mid \leqslant \epsilon \varphi(\mathfrak{t}), \; \epsilon \gt 0, \; \mathfrak{t}\in J. \end{equation}$

(4.2)

Operating by $\; \mathfrak{I}^{\alpha-1}\;$ on both sides of (4.2) and then integrating, we get

$\begin{eqnarray*} \bigg| \mathfrak{y}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))&-&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))d\theta\right)d\mathfrak{s}\bigg| \\ &\leqslant& \frac{\epsilon}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\varphi(\mathfrak{s}) d\mathfrak{s}, \\ &\leq& \epsilon\; \mathfrak{I}^{\alpha} \varphi (\mathfrak{t})\\ &\leq & \epsilon\; \lambda_{\varphi} \varphi(\mathfrak{t}). \end{eqnarray*}$

Also, we have

$\begin{eqnarray*} |\mathfrak{y}(\mathfrak{t})-\mathfrak{x}(\mathfrak{t})|& = & \bigg|\mathfrak{y}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))-\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ & = &\bigg|\mathfrak{y}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))-\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\\ &+& \mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{ y}(\theta))d\theta\right)d\mathfrak{s}\\ &-& \mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))-\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ &\leq& \bigg|\mathfrak{y}(\mathfrak{t})-\mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))-\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ &+& | \mathfrak{h}(\mathfrak{t}, \mathfrak{y}(\mathfrak{t}))-\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))| \\ &+&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg |\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{y}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))d\theta\right)\\ &-&\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)\bigg|d\mathfrak{s}\\ &\leq& \epsilon \; \lambda_{\varphi} \varphi(\mathfrak{t}) +K_1|\mathfrak{y}(\mathfrak{t})-\mathfrak{x}(\mathfrak{t})|\\ &+&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1} \left[K|\mathfrak{y}(\mathfrak{t})-\mathfrak{x}(\mathfrak{t})|+\frac{H}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}- \theta)^{\beta-1}|\mathfrak{g}_2(\theta, \mathfrak{y}(\theta))-\mathfrak{g}_2(\theta, \mathfrak{ x}(\theta))|d\theta\right]d\mathfrak{s} \end{eqnarray*}$

$\begin{eqnarray*} ||\mathfrak{y}-\mathfrak{x}||&\leq& \epsilon\; \lambda_{\varphi} \varphi(\mathfrak{t})+K_1||\mathfrak{y}-\mathfrak{x}||\\ &+&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1} \left[K||\mathfrak{y}(\mathfrak{t})-\mathfrak{x}(\mathfrak{t})||+\frac{H\; K_2||\mathfrak{x}-\mathfrak{y}||T^\beta}{\Gamma(\beta+1)}\right]d\mathfrak{s}\\ &\leq&\epsilon \lambda_{\varphi} \varphi(\mathfrak{t})+K_1||\mathfrak{y}-\mathfrak{x}||\\ &+&\frac{K\; T\; ^\alpha||\mathfrak{y}-\mathfrak{x}||}{\Gamma(\alpha+1)} +\frac{H\; K_2\; T^{\beta+\alpha} ||\mathfrak{x}-\mathfrak{y}||}{\Gamma(\beta+1)\Gamma(\alpha+1)}. \end{eqnarray*}$

Then

$||\mathfrak{y}-\mathfrak{x}||\leq \epsilon\; \lambda_{\varphi}\; \varphi(\mathfrak{t})\left[1- \left( K_1+\frac{K\; T^\alpha}{\Gamma(\alpha+1)} +\frac{H\; K_2\; T^{\beta+\alpha} }{\Gamma(\beta+1)\Gamma(\alpha+1)}\right)\right]^{-1} = \mathfrak{c}\; \epsilon \varphi(\mathfrak{t}),$

then the initial problem (1.1) is Ulam-Heyers-Rassias stable, and hence the proof is completed.

5. Existence and attractivity of solutions on half line

In this section, we prove some results on the existence of mild solutions and attractivity for the neutral fractional differential equation (1.1) by applying Schauder fixed point theorem. Denote $\mathfrak{B}C = \mathfrak{ B}C(J), \; J = [0, +\infty)\;$ and consider the following assumptions:

(I) $\mathfrak{h}:J \times R \longrightarrow R$ is a continuous function and there exists a continuous function $\; K_\mathfrak{h}(\mathfrak{t})\;$ such that:

$\mid \mathfrak{h}(\mathfrak{t}, \mathfrak{x})- \mathfrak{h}(\mathfrak{t}, \mathfrak{y})\mid \leqslant K_\mathfrak{h}(\mathfrak{t}) \mid \mathfrak{x}-\mathfrak{y}\mid \; \; \mbox{ for each }\; \; \mathfrak{t}\in J \; \mbox{ and } \; \mathfrak{x}, \mathfrak{y}\in R,$

where $\; K_\mathfrak{h}^* = \sup\limits_{\mathfrak{t}\geq 0}K_\mathfrak{h}(\mathfrak{t}) < 1, \; \lim\limits_{\mathfrak{\mathfrak{t}}\rightarrow \infty} K_\mathfrak{h}(\mathfrak{t}) = 0, \;$ and $\; \lim\limits_{\mathfrak{t}\rightarrow \infty}\mathfrak{h}(\mathfrak{t}, 0) = 0.$

(II) $\mathfrak{g} _{1}: J \times R \times R\longrightarrow R$ satisfies Carathéodory condition and there exist an integrable function $\mathfrak{a}_1:R_+\longrightarrow R_+$ and a positive constant $\mathfrak{b}$ such that:

$\mid \mathfrak{g} _{1}(\mathfrak{t}, \mathfrak{x}, \mathfrak{y})\mid \leq \frac{\mathfrak{a}_1(\mathfrak{t})}{1+|\mathfrak{x}|}+\mathfrak{b}|\mathfrak{y}|\; \; \; \mbox{ for each}\; \; \mathfrak{t}\in J \; \mbox{ and }\; \mathfrak{x}, \; \mathfrak{y} \in R .$

(III) $\mathfrak{g} _{2}:J \times R \longrightarrow R$ satisfies Carathéodory condition and there exists an integrable function

$\mathfrak{a}_2:R_+\longrightarrow R_+$ such that:

$\mid \mathfrak{g} _{2}(\mathfrak{t}, \mathfrak{x})\mid \leq \frac{\mathfrak{a}_2(\mathfrak{t})}{1+|\mathfrak{x}|} \; \; \; \mbox{ for each}\; \; \mathfrak{t}\in J \; \mbox{ and } \; \mathfrak{x} \in R.$

(IV) Let

$\begin{eqnarray*} \lim\limits_{\mathfrak{t} \rightarrow \infty} \int_0^\mathfrak{t} \frac{(\mathfrak{t}-\mathfrak{s})^{\alpha-1}}{\Gamma(\alpha)}\mathfrak{a}_1(\mathfrak{s})\; d\mathfrak{s} = 0\\ \\ \mathfrak{a}_1^* = \sup\limits_{\mathfrak{t}\in J} \int_0^\mathfrak{t} \frac{(\mathfrak{t}-\mathfrak{s})^{\alpha-1}}{\Gamma(\alpha)}\mathfrak{a}_1(\mathfrak{s})\; d\mathfrak{s}\\ \\ \lim\limits_{\mathfrak{t} \rightarrow \infty} \int_0^\mathfrak{t} \frac{(\mathfrak{t}-\mathfrak{s})^{\alpha+\beta-1}}{\Gamma(\alpha+\beta)}\mathfrak{a}_2(\mathfrak{s})\; d\mathfrak{s} = 0\\ \\ \mathfrak{a}_2^* = \sup\limits_{\mathfrak{t}\in J} \int_0^\mathfrak{t} \frac{(\mathfrak{t}-\mathfrak{s})^{\alpha+\beta-1}}{\Gamma(\alpha+\beta)}\mathfrak{a}_2(\mathfrak{s})\; d\mathfrak{s} \end{eqnarray*}$

By a mild solution of the Eq (1.1) we mean that a function $x \in C(J, R)$ such that $\mathfrak{x}$ satisfies the equation in (3.1).

Theorem 4. Let assumptions (I)-(IV) be satisfied. Then there exists at least one mild solution for the nonlinear implicit neutral differential equation of fractional order (1.1). Moreover, mild solutions of IVP (1.1) are locally attractive.

Proof. For any $\mathfrak{x}\in \mathfrak{B}C, \;$ define the operator $\; \mathbb{A} \;$ by

$\mathbb{A}\mathfrak{x}(\mathfrak{t}) = \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{ x}(\theta))d\theta\right)d\mathfrak{s}.$

The operator $\; \mathbb{A} \;$ is well defined and maps $\; \mathfrak{B}C\;$ into $\; \mathfrak{B}C.\;$ Obviously, the map $\; \mathbb{A}(\mathfrak{x}) \;$ is continuous on $\; J\;$ for any $\; \mathfrak{x} \in \mathfrak{B}C\;$ and for each $\; \mathfrak{t} \in J, \;$ we have

$\begin{eqnarray*} |\mathbb{A}\mathfrak{x}(\mathfrak{t})|&\leq&|\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))|+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg|\mathfrak{g}_1\bigg(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\bigg)\bigg|d\mathfrak{s}\\ &\leq& |\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))-\mathfrak{h}(\mathfrak{t}, 0)|+|\mathfrak{h}(\mathfrak{t}, 0)|\\ &+&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg[\frac{\mathfrak{a}_1(\mathfrak{s})} {1+|\mathfrak{x}(\mathfrak{s})|}+\mathfrak{b} \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}|\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))|d\theta\bigg]d\mathfrak{s}\\ &\leq& K_\mathfrak{h}(\mathfrak{t})|\mathfrak{x}(\mathfrak{t})|+|\mathfrak{h}(\mathfrak{t}, 0)|+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg[\mathfrak{a}_1(\mathfrak{s})+\mathfrak{b} \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\frac{\mathfrak{a}_2(\theta)}{1+|\mathfrak{x}(\theta)|}d\theta\bigg]d\mathfrak{s}\\ &\leq& K_\mathfrak{h}^*|\mathfrak{x}(\mathfrak{t})|+|\mathfrak{h}(\mathfrak{t}, 0)|+\mathfrak{a}_1^*+\mathfrak{b}\frac{1}{\Gamma(\alpha+\beta)} \int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha+\beta-1}\frac{\mathfrak{a}_2(\theta)}{1+|\mathfrak{x}(\theta)|}d\theta\\ &\leq& K_\mathfrak{h}^*M+|\mathfrak{h}(\mathfrak{t}, 0)|+\mathfrak{a}_1^*+\mathfrak{b}\frac{1}{\Gamma(\alpha+\beta)} \int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha+\beta-1}\mathfrak{a}_2(\theta)d\theta\\ &\leq& K_\mathfrak{h}^*M+|\mathfrak{h}(\mathfrak{t}, 0)|+\mathfrak{a}_1^*+\mathfrak{b} \mathfrak{a}_2^*\leq M. \end{eqnarray*}$

Then

$\begin{equation} ||\mathbb{A}\mathfrak{x}(\mathfrak{t})||_{\mathfrak{B}C} \leq M, \; \; M = \big(|\mathfrak{h}(\mathfrak{t}, 0)|+\mathfrak{a}_1^*+\mathfrak{b} \mathfrak{a}_2^* \big)\big(1-K_\mathfrak{h}^*\big)^{-1}. \end{equation}$

(5.1)

Thus $\; \mathbb{A}(\mathfrak{x}) \in \mathfrak{B}C.\;$ This clarifies that operator $\; \mathbb{A} \;$ maps $\; \mathfrak{B}C\;$ into itself.

Finding the solutions of IVP (1.1) is reduced to find solutions of the operator equation $\; \mathbb{A}(\mathfrak{x}) = \mathfrak{x}.\;$ Eq (5.1) implies that $\; \mathbb{A}\;$ maps the ball $\; \mathbb{B}_M : = \mathbb{B}(0, M) = \{\mathfrak{x}\in \mathfrak{B}C : ||\mathfrak{x}(\mathfrak{t})||_{\mathfrak{B}C} \leq M\}\;$ into itself. Now, our proof will be established in the following steps:

Step 1: $\; \mathbb{A} \;$ is continuous.

Let $\; \{ \mathfrak{x}_n\}_{n\in N}\;$ be a sequence such that $\; \mathfrak{x}_n \rightarrow \mathfrak{x}\;$ in $\; \mathbb{B}_M.\;$ Then, for each $\; \mathfrak{t} \in J, \;$ we have

$\begin{eqnarray*} \mid \mathbb{A}\mathfrak{x}_n(\mathfrak{t})-\mathbb{A}{\mathfrak{x}}(\mathfrak{t})\mid& = & \bigg| \mathfrak{h}(\mathfrak{t}, \mathfrak{x}_n(\mathfrak{t}))\\ &+&\dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}_n(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}_n(\theta))d\theta\bigg)d\mathfrak{s}\\ \; \; \; \; \; \; &-&\mathfrak{h}(\mathfrak{t}, {\mathfrak{x}}(\mathfrak{t}))\\ &-&\dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_{1}\bigg(\mathfrak{s}, {\mathfrak{x}}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, {\mathfrak{x}}(\theta))d\theta\bigg)d\mathfrak{s}\bigg|\\ &\leqslant& K_\mathfrak{h}(\mathfrak{t})|\mathfrak{x}_n(\mathfrak{t})-{\mathfrak{x}}(\mathfrak{t})|\\ &+& \dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg| \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}_n(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}_n(\theta))d\theta\bigg)\\ \; \; \; \; \; \; &-& \mathfrak{g}_{1}\bigg(\mathfrak{s}, {\mathfrak{x}}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, {\mathfrak{x}}(\theta))d\theta\bigg)\bigg| d\mathfrak{s}\\ &\leqslant& K_\mathfrak{h}^*|\mathfrak{x}_n(\mathfrak{t})-{\mathfrak{x}}(\mathfrak{t})|\\ &+& \dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg| \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}_n(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}_n(\theta))d\theta\bigg)\\ \; \; \; \; \; \; &-& \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}(\theta))d\theta\bigg)\bigg| d\mathfrak{s} \end{eqnarray*}$

Assumptions (II) and (III) implies that:

$\mathfrak{g}_1\bigg(\mathfrak{t}, \mathfrak{x}_n, \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x}_n)\bigg)\rightarrow \mathfrak{g}_1\bigg(\mathfrak{t}, \mathfrak{x}, \mathfrak{I}^\beta \mathfrak{g}_2(\mathfrak{t}, \mathfrak{x})\bigg) \mbox{\; as \; } n\rightarrow \infty.$

Using Lebesgue dominated convergence theorem, we have

$|| \mathbb{A}\mathfrak{x}_n(\mathfrak{t})-\mathbb{A}\mathfrak{x}(\mathfrak{t})||_{\mathfrak{B}C} \rightarrow 0 \; \; \mbox{ as}\; \; n\rightarrow \infty.$

Step 2: $\; \mathbb{A}(\mathbb{B}_M)\;$ is uniformly bounded.

It is obvious since $\; \mathbb{A}(\mathbb{B}_M) \subset \mathbb{B}_M\;$ and $\; \mathbb{B}_M\;$ is bounded.

Step 3: $\; \mathbb{A}(\mathbb{B}_M)\;$ is equicontinuous on every compact subset $\; [0, T]\;$ of $\; J, \; T > 0\;$ and $\mathfrak{t} _{1}, \mathfrak{t} _{2} \in [0, T], \; \mathfrak{t}_{2} > \mathfrak{t}_{1}$ (without loss of generality), we get

$\mid \mathbb{A} \mathfrak{x} (\mathfrak{t}_{2} )- \mathbb{A}\mathfrak{x}(\mathfrak{t}_{1}) \mid \leq \bigg| \mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_2))+\frac{1}{\Gamma(\alpha)}\int_0^{\mathfrak{t}_2} (\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\\ - \mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))+\frac{1}{\Gamma(\alpha)}\int_0^{\mathfrak{t}_1}(\mathfrak{t}_{1}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ \leq \mid \mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_2))- \mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))|$

$+ \frac{1}{\Gamma(\alpha)}\bigg| \int_0^{\mathfrak{t}_2}(\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s} \\ - \int_0^{\mathfrak{t}_1}(\mathfrak{t}_{1}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ \leq \mid \mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_2))- \mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))+\mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_1))- \mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_1))| \\ + \frac{1}{\Gamma(\alpha)}\bigg| \int_0^{\mathfrak{t}_1}(\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\\ + \frac{1}{\Gamma(\alpha)} \int_{\mathfrak{t}_1}^{\mathfrak{t}_2}(\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s} \\ - \int_0^{\mathfrak{t}_1}(\mathfrak{t}_{1}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ \leq K_\mathfrak{h}(\mathfrak{t})\mid \mathfrak{x}(\mathfrak{t}_2)-\mathfrak{x}(\mathfrak{t}_1)\mid+|\mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_1))-\mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))| \\ + \frac{1}{\Gamma(\alpha)}\bigg| \int_0^{\mathfrak{t}_1}(\mathfrak{t}_{1}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\\ + \frac{1}{\Gamma(\alpha)} \int_{\mathfrak{t}_1}^{\mathfrak{t}_2}(\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\mathfrak{\mathfrak{g}}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s} \\ - \int_0^{\mathfrak{t}_1}(\mathfrak{t}_{1}-\mathfrak{s})^{\alpha-1}\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)d\mathfrak{s}\bigg|\\ \leq K_\mathfrak{h}(\mathfrak{t})\mid \mathfrak{x}(\mathfrak{t}_2)-\mathfrak{x}(\mathfrak{t}_1)\mid+|\mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_1))-\mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))| \\ + \frac{1}{\Gamma(\alpha)} \int_{\mathfrak{t}_1}^{\mathfrak{t}_2}(\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\bigg|\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)\bigg|d\mathfrak{s}\\ \leq K_\mathfrak{h}(\mathfrak{t})\mid \mathfrak{x}(\mathfrak{t}_2)-\mathfrak{x}(\mathfrak{t}_1)\mid+|\mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_1))-\mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))| \\ + \frac{1}{\Gamma(\alpha)} \int_{\mathfrak{t}_1}^{\mathfrak{t}_2}(\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\bigg[\frac{\mathfrak{a}_1(\mathfrak{s})}{1+|\mathfrak{x}(\mathfrak{s})|}+\mathfrak{b} \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}|\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))|d\theta \bigg]d\mathfrak{s} \\ \leq K_\mathfrak{h}^*\mid \mathfrak{x}(\mathfrak{t}_2)-\mathfrak{x}(\mathfrak{t}_1)\mid+|\mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_1))-\mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))| \\ + \frac{1}{\Gamma(\alpha)} \int_{\mathfrak{t}_1}^{\mathfrak{t}_2}(\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\bigg[\mathfrak{a}_1(\mathfrak{s})+\mathfrak{b} \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\frac{\mathfrak{a}_2(\theta)}{1+|\mathfrak{x}(\theta)|}d\theta \bigg]d\mathfrak{s} \\ \leq K_\mathfrak{h}^*\mid \mathfrak{x}(\mathfrak{t}_2)-\mathfrak{x}(\mathfrak{t}_1)\mid+|\mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_1))-\mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))| \\ + \frac{1}{\Gamma(\alpha)} \int_{\mathfrak{t}_1}^{\mathfrak{t}_2}(\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\bigg[\mathfrak{a}_1(\mathfrak{s})+\mathfrak{b} \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{a}_2(\theta)d\theta \bigg]d\mathfrak{s}.$

Thus, for $\; \mathbf{a}_\mathfrak{i} = \sup\limits_{ \mathfrak{t}\in [0, T]}\mathfrak{a}_\mathfrak{i}, \; \mathfrak{i} = 1, 2\;$ and from the continuity of the functions $\; \mathfrak{a}_\mathfrak{i}\;$ we obtain

$\begin{eqnarray*} \mid \mathbb{A} \mathfrak{x} (\mathfrak{t}_{2} )- \mathbb{A}\mathfrak{x}(\mathfrak{t}_{1}) \mid&\leq& K_\mathfrak{h}^*\mid \mathfrak{x}(\mathfrak{t}_2)-\mathfrak{x}(\mathfrak{t}_1)\mid+|\mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_1)) -\mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))| \\ &+&\frac{1}{\Gamma(\alpha)} \int_{\mathfrak{t}_1}^{\mathfrak{t}_2}(\mathfrak{t}_{2}-\mathfrak{s})^{\alpha-1}\bigg[\mathfrak{a}_1(\mathfrak{s})+ \frac{\mathfrak{b}\mathbf{a}_2}{\Gamma(\beta+1)}\mathfrak{s}^\beta \bigg]d\mathfrak{s}.\\ &\leq& K_\mathfrak{h}^*\mid \mathfrak{x}(\mathfrak{t}_2)-\mathfrak{x}(\mathfrak{t}_1)\mid+|\mathfrak{h}(\mathfrak{t}_{2}, \mathfrak{x}(\mathfrak{t}_1))-\mathfrak{h}(\mathfrak{t}_{1}, \mathfrak{x}(\mathfrak{t}_1))|\\ &+&\frac{\mathbf{a}_1}{\Gamma(\alpha+1)} (\mathfrak{t}_{2}-\mathfrak{t}_1)^{\alpha}+ \frac{\mathfrak{b}\mathbf{a}_2}{\Gamma(\alpha+\beta+1)}(\mathfrak{t}_{2}-\mathfrak{t}_1)^{\alpha+\beta} . \end{eqnarray*}$

Continuity of $\; \mathfrak{h}\;$ implies that

$|\; (\mathbb{A}\mathfrak{x})(\mathfrak{t}_2) \; -\; (\mathbb{A}\mathfrak{x})(\mathfrak{t}_1)\; |\; \rightarrow 0\; \; \; \; \mbox{as}\; \; \; \mathfrak{t}_2\rightarrow \mathfrak{t}_1.$

Step 4: $\; \mathbb{A}(\mathbb{B}_M)\;$ is equiconvergent.

Let $\; \mathfrak{t} \in J\;$ and $\; \mathfrak{x} \in \mathbb{B}_M\;$ then we have

$\begin{eqnarray*} |\mathbb{A}\mathfrak{x}(\mathfrak{t})|&\leq&|\mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))-\mathfrak{h}(\mathfrak{t}, 0)|+|\mathfrak{h}(\mathfrak{t}, 0)|\\ &+&\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg|\mathfrak{g}_1\left(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))d\theta\right)\bigg|d\mathfrak{s}\\ &\leq& K_\mathfrak{h}(\mathfrak{t})|\mathfrak{x}(\mathfrak{t})|+ |\mathfrak{h}(\mathfrak{t}, 0)|+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg[\frac{\mathfrak{a}_1(\mathfrak{s})}{1+|\mathfrak{x}(\mathfrak{s})|}+\mathfrak{b} \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}|\mathfrak{g}_2(\theta, \mathfrak{x}(\theta))|d\theta\bigg]d\mathfrak{s}\\ &\leq&K_\mathfrak{h}(\mathfrak{t})|\mathfrak{x}(\mathfrak{t})|+ |\mathfrak{h}(\mathfrak{t}, 0)|+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg[\mathfrak{a}_1(\mathfrak{s})+\mathfrak{b} \frac{1}{\Gamma(\beta)}\int_0^\mathfrak{s}(\mathfrak{s}-\theta)^{\beta-1}\frac{\mathfrak{a}_2(\theta)}{1+|\mathfrak{x}(\theta)|}d\theta\bigg]d\mathfrak{s}\\ &\leq& K_\mathfrak{h}(\mathfrak{t})|\mathfrak{x}(\mathfrak{t})|+ |\mathfrak{h}(\mathfrak{t}, 0)|+\frac{1}{\Gamma(\alpha)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha-1}\mathfrak{a}_1(\mathfrak{s})d\mathfrak{s}+\mathfrak{b} \frac{1}{\Gamma(\alpha+\beta)}\int_0^\mathfrak{t}(\mathfrak{t}-\mathfrak{s})^{\alpha+\beta-1} \mathfrak{a}_2(\mathfrak{s})d\mathfrak{s}. \end{eqnarray*}$

In view of assumptions (I) and (IV), we obtain

$|\mathbb{A}\mathfrak{x}(\mathfrak{t})|\rightarrow 0 \; \; \mbox{ as }\; \; \mathfrak{t} \rightarrow \infty.$

Then $\; \mathbb{A}\;$ has a fixed point $\; \mathfrak{x}\;$ which is a solution of IVP (1.1) on $\; J.\;$

Step 5: Local attactivity of mild solutions. Let $\; \mathfrak{x}^*\;$ be a mild solution of IVP (1.1). Taking $\; \mathfrak{x}\in \mathbb{B}(\mathfrak{x}^*, 2M), \;$ we have

$\begin{eqnarray*} \big|\mathbb{A}\mathfrak{x}(\mathfrak{t})-\mathfrak{x}^*\; (\mathfrak{t}) \big| & = &\big|\mathbb{A}\mathfrak{x}(\mathfrak{t})-\mathbb{A}\mathfrak{x}^*\; (\mathfrak{t}) \big|\\ &\leq& \mid \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}^*(t))\mid\\ &+& \dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg| \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}(\theta))d\theta\bigg)\\ &-& \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}^*(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}^*(\theta))d\theta\bigg) \bigg| d\mathfrak{s}\\ &\leq& K_\mathfrak{h}(\mathfrak{t})\mid \mathfrak{x}(\mathfrak{t})- \mathfrak{x}^*(\mathfrak{t})\mid\\ &+& \dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg[\bigg| \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}(\theta))d\theta\bigg)\bigg|\\ &+& \bigg| \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}^*(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}^*(\theta))d\theta\bigg) \bigg|\bigg] d\mathfrak{s}\\ &\leq& K_\mathfrak{h}^*\mid \mathfrak{x}(\mathfrak{t})- \mathfrak{x}^*(\mathfrak{t})\mid\\ &+& \dfrac{2}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg| \mathfrak{a}_1(\mathfrak{s})+\dfrac{\mathfrak{b}}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{a}_{2}(\theta)d\theta\bigg|d\mathfrak{s}\\ &\leq& K_\mathfrak{h}^*\mid \mathfrak{x}(\mathfrak{t})- \mathfrak{x}^*(\mathfrak{t})\mid + 2\mathfrak{a}_1^*+2\mathfrak{b}\dfrac{1}{\Gamma(\beta+\alpha)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\alpha+\beta-1}\mathfrak{a}_{2}(\theta)d\theta \\ &\leq&2\big( K_\mathfrak{h}^*\mid \mathfrak{x}(\mathfrak{t})\mid+ |\mathfrak{h}(\mathfrak{t}, 0)| + \mathfrak{a}_1^*+\mathfrak{b} \mathfrak{a}_2^*\big)\\ &\leq& 2\big( K_\mathfrak{h}^*M+| \mathfrak{h}(t, 0)| + \mathfrak{a}_1^*+\mathfrak{b} \mathfrak{a}_2^*\big)\leq 2M. \end{eqnarray*}$

We have

$||\mathbb{A}\mathfrak{x}(\mathfrak{t})-\mathfrak{x}^*(\mathfrak{t})||_{\mathfrak{B}C}\leq 2M.$

Hence $\; \mathbb{A}\;$ is a continuous function such that $\; \mathbb{A}(\mathbb{B}(\mathfrak{x}^*, 2M))\subset \mathbb{B}(\mathfrak{x}^*, 2M).\;$

Moreover, if $\; x\;$ is a mild solution of IVP (1.1), then

$\begin{eqnarray*} \big|\mathfrak{x}(\mathfrak{t})-\mathfrak{x}^*(\mathfrak{t}) \big| & = & \big|\mathbb{A}\mathfrak{x}(\mathfrak{t})-\mathbb{A}\mathfrak{x}^*(\mathfrak{t}) \big|\\ &\leq& \mid \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}^*(\mathfrak{t}))\mid\\ &+& \dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg| \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}(\theta))d\theta\bigg)\\ &-& \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}^*(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}^*(\theta))d\theta\bigg) \bigg| d\mathfrak{s}\\ &\leq& K_\mathfrak{h}(\mathfrak{t})\mid \mathfrak{x}(\mathfrak{t})- \mathfrak{x}^*(\mathfrak{t})\mid\\ &+& \dfrac{1}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1}\bigg[\bigg| \mathfrak{g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}(\theta))d\theta\bigg)\bigg|\\ &+& \bigg|\mathfrak{ g}_{1}\bigg(\mathfrak{s}, \mathfrak{x}^*(\mathfrak{s}), \dfrac{1}{\Gamma(\beta)}\int_{0}^{\mathfrak{s}} (\mathfrak{s}-\theta)^{\beta-1}\mathfrak{g}_{2}(\theta, \mathfrak{x}^*(\theta))d\theta\bigg) \bigg|\bigg] d\mathfrak{s}\\ &\leq& K_\mathfrak{h}^*\mid \mathfrak{x}(\mathfrak{t})-\mathfrak{ x}^*(\mathfrak{t})\mid\\ &+& \dfrac{2}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1} \mathfrak{a}_1(\mathfrak{s})d\mathfrak{s}+\dfrac{2\mathfrak{b}}{\Gamma(\alpha+\beta)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\theta)^{\alpha+\beta-1}\mathfrak{a}_{2}(\theta)d\theta. \end{eqnarray*}$

Then

$\begin{equation} |\mathfrak{x}(\mathfrak{t})-\mathfrak{x}^*(\mathfrak{t})|\leq (1-K_\mathfrak{h}^*)^{-1} \bigg[ \dfrac{2}{\Gamma(\alpha)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\mathfrak{s})^{\alpha-1} \mathfrak{a}_1(\mathfrak{s})d\mathfrak{s}+\dfrac{2b}{\Gamma(\alpha+\beta)}\int_{0}^{\mathfrak{t}} (\mathfrak{t}-\theta)^{\alpha+\beta-1}\mathfrak{a}_{2}(\theta)d\theta \bigg]. \end{equation}$

(5.2)

In view of assumption of (IV) and estimation (5.2), we get

$\lim\limits_{t\rightarrow \infty} |\mathfrak{x}(\mathfrak{t})-\mathfrak{x}^*(\mathfrak{t})| = 0.$

Then, all mild solutions of IVP (1.1) are locally attractive.

6. Applications

As particular cases of the IVP (1.1), we have

● Taking $\; \mathfrak{g}_1(\mathfrak{t}, \mathfrak{x}, \mathfrak{y}) = \mathfrak{g}_1(\mathfrak{t}, \mathfrak{x}), \;$ we obtain the initial value problem

$\begin{align*} \left\{ \begin{array}{c} ^C\mathfrak{D}^\alpha \left[ \mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\right] = \mathfrak{g}_1( \mathfrak{t}, \mathfrak{x}(\mathfrak{t})) \; \; \; \mathfrak{t}\in J, \; 1 \lt \alpha \leq 2 , \\\\ \bigg(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\bigg)\bigg|_{\mathfrak{t} = 0} = 0\; \; \; \; \mbox{and}\; \; \; \; \; \frac{d}{d\mathfrak{t}}\big[\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big]_{\mathfrak{t} = 0} = 0\end{array} \right. \end{align*}$

● Letting $\; \alpha \rightarrow 2, \; \beta \rightarrow 1, \;$ as a particular case of Theorem 1 we can deduce an existence result for the initial value problem for implicit second-order differe-integral equation

$\begin{equation*} \left\{ \begin{array}{c} \frac{d^2}{d\mathfrak{t}^2} \left( \mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\right) = \mathfrak{g}_1\bigg( \mathfrak{t}, \mathfrak{x}(\mathfrak{t}), \int_0^\mathfrak{t } \mathfrak{g}_2(\mathfrak{s}, \mathfrak{x}(\mathfrak{s}))d\mathfrak{s}\bigg) \; \; \; \mathfrak{t}\in J, \\\\ \big(\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\big)\bigg|_{\mathfrak{t} = 0} = 0\; \; \; \; \mbox{and}\; \; \; \; \; \frac{d}{d\mathfrak{t}}\bigg[\mathfrak{x}(\mathfrak{t})- \mathfrak{h}(\mathfrak{t}, \mathfrak{x}(\mathfrak{t}))\bigg]_{\mathfrak{t} = 0} = 0\end{array} \right. \end{equation*}$

As particular cases we can deduce existence results for some initial value problem of second order differential equations (when $\; \mathfrak{h} = 0\;$ ) and $\; \alpha \rightarrow 2, \;$ we get:

● Taking $\; \mathfrak{g}_1(\mathfrak{t}, \mathfrak{x}, \mathfrak{y}) = -\lambda^2 \mathfrak{x}(\mathfrak{t}), \; \lambda \in R^+,$ then we obtain a second order differential equation of simple harmonic oscillator

$\begin{equation*} \left\{ \begin{array}{c} \frac{d^2\mathfrak{x}(\mathfrak{t})}{d\mathfrak{t}^2} = -\lambda^2 \mathfrak{x}(\mathfrak{t}) \; \; \; \mathfrak{t}\in J, \\\\ \mathfrak{x}(0) = 0\; \; \; \; \mbox{and}\; \; \; \; \; \mathfrak{x}'(0) = 0\end{array} \right. \end{equation*}$

● Taking $\; \mathfrak{g}_1(\mathfrak{t}, \mathfrak{x}, \mathfrak{y}) = \left(\frac{\mathfrak{t}^2-k}{\mathfrak{t}^2}\right)\mathfrak{x}+\mathfrak{q}(\mathfrak{x}), \; k \in R\;$ where $\; \mathfrak{q}(\mathfrak{x})\;$ is continuous function, then we obtain Riccati differential equation of second order

$\begin{equation*} \left\{ \begin{array}{c} \mathfrak{t}^2\frac{d^2\mathfrak{x}(\mathfrak{t})}{d\mathfrak{t}^2} -(\mathfrak{t}^2-k)\mathfrak{x}(\mathfrak{t}) = \mathfrak{t}^2 \mathfrak{q}( \mathfrak{x}(\mathfrak{t})) \; \; \; \mathfrak{t}\in J, \\\\ \mathfrak{x}(0) = 0\; \; \; \; \mbox{and}\; \; \; \; \; \mathfrak{x}'(0) = 0\end{array} \right. \end{equation*}$

● Taking $\; \mathfrak{g}_1(\mathfrak{t}, \mathfrak{x}, \mathfrak{y}) = -(\mathfrak{t}^2-2l\mathfrak{t}-k)\mathfrak{x}+\mathfrak{q}(\mathfrak{x}), \; k\in R\;$ where $\; \mathfrak{q}(\mathfrak{x})\;$ is continuous function and $\; l\;$ is fixed, then we obtain Coulomb wave differential equation of second order

$\begin{equation*} \left\{ \begin{array}{c} \frac{d^2\mathfrak{x}(\mathfrak{t})}{d\mathfrak{t}^2}+ ( \mathfrak{t}^2-2l\mathfrak{t}-k)\mathfrak{x} = \mathfrak{q}( \mathfrak{x}(\mathfrak{t})) \; \; \; \mathfrak{t}\in J, \\\\ \mathfrak{x}(0) = 0\; \; \; \; \mbox{and}\; \; \; \; \; \mathfrak{x}'(0) = 0\end{array} \right. \end{equation*}$

● Taking $\; \mathfrak{g}_1(\mathfrak{t}, \mathfrak{x}, \mathfrak{y}) = \left(\frac{-8\pi^2m}{\hbar^2} \right)(E \mathfrak{x}-\frac{k\mathfrak{t}^2}{2} \mathfrak{x})+\mathfrak{q}(\mathfrak{x}), \; k \in R\;$ where $\; \mathfrak{q}(\mathfrak{x})\;$ is continuous function and $\; \hbar\;$ is the Planket's constant and $\; E, \; k\;$ are positive real numbers, then we obtain of Schrödinger wave differential equation for simple harmonic oscillator

$\begin{equation*} \left\{ \begin{array}{c} \frac{d^2\mathfrak{x}(\mathfrak{t})}{d\mathfrak{t}^2} = \left(\frac{-8\pi^2m}{\hbar^2} \right)( E \mathfrak{x}(\mathfrak{t})- \frac{k\mathfrak{t}^2}{2} \mathfrak{x}(\mathfrak{t}))+ \mathfrak{q}( \mathfrak{x}(\mathfrak{t})) \; \; \; \mathfrak{t}\in J, \\\\ \mathfrak{x}(0) = 0\; \; \; \; \mbox{and}\; \; \; \; \; \mathfrak{x}'(0) = 0.\end{array} \right. \end{equation*}$

7. Conclusions

Sufficient conditions for the existence of solutions for a class of neutral integro-differential equations of fractional order (1.1) are discussed which involved many key functional differential equations that appear in applications of nonlinear analysis. Also, some types of Ulam stability for this class of implicit fractional differential equation are established. Some applications and particular cases are presented. Finally, the existence of at least one mild solution for this class of equations on an infinite interval by applying Schauder fixed point theorem and the local attractivity of solutions are proved.

Acknowledgements

The authors express their thanks to the anonymous referees for their valuable comments and remarks.

Conflict of interest

The authors declare that they have no competing interests.

Acknowledgments

The authors acknowledge the support of the National Science Foundation Grant Number 1828793. Portions of this work were supported under a contract from the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration. We thank Louis Sumrall for the S. putrefaciens dataset.

Conflict of interest

The authors declare no conflict of interest.

Author contributions:

Acres—processed DHM data, performed tracking, prepared figures, wrote paper.
Nadeau—collected and processed DHM data, wrote paper.

References

[1]	Armitage JP (1997) Three hundred years of bacterial motility. Found Mod Biochem 3: 107-171. doi: 10.1016/S1874-5660(97)80007-X
[2]	Belas R, Zhulin IB, Yang Z (2008) Bacterial signaling and motility: sure bets. J Bacteriol 190: 1849-1856. doi: 10.1128/JB.01943-07
[3]	Holscher T, Bartels B, Lin YC, et al. (2015) Motility, chemotaxis and aerotaxis contribute to competitiveness during bacterial pellicle biofilm development. J Mol Biol 427: 3695-3708. doi: 10.1016/j.jmb.2015.06.014
[4]	Mitchell JG, Kogure K (2006) Bacterial motility: links to the environment and a driving force for microbial physics. FEMS Microbiol Ecol 55: 3-16. doi: 10.1111/j.1574-6941.2005.00003.x
[5]	Velho Rodrigues MF, Lisicki M, Lauga E (2021) The bank of swimming organisms at the micron scale (BOSO-Micro). Plos One 16: e0252291. doi: 10.1371/journal.pone.0252291
[6]	Vaituzis Z, Doetsch RN (1969) Motility tracks: technique for quantitative study of bacterial movement. Appl Microbiol 17: 584-588. doi: 10.1128/am.17.4.584-588.1969
[7]	Barnkob R, Rossi M (2020) General defocusing particle tracking: fundamentals and uncertainty assessment. Exp Fluids 61: 1-14. doi: 10.1007/s00348-020-2937-5
[8]	de Anda J, Lee EY, Lee CK, et al. (2017) High-speed “4D” computational microscopy of bacterial surface motility. Acs Nano 11: 9340-9351. doi: 10.1021/acsnano.7b04738
[9]	Lofroth M, Avci E (2018) Auto-focusing approach on multiple micro objects using the prewitt operator. Int J Intell Robot Appl 2: 413-424. doi: 10.1007/s41315-018-0070-x
[10]	Taute KM, Gude S, Tans SJ, et al. (2015) High-throughput 3D tracking of bacteria on a standard phase contrast microscope. Nat Commun 6: 8776. doi: 10.1038/ncomms9776
[11]	Qi M, Gong X, Wu B, et al. (2017) Landing dynamics of swimming bacteria on a polymeric surface: Effect of surface properties. Langmuir 33: 3525-3533. doi: 10.1021/acs.langmuir.7b00439
[12]	Yuan S, Qi M, Peng Q, et al. (2021) Adaptive behaviors of planktonic Pseudomonas aeruginosa in response to the surface-deposited dead siblings. Colloid Surface B 197: 111408. doi: 10.1016/j.colsurfb.2020.111408
[13]	Cheong FC, Wong CC, Gao YF, et al. (2015) Rapid, high-throughput tracking of bacterial motility in 3D via phase-contrast holographic video microscopy. Biophys J 108: 1248-1256. doi: 10.1016/j.bpj.2015.01.018
[14]	Farthing NE, Findlay RC, Jikeli JF, et al. (2017) Simultaneous two-color imaging in digital holographic microscopy. Opt Express 25: 28489-28500. doi: 10.1364/OE.25.028489
[15]	Flewellen JL, Zaid IM, Berry RM (2019) A multi-mode digital holographic microscope. Rev Sci Instrum 90: 023705. doi: 10.1063/1.5066556
[16]	Hook AL, Flewellen JL, Dubern JF, et al. (2019) Simultaneous tracking of Pseudomonas aeruginosa motility in liquid and at the solid-liquid interface reveals differential roles for the flagellar stators. Msystems 4: e00390-19.
[17]	Huang G, Tian WZ, Qi M, et al. (2018) Improving axial resolution for holographic tracking of colloids and bacteria over a wide depth of field by optimizing different factors. Opt Express 26: 9920-9930. doi: 10.1364/OE.26.009920
[18]	Rieu M, Vieille T, Radou G, et al. (2021) Parallel, linear, and subnanometric 3D tracking of microparticles with Stereo Darkfield Interferometry. Sci Adv 7: eabe3902. doi: 10.1126/sciadv.abe3902
[19]	Yamato K, Chiba H, Oku H (2020) High speed three dimensional tracking of swimming cell by synchronous modulation between tece camera and tag lens. IEEE Robot Autom Lett 5: 1907-1914. doi: 10.1109/LRA.2020.2969911
[20]	Bedrossian M, El-Kholy M, Neamati D, et al. (2018) A machine learning algorithm for identifying and tracking bacteria in three dimensions using Digital Holographic Microscopy. AIMS Biophys 5: 36-49. doi: 10.3934/biophy.2018.1.36
[21]	Cerda M, Navarro CA, Silva J, et al. (2018) A high-speed tracking algorithm for dense granular media. Comput Phys Commun 227: 8-16. doi: 10.1016/j.cpc.2018.02.010
[22]	Shao S, Mallery K, Kumar SS, et al. (2020) Machine learning holography for 3D particle field imaging. Opt Express 28: 2987-2999. doi: 10.1364/OE.379480
[23]	Wang GC, Huang G, Gong XJ, et al. (2020) Method for 3D tracking behaviors of interplaying bacteria individuals. Opt Express 28: 28060-28071. doi: 10.1364/OE.401032
[24]	Wang A, Garmann RF, Manoharan VN (2016) Tracking E-coli runs and tumbles with scattering solutions and digital holographic microscopy. Opt Express 24: 23719-23725. doi: 10.1364/OE.24.023719
[25]	Sheng J, Malkiel E, Katz J, et al. (2007) Digital holographic microscopy reveals prey-induced changes in swimming behavior of predatory dinoflagellates. P Natl Acad Sci USA 104: 17512-17517. doi: 10.1073/pnas.0704658104
[26]	Sheng J, Malkiel E, Katz J, et al. (2010) A dinoflagellate exploits toxins to immobilize prey prior to ingestion. P Natl Acad Sci USA 107: 2082-2087. doi: 10.1073/pnas.0912254107
[27]	Molaei M, Sheng J (2014) Imaging bacterial 3D motion using digital in-line holographic microscopy and correlation-based de-noising algorithm. Opt Express 22: 32119-32137. doi: 10.1364/OE.22.032119
[28]	Wang A, Garmann RF, Manoharan VN (2016) Tracking E. coli runs and tumbles with scattering solutions and digital holographic microscopy. Opt Express 24: 23719-23725. doi: 10.1364/OE.24.023719
[29]	Yu X, Hong J, Liu C, et al. (2014) Review of digital holographic microscopy for three-dimensional profiling and tracking. Opt Eng 53: 112306. doi: 10.1117/1.OE.53.11.112306
[30]	Stocker R (2011) Reverse and flick: Hybrid locomotion in bacteria. P Natl Acad Sci USA 108: 2635-2636. doi: 10.1073/pnas.1019199108
[31]	Bubendorfer S, Held S, Windel N, et al. (2012) Specificity of motor components in the dual flagellar system of Shewanella putrefaciens CN-32. Mol Microbiol 83: 335-350. doi: 10.1111/j.1365-2958.2011.07934.x
[32]	Butler SM, Camilli A (2004) Both chemotaxis and net motility greatly influence the infectivity of Vibrio cholerae. P Natl Acad Sci USA 101: 5018-5023. doi: 10.1073/pnas.0308052101
[33]	Butler SM, Camilli A (2005) Going against the grain: chemotaxis and infection in Vibrio cholerae. Nat Rev Microbiol 3: 611-620. doi: 10.1038/nrmicro1207
[34]	Paździor E, Pękala-Safińska A, Wasyl D (2019) Phenotypic diversity and potential virulence factors of the Shewanella putrefaciens group isolated from freshwater fish. J Vet Res 63: 321-332. doi: 10.2478/jvetres-2019-0046
[35]	Pękala A, Kozińska A, Paździor E, et al. (2015) Phenotypical and genotypical characterization of Shewanella putrefaciens strains isolated from diseased freshwater fish. J Fish Dis 38: 283-293. doi: 10.1111/jfd.12231
[36]	Vignier N, Barreau M, Olive C, et al. (2013) Human infection with Shewanella putrefaciens and S. algae: report of 16 cases in Martinique and review of the literature. Am J Trop Med Hyg 89: 151-156. doi: 10.4269/ajtmh.13-0055
[37]	Kühn J, Niraula B, Liewer K, et al. (2014) A Mach-Zender digital holographic microscope with sub-micrometer resolution for imaging and tracking of marine micro-organisms. Rev Sci Instrum 85: 123113. doi: 10.1063/1.4904449
[38]	Nadeau JL, Cho YB, Kühn J, et al. (2016) Improved tracking and resolution of bacteria in holographic microscopy using dye and fluorescent protein labeling. Front Chem 4: 17. doi: 10.3389/fchem.2016.00017
[39]	Fregoso SF, Lima F, Wallace JK, et al. (2020) DHMx (Digital holographic microscope experience) software tool, Digital Holography and Three-Dimensional Imaging. Opt Soc Am HF2G.7.
[40]	Mann CJ, Kim MK (2006) Quantitative phase-contrast microscopy by angular spectrum digital holography, Three-Dimensional and Multidimensional Microscopy: Image Acquisition and Processing XIII. International Society for Optics and Photonics 6090: 60900B.
[41]	Cohoe D, Hanczarek I, Wallace JK, et al. (2019) Multiwavelength imaging and unwrapping of protozoa in amplitude and phase using custom Fiji plug-ins. Frontiers Phys 7: 94. doi: 10.3389/fphy.2019.00094
[42]	Colomb T, Kühn J, Charrière F, et al. (2006) Total aberrations compensation in digital holographic microscopy with a reference conjugated hologram. Opt Express 14: 4300-4306. doi: 10.1364/OE.14.004300
[43]	Schindelin J, Arganda-Carreras I, Frise E, et al. (2012) Fiji: an open-source platform for biological-image analysis. Nat Methods 9: 676-682. doi: 10.1038/nmeth.2019
[44]	Taute KM, Gude S, Tans SJ, et al. (2015) High-throughput 3D tracking of bacteria on a standard phase contrast microscope. Nat Commun 6: 8776. doi: 10.1038/ncomms9776
[45]	Vater SM, Weiße S, Maleschlijski S, et al. (2014) Swimming behavior of Pseudomonas aeruginosa studied by holographic 3D tracking. PLoS One 9: e87765. doi: 10.1371/journal.pone.0087765
[46]	Heydt M, Divós P, Grunze M, et al. (2009) Analysis of holographic microscopy data to quantitatively investigate three-dimensional settlement dynamics of algal zoospores in the vicinity of surfaces. Eur Phys J E 30: 141-148. doi: 10.1140/epje/i2009-10459-9
[47]	Dixon L, Cheong FC, Grier DG (2011) Holographic deconvolution microscopy for high-resolution particle tracking. Opt Express 19: 16410-16417. doi: 10.1364/OE.19.016410
[48]	Latychevskaia T, Fink HW (2014) Holographic time-resolved particle tracking by means of three-dimensional volumetric deconvolution. Opt Express 22: 20994-21003. doi: 10.1364/OE.22.020994
[49]	Goto T, Nakata K, Baba K, et al. (2005) A fluid-dynamic interpretation of the asymmetric motion of singly flagellated bacteria swimming close to a boundary. Biophys J 89: 3771-3779. doi: 10.1529/biophysj.105.067553
[50]	Lauga E, Diluzio WR, Whitesides GM, et al. (2006) Swimming in circles: Motion of bacteria near solid boundaries. Biophys J 90: 400-412. doi: 10.1529/biophysj.105.069401
[51]	Das D, Lauga E (2019) Transition to bound states for bacteria swimming near surfaces. Phys Rev E 100: 043117. doi: 10.1103/PhysRevE.100.043117
[52]	Kudo S, Imai N, Nishitoba M, et al. (2005) Asymmetric swimming pattern of Vibrio alginolyticuscells with single polar flagella. FEMS Microbiol Lett 242: 221-225. doi: 10.1016/j.femsle.2004.11.007
[53]	Chen M, Zhao Z, Yang J, et al. (2017) Length-dependent flagellar growth of Vibrio alginolyticus revealed by real time fluorescent imaging. ELife 6: e22140. doi: 10.7554/eLife.22140
[54]	Molaei M, Barry M, Stocker R, et al. (2014) Failed escape: Solid surfaces prevent tumbling of Escherichia coli. Phys Rev Lett 113: 068103. doi: 10.1103/PhysRevLett.113.068103
[55]	Bianchi S, Saglimbeni F, Frangipane G, et al. (2019) 3D dynamics of bacteria wall entrapment at a water–air interface. Soft Matter 15: 3397-3406. doi: 10.1039/C9SM00077A
[56]	Berg HC (1971) How to track bacteria. Rev Sci Instrum 42: 868-871. doi: 10.1063/1.1685246
[57]	Frymier PD, Ford RM, Berg HC, et al. (1995) Three-dimensional tracking of motile bacteria near a solid planar surface. P Natl Acad Sci USA 92: 6195-6199. doi: 10.1073/pnas.92.13.6195
[58]	Bianchi S, Saglimbeni F, Di Leonardo R (2017) Holographic imaging reveals the mechanism of wall entrapment in swimming bacteria. Phys Rev X 7: 011010.
[59]	Tokárová V, Perumal AS, Nayak M, et al. (2021) Patterns of bacterial motility in microfluidics-confining environments. P Natl Acad Sci USA 118: e2013925118. doi: 10.1073/pnas.2013925118
[60]	Conrad JC, Poling-Skutvik R (2018) Confined flow: Consequences and implications for bacteria and biofilms. Annu Rev Chem Biomol Eng 9: 175-200. doi: 10.1146/annurev-chembioeng-060817-084006
[61]	Drescher K, Dunkel J, Cisneros LH, et al. (2011) Fluid dynamics and noise in bacterial cell-cell and cell-surface scattering. P Natl Acad Sci USA 108: 10940-10945. doi: 10.1073/pnas.1019079108
[62]	Morse M, Huang A, Li G, et al. (2013) Molecular adsorption steers bacterial swimming at the air/water interface. Biophys J 105: 21-28. doi: 10.1016/j.bpj.2013.05.026
[63]	Mitchell JG, Kogure K (2006) Bacterial motility: links to the environment and a driving force for microbial physics. FEMS Microbiol Ecol 55: 3-16. doi: 10.1111/j.1574-6941.2005.00003.x
[64]	Donlan RM (2002) Biofilms: Microbial life on surfaces. Emerg Infect Dis J 8: 881. doi: 10.3201/eid0809.020063
[65]	Kearns DB (2010) A field guide to bacterial swarming motility. Nat Rev Microbiol 8: 634-644. doi: 10.1038/nrmicro2405
[66]	Guttenplan SB, Kearns DB (2013) Regulation of flagellar motility during biofilm formation. FEMS Microbiol Rev 37: 849-871. doi: 10.1111/1574-6976.12018

biophy-08-04-030-video1.avi

biophy-08-04-030-video2.avi

This article has been cited by:

1.	Xuming Chen, Jianfa Zhu, Liangxiao Li, Chengwen Long, Uniqueness of system integration scheme of artificial intelligence technology in fractional differential mathematical equation, 2022, 0, 2444-8656, 10.2478/amns.2022.2.0104
2.	Abha Singh, Abdul Hamid Ganie, Mashael M. Albaidani, Antonio Scarfone, Some New Inequalities Using Nonintegral Notion of Variables, 2021, 2021, 1687-9139, 1, 10.1155/2021/8045406
3.	Liang Song, Shaodong Chen, Guoxin Wang, Oscillation Analysis Algorithm for Nonlinear Second-Order Neutral Differential Equations, 2023, 11, 2227-7390, 3478, 10.3390/math11163478

Reader Comments

Your name:*

Email:*
© 2021 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Biophysics

1.0 2.2

Metrics

Article views(3196) PDF downloads(144) Cited by(4)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Biophysics

2D vs 3D tracking in bacterial motility analysis

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Main results

4. Stability of solutions of the IVP (1.1)

4.1. Ulam-Hyers stability

4.2. Ulam-Hyers-Rassias stability

5. Existence and attractivity of solutions on half line

6. Applications

7. Conclusions

Acknowledgements

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Abstract

1. Introduction

2. Preliminaries

3. Main results

4. Stability of solutions of the IVP (1.1)

4.1. Ulam-Hyers stability

4.2. Ulam-Hyers-Rassias stability

5. Existence and attractivity of solutions on half line

6. Applications

7. Conclusions

Acknowledgements

Conflict of interest

References

AIMS Biophysics

2D vs 3D tracking in bacterial motility analysis

Related Papers:

Abstract

1. Introduction

2. Preliminaries

3. Main results

4. Stability of solutions of the IVP (1.1)

4.1. Ulam-Hyers stability

4.2. Ulam-Hyers-Rassias stability

5. Existence and attractivity of solutions on half line

6. Applications

7. Conclusions

Acknowledgements

Conflict of interest

References

This article has been cited by:

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

Abstract

1. Introduction

2. Preliminaries

3. Main results

4. Stability of solutions of the IVP (1.1)

4.1. Ulam-Hyers stability

4.2. Ulam-Hyers-Rassias stability

5. Existence and attractivity of solutions on half line

6. Applications

7. Conclusions

Acknowledgements

Conflict of interest

References