Loading [MathJax]/jax/output/SVG/jax.js
Research article Special Issues

Least energy sign-changing solutions for Kirchhoff-Schrödinger-Poisson system on bounded domains

  • Received: 28 October 2022 Revised: 16 February 2023 Accepted: 20 February 2023 Published: 20 March 2023
  • We investigate the following nonlinear system

    {(a+bΩ|u|2dx)Δu+ϕu=λu+μ|u|2u, xΩ,Δϕ=u2, xΩ,u=ϕ=0, xΩ,

    with a,b>0, λ,μR, and ΩR3 is bounded with smooth boundary. Let λ1>0 be the first eigenvalue of (Δu,H10(Ω)). We get that for certain ˜μ>0 there exists at least one least energy sign-changing solution for the above system if λ<aλ1 and μ>˜μ. In addition, we remark that the nonlinearity λu+μ|u|2u does not satisfy the growth conditions.

    Citation: Xia Su, Wen Guan, Xia Li. Least energy sign-changing solutions for Kirchhoff-Schrödinger-Poisson system on bounded domains[J]. Electronic Research Archive, 2023, 31(5): 2959-2973. doi: 10.3934/era.2023149

    Related Papers:

    [1] Gilberto González-Parra, Cristina-Luisovna Pérez, Marcos Llamazares, Rafael-J. Villanueva, Jesus Villegas-Villanueva . Challenges in the mathematical modeling of the spatial diffusion of SARS-CoV-2 in Chile. Mathematical Biosciences and Engineering, 2025, 22(7): 1680-1721. doi: 10.3934/mbe.2025062
    [2] Fang Wang, Lianying Cao, Xiaoji Song . Mathematical modeling of mutated COVID-19 transmission with quarantine, isolation and vaccination. Mathematical Biosciences and Engineering, 2022, 19(8): 8035-8056. doi: 10.3934/mbe.2022376
    [3] Yangyang Yu, Yuan Liu, Shi Zhao, Daihai He . A simple model to estimate the transmissibility of the Beta, Delta, and Omicron variants of SARS-COV-2 in South Africa. Mathematical Biosciences and Engineering, 2022, 19(10): 10361-10373. doi: 10.3934/mbe.2022485
    [4] Sarafa A. Iyaniwura, Rabiu Musa, Jude D. Kong . A generalized distributed delay model of COVID-19: An endemic model with immunity waning. Mathematical Biosciences and Engineering, 2023, 20(3): 5379-5412. doi: 10.3934/mbe.2023249
    [5] Anichur Rahman, Muaz Rahman, Dipanjali Kundu, Md Razaul Karim, Shahab S. Band, Mehdi Sookhak . Study on IoT for SARS-CoV-2 with healthcare: present and future perspective. Mathematical Biosciences and Engineering, 2021, 18(6): 9697-9726. doi: 10.3934/mbe.2021475
    [6] Matthew Hayden, Bryce Morrow, Wesley Yang, Jin Wang . Quantifying the role of airborne transmission in the spread of COVID-19. Mathematical Biosciences and Engineering, 2023, 20(1): 587-612. doi: 10.3934/mbe.2023027
    [7] A. M. Elaiw, Raghad S. Alsulami, A. D. Hobiny . Global dynamics of IAV/SARS-CoV-2 coinfection model with eclipse phase and antibody immunity. Mathematical Biosciences and Engineering, 2023, 20(2): 3873-3917. doi: 10.3934/mbe.2023182
    [8] Rahat Zarin, Usa Wannasingha Humphries, Amir Khan, Aeshah A. Raezah . Computational modeling of fractional COVID-19 model by Haar wavelet collocation Methods with real data. Mathematical Biosciences and Engineering, 2023, 20(6): 11281-11312. doi: 10.3934/mbe.2023500
    [9] Salman Safdar, Calistus N. Ngonghala, Abba B. Gumel . Mathematical assessment of the role of waning and boosting immunity against the BA.1 Omicron variant in the United States. Mathematical Biosciences and Engineering, 2023, 20(1): 179-212. doi: 10.3934/mbe.2023009
    [10] Hiroshi Nishiura, Seiko Fujiwara, Akifumi Imamura, Takuma Shirasaka . HIV incidence before and during the COVID-19 pandemic in Japan. Mathematical Biosciences and Engineering, 2024, 21(4): 4874-4885. doi: 10.3934/mbe.2024215
  • We investigate the following nonlinear system

    {(a+bΩ|u|2dx)Δu+ϕu=λu+μ|u|2u, xΩ,Δϕ=u2, xΩ,u=ϕ=0, xΩ,

    with a,b>0, λ,μR, and ΩR3 is bounded with smooth boundary. Let λ1>0 be the first eigenvalue of (Δu,H10(Ω)). We get that for certain ˜μ>0 there exists at least one least energy sign-changing solution for the above system if λ<aλ1 and μ>˜μ. In addition, we remark that the nonlinearity λu+μ|u|2u does not satisfy the growth conditions.



    In this paper, we are concerned with the existence and multiplicity of weak solutions for the damped-like fractional differential system

    {ddt(p(t)(12 0Dξt(u(t))+12 tDξT(u(t))))+r(t)(12 0Dξt(u(t))+12 tDξT(u(t)))+q(t)u(t)=λF(t,u(t)),  a.e. t[0,T],u(0)=u(T)=0, (1.1)

    where 0Dξt and tDξT are the left and right Riemann-Liouville fractional integrals of order 0ξ<1, respectively, p,r,qC([0,T],R), L(t):=t0(r(s)/p(s))ds, 0<meL(t)p(t)M and q(t)p(t)0 for a.e. t[0,T], u(t)=(u1(t),u2(t),un(t))T, ()T denotes the transpose of a vector, n1 is a given positive integer, λ>0 is a parameter, F(t,x) is the gradient of F with respect to x=(x1,,xn)Rn, that is, F(t,x)=(Fx1,,Fxn)T, and there exists a constant δ(0,1) such that F:[0,T]ׯBδ0R (where ¯Bδ0 is a closed ball in RN with center at 0 and radius δ) satisfies the following condition

    (F0) F(t,x) is continuously differentiable in ¯Bδ0 for a.e. t[0,T], measurable in t for every x¯Bδ0, and there are aC(¯Bδ0,R+) and bL1([0,T];R+) such that

    |F(t,x)|a(|x|)b(t)

    and

    |F(t,x)|a(|x|)b(t)

    for all x¯Bδ0 and a.e. t[0,T].

    In recent years, critical point theory has been extensively applied to investigate the existence and multiplicity of fractional differential equations. An successful application to ordinary fractional differential equations with Riemann-Liouville fractional integrals was first given by [1], in which they considered the system

    {ddt((12 0Dξt(u(t))+12 tDξT(u(t))))=F(t,u(t)),  a.e. t[0,T],u(0)=u(T)=0. (1.2)

    They established the variational structure and then obtained some existence results for system (1.2). Subsequently, this topic attracted lots of attention and a series of existence and multiplicity results are established (for example, see [2,3,4,5,6,7,8,9,10,11,12] and reference therein). It is obvious that system (1.1) is more complicated than system (1.2) because of the appearance of damped-like term

    r(t)(12 0Dξt(u(t))+12 tDξT(u(t))).

    In [13], the variational functional for system (1.1) with λ=1 and N=1 has been established, and in [14], they investigated system (1.1) with λ=1, N=1 and an additional perturbation term. By mountain pass theorem and symmetric mountain pass theorem in [15] and a local minimum theorem in [16], they obtained some existence and multiplicity results when F satisfies superquadratic growth at infinity and some other reasonable conditions at origin.

    In this paper, motivated by the idea in [17,18], being different from those in [13,14], we consider the case that F has subquadratic growth only near the origin and no any growth condition at infinity. Our main tools are Ekeland's variational principle in [19], a variant of Clark's theorem in [17] and a cut-off technique in [18]. We obtain that system (1.1) has a ground state weak solution uλ if λ is in some given interval and then some estimates of uλ are given, and when F(t,x) is also even about x near the origin for a.e. t[0,T], for each given λ>0, system (1.1) has infinitely many weak solutions {uλn} with uλn0 as n. Next, we make some assumptions and state our main results.

    (f0) There exist constants M1>0 and 0<p1<2 such that

    F(t,x)M1|x|p1 (1.3)

    for all x¯Bδ0 and a.e. t[0,T].

    (f1) There exist constants M2>0 and 0<p2<p1<2 such that

    F(t,x)M2|x|p2 (1.4)

    for all x¯Bδ0 and a.e. t[0,T].

    (f0) There exist constants M1>0 and 0<p1<1 such that (1.3) holds.

    (f1) There exist constants M2>0 and 0<p2<p1<1 such that (1.4) holds.

    (f2) There exists a constant η(0,2) such that

    (F(t,x),x)ηF(t,x)

    for all x¯Bδ0 and a.e. t[0,T].

    (f3) F(t,x)=F(t,x) for all x¯Bδ0 and a.e. t[0,T].

    Theorem 1.1. Suppose that (F0), (f0), (f1) and (f2) hold. If

    0<λmin{|cos(πα)|2C,(1B)2p2(δ2)2p2|cos(πα)|2C},

    then system (1.1) has a ground state weak solution uλ satisfying

    uλ2p2min{1,(1B)2p2(δ2)2p2},uλ2p2{B2p2,(δ2)2p2}.

    where

    B=T2α12mΓ(α)(2α1)12,C=max{p,η}max{M1,M2}Tmaxt[0,T]eL(t)max{Bp1,Bp2}.

    If (f0) and (f1) are replaced by the stronger conditions (f0) and (f1), then (f2) is not necessary in Theorem 1.1. So we have the following result.

    Theorem 1.2. Suppose that (F0), (f0) and (f1) hold. If

    0<λmin{|cos(πα)|3C,(1B)2p1(δ2)2p1|cos(πα)|3C},

    then system (1.1) has a ground state weak solution uλ satisfying

    uλ2p1min{1,(1B)2p1(δ2)2p1},uλ2p1{B2p1,(δ2)2p1},

    where C=maxt[0,T]eL(t)max{(1+ρ0)a0BT0b(t)dt,M1p1TBp1,ρ0M1TBp1+1}, a0=maxs[0,δ]a(s) and ρ0=maxs[δ2,δ]|ρ(s)| and ρ(s)C1(R,[0,1]) is any given even cut-off function satisfying

    ρ(s)={1,if |s|δ/2,0,if |s|>δ. (1.5)

    Theorem 1.3. Suppose that (F0), (f0), (f1) and (f3) hold. Then for each λ>0, system (1.1) has a sequence of weak solutions {uλn} satisfying {uλn}0, as n.

    Remark 1.1. Theorem 1.1-Theorem 1.3 still hold even if r(t)0 for all t[0,T], that is, the damped-like term disappears, which are different from those in [2,3,4,5,6,7,8,9,10,11,12] because all those assumptions with respect to x in our theorems are made only near origin without any assumption near infinity.

    The paper is organized as follows. In section 2, we give some preliminary facts. In section 3, we prove Theorem 1.1–Theorem 1.3.

    In this section, we introduce some definitions and lemmas in fractional calculus theory. We refer the readers to [1,9,20,21,22]. We also recall Ekeland's variational principle in [19] and the variant of Clark's theorem in [17].

    Definition 2.1. (Left and Right Riemann-Liouville Fractional Integrals [22]) Let f be a function defined on [a,b]. The left and right Riemann-Liouville fractional integrals of order γ for function f denoted by aDγtf(t) and tDγbf(t), respectively, are defined by

    aDγtf(t)=1Γ(γ)ta(ts)γ1f(s)ds,t[a,b],γ>0,tDγbf(t)=1Γ(γ)bt(st)γ1f(s)ds,t[a,b],γ>0.

    provided the right-hand sides are pointwise defined on [a,b], where Γ>0 is the Gamma function.

    Definition 2.2. ([22]) For nN, if γ=n, Definition 2.1 coincides with nth integrals of the form

    aDntf(t)=1(n1)!ta(ts)n1f(s)ds,t[a,b],nN,tDnbf(t)=1(n1)!bt(ts)n1f(s)ds,t[a,b],nN.

    Definition 2.3. (Left and Right Riemann-Liouville Fractional Derivatives [22]) Let f be a function defined on [a,b]. The left and right Riemann-Liouville fractional derivatives of order γ for function f denoted by aDγtf(t) and tDγbf(t), respectively, are defined by

    aDγtf(t)=dndtnaDγntf(t)=1Γ(nγ)dndtn(ta(ts)nγ1f(s)ds),tDγbf(t)=(1)ndndtntDγnbf(t)=(1)nΓ(nγ)dndtn(bt(st)nγ1f(s)ds).

    where t[a,b],n1γ<n and nN. In particular, if 0γ<1, then

    aDγtf(t)=ddtaDγ1tf(t)=1Γ(1γ)ddt(ta(ts)γf(s)ds),t[a,b],tDγbf(t)=ddttDγ1bf(t)=1Γ(1γ)ddt(bt(st)γf(s)ds),t[a,b].

    Remark 2.1. ([9,13]) The left and right Caputo fractional derivatives are defined by the above-mentioned Riemann-Liuville fractional derivative. In particular, they are defined for function belonging to the space of absolutely continuous functions, which we denote by AC([a,b],RN). ACk([a,b],RN)(k=0,1,...) are the space of the function f such that fCk([a,b],RN). In particular, AC([a,b],RN)=AC1([a,b],RN).

    Definition 2.4. (Left and Right Caputo Fractional Derivatives [22]) Let γ0 and nN.

    (ⅰ) If γ(n1,n) and fACn([a,b],RN), then the left and right Caputo fractional derivatives of order γ for function f denoted by caDγtf(t) and ctDγbf(t), respectively, exist almost everywhere on [a,b]. caDγtf(t) and ctDγbf(t) are represented by

    caDγtf(t)=aDγntfn(t)=1Γ(nγ)(ta(ts)nγ1f(n)(s)ds),ctDγbf(t)=(1)n tDγnbfn(t)=(1)nΓ(nγ)(bt(st)nγ1f(n)(s)ds),

    respectively, where t[a,b]. In particular, if 0<γ<1, then

    caDγtf(t)=aDγ1tf(t)=1Γ(1γ)(ta(ts)γf(s)ds),t[a,b],ctDγbf(t)=tDγ1bf(t)=1Γ(1γ)(bt(st)γf(s)ds),t[a,b].

    (ⅱ) If γ=n1 and fACn([a,b],RN), then caDγtf(t) and ctDγbf(t) are represented by

    caDn1tf(t)=f(n1)(t),t[a,b],ctDn1bf(t)=(1)n1f(n1)(t),t[a,b].

    In particular, caD0tf(t)= ctD0bf(t)=f(t), t[a,b].

    Lemma 2.1. ([22]) The left and right Riemann-Liouville fractional integral operators have the property of a semigroup, i.e.

    aDγ1t(aDγ2tf(t))=aDγ1γ2tf(t),tDγ1b(tDγ2bf(t))=tDγ1γ2bf(t),γ1,γ2>0,

    in any point t[a,b] for continuous function f and for almost every point in [a,b] if the function fL1([a,b],RN).

    For 1r<, define

    uLr=(T0|u(t)|rdt)1r (2.1)

    and

    u=maxt[0,T]|u(t)|. (2.2)

    Definition 2.5. ([1]) Let 0<α1 and 1<p<. The fractional derivative space Eα,p0 is defined by closure of C0([0,T],RN) with respect to the norm

    uα,p=(T0|u(t)|pdt+T0|c0Dαtu(t)|pdt)1p. (2.3)

    Remark 2.2. ([9]) Eα,p0 is the space of functions uLp([0,T],RN) having an α-order Caputo fractional derivative c0Dαtu(t)Lp([0,T],RN) and u(0)=u(T)=0.

    Lemma 2.2. ([1]) Let 0<α1 and 1<p<. Eα,p0 is a reflexive and separable Banach space.

    Lemma 2.3. ([1]) Assume that 1<p< and α>1p. Then Eα,p0 compactly embedding in C([0,T],RN).

    Lemma 2.4. ([1]) Let 0<α1 and 1<p<. For all uEα,p0, we have

    uLpTαΓ(α+1)c0DαtuLp. (2.4)

    Moreover, if α>1p and 1p+1q=1, then

    uTα1pΓ(α)((α1)q+1)1qc0DαtuLp. (2.5)

    Definition 2.6. ([13]) Assume that X is a Banach space. An operator A:XX is of type (S)+ if, for any sequence {un} in X, unu and lim supn+A(un),unu0 imply unu.

    Let φ:XR. A sequence {un}X is called (PS) sequence if the sequence {un} satisfies

    φ(un) is bounded, φ(un)0.

    Furthermore, if every (PS) sequence {un} has a convergent subsequence in X, then one call that φ satisfies (PS) condition.

    Lemma 2.5. ([19]) Assume that X is a Banach space and φ:XR is Gˆateaux differentiable, lower semi-continuous and bounded from below. Then there exists a sequence {xn} such that

    φ(xn)infXφ,φ(xn)0.

    Lemma 2.6. ([17]) Let X be a Banach space, φC1(X,R). Assume φ satisfies the (PS) condition, is even and bounded below, and φ(0)=0. If for any kN, there exist a k-dimensional subspace Xk of X and ρk>0 such that supXkSpkφ<0, where Sρ={uX|u=ρ}, then at least one of the following conclusions holds.

    (ⅰ) There exist a sequence of critical points {uk} satisfying φ(uk)<0 for all k and uk0 as k.

    (ⅱ) There exists a constant r>0 such that for any 0<a<r there exists a critical point u such that u=a and φ(u)=0.

    Remark 2.3. ([17]) Lemma 2.6 implies that there exist a sequence of critical points uk0 such that φ(uk)0, φ(uk)0 and uk0 as k.

    Now, we establish the variational functional defined on the space Eα,20 with 12<α1. We follow the same argument as in [13] where the one-dimensional case N=1 and λ=1 for system (1.1) was investigated. For reader's convenience, we also present the details here. Note that L(t):=t0(r(s)/p(s))ds,0<meL(t)p(t)M and q(t)p(t)0 for a.e. t[0,T]. Then system (1.1) is equivalent to the system

    {ddt(eL(t)p(t)(12 0Dξt(u(t))+12 tDξT(u(t))))+eL(t)q(t)u(t)=λeL(t)F(t,u),  a.e. t[0,T],u(0)=u(T)=0. (2.6)

    By Lemma 2.1, for every uAC([0,T],RN), it is easy to see that system (2.6) is equivalent to the system

    {ddt[eL(t)p(t)(12 0Dξ2t(0Dξ2tu(t))+12 tDξ2T(tDξ2Tu(t)))]+eL(t)q(t)u(t)=λeL(t)F(t,u),  a.e. t[0,T],u(0)=u(T)=0, (2.7)

    where ξ[0,1).

    By Definition 2.4, we obtain that uAC([0,T],RN) is a solution of problem (2.7) if and only if u is a solution of the following system

    {ddt(eL(t)p(t)(12 0Dα1t(c0Dαtu(t))12 tDα1T(ctDαTu(t))))+eL(t)q(t)u(t)=λeL(t)F(t,u),u(0)=u(T)=0, (2.8)

    for a.e. t[0,T], where α=1ξ2(12,1]. Hence, the solutions of system (2.8) correspond to the solutions of system (1.1) if uAC([0,T],RN).

    In this paper, we investigate system (2.8) in the Hilbert space Eα,20 with the corresponding norm

    u=(T0eL(t)p(t)(|c0Dαtu(t)|2+|u(t)|2)dt)12.

    It is easy to see that u is equivalent to uα,2 and

    mT0|c0Dαtu(t)|2dtT0eL(t)p(t)|c0Dαtu(t)|2dtMT0|c0Dαtu(t)|2dt.

    So

    uL2TαmΓ(α+1)(T0eL(t)p(t)|c0Dαtu(t)|2dt)12,

    and

    uBu, (2.9)

    where

    B=T2α12mΓ(α)(2α1)12>0.

    (see [13]).

    Lemma 2.7. ([13]) If 12<α1, then for every uEα,20, we have

    |cos(πα)|u2T0eL(t)p(t)(c0Dαtu(t),ctDαTu(t))dt+T0eL(t)p(t)|u(t)|2dtmax{Mm|cos(πα)|,1}u2. (2.10)

    We follow the idea in [17] and [18]. We first modify and extend F to an appropriate ˜F defined by

    ˜F(t,x)=ρ(|x|)F(t,x)+(1ρ(|x|))M1|x|p1,  for all xRN,

    where ρ is defined by (1.5).

    Lemma 3.1. Let (F0), (f0), (f1) (or (f0), (f1)), (f2) and (f3) be satisfied. Then

    (˜F0) ˜F(t,x) is continuously differentiable in xRN for a.e. t[0,T], measurable in t for every xRN, and there exists bL1([0,T];R+) such that

    |˜F(t,x)|a0b(t)+M1|x|p1,|˜F(t,x)|(1+ρ0)a0b(t)+M1p1|x|p11+ρ0M1|x|p1

    for all xRN and a.e. t[0,T];

    (˜f0) ˜F(t,x)M1|x|p1 for all xRN and a.e. t[0,T];

    (˜f1) ˜F(t,x)max{M1,M2}(|x|p1+|x|p2) for all xRN and a.e. t[0,T];

    (˜f2) (˜F(t,x),x)θ˜F(t,x) for all xRN and a.e. t[0,T], where θ=max{p1,η};

    (˜f3) ˜F(t,x)=˜F(t,x) for all xRN and a.e. t[0,T].

    Proof. We only prove (˜f0), (˜f1) and (˜f2). (˜F0) can be proved by a similar argument by (F0). By the definition of ˜F(t,x), (f0) and (f1) (or (f0) and (f1)), we have

    M1|x|p1˜F(t,x)=F(t,x)M2|x|p2, if |x|δ/2,
    ˜F(t,x)=M1|x|p1, if |x|>δ,
    ˜F(t,x)F(t,x)+M1|x|p1M1|x|p1+M2|x|p2, if δ/2<|x|δ

    and

    ˜F(t,x)ρ(|x|)M1|x|p1+(1ρ(|x|))M1|x|p1=M1|x|p1, if δ/2<|x|δ.

    Hence, (˜f1) holds. Note that

    θ˜F(t,x)(˜F(t,x),x)=ρ(|x|)(θF(t,x)(F(t,x),x))+(θp1)(1ρ(|x|))M1|x|p1|x|ρ(|x|)(F(t,x)M1|x|p1).

    It is obvious that the conclusion holds for 0|x|δ/2 and |x|>δ. If δ/2<|x|δ, by using θp1, (f2), (˜f1) and the fact sρ(s)0 for all sR, we can get the conclusion (˜f2). Finally, since ρ(|x|) is even for all xRN, by (f3) and the definition of ˜F(t,x), it is easy to get (˜f3).

    Remark 3.1. From the proof of Lemma 3.1, it is easy to see that (F0), (f0) (or (f0)) and (f1) (or (f1)) independently imply (˜F0), (˜f0) and (˜f1), respectively.

    Consider the modified system

    {ddt(eL(t)p(t)(12 0Dα1t(c0Dαtu(t))12 tDα1T(ctDαTu(t))))+eL(t)q(t)u(t)=λeL(t)˜F(t,u),u(0)=u(T)=0, (3.1)

    for a.e. t[0,T], where α=1ξ2(12,1].

    If the equality

    T0eL(t)[12p(t)((c0Dαtu(t),ctDαTv(t))+(ctDαTu(t),c0Dαtv(t)))+p(t)(u(t),v(t))+(q(t)p(t))(u(t),v(t))λ(˜F(t,u(t)),v(t))]dt=0

    holds for every vEα,20, then we call uEα,20 is a weak solution of system (3.1).

    Define the functional ˜J:Eα,20R by

    ˜J(u)=T0eL(t)[12p(t)((c0Dαtu(t),ctDαTu(t))+|u(t)|2)+12(q(t)p(t))|u(t)|2λ˜F(t,u(t))]dt, for all uEα,20.

    Then (˜F0) and Theorem 6.1 in [9] imply that ˜JC1(Eα,20,R), and for every u,vEα,20, we have

    ˜J(u),v=T0eL(t)[12p(t)((c0Dαtu(t),ctDαTv(t))+(ctDαTu(t),c0Dαtv(t)))+p(t)(u(t),v(t))+(q(t)p(t))(u(t),v(t))λ(˜F(t,u(t)),v(t))]dt.

    Hence, a critical point of ˜J(u) corresponds to a weak solution of problem (3.1).

    Let

    Au,v:=T0eL(t)[12p(t)((c0Dαtu(t),ctDαTv(t))+(ctDαTu(t),c0Dαtv(t)))+p(t)(u(t),v(t))+(q(t)p(t))(u(t),v(t))]dt.

    Lemma 3.2. ([13])

    γ1u2Au,uγ2u2,for all uEα,20, (3.2)

    where γ1=|cos(πα)| and γ2=(max{Mm|cosπα|,1}+maxt[0,T](q(t)p(t))).

    Lemma 3.3. Assume that (F0), (f0) and (f1) (or (f0) and (f1)) hold. Then for each λ>0, ˜J is bounded from below on Eα,20 and satisfies (PS) condition.

    Proof. By (˜f1), (2.9) and (3.2), we have

    ˜J(u)=12Au,uλT0eL(t)˜F(t,u(t))dtγ12u2λmax{M1,M2}T0eL(t)(|u(t)|p1+|u(t)|p2)dtγ12u2λmax{M1,M2}Tmaxt[0,T]eL(t)(up1+up2)γ12u2λmax{M1,M2}Tmaxt[0,T]eL(t)[Bp1up1+Bp2up2].

    It follows from 0<p2<p1<2 that

    ˜J(u)+, as u.

    Hence, ˜J is coercive and then is bounded from below. Now we prove that ˜J satisfies the (PS) condition. Assume that {un} is a (PS) sequence of ˜J, that is,

    ˜J(un) is bounded, ˜J(un)0. (3.3)

    Then by the coercivity of ˜J and (3.3), there exists C0>0 such that unC0 and then by Lemma 2.3, there exists a subsequence (denoted again by {un}) such that

    unu, weakly in Eα,20, (3.4)
    unu, a.e.  in C([0,T],R). (3.5)

    Therefore, the boundness of {un} and (3.3) imply that

    |˜J(un),unu|˜J(un)(Eα,20)unu,˜J(un)(Eα,20)(un+u)0, (3.6)

    where (Eα,20) is the dual space of Eα,20, and (˜F0), (2.9) together with (3.5) imply that

    |λT0(˜F(t,un(t)),un(t)u(t))dt|λT0|˜F(t,un(t))||(un(t)u(t))|dtλunuT0[(1+ρ0)a0b(t)+M1p1|un(t)|p11+ρ0M1|un(t)|p1]dtλunu[(1+ρ0)a0T0b(t)dt+M1p1TBp11Cp110+M1Tρ0Bp1Cp10]0. (3.7)

    Note that

    ˜J(un),unu=Aun,unuλT0(˜F(t,un(t)),un(t)u(t))dt.

    Then (3.6) and (3.7) imply that limnAun,unu=0. Moreover, by (3.4), we have

    limnAu,unu=0.

    Therefore

    limnAunAu,unu=0.

    Since A is of type (S)+ (see [13]), by Definition 2.6, we obtain unu in Eα,20.

    Define a Nehari manifold by

    Nλ={uEα,20/{0}|˜Jλ(u),u=0}.

    Lemma 3.4. Assume that (F0) and (f0) (or (f0)) hold. For each λ>0, ˜Jλ has a nontrivial least energy (ground state) weak solution uλ, that is, uλNλ and ˜Jλ(uλ)=infNλ˜Jλ. Moreover, the least energy can be estimated as follows

    ˜Jλ(uλ)Gλ:=(p1/γ2)p12p1[λM1mint[0,T]eL(t)T0|w0(t)|p1dt]22p1(p12)2.

    where w0=ww, and w=(TπsinπtT,0,,0)Eα,20.

    Proof. By Lemma 3.3 and ˜JC1(Eα,20,R), for each λ>0, Lemma 2.5 implies that there exists some uλEα,20 such that

    ˜J(uλ)=infvEα,20˜J(v)and ˜J(uλ)=0. (3.8)

    By (3.2) and (˜f0), we have

    ˜Jλ(sw0)=12A(sw0),sw0λT0eL(t)˜F(t,sw0(t))dtγ22s2w02λT0eL(t)M1|sw0(t)|p1dtγ22s2λM1mint[0,T]eL(t)sp1T0|w0(t)|p1dt. (3.9)

    for all s[0,). Define g:[0,+)R by

    g(s)=γ22s2λM1mint[0,T]eL(t)sp1T0|w0(t)|p1dt.

    Then g(s) achieves its minimum at

    s0,λ=(p1λM1mint[0,T]eL(t)T0|w0(t)|p1dtγ2)12p1

    and

    g(s0,λ)=(p1/γ2)p12p1[λM1mint[0,T]eL(t)T0|w0(t)|p1dt]22p1(p12)2.

    Note that p1<2. So g(s0,λ)<0. Hence, (3.9) implies that

    ˜Jλ(uλ)=infvEα,20˜Jλ(v)˜Jλ(s0,λw0)g(s0,λ)<0=˜Jλ(0)

    and then uλ0 which together with (3.8) implies that uλNλ and ˜Jλ(uλ)=infNλ˜Jλ.

    Lemma 3.5. Assume that (F0), (f1) and (f2) hold. If 0<λ|cos(πα)|2C, then the following estimates hold

    uλ2p22λC|cos(πα)|,uλ2p22λCB2p2|cos(πα)|.

    Proof. It follows from Lemma 3.1, (2.9) and ˜J(uλ),uλ=0 that

    T0eL(t)[p(t)(c0Dαtuλ(t),ctDαTuλ(t))+p(t)(uλ(t),uλ(t))+(q(t)p(t))(uλ(t),uλ(t))]dt=λT0eL(t)(˜F(t,uλ(t)),uλ(t))dtλθT0eL(t)˜F(t,uλ(t))dtλθmax{M1,M2}maxt[0,T]eL(t)T0(|uλ(t)|p1+|uλ(t)|p2)dtλθmax{M1,M2}Tmaxt[0,T]eL(t)(uλp1+uλp2)λθmax{M1,M2}Tmaxt[0,T]eL(t)[Bp1uλp1+Bp2uλp2]λC(uλp1+uλp2). (3.10)

    We claim that uλ1 uniformly for all 0<λ|cos(πα)|2C. Otherwise, we have a sequence of {λn|cos(πα)|2C} such that uλn>1. Thus uλnp2<uλnp1 since p2<p1<2. By (2.10) and (3.10), we obtain

    T0eL(t)[p(t)(c0Dαtuλn(t),ctDαTuλn(t))+p(t)(uλn(t),uλn(t))+(q(t)p(t))(uλn(t),uλn(t))]dt|cos(πα)|uλn2+T0eL(t)(q(t)p(t))|uλn(t)|2dt. (3.11)

    By (3.10) and (3.11), we obtain

    |cos(πα)|uλn2+T0eL(t)(q(t)p(t))|uλn(t)|2dtλnC(uλnp1+uλnp2).

    Since q(t)p(t)>0,

    |cos(πα)|uλn2λnC(uλnp1+uλnp2)2λnCuλnp1.

    Then

    uλn2p12λnC|cos(πα)|1,

    which contradicts with the assumption uλn>1. Now, from (3.10) we can get

    |cos(πα)|uλ2λC(uλp1+uλp2)2λCuλp2.

    So

    uλ2p22λC|cos(πα)|.

    By (2.9), we can obtain

    uλBuλB(2λC|cos(πα)|)12p2.

    Observe that, in the proof of Lemma 3.5, (˜f2) is used only in (3.10). If we directly use (˜F0) to rescale (˜F(t,uλ(t)),uλ(t)) in (3.10). Then the assumption (f2) is not necessary but we have to pay the price that p(0,1). To be precise, we have the following lemma.

    Lemma 3.6. Assume that (F0) and (f0) hold. If 0<λ|cos(πα)|3C, then the following estimates hold

    uλ2p13λC|cos(πα)|,uλ2p13λCB2p1|cos(πα)|.

    Proof. It follows from (F0), Lemma 3.1, Remark 3.1, (2.9) and ˜J(uλ),uλ=0 that

    T0eL(t)[p(t)(c0Dαtuλ(t),ctDαTuλ(t))+p(t)(uλ(t),uλ(t))+(q(t)p(t))(uλ(t),uλ(t))]dt=λT0eL(t)(˜F(t,uλ(t)),uλ(t))dtλmaxt[0,T]eL(t)T0|˜F(t,uλ(t))||uλ(t)|dtλmaxt[0,T]eL(t)T0[(1+ρ0)a0b(t)|uλ(t)|+M1p1|uλ(t)|p1+ρ0M1|uλ(t)|p1+1]dtλmaxt[0,T]eL(t)[(1+ρ0)a0uλT0b(t)dt+M1p1uλp1+ρ0M1Tuλp1+1]λmaxt[0,T]eL(t)[(1+ρ0)a0BuλT0b(t)dt+M1p1TBp1uλp1+ρ0M1TBp1+1uλp1+1]λC(uλ+uλp1+uλp1+1). (3.12)

    We claim that uλ1 uniformly for all 0<λ|cos(πα)|3C. Otherwise, we have a sequence of {λn|cos(πα)|3C} such that uλn>1. Thus uλnp1<uλn<uλnp1+1 since p1<1. By (2.10) and (3.12), we obtain

    T0eL(t)[p(t)(c0Dαtuλn(t),ctDαTuλn(t))+p(t)(uλn(t),uλn(t))+(q(t)p(t))(uλn(t),uλn(t))]dt|cos(πα)|uλn2+T0eL(t)(q(t)p(t))|uλn(t)|2dt. (3.13)

    By (3.12) and (3.13), we obtain

    |cos(πα)|uλn2+T0eL(t)(q(t)p(t))|uλn|2dtλnC(uλ+uλp1+uλp1+1).

    Since q(t)p(t)>0,

    |cos(πα)|uλn2λnC(uλ+uλp1+uλp1+1)3λnCuλnp1+1.

    Then

    uλn1p13λnC|cos(πα)|1,

    which contradicts with the assumption uλn>1. Now, we can get from (3.12) that

    |cos(πα)|uλ2λC(uλ+uλp1+uλp1+1)3λCuλp1.

    So

    uλ2p13λC|cos(πα)|.

    By (2.9), we can obtain

    uλBuλB(3λC|cos(πα)|)12p1.

    Proof of Theorem 1.1. Since 0<λmin{|cos(πα)|2C,(1B)2p2(δ2)2p2|cos(πα)|2C}, Lemma 3.5 implies that

    uλδ2.

    Therefore, for all 0<λmin{|cos(πα)|2C,(1B)2p2(δ2)2p2|cos(πα)|2C}, we have ˜F(t,u(t))=F(t,u(t)) and then uλ is a nontrivial weak solution of the original problem (1.1). Moreover, Lemma 3.5 implies that limλ0uλ=0 as λ0 and

    uλ2p2min{1,(1B)2p2(δ2)2p2},uλ2p2B2p2{1,(1B)2p2(δ2)2p2}.

    Proof of Theorem 1.2. Note that 0<λmin{|cos(πα)|3C,(1B)2p1(δ2)2p1|cos(πα)|3C}. Similar to the proof of Theorem 1.1, by Lemma 3.6, it is easy to complete the proof.

    Proof of Theorem 1.3. By Lemma 3.1 and Lemma 3.3, we obtain that ˜J satisfies (PS) condition and is even and bounded from below, and ˜J(0)=0. Next, we prove that for any kN, there exists a subspace k-dimensional subspace XkEα,20 and ρk>0 such that

    supuXkSρk˜Jλ(u)<0.

    In fact, for any kN, assume that Xk is any subspace with dimension k in Eα,20. Then by (2.10) and Lemma 3.1, there exist constants C1,C2>0 such that

    ˜J(u)max{Mm|cos(πα)|,1}u2+12T0eL(t)(q(t)p(t))|u(t)|2dtλT0eL(t)˜F(t,u(t))dtmax{Mm|cos(πα)|,1}u2+C12u2λC2T0˜F(t,u(t))dtmax{Mm|cos(πα)|,1}u2+C1B22u2λC2M1T0|u(t)|p1dt[max{Mm|cos(πα)|,1}+C1B22]u2λC2M1up1Lp1.

    Since all norms on Xk are equivalent and p1<2, for each fixed λ>0, we can choose ρk>0 small enough such that

    supuXkSρk˜Jλ(u)<0.

    Thus, by Lemma 2.6 and Remark 2.3. ˜Jλ has a sequence of nonzero critical points {uλn}Eα,20 converging to 0 and ˜Jλ(uλn)0. Hence, for each fixed λ>0, (3.1) has a sequence of weak solutions {uλn}Eα,20 with uλn0, as n. Furthermore, there exists n0 large enough such that uλnδ2B for all nn0 and then (2.9) implies that uλnδ2 for all nn0. Thus, ˜F(t,u(t))=F(t,u(t)) and then {uλn}n0 is a sequence of weak solutions of the original problem (1.1) for each fixed λ>0.

    When the nonlinear term F(t,x) is local subquadratic only near the origin with respect to x, system (1.1) with λ in some given interval has a ground state weak solution uλ. If the nonlinear term F(t,x) is also locally even near the origin with respect to x, system (1.1) with λ>0 has infinitely many weak solutions {uλn}.

    This project is supported by Yunnan Ten Thousand Talents Plan Young & Elite Talents Project and Candidate Talents Training Fund of Yunnan Province (No: 2017HB016).

    The authors declare that they have no conflicts of interest.



    [1] G. Kirchhoff, Mechanik, Teubner, Leipzig, 1883.
    [2] D. Cassani, Z. Liu, C. Tarsi, J. Zhang, Multiplicity of sign-changing solutions for Kirchhoff-type equations, Nonlinear Anal., 186 (2019), 145–161. https://doi.org/10.1016/j.na.2019.01.025 doi: 10.1016/j.na.2019.01.025
    [3] B. Cheng, X. H. Tang, Ground state sign-changing solutions for asymptotically 3-linear Kirchhoff-type problems, Complex Var. Elliptic Equations, 62 (2017), 1093–1116. https://doi.org/10.1080/17476933.2016.1270272 doi: 10.1080/17476933.2016.1270272
    [4] Y. B. Deng, S. J. Peng, W. Shuai, Existence and asymptotic behavior of nodal solutions for the Kirchhoff-type problems in R3, J. Funct. Anal., 269 (2015), 3500–3527. https://doi.org/10.1016/j.jfa.2015.09.012 doi: 10.1016/j.jfa.2015.09.012
    [5] G. M. Figueiredo, J. R. Santos Júnior, Existence of a least energy nodal solution for a Schrödinger-Kirchhoff equation with potential vanishing at infinity, J. Math. Phys., 56 (2015), 051506. https://doi.org/10.1063/1.4921639 doi: 10.1063/1.4921639
    [6] F. Y. Li, C. Gao, X. Zhu, Existence and concentration of sign-changing solutions to Kirchhoff-type system with Hartree-type nonlinearity, J. Math. Anal. Appl., 448 (2017), 60–80. https://doi.org/10.1016/j.jmaa.2016.10.069 doi: 10.1016/j.jmaa.2016.10.069
    [7] Q. Li, X. Du, Z. Zhao, Existence of sign-changing solutions for nonlocal Kirchhoff-Schrödinger-type equations in R3, J. Math. Anal. Appl., 477 (2019), 174–186. https://doi.org/10.1016/j.jmaa.2019.04.025 doi: 10.1016/j.jmaa.2019.04.025
    [8] S. Lu, Signed and sign-changing solutions for a Kirchhoff-type equation in bounded domains, J. Math. Anal. Appl., 432 (2015), 965–982. https://doi.org/10.1016/j.jmaa.2015.07.033 doi: 10.1016/j.jmaa.2015.07.033
    [9] A. M. Mao, S. X. Luan, Sign-changing solutions of a class of nonlocal quasilinear elliptic boundary value problems, J. Math. Anal. Appl., 383 (2011), 239–243. https://doi.org/10.1016/j.jmaa.2011.05.021 doi: 10.1016/j.jmaa.2011.05.021
    [10] A. M. Mao, Z. T. Zhang, Sign-changing and multiple solutions of Kirchhoff type problems without the P.S. condition, Nonlinear Anal., 70 (2009), 1275–1287. https://doi.org/10.1016/j.na.2008.02.011 doi: 10.1016/j.na.2008.02.011
    [11] M. Shao, A. Mao, Signed and sign-changing solutions of Kirchhoff type problems, J. Fixed. Point Theory Appl., 20 (2018), 1–20. https://doi.org/10.1007/s11784-018-0486-9 doi: 10.1007/s11784-018-0486-9
    [12] W. Shuai, Sign-changing solutions for a class of Kirchhoff-type problem in bounded domains, J. Differ. Equations, 259 (2015), 1256–1274. https://doi.org/10.1016/j.jde.2015.02.040 doi: 10.1016/j.jde.2015.02.040
    [13] J. Sun, L. Li, M. Cencelj, B. Gabrovšek, Infinitely many sign-changing solutions for Kirchhoff type problems in R3, Nonlinear Anal., 186 (2019), 33–54. https://doi.org/10.1016/j.na.2018.10.007 doi: 10.1016/j.na.2018.10.007
    [14] X. H. Tang, B. Cheng, Ground state sign-changing solutions for Kirchhoff type problems in bounded domains, J. Differ. Equations, 261 (2016), 2384–2402. https://doi.org/10.1016/j.jde.2016.04.032 doi: 10.1016/j.jde.2016.04.032
    [15] D. B. Wang, Least energy sign-changing solutions of Kirchhoff-type equation with critical growth, J. Math. Phys., 61 (2020), 011501. https://doi.org/10.1063/1.5074163 doi: 10.1063/1.5074163
    [16] L. Wang, B. L. Zhang, K. Cheng, Ground state sign-changing solutions for the Schrödinger-Kirchhoff equation in R3, J. Math. Anal. Appl., 466 (2018), 1545–1569. https://doi.org/10.1016/j.jmaa.2018.06.071 doi: 10.1016/j.jmaa.2018.06.071
    [17] H. Ye, The existence of least energy nodal solutions for some class of Kirchhoff equations and Choquard equations in RN, J. Math. Anal. Appl., 431 (2015), 935–954. https://doi.org/10.1016/j.jmaa.2015.06.012 doi: 10.1016/j.jmaa.2015.06.012
    [18] Z. T. Zhang, K. Perera, Sign changing solutions of Kirchhoff type problems via invariant sets of descentow, J. Math. Anal. Appl., 317 (2006), 456–463. https://doi.org/10.1016/j.jmaa.2005.06.102 doi: 10.1016/j.jmaa.2005.06.102
    [19] X. Zhong, C. L. Tang, The existence and nonexistence results of ground state nodal solutions for a Kirchhoff type problem, Commun. Pure Appl. Anal., 16 (2017), 611–628. https://doi.org/10.3934/cpaa.2017030 doi: 10.3934/cpaa.2017030
    [20] C. O. Alves, M. A. S. Souto, Existence of least energy nodal solution for a Schrödinger-Poisson system in bounded domains, Z. Angew. Math. Phys., 65 (2014), 1153–1166. https://doi.org/10.1007/s00033-013-0376-3 doi: 10.1007/s00033-013-0376-3
    [21] S. Chen, X. Tang, Ground state sign-changing solutions for a class of Schrödinger-Poisson type problems in R3, Z. Angew. Math. Phys., 67 (2016). https://doi.org/10.1007/s00033-016-0695-2 doi: 10.1007/s00033-016-0695-2
    [22] S. Khoutir, Least energy sign-changing solutions for a class of Schrödinger-Poisson system on bounded domains, J. Math. Phys., 62 (2021), 031509. https://doi.org/10.1063/5.0040741 doi: 10.1063/5.0040741
    [23] S. Kim, J. Seok, On nodal solutions of the nonlinear Schrödinger-Poisson equations, Commun. Contemp. Math., 14 (2012), 1250041. https://doi.org/10.1142/S0219199712500411 doi: 10.1142/S0219199712500411
    [24] Z. Liang, J. Xu, X. Zhu, Revisit to sign-changing solutions for the nonlinear Schrödinger-Poisson system in R3, J. Math. Anal. Appl., 435 (2016), 783–799. https://doi.org/10.1016/j.jmaa.2015.10.076 doi: 10.1016/j.jmaa.2015.10.076
    [25] Z. Liu, Z. Wang, J. Zhang, Infinitely many sign-changing solutions for the nonlinear Schrödinger-Poisson system, Ann. Mat. Pura Appl., 4 (2016), 775–794. https://doi.org/10.1007/s10231-015-0489-8 doi: 10.1007/s10231-015-0489-8
    [26] W. Shuai, Q. Wang, Existence and asymptotic behavior of sign-changing solutions for the nonlinear Schrödinger-Poisson system in R3, Z. Angew. Math. Phys., 66 (2015), 3267–3282. https://doi.org/10.1007/s00033-015-0571-5 doi: 10.1007/s00033-015-0571-5
    [27] D. B. Wang, H. Zhang, W. Guan, Existence of least-energy sign-changing solutions for Schrödinger-Poisson system with critical growth, J. Math. Anal. Appl., 479 (2019), 2284–2301. https://doi.org/10.1016/j.jmaa.2019.07.052 doi: 10.1016/j.jmaa.2019.07.052
    [28] Z. Wang, H. Zhou, Sign-changing solutions for the nonlinear Schrödinger-Poisson system in R3, Calc. Var. Partial Differ. Equations, 52 (2015), 927–943. https://doi.org/10.1007/s00526-014-0738-5 doi: 10.1007/s00526-014-0738-5
    [29] X. Zhong, C. L. Tang, Ground state sign-changing solutions for a Schrödinger-Poisson system with a critical nonlinearity in R3, Nonlinear Anal. Real World Appl., 39 (2018), 166–184. https://doi.org/10.1016/j.nonrwa.2017.06.014 doi: 10.1016/j.nonrwa.2017.06.014
    [30] V. Bobkov, Least energy nodal solutions for elliptic equations with indefinite nonlinearity, Electron. J. Qual. Theory Differ. Equations, 2014 (2014), 1–15. https://doi.org/10.14232/ejqtde.2014.1.56 doi: 10.14232/ejqtde.2014.1.56
    [31] D. Ruiz, The Schrödinger-Poisson equation under the effect of a nonlinear local term, J. Funct. Anal., 237 (2006), 655–674. https://doi.org/10.1016/j.jfa.2006.04.005 doi: 10.1016/j.jfa.2006.04.005
    [32] M. Willem, Minimax Theorems, Birkhäuser, Bosten, 1996.
  • This article has been cited by:

    1. Ji Hyun Yun, Eunsol Lee, Jeong Hwa Choi, Hyun Kyun Ki, Jiho Park, Seroprevalence of Varicella-Zoster Virus and Measles among Healthcare Workers in a Tertiary Medical Center in Korea, 2022, 10, 2076-393X, 1956, 10.3390/vaccines10111956
    2. Ayako Suzuki, Hiroshi Nishiura, Seasonal transmission dynamics of varicella in Japan: The role of temperature and school holidays, 2022, 20, 1551-0018, 4069, 10.3934/mbe.2023190
    3. Lingling Wang, Xu Yang, Xuwen Wang, Ping Shi, Xuhui Zhang, Assessing vaccine effectiveness for varicella in Wuxi, China: a time-series analysis, 2023, 55, 0785-3890, 10.1080/07853890.2023.2246369
    4. Cheng-Cheng Zhu, Jiang Zhu, Jie Shao, Epidemiological Investigation: Important Measures for the Prevention and Control of COVID-19 Epidemic in China, 2023, 11, 2227-7390, 3027, 10.3390/math11133027
    5. S. Hariharan, L. Shangerganesh, A. Debbouche, V. Antonov, Stability analysis of spatiotemporal reaction–diffusion mathematical model incorporating the varicella virus transmission, 2023, 138, 2190-5444, 10.1140/epjp/s13360-023-04777-6
    6. Ulrikka Nygaard, Mette Holm, Helena Rabie, Maren Rytter, The pattern of childhood infections during and after the COVID-19 pandemic, 2024, 8, 23524642, 910, 10.1016/S2352-4642(24)00236-0
    7. Ana Clara Poncheke Ratto, Lara De Camargo Stein, Letícia Dardes Barbosa , Marina Fernandes Xavier , Milena Roberta Freire da Silva, The Impacts of Social Isolation and Vaccination on cases of Varicella in Brazil: An Ecological Study, 2025, 9, 2675-1542, 10.53843/bms.v9i13.760
    8. Chunmei Duan, Yan Zhang, Qian Zhang, Shuzhi Zhang, Peisong Zhong, Gang Gong, Yefan Zhu, Jie Fei, Jingjing Zhao, Yanling Sun, Yaqing Jin, Yunjie Ren, Yihan Lu, Ye Yao, Hongjie Yu, Impact of natural and socio-economic factors on varicella incidence in children in Shanghai, 2013-2022, 2025, 13, 2296-2565, 10.3389/fpubh.2025.1565717
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

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

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1731) PDF downloads(81) Cited by(0)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog