Hilfer–Hadamard Fractional Boundary Value Problems with Nonlocal Mixed Boundary Conditions

Abstract

This paper is concerned with the existence and uniqueness of solutions for a Hilfer–Hadamard fractional differential equation, supplemented with mixed nonlocal (multi-point, fractional integral multi-order and fractional derivative multi-order) boundary conditions. The existence of a unique solution is obtained via Banach contraction mapping principle, while the existence results are established by applying the fixed point theorems due to Krasnoselskiĭ and Schaefer and Leray–Schauder nonlinear alternatives. We demonstrate the application of the main results by presenting numerical examples. We also derive the existence results for the cases of convex and non-convex multifunctions involved in the multi-valued analogue of the problem at hand.

1. Introduction

Fractional calculus is regarded as the generalization of the integer-order integration and differentiation in the sense that it deals with derivative and integral operators of an arbitrary real or complex order. This branch of mathematical analysis gained much importance during the last few decades owing to its widespread applications in a variety of disciplines, such as mechanical engineering, bioengineering, biology, physics, chemistry, economics, viscoelasticity, acoustics, optics, robotics, control theory, electronics, etc. The main reason for the popularity of fractional calculus is that mathematical models based on fractional-order operators are considered to be more realistic than the ones relying on classical calculus as such operators are nonlocal in nature and can trace the history of the phenomena under consideration. For the theoretical development of the subject, we refer the reader to the monographs [1,2,3,4,5,6,7,8,9] and the references therein. For some recent applications of fractional calculus concerning structural mechanics and, more specifically, nonlocal elasticity, see [10,11,12].

Fractional-order boundary value problems constitute an important and interesting area of research. It reflects from the literature on the topic that a good deal of work on fractional differential equations involve either Caputo or Riemann–Liouville fractional derivatives. However, these derivatives are found to be inappropriate in the study of some engineering problems. In order to tackle such inaccuracies, some new fractional-order derivative operators such as Hadamard, Erdeyl–Kober, Katugampola, etc., were proposed. In [13], Hilfer introduced a new derivative, which is known as the Hilfer fractional derivative and can generalize both Riemann–Liouville and Caputo derivatives. For some applications of this derivative, we refer the interested reader to the investigations [14,15]. For some recent results on initial and boundary value problems involving Hilfer fractional derivative, for instance, see [16,17,18,19,20,21,22] and the references cited therein.

The fractional derivative presented by Hadamard in 1892 [23] differs from the well known Caputo derivative in two significant ways: (i) its kernel involves a logarithmic function with an arbitrary exponent and (ii) the Hadamard derivative of a constant is not zero. One can find applications of the Hadamard derivative and integral operators in the paper [24] and the monograph [2]. The Hadamard calculus can be obtained by changing $d/dt\to td/dt$, ${(t-s)}^{(·)}\to {({log}_{e}t-{log}_{e}s)}^{(·)}$ and $ds\to (1/s)ds$ in Riemann–Liouville and Caputo fractional derivatives. Later, the modification of Hilfer fractional derivative resulted in the concept of the Hilfer–Hadamard derivative.

Existence results for Hilfer–Hadamard fractional differential equations of order in $(0,1]$ were studied by several researchers, for instance, see [25,26,27]. To the best of our knowledge, only a few results are available in the literature concerning boundary value problems for Hilfer–Hadamard fractional differential equations of order in $(1,2].$ Recently, in [28], the authors applied the tools of the fixed-point theory to study the existence and uniqueness of solutions for a boundary value problem of Hilfer–Hadamard fractional differential equations with nonlocal integro-multi-point boundary conditions:

$\left\{\begin{array}{cc}\hfill & {}^{HH}{D}_1^{\alpha ,\beta }x\left(t\right)=f(t,x\left(t\right)),t\in [1,T],\hfill \\ \hfill & {\displaystyle x\left(1\right)=0,\sum _{i=1}^m{\theta }_{i}x\left({\xi }_i\right)=\lambda {}^H{I}_1^{\delta }x\left(\eta \right),}\hfill \end{array}\right.$

where ${}^{HH}{D}_1^{\alpha ,\beta }$ denote the Hilfer–Hadamard fractional derivative operator of order $\alpha \in (1,2]$ and type $\beta \in [0,1]$, ${\theta }_i,\lambda \in \mathbb{R}$ and $i=1,2,\dots ,m$, are given constants, $f:[1,T]\times \mathbb{R}\to \mathbb{R}$ is a given continuous function and ${}^H{I}^{\delta }$ is the Hadamard fractional integral of order $\delta >0$ and $\eta ,{\xi }_i\in (1,T)$, $i=1,2,\dots ,m.$

In [29], the existence of solutions for the following system of sequential fractional differential equations involving Hilfer–Hadamard type differential operators ${}_H{\mathfrak{D}}^{.,.}$ of different orders was discussed:

$\left\{\begin{array}{c}{(}_H{\mathfrak{D}}_{1^{+}}^{{\alpha }_1,{\beta }_1}+{\lambda }_1{}_H{\mathfrak{D}}_{1^{+}}^{{\alpha }_1-1,{\beta }_1})u\left(t\right)=f(t,u\left(t\right),v\left(t\right)),1<{\alpha }_1\le 2,t\in [1,e],\hfill \\ {(}_H{\mathfrak{D}}_{1^{+}}^{{\alpha }_2,{\beta }_2}+{\lambda }_2{}_H{\mathfrak{D}}_{1^{+}}^{{\alpha }_2-1,{\beta }_2})v\left(t\right)=g(t,u\left(t\right),v\left(t\right)),1<{\alpha }_2\le 2,t\in [1,e],\hfill \\ u\left(1\right)=0,u\left(e\right)=A_1,A_1\in {\mathbb{R}}_{+},\hfill \\ v\left(1\right)=0,v\left(e\right)=A_2,A_2\in {\mathbb{R}}_{+},\hfill \end{array}\right.$

where ${\lambda }_1,{\lambda }_2\in {\mathbb{R}}_{+}$ and f $g:[1,e]\times {\mathbb{R}}^2\to \mathbb{R}$ are given continuous functions.

Motivated by the aforementioned work, our goal in this paper is to enrich the literature on boundary value problems of Hilfer–Hadamard fractional differential equations of order in $(1,2].$ In precise terms, we introduce and study a nonlocal mixed Hilfer–Hadamard boundary value problem of the following form:

(1)$\left\{\begin{array}{cc}\hfill & {}^{HH}{D}_1^{\alpha ,\beta }x\left(t\right)=f(t,x\left(t\right)),t\in [1,T],\hfill \\ \hfill & {\displaystyle x\left(1\right)=0,x\left(T\right)=\sum _{j=1}^m{\eta }_{j}x\left({\xi }_j\right)+\sum _{i=1}^n{\zeta }_i{}^H{I}_1^{{\varphi }_i}x\left({\theta }_i\right)+\sum _{k=1}^r{\lambda }_k{}_H{D}_1^{{\omega }_k}x\left({\mu }_k\right),}\hfill \end{array}\right.$

where ${}^{HH}{D}_1^{\alpha ,\beta }$ denotes the Hilfer–Hadamard fractional derivative operator of order $\alpha \in (1,2]$ and type $\beta \in [0,1]$ and ${\eta }_j,{\zeta }_i,{\lambda }_k\in \mathbb{R}$ are given constants, $f:[1,T]\times \mathbb{R}\to \mathbb{R}$ is a given continuous function, ${}^H{I}^{{\varphi }_i}$ is the Hadamard fractional integral operator of order ${\varphi }_i>0$ and ${\xi }_j,{\theta }_i,{\mu }_k\in (1,T)$, $j=1,2,\dots ,m$, $i=1,2,\dots ,n,$$k=1,2,\dots ,r.$ We also study the multi-valued analogue of the problem (1).

Concerning the significance of problem (1), we recall that the Hilfer fractional derivative interpolates between the Riemann–Liouville and Caputo derivatives [13]. Analogously, the Hilfer–Hadamard type fractional derivative covers the cases of the Riemann–Liouville–Hadamard and Caputo–Hadamard fractional derivatives. Thus, the present study will be useful for improving the works related to glass forming materials [14], Turbulent Flow Model [30], etc. Furthermore, several results involving the Caputo–Hadamard fractional derivative [31,32,33,34] can be extended to the framework of Hilfer–Hadamard fractional derivative.

It is well known that the nonlocal condition is more appropriate than the local condition (initial and/or boundary) with respect to describing certain features of applied mathematics and physics correctly (see the survey paper [35]). More specifically, the boundary conditions arising in the study of boundary value problems of nonlocal elasticity are always nonlocal in nature. This is due to the fact that the long-range interactions within nonlocal solids give rise to nonlocal traction (force) boundary conditions.

Here, we remark that there are only two articles [28,29] in the literature (to the best of our knowledge) concerning boundary value problems for Hilfer–Hadamard fractional differential equations of the order in $(1,2].$ Much of the known studies in the literature deals with initial value problems of Hilfer–Hadamard fractional differential equations of the order in $(0,1].$ The two classes of problems are entirely different. The methodology employed to study the Hilfer–Hadamard fractional differential equations of the order in $(0,1]$ is different from the one applied to such equations of the order in $(1,2].$ Thus, our main objective in this paper is to enrich the new research area on Hilfer–Hadamard fractional differential equations of the order in $(1,2].$ Moreover, the mixed boundary conditions introduced in the present study are of a more general type and include multi-point, fractional integral multi-order and fractional derivative multi-order contributions.

One can notice that the boundary conditions considered in problem (1) reduce to several special cases such as (i) nonlocal multi-point boundary conditions if we choose all ${\zeta }_i=0,i=1,2,\dots ,n$ and ${\lambda }_k=0,k=1,2,\dots ,r$; (ii) nonlocal Hadamard fractional integral boundary conditions when all ${\eta }_j=0,j=1,2,\dots ,m$ and ${\lambda }_k=0,k=1,2,\dots ,r$; and (iii) nonlocal Hadamard fractional boundary conditions if we take all ${\eta }_j=0,j=1,2,\dots ,m$ and ${\zeta }_i=0,i=0,k=1,2,\dots ,n.$ Likewise, we can consider the combination of nonlocal multipoint and Hadamard fractional integral conditions when we fix all ${\lambda }_k=0,$$k=1,2,\dots ,r$ and so on. Thus, the results presented in this paper are significant as they specialize to the ones associated with several interesting boundary conditions. Another novelty in the present work is concerned with the derivation of the existence results for the Hilfer–Hadamard fractional differential inclusions of the order in $(1,2]$ supplemented with the mixed boundary conditions. Thus, the investigation of single-valued and multi-valued nonlocal nonlinear Hilfer–Hadamard fractional boundary value problems of the order in $(1,2]$ enhances the scope of the literature on the topic.

The remaining part of this manuscript is arranged as follows. Section 2 contains some basic notions and known results of fractional differential calculus. In Section 3, we first prove an auxiliary result that plays a key role in transforming the given problem into a fixed point problem. Then, based on Banach’s contraction mapping principle, we establish the existence of a unique solution for the problem (1). By using the fixed point theorems due to Krasnoselskiĭ and Schaefer and nonlinear alternative of Leray–Schauder type, we prove some existence results for problem (1). Examples illustrating the applicability of the main results are also presented in this section. The existence results for the multi-valued analogue of the problem (1) are obtained in Section 4. Some interesting observations are presented in the last section of the paper.

2. Preliminaries

In this section, we recall some basic concepts.

Definition 1.

(Hadamard fractional integral [2]). Let $f:[a,\infty )\to \mathbb{R}.$ Then, the Hadamard fractional integral of order $\alpha >0$ is defined as follows:

(2) ${}^H{I}_a^{\alpha }f\left(t\right)=\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_a^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{f\left(z\right)}{z}dz,t>a,$

provided that the integral exists, where $log(.)={log}_e(.).$

Definition 2.

(Hadamard fractional derivative [2]). For a function $f:[a,\infty )\to \mathbb{R}$, the Hadamard fractional derivative of order $\alpha >0$ is defined as follows:

(3) ${}_H{D}_a^{\alpha }f\left(t\right)={\delta }^n\left({}^H{I}_a^{n-\alpha }f\right)\left(t\right),n=\left[\alpha \right]+1,$

where ${\delta }^n={\left(t\frac{d}{dt}\right)}^n$ and $\left[\alpha \right]$ denote the integer part of the real number α.

Lemma 1.

[2] If $\alpha >0,\beta >0$ and $0<a<b<\infty ,$ then

•  (i) 

  ${\displaystyle \left({}^H{I}_{a^{+}}^{\alpha }{\left(log\frac{t}{a}\right)}^{\beta -1}\right)\left(x\right)=\frac{\mathsf{\Gamma }\left(\beta \right)}{\mathsf{\Gamma }(\beta +\alpha )}{\left(log\frac{x}{a}\right)}^{\beta +\alpha -1};}$

•  (ii) 

  ${\displaystyle \left({}_H{D}_{a^{+}}^{\alpha }{\left(log\frac{t}{a}\right)}^{\beta -1}\right)\left(x\right)=\frac{\mathsf{\Gamma }\left(\beta \right)}{\mathsf{\Gamma }(\beta -\alpha )}{\left(log\frac{x}{a}\right)}^{\beta -\alpha -1}.}$

In particular, if $\beta =1$ , then the following is the case:

$\left({}_H{D}_{a^{+}}^{\alpha }\right)\left(1\right)=\frac{1}{\mathsf{\Gamma }(1-\alpha )}{\left(log\frac{x}{a}\right)}^{-\alpha }\ne 0,0<\alpha <1.$

Definition 3.

(Hilfer–Hadamard fractional derivative [15]). Let $f\in L^1(a,b)$ and $n-1<\alpha <n,$$0\le \beta \le 1.$ We define the Hilfer–Hadamard fractional derivative of order α and type β for f as follows:

$\begin{array}{ccc}\hfill {(}^{HH}{D}_a^{\alpha ,\beta }f)\left(t\right)& =& {(}^H{I}_a^{\beta (n-\alpha )}{\delta }^n{}^H{I}_a^{(n-\alpha )(1-\beta )}f)\left(t\right)\hfill \\ & =& {(}^H{I}_a^{\beta (n-\alpha )}{\delta }^n{}^H{I}_a^{(n-\gamma )}f)\left(t\right)\hfill \\ & =& {(}^H{I}_a^{\beta (n-\alpha )}{}_H{D}_a^{\gamma }f)\left(t\right),\gamma =\alpha +n\beta -\alpha \beta ,\hfill \end{array}$

where ${}^H{I}_a^{(.)}$ and ${}_H{D}_a^{(.)}$ are defined by (2) and (3), respectively.

Here, we remark that the Hilfer–Hadamard fractional derivative reduces to the Hadamard fractional derivative for $\beta =0$ and corresponds to the Caputo–Hadamard derivative for $\beta =1$ given in the following equation:

${}_H^C{D}_a^{\alpha }f\left(t\right)=\left({}^H{I}_a^{n-\alpha }{\delta }^{n}f\right)\left(t\right),n=\left[\alpha \right]+1.$

Next, we recall the following known theorem that will be used in the sequel.

Theorem 1.

([5]). Let $\alpha >0$, $0\le \beta \le 1$, $\gamma =\alpha +n\beta -\alpha \beta$, $n=\left[\alpha \right]+1$ and $0<a<b<\infty$. If $f\in L^1(a,b)$ and ${(}^H{I}_a^{n-\gamma }f)\left(t\right)\in A{C}_{\delta }^n[a,b]$, then the following is the case:

$\begin{array}{ccc}\hfill {}^H{I}_a^{\alpha }{(}^{HH}{D}_a^{\alpha ,\beta }f)\left(t\right)& =& {}^H{I}_a^{\gamma }{(}^{HH}{D}_a^{\gamma }f)\left(t\right)\hfill \\ & =& f\left(t\right)-\sum _{j=0}^{n-1}\frac{\left({\delta }^{(n-j-1)}{(}^H{I}_a^{n-\gamma }f)\right)\left(a\right)}{\mathsf{\Gamma }(\gamma -j)}{\left(log\frac{t}{a}\right)}^{\gamma -j-1}.\hfill \end{array}$

Observe that $\mathsf{\Gamma }(\gamma -j)$ exists for all $j=1,2,\dots ,n-1$ for $\gamma \in [\alpha ,n].$

3. Main Results

This section is concerned with the existence and uniqueness of solutions for the nonlinear Hilfer–Hadamard fractional boundary value problem (1). First of all, we prove an auxiliary lemma dealing with the linear variant of the boundary value problem (1), which will be used to transform the problem at hand into an equivalent fixed point problem. In the case $n=\left[\alpha \right]+1=2$, we have $\gamma =\alpha +(2-\alpha )\beta$.

Lemma 2.

Let $h\in C\left(\right[1,T],\mathbb{R})$ and that

(4) $\begin{array}{ccc}\hfill \mathsf{\Lambda }& =& {(logT)}^{\gamma -1}-\sum _{j=1}^m{\eta }_j{(log{\xi }_j)}^{\gamma -1}-\sum _{i=1}^n{\zeta }_i\frac{\mathsf{\Gamma }\left(\gamma \right)}{\mathsf{\Gamma }(\gamma +{\varphi }_i)}{\left(log{\theta }_i\right)}^{\gamma +{\varphi }_i-1}\hfill \\ & & -\sum _{k=1}^r{\lambda }_k\frac{\mathsf{\Gamma }\left(\gamma \right)}{\mathsf{\Gamma }(\gamma -{\omega }_k)}{(log{\mu }_k)}^{\gamma -{\omega }_k-1}\ne 0.\hfill \end{array}$

Then, x is a solution of the following linear Hilfer–Hadamard fractional boundary value problem:

(5) $\left\{\begin{array}{cc}\hfill & {}^{HH}{D}_1^{\alpha ,\beta }x\left(t\right)=h(t,),t\in [1,T],\hfill \\ \hfill & {\displaystyle x\left(1\right)=0,x\left(T\right)=\sum _{j=1}^m{\eta }_{j}x\left({\xi }_j\right)+\sum _{i=1}^n{\zeta }_i{}^H{I}_1^{{\varphi }_i}x\left({\theta }_i\right)+\sum _{k=1}^r{\lambda }_k{}_H{D}_1^{{\omega }_k}x\left({\mu }_k\right),}\hfill \end{array}\right.$

if and only if it satisfies the integral equation:

(6) $\begin{array}{ccc}\hfill x\left(t\right)& =& {}^H{I}_1^{\alpha }h\left(t\right)+\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m{}^H{I}^{\alpha }h\left({\xi }_j\right)+\sum _{i=1}^n{\zeta }_i{}^H{I}^{\alpha +{\varphi }_i}h\left({\theta }_i\right)\hfill \\ & & +\sum _{k=1}^r{\lambda }_k{}^H{I}^{\alpha -{\omega }_k}h\left({\mu }_k\right)-{}^H{I}^{\alpha }h\left(T\right)\},t\in [1,T].\hfill \end{array}$

Proof. 

Applying the Hadamard fractional integral operator of order α from 1 to t on both sides of Hilfer–Hadamard fractional differential equation in (5) and using Theorem 1, we find that

(7)$x\left(t\right)-\frac{\delta {(}_H{I}_{1^{+}}^{2-\gamma }x)\left(1\right)}{\mathsf{\Gamma }\left(\gamma \right)}{(logt)}^{\gamma -1}-\frac{{(}_H{I}_{1^{+}}^{2-\gamma }x)\left(1\right)}{\mathsf{\Gamma }(\gamma -1)}{(logt)}^{\gamma -2}={}^H{I}_1^{\alpha }h\left(t\right),$

which can be rewritten as follows:

(8)$\begin{array}{cc}\hfill x\left(t\right)=& c_0{(logt)}^{\gamma -1}+c_1{(logt)}^{\gamma -2}+\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t\frac{h\left(s\right)}{s}{\left(log\frac{t}{s}\right)}^{\alpha -1}ds,\hfill \end{array}$

where $c_0$ and $c_1$ are arbitrary constants. Using the first boundary condition ($x\left(1\right)=0$) in (8) yields $c_1=0$, since $\gamma \in [\alpha ,2]$. In consequence, (8) takes the following form:

(9)$\begin{array}{c}\hfill x\left(t\right)=c_0{(logt)}^{\gamma -1}+\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{s}\right)}^{\alpha -1}\frac{h\left(s\right)}{s}ds.\end{array}$

Now, inserting (9) into the second boundary condition:

$x\left(T\right)=\sum _{j=1}^m{\eta }_{j}x\left({\xi }_j\right)+\sum _{i=1}^n{\zeta }_i{}^H{I}_1^{{\varphi }_i}x\left({\theta }_i\right)+\sum _{k=1}^r{\lambda }_k{}_H{D}_1^{{\omega }_k}x\left({\mu }_k\right),$

and using notation (4), we obtain the following:

$c_0=\frac{1}{\mathsf{\Lambda }}\left\{\sum _{j=1}^m{}^H{I}^{\alpha }h\left({\xi }_j\right)+\sum _{i=1}^n{\zeta }_i{}^H{I}^{\alpha +{\varphi }_i}h\left({\theta }_i\right)+\sum _{k=1}^r{\lambda }_k{}^H{I}^{\alpha -{\omega }_k}h\left({\mu }_k\right)-{}^H{I}^{\alpha }h\left(T\right)\right\}.$

Substituting the value of $c_0$ in (9) results in Equation (6) as desired. By direct computation, one can obtain the converse of the lemma. The proof is completed. □

Let $X=C\left([1,T],\mathbb{R}\right)$ be the Banach space endowed with the norm

$\parallel x\parallel :=\underset{t\in [1,T]}{max}\left|x\left(t\right)\right|.$

In view of Lemma 2 and Definition 1, we introduce an operator $\mathcal{F}:X\to X$ associated with the problem (1) as follows:

(10)$\begin{array}{ccc}\hfill \mathcal{F}\left(x\right)\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{f(z,x(z\left)\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{f(z,x(z\left)\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{f(z,x(z\left)\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{f(z,x(z\left)\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{f(z,x(z\left)\right)}{z}dz\},t\in [1,T].\hfill \end{array}$

In the sequel, we use the following notation:

(11)$\begin{array}{ccc}\hfill \mathsf{\Omega }& =& \frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}.\hfill \end{array}$

3.1. Uniqueness Result

Here, by applying Banach’s contraction mapping principle [36], we prove the existence of a unique solution for problem (1).

Theorem 2.

Suppose that the following condition holds:

• $\left(H_1\right)$

  There exists a constant $l>0$ such that for all $t\in [1,T]$ and $u_i\in \mathbb{R},i=1,2,$.

  $|f(t,u_1)-f(t,u_2)|\le l|u_1-u_2|.$

Then, the nonlinear Hilfer–Hadamard fractional boundary value problem (1) has a unique solution on$[1,T]$ if $l\mathsf{\Omega }<1,$ where Ω is defined by (11).

Proof. 

We will verify that the operator $\mathcal{F}$ defined by (10) satisfies the hypotheses of Banach’s contraction mapping principle. Fixing $N={max}_{t\in [1,T]}\left|f(t,0)\right|<\infty$ and using the assumption $\left(H_1\right)$, we obtain the following:

(12)$\begin{array}{c}\hfill \left|f\right(t,x\left(t\right)\left)\right|\le l\left|x\right(t\left)\right|+\left|f\right(t,0\left)\right|\le l\parallel x\parallel +N.\end{array}$

The proof is divided into two steps.

Step I: We show that $\mathcal{F}\left(B_r\right)\subset B_r$, where $B_r=\{x\in X:\parallel x\parallel <r\}$ with $r\ge N\mathsf{\Omega }/(1-l\mathsf{\Omega }).$ Let $x\in B_r.$ Then, we have the following:

$\begin{array}{ccc}\hfill \left|\mathcal{F}\right(x\left)\right(t\left)\right|& \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{|{\eta }_j|}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{|{\zeta }_i|}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\}\hfill \\ & \le & \frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}(l\parallel x\parallel +N)+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}(l\parallel x\parallel +N)\hfill \\ & \le & [\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\left\}\right](lr+N)\hfill \\ & =& \mathsf{\Omega }(lr+N)\le r.\hfill \end{array}$

Thus, the following is the case:

$\parallel \mathcal{F}\left(x\right)\parallel =\underset{t\in [1,T]}{max}\left|\mathcal{F}\left(u\right)\left(t\right)\right|\le r,$

which means that $\mathcal{F}\left(B_r\right)\subset B_r.$

Step II: To show that the operator $\mathcal{F}$ is a contraction, let $x_1,x_2\in X.$ Then, for any $t\in [1,T]$, we have the following:

$\begin{array}{ccc}& & |\mathcal{F}\left(x_2\right)\left(z\right)-\mathcal{F}\left(x_1\right)\left(z\right)|\hfill \\ & \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{|f(z,x_2\left(z\right))-f(z,x_1\left(z\right))|}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{j=1}^m\frac{|{\eta }_j|}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{|f(z,x_2\left(z\right))-f(z,x_1\left(z\right))|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{|{\zeta }_i|}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{|f(z,x_2\left(z\right))-f(z,x_1\left(z\right))|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{|f(z,x_2\left(z\right))-f(z,x_1\left(z\right))|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{|f(z,x_2\left(z\right))-f(z,x_1\left(z\right))|}{z}dz\}\hfill \\ & \le & [\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\left\}\right]l\parallel x_2-x_1\parallel .\hfill \end{array}$

Thus, the following is the case:

$\parallel \mathcal{F}\left(x_2\right)-\mathcal{F}\left(x_1\right)\parallel =\underset{t\in [1,T]}{max}|\mathcal{F}\left(x_2\right)\left(t\right)-\mathcal{F}\left(x_1\right)\left(t\right)|\le l\mathsf{\Omega }\parallel x_2-x_1\parallel ,$

which, in view of $l\mathsf{\Omega }<1$, shows that the operator $\mathcal{F}$ is a contraction. Hence, the operator $\mathcal{F}$ has a unique fixed point by Banach’s contraction mapping principle. Therefore, the problem (1) has a unique solution on $[1,T].$ The proof is completed. □

3.2. Existence Results

In this subsection, we present different criteria for the existence of solutions for the problem (1). First, we prove an existence result based on Krasnoselskiĭ’s fixed point theorem [37].

Theorem 3.

Let $f:[1,T]\times \mathbb{R}\to \mathbb{R}$ be a continuous function satisfying $\left(H_1\right).$ In addition, we assume that the following condition is satisfied:

• $\left(H_2\right)$

  There exists a continuous function $\varphi \in C([1,T],{\mathbb{R}}^{+})$ such that

  $\left|f\right(t,x\left)\right|\le \varphi \left(t\right),for each(t,u)\in [1,T]\times \mathbb{R}.$

Then, the nonlinear Hilfer–Hadamard fractional boundary value problem (1) has at least one solution on $[1,T],$ provided that the following condition holds:

(13) $\begin{array}{ccc}\hfill & & \frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}l<1.\hfill \end{array}$

Proof. 

By assumption $\left(H_2\right)$, we can fix $\rho \ge \mathsf{\Omega }\parallel \varphi \parallel$ and consider a closed ball $B_{\rho }=\{x\in C([1,T],\mathbb{R}):\parallel x\parallel \le \rho \},$ where $\parallel \varphi \parallel ={sup}_{t\in [1,T]}\left|\varphi \left(t\right)\right|$ and $\mathsf{\Omega }$ is given by (11). We verify the hypotheses of Krasnoselskiĭ’s fixed point theorem [37] by splitting the operator $\mathcal{F}$ defined by (10) on $B_{\rho }$ to $C\left(\right[1,T],\mathbb{R})$ as $\mathcal{F}={\mathcal{F}}_1+{\mathcal{F}}_2,$ where ${\mathcal{F}}_1$ and ${\mathcal{F}}_2$ are defined by the following:

$\begin{array}{ccc}\hfill \left({\mathcal{F}}_{1}x\right)\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{f(z,x(z\left)\right)}{z}dz,t\in [1,T],\hfill \end{array}$

$\begin{array}{ccc}\hfill \left({\mathcal{F}}_{2}x\right)\left(t\right)& =& \frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{f(z,x(z\left)\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{f(z,x(z\left)\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{f(z,x(z\left)\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{f(z,x(z\left)\right)}{z}dz\},t\in [1,T].\hfill \end{array}$

For any $x,y\in B_{\rho },$ we have the following:

$\begin{array}{ccc}\hfill |\left({\mathcal{F}}_{1}x\right)\left(t\right)+\left({\mathcal{F}}_{2}y\right)\left(t\right)|& \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{j=1}^m\frac{|{\eta }_j|}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{|{\zeta }_i|}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\}\hfill \\ & \le & [\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\left\}\right]\parallel \varphi \parallel \hfill \\ \hfill & =& \mathsf{\Omega }\parallel \varphi \parallel \le \rho .\hfill \end{array}$

Hence, $\parallel {\mathcal{F}}_{1}x+{\mathcal{F}}_{2}y\parallel \le \rho ,$ which shows that ${\mathcal{F}}_{1}x+{\mathcal{F}}_{2}y\in B_{\rho }.$ By condition (13), it is easy to prove that the operator ${\mathcal{F}}_2$ is a contraction mapping. The operator ${\mathcal{F}}_1$ is continuous by the continuity of $f.$ Moreover, ${\mathcal{F}}_1$ is uniformly bounded on $B_{\rho },$ since

$\parallel {\mathcal{F}}_{1}x\parallel \le \frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\parallel \varphi \parallel .$

Finally, we prove the compactness of the operator ${\mathcal{F}}_1.$ For $t_1,t_2\in [1,T],t_1<t_2,$ we have the following case:

$\begin{array}{ccc}& & |{\mathcal{F}}_{1}x\left(t_2\right)-{\mathcal{F}}_{1}x\left(t_1\right)|\hfill \\ & \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{t_1}\left[{\left(log\frac{t_2}{z}\right)}^{\alpha -1}-{\left(log\frac{t_1}{z}\right)}^{\alpha -1}\right]\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_{t_1}^{t_2}{\left(log\frac{t_2}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & \le & \frac{\parallel \varphi \parallel }{\mathsf{\Gamma }(\alpha +1)}\left[2{(log{t}_2-log{t}_1)}^{\alpha }+|{(log{t}_2)}^{\alpha }-{(log{t}_1)}^{\alpha }|\right],\hfill \end{array}$

which tends to zero independently of $x\in B_{\rho },$ as $t_1\to t_2.$ Thus, ${\mathcal{F}}_1$ is equicontinuous. By the application of the Arzelá–Ascoli theorem, we deduce that operator ${\mathcal{F}}_1$ is compact on $B_{\rho }.$ Thus, the hypotheses of Krasnoselskiĭ’s fixed point theorem [37] hold true. In consequence, there exists at least one solution for the nonlinear Hilfer–Hadamard fractional boundary value problem (1) on $[1,T],$ which completes the proof. □

Our next existence result is based on Schaefer’s fixed point theorem [38].

Theorem 4.

Let $f:[1,T]\times \mathbb{R}\to \mathbb{R}$ be a continuous function satisfying the following assumption:

• $\left(H_3\right)$

  There exists a real constant $M>0$ such that for all $t\in [1,T],u\in \mathbb{R},$

  $\left|f\right(t,u\left)\right|\le M.$

Then, there exists at least one solution for the nonlinear Hilfer–Hadamard fractional boundary value problem (1) on $[1,T].$

Proof. 

We will prove that the operator $\mathcal{F}$, defined by (10), has a fixed point by using Schaefer’s fixed point theorem [38]. The proof is given in two steps.

Step I.We show that the operator $\mathcal{F}:X\to X$ is completely continuous.

Let us first establish that $\mathcal{F}$ is continuous. Let $\left\{x_n\right\}$ be a sequence such that $x_n\to x$ in X. Then, for each $t\in [1,T]$, we have the following:

$\begin{array}{ccc}& & |\mathcal{F}\left(x_n\right)\left(t\right)-\mathcal{F}\left(x\right)\left(t\right)|\hfill \\ & \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{|f(z,x_n\left(z\right))-f(z,x\left(z\right))|}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{j=1}^m\frac{|{\eta }_j|}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{|f(x,x_n\left(z\right))-f(z,x\left(z\right))|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{|{\zeta }_i|}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{|f(z,x_n\left(z\right))-f(z,x\left(z\right))|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{|f(z,x_n\left(z\right))-f(z,x\left(z\right))|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{|f(z,x_n\left(z\right))-f(z,x\left(z\right))|}{z}dz\}.\hfill \end{array}$

Taking into account the fact that f is continuous, that is, $|f(s,x_n\left(s\right))-f(s,x\left(s\right))|\to 0as{x}_n\to x,$ we obtain from the foregoing inequality that the following is the case:

$\parallel \mathcal{F}\left(x_n\right)-\mathcal{F}\left(x\right)\parallel \to 0as{x}_n\to x.$

Hence, $\mathcal{F}$ is continuous.

Now we show that the operator $\mathcal{F}$ which maps bounded sets into bounded sets in $X.$ For $R>0$, let $B_R=\{x\in X:\parallel x\parallel \le R\}.$ Then, for $t\in [1,T]$, we have the following case:

$\begin{array}{ccc}\hfill \left|\mathcal{F}\right(x\left)\right(t\left)\right|& \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{j=1}^m\frac{|{\eta }_j|}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{|{\zeta }_i|}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\}\hfill \\ & \le & \frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}M+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}M,\hfill \end{array}$

which, after taking the norm for $t\in [1,T],$ results in the following inequality:

$\begin{array}{ccc}\hfill \parallel \mathcal{F}\left(x\right)\parallel & \le & \frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}M+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}M.\hfill \end{array}$

Next, we show that bounded sets are mapped into equicontinuous sets by $\mathcal{F}.$ For $t_1,t_2\in [1,T],$$t_1<t_2$ and $u\in B_R,$ we obtain the following:

$\begin{array}{ccc}& & |\mathcal{F}\left(x\right)\left(t_2\right)-\mathcal{F}\left(x\right)\left(t_1\right)|\hfill \\ & \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{t_1}\left[{\left(log\frac{t_2}{z}\right)}^{\alpha -1}-{\left(log\frac{t_1}{z}\right)}^{\alpha -1}\right]\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_{t_1}^{t_2}{\left(log\frac{t_2}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{|{(log{t}_2)}^{\gamma -1}-{(log{t}_1)}^{\gamma -1}|}{\left|\mathsf{\Lambda }\right|}\{\sum _{j=1}^m\frac{|{\eta }_j|}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{|{\zeta }_i|}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\}\hfill \\ & \le & \frac{M}{\mathsf{\Gamma }(\alpha +1)}\left[2{(log{t}_2-log{t}_1)}^{\alpha }+|{(log{t}_2)}^{\alpha }-{(log{t}_1)}^{\alpha }|\right]\hfill \\ & & +\frac{|{(log{t}_2)}^{\gamma -1}-{(log{t}_1)}^{\gamma -1}|}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}M,\hfill \end{array}$

which tends to zero independently of $x\in B_R,$ as $t_1\to t_2.$ Thus, the operator $\mathcal{F}:X\to X$ is completely continuous by applying the Arzelá–Ascoli theorem.

Step II.: We show that the set $\mathcal{E}=\{x\in X|x=\nu \mathcal{F}\left(x\right),0\le \nu \le 1\}$ is bounded. Let $x\in \mathcal{E}$, then $x=\nu \mathcal{F}\left(x\right)$. For any $t\in [1,T]$, we have $x\left(t\right)=\nu \mathcal{F}\left(x\right)\left(t\right)$. Then, in view of the hypothesis $\left(H_3\right)$, as in Step I, we obtain

$\begin{array}{ccc}\hfill \left|x\right(t\left)\right|& \le & \frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}M+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}M.\hfill \end{array}$

Thus, the following is the case:

$\begin{array}{ccc}\hfill \parallel x\parallel & \le & \frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}M+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}M,\hfill \end{array}$

which shows that the set $\mathcal{E}$ is bounded. Thus, it follows by Schaefer’s fixed point theorem [38] that the operator $\mathcal{F}$ has at least one fixed point. Therefore, there exists at least one solution for the nonlinear Hilfer–Hadamard fractional boundary value problem (1) on $[1,T].$ This completes the proof. □

We apply the Leray–Schauder nonlinear alternative [39] to prove our last existence result.

Theorem 5.

Let $f\in C\left(\right[1,T]\times \mathbb{R},\mathbb{R}).$ In addition, it is assumed that the following conditions are satisfied:

• $\left(H_4\right)$

  There exist $p\in C([1,T],{\mathbb{R}}^{+})$ and a continuous nondecreasing function $\psi :{\mathbb{R}}^{+}$$\to {\mathbb{R}}^{+}$ such that $\left|f\right(t,u\left)\right|\le p\left(t\right)\psi (\parallel u\parallel )for each(t,u)\in [1,T]\times \mathbb{R};$

• $\left(H_5\right)$

  There exists a constant $K>0$ such that

  $\frac{K}{\mathsf{\Omega }\parallel p\parallel \psi \left(K\right)}>1,$

  where Ω is defined by (11).

Then, the nonlinear Hilfer–Hadamard fractional boundary value problem (1) has at least one solution on $[1,T].$

Proof. 

As argued in Theorem 4, one can obtain that the operator $\mathcal{F}$ is completely continuous. Next, we establish that we can find an open set $U\subseteq C\left(\right[1,T],\mathbb{R})$ with $x\ne \mu \mathcal{F}\left(x\right)$ for $\mu \in (0,1)$ and $x\in \partial U.$

Let $x\in C\left(\right[1,T],\mathbb{R})$ be such that $x=\mu \mathcal{F}\left(x\right)$ for some $0<\mu <1.$ Then, for each $t\in [1,T],$ we have the following case:

$\begin{array}{ccc}\hfill \left|x\right(t\left)\right|& \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{j=1}^m\frac{|{\eta }_j|}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{|{\zeta }_i|}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\}\hfill \\ & \le & [\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\left\}\right]\parallel p\parallel \psi (\parallel x\parallel ).\hfill \end{array}$

Consequently, we obtain

$\frac{\parallel x\parallel }{\mathsf{\Omega }\parallel p\parallel \psi (\parallel x\parallel )}\le 1.$

In view of $\left(H_5\right)$, there is no solution x such that $\parallel x\parallel \ne K$. Let us set the following:

$U=\{x\in C([1,T],\mathbb{R}):\parallel x\parallel <K\}.$

The operator $\mathcal{F}$$:\overline{U}\to C([1,T],\mathbb{R})$ is continuous and completely continuous. Note that there is no $u\in \partial U$ such that $x=\mu \mathcal{F}\left(x\right)$ for some $\mu \in (0,1),$ by the choice of $U.$ Thus, it follows by the Leray–Schauder nonlinear alternative [39] that $\mathcal{F}$ has a fixed point $x\in \overline{U}$ which is a solution of the nonlinear Hilfer–Hadamard fractional boundary value problem (1). This ends the proof. □

3.3. Examples

Consider the following Hilfer–Hadamard fractional boundary value problem:

(14)$\left\{\begin{array}{cc}\hfill & {}^{HH}{D}_1^{\alpha ,\beta }x\left(t\right)=f(t,x\left(t\right)),t\in [1,T],\hfill \\ \hfill & {\displaystyle x\left(1\right)=0,x\left(T\right)=\sum _{j=1}^m{\eta }_{j}x\left({\xi }_j\right)+\sum _{i=1}^n{\zeta }_i{}^H{I}_1^{{\varphi }_i}x\left({\theta }_i\right)+\sum _{k=1}^r{\lambda }_k{}_H{D}_1^{{\omega }_k}x\left({\mu }_k\right),}\hfill \end{array}\right.$

with $\alpha =5/3,\beta =3/4,T=5,m=4,n=3,r=2,{\eta }_1=1/15,{\eta }_2=1/10,{\eta }_3=2/15,{\eta }_4=1/6,{\xi }_1=5/4,{\xi }_2=3/2,{\xi }_3=7/4,{\xi }_4=2,{\zeta }_i=1/18,{\zeta }_2=1/9,{\zeta }_3=1/6,{\varphi }_1=1/2,{\varphi }_2=1,{\varphi }_3=3/2,{\theta }_1=5/2,{\theta }_2=3,{\theta }_3=7/2,{\lambda }_1=1/28,{\lambda }_2=1/14,{\omega }_1=1/4,{\omega }_2=2/3,{\mu }_1=4,{\mu }_2=9/2$ and $f(t,x(t\left)\right)$ to be fixed later. Using the given data, it is found that $\gamma =23/12,\mathsf{\Lambda }\approx 0.662923$ ($\mathsf{\Lambda }$ is given by (4)), $\mathsf{\Omega }\approx 39.388095$ ($\mathsf{\Omega }$ is given by (11)).

(a). For illustrating Theorem 2, we take

(15)$f(t,x)=\frac{a{e}^{-{(t-1)}^2}}{15}arctanx+\frac{sint+1}{\sqrt{t^3+1}},t\in [1,5],$

where a is a positive real constant. Obviously, the nonlinear function $f(t,x)$ satisfies the assumption $\left(H_1\right)$ with $l=a/15$ and the condition $l\mathsf{\Omega }<1$ holds for $a<0.3808257.$ Thus, the hypothesis of Theorem 2 is satisfied; hence, the problem (14) with $f(t,x)$ given by (15) has a unique solution on $[1,5].$

(b). In order to illustrate Theorem 3, we take the following function:

(16)$f(t,x)=\frac{\left|x\right|}{[{(t-1)}^2+100](1+|x\left|\right)}+cost,t\in [1,5].$

It is easy to verify that the nonlinear function $f(t,x)$ satisfies the assumption $\left(H_1\right)$ with $l=1/100$ and that

$\begin{array}{ccc}& & \frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}l\approx 0.379190<1.\hfill \end{array}$

Thus, the assumptions of Theorem 3 hold true. Therefore, by the conclusion of Theorem 3, there exists at least one solution for the problem (14) with $f(t,x)$ given by (16) on $[1,5].$

(c). Now, we illustrate Theorem 4 with the aid of the following nonlinear function:

(17)$f(t,x)=\frac{e^{1-t}}{t^2+4}{cos}^4\left(\frac{3+2\left|x\right|}{1+{\left|x\right|}^2}\right)+\frac{1}{10}\sqrt{(t^2+11)},t\in [1,5].$

It is easy to obtain that $\left|f\right(t,x\left)\right|\le 4/5.$ Thus, by the conclusion of Theorem 4, the problem (14) with $f(t,x)$ given by (17) has at least one solutions on $[1,5].$

(d). Finally, we demonstrate the application of Theorem 5 by considering the nonlinear function:

(18)$f(t,x)=\frac{2}{[{(t-1)}^2+100]}\left(\left|x\right|cos(5+2\left|x\right|)+\frac{1}{2}\right),t\in [1,5].$

Note that $\left|f\right(t,x\left)\right|\le p\left(t\right)\psi (\parallel x\parallel ),$ where $p\left(t\right)=2{[{(t-1)}^2+100]}^{-1}$ ($\parallel p\parallel =1/50$) and $\psi (\parallel x\parallel )=\parallel x\left|\right|+1/2.$ By the condition $\left(H_5\right),$ we find that $K>1.85584.$ Thus, all the assumptions of Theorem 5 hold true; hence, its conclusion ensures the existence of at least one solution for the problem (14) with $f(t,x)$ given by (18) on $[1,5].$

4. Multi-Valued Case

This section is devoted to the study of the multi-valued case of the boundary value problem (1) given as follows:

(19)$\left\{\begin{array}{cc}\hfill & {}^{HH}{D}_1^{\alpha ,\beta }x\left(t\right)\in F(t,x\left(t\right)),t\in [1,T],\hfill \\ \hfill & {\displaystyle x\left(1\right)=0,x\left(T\right)=\sum _{j=1}^m{\eta }_{j}x\left({\xi }_j\right)+\sum _{i=1}^n{\zeta }_i{}^H{I}_1^{{\varphi }_i}x\left({\theta }_i\right)+\sum _{k=1}^r{\lambda }_k{}_H{D}_1^{{\omega }_k}x\left({\mu }_k\right),}\hfill \end{array}\right.$

where the symbols are the same as defined in problem (1) and $F:J\times \mathbb{R}\to \mathcal{P}\left(\mathbb{R}\right)$ is a multi-valued map. By $\mathcal{P}\left(\mathbb{R}\right)$, we denote the family of all nonempty subsets of $\mathbb{R}$.

Next, we define

${\mathcal{P}}_p=\{Y\in \mathcal{P}\left(X\right):Y\ne \varnothing \mathrm{and} \mathrm{has} \mathrm{the} \mathrm{property}p\},$

where $(X,\parallel ·\parallel )$ is a normed space. Thus, ${\mathcal{P}}_{cl},$${\mathcal{P}}_b,$${\mathcal{P}}_{cp}$ and ${\mathcal{P}}_{cp,c}$ respectively, denote the classes of all closed, bounded, compact and compact and convex sets in $X.$

Define the set of selections of F for each $\omega \in C\left(\right[1,T],\mathbb{R})$ as

$S_{F,\omega }:=\{z\in L^1([1,T],\mathbb{R}):z\left(t\right)\in F(t,\omega \left(t\right))fora.e.t\in [1,T]\}.$

4.1. Existence Results for the Problem (19)

Let us first define the solution for Hilfer–Hadamard inclusions fractional boundary value problem (19).

Definition 4.

A function $x\in C\left(\right[1,T],\mathbb{R})$ is called a solution of the multi-valued problem (19) if we can find a function $v\in L^1([1,T],\mathbb{R})$ with $v\left(t\right)\in F(t,x)$ almost everywhere on $[1,T]$ such that

$\begin{array}{ccc}\hfill x\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz+\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\},\hfill \end{array}$

where Λ is given by (4).

4.1.1. Case 1: Convex-Valued Multifunctions

In the first theorem, dealing with convex-valued multifunctions, we assume that the multifunction F is $L^1$-Carathéodory and apply the nonlinear alternative for Kakutani maps [39] and closed graph operator theorem [40] to prove it.

Theorem 6.

Assume that the following conditions are satisfied:

• $\left(A_1\right)$

  The multifunction $F:[1,T]\times \mathbb{R}\to {\mathcal{P}}_{cp,c}\left(\mathbb{R}\right)$ is $L^1$-Carathéodory;

• $\left(A_2\right)$

  There exist a nondecreasing function $\chi \in C\left(\right[0,\infty \left)\right(0,\infty \left)\right)$ and a continuous function $q:[1,T]\to {\mathbb{R}}^{+}$ such that

  ${\parallel F(t,\omega )\parallel }_{\mathcal{P}}:=sup\left\{\right|z|:z\in F(t,\omega )\}\le q\left(t\right)\chi (\parallel \omega \parallel )for each(t,\omega )$ $\in [1,T]\times \mathbb{R};$

• $\left(A_3\right)$

  There exists $M>0$ satisfying the following inequality:

  ${\displaystyle \frac{M}{{\displaystyle \chi \left(M\right)\parallel q\parallel \mathsf{\Omega }}}>1,}$

  where Ω is given by (11).

Then, there exists at least one solution for the inclusions problem (19) on $[1,T].$

Proof. 

Associated with the Hilfer–Hadamard fractional inclusions boundary value problem (19), we introduce a multi-valued operator, $N:C\left(\right[1,T],\mathbb{R})\to \mathcal{P}\left(C\right([1,T],\mathbb{R}\left)\right)$, as follows:

$N\left(x\right)=\left\{\begin{array}{cc}& h\in C\left(\right[1,T],\mathbb{R}):\hfill \\ & h\left(t\right)=\left\{\begin{array}{c}{\displaystyle \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz}\hfill \\ {\displaystyle +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz}\hfill \\ {\displaystyle +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v\left(z\right)}{z}dz}\hfill \\ {\displaystyle +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v\left(z\right)}{z}dz}\hfill \\ {\displaystyle -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\},v\in S_{F,x}.}\hfill \end{array}\right.\hfill \end{array}\right\}$

In what follows, we will prove in several steps that the operator N satisfies the hypotheses of the Leray–Schauder nonlinear alternative for Kakutani maps [39].

Step 1. N is bounded on bounded sets of $C\left(\right[1,T],\mathbb{R}).$

For a fixed $r>0,$ let $B_r=\{x\in C[1,T],\mathbb{R}):\parallel x\parallel \le r\}$ be a bounded set in $C\left(\right[1,T],\mathbb{R})$. For each $h\in N\left(x\right)$ and $x\in B_r,$ there exists $v\in S_{F,x}$ such that

$\begin{array}{ccc}\hfill h\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz+\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\},t\in [1,T].\hfill \end{array}$

For $t\in [1,T],$ using the assumption $\left(A_2\right)$, we obtain the following inequality:

$\begin{array}{ccc}\hfill \left|h\right(t\left)\right|& \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|v\right(z\left)\right|}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{|{\eta }_j|}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{\left|v\right(z\left)\right|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{|{\zeta }_i|}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{\left|v\right(z\left)\right|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{\left|v\right(z\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|v\right(z\left)\right|}{z}dz\}\hfill \\ & \le & [\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\left\}\right]\parallel p\parallel \chi (\parallel x\parallel ),\hfill \end{array}$

which yields

$\parallel h\parallel \le \parallel p\parallel \chi \left(r\right)\mathsf{\Omega }.$

Step 2. Bounded sets are mapped by N into equicontinuous sets of $C\left(\right[1,T],\mathbb{R}).$

Let $x\in B_r$ and $h\in N\left(x\right).$ Then, there exists $v\in S_{F,x}$ such that

$\begin{array}{ccc}\hfill h\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz+\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\},t\in [1,T].\hfill \end{array}$

Let $t_1,t_2\in [1,T],t_1<t_2$. Then, we obtain

$\begin{array}{ccc}& & |\mathcal{F}\left(x\right)\left(t_2\right)-\mathcal{F}\left(x\right)\left(t_1\right)|\hfill \\ & \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{t_1}\left[{\left(log\frac{t_2}{z}\right)}^{\alpha -1}-{\left(log\frac{t_1}{z}\right)}^{\alpha -1}\right]\frac{\left|f\right(z,x\left(z\right)\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_{t_1}^{t_2}{\left(log\frac{t_2}{z}\right)}^{\alpha -1}\frac{\left|v\right(z\left)\right|}{z}dz\hfill \\ & & +\frac{|{(log{t}_2)}^{\gamma -1}-{(log{t}_1)}^{\gamma -1}|}{\left|\mathsf{\Lambda }\right|}\{\sum _{j=1}^m\frac{|{\eta }_j|}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{\left|v\right(z\left)\right|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{|{\zeta }_i|}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{\left|v\right(z\left)\right|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{\left|v\right(z\left)\right|}{z}dz\hfill \\ & & +\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{\left|v\right(z\left)\right|}{z}dz\}\hfill \\ & \le & \frac{\parallel p\parallel \chi \left(r\right)}{\mathsf{\Gamma }(\alpha +1)}\left[2{(log{t}_2-log{t}_1)}^{\alpha }+|{(log{t}_2)}^{\alpha }-{(log{t}_1)}^{\alpha }|\right]\hfill \\ & & +\frac{|{(log{t}_2)}^{\gamma -1}-{(log{t}_1)}^{\gamma -1}|}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\}\parallel p\parallel \chi \left(r\right)\to 0,\hfill \end{array}$

as $t_1\to t_2$ independently of $x\in B_r.$ Hence, $N:C\left(\right[1,T],\mathbb{R})\to \mathcal{P}\left(C\right[1,T],\mathbb{R}))$ is completely continuous by Arzelá–Ascoli theorem.

Step 3. For each $x\in C\left(\right[1,T],\mathbb{R}),$$N\left(x\right)$ is convex.

For $h_1,h_2\in N\left(x\right),$ there exist $v_1,v_2\in S_{F,x}$ such that, for each $t\in [1,T]$, we have

$\begin{array}{ccc}\hfill h_{\varpi }\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v_{\varpi }\left(z\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v_{\varpi }\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v_{\varpi }\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v_{\varpi }\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v_{\varpi }\left(z\right)}{z}dz\},\varpi =1,2.\hfill \end{array}$

Let $0\le \sigma \le 1$. Then, for each $t\in [1,T]$, we have the following:

$\begin{array}{ccc}& & [\sigma h_1+(1-\sigma )h_2]\left(t\right)\hfill \\ & =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{[\sigma v_1\left(z\right)+(1-\sigma )v_2\left(z\right)]}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{[\sigma v_1\left(z\right)+(1-\sigma )v_2\left(z\right)]}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\sigma }_i}{\left(log\frac{{\sigma }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{[\sigma v_1\left(z\right)+(1-\sigma )v_2\left(z\right)]}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{[\sigma v_1\left(z\right)+(1-\sigma )v_2\left(z\right)]}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{[\sigma v_1\left(z\right)+(1-\sigma )v_2\left(z\right)]}{z}dz\}.\hfill \end{array}$

Since $S_{F,x}$ is convex (F has convex values), $\sigma h_1+(1-\sigma )h_2\in N\left(x\right).$ In consequence, N is convex-valued.

Next, it will be shown that the operator N is upper semicontinuous. By using the fact that a completely continuous operator, which has a closed graph, is upper semicontinuous ([41] (Proposition 1.2)), it is enough to prove that the operator N has a closed graph. This will be established in the following step.

Step 4. The graph of N is closed.

Let $x_n\to x_{*}$, $h_n\in N\left(x_n\right)$ and $h_n\to h_{*}$. Then, we show that $h_{*}\in N\left(x_{*}\right)$. Observe that $h_n\in N\left(x_n\right)$ implies that there exists $v_n\in S_{F,x_n}$ such that, for each $t\in [1,T]$, we have

$\begin{array}{ccc}\hfill h_n\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v_n\left(z\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v_n\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v_n\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v_n\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v_n\left(z\right)}{z}dz\}.\hfill \end{array}$

For each $t\in [1,T]$, we must have $v_{*}\in S_{F,x_{*}}$ and the following expression:

$\begin{array}{ccc}\hfill h_{*}\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v_{*}\left(z\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v_{*}\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v_{*}\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v_{*}\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v_{*}\left(z\right)}{z}dz\}.\hfill \end{array}$

Consider the continuous linear operator $\mathsf{\Phi }:L^1([1,T],\mathbb{R})\to C([1,T],\mathbb{R})$ defined as follows:

$\begin{array}{ccc}\hfill v\to \mathsf{\Phi }\left(v\right)\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\}.\hfill \end{array}$

Clearly $\parallel h_n-h_{*}\parallel \to 0$ as $n\to \infty ,$ and consequently, by the closed graph operator theorem [40], $\mathsf{\Phi }\circ S_{F,x}$ is a closed graph operator. Moreover, we have $h_n\in \mathsf{\Phi }\left(S_{F,x_n}\right)$ and the following:

$\begin{array}{ccc}\hfill h_{*}\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v_{*}\left(z\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v_{*}\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v_{*}\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v_{*}\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v_{*}\left(z\right)}{z}dz\},\hfill \end{array}$

for some $v_{*}\in S_{F,x_{*}}.$ Thus, N has a closed graph, which implies that the operator N is upper semicontinuous.

Step 5. There exists an open set $U\subseteq C\left(\right[1,T],\mathbb{R})$ such that, for any $k\in (0,1)$ and all $x\in \partial U,$ $x\notin kN\left(x\right).$

Let $x\in kN\left(x\right),$$k\in (0,1).$ Then, there exists $v\in L^1([1,T],\mathbb{R})$ with $v\in S_{F,x}$ such that, for $t\in [1,T]$, we have the following case.

$\begin{array}{ccc}\hfill x\left(t\right)& =& k\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +k\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\}.\hfill \end{array}$

Following the computation as in Step 2, for each $t\in [1,T],$ we have the following inequality:

$\begin{array}{ccc}\hfill \left|x\right(t\left)\right|& \le & [\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\left\}\right]\parallel p\parallel \chi (\parallel x\parallel )\hfill \\ & =& \parallel p\parallel \chi (\parallel x\parallel )\mathsf{\Omega }.\hfill \end{array}$

In consequence, we obtain

${\displaystyle \frac{\parallel x\parallel }{{\displaystyle \chi (\parallel x\parallel )\parallel p\parallel \mathsf{\Omega }}}\le 1.}$

By assumption $\left(A_3\right)$, we can find M such that $\parallel x\parallel \ne M$. Let us set the following:

$U=\{x\in C([1,T],\mathbb{R}):\parallel x\parallel <M\}.$

Notice that $N:\overline{U}\to \mathcal{P}\left(C\right([1,T],$$\mathbb{R}\left)\right)$ is a compact multi-valued map with convex closed values, which is upper semicontinuous; moreover, from the choice of $U,$ there is no $x\in \partial U$ such that $x\in kN\left(x\right)$ for some $k\in (0,1).$ Hence, we deduce by the Leray–Schauder nonlinear alternative for Kakutani maps [39] that N has a fixed point $x\in \overline{U}.$ This implies the existence of at least one solution for the inclusions problem (19) on $[1,T].$ The proof is complete. □

4.1.2. Case 2: Nonconvex Valued Multifunctions

Here, we prove an existence result for the Hilfer–Hadamard fractional inclusions boundary value problem (19) with a non-convex valued multi-valued map via a fixed point theorem for multivalued maps due to Covitz and Nadler [42].

Definition 5. ([43]) Let $(X,d)$ be a metric space induced from the normed space $(X;\parallel ·\parallel )$ and $H_d:\mathcal{P}\left(X\right)\times \mathcal{P}\left(X\right)\to \mathbb{R}\cup \{\infty \}$ be defined as follows:

$H_d(A,B)=max\{\underset{a\in A}{sup}d(a,B),\underset{b\in B}{sup}d(A,b)\},$

where $d(A,b)={inf}_{a\in A}d(a,b)$ and $d(a,B)={inf}_{b\in B}d(a,b)$.

Theorem 7.

Let the following conditions hold:

• $\left(B_1\right)$

  $F:[1,T]\times \mathbb{R}\to {\mathcal{P}}_{cp}\left(\mathbb{R}\right)$ is such that $F(·,x):[1,T]\to {\mathcal{P}}_{cp}\left(\mathbb{R}\right)$ is measurable for each $x\in \mathbb{R}$;

• $\left(B_2\right)$

  $H_d(F(t,x),F(t,\overline{x}))\le \varrho \left(t\right)|x-\overline{x}|$ for almost all $t\in [1,T]$ and $x,\overline{x}\in \mathbb{R}$ with $\varrho \in C([1,T],{\mathbb{R}}^{+})$ and $d(0,F(t,0\left)\right)\le \varrho \left(t\right)$ for almost all $t\in [1,T]$.

Then, the Hilfer–Hadamard inclusions fractional boundary value problem (19) has at least one solution on $[1,T]$ if

$\mathsf{\Omega }\parallel \varrho \parallel <1,$

where Ω is given by (11).

Proof. 

We verify that the operator $N:C\left(\right[1,T],\mathbb{R})\to \mathcal{P}\left(C\right([1,T],\mathbb{R}\left)\right),$ defined at the beginning of the proof of Theorem 6, satisfies the hypotheses of the fixed point theorem for multivalued maps due to Covitz and Nadler [42].

Step I.N is nonempty and closed for every $v\in S_{F,x}.$

It follows by the measurable selection theorem ([44], (Theorem III.6) that the set-valued map $F(·,x(·\left)\right)$ is measurable; hence, it admits a measurable selection $v:[1,T]\to \mathbb{R}$. In view of the assumption $\left(B_2\right),$ we obtain $\left|v\right(t\left)\right|\le \varrho \left(t\right)(1+|x\left(t\right)\left|\right),$ that is, $v\in L^1([1,T],\mathbb{R})$; hence, F is integrably bounded. In consequence, we deduce that $S_{F,x}\ne \varnothing$.

Now, we show that $N\left(x\right)\in {\mathcal{P}}_{cl}\left(C\right([1,T],$$\mathbb{R}\left)\right)$ for each $x\in C\left(\right[1,T],\mathbb{R})$. For that, let ${\left\{u_n\right\}}_{n\ge 0}\in N\left(x\right)$ with $u_n\to u(n\to \infty )$ in $C\left(\right[1,T],\mathbb{R}).$ Then, $u\in C\left(\right[1,T],\mathbb{R})$, and we can find $v_n\in S_{F,x_n}$ satisfying the following equation:

$\begin{array}{ccc}\hfill u_n\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v_n\left(z\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v_n\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v_n\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v_n\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v_n\left(z\right)}{z}dz\},\hfill \end{array}$

for each $t\in [1,T].$ Then, we can obtain a sub-sequence (if necessary) $v_n$ converging to v in $L^1([1,T],\mathbb{R}),$ as F has compact values. Thus, $v\in S_{F,x}$ and for each $t\in [1,T]$, we have the following:

$\begin{array}{ccc}\hfill u_n\left(t\right)\to v\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v\left(z\right)}{z}dz\}.\hfill \end{array}$

Thus, $u\in N\left(x\right).$

Step II. In this step, it will be shown that there exists $0<m_0<1$$(m_0=\mathsf{\Omega }\parallel \varrho \parallel )$ such that the following is the case.

$H_d(N\left(x\right),N\left(\overline{x}\right))\le m_0\parallel x-\overline{x}\parallel for eachx,\overline{x}\in C([1,T],\mathbb{R}).$

Let $x,\overline{x}\in C([1,T],\mathbb{R})$ and $h_1\in N\left(x\right)$. Then, there exists $v_1\left(t\right)\in F(t,x\left(t\right))$ such that, for each $t\in [1,T]$, we have

$\begin{array}{ccc}\hfill h_1\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v_1\left(z\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v_1\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v_1\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v_1\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v_1\left(z\right)}{z}dz\}.\hfill \end{array}$

By $\left(B_2\right)$, we have

$H_d(F(t,x),F(t,\overline{x}))\le \varrho \left(t\right)|x\left(t\right)-\overline{x}\left(t\right)|.$

Thus, there exists $\vartheta \left(t\right)\in F(t,\overline{x}\left(t\right))$ such that

$|v_1\left(t\right)-\vartheta |\le \varrho \left(t\right)|x\left(t\right)-\overline{x}\left(t\right)|,t\in [1,T].$

Let us define $\mathcal{V}:[1,T]\to \mathcal{P}\left(\mathbb{R}\right)$ by

$\mathcal{V}\left(t\right)=\{\vartheta \in \mathbb{R}:|v_1\left(t\right)-\vartheta |\le \varrho \left(t\right)|x\left(t\right)-\overline{x}\left(t\right)\left|\right\}.$

Then, there exists a function $v_2\left(t\right)$ that is a measurable selection of $\mathcal{V},$ since the multivalued operator $\mathcal{V}\left(t\right)\cap F(t,\overline{x}\left(t\right))$ is measurable (Proposition III.4 [44]). Hence, $v_2\left(t\right)\in F(t,\overline{x}\left(t\right))$ and for each $t\in [1,T]$, we have $|v_1\left(t\right)-v_2\left(t\right)|\le \varrho \left(t\right)|x\left(t\right)-\overline{x}\left(t\right)|$. Thus, for each $t\in [1,T]$, we have

$\begin{array}{ccc}\hfill h_2\left(t\right)& =& \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{v_2\left(z\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{v_2\left(z\right)}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{v_2\left(z\right)}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{v_2\left(z\right)}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{v_2\left(z\right)}{z}dz\}.\hfill \end{array}$

Hence, we have the following:

$\begin{array}{ccc}\hfill \left|h_1\left(t\right)-h_2\left(t\right)\right|& \le & \frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^t{\left(log\frac{t}{z}\right)}^{\alpha -1}\frac{|v_2\left(z\right)-v_1\left(z\right)\left|\right)}{z}dz\hfill \\ & & +\frac{{(logt)}^{\gamma -1}}{\mathsf{\Lambda }}\{\sum _{j=1}^m\frac{{\eta }_j}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^{{\xi }_j}{\left(log\frac{{\xi }_j}{z}\right)}^{\alpha -1}\frac{|v_2\left(z\right)-v_1\left(z\right)|}{z}dz\hfill \\ & & +\sum _{i=1}^n\frac{{\zeta }_i}{\mathsf{\Gamma }(\alpha +{\varphi }_i)}{\int }_1^{{\theta }_i}{\left(log\frac{{\theta }_i}{z}\right)}^{\alpha +{\varphi }_i-1}\frac{|v_2\left(z\right)-v_1\left(z\right)|}{z}dz\hfill \\ & & +\sum _{k=1}^r\frac{{\lambda }_k}{\mathsf{\Gamma }(\alpha -{\omega }_k)}{\int }_1^{{\mu }_k}{\left(log\frac{{\mu }_k}{z}\right)}^{\alpha -{\omega }_k-1}\frac{|v_2\left(z\right)-v_1\left(z\right)|}{z}dz\hfill \\ & & -\frac{1}{\mathsf{\Gamma }\left(\alpha \right)}{\int }_1^T{\left(log\frac{T}{z}\right)}^{\alpha -1}\frac{|v_2\left(z\right)-v_1\left(z\right)|}{z}dz\}\hfill \\ & \le & [\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\frac{{(logT)}^{\gamma -1}}{\left|\mathsf{\Lambda }\right|}\{\sum _{i=j}^m\frac{|{\eta }_j|{(log{\xi }_j)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}+\sum _{i=1}^n\frac{|{\zeta }_i|{(log{\theta }_i)}^{\alpha +{\varphi }_i}}{\mathsf{\Gamma }(\alpha +{\varphi }_i+1)}\hfill \\ & & +\sum _{k=1}^r\frac{|{\lambda }_k|{(log{\mu }_k)}^{\alpha -{\omega }_k}}{\mathsf{\Gamma }(\alpha -{\omega }_k+1)}+\frac{{(logT)}^{\alpha }}{\mathsf{\Gamma }(\alpha +1)}\left\}\right]\parallel \varrho \parallel \parallel x-\overline{x}\parallel ,\hfill \end{array}$

which yields

$\parallel h_1-h_2\parallel \le \mathsf{\Omega }\parallel \varrho \parallel \parallel x-\overline{x}\parallel .$

On switching the roles of x and $\overline{x}$, we obtain the following case:

$H_d(N\left(x\right),N\left(\overline{x}\right))\le \mathsf{\Omega }\parallel \varrho \parallel \parallel x-\overline{x}\parallel ,$

which shows that N is a contraction. Consequently, the conclusion of the fixed point theorem for multivalued maps due to Covitz and Nadler [42] applies; hence, the operator N has a fixed point x which corresponds to a solution of the Hilfer–Hadamard inclusions fractional boundary value problem (19). The proof is complete. □

4.2. Examples

Let us consider the Hilfer–Hadamard fractional inclusions boundary value problem:

(20)$\left\{\begin{array}{cc}\hfill & {}^{HH}{D}_1^{\alpha ,\beta }x\left(t\right)\in F(t,x\left(t\right)),t\in [1,T],\hfill \\ \hfill & {\displaystyle x\left(1\right)=0,x\left(T\right)=\sum _{j=1}^m{\eta }_{j}x\left({\xi }_j\right)+\sum _{i=1}^n{\zeta }_i{}^H{I}_1^{{\varphi }_i}x\left({\theta }_i\right)+\sum _{k=1}^r{\lambda }_k{}_H{D}_1^{{\omega }_k}x\left({\mu }_k\right),}\hfill \end{array}\right.$

with the values of the parameters taken in the problem (14), while the multi-valued map $F:J\times \mathbb{R}\to \mathcal{P}\left(\mathbb{R}\right)$ will be defined below.

We first illustrate Theorem 6 by taking the following multi-valued map:

(21)$F(t,x\left(t\right))=\left[\frac{e^{-{(1-t)}^2}}{\sqrt{{(1-t)}^4}+75)}\left(\frac{{\left|x\right|}^2}{(1+|x\left|\right)}+\frac{1}{5}\right),\frac{1}{100}\left(x+\frac{cosx}{9}\right)\right].$

It is easy to check that ${\parallel F(t,x)\parallel }_{\mathcal{P}}\le q\left(t\right)\chi (\parallel x\parallel ),$ where $q\left(t\right)$ and $\chi (\parallel x\parallel )$ are given as follows:

$q\left(t\right)=\frac{e^{-{(1-t)}^2}}{\sqrt{{(1-t)}^4}+75)},\chi (\parallel x\parallel )=\parallel x\parallel +\frac{1}{5}.$

Using the values $\parallel q\parallel =1/75,\chi (\parallel x\parallel )=\parallel x\parallel +\frac{1}{5}$ and $\mathsf{\Omega }\approx 39.388095$ (see Section 3.3) in the condition $\left(A_3\right)$, that is, ${\displaystyle \frac{M}{{\displaystyle \chi \left(M\right)\parallel q\parallel \mathsf{\Omega }}}>1,}$ we find that $M>0.221207.$ Thus, all the assumptions of Theorem 6 are satisfied; hence, its conclusion implies that problem (20) with $F(t,x(t\left)\right)$ given by (21) has a solution on $[1,5].$

Next, we demonstrate the application of Theorem 7 by choosing the following multi-valued map:

(22)$F(t,x\left(t\right))=\left[\frac{1+arctanx}{{(20+t)}^2},\left(\frac{\sqrt{t^2+11}}{360}\right)\left(\frac{1+5\left|x\right|}{1+4\left|x\right|}\right)\right].$

Clearly, F is measurable for all $x\in \mathbb{R},$ satisfying the following inequality:

$H_d(F(t,x),F(t,\overline{x}))\le \left(\frac{\sqrt{t^2+11}}{360}\right)|x-\overline{x}|,x,\overline{x}\in \mathbb{R},t\in [1,5].$

Fixing $\varrho \left(t\right)=\sqrt{t^2+11}/360$, we have $\parallel \varrho \parallel =1/60$ and $d(0,\mathcal{F}(t,0\left)\right)\le \varrho \left(t\right)$, $t\in [1,5]$. Furthermore, $\parallel \varrho \parallel \mathsf{\Omega }\approx 0.656468<1.$ Clearly, the hypothesis of Theorem 7 holds true; consequently, we deduce by its conclusion that there exists a solution for the problem (20) with $F(t,x(t\left)\right)$ given by (22) on $[1,5].$

5. Conclusions

In this paper, we have presented the existence and uniqueness criteria for the solutions of a Hilfer–Hadamard fractional differential equation complemented with mixed nonlocal (multi-point, fractional integral multi-order and fractional derivative multi-order) boundary conditions. Firstly, we have converted the given nonlinear problem into a fixed point problem. Once the fixed point operator is available, we can make use of the Banach contraction mapping principle to obtain the uniqueness result. The first two existence results are proved by applying the fixed point theorems due to Schaefer and Krasnoselskiĭ, while the third existence result is based on Leray–Schauder nonlinear alternative. Next, we present the existence results for the corresponding inclusions problem. The first result for the inclusions problem deals with the convex-valued multivalued map and is obtained by applying the Leray–Schauder alternative for multivalued maps, while the non-convex valued multivalued map case relies on the Covitz–Nadler fixed point theorem for contractive multivalued maps. All the results obtained for single and multivalued maps are well illustrated by numerical examples. It is imperative to mention that our results are new in the given configuration and enrich the literature on boundary value problems involving Hilfer–Hadamard fractional differential equations and inclusions of order in $(1,2].$

Author Contributions

Conceptualization, B.A. and S.K.N.; methodology, B.A. and S.K.N.; validation, B.A. and S.K.N.; formal analysis, B.A. and S.K.N.; writing—original draft preparation, B.A. and S.K.N. All authors have read and agreed to the published version of the manuscript.

Funding

This project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, under grant no. (RG-43-130-41).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

This project was funded by the Deanship of Scientific Research (DSR) at King Abdulaziz University, Jeddah, under grant no. (RG-43-130-41). The authors, therefore, acknowledge with thanks DSR technical and financial support. The authors also thank the reviewers for their constructive remarks on our work.

Conflicts of Interest

The authors declare no conflict of interest.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Diethelm, K. The Analysis of Fractional Differential Equations; Lecture Notes in Mathematics; Springer: New York, NY, USA, 2010.

2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of the Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204.

3. Lakshmikantham, V.; Leela, S.; Devi, J.V. Theory of Fractional Dynamic Systems; Cambridge Scientific Publishers: Cambridge, UK, 2009.

4. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Differential Equations; John Wiley: NewYork, NY, USA, 1993.

5. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and Breach Science: Yverdon, Switzerland, 1993.

6. Podlubny, I. Fractional Differential Equations; Academic Press: New York, NY, USA, 1999.

7. Ahmad, B.; Alsaedi, A.; Ntouyas, S.K.; Tariboon, J. Hadamard-Type Fractional Differential Equations, Inclusions and Inequalities; Springer: Cham, Switzerland, 2017.

8. Zhou, Y. Basic Theory of Fractional Differential Equations; World Scientific: Singapore, 2014.

9. Ahmad, B.; Ntouyas, S.K. Nonlocal Nonlinear Fractional-Order Boundary Value Problems; World Scientific: Singapore, 2021.

10. Patnaik, S.; Sidhardh, S.; Semperlotti, F. Towards a unified approach to nonlocal elasticity via fractional-order mechanics. Int. J. Mech. Sci.; 2021; 189, 105992. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2020.105992]

11. Patnaik, S.; Sidhardh, S.; Semperlotti, F. Geometrically nonlinear analysis of nonlocal plates using fractional calculus. Int. J. Mech. Sci.; 2020; 179, 105710. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2020.105710]

12. Stempin, P.; Sumelka, W. Space-fractional Euler-Bernoulli beam model-Theory and identification for silver nanobeam bending. Int. J. Mech. Sci.; 2020; 186, 105902. [DOI: https://dx.doi.org/10.1016/j.ijmecsci.2020.105902]

13. Applications of Fractional Calculus in Physics; Hilfer, R. World Scientific: Singapore, 2000.

14. Hilfer, R. Experimental evidence for fractional time evolution in glass forming materials. J. Chem. Phys.; 2002; 284, pp. 399-408. [DOI: https://dx.doi.org/10.1016/S0301-0104(02)00670-5]

15. Hilfer, R.; Luchko, Y.; Tomovski, Z. Operational method for the solution of fractional differential equations with generalized Riemann-Liouville fractional derivatives. Frac. Calc. Appl. Anal.; 2009; 12, pp. 299-318.

16. Furati, K.M.; Kassim, N.D.; Tatar, N.E. Existence and uniqueness for a problem involving Hilfer fractional derivative. Comput. Math. Appl.; 2012; 64, pp. 1616-1626. [DOI: https://dx.doi.org/10.1016/j.camwa.2012.01.009]

17. Gu, H.; Trujillo, J.J. Existence of mild solution for evolution equation with Hilfer fractional derivative. Appl. Math. Comput.; 2015; 257, pp. 344-354. [DOI: https://dx.doi.org/10.1016/j.amc.2014.10.083]

18. Wang, J.; Zhang, Y. Nonlocal initial value problems for differential equations with Hilfer fractional derivative. Appl. Math. Comput.; 2015; 266, pp. 850-859. [DOI: https://dx.doi.org/10.1016/j.amc.2015.05.144]

19. Sousa, J.V.d.; de Oliveira, E.C. On the ψ-Hilfer fractional derivative. Commun. Nonlinear Sci. Numer. Simul.; 2018; 60, pp. 72-91. [DOI: https://dx.doi.org/10.1016/j.cnsns.2018.01.005]

20. Asawasamrit, S.; Kijjathanakorn, A.; Ntouyas, S.K.; Tariboon, J. Nonlocal boundary value problems for Hilfer fractional differential equations. Bull. Korean Math. Soc.; 2018; 55, pp. 1639-1657.

21. Wongcharoen, A.; Ahmad, B.; Ntouyas, S.K.; Tariboon, J. Three-point boundary value problems for Langevin equation with Hilfer fractional derivative. Adv. Math. Phys.; 2020; 2020, 9606428. [DOI: https://dx.doi.org/10.1155/2020/9606428]

22. Ntouyas, S.K. A survey on existence results for boundary value problems of Hilfer fractional differential equations and inclusions. Foundations; 2021; 1, pp. 63-98. [DOI: https://dx.doi.org/10.3390/foundations1010007]

23. Hadamard, J. Essai sur l’etude des fonctions donnees par leur developpment de Taylor. J. Mat. Pure Appl. Ser.; 1892; 8, pp. 101-186.

24. Butzer, P.L.; Kilbas, A.A.; Trujillo, J.J. Mellin transform analysis and integration by parts for Hadamard-type fractional integrals. J. Math. Anal. Appl.; 2002; 270, pp. 1-15. [DOI: https://dx.doi.org/10.1016/S0022-247X(02)00066-5]

25. Qassim, M.D.; Furati, K.M.; Tatar, N.-E. On a differential equation involving Hilfer-Hadamard fractional derivative. Abstr. Appl. Anal.; 2012; 2012, 391062. [DOI: https://dx.doi.org/10.1155/2012/391062]

26. Vivek, D.; Shahb, K.; Kanagarajan, K. Dynamical analysis of Hilfer-Hadamard type fractional pantograph equations via successive approximation. J. Taibah Univ. Sci.; 2019; 13, pp. 225-230. [DOI: https://dx.doi.org/10.1080/16583655.2018.1558613]

27. Bachira, F.S.; Abbas, S.; Benbachir, M.; Benchohra, M. Hilfer-Hadamard fractional differential equations; Existence and Attractivity, Adv. Theory Nonl. Anal. Appl.; 2021; 5, pp. 49-57.

28. Promsakon, C.; Ntouyas, S.K.; Tariboon, J. Hilfer-Hardamard nonlocal integro-multi-point fractional boundary value problems. J. Funct. Spaces; submitted for publication

29. Saengthong, W.; Thailert, E.; Ntouyas, S.K. Existence and uniqueness of solutions for system of Hilfer-Hadamard sequential fractional differential equations with two point boundary conditions. Adv. Differ. Equ.; 2019; 2019, 525. [DOI: https://dx.doi.org/10.1186/s13662-019-2459-8]

30. Wang, G.; Ren, X.; Zhang, L.; Ahmad, B. Explicit iteration and unique positive solution for a Caputo-Hadamard fractional turbulent flow model. IEEE Access; 2019; 7, pp. 109833-109839. [DOI: https://dx.doi.org/10.1109/ACCESS.2019.2933865]

31. Graef, J.R.; Grace, S.R.; Tunc, E. Asymptotic behavior of solutions of nonlinear fractional equations with Caputo-type Hadamard derivatives. Fract. Calc. Appl. Anal.; 2017; 20, pp. 71-87. [DOI: https://dx.doi.org/10.1515/fca-2017-0004]

32. Abbas, S.; Benchohra, M.; Hamidi, N.; Henderson, J. Caputo-Hadamard fractional differential equations in Banach spaces. Fract. Calc. Appl. Anal.; 2018; 21, pp. 1027-1045. [DOI: https://dx.doi.org/10.1515/fca-2018-0056]

33. Ma, L. Comparison theorems for Caputo-Hadamard fractional differential equations. Fractals; 2019; 27, 1950036. [DOI: https://dx.doi.org/10.1142/S0218348X19500361]

34. Li, C.; Li, Z.; Wang, Z. Mathematical analysis and the local discontinuous Galerkin method for Caputo-Hadamard fractional partial differential equation. J. Sci. Comput.; 2020; 85, 41. [DOI: https://dx.doi.org/10.1007/s10915-020-01353-3]

35. Ntouyas, S.K. Nonlocal Initial and Boundary Value Problems: A survey. Handbook on Differential Equations: Ordinary Differential Equations; Canada, A.; Drabek, P.; Fonda, A. Elsevier Science B.V.: Amsterdam, The Netherlands, 2005; pp. 459-555.

36. Deimling, K. Nonlinear Functional Analysis; Springer: New York, NY, USA, 1985.

37. Krasnoselskiĭ, M.A. Two remarks on the method of successive approximations. Uspekhi Mat. Nauk; 1955; 10, pp. 123-127.

38. Smart, D.R. Fixed Point Theory; Cambridge University Press: Cambridge, UK, 1974.

39. Granas, A.; Dugundji, J. Fixed Point Theory; Springer: New York, NY, USA, 2003.

40. Lasota, A.; Opial, Z. An application of the Kakutani-Ky Fan theorem in the theory of ordinary differential equations, Bull. Acad. Pol. Sci. Ser. Sci. Math. Astronom. Phys.; 1965; 13, pp. 781-786.

41. Deimling, K. Multivalued Differential Equations; De Gruyter: Berlin, Germany, 1992.

42. Covitz, H.; Nadler, S.B., Jr. Multivalued contraction mappings in generalized metric spaces. Israel J. Math.; 1970; 8, pp. 5-11. [DOI: https://dx.doi.org/10.1007/BF02771543]

43. Kisielewicz, M. Differential Inclusions and Optimal Control; Kluwer: Dordrecht, The Netherlands, 1991.

44. Castaing, C.; Valadier, M. Convex Analysis and Measurable Multifunctions; Lecture Notes in Mathematics 580 Springer: Berlin/Heidelberg, Germany, New York, NY, USA, 1977.