\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2013 (2013), No. 136, pp. 1--12.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu
\newline ftp ejde.math.txstate.edu}
\thanks{\copyright 2013 Texas State University - San Marcos.}
\vspace{9mm}}

\begin{document}
\title[\hfilneg EJDE-2013/136\hfil Infinitely many solutions]
{Infinitely many solutions for a  perturbed nonlinear fractional
boundary-value problem}

\author[C. Bai \hfil EJDE-2013/136\hfilneg]
{Chuanzhi Bai}  % in alphabetical order

\address{Chuanzhi Bai \newline
Department of Mathematics,
Huaiyin Normal University\\
Huaian, Jiangsu 223300, China}
\email{czbai8@sohu.com}


\thanks{Submitted January 18, 2013. Published June 20, 2013.}
\subjclass[2000]{58E05, 34B15, 26A33}
\keywords{Infinitely many solutions; fractional boundary value problem;
\hfill\break\indent critical point}

\begin{abstract}
 Using variational methods, we prove the existence of
 infinitely many  solutions for a class of nonlinear fractional
 boundary-value problems depending on two parameters.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{example}[theorem]{Example}
\allowdisplaybreaks

\section{Introduction}

In recent years, some fixed point theorems and monotone iterative
methods have been  applied successfully to investigate the
existence of solutions for nonlinear fractional boundary-value
problems, see for example,
\cite{ABH,AN,B1,B2,B4,BL,BHN,GK,K,LZL,WZWH,WZ,Z} and the references
therein. But till now, there are few results on the solutions to
fractional boundary value problems which  are established by the
variational methods. It is often very difficult to establish a
suitable space and variational functional for fractional boundary
value problem for several reasons. First and foremost, the
composition rule  in general fails to be satisfied by fractional
integral and fractional derivative operators. Furthermore, the
fractional integral is a singular integral operator and fractional
derivative operator is non-local. Besides, the adjoint of a
fractional differential operator is not the negative of itself.
Recently, by using critical point theory, Jiao and Zhou \cite{JZ}
studied the  fractional BVP
\begin{gather*}
 _t D_T^{\alpha}(_0 D_t^{\alpha} u(t))  =  \nabla F(t, u(t)), \quad
 \text{a.e. } t \in [0, T], \\
   u(0) =0, \quad  u(T) = 0,
  \end{gather*}
 where $_0 D_t^{\alpha}$ and $_t D_T^{\alpha}$ are the left and
right Riemann-Liouville fractional derivatives of order
$0 < \alpha \le 1$ respectively,  $F : [0, T] \times \mathbb{R}^N \to
\mathbb{R}$ is a given function satisfying some assumptions and
$\nabla F(t, x)$  is the gradient of $F$ at $x$.


In \cite{B3},  by using a local minimum theorem established by
 Bonanno in \cite{B5}, we  provided  a new approach to investigated the
existence of solutions for the following  fractional boundary-value problem
\begin{gather*}
\frac{d}{dt} \Big({}_0 D_t^{\alpha-1}({}_0^c D_t^{\alpha} u(t)) -
{}_t D_T^{\alpha-1}({}_t^c D_T^{\alpha} u(t))\Big)
+ \lambda f(u(t)) = 0, \quad \text{a.e. }  t \in [0, T], \\
   u(0) =0,\quad  u(T) = 0,
\end{gather*}
  where $\alpha \in (1/2, 1]$, ${}_0 D_t^{\alpha-1}$ and
${}_t D_T^{\alpha-1}$ are the left and
 right  Riemann-Liouville fractional integrals of order $1-\alpha$ respectively,
$_0^c D_t^{\alpha}$ and $_t^c D_T^{\alpha}$ are the left and right
Caputo  fractional derivatives of order $0< \alpha \leq 1$
respectively, $\lambda$ is a positive real parameter,  and
$f :\mathbb{R} \to \mathbb{R}$ is a continuous function.

The purpose of this article is to establish the existence of
infinitely many  solutions for the following perturbed nonlinear
fractional boundary value problem
\begin{equation}\label{eq:1.1}
\begin{gathered}
 _t D_T^{\alpha}(_0 D_t^{\alpha} u(t)) = \lambda a(t) f(u(t))
+ \mu g(t, u(t)), \quad\text{a.e. } t \in [0, T], \\
    u(0) = 0,\quad  u(T) = 0,
  \end{gathered}
\end{equation}
  where $_0 D_t^{\alpha}$ and $_t D_T^{\alpha}$ are the left and
right Riemann-Liouville fractional derivatives of order $0 <
\alpha \le 1$ respectively, $\lambda$ and $\mu$ are non-negative
parameters, $a : [0,  T] \to \mathbb{R}$, $f : \mathbb{R} \to
\mathbb{R}$ and  $g : [0,  T] \times \mathbb{R} \to \mathbb{R}$
are three  given continuous functions.


Precisely, under appropriate hypotheses on the nonlinear term $f,
g$, the existence of two intervals $\Lambda$ and $J$ such that, for
each $\lambda \in \Lambda$ and $\mu \in J$,  BVP \eqref{eq:1.1}
admits a sequence of pairwise distinct solutions is proved.   Our
analysis is mainly based on a recent critical point theorem of
Bonanno and Molica Bisci \cite{BM1}, which is a more precise version
of Ricceri's Variational Principle \cite{R}. This theorem  and its
variations have been used in several works in order to obtain
existence of infinitely many solutions for different kinds of
problems (see, for instance,
\cite{AH,BM1,BB,BM2,BMD,BMO,BT} and references
therein). The technique used in this paper in order to approach
perturbed nonlinearity depending on two parameters has been adopted
first in the paper \cite{BC}. Moreover, among authors that follow
this technique, we recall the papers
\cite{BMW,CD,DS,DM,HT}.


This article is organized as follows. In section 2, we present some
necessary preliminary facts that will be needed in the paper. In
section 3,  we establish our main two existence results and give
an example to show the effectiveness of the our results.


\section{Preliminaries}


 In this section, we first introduce some necessary definitions and properties
of the fractional calculus which are used in this paper.

\begin{definition}\label{d2.1}\rm
  Let $f$ be a function defined on
 $[a , b]$. The left and right Riemann-Liouville fractional
 integrals  of order $\alpha$ for function $f$ denoted by
 $_a D_t^{-\alpha}f(t)$ and   $_t D_b^{-\alpha}f(t)$,  respectively,
are defined by
\begin{gather*}
_a D_t^{-\alpha}f(t) = \frac{1}{\Gamma(\alpha)} \int_a^t
(t-s)^{\alpha-1} f(s) ds, \quad t \in [a, b], \; \alpha > 0, \\
 _t D_b^{-\alpha}f(t) = \frac{1}{\Gamma(\alpha)} \int_t^b
(s-t)^{\alpha-1} f(s) ds, \quad t \in [a, b], \; \alpha > 0,
\end{gather*}
 provided the right-hand sides are pointwise defined on $[a, b]$,
where $\Gamma(\alpha)$ is the gamma function.
\end{definition}


\begin{definition}\label{d2.2}\rm
  Let $f$ be a function defined on
 $[a , b]$.  For $n-1 \le \gamma < n$($n \in \mathbb{N}$), the left and right Riemann-Liouville fractional
  derivatives  of order $\gamma$ for function $f$ denoted by  $_a D_t^{\gamma}f(t)$ and
  $_t D_b^{\gamma}f(t)$,  respectively, are defined by
$$  _a D_t^{\gamma}f(t) = \frac{d^n}{d t^n}
{\empty}_a D_t^{\gamma-n}f(t) =  \frac{1}{\Gamma(n-\gamma)}
\frac{d^n}{d t^n} \int_a^t (t-s)^{n-\gamma-1} f(s) ds, \quad t \in
[a, b]
$$
 and
$$ _t D_b^{\gamma}f(t) = (-1)^n \frac{d^n}{d t^n} {\empty}_t
D_b^{\gamma-n}f(t) = \frac{(-1)^n}{\Gamma(n-\gamma)} \frac{d^n}{d
t^n} \int_t^b (s-t)^{n-\gamma-1} f(s) ds, \quad t \in [a, b].
$$
\end{definition}


According to \cite{KST,P}, if $f \in AC^n([a, b],
\mathbb{R}^N)$, then by performing repeatedly integration by parts
and differentiation, for $n-1 \le \gamma < n$, we have
\begin{equation}\label{eq:2.1}
 _a D_t^{\gamma}f(t) = \frac{1}{\Gamma(n-\gamma)}  \int_a^t
\frac{f^{(n)}(s)}{(t-s)^{\gamma+1-n}} ds + \sum
_{j=0}^{n-1} \frac{f^j(a)}{\Gamma(j-\gamma+1)}
(t-a)^{j-\gamma},
\end{equation}
 and
\begin{equation}\label{eq:2.2}
_t D_b^{\gamma}f(t) = \frac{(-1)^n}{\Gamma(n-\gamma)}  \int_t^b
\frac{f^{(n)}(s)}{(s-t)^{\gamma+1-n}} ds + \sum
_{j=0}^{n-1} \frac{(-1)^j f^j(b)}{\Gamma(j-\gamma+1)}
(b-t)^{j-\gamma},
\end{equation}
where $t \in [a, b]$.

From \cite{KST}, \cite{SKM}, we have the following property of
fractional integration.

\begin{proposition}\label{pro:2.1}
  If  $f \in L^p([a, b],
\mathbb{R}^N)$, $g \in L^q([a, b], \mathbb{R}^N)$ and $p \ge 1$,
$q \ge 1$, $1/p + 1/q \le 1 + \gamma$ or $p \not= 1$, $q \not= 1$,
$1/p + 1/q = 1 + \gamma$, then
\begin{equation}\label{eq:2.3}
 \int_a^b [_a D_t^{-\gamma}f(t)] g(t) dt
= \int_a^b[_t D_b^{-\gamma}g(t)] f(t) dt, \quad \gamma > 0.
\end{equation}
\end{proposition}

By properties \eqref{eq:2.1}--\eqref{eq:2.3}, one has (see \cite{JZ})

\begin{proposition}\label{pro:2.2}
 If $f(a) = f(b) = 0$, $ f' \in L^{\infty}([a, b], \mathbb{R}^N)$ and $g \in
L^1([a, b], \mathbb{R}^N)$, or $g(a) = g(b) = 0$, $ g' \in
L^{\infty}([a, b], \mathbb{R}^N)$ and $f \in L^1([a, b],
\mathbb{R}^N)$, then
\begin{equation}\label{eq:2.4}
 \int_a^b \left[_a D_t^{\alpha}f(t)\right] g(t) dt = \int_a^b
\left[_t D_b^{\alpha}g(t)\right] f(t) dt, \quad  0 < \alpha \le 1.
\end{equation}
\end{proposition}

To establish a variational structure for BVP
\eqref{eq:1.1}, it is necessary to construct appropriate function
spaces.

 \begin{definition}\label{d2.3}\rm
 Let $0 < \alpha \le 1$. The fractional derivative space
$E_0^{\alpha}$  is defined
 by the closure of $C_0^{\infty}([0, T], \mathbb{R})$ with respect
 to the norm
$$
\|u\|_{\alpha} = \Big(\int_0^T |_0 D_t^{\alpha} u(t)|^2 dt +
\int_0^T |u(t)|^2 dt   \Big)^{1/2},
   \quad \forall u \in E_0^{\alpha}.
$$
\end{definition}

It is obvious that the fractional derivative space $E_0^{\alpha}$
is the space of functions $u \in L^2[0, T]$ having an
$\alpha$-order  fractional derivative $_0 D_t^{\alpha} u \in
L^2[0, T]$   and $u(0) = u(T) = 0$.

\begin{proposition}[\cite{JZ}] \label{pro:2.3}
 Let $0  < \alpha \le 1$. The fractional derivative space $E_0^{\alpha}$ is
reflexive and separable Banach space.
\end{proposition}

\begin{proposition}[\cite{JZ}] \label{pro:2.4}
   Let $ 1/2 < \alpha \le 1$. For all $u \in E_0^{\alpha}$, we have
\begin{itemize}
\item[(i)]
\begin{equation}\label{eq:2.5}
\|u\|_{L^2} \le \frac{T^{\alpha}}{\Gamma(\alpha+1)} \|\empty_0
D_t^{\alpha}u\|_{L^2}.
\end{equation}


\item[(ii)]  \begin{equation}\label{eq:2.6}
  \|u\|_{\infty} = \max _{t \in [0, T]}
|u(t)| \le  \frac{T^{\alpha -
1/2}}{\Gamma(\alpha)(2(\alpha-1)+1)^{1/2}} \|_0 D_t^{\alpha}
u\|_{L^2}.  \end{equation}
\end{itemize}
\end{proposition}

 By \eqref{eq:2.5}, we can consider $E_0^{\alpha}$ with respect to the norm
 \begin{equation}\label{eq:2.7}
 \|u\|_{\alpha} = \Big(\int_0^T |_0 D_t^{\alpha} u(t)|^2
dt\Big)^{1/2}  = \|_0 D_t^{\alpha} u\|_{L^2},    \quad \forall u
\in E_0^{\alpha}
 \end{equation}
 in the following analysis.


We are now in a position to give the definition for the solution
of BVP \eqref{eq:1.1}.

\begin{definition}\label{d2.4}\rm
    A function $ u : [0, T] \to
\mathbb{R}$ is called a solution of BVP \eqref{eq:1.1} if
\begin{itemize}
\item[(i)]  $_t D_T^{\alpha-1} (_0 D_t^{\alpha} u(t))$ and $_0
D_t^{\alpha-1} u(t)$  exist for almost every $t \in [0, T]$, and

\item[(ii)]  $u$ satisfies \eqref{eq:1.1}.
\end{itemize}
\end{definition}


By using  \eqref{eq:2.4} and Definition \ref{d2.4},  we can give
the definition of weak solution for BVP \eqref{eq:1.1} as follows.

\begin{definition}\label{d2.5}\rm
 By the weak solution of BVP \eqref{eq:1.1}, we mean that the function
  $u \in E_0^{\alpha}$ such that
$a(\cdot) f(u(\cdot)) \in L^1[0, T]$, $g(\cdot, u(\cdot)) \in
L^1[0, T]$  and satisfies
 $$  \int_0^T [_0 D_t^{\alpha} u(t) \cdot {\empty}_0 D_t^{\alpha}
v(t) - \lambda a(t) f(u(t)) \cdot v(t) - \mu g(t, u(t))\cdot
v(t)]dt = 0  $$
 for all $v \in C_0^{\infty}([0, T], \mathbb{R})$.
\end{definition}

By \cite[Theorem 5.1]{JZ}, we have

\begin{theorem}\label{thm:2.1}
 Let $0 < \alpha \le 1$ and $u \in E_0^{\alpha}$.
If $u$ is a non-trivial weak solution of 
 \eqref{eq:1.1}, then $ u$ is also a non-trivial solution
 of  \eqref{eq:1.1}.
\end{theorem}


Our main tools is an infinitely many critical points theorem
\cite{BM1} which is recalled below.

\begin{theorem}\label{thm:2.2}
  Let $X$ be a reflexive real
Banach space; $\Phi, \Psi  : X \to \mathbb{R}$ be two
 G\^{a}teaux   differentiable functionals such that $\Phi$ is  sequentially weakly
 lower semicontinuous,  strongly continuous, and coercive and $\Psi$ is
 sequentially weakly upper
 semicontinuous. For every $r > \inf_X \Phi$, let us put
$$
\varphi(r) = \inf _{u \in \Phi^{-1}(]-\infty, r[)}
\frac{\sup _{u \in \Phi^{-1}(]-\infty, r[)} \Psi(v) -
\Psi(u)}{r - \Phi(u)}
 $$
 and
 $$ \gamma = \liminf _{r \to + \infty} \varphi(r), \quad
 \delta = \liminf _{r \to (\inf_X \Phi)^+} \varphi(r).
$$

(1)   If $\gamma < + \infty$ then, for each $\lambda \in ]0, \frac{1}{\gamma}[$,
 the following alternative holds: either the functional $\Phi - \lambda \Psi$
 has a global minimum,
 or there exists a sequence $\{u_n\}$ of critical points (local minima)
 of $\Phi - \lambda \Psi$ such that
 $\lim_{n \to + \infty} \Phi(u_n) = + \infty$.

(2)  If $\delta < + \infty$ then, for each $\lambda \in ]0,
\frac{1}{\delta}[$,  the following alternative holds: either there
exists a global minimum of $\Phi$ which is a local minimum of
 $\Phi - \lambda \Psi$, or there exists a sequence $\{u_n\}$ of pairwise
 distinct critical points (local minima)
 of $\Phi - \lambda \Psi$, with  $\lim_{n \to + \infty} \Phi(u_n) = \inf_X \Phi$,
 which weakly converges  to a global minimum of $\Phi$.

\end{theorem}


\section{Main results}

We define the functional $\Phi, \Psi : E_0^{\alpha} \to R$ by
setting, for every $u \in E_0^{\alpha}$,
\begin{gather}\label{eq:3.1}
 \Phi(u) := \frac{1}{2} \|u\|_{\alpha}^2,
\\
\label{eq:3.2}
 \Psi(u) := \int_0^T \left[a(t)
F(u(t)) + \frac{\mu}{\lambda} G(t, u(t))\right] dt,
\end{gather}
 where $F(u) = \int_0^u f(s)ds$ and $G(t, u) = \int_0^u g(t, x)
dx$.  Clearly, $\Phi$ and $\Psi$ are G\^{a}teaux  differentiable
functional whose G\^{a}teaux  derivative at the point
$u \in E_0^{\alpha}$ are given by
\begin{gather*}
\Phi'(u)v =  \int_0^T
{\empty}_0 D_t^{\alpha} u(t) \cdot {\empty}_0 D_t^{\alpha} v(t) dt,
\\
\Psi'(u)v =  \int_0^T
\left(a(t) f(u(t)) + \frac{\mu}{\lambda} g(t, u(t))\right) v(t) dt
\end{gather*}
 for every $v \in E_0^{\alpha}$.   Hence, a critical point of
$I_{\lambda} = \Phi - \lambda \Psi$,   gives us a weak solution of
\eqref{eq:1.1}, which is also a solution of \eqref{eq:1.1}
by Theorem \ref{thm:2.1}.


  If $\alpha > 1/2$, by Proposition \ref{d2.4} and \eqref{eq:2.7},
 one has
\begin{equation}\label{eq:3.3}
  \|u\|_{\infty} \le  M \Big(\int_0^T |_0 D_t^{\alpha} u(t)|^2 dt\Big)^{1/2}
  = M \|u\|_{\alpha}, \quad  u \in E_0^{\alpha},
\end{equation}
 where
\begin{equation}\label{eq:3.4}
M =  \frac{T^{\alpha-\frac{1}{2}}}{\Gamma(\alpha) \sqrt{2(\alpha-1)
+ 1}}.
\end{equation}

Given a constant $0 < h  < 1/2$, put
\begin{align*}
A(\alpha, h)
&:= \frac{1}{ 2 h^2 T^2} \Big[ \frac{1 + h^{3-2\alpha}
+ (1-h)^{3-2\alpha}}{3-2\alpha} T^{3-2\alpha}  \\
& \quad  - 2 \int_{(1-h)T}^T t^{1-\alpha} (t-(1-h)T)^{1-\alpha} dt 
  -  2 \int_{hT}^T t^{1-\alpha} (t-hT)^{1-\alpha} dt \\
&  \quad  + 2 \int_{(1-h)T}^T (t - hT)^{1-\alpha}
(t-(1-h)T)^{1-\alpha} dt \Big],
\end{align*}
\begin{equation}\label{eq:3.5}
 K := \frac{\Gamma^2(2-\alpha)
\int_{hT}^{(1-h)T} a(t) dt}{2 M^2 A(\alpha, h) \int_0^T a(t) dt},
\end{equation}
\begin{equation}\label{eq:3.6}
\lambda_1 = \begin{cases}
\frac{A(\alpha, h)}{\Gamma^2(2-\alpha)
 \int_{hT}^{(1-h)T} a(t) dt \cdot \limsup_{\xi \to + \infty} \frac{F(\xi)}{\xi^2}}, &
\text{if }  \limsup _{\xi \to + \infty} \frac{F(\xi)}{\xi^2} < + \infty, \\
0, &  \text{if }  \lim\sup _{\xi \to + \infty}
\frac{F(\xi)}{\xi^2} = + \infty
\end{cases}
\end{equation}
\begin{equation}\label{eq:3.7}
\lambda_2 = \frac{1}{2 M^2 \int_0^T a(t) dt \cdot \lim\inf
_{\xi \to + \infty} \frac{F(\xi)}{\xi^2} },
\end{equation}
where $M$ as in \eqref{eq:3.4}.

\begin{theorem}\label{thm:3.1}
  Let $1/2 < \alpha \le 1$, $0 < h < 1/2$,
 $a : [0, T] \to \mathbb{R}$ and $f : \mathbb{R} \to \mathbb{R}$ be
 two nonnegative continuous functions, and assume that
% \item[(H1)]
\begin{equation}\label{eq:3.8}
   0 <  \liminf _{\xi \to
+ \infty} \frac{F(\xi)}{\xi^2} < K  \limsup _{\xi \to +
\infty} \frac{F(\xi)}{\xi^2},
\end{equation}
  where $K$ is given in \eqref{eq:3.5}.
For every $\lambda \in \Lambda
:=  ]\lambda_1, \lambda_2[$ ($\lambda_1$ and $\lambda_2$ are given
in \eqref{eq:3.6} and \eqref{eq:3.7} respectively)
and for every $g \in C([0, T] \times \mathbb{R})$ such that
\begin{equation} \label{H2}
\begin{gathered}
 G(t, u) \ge 0, \quad \forall (t,u) \in [0, T]
\times [0, + \infty[,  \\
 G_{\infty} := \limsup_{\xi \to
+\infty} \frac{\int_0^T \max_{|x| \le \xi} G(t, x) dt}{\xi^2} < +
\infty,
\end{gathered}
\end{equation}
if we put
$$
\mu_* := \frac{1}{2 M^2  G_{\infty}}\Big(1 - 2 M^2\lambda \int_0^T a(t) dt
\cdot \lim\inf _{\xi \to + \infty}
 \frac{F(\xi)}{\xi^2}\Big),
$$
 with $\mu_* = + \infty$ when $G_{\infty} = 0$,
then \eqref{eq:1.1}  possesses an unbounded sequence of solutions in
$E_0^{\alpha}$ for every $\mu \in J:= [0, \mu_*[$.
\end{theorem}

\begin{proof}
Our aim is to apply part (1) of  Theorem \ref{thm:2.2}.
  Let $\Phi, \Psi$ be the functionals defined in \eqref{eq:3.1} and
\eqref{eq:3.2} respectively. By the Lemma 5.1 in \cite{JZ}, $\Phi$
is continuous and convex, so it is weakly sequentially lower
semicontinuous, moreover, $\Phi$ is continuously G\^{a}teaux
differentiable and coercive.  The functional $\Psi$ is well defined,
continuously G\^{a}teaux  differentiable and with compact
derivative, hence it is sequentially weakly upper semicontinuous. It
is well known that the critical point of the functional $\Phi -
\lambda \Psi$ in $E_0^{\alpha}$ is exactly the solution of
\eqref{eq:1.1}.


  Let $\rho_n$ be a sequence of positive numbers  such that
   $\lim_{n \to \infty} \rho_n = + \infty$ and
$$
\lim _{n \to \infty} \frac{F(\rho_n)}{\rho_n^2}
 = \liminf _{\xi \to + \infty} \frac{F(\xi)}{\xi^2}.
$$
Let $r_n = \frac{\rho_n^2}{2 M^2}$ for all $n \in \mathbb{N}$. By
\eqref{eq:3.3}, for all $v \in E_0^{\alpha}$ such that $\Phi(v)
\le r_n$, one has  $\|v\|_{\infty} \le \rho_n$. Thus,
\begin{align*}
 \varphi(r_n) & = \inf _{u \in \Phi^{-1}(]-\infty, r_n[)}
\frac{\sup _{v \in \Phi^{-1}(]-\infty, r_n[)} \Psi(v) -
\Psi(u)}{r_n - \Phi(u)} \\
 & \le \frac{\sup _{v \in \Phi^{-1}(]-\infty, r_n[)}
 \Psi(v)}{r_n} \\
 & \le  2 M^2 \int_0^T a(t) dt \cdot  \frac{F(\rho_n)}{\rho_n^2}
 + \frac{2 M^2 \mu}{\lambda}
\frac{\int_0^T \max_{|\xi| \le \rho_n}G(t, \xi) dt}{\rho_n^2}.
\end{align*}
 So,
 $$
\gamma \le \lim \inf _{n \to + \infty} \varphi(r_n)
\le 2 M^2 \int_0^T a(t) dt \cdot
 \lim \inf _{\xi \to + \infty} \frac{F(\xi)}{\xi^2} + \frac{2 M^2
 \mu}{\lambda} G_{\infty}  < + \infty.
$$
  Thus, it is easy to verify that when $G_{\infty} > 0$,
for every $\mu \in J$,
  $$
\gamma <  2 M^2 \int_0^T a(t) dt \cdot  \liminf _{\xi \to
+ \infty} \frac{F(\xi)}{\xi^2} + \frac{2 M^2  \mu_*}{\lambda}
G_{\infty} = \frac{1}{\lambda},
$$
 while, when $G_{\infty} = 0$, we
have by $\lambda \in ]\lambda_1, \lambda_2[$ that,
$$
\gamma <  2 M^2 \int_0^T a(t) dt
\cdot  \liminf _{\xi \to + \infty} \frac{F(\xi)}{\xi^2} <
\frac{1}{\lambda}.
$$
 Thus, we conclude that
$$ \Lambda \subset ]0,
\frac{1}{\gamma}[,
$$
by the definition of $\Lambda$.   Now, we claim that the
functional $\Phi - \lambda \Psi$ is unbounded from below. Let
$\{\eta_n\}$ be a positive real sequence such that
 $\lim_{n \to \infty} \eta_n = + \infty$. We consider a function
 $v_n$ defined by setting
\begin{equation}\label{eq:3.9}
 v_n(t) = \begin{cases}
\frac{\Gamma(2-\alpha) \eta_n}{hT} t, &  t \in [0, h T[,\\
 \Gamma(2-\alpha) \eta_n, &  t \in [hT, (1-h)T],\\
 \frac{\Gamma(2-\alpha) \eta_n}{hT}(T - t), &  t \in ](1-h)T, T],
\end{cases}
\end{equation}
where $0 < h < 1/2$.  Clearly  $v_n(0) = v_n(T) = 0$ and
$v_n \in L^2[0, T]$. A direct calculation shows that
\begin{equation*}
 {\empty}_0 D_t^{\alpha} v_n(t)
= \begin{cases}
\frac{\eta_n}{hT} t^{1-\alpha}, &  t \in [0, hT[,\\
\frac{\eta_n}{hT} \left(t^{1-\alpha} - (t-hT)^{1-\alpha}\right),
   &   t \in [hT, (1-h)T],\\
\frac{\eta_n}{hT} \left(t^{1-\alpha} - (t-hT)^{1-\alpha} -
(t-(1-h)T)^{1-\alpha}\right), &   t \in ](1-h)T, T]
 \end{cases}
\end{equation*}
 and
\begin{align*}
&\int_0^T (_0 D_t^{\alpha} v_n(t))^2 dt \\
& = \int_0^{hT} + \int_{hT}^{(1-h)T}
 + \int_{(1-h)T}^T (_0 D_t^{\alpha} v_n(t)^2 dt
\\
& = \frac{\eta_n^2}{h^2 T^2}\Big[\int_0^T t^{2(1-\alpha)} dt +
\int_{hT}^T  (t - hT)^{2(1-\alpha)} dt  \\
& \quad    + \int_{(1-h)T}^T  (t -(1-h)T)^{2(1-\alpha)} dt  \\
& \quad  - 2 \int_{hT}^T t^{1-\alpha} (t-hT)^{1-\alpha} dt -
2 \int_{(1-h)T}^T  t^{1-\alpha} (t-(1-h)T)^{1-\alpha} dt \\
&  \quad   + 2 \int_{(1-h)T}^T (t - hT)^{1-\alpha}
(t-(1-h)T)^{1-\alpha} dt \Big] \\
& = \frac{\eta_n^2 }{ h^2 T^2}\Big[\frac{1 + h^{3-2\alpha} +
(1-h)^{3-2\alpha}}{3-2\alpha} T^{3-2\alpha}  \\
& \quad  - 2 \int_{(1-h)T}^T t^{1-\alpha}
(t-(1-h)T)^{1-\alpha} dt  -  2 \int_{hT}^T t^{1-\alpha}
(t-hT)^{1-\alpha} dt  \\
 & \quad   + 2 \int_{(1-h)T}^T (t - hT)^{1-\alpha}
(t-(1-h)T)^{1-\alpha} dt \Big] \\
& = 2 A(\alpha, h) \eta_n^2,
\end{align*}
for each $n \in \mathbb{N}$. Thus,  $v_n \in E_0^{\alpha}$, and
\begin{equation}\label{eq:3.10}
 \Phi(v_n) = \frac{1}{2} \|v_n\|_{\alpha}^2 = A(\alpha, h)\eta_n^2.
\end{equation}
 Putting together \eqref{H2} and the nonnegative of $f$,  one has
\begin{equation}\label{eq:3.11}
\begin{aligned}
 \Psi(v_n) & = \int_0^T \big[a(t)
F(v_n(t)) + \frac{\mu}{\lambda} G(t, v_n(t))\big] dt \\
 &  \ge \int_0^T a(t) F(v_n(t))dt \\
&  = \int_0^{hT} a(t) F\Big(\frac{ \Gamma(2-\alpha) \eta_n}{hT}
t\Big) dt + \int_{hT}^{(1-h)T} a(t) F\Big(\Gamma(2-\alpha)
\eta_n\Big) dt  \\
& \quad + \int_{(1-h)T}^T  a(t) F\Big(\frac{ \Gamma(2-\alpha)
\eta_n}{hT} (T-t)\Big) dt \\
&  \ge F(\Gamma(2-\alpha) \eta_n) \int_{hT}^{(1-h)T}a(t) dt.
\end{aligned}
\end{equation}
 Therefore,  from \eqref{eq:3.10} and \eqref{eq:3.11} we achieve
 $$
 \Phi(v_n) - \lambda \Psi(v_n) \le A(\alpha, h) \eta_n^2 - \lambda
F\left(\Gamma(2-\alpha) \eta_n\right)
 \int_{hT}^{(1-h)T} a(t) dt.
$$
 From  \eqref{eq:3.8}, we know that $\lambda _1 < \lambda_2$.
Let
\begin{equation}\label{eq:3.12}
 B = \limsup _{\xi \to + \infty} \frac{F(\xi)}{\xi^2}.
\end{equation}

If  $B < + \infty$, we set $\epsilon \in \big]0, B -
\frac{A(\alpha, h)}{\lambda  \Gamma^2(2-\alpha)
 \int_{hT}^{(1-h)T} a(t) dt} \big[$,  then from \eqref{eq:3.12}
  there exists $N_1$ such that
 $$
F\left(\Gamma(2-\alpha) \eta_n\right)
> (B - \epsilon) \Gamma^2(2-\alpha) \eta_n^2, \quad \forall n >N_1.
$$
 Hence,
\begin{equation}\label{eq:3.13}
\begin{aligned}
 \Phi(v_n) - \lambda \Psi(v_n)
& < A(\alpha, h) \eta_n^2 - \lambda (B
- \epsilon) \Gamma^2(2-\alpha) \eta_n^2 \int_{hT}^{(1-h)T} a(t) dt
\\
&  = \eta_n^2 \Big(A(\alpha, h) - \lambda (B - \epsilon)
\Gamma^2(2-\alpha) \int_{hT}^{(1-h)T} a(t) dt\Big),
\end{aligned}
\end{equation}
for  $n >N_1$.   From the choice of $\epsilon$, we have
 $$
 \lim _{n \to + \infty}(\Phi(v_n) - \lambda \Psi(v_n)) = - \infty.
$$

On the other hand, if  $B = + \infty$, we fix
$\Theta > \frac{A(\alpha, h)}{\lambda  \Gamma^2(2-\alpha)
  \int_{hT}^{(1-h)T} a(t) dt}$, then from \eqref{eq:3.12} there exists $N_{\Theta}$
  such that
$$
F(\Gamma(2-\alpha) \eta_n)> \Theta \Gamma^2(2-\alpha) \eta_n^2, \quad \forall n >
N_{\Theta}.
$$
 Therefore,
\begin{align*}
 \Phi(v_n) - \lambda \Psi(v_n)
& \le A(\alpha, h) \eta_n^2 - \lambda
F\left(\Gamma(2-\alpha) \eta_n\right) \int_{hT}^{(1-h)T} a(t) dt
\\
&  < A(\alpha, h) \eta_n^2  - \lambda \Theta \Gamma^2(2-\alpha)
\eta_n^2 \int_{hT}^{(1-h)T} a(t) dt \\
 & = \eta_n^2 \Big(A(\alpha,h)  - \lambda \Theta \Gamma^2(2-\alpha)
 \int_{hT}^{(1-h)T} a(t) dt\Big), \quad
 n > N_{\Theta}.
\end{align*}
 Taking into account the choice of $\Theta$, one has
$$
 \lim _{n \to + \infty}(\Phi(v_n) - \lambda \Psi(v_n)) = - \infty.
$$
 By Theorem \ref{thm:2.2}, the functional $\Phi - \lambda \Psi$
admits a sequence $u_n$ of critical points such that $\lim_{n \to
+ \infty} \Phi(u_n) = + \infty$. It follows from \eqref{eq:3.1} that
 $$
\|u_n\|_{\alpha} = \sqrt{2 \Phi(u_n)},
$$
 which implies  $\lim_{n \to + \infty} \|u_n\|_{\alpha} = + \infty $.
This completes the proof
in view of the relation between the critical points of $\Phi -
\lambda \Psi$ and the solutions of problem (1.1) pointed out in
Theorem \ref{thm:2.1}.
 \end{proof}


In the following, arguing in a similar way, but applying case (2)
of Theorem \ref{thm:2.2}, we can establish the  existence of
infinitely many solutions of \eqref{eq:1.1} converging at
zero.
For convenience, let
\begin{equation}\label{eq:3.14}
\lambda_3 = \begin{cases}
\frac{A(\alpha, h)}{\Gamma^2(2-\alpha) \int_{hT}^{(1-h)T} a(t) dt
\cdot \lim\sup _{\xi \to 0^+} \frac{F(\xi)}{\xi^2}},
 &  \text{if }   \limsup _{\xi \to 0^+} \frac{F(\xi)}{\xi^2} < + \infty, \\
0, & \text{if }   \limsup _{\xi \to 0^+}
\frac{F(\xi)}{\xi^2} = + \infty,
\end{cases}
\end{equation}
\begin{equation}\label{eq:3.15}
\lambda_4 = \frac{1}{2 M^2 \int_0^T a(t) dt \cdot \lim\inf
_{\xi \to 0^+} \frac{F(\xi)}{\xi^2} }.
\end{equation}

\begin{theorem}\label{thm:3.2}
  Let $1/2 < \alpha \le 1$, $0 < h < 1/2$,
 $a : [0, T] \to \mathbb{R}$ and $f : \mathbb{R} \to \mathbb{R}$ be
 two nonnegative continuous functions, and
assume that % \label{H3}
\begin{equation}\label{eq:3.16}
0 <  \liminf _{\xi \to 0^+} \frac{F(\xi)}{\xi^2} < K
\limsup _{\xi \to 0^+} \frac{F(\xi)}{\xi^2},
\end{equation}
 where $K$ is given in \eqref{eq:3.5}.
For every $\lambda \in
\Lambda_1 :=  ]\lambda_3, \lambda_4[$ ($\lambda_3$ and $\lambda_4$
are given in \eqref{eq:3.14} and \eqref{eq:3.15} respectively)
 and for every $g \in C([0, T] \times \mathbb{R})$ such that
\begin{equation} \label{H4}
\begin{gathered}
G(t, u) \ge 0\quad\text{for all }(t, u) \in [0, T] \times [0, \tau],
 \text{for  some }\tau > 0,
\\
G_0 := \limsup_{\xi \to   0^+}
\frac{\int_0^T \max_{|x|\le \xi} G(t, x) dt}{\xi^2} < + \infty,
\end{gathered}
\end{equation}
if we put
$$
\mu_* := \frac{1}{2 M^2
 G_0}\Big(1 - 2 M^2\lambda \int_0^T a(t) dt \cdot \lim\inf _{\xi \to 0^+}
 \frac{F(\xi)}{\xi^2}\Big),
$$
 with $\mu_* = + \infty$ when $G_0 = 0$,  then \eqref{eq:1.1} admits a
 sequence $\{u_n\}$  of solutions such that $u_n \to 0$ strongly  in
$E_0^{\alpha}$ for every $\mu \in J:= [0, \mu_*[$.
\end{theorem}

\begin{proof}  Fix $ \lambda \in \Lambda_1$, and pick $\mu \in [0, \mu_*[$. We
want to apply Theorem \ref{thm:2.2}(2), with $X = E_0^{\alpha}$,
and $\Phi, \Psi$ be the functionals defined in \eqref{eq:3.1} and
\eqref{eq:3.2} respectively. Let $k_n$ be a sequence of positive
numbers such that $\lim_{n \to \infty} k_n = 0$ and
 $$  \lim _{n \to \infty} \frac{F(k_n)}{k_n^2}
 = \lim \inf _{\xi \to 0^+} \frac{F(\xi)}{\xi^2}. $$

Putting  $r_n = \frac{k_n^2}{2 M^2}$ for all $n \in \mathbb{N}$
and working as in the proof of Theorem \ref{thm:3.1}, it follows
that $\delta < + \infty$, where $\delta$ is as defined in Theorem
\ref{thm:2.2}, and also $\Lambda_1 \subset \left]0,
\frac{1}{\delta}\right[$. Now we claim that
\begin{equation}\label{eq:3.17}
  \Phi - \lambda \Psi \quad  \text{  does  not  have  a local minimum
 at  zero}.
\end{equation}
Let $\{\eta_n\}$ be a sequence of positive numbers in $]0, \eta[$ ($\eta > 0$)
 such that $\eta_n \to
0$, and $\{v_n\}$ be the sequence in $E_0^{\alpha}$ defined in
\eqref{eq:3.9}. From \eqref{H4} and the nonnegative of $f$ one has that
\eqref{eq:3.11} holds.

By  condition \eqref{eq:3.16}, we know that $\lambda _3 < \lambda_4$. Let
\begin{equation}\label{eq:3.18}
  B_1 = \limsup _{\xi \to + 0^+} \frac{F(\xi)}{\xi^2}.
\end{equation}
    For  the case  :  $B_1 < + \infty$, one has \eqref{eq:3.13} holds. From
the choice of $\varepsilon$, we have
$$
\lim _{n \to + \infty} (\Phi(v_n) - \lambda \Psi(v_n)) < 0
= \Phi(0) - \lambda \Psi(0).
$$
 for each $n \in \mathbb{N}$ large enough, which implies \eqref{eq:3.17} holds
  in view of fact that $\|v_n\|\to 0$. Similarly, for the case
 $B_1 = + \infty$, one has \eqref{eq:3.17} holds.

 Observing that $\min_X \Phi = \Phi(0)$, the conclusion follows
  from Theorem \ref{thm:2.2} case (2).
\end{proof}


Finally,  we give an example to show the effectiveness of the
results obtained here.

\begin{example} \label{examp3.1} \rm
  Let $\alpha = 0.8$ and $T = 1$.
Consider the following  boundary-value problem

\begin{equation}\label{eq:3.19}
\begin{gathered}
_t D_1^{0.8}(_0 D_t^{0.8} u(t)) =
 \lambda a(t) f(u(t)) + \mu g(t, u(t)), \quad\text{a.e. } t \in [0, 1], \\
   u(0) = u(1) = 0,
  \end{gathered}
\end{equation}
where $a : [0, 1] \to \mathbb{R}$, $f: \mathbb{R} \to
\mathbb{R}$  and $g : [0, 1] \times \mathbb{R} \to \mathbb{R}$ are
the nonnegative  and continuous functions defined as follows
\begin{gather*}
a(t) = \begin{cases}
 \frac{t}{3}, &  t \in [0, \frac{1}{3}[\\
-12\big(t-\frac{1}{2}\big)^2 + \frac{4}{9}, & t \in [\frac{1}{3}, \frac{2}{3}[\\
 \frac{1-t}{3}, &  t \in [\frac{2}{3}, 1],
\end{cases}
\\
f(x) = \begin{cases}
|x|(3 + 2 \cos(\ln(|x|)) - \sin(\ln(|x|))), &   x \neq  0 \\
0, &  x = 0,
\end{cases}
\\
 g(t, x) = \begin{cases}
\frac{4}{5} + t \sqrt[3]{x}, &  x \le 1 \\
t + \frac{x}{5} (3 + 2\sin(\ln x) + \cos(\ln x)), &  x > 1.
\end{cases}
\end{gather*}
 Then, for every  $\lambda \in ]5.4503, 5.4911[$ and
$\mu \in [0, 0.8132(1-01821 \lambda)[$,   BVP \eqref{eq:3.19}
admits an unbounded sequence of solutions in $E_0^{0.8}$.
 In fact, we have
\begin{equation*}
 F(x) = \begin{cases}
x^2(\frac{3}{2} +  \cos(\ln x)), &  x > 0 \\
0, &  x = 0 \\
- x^2(\frac{3}{2} +  \cos(\ln(|x|))), &  x < 0,
\end{cases}
\end{equation*}
 and
\begin{equation*}
 G(t, x) = \begin{cases}
 \frac{4}{5}x + \frac{3}{4} t x^{\frac{4}{3}}, &  x \le 1 \\
\frac{1}{2} - \frac{t}{4} + tx + \frac{x^2}{5} (\frac{3}{2} +
\sin(\ln x)),  & x > 1.
\end{cases}
\end{equation*}
 It is easy to verify that
\begin{gather*}
\liminf _{\xi \to + \infty} \frac{F(\xi)}{\xi^2} =
\frac{1}{2}, \quad \limsup _{\xi \to + \infty}
\frac{F(\xi)}{\xi^2} = \frac{5}{2},
\\
 G_{\infty} = \limsup_{\xi \to +\infty}
\frac{\int_0^1 \max_{|x| \le \xi} G(t, x) dt}{\xi^2} =
\frac{1}{2}.
\end{gather*}
  Moreover,  it is easy to calculate that $M = 1.1089$,
$A(\alpha, h) = A(0.8, 1/3) = 1.2762$ (here $h= 1/3$)  and
 $$
\frac{\liminf _{\xi \to + \infty} F(\xi)/\xi^2}
{\limsup _{\xi \to + \infty} F(\xi)/\xi^2}
= 0.2 < 0.2015
= \frac{\Gamma^2(1.2) \int_{1/3}^{2/3} a(t) dt}
{2 M^2 A(0.8,1/3) \int_0^1 a(t) dt} = K,
$$
  which implies that condition \eqref{eq:3.8} holds. Obviously,
condition \eqref{H2} holds.  Thus, by Theorem \ref{thm:3.1}, for each
$\lambda \in ]\lambda_1, \lambda_2[ = ]5.4503, 5.4911[$ and
$\mu \in [0, 0.8132(1-0.1821 \lambda)[$, the problem \eqref{eq:3.19}
has an unbounded sequence of solutions in $E_0^{0.8}$.
\end{example}


\subsection*{Acknowledgments}
The author want to thank the anonymous referees for their valuable
comments and suggestions.  This work is supported by grant BK2011407
from the Natural Science Foundation of Jiangsu Province,
 and by grant 11271364 from the Natural Science Foundation of China.

\begin{thebibliography}{00}

\bibitem{AH}   G. A. Afrouzi,  A. Hadjian;  Infinitely many solutions
for a class of Dirichlet quasilinear elliptic systems,\emph{ J.
Math. Anal. Appl.},  393 (2012), 265--272.

\bibitem{ABH}  R. P. Agarwal, M. Benchohra, S. Hamani;  A survey on
existence results for boundary value problems of nonlinear
fractional differential equations and inclusions,  \emph{Acta Appl.
Math.}, 109 (2010), 973--1033.

\bibitem{AN}   B. Ahmad, J. J. Nieto;  Existence results for a
coupled system of nonlinear fractional differential equations with
three-point boundary conditions, \emph{Comput. Math. Appl.}, 58
(2009), 1838--1843.

\bibitem{B1}  C. Bai;  Positive solutions for nonlinear fractional
differential equations with coefficient that changes sign,
\emph{Nonlinear  Anal.},   64 (2006), 677--685.

\bibitem{B2}  C. Bai; Impulsive periodic boundary value problems
for fractional differential equation involving Riemann-Liouville
sequential fractional derivative,  \emph{J. Math. Anal. Appl.}, 384
(2011), 211--231.

\bibitem{B3}  C. Bai;  Existence of solutions for a nonlinear fractional
boundary value problem via a local minimum theorem, \emph{ Electron.
J. Diff. Equ.}, 2012 (2012), No. 176, pp. 1--9.

\bibitem{B4}  C. Bai;  Existence of positive solutions for boundary
value problems of fractional functional differential equations,
 \emph{Electron. J. Qual. Theory Differ. Equ.},  2010 (2010),  No. 30, 1--14.

\bibitem{BL}  Z. Bai,  H. Lu;  Positive solutions for boundary
value problem of nonlinear fractional differential equation,
\emph{J. Math. Anal. Appl.},   311 (2005), 495--505.

\bibitem{BHN}   M. Benchohra, S. Hamani, S. K. Ntouyas;  Boundary
value problems for differential equations with fractional order and
nonlocal conditions, \emph{Nonlinear Anal.},   71 (2009),
2391--2396.

\bibitem{BM1}  G. Bonanno, G. Molica Bisci;  Infinitely many
solutions for a boundary value problem with discontinuous
nonlinearities, \emph{Bound. Value Probl.},  2009 (2009), 1--20.

\bibitem{BB}    G. Bonanno, B. Di Bella;  Infinitely many solutions
for a fourth-order elastic beam equation, \emph{NoDEA Nonlinear
Differential Equations Appl.},  18 (2011), 357--368.

\bibitem{BM2}   G. Bonanno, G. Molica Bisci;  Infinitely many
solutions for a Dirichlet problem involving the p-Laplacian, \emph{
Proc. Roy. Soc. Edinburgh Sect. A},  140 (2010), 737--752.

\bibitem{BMD}  G. Bonanno, G. Molica Bisci, V. R\v{a}dulescu;
Infinitely many solutions for a class of nonlinear eigenvalue
problem in Orlicz-Sobolev spaces, \emph{Comptes Rendus
Math\'{e}matique}, 349 (2011), 263--268.

\bibitem{BMO}   G. Bonanno, G. Molica Bisci, D. O'Regan;
Infinitely many weak solutions for a class of quasilinear elliptic
systems, \emph{ Math. Comput. Modelling},  52 (2010), 152--160.

\bibitem{BT}  G. Bonanno, E. Tornatore;  Infinitely many
solutions for a mixed boundary value problem, \emph{ Ann. Polon.
Math.},  99 (2010), 285--293.

\bibitem{B5}  G. Bonanno;   A critical point theorem via the
Ekeland variational principle, \emph{ Nonlinear Anal.},  75 (2012),
2992--3007.

\bibitem{BC} G. Bonanno, A. Chinn\'{i}; Existence of three solutions for a perturbed
 two-point boundary value problem,  \emph{ Appl. Math.
Lett.},  23 (2010),  807--811.


\bibitem{BMW} G. Bonanno, D. Motreanu, P. Winkert;  Variational  hemivariational
 inequalities with small perturbations of nonhomogeneous
Neumann boundary conditions,  \emph{ J. Math. Anal. Appl.},  381
(2011), 627--637.

\bibitem{CD}  P. Candito, G. D'Agu\'{i};  Three solutions to a perturbed non-linear
discrete Dirichlet problem,  \emph{ J. Math. Anal. Appl.},  375
(2011), 594--601.

\bibitem{DS} G. D'Agu\'{i}, A. Sciammetta;  Infinitely many solutions to elliptic
problems with variable exponent and nonhomogeneous Neumann condi-
tions,  \emph{ Nonlinear Anal.},  75 (2012), 5612--5619.


\bibitem{DM}  G. D'Agu\'{i}, G. Molica Molica;  Infinitely many solutions for perturbed
hemivariational inequalities,  \emph{ Bound. Value Probl.},  2010,
Art. ID 363518, 19 pp.

\bibitem{HT} S. Heidarkhani, Y. Tian;  Multiplicity results for a class of
gradient systems depending on two parameters, \emph{ Nonlinear
Anal.},  73 (2010), 547--554.

\bibitem{GK}  J. Graef, L. Kong;  Positive solutions for a class
of higher order boundary value problems with fractional
q-derivatives, \emph{ Appl. Math. Comput.},  218 (2012), 9682--9689.

\bibitem{JZ}  F. Jiao, Y. Zhou;  Existence results for fractional
boundary value problem via critical point theory,
\emph{International Journal of Bifurcation and Chaos}, 22 (2012)
1250086.

\bibitem{KST}  A. A. Kilbas, H. M. Srivastava, J. J. Trujillo;
Theory and Applications of Fractional Differential Equations,
Elsevier, Amsterdam, 2006.

\bibitem{K}   N. Kosmatov;  Integral equations and initial value
problems for nonlinear differential equations of fractional order,
\emph{Nonlinear Anal.},  70 (2009), 2521--2529.

\bibitem{LZL}  Y. Liu,  W. Zhang,  X. Liu;  A sufficient
condition for the existence of a positive solution for a nonlinear
fractional differential equation with the Riemann-Liouville
derivative, \emph{Appl. Math. Lett.},  25 (2012), 1986--1992.

\bibitem{P}  I. Podlubny;  Fractional Differential Equations,
Academic Press, San Diego, 1999.

\bibitem{R} B. Ricceri;  A general variational principle and
some of its applications, \emph{J. Comput. Appl. Math.},  113
(2000), 401--410.

\bibitem{SKM}  S. G. Samko, A. A. Kilbas, O. I. Marichev;  Fractional
Integral and Derivatives : Theory and Applications, Gordon and
Breach, Longhorne, PA, 1993.

\bibitem{WZWH}   J. Wang, Y. Zhou, W. Wei, H. Xu;  Nonlocal problems
for fractional integrodifferential equations via fractional
operators and optimal controls, \emph{Comput. Math. Appl.},  62
(2011), 1427--1441.

\bibitem{WZ}  J. Wang, Y. Zhou;  A class of fractional evolution
equations and optimal controls, \emph{Nonlinear Anal. RWA},  12
(2011), 262--272.

\bibitem{Z} S. Zhang;  Positive solutions to singular boundary
value problem for nonlinear fractional differential equation,
\emph{Comput. Math. Appl.},  59 (2010), 1300--1309.

\end{thebibliography}

\end{document}
