Positive Solutions for a Class of Quasilinear Schrödinger Equations with Nonlocal Term

Peng Liao, Rui Ping, Shaoxiong Chen^*

School of Mathematics, Yunnan Normal University, Kunming, China.

DOI: 10.4236/jamp.2022.102027                

Abstract

This paper is considered the existence of positive solutions for a class of generalized quasilinear Schrödinger equations with nonlocal term in R^N which have appeared from plasma physics, as well as high-power ultrashort laser in matter. We use a charge of variables and obtain the existence of solutions via minimization argument.

Keywords

Quasilinear Schrödinger Equation, Minimization, Implicit Function Theorem

1. Introduction

In this paper, we consider investigating the existence of solutions for the following generalized quasilinear Schrödinger equation with nonlocal term

$\begin{array}{l}-div\left(g^2\left(u\right)\nabla u\right)+g\left(u\right)g^{\prime }\left(u\right){\left|\nabla u\right|}^2+V\left(x\right)u\\ =\lambda \left[{\left|x\right|}^{-\mu }\ast {\left|u\right|}^p\right]{\left|u\right|}^{p-2}u+\beta {\left|u\right|}^{q-1}u\mathrm{,}\end{array}$ (1.1)

where $N\ge 3$ , $0<\mu <N$ , $\beta <\frac{N-2}{2N}$ , $1\le q\le \frac{N}{N-2}$ , $\frac{2N-\mu }{N}\le p<\frac{2N-\mu }{N-2}$ , the function $V\in C\left({ℝ}^N\mathrm{,}{ℝ}^{+}\right)$ , g is a ${\mathcal{C}}^1$ function with $g^{\prime }\left(t\right)\le 0$ for all $t>0$ , $g\left(0\right)=0$ , ${\lim }_{t\to +\infty }g\left(t\right)=a$ .

When $g\left(u\right)=1$ , (1.1) boils down to the so called nonlinear Choquard or Choquard-Pekar equation

$-\Delta u+V\left(x\right)u=\lambda \left[{\left|x\right|}^{-\mu }\ast {\left|u\right|}^p\right]{\left|u\right|}^{p-2}u+\beta {\left|u\right|}^{q-1}u$ (1.2)

Such like equation has several physical origins. The problem

$-\Delta u+u=\left[{\left|x\right|}^{-1}\ast {\left|u\right|}^2\right]u\mathrm{,}$ (1.3)

appeared at least as early as in 1954, in a work by Pekar describing the quantum mechanics of a polaron at rest [1]. In 1976, Choquard used (1.3) to describe an electron trapped in its own hole and in a certain approximation to Hartree-Fock theory of one component plasma [2]. In 1996, Penrose proposed (1.3) as a model of self-gravitating matter, in a program in which quantum state reduction is understood as a gravitational phenomenon [3]. In this context, equation of type (1.3) is usually called the nonlinear Schrödinger-Newtonequation. The first investigations for the existence and symmetry of the solutions to (1.3) go back to the works of Lieb [2] and Lions [4]. In [2], by using symmetric decreasing rearrangement inequalities, Lieb proved that the ground state solution of Equation (1.3) is radial and unique up to translations. Lions [4] showed the existence of a sequence of radially symmetric solutions. Ma and Zhao [5] considered the generalized Choquard equation

$-\Delta u+u=\left[{\left|x\right|}^{-\mu }\ast {\left|u\right|}^q\right]{\left|u\right|}^{q-2}u\left(q\ge 2\right)\mathrm{,}$ (1.4)

and proved that every positive solution of it is radially symmetric and monotone decreasing about some fixed point, under the assumption that a certain set of real numbers, defined in terms of N, and q, is nonempty. Under the same assumption, Cingolani, Clapp, and Secchi [6] gave some existence and multiplicity results in the electromagnetic case and established the regularity and some decay asymptotically at infinity of the ground states. In [7], Moroz and Van Schaftingen eliminated this restriction and showed the regularity, positivity and radial symmetry of the ground states for the optimal range of parameters and derived decay asymptotically at infinity for them as well. Moreover, they [8] also obtained a similar conclusion under the assumption of Berestycki-Lions type nonlinearity. We point out that the existence, multiplicity, and concentration of such like equation have been established by many authors. We refer the readers to [9] [10] for the existence of sign-changing solutions, [11] [12] for the existence and concentration behavior of the semiclassical solutions and [13] for the critical nonlocal part with respect to the Hardy-Littlewood-Sobolev inequality. For more details associated with the Choquard equation, please refer to [14] [15] [16] and the references in. Li, Teng, Zhang, Nie [17] investigate the existence of solutions for the following generalized quasilinear Schrödinger equation with nonlocal term

$\begin{array}{l}-div\left(g^2\left(u\right)\nabla u\right)+g\left(u\right)g^{\prime }\left(u\right)+{\left|\nabla u\right|}^2+V\left(x\right)u\\ =\lambda \left[{\left|x\right|}^{-\mu }\ast {\left|u\right|}^p\right]{\left|u\right|}^{p-2}u\mathrm{,}x\in {ℝ}^N\end{array}$ (1.5)

and prove the existence of solution.

In this paper, our main ideas come from [18] and the assumption of g from [19]. Our purpose is to search for the existence of nontrivial solutions of (1.1) by implicit function theorem. For convenience, we introduce several notations: C denotes a positive (possibly different) constant, $L^p\left({ℝ}^N\right)$ denotes the usual Lebesgue space with norms ${\Vert u\Vert }_{L^p\left({ℝ}^N\right)}={\left({\displaystyle {\int }_{{ℝ}^N}}{\left|u\right|}^p\text{d}x\right)}^{\frac{1}{p}}$ , $1\le p<\infty$ , $C_0^{\infty }\left({ℝ}^N\right)$ be the collection of smooth functions with compact support. Next, we introduce the energy functional of Equation (1.1)

$\begin{array}{c}I\left(u\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left[g^2\left(u\right){\left|\nabla u\right|}^2+V\left(x\right)u^2\right]\text{d}x\\ -\frac{\lambda }{2p}{\displaystyle {\int }_{{ℝ}^N}}\left[{\left|x\right|}^{-\mu }\ast {\left|u\right|}^p\right]{\left|u\right|}^p\text{d}x-\frac{\beta }{q+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|u\right|}^{q+1}\text{d}x\end{array}$ (1.6)

however, J is not well defined in $H^1\left({ℝ}^N\right)$ because of the term ${\displaystyle {\int }_{{ℝ}^N}}{g}^2\left(u\right){\left|\nabla u\right|}^2\text{d}x$ . To overcome this difficulty, we make a change of variable constructed by Shen and Wang in [20]: $v=G\left(u\right)={\displaystyle {\int }_0^u}g\left(s\right)\text{d}s$ , then,

$\begin{array}{c}J\left(v\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left[{\left|\nabla v\right|}^2+V\left(x\right){\left|G^{-1}\left(v\right)\right|}^2\right]\text{d}x\\ -\frac{\lambda }{2p}{\displaystyle {\int }_{{ℝ}^N}}\left(\left[{\left|x\right|}^{-\mu }\ast {\left|G^{-1}\left(v\right)\right|}^p\right]{\left|G^{-1}\left(v\right)\right|}^p\right)\text{d}x+{\displaystyle {\int }_{{ℝ}^N}}\left[\frac{\beta }{q+1}{\left|G^{-1}\left(v\right)\right|}^{q+1}\right]\text{d}x\\ =\frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left[{\left|\nabla v\right|}^2+V\left(x\right){\left|G^{-1}\left(v\right)\right|}^2\right]\text{d}x\\ -\frac{\lambda }{2p}{\displaystyle {\int }_{{ℝ}^{2N}}}\left(\frac{{\left|G^{-1}\left(v\left(y\right)\right)\right|}^p{\left|G^{-1}\left(v\left(x\right)\right)\right|}^p}{{\left|x-y\right|}^{\mu }}\right)\text{d}x\text{d}y-{\displaystyle {\int }_{{ℝ}^N}}\left[\frac{\beta }{q+1}{\left|G^{-1}\left(v\right)\right|}^{q+1}\right]\text{d}x\end{array}$ (1.7)

We say that u is a weak solution of (1.1), if

$\begin{array}{c}\langle I^{\prime }\left(u\right),\phi \rangle ={\displaystyle {\int }_{R^N}}\left(g^2\left(u\right)\nabla u\nabla \phi +g\left(u\right)g^{\prime }\left(u\right){\left|\nabla u\right|}^2\phi +V\left(x\right)u\phi \right)\text{d}x\\ -\lambda {\displaystyle {\int }_{{ℝ}^N}}\left[{\left|x\right|}^{-\mu }\ast {\left|u\right|}^p\right]{\left|u\right|}^{p-2}u\phi -\beta {\displaystyle {\int }_{{ℝ}^N}}{\left|u\right|}^{q-1}u\phi \text{d}x\end{array}$ (1.8)

for all $\phi \in C_0^{\infty }\left({ℝ}^N\right)$ . Let $\phi =\left(\frac{1}{g\left(u\right)}\psi \right)$ , by [20], we know that the above formula is equivalent to

$\begin{array}{c}\langle J^{\prime }\left(v\right),\psi \rangle ={\displaystyle {\int }_{{ℝ}^N}}\nabla v\nabla \psi +V\left(x\right)\frac{G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\psi \text{d}x\\ -\lambda {\displaystyle {\int }_{{ℝ}^N}}\frac{\left[{\left|x\right|}^{-2}\ast {\left|G^{-1}\left(v\right)\right|}^p\right]{\left|G^{-1}\left(v\right)\right|}^{p-2}{G}^{-1}\left(v\right)\psi }{g\left(G^{-1}\left(v\right)\right)}\text{d}x\\ -\beta {\displaystyle {\int }_{{ℝ}^N}}\frac{{\left|G^{-1}\left(v\right)\right|}^{q-1}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\psi \text{d}x\end{array}$ (1.9)

for all $\psi \in C_0^{\infty }\left({ℝ}^N\right)$ . Therefore, in order to find the solution of (1.1), it suffices to study the solution of following equation:

$\begin{array}{l}-\Delta v+V\left(x\right)\frac{G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}-\lambda \frac{\left[{\left|x\right|}^{-\mu }\ast {\left|G^{-1}\left(v\right)\right|}^p\right]{\left|G^{-1}\left(v\right)\right|}^{p-2}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\\ -\beta \frac{{\left|G^{-1}\left(v\right)\right|}^{q-1}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}=0\end{array}$ (1.10)

J is defined on the space

$H_V^1\left({ℝ}^N\right)=\left\{v\in H^1\left({ℝ}^N\right):{\displaystyle {\int }_{{ℝ}^N}}{\left|\nabla v\right|}^2+V\left(x\right)v^2\text{d}x<+\infty \right\}$

we can define the norm on $H_V^1\left({ℝ}^N\right)$ by

${\Vert v\Vert }_{H_V^1\left({ℝ}^N\right)}^2={\displaystyle {\int }_{{ℝ}^N}}{\left|\nabla v\right|}^2+V\left(x\right)v^2\text{d}x$

then, $H_V^1\left({ℝ}^N\right)$ is a Banach space. In the following, we always assume $V\in C\left({ℝ}^N\mathrm{,}{ℝ}^{+}\right)$ and ${\inf }_{{ℝ}^N}V\left(x\right)\ge 1$ . Let us consider the following assumptions of potential function $V\left(x\right)$ :

(V₁) ${\lim }_{\left|x\right|\to \infty }V\left(x\right)=+\infty$ ;

(V₂) $V\left(x\right)$ is radially symmetric.

Next, we will introduce the properties of some functions.

Lemma 1.1. [19] The function $g\left(t\right)\mathrm{,}{G}^{-1}\left(t\right)\mathrm{,}G\left(t\right)$ enjoys the following properties.

(g₁) the function $G\left(t\right)$ and $G^{-1}\left(t\right)$ are strictly increasing and odd;

(g₂) $\left|t\right|\le \left|G^{-1}\left(t\right)\right|\le \left|t\right|/a$ for all $t\in ℝ$ ;

(g₃) $G^{-1}\left(t\right)/t$ is nondecreasing for all $t\in ℝ$ and ${\lim }_{t\to 0}{G}^{-1}\left(t\right)/t=1$ , ${\lim }_{t\to \infty }{G}^{-1}\left(t\right)/t=1/a$ ;

(g₄) $t^2\le \left(t/g\left(t\right)\right)G\left(t\right)\le t^2/a$ for all $t\in ℝ$ .

Next, we set forth some preliminary results.

2. Preliminary Results

To begin with, we prove some functions are continuous, more detailed see [21].

Lemma 2.1 If ${\Vert v_n-v\Vert }_{H_V^1\left({ℝ}^N\right)}\to 0$ , then ${\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)G^{-1}{\left(v_n\right)}^2-V\left(x\right)G^{-1}{\left(v\right)}^2\text{d}x\to 0$ .

Proof: By sobolev imbedding inequality, Lemma1.1 and definition of g, we have

$\begin{array}{l}{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)G^{-1}{\left(v_n\right)}^2-V\left(x\right)G^{-1}{\left(v\right)}^2\text{d}x\\ =2{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)\frac{G^{-1}\left(v+\theta \left(v_n-v\right)\right)}{g\left(G^{-1}\left(v+\theta \left(v_n-v\right)\right)\right)}\left(v_n-v\right)\text{d}x\\ \le \frac{2}{a^2}{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)\left|v+\theta \left(v_n-v\right)\right|\left|v_n-v\right|\text{d}x\\ \le \frac{2}{a^2}{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right){\left|v_n-v\right|}^2\text{d}x{\displaystyle {\int }_{R^N}}V\left(x\right){\left|v+\theta \left(v_n-v\right)\right|}^2\text{d}x\\ \le \frac{2}{a^2}{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right){\left|v_n-v\right|}^2\text{d}x\left({\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)v^2\text{d}x+{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right){\left|v_n-v\right|}^2\text{d}x\right)\\ \le \frac{2}{a^2}{\Vert v_n-v\Vert }_{H_V^1\left({ℝ}^N\right)}\left(C+{\Vert v_n-v\Vert }_{H_V^1\left({ℝ}^N\right)}\right)\to 0,n\to +\infty \end{array}$

where $\theta \in \left(\mathrm{0,1}\right)$ . $■$

Lemma 2.2. The map: $v\to G^{-1}\left(v\right)$ from $H_V^1\left({ℝ}^N\right)$ into $L^r\left({ℝ}^N\right)$ is continuous for $2\le r<2^{\ast }$ .

Proof: By the definition of g, we have

$\begin{array}{c}{\displaystyle {\int }_{{ℝ}^N}}{\left|G^{-1}\left(v_n\right)-G^{-1}\left(v\right)\right|}^r\text{d}x\le {\displaystyle {\int }_{{ℝ}^N}}{\left|G^{-1}\left(v_n\right)-G^{-1}\left(v\right)\right|}^r\text{d}x\\ \le \frac{1}{a^r}{\displaystyle {\int }_{{ℝ}^N}}\left({\left|v_n\right|}^r+{\left|v\right|}^r\right)\text{d}x\mathrm{.}\end{array}$

Assume $v_n\to v$ in $H_V^1\left({ℝ}^N\right)$ , moreover, the imbedding from $H_V^1\left({ℝ}^N\right)$ into $L^r\left({ℝ}^N\right)$ is compact where $r\in \left[{\mathrm{2,2}}^{\ast }\right)$ , from Lemma 3.4 [22], we get the result. $■$

Next, we introduce some minimization with corresponding energy functional and define

$m_b=\underset{u\in M_b}{\inf }E(\; u\; )$

where

$M_b=\left\{u\in H^1\left({ℝ}^N\right):{\Vert u\Vert }_{L^{p+1}}=b\right\},a>0$

and

$E\left(u\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left[g^2\left(u\right){\left|\nabla u\right|}^2+V\left(x\right)u^2\right]dx-\frac{\beta }{q+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|u\right|}^{q+1}\text{d}x.$

We also define

${\omega }_b=\underset{v\in W_b}{\inf }F(\; v\; )$

where

$W_b=\left\{v\in H_V^1\left({ℝ}^N\right):{\Vert G^{-1}\left(v\right)\Vert }_{L^{p+1}}=b\right\},b>0$

and

$F\left(v\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left({\left|\nabla v\right|}^2+V\left(x\right)G^{-1}{\left(v\right)}^2\right)\text{d}x-\frac{\beta }{q+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|G^{-1}\left(v\right)\right|}^{q+1}\text{d}x$

Therefore, we have following fact.

Lemma 2.3. $m_b={\omega }_b$ for every $a>0$ .

Proof: For any $v\in W_b$ , let $u=G^{-1}\left(v\right)$ , from the definition of g, we get

$\begin{array}{l}{\displaystyle {\int }_{{ℝ}^N}}{\left|\nabla u\right|}^2\text{d}x={\displaystyle {\int }_{{ℝ}^N}}\frac{{\left|\nabla v\right|}^2}{g^2\left(G^{-1}\left(v\right)\right)}\text{d}x\le \frac{1}{a^2}{\displaystyle {\int }_{{ℝ}^N}}{\left|\nabla v\right|}^2\text{d}x<+\infty ,\\ {\displaystyle {\int }_{R^N}}{u}^2\text{d}x\le {\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)G^{-1}{\left(v\right)}^2\text{d}x<+\infty ,\end{array}$

so $u\in M_b$ . It follows that $F\left(v\right)=E\left(G^{-1}\left(v\right)\right)=E\left(u\right)\ge m_b$ , hence ${\omega }_b\ge m_b$ , moreover, for any $u\in M_b$ , let $v=G\left(u\right)$ , then $u=G^{-1}\left(v\right)$ . We assume $E\left(u\right)<+\infty$ , since $u\in H_V^1\left({ℝ}^N\right)$ , $2<q+1<2^{\ast }$ , then $u\in L^{q+1}\left({ℝ}^N\right)$ . We have

$\frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left[g^2\left(u\right){\left|\nabla u\right|}^2+V\left(x\right)u^2\right]\text{d}x=E\left(u\right)+\frac{\beta }{q+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|u\right|}^{q+1}\text{d}x<+\infty$ .

Then ${\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)G^{-1}{\left(v\right)}^2={\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)u^2\text{d}x<+\infty$ . It shows that $v\in W_a$ , which implies that $E\left(u\right)=E\left(G^{-1}\left(v\right)\right)=F\left(v\right)\ge {\omega }_b$ , hence $m_b\ge {\omega }_b$ , this completes the proof. $■$

Lemma 2.4. 1) $F\left(v\right)$ iswell defined and continuous for $2\le r\le 2^{\ast }$ .

2) $F\left(v\right)$ is Gateaux-differentiable. For $v\in H_V^1\left({ℝ}^N\right)$ ,the G-derivative $F^{\prime }\left(v\right)$ is a continuous function, and $F^{\prime }\left(v\right)$ is continuous in v in the strongly-weak topology, that is,if $v_n\to v$ strongly in $H_V^1\left({ℝ}^N\right)$ ,then $F^{\prime }\left(v_n\right)\rightharpoonup F^{\prime }\left(v\right)$ weakly.

Proof: (1) For any $v\in H_V^1\left({ℝ}^N\right)$ , we have ${\displaystyle {\int }_{{ℝ}^N}}{\left|G^{-1}\left(v\right)\right|}^{q+1}\text{d}x\le C{\displaystyle {\int }_{{ℝ}^N}}{\left|v\right|}^{q+1}\text{d}x<+\infty$ , then

$\begin{array}{c}F\left(v\right)=\frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}{\left|\nabla v\right|}^2+V\left(x\right)G^{-1}{\left(v\right)}^2\text{d}x-\frac{\beta }{q+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|G^{-1}\left(v\right)\right|}^{q+1}\text{d}x\\ \le \frac{1}{2a^2}{\displaystyle {\int }_{{ℝ}^N}}{\left|\nabla v\right|}^2+V\left(x\right)v^2\text{d}x+\left|\frac{\beta }{q+1}\right|{\displaystyle {\int }_{{ℝ}^N}}{\left|G^{-1}\left(v\right)\right|}^{q+1}\text{d}x\\ \le \frac{1}{2a^2}{\displaystyle {\int }_{{ℝ}^N}}{\left|\nabla v\right|}^2+V\left(x\right)v^2\text{d}x+\frac{\beta }{\left(q+1\right)a^{q+1}}{\displaystyle {\int }_{{ℝ}^N}}{\left|v\right|}^{q+1}\text{d}x<+\infty .\end{array}$

with the proof of continuity, note that J consists of three terms. By Lemma 1.1, we need to check the superlinear term only.

$\begin{array}{l}\left|\frac{1}{q+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|G^{-1}\left(v_n\right)\right|}^{q+1}\text{d}x-\frac{1}{q+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|G^{-1}\left(v\right)\right|}^{q+1}\text{d}x\right|\\ =\left|{\displaystyle {\int }_0^1}\text{d}t{\displaystyle {\int }_{{ℝ}^N}}\frac{{\left|G^{-1}\left(v+t\left(v_n-v\right)\right)\right|}^{q-1}{G}^{-1}\left(v+t\left(v_n-v\right)\right)}{g\left(G^{-1}\left(v+t\left(v_n-v\right)\right)\right)}\left(v_n-v\right)\text{d}x\right|\\ \le C{\displaystyle {\int }_0^1}{\displaystyle {\int }_{{ℝ}^N}}{\left|v+t\left(v_n-v\right)\right|}^q\left|v_n-v\right|\text{d}t\text{d}x\\ \le {\Vert {\left|v+t\left(v_n-v\right)\right|}^q\Vert }_{L^2\left({ℝ}^N\right)}{\Vert v_n-v\Vert }_{L^2\left({ℝ}^N\right)}\\ \le C{\Vert v_n-v\Vert }_{H_V^1}\end{array}$

where $1\le q\le \frac{N}{N-2}$ .

For (2) we consider the second and the third terms of the functional J, we see for $\varphi \in H_V^1\left({ℝ}^N\right)$ , using Hölder inequality, we get

$\begin{array}{l}\left|\frac{1}{2t}{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)\left(G^{-1}{\left(v+t\varphi \right)}^2-G^{-1}{\left(v\right)}^2\right)\text{d}x-{\displaystyle {\int }_{{ℝ}^N}}\frac{V\left(x\right)G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\varphi \text{d}x\right|\\ =\left|{\displaystyle {\int }_0^1}\text{d}s{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)\left(\frac{G^{-1}\left(v+ts\varphi \right)}{g\left(G^{-1}\left(v+ts\varphi \right)\right)}-\frac{G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\right)\varphi \text{d}x\right|\\ \le {\displaystyle {\int }_0^1}\text{d}s{\left({\displaystyle {\int }_{{ℝ}^N}}V\left(x\right){\left|\frac{G^{-1}\left(v+ts\varphi \right)}{g\left(G^{-1}\left(v+ts\varphi \right)\right)}-\frac{G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\right|}^2\text{d}x\right)}^{\frac{1}{2}}{\displaystyle {\int }_0^1}\text{d}s{\left({\displaystyle {\int }_{{ℝ}^N}}V\left(x\right){\varphi }^2\text{d}x\right)}^{\frac{1}{2}}\end{array}$ (2.1)

Using the definition of g and Lemma 1.1, we know

$\begin{array}{c}{\left|\frac{G^{-1}\left(v+ts\varphi \right)}{g\left(G^{-1}\left(v+ts\varphi \right)\right)}-\frac{G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\right|}^2\le {\left|G^{-1}\left(v+ts\varphi \right)+G^{-1}\left(v\right)\right|}^2\\ \le C\left({\left|G^{-1}\left(v+ts\varphi \right)\right|}^2+{\left|G^{-1}\left(v\right)\right|}^2\right)\\ \le C\left({\left|v+ts\varphi \right|}^2+{\left|v\right|}^2\right)\\ \le C\left({\left|v\right|}^2+{\left|\varphi \right|}^2\right)\mathrm{.}\end{array}$

By the dominated convergence theorem

$\left|{\displaystyle {\int }_0^1}\text{d}s{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)\left(\frac{G^{-1}\left(v+ts\varphi \right)}{g\left(G^{-1}\left(v+ts\varphi \right)\right)}-\frac{G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\right)\varphi \text{d}x\right|\to 0\left(t\to 0\right)$

For the third term, we have

$\begin{array}{l}\left|\frac{1}{t}{\displaystyle {\int }_{{ℝ}^N}}\frac{1}{q+1}\left({\left|G^{-1}\left(v+t\varphi \right)\right|}^{q+1}-{\left|G^{-1}\left(v\right)\right|}^{q+1}\right)\text{d}x-{\displaystyle {\int }_{{ℝ}^N}}\frac{G^{-1}{\left(v\right)}^{q-1}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\varphi \text{d}x\right|\\ =\left|\frac{1}{t}{\displaystyle {\int }_0^1}\text{d}s{\displaystyle {\int }_{{ℝ}^N}}\left(\frac{{\left|G^{-1}\left(v+ts\varphi \right)\right|}^{q-1}{G}^{-1}\left(v+ts\varphi \right)}{g\left(G^{-1}\left(v+ts\varphi \right)\right)}-\frac{{\left|G^{-1}\left(v\right)\right|}^{q-1}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\right)\varphi \text{d}x\right|\end{array}$

Similarly to above, by the dominated convergence theorem

$\underset{t\to 0}{\lim }\left|\frac{1}{t}{\displaystyle {\int }_0^1}\text{d}s{\displaystyle {\int }_{{ℝ}^N}}\left(\frac{{\left|G^{-1}\right|}^{q-1}\left(v+ts\varphi \right)}{g\left(G^{-1}\left(v+ts\varphi \right)\right)}-\frac{{\left|G^{-1}\right|}^{q-1}\left(v\right)G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\right)\varphi \text{d}x\right|=0$

The Gateaux derivative $J^{\prime }\left(v\right)$ has the form

$\begin{array}{c}\langle J^{\prime }\left(v\right),\varphi \rangle ={\displaystyle {\int }_{{ℝ}^N}}\nabla v\nabla \varphi \text{d}x+{\displaystyle {\int }_{{ℝ}^N}}\frac{V\left(x\right)G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\varphi \text{d}x-{\displaystyle {\int }_{{ℝ}^N}}\frac{{\left|G^{-1}\left(v\right)\right|}^{q-1}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\varphi \text{d}x\\ \le C{\Vert v\Vert }_{H_V^1\left({ℝ}^N\right)}{\Vert \varphi \Vert }_{H_V^1\left({ℝ}^N\right)}+C{\Vert G^{-1}\left(v\right)\Vert }_{L^{2q}\left({ℝ}^N\right)}{\Vert \varphi \Vert }_{L^2(\; ℝ\; N\; )}\end{array}$

from Sobolev imbedding theorem, we get $J^{\prime }\left(v\right)$ is a continuous linear functional on $H_V^1\left({ℝ}^N\right)$ .

Finally, the continuity with strong-weak topology is easy to check, as $v_n\to v$ in $H_V^1\left({ℝ}^N\right)$ , for any $\varphi \in H_V^1\left({ℝ}^N\right)$ ,

${\displaystyle {\int }_{{ℝ}^N}}V\left(x\right)\frac{G^{-1}\left(v_n\right)}{g\left(G^{-1}\left(v_n\right)\right)}-V\left(x\right)\frac{G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\text{d}x\to 0.$ $■$

Remark 2.5. Lemma 2.4 does not show that $F\left(v\right)$ is $C^{\infty }$ , so we cannot use the Lagrange multiplier theorem. But we can get our conclusion we want exactly by a similar argument for the Lagrange multiplier theorem. Next, we state our main conclusion. The idea of our proof is based on the work in [18] [22] [23].

3. Main Conclusion

Theorem 3.1. Let $N\ge 3$ , $1\le q<\frac{N}{N-2}$ , $\frac{2N-\mu }{N}\le p<\frac{2N-\mu }{N-2}$ , $\beta \in ℝ$ and

$q<p$ . Assume (V₁) or (V₂) holds. Then for every $b>0$ , there exists $\lambda \left(b\right)\in ℝ$ such that Equation (1.1) with $\lambda =\lambda \left(b\right)$ has a positive weak solution $u\in M_b$ .

Remark 3.2. From the assumption of V, we know $H_V^1\left({ℝ}^N\right)$ embedding into $L^p\left({ℝ}^N\right)$ is compact. In the process of the proof of theorem 3.1, it is important for us to construct auxiliary function, then by implicit function theorem to prove it and lemma 3.4 [22] play a great role in this paper. Moreover, when $\frac{N-2}{N}<q\le 2^{*}$ is a open question for Equation (1.1), someone could do it if they are interested.

Proof of Theorem 3.1: Step 1: By the assumptions of (V₁) or (V₂), ${\omega }_b$ is achieved at some $0\le v_b\le W_b$ with $v_b\ne 0$ . Let $\left\{v_n\right\}\in W_b$ be a minimizing sequence for ${\omega }_b$ . Set $u_n=G^{-1}\left(v_n\right)$ . Then $\left\{u_n\right\}\in M_b$ is a minimizing sequence for $m_b$ . We can assume $u_n\ge 0$ . It shows that $E\left(u_n\right)\to m_b$ , so there exists $C>0$ such that

$\begin{array}{c}C\ge E\left(u_n\right)\\ =\frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left[g^2\left(u_n\right){\left|\nabla u_n\right|}^2+V\left(x\right)u_n^2\right]\text{d}x-\frac{\beta }{q+1}{\displaystyle {\int }_{R^N}}{\left|u_n\right|}^{q+1}\text{d}x\\ \ge \frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left[{\left|\nabla u\right|}^2+{\left|u_n\right|}^2\right]\text{d}x-\frac{\beta }{q+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|u_n\right|}^{q+1}\text{d}x.\end{array}$

By Hölder inequality,

$\begin{array}{c}{\displaystyle {\int }_{{ℝ}^N}}{\left|u_n\right|}^{q+1}\text{d}x\le {\left({\displaystyle {\int }_{{ℝ}^N}}{\left|u_n\right|}^2\text{d}x\right)}^{\frac{\lambda \left(q+1\right)}{2}}{\left({\displaystyle {\int }_{{ℝ}^N}}{\left|u_n\right|}^{p+1}\text{d}x\right)}^{\frac{\left(1-\lambda \right)\left(q+1\right)}{p+1}}\\ \le \frac{\lambda \left(q+1\right)}{2}{\displaystyle {\int }_{{ℝ}^N}}{\left|u_n\right|}^2\text{d}x+\frac{\left(1-\lambda \right)\left(q+1\right)}{p+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|u_n\right|}^{p+1}\text{d}x\\ =\frac{\lambda \left(q+1\right)}{2}{\displaystyle {\int }_{{ℝ}^N}}{\left|u_n\right|}^2\text{d}x+\frac{\left(1-\lambda \right)\left(q+1\right)}{p+1}{a}^{p+1}\end{array}$

where $\lambda =\frac{2\left(p-q\right)}{\left(q+1\right)\left(p-1\right)}$ . Then

$\begin{array}{c}C\ge E\left(u_n\right)\\ \ge \frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left[{\left|\nabla u_n\right|}^2+V\left(x\right)u_n^2\right]\text{d}x\\ -\frac{\beta }{q+1}\left(\frac{\lambda \left(q+1\right)}{2}{\displaystyle {\int }_{{ℝ}^N}}V\left(x\right){\left|u_n\right|}^2\text{d}x+\frac{\left(1-\lambda \right)\left(q+1\right)}{p+1}{a}^{p+1}\right)\\ \ge \left(\frac{1}{2}-\frac{\beta \left(p-q\right)}{\left(q+1\right)\left(p-1\right)}\right)\left({\displaystyle {\int }_{{ℝ}^N}}\left[{\left|\nabla u_n\right|}^2+V\left(x\right){\left|u_n\right|}^2\right]\text{d}x\right)-\frac{\beta \left(q-1\right)}{\left(q+1\right)\left(p-1\right)}{a}^{p+1}\end{array}$

Because of $\beta <\frac{N-2}{2N}$ , $\frac{1}{2}-\frac{\beta \left(p-q\right)}{\left(q+1\right)\left(p-1\right)}>0$ . It implies that $u_n\left(x\right)$ is bounded in $H_V^1\left({ℝ}^N\right)$ . By the compact embedding result from $H_V^1\left({ℝ}^N\right)$ into $L^r\left({ℝ}^N\right)$ for $2\le r<2^{\ast }$ . We may assume that $u_n\rightharpoonup u_b$ in $H_V^1\left({ℝ}^N\right)$ , $u_n\to u_b$ in $L^r\left({ℝ}^N\right)$ for $2\le r<2^{\ast }$ and $u_n\left(x\right)\to u_b\left(x\right)$ a.e $x\in {ℝ}^N$ . Hence $u_b\in M_b$ , since $u_n\ge \mathrm{0,}{u}_b\ge 0$ and $u_b\ne 0$ .

Using the same argument as the process of the proof of Lemma 2.1 in [20] and noting that $u_n\to u_b$ in $L^{q+1}\left({ℝ}^N\right)$ . We have

$\begin{array}{c}{m}_b=\underset{n\to \infty }{\lim }E\left(u_n\right)\\ \ge \underset{n\to \infty }{\lim \inf }\left\{\frac{1}{2}{\displaystyle {\int }_{{ℝ}^N}}\left[g^2\left(u_n\right){\left|\nabla u_n\right|}^2+V\left(x\right)u_n^2\right]\text{d}x-\frac{\beta }{q+1}{\displaystyle {\int }_{{ℝ}^N}}{\left|u_n\right|}^{q+1}\text{d}x\right\}\\ \ge E(\; u\; b\; )\end{array}$

Hence $m_b$ is achieved at $u_b$ and

$v_b=G\left(u_b\right)\in W_b,F\left(v_b\right)=E\left(G^{-1}\left(v_b\right)\right)=E\left(u_b\right)={\omega }_b$

and the property of g implies $v_b\ge 0$ and $v_b\ne 0$ .

Step 2: Set $h_p=\frac{1}{2p}{\displaystyle {\int }_{{ℝ}^{2N}}}\frac{{\left|G^{-1}\left(v\left(x\right)\right)\right|}^p{\left|G^{-1}\left(v\left(y\right)\right)\right|}^p}{{\left|x-y\right|}^{\mu }}\text{d}x\text{d}y$ for $\frac{2N-\mu }{N}<p<\frac{2N-\mu }{N-2}$ . Then $h_p\left(v\right)\in C^1\left(H_V^1\left({ℝ}^N\right)\mathrm{,}ℝ\right)$ . Actually, for any $\phi \in H_V^1\left({ℝ}^N\right)$ , by Sobolev inequality and Hölder inequality, we get

$\begin{array}{c}\left|\langle {h^{\prime }}_p\left(v\right),\phi \rangle \right|=\left|{\displaystyle {\int }_{{ℝ}^N}}\frac{\left[{\left|x\right|}^{-\mu }\ast {\left|G^{-1}\left(v\right)\right|}^p\right]{\left|G^{-1}\left(v\left(x\right)\right)\right|}^{r-2}{G}^{-1}\left(v\left(x\right)\right)}{g\left(G^{-1}\left(v\left(x\right)\right)\right)}\phi \right|\text{d}x\\ \le {\displaystyle {\int }_{{ℝ}^N}}{\displaystyle {\int }_{{ℝ}^N}}\frac{{\left|G^{-1}\left(v\left(y\right)\right)\right|}^p{\left|G^{-1}\left(v\left(x\right)\right)\right|}^{p-1}}{{\left|x-y\right|}^{\mu }}\left|\phi \right|\text{d}y\text{d}x\\ \le C{\Vert G^{-1}\left(v\left(y\right)\right)\Vert }_{L^{pr}\left({ℝ}^N\right)}^p{\Vert G^{-1}\left(v\left(x\right)\right)\Vert }_{L^{pr}\left({ℝ}^N\right)}^{p-1}{\Vert \phi \Vert }_{L^{pr}\left({ℝ}^N\right)}\\ \le C{\Vert \phi \Vert }_{L^{pr}\left({ℝ}^N\right)}\\ \le C{\Vert \phi \Vert }_{H_V^1(\; ℝ\; N\; )}\end{array}$

where $r=\frac{2N}{2N-\mu }$ and $2\le pr<2^{\ast }$ , so ${h^{\prime }}_p\left(v\right)\in {\left(H_V^1\left({ℝ}^N\right)\right)}^{\ast }$ .

Let $v_n\to v$ in $H_V^1\left({ℝ}^N\right)$ . Up to a subsequence, we can assume $v_n\to v$ a.e. in $ℝ$ and $v_n\to v$ in $L^r\left({ℝ}^N\right)$ for $2\le r<2^{\ast }$ . Hence

$\begin{array}{l}\left|\langle {h^{\prime }}_p\left(v_n\right)-{h^{\prime }}_p\left(v\right),\phi \rangle \right|\\ =|{\displaystyle {\int }_{{ℝ}^N}}(\frac{\left[{\left|x\right|}^{-\mu }\ast {\left|G^{-1}\left(v_n\right)\right|}^p\right]{\left|G^{-1}\left(v_n\right)\right|}^{p-2}{G}^{-1}\left(v_n\right)}{g\left(G^{-1}\left(v_n\right)\right)}\\ -\frac{\left[{\left|x\right|}^{-\mu }\ast {\left|G^{-1}\left(v\right)\right|}^p\right]{\left|G^{-1}\left(v\right)\right|}^{p-2}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)})\phi \text{d}x|\end{array}$

$\begin{array}{l}\le \left|{\displaystyle {\int }_{{ℝ}^N}}\frac{\left[{\left|x\right|}^{-\mu }\ast \left({\left|G^{-1}\left(v_n\right)\right|}^p-{\left|G^{-1}\left(v\right)\right|}^p\right)\right]{\left|G^{-1}\left(v_n\right)\right|}^{p-2}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v_n\right)\right)}\phi \right|\\ +\left|{\displaystyle {\int }_{{ℝ}^N}}\left[{\left|x\right|}^{-\mu }\ast {\left|G^{-1}\left(v\right)\right|}^p\right]\left(\frac{{\left|G^{-1}\left(v_n\right)\right|}^{p-1}{G}^{-1}\left(v_n\right)}{g\left(G^{-1}\left(v_n\right)\right)}-\frac{{\left|G^{-1}\left(v\right)\right|}^{p-1}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\right)\phi \text{d}x\right|\\ \le C{\Vert {\left|G^{-1}\left(v_n\right)\right|}^p-{\left|G^{-1}\left(v\right)\right|}^p\Vert }_{L^r\left({ℝ}^N\right)}{\Vert G^{-1}\left(v_n\right)\Vert }_{L^{pr}\left({ℝ}^N\right)}^{p-1}{\Vert \phi \Vert }_{L^{pr}\left({ℝ}^N\right)}\\ +C{\Vert G^{-1}\left(v\right)\Vert }_{L^{pr}\left({ℝ}^N\right)}^p{\Vert \left(\frac{{\left|G^{-1}\left(v_n\right)\right|}^{p-1}{G}^{-1}\left(v_n\right)}{g\left(G^{-1}\left(v_n\right)\right)}-\frac{{\left|G^{-1}\left(v\right)\right|}^{p-1}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\right)\phi \Vert }_{L^r\left({ℝ}^N\right)}\end{array}$ (3.1)

Since $v_n\to v$ in $H_V^1\left({ℝ}^N\right)$ , we get ${\Vert v_n\Vert }_{H_V^1\left({ℝ}^N\right)}\le C$ , $2\le pr<2^{*}$ and by definition g, then

$\begin{array}{l}{\Vert {\left|G^{-1}\left(v_n\right)\right|}^p-{\left|G^{-1}\left(v\right)\right|}^p\Vert }_{L^r\left({ℝ}^N\right)}{\Vert G^{-1}\left(v_n\right)\Vert }_{L^{pr}\left({ℝ}^N\right)}^{p-1}{\Vert \phi \Vert }_{L^{pr}\left({ℝ}^N\right)}\\ \le C\left({\Vert v_n\Vert }_{L^{pr}\left({ℝ}^N\right)}^p+{\Vert v\Vert }_{L^{pr}\left({ℝ}^N\right)}^p\right),\end{array}$

then by Lemma 3.4 [22] and assumption, we get

${\Vert {\left|G^{-1}\left(v_n\right)\right|}^p-{\left|G^{-1}\left(v\right)\right|}^p\Vert }_{L^r\left({ℝ}^N\right)}{\Vert G^{-1}\left(v_n\right)\Vert }_{L^{pr}\left({ℝ}^N\right)}^{p-1}{\Vert \phi \Vert }_{L^{pr}\left({ℝ}^N\right)}\to \mathrm{0,}n\to +\infty \mathrm{.}$

Similarly,

${\Vert G^{-1}\left(v\right)\Vert }_{L^{pr}\left({ℝ}^N\right)}^p{\Vert \left(\frac{{\left|G^{-1}\left(v_n\right)\right|}^{p-1}{G}^{-1}\left(v_n\right)}{g\left(G^{-1}\left(v_n\right)\right)}-\frac{{\left|G^{-1}\left(v\right)\right|}^{p-1}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\right)\phi \Vert }_{L^r\left({ℝ}^N\right)}\to \mathrm{0,}n\to +\infty \mathrm{.}$

Hence, follows that $h_p\left(v\right)\in C^1\left(H_V^1\left({ℝ}^N\right)\mathrm{,}ℝ\right)$ for $\frac{2N-\mu }{N}<p<\frac{2N-\mu }{N-2}$ .

Step 3: For any $b\ge 0$ , there exists $\lambda \left(b\right)\in ℝ$ such that $0<u_b=G^{-1}\left(v_b\right)\in M_b$ is a weak solution of Equation (1.1) with $\lambda =\lambda \left(b\right)$ . In fact, by Lemma 2.4,

$\langle F\left(v\right),\phi \rangle ={\displaystyle {\int }_{{ℝ}^N}}\nabla v\nabla \phi +{\displaystyle {\int }_{{ℝ}^N}}\frac{V\left(x\right)G^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\phi \text{d}x-\beta {\displaystyle {\int }_{{ℝ}^N}}\frac{{\left|G^{-1}\left(v\right)\right|}^{q-1}{G}^{-1}\left(v\right)}{g\left(G^{-1}\left(v\right)\right)}\phi \text{d}x$

and $F\left(v\right)\in {\left(H_V^1\left({ℝ}^N\right)\right)}^{\ast }$ for all $v\in H_V^1\left({ℝ}^N\right)$ . Since ${h^{\prime }}_p\left(v\right)\in C\left(H_V^1\left({ℝ}^N\right)\mathrm{,}ℝ\right)$ and $v_b\in W_b$ , the implicit function theorem implies that for all $v\in \mathcal{N}\left({h^{\prime }}_p\left(v_b\right)\right)$ (the null space of ${h^{\prime }}_p\left(v_b\right)$ ), there exist a $C^1$ -map $g\mathrm{<mo>\; :\; </mo>}\left[\mathrm{0,1}\right]\to W_b$ such that $f\left(0\right)=v_a$ and $g^{\prime }\left(0\right)=v$ . Now, we prove $\langle F\left(v_b\right)\mathrm{,}v\rangle =0$ for all $v\in \mathcal{N}\left({h^{\prime }}_p\left(v_b\right)\right)$ . Indeed, for every $t>0$ , $f\left(t\right)=v_b+v+o\left(t\right)\in W_b$ , where $\frac{o\left(t\right)}{t}\to 0$ as $t\to 0$ . Let

$\Phi \left(r\right)=F\left(v_b+r\left(tv+o\left(t\right)\right)\right)$

By Lemma 2.4

$\begin{array}{c}{\Phi }^{\prime }\left(r\right)=\underset{\delta \to 0}{\lim }\frac{\Phi \left(r+\delta \right)-\Phi \left(r\right)}{\delta }\\ =\underset{\delta \to 0}{\lim }\frac{F\left(v_a+\left(r+\delta \right)\left(tv+o\left(t\right)\right)\right)-F\left(v_a+r\left(tv+o\left(t\right)\right)\right)}{\delta }\\ =\langle F^{\prime }\left(v_b+r\left(tv+o\left(t\right)\right)\right)\mathrm{,}tv+o\left(t\right)\rangle \mathrm{.}\end{array}$

Hence there exist an $\theta \in \left(\mathrm{0,1}\right)$ such that

$\begin{array}{l}F\left(v_b+tv+o\left(t\right)\right)-F\left(v_b\right)=\langle F^{\prime }\left(v_b+\theta \left(tv+o\left(t\right)\right)\right),tv+o\left(t\right)\rangle \\ =t\langle F^{\prime }\left(v_b+\theta \left(tv+o\left(t\right)\right)\right),v\rangle +t\langle F^{\prime }\left(v_a+\theta \left(tv+o\left(t\right)\right)\right),\frac{o\left(t\right)}{t}\rangle .\end{array}$

Take limit $t\to 0$ , by Lemma 2.4, one has $F^{\prime }\left(v_b+\theta \left(tv+o\left(t\right)\right)\right)\rightharpoonup F^{\prime }\left(v_b\right)$ weakly. It follows that $\langle F^{\prime }\left(v_b+\theta \left(tv+o\left(t\right)\right)\right)\mathrm{,}v\rangle \to \langle F^{\prime }\left(v_a\right)\mathrm{,}v\rangle$ and

$\left\{F^{\prime }\left(v_b+\theta \left(tv+o\left(t\right)\right)\right)\right\}$ is bounded. Since $\frac{o\left(t\right)}{t}\to 0$ as $t\to 0$ , we have $\langle {F^{\prime }}_b\left(v_b+\theta \left(tv+o\left(t\right)\right)\right)\mathrm{,}\frac{o\left(t\right)}{t}\rangle \to 0$ .

Since $F\left(v_b\right)={\omega }_b$ , one has

$\begin{array}{l}0\le F\left(v_b+tv+o\left(t\right)\right)-F\left(v_b\right)\\ =t\langle F^{\prime }\left(v_b+\theta \left(tv+o\left(t\right)\right)\right)\mathrm{,}v\rangle +t\langle F^{\prime }\left(v_b+\theta \left(tv+o\left(t\right)\right)\right)\mathrm{,}\frac{o\left(t\right)}{t}\rangle \mathrm{.}\end{array}$

Hence

$0\le \langle {F^{\prime }}_b\left(v_a+\theta \left(tv+o\left(t\right)\right)\right)\mathrm{,}v\rangle +\langle F^{\prime }\left(v_a+\theta \left(tv+o\left(t\right)\right)\right)\mathrm{,}\frac{o\left(t\right)}{t}\rangle \mathrm{.}$

Take limit $t\to 0$ , we get $\langle F^{\prime }\left(v_b\right)\mathrm{,}v\rangle \ge 0$ . By arbitrariness of v, one has $\langle F^{\prime }\left(v_b\right)\mathrm{,}-v\rangle \ge 0$ . It follows that $\langle F^{\prime }\left(v_b\right)\mathrm{,}v\rangle =0$ , for every $v\in \mathcal{N}\left({h^{\prime }}_p\left(v_b\right)\right)$ . Set $v^{\prime }\in H_V^1\left({ℝ}^N\right)$ be such that $\langle {h^{\prime }}_p\left(v_a\right),v^{\prime }\rangle =1$ , for every $\phi \in H_V^1\left({ℝ}^N\right)$ , let

$\psi =\phi -\langle {h^{\prime }}_p\left(v_b\right)\mathrm{,}\phi \rangle v^{\prime }\mathrm{.}$

Then $\psi \in \mathcal{N}\left({h^{\prime }}_p\left(v_b\right)\right)$ . It means $\langle F^{\prime }\left(v_b\right)\mathrm{,}\psi \rangle =0$ , i.e.

$\langle F^{\prime }\left(v_b\right)\mathrm{,}\phi \rangle =\langle F^{\prime }\left(v_b\right)\mathrm{,}{v}^{\prime }\rangle \langle {h^{\prime }}_p\left(v_a\right)\mathrm{,}\phi \rangle \mathrm{.}$

Put $\lambda =\lambda \left(b\right)=\langle F^{\prime }\left(v_b\right)\mathrm{,}{v}^{\prime }\rangle$ , we have

$\langle F^{\prime }\left(v_b\right)\mathrm{,}\phi \rangle =\lambda \langle {h^{\prime }}_p\left(v_b\right)\mathrm{,}\phi \rangle \mathrm{,}$

namely,

$\begin{array}{l}{\displaystyle {\int }_{{ℝ}^N}}\nabla v_b\nabla \phi \text{d}x+{\displaystyle {\int }_{{ℝ}^N}}\frac{V\left(x\right)G^{-1}\left(v_b\right)}{g\left(G^{-1}\left(v_b\right)\right)}\phi \text{d}x\\ =\lambda {\displaystyle {\int }_{{ℝ}^N}}\left[{\left|x\right|}^{-\mu }\ast {\left|G^{-1}\left(v_b\right)\right|}^p\right]\frac{{\left|G^{-1}\left(v_b\right)\right|}^{p-2}{G}^{-1}\left(v_b\right)}{g\left(G^{-1}\left(v_b\right)\right)}\phi \text{d}x\\ +\beta {\displaystyle {\int }_{{ℝ}^N}}\frac{{\left|G^{-1}\left(v_b\right)\right|}^{q-1}{G}^{-1}\left(v_b\right)}{g\left(G^{-1}\left(v_b\right)\right)}\phi \text{d}x.\end{array}$

It implies that $u_b=G^{-1}\left(v_b\right)\in M_b$ is a weak solution of Equation (1.1). Moreover, the maximum principle implies $u_b>0$ . $■$

Acknowledgements

This work is supported by the Natural Science Foundation of China (11961081).

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Pekar, S. (1945) Untersuchungen über die Elektronentheorie der Kristale. Akademie Verlag, Berlin.
[2]	Lieb, E.H. (1977) Existence and Uniqueness of the Minimizing Solution of Choquard’s Nonlinear Equation. Studies in Applied Mathematics, 57, 93-105. https://doi.org/10.1002/sapm197757293
[3]	Penrose, R. (1996) On Gravity’s Role in Quantum State Reductions. General Relativity and Gravitation, 28, 581-600. https://doi.org/10.1007/BF02105068
[4]	Lions, P.L. (1980) The Choquard Equation and Related Questions. Nonlinear Analysis: Theory, Methods & Applications, 4, 1063-1072. https://doi.org/10.1016/0362-546X(80)90016-4
[5]	Ma, L. and Zhao, L. (2010) Classification of Positive Solitary Solutions of the Nonlinear Choquard Equation. Archive for Rational Mechanics and Analysis, 195, 455-467. https://doi.org/10.1007/s00205-008-0208-3
[6]	Cingolani, S., Clapp, M. and Secchi, S. (2012) Multiple Solutions to a Magnetic Nonlinear Choquard Equation. Zeitschrift für Angewandte Mathematik und Physik, 63, 233-248. https://doi.org/10.1007/s00033-011-0166-8
[7]	Moroz, V. and Van Schaftingen, J. (2013) Groundstates of Nonlinear Choquard Equations: Existence, Qualitative Properties and Decay Asymptotics. Journal of Functional Analysis, 265, 153-184. https://doi.org/10.1016/j.jfa.2013.04.007
[8]	Moroz, V. and Van Schaftingen, J. (2015) Existence of Groundstates for a Class of Nonlinear Choquard Equations. Transactions of the American Mathematical Society, 367, 6557-6579. https://doi.org/10.1090/S0002-9947-2014-06289-2
[9]	Clapp, M. and Salazar, D. (2013) Positive and Sign Changing Solutions to a Nonlinear Choquard Equation. Journal of Mathematical Analysis and Applications, 407, 1-15. https://doi.org/10.1016/j.jmaa.2013.04.081
[10]	Ghimenti, M. and Van Schaftingen, J. (2016) Nodal Solutions for the Choquard Equation. Journal of Functional Analysis, 271, 107-135. https://doi.org/10.1016/j.jfa.2016.04.019
[11]	Wei, J.C. and Winter, M. (2009) Strongly Interacting Bumps for the Schröding-Newton Equation. Journal of Mathematical Physics, 50, Article ID: 012905. https://doi.org/10.1063/1.3060169
[12]	Moroz, V. and Van Schaftingen, J. (2015) Semi-Classical States for the Choquard Equation. Calculus of Variations and Partial Differential Equations, 52, 199-235. https://doi.org/10.1007/s00526-014-0709-x
[13]	Moroz, V. and Van Schaftingen, J. (2015) Groundstates of Nonlinear Choquard Equations: Hardy-Littlewood-Sobolev Critical Exponent. Communications in Contemporary Mathematics, 17, Article ID: 1550005. https://doi.org/10.1142/S0219199715500054
[14]	Cassani, D. and Zhang, J. (2018) Choquard-Type Equations with Hardy-Littlewood-Sobolev Upper-Critical Growth. Advances in Nonlinear Analysis, 8, 1184-1212. https://doi.org/10.1515/anona-2018-0019
[15]	Li, Q.Q., Teng, K. and Zhang, J. (2020) Ground State Solutions for Fractional Choquard Equations Involving Upper Critical Exponent. Nonlinear Analysis, 197, Article ID: 111846. https://doi.org/10.1016/j.na.2020.111846
[16]	Seok, J. (2018) Limit Profiles and Uniqueness of Ground States to the Nonlinear Choquard Equations. Advances in Nonlinear Analysis, 8, 1083-1098. https://doi.org/10.1515/anona-2017-0182
[17]	Li, Q.Q., Teng, K.M., Zhang, J. and Nie, J.J. (2020) An Existence Result for a Generalized Quasilinear Schrödinger Equation with Nonlocal Term. Journal of Functional Spaces, 2020, Article ID: 6430104. https://doi.org/10.1155/2020/6430104
[18]	Chen, S.X. (2013) Existence of Positive Solutions for a Class of Quasilinear Schrödinger Equations on . Journal of Mathematical Analysis and Applications, 405, 595-607. https://doi.org/10.1016/j.jmaa.2013.04.031
[19]	Li, Q.Q. and Wu, X. (2016) Multiple Solutions for Genneralized Quasilinear Schrödinger Equations. Mathematical Methods in the Applied Science, 40, 1359-1366. https://doi.org/10.1002/mma.4050
[20]	Shen, Y.T. and Wang, Y.J. (2013) Soliton Solutions for Genneralized Quasilinear Schrödinger Equations. Nonlinear Analysis, 80, 194-201. https://doi.org/10.1016/j.na.2012.10.005
[21]	Liu, J.Q., Wang, Y.Q. and Wang, Z.Q. (2003) Soliton Solutions for Quasilinear Schrödinger Equations: II. Journal of Differential Equations, 187, 473-493. https://doi.org/10.1016/S0022-0396(02)00064-5
[22]	Chen, S.X. and Wu, X. (2019) Existence of Positive Solutions for a Class of Quasilinear Schrödinger Equations of Choquard Type. Journal of Mathematical Analysis and Applications, 475, 1754-1777. https://doi.org/10.1016/j.jmaa.2019.03.051
[23]	Liu, J., Wang, Y. and Wang, Z.Q. (2003) Soliton Solutions for Quasilinear Schrödinger Equations: I. Proceedings of the American Mathematical Society, 131, 441-448. https://doi.org/10.1090/S0002-9939-02-06783-7