Existence of Solutions for Klein-Gordon-Maxwell Equations Involving Hardy-Sobolev Critical Exponents

Xin Zhang

School of Mathematics, Liaoning Normal University, Dalian, China.

DOI: 10.4236/oalib.1107709                

Abstract

We investigate the existence of solutions for Klein-Gordon-Maxwell equations involving Hardy-Sobolev critical exponents. By means of the Ekeland’s variational principle and the Mountain Pass Theorem, we obtain that there is at least a nontrivial solution for the subcritical system. Then we prove that there are at least two different solutions for the critical system.

Keywords

Existence of Solutions, Klein-Gordon-Maxwell Equations, Hardy-Sobolev Critical Exponents

1. Introduction and Preliminaries

In recent years, great attention has been given to problems driven by the Laplacian. One of the reasons for this comes from the fact that this operator appears in several applications in different subjects, such as flame propagation, free boundary obstacle problems, and ultrarelativistic limits of quantum mechanics. In particular, from a probabilistic point of view, the Laplace operator is the infinitesimal generator of a Lévy process. For more details and applications, see other references [1] [2].

Problems with two nonlinearities recently have been studied by several authors. In particular, such problems were considered in [3] [4] for the Laplacian, the p-Laplacian, the Biharmonic operator and the fractional Laplacian. In [5], Ghoussoub, Robert and Shakerian investigated problems with doubly critical nonlinear terms, with either critical Sobolev term or critical Hardy-Sobolev term, for the Laplacian and the fractional Laplacian.

Solutions of critical Sobolev problems were found in [6], as critical points of a suitable functional, by the Mountain-Pass lemma without the ${\left(PS\right)}_c$ condition. In this case, the ${\left(PS\right)}_c$ condition only holds true for c in certain intervals related to the best Sobolev constant. In the control of the Mountain-Pass level, the extremal function of the best Sobolev constant plays an important role.

For example, Jannelli (see [7] ) considered that the problem

$(\begin{array}{c}\begin{array}{l}-\Delta u-\lambda u=\mu \frac{u}{{\left|x\right|}^2}+\frac{{\left|u\right|}^{2^{\ast }\left(s\right)-2}}{{\left|x\right|}^s}u,x\in \Omega ,\\ u=0,x\in \partial \Omega ,\end{array}\end{array}$ (1)

they proved the existence of nontrivial solutions for the preceding equations involving a critical Hardy-Sobolev exponent in a bounded domain.

Then Kang and Peng in [8] considered the following problem based on the above equations

$(\begin{array}{c}\begin{array}{l}-\Delta u-\mu \frac{u}{{\left|x\right|}^2}=\frac{u^{2^{\ast }\left(s\right)-2}}{{\left|x\right|}^s}u+\lambda {\left|u\right|}^{r-2}u,x\in \Omega ,\\ u=0,x\in \partial \Omega .\end{array}\end{array}$ (2)

They established the problem above has at least a pair of sign-changing solutions with $\lambda >0$ and $2<r<2^{\ast }$ in this reference [8]. Kang and Peng in [9] proved that problem above has at least one positive solution under some conditions for $r\mathrm{,}\mu$ and $\lambda$. They also proved in [10] that as $r=2$, the problem above has at least a pair of sign-changing solutions for $N\ge 7$, $\mu \in \left[\mathrm{0,}\bar{\mu }-4\right)$ and $0<\lambda <{\lambda }_1\left(u\right)$.

Motivated by the study of solitary waves of the nonlinear Klein-Gordon equation interacting with an electromagnetic field, Benci and Fortunato derived in [11] a model that is described by the following elliptic system

$(\begin{array}{c}\begin{array}{l}-\Delta u+\left[m_0^2-{\left(\omega +\varphi \right)}^2\right]u={\left|u\right|}^{p-2}u\mathrm{,}\\ \Delta \varphi =\left(\omega +\varphi \right)u^2\mathrm{,}\end{array}\end{array}$ (3)

where $m_0$ and $\omega$ are real constants. They proved existence of infinitely many radially symmetric solutions $\left(u\mathrm{,}\varphi \right)\in H^1\left({ℝ}^3\right)\times D^{\mathrm{1,2}}\left({ℝ}^3\right)$ for the above system when $\left|m_0\right|>\left|\omega \right|$ and for sub-critical exponents p satisfying $4<p<2^{\ast }$. As in [12], they derived a variational identity to prove the non-existence of nontrivial weak solutions for the system above. In [13], Cassani investigated the critical case i.e. $p=2^{\ast }=\frac{2N}{N-2}$, the critical Sobolev exponent. Moreover, Cassani used a Pohozaev-type argument, which points out an invariance property for the problem (3), to prove non-existence of solutions with a suitable decay at infinity and in particular it turns out to be the case of radially symmetric solutions. In [13], Cassani replaced the first equation of the System (3), adding a lower order perturbation, by the following

$-\Delta u+\left[m_0^2-{\left(\omega +\varphi \right)}^2\right]u=\mu {\left|u\right|}^{p-2}u+{\left|u\right|}^{2^{\ast }-2}u\mathrm{,}$ (4)

where $\mu >0$ and $4\le p<6=2^{\ast }$. In this case, they recovered a Mountain-Pass type solution for Equation (4) and the second equation of System (3).

The above-mentioned equations with Hardy-Sobolev critical exponents are restricted to bounded regions, and the two Klein-Gordon-Maxwell systems mentioned above involve Sobolev critical exponents, not Hardy-Sobolev critical exponents. Moreover, some parts of the certification process did not give a specific certification process. We are inspired by the proof methods of a nontrivial solution and infinitely many solutions in the above-mentioned literatures, and investigate the existence of solutions for subcritical equations and critical equations with Hardy-Sobolev critical exponents in ${ℝ}^3$.

In this paper, firstly we study existence of solutions for the following Klein-Gordon-Maxwell equations involving Hardy-Sobolev critical exponents

$(\begin{array}{c}\begin{array}{ll}-\Delta u+\left[m^2-{\left(\omega +\varphi \right)}^2\right]u=\frac{{\left|u\right|}^{p-2}u}{{\left|x\right|}^s}\mathrm{,}\hfill & x\in {ℝ}^3\hfill \\ -\Delta \varphi +u^2\varphi =-\omega u^2\mathrm{,}\hfill & x\in {ℝ}^3\hfill \end{array}\end{array}$ (5)

where m and $\omega$ are real constants, $2<p<2^{\ast }\left(s\right)=\frac{2\left(N-s\right)}{N-2}$, $0<s<2$.

Then we study existence of solutions for the following equations

$(\begin{array}{c}\begin{array}{ll}-\Delta u+\left[m^2-{\left(\omega +\varphi \right)}^2\right]u=\mu \frac{{\left|u\right|}^{q-2}u}{{\left|x\right|}^s}+\frac{{\left|u\right|}^{p-2}u}{{\left|x\right|}^s}+\frac{{\left|u\right|}^{2^{\ast }\left(s\right)-2}u}{{\left|x\right|}^s}\mathrm{,}\hfill & x\in {ℝ}^3\hfill \\ -\Delta \varphi +u^2\varphi =-\omega u^2\mathrm{,}\hfill & x\in {ℝ}^3\hfill \end{array}\end{array}$ (6)

where m, $\mu$ and $\omega$ are real constants, $1<q<2$, $2<p<2^{\ast }\left(s\right)=\frac{2\left(N-s\right)}{N-2}$, $0<s<2$.

In Section 2, we first give main results for the Systems (5) and (6). In Section 3, we prove the existence of solutions of System (5). In Section 4, we establish existence of solutions of System (6).

2. Main Results

Throughout this paper, we denote the $L^p\left({ℝ}^3\mathrm{,}{\left|x\right|}^{-s}\text{d}x\right)$ norm by

${\left|u\right|}_{p\mathrm{,}s}={\left({\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x\right)}^{\frac{1}{p}}$

and $L^p\left({ℝ}^3\mathrm{,}\text{d}x\right)$ norm by ${\left|u\right|}_p={\left({\displaystyle {\int }_{{ℝ}^3}}{\left|u\right|}^p\text{d}x\right)}^{\frac{1}{p}}$. For simplicity, set ${\left|u\right|}_{2^{\ast }\left(0\right),0}={\left|u\right|}_{2^{\ast }}={\left|u\right|}_6$, since $2^{\ast }=\frac{2N}{N-2}$.

$D^{\mathrm{1,2}}\left({ℝ}^3\right)=\left\{u\in L^{2^{\ast }}\left({ℝ}^3\right)\mathrm{<mo>\; :\; </mo>}\left|\nabla u\right|\in L^2\left({ℝ}^3\right)\right\}$ is a Sobolev space with norm

${\Vert u\Vert }_{D^{\mathrm{1,2}}}={\left({\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla u\right|}^2\text{d}x\right)}^{\frac{1}{2}}\mathrm{.}$

Moreover, $H^1\left({ℝ}^3\right)$ is the usual Sobolev space with norm

${\Vert u\Vert }_{H^1}={\left({\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+u^2\right)\text{d}x\right)}^{\frac{1}{2}}$

continuously embedded in $L^p\left({ℝ}^3\right)$ for $2\le p\le 2^{\ast }\left(s\right)$. Define Hardy-Sobolev best constant as follows

$S_s=\underset{u\in D^{\mathrm{1,2}}\left({ℝ}^3\right)\mathrm{,}u\ne 0}{inf}\frac{{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla u\right|}^2\text{d}x}{{\left({\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\right)}^{\frac{2}{2^{\ast }\left(s\right)}}}\mathrm{.}$

Theorem 2.1. If one of the following conditions is satisfied:

1) $2<p<4$ and $\left(\frac{p}{2}-1\right)m^2>{\omega }^2$ , or

2) $4\le p<2^{\ast }\left(s\right)$ and $\left|m\right|>\left|\omega \right|$.

Then System (5) admits at least a nontrivial solution.

Theorem 2.2. If one of the following conditions is satisfied:

1) $4\le p<2^{\ast }\left(s\right)$ and $\left|m\right|>\left|\omega \right|$ , or

2) $2<p<4$ and $\left(\frac{p}{2}-1\right)m^2>{\omega }^2$.

Then there exists a constant $\tilde{m}>0$ such that System (6) admits at least two different solutions $\left(u\mathrm{,}\varphi \right)$ satisfying ${\Vert u\Vert }_{H^1}<+\infty$, ${\Vert \varphi \Vert }_{D^{1,2}}<+\infty$ when $0<\mu <\tilde{m}$.

We define the functional of System (5)

$\begin{array}{l}F\left(u,\varphi \right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2-{\left|\nabla \varphi \right|}^2+\left[m^2-{\left(\omega +\varphi \right)}^2\right]u^2\right)\text{d}x-\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x.\hfill \end{array}$ (7)

Define the functional of System (6) as follows

$\begin{array}{c}\tilde{F}\left(u,\varphi \right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2-{\left|\nabla \varphi \right|}^2+\left[m^2-{\left(\omega +\varphi \right)}^2\right]u^2\right)\text{d}x\\ -\frac{\mu }{q}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^q}{{\left|x\right|}^s}\text{d}x-\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x-\frac{1}{2^{\ast }\left(s\right)}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\mathrm{.}\end{array}$ (8)

Remark 2.1. The functional F and $\tilde{F}$ are strongly indefinite i.e. unbounded from below and from above on infinite dimensional subspaces. In order to avoid this indefiniteness, which rules out many of the usual tools of critical point theory, a reduction method is performed in [11] which we now recall. For u and $\varphi$ defined above, we have the following lemmas.

3. The Proof of Theorem 2.1

Lemma 3.1. For every $u\in H^1\left({ℝ}^3\right)$ ,

1) there exists a unique function $\varphi =\Phi \left(u\right)\in D^{\mathrm{1,2}}\left({ℝ}^3\right)$ that solves the second equation of System (5);

2) if u is radially symmetric, then $\Phi \left(u\right)$ is radial too;

3) $\Phi \left(u\right)\left(x\right)\le 0$, moreover, $\Phi \left(u\right)\left(x\right)\ge -\omega$, if $u\left(x\right)\ne 0$ and $\omega >0$.

Proof. The first result is proved in Lemma 3 of [14]. While the second one, though not explicitly stated, is proved in Lemma 5 of [14]. The third result can be found in Lemma 2.3 of [15]. □

Lemma 3.2. The map $\Phi$ is $C^1$ and $G_{\varphi }=\left\{\left(u,\varphi \right)\in H^1\left({ℝ}^3\right)\times D^{\mathrm{1,2}}\left({ℝ}^3\right)\mathrm{<mo>\; |\; </mo>}{F^{\prime }}_{\varphi }\left(u\mathrm{,}\varphi \right)=0\right\}$.

Proof. Noticing that $\Phi \left(u\right)$ is a solution of the second equation in System (5), we have

$-{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(u\right)\right|}^2\text{d}x={\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(u\right)u^2\text{d}x+{\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u\right)u^2\text{d}x\mathrm{.}$ (9)

In addition,

$\begin{array}{c}F\left(u\mathrm{,}\Phi \left(u\right)\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2-{\left|\nabla \Phi \left(u\right)\right|}^2+\left[m^2-{\left(\omega +\Phi \left(u\right)\right)}^2\right]u^2\right)\text{d}x\\ -\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x\mathrm{,}\end{array}$

${F^{\prime }}_{\varphi }\left(u\mathrm{,}\Phi \left(u\right)\right)=-{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(u\right)\right|}^2\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(u\right)u^2\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u\right)u^2\text{d}x\mathrm{.}$

According to (9), one gets ${F^{\prime }}_{\varphi }\left(u\mathrm{,}\Phi \left(u\right)\right)=0$ for any $\left(u,\varphi \right)\in H^1\left({ℝ}^3\right)\times D^{\mathrm{1,2}}\left({ℝ}^3\right)$. Thus $F^{\prime }\left(u\mathrm{,}\Phi \left(u\right)\right)={F^{\prime }}_u\left(u\mathrm{,}\Phi \left(u\right)\right)+{F^{\prime }}_{\varphi }\left(u\mathrm{,}\Phi \left(u\right)\right){\Phi }^{\prime }\left(u\right)={F^{\prime }}_u\left(u\mathrm{,}\Phi \left(u\right)\right)$. □

Define $I\left(u\right)=F\left(u,\varphi \right)$. If $u\mathrm{,}v\in H^1\left({ℝ}^3\right)$, then one has

$I^{\prime }\left(u\right)v={\displaystyle {\int }_{{ℝ}^3}}\left(\nabla u\cdot \nabla v+\left[m^2-{\left(\omega +\varphi \right)}^2\right]uv\right)\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{p-2}}{{\left|x\right|}^s}uv\text{d}x\mathrm{.}$ (10)

Lemma 3.3. The following statements are equivalent:

1) $\left(u,\varphi \right)\in H^1\left({ℝ}^3\right)\times D^{\mathrm{1,2}}\left({ℝ}^3\right)$ is a solution of System (5);

2) u is a critical point for I and $\varphi =\Phi \left(u\right)$.

Proof. 2) $\Rightarrow$ 1) Obviously.

1) $\Rightarrow$ 2) Suppose ${F^{\prime }}_u\left(u,\varphi \right)$ and ${F^{\prime }}_{\varphi }\left(u,\varphi \right)$ denote the partial derivatives of F at $\left(u,\varphi \right)\in H^1\left({ℝ}^3\right)\times D^{\mathrm{1,2}}\left({ℝ}^3\right)$. Then for every $v\in H^1\left({ℝ}^3\right)$ and $\psi \in D^{\mathrm{1,2}}\left({ℝ}^3\right)$, one gets

$\begin{array}{l}{F^{\prime }}_u\left(u,\varphi \right)\left[v\right]={\displaystyle {\int }_{{ℝ}^3}}\left(\nabla u\cdot \nabla v+\left[m^2-{\left(\omega +\varphi \right)}^2\right]uv\right)\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{p-2}}{{\left|x\right|}^s}uv\text{d}x\mathrm{,}\hfill \end{array}$ (11)

${F^{\prime }}_{\varphi }\left(u,\varphi \right)\left[\psi \right]=-{\displaystyle {\int }_{{ℝ}^3}}\nabla \varphi \cdot \nabla \psi \text{d}x-{\displaystyle {\int }_{{ℝ}^3}}\left(\omega +\varphi \right)\psi u^2\text{d}x\mathrm{.}$ (12)

By the standard computations, we can prove that ${F^{\prime }}_u\left(u,\varphi \right)$ and ${F^{\prime }}_{\varphi }\left(u,\varphi \right)$ are continuous. From (11) and (12), it is easy to obtain that its critical points are solutions of System (5), by 1) of Lemma 3.1, one has $\varphi =\Phi \left(u\right)$. □

Lemma 3.4. For $u\in H^1\left({ℝ}^3\right)$ , if $\left|m\right|>\left|\omega \right|$ , then there exist some constants ${\rho }_1,{\alpha }_1>0$ such that ${I\left(u\right)|}_{{\Vert u\Vert }_{H^1}={\rho }_1}\ge {\alpha }_1>0$ .

Proof. From (7), one obtains

$\begin{array}{c}I\left(u\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2\right)\text{d}x-\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \varphi \right|}^2\text{d}x\\ -{\displaystyle {\int }_{{ℝ}^3}}\omega \varphi u^2\text{d}x-\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\varphi }^2{u}^2\text{d}x-\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x\mathrm{.}\end{array}$ (13)

Substituting (9) into (13), we have

$\begin{array}{c}I\left(u\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2\right)\text{d}x+\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \varphi \right|}^2\text{d}x\\ +\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\varphi }^2{u}^2\text{d}x-\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x\\ \ge \frac{min\left\{\mathrm{1,}{m}^2-{\omega }^2\right\}}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+u^2\right)\text{d}x-\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x\\ \ge C_1{\Vert u\Vert }_{H^1}^2-C_2{\Vert u\Vert }_{H^1}^p\ge {\alpha }_1\mathrm{,}\text{for}u\in H^1\left({ℝ}^3\right)\mathrm{,}{\Vert u\Vert }_{H^1}={\rho }_1\mathrm{.}\end{array}$

Thus

${I\left(u\right)|}_{{\Vert u\Vert }_{H^1}={\rho }_1}\ge {\alpha }_1>0$

and the proof is completed. □

Lemma 3.5. Under the assumptions of Theorem 2.1, there exists a function ${\eta }_1\in H^1\left({ℝ}^3\right)$ with ${\Vert {\eta }_1\Vert }_{H^1}>{\rho }_1$ such that $I\left({\eta }_1\right)<0$ .

Proof. It is easy to obtain

$\begin{array}{c}\underset{t\to +\infty }{lim}I\left(tu\right)=\frac{t^2}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2\right)\text{d}x-t^2{\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(tu\right)u^2\text{d}x\\ -\frac{t^2}{2}{\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(tu\right)u^2\text{d}x-\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(tu\right)\right|}^2\text{d}x-\frac{t^p}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x\\ \le \frac{t^2}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2\right)\text{d}x-\omega t^2{\displaystyle {\int }_{{ℝ}^3}}\Phi \left(tu\right)u^2\text{d}x-\frac{t^p}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x\\ \le \frac{t^2}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2+2{\omega }^2{u}^2\right)\text{d}x-\frac{t^p}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x\\ \to -\infty \mathrm{,}\end{array}$

which implies that $I\left(u\right)\to -\infty$, as ${\Vert u\Vert }_{H^1}\to \infty$.

The lemma is proved by taking ${\eta }_1=tu$ with $t>0$ large enough and $u\ne 0$. Therefore we know that there exists ${\eta }_1\in H^1\left({ℝ}^3\right)$, ${\Vert {\eta }_1\Vert }_{H^1}>{\rho }_1$ such that $I\left({\eta }_1\right)<0$. □

Therefore, there exists a sequence $\left\{u_n\right\}\in H^1\left({ℝ}^3\right)$, so-called Palais-Smale sequence, such that

$I\left(u_n\right)\to c>\mathrm{0,}\text{and}{I}^{\prime }\left(u_n\right)\to \mathrm{0,}n\to \infty \mathrm{,}$ (14)

where

$c=\underset{\gamma \in \Gamma }{inf}\underset{t\in \left[\mathrm{0,1}\right]}{max}I\left(\gamma (\; t\; )\right)$

with

$\Gamma =\left\{\gamma \in C\left(\left[0,1\right],H^1\left({ℝ}^3\right)\right):\gamma \left(0\right)=0,\gamma \left(1\right)={\eta }_1\right\}.$

Since System (5) is set on ${ℝ}^3$, it is well known that the Sobolev embedding $H^1\left({ℝ}^3\right)$ ↪ $L^p\left({ℝ}^3\mathrm{,}{\left|x\right|}^{-s}\text{d}x\right)\left(2\le p\le 2^{\ast }\left(s\right)\right)$ is not compact and then it is usually difficult to prove that a Palais-Smale sequence is strongly convergent when we seek solutions of System (5) by variational methods. A standard tool to overcome the problem is to restrict ourselves to radial functions, namely we look at the functional I on the subspace $H_r^1\left({ℝ}^3\right)=\left\{u\in H^1\left({ℝ}^3\right)\mathrm{<mo>\; |\; </mo>}u\left(x\right)=u\left(\left|x\right|\right)\right\}$ and $D_r^{\mathrm{1,2}}\left({ℝ}^3\right)=\left\{u\in D^{\mathrm{1,2}}\left({ℝ}^3\right)\mathrm{<mo>\; |\; </mo>}u\left(x\right)=u\left(\left|x\right|\right)\right\}$ compactly embedded in $L_r^p\left({ℝ}^3\mathrm{,}{\left|x\right|}^{-s}\text{d}x\right)$ for $2<p<2^{\ast }\left(s\right)$ and $L_r^p\left({ℝ}^3\mathrm{,}\text{d}x\right)$ for $2<p<2^{\ast }=\frac{2N}{N-2}=6$. By standard arguments, one sees that if a critical point $u\in H_r^1\left({ℝ}^3\right)$ for the functional ${I|}_{H_r^1\left({ℝ}^3\right)}$ is also a critical point of I.

Lemma 3.6. The PS sequence $\left\{u_n\right\}\in H_r^1\left({ℝ}^3\right)$ given in (14) is bounded. Moreover, $\Phi \left(u_n\right)\in D_r^{\mathrm{1,2}}\left({ℝ}^3\right)$ is bounded, too.

Proof. Case 1. $2<p<4$. There exists a constant $C_3>0$, then by (7), (9) and (10), we get

$\begin{array}{l}{C}_3+o\left(1\right){\Vert u_n\Vert }_{H_r^1}\\ \ge pI\left(u_n\right)-\langle I^{\prime }\left(u_n\right)\mathrm{,}{u}_n\rangle \\ =\left(\frac{p}{2}-1\right){\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u_n\right|}^2+\left(m^2-{\omega }^2\right)u_n^2\right)\text{d}x-\left(\frac{p}{2}-2\right){\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(u_n\right)u_n^2\text{d}x\\ +{\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u_n\right)u_n^2\text{d}x\\ \ge \left(\frac{p}{2}-1\right){\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla u_n\right|}^2\text{d}x+{\displaystyle {\int }_{{ℝ}^3}}\left[\left(\frac{p}{2}-1\right)\left(m^2-{\omega }^2\right)+\left(\frac{p}{2}-2\right){\omega }^2\right]u_n^2\text{d}x\\ \ge C_4{\Vert u_n\Vert }_{H_r^1}^2\end{array}$

for n large enough. Therefore, it follows that $\left\{u_n\right\}$ is bounded in $H_r^1\left({ℝ}^3\right)$.

Case 2. $4\le p<2^{\ast }\left(s\right)$. We have that there is a positive constant ${\tilde{C}}_3$ such that

${\tilde{C}}_3+{\Vert u_n\Vert }_{H_r^1}\ge I\left(u_n\right)-\frac{1}{p}\langle I^{\prime }\left(u_n\right)\mathrm{,}{u}_n\rangle \mathrm{.}$

According to (9), (10) and (13), $\Phi \left(u_n\right)\ge -\omega$, one has

$\begin{array}{l}{\tilde{C}}_3+{\Vert u_n\Vert }_{H_r^1}\\ \ge I\left(u_n\right)-\frac{1}{p}\langle I^{\prime }\left(u_n\right)\mathrm{,}{u}_n\rangle \\ =\left(\frac{1}{2}-\frac{1}{p}\right){\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u_n\right|}^2+\left(m^2-{\omega }^2\right)u_n^2\right)\text{d}x+\left(\frac{1}{2}+\frac{1}{p}\right){\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u_n\right)u_n^2\text{d}x\\ +\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(u_n\right)\right|}^2\text{d}x+\frac{2}{p}{\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(u_n\right)u_n^2\text{d}x\end{array}$

$\begin{array}{l}=\left(\frac{1}{2}-\frac{1}{p}\right){\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u_n\right|}^2+\left(m^2-{\omega }^2\right)u_n^2\right)\text{d}x+\left(\frac{1}{2}-\frac{1}{p}\right){\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u_n\right)u_n^2\text{d}x\\ +\left(\frac{1}{2}-\frac{2}{p}\right){\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(u_n\right)\right|}^2\text{d}x\\ \ge C_5{\Vert u_n\Vert }_{H_r^1}^2\end{array}$

for n large enough. It follows that $\left\{u_n\right\}$ is bounded in $H_r^1\left({ℝ}^3\right)$.

According to Equation (9), one has

${\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(u_n\right)\right|}^2\text{d}x=-{\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(u_n\right)u_n^2\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u_n\right)u_n^2\text{d}x\mathrm{,}$

then by Hölder inequality

$\begin{array}{l}{\displaystyle {\int }_{{ℝ}^3}}\left|f\left(x\right)\cdot g\left(x\right)\right|\text{d}x\le {\left({\displaystyle {\int }_{{ℝ}^3}}{\left|f\left(x\right)\right|}^p\text{d}x\right)}^{\frac{1}{p}}\cdot {\left({\displaystyle {\int }_{{ℝ}^3}}{\left|g\left(x\right)\right|}^q\text{d}x\right)}^{\frac{1}{q}}\\ ={\left|f\left(x\right)\right|}_p\cdot {\left|g\left(x\right)\right|}_q\mathrm{,}\left(p>1,q>1,\frac{1}{p}+\frac{1}{q}=1\right)\end{array}$

and Sobolev inequality

${\left|u\right|}_q\le C{\Vert u\Vert }_{H_r^1}\mathrm{.}$

One obtains

${\Vert \Phi \left(u_n\right)\Vert }_{D_r^{\mathrm{1,2}}}^2\le C_6{\Vert \Phi \left(u_n\right)\Vert }_{D_r^{\mathrm{1,2}}}{\Vert u_n\Vert }_{H_r^1}^2\mathrm{.}$

Thus $\left\{\Phi \left(u_n\right)\right\}$ is bounded in $D_r^{\mathrm{1,2}}\left({ℝ}^3\right)$ by the boundedness of $\left\{u_n\right\}\subset H_r^1\left({ℝ}^3\right)$. □

Up to subsequence, we may assume that there exists $u\in H_r^1\left({ℝ}^3\right)$ and $\phi \in D_r^{\mathrm{1,2}}\left({ℝ}^3\right)$ such that

$u_n\rightharpoonup u\text{in}{H}_r^1\left({ℝ}^3\right)\mathrm{,}$ (15)

$u_n\to u\text{in}{L}_r^p\left({ℝ}^3\mathrm{,}\text{d}x\right)\text{for}2<p<\mathrm{6,}$ (16)

$\Phi \left(u_n\right)\rightharpoonup \phi \text{in}{D}_r^{\mathrm{1,2}}\left({ℝ}^3\right)\mathrm{.}$ (17)

Lemma 3.7. $\phi =\Phi \left(u\right)$ and $\Phi \left(u_n\right)\to \Phi \left(u\right)$ in $D_r^{\mathrm{1,2}}\left({ℝ}^3\right)$ .

Proof. First we prove the uniqueness. For every fixed $u\in H_r^1\left({ℝ}^3\right)$, we consider the following minimizing problem $\underset{\varphi \in D_r^{\mathrm{1,2}}\left({ℝ}^3\right)}{inf}{E}_u\left(\varphi \right)$, where $E_u\mathrm{<mo>\; :\; </mo>}{D}_r^{\mathrm{1,2}}\left({ℝ}^3\right)\to ℝ$ defined as energy functional of the second equation in System (5).

$\begin{array}{l}{E}_u\left(\varphi \right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \varphi \right|}^2\text{d}x+\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\varphi }^2{u}^2\text{d}x+{\displaystyle {\int }_{{ℝ}^3}}\omega \varphi u^2\text{d}x\mathrm{.}\hfill \end{array}$

In fact, by the proof of Lemma 2.1 in [16], one can know

$\Phi \left(u_n\right)\to \phi \mathrm{,}\text{locallyuniformlyin}{ℝ}^3\mathrm{,}$

so we obtain

${\displaystyle {\int }_{{ℝ}^3}}\Phi \left(u_n\right)u_n^2\text{d}x\to {\displaystyle {\int }_{{ℝ}^3}}\phi u^2\text{d}x\mathrm{,}{\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u_n\right)u_n^2\text{d}x\to {\displaystyle {\int }_{{ℝ}^3}}{\phi }^2{u}^2\text{d}x\mathrm{.}$

From the weak lower semicontinuity of the norm in $D_r^{\mathrm{1,2}}\left({ℝ}^3\right)$ and the convergence above, one has

$E_u\left(\phi \right)\le \underset{n\to \infty }{\lim \inf }{E}_{u_n}\left(\Phi \left(u_n\right)\right)\le \underset{n\to \infty }{\lim \inf }{E}_{u_n}\left(\Phi \left(u\right)\right)=E_u\left(\Phi \left(u\right)\right),$

then by 1) of Lemma 3.1, $\phi =\Phi \left(u\right)$.

Next, we prove that $\left\{\Phi \left(u_n\right)\right\}$ converges strongly in $D_r^{\mathrm{1,2}}\left({ℝ}^3\right)$. Since $\Phi \left(u_n\right)$ satisfies the following equation

${\displaystyle {\int }_{{ℝ}^3}}\nabla \Phi \left(u_n\right)\nabla \psi \text{d}x+{\displaystyle {\int }_{{ℝ}^3}}{u}_n^2\Phi \left(u_n\right)\psi \text{d}x=-\omega {\displaystyle {\int }_{{ℝ}^3}}{u}_n^2\psi \text{d}x\mathrm{,}\psi \in D_r^{\mathrm{1,2}}\left({ℝ}^3\right)\mathrm{.}$ (18)

Let us take the difference between (18) and the corresponding equation for $\Phi \left(u\right)$ to have

$\begin{array}{l}{\displaystyle {\int }_{{ℝ}^3}}\nabla \left(\Phi \left(u_n\right)-\Phi \left(u\right)\right)\nabla \psi +u_n^2\left(\Phi \left(u_n\right)-\Phi \left(u\right)\right)\psi +\left(u_n^2-u^2\right)\Phi \left(u\right)\psi \text{d}x\\ =-\omega {\displaystyle {\int }_{{ℝ}^3}}\left(u_n^2-u^2\right)\psi \text{d}x\mathrm{,}\psi \in D_r^{\mathrm{1,2}}\left({ℝ}^3\right)\mathrm{.}\end{array}$ (19)

Testing with $\psi =\left(\Phi \left(u_n\right)-\Phi \left(u\right)\right)$, by the Hölder inequality, the following holds

$\begin{array}{l}{\Vert \Phi \left(u_n\right)-\Phi \left(u\right)\Vert }_{D_r^{\mathrm{1,2}}}^2\\ \le \left|\omega \right|{\displaystyle {\int }_{{ℝ}^3}}\left|u_n^2-u^2\right|\left|\Phi \left(u_n\right)-\Phi \left(u\right)\right|\text{d}x+{\displaystyle {\int }_{{ℝ}^3}}\left|u_n^2-u^2\right|\left|\Phi \left(u\right)\right|\left|\Phi \left(u_n\right)-\Phi \left(u\right)\right|\text{d}x\\ \le \left|\omega \right|{\Vert \Phi \left(u_n\right)-\Phi \left(u\right)\Vert }_6{\left|u_n^2-u^2\right|}_{\frac{6}{5}}+{\left|\Phi \left(u\right)\right|}_6{\left|\Phi \left(u_n\right)-\Phi \left(u\right)\right|}_6{\left|u_n^2-u^2\right|}_{\frac{3}{2}}\\ \le C_7{\left|u_n^{}-u\right|}_{\frac{12}{5}}+C_8{\left|u_n^{}-u\right|}_3\mathrm{,}\end{array}$ (20)

according to (16), one has $\Phi \left(u_n\right)\to \Phi \left(u\right)$ strongly in $D_r^{\mathrm{1,2}}\left({ℝ}^3\right)$. □

Lemma 3.8. $\left\{u_n\right\}\in H_r^1\left({ℝ}^3\right)$ has a strongly convergent subsequence in $H_r^1\left({ℝ}^3\right)$ .

Proof. Consider a sequence $\left\{u_n\right\}$ in $H_r^1\left({ℝ}^3\right)$, which satisfies $I\left(u_n\right)\to c$, $I^{\prime }\left(u_n\right)\to 0$, and $sup{\Vert u_n\Vert }_{H_r^1}<+\infty$. Going if necessary to a subsequence, since the embedding $H_r^1\left({ℝ}^3\right)$ ↪ $L_r^p\left({ℝ}^3\mathrm{,}{\left|x\right|}^{-s}\text{d}x\right)$ is compact for any $p\in \left({\mathrm{2,2}}^{\ast }\left(s\right)\right)$, we have

$u_n\to u\text{in}{L}_r^p\left({ℝ}^3\mathrm{,}{\left|x\right|}^{-s}\text{d}x\right)\mathrm{.}$ (21)

According to (10), one obtains

$\begin{array}{l}{I}^{\prime }\left(u_n\right)\left(u_n-u\right)\\ ={\displaystyle {\int }_{{ℝ}^3}}\left(\nabla u_n\cdot \nabla \left(u_n-u\right)+\left[m^2-{\left(\omega +\Phi \left(u_n\right)\right)}^2\right]u_n\left(u_n-u\right)\right)\text{d}x\\ -{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^{p-2}}{{\left|x\right|}^s}{u}_n\left(u_n-u\right)\text{d}x\mathrm{.}\end{array}$

Similarly, one gets

$\begin{array}{l}{I}^{\prime }\left(u\right)\left(u_n-u\right)\\ ={\displaystyle {\int }_{{ℝ}^3}}\left(\nabla u\cdot \nabla \left(u_n-u\right)+\left[m^2-{\left(\omega +\Phi \left(u\right)\right)}^2\right]u\left(u_n-u\right)\right)\text{d}x\\ -{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{p-2}}{{\left|x\right|}^s}u\left(u_n-u\right)\text{d}x\mathrm{.}\end{array}$

By (10), we easily get that

$\begin{array}{l}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla \left(u_n-u\right)\right|}^2+\left(m^2-{\omega }^2\right){\left(u_n-u\right)}^2\right)\text{d}x\\ =\langle I^{\prime }\left(u_n\right)-I^{\prime }\left(u\right)\mathrm{,}{u}_n-u\rangle +{\displaystyle {\int }_{{ℝ}^3}}\left({\Phi }^2\left(u_n\right)u_n-{\Phi }^2\left(u\right)u\right)\left(u_n-u\right)\text{d}x\\ +{\displaystyle {\int }_{{ℝ}^3}}\left(\frac{{\left|u_n\right|}^{p-2}}{{\left|x\right|}^s}{u}_n-\frac{{\left|u\right|}^{p-2}}{{\left|x\right|}^s}u\right)\left(u_n-u\right)\text{d}x\\ +2\omega {\displaystyle {\int }_{{ℝ}^3}}\left(\Phi \left(u_n\right)u_n-\Phi \left(u\right)u\right)\left(u_n-u\right)\text{d}x\mathrm{.}\end{array}$ (22)

It is clear that

$\begin{array}{l}\langle I^{\prime }\left(u_n\right)-I^{\prime }\left(u\right)\mathrm{,}{u}_n-u\rangle \to \mathrm{0,}\text{as}n\to \infty \mathrm{.}\hfill \end{array}$ (23)

Furthermore, in view of (21), we have

$\begin{array}{c}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^{p-1}}{{\left|x\right|}^s}\left(u_n-u\right)\text{d}x\le {\left({\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^p}{{\left|x\right|}^s}\text{d}x\right)}^{\frac{p-1}{p}}{\left({\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n-u\right|}^p}{{\left|x\right|}^s}\text{d}x\right)}^{\frac{1}{p}}\\ ={\left|u_n\right|}_{p\mathrm{,}s}^{p-1}{\left|u_n-u\right|}_{p\mathrm{,}s}\to \mathrm{0,}\text{as}n\to \infty \mathrm{.}\end{array}$ (24)

Similarly, we also obtain that

${\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{p-1}}{{\left|x\right|}^s}\left(u_n-u\right)\text{d}x\to \mathrm{0,}\text{as}n\to \infty \mathrm{.}$ (25)

Thus combining (24) and (25), one gets that

$\begin{array}{l}{\displaystyle {\int }_{{ℝ}^3}}\left(\frac{{\left|u_n\right|}^{p-2}}{{\left|x\right|}^s}{u}_n-\frac{{\left|u\right|}^{p-2}}{{\left|x\right|}^s}u\right)\left(u_n-u\right)\text{d}x\\ ={\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^{p-1}}{{\left|x\right|}^s}\left(u_n-u\right)\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{p-1}}{{\left|x\right|}^s}\left(u_n-u\right)\text{d}x\\ \to \mathrm{0,}\text{as}n\to \infty \mathrm{.}\end{array}$ (26)

By Hölder inequality and Sobolev inequality, one has

$\begin{array}{l}\left|{\displaystyle {\int }_{{ℝ}^3}}\left(\Phi \left(u_n\right)-\Phi \left(u\right)\right)u_n\left(u_n-u\right)\text{d}x\right|\\ \le {\left|\left(\Phi \left(u_n\right)-\Phi \left(u\right)\right)\left(u_n-u\right)\right|}_2{\left|u_n\right|}_2\\ \le {\left|\Phi \left(u_n\right)-\Phi \left(u\right)\right|}_6{\left|u_n-u\right|}_3{\left|u_n\right|}_2\\ \le C_9{\Vert \Phi \left(u_n\right)-\Phi \left(u\right)\Vert }_{D_r^{\mathrm{1,2}}}{\left|u_n-u\right|}_3{\Vert u_n\Vert }_{H_r^1}\mathrm{.}\end{array}$

According to (16), one gets ${\displaystyle {\int }_{{ℝ}^3}}\left(\Phi \left(u_n\right)-\Phi \left(u\right)\right)u_n\left(u_n-u\right)\text{d}x\to 0$, as $n\to \infty$.

And

$\begin{array}{l}\left|{\displaystyle {\int }_{{ℝ}^3}}\Phi \left(u\right){\left(u_n-u\right)}^2\text{d}x\right|\le {\left|\Phi \left(u\right)\right|}_6{\left|u_n-u\right|}_3{\left|u_n-u\right|}_2\\ \le C_{10}{\Vert \Phi \left(u\right)\Vert }_{D_r^{\mathrm{1,2}}}{\left|u_n-u\right|}_3{\Vert u_n-u\Vert }_{H_r^1}\to \mathrm{0,}\text{as}n\to \infty \mathrm{.}\end{array}$

Thus we get that

$\begin{array}{l}{\displaystyle {\int }_{{ℝ}^3}}\left(\Phi \left(u_n\right)u_n-\Phi \left(u\right)u\right)\left(u_n-u\right)\text{d}x\\ ={\displaystyle {\int }_{{ℝ}^3}}\left(\Phi \left(u_n\right)-\Phi \left(u\right)\right)u_n\left(u_n-u\right)\text{d}x+{\displaystyle {\int }_{{ℝ}^3}}\Phi \left(u\right){\left(u_n-u\right)}^2\text{d}x\\ \to \mathrm{0,}\text{as}n\to \infty \mathrm{.}\end{array}$ (27)

We observe that the sequence $\left\{{\Phi }^2\left(u_n\right)u_n\right\}$ is bounded in $L^{\frac{3}{2}}\left({ℝ}^3\right)$, since

${\left|{\Phi }^2\left(u_n\right)u_n\right|}_{\frac{3}{2}}\le {\left|\Phi \left(u_n\right)\right|}_6^2{\left|u_n\right|}_3\mathrm{,}$

so that

$\begin{array}{l}\left|{\displaystyle {\int }_{{ℝ}^3}}\left({\Phi }^2\left(u_n\right)u_n-{\Phi }^2\left(u\right)u\right)\left(u_n-u\right)\text{d}x\right|\\ \le {\left|{\Phi }^2\left(u_n\right)u_n-{\Phi }^2\left(u\right)u\right|}_{\frac{3}{2}}{\left|u_n-u\right|}_3\\ \le \left({\left|{\Phi }^2\left(u_n\right)u_n\right|}_{\frac{3}{2}}+{\left|{\Phi }^2\left(u\right)u\right|}_{\frac{3}{2}}\right){\left|u_n-u\right|}_3\mathrm{.}\end{array}$

By (16), one has

${\displaystyle {\int }_{{ℝ}^3}}\left({\Phi }^2\left(u_n\right)u_n-{\Phi }^2\left(u\right)u\right)\left(u_n-u\right)\text{d}x\to \mathrm{0,}\text{as}n\to \infty \mathrm{.}$ (28)

Therefore according to (23)-(28) and $\left|m\right|>\left|\omega \right|$, we obtain that

$\begin{array}{c}{\Vert u_n-u\Vert }_{H_r^1}\le max\left\{\mathrm{1,}{m}^2-{\omega }^2\right\}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla \left(u_n-u\right)\right|}^2+{\left(u_n-u\right)}^2\right)\text{d}x\\ \to \mathrm{0,}\text{as}n\to \infty \mathrm{.}\end{array}$

Thus $\left\{u_n\right\}$ has a strongly convergent subsequence in $H_r^1\left({ℝ}^3\right)$.

Consequently, we conclude that

$I\left(u\right)=c,I^{\prime }\left(u\right)=0.$ □

Next we begin to prove Theorem 2.1.

Proof. We only need to prove that $u\ne 0$. Suppose by contradiction that $u=0$, and hence $\Phi \left(u\right)=0$. Since as $n\to \infty$, $\langle I^{\prime }\left(u_n\right)\mathrm{,}{u}_n\rangle \to 0$, $u_n\to 0$ in $L_r^p\left({ℝ}^3\mathrm{,}{\left|x\right|}^{-s}\text{d}x\right)\left(2<p<2^{\ast }\left(s\right)\right)$ and $L_r^p\left({ℝ}^3\mathrm{,}\text{d}x\right)\left(2<p<2^{\ast }\right)$. Thus we get

${\displaystyle {\int }_{{ℝ}^3}}{u}_n^2{\Phi }^2\left(u_n\right)\text{d}x+2\omega {\displaystyle {\int }_{{ℝ}^3}}\Phi \left(u_n\right)u_n^2\text{d}x\to \mathrm{0,}$

${\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^p}{{\left|x\right|}^s}\text{d}x\to 0.$

We may assume

${\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u_n\right|}^2+\left(m^2-{\omega }^2\right)u_n^2\right)\text{d}x=\tilde{C}{\Vert u_n\Vert }_{H_r^1}^2\mathrm{,}\tilde{C}>0.$

Set

${\Vert u_n\Vert }_{H_r^1}^2\to l\mathrm{,}l\ge \mathrm{0,}$

obviously, $l=0\iff u_n\to 0$ in $H_r^1\left({ℝ}^3\right)$. As a consequence we obtain that

$I\left(u_n\right)\to \frac{\tilde{C}l}{2}\mathrm{,}n\to \infty \mathrm{.}$

According to $I\left(u_n\right)\to c>0$, we get

$c=\frac{\tilde{C}l}{2}>0,$

which implies that $l=0$ is impossible, i.e., which contradicts with $u=0$. Therefore, u is a nontrivial solution of System (5).

This theorem is mainly based on satisfying the conditions of the Mountain Pass Theorem, and then there is a (PS) sequence, proving that the (PS) sequence is bounded and $u_n$ has a strongly convergent sub-sequence in $H_r^1\left({ℝ}^3\right)$, so as to prove that the system of Equation (5) has at least a nontrivial solution. □

4. The Proof of Theorem 2.2

Similarly, we also have the following lemmas.

Lemma 4.1. For every $u\in H^1\left({ℝ}^3\right)$ ,

1) there exists a unique function $\varphi =\Phi \left(u\right)\in D^{\mathrm{1,2}}\left({ℝ}^3\right)$ that solves the second equation of System (6);

2) if u is radially symmetric, then $\Phi \left(u\right)$ is radial too;

3) $\Phi \left(u\right)\left(x\right)\le 0$, moreover, $\Phi \left(u\right)\left(x\right)\ge -\omega$, if $u\left(x\right)\ne 0$ and $\omega >0$.

Lemma 4.2. The map $\Phi$ is $C^1$ and ${\tilde{G}}_{\varphi }=\left\{\left(u\mathrm{,}\varphi \right)\in H^1\left({ℝ}^3\right)\times D^{\mathrm{1,2}}\left({ℝ}^3\right)\mathrm{<mo>\; |\; </mo>}{{\tilde{F}}^{\prime }}_{\varphi }\left(u\mathrm{,}\varphi \right)=0\right\}$.

Likewise, define $J\left(u\right)=\tilde{F}\left(u\mathrm{,}\varphi \right)$.

Lemma 4.3. The pair $\left(u\mathrm{,}\varphi \right)$ is a weak solution of System (6) if and only if it is a critical point of J in $H^1\left({ℝ}^3\right)\times D^{\mathrm{1,2}}\left({ℝ}^3\right)$ .

Lemma 4.4 For $u\in H^1\left({ℝ}^3\right)$ , if $\left|m\right|>\left|\omega \right|$ , then there exist some constants ${\rho }_2\mathrm{,}{\alpha }_2\mathrm{,}\tilde{m}>0$ such that ${J\left(u\right)|}_{{\Vert u\Vert }_{H^1}={\rho }_2}\ge {\alpha }_2$ for all $\mu$ satisfying $0<\mu <\tilde{m}$ .

Proof. From (8), one obtains

$\begin{array}{c}J\left(u\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2\right)\text{d}x-\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \varphi \right|}^2\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}\omega \varphi u^2\text{d}x\\ -\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\varphi }^2{u}^2\text{d}x-\frac{\mu }{q}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^q}{{\left|x\right|}^s}\text{d}x-\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x-\frac{1}{2^{\ast }\left(s\right)}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\mathrm{.}\end{array}$ (29)

Substituting (9) into (29), we have

$\begin{array}{c}J\left(u\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2\right)\text{d}x+\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \varphi \right|}^2\text{d}x+\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\varphi }^2{u}^2\text{d}x\\ -\frac{\mu }{q}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^q}{{\left|x\right|}^s}\text{d}x-\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x-\frac{1}{2^{\ast }\left(s\right)}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\\ \ge \frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2\right)\text{d}x-\frac{\mu }{q}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^q}{\mathrm{<mo>\; |\; </mo>}x{\mathrm{<mo>\; |\; </mo>}}^s}\text{d}x\\ -\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x-\frac{1}{2^{\ast }\left(s\right)}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\end{array}$

$\begin{array}{l}\ge \frac{min\left\{\mathrm{1,}{m}^2-{\omega }^2\right\}}{2}{\Vert u\Vert }_{H^1}^2-\frac{C_{11}\mu }{q}{\Vert u\Vert }_{H^1}^q-\frac{C_{12}}{p}{\Vert u\Vert }_{H^1}^p-\frac{1}{2^{\ast }\left(s\right)}{S}_s^{-\frac{2^{\ast }\left(s\right)}{2}}{\Vert u\Vert }_{H^1}^{2^{\ast }\left(s\right)}\\ ={\Vert u\Vert }_{H^1}^q(\frac{min\left\{\mathrm{1,}{m}^2-{\omega }^2\right\}}{2}{\Vert u\Vert }_{H^1}^{2-q}-\frac{C_{12}}{p}{\Vert u\Vert }_{H^1}^{p-q}-\frac{1}{2^{\ast }\left(s\right)}{S}_s^{-\frac{2^{\ast }\left(s\right)}{2}}{\Vert u\Vert }_{H^1}^{2^{\ast }\left(s\right)-q}-\frac{C_{11}\mu }{q})\mathrm{.}\end{array}$

Set $h_1\left(\rho \right)=\frac{\min \left\{1,m^2-{\omega }^2\right\}}{2}{\rho }^{2-q}-\frac{C_{12}}{p}{\rho }^{p-q}-\frac{1}{2^{\ast }\left(s\right)}{S}_s^{-\frac{2^{\ast }\left(s\right)}{2}}{\rho }^{2^{\ast }\left(s\right)-q},\rho \ge 0$.

Since

$\begin{array}{c}{h^{\prime }}_1\left(\rho \right)=\frac{\left(min\left\{\mathrm{1,}{m}^2-{\omega }^2\right\}\right)\left(2-q\right)}{2}{\rho }^{1-q}-\frac{C_{12}\left(p-q\right)}{p}{\rho }^{p-q-1}\\ -\frac{2^{\ast }\left(s\right)-q}{2^{\ast }\left(s\right)}{S}_s^{-\frac{2^{\ast }\left(s\right)}{2}}{\rho }^{2^{\ast }\left(s\right)-q-1}\\ ={\rho }^{1-q}(\frac{\left(min\left\{\mathrm{1,}{m}^2-{\omega }^2\right\}\right)\left(2-q\right)}{2}-\frac{C_{12}\left(p-q\right)}{p}{\rho }^{p-2}\\ -\frac{2^{\ast }\left(s\right)-q}{2^{\ast }\left(s\right)}{S}_s^{-\frac{2^{\ast }\left(s\right)}{2}}{\rho }^{2^{\ast }\left(s\right)-2})\mathrm{.}\end{array}$

Evidently, when $\rho \ge 0$ is small enough, ${h^{\prime }}_1\left(\rho \right)$ is greater than 0, and $h_1\left(\rho \right)$ increases monotonically. When $\rho \ge 0$ is large enough, ${h^{\prime }}_1\left(\rho \right)$ is less than 0, and $h_1\left(\rho \right)$ decreases monotonically. Therefore there exists a maximum point $\tilde{\xi }$ such that

${h^{\prime }}_1\left(\tilde{\xi }\right)=0.$

Obviously, we observe that $\tilde{\xi }=0$ or $\tilde{\xi }>0$, while $\tilde{\xi }=0$ is impossible. Then $\tilde{\xi }>0$, next we obtain $\underset{\rho \ge 0}{\max }{h}_1\left(\rho \right)=h_1\left(\tilde{\xi }\right)>0$. Choosing that ${\Vert u\Vert }_{H^1}=\tilde{\xi }={\rho }_2$, we deduce, for all $\mu$ satisfying $0<\mu <\frac{q{h}_1\left({\rho }_2\right)}{2C_{11}}=\tilde{m}$,

${J\left(u\right)|}_{{\Vert u\Vert }_{H^1}={\rho }_2}\ge {\rho }_2^q\cdot \left(h_1\left({\rho }_2\right)-\frac{C_{11}\mu }{q}\right)\ge \frac{{\rho }_2^q\cdot h_1\left({\rho }_2\right)}{2}={\alpha }_2>0$

and the proof is completed. □

Lemma 4.5. Under the assumptions of Theorem 2.2, there exists a function ${\eta }_2\in H^1\left({ℝ}^3\right)$ with ${\Vert {\eta }_2\Vert }_{H^1}>{\rho }_2$ such that $J\left({\eta }_2\right)<0$ .

Similarly, a standard tool is to restrict ourselves to radial functions, namely $H_r^1\left({ℝ}^3\right)=\left\{u\in H^1\left({ℝ}^3\right)\mathrm{<mo>\; |\; </mo>}u\left(x\right)=u\left(\left|x\right|\right)\right\}$ and $D_r^{\mathrm{1,2}}\left({ℝ}^3\right)=\left\{u\in D^{\mathrm{1,2}}\left({ℝ}^3\right)\mathrm{<mo>\; |\; </mo>}u\left(x\right)=u\left(\left|x\right|\right)\right\}$ compactly embedded in $L_r^p\left({ℝ}^3\mathrm{,}{\left|x\right|}^{-s}\text{d}x\right)$ for $2<p<2^{\ast }\left(s\right)$ and $L_r^p\left({ℝ}^3\mathrm{,}\text{d}x\right)$ for $2<p<2^{\ast }=\frac{2N}{N-2}=6$. Moreover, one sees that if a critical point $u\in H_r^1\left({ℝ}^3\right)$ for the functional ${J|}_{H_r^1\left({ℝ}^3\right)}$ is also a critical point of J.

Lemma 4.6. Under the assumptions of Theorem 2.2, if $\left\{u_n\right\}\in H_r^1\left({ℝ}^3\right)$ is a bounded Palais-Smale sequence of J, then $\left\{u_n\right\}$ has a strongly convergent subsequence in $H_r^1\left({ℝ}^3\right)$ .

Proof. Consider a sequence $\left\{u_n\right\}$ in $H_r^1\left({ℝ}^3\right)$, which satisfies $J\left(u_n\right)\to c$, $J^{\prime }\left(u_n\right)\to 0$, and $\sup {\Vert u_n\Vert }_{H_r^1}<+\infty$. Going if necessary to a subsequence, we assume

$u_n\rightharpoonup u\text{in}{H}_r^1\left({ℝ}^3\right)\mathrm{.}$

Since the embedding $H_r^1\left({ℝ}^3\right)$ ↪ $L_r^p\left({ℝ}^3\mathrm{,}{\left|x\right|}^{-s}\text{d}x\right)$ is compact for any $p\in \left({\mathrm{2,2}}^{\ast }\left(s\right)\right)$, we have

$u_n\to u\text{in}{L}_r^p\left({ℝ}^3\mathrm{,}{\left|x\right|}^{-s}\text{d}x\right)\mathrm{.}$ (30)

Moreover, likewise, for $p\in \left({\mathrm{2,2}}^{\ast }\right)$, we also get

$u_n\to u\text{in}{L}_r^p\left({ℝ}^3\mathrm{,}\text{d}x\right)\mathrm{.}$ (31)

According to

$\begin{array}{c}{J}^{\prime }\left(u\right)v={\displaystyle {\int }_{{ℝ}^3}}\left(\nabla u\cdot \nabla v+\left[m^2-{\left(\omega +\varphi \right)}^2\right]uv\right)\text{d}x-\mu {\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{q-2}}{{\left|x\right|}^s}uv\text{d}x\\ -{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{p-2}}{{\left|x\right|}^s}uv\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{2^{\ast }\left(s\right)-2}}{{\left|x\right|}^s}uv\text{d}x\mathrm{,}\end{array}$ (32)

one obtains

$\begin{array}{l}{J}^{\prime }\left(u_n\right)\left(u_n-u\right)\\ ={\displaystyle {\int }_{{ℝ}^3}}\left(\nabla u_n\cdot \nabla \left(u_n-u\right)+\left[m^2-{\left(\omega +\Phi \left(u_n\right)\right)}^2\right]u_n\left(u_n-u\right)\right)\text{d}x\\ -\mu {\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^{q-2}}{{\left|x\right|}^s}{u}_n\left(u_n-u\right)\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^{p-2}}{{\left|x\right|}^s}{u}_n\left(u_n-u\right)\text{d}x\\ -{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^{2^{\ast }\left(s\right)-2}}{{\left|x\right|}^s}{u}_n\left(u_n-u\right)\text{d}x\mathrm{.}\end{array}$

Similarly, one gets

$\begin{array}{l}{J}^{\prime }\left(u\right)\left(u_n-u\right)\\ ={\displaystyle {\int }_{{ℝ}^3}}\left(\nabla u\cdot \nabla \left(u_n-u\right)+\left[m^2-{\left(\omega +\Phi \left(u\right)\right)}^2\right]u\left(u_n-u\right)\right)\text{d}x\\ -\mu {\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{q-2}}{{\left|x\right|}^s}u\left(u_n-u\right)\text{d}x-{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{p-2}}{{\left|x\right|}^s}u\left(u_n-u\right)\text{d}x\\ -{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{2^{\ast }\left(s\right)-2}}{{\left|x\right|}^s}u\left(u_n-u\right)\text{d}x\mathrm{.}\end{array}$

By (32), we easily get that

$\begin{array}{l}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla \left(u_n-u\right)\right|}^2+\left(m^2-{\omega }^2\right){\left(u_n-u\right)}^2\right)\text{d}x\\ =\langle J^{\prime }\left(u_n\right)-J^{\prime }\left(u\right)\mathrm{,}{u}_n-u\rangle +\mu {\displaystyle {\int }_{{ℝ}^3}}\left(\frac{{\left|u_n\right|}^{q-2}}{{\left|x\right|}^s}{u}_n-\frac{{\left|u\right|}^{q-2}}{{\left|x\right|}^s}u\right)\left(u_n-u\right)\text{d}x\\ +{\displaystyle {\int }_{{ℝ}^3}}\left(\frac{{\left|u_n\right|}^{p-2}}{{\left|x\right|}^s}{u}_n-\frac{{\left|u\right|}^{p-2}}{{\left|x\right|}^s}u\right)\left(u_n-u\right)\text{d}x\\ +2\omega {\displaystyle {\int }_{{ℝ}^3}}\left(\Phi \left(u_n\right)u_n-\Phi \left(u\right)u\right)\left(u_n-u\right)\text{d}x\\ +{\displaystyle {\int }_{{ℝ}^3}}\left(\frac{{\left|u_n\right|}^{2^{\ast }\left(s\right)-2}}{{\left|x\right|}^s}{u}_n-\frac{{\left|u\right|}^{2^{\ast }\left(s\right)-2}}{{\left|x\right|}^s}u\right)\left(u_n-u\right)\text{d}x\\ +{\displaystyle {\int }_{{ℝ}^3}}\left({\Phi }^2\left(u_n\right)u_n-{\Phi }^2\left(u\right)u\right)\left(u_n-u\right)\text{d}x\mathrm{.}\end{array}$ (33)

Then we use similar method in the proof of Lemma 3.8, we obtain that $\left\{u_n\right\}$ has a strongly convergent subsequence in $H_r^1\left({ℝ}^3\right)$. □

Next we begin to prove Theorem 2.2.

Proof. The proof is divided into two steps.

Step 1. There exists $u_1\in H_r^1\left({ℝ}^3\right)$ such that $J^{\prime }\left(u_1\right)=0$ and $J\left(u_1\right)<0$.

We choose a function $v\in H_r^1\left({ℝ}^3\right)$. Since $\Phi \left(u\right)\ge -\omega$, one has

$\begin{array}{c}J\left(tv\right)=\frac{t^2}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla v\right|}^2+\left(m^2-{\omega }^2\right)v^2\right)\text{d}x-t^2{\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(tv\right)v^2\text{d}x\\ -\frac{t^2}{2}{\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(tv\right)v^2\text{d}x-\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(tv\right)\right|}^2\text{d}x\\ -\frac{\mu t^q}{q}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|v\right|}^q}{{\left|x\right|}^s}\text{d}x-\frac{t^p}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|v\right|}^p}{{\left|x\right|}^s}\text{d}x-\frac{t^{2^{\ast }\left(s\right)}}{2^{\ast }\left(s\right)}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|v\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\\ \le \frac{t^2}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla v\right|}^2+\left(m^2-{\omega }^2\right)v^2+2{\omega }^2{v}^2\right)\text{d}x-\frac{\mu t^q}{q}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|v\right|}^q}{{\left|x\right|}^s}\text{d}x\\ <0\end{array}$ (34)

for $t>0$ small enough. Thus we have $c_1=\inf \left\{J\left(u\right):u\in {\bar{B}}_{{\rho }_2}\right\}<0$, where ${\rho }_2>0$ is given by Lemma 4.4, $B_{{\rho }_2}=\left\{u\in H_r^1\left({ℝ}^3\right):{\Vert u\Vert }_{H_r^1}<{\rho }_2\right\}$. By the Ekeland’s variational principle [17], Let $\epsilon =\frac{1}{n}>0$, $\delta ={\Vert \zeta -u_n\Vert }_{H_r^1}>0$ for all $\zeta \in {\bar{B}}_{{\rho }_2}$, then there exists a sequence $\left\{u_n\right\}\subset {\bar{B}}_{{\rho }_2}$ such that

$c_1\le J\left(u_n\right)\le c_1+\epsilon$

and

$c_1+\frac{\epsilon }{\delta }{\Vert \zeta -u_n\Vert }_{H_r^1}\le J\left(\zeta \right)+\frac{\epsilon }{\delta }{\Vert \zeta -u_n\Vert }_{H_r^1}\mathrm{.}$

Then we obtain that

$J\left(\zeta \right)\ge J\left(u_n\right)-\frac{\epsilon }{\delta }{\Vert \zeta -u_n\Vert }_{H_r^1}\mathrm{.}$

Obviously, in view of Lemma 4.4, $u_n\in B_{{\rho }_2}$, for n large enough. Thus for any $\varphi \in H_r^1\left({ℝ}^3\right)$ with ${\Vert \varphi \Vert }_{H_r^1}=1$, we can take $t>0$ such that $\left(u_n+t\varphi \right)\in {\bar{B}}_{{\rho }_2}$ for n large enough (see [18] ). Then we have

$\frac{J\left(u_n+t\varphi \right)-J\left(u_n\right)}{t}\ge -\epsilon \mathrm{.}$

Letting $t\to 0$, we get

$\langle J^{\prime }\left(u_n\right)\mathrm{,}\varphi \rangle \ge -\epsilon \mathrm{.}$

We replace $\varphi$ by $-\varphi$ in the above inequality, then it follows that

$\langle J^{\prime }\left(u_n\right)\mathrm{,}-\varphi \rangle \ge -\epsilon \mathrm{,}$

i.e.,

$\langle J^{\prime }\left(u_n\right)\mathrm{,}\varphi \rangle \le \epsilon \mathrm{.}$

Thus one obtains

$\left|\langle J^{\prime }\left(u_n\right)\mathrm{,}\varphi \rangle \right|\le \epsilon \mathrm{,}$

which implies $J^{\prime }\left(u_n\right)\to 0$ as $n\to \infty$.

Hence we conclude that $\left\{u_n\right\}$ is bounded PS sequence of J for $c_1$. Therefore, by Lemma 4.6, we get that there exists a function $u_1\in H_r^1\left({ℝ}^3\right)$ such that $J^{\prime }\left(u_1\right)=0$ and $J\left(u_1\right)=c_1<0$.

Step 2. There exists $u_2\in H_r^1\left({ℝ}^3\right)$ such that $J^{\prime }\left(u_2\right)=0$ and $J\left(u_2\right)>0$.

From Lemma 4.4, Lemma 4.5 and the Mountain Pass Theorem [19], there is a sequence $\left\{u_n\right\}\subset H_r^1\left({ℝ}^3\right)$ such that

$J\left(u_n\right)\to c_2=\underset{\gamma \in \Gamma }{\inf }\underset{t\in \left[0,1\right]}{\max }J\left(\gamma \left(t\right)\right)\ge {\alpha }_2>0$

and

$J^{\prime }\left(u_n\right)\to \mathrm{0,}$

where $\Gamma =\left\{\gamma \in C\left(\left[0,1\right],H_r^1\left({ℝ}^3\right)\right):\gamma \left(0\right)=0,\gamma \left(1\right)={\eta }_2\right\}$. From Lemma 4.6, we only need to prove that $\left\{u_n\right\}$ is bounded in $H_r^1\left({ℝ}^3\right)$.

Case 1): $4\le p<2^{\ast }\left(s\right)$. From (8) and (9), one has

$\begin{array}{c}J\left(u\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2\right)\text{d}x+\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u\right)u^2\text{d}x+\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(u\right)\right|}^2\text{d}x\\ -\frac{\mu }{q}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^q}{{\left|x\right|}^s}\text{d}x-\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x-\frac{1}{2^{\ast }\left(s\right)}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\mathrm{.}\end{array}$

Then by (32), one obtains

$\begin{array}{l}\left(c_2+1\right)+o\left(1\right){\Vert u_n\Vert }_{H_r^1}\\ \ge J\left(u_n\right)-\frac{1}{p}\langle J^{\prime }\left(u_n\right)\mathrm{,}{u}_n\rangle \\ =\left(\frac{1}{2}-\frac{1}{p}\right){\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u_n\right|}^2+\left(m^2-{\omega }^2\right)u_n^2\right)\text{d}x+\left(\frac{1}{2}+\frac{1}{p}\right){\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u_n\right)u_n^2\text{d}x\\ +\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(u_n\right)\right|}^2\text{d}x+\mu \left(\frac{1}{p}-\frac{1}{q}\right){\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^q}{{\left|x\right|}^s}\text{d}x\\ +\frac{2}{p}{\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(u_n\right)u_n^2\text{d}x+\left(\frac{1}{p}-\frac{1}{2^{\ast }\left(s\right)}\right){\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\end{array}$

$\begin{array}{l}=\left(\frac{1}{2}-\frac{1}{p}\right){\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u_n\right|}^2+\left(m^2-{\omega }^2\right)u_n^2\right)\text{d}x+\left(\frac{1}{2}-\frac{1}{p}\right){\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u_n\right)u_n^2\text{d}x\\ +\left(\frac{1}{2}-\frac{2}{p}\right){\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla \Phi \left(u_n\right)\right|}^2\text{d}x+\mu \left(\frac{1}{p}-\frac{1}{q}\right){\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^q}{{\left|x\right|}^s}\text{d}x\\ +\left(\frac{1}{p}-\frac{1}{2^{\ast }\left(s\right)}\right){\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\\ \ge C_{13}{\Vert u_n\Vert }_{H_r^1}^2+\mu \left(\frac{1}{p}-\frac{1}{q}\right)C_{14}{\Vert u_n\Vert }_{H_r^1}^q\end{array}$

for n large enough. Therefore, it follows that $\left\{u_n\right\}$ is bounded in $H_r^1\left({ℝ}^3\right)$.

Case 2): $2<p<4$. From (8) and (9), one has

$\begin{array}{c}J\left(u\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u\right|}^2+\left(m^2-{\omega }^2\right)u^2\right)\text{d}x-\frac{1}{2}{\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(u\right)u^2\text{d}x-\frac{\mu }{q}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^q}{{\left|x\right|}^s}\text{d}x\\ -\frac{1}{p}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^p}{{\left|x\right|}^s}\text{d}x-\frac{1}{2^{\ast }\left(s\right)}{\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\mathrm{.}\end{array}$

Then by (32) and $\left(\frac{p}{2}-1\right)m^2>{\omega }^2$ and $\Phi \left(u_n\right)\ge -\omega$, $1<q<2$, we get

$\begin{array}{l}p\left(c_2+1\right)+{\Vert u_n\Vert }_{H_r^1}\\ \ge pJ\left(u_n\right)-\langle J^{\prime }\left(u_n\right)\mathrm{,}{u}_n\rangle \\ =\left(\frac{p}{2}-1\right){\displaystyle {\int }_{{ℝ}^3}}\left({\left|\nabla u_n\right|}^2+\left(m^2-{\omega }^2\right)u_n^2\right)\text{d}x-\left(\frac{p}{2}-2\right){\displaystyle {\int }_{{ℝ}^3}}\omega \Phi \left(u_n\right)u_n^2\text{d}x\\ +{\displaystyle {\int }_{{ℝ}^3}}{\Phi }^2\left(u_n\right)u_n^2\text{d}x+\mu \left(1-\frac{p}{q}\right){\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^q}{{\left|x\right|}^s}\text{d}x+\left(1-\frac{p}{2^{\ast }\left(s\right)}\right){\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^{2^{\ast }\left(s\right)}}{{\left|x\right|}^s}\text{d}x\end{array}$

$\begin{array}{l}\ge \left(\frac{p}{2}-1\right){\displaystyle {\int }_{{ℝ}^3}}{\left|\nabla u_n\right|}^2\text{d}x+{\displaystyle {\int }_{{ℝ}^3}}\left[\left(\frac{p}{2}-1\right)\left(m^2-{\omega }^2\right)+\left(\frac{p}{2}-2\right){\omega }^2\right]u_n^2\text{d}x\\ +\mu \left(1-\frac{p}{q}\right){\displaystyle {\int }_{{ℝ}^3}}\frac{{\left|u_n\right|}^q}{{\left|x\right|}^s}\text{d}x\\ \ge C_{15}{\Vert u_n\Vert }_{H_r^1}^2+\mu \left(1-\frac{p}{q}\right)C_{16}{\Vert u_n\Vert }_{H_r^1}^q\end{array}$

for n large enough. It follows that $\left\{u_n\right\}$ is bounded in $H_r^1\left({ℝ}^3\right)$.

This theorem is mainly based on the Mountain Pass Theorem and the Ekeland variational principle to prove that the System (6) has at least two different solutions. □

Remark 4.1. Here a simple infimum definition is given as follows. Suppose S is a set of numbers in $ℝ$ . If the number ${\xi }_0$ satisfies:

1) $\forall x\in S\mathrm{,}x\ge {\xi }_0$, that is, ${\xi }_0$ is a lower bound of S.

2) $\forall {\beta }_0>{\xi }_0$, there exists a $x_0\in S$ that satisfies $x_0<{\beta }_0$, that is, ${\xi }_0$ is maximum lower bound of S, then the number ${\xi }_0$ is the infimum of the number set. Referred to as: ${\xi }_0=infS$.

Conflicts of Interest

The author declares no conflicts of interest.

References

[1]	Bertoin, J. (1996) Lévy Processes. Volume 121 of Cambridge Tracts in Mathematics. Cambridge University Press, Cambridge.
[2]	Bucur, C. and Valdinoci, E. (2016) Nonlocal Diffusion and Applications. Volume 20 of Lecture Notes of the Unione Matematica Italiana. Springer, Cham. https://doi.org/10.1007/978-3-319-28739-3
[3]	Bhakta, M. (2015) Semilinear Elliptic Equation with Biharmonic Operator and Multiple Critical Nonlinearities. Advanced Nonlinear Studies, 15, 835-848. https://doi.org/10.1515/ans-2015-0405
[4]	Filippucci, R., Pucci, P. and Robert, F. (2009) On a p-Laplace Equation with Multiple Critical Nonlinearities. Journal de Mathématiques Pures et Appliquées, 91, 156-177. https://doi.org/10.1016/j.matpur.2008.09.008
[5]	Ghoussoub, N. and Robert, F. (2006) Concentration Estimates for Emden-Fowler Equations with Boundary Singularities and Critical Growth. International Mathematics Research Papers, 2006, Article ID: 21867.
[6]	Brézis, H. and Nirenberg, L. (1983) Positive Solutions of Nonlinear Elliptic Equations Involving Critical Sobolev Exponents. Communications on Pure and Applied Mathematics, 36, 437-477. https://doi.org/10.1002/cpa.3160360405
[7]	Jannelli, E. (1999) The Role Played by Space Dimension in Elliptic Critical Problems. Journal of Differential Equations, 156, 407-426. https://doi.org/10.1006/jdeq.1998.3589
[8]	Kang, D. and Peng, S. (2004) Existence of Solutions for Elliptic Equations with Critical Sobolev-Hardy Exponents. Nonlinear Analysis, 56, 1151-1164. https://doi.org/10.1016/j.na.2003.11.008
[9]	Kang, D. and Peng, S. (2005) Solutions for Semilinear Elliptic Problems with Critical Sobolev-Hardy Exponents and Hardy Potential. Applied Mathematics Letters, 18, 1094-1100. https://doi.org/10.1016/j.aml.2004.09.016
[10]	Kang, D. and Peng, S. (2004) Sign-Changing Solutions for Elliptic Problems with Critical Sobolev-Hardy Exponents. Journal of Mathematical Analysis and Applications, 291, 488-499. https://doi.org/10.1016/j.jmaa.2003.11.012
[11]	Benci, V. and Fortunato, D. (2002) Solitary Waves of the Nonlinear Klein-Gordon Equation Coupled with the Maxwell Equations. Reviews in Mathematical Physics, 14, 409-420. https://doi.org/10.1142/S0129055X02001168
[12]	Berestycki, H. and Lions, P.-L. (1983) Nonlinear Scalar Field Equations, I—Existence of a Ground State. Archive for Rational Mechanics and Analysis, 82, 313-345. https://doi.org/10.1007/BF00250555
[13]	Cassani, D. (2004) Existence and Non-Existence of Solitary Waves for the Critical Klein-Gordon Equation Coupled with Maxwell’s Equations. Nonlinear Analysis, 58, 733-747. https://doi.org/10.1016/j.na.2003.05.001
[14]	d’Avenia, P. and Pisani, L. (2002) Nonlinear Klein-Gordon Equations Coupled with Born-Infeld Type Equations. Electronic Journal of Differential Equations, 2002, 1-13.
[15]	Mugnai, D. (2004) Coupled Klein-Gordon and Born-Infeld-Type Equations: Looking for Solitary Waves. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 460, 1519-1527. https://doi.org/10.1098/rspa.2003.1267
[16]	Yu, Y. (2010) Solitary Waves for Nonlinear Klein-Gordon Equations Coupled with Born-Infeld Theory. Annales de l’Institut Henri Poincaré, Analyse Non Linéaire, 27, 351-376. https://doi.org/10.1016/j.anihpc.2009.11.001
[17]	Ekeland, I. (1974) On the Variational Principle. Journal of Mathematical Analysis and Applications, 47, 324-353. https://doi.org/10.1016/0022-247X(74)90025-0
[18]	Chen, S.-J. and Song, S.-Z. (2017) The Existence of Multiple Solutions for the Klein-Gordon Equation with Concave and Convex Nonlinearities Coupled with Born-Infeld Theory on R3. Nonlinear Analysis: Real World Applications, 38, 78-95. https://doi.org/10.1016/j.nonrwa.2017.04.008
[19]	Willem, M. (1996) Minimax Theorems (Progress in Nonlinear Differential Equations and Their Applications). Vol. 24, Birkhäuser Boston, Inc., Boston.