JFS Journal of Function Spaces 2314-8888 2314-8896 Hindawi 10.1155/2019/8981528 8981528 Research Article Multiple Solutions for a Class of Fractional Schrödinger-Poisson System https://orcid.org/0000-0002-5169-0094 Chen Lizhen 1 https://orcid.org/0000-0002-1818-5862 Li Anran 2 Wei Chongqing 2 Avery Richard I. 1 School of Applied Mathematics Shanxi University of Finance and Economics Taiyuan 030006 China sxufe.edu.cn 2 School of Mathematical Sciences Shanxi University Taiyuan 030006 China sxu.edu.cn 2019 3172019 2019 08 06 2019 16 07 2019 3172019 2019 Copyright © 2019 Lizhen Chen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

We investigate a class of fractional Schrödinger-Poisson system via variational methods. By using symmetric mountain pass theorem, we prove the existence of multiple solutions. Moreover, by using dual fountain theorem, we prove the above system has a sequence of negative energy solutions, and the corresponding energy values tend to 0. These results extend some known results in previous papers.

National Natural Science Foundation of China 11701346 11571209
1. Introduction

We consider the following system via variational methods:SP-Δαw+Vxw+ϕw=gx,w+λhxwq-2w,xR3,-Δβϕ=w2,xR3,where λ>0, 1<q<2, α,β(0,1], 2β+4α>3. (-Δ)α and (-Δ)β represent the Laplace operator of the fractional order. If α=β=1, then the system SP degenerates into the standard Schrödinger-Poisson system, which describes the interaction between the same charged particles when the magnetic effect can be ignored . In recent years, the existence, multiplicity, and centralization of solutions for the Schrödinger-Poisson system have been deeply studied via variational methods, and a great number of works have been obtained, see, for example, . On the other hand, (-Δ)α is a class of nonlocal pseudo-differential operators. Since nonlocal differential equations can better and more fully describe the physical experimental phenomena than classical local differential operators, the study of nonlinear fractional Laplace equation has become one of the most popular research fields in nonlinear analysis.

In the literature , Wei considered the following system:(1)-Δαw+Vxw+ϕw=gx,w,xR3,-Δαϕ=γαw2,xR3.By using the critical point theory, the author obtained infinitely many solutions when α0,1. In the literature , Teng studied a system of the form(2)-Δαw+Vxw+ϕw=θwq-1w+w2α-2w,xR3,-Δβϕ=w2,xR3,where q(1,3+2α/3-2α), α,β(0,1), 2β+2α>3. In , Zhang, Marcos, and Squassina used perturbation approach to obtain the existence of solutions for the following system when the nonlinear term is subcritical or critical(3)-Δαw+λϕw=gw,xR3,-Δβϕ=λw2,xR3,where λ>0, α,β[0,1]. In , Duarte and Souto investigated the following system via variational methods(4)-Δαw+Vxw+ϕw=gw,xR3,-Δβϕ=w2,xR3,where α(3/4,1), β(0,1), V:R3R is a periodic potential. A positive solution and a ground state solution were got in . In , Li studied a system of the following form:(5)-Δαw+Vxw+ϕw=gx,w,xR3,-Δβϕ=w2,xR3,where α,β(0,1], 2β+4α>3. Combining the perturbation method with mountain pass theorem, the existence of nontrivial solutions was obtained in . In , Yu, Zhao, and Zhao studied the following fractional Schrödinger–Poisson system with critical growth via variational methods(6)ε2α-Δαw+Vxw+ϕw=w2α-2w+gw,xR3,ε2α-Δαϕ=w2,xR3,where α(3/4,1),2α=6/3-2α, the potential V is continuous with positive infimum, g is continuous and subcritical at infinity. Under some Monotone hypothesis on g, the existence of positive ground state solution is got in . For small ε>0, a multiple result is also got in .

Inspired by , in this paper, we prove the existence of multiple solutions for system SP by symmetric mountain pass theorem. Moreover, we prove the system SP has a sequence of negative energy solutions by dual fountain theorem. The assumptions on V and nonlinearity g in this paper are given below:

(V) VC(R3,R), V0infxR3V(x)>0 and lim|x|+V(x)=+;

(H1) gC(R3×R,R), and there exists C1>0 such that |g(x,w)|C1(|w|+|w|p-1), where p(4,2α), 2α=6/3-2α;

(H2) there exist κ>4 and r>0 such that 0<κG(x,w)κ0wg(x,t)dtg(x,w)w, for |w|r. Moreover, infxR3,|w|=rG(x,w)>0;

(H3) g(x,-w)=-g(x,w), xR3, wR;

(H4) h:R3R+, and hL2/2-q(R3).

2. Preliminaries

For 1ν<, Lν(R3) denotes the usual Lebesgue space with norm wν=(R3|w|νdx)1/ν. Fix α(0,1), fractional Sobolev space Hα(R3) denoted as(7)HαR3wL2R3:R31+l2αFwl2dl<,equipped with the norm(8)wHα=R3Fwl2+l2αFwl2dl1/2,where Fw denotes the Fourier transform of function w. Let gC0(R3); the fractional Laplacian operator (-Δ)α:C0(R3)C0R3 is defined by(9)-Δαg=F-1l2αFg,lR3.According to Plancherel theorem , one has Fw2=w2, |l|αFw2=(-Δ)α/2w2. By (8), we define the equivalent norm(10)wHα=R3-Δα/2wx2+wx2dx1/2.Dα,2(R3) is denoted as(11)Dα,2R3=wL2αR3:lαFwlL2R3.In particular, Dα,2(R3) is the completion of C0(R3), with respect to the norm(12)wDα,2=R3-Δα/2w2dx1/2=R3l2αFwl2dl1/2.Lq(R3,h), q(1,2) represents a weighted Lebesgue space, that is,(13)LqR3,h=w:R3RismeasurableandR3hxwqdx<+,equipped with the norm(14)wLqR3,h=R3hxwqdx1/q.

For convenience, we use C to represent any positive constants which may change from line to line. According to , the embedding Hα(R3)Lν(R3) is continuous for all ν[2,2α], i.e., there exists Mν>0 satisfying(15)wνMνwHα,wHαR3.So, the embedding Hα(R3)Lν(R3) is continuous when 2β+4α>3. Fix wHα(R3); we define the nonlinear operator Lw:Dβ,2(R3)R by(16)Lwv=R3w2vdx.Hence(17)LwvR3wx12/3+2βdx3+2β/6R3vx2βdx1/2βCR3wx12/3+2βdx3+2β/6vDβ,2Cw2vDβ,2.By Lax-Milgram theorem, we can find a unique ϕwβDβ,2(R3) such that(18)R3-Δβ/2ϕwβ-Δβ/2vdx=R3w2vdx,vDβ,2R3,and ϕwβ is expressed as(19)ϕwβx=cβR3w2yx-y3-2βdy,wherecβ=Γ3/2-2βπ3/222βΓβ.According to (19), ϕwβ0 for all xR3. Since β(0,1], 2β+4α>3, we can also get 12/3+2β(2,2α). Together with (17) and (18),(20)ϕwβDβ,2=R3-Δβ/2ϕwβ2dxCR3wx12/3+2βdx3+2β/6Cw2.From Hölder’s inequality and (19),(21)R3ϕwβw2dxCR3wx12/3+2αdx3+2α/6ϕwβDβ,2Cw2ϕwβDβ,2Cw4.Evidently,(22)R3ϕwβw2dxCw4.Substituting (18) into system SP, system SP is equivalent toS-Δαw+Vxw+ϕwβw=gx,w+λhxwq-2w.

For the equation S, we define the work space E as(23)EwHαR3:R3Vxw2dx<,with(24)w=R3-Δα/2w2+Vxw2dx1/2.

Lemma 1 (see [<xref ref-type="bibr" rid="B19">19</xref>]).

Assume condition (V) holds, α(0,1), ν[2,2α); then the embedding ELν(R3) is compact.

From (V) and (H1)-(H4), I:ER(25)Iw=12R3-Δα/2w2+Vxw2dx+14R3ϕwβw2dx-R3Gx,wdx-λqR3hxwqdx.is well defined. Moreover, by Lemma 1, IC1(E,R) with(26)Iw,v=R3-Δα/2w-Δα/2v+Vxwvdx+R3ϕwβwvdx-R3gx,wvdx-λR3hxwq-2wvdx,vE

Proposition 2.

(i) If equation S has a solution wE, then system SP has a solution (w,ϕ)E×Dβ,2(R3).

(ii) If for every vE, the following equation(27)R3-Δα/2w-Δα/2v+Vxwvdx+R3ϕwβwvdx-R3gx,wvdx-λR3hxwq-2wvdx=0holds, then wE is a solution of S.

Set {ei}i=1 as a set of normalized orthogonal basis of E, Xi=Rei. Yk=i=1kXi,Zk=i=k+1Xi¯,kN. Obviously, E=YkZk.

Definition 3 (see [<xref ref-type="bibr" rid="B20">20</xref>, <xref ref-type="bibr" rid="B21">21</xref>]).

Set IC1(E,R), cR. If any sequence {wn}E satisfying(28)Iwnc,Iwn0,n,has a convergent subsequence, then I satisfies the (PS)c condition. Any sequence satisfying (28) is called the (PS)c sequence.

Definition 4 (see [<xref ref-type="bibr" rid="B21">21</xref>, <xref ref-type="bibr" rid="B22">22</xref>]).

Set IC1(E,R), cR. If every sequence {wnj}E, satisfying(29)wnjYnj,Iwnjc,IYnj0,n,has a convergent subsequence, then I satisfies the (PS)c condition.

Proposition 5 (see [<xref ref-type="bibr" rid="B23">23</xref>]).

Let E=YZ be a Banach space with dimY<. Assume IC1(E,R) is an even functional and satisfies the (PS)c condition and

(A1) there exist ω, μ>0 satisfying IBωZ=infwZ,w=ωI(w)μ;

(A2) for every linear subspace UE with dimU<, there exists a constant L=L(U) such that maxwU,wLI(w)<0,

then I has a list of unbounded critical points.

Proposition 6 (see [<xref ref-type="bibr" rid="B22">22</xref>]).

Assume that IC1(E,R) is an even functional, k0N. If for any k>k0, there exist rk>γk>0 such that

bk=infIw:wZk,w=rk0;

ak=max{I(w):wYk,w=γk}<0;

ck=inf{I(w):wZk,wrk}0, k;

for every c[ck0,0], I satisfies the (PS)c condition,

then I has a sequence of negative critical points that converge to 0.

3. Main Results Lemma 7.

If hypotheses (V) and (H1)-(H4) hold, then for any cR, I satisfies the (PS)c condition.

Proof.

First, we prove the (PS)c sequence {wn} of I is bounded. According to (H4), it is easy to get that(30)R3hxwnqdxR3hx2/2-qdx2-q/2R3wn2dxq/2h2/2-qw2qM2qh2/2-qwnq.By condition (H2), there exists r>0 such that(31)gx,wwκGx,w,wr.Moreover, for any given C0(0,1/16V0), we can choose a constant δ>0 such that(32)1κgx,ww-Gx,wC0w2,forwδ.From condition (H1), when δ|u|r, one has(33)1κgx,ww-Gx,wC1δ+rp-2w2.So, for any |w|r,(34)1κgx,ww-Gx,wC0w2+C1δ+rp-2w2.For lim|x|+V(x)=+, there exists L>r>0 such that(35)116Vx>C1δ+rp-2,xL.Now (34) implies(36)14R3Vxwn2dx+wnxr1κgx,wnwn-Gx,wndx14R3Vxwn2dx-wnxrC0wn2+C1δ+rp-2wn2dx=18R3Vxwn2dx-wnxrC0wn2dx+18R3Vxwn2dx-wnxrC1δ+rp-2wn2dxwnxr116V0-C0wn2dx+18R3Vxwn2dx-C1δ+rp-2r2·measxR3xL18R3Vxwn2dx-C1δ+rp-2r2·measxR3xL.Since {wn} is the (PS)c sequence, when n is large enough,(37)c+wnIwn-1κIwn,wn=12-1κR3wn2+Vxwn2dx+14-1κR3ϕwn2wn2dx+R31κgx,wnwn-Gx,wndx+1κ-1qλR3hxwnqdx14R3wn2+Vxwn2dx+R31κgx,wnwn-Gx,wndx-1q-1κλM2qh2/2-qwnq14R3wn2dx+14R3Vxwn2dx+wnxr1κgx,wnwn-Gx,wndx-1q-1κλM2qh2/2-qwnq.Thus, according to (36), when n is large enough,(38)c+wn14R3wn2dx+14R3Vxwn2dx+wnxr1κgx,wnwn-Gx,wndx-1q-1κλM2qh2/2-qwnq18wn2-C1δ+rp-2r2·measxR3xL-1q-1κλM2qh2/2-qwnq.Therefore, {wn} is a bounded sequence in E. By Lemma 1, there exists w^E such that(39)wnw^inE,wnw^inLνR3,wnxw^xa.e.  onR3.

Next, we define the linear operator Bφ:ER as(40)Bφv=-Δα/2φ-Δα/2v.From Hölder’s inequality, we obtain(41)Bφvφv,vE.Now by Lemma 1 and (22),(42)R3ϕwnβwnwn-w^dxϕwnβ2βwn12/3+2βwn-w^12/3+2βCϕwnβDβ,2wn12/3+2βwn-w^12/3+2βCwn12/3+2β3wn-w^12/3+2βCwn3wn-w^12/3+2β.Similarly, we can also prove(43)R3ϕw^βw^wn-w^dxCw^3wn-w^12/3+2β.Since wnw^ in Lν(R3)(ν[2,2α)), limnR3(ϕwnβwn-ϕw^βw^)(wn-w^)dx=0.

At last, combining Hölder’s inequality with (H1) and (H4), we can easily get(44)limnR3gx,wn-gx,w^wn-w^dx=0,limnR3hxwnq-2wn-hxw^q-2w^wn-w^dx=0.Thus(45)o1=Iwn-Iw^,wn-w^=Bwnwn-w^-Bwwn-w^+R3Vxwnwn-w^-Vxw^wn-w^dx+R3gx,wn-gx,w^wn-w^dx+λR3hxwnq-2wn-w^q-2wwn-w^dx+R3ϕwnβwn-ϕw^βw^wn-w^dx=Bwnwn-w^-Bwwn-w^+R3Vxwnwn-w^-Vxw^wn-w^dx+o1,that is,(46)wn-w^2=Bwnwn-w^-Bwwn-w^+R3Vxwn-w^2dx0.

Lemma 8.

If hypotheses (V) and (H1)-(H4) hold, then I satisfies (PS)c condition for all cR.

Proof.

By Definition 4, we just prove the following fact: if for any cR, {wnj}E, and wnjYnj,I(wnj)c,IYnj0, as nj, then {wnj} has a convergence subsequence. The proof method is similar to Lemma 7.

Lemma 9 (see [<xref ref-type="bibr" rid="B24">24</xref>]).

For 2ν<2α, kN, set(47)βνksupwν:wZk,w=1,and then βν(k)0, k.

Theorem 10.

If hypotheses (V) and (H1)-(H4) hold, then we can find λ0>0, such that system (SP) has multiple solutions for every λ<λ0. Moreover, the corresponding energy values tend to infinity.

Proof.

According to Lemma 7, I satisfies (PS)c condition. We only need to prove that I satisfies (A1) and (A2). By virtue of (H1),(48)Gx,wC12w2+C1pwp,x,wR3×R.From Lemma 9, we can get(49)Iw=12R3-Δα/2w2+Vxw2dx+14R3ϕwβw2dx-R3Gx,wdx-λqR3hxwqdx12w2-C12R3w2dx-C1pR3wpdx-M2qλh2/2-qwq12w2-C12β22kw2-C1pMppwp-M2qλh2/2-qwqw212-C12β22k-C1Mppwp-2-M2qλh2/2-qwq-2.Take a sufficiently large k such that β22(k)<1/2C1. Combining the above inequality, we obtain(50)Iww214-C1Mppwp-2-M2qλh2/2-qwq-2.Set(51)ηt=14-C1Mpptp-2-M2qλh2/2-qtq-2,t>0.Since 1<q<2<p, there exists(52)ωλλ2-qM2qh2/2-qC1Mppp-21/p-q>0,such that maxtR+η(t)=η(ωλ). Therefore, foreveryλ<λ0(2-q/4C1Mpp(p-q))p-q/p-2·C1(p-2)Mpp/M2q(2-q)h2/2-q,(53)Iwωλ2ηωλμ>0,withw=ωλ.

On the other hand, by conditions (H1) and (H2), there exist positive constants C2,C3 such that(54)Gx,wC2wκ-C3w2,x,wR3×R.Since all the norms are equivalent in every finite linear subspace UE, then for wU(55)Iw=12R3-Δα/2w2+Vxw2dx+14R3ϕwβw2dx-R3Gx,wdx-λqR3hxwqdx12w2+Cw4-C2wκκ-C3w22-λqwLqR3,hq.For q<2<4<κ, I(w)- as w. Then there exists L=L(U)>0 such that maxwU,wLI(w)<0. Thus, according to Proposition 5, the system (SP) has a list of solutions {(wn,ϕn)}E×Dβ,2(R3), and the corresponding energy values tend to infinity.

Theorem 11.

If hypotheses (V) and (H1)-(H4) hold, then the system (SP) has a sequence of negative energy solutions for all λ>0, and the energy values tend to 0.

Proof.

By Lemma 8, for all cR, I satisfies the (PS)c condition. It now remains to show that (C1)-(C3) are satisfied. According to Lemma 9, for every ν[2,2α), βν(k)0, as k. Thus there exists k1>0 such that β2(k)1/2C1 for k>k1. For 4<p<2α, there exists L(0,1) such that(56)18w2C1pMppwp,withwL.

Hence, for wZk with wL, it follows that(57)Iw=12R3-Δα/2w2+Vxw2dx+14R3ϕwβw2dx-R3Gx,wdx-λqR3hxwqdx12w2-C12R3w2dx-C1pR3wpdx-λqR3hxwqdx12w2-C12β22ku2-C1pMppwp-λqR3hxwqdx14w2-C1pMppwp-λqh2/2-qβ2qkwq18w2-λqh2/2-qβ2qkwq.For every k>k1, let rk=(8/qλh2/2-qβ2q(k))1/(2-q). By Lemma 9, rk0, as k. Thus, there exists k0>k1, such that for every kk0, I(w)0, for wZk with w=rk.

Secondly, since for every fixed kN, the norms are equivalent in Yk, when k is sufficiently large, there exists a small enough γk such that 0<γk<rk and I(w)<0 for wYk with w=γk.

Finally, according to (C3), when kk0, for uZk, with wrk, one has(58)Iw-λqh2/2-qβ2qkwq-λqh2/2-qβ2qkrkq.Since β2(k)0, rk0, as k, therefore (C3) holds. By Proposition 6, I has a list of solutions {(wn,ϕn)}E×Dβ,2(R3) such that(59)12R3-Δα/2wn2+Vxwn2dx+14R3ϕnβwn2dx-R3Gx,wndx-λqR3hxwnqdx0.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that they have no competing interests.

Authors’ Contributions

All authors read and approved the final manuscript.

Acknowledgments

This work is supported by the National Natural Science Foundation of China (Grant No. 11701346, 11571209).

Benci V. Fortunato D. An eigenvalue problem for the Schrödinger-Maxwell equations Topological Methods in Nonlinear Analysis 1998 11 2 283 293 10.12775/TMNA.1998.019 MR1659454 Zbl0926.35125 D'Aprile T. Mugnai D. Solitary waves for nonlinear Klein-Gordon-Maxwell and Schrödinger-Maxwell equations Proceedings of the Royal Society of Edinburgh, Section: A Mathematics 2004 134 5 893 906 10.1017/s030821050000353x MR2099569 Ruiz D. The Schrödinger-Poisson equation under the effect of a nonlinear local term Journal of Functional Analysis 2006 237 2 655 674 10.1016/j.jfa.2006.04.005 MR2230354 Zbl1136.35037 2-s2.0-33646917496 Mao A. Yang L. Qian A. Luan S. Existence and concentration of solutions of Schrödinger-Possion system Applied Mathematics Letters 2017 68 8 12 10.1016/j.aml.2016.12.014 MR3614271 Zhao L. Zhao F. Positive solutions for Schrödinger–Poisson equations with a critical exponent Nonlinear Analysis: Theory, Methods & Applications 2009 70 6 2150 2164 10.1016/j.na.2008.02.116 Zbl1156.35374 Sun M. Su J. Zhao L. Solutions of a Schrödinger-Poisson system with combined nonlinearities Journal of Mathematical Analysis and Applications 2016 442 2 385 403 10.1016/j.jmaa.2016.04.053 MR3504005 Zbl1343.35097 Shao M. Mao A. Multiplicity of solutions to Schrödinger-Poisson system with concave-convex nonlinearities Applied Mathematics Letters 2018 83 212 218 10.1016/j.aml.2018.04.005 MR3795694 Zbl06892582 Wang L. Ma S. Wang X. On the existence of solutions for nonhomogeneous Schrödinger-Poisson system Boundary Value Problems 2016 2016 1, article no. 76 10.1186/s13661-016-0584-9 Zbl1341.35043 Wei Z. L. Existence of Infinitely Many Solutions for the Fractional Schrödinger-Maxwell Eqautions 2015 https://arxiv.org/abs/1508.03088 Teng K. Existence of ground state solutions for the nonlinear fractional Schrödinger–Poisson system with critical Sobolev exponent Journal of Differential Equations 2016 261 6 3061 3106 10.1016/j.jde.2016.05.022 Zbl1386.35458 Zhang J. do Ó J. M. Squassina M. Fractional Schrödinger-Poisson systems with a general subcritical or critical nonlinearity Advanced Nonlinear Studies 2016 16 1 15 30 10.1515/ans-2015-5024 MR3456743 Duarte R. C. Souto M. A. Fractional Schrödinger-Poisson equations with general nonlinearities Electronic Journal of Differential Equations 2016 319 1 19 Li K. Existence of non-trivial solutions for nonlinear fractional Schrödinger-Poisson equations Applied Mathematics Letters 2017 72 1 9 10.1016/j.aml.2017.03.023 MR3652422 Zbl1398.35272 Yu Y. Zhao F. Zhao L. The existence and multiplicity of solutions of a fractional Schrodinger-Poisson with critical exponent Science China Mathematics 2018 61 6 1039 1062 10.1007/s11425-016-9074-6 MR3802758 Xiang M. Zhang B. Guo X. Infinitely many solutions for a fractional Kirchhoff type problem via fountain theorem Nonlinear Analysis. Theory, Methods & Applications. An International Multidisciplinary Journal 2015 120 299 313 10.1016/j.na.2015.03.015 MR3348060 Zbl1328.35287 Li Q. Q. Wu X. A new result on high energy solutions for Schrödinger-Kirchhoff type equations in RN Applied Mathematics Letters 2014 30 24 27 10.1016/j.aml.2013.12.002 Di Nezza E. Palatucci G. Valdinoci E. Hitchhiker's guide to the fractional Sobolev spaces Bulletin des Sciences Mathématiques 2012 136 5 521 573 10.1016/j.bulsci.2011.12.004 MR2944369 Pucci P. Xiang M. Zhang B. Multiple solutions for nonhomogeneous Schrödinger-Kirchhoff type equations involving the fractional p-Laplacian in RN Calculus of Variations and Partial Differential Equations 2015 54 3 2785 2806 10.1007/s00526-015-0883-5 MR3412392 Teng K. Multiple solutions for a class of fractional Schrodinger equations in RN Nonlinear Analysis: Real World Applications 2015 21 76 86 10.1016/j.nonrwa.2014.06.008 MR3261580 Brézis H. Coron J. Nirenberg L. Free vibrations for a nonlinear wave equation and a theorem of P. Rabinowitz Communications on Pure and Applied Mathematics 1980 33 5 667 684 10.1002/cpa.3160330507 Zbl0484.35057 Willem M. Minimax Theorems 1996 Boston, Mass, USA Birkhäuser 10.1007/978-1-4612-4146-1 MR1400007 Zbl0856.49001 Bartsch T. Willem M. On an elliptic equation with concave and convex nonlinearities Proceedings of the American Mathematical Society 1995 123 11 3555 3561 10.2307/2161107 MR1301008 Zbl0848.35039 Bartolo P. Benci V. Fortunato D. Abstract critical point theorems and applications to some nonlinear problems with ‘strong’ resonance at infinity Nonlinear Analysis. Theory, Methods & Applications 1983 7 9 981 1012 10.1016/0362-546x(83)90115-3 MR713209 Zbl0522.58012 2-s2.0-0009179578 Binlin Z. Molica Bisci G. Servadei R. Superlinear nonlocal fractional problems with infinitely many solutions Nonlinearity 2015 28 7 2247 2264 MR3366642 10.1088/0951-7715/28/7/2247 Zbl1322.35158 2-s2.0-84933073626