Three Weak Solutions for a Critical Non-Local Problem with Strong Singularity in High Dimension

Abstract

In this paper, we deal with a strongly singular problem involving a non-local operator, a critical nonlinearity, and a subcritical perturbation. We apply techniques from non-smooth analysis to the energy functional, in combination with the study of the topological properties of the sublevels of its smooth part, to prove the existence of three weak solutions: two points of local minimum and a third one as a mountain pass critical point.

1. Introduction

In the present paper, we study a non-local two-parameter problem of the following type:

$$\left(\mathcal{P}\right)\left\{\begin{array}{cc}-\left(a+b{\int }_{\Omega }{|\nabla u|}^{2}dx\right)\Delta u=\mu \frac{s\left(x\right)}{u^{\gamma }}+{\left|u\right|}^{2^{★}-2}u+\lambda f(x,u)\hfill & \mathrm{in}\Omega \hfill \\ u>0\hfill & \mathrm{in}\Omega \hfill \\ u=0\hfill & \mathrm{on}\partial \Omega .\hfill \end{array}\right.$$

We assume that $\Omega$ is a bounded domain in ${\mathbb{R}}^N$ with smooth boundary $\partial \Omega$ ($N\ge 4$), $s:\Omega \mapsto [0,+\infty ]$ is in ${L^{\left(2^{★}\right)}}^{\prime }(\Omega )$, and $f:\Omega \times \mathbb{R}\mapsto [0,+\infty ]$ is a Carathèodory function not identically zero and satisfying $f(x,0)=0$ a.e. in $\Omega$. Moreover, $a,b>0$, $\gamma >1$, and $\lambda ,\mu >0$ are real parameters. As usual, we denote by $2^{★}:=\frac{2N}{N-2}$ the Sobolev critical exponent and by ${\left(2^{★}\right)}^{\prime }=\frac{2^{★}}{2^{★}-1}$ its conjugate.

The peculiarity of the above problem is the combination of a Kirchhoff-type operator, which is responsible for the non-local nature of the problem, of a strong singular term, a critical nonlinearity, and a subcritical perturbation. The existence or multiplicity of solutions to Kirchhoff-type problems with critical terms is frequently studied by mountain pass arguments combined with the Lions concentration-compactness principle, both when $N<4$ (see [1,2,3]) and in the higher dimensions ($N\ge 4$) (see [4,5,6,7]). Note that, in order to employ the concentration-compactness principle, a and b need to satisfy suitable constraints. Moreover, in high dimensions, the effect of a non-local operator, combined with a critical nonlinearity, forces the energy functional to be coercive, and the interplay between the Kirchhoff and the critical term allows us to establish some variational properties that will be crucial to our arguments (see [8]).

After the pioneering work of [9], the interest in singular problems has been increasing over the years; existence and multiplicity results have been obtained, both for the weak singular case (i.e., $\gamma <1$) and for the strong singular case (i.e., $\gamma \ge 1$). Indeed, due to the presence of the singular term, when $\gamma <1$, problem $\left(\mathcal{P}\right)$ is variational, and it is possible to associate with $\left(\mathcal{P}\right)$ an energy functional $E_{\lambda ,\mu }$ which, although not differentiable, has useful properties in the natural Sobolev space $H_0^1(\Omega )$. On the contrary, in the case $\gamma \ge 1$, in general, the energy functional $E_{\lambda ,\mu }$ is no longer defined on the whole $H_0^1(\Omega )$. In the seminal work of [10], it is proved that, even in the semilinear case, the pure singular problem is not expected to possess solutions in $H_0^1(\Omega )$ for $\gamma \ge 3$. However, if the singular term is multiplied by a suitable positive weight, one can still obtain solutions in $H_0^1(\Omega )$ via sub- and super-solution techniques (see [11]).

The presence of a non-local term, in the context of singular problems, makes the analysis more challenging, since it becomes quite difficult to establish any sort of comparison principle (which is one of the main tools to produce solutions in the presence of a strong singularity).

In the present work, we study the interaction between a non-local, a critical, and a singular term to investigate the existence of multiple solutions for $\left(\mathcal{P}\right)$ in $H_0^1(\Omega )$. Note that, if a and b are big enough and $\lambda =0$, problem $\left(\mathcal{P}\right)$ has precisely one solution, which corresponds to the unique global minimizer of $E_{0,\mu }$ (see Remark 3). The presence of a perturbation breaks this uniqueness property, and the existence of three weak solutions is obtained for big $\lambda$s and small $\mu$s. The study of existence of three solutions for singular problems driven by the p-Laplace operator has been developed in [12] in the low-dimensional case and in [13] for any dimension. Later, in [14], an equivalence-type result was obtained in the setting of Orlicz spaces for the non-local case.

In the present work, we extend the results of [13] to a more general problem: with respect to [13], the presence of a non-local operator and a critical term requires some extra variational properties. Thus, even if the underlying idea is the same, our multiplicity result is not straightforward. Employing the results of [8], we are still able to prove that the energy functional $E_{\lambda ,\mu }$ associated with $\left(\mathcal{P}\right)$ is sequentially weakly lower semicontinuous and that its derivative satisfies some form of compactness property.

We will prove the existence of two local minimizers via topological arguments. The idea is to show that the functional $E_{\lambda ,0}$ for $\lambda >0$ is large enough, has two local minimizers, and that the topology of $E_{\lambda ,0}$, after being perturbed by the singular term, changes little enough so that $E_{\lambda ,\mu }$ still has two local minimizers in $H_0^1(\Omega )$, as long as the parameter $\mu >0$ is small enough. The existence of a third solution follows at once by employing a suitable version of the mountain pass theorem for Szulkin functionals.

Let us first present the definition of weak solution in our framework (see [11]).

Definition 1.

A weak solution for $\left(\mathcal{P}\right)$ is a function $u\in H_0^1(\Omega )$ such that

(i) 

$u>0$ almost everywhere in Ω,

(ii) 

$s{u}^{-\gamma }\varphi \in L^1(\Omega )\forall \varphi \in H_0^1(\Omega )$,

(iii) 

$\forall \varphi \in H_0^1(\Omega )$; there holds

$$\left(a+b{\int }_{\Omega }{|\nabla u|}^{2}dx\right){\int }_{\Omega }\nabla u\nabla \varphi dx-\mu {\int }_{\Omega }s{u}^{-\gamma }\varphi dx-{\int }_{\Omega }{u}^{2^{★}-1}\varphi dx-\lambda {\int }_{\Omega }f(x,u)\varphi dx=0.$$

In the subsequent work, we will assume the following:

(f₁) 

there exist $c_1,c_2>0$, $q\in (2,2^{★})$, such that $0\le f(x,t)\le c_1+c_2{t}^{q-1}$ for a.e. $x\in \Omega$, $t\ge 0$;

(f₂) 

${lim}_{t\to 0^{+}}\frac{f(x,t)}{t}=0$ uniformly in $\Omega$;

(s₁) 

there exists $\overline{u}\in C_0^1\left(\overline{\Omega }\right)$, with $\overline{u}>0$, such that $s{\overline{u}}^{-\gamma }\in {L\left(2^{★}\right)}^{\prime }(\Omega )$.

Let $F:\Omega \times [0,+\infty ]\to \mathbb{R}$ be a primitive of f, i.e.,

$$F(x,t)={\int }_0^{t}f(x,s)ds.$$

Also put $f(x,t)=0$ for all $t\le 0$.

In order to state our main result, let us introduce some useful notation. We endow $H_0^1(\Omega )$ with the classical norm $\parallel ·\parallel :={\left({\int }_{\Omega }{|\nabla ·|}^{2}dx\right)}^{\frac{1}{2}}$ and $L^p(\Omega )$ with the standard Lebesgue norm ${\parallel ·\parallel }_p:={\left({\int }_{\Omega }{|·|}^p\right)}^{\frac{1}{p}}$, $p\le 2^{★}$. Let $u_{+}:=max\{u,0\}$ and $u_{-}:=max\{-u,0\}$.

Denote by $S_0$ the embedding constant of $H_0^1(\Omega )\hookrightarrow L^{2^{★}}(\Omega )$, i.e.,

$$S_0=\underset{u\in H_0^1(\Omega )\setminus \left\{0\right\}}{sup}\frac{{\parallel u\parallel }_{2^{★}}^{2^{★}}}{{\parallel u\parallel }^{2^{★}}}.$$

Now, define the constants

$$C_1\left(N\right):=\left\{\begin{array}{ccc}\frac{{4(N-4)}^{\frac{N-4}{2}}}{S_0^{-\frac{N-2}{2}}{N}^{\frac{N-2}{2}}}\hfill & \hfill \mathrm{if}& \hfill N>4\\ S_0\hfill & \hfill \mathrm{if}& \hfill N=4\end{array}\right.,$$

$$C_2\left(N\right):=\left\{\begin{array}{ccc}\frac{2{(N-4)}^{\frac{N-4}{2}}}{S_0^{-\frac{N-2}{2}}}\frac{1}{{\left[\frac{{(N-2)}^2}{N+2}\right]}^{\frac{N-2}{2}}}\hfill & \hfill \mathrm{if}& \hfill N>4\\ 3S_0\hfill & \hfill \mathrm{if}& \hfill N=4.\end{array}\right.$$

The above constants were introduced in [8], where some useful variational properties of the energy functional involving a Kirchhoff and a critical term were proved.

The following result is the main purpose of this work.

Theorem 1.

Assume conditions $\left(f_1\right),\left(f_2\right)$, and $\left(s_1\right)$, and put

$${\lambda }^{★}=inf\left\{\frac{\frac{a}{2}{\parallel u\parallel }^2+\frac{b}{4}{\parallel u\parallel }^4-\frac{1}{2^{★}}{\parallel u_{+}\parallel }_{2^{★}}^{2^{★}}}{{\int }_{\Omega }F(x,u_{+}\left(x\right))dx}:u\in H_0^1(\Omega ),{\int }_{\Omega }F(x,u_{+}\left(x\right))dx>0\right\}.$$

(1)

Then, the following holds true.

(i) 

If $a^{\frac{N-4}{2}}b>C_1\left(N\right)$ then, for each $\lambda >{\lambda }^{★}$, there exists ${\mu }^{★}={\mu }^{★}\left(\lambda \right)>0$ such that, for all $0<\mu <{\mu }^{★}$, problem $\left(\mathcal{P}\right)$ has at least two weak solutions.

(ii) 

If $a^{\frac{N-4}{2}}b>C_2\left(N\right)$ then, for each $\lambda >{\lambda }^{★}$, there exists ${\mu }^{★}={\mu }^{★}\left(\lambda \right)>0$ such that, for all $0<\mu <{\mu }^{★}$, problem $\left(\mathcal{P}\right)$ has at least three weak solutions.

Remark 1.

Note that $C_2\left(N\right)>C_1\left(N\right)$ for $N\ge 4$.

2. Abstract Results

We now state the preparatory results upon which the proof of Theorem 1 is based.

Define $S:\Omega \times \mathbb{R}\mapsto [-\infty ,+\infty ]$ by

$$\begin{array}{c}S(x,t)=\left\{\begin{array}{ccc}\frac{1}{\gamma -1}s\left(x\right)t^{1-\gamma },\hfill & \mathrm{if}\hfill & x\in \Omega \mathrm{and}t>0,\hfill \\ +\infty ,\hfill & \mathrm{if}\hfill & x\in \Omega \mathrm{and}t\le 0.\hfill \end{array}\right.\hfill \end{array}$$

The singular operator $\Psi :H_0^1(\Omega )\to [-\infty ,+\infty ]$ is given by

$$\Psi \left(u\right)=\left\{\begin{array}{ccc}{\int }_{\Omega }S(x,u)dx,\hfill & \mathrm{if}\hfill & S(x,u)\in L^1(\Omega ),\hfill \\ \\ +\infty ,\hfill & \mathrm{if}\hfill & S(x,u)\notin L^1(\Omega ).\hfill \end{array}\right.$$

(2)

Also, we define a functional ${\Phi }_{\lambda }:H_0^1(\Omega )\mapsto \mathbb{R}$ by

$${\Phi }_{\lambda }\left(u\right):=\frac{1}{2}{a\parallel u\parallel }^2+b\frac{1}{4}{\parallel u\parallel }^4-\frac{1}{2^{★}}{\parallel u_{+}\parallel }_{2^{★}}^{2^{★}}-\lambda {\int }_{\Omega }F(x,u_{+}\left(x\right))dx.$$

Finally, we can introduce the extended energy functional $E_{\lambda ,\mu }:H_0^1(\Omega )\mapsto [-\infty ,+\infty ]$ by

$$E_{\lambda ,\mu }\left(u\right):={\Phi }_{\lambda }\left(u\right)+\mu \Psi \left(u\right).$$

Lemma 1.

Assume $a^{\frac{N-4}{2}}b\ge C_1\left(N\right)$. Then, the functional ${\Phi }_0$ is sequentially weakly lower semicontinuous.

The proof is the same as in Lemma 2.1 of [8]. Here, we indicate the main steps for completeness.

Proof. 

Let $\left\{u_n\right\}\subset H_0^1(\Omega )$ such that $u_n\rightharpoonup u\in H_0^1(\Omega )$ (thus, $\left\{u_n\right\}$ is bounded). Let us define an auxiliary functional $\Gamma :H_0^1(\Omega )\mapsto \mathbb{R}$ by

$$\Gamma \left(u\right)=\frac{1}{2}{a\parallel u\parallel }^2+\frac{1}{4}{b\parallel u\parallel }^4-\frac{1}{2^{★}}{\parallel u\parallel }_{2^{★}}^{2^{★}}.$$

Simple computations and the Brezis–Lieb lemma give us the relations

$$\begin{array}{cc}\hfill \parallel u_n{\parallel }^2-{\parallel u\parallel }^2& =\parallel u_n{-u\parallel }^2+o_n\left(1\right),\hfill \\ \hfill \parallel u_n{\parallel }^4-{\parallel u\parallel }^4& \ge \parallel u_n{-u\parallel }^4+o_n\left(1\right),\hfill \\ \hfill \parallel u_n{\parallel }_{2^{★}}^{2^{★}}-{\parallel u\parallel }_{2^{★}}^{2^{★}}& =\parallel u_n{-u\parallel }_{2^{★}}^{2^{★}}+o_n\left(1\right).\hfill \end{array}$$

Combining the relations above, we obtain

$$\begin{array}{ccc}\hfill \Gamma \left(u_n\right)-\Gamma \left(u\right)& \ge \hfill & \parallel u_n{-u\parallel }^2\left(\frac{1}{2}a+\frac{1}{4}b\parallel u_n{-u\parallel }^2-\frac{S_0}{2^{★}}{\parallel u_n-u\parallel }^{2^{★}-2}\right)+o_n\left(1\right).\hfill \end{array}$$

Define, for $N\ge 4$, the convex function $g:[0,+\infty [\mapsto \mathbb{R}$ given by

$$g\left(t\right)\equiv \left(\frac{a}{2}+\frac{b}{4}{t}^2-\frac{S_0}{2^{★}}{t}^{2^{★}-2}\right).$$

For $N>4$, g attains its minimum at

$$t_0={\left(\frac{2(2^{★}-2)}{2^{★}}\frac{S_0}{b}\right)}^{\frac{1}{4-2^{★}}},$$

and simple computations show that

$$\begin{array}{ccc}\hfill g\left(t_0\right)\ge 0& \iff & a^{\frac{N-4}{2}}b\ge C_1\left(N\right).\hfill \end{array}$$

For $N=4$ (in this case, $2^{★}=4$), we rewrite $g\left(t\right)\equiv \left(\frac{a}{2}+\left(b-S_0\right)\frac{t^2}{4}\right)$, which attains its minimum value $\frac{a}{2}$ at $t_0=0$ for $b\ge S_0$.

Therefore, for $N\ge 4$, if $a^{\frac{N-4}{2}}b\ge C_1\left(N\right)$, it holds that

$$\underset{n\to +\infty }{lim\; inf}\left[\Gamma \left(u_n\right)-\Gamma \left(u\right)\right]\ge \underset{n\to +\infty }{lim\; inf}\left[\parallel u_n{-u\parallel }^{2}g\left(t_0\right)\right]\ge 0.$$

To finish the proof, it is enough to note that

$${\Phi }_0\left(u\right)\equiv \Gamma \left(u\right)+\frac{1}{2^{★}}{\parallel u_{-}\parallel }_{2^{★}}^{2^{★}}.$$

□

Corollary 1.

Assume $a^{\frac{N-4}{2}}b\ge C_1\left(N\right)$ and $\left(f_1\right)$. Then, for each $\lambda >0$, the functional ${\Phi }_{\lambda }$ is sequentially weakly lower semicontinuous.

Lemma 2.

Assume the following: either $N=4$ and $a^{\frac{N-4}{2}}b\ge C_1\left(N\right)$, or $N>4$ and $a,b>0$ are arbitrary. Then, the functional ${\Phi }_0$ is coercive.

Proof. 

Indeed,

$$\begin{array}{ccc}\hfill {\Phi }_0\left(u\right)& \ge & \frac{1}{2}{a\parallel u\parallel }^2+b\frac{1}{4}{\parallel u\parallel }^4-\frac{1}{2^{★}}{\parallel u\parallel }_{2^{★}}^{2^{★}}\hfill \\ & \ge & \frac{1}{2}{a\parallel u\parallel }^2+b\frac{1}{4}{\parallel u\parallel }^4-S_0\frac{1}{2^{★}}{\parallel u\parallel }^{2^{★}}.\hfill \end{array}$$

Therefore, for $2^{★}<4$ (when $N>4$), the claim is true for all $a,b>0$; if $2^{★}=4$ (when $N=4$), the conclusion follows only for $b\ge S_0$. □

Corollary 2.

Under the same conditions as in Lemma 2, the energy functional $E_{\lambda ,\mu }$ is coercive for all $\lambda ,\mu \ge 0$.

We denote by $\mathcal{W}$ the class of functionals $\Phi :X\to \mathbb{R}$ having the following property:

$\left(\mathcal{W}\right)$

if ${\left\{u_n\right\}}_{n\in \mathbb{N}}$ is a sequence in X such that $u_n\rightharpoonup u$ (weakly) and

$$\underset{n}{lim\; inf}\Phi \left(u_n\right)\le \Phi \left(u\right);$$

then, it has a sub-sequence strongly converging to u.

Proposition 1.

Assume $a^{\frac{N-4}{2}}b>C_1\left(N\right)$. Then, ${\Phi }_0\in \mathcal{W}$.

Proof. 

Let ${\left\{u_n\right\}}_{n\in \mathbb{N}}$ be a sequence in $H_0^1(\Omega )$ such that $u_n\rightharpoonup u$ and ${lim\; inf}_n{\Phi }_0\left(u_n\right)\le {\Phi }_0\left(u\right).$ By Lemma 1, ${lim\; inf}_n{\Phi }_0\left(u_n\right)={\Phi }_0\left(u\right).$ Thus, recalling that

$$\Gamma \left(u\right)={\Phi }_0\left(u\right)-\frac{1}{2^{★}}{\parallel u_{-}\parallel }_{2^{★}}^{2^{★}},$$

by the weak lower semicontinuity of the norm, we deduce that

$$\underset{n\to \infty }{lim\; inf}\Gamma \left(u_n\right)\le {\Phi }_0\left(u\right)-\frac{1}{2^{★}}{\parallel u_{-}\parallel }_{2^{★}}^{2^{★}}=\Gamma \left(u\right).$$

(3)

From the proof of Lemma 1, we observe that

$$\Gamma \left(u_n\right)-\Gamma \left(u\right)\ge {\parallel u_n-u\parallel }^{2}g\left(t_0\right)+o_n\left(1\right),$$

and thus

$$\underset{n\to +\infty }{lim\; inf}\left[\Gamma \left(u_n\right)-\Gamma \left(u\right)\right]\ge \underset{n\to +\infty }{lim\; inf}\left[\parallel u_n{-u\parallel }^{2}g\left(t_0\right)\right].$$

Now, if $a^{\frac{N-4}{2}}b>C_1\left(N\right)$, one has $g\left(t_0\right)>0$, and the claim follows by (3). □

The following result will be crucial in the subsequent work.

Theorem 2

([15]). Let X be a reflexive and separable real Banach space, and let $I,J:X\to \mathbb{R}$ be two sequentially weakly lower semicontinuous functionals with $I\in \mathcal{W}$. Assume that

$$\underset{\parallel u\parallel \to \infty }{lim}I\left(u\right)+J\left(u\right)=+\infty .$$

Then, any strict local minimizer of $I+J$ in the strong topology is the same in the weak topology.

Corollary 3.

Assume $a^{\frac{N-4}{2}}b>C_1\left(N\right)$ and $\left(f_1\right),\left(f_2\right)$. Then, 0 is a local minimizer of ${\Phi }_{\lambda }$ in the weak topology.

Proof. 

We apply Theorem 2 choosing $I={\Phi }_0$, and J is defined by $J\left(u\right)=-\lambda {\int }_{\Omega }F(x,u_{+}\left(x\right))dx$. By Lemma 1 and assumption $\left(f_1\right)$, I and J are sequentially weakly lower semicontinuous, and Proposition 1 ensures that $I\in \mathcal{W}$. The coercivity of $I+J={\Phi }_{\lambda }$ follows by Lemma 2 (for $\mu =0$). By $\left(f_2\right)$, 0 is a local minimizer of ${\Phi }_{\lambda }$ in the strong topology, and the claim follows. □

The singular term prevents the application of the classical critical point theory for $C^1$ functionals. We will need the following.

Definition 2

([16]). Let X be a real Banach space, $\Phi \in C^1\left(X\right)$, and $\Psi :X\mapsto \mathbb{R}\cup \{+\infty \}$ proper, convex and lower semicontinuous. Then, $I=\Phi +\Psi$ is referred to as the Szulkin functional. Moreover, a point $u\in X$ is said to be critical for the Szulkin functional I if $u\in \mathrm{dom}\Psi \left(u\right)=\{u\in X:\Psi (u)<+\infty \}$ and

$$\langle {\Phi }^{\prime }\left(u\right),v-u\rangle +\Psi \left(v\right)-\Psi \left(u\right)\ge 0\mathit{for}\mathit{all}v\in X.$$

It is well-known that a local minimum of I is a critical point of I.

Lemma 3

([13], Lemma 3.1). Assume $\left(s_1\right)$. Then,

$$\mathrm{int}\left(C_0^1{\left(\overline{\Omega }\right)}_{+}\right)\subset \mathrm{dom}(\Psi )$$

where $\mathrm{int}\left(C_0^1{\left(\overline{\Omega }\right)}_{+}\right)$ denotes the interior in the ordered Banach space $C_0^1\left(\overline{\Omega }\right)$ of the positive cone

$$C_0^1{\left(\overline{\Omega }\right)}_{+}=\{u\in C_0^1\left(\overline{\Omega }\right):u\left(x\right)\ge 0\forall x\in \Omega \}.$$

Remark 2.

The energy functional $E_{\lambda ,\mu }$ associated with problem $\left(\mathcal{P}\right)$ is a Szulkin functional. Indeed, ${\Phi }_{\lambda }$ is of class $C^1\left(H_0^1(\Omega )\right)$, while Ψ, defined as in (2), is a convex and lower semicontinuous (l.s.c.) functional. Moreover, it is proper from the previous result.

Definition 3.

We say that the operator $T:X\mapsto X^{★}$ satisfies the ${\left(S\right)}_{+}$ condition if the following is true: let ${\left\{u_n\right\}}_{n\ge 1}\subset X$ such that $u_n\rightharpoonup u\in X$ and

$$\underset{n\to +\infty }{lim\; sup}〈T\left(u_n\right),u_n-u〉\le 0.$$

Then, $u_n\mapsto u\in X$.

Proposition 2.

Assume $a^{\frac{N-4}{2}}b>C_2\left(N\right)$. Then, the operator ${\Phi }_0^{\prime }$ satisfies the ${\left(S\right)}_{+}$ condition.

Proof. 

Let $n\ge 1$ and $t\in [0,1]$. A straightforward calculation shows that

$$\begin{array}{cc}\hfill & {\Phi }_0^{″}(tu+(1-t)u_n)(u_n-u)(u_n-u)\ge \hfill \\ \hfill & a\parallel u_n{-u\parallel }^2+b\parallel tu+(1-t)u_n{\parallel }^2\parallel u_n{-u\parallel }^2-(2^{★}-1)S_0\parallel {\left(tu+(1-t)u_n\right)}_{+}{\parallel }^{\frac{4}{N-2}}{\parallel u_n-u\parallel }^2\ge \hfill \\ \hfill & a\parallel u_n{-u\parallel }^2+b\parallel tu+(1-t)u_n{\parallel }^2\parallel u_n{-u\parallel }^2-(2^{★}-1)S_0\parallel tu+(1-t)u_n{\parallel }^{\frac{4}{N-2}}{\parallel u_n-u\parallel }^2=\hfill \\ \hfill & \left[a+b\parallel tu+(1-t)u_n{\parallel }^2-(2^{★}-1)S_0{\parallel tu+(1-t)u_n\parallel }^{\frac{4}{N-2}}\right]{\parallel u_n-u\parallel }^2.\hfill \end{array}$$

Now, define the auxiliary function

$$g\left(t\right)\equiv a+b{t}^2-\frac{N+2}{N-2}{S}_0{t}^{\frac{4}{N-2}}.$$

If $N>4$, g attains its minimum at

$$t_0={\left[S_0\frac{2(N+2)}{b{(N-2)}^2}\right]}^{\frac{N-2}{2(N-4)}},$$

and $g\left(t_0\right)>0$ if and only if

$$a{b}^{\frac{2}{N-4}}>\left[\frac{N-4}{2}\right]{\left[\frac{2(N+2)}{{(N-2)}^2}{S}_0\right]}^{\frac{N-2}{N-4}}.$$

If $N=4$,

$$g\left(t\right)=a+(b-3S_0)t^2,$$

and if $b>3S_0$, $g\left(t\right)\ge a$. Subsequently, if $a^{\frac{N-4}{2}}b>C_2\left(N\right)$, there exists a constant $C>0$ such that

$${\Phi }_0^{″}(tu+(1-t)u_n)(u_n-u)(u_n-u)\ge C{\parallel u_n-u\parallel }^2.$$

Let $\left\{u_n\right\}$ be a sequence such that $u_n\rightharpoonup u$ and

$$\underset{n\to +\infty }{lim\; sup}〈{{\Phi }_0}^{\prime }\left(u_n\right),u_n-u〉\le 0.$$

Then,

$$\underset{n\to +\infty }{lim\; sup}〈{{\Phi }_0}^{\prime }\left(u\right),u_n-u〉=0,$$

and so

$$\begin{array}{cc}\hfill 0& \ge \underset{n\to +\infty }{lim\; sup}〈{{\Phi }_0}^{\prime }\left(u_n\right),u_n-u〉=\underset{n\to +\infty }{lim\; sup}〈{\Phi }_0^{\prime }\left(u_n\right)-{\Phi }_0^{\prime }\left(u\right),u_n-u〉\hfill \\ & =\underset{n\to +\infty }{lim\; sup}{\int }_0^1{\Phi }_0^{″}(tu+(1-t)u_n)(u_n-u)(u_n-u)dt\hfill \\ & \ge \underset{n\to +\infty }{lim\; sup}C{\parallel u_n-u\parallel }^2.\hfill \end{array}$$

Thus, we conclude that

$$\underset{n\to +\infty }{lim}{\parallel u_n-u\parallel }^2=0;$$

that is our claim. □

From Proposition 2, it follows that:

Corollary 4.

Assume $a^{\frac{N-4}{2}}b>C_2\left(N\right)$ and $\left(f_1\right)$. Then, the operator ${\Phi }_{\lambda }^{\prime }$ satisfies the ${\left(S\right)}_{+}$ condition.

Remark 3.

Proposition 2 also shows that ${\Phi }_0^{\prime }$ is strongly monotone; thus, in particular, ${\Phi }_0$ is strictly convex.

Definition 4.

Let $I=\Phi +\Psi$ be a Szulkin functional defined on a Banach space X. We say that I satisfies the Palais–Smale condition if, for any sequence ${\left\{u_n\right\}}_{n\ge 1}\subset X$ and ${\left\{{ϵ}_n\right\}}_{n\ge 1}\subset {\mathbb{R}}_{+}$ such that $I\left(u_n\right)\mapsto c\in \mathbb{R}$, ${ϵ}_n\mapsto 0$, and

$$〈{\Phi }^{\prime }\left(u_n\right),v-u_n〉+\Psi \left(v\right)-\Psi \left(u_n\right)\ge -{ϵ}_n\parallel v-u_n\parallel$$

for all $n\ge 1$ and for all $v\in X$, ${\left\{u_n\right\}}_{n\ge 1}$ possesses a strongly convergent sub-sequence.

Proposition 3.

Assume $a^{\frac{N-4}{2}}b>C_2\left(N\right)$, $\left(f_1\right)$. Then, the energy functional $E_{\lambda ,\mu }$ satisfies the Palais–Smale condition.

Proof. 

Let ${\left\{u_n\right\}}_{n\ge 1}$ be a sequence such that $\left\{E_{\lambda ,\mu }\left(u_n\right)\right\}$ converges to some $c\in \mathbb{R}$, and let ${\left\{{ϵ}_n\right\}}_{n\ge 1}$ be a sequence of positive real numbers such that ${ϵ}_n\mapsto 0$. Since $E_{\lambda ,\mu }$ is coercive by Corollary 2, $u_n\rightharpoonup u\in H_0^1(\Omega )$ up to a sub-sequence. Moreover, we note that $E_{\lambda ,\mu }$ is sequentially weakly lower semicontinuous, and thus $E_{\lambda ,\mu }\left(u\right)<+\infty$. Subsequently, $u\in \mathrm{dom}(\Psi )$. Setting $v=u$ in the inequality in Definition 4, we obtain

$$〈{\Phi }_{\lambda }^{\prime }\left(u_n\right),u_n-u〉\le \mu \Psi \left(u\right)-\mu \Psi \left(u_n\right)+{ϵ}_n\parallel u-u_n\parallel .$$

Now, since $\Psi$ is sequentially weakly lower semicontinuous, there holds

$$\begin{array}{ccc}\hfill \underset{n\to +\infty }{lim\; sup}〈{\Phi }_{\lambda }^{\prime }\left(u_n\right),u_n-u〉& \le & \underset{n\to +\infty }{lim\; sup}\left(\mu \Psi \left(u\right)-\mu \Psi \left(u_n\right)+{ϵ}_n\parallel u-u_n\parallel \right)\hfill \\ & =& \mu \underset{n\to +\infty }{lim\; sup}\left(\Psi \left(u\right)-\Psi \left(u_n\right)\right)\hfill \\ & \le & 0.\hfill \end{array}$$

From Corollary 4, ${\Phi }_{\lambda }^{\prime }$ is of type ${\left(S\right)}_{+}$, and thus $u_n\mapsto u$ strongly in $H_0^1(\Omega )$. □

Proposition 4.

Any critical point $u\in H_0^1(\Omega )$ of $E_{\lambda ,\mu }$ (in the sense of Szulkin) is a weak solution for problem $\left(\mathcal{P}\right)$.

Proof. 

By Definition 2, for all $v\in H_0^1(\Omega )$, there holds

$$〈{\Phi }_{\lambda }^{\prime }\left(u\right),v-u〉+\mu \left(\Psi \left(v\right)-\Psi \left(u\right)\right)\ge 0.$$

(4)

Since $u\in \mathrm{dom}(\Psi )$, by the definition of S, it is clear that $u>0$ almost everywhere in $\Omega$. Let $\phi \in H_0^1(\Omega )$, $\phi \ge 0$, and $v:=u+ϵ\phi$ in (4).

Then,

$$\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla \phi 〉dx-{\int }_{\Omega }{u}^{2^{★}-1}\phi dx-\lambda {\int }_{\Omega }f(x,u)\phi dx\ge \mu {\int }_{\Omega }\frac{S(x,u)-S(x,u+ϵ\phi )}{ϵ}dx.$$

Taking $ϵ\mapsto 0$ and applying Fatou’s Lemma, we obtain that

$$s\left(x\right)u^{-\gamma }\phi \in L^1(\Omega )$$

and

$$\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla \phi 〉dx-{\int }_{\Omega }{u}^{2^{★}-1}\phi dx-\lambda {\int }_{\Omega }f(x,u)\phi dx\ge \mu {\int }_{\Omega }s\left(x\right)u^{-\gamma }\phi dx.$$

(5)

Let $v:=(1-ϵ)u$, with $ϵ\in (0,1)$, in (4).

Then,

$$\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla u〉dx-{\int }_{\Omega }{u}^{2^{★}}dx-\lambda {\int }_{\Omega }f(x,u)udx-\mu {\int }_{\Omega }\frac{S(x,(1-ϵ\left)u\right)-S(x,u)}{ϵ}dx\le 0.$$

(6)

Applying the mean value theorem for functions of a real variable, there exists $\tau =\tau \left(ϵ\right)\in (0,ϵ)$ such that (6) can be rewritten as

$$\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla u〉dx-{\int }_{\Omega }{u}^{2^{★}-1}udx-\lambda {\int }_{\Omega }f(x,u)udx-\mu {(1-\tau \left(ϵ\right))}^{-\gamma }{\int }_{\Omega }s\left(x\right)u^{1-\gamma }dx\le 0.$$

Taking $ϵ\mapsto 0$, we obtain

$$\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla u〉dx-{\int }_{\Omega }{u}^{2^{★}}dx-\lambda {\int }_{\Omega }f(x,u)udx-\mu {\int }_{\Omega }s\left(x\right)u^{1-\gamma }dx\le 0.$$

Therefore, by (5), we deduce that

$$\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla u〉dx-{\int }_{\Omega }{u}^{2^{★}}dx-\lambda {\int }_{\Omega }f(x,u)udx-\mu {\int }_{\Omega }s\left(x\right)u^{1-\gamma }dx=0.$$

Now, take any $\phi \in H_0^1(\Omega )$ to test (5) with the function

$$v:={(u+ϵ\phi )}^{+}$$

and obtain

$$\begin{array}{cc}\hfill & \left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla (u+ϵ\phi )〉dx-{\int }_{\Omega }{u}^{2^{★}-1}(u+ϵ\phi )dx-\lambda {\int }_{\Omega }f(x,u)(u+ϵ\phi )dx\hfill \\ & -\mu {\int }_{\Omega }s\left(x\right)u^{-\gamma }(u+ϵ\phi )dx-\left(a+{b\parallel u\parallel }^2\right){\int }_{\{u+ϵ\phi <0\}}〈\nabla u,\nabla (u+ϵ\phi )〉dx\hfill \\ & +{\int }_{\{u+ϵ\phi <0\}}{u}^{2^{★}-1}(u+ϵ\phi )dx+\lambda {\int }_{\{u+ϵ\phi <0\}}f(x,u)(u+ϵ\phi )dx+\mu {\int }_{\{u+ϵ\phi <0\}}s\left(x\right)u^{-\gamma }(u+ϵ\phi )dx\ge 0.\hfill \end{array}$$

From the inequality above, we have

$$\begin{array}{cc}\hfill & \stackrel{=0}{\overbrace{\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }{|\nabla u|}^{2}dx-{\int }_{\Omega }{u}^{2^{★}}dx-\lambda {\int }_{\Omega }f(x,u)udx-\mu {\int }_{\Omega }s\left(x\right)u^{1-\gamma }dx}}+\hfill \\ & ϵ\left[\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla \phi 〉dx-{\int }_{\Omega }{u}^{2^{★}-1}\phi dx-\lambda {\int }_{\Omega }f(x,u)\phi dx-\mu {\int }_{\Omega }s\left(x\right)u^{-\gamma }\phi dx\right]+\hfill \\ & \stackrel{\le 0}{\overbrace{-\left(a+{b\parallel u\parallel }^2\right){\int }_{\{u+ϵ\phi <0\}}{|\nabla u|}^{2}dx}}-ϵ\left(a+{b\parallel u\parallel }^2\right){\int }_{\{u+ϵ\phi <0\}}〈\nabla u,\nabla \phi 〉dx+\hfill \\ & \stackrel{\le 0}{\overbrace{{\int }_{\{u+ϵ\phi <0\}}{u}^{2^{★}-1}(u+ϵ\phi )dx+\lambda {\int }_{\{u+ϵ\phi <0\}}f(x,u)(u+ϵ\phi )dx+\mu {\int }_{\{u+ϵ\phi <0\}}s\left(x\right)u^{-\gamma }(u+ϵ\phi )dx}}\ge 0.\hfill \end{array}$$

In other words,

$$\begin{array}{cc}\hfill & \left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla \phi 〉dx-{\int }_{\Omega }{u}^{2^{★}-1}\phi dx-\lambda {\int }_{\Omega }f(x,u)\phi dx-\mu {\int }_{\Omega }s\left(x\right)u^{-\gamma }\phi dx\hfill \\ & -\left(a+{b\parallel u\parallel }^2\right){\int }_{\{u+ϵ\phi <0\}}〈\nabla u,\nabla \phi 〉dx\ge 0.\hfill \end{array}$$

Taking the limit as $ϵ\mapsto 0$, we conclude that

$$\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla \phi 〉dx-{\int }_{\Omega }{u}^{2^{★}-1}\phi dx-\lambda {\int }_{\Omega }f(x,u)\phi dx-\mu {\int }_{\Omega }s\left(x\right)u^{-\gamma }\phi dx\ge 0,$$

and, by the arbitrariness of $\phi$, we conclude that

$$\left(a+{b\parallel u\parallel }^2\right){\int }_{\Omega }〈\nabla u,\nabla \phi 〉dx-{\int }_{\Omega }{u}^{2^{★}-1}\phi dx-\lambda {\int }_{\Omega }f(x,u)\phi dx-\mu {\int }_{\Omega }s\left(x\right)u^{-\gamma }\phi dx=0;$$

that is, u is a weak solution to $\left(\mathcal{P}\right)$. □

Lemma 4.

Under condition $\left(f_1\right)$, all critical points $u\in H_0^1(\Omega )$ of ${\Phi }_{\lambda }$ belong to $\mathrm{int}\left(C_0^1{\left(\overline{\Omega }\right)}_{+}\right)$.

Proof. 

Assume u as in the statement and note first that $u\ge 0$. Put $k=\left(a+{b\parallel u\parallel }^2\right)$. Thus, we can write

$$\left\{\begin{array}{cc}-\Delta u=a\left(x\right)(1+u(x\left)\right)\hfill & in\Omega \hfill \\ u=0\hfill & in\partial \Omega ,\hfill \end{array}\right.$$

(7)

where

$$a\left(x\right):=\frac{1}{k}·\frac{u{\left(x\right)}^{2^{★}-1}+\lambda f(x,\left(u\left(x\right)\right)}{1+u\left(x\right)}.$$

Let us prove that $a\left(x\right)\in L^{\frac{N}{2}}(\Omega )$.

In fact, by $\left(f_1\right)$, denoting by c a constant whose value may vary from line to line, there holds

$$\begin{array}{ccc}\hfill {\left(\frac{1}{k}·\frac{u{\left(x\right)}^{2^{★}-1}+\lambda f(x,\left(u\left(x\right)\right)}{1+u\left(x\right)}\right)}^{\frac{N}{2}}& \le & c{\left(\frac{u{\left(x\right)}^{2^{★}-1}+c(1+u{\left(x\right)}^{q-1})}{1+u\left(x\right)}\right)}^{\frac{N}{2}}\hfill \\ & \le & c(1+u{\left(x\right)}^{2^{★}})\in L^1(\Omega ).\hfill \end{array}$$

Therefore, by the Brezis–Kato arguments (see, for instance, Lemma B.3 of [17]), $u\in L^p(\Omega )$ for all $p>1$. This proves that the right-hand side of (7) belongs to $L^p(\Omega )$ for all $1<p<+\infty$, which in turn allows us to use the Calderón–Zigmund inequalities to prove that $u\in W_0^{2,p}(\Omega )$ for all $1<p<+\infty$. Finally, the Sobolev immersions imply that $u\in C^{1,\alpha }\left(\overline{\Omega }\right)$ for some $\alpha \in (0,1)$. From Theorem 1 in [18], we conclude that $u>0$ a.e. in $\Omega$. □

We conclude this section recalling the following topological result, which will be useful in the subsequent work. Its proof will be used in the proof of our main result and, for this reason, we provide it here.

Proposition 5

([13], Proposition 2.2). Let X be a Hausdorff topological space. Let ${\left\{K_n\right\}}_{n\ge 1}\subset X$ be a decreasing sequence of compact subsets such that ${\bigcap }_{n\ge 1}{K}_n$ is the disjoint union of non-empty compact sets ${\mathcal{S}}_i$ ($i=1,2$). Then, for some $\overline{n}\ge 1$, the set $K_{\overline{n}}$ is the disjoint union of non-empty compact sets ${\mathcal{C}}_i$ ($i=1,2$), where ${\mathcal{S}}_i\subset {\mathcal{C}}_i$ ($i=1,2$).

Proof. 

Since X is Hausdorff and ${\mathcal{S}}_i$ is compact, there are disjoint open sets ${\mathcal{O}}_i$ ($i=1,2$) so that ${\mathcal{S}}_i\subset {\mathcal{O}}_i$ ($i=1,2$). Moreover,

$$\begin{array}{ccc}\hfill K_1& \subset & \left[X\setminus \bigcap _{n=1}^{\infty }{K}_n\right]\bigcup \left[\bigcap _{n=1}^{\infty }{K}_n\right]\hfill \\ & =& \left[\bigcup _{n=1}^{\infty }\left(X\setminus K_n\right)\right]\bigcup \left[\bigcup _{i=1}^2{\mathcal{S}}_i\right]\hfill \\ & \subseteq & \left[\bigcup _{n=1}^{\infty }\left(X\setminus K_n\right)\right]\bigcup _{i=1}^2{\mathcal{O}}_i.\hfill \end{array}$$

By the compactness of $K_1$, there is a finite index set J such that

$$K_1\subset \left[\bigcup _{n\in J}\left(X\setminus K_n\right)\right]\bigcup _{i=1}^2{\mathcal{O}}_i.$$

Since the set sequence ${\{X\setminus K_n\}}_{n\ge 1}$ is increasing, for $\overline{n}=maxJ$, there holds

$$\begin{array}{c}\hfill K_{\overline{n}}\subset K_1\subset \left(X\setminus K_{\overline{n}}\right)\bigcup _{i=1}^2{\mathcal{O}}_i.\end{array}$$

In other words,

$$K_{\overline{n}}\subset \bigcup _{i=1}^2{\mathcal{O}}_i.$$

(8)

Define

$$\begin{array}{ccc}\hfill C_i:=K_{\overline{n}}\cap {\mathcal{O}}_i& for& i=1,2.\hfill \end{array}$$

(9)

Relations (8) and (9) tell us that ${\mathcal{C}}_i$ ($i=1,2$) are disjoint and compact. It is also clear that ${\mathcal{S}}_i\subset {\mathcal{C}}_i$ ($i=1,2$). □

3. Main Result

Now we are in a position to prove Theorem 1. First, we prove the existence of two local minimizers for the energy functional associated with problem $\left(\mathcal{P}\right)$. The proof follows as in [13], Theorem 1.1 (see also [19]).

Lemma 5.

Assume $a^{\frac{N-4}{2}}b>C_1\left(N\right)$. Under conditions $\left(f_1\right)$, $\left(f_2\right)$, and $\left(s_1\right)$, the energy functional $E_{\lambda ,\mu }$ has at least two local minimizers.

Proof. 

From assumption $\left(f_2\right)$, it follows that $0\in H_0^1(\Omega )$ is a strict local minimizer of ${\Phi }_{\lambda }$ in the strong topology. Consequently, 0 is also a strict local minimizer in the weak topology from Corollary 3.

In other words, there is a weak neighborhood of zero $\mathcal{O}\subset H_0^1(\Omega )$ such that ${\Phi }_{\lambda }\left(u\right)>0\forall u\in \mathcal{O}\setminus \left\{0\right\}$. Thus,

$${\Phi }_{\lambda }^{-1}\left(\left[-\infty ,0\right]\right)=\bigcap _{n\ge 1}^{\infty }{\Phi }_{\lambda }^{-1}\left(\left[-\infty ,\frac{1}{n}\right]\right)=\left\{0\right\}\cup \left({\Phi }_{\lambda }^{-1}\left(\right[-\infty ,0\left]\right)\setminus \mathcal{O}\right).$$

By the definition of ${\lambda }^{★}$, for $\lambda >{\lambda }^{★}$, ${inf}_{u\in H_0^1(\Omega )}{\Phi }_{\lambda }<0$, and thus ${\Phi }_{\lambda }^{-1}\left(\right[-\infty ,0\left]\right)\setminus \mathcal{O}\ne \varnothing$. Note that the sublevel sets of ${\Phi }_{\lambda }$ are weakly compact and metrizable (thus sequentially weakly compact) with respect to the weak topology. By Proposition 5, for some $\overline{n}\ge 1$, there exist weakly compact and sequentially weakly compact disjoint sets ${\mathcal{C}}_i$ ($i=1,2$) with $\left\{0\right\}\subset {\mathcal{C}}_1$ and ${\Phi }_{\lambda }^{-1}\left(\right[-\infty ,0\left]\right)\setminus \mathcal{O}\subset {\mathcal{C}}_2$, such that

$${\Phi }_{\lambda }^{-1}\left(\left[-\infty ,\frac{1}{\overline{n}}\right]\right)={\mathcal{C}}_1\cup {\mathcal{C}}_2.$$

Since ${\Phi }_{\lambda }+\mu \Psi$ is sequentially weakly lower semicontinuous and ${\mathcal{C}}_i$ ($i=1,2$) are sequentially weakly compact, we obtain the existence of $u_i\in {\mathcal{C}}_i$ ($i=1,2$) such that

$$\underset{{\mathcal{C}}_i}{inf}({\Phi }_{\lambda }+\mu \Psi )={\Phi }_{\lambda }\left(u_i\right)+\mu \Psi \left(u_i\right).$$

We wish to prove that, for $i=1,2$, there holds

$$u_i\in {\mathcal{O}}_i^{\prime }:=\{u\in {\mathcal{O}}_i:{\Phi }_{\lambda }\left(u\right)<\frac{1}{\overline{n}}\},$$

where ${\mathcal{O}}_i$ ($i=1,2$) come from the proof of Proposition 5. Notice also that ${\mathcal{O}}_i^{\prime }\subset {\mathcal{C}}_i$ ($i=1,2$).

The ${\mathcal{O}}_i$ values are weakly open and therefore strongly open, and ${\Phi }_{\lambda }$ is strongly continuous. Therefore, the ${\mathcal{O}}_i^{\prime }$ values are strongly open. Let us prove that

$$\mathrm{dom}(\Psi )\cap {\mathcal{O}}_i^{\prime }\ne \varnothing .$$

Indeed, by Lemma 3, it is enough to prove that

$$\mathrm{int}\left(C_0^1\left(\overline{\Omega }\right)\right)\cap {\mathcal{O}}_i^{\prime }\ne \varnothing .$$

Since $0\in {\mathcal{O}}_1^{\prime }$, then there exists a positive number $\epsilon$ such that $B(0,\epsilon )\subset {\mathcal{O}}_1^{\prime }$. Then, $\mathrm{dom}(\Psi )\cap {\mathcal{O}}_1^{\prime }\ne \varnothing$.

On the other hand, by the coercivity and the sequential weak lower semicontinuity of ${\Phi }_{\lambda }$, there exists $u_0\in {\mathcal{O}}_2^{\prime }$ as a global minimizer of ${\Phi }_{\lambda }$. Thus, $u_0$ is a critical point of ${\Phi }_{\lambda }$, and by Lemma 4, $u_0\in \mathrm{int}\left(C_0^1{\left(\overline{\Omega }\right)}_{+}\right)$; thus, $\mathrm{dom}(\Psi )\cap {\mathcal{O}}_2^{\prime }\ne \varnothing$.

Then, we may define in $\mathbb{R}$ the constant

$$m_i:=\underset{\mathrm{dom}(\Psi )\cap {\mathcal{O}}_i^{\prime }}{inf}\frac{\Psi -{min}_{{\mathcal{C}}_i}\Psi }{\frac{1}{\overline{n}}-{\Phi }_{\lambda }}i=1,2.$$

Take ${\mu }^{★}>0$ such that

$$\frac{1}{{\mu }^{★}}>max\{m_i:i=1,2\}.$$

Therefore, for all $0<\mu <{\mu }^{★}$, there are $y_i\in \mathrm{dom}(\Psi )\cap {\mathcal{O}}_i^{\prime }$ such that

$$\mu \Psi \left(y_i\right)+{\Phi }_{\lambda }\left(y_i\right)<\mu \underset{{\mathcal{C}}_i}{min}\Psi +\frac{1}{\overline{n}}.$$

Thus, we may find $u_1\in \mathrm{dom}(\Psi )\cap C_1,u_2\in \mathrm{dom}(\Psi )\cap C_2$ s.t.:

$$\underset{C_i}{min}({\Phi }_{\lambda }+\mu \Psi )={\Phi }_{\lambda }\left(u_i\right)+\mu \Psi \left(u_i\right),i=1,2.$$

We wish to prove that $u_i\in {\mathcal{O}}_i^{\prime }$, $i=1,2$. Assume by contradiction that, for some $i\in \{1,2\}$, it holds that $u_i\notin {\mathcal{O}}_i^{\prime }$. This means that ${\Phi }_{\lambda }\left(u_i\right)\ge \frac{1}{\overline{n}}$.

Then,

$$\begin{array}{ccc}\hfill \underset{{\mathcal{O}}_i^{\prime }}{inf}({\Phi }_{\lambda }+\mu \Psi )& \ge & \underset{{\mathcal{C}}_i}{inf}({\Phi }_{\lambda }+\mu \Psi )\hfill \\ & =& {\Phi }_{\lambda }\left(u_i\right)+\mu \Psi \left(u_i\right)\hfill \\ & \ge & \frac{1}{\overline{n}}+\mu \underset{{\mathcal{C}}_i}{min}\Psi \hfill \\ & >& \mu \Psi \left(y_i\right)+{\Phi }_{\lambda }\left(y_i\right),\hfill \end{array}$$

which is a contradiction. □

Proof. 

$\left(1\right)$ By Lemma 5, the energy functional $E_{\lambda ,\mu }$ has at least two local minimizers. Proposition 1.1 of [16] guarantees that local minimizers are critical points in the sense of Szulkin, and from Proposition 4, the claim follows.

$\left(2\right)$ From Proposition (3), the energy functional $E_{\lambda ,\mu }$ satisfies the Palais–Smale condition. Since $C_2\left(N\right)>C_1\left(N\right)$, $E_{\lambda ,\mu }$ has two local minimizers and given Corollary 3.3 of [16], the thesis is proved. □