Research on the degenerate versions of special polynomials provides a new area, introducing the λ-analogue of special polynomials and numbers, such as λ-Sheffer polynomials. In this paper, we propose a new variant of type 2 Bernoulli polynomials and numbers by modifying a generating function. Then we derive explicit expressions and their representations that provide connections among existing λ-Sheffer polynomials. Also, we provide the explicit representations of the proposed polynomials in terms of the degenerate Lah-Bell polynomials and the higher-order degenerate derangement polynomials to confirm the presented identities.
Citation: Jongkyum Kwon, Patcharee Wongsason, Yunjae Kim, Dojin Kim. Representations of modified type 2 degenerate poly-Bernoulli polynomials[J]. AIMS Mathematics, 2022, 7(6): 11443-11463. doi: 10.3934/math.2022638
[1] | Saak Gabriyelyan . Limited type subsets of locally convex spaces. AIMS Mathematics, 2024, 9(11): 31414-31443. doi: 10.3934/math.20241513 |
[2] | Mahmut Karakuş . Spaces of multiplier σ-convergent vector valued sequences and uniform σ-summability. AIMS Mathematics, 2025, 10(3): 5095-5109. doi: 10.3934/math.2025233 |
[3] | Xinwei Su, Shuqin Zhang, Lixin Zhang . Periodic boundary value problem involving sequential fractional derivatives in Banach space. AIMS Mathematics, 2020, 5(6): 7510-7530. doi: 10.3934/math.2020481 |
[4] | Imran Ali, Faizan Ahmad Khan, Haider Abbas Rizvi, Rais Ahmad, Arvind Kumar Rajpoot . Second order evolutionary partial differential variational-like inequalities. AIMS Mathematics, 2022, 7(9): 16832-16850. doi: 10.3934/math.2022924 |
[5] | Heng Yang, Jiang Zhou . Compactness of commutators of fractional integral operators on ball Banach function spaces. AIMS Mathematics, 2024, 9(2): 3126-3149. doi: 10.3934/math.2024152 |
[6] | Narjes Alabkary . Hyper-instability of Banach algebras. AIMS Mathematics, 2024, 9(6): 14012-14025. doi: 10.3934/math.2024681 |
[7] | Hasanen A. Hammad, Hassen Aydi, Maryam G. Alshehri . Solving hybrid functional-fractional equations originating in biological population dynamics with an effect on infectious diseases. AIMS Mathematics, 2024, 9(6): 14574-14593. doi: 10.3934/math.2024709 |
[8] | Muneerah Al Nuwairan, Ahmed Gamal Ibrahim . Nonlocal impulsive differential equations and inclusions involving Atangana-Baleanu fractional derivative in infinite dimensional spaces. AIMS Mathematics, 2023, 8(5): 11752-11780. doi: 10.3934/math.2023595 |
[9] | Sara Smail, Chafika Belabbaci . A characterization of Wolf and Schechter essential pseudospectra. AIMS Mathematics, 2024, 9(7): 17146-17153. doi: 10.3934/math.2024832 |
[10] | Ateq Alsaadi, Manochehr Kazemi, Mohamed M. A. Metwali . On generalization of Petryshyn's fixed point theorem and its application to the product of n-nonlinear integral equations. AIMS Mathematics, 2023, 8(12): 30562-30573. doi: 10.3934/math.20231562 |
Research on the degenerate versions of special polynomials provides a new area, introducing the λ-analogue of special polynomials and numbers, such as λ-Sheffer polynomials. In this paper, we propose a new variant of type 2 Bernoulli polynomials and numbers by modifying a generating function. Then we derive explicit expressions and their representations that provide connections among existing λ-Sheffer polynomials. Also, we provide the explicit representations of the proposed polynomials in terms of the degenerate Lah-Bell polynomials and the higher-order degenerate derangement polynomials to confirm the presented identities.
In this paper, we are interested in establishing the existence and nonexistence results of nontrivial solutions for the coupled fractional Schrödinger systems of Choquard type
{(−Δ)su+λ1u=(Iα∗|u|p)|u|p−2u+βvin RN,(−Δ)sv+λ2v=(Iα∗|v|p)|v|p−2v+βuin RN, | (1.1) |
where s∈(0,1), N≥3, α∈(0,N), p>1, λi>0 are constants for i=1, 2, β>0 is a parameter, and Iα(x) is the Riesz Potential defined as
Iα(x)=Γ(N−α2)Γ(α2)πN22α|x|N−α,x∈RN∖{0}, |
where Γ is the Gamma function.
Here, the nonlocal Laplacian operator (−Δ)s with s∈(0,1) of a function u:RN→R is expressed by the formula
(−Δ)su(x)=C(N,s)P.V.∫RNu(x)−u(z)|x−z|N+2sdz, |
where P.V. stand for the Cauchy principal value on the integral, and C(N,s) is some positive normalization constant (see [1] for details).
It can also be defined as a pseudo-differential operator
F((−Δ)sf)(ξ)=|ξ|2sF(f)(ξ)=|ξ|2sˆf(ξ), |
where F is the Fourier transform.
The problem (1.1) presents nonlocal characteristics in the nonlinearity as well as in the (fractional) diffusion because of the appearance of the terms (Iα∗|u|p)|u|p−2u and (Iα∗|v|p)|v|p−2v. This phenomenon raises some mathematical puzzles that make the study of such problems particularly interesting. We point out that when s=1, λ1=1, p=2, N=3, α=2 and β=0, (1.1) reduces to the Choquard-Pekar equation
−Δu+u=(I2∗|u|2)u,in R3, | (1.2) |
which appeared in 1954 by Pekar [2] describing a polaron at rest in the quantum theory. In 1976, Choquard [3] used this equation to model an electron trapped in its own hole and considered it as an approximation to Hartree-Fock theory of one-component plasma. Subsequently, in 1996 Penrose [4] investigated it as a model for the self-gravitating collapse of a quantum mechanical wave function; see also [5]. The first investigations for existence and uniqueness of ground state solutions of (1.2) go back to the work of Lieb [6]. Lions [7] generalized the result in [6] and proved the existence and multiplicity of positive solutions of (1.2). In addition, the existence and qualitative results of solutions of power type nonlinearities |u|p−2u and for more generic values of α∈(0,N) are discussed by variational method, where N≥3, see [8,9,10,11,12]. Under almost necessary conditions on the nonlinearity F in the spirit of H. Berestycki and P. L. Lions [13], Moroz and Schaftingen [14] considered the existence of a ground state solution u∈H1(RN) to the nonlinear Choquard equation
−Δu+u=(Iα∗F(u))F′(u),in RN. |
When s∈(0,1), Laskin [15] introduced the fractional power of the Laplace operator in (1.1) as an extension of the classical local Laplace operator in the study of nonlinear Schrödinger equations, replacing the path integral over Brownian motions with Lévy flights [16]. This operator has concrete applications in a wide range of fields, see [1,17] and the references therein. Equations involving the fractional Laplacian together with local nonlinearities and the system of weakly coupled equations has been investigated extensively in recent years, and some research results can be found in [18,19,20,21].
When β=0, the system (1.1) can be reduced to two single Choquard equations
(−Δ)su+λ1u=(Iα∗|u|p)|u|p−2uin RN | (1.3) |
and
(−Δ)sv+λ2v=(Iα∗|v|p)|v|p−2vin RN. | (1.4) |
Equations (1.3) and (1.4) arise from the search for standing wave solutions of the following time-dependent fractional Choquard equation:
i∂Ψ∂t=(−Δ)sΨ+λΨ−(Iα∗|Ψ|p)|Ψ|p−2Ψ,(t,x)∈R+×RN, |
where i denotes the imaginary unit.
In [22], by minimizing
S(u)=‖(−Δ)s2u‖22+λ1‖u‖22(∫RN(Iα∗|u|p)|u|p)1p |
on Hs(RN)∖{0}, the authors obtained the existence of ground state solution of (1.3) with p∈(1+αN,N+αN−2s) (see [22, Theorem 4.2]).
Of course, scalar problems can be extended to systems. It is easy to see that the system (1.1) can be regarded as a counterpart of the following systems with standard Laplace operator
{−Δu+u=(Iα∗|u|p)|u|p−2u+λvin RN,−Δv+v=(Iα∗|v|p)|v|p−2v+λuin RN. |
In [23], Chen and Liu studied the systems of Choquard type, when p∈(1+αN,N+αN−2), they obtained the existence of ground state solutions of the systems. Yang et al. [24] considered the corresponding critical case.
Motivated by the above mentioned works, in this paper, we aim to study the existence of positive ground state solutions of the systems (1.1). This class of systems has two new characteristics: One is the presence of the fractional Laplace and the Choquard type functions which are nonlocal, the other is its lack of compactness inherent to problems defined on unbounded domains. In order to overcome such difficulties, next we introduce a special space where we are able to recover some compactness.
First we use ‖⋅‖p denote the norm of Lp(RN) for any 1≤p<∞. The Hilbert space Hs(RN) is defined by
Hs(RN):={u∈L2(RN):∫RN∫RN|u(x)−u(z)|2|x−z|N+2sdxdz<+∞} |
with the scalar product and norm given by
⟨u,v⟩:=∫RN(−Δ)s2u(−Δ)s2vdx+∫RNuvdx,‖u‖:=(‖(−Δ)s2u‖22+‖u‖22)12, |
where
‖(−Δ)s2u‖22:=C(N,s)2∫RN∫RN|u(x)−u(z)|2|x−z|N+2sdxdz. |
The radial space Hsr(RN) of Hs(RN) is defined as
Hsr(RN):={u∈Hs(RN)|u(x)=u(|x|)} |
with the Hs(RN) norm.
Let
‖u‖2λi:=‖(−Δ)s2u‖22+λi‖u‖22,i=1,2 |
for convenience. It is easy to obtain that ‖⋅‖λi and ‖⋅‖ are equivalent norms in Hs(RN). Denote H:=Hs(RN)×Hs(RN) and Hr:=Hsr(RN)×Hsr(RN). The norm of H is given by
‖(u,v)‖2H=‖u‖2λ1+‖v‖2λ2,for all (u,v)∈H. |
The energy functional Eβ associated to (1.1) is
Eβ(u,v)=12∫RN[|(−Δ)s2u|2+|(−Δ)s2v|2+λ1|u|2+λ2|v|2]dx−12p∫RN(Iα∗|u|p)|u|pdx−12p∫RN(Iα∗|v|p)|v|pdx−β∫RNuvdx,for all (u,v)∈H. | (1.5) |
It is easy to obtain that Eβ∈C1(H,R) and
⟨E′β(u,v),(φ,ψ)⟩=∫RN[(−Δ)s2u(−Δ)s2φ+(−Δ)s2v(−Δ)s2ψ+λ1uφ+λ2vψ]dx−∫RN(Iα∗|u|p)|u|p−2uφdx−∫RN(Iα∗|v|p)|v|p−2vψdx−β∫RN(vφ+uψ)dx | (1.6) |
for all (φ,ψ)∈H.
(u,v) is called a nontrivial solution of (1.1) if uβ≢0, vβ≢0 and (u,v)∈H solves (1.1). A positive ground state solution (u,v) of (1.1) is a nontrivial solution of (1.1) such that u>0, v>0 which has minimal energy among all nontrivial solutions. In order to find positive ground state solutions of (1.1), we need to investigate the existence of the minimum value of Eβ, defined in (1.5) under the Nehari manifold constraint
Nβ={(u,v)∈H∖{(0,0)}:⟨E′β(u,v),(u,v)⟩=0}. | (1.7) |
Define
mβ=inf{Eβ(u,v):(u,v)∈Nβ}. |
Furthermore, define E0,i:Hs(RN)→R by
E0,i(u)=12∫RN|(−Δ)s2u|2dx+λi2∫RNu2dx−12p∫RN(Iα∗|u|p)|u|pdx, i=1,2. | (1.8) |
We introduce the Nehari manifolds
N0,i:={u∈Hs(RN)∖{0}:‖(−Δ)s2u‖22+λi‖u‖22−∫RN(Iα∗|u|p)|u|pdx=0}, i=1,2. | (1.9) |
A ground state solution of (1.3) (or (1.4)) is a solution with minimal energy E0,1 (or E0,2) and can be characterized as
minu∈N0,1E0,1(u) (orminu∈N0,2E0,2(u)). |
The main results of our paper are the following.
Theorem 1.1. Suppose s∈(0,1), N≥3, α∈(0,N) and p∈(1+αN,α+NN−2s), then the system (1.1) possesses a positive radial ground state solution (uβ,vβ)∈Nβ with Eβ(uβ,vβ)=mβ>0 for any 0<β<√λ1λ2. Moreover, (uβ,vβ)→(u0,v0) in H as β→0+, where (u0,v0) is a positive radial ground state solution for the system (1.1) with β=0, namely, u0 and v0 are positive radial ground state solutions to problems (1.3) and (1.4), respectively.
Remark 1.1. In comparison with [19], this paper has several new features. Firstly, the system (1.1) contains the Choquard type terms which are more difficult to deal with. Secondly, Lemma 3.11 in [19] shows that (uβ,vβ)→(u0,v0) in H as β→0+, where either v0≡0 and u0 is a ground state solution to one single equation, or u0≡0 and v0 is a ground state solution to the other single equation. While we prove that (u0,v0) is a positive radial ground state solution for the system (1.1) with β=0. Finally, the difference in asymptotic behavior is that it is obtained in this paper that u0>0 and v0>0 are positive radial ground state solutions to problems (1.3) and (1.4), respectively (see Theorem 1.3 in [19]).
Finally, by using the Pohožaev identity (4.1) of the system (1.1), we have the following non-existence result.
Theorem 1.2. Suppose p≥α+NN−2s or p≤1+αN, then the system (1.1) does not admit non-trivial solutions.
Remark 1.2. According to Theorem 1.2, we can know that the range of p∈(1+αN,α+NN−2s) is optimal for the existence of nontrivial solutions to the system (1.1).
The rest of this paper is as following. In Section 2, we introduce some preliminary results and notions. In Section 3, we obtain the existence of ground state solutions of the system (1.1) and we also investigate their asymptotic behaviour. In Section 4, we get the nonexistence result.
Throughout this paper, we use "→" and "⇀" to denote the strong convergence and weak convergence in the correlation function space, respectively. on(1) denotes a sequence which converges to 0 as n→∞. C will always denote a positive constants, which may vary from line to line.
It is well known that the following properties which follow from the fractional Sobolev embedding
Hs(RN)↪Lq(RN),q∈[2,2∗s], where 2∗s:=2NN−2s. |
If 1+αN<p<α+NN−2s, we have that 2<2NpN+α<2∗s, the space Hsr(RN) compactly embedded into L2NpN+α(RN).
First of all, let us recall the Hardy-Littlewood-Sobolev inequality.
Lemma 2.1. (Hardy-Littlewood-Sobolev inequality [23]) Let 0<α<N, r, q>1 and 1≤s<t<∞ be such that
1r+1q=1+αN,1s−1t=αN. |
(i) For any u∈Lr(RN) and v∈Lq(RN), we have
|∫RN(Iα∗u)v|≤C(N,α,q)‖u‖r‖v‖q. | (2.1) |
If p∈(1+αN,α+NN−2s) and r=q=2NN+α, then
|∫RN(Iα∗|u|p)|u|p|≤C(N,α,p)‖u‖2p2NpN+α, | (2.2) |
where the sharp constant C(N,α,p) is
C(N,α,p)=Cα(N)=πN−α2Γ(α2)Γ(N+α2){Γ(N2)Γ(N)}−αN. |
(ii) For any u∈Ls(RN), we have
‖Iα∗u‖t≤C(N,α,s)‖u‖s. | (2.3) |
Here, C(N,α,s) is a positive constant which depends only on N, α and s, and satisfies
lim supα→0αC(N,α,s)≤2s(s−1)ωN−1, |
where ωN−1 denotes the surface area of the N−1 dimensional unit sphere.
Next, the following result is crucial in the proof of the Theorem 1.1.
Lemma 2.2. Assumption N∈N, 0<α<N and p∈(1+αN,α+NN−2s). Let {un}⊂Hs(RN) be a sequence satisfying that un⇀u weakly in Hs(RN) as n→∞, then
limn→∞∫RN(Iα∗|un|p)|un|p−∫RN(Iα∗|un−u|p)|un−u|p=∫RN(Iα∗|u|p)|u|p. | (2.4) |
To show Lemma 2.2, we state the classical Brezis-Lieb lemma [25].
Lemma 2.3. Let Ω⊆RN be an open subset and 1≤r<∞. If
(i) {un}n∈N is bounded in Lr(Ω).
(ii) un→u almost everywhere on Ω as n→∞, then for every q∈[1,r],
limn→∞∫Ω||un|q−|un−u|q−|u|q|rq=0. | (2.5) |
Here we also need to mention sufficient conditions for weak convergence (see for example [25, Proposition 4.7.12]).
Lemma 2.4. Assume Ω be an open subset of RN, 1<q<∞ and the sequence {un}n∈N is bounded in Lq(Ω). If un→u almost everywhere on Ω as n→∞, we have that un⇀u weakly in Lq(Ω).
In view of Lemmas 2.3 and 2.4 we have the following proof.
Proof of Lemma 2.2. For every n∈N. We have that
∫RN(Iα∗|un|p)|un|p−∫RN(Iα∗|un−u|p)|un−u|p=∫RN(Iα∗(|un|p−|un−u|p))(|un|p−|un−u|p)+2∫RN(Iα∗(|un|p−|un−u|p))|un−u|p. |
Since 1+αN<p<α+NN−2s, we have that 2<2NpN+α<2∗s, then the space Hs(RN) is embedded continuously in L2NpN+α(RN). Moreover, un⇀u weakly in Hs(RN) as n→∞. Thus, the sequence {un}n∈N is bounded in L2NpN+α(RN). By (2.5) with q=p and r=2NpN+α, we have that
|un|p−|un−u|p→|u|p |
strongly in L2NN+α(RN) as n→∞. By (2.3), we have that Iα defines a linear continuous map from L2NN+α(RN) to L2NN−α(RN), then
Iα∗(|un|p−|un−u|p)→Iα∗|u|p |
in L2NN−α(RN) as n→∞. By (2.2), we have
∫RN(Iα∗(|un|p−|un−u|p))(|un|p−|un−u|p)=∫RN(Iα∗|u|p)|u|p+on(1). |
In view of Lemma 2.4, we get |un−u|p⇀0 weakly in L2NN+α(RN) as n→∞. Thus,
∫RN(Iα∗(|un|p−|un−u|p))|un−u|p=on(1). |
The proof is thereby complete.
Lemma 2.5. Let 0<α<N, p∈(1+αN,α+NN−2s) and the sequence {un}n∈N⊂Hs(RN) be such that un⇀u∈Hs(RN) weakly in Hs(RN) as n→∞. Let ϕ∈Hs(RN), we have
limn→∞∫RN(Iα∗|un|p)|un|p−2unϕ=∫RN(Iα∗|u|p)|u|p−2uϕ. | (2.6) |
Proof. Since un⇀u weakly in Hs(RN) as n→∞, then un→u a.e. in RN. By the fractional Sobolev embedding Hs(RN)↪Lq(RN) with q∈[2,2∗s], we see that {un}n∈N is bounded in L2(RN)∩L2∗s(RN). Since 2<2NpN+α<2∗s, then {|un|p} and {|un|q−2un} are bounded in L2NN+α(RN) and Lqq−1(RN) with q∈[2,2∗s], respectively, up to a subsequence, we get
|un|q−2un⇀|u|q−2u weakly in Lqq−1(RN), |
|un|p⇀|u|p weakly in L2NN+α(RN). | (2.7) |
In view of the Rellich theorem, un→u in Ltloc(RN) for t∈[1,2∗s) and |un|p−2un→|u|p−2u in L2Np(p−1)(N+α)loc(RN) (see [26, Theorem A.2]), then we have that |un|p−2unϕ→|u|p−2uϕ in L2NN+α(RN) for any ϕ∈C∞0(RN), where C∞0(RN) denotes the space of the functions infinitely differentiable with compact support in RN. By (2.3), we get
Iα∗(|un|p−2unϕ)→Iα∗(|u|p−2uϕ) | (2.8) |
in L2NN−α(RN). Therefore, by (2.7) and (2.8) we get
∫RN(Iα∗|un|p)|un|p−2unϕ−∫RN(Iα∗|u|p)|u|p−2uϕ=∫RN(Iα∗(|un|p−2unϕ))|un|p−∫RN(Iα∗(|u|p−2uϕ))|u|p=∫RN[Iα∗(|un|p−2unϕ)−Iα∗(|u|p−2uϕ)]|un|p+∫RN(Iα∗(|u|p−2uϕ))(|un|p−|u|p)→ 0 |
as n→∞. Since C∞0(RN) is dense in Hs(RN), we reach the conclusion.
Lemma 2.6. (see [27, Theorem 3.7]) Let f, g and h be three non-negative Lebesgue measurable functions on RN. Let
W(f,g,h):=∫RN∫RNf(x)g(y)h(x−y)dxdy, |
we get
W(f∗,g∗,h∗)≥W(f,g,h), |
where f∗, g∗ and h∗ denote the symmetric radial decreasing rearrangement of f, g and h.
Lemma 2.7. (see [22, Theorem 1.1]) Under the assumptions of Theorem 1.1, there exists a ground state solution u∈Hs(RN) (v∈Hs(RN)) to problem (1.3) (1.4) which is positive, radially symmetric. Moreover, the minima of the energy functional E0,1 (E0,2) on the Nehari manifold N0,1 (N0,2) defined in (1.9) satisfies minu∈N0,1E0,1(u)>0 (minu∈N0,2E0,2(u)>0).
For any (u,v)∈Nβ, we have
Eβ(u,v)=(12−12p)(‖(u,v)‖2H−2β∫RNuvdx)=(12−12p)(∫RN(Iα∗|u|p)|u|pdx+∫RN(Iα∗|v|p)|v|pdx). |
This shows that Eβ is coercive on Nβ. Next we show, through a series of lemmas, that mβ is attained by some (u,v)∈Nβ which is a critical point of Eβ considered on the whole space H, and therefore a ground state solution to (1.1).
We begin with some basic properties of Eβ and Nβ.
Lemma 3.1. For every (u,v)∈H∖{(0,0)}, there exists some t>0 such that (tu,tv)∈Nβ.
Proof. Indeed, (tu,tv)∈Nβ is equivalent to
‖(tu,tv)‖2H=∫RN(Iα∗|tu|p)|tu|p+∫RN(Iα∗|tv|p)|tv|p+2βt2∫RNuv, |
which is solved by
t=(‖(u,v)‖2H−2β∫RNuv∫RN(Iα∗|u|p)|u|p+∫RN(Iα∗|v|p)|v|p)12p−2. | (3.1) |
By inequality
2β∫RNuv<2√λ1λ2∫RNuv≤∫RNλ1u2+λ2v2≤‖u‖2λ1+‖v‖2λ1=‖(u,v)‖2H, |
we have that
‖(u,v)‖2H−2β∫RNuv>‖(u,v)‖2H−‖(u,v)‖2H=0. |
Therefore we get t>0.
Lemma 3.2. The following assertions hold:
(i) There exists c>0 such that ‖(u,v)‖H≥c for any (u,v)∈Nβ.
(ii) mβ=inf(u,v)∈NβEβ(u,v)>0 for all fixed 0<β<√λ1λ2.
(iii) Let u1, v1 are positive solutions of (1.3) and (1.4) respectively, and let t>0 be such that (tu1,tv1)∈Nβ, then 0<t<1.
Proof. (i) In view of the definition of Nβ, by the Hardy-Littlewood-Sobolev inequality (2.2), for any (u,v)∈Nβ, we have
‖u‖2λ1+‖v‖2λ2= ∫RN(Iα∗|u|p)|u|p+∫RN(Iα∗|v|p)|v|p+2β∫RNuv≤C(N,α,p)(‖u‖2p2NpN+α+‖v‖2p2NpN+α)+β√λ1λ2(2√λ1λ2∫RNuv)≤C1C(N,α,p)(‖u‖2pλ1+‖v‖2pλ2)+β√λ1λ2(∫RNλ1u2+λ2v2)≤C1C(N,α,p)(‖u‖2λ1+‖v‖2λ2)p+β√λ1λ2(‖u‖2λ1+‖v‖2λ2), |
where C1>0 denotes the fractional Sobolev embedding constant and C1 does not depend on u and v. This means that
(1−β√λ1λ2)‖(u,v)‖2H≤C1C(N,α,p)‖(u,v)‖2pH. |
Since 0<β<√λ1λ2, we have ‖(u,v)‖H≥c, where
c=(√λ1λ2−βC1C(N,α,p)√λ1λ2)12p−2>0. | (3.2) |
(ii) For any (u,v)∈Nβ, we have
Eβ(u,v)=(12−12p)(‖(u,v)‖2H−2β∫RNuv)≥(12−12p)(‖(u,v)‖2H−β√λ1λ2(‖u‖2λ1+‖v‖2λ2))≥(12−12p)(1−β√λ1λ2)‖(u,v)‖2H. | (3.3) |
Since p>1, we obtain mβ≥(12−12p)(1−β√λ1λ2)c2>0.
(iii) Since u1, v1 are positive solutions of (1.3) and (1.4) respectively, and (tu1,tv1)∈Nβ, we have
‖u1‖2λ1+‖v1‖2λ2=∫RN(Iα∗|u1|p)|u1|p+∫RN(Iα∗|v1|p)|v1|p | (3.4) |
and
t2(‖u1‖2λ1+‖v1‖2λ2−2β∫RNu1v1)=t2p(∫RN(Iα∗|u1|p)|u1|p+∫RN(Iα∗|v1|p)|v1|p). | (3.5) |
Combining (3.4) and (3.5), we have
t2p−2=‖u1‖2λ1+‖v1‖2λ2−2β∫RNu1v1‖u1‖2λ1+‖v1‖2λ2<1. |
The proof is complete.
Proof of Theorem 1.1. Let (un,vn)∈Nβ be a minimizing sequence for Eβ, namely such that Eβ(un,vn)→mβ. By (3.3), we know that {(un,vn)}n∈N is bounded in H. In view of Lemma 3.1, there exists tn>0 such that (tn|un|,tn|vn|)∈Nβ. Then
t2p−2n=‖(|un|,|vn|)‖2H−2β∫RN|un||vn|∫RN(Iα∗|un|p)|un|p+∫RN(Iα∗|vn|p)|vn|p≤‖(un,vn)‖2H−2β∫RNunvn∫RN(Iα∗|un|p)|un|p+∫RN(Iα∗|vn|p)|vn|p=1. |
Hence, we have that 0<tn≤1. Since
Eβ(tn|un|,tn|vn|)=(12−12p)t2pn(∫RN(Iα∗|un|p)|un|p+∫RN(Iα∗|vn|p)|vn|p)≤(12−12p)(∫RN(Iα∗|un|p)|un|p+∫RN(Iα∗|vn|p)|vn|p)=Eβ(un,vn). |
For this reason we can assume that un≥0 and vn≥0. Let u∗n and v∗n denote the symmetric decreasing rearrangement of un, respectively vn. By Lemma 2.6 with f(x)=|un(x)|p, g(y)=|un(y)|p, h(x−y)=|x−y|α−N, we have
∫RN(Iα∗|u∗n|p)|u∗n|p≥∫RN(Iα∗|un|p)|un|p. | (3.6) |
In addition, it is well known that
∫RN|(−Δ)s2u∗n|2≤∫RN|(−Δ)s2un|2and∫RN|u∗n|2=∫RN|un|2 | (3.7) |
(see [28, Theorem 3]). By Hardy-Littlewood inequality and Riesz rearrangement inequality (see [28]),
∫RNu∗nv∗n≥∫RNunvn. | (3.8) |
By (3.6)–(3.8) we have
Eβ(u∗n,v∗n)=12(‖u∗n‖2λ1+‖v∗n‖2λ2)−12p∫RN(Iα∗|u∗n|p)|u∗n|p−12p∫RN(Iα∗|v∗n|p)|v∗n|p−β∫RNu∗nv∗n≤12(‖un‖2λ1+‖vn‖2λ2)−12p∫RN(Iα∗|un|p)|un|p−12p∫RN(Iα∗|vn|p)|vn|p−β∫RNunvn=Eβ(un,vn). |
Therefore, we can further assume that (un,vn)∈Hr. By (3.3), we have that {(un,vn)} is bounded in H, there exists (uβ,vβ)∈H and uβ≥0, vβ≥0 such that up to subsequences, (un,vn)⇀(uβ,vβ) weakly in H. Moreover, we also can assume that un→uβ, vn→vβ a.e. in RN and (uβ,vβ)∈Hr. Since {(un,vn)}n∈N⊂Nβ, we have
∫RN(Iα∗|un|p)|un|p+∫RN(Iα∗|vn|p)|vn|p= ‖un‖2λ1+‖vn‖2λ2−2β∫RNunvn≥(1−β√λ1λ2)‖(un,vn)‖2H≥(1−β√λ1λ2)c2. |
By (2.4), we obtain
∫RN(Iα∗|uβ|p)|uβ|p+∫RN(Iα∗|vβ|p)|vβ|p≥(1−β√λ1λ2)c2>0, |
which means uβ≢0 or vβ≢0.
By (2.4) and Fatou's lemma, we have
‖uβ‖2λ1+‖vβ‖2λ2−2β∫RNuβvβ≤∫RN(Iα∗|uβ|p)|uβ|p+∫RN(Iα∗|vβ|p)|vβ|p. |
Let t>0 such that (tuβ,tvβ)∈Nβ, we have
t=(‖(uβ,vβ)‖2H−2β∫RNuβvβ∫RN(Iα∗|uβ|p)|uβ|p+∫RN(Iα∗|vβ|p)|vβ|p)12p−2≤1. |
Hence,
mβ≤Eβ(tuβ,tvβ)= (12−12p)t2p(∫RN(Iα∗|u|p)|uβ|p+∫RN(Iα∗|vβ|p)|vβ|p)≤(12−12p)(∫RN(Iα∗|uβ|p)|uβ|p+∫RN(Iα∗|vβ|p)|vβ|p)=limn→∞Eβ(un,vn)=mβ. |
Thus, we can deduce that t=1 and mβ is achieved by (uβ,vβ)∈Nβ with uβ≥0, vβ≥0. Now we know that (uβ,vβ) be non-negative and radial ground state solution of (1.1). Since (1.1) has no semitrivial solution, namely (uβ,0) and (0,vβ) are no solutions of (1.1), we infer that uβ≢0 and vβ≢0. By the strong maximum principle, we get uβ>0 and vβ>0, then (uβ,vβ) be positive and radial ground state solution of (1.1).
Next we consider the asymptotic behavior of the ground state solution.
Suppose {βn} be a sequence which satisfies βn∈(0,min{12,√λ1λ2}) and βn→0 as n→∞. Let (uβn,vβn) be the positive radial ground state solution of (1.1) obtained above, we claim {(uβn,vβn)} is bounded in H. Indeed, let ϕ, ψ are the positive solutions of (1.3) and (1.4) respectively. By (iii) of Lemma 3.2, we have that (tnϕ,tnψ)∈Nβn, where 0<tn<1. Hence, by (1.5) and (1.6), we have
Eβn(uβn,vβn)≤Eβn(tnϕ,tnψ)=Eβn(tnϕ,tnψ)−12p⟨E′βn(tnϕ,tnψ),(tnϕ,tnψ)⟩=(12−12p)(‖(tnϕ,tnψ)‖2H−2βnt2n∫RNϕψ)<(12−12p)‖(ϕ,ψ)‖2H:=D. |
Therefore, let c0=min{12,√λ1λ2}, for n large enough, we have
D>Eβn(uβn,vβn)=Eβn(uβn,vβn)−12p⟨E′βn(uβn,vβn),(uβn,vβn)⟩≥(12−12p)(1−βn)‖(uβn,vβn)‖2H>c0(12−12p)‖(uβn,vβn)‖2H, |
from which we deduce that {(uβn,vβn)} is bounded in H. Thus, there exists (u0,v0)∈H such that, up to a subsequences, (uβn,vβn)⇀(u0,v0) in H as n→∞ and u0≥0,v0≥0. Moreover by (3.2) we have that
cn=(√λ1λ2−βnC1C(N,α,p)√λ1λ2)12p−2 |
is an increasing sequence and ‖(uβn,vβn)‖2H>c1>0, hence we have that u0≢0 or v0≢0. It is easy to observe that E′0(u0,v0)=0, thus u0, v0 are the solutions of (1.3) and (1.4), respectively. Since
‖(uβn,vβn)−(u0,v0)‖2H= ⟨E′βn(uβn,vβn)−E′0(u0,v0),(uβn,vβn)−(u0,v0)⟩+∫RN(Iα∗|uβn|p)|uβn|p+∫RN(Iα∗|vβn|p)|vβn|p−∫RN(Iα∗|uβn|p)|uβn|p−2uβnu0−∫RN(Iα∗|vβn|p)|vβn|p−2vβnv0+∫RN(Iα∗|u0|p)(|u0|p−|u0|p−2u0uβn)+∫RN(Iα∗|v0|p)(|v0|p−|v0|p−2v0vβn)+βn∫RN(2uβnvβn−uβnv0−vβnu0), | (3.9) |
by Lemmas 2.1, 2.2, 2.5 and above equality (3.9), we can conclude that (uβn,vβn)→(u0,v0) in H as n→∞.
In view of Lemma 2.7, we can assume that u1, v1 are positive ground state solutions to (1.3) and (1.4) respectively, and let tn>0 such that (tnu1,tnv1)∈Nβn. In view of (iii) of Lemma 3.2, we know that 0<tn<1. Furthermore, by (3.1) we have that
tn=(‖(u,v)‖2H−2βn∫RNuv∫RN(Iα∗|u|p)|u|p+∫RN(Iα∗|v|p)|v|p)12p−2 |
is an increasing sequence and tn>t1>0, then we know that tn→1. Consequently, we have
E0(u1,v1)≤E0(u0,v0)=limn→∞Eβn(uβn,vβn)≤limn→∞Eβn(tnu1,tnv1)=E0(u1,v1). | (3.10) |
Obviously E0(u0,v0) is the sum of the energy of u0 and v0 for the single equation (1.3) and (1.4) respectively, namely
E0(u0,v0)=E0,1(u0)+E0,2(v0), |
where E0,2:Hs(RN)→R is the energy functional of (1.4), which is defined similarly to E0,1, and E0(u1,v1) is the sum of the energy of u1 and v1 for the single equation (1.3) and (1.4), respectively, namely
E0(u1,v1)=E0,1(u1)+E0,2(v1). |
Since u1, v1 are positive ground state solutions to (1.3) and (1.4) respectively, we have
E0,1(u0)≥E0,1(u1)andE0,2(v0)≥E0,2(v1). |
By (3.10), we get E0,1(u0)=E0,1(u1) and E0,2(v0)=E0,2(v1). By Lemma 2.7, we know that u0, v0 are positive ground state solutions of (1.3) and (1.4) respectively.
Let u∗0 and v∗0 denote the symmetric decreasing rearrangement of u0 and v0 respectively. By Lemma 2.6 with f(x)=|u0(x)|p, g(y)=|u0(y)|p, h(x−y)=|x−y|α−N, we have
∫RN(Iα∗|u∗0|p)|u∗0|p≥∫RN(Iα∗|u0|p)|u0|p. | (3.11) |
In addition, we know that
∫RN|(−Δ)s2u∗0|2≤∫RN|(−Δ)s2u0|2and∫RN|u∗0|2=∫RN|u0|2 | (3.12) |
(see [28, Theorem 3]). By (3.11) and (3.12) we have
E0(u∗0,v∗0)=12(‖u∗0‖2λ1+‖v∗0‖2λ2)−12p∫RN(Iα∗|u∗0|p)|u∗0|p−12p∫RN(Iα∗|v∗0|p)|v∗0|p≤12(‖u0‖2λ1+‖v0‖2λ2)−12p∫RN(Iα∗|u0|p)|u0|p−12p∫RN(Iα∗|v0|p)|v0|p=E0(u0,v0). |
Therefore, we can further assume that (u0,v0)∈Hr. This completes the proof of Theorem 1.1.
In this section, in order to prove the nonexistence of nontrivial solutions, we need to use the following Pohožaev identity type:
Lemma 4.1. Let N≥3 and (u,v)∈H be any solution of (1.1). Then, (u,v) satisfies the Pohožaev identity
N−2s2∫[|(−Δ)s2u|2+ |(−Δ)s2v|2]dx+N2∫(λ1|u|2+λ2|v|2)dx=N+α2p(∫(Iα∗|u|p)|u|pdx+∫(Iα∗|v|p)|v|pdx)+Nβ∫uvdx. | (4.1) |
Proof. The proof is similar to the argument of Theorem 1.13 in [22].
Proof of Theorem 1.2. Let ⟨E′β(u,v),(u,v)⟩=0, by (1.6), we have
∫[|(−Δ)s2u|2+|(−Δ)s2v|2+λ1|u|2+λ2|v|2]dx=∫(Iα∗|u|p)|u|pdx+∫(Iα∗|v|p)|v|pdx+2β∫uvdx | (4.2) |
for all (u,v)∈H.
Combining the Pohožaev identity (4.1) and (4.2), we can see that
0= (N−2s−N+αp)∫[|(−Δ)s2u|2+|(−Δ)s2v|2]dx+(N−N+αp)∫(λ1|u|2+λ2|v|2)dx+(N+αp−N)∫2βuvdx.= (N−2s−N+αp)∫[|(−Δ)s2u|2+|(−Δ)s2v|2]dx+(N−N+αp)∫(λ1|u|2+λ2|v|2−2βuv)dx. | (4.3) |
Since λ1>0, λ2>0 and 0<β<√λ1λ2, we have
λ1|u|2+λ2|v|2≥2√λ1λ2uv>2βuv. |
Thus, if both the coefficients are non-positive, that is
N−2s−N+αp≤0andN−N+αp≤0, |
then we get p≤1+αN, which jointly with (4.3) leads us to a contradiction. Therefore, the solution of (1.1) is the trivial one. Similarly, if they are nonnegative, that is p≥N+αN−2s, we get that nontrivial solutions of (1.1) cannot exist. Therefore, the range of 1+αN<p<N+αN−2s is optimal for the existence of nontrival solutions of the problem (1.1). This completes the proof.
In this present paper, we combine the critical point theory and variational method to investigate a class of coupled fractional systems of Choquard type. By using constrained minimization method and Hardy-Littlewood-Sobolev inequality, we establish the existence and asymptotic behaviour of positive ground state solutions of the systems. Furthermore, nonexistence of nontrivial solutions is also obtained. In the next work, we will focus on the research of normalized solutions to fractional couple Choquard systems.
This research was funded by the National Natural Science Foundation of China (61803236) and Natural Science Foundation of Shandong Province (ZR2018MA022).
The authors declare that they have no conflicts of interest.
[1] |
S. Araci, Degenerate poly type 2-Bernoulli polynomials, Math. Sci. Appl. E-Notes, 9 (2021), 1–8. https://doi.org/10.36753/mathenot.839111 doi: 10.36753/mathenot.839111
![]() |
[2] | A. Bayad, D. S. Kim, Some identities for the Bernoulli, the Euler and the Genocchi numbers and polynomials, Adv. Stud. Contemp. Math., 20 (2010), 23–28. |
[3] |
A. Bayad, M. Hajli, On the multidimensional zeta functions associated with theta functions, and the multidimensional Appell polynomials, Math. Method. Appl. Sci., 43 (2020), 2679–2694 https://doi.org/10.1002/mma.6075 doi: 10.1002/mma.6075
![]() |
[4] |
A. Bayad, T. Kim, Results on values of Barnes polynomials, Rocky Mountain J. Math., 43 (2013), 1857–1869. https://doi.org/10.1216/RMJ-2013-43-6-1857 doi: 10.1216/RMJ-2013-43-6-1857
![]() |
[5] |
L. Carlitz, A degenerate Staudt-Clausen theorem, Arch. Math., 7 (1956), 28–33. https://doi.org/10.1007/BF01900520 doi: 10.1007/BF01900520
![]() |
[6] | L. Carlitz, Degenerate Stirling, Bernoulli and Eulerian numbers, Util. Math., 15 (1979), 51–88. |
[7] |
U. Duran, M. Acikgoz, On generalized degenerate Gould-Hopper based fully degenerate Bell polynomials, JMCS, 21 (2020), 243–257. http://doi.org/10.22436/jmcs.021.03.07 doi: 10.22436/jmcs.021.03.07
![]() |
[8] |
M. Hajli, The spectral properties of a continuous family of zeta functions, Int. J. Number Theory, 16 (2020), 693–717. https://doi.org/10.1142/S1793042120500359 doi: 10.1142/S1793042120500359
![]() |
[9] |
H. Haroon, W. A. Khan, Degenerate Bernoulli numbers and polynomials associated with degenerate Hermite polynomials, Commun. Korean Math. Soc., 33 (2018), 651–669. https://doi.org/10.4134/CKMS.c170217 doi: 10.4134/CKMS.c170217
![]() |
[10] |
M. Hajli, On a formula for the regularized determinant of zeta functions with application to some Dirichlet series, Q. J. Math., 71 (2020), 843–865. https://doi.org/10.1093/qmathj/haaa006 doi: 10.1093/qmathj/haaa006
![]() |
[11] | G. W. Jang, T. Kim, A note on type 2 degenerate Euler and Bernoulli polynomials, Adv. Stud. Contemp. Math., 29 (2019), 147–159. |
[12] |
L. C. Jang, D. S. Kim, T. Kim, H. Lee, p-adic integral on Zp associated with degenerate Bernoulli polynomials of the second kind, Adv. Differ. Equ., 2020 (2020), 278. https://doi.org/10.1186/s13662-020-02746-2 doi: 10.1186/s13662-020-02746-2
![]() |
[13] | M. Kaneko, Poly-Bernoulli numbers, J. Theor. Nombr. Bordx., 9 (1997), 221–228. |
[14] |
W. A. Khan, R. Ali, K. A. H. Alzobyd, N. Ahmed, A new family of degenerate poly-Genocchi polynomials with its certain properties, J. Funct. Spaces, 2021 (2021), 6660517. https://doi.org/10.1155/2021/6660517 doi: 10.1155/2021/6660517
![]() |
[15] | D. S. Kim, T. Kim, Lah-Bell numbers and polynomials, Proc. Jangjeon Math. Soc., 23 (2020), 577–586. |
[16] |
D. S. Kim, T. Kim, A note on a new type of degenerate Bernoulli numbers, Russ. J. Math. Phys., 27 (2020), 227–235. https://doi.org/10.1134/S1061920820020090 doi: 10.1134/S1061920820020090
![]() |
[17] |
D. S. Kim, T. Kim, Degenerate Sheffer sequences and λ-Sheffer sequences, J. Math. Anal. Appl., 493 (2021), 124521. https://doi.org/10.1016/j.jmaa.2020.124521 doi: 10.1016/j.jmaa.2020.124521
![]() |
[18] |
T. Kim, D. S. Kim, Degenerate Laplace transform and degenerate gamma function, Russ. J. Math. Phys., 24 (2017), 241–248. https://doi.org/10.1134/S1061920817020091 doi: 10.1134/S1061920817020091
![]() |
[19] |
T. Kim, D. S. Kim, A note on type 2 Changhee and Daehee polynomials, RACSAM, 113 (2019), 2763–2771. https://doi.org/10.1007/s13398-019-00656-x doi: 10.1007/s13398-019-00656-x
![]() |
[20] |
T. Kim, D. S. Kim, Degenerate polyexponential functions and degenerate Bell polynomials, J. Math. Anal. Appl., 487 (2020), 124017. https://doi.org/10.1016/j.jmaa.2020.124017 doi: 10.1016/j.jmaa.2020.124017
![]() |
[21] |
T. Kim, D. S. Kim, J. Kwon, H. Lee, Degenerate polyexponential functions and type 2 degenerate poly-Bernoulli numbers and polynomials, Adv. Differ. Equ., 2020 (2020), 168. https://doi.org/10.1186/s13662-020-02636-7 doi: 10.1186/s13662-020-02636-7
![]() |
[22] |
T. Kim, D. S. Kim, H. Lee, L. Jang, A note on degenerate derangement polynomials and numbers, AIMS Mathematics, 6 (2021), 6469–6481. https://doi.org/10.3934/math.2021380 doi: 10.3934/math.2021380
![]() |
[23] | T. Kim, A note on degenerate Stirling polynomials of the second kind, Proc. Jangjeon Math. Soc., 20 (2017), 319–331. |
[24] | T. K. Kim, D. S. Kim, H. I. Kwon, A note on degenerate Stirling numbers and their applications, Proc. Jangjeon Math. Soc., 21 (2018), 195–203. |
[25] |
F. Qi, J. L. Wang, B. N. Guo, Simplifying differential eqpuations concerning degenerate Bernoulli and Euler numbers, T. A. Razmadze Math. In., 172 (2018), 90–94. https://doi.org/10.1016/j.trmi.2017.08.001 doi: 10.1016/j.trmi.2017.08.001
![]() |
[26] | S. Roman, The umbral calculus (Pure and Applied Mathematics), London: Academic Press, 1984. |
1. | Hamza El-Houari, Hicham Moussa, On a class of generalized Choquard system in fractional Orlicz-Sobolev spaces, 2024, 540, 0022247X, 128563, 10.1016/j.jmaa.2024.128563 | |
2. | Huiqin Lu, Kexin Ouyang, Existence of Positive Ground State Solutions for Fractional Choquard Systems in Subcritical and Critical Cases, 2023, 11, 2227-7390, 2938, 10.3390/math11132938 |