\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2016 (2016), No. 198, pp. 1--18.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu}
\thanks{\copyright 2016 Texas State University.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2016/198\hfil Infinitely many positive solutions]
{Infinitely many positive solutions for fractional differential inclusions}

\author[G. Bin, Y.-X. Cui, J.-C. Zhang \hfil EJDE-2016/198\hfilneg]
{Ge Bin, Ying-Xin Cui, Ji-Chun Zhang}

\address{Ge Bin (corresponding author)\newline
Department of Applied Mathematics,
Harbin Engineering University,
Harbin 150001, China}
\email{gebin791025@hrbeu.edu.cn}

\address{Ying-Xin Cui \newline
Department of Applied Mathematics,
Harbin Engineering University,
Harbin 150001, China}
\email{605495064@qq.com}


\address{Ji-Chun Zhang \newline
Department of Applied Mathematics,
Harbin Engineering University,
Harbin 150001, China}
\email{986799294@qq.com}

\thanks{Submitted March 2, 2016. Published July 24, 2016.}
\subjclass[2010]{35A15, 34B15, 58E05, 26A33}
\keywords{Fractional differential inclusions; oscillatory nonlinearities; 
\hfill\break\indent infinitely many solutions;
variational methods; nonsmooth critical point theory}

\begin{abstract}
 In this article, we study a class of fractional differential
 inclusions problem. By nonsmooth variational methods and the theory
 of the fractional derivative spaces, we establish the existence of
 infinitely many positive solutions of the problem under suitable
 oscillatory assumptions on the potential $F$ at zero or at infinity.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{remark}[theorem]{Remark}
\newtheorem{definition}[theorem]{Definition}
\allowdisplaybreaks

\section{Introduction}

In this article, we  consider the existence and
multiplicity of solutions for the fractional
differential inclusion
\begin{equation}
\begin{gathered}
\frac{d}{dt}\Big(\frac{1}{2}\,{_0}D_t^{-\beta}(u'(t))+\frac{1}{2}
\,{_0}D_T^{-\beta}(u'(t))\Big) \in\partial
F(t,u(t)), \quad \text{a.a. }t\in[0,T], \\
 u(0)=u(T)=0,
 \end{gathered} \label{eP}
\end{equation}
where ${_0}D_t^{-\beta}$ and ${_0}D_T^{-\beta}$
are the left and right Riemann-Liouville fractional integrals of
order $0\leq\beta<1$, respectively, $F:[0,T]\times
{\mathbb{R}^{N}}\to \mathbb{R}$ is locally Lipschitz
function in the $t$-variable integrand (in general it can be
nonsmooth), and $\partial F(t,x)$ is the subdifferential with
respect to the $t$-variable in the sense of Clarke \cite{1}.

Fractional differential equations and inclusions have been proved
that they are very valued tools in the modeling of many phenomena in
various fields of science and engineering, such as, viscoelasticity,
electrochemistry, electromagnetism, economics, optimal control,
porous media, etc. In consequence, the subject of fractional
differential equations and inclusions is gaining much importance and
attention. For details and examples, see \cite{31,32,4,35,10},
and the references therein.

Recently, variational methods have turned out to be a very effective
analytical tool in the study of nonlinear problems. The classical
point theory for $C^1$ functional was developed in the sixties and
seventies, see \cite{13,14,11,12}. The need of specific applications
(such as nonsmooth mechanics, nonsmooth gradient systems, etc.) and
the impressive progress in nonsmooth analysis and multivalued
analysis led to extensions of the critical point theory to
nondifferentiable functions, locally Lipschitz functions in
particular. The nonsmooth critical point theory for locally
Lipschitz functions started with the work of Chang \cite{14}. Chang
proposed a generalization of the well-known Palais-Smale condition
and obtained various minimax principles concerning the existence and
characterization of critical points for locally Lipschitz functions.
Chang used his theory to study semilinear elliptic boundary value
problem with a discontinuous nonlinearity.

There are some papers which are devoted to the boundary value
problems for fractional differential inclusion, see
\cite{16,15,18,17}. And the main tools they use are fixed point
theory for multi-valued contractions. In particular, if
$F(x,\cdot)\in C^1(\mathbb{R}^N)$ for a.a. $x \in \mathbb{R}^N$,
then problem \eqref{eP} becomes
\begin{equation}
 \begin{gathered}
\frac{d}{dt}\Big(
\frac{1}{2}\,{_0}D_t^{-\beta}(u'(t))+\frac{1}{2}\,{_0}D_T^{-\beta}(u'(t))\Big)
=\nabla F(t,u(t)), \quad \text{a.a. }t\in[0,T], \\
 u(0)=u(T)=0.
 \end{gathered} \label{eP1}
\end{equation}
Thus a solution $u$ of \eqref{eP} is a weak solution to the problem \eqref{eP1}. So,
in some sense, the solutions of \eqref{eP} can be considered as generalized
solutions of \eqref{eP1}, thus, the formulation of
\eqref{eP} is completely justified.

In the past decade, there are many papers dealing with the existence of
multiple solutions of fractional boundary value
problems \cite{20,23,33,34,19,22,36,21} and the references therein.
 For example, Jiao and Zhou  \cite{19} got one nontrivial solutions for
problem \eqref{eP1} using the mountain pass theorem. Chen and Tang \cite{20}
studied the existence and multiplicity of solutions for the system \eqref{eP1}
when the nonlinearity $F(t,\cdot)$ are superquadratic,
asymptotically quadratic, and subquadratic, respectively.
In \cite{23}, by using the minmax methods in critical point theory,
the authors proved the existence of infinitely many solutions under
suitable conditions. Inspired by the above-mentioned papers,
we study problem \eqref{eP} from a more extensive viewpoint.
So we deal with the existence of infinitely many
solutions for problem \eqref{eP} with the potential $F(x,t)$
exhibits an oscillation at the origin or at infinity.
Indeed, our main results (see Theorems \ref{thm3.1} and \ref{thm3.2} below) give
sufficient conditions on the oscillatory terms such that problem
\eqref{eP} has infinitely many positive solutions. As a byproduct, these
solutions can be constructed in such a way that their norms in
$E^{\alpha}$ tend to zero (to infinity, respectively) whenever the
nonlinearity oscillates at zero (at infinity, respectively).

This article is organized as follows.
In section 2, we present some necessary
preliminary knowledge on the fractional derivative space $E_0^{\alpha,p}$
and generalized gradient of the locally Lipschitz function.
In section 3, we give the main results of this paper.


\section{Preliminaries}

 In this part, we recall some definitions and display
the variational setting which has been established for our problem.

\begin{definition}[\cite{15}] \label{def2.1}
 Let $f(t)$ be a function defined on $[a,b]$ and  $\tau>0$.
The left and right Riemann-Liouville fractional integrals
 of order $\tau$ for function $f(t)$ denoted by
$_{a}D_t^{-\tau}f(t)$ and $_{t}D_b^{-\tau}f(t)$, respectively, are
 defined by
\begin{equation}\label{2.1}
\begin{gathered}
_{a}D_t^{-\tau}f(t)=\frac{1}{\Gamma(\tau)}\int_a^t(t-s)^{\tau-1}f(s)ds,\;t\in[a,b] ,\\
_{t}D_b^{-\tau}f(t)=\frac{1}{\Gamma(\tau)}\int_t^b(t-s)^{\tau-1}f(s)ds,\;t\in[a,b],
\end{gathered}
\end{equation}
provided the right-hand sides are pointwise defined on $[a,b]$, where
$\Gamma$ is the gamma function.
\end{definition}

\begin{definition}[\cite{15}] \label{def2.2} \rm
 Let $f(t)$ be a function defined on  $[a,b]$. The left and right Riemann-Liouville
fractional derivatives of order $\tau$ for function $f(t)$ denoted by
$_{a}D_t^{\tau}f(t)$ and $_{t}D_b^{\tau}f(t)$, respectively, are defined
by
\begin{equation}\label{2.2}
\begin{gathered}
 {}_{a}D_t^{\tau}f(t)=\frac{d^n}{dt^n}\,_{a}
D_t^{\tau-n}f(t)=\frac{1}{\Gamma(n-\tau)}\frac{d^n}{dt^n}
\Big(\int_a^t(t-s)^{n-\tau-1}f(s)ds\Big),\\
{}_{t} D_b^{\tau}f(t)=(-1)^n\frac{d^n}{dt^n}\, _{t}D_b^{\tau-n}f(t)
=\frac{1}{\Gamma(n-\tau)}\frac{d^n}{dt^n}\Big(\int_t^b(t-s)^{n-\tau-1}f(s)ds\Big),
\end{gathered}
\end{equation}
where $t\in[a,b]$, $n-1\leq\tau<n$ and $n\in\mathbb{N}$.
\end{definition}

The left and the right Caputo fractional derivatives are defined via the
above Riemann-Liouville fractional derivatives. In particular, they
are defined for the function belong-ing to the space of absolutely
continuous functions, which we denote by $AC([a,b],\mathbb{R}^N)$.
$AC^k([a,b],\mathbb{R}^N)(k=1,2,\cdots)$ is the space of functions
$f$ such that $f\in C^k([a,b],\mathbb{R}^N)$. In particular,
$AC([a,b],\mathbb{R}^N)=AC^1([a,b],\mathbb{R}^N)$.

\begin{definition}[\cite{15}] \label{def2.3} \rm
 Let $\tau\geq 0$ and $n\in \mathbb{N}$.
If $\tau\in[n-1,n)$ and $f(t)\in AC^n([a,b],\mathbb{R}^N)$, then the
left and right Caputo fractional derivative of order $\tau$ for
function $f(t)$ denoted by $_{a}^c D_t^{\tau}f(t)$ and $_{t}^c
D_b^{\tau}f(t)$, respectively, exist almost everywhere on $[a,b]$.
$_{a}^c D_t^{\tau}f(t)$ and $_{t}^c D_b^{\tau}f(t)$ are represented
by
\begin{equation}\label{2.3}
\begin{gathered}
_{a}^c D_t^{\tau}f(t)=\,_{a}D_t^{\tau-n}f^{(n)}(t)
=\frac{1}{\Gamma(n-\tau)}\Big(\int_a^t(t-s)^{n-\tau-1}f^{(n)}(s)ds\Big) ,\\
_{t}^c D_b^{\tau}f(t)=(-1)\,^n_{t}
D_b^{\tau-n}f^{(n)}(t)
=\frac{1}{\Gamma(n-\tau)}\Big(\int_t^b(t-s)^{n-\tau-1}f^{(n)}(s)ds\Big),
\end{gathered}
\end{equation}
respectively, where $t\in[a,b]$.
\end{definition}

\begin{definition}[\cite{16}] \label{def2.4} \rm
 Define $0<\alpha\leq 1$ and
$1<p<\infty$. The fractional derivative space $E_0^{\alpha,p}$ is
defined by the closure of $C_0^\infty ([0,T],\mathbb{R}^N)$ with
respect to the norm
\begin{equation}\label{2.4}
\begin{aligned}
\|u\|_{\alpha,p}=\Big(\int_0^T|u(t)|^pdt+\int_0^T|_0^cD_t^\alpha
u(t)|^pdt\Big)^{1/p},\quad \forall u\in E_0^{\alpha,p},
\end{aligned}\end{equation}
 where $C_0^\infty([0,T],\mathbb{R}^N)$ denotes the set of all functions
$u\in C^\infty([0,T],\mathbb{R}^N)$ with $u(0)=u(T)=0$.
It is obvious that the fractional
derivative space $E_0^{\alpha,p}$ is the space of functions $u\in
L^p([0,T],\mathbb{R}^N)$ having an $\alpha$-order Caputo fractional
derivative $_{0}^c D_t^{\alpha}u\in L^p([0,T],\mathbb{R}^N)$  and
$u(0)=u(T) = 0$.
\end{definition}

\begin{proposition}[\cite{16}] \label{prop2.1}
Let $0<\alpha\leq 1$ and $1<p<\infty$.
The fractional derivative space $E_0^{\alpha,p}$ is a reflexive and
separable space.
\end{proposition}

\begin{proposition}[\cite{16}] \label{prop2.2} \rm
Let $0<\alpha\leq 1$ and $1<p<\infty$.
For all $u\in E_0^{\alpha,p}$, we have
\begin{equation}\label{2.1b}
\|u\|_{L^p}\leq\frac{T^\alpha}{\Gamma(\alpha+1)}\|_0^c D_t^\alpha
u\|_{L^p}.
\end{equation}
Moreover, if $\alpha>\frac{1}{p}$ and $\frac{1}{p}+\frac{1}{q}=1$,
then
\begin{equation}\label{2.2b}
\|u\|_{\infty}\leq\frac{T^{\frac{\alpha-1}{p}}}
{\Gamma(\alpha)((\alpha-1)q+1)^{1/q}}\|_0^c
D_t^\alpha u\|_{L^p}.
\end{equation}
According to \cite{16}, we can consider $E_0^{\alpha,p}$ with respect to
the norm
\begin{equation}\label{2.2c}
\|u\|_{\alpha,p}=\|_0^c D_t^\alpha u\|_{L^p}=\Big(\int_0^T|_0^c
D_t^\alpha u|^p dt\Big)^{\frac{1}{p}}.
\end{equation}
\end{proposition}

\begin{proposition}[\cite{16}] \label{prop2.3}
 Define $0<\alpha\leq 1$ and $1<p<\infty$. Assume that $\alpha>\frac{1}{p}$ and
 the sequence ${u_k}$ converges weakly to $u\in E_0^{\alpha,p}$, i.e.
$u_k\rightharpoonup u$.
 Then $u_k\to u$ in $C([0,T],\mathbb{R}^N)$, i.e. $\|u_k-u\|_\infty\to0$,
as $k\to\infty$.
\end{proposition}

Using Definition 2.3, for any $u\in
AC([0,T],\mathbb{R}^N)$, problem \eqref{eP} is equivalent to
the problem
\begin{equation}
  \begin{gathered}
\frac{d}{dt}\Big( \frac{1}{2}\,{_0}D_t^{\alpha-1}(_0^c D_t^\alpha
u(t))-\frac{1}{2}\,{_t}D_T^{\alpha-1}(_t^c D_T^\alpha
u(t))\Big)\in\partial F(t,u(t)),\quad \text{a.e. }t\in[0,T],\\
u(0)=u(T)=0,
\end{gathered} \label{eP1b}
\end{equation}
 where $ \alpha=1-\beta\in(\frac{1}{2},1]$. In
the following, we will treat problem \eqref{eP1} in the Hilbert space
$E^\alpha=E_0^{\alpha,2}$ with the corresponding norm
$\|u\|_\alpha=\|u\|_{\alpha,2}$.

\begin{definition}[\cite{16}] \label{def2.5}
 A function $u\in AC([0,T],\mathbb{R}^N)$ is called a solution of \eqref{eP} if
\begin{itemize}
\item[(i)] $D^\alpha(u(t))$ is derivative for almost every $t\in [0,T]$,
and

\item[(ii)] $u$ satisfies \eqref{eP},

\end{itemize}
where $D^\alpha(u(t)):=\frac{1}{2}_0D_t^{\alpha-1}(_0^cD_t^\alpha
u(t))-\frac{1}{2}_tD_T^{\alpha-1}(_t^c D_T^\alpha u(t))$.
\end{definition}

\begin{proposition}[\cite{16}] \label{prop2.4}
If $\frac{1}{2}<\alpha\leq 1$, then for any $u\in E^\alpha$, we have
\begin{equation}\label{2.3b}
|\cos (\pi\alpha)|\|u\|_\alpha^2\leq-\int_0^T\Big(\,_0^cD_t^\alpha
u(t),\,_t^cD_T^\alpha u(t)\Big)dt \leq\frac{1}{|\cos (\pi\alpha)|}\|u\|_\alpha^2.
\end{equation}
\end{proposition}

\begin{proposition}[\cite{16}] \label{prop2.5}
 Let ${1/2}< \alpha \leq 1$ be satisfied. If $u\in E^\alpha$, then the functional
$J: E^\alpha\to \mathbb{R}$ defined by
$$
J(u)=-\frac{1}{2}\int_0^T(\,_0^cD_t^\alpha u(t),\, _t^cD_T^\alpha u(t))dt
$$
is convex and continuous on $E^\alpha$.
\end{proposition}

Let $X$ be a Banach space and $X^{*}$ be its topological dual space
and we denote $\langle\cdot, \cdot\rangle$ as the duality bracket
for pair $(X^{*}, X)$. A function $\varphi : X\mapsto\mathbb{R}$ is
said to be locally Lipschitz, if for every $x\in X$, we can find a
neighbourhood $U$ of $x$ and a constant $k>0$(depending on $U$),
such that $|\varphi(y)-\varphi(z)|\leq k\|y-z\|, \forall y,z\in U$.

 For a locally Lipschitz function $\varphi : X\mapsto\mathbb{R}$ we define
\[
\varphi^{0}(x; h)=\limsup_{x'\to x; \lambda\downarrow
0}\frac{\varphi(x'+\lambda h)-\varphi(x')}{\lambda}.
\]
It is obvious that the function $h \mapsto \varphi^{0}(x; h)$ is
sublinear, continuous and so is the support function of a nonempty,
convex and $w^{*}-$compact set $\partial \varphi(x)\subseteq X^{*}$,
defined by
\[
\partial \varphi(x)=\{x^{*}\in X^{*}; \langle
x^{*}, h\rangle\leq\varphi^{0}(x; h),~\forall h\in X\}.
\]
The multifunction $\partial\varphi: X\mapsto 2^{X^{*}}$ is called
the generalized subdifferential of $\varphi$.

If $\varphi$ is also convex, then $\partial \varphi(x)$ coincides
with subdifferential in the sense of convex analysis, defined by
\[
\partial_{C}\varphi(x)=\{x^{*}\in X^{*}: \langle x^{*}, h\rangle\leq \varphi(x+h)-\varphi(x)\;\;{\rm for }\;\;h\in X\}.
\]
If $\varphi\in C^{1}(X)$, then $\partial\varphi(x)=\{\varphi'(x)\}$.

A point $x\in X$ is a critical point of $\varphi$, if
$0\in\partial\varphi(x)$. It is easily seen that, if $x\in X$ is a
local minimum of $\varphi$, then $0\in\partial\varphi(x)$.

\begin{lemma} \label{lem2.1}
The functional
\begin{equation}\label{2.2d}
\varphi(u)=\int_0^T \big[-\frac{1}{2}(\,_0^c D_t^\alpha u(t),\,_t^c
D_T^\alpha u(t))\big]dt-\int_0^TF(t,u(t))dt
\end{equation}
is locally Lipschitz on $E^\alpha$. Moreover, for $u,v\in E^\alpha$,
we have
\begin{equation}\label{2.2e}
\begin{aligned}
\langle\zeta,v\rangle
&=-\int_0^T\frac{1}{2}\big[(\,_0^c D_t^\alpha
u(t),\,_t^c D_T^\alpha v(t)) +(\,_t^c D_T^\alpha u(t),\,_0^c
D_t^\alpha v(t))\big]dt\\
&\quad -\int_0^T (q(t),v(t)) dt,
\end{aligned}
\end{equation}
where $\zeta \in\partial\varphi(u)$ and
$q(t)\in\partial(F(t,u(t)))$.
\end{lemma}

\begin{proof}
Let $I(u)=\int_0^TF(t,u(t))dt$, then
 $\varphi(u)=J(u)-I(u)$. Obviously, $J(u)$ is locally Lipschitz.
 For $\varepsilon$ is smaller enough, there existent
 $B_\varepsilon(0)\subset\mathbb{N}$.
For any $  {u_1(t), u_2(t)}\in B_\varepsilon(0)$ we have
$$
F(t,u_1(t))-F(t,u_2(t))=\langle\partial F(t,\bar{u}(t)),
 u_1(t)-u_2(t)\rangle,
$$
where $\bar{u}(t)=\lambda  u_1(t)+(1-\lambda)u_2(t)$,
for $\lambda\in(0,1)$.
Furthermore,
$$
\|\bar{u}\|_{E^\alpha}=\|\lambda
 u_1+(1-\lambda)u_2\|_{E^\alpha}\leq\|\lambda
 u_1\|_\alpha+\|(1-\lambda)u_2\|_\alpha
\leq\|u_1\|_\alpha+\|u_2\|_\alpha\leq 2\varepsilon.
$$
Thus, we obtain
\begin{align*}
|I(u_1)-I(u_2)|
&\leq\int_0^T c(1+|\bar{u}(t)|^{\alpha(t)-1})|u_1(t)-u_2(t)|dt\\
&\leq c\int_0^T|u_1(t)-u_2(t)|dt+c\int_0^T||\bar{u}(t)|^{\alpha_1-1}|u_1(t)-u_2(t)|dt\\
&\leq c_1\|u_1-u_2\|_{E^\alpha}+c_2\|\bar{u}\|_{E^\alpha}^{\alpha_1-1}\|u_1-u_2\|_{E^\alpha}\\
&\leq c_1\|u_1-u_2\|_{E^\alpha}+c_2(2\varepsilon)^{\alpha_1-1}\|u_1-u_2\|_{E^\alpha}\\
&\leq L\|u_1-u_2\|_{E^\alpha},
\end{align*}
where $\alpha_1=\min_{t\in[0,T]}\alpha(t)$, and $c_1, c_2$
are positive contents.
\end{proof}

\begin{proposition}[\cite{1}] \label{prop2.6}
 Let $x$ and $y$ be point in Banach space $X$, and suppose that $f$ is
Lipschitz on an open set containing
the line segment $[x,y]$. Then there exists a point $u$ in $(x,y)$
such that
\[
f(y)-f(x)\in\langle \partial f(u),y-x \rangle.
\]
\end{proposition}

\section{Main results and their proofs}

Now we are in a position to state our first main result which deals with
the case when the nonlinearity $F(x,t)$
exhibits an oscillation at the origin.
Our hypotheses on nonsmooth potential $F(x,t)$ are listed as follows.
\smallskip

\noindent\textbf{(H1)}  $F: [0,T]\times\mathbb{R}^{N}\to\mathbb{R}$
is a function, $F(t, 0)=0$ for almost all $t\in[0,T]$ and satisfies
the following facts:
\begin{itemize}
\item[(1)] For all $x\in \mathbb{R}^{N}$, $t\mapsto F(t, x)$ is measurable;

\item[(2)] For almost all $t\in[0,T]$, $x\mapsto F(t,x)$ is locally
Lipschitz;

\item[(3)] There exist a positive constant $c$ such that for almost all
$x\in\mathbb{R}^N$, all $t\in[0,T]$ and $\omega\in\partial F(t,x)$
$$
|\omega|\leq c(1+|x|^{\alpha(t)-1})
$$
where $1<\alpha(t)<+\infty$;

\item[(4)] $-\infty<\liminf_{|x|\to
0^+}\frac{F(t,x)}{|x|^2}\leq\limsup_{|x|\to
0^+}\frac{F(t,x)}{|x|^{2}}=+\infty$ uniformly for a.e. $t\in[0,T]$;

\item[(5)] For every $k\in \mathbb N$, there exists $e_k\in\mathbb{R}^{N}$
with $|e_k|=1$ and there are two sequences $\{a_k\}$ and $\{b_k\}$
in $(0,+\infty)$ with $a_k<b_k$,
$\lim_{k\to+\infty}b_k=0$ such that
$$
\sup\{\omega\cdot e_k:\omega\in\partial F(t,x),
\quad\text{a.e. }t\in[0,T], \;x\in[a_k,b_k]e_k\}\geq 0.
$$
\end{itemize}

\begin{remark} \label{rmk3.1}\rm
 Hypotheses {\rm (H1)(4)} and  {\rm (H1)(5)} imply an oscillatory behaviour
of $F$ near the origin.
\end{remark}

\begin{remark} \label{rmk3.2}  \rm
A simple example of a nonsmooth potential function satisfying
\[
F(t,x)=\begin{cases}
 0, & \text{if } |x|=0\text{ or } |x|\in [\frac{1}{2\pi},+\infty),\\
 |x|^{\beta(t)}\sin \frac{1}{|x|}, & \text{if }
|x|\in[\frac{1}{(2k+1)\pi},\frac{1}{2k\pi}),\\
|x|^{\alpha(t)}\sin \frac{1}{|x|}, & \text{if }
|x|\in[\frac{1}{(2k+2)\pi},\frac{1}{(2k+1)\pi}],
\end{cases}
\]
where $k\in N$ with $k\geq 1$, $1<\beta(t)<2<\alpha(t)$.
\end{remark}

\begin{proof}
Obviously, (H)(1) and (H1)(2) are satisfied.
It is also obvious that $x \mapsto F(t,x)$ is locally Lipschitz.
Then
\[
\partial F(t,x)= \begin{cases}
 0, &\text{if } |x|=0 \text{ or } |x|>\frac{1}{2\pi} ,\\
 \alpha(t)|x|^{\beta(t)-2}x\sin \frac{1}{|x|}-|x|^{\beta(t)-3}x\cos \frac{1}{|x|},
 &\text{if }|x|\in \big(\frac{1}{(2k+1)\pi},\frac{1}{2k\pi}\big),\\
\beta(t)|x|^{\alpha(t)-2}x\sin \frac{1}{|x|}-|x|^{\alpha(t)-3}x\cos \frac{1}{|x|}, 
&\text{if }|x|\in \big(\frac{1}{(2k+2)\pi},\frac{1}{(2k+1)\pi}\big),\\
[|x|^{\beta(t)-3}x,|x|^{\alpha(t)-3}x], &\text{if }|x|=\frac{1}{(2k+1)\pi},\\
[-|x|^{\beta(t)-3}x,-|x|^{\alpha(t)-3}x], &\text{if }|x|=\frac{1}{(2k+2)\pi},\\
[-|x|^{\beta(t)-3}x,0], &\text{if }|x|=\frac{1}{2\pi},
\end{cases}
\]
Hence, there exists a constant $c>0$ such that
$$
|w|\leq c(1+|x|^{\alpha(t)-1})\quad \text{for all }w\in \partial
F(t,x).
$$ 
So condition (H1)(3) holds. Then, for any $1\leq k
\in N$, we can choose
\[
a_k:=\frac{1}{(2k+2)\pi},\quad b_k:=\frac{1}{(2k+\frac{3}{2})\pi},
\]
which means $a_k< b_k$, $\lim_{k\to+\infty}b_k=0$ and 
\[
\sup \{w\cdot e_k:w\in\partial F(t,x),\text{ a.e. $t\in[0,T]$ and }
 x\in[a_k,b_k]e_k\}\leq 0.
\]
So condition (H1)(5) is satisfied.

On the other hand, for any $1\leq k \in N$, we can choose
$c_k:=\frac{1}{(2k+\frac{1}{2})\pi}$, which implies
$\lim_{k\to +\infty}c_k=0$,
\begin{gather*}
\limsup_{k\to +\infty}\frac{F(t,c_ke_k)}{|c_ke_k|^{2}}
=\limsup_{k\to +\infty}\frac{|c_ke_k|^{\beta(t)}
 \sin\frac{1}{|c_ke_k|}}{|c_ke_k|^2}
=\limsup_{k\to +\infty}\frac{1}{|c_ke_k|^{2-\beta(t)}}=+\infty,\\
-\infty<-1\leq\liminf_{|x|\to 0^+}\frac{F(t,x)}{|x|^2}
=\liminf_{|x|\to 0^+}\frac{|x|^{\alpha(t)}\sin \frac{1}{|x|}}{|x|^2}
=\liminf_{|x|\to 0^+}|x|^{\alpha(t)-2}\sin \frac{1}{|x|}\leq0
\end{gather*}
uniformly for a.e. $t\in[0,T]$. So condition (H1)(4) holds.
\end{proof}

\begin{theorem} \label{thm3.1} 
 Suppose that {\rm (H1)} holds. Then there
exists a sequence $\{u_n\}\subset E^\alpha $ of distinct positive
solution of problem \eqref{eP} such that
$$
\lim_{n\to+\infty}\|u_n\|_{\alpha}=\lim_{n\to+\infty}|u_n|_{\infty}=0.
$$
\end{theorem}

\begin{proof} For every fixed $k\in \mathbb{N}$, consider the set
$$
S_k=\{u\in E^{\alpha}: u(t)\neq 0\text{ and } u(t)\in [0, b_k]e_k\quad
\text{a.e. } t\in[0,T]\},
$$ 
where $b_k$ is from (H1)(5).
The proof is divided into four steps as follows.
\smallskip

\noindent\textbf{Step 1.} 
We  claim that $\varphi$ is bounded from below on
$S_k$ and its infimum $m_k$ on $S_k$ is attained at $u_k\in S_k$.

On account of (H1)(3) and Proposition \ref{prop2.6}, for every
 $u\in S_k$, we have
$$
F(t,x)-F(t,0)\in \langle\partial F(t,\xi),x\rangle,
$$
where $\xi=\lambda x$, and $\lambda\in(0,1)$. Furthermore, we have
\begin{equation}\label{2.3c}
|\omega|\leq c(1+|\xi|^{\alpha(t)-1})
=c(1+|\lambda|^{\alpha(t)-1}|x|^{\alpha(t)-1})
\leq c(1+|x|^{\alpha(t)-1}).
\end{equation}
Applying the Mean Value Theorem and \eqref{2.3}, for any 
$ \omega \in\partial F(t,\xi)$, we have
$$
|F(t,x)-F(t,0)|=|\langle\omega, x\rangle|
\leq|\omega|\cdot|x|\leq c(|x|+|x|^{\alpha(t)}),
$$ 
 That is,
\begin{equation}\label{2.4b}
|F(t,x)|\leq c(|x|+|x|^{\alpha(t)})
\leq c(1+|x|^{\alpha(t)}).
\end{equation}
Thus,
\begin{equation}\label{2.5}
\begin{aligned}
\varphi(u)&= \int_0^T\big[-\frac{1}{2}(\,_0^c D_t^\alpha u(t),\,_t^c
D_T^\alpha u(t))\big]dt-\int_0^TF(t,u(t))dt\\
&\geq \frac{|\cos(\pi\alpha)|}{2}\|u\|_\alpha^2-\int_0^Tc(1+|u(t)|^{\alpha(t)})dt\\
&\geq \frac{|\cos(\pi\alpha)|}{2}\|u\|_\alpha^2-\int_0^Tc(1+|u(t)|^{\alpha_0})dt\\
&\geq \frac{|\cos(\pi\alpha)|}{2}\|u\|_\alpha^2-cT-c\int_0^T|u(t)|^{\alpha_0}dt\\
&\geq \frac{|\cos(\pi\alpha)|}{2}\|u\|_\alpha^2-cT-cT|b_k|^{\alpha_0}\\
&\geq  -cT-cT|b_k|^{\alpha_0},
\end{aligned}
\end{equation}
where $\alpha_0=\inf_{t\in[0,T]}\alpha(t)$. It is clear that
$S_k$ is convex and closed, thus weakly closed in $E^\alpha$. Let
$m_k=\inf_{S_k}\varphi$, and $\{u_k^n\}_{n=1}^\infty$ be a
sequence in $S_k$ such that $m_k\leq\varphi(u_k^n)\leq
m_k+\frac{1}{n}$ for all $n\in \mathbb{N}$. Then
\begin{equation}\label{2.6}
\begin{aligned}
m_k+\frac{1}{n}
&\geq \varphi(u_k^n)\\
&= \int_o^T\big[-\frac{1}{2}(\,_0^c
D_t^\alpha u_k^n(t),\,_t^c D_T^\alpha
u_k^n(t))\big]dt-\int_0^TF(t,u_k^n(t))dt ,
\end{aligned}
\end{equation}
which implies 
\begin{equation}\label{2.7}
\begin{aligned}
\frac{|\cos(\pi\alpha)|}{2}\|u_k^n\|_\alpha^2
&\leq \int_0^T\big[-\frac{1}{2}(_0^cD_t^\alpha
u_k^n(t),_t^cD_T^\alpha u_k^n(t))\big]dt\\
&\leq m_k+\frac{1}{n}+\int_0^TF(t,u_k^n(t))dt\\
&\leq m_k+\frac{1}{n}+\int_0^Tc(1+|u_k^n(t)|^{\alpha_0})dt\\
&\leq  m_k+\frac{1}{n}+cT+cT|b_k|^{\alpha_0} ,
\end{aligned}
\end{equation} 
for all $n\in \mathbb{N}$, thus $\{u_k^n(t)\}_{n=1}^{\infty}$ is bounded in
$E^\alpha$.

By Proposition \ref{prop2.1}, one can easily see that there exists  
$\{u_k^n\}_{n=1}^\infty\in E^\alpha$ such that 
$u_k^n \rightharpoonup u_k$ in $E^\alpha$. We will show that $\varphi$ 
is weak lower semicontinuous.  Let $u_k^n\rightharpoonup u_k $ weakly in
$E^\alpha$, and by Proposition \ref{prop2.3}, we obtain the following results:
\begin{gather*}
E^\alpha \hookrightarrow L^{p}(\mathbb{R}^N),\\
u_{k}^n(t)\to u_k(t)\;{\rm a.e.}\;t\in[0,T], \\
F(t,u_{k}^n(t))\to F(t,u_k(t))\;{\rm a.e.}\;t\in[0,T].
\end{gather*}
By Fatou's lemma, 
$$
\limsup_{n\to\infty}\int_0^T F(t,u_{k}^n(t))dt\leq \int_0^T F(t,u_k(t))dt.
$$ 
On the other hand, by Proposition \ref{prop2.5}, we have 
$\lim_{n\to\infty} J(u_k^n)=J(u_k)$; that is,
$$
\lim_{n\to\infty}\int_0^T[-\frac{1}{2}(\;_0^cD_t^\alpha
u_k^n(t),\, _t^cD_T^\alpha u_k^n(t))]dt
=\int_0^T[-\frac{1}{2}(_0^cD_t^\alpha u_k(t),\,
_t^cD_T^\alpha u_k(t))]dt.
$$ 
Thus,
\begin{equation}\label{2.8}
\begin{aligned}
\liminf_{n\to\infty}\varphi(u_k^n)
&= \liminf_{n\to\infty}\int_0^T[-\frac{1}{2}(\;_0^cD_t^\alpha
u_k^n(t),\, _t^cD_T^\alpha u_k^n(t))]dt\\
&\quad -\limsup_{n\to\infty}\lambda\int_0^T
F(t,u_k^n(t))dt\\ 
&\geq \int_0^T[-\frac{1}{2}(_0^cD_t^\alpha u_k(t),\,
_t^cD_T^\alpha u_k(t))]dt-\lambda\int_0^T F(t,u_k^n(t))dt\\
&= \varphi(u_k).
\end{aligned}
\end{equation}
Then $\varphi$ is weak lower semicontinuous, and
$$
m_k\leq\varphi(u_k)\leq \lim_{\overline{n\to
+\infty}} \varphi(u_k^n)\leq m_k+\frac{1}{n},
$$ 
which implies 
$\varphi(u_k)=m_k$. Hence, $u_k$ is a minimum point of $\varphi$
over $S_k$.
\smallskip

\noindent\textbf{Step 2.} 
We show that $u_k(t)\in[0,a_k]e_k$ a.e. $t\in[0,T]$.
Let $A=\{t\in[0,T]:u_k(t)\not\in
[0,a_k]e_k\}=\{t\in[0,T]:u_k(t)\in[a_k,b_k]e_k\}$. We will prove that
$\operatorname{meas}(A)=0$. Define the function 
$h:[0,+\infty)e_k\to[0,+\infty)e_k$ by
\[
h(s)=\begin{cases}
a_ke_k, & \text{if } s\in[a_k,+\infty]e_k, \\
s, &\text{if } s\in[0,a_k]e_k.
\end{cases}
\]
Now, we set $v_k=h\circ u_k$. Since $h$ is a Lipschitz function and
$h(0)=0$, the theorem of Marcus-Mizel \cite{19} shows that $v_k\in
E^\alpha$. Moreover, $v_k(t)\in [0,a_k]e_k$ for a.e. $t\in [0,T]$.
Consequently, $v_k\in S_k$ and
\[
v_k(t)=\begin{cases}
 u_k(t), & \text{if } t\in[0,T]\backslash A, \\
 a_ke_k, & \text{if } t\in A.
\end{cases}
\]
By straightforward computations, we obtain
\begin{equation}\label{2.9}
\begin{aligned}
&\varphi(v_k)-\varphi(u_k)\\
&= \int_{[0,T]}\big[-\frac{1}{2}(_0^cD_t^\alpha
v_k(t),\;_t^cD_T^\alpha v_k(t))\big]dt
-\int_{[0,T]}F(t,v_k(t))dt\\
&\quad- \int_{[0,T]}\big[-\frac{1}{2}(_0^cD_t^\alpha
u_k(t),\;_t^cD_T^\alpha u_k(t))\big]dt+\int_{[0,T]}F(t,u_k(t))dt\\
&= \int_{[0,T]\backslash A}\big[-\frac{1}{2}(_0^cD_t^\alpha
u_k(t),\;_t^cD_T^\alpha
u_k(t))\big]dt\\
&\quad +\int_A\big[-\frac{1}{2}(_0^cD_t^\alpha a_ke_k,\;
_t^cD_T^\alpha a_ke_k)\big]dt
 -\int_{[0,T]\backslash A}F(t,u_k(t))dt\\
&\quad -\int_{A}F(t,a_ke_k)dt
-\int_{[0,T]\backslash A}\big[-\frac{1}{2}(_0^cD_t^\alpha u_k(t),\;
_t^cD_T^\alpha
u_k(t))\big]dt\\
&\quad -\int_A\big[-\frac{1}{2}(_0^cD_t^\alpha u_k(t),\;_t^cD_T^\alpha
u_k(t))\big]dt
+\int_{[0,T]\backslash A} F(t,u_k(t))\\
&\quad +\int_A F(t,u_k(t))dt\\
&= -\int_A\big[-\frac{1}{2}(_0^cD_t^\alpha u_k(t),\;_t^cD_T^\alpha
u_k(t))\big]dt -\int_A[F(t,a_ke_k)-F(t,u_k(t))]dt.
\end{aligned}
\end{equation}

For every $t\in A$, $u_k(t)\in[a_k,b_k]e_k$, there exists a map
$\lambda: A\to[0,1]$ such that
$u_k(t)=a_ke_k+\lambda(t)(b_k-a_k)e_k$.

By the Mean Value Theorem, it holds
\begin{equation}\label{2.10}
\begin{aligned}
&\int_{A}[F(t,a_ke_k)-F(t,u_k(t))]dt\\&= \int_{A} \xi_k(t)
\cdot(a_ke_k-u_k(t))dt\\
&= \int_{A} \xi_k(t)\cdot[a_ke_k-a_ke_k-\lambda(t)(b_k-a_k)e_k]  dt\\
&= \int_{A} \xi_k(t)\cdot\lambda(t)(a_k-b_k)e_k  dt,
\end{aligned}
\end{equation}
where $\xi_k(t)\in\partial F(t,\tau_k(t))$ for some 
$\tau_k(t) \in [a_ke_k, u_k(t)]\subseteq [a_k, b_k]e_k$ for a.e. $t\in A$.

  By (H1)(5), we have $\xi_k(t)\cdot e_k \leq 0$ for
a.e. $t\in A$. Consequently,
\begin{equation}\label{2.11}
\int_{A}[F(t,a_ke_k)-F(t,u_k(t))]dt\geq 0.
\end{equation}
In conclusion, every term of the expression
$\varphi(v_k)-\varphi(u_k)\leq0$. On the other hand, since 
$v_k\in S_k$, then $\varphi(v_k)\geq\varphi(u_k)=\inf_{S_k}\varphi$.
So, $\varphi(v_k)-\varphi(u_k)=0$. Namely,
\begin{equation}\label{2.12}
-\int_A\big[-\frac{1}{2}(_0^cD_t^\alpha u_k(t),_t^cD_T^\alpha
u_k(t))\big]dt-\int_A[F(t,a_ke_k)-F(t,u_k(t))]dt=0,
\end{equation}
which implies that $\operatorname{meas}(A)=0$.
\smallskip

\noindent\textbf{Step 3.} 
We  show that $u_k$ is a local minimum point in
$E^\alpha$ for every $k\in\mathbb{N}$.
Let $A'=\{t\in[0,T]:u(t)\not\in [0,a_k]e_k\}=\{t\in[0,T]:
u(t)\in(a_k,b_k]e_k\}$. Set $v=h\circ u$, then we have
\begin{equation}\label{2.13}
\begin{aligned}
&\varphi(u)-\varphi(v)\\
&= \int_{[0,T]}\big[-\frac{1}{2}(_0^cD_t^\alpha
u(t),\;_t^cD_T^\alpha u(t))\big]dt
 -\int_{[0,T]}F(t,u(t))dt\\
&\quad -\int_{[0,T]}[-\frac{1}{2}(_0^cD_t^\alpha
v(t),\;_t^cD_T^\alpha v(t))]dt+\int_{[0,T]}F(t,v(t))dt\\
&= \int_{[0,T]\backslash A'}\big[-\frac{1}{2}(_0^cD_t^\alpha
u(t),\;_t^cD_T^\alpha u(t))\big]dt \\
&\quad +\int_{A'}\big[-\frac{1}{2}(_0^cD_t^\alpha
u(t),\;_t^cD_T^\alpha u(t))\big]dt\\
&\quad -\int_{[0,T]\backslash A'}F(t,u(t))dt-\int_{A'}
F(t,u(t))dt \\
&\quad -\int_{[0,T]\backslash
A'}\big[-\frac{1}{2}(_0^cD_t^\alpha u(t),\;_t^cD_T^\alpha
u(t))\big]dt\\
&\quad -\int_{A'}\big[-\frac{1}{2}(_0^cD_t^\alpha
a_ke_k,\;_t^cD_T^\alpha a_ke_k)\big]dt \\
&\quad +\int_{[0,T]\backslash A'} F(t,u(t))+\int_{A'} F(t,a_ke_k)dt\\
&= \int_{A'}\big[-\frac{1}{2}(_0^cD_t^\alpha
u(t),\;_t^cD_T^\alpha u(t))\big]dt
 +\int_{A'}[F(t,a_ke_k)-F(t,u(t))]dt.
\end{aligned}
\end{equation}
  From assumption (H1)(5), we have
\begin{equation}\label{2.14}
\int_{A'}[F(t,a_ke_k)-F(t,u(t))]dt
=\int_{A'}\xi_k(t)\cdot(a_ke_k-u(t))dt \geq0,
\end{equation}
for a.e. $t\in A'$, where $\xi_k(t)\in\partial F(t,\tau(t))$,
$\tau(t)\in[a_ke_k,u(t)]\subseteq[a_k,b_k]e_k$, a.e. $t\in A'$.
Consequently,
\begin{equation}\label{2.15}
\varphi(u)-\varphi(v)\geq 0.
\end{equation}
On the other hand, by $v\in S_k$, we have
\begin{equation}\label{3.5}
\varphi(v)\geq\varphi(u_k).
\end{equation}
In view of \eqref{2.13}, we derive
\begin{equation}\label{2.16}
\varphi(u)-\varphi(v)\geq\int_{A'}\big[-\frac{1}{2}(_0^cD_t^\alpha
u(t),\;_t^cD_T^\alpha u(t))\big]dt .
\end{equation}
Moreover, we have
\begin{equation}\label{2.17}
\begin{aligned}
\varphi(u)
&\geq \varphi(v)+\int_{A'}\big[-\frac{1}{2}(_0^cD_t^\alpha
u(t),\;_t^cD_T^\alpha u(t))\big]dt\\
 &\geq \varphi(u_k)+\int_{A'}\big[-\frac{1}{2}(_0^cD_t^\alpha
u(t),\;_t^cD_T^\alpha u(t))\big]dt\\
 &\geq \varphi(u_k)+\int_{[0,T]}\big[-\frac{1}{2}(_0^cD_t^\alpha
u(t),\;_t^cD_T^\alpha u(t))\big]dt\\
&\quad -\int_{[0,T]\backslash
A'}\big[-\frac{1}{2}(_0^cD_t^\alpha u(t),\;_t^cD_T^\alpha
u(t))\big]dt\\
 &\geq \varphi(u_k)+\int_{[0,T]}\big[-\frac{1}{2}(_0^cD_t^\alpha
(u(t)-v(t)),\;_t^cD_T^\alpha (u(t)-v(t))\big]dt\\
 &\geq \varphi(u_k)+\frac{|\cos(\pi\alpha)|}{2}\|u-v\|_\alpha^2.
\end{aligned}\end{equation}

Since $h$ is continuous, there exists $\delta>0$ such that, for
every $u\in E^\alpha$ with $\|u-v\|_\alpha < \delta$, which
implies that $u_k$ is a local minimum of $\varphi$.
\smallskip

\noindent\textbf{Step 4.} We prove that
$m_k=\inf_{S_k}\varphi<0$ and
$\lim_{k\to+\infty}m_k=0$.
Let $B_{r_0}(t_0)\subset[0,T]$ be the  ball with radius
$r_0\in(0,1)$ and center $t_0\in[0,T]$. For $\xi\in\mathbb R^{N}$,
define
\begin{equation}\label{2.18}
\eta_\xi(t)= \begin{cases}
 0, & \text{if } t\in[0,T]\backslash B_{r_0}(t_0), \\
\xi, & \text{if } t\in B_{\frac{r_0}{2}}(t_0), \\
\frac{2\xi}{r_0}(r_0-|t-t_0|), &\text{if } 
t\in B_{r_0}(t_0)\backslash B_{\frac{r_0}{2}}(t_0).
\end{cases}
\end{equation}

It is clear that $\eta_\xi\in E^\alpha$ and
\begin{gather}\label{2.19}
|\eta_\xi(t)|\leq\frac{2|\xi|}{r_0}, \\
\label{2.20}
\begin{aligned}
|_0^cD_t^\alpha\eta_\xi(t)|
&= \big|\frac{1}{\Gamma(1-\alpha)} \int_0^t(t-s)^{-\alpha}\eta_\xi' ds \big|\\
&\leq \frac{1}{\Gamma(1-\alpha)}\Big(\int_0^t(t-s)^{-\alpha}|\eta_\xi'|ds\Big)\\
&\leq \frac{1}{\Gamma(1-\alpha)}\frac{2|\xi|}{r_0}\frac{t^{1-\alpha}}{1-\alpha}ds,
\end{aligned} \\
\label{2.21}
\begin{aligned}
\|\eta_\xi\|_\alpha^2&= \int_0^T|_0^cD_t^\alpha\eta_\xi(t)|^2 dt\\
&\leq \int_0^T\frac{1}{\Gamma^2(1-\alpha)}\frac{4|\xi|^2}{r_0^2}\frac{t^{2-2\alpha}}{(1-\alpha)^2}dt\\
&\leq \frac{1}{\Gamma^2(1-\alpha)}\frac{4\xi^2}{r_0^2}\frac{1}{(1-\alpha)^2}\int_0^T t^{2-2\alpha}dt\\
&\leq \frac{4|\xi|^2}{\Gamma^2(1-\alpha)r_0^2(1-\alpha)^2(3-2\alpha)}T^{3-2\alpha}.
\end{aligned}
\end{gather}

From the left part of (H1)(4) we deduce that the existence
of some $l_0>0$ and $\lambda_0\in[0,a_k]e_k$, such that
\begin{equation}\label{2.22}
\operatorname{ess\,inf}_{t\in[0,T]}F(t,x)\geq -l_0 |x|^2\quad
\text{for all } x\in[0,\lambda_0]e_k.
\end{equation}

There exist $L_0>0$ large enough to enable
\begin{equation}\label{2.23}
C(r_0,\alpha,T)+l_0T<\frac{1}{3}L_0r_0, \quad
C(r_0,\alpha,T)=\frac{1}{2|\cos(\pi\alpha)|}
\frac{4T^{3-2\alpha}}{\Gamma^2(1-\alpha)r_0^2(3-2\alpha)}.
\end{equation}
Taking into account the right part of (H1)(4), there is a
sequence $\{\xi_k\}\in[0,\lambda_0]$ such that
$\{\xi_k\}\in[0,a_k]e_k$ and
\begin{equation}\label{2.24}
\operatorname{ess\,sup}_{t\in[0,T]}F(t,\xi_k)>L_0|\xi_k|^2\quad
\text{for all }k\in N.
\end{equation}
Note that
$\frac{2\xi_k}{r_0}(r_0-|t-t_0|)\in[0,\xi_k]\subset[0,\lambda_0]e_k
$, for every $t\in B_{r_0}(t_0)\backslash B_{\frac{r_0}{2}}(t_0)$,
because of $|t-t_0|\in (\frac{r_0}{2},r_0)$ and $ r_0-|t-t_0|\in
(0,\frac{r_0}{2})$, $\forall t\in B_{r_0}(t_0)\backslash
B_{\frac{r_0}{2}}(t_0)$.

In view of proposition \ref{prop2.4} and \eqref{2.21}, we deduce
\begin{equation}\label{2.25}
\begin{aligned}
&\int_0^T\big[-\frac{1}{2}(_0^cD_t^\alpha
\eta_{\xi_k}(t),\;_t^cD_T^\alpha \eta_{\xi_k}(t))\big]dt\\
&\leq \frac{1}{2|\cos(\pi\alpha)|}\|\eta_{\xi_k}(t)\|_\alpha^2\\
 &\leq \frac{1}{2|\cos(\pi\alpha)|}\frac{4T^{3-2\alpha}}{\Gamma^2(1-\alpha)r_0^2(3-2\alpha)}|\xi_k|^2\\
 &= C(r_0,\alpha,T)|\xi_k|^2,
\end{aligned}
\end{equation}

And combining \eqref{2.22} with \eqref{2.24}, we obtain
\begin{equation}\label{2.26}
\begin{aligned}
&\int_0^T
F(t,\eta_\xi(t))dt\\
&=  \int_{B_{\frac{r_0}{2}}(t_0)}F(t,\eta_{\xi_k}(t))dt+\int_{B_{r_0}(t_0)\backslash
B_{\frac{r_0}{2}}(t_0)}F(t,\eta_{\xi_k}(t))dt\\
&\geq \int_{B_{\frac{r_0}{2}}(t_0)}F(t,\xi_k
(t))dt+\int_{B_{r_0}(t_0)\backslash
B_{\frac{r_0}{2}}(t_0)}F(t,\frac{2\xi_k}{r_0}(r_0-|t-t_0|))dt\\
&\geq \int_{B_{\frac{r_0}{2}}(t_0)}-l_0
|\xi_k|^2dt+\int_{B_{r_0}(t_0)\backslash
B_{\frac{r_0}{2}}(t_0)}L_0|\frac{2\xi_k}{r_0}(r_0-|t-t_0|)|^2dt\\
&=
L_0\frac{4|\xi_k|^2}{r_0^2}[\int_{t_0-r_0}^{t_0-\frac{r_0}{2}}(r_0-|t-t_0|)^2dt+\int_{t_0+\frac{r_0}{2}}^{t_0+r_0}(r_0-|t-t_0|)^2] -l_0r_0|\xi_k|^2\\
&= L_0\frac{4|\xi_k|^2}{r_0^2}[\int_{t_0-r_0}^{t_0-\frac{r_0}{2}}(r_0+t-t_0)^2dt+\int_{t_0+\frac{r_0}{2}}^{t_0+r_0}(r_0-t+t_0)^2] -l_0r_0|\xi_k|^2\\
&\geq  \frac{1}{3}L_0r_0|\xi_k|^2-l_0T|\xi_k|^2.
\end{aligned}
\end{equation}

Let $k\in \mathbb N$ be a fixed number and let $\eta_{\xi_k}\in
E^\alpha$ be the function from \eqref{2.18} corresponding to the
value $|\xi_k|>0$. Then $\eta_{\xi_k}\in S_k$, and on account of
\eqref{2.23}, \eqref{2.25} and \eqref{2.26}, one has
\begin{equation}\label{2.27}
\begin{aligned}
\varphi(\eta_{\xi_k})
&= \int_0^T\big[-\frac{1}{2}(_0^cD_t^\alpha
\eta_{\xi_k}(t),\;_t^cD_T^\alpha \eta_{\xi_k}(t))\big]dt-\int_0^T
F(t,\eta_\xi(t))dt\\
&\leq C(r_0,\alpha,T)|\xi_k|^2-\frac{1}{3}L_0r_0|\xi_k|^2+l_0T|\xi_k|^2 \\
&\leq (C(r_0,\alpha,T)+l_0T-\frac{1}{3}L_0r_0)|\xi_k|^2
<0.
\end{aligned}\end{equation}

From Step 3 and \eqref{2.27}, we deduce
\begin{equation}\label{2.28}
m_k=\varphi(u_k)=\inf_{S_k}\varphi\leq\varphi(\eta_{\xi_k})<0.
\end{equation}

Now we prove that $\lim_{k\to+\infty}m_k = 0$.
Observe that by assumption (H1)(3), one can find a
positive constant $c$ and $\omega\in\partial F(t,x)$ such that
\begin{equation}\label{2.29}
|\omega|\leq c(1+|x|^{\alpha_0}), \quad \forall t\in [0,T], x\in
\mathbb{R}^N.
\end{equation}
where $\alpha_1=\max_{t\in[0,T]}\alpha(t)$.

Applying the Mean Value Theorem and Step 1, for every $x\in[0,a_k]e_k$ and all
$t\in[0,T]$, there exists a constant $c>0$ such that
\begin{equation}\label{2.30}
|F(t,x)|=|F(t,x)-F(t,0)| \leq  c(1+|x|^{\alpha_1}).
\end{equation}
Therefore
\begin{equation}\label{2.31}
\begin{aligned}
m_k&=\varphi(u_k)\\
&=\int_0^T\big[-\frac{1}{2}(_0^cD_t^\alpha u_k(t), _t^cD_T^\alpha
u_k(t))\big]dt-\int_0^T
F(t,u_k(t))dt\\
&\geq\frac{|\cos(\pi\alpha)|}{2}\|u_k\|_\alpha^2-\int_0^T
F(t,u_k(t))dt\\
&\geq-\int_0^T F(t,u_k(t))dt\\
&\geq-\int_0^T
\big[c|u_k(t)|+c|u_k(t)|^{\alpha_1}\big]dt\\
&\geq-c T(|b_k|+|b_k|^{\alpha_1}).
\end{aligned}
\end{equation}

Since $\lim_{k\to+\infty}b_k=0$, we have
$\lim_{k\to+\infty}m_k\geq 0$. Note that $m_k<0$,
hence $\lim_{k\to+\infty}m_k=0$.

Finally, since $u_k$ are local minima of $\varphi$, they are
critical points of $\varphi$, thus weak solutions of \eqref{eP}. Due to
Step 2, there are infinitely many distinct $u_k$ with
$\lim_{k\to+\infty}|u_k|_{\infty}=0$. Moreover, we
have
\begin{equation}\label{2.32}
\begin{aligned}
\frac{|\cos(\pi\alpha)|}{2}\|u_k\|_\alpha^2&\leq \int_0^T\big[-\frac{1}{2}(_0^cD_t^\alpha
u_k(t),\;_t^cD_T^\alpha u_k(t))\big]dt\\ &= m_k+\int_0^T F(t,u_k(t))dt\\
&\leq  m_k+cT(|b_k|+|b_k|^{\alpha_1}),
\end{aligned}\end{equation}
 which means that
$\lim_{k\to+\infty}\|u_k\|_\alpha=0$.
\end{proof}

Next, we will state the counterpart of Theorem \ref{thm3.1} when the nonlinearity
 oscillates at infinity. The hypotheses
on the nonsmooth potential $F(x,t)$ are the following:

Our hypotheses on nonsmooth potential $F(x,t)$ are as follows.
\smallskip

\noindent\textbf{(H2)} $F: [0,T]\times \mathbb{R}^{N}\to\mathbb{R}$
is a function, $F(t, 0)=0$ for almost all $t\in[0,T]$ and satisfies
the following facts:
\begin{itemize}
\item[(1)] For all $x\in\mathbb{R}^{N}$, $t\mapsto F(t, x)$ is measurable;

\item[(2)] For almost all $t\in[0,T]$, $x\mapsto F(t,x)$ is locally
Lipschitz;

\item[(3)] There exist a positive constant $c$ such that for almost all
$x\in\mathbb{R}^N$, all $t\in[0,T]$ and $\omega\in\partial F(t,x)$
$$|\omega|\leq c(1+|x|^{\alpha(t)-1})$$
where $1<\alpha(t)<+\infty$;

\item[(4)] 
\[
-\infty<\liminf_{|x|\to
+\infty}\frac{F(t,x)}{|x|^2}\leq\limsup_{|x|\to
+\infty}\frac{F(t,x)}{|x|^2}=+\infty
\]
 uniformly for a.e. $x\in \mathbb{R}^{N}$;

\item[(5)] For every $k\in \mathbb N$, there exists $e_k\in\mathbb{R}^{N}$
with $|e_k|=1$ and there are two sequences $\{a_k\}$ and $\{b_k\}$
in $(0,+\infty)$ with $a_k<b_k$,
$\lim_{k\to+\infty}b_k=0$ such that
$$ 
\sup\{\omega\cdot e_k:\omega\in\partial F(t,x), \text{ a.e. }t\in[0,T], 
\;x\in[a_k,b_k]e_k\}\leq 0.
$$
\end{itemize}

\begin{remark} \label{rmk3.3} \rm
 Hypotheses (H2)(4) and (H2)(5) imply an oscillatory behaviour of $F$ 
near the infinity.
\end{remark}

\begin{remark} \label{rmk3.4} \rm
A simple example of a nonsmooth potential function satisfying
 (H2) is
\[
F(t,x)=\begin{cases}
 |x|^{\alpha(t)}\sin |x|, & \text{if } |x|\in\big[ 2k \pi, (2k+1)\pi\big), \\
 |x|^{\beta(t)}\sin |x|, &\text{if } |x|\in\big[(2k+1)\pi,(2k+2)\pi\big],
\end{cases}
\]
where $k\in N$ with $k\geq 1$, $1<\beta(t)<2<\alpha(t)<\infty$.
\end{remark}

\begin{proof}  Obviously, Hypothesis (H2)(1) and (H2)(2) 
 are satisfied.
Clearly, $x \mapsto F(t,x)$ is locally Lipschitz. Then for any
$1\leq k \in N$,
\begin{align*}
&\partial F(t,x) \\
&= \begin{cases}
\alpha(t)|x|^{\alpha(t)-2}x\sin |x|+|x|^{\alpha(t)-1}x\cos |x|, 
& \text{if }|x|\in\big(2k\pi, (2k+1)\pi\big),\\
\beta(t)|x|^{\beta(t)-2}x\sin |x|+|x|^{\beta(t)-1}x\cos |x|, 
&\text{if }|x|\in\big((2k+1)\pi,(2k+2)\pi\big),\\
[ -x|x|^{\alpha(t)-1}, -x|x|^{\beta(t)-1}\}, 
&\text{if }|x|=(2k+1)\pi,\\
[x|x|^{\alpha(t)-1},x|x|^{\beta(t)-1}\}, 
&\text{if }|x|=2k\pi,
\end{cases}
\end{align*}
where $\{\gamma,\delta\}=\{\xi:\xi=\lambda\gamma+(1-\lambda)\delta,
\lambda\in[0,1]\}$. Then, there exists a constant $c>0$ and
$\theta(t)=\alpha(t)+1$ such that
$$
|w|\leq c(1+|x|^{\theta(t)-1})\quad \text{for all } w\in \partial F(t,x).
$$ 
So condition (H2)(3) holds. Then, for any $1\leq k
\in N$, we can choose
\[
\quad a_k:=(2k+1)\pi,\quad  b_k:=(2k+\frac{3}{2})\pi,
\]
which implies $a_k< b_k$, $ \lim_{k\to+\infty}a_k=+\infty$ and
\[
\sup \{w\cdot e_k:w\in\partial F(x,t),\text{ a.e. $t\in[0,T]$  and }
x\in[a_k,b_k]e_k\}\leq 0.
\]
So condition (H2)(5) is satisfied.

On the other hand, for any $1\leq k \in N$, we can choose
$c_k:=(2k+\frac{1}{2})\pi$, which means
$\lim_{k\to +\infty}c_k=+\infty$,
\begin{gather*}
\limsup_{k\to
+\infty}\frac{F(t,c_ke_k)}{|c_k|^2}=\limsup_{k\to
+\infty}|c_k|^{\alpha(t)-2}\sin |c_k|=\limsup_{k\to
+\infty}|c_k|^{\alpha(t)-2}=+\infty,\\
-\infty<1\leq\liminf_{|x|\to
+\infty}\frac{F(t,x)}{|x|^2}=\liminf_{|x|\to
+\infty}\frac{|x|^{\beta(t)}\sin |x|}{|x|^2}=\liminf_{|x|\to
+\infty}|x|^{\beta(t)-2}\sin |x|\leq 0
\end{gather*}
uniformly for a.e. $t\in[0,T]$. So condition (H2)(4) holds.
\end{proof}

\begin{theorem} \label{thm3.2} 
 Suppose that $(H2)$ holds. Then there
exists a sequence $\{u_n\}\subset E^\alpha $ of distinct positive
solution of problem \eqref{eP} such that
$$
\lim_{n\to+\infty}\|u_n\|_{\alpha}=\lim_{n\to+\infty}|u_n|_{\infty}=+\infty.
$$
\end{theorem}

\begin{proof} For every fixed $k\in
\mathbb N$, consider the set
$$
T_k=\{u\in E^\alpha: u(x)\neq 0\;{\rm and}\; u(x)\in [0,b_k]e_k\;{\rm a.e.}\;
x\in \mathbb R^N \},
$$ 
where $b_k$ is from  (H2)(5). The
first part of the proof is similar to that of Theorem \ref{thm3.1}. Indeed,
we can prove that the functional $\varphi$ is bounded from below on
$T_k$ and its infimum on $T_k$ is attained (see Step 1 of Theorem
\ref{thm3.1}). Moreover, if $u_k\in T_k$ is chosen such that
$\varphi(u_k)=\inf_{T_k}$, then $u_k(t)\in[0,a_k]e_k$ a.e.
$t\in[0,T]$ (see Step 2 of Theorem \ref{thm3.1}), and $u_k$ is a local
minimum point of $\varphi$ in $E^\alpha$ 
(see Step 3 of Theorem \ref{thm3.1}). Instead of Step 4, we prove
\smallskip

\noindent\textbf{Step 4.} Let
$\vartheta_k=\inf_{T_k}\varphi=\varphi(u_k)$, then
$\lim_{k\to+\infty}\vartheta_k=-\infty$.
From (H2)(4), we deduce that there exist $l_{\infty}>0$
and $\lambda_{\infty}>0$ such that
\begin{equation}\label{2.33}
\operatorname{ess\,inf}_{t\in[0,T]}F(t,x)\geq -l_\infty |x|^2\quad\text{for all }
|x|>\lambda_\infty.
\end{equation}
There exist $L_\infty>0$ be large enough to enable
\begin{equation}\label{2.34}
C(r_0,\alpha,T)+l_\infty T<L_\infty r_0.
\end{equation}
From the right hand side of (H2)(4), we deduce that there
is a sequence $\{\xi_k\}\subset \mathbb R^N$ such that
$\lim_{k\to+\infty}|\xi_k|=+\infty$, and
\begin{equation}\label{2.35}
\operatorname{ess\,inf}_{t\in[0,T]}F(t,\xi_k)>L_\infty |\xi_k|^2
\quad \text{for all }k\in \mathbb N.
\end{equation}
It is easy to see that
\begin{equation}\label{2.36}
|\eta_{\xi_k}(t)|\leq|\xi_k|,\quad \forall t\in B_{r_0}(t_0)
\backslash B_{\frac{r_0}{2}}(t_0),
\end{equation}
since
\[
\eta_{\xi_k}(t)=\frac{2\xi_k}{r_0}(r_0-|t-t_0|),\;{\forall}\;t\in
B_{r_0}(t_0)\backslash B_{\frac{r_0}{2}}(t_0).
\]

Let $k\in \mathbb N$ be fixed and let $\eta_{\xi_k}\in E^\alpha$ be
the function from \eqref{2.18} corresponding to the value
$\xi_k\in\mathbb R^N$. Then $\eta_{\xi_k}\in T_{b_k}$, and on
account of \eqref{2.33} and \eqref{2.35}, we have
\begin{equation}\label{2.37}
\begin{aligned}
\varphi(\eta_{\xi_k})
&= \int_0^T\big[-\frac{1}{2}(_0^cD_t^\alpha
\eta_{\xi_k}(t),\;_t^cD_T^\alpha \eta_{\xi_k}(t))\big]dt-\int_0^T
F(t,\eta_\xi(t))dt\\
&\leq  \frac{1}{2|\cos(\pi\alpha)|}\|\eta_{\xi_k}\|_\alpha^2
 -\int_{B_{\frac{r_0}{2}}(t_0)}F(t,\eta_{\xi_k}(t))dt\\
&\quad-\int_{(B_{r_0}(t_0)\backslash
B_{\frac{r_0}{2}}(t_0))\bigcap\{t:|\eta_{\xi_k}(t)|>\lambda_\infty\}}
 F(t,\eta_{\xi_k}(t))dt\\
&\quad -\int_{(B_{r_0}(t_0)\backslash
B_{\frac{r_0}{2}}(t_0))\bigcap\{t:|\eta_{\xi_k}(t)|\leq\lambda_\infty\}}
 F(t,\eta_{\xi_k}(t))dt\\
&\leq \frac{1}{2|\cos(\pi\alpha)|}\frac{4T^(3-2\alpha)}{\Gamma^2(1-\alpha)
 r_0^2(3-2\alpha)}|\xi_k|^2-\frac{1}{3}L_\infty r_0|\xi_k|^2\\
&\quad +l_\infty T|\xi_k|^2-cT(1+\lambda_\infty^{\alpha_0})\\
&= (C(r_0,\alpha,T)-L_\infty r_0+l_\infty T)|\xi_k|^2+cT\lambda_\infty^{\alpha_1}.
\end{aligned}
\end{equation}
From \eqref{2.34}, \eqref{2.37} and
$\lim_{k\to+\infty}|\xi_k|=+\infty$, we conclude that
\begin{equation}\label{2.38}
\lim_{k\to+\infty}\varphi(\eta_{\xi_k})=-\infty.
\end{equation}
On the other hand, from
$\varphi(u_{m_k})=\min_{T_{b_{m_k}}}\varphi$, we have
\[
\varphi(u_{m_k})\leq\varphi(\eta_{\xi_k}(t)).
\]
On account of \eqref{2.38}, we have
\begin{equation}\label{2.39}
\lim_{k\to+\infty}\varphi(u_{m_k})=-\infty.
\end{equation}
Since the sequence  $\{\varphi(u_k)\}$ is non-increasing, so, we have
$$
\lim_{k\to+\infty}\vartheta_k=\lim_{k\to+\infty}\varphi(u_k)=-\infty.
$$
\smallskip

\noindent\textbf{Step 5.} We prove that
\[
\lim_{k\to+\infty}|u_{k}|_\infty=\lim_{k\to+\infty}\|u_{k}\|_{\alpha}=+\infty\,.
\]
Arguing by contradiction, assume that there exists a subsequence
$\{u_{n_k}\}$ of $\{u_{k}\}$ such that $|u_{n_k}|_\infty\leq M$ for
some $M>0$. In particular, $\{u_{n_k}\}\subset T_{b_l}$ for some
$l\in \mathbb N$. Thus, for every $n_k>l$, we have
\begin{equation}\label{2.40}
\vartheta_l\geq\vartheta_{n_k}=\inf_{T_{n_k}}\varphi=\varphi(u_{n_k})\geq
\inf_{T_{l}}\varphi=\vartheta_l.
\end{equation}
So, $\vartheta_{n_k}=\vartheta_l$ for every $n_k>l$. This fact
contradicts with \eqref{2.39} which completes the first part of the
proof.

Next, we prove that $\lim_{k\to+\infty}\|u_{k}\|_\alpha=+\infty$.
Note that $1 <\alpha_1<+\infty$, then by Proposition \ref{prop2.3}, we have
$E^\alpha\hookrightarrow C([0,T],\mathbb R^N)$ (compact embedding).
Furthermore, there exists $c_1>0$ such that 
$|u_k|_\infty\leq c_1\|u_k\|_{\alpha}$. Hence, there exists a constant $c_2>0$, 
such that
\begin{equation}\label{2.41}
\begin{aligned}
\int_0^T F(t,u_k(t))dt
&\leq \int_0^T c(1+|u_k(t)|^{\alpha_1})dt \\
&\leq c T+c|u_k(t)|_\infty^{\alpha_1}T\\
&\leq c T+c c_1^{\alpha_1}\|u_k\|_{\alpha}^{\alpha_1}T\\
&\leq cT+c_2\|u_k\|_{\alpha}^{\alpha_1}T.
\end{aligned}
\end{equation}

Let us assume that there exists a subsequence $\{u_{n_k}\}$ of
$\{u_{k}\}$ such that for some $M>0$, we have
$\|u_{n_k}\|_\alpha\leq M$. In particular, by the above
inequality,
\begin{equation}\label{2.42}
\begin{aligned}
|\varphi(u_{n_k})|
&=\Big|\int_0^T\big[-\frac{1}{2}(_0^cD_t^\alpha
u_{n_k}(t),\;_t^cD_T^\alpha u_{n_k}(t))\big]dt-\int_0^T
F(t,u_{n_k}(t))dt\Big|\\
&\leq \int_0^T\big[-\frac{1}{2}(_0^cD_t^\alpha
u_{n_k}(t),\;_t^cD_T^\alpha u_{n_k}(t))\big]dt +\big|\int_0^T
F(t,u_{n_k}(t))dt\big|\\
&\leq \frac{1}{2|\cos(\pi\alpha)|}\|u_{n_k}\|_\alpha^2+cT+c_2\|u_k\|_{\alpha}T
\end{aligned}
\end{equation}
is bounded. Hence $\vartheta_{n_k}=\varphi(u_{n_k})$ is also
bounded. This fact contradicts with
$\lim_{k\to+\infty}\vartheta_k=-\infty$.
\end{proof}

\subsection*{Acknowledgments}
This research was supported by the NNSF of China (No. 11201095),
the Youth Scholar Backbone Supporting Plan Project of Harbin Engineering University,
the Fundamental Research Funds for the Central Universities(No. 2016),
Postdoctoral research startup foundation of Heilongjiang (No. LBH-Q14044),
the Science Research Funds for Overseas Returned Chinese Scholars
of Heilongjiang Province (No. LC201502).



\begin{thebibliography}{10}

\bibitem{13} A. Ambrosetti, P. H. Rabinowitz;
\emph{Dual variational methods in  critical point theory and applications}. 
J. Func. Anal., 14 (1973), 349-381.

\bibitem{31} D. Anderson, R. Avery;
\emph{Fractional-order boundary value problem with
Sturm-Liouville boundary conditions}. Electron. J. Differential
Equations, 29 (2015) 1-10.

\bibitem{32} I. Bachar, H. Maagli, V. Radulescu;
\emph{Fractional Navier boundary value problems}. 
Bound. Value Probl., 79 (2016), 1-14.

\bibitem{1} F. Clarke;
\emph{Optimization and Nonsmooth Analysis}. John Wiley and Sons, New York, 1983.

\bibitem{14} K. Chang;
\emph{Variational methods for non-differentiable functionals and
 their applications to partial differential equations}. J. Math.
 Anal. Appl., 80 (1981), 102-129.

\bibitem{16} Y. K. Chang, J. J. Nieto;
\emph{Some new existence results for fractional
differential inclusions with boundary conditions}. Math. Comput.
Model., 49 (2009), 605-609.

\bibitem{20} J. Chen, X. H. Tang;
\emph{Existence and multiplicity of solutions for some fractional
boundary value problem via critical point theory}. Abstr. Appl. Anal., ID: 648635,
2012: 1-21.

\bibitem{23} J. Chen, X. H. Tang;
\emph{Infinitely many solutions for a class of fractional boundary value
problem}. Bull. Malays. Math. Sci. Soc., 36 (2013), 1083-1097.

\bibitem{33} L. Del Pezzo, J. Rossi, N. Saintier, A. Salort;
\emph{An optimal mass transport approach for limits of eigenvalue problems
 for the fractional $p$-Laplacian}. Adv. Nonlinear Anal., 4 (2015), 235-249.

\bibitem{34} S. Goyal, K. Sreenadh;
\emph{Existence of multiple solutions of
$p$-fractional Laplace operator with sign-changing weight function}.
Adv. Nonlinear Anal., 4 (2015),  37-58.

\bibitem{19} F. Jiao, Y. Zhou;
\emph{Existence of solutions for a class of fractional
boundary value problem via critical point theory}. Comput. Math.
Appl., 62 (2011), 1181-1191.

\bibitem{22} F. Jiao, Y. Zhou;
\emph{Existence results for fractional boundary value problem via critical point
theory}. J. Bifur. Chaos. App. Sci. Engrg., 22 (2012), 1250086.


\bibitem{4} F. Mainardi;
\emph{Fractional calculus: Some basic problems in
  Continuum and Statistical mechanics, Fractals and Fractional
  Calculus in Continuum Mechanics (Udine, 1996)}, 291-348, CISM Course and Lectures,
  378, Springer, Vienna, 1997.

\bibitem{35} G. Molica Bisci, V. Radulescu, R. Servadei;
\emph{Variational methods for nonlocal fractional problems}, 
Encyclopedia of Mathematics and its
Applications, 162. Cambridge University Press, Cambridge, 2016.

\bibitem{36} G. Molica Bisci, D. Repovs;
\emph{On doubly nonlocal fractional elliptic equations}. 
Atti Accad. Naz. Lincei Rend. Lincei Mat. Appl., 26 (2015), 161-176.

\bibitem{11} J. Mawhin, M. Willem;
\emph{Critical point theory and hamiltonian  systems}. Spinger, New York, 1989.

\bibitem{15} A. Ouahab;
\emph{Some results for fractional boundary value problem of
 differential inclusions}. Nonlinear Anal., 69 (2008), 3877-3896.

\bibitem{12} P. H. Rabinowitz;
\emph{Minimax methods in critical point theory with
applications to differential equations}, CBMS Reg. Conf. Ser. Math.
65, Amer. Math. Soc., Providence, RI, 1986.

\bibitem{21} H. R. Sun, Q. G. Zhang;
\emph{Existence of solutions for a fractional boundary value problem via
the Mountain Pass method and an iterative technique}. Computers Math. Appl., 64 (2012),
3436-3443.

\bibitem{18} K. M. Teng, H. G. Jia, H. F. Zhang;
\emph{Existence and multiplicity results  for fractional differential 
inclusions with Dirichlet boundary  conditions}. 
Appl. Math. Comput., 220 (2013), 792-801.

\bibitem{10} Z. L. Wei, W. Dong, J. L. Che;
\emph{Periodic boundary value problems for
 fractional differential equation involving a Riemann-Liouville fractrional 
derivative}. Nonlinear  Anal., 73 (2010), 3232-3238.

\bibitem{17} J. R. Wang, Y. Zhou;
\emph{Existence and controllability results for
fractional semilinear differential inclusions}. Nonlinear Anal. Real
World Appl., 12 (2012), 3642-3653.

\end{thebibliography}

\end{document}

