Hostname: page-component-745bb68f8f-cphqk Total loading time: 0 Render date: 2025-02-06T09:35:44.241Z Has data issue: false hasContentIssue false

On a planar Schrödinger–Poisson system involving a non-symmetric potential

Published online by Cambridge University Press:  05 December 2022

Riccardo Molle
Affiliation:
Dipartimento di Matematica, Università di Roma “Tor Vergata”, Via della Ricerca Scientifica n. 1, 00133 Roma, Italy (molle@mat.uniroma2.it; andreasardilli6@gmail.com)
Andrea Sardilli
Affiliation:
Dipartimento di Matematica, Università di Roma “Tor Vergata”, Via della Ricerca Scientifica n. 1, 00133 Roma, Italy (molle@mat.uniroma2.it; andreasardilli6@gmail.com)
Rights & Permissions [Opens in a new window]

Abstract

We prove the existence of a ground state positive solution of Schrödinger–Poisson systems in the plane of the form

\[ -\Delta u + V(x)u + \frac{\gamma}{2\pi} \left(\log|\cdot| \ast u^2 \right)u = b |u|^{p-2}u \quad\text{in}\ \mathbb{R}^2, \]
where $p\ge 4$, $\gamma,b>0$ and the potential $V$ is assumed to be positive and unbounded at infinity. On the potential we do not require any symmetry or periodicity assumption, and it is not supposed it has a limit at infinity. We approach the problem by variational methods, using a variant of the mountain pass theorem and the Cerami compactness condition.

Type
Research Article
Copyright
© The Author(s), 2022. Published by Cambridge University Press on Behalf of The Edinburgh Mathematical Society

1. Introduction

This work deals with the existence of positive solutions of Schrödinger–Poisson systems of the form

(S)\begin{equation} \left\{ \begin{array}{@{}ll} - \Delta u + V(x)u + \gamma\phi(x)u = b\, |u|^{p-2}u & \text{in } \mathbb{R}^2\\ \Delta \phi = u^2 & \text{in } \mathbb{R}^2, \end{array}\right. \end{equation}

where we assume $b>0$, the nonlinearity $p\ge 4$, the coupling coefficient $\gamma >0$ and the potential $V$ is a positive measurable function.

For their relevance in physics, Schrödinger–Poisson systems have been extensively studied in the three-dimensional case, see, for example, the survey paper [Reference Ambrosetti1], the recent paper [Reference Dutko, Mercuri and Tyler15] or also [Reference Azzollini and Pomponio2, Reference Bellazzini, Jeanjean and Luo3, Reference Cerami and Molle8, Reference Cerami and Molle9, Reference Ruiz25] and the references therein. On the contrary, in dimension two they are much less studied, also due to some difficulties intrinsic to the planar case (see [Reference Bonheure, Cingolani and Secchi4, Reference Cingolani and Weth10, Reference Chen and Pan11, Reference Choquard, Stubbe and Vuffray13, Reference Du and Weth14, Reference Masaki20, Reference Masaki21, Reference Stubbe27] and references therein).

Concerning the semiclassical analysis related to (S), very little is known. The autonomous case has been treated first in [Reference Masaki20], while the case where a potential is involved has been studied only recently, in [Reference Bonheure, Cingolani and Secchi4], where the problem is fronted by perturbation methods.

We observe that in the two-dimensional Choquard problem with logarithmic kernel, analogous difficulties come out. We refer the reader to the recent paper [Reference Cassani and Tarsi7] for more details and references on this subject. In [Reference Cassani and Tarsi7], the authors considered non homogeneous nonlinearities and the positive potential is assumed to be continuous and one-periodic.

Since the Poisson equation in (S) gives

(1.1)\begin{equation} \phi=\phi_u(x) = \displaystyle\frac{1}{2\pi}\displaystyle\int_{\mathbb{R}^2} \log|x-y|\, u^2(y)\, {\rm d}y, \end{equation}

the system (S) can be reduced to the single nonlocal equation

(E)\begin{equation} {-}\Delta u + V(x)u + \displaystyle\frac{\gamma}{2\pi} (\log(|\cdot|) \ast u^2 )u = b |u|^{p-2}u \quad\text{in}\ \mathbb{R}^2. \end{equation}

It is natural to expect that the solutions of equation (E) correspond to the critical points of the following action functional, where we have normalized the constants:

\begin{align*} I(u)& =\frac 12 \int_{\mathbb{R}^2}(|\nabla u|^2+V(x)u^2)\, {\rm d}x\\ & \quad+ \frac 14 \iint_{\mathbb{R}^2\times\mathbb{R}^2}\log(|x-y|)u^2(x)u^2(y)\,{\rm d}x\, {\rm d}y-\frac 1p\int_{\mathbb{R}^2}|u|^p\,{\rm d}x. \end{align*}

However, a first difficulty is that $I$ is not well defined on the usual Sobolev space $H^1(\mathbb {R}^2)$, for the presence of the sign-changing and unbounded logarithmic kernel. Moreover, in the framework considered in this paper, another point to be fixed is related to the potential, that is assumed to be unbounded and possibly not coercive at infinity. To overcome the difficulty related to the Newtonian kernel, Stubbe in [Reference Stubbe27] introduced the following Hilbert space, to study (S) in the autonomous case $V\equiv 0$, $b=0$,

(1.2)\begin{equation} \widetilde X:=\left\{u\in H^1(\mathbb{R}^2)\ :\ \displaystyle\int_{\mathbb{R}^2}[\log (1+|x|)]u^2(x)\, {\rm d}x<\infty\right\}. \end{equation}

The variational framework introduced in [Reference Stubbe27] was then developed in detail by Cingolani and Weth in [Reference Cingolani and Weth10] to front non-autonomous problems of the form (S), with $b>0$, when the potential is assumed to be $\mathbb {Z}^2$-periodic or constant. The Cerami condition turns out to be the more convenient compactness condition to tackle the two-dimensional case in a variational way, indeed it is not evident how to show the boundedness of the Palais–Smale sequences. However, since the logarithmic convolution term does not have a definite sign on $\widetilde X$, also the boundedness of Cerami sequences becomes a major difficulty. The idea, developed in [Reference Stubbe27], to handle the different signs of the kernel is to take advantage of the relation $\log r=\log (1+r)-\log (1+\frac 1r)$, for all $r>0$, to work in $\widetilde X$ and decompose the action of the nonlocal part in its positive and negative contribution. Once obtained the boundedness of the Cerami sequences, another difficulty is to prove that the weak limit actually is a strong limit for the Cerami sequences. To this aim, some more assumptions on the potential has been done; for example, its $\mathbb {Z}^2$ periodicity considered in [Reference Cingolani and Weth10] (see also [Reference Cao, Dai and Zhang6] and references therein), or its axial symmetry, assumed for example in [Reference Chen and Pan11, Reference Chen and Tang12].

To the best of our knowledge, up to the semiclassical analysis performed in [Reference Bonheure, Cingolani and Secchi4], there are not works concerning problem (S) without any periodicity or symmetry assumption on the potential. Our contribution in this paper is to consider a class of problems of the form (S), where $V$ does not enjoy those assumptions. On $V$, we assume the following conditions

(V)\begin{equation} \begin{array}{cl} (a) & V\in L^1_{\mathop{\rm loc}\nolimits}(\mathbb{R}^2)\\ (b) & \inf_{\mathbb{R}^2}V>0\\ (c) & \left|\{x \in \mathbb{R}^2: V(x)\leq M \}\right| <{+} \infty,\ \forall\ M>0. \end{array} \end{equation}

If the potential verifies (V), the functional $I$ could be not well defined on $\widetilde X$, so that its natural domain turns out to be the weighted Hilbert space

(1.3)\begin{equation} X:=\left\{u\in\widetilde X\ :\ \displaystyle\int_{\mathbb{R}^2}V(x)u^2\, {\rm d}x<\infty\right\} \end{equation}

endowed with the norm

\[ \|u\|_X^2:=\int_{\mathbb{R}^2}\left(|\nabla u|^2+V(x)u^2+\log(1+|x|)u^2\right){\rm d}x\quad u\in X. \]

We assume (V)$(a)$ to avoid pathological potentials such that $X=\emptyset$ and to guarantee that ${\mathcal {C}}^\infty _c(\mathbb {R}^2)\subseteq X$, so that the weak formulation of problem (E) is meaningful.

In Proposition 3.2, we employ assumption (V)$(c)$ to get the compactness condition. In this regard, it is worth to observe that even if the coercivity of the logarithmic weight in the definition of $\widetilde X$ guarantees the compactness of the embedding of $\widetilde X$ in the Lebesgue spaces $L^p$, for $p\in [2,+\infty )$, the logarithmic weight in the functional $I$ does not work in the compactness of the functional, because its contribution is invariant by translation.

We prove the following result.

Theorem 1.1 If $p\ge 4,$ and $V$ verifies (V), then equation (E) has at least a non-trivial weak solution $\bar u\in X$, with $\bar u\ge 0$. Moreover, if $p>4$ then the solution $\bar u$ is a ground state solution:

\[ I(\bar u)=\inf\{I(u)\ :\ u\in X\setminus\{0\},\ u\mbox{ is a solution of (E)}\}. \]

If $V\in L^q_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$ for some $q>1,$ then $\bar u\in {\mathcal {C}}^{0,\alpha }_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$ and $\bar u(x)>0,$ for all $x\in \mathbb {R}^2$.

If $V\in L^q_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$ for some $q>2,$ then $\bar u\in {\mathcal {C}}^{1,\alpha }_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$.

If $V\in {\mathcal {C}}^{0,\alpha }_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2),$ then $\bar u\in {\mathcal {C}}^2(\mathbb {R}^2)$ is a classical solution.

To prove Theorem 1.1, we use a variant of the mountain pass Theorem (see Theorem 2.6). Indeed, it is not convenient to try a direct application of the mountain pass Theorem because it is not straightforward to see that $I$ is a positive quadratic form near 0 with respect to the $X$-norm. On the contrary, it is easy to see that $I$ is quadratic in zero with respect to the usual $H^1$-norm. Let us remark that for $p=4$, the analysis of I on the fiber  $t\mapsto t u(\cdot )$ is not sufficient to guarantee the mountain pass structure, so for $p=4$ a different fiber $t\mapsto u_t$ has to be considered. Moreover, the assumption $p>4$ is necessary to apply the fibering techniques that guarantee that the mountain pass solution is a ground state. To apply Theorem 2.6, a basic tool is the Cerami compactness condition, that can be verified because of (V) and $p\ge 4$.

The paper is organized as follows: in § 2, we introduce the main notation and some preliminary results, in § 3, we prove the compactness condition and in § 4, we prove Theorem 1.1.

2. Variational framework and preliminary facts

Since in the present work we focus on (E) in the case $\gamma,b,\inf V > 0$, to simplify the notation, we assume

\[ \gamma = 2\pi,\quad b=1,\ \inf V =1. \]

We use the following notation:

  • $\int =\int _{\mathbb {R}^2}$.

  • $B_r(y)$ denotes the open ball of radius $r>0$ and centre $y\in \mathbb {R}^2$.

  • $L^p$, $1\leq p<+\infty$, and $H^1$ are the usual Lebesgue and Sobolev spaces in $\mathbb {R}^2$, with norm $|u|_{p}=(\int |u|^p{\rm d}x)^{1/p}$ and $\|u\|=(\int (|\nabla u|^2+u^2)\,{\rm d}x)^{1/2}$, respectively.

  • We denote by $C, c_1, c_2,\ldots$ various constants that can also vary from one line to another.

First, in Lemma 2.1, we recall from [Reference Cingolani and Weth10, Reference Stubbe27] that the nonlocal term is well defined and regular in $\widetilde X$. To state this result, let us introduce the following bilinear forms:

\begin{align*} & (u, v) \mapsto B_1(u, v):= \iint \log(1+ |x-y|)u(x)v(y)\, {\rm d}x\,{\rm d}y,\\ & (u, v)\mapsto B_2(u,v):= \iint \log\left(1+\frac{1}{|x-y|}\right)u(x) v(y)\, {\rm d}x\,{\rm d}y,\\ & (u, v)\mapsto B_0(u, v):= B_1(u, v) - B_2(u, v) = \iint \log(|x-y|)u(x) v(y)\, {\rm d}x\, {\rm d}y, \end{align*}

where $u, v:\mathbb {R}^2 \to \mathbb {R}$ are measurable functions such that the integrals are well defined. Correspondingly, let

\[ V_1(u) := B_1(u^2, u^2),\quad V_2(u):= B_2(u^2, u^2),\quad V_0(u):= B_0(u^2, u^2) = V_1(u) - V_2(u). \]

With this notation, the functional related to problem (E) takes the form

\[ I(u) = \frac12 \int (|\nabla u|^2 + V(x)u^2)\ {\rm d}x + \frac14 V_0(u) - \frac 1p |u|_p^p. \]

Lemma 2.1 The functional $I$ is a well-defined $C^2$ functional on the Hilbert space X. Moreover, $V_2$ is of class $C^1$ on $L^\frac 83 (\mathbb {R}^2)$ and the following formulas hold

(2.1)\begin{equation} V'_i(u)[v] = 4B_i(u^2, uv), \quad \forall u, v \in X,\ i=0, 1, 2. \end{equation}

In particular,

\[ V_i'(u)[u]=4V_i(u), \quad \forall u \in X,\ i=0, 1, 2. \]

For the proof, see [Reference Cingolani and Weth10, Lemma 2.2] and [Reference Sardilli26, Lemma 3.1.5].

Moreover, let us recall the Hardy–Littlewood–Sobolev inequality (see [Reference Lieb and Loss18, Theorem 4.3]):

(2.2)\begin{equation} \left|\iint_{\mathbb{R}^{N}\times \mathbb{R}^N} \displaystyle\frac {f(x) h(y)}{|x-y|^\lambda}\, {\rm d}x{\rm d}y\right|\leq C(\lambda, p )|f|_p|h|_{r}, \quad \forall f \in L^{p}(\mathbb{R}^N),\ \forall g \in L^r(\mathbb{R}^N), \end{equation}

where $p,r>1$ and $0 <\lambda < N$ are such that $\frac {1}{p} + \frac {\lambda }{N} + \frac {1}{r} = 2$. Then, by (2.2) and $\log (1+r)\le r$, $\forall r>0$, it is readily seen that

(2.3)\begin{equation} V_2(u)\le \iint \displaystyle\frac{u^2(x)u^2(y)}{|x-y|} \, {\rm d}x\, {\rm d}y\le c\, |u^2|_{4/3} |u^2|_{4/3}=c|u|^4_{8/3}. \end{equation}

Next, we introduce the Cerami compactness condition.

Definition 2.2 Let $X$ be a Banach space, $I\in C^1(X)$ and $c \in \mathbb {R}$. A sequence $\{x_n\}_n$ is called a $(C)$ sequence for $I$ at the level $c$ if

(C)\begin{equation} {I(x_n) \to c\ \text{as}\ n\to\infty\ \text{and} \ } {\|I'(x_n)\|_{\mathcal{L}}(1 + \|x_n\|_X) \to 0,\ \text{as}\ n\to\infty.} \end{equation}

We say that $I$ satisfies the Palais–Smale–Cerami condition (at the level $c$), $(C)$ condition for short, if any $(C)$ sequence (at the level $c$) possesses a converging subsequence.

To analyse the compactness condition, we report some known facts: Lemma 2.3 is Lemma 2.6 in [Reference Cingolani and Weth10] and states a continuity property of the bilinear form $B_1$, Lemma 2.4 provides the compactness of the embedding of $X$ in the Lebesgue spaces.

Lemma 2.3 Let $\{u_n\}_n, \{v_n\}_n, \{w_n\}_n$ be bounded sequences in X such that $u_n \to u$ weakly in X. Then, for every $z\in X$, we have $B_1(v_nw_n, z(u_n - u)) \to 0$ as $n\to \infty$.

Lemma 2.4 The space $X$ is compactly embedded in $L^q(\mathbb {R}^2)$ for all $q \in [2, +\infty )$.

Since bounded sets in $X$ are also bounded sets in $\widetilde X$ (see (1.2), (1.3)), to verify Lemma 2.4 it is sufficient to exploit the coercivity of the map $|x|\mapsto \log (1+|x|)$ and apply, for example, Theorem XIII.65 in [Reference Reed and Simon24] (see also [Reference Rabinowitz23]).

Finally, let us state the topological tools we will use in the proof of Theorem 1.1.

Lemma 2.5 see [Reference Li and Wang17, Lemma 2.6]

Let $X$ be a Banach space, $I \in C^{1}(X),$ $c \in \mathbb {R},$ $\varepsilon >0$ and suppose that there exists $\alpha >0$ such that

(2.4)\begin{equation} (1 + \|x\|_X) \|I'(x)\|_{\mathcal{L}} \geq \alpha > 0 \quad \forall x \in I^{c+2\epsilon}_{c-2\epsilon}. \end{equation}

Then there exists a continuous map $\eta :X \to X$, such that:

  1. 1. $\eta (I^{c+\epsilon }) \subset (I^{c-\epsilon });$

  2. 2. $\eta (x) = x,$ for all $x \not \in I^{c+2\epsilon }_{c-2\epsilon },$

where $I^a=\{x\in X\ :\ I(x)\le a\}$ and $I^b_a=\{x\in X\ :\ a\le I(x)\le b\}$, for all $a< b$.

From Lemma 2.5 there follows the following variant of the Mountain Pass Theorem.

Theorem 2.6 Let $I \in C^{1}(X)$, X a Banach space. We suppose that $I(0) = 0$ and there exists a closed subspace $S \subset X$ that disconnects $X$ in two path-connected components $X_1$ and $X_2$. We assume that $0 \in X_1$ and there exists $A>0$, $\bar u \in X_2$ such that

  • $I\big |_S\geq A;$

  • $I(\bar u)\leq 0$.

Set

(2.5)\begin{equation} c = \inf_{\gamma \in \Gamma} \max_{t \in [0,1]} I(\gamma(t)) \end{equation}

where $\Gamma :=\{\gamma : [0,1] \to X: \gamma\mbox{ is continuous,}\ \gamma (0)=0,\ \gamma (1) \in X_2,\ I(\gamma (1))\leq 0 \}$. If $I$ verifies the (C) condition at $c,$ then $c$ is a critical value for $I$.

The proof of Theorem 2.6 is standard so we give here only a sketch of it, for the sake of completeness.

If, by contradiction, $c$ is not a critical value for $I$, then (2.4) has to be verified for suitable $\varepsilon,\alpha >0$, because $I$ verifies the (C) condition at the level $c$. We can choose $\varepsilon < A$. Now, let $\gamma _\varepsilon \in \Gamma$ be such that $\max _{[0,1]}I(\gamma _\varepsilon (t))\le c+\varepsilon$. If we consider $\tilde \gamma _\varepsilon :=\eta \circ \gamma _\varepsilon$, where $\eta$ is the deformation provided by Lemma 2.5, it turns out that

\[ \tilde\gamma_\varepsilon\in\Gamma,\quad \max_{[0,1]}I(\tilde \gamma_\varepsilon(t))\le c-\varepsilon, \]

contrary to the definition of $c$.

3. Compactness

This section is devoted to prove that the (C) condition holds at positive values, that is the range where we are looking for critical levels:

Remark 3.1 Let $u$ be a nontrivial critical point of $I$, then $I(u)>0$.

Indeed, by $p\ge 4$ and assumption (V)$(b)$,

\begin{align*} I(u) & =I(u) -\displaystyle\frac{1}{4}I'(u)[u] = \frac{1}{4}\displaystyle\int (|\nabla u|^2{\rm d}x + V(x)u^2)\, {\rm d}x +\left(\frac 14- \frac{1}{p}\right)|u|^p_p \\ & \geq\displaystyle\frac{1}{4} \displaystyle\int ( |\nabla u|^2 + V(x)u^2)\,{\rm d}x \geq C\|u\|^2>0. \end{align*}

Proposition 3.2 Assume that $p\ge 4$ and $V$ satisfies (V), and let $\{u_n\}_n$ be a $(C)$ sequence at the level $c$. If $c>0$ then $\{u_n\}_n$ is relatively compact.

A key point to prove the proposition is the following lemma that states the boundedness of the $(C)$ sequences at the positive levels.

Lemma 3.3 Under the assumptions of Proposition 3.2, the sequence $\{u_n\}_n$ is bounded in $X$. Moreover, there exist a ball $B_2(\bar x)$ and a constant $d>0$ such that, up to a subsequence,

(3.1)\begin{equation} \displaystyle\int_{B_2(\bar x) }u_n^2(x)\, {\rm d}x \ge d \quad \forall n\in\mathbb{N}. \end{equation}

Proof. Since $\{u_n\}_n$ is a (C) sequence, it turns out that

(3.2)\begin{equation} I'(u_n)[u_n] = o(1) \end{equation}

and so we have

\begin{align*} c+o(1) & =I(u_n) -\displaystyle\frac{1}{4}I'(u_n)[u_n] \\ & =\displaystyle\frac{1}{4} \displaystyle\int ( |\nabla u_n|^2 + V(x)u_n^2)\,{\rm d}x +\left(\frac 14-\frac 1p\right)\, |u_n|_p^p\\ & \geq C \|u_n\|^2. \end{align*}

Therefore, $\{u_n\}_n$ is bounded in $H^1(\mathbb {R}^2)$, hence in $L^q(\mathbb {R}^2)$ for all $q\in [2, +\infty )$, and

(3.3)\begin{equation} \displaystyle\int (|\nabla u_n|^2 + V(x)u_n^2) \,{\rm d}x\le C\quad\forall n\in\mathbb{N}. \end{equation}

So, we are left to prove that

\[ |u_n|^2_{{\ast}}:=\int \log(1 + |y|)u^2_n(x)\, {\rm d}x< C,\quad \forall n\in\mathbb{N}. \]

Let us define

\[ \delta := \limsup_{n\to\infty} \left(\max_{y\in\mathbb{R}^2} \int_{B_1(y)}|u_n|^2{\rm d}x\right). \]

We claim that $\delta >0$. Assume, by contradiction, that $\delta = 0$. Then by [Reference Lions19, Lemma I.1] there follows

(3.4)\begin{equation} u_n \to 0\text{ in }\ L^{q}(\mathbb{R}^2),\quad\forall q\in (2,+\infty) \end{equation}

and so, by (2.3), we get

(3.5)\begin{equation} V_2(u_n) \to 0 \quad \text{ as }\quad n\to+\infty. \end{equation}

From (3.2), (3.4) and (3.5), we infer

\[ \int \left(|\nabla u_n|^2 + V(x)u_n^2\right){\rm d}x + V_1(u_n) = I'(u_n)[u_n] + V_2(u_n) + |u_n|^p_p=o(1), \]

that, by the positivity of $V_1$, implies

(3.6)\begin{equation} \displaystyle\int (|\nabla u_n|^2 + V(x)u_n^2)\,{\rm d}x=o(1)\quad \mbox{ and }\quad V_1(u_n)=o(1). \end{equation}

Summarizing, from (3.4), (3.5) and (3.6), we infer

\[ I(u_n) = \frac{1}{2}\int (|\nabla u_n|^2 + V(x)u^2_n)\,{\rm d}x + \frac{1}{4}V_1(u_n) - \frac{1}{4} V_2(u_n) - \frac{1}{p}|u_n|^p_p =o(1), \]

contrary to $I(u_n)\to c>0$, and the claim follows.

Then, $\delta > 0$ and there exists a sequence $\{x_n\}_n \subset \mathbb {R}^2$ such that

(3.7)\begin{equation} \displaystyle\int_{B_1(x_n)} u^2_n(x)\, {\rm d}x > \displaystyle\frac{\delta}{2}, \end{equation}

for n large enough.

The task is now to prove that $\{x_n\}$ is bounded. Assume by contradiction that $|x_n| \to \infty$, up to a subsequence.

Let us fix an arbitrary $M>0$ and call

\[ A_n := \{x \in B_1(x_n)\ : \ V(x)\leq M\}, \quad B_n := \{x \in B_1(x_n)\ : \ V(x)>M\}. \]

By (V) $(c)$ and $|x_n| \to \infty$, we obtain

(3.8)\begin{equation} \left|A_n\right| \leq \left|\{x \in \mathbb{R}^2\ :\ |x|\geq|x_n|-2, \ V(x)\leq M\}\right| \mathop{\longrightarrow} 0,\quad\mbox{ as }n\to\infty. \end{equation}

Observe that from (3.8) we infer, for any fixed $q>2$,

\[ \int_{A_n}u_n^2(x)\, {\rm d}x\le |u_n|_q^2|A_n|^{1-\frac 2q}=o(1). \]

As a consequence

\[ \int_{B_n}u_n^2(x)\, {\rm d}x =\int_{B_1(x_n)}u_n^2(x)\,{\rm d}x-\int_{A_n}u_n^2(x)\,{\rm d}x\ge \frac\delta 2+o(1). \]

Then, taking into account (2.3), we have

(3.9)\begin{equation} \begin{aligned} c+o(1)= I(u_n) & \geq \displaystyle\frac{1}{2} \displaystyle\int_{B_n}V(x)u_n^2{\rm d}x - \frac{1}{4}V_2(u_n) - C_1\frac{1}{p}\|u_n\|^p \\ & \geq \frac\delta 4\, M - C_2, \end{aligned} \end{equation}

where $C_2$ is a constant independent of $M$. Letting $M\to \infty$ in (3.9), we get a contradiction, hence $\{x_n\}$ has to be bounded.

Therefore, there exists $\bar {x}\in \mathbb {R}^2$ such that $x_n\to \bar {x}$, up to a subsequence. Taking into account (3.7), we can also assume that

(3.10)\begin{equation} \displaystyle\int_{B_2(\bar{x})}u_n^2{\rm d}x \geq \displaystyle\frac{\delta}{2} > 0, \end{equation}

that proves (3.1).

Now, observe that for every $R>0$

\[ 1 + |x - y|\geq 1 +\frac{|y|}{2}\geq \sqrt{1 + |y|} \quad \forall y\in \mathbb{R}^2\setminus B_{2R}(0),\ \forall x\in B_R(0). \]

Then, fixing $\bar R>|\bar x|+2$ and taking into account (3.10), we get

(3.11)\begin{equation} \begin{aligned} V_1(u_n) =B_1(u_n^2,u_n^2) & \ge \displaystyle\int_{\mathbb{R}^2\setminus B_{2\bar R}(0)}\left(\int_{B_{\bar R}(0)} \log(1 + |x - y|)u^2_n(x)\, {\rm d}x\right)u^2_n(y)\, {\rm d}y \\ & \geq \left(\displaystyle\int_{B_2(\bar x)}u_n^2(x){\rm d}x\right)\cdot \left(\int_{\mathbb{R}^2\setminus B_{2\bar R}(0)} \log\left(1 + \displaystyle\frac{|y|}{2}\right)u^2_n(y)\, {\rm d}y\right) \\ & \geq \frac{\delta}{4}\int_{\mathbb{R}^2\setminus B_{2\bar R}(0)} \log(1 + |y|)u^2_n(y)\, {\rm d}y\\ & = \frac{\delta }{4}\left(|u_n|^2_{{\ast}} - \int_{B_{2\bar R}(0)} \log(1 + |y|)\, u^2_n(y)\, {\rm d}y\right) \\ & \geq \frac{\delta}{4}(|u_n|^2_{{\ast}} - \log (1+2\bar R) |u_n|^2_2). \end{aligned} \end{equation}

Thus, we conclude that

(3.12)\begin{equation} |u_n|^2_{{\ast}} \leq \displaystyle\frac{4}{\delta}V_1(u_n) + C |u_n|^2_2. \end{equation}

Since $\{|u_n|_2\}_n$ is bounded, if we prove that $\{V_1(u_n)\}_n$ is bounded, we are done.

Observe that

(3.13)\begin{equation} {-}\displaystyle\frac{1}{4}V_1(u_n) + \frac{1}{4} V_2(u_n) + \left(\frac{1}{2} - \frac{1}{p}\right)|u|^p_p = I(u_n) -\frac{1}{2}I'(u_n)[u_n] = c + o(1), \end{equation}

from which we infer, taking into account (2.3),

(3.14)\begin{equation} \displaystyle\frac{1}{4}V_1(u_n) = \frac{1}{4}V_2(u_n) + \left(\frac{1}{2} - \frac{1}{p}\right)|u_n|^p_p -c +o(1) \leq C|u_n|^4_{ {8}/{3}} + C|u_n|^p_p -c+o(1)\leq C, \end{equation}

that is the desired result.

Proof of Proposition 3.2. By Lemma 3.3, $\{u_n\}_n$ is bounded in the Hilbert space $X$, so, taking also into account Lemma 2.4, there exists $\bar u\in X$ such that, up to a subsequence,

(3.15)\begin{equation} u_n \to \bar{u}\text{ weakly in } X,\mbox{ in }\ L^s(\mathbb{R}^2), \text{ for every }s\in[2,+\infty), \mbox{ and a.e. in }\mathbb{R}^2. \end{equation}

Clearly, the sequence $\{u_n\}$ is bounded also in the weighted Hilbert space

\[ H_V:=\left\{u\in H^1(\mathbb{R}^2)\ :\ \int V(x)u^2(x)\, {\rm d}x<\infty\right\},\quad\|u\|_V^2:=\int(|\nabla u|^2+V(x)u^2)\, {\rm d}x. \]

Hence, we can assume that $\{u_n\}$ weakly converges in $H_V$ to a function $\tilde u\in H_V$. Observe that $\tilde u=\bar u$. Indeed for every fixed $\varphi \in {\mathcal {C}}^\infty _0(\mathbb {R}^2)$ the map

\[ u\mapsto \int u\varphi\, {\rm d}x \]

is a continuous linear form on $X$ and on $H_V$, so that

\[ \lim_{n\to\infty}\int u_n\varphi\, {\rm d}x=\int\bar u\varphi\, {\rm d}x \quad\mbox{ and }\quad \lim_{n\to\infty}\int u_n\varphi\, {\rm d}x=\int\tilde u\varphi\, {\rm d}x. \]

Then we can conclude that

\[ \int (\bar u-\tilde u)\varphi\, {\rm d}x =0\quad\forall \varphi\in{\mathcal{C}}^\infty_0(\mathbb{R}^2) \]

and the assertion follows by the fundamental lemma of Calculus of Variations (see e.g. [Reference Brezis5, Corollary 4.24]).

Now, from the boundedness in $X$ of the (C) sequence $\{u_n\}_n$, we get

(3.16)\begin{equation} \begin{aligned} o(1) & = I'(u_n)[u_n - \bar{u}]= \displaystyle\int (|\nabla u_n|^2 + V(x)u_n^2)\, {\rm d}x - \int ( \nabla u_n \cdot \nabla \bar{u}+ V(x)u_n\bar{u})\,{\rm d}x \\ & \quad+ V'_1(u_n)(u_n - \bar{u}) - V'_2(u_n)(u_n - \bar{u}) - \displaystyle\int |u_n|^{p-2}u_n(u_n - \bar{u})\, {\rm d}x.\\ \end{aligned} \end{equation}

The weak convergence in $H_V$ yields

(3.17)\begin{equation} \liminf_{n\to\infty} \displaystyle\int (|\nabla u_n|^2 + V(x)u_n^2)\, {\rm d}x \geq \int (|\nabla \bar{u}|^2 + V(x)\bar{u}^2)\, {\rm d}x. \end{equation}

By Lemmas 2.1 and 2.3,

(3.18)\begin{equation} \begin{aligned} V'_1(u_n)[u_n - \bar{u}] & = 4B_1(u_n^2, u_n(u_n - \bar{u})) = 4B_1(u_n^2, (u_n - \bar{u})^2) + 4B_1(u_n^2, \bar{u}(u_n - \bar{u}))\\ & =4B_1(u_n^2, (u_n - \bar{u})^2)+o(1). \end{aligned} \end{equation}

By (2.2) and (3.15)

(3.19)\begin{equation} \begin{aligned} |V'_2(u_n)[u_n - \bar{u}]| & = |4B_2(u_n^2, u_n(u_n - \bar{u}))| \\ & = 4\left| \iint \log\left(1 + \displaystyle\frac{1}{|x - y|}\right)u^2_n(x)u_n(y)(u_n - \bar{u})(y)\,{\rm d}x\,{\rm d}y\right|\\ & \leq 4\iint \frac{u_n^2(x)\,|u_n(y)|\, |(u_n - \bar{u})(y)|}{|x-y|}\, {\rm d}x\,{\rm d}y \\ & \leq 4 |u_n|^3_{\frac{8}{3}}|u_n - \bar{u}|_{\frac{8}{3}}=o(1). \end{aligned} \end{equation}

By (3.15),

(3.20)\begin{equation} \left|\displaystyle\int |u_n|^{p-2}u_n(u_n - \bar{u})\, {\rm d}x\right|\le |u_n|_p^{p-1}|u_n-u|_p=o(1). \end{equation}

From (3.16), taking into account (3.17), (3.18), (3.19) and (3.20), we infer

\begin{align*} o(1) & = \int \left(|\nabla u_n|^2 + V(x)u_n^2\right){\rm d}x - \int \left(|\nabla \bar u |^2 + V(x)\bar u^2\right){\rm d}x + 4B_1(u_n^2, (u_n - \bar{u})^2) + o(1)\\ & \geq o(1). \end{align*}

Thus, it follows that

(3.21)\begin{equation} \displaystyle\int(|\nabla u_n|^2 + V(x)u_n^2)\,{\rm d}x \to \int(|\nabla \bar{u}|^2 + V(x)\bar{u}^2)\,{\rm d}x, \end{equation}

and

(3.22)\begin{equation} B_1(u_n^2, (u_n - \bar{u})^2) \to 0, \text{ as } n\to+\infty. \end{equation}

Taking into account (3.1) and arguing exactly as in (3.11), we get

(3.23)\begin{equation} B_1(u_n^2, (u_n - \bar{u})^2) \ge\displaystyle\frac{d}{2}(|u_n - \bar{u}|^2_*-C|u_n - \bar{u}|_2^2). \end{equation}

From (3.23), (3.22) and $u_n\to \bar u$ in $L^2$, there follows

\[ |u_n-\bar u|_*\to 0, \]

that, together with (3.21), gives $\|u_n-\bar u\|_X\to 0$.

4. Proof of Theorem 1.1

To prove Theorem 1.1, we are going to apply Theorem 2.6. Since condition $(C)$ holds on $(0,+\infty )$ by Proposition 3.2, we are left to show that also the geometric conditions hold.

Observe that, by (2.3) and by the Sobolev inequality, for every $u \in X \setminus \{0\}$

(4.1)\begin{equation} \begin{aligned} I(u) & \ge \displaystyle\frac{1}{2} \displaystyle\int_{\mathbb{R}^2} (|\nabla u|^2 + V(x)u^2)\, {\rm d}x - \frac{1}{4}\iint \log\left(1 + \frac{1}{|x-y|}\right)u^2(x)u^2(y)\, {\rm d}x\,{\rm d}y - \frac{1}{p}|u|^p_p\\ & \geq c_1 \|u\|^2 - c_2|u|^4_{\frac{8}{3}} - c_3\|u\|^p \\ & \geq c_1\|u\|^2 - c_4\|u\|^4 - c_3\|u\|^p. \end{aligned} \end{equation}

Since $p\ge 4$, there exist $A >0$ and $\bar {\rho }>0$ such that $I_{|_{S }}\geq A$, where $S :=\{u \in X: \|u\| = \bar \rho \}$. The closed set $S$ disconnects $X$ in the two path-connected components

\[ X_1:=\{u \in X: \|u\|<\bar{\rho} \}\quad\mbox{ and }\quad X_2:= \{u \in X: \|u\|>\bar{\rho} \}. \]

Clearly, $0 \in X_1$ and $I(0)=0$.

Now, we are going to show that there exists $\bar u\in X_2$ such that $I(\bar u)<0$. If $p>4$, let us fix $\tilde u\in X\setminus \{0\}$ and observe that $\lim \limits _{t\to \infty }I(t\tilde u)=-\infty$, so the claim follows.

If $p=4$, then let us fix $\tilde u\in {\mathcal {C}} (\mathbb {R}^2)$ with compact support and define $u_t:=t^2\tilde u(t\,\cdot )$, for $t>0$. We can assume that $0$ is a Lebesgue point for $V$, then, by standard computations,

(4.2)\begin{equation} \begin{aligned} I(u_t) & =\left[\displaystyle\frac 12 |\nabla \tilde u|_2^2\right]\, t^4+\left[\frac {V(0)+o(1)}{2}\, |\tilde u|_2^2\right]\, t^2 \\ & \quad+\left[\frac 14 \iint \log(|x-y|)\tilde u^2(x)\tilde u^2(y)\,{\rm d}x\,{\rm d}y\right]\, t^4\\ & \quad-\left[\frac 14 |\tilde u|_2^4\right]\, t^4\log t -\left[\frac 14 |\tilde u|_4^4\right]\, t^6, \end{aligned} \end{equation}

where $o(1)\to 0$ as $t\to \infty$. Since (4.2) implies $\lim \limits _{t\to \infty }I( u_t)=-\infty$, and $\|u_t\|\to \infty$ as $t\to \infty$, the existence of $\bar u\in X_2$ such that $I(\bar u)<0$ follows.

Then, the value

(4.3)\begin{equation} c:=\inf_{\gamma \in \Gamma}\max_{t\in[0,1]} I(\gamma(t)), \end{equation}

where $\Gamma := \{\gamma :[0,1] \to X: \gamma\mbox{ is continuous,}\ \gamma (0)=0, \gamma (1) \in X_2, I(\gamma (1))\leq 0 \}$, is a critical value for $I$, according to Theorem 2.6, with $c\ge A>0$.

Now, let $p>4$. To verify that $c$ is a ground state level, let us fix any nontrivial solution $\tilde u$ and consider the path $t\mapsto t\tilde u$. A direct computation shows that the real function $t\mapsto I( t\tilde u)$ has a unique critical point $t_{\tilde u}$, that corresponds to a maximum. From $\frac {d\,}{dt}I(t\tilde u)=I'(t\tilde u)[\tilde u]$ we infer $t_{\tilde u}=1$. Then, taking into account that $\lim \limits _{t\to \infty }I(t\tilde u)=-\infty$ and (4.3), it is readily seen that

\[ c\le \max_{t\ge 0}I(t\tilde u)=I(\tilde u), \]

that is our claim.

Our next goal is to show that there exists a constant sign solution. Let $u$ be a critical point for $I$ at the level $c$ and consider the path $t\mapsto t|u|$. Since the functional is even, $c=\max _{t\ge 0}I(t|u|)=I(|u|)$ from which we infer that $\bar u:=|u|$ is a critical point for $I$, too. So, $\bar u$ is the weak solution we are looking for.

Now, we consider the regularity issue. For classical regularity results, we refer the reader to [Reference Lieb and Loss18, §10] (in particular Theorems 10.2 and 10.3), or to [Reference Gilbarg and Trudinger16, §8].

Observe that

(4.4)\begin{equation} -\Delta \bar u=\varphi(x)\quad\mbox{ for }\varphi(x)={-}V(x)\bar u+\phi_u(x)\bar u+\bar u^{p-1}, \end{equation}

(see (1.1)).

Assume $V\in L^q_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$ for some $q>1$. Since $\bar u\in L^r(\mathbb {R}^2)$ for every $r\in [2,\infty )$, and $\phi _u\in {\mathcal {C}}^3(\mathbb {R}^2)$ (see [Reference Cingolani and Weth10, Proposition 2.3]), we get $\varphi \in L^s_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$ for every $s\in [1,q)$. So we are in position to apply classical regularity results and obtain $\bar u\in {\mathcal {C}}^{0,\alpha }_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$. Moreover, we can also write

\begin{align*} -\Delta \bar u=\psi(x)\bar u\quad\mbox{ for }\psi(x)={-}V(x) +\phi_u(x) +\bar u^{p-2} \end{align*}

and then, since $\psi \in L^q_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$, by the Harnack inequality, we can conclude $\bar u>0$ (see [Reference Pucci and Serrin22, Theorem 7.2.1]).

If $V\in L^q_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$ for some $q>2$, then we can argue as before and conclude $\varphi \in L^s_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$ for every $s\in [1,q)$, so that $\bar u\in {\mathcal {C}}_{\mathop {\rm loc}\nolimits }^{1,\alpha }(\mathbb {R}^2)$.

If $V\in {\mathcal {C}}^{0,\alpha }_{\mathop {\rm loc}\nolimits }(\mathbb {R}^2)$, then the RHS in (4.4) is locally Hölder continuous, so $u\in {\mathcal {C}}^2(\mathbb {R}^2)$ by standard elliptic regularity (see, for example, [Reference Lieb and Loss18, §10]).

Acknowledgements

R.M. has been supported by the INdAM-GNAMPA group. He acknowledges also the MIUR Excellence Department Project awarded to the Department of Mathematics, University of Rome Tor Vergata, CUP E83C18000100006.

References

Ambrosetti, A., Schrödinger-Poisson systems, Milan J. Math. 76 (2008), 257274.CrossRefGoogle Scholar
Azzollini, A. and Pomponio, A., Ground state solutions for the nonlinear Schrödinger-Maxwell equations, J. Math. Anal. Appl. 345(1) (2008), 90108.CrossRefGoogle Scholar
Bellazzini, J., Jeanjean, L. and Luo, T., Existence and instability of standing waves with prescribed norm for a class of Schrödinger–Poisson equations, Proc. Lond. Math. Soc. 107(3) (2013), 303339.CrossRefGoogle Scholar
Bonheure, D., Cingolani, S. and Secchi, S., Concentration phenomena for the Schrödinger-Poisson system in ℝ2, Discrete Contin. Dyn. Syst. Ser. S 14(5) (2021), 16311648.Google Scholar
Brezis, H., Functional analysis, Sobolev spaces and partial differential equations (New York, Springer, 2011).CrossRefGoogle Scholar
Cao, D., Dai, W. and Zhang, Y., Existence and symmetry of solutions to 2-D Schrödinger-Newton equations, Dyn. Partial Differ. Equ. 18(2) (2021), 113156.CrossRefGoogle Scholar
Cassani, D. and Tarsi, C., Schrödinger-Newton equations in dimension two via a Pohozaev-Trudinger log-weighted inequality, Calc. Var. Partial Differ. Equ. 60(5) (2021), 197.CrossRefGoogle Scholar
Cerami, G. and Molle, R., Positive bound state solutions for some Schrödinger-Poisson systems, Nonlinearity 29(10) (2016), 31033119.CrossRefGoogle Scholar
Cerami, G. and Molle, R., Multiple positive bound states for critical Schrödinger-Poisson systems, ESAIM Control Optim. Calc. Var. 25 (2019), 7–3.CrossRefGoogle Scholar
Cingolani, S. and Weth, T., On the planar Schrödinger-Poisson system, Ann. Inst. H. Poincaré C Anal. Non Linéaire 33(1) (2016), 169197.CrossRefGoogle Scholar
Chen, W. and Pan, H., On the planar axially symmetric Schrödinger-Poisson systems with Choquard nonlinearity, J. Math. Anal. Appl. 504(1) (2021), 125378.CrossRefGoogle Scholar
Chen, S. and Tang, X., Axially symmetric solutions for the planar Schrödinger-Poisson system with critical exponential growth, J. Differ. Equ. 269(11) (2020), 91449174.CrossRefGoogle Scholar
Choquard, P., Stubbe, J. and Vuffray, M., Stationary solutions of the Schrödinger-Newton model – an ODE approach, Differ. Int. Equ. 21(7–8) (2008), 665679.Google Scholar
Du, M., Weth, Tobias T., Ground states and high energy solutions of the planar Schrödinger-Poisson system, Nonlinearity 30(9) (2017), 34923515.CrossRefGoogle Scholar
Dutko, T., Mercuri, C. and Tyler, T. M., Groundstates and infinitely many high energy solutions to a class of nonlinear Schrödinger-Poisson systems, Calc. Var. Partial Differ. Equ. 60(5) (2021), 174.Google Scholar
Gilbarg, D. and Trudinger, N. S., Elliptic partial differential equations of second order (Berlin, Springer-Verlag, 2001).Google Scholar
Li, G. and Wang, C., The existence of a nontrivial solution to a nonlinear elliptic problem of linking type without the Ambrosetti-Rabinowitz condition, Ann. Acad. Sci. Fenn. Math. 36(2) (2011), 461480.CrossRefGoogle Scholar
Lieb, E. H. and Loss, M., Analysis, 2nd edn. Graduate Studies in Mathematics, Volume 14 (American Mathematical Society, Providence, RI, 2001).Google Scholar
Lions, P.-L., The concentration-compactness principle in the calculus of variations. The locally compact case, part 2, Ann. Inst. H. Poincaré Anal. Non Linéaire 1(4) (1984), 223283.CrossRefGoogle Scholar
Masaki, S., Local existence and WKB approximation of solutions to Schrödinger-Poisson system in the two-dimensional whole space, Comm. Partial Differ. Equ. 35(12) (2010), 22532278.CrossRefGoogle Scholar
Masaki, S., Energy solution to a Schrödinger–Poisson system in the two-dimensional whole space, SIAM J. Math. Anal. 43(6) (2011), 27192731.CrossRefGoogle Scholar
Pucci, P. and Serrin, J., The maximum principle, Progress in Nonlinear Differential Equations and their Applications, Volume 73 (Birkhäuser Verlag, Basel, 2007).CrossRefGoogle Scholar
Rabinowitz, P. H., On a class of nonlinear Schrödinger equations, Z. Angew. Math. Phys. 43(2) (1992), 270291.CrossRefGoogle Scholar
Reed, M. and Simon, B., Methods of modern mathematical physics. IV. Analysis of operators (New York-London: Academic Press, 1978).Google Scholar
Ruiz, D., The Schrödinger-Poisson equation under the effect of a nonlinear local term, J. Funct. Anal. 237(2) (2006), 655674.CrossRefGoogle Scholar
Sardilli, A., Master's degree thesis, Roma “Tor Vergata”, 2021.Google Scholar
Stubbe, J., Bound states of two-dimensional Schrödinger–Newton equations, e-print arXiv:0807.4059v1, 2008.Google Scholar