\documentclass[reqno]{amsart}
\usepackage{hyperref}

\AtBeginDocument{{\noindent\small
\emph{Electronic Journal of Differential Equations},
Vol. 2018 (2018), No. 174, pp. 1--21.\newline
ISSN: 1072-6691. URL: http://ejde.math.txstate.edu or http://ejde.math.unt.edu}
\thanks{\copyright 2018 Texas State University.}
\vspace{8mm}}

\begin{document}
\title[\hfilneg EJDE-2018/174\hfil 
 Solutions for magnetic Schr\"odinger equations]
{Nontrivial complex solutions for magnetic Schr\"odinger equations
with critical  nonlinearities}

\author[S. Barile, G. M. Figueiredo \hfil EJDE-2018/174\hfilneg]
{Sara Barile, Giovany M. Figueiredo}


\address{Sara Barile (corresponding author) \newline
Dipartimento di Matematica,
Universit\`a degli Studi di Bari Aldo Moro,
Via E. Orabona 4, 70125 Bari, Italy}
\email{sara.barile@uniba.it}

\address{Giovany M. Figueiredo \newline
Universidade de Brasilia - UNB,
Departamento de Matem\'atica,
Campus Universit\'ario Darcy Ribeiro,
Brasilia - DF, CEP 70.910-900, Brazil}
\email{giovany@unb.br}

\dedicatory{Communicated by Raffaella Servadei}

\thanks{Submitted September 14, 2017. Published October 22, 2018.}
\subjclass[2010]{35B33, 35J20, 35Q55}
\keywords{Magnetic Schr\"odinger equations; critical nonlinearities;
\hfill\break\indent  minimization problem; concentration-compactness methods;
 Pohozaev manifold}

\begin{abstract}
 Using minimization arguments we establish the existence of a complex  
 solution to the magnetic  Schr\"odinger equation
 $$
 - (\nabla + i A(x) )^2  u + u = f(|u|^2) u \quad \text{in }\mathbb{R}^N,
 $$
 where $N \geq 3$, $A{:}\mathbb{R}^N \to \mathbb{R}^N$ is the magnetic
 potential and $f$ satisfies some critical growth assumptions.
 First we obtain bounds from a real Pohozaev manifold. Then relate them
 to Sobolev imbedding constants and to the least energy  level  associated
 with the real equation in absence of the magnetic field (i.e., with $A(x)=0$).
 We also apply the Lions Concentration Compactness Principle to the modula
 of the  minimizing sequences involved.
\end{abstract}

\maketitle
\numberwithin{equation}{section}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{remark}[theorem]{Remark}
\allowdisplaybreaks

\section{Introduction}

The aim of this article is to study  the  magnetic Schr\"odinger equation
\begin{equation}\label{LE1}
- (\nabla + i A(x) )^2  u + u =f(|u|^2) u  \quad \text{in } \mathbb{R}^N,
 \end{equation}
where  $u{:} \mathbb{R}^N \to \mathbb{C}$, $N \geq 3$,   $i$ is the imaginary unit,
$A=(A_1,\dots,A_N){:} \mathbb{R}^N \to \mathbb{R}^N$ is the magnetic (or vector)
 potential and the nonlinear term $f{:} \mathbb{R}^+ \to \mathbb{R}$ is a regular
function satisfying suitable  assumptions and  having critical growth
at infinity  with critical Sobolev exponent $2^* = \frac{2N}{N-2}$
for $N \geq 3$.   In the recent literature,  magnetic
Schr\"odinger equations  have been studied in the critical case from
different points of view but in  few papers which we recall here in the following.

 Esteban and  Lions \cite{EstLions} (1989), found  solutions to
$$
(-i \nabla +  A(x) )^2  u + \lambda u = |u|^4 u \quad \text{in }\mathbb{R}^3
 $$
with $\lambda \in \mathbb{R}$  by solving constrained minimization problems with
Concentration-Compactness methods.  Arioli and Szulkin \cite{ariolisz} (2003),
found  non-trivial solutions to
$$
(-i \nabla +  A(x) )^2  u +  V(x) u  =|u|^{2^* - 2} u  \quad
\text{in }\mathbb{R}^N,
$$
with $A$ and $V$ locally Lebesgue measurable  by means of constrained
minimization and Concentration-Compactness arguments under suitable
assumptions on the spectrum of the operator $(-i \nabla + A(x))^2 + V$.
Chabrowski-Szulkin \cite{chsz} (2005)  proved  the existence of a non
trivial solution to
$$
(-i \nabla +   A(x) )^2  u +  V(x) u  =Q(x) |u|^{2^* - 2} u  \quad
\text{in } \mathbb{R}^N,
$$
when the electric potential $V$ changes sign  by a min-max type argument
based on a topological linking.  Certain regularity properties of solutions
for a rather general class of equations involving the operator
$(-i \nabla +  A(x) )^2$ are also established.
 Barile, Cingolani, Secchi \cite{Barcs} (2006) established
 existence results  by abstract perturbation techniques to
$$
(i \nabla + \varepsilon  A(x) )^2  u + \varepsilon^{\alpha} V(x) u
=|u|^{2^* - 2} u  \quad \text{in } \mathbb{R}^N,
$$
where  $\varepsilon \in (0, \varepsilon_0)$,  $\alpha \in [1, 2]$, $N > 4$
and the potentials $A$ and $V$ are bounded continuous and Lebesgue measurable.
Han \cite{pigong} (2006) showed the existence of a non-trivial complex solution
to
$$  (-i \nabla +  A(x) )^2  u  - V(x)  u  =|u|^{2^* - 2} u  \quad
\text{in }\mathbb{R}^N,
$$
with $N \geq 3$, have been established by a Mountain Pass Theorem  under
suitable assumptions on the integrability and the  behaviour
of the magnetic potential $A$ and the electric potential $V$.
Wang \cite{wangF} (2008) established  existence results for a nontrivial
solution to
\begin{gather*}
 (-i \nabla +  A(x) )^2  u  + \lambda  V(x)  u  =K(x)  |u|^{2^* - 2} u  \quad
\text{in } \mathbb{R}^N,  \\
 \lim_{|x| \to + \infty } u(x) =0
\end{gather*}
with $N \geq 3$, by means of Linking Theorem applied twice when $\lambda  > 0$,
the magnetic potential  $A \in L^2_{\rm loc}(\mathbb{R}^N)$, the electric
potential $V(x)$ is sign-changing, $K(x)$ is positive bounded and continuous
and $V, K$ satisfying suitable local assumptions.
 Liang and Zhang \cite{liangzhang} (2011) studied standing waves
solutions $\psi(x, t) = e^{-\frac{i E t}{\hbar}} u(x)$,
$(t, x) \in \mathbb{R} \times \mathbb{R}^N$, $N \geq 3$,  to
$$
i \hbar \frac{\partial \psi}{\partial t}
= - \frac{\hbar^2}{2m} ( \nabla + i  A(x) )^2 \psi + W(x) \psi
 - K(x)  |\psi|^{2^* - 2} \psi - h(x, |\psi|^2) \psi,
$$
thus establishing  the existence of at least one solution and,
for any $m \in \mathbb{N}$, the existence of at least $m$ pairs of solutions
under suitable  assumptions.
 Ding and Liu \cite{dingliu} (2013) proved then  existence and have also described
 concentration phenomena  of  (ground states) solutions to
$$
(-i \varepsilon \nabla +  A(x) )^2  u  + V(x) u
 = W(x) (g(|u|) + |u|^{2^* - 2}) u  \quad \text{in } \mathbb{R}^N
$$
in the semiclassical limit  (i.e.\ as $\varepsilon \to 0$) when
$A \in C^1(\mathbb{R}^N, \mathbb{R}^N) \cap L^{\infty}(\mathbb{R}^N, \mathbb{R}^N)$,
$V, W$ are positive and satisfy proper boundedness assumptions  and
$g(|u|) u$ is superlinear and subcritical.
 Liang and Song \cite{liangsongBVP} (2014) treated
$$
- \varepsilon^2  (\nabla + i  A(x) )^2  u + V(x) u
= |u|^{2^* - 2} u + h(x, |u|^2) u \quad \text{in } \mathbb{R}^N
$$
where $N \geq 3$ and $V(x)$ is a nonnegative potential  by establishing
for $\varepsilon > 0$ sufficiently small  the existence of both one solution
and $m$ pairs of solutions for every $m \in \mathbb{N}$ by means of Lions' second
Concentration-Compactness method and Concentration Compactness principle
at infinity in order to recover a $(PS)_c$ condition.
   Alves and Figueiredo \cite{alvesfig} (2014) studied the multiplicity
of nontrivial solutions to
\[
\begin{gathered}
(- i \nabla -  A(x) )^2  u = \mu |u|^{q-2} u + |u|^{2^* - 2} u  \quad
  \text{in } \Omega, \\
u=0  \quad  \text{in } \partial \Omega,
\end{gathered}
\]
where $\Omega$ is a smooth bounded domain  of $\mathbb{R}^N$ with $N \geq 4$,
 $A$ is continuous on $\overline \Omega$, $2 \leq q < 2^*$  thus relating
the number of solutions with the topology of $\Omega$ by
Ljusternik-Schnirelmann theory.
Tang and Wang \cite{tangwang} (2015)  studied
$$
- ( \nabla + i   A(x) )^2  u + (\lambda a(x) - \delta) u
= |u|^{2^* - 1} u \quad \text{in } \mathbb{R}^N
$$
with $\delta > 0$ and the electric potential can be negative in some  domain;
specifically, by variational and Nehari methods they have established
the existence of a least energy solution $u_{\lambda}$ which localizes
at the bottom of the potential well as $\lambda \to + \infty$.

 Throughout this article, we use the following assumptions:
\begin{itemize}
\item[(A1)]  $A \in L^2_{\rm loc}(\mathbb{R}^N, \mathbb{R}^N)$ and
there exists $x_0 \in \mathbb{R}^N$ such that  $A$ is continuous at $x_0$;

\item[(A2)]  $ A(\frac{x}{\sigma}) = \sigma   A(x)$  for every
$x \in \mathbb{R}^N$ and for every $\sigma > 0$;

\item[(A3)]  $A \in L^2_{rad}(\mathbb{R}^N, \mathbb{R}^N)$;

\item[(A4)] $f \in C^1(\mathbb{R}^+, \mathbb{R})$  and $\lim_{s \to 0^+} f(s) =0$;

\item[(A5)] $0< \limsup_{s \to + \infty} f(s)/s^{(2^* -2)/2} \leq 1$;

\item[(A6)] $ f(s) s - F(s) \geq 0$ for every   $s \in \mathbb{R}$ with $s \geq 0$
where $F(s)= \int_0^s f( t) \, dt$;

\item[(A7)] there exist $\lambda > 0$ and $q \in (2, 2^*)$ such that
$$
f(s) \geq \lambda s^{(q-2)/2} \quad \text{for every  $s \in \mathbb{R}$ with $s \geq 0$}.
$$
\end{itemize}
A typical example of a function satisfying conditions (A1)--(A3) is
$A(x)=A/|x|$,  where $A$ is a constant vector.

 Note that by hypothesis (A2), the magnetic potential $A$ is homogeneous
of degree $-1$. This hypothesis is just used in Section 4 when we show
that the limit of minimizing sequence is a nontrivial complex solution
of Problem \ref{LE1}.

We shall prove the following  result in the case $N \geq 3$.
For the definition of the Sobolev embeddings constants $S$ and $c_A$
which appear in the next theorem see Section \ref{tools}.


\begin{theorem}\label{groundN3}
Suppose  that  $N \geq 3$  and {\rm (A1)--(A7)} hold.
Then, there exists $\lambda^* > 0$ such that for all $\lambda > \lambda^*$
problem  \eqref{LE1} admits a nontrivial complex solution.
Precisely, the constant is
$$
\lambda^* =  \Big[ 2^{(2-N)/2} S^{-N/2} N (\frac{2N}{N-2})^{(N-2)/2}\Big]^{(q-2)/2}
c_{A}^{(q-2)/2}.
$$
\end{theorem}

From our point of view, the main contributions of this article are as follows:

\noindent(1)  Inspired by  recent results obtained by Alves, Souto and Montenegro
\cite{AlvesSM} (see also Zhang and Zou \cite{ZZou}) for equation \eqref{LE1}
 when $A(x) = 0$, we aim to establish the existence of a complex solution
to the magnetic equation  \eqref{LE1} by means of Concentration Compactness
Principle of Lions \cite{Lionsibero, Lionsibero2} in the case $N \geq 3$.
Concerning  the case $N =2$, we mention the  paper by Barile and
 Figueiredo \cite{N2magnetic}.

\noindent(2) Since we do not know the Pohozaev identity associated with
the problem (P), we use Pohozaev's identity of real problem, causing
a modification in the arguments that can be found in \cite{AlvesSM}
 (see also Zhang and Zou \cite{ZZou}).

\begin{remark} \rm
Condition (A5) can be replaced by
 $$
0< \limsup_{s \to + \infty} \frac{f(s)}{s^{(2^* -2)/2}} \leq \mu
$$
with $\mu > 0$ or more in general by
 $$
0< \limsup_{s \to + \infty} \frac{f(s)}{s^{(2^* -2)/2}} < + \infty.
$$
For simplicity but  without loosing in generality, we assume (A5)
thus studying the case $\mu =1$.
\end{remark}

\begin{remark}\label{rmkfF} \rm
From (A4) and (A5), for any $\varepsilon > 0$ there exists $C_{\varepsilon} > 0$
such that
$$
|f(s) | \leq \varepsilon + C_{\varepsilon} |s|^{(2^* -2)/2}
\quad \text{for all } s \geq 0
$$
and by integration
$$
|F(s) | \leq \varepsilon |s| + \overline{C}_{\varepsilon} |s|^{2^* /2} \quad
\text{for all } s \geq 0,
$$
where $\overline{C}_{\varepsilon} = (2 C_{\varepsilon})/2^*$.
\end{remark}

\begin{remark}\label{rmkfFq} \rm
From  hypothesis  (A7), by integration we obtain
$$
F(s) \geq \frac{2}{q} \lambda  |s|^{q/2} \quad \text{for every  $s \in \mathbb{R}$ with
$s \geq 0$}.
$$
\end{remark}

This article is organized as follows.
In Section \ref{tools}, we fix notation and variational tools.
In Section \ref{prelimN3} we establish some preliminary results
which will be useful in  Section \ref{sectmainres} for proving
 Theorem \ref{groundN3}.

\section{Notation and variational tools}\label{tools}

 To introduce the variational structure of the problem, we set
$$
H^1_{A}(\mathbb{R}^N, \mathbb{C})
= \big\{u \in L^2(\mathbb{R}^N, \mathbb{C}):  \int_{\mathbb{R}^N} |\nabla_A u|^2 dx
 < + \infty   \big\}
$$
with $\nabla_A u = (\nabla + i A(x) ) u$. The space
$ H^1_{A}(\mathbb{R}^N, \mathbb{C})$ is an Hilbert space endowed with the scalar product
$$
(u, v)_A = \operatorname{Re} \int_{\mathbb{R}^N} \left(\nabla_A u
\cdot \overline{\nabla_A v } + u \overline v \right) \, dx \quad
\text{for any } u, v \in H^1_{A}(\mathbb{R}^N, \mathbb{C})
$$
where $\operatorname{Re}$  and the bar denote the real part of a complex number
and  the complex conjugation respectively.
The norm induced by this inner product is
$$
\| u \|_{A}  =  \Big(\int_{\mathbb{R}^N} ( |\nabla_A u|^2 + |u|^2 ) \, dx
\Big)^{1/2} \quad \text{for  }u \in H^1_{A}(\mathbb{R}^N, \mathbb{C})
$$
and $C_0^{\infty}(\mathbb{R}^N, \mathbb{C})$ is dense in $H^1_{A}(\mathbb{R}^N, \mathbb{C})$
with respect to the norm $\|  \cdot \|_{A} $
(see \cite[Section 2]{EstLions} and \cite[Theorem 7.22]{LL}).
 We denote by $H^{-1}_{A}(\mathbb{R}^N, \mathbb{C})$ the dual space of
$H^1_{A}(\mathbb{R}^N, \mathbb{C})$. Recall that for every
$u \in H^1_{A}(\mathbb{R}^N, \mathbb{C})$ one has
$$
\int_{\mathbb{R}^N}  |\nabla_A u|^2 dx
= \int_{\mathbb{R}^N}  |\nabla u|^2 dx
+ \int_{\mathbb{R}^N}  |A(x)|^2 |u|^2 dx
- 2 \operatorname{Re} \int_{\mathbb{R}^N} \nabla u \cdot i A(x) \overline u \, dx.
$$
Since there is no relation between   $H^1_{A}(\mathbb{R}^N, \mathbb{C})$  and
$H^1(\mathbb{R}^N, \mathbb{C})$; that is,
   $H^1_{A}(\mathbb{R}^N, \mathbb{C})  \not \subset H^1(\mathbb{R}^N, \mathbb{R})$ and
$H^1(\mathbb{R}^N, \mathbb{C}) \not \subset    H^1_{A}(\mathbb{R}^N, \mathbb{C})$,
we will frequently use in this paper the following diamagnetic inequality
(see \cite[Theorem 7.21]{LL})
\begin{equation}\label{diam}
|\nabla |u|(x)| \leq |\nabla_A u(x)| \quad \text{for almost every  }
x \in \mathbb{R}^N\,.
\end{equation}
This implies that, if  $u \in H^1_{A}(\mathbb{R}^N, \mathbb{C})$ then
 $|u| \in H^1(\mathbb{R}^N, \mathbb{R})$. Therefore,
$u \in L^p(\mathbb{R}^N, \mathbb{C})$ for any $p \in [2, 2^*]$.

 By adapting standard variational arguments exploited in existing literature
and by exploiting radial assumptions it is not difficult to prove that there
exists $\phi$ a solution to
\begin{equation}\label{Subcritical}
- (\nabla + i A(x) )^2 \phi + \phi   = |\phi|^{q-2}\phi \quad
 \text{in } \mathbb{R}^{N},\;
\phi   \in H^1_{A}(\mathbb{R}^N, \mathbb{C}).
\end{equation}
Note that if $I_{q}$ is the functional associated with problem \ref{Subcritical},
then $ I_{q}(\phi)=c_{A}$, where
\begin{gather}
 c_A = \inf_{\gamma \in \Gamma_A} \max_{t \in [0, 1]} I_q(\gamma(t)) > 0,
\nonumber \\
\Gamma_A = \big\{\gamma \in C\left([0, 1], H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C}) \right):
\gamma(0) =0  \text{ and }  I_q(\gamma(1)) < 0  \big\},\nonumber\\
\label{Sara1}
\|\phi\|^{2}_{A}= \int _{\mathbb{R}^{N}}|\phi|^{q} dx,\\
\label{Sara2}
\|\phi\|^{2}_{A} =\frac{2q}{q-2} c_A.
\end{gather}

 Moreover, we consider the space
$$
\mathcal{D}^{1, 2}_A(\mathbb{R}^N, \mathbb{C})
= \{ u \in L^{2^*}(\mathbb{R}^N, \mathbb{C}):
    \int_{\mathbb{R}^N} |\nabla_A u|^2 dx < + \infty   \},
$$
which is the closure of $C_0^{\infty}(\mathbb{R}^N, \mathbb{C})$ with respect to the norm
$$
\| u \|_{\mathcal{D}^{1, 2}_A(\mathbb{R}^N, \mathbb{C})}
 =  \Big(\int_{\mathbb{R}^N}  |\nabla_A u|^2 dx   \Big)^{1/2}  \quad
\text{for }u \in \mathcal{D}^{1, 2}_A(\mathbb{R}^N, \mathbb{C})
$$
corresponding to the inner product
$$
(u, v)_{\mathcal{D}^{1, 2}_A(\mathbb{R}^N, \mathbb{C})}
= \operatorname{Re} \int_{\mathbb{R}^N} \nabla_A u  \cdot
\overline{\nabla_A v}   \, dx  \quad \text{for  }
u, v \in \mathcal{D}^{1, 2}_A(\mathbb{R}^N, \mathbb{C}).
$$
Recall that   $\mathcal{D}^{1, 2}_A(\mathbb{R}^N, \mathbb{C}) \hookrightarrow
 L^{2^*}(\mathbb{R}^N, \mathbb{C})$. It is  also useful to  define
$$
\mathcal{D}^{1, 2}(\mathbb{R}^N, \mathbb{R}) = \big\{ u \in L^{2^*}(\mathbb{R}^N, \mathbb{R}):
   \int_{\mathbb{R}^N} |\nabla u|^2 dx < + \infty   \big\}
$$
which is the closure of $C_0^{\infty}(\mathbb{R}^N, \mathbb{R})$ with respect to the norm
$$
\| u \|_{\mathcal{D}^{1, 2}(\mathbb{R}^N, \mathbb{R})}
=  \Big(\int_{\mathbb{R}^N}  |\nabla u|^2 dx   \Big)^{1/2} \quad
\text{for  } u \in \mathcal{D}^{1, 2}(\mathbb{R}^N, \mathbb{R})
 $$
corresponding to the inner product
$$
(u, v)_{\mathcal{D}^{1, 2}(\mathbb{R}^N, \mathbb{R})}
 = \int_{\mathbb{R}^N} \nabla u  \cdot \nabla v   \, dx \quad
\text{for  }u, v \in \mathcal{D}^{1, 2}(\mathbb{R}^N, \mathbb{R}).
$$
Recall that  $\mathcal{D}^{1, 2}(\mathbb{R}^N, \mathbb{R}) \hookrightarrow
L^{2^*}(\mathbb{R}^N, \mathbb{R})$ and we denote by $S_0 >0 $ the best constant
of Sobolev embedding
$\mathcal{D}^{1, 2}(\mathbb{R}^N, \mathbb{R}) \hookrightarrow L^{2^*}(\mathbb{R}^N, \mathbb{R})$;
that is,
$$
S_0  \Big( \int_{\mathbb{R}^N} |u|^{2^*} \, dx   \Big)^{2/2^*}
\leq  \int_{\mathbb{R}^N} |\nabla u|^{2} \, dx  \quad
\text{for all } u \in \mathcal{D}^{1, 2}(\mathbb{R}^N, \mathbb{R}).
$$
If $S$ denotes the best constant of the imbedding
$D^{1,2}_{A}(\mathbb{R}^N,\mathbb{C}) \to L^{2^{*}}(\mathbb{R}^N, \mathbb{C})$,
that is,
$$
S=\inf_{u \in D^{1,2}_{A}(\mathbb{R}^N,\mathbb{C})}\frac{
\int_{\mathbb{R}^N}|\nabla_{A} u|^{2} dx }{( \int_{\mathbb{R}^N}|
u|^{2^{*}} dx )^{2/2^{*}}},
$$
we have that $S=S_0$, for details see \cite[Theorem 1.1]{ariolisz}.

 The energy functional $I_A{:} H^1_{A}(\mathbb{R}^N, \mathbb{C}) \to \mathbb{R}$
associated with \eqref{LE1} is defined as
$$
I_A(u) = \frac 12 \| u \|^2_{A} - \frac 12 \int_{\mathbb{R}^N} F(|u|^2) \, dx
\quad \text{for  } u \in H^1_{A}(\mathbb{R}^N, \mathbb{C}).
$$
Under assumptions  (A4), (A5) (see in particular Remark \ref{rmkfF}
and a direct application of \cite[Theorem 1.22]{Schectzou}),
 we obtain $I_A \in C^1(H^1_{A}(\mathbb{R}^N, \mathbb{C}), \mathbb{R})$  with
G\^{a}teaux differential
\[
I'_A(u)v = \operatorname{Re} \int_{\mathbb{R}^N}
(\nabla_A u  \cdot \overline{\nabla_A v } + u \overline v ) \, dx
 - \operatorname{Re} \int_{\mathbb{R}^N} f(|u|^2) u \overline v \, dx
\]
for all $u, v \in H^1_{A}(\mathbb{R}^N, \mathbb{C})$,
and its critical points are  the weak solutions to \eqref{LE1}.

We denote by $H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$ and $H^1_{rad}(\mathbb{R}^N, \mathbb{R})$
the subspaces of $H^1_{A}(\mathbb{R}^N, \mathbb{C})$ and $H^1(\mathbb{R}^N, \mathbb{R})$
formed by the radial functions, that is
\begin{gather*}
H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C}) = \{u \in H^1_{A}(\mathbb{R}^N, \mathbb{C}):
 u(x)= u(|x|)
\text{ for  }x \in \mathbb{R}^N\}, \\
H^1_{ rad}(\mathbb{R}^N, \mathbb{R}) = \{u \in H^1(\mathbb{R}^N, \mathbb{R}):
 u(x)= u(|x|) \text{ for  } x \in \mathbb{R}^N  \}.
\end{gather*}
Now, for finding  a nontrivial complex solution to \eqref{LE1} let
\begin{equation}\label{DefD}
D_A =  \inf \big\{\frac 12 \int_{\mathbb{R}^N}  |\nabla_A u|^2 dx
: u \in \mathcal{M}_A \big\}
\end{equation}
where
$$
\mathcal{M}_A= \big\{ u \in H^1_{A}(\mathbb{R}^N, \mathbb{C}) \setminus \{0 \}
: \int_{\mathbb{R}^N} G(u) \, dx =1 \big\}  \quad \text{with } N \geq 3
$$
with  $g(u) = f(|u|^2) u -u$ and
$G(u) = \frac 12 \left(F(|u|^2)  - |u|^2 \right)$ for
 $u \in H^1_{A}(\mathbb{R}^N, \mathbb{C})$. Hereafter, we can denote by
$$
D_A = \inf_{u \in \mathcal{M}_A} T_A(u)
$$
where for simplicity of notation
$$
T_A(u) = \frac 12 \int_{\mathbb{R}^N}  |\nabla_A u |^2 dx \quad
\text{for }u \in H^1_{A}(\mathbb{R}^N, \mathbb{C}).
$$
Moreover,   we set
\begin{gather*}
J(u) = \int_{\mathbb{R}^N} G(u) \, dx
= \frac 12 \int_{\mathbb{R}^N}  \left(F(|u|^2) - |u|^2 \right) \, dx \quad
\text{for   }u \in H^1_{A}(\mathbb{R}^N, \mathbb{C}),\\
 J'(u)v = \operatorname{Re} \int_{\mathbb{R}^N} g(u) \overline v  \, dx
= \operatorname{Re} \int_{\mathbb{R}^N}  \big(f(|u|^2)u - u \big)
\overline v  \, dx \quad \text{for  }u, v \in H^1_{A}(\mathbb{R}^N, \mathbb{C}).
\end{gather*}

\begin{lemma}\label{mountpassIA}
Suppose  {\rm (A1), (A4), (A5), (A7)} hold.
Then,  the functional $I_A$ has a mountain pass geometry, that is
\begin{itemize}
\item[$(i)$] $I_A(0) =0$,
\item[$(ii)$]  there exist $\rho_0, \delta_0  > 0$ such that
  $I_A(u) \geq \delta_0$ for all $u \in H^1_{A}(\mathbb{R}^N, \mathbb{C})$ with
 $\|u \|_A = \rho_0$;
\item[$(iii)$] there exists $u_0 \in H^1_{A}(\mathbb{R}^N, \mathbb{C})$ such that
 $\|u_0 \|_A > \rho_0$ and $I_A(u_0) \leq 0$.
\end{itemize}
\end{lemma}

\begin{proof}
(i) follows easily by Remark \ref{rmkfF}.
(ii) By exploiting  Remark \ref{rmkfF} again and  Sobolev embeddings,
for any $\varepsilon > 0$ there exists $\overline{C}_{\varepsilon} > 0$
such that
\begin{align*}
I_A(u) &=  \frac 12 \| u \|^2_{A} - \frac 12 \int_{\mathbb{R}^N} F(|u|^2) \, dx \\
&\geq  \frac 12 \| u \|^2_{A}  - \frac 12  \varepsilon  \int_{\mathbb{R}^N} |u|^2 dx - \frac 12 \, \overline{C}_{\varepsilon} \int_{\mathbb{R}^N} |u|^{2^*} \, dx \\
&\geq  \frac 12 (1- \varepsilon) \| u \|^2_{A} -  \overline{C}'_{\varepsilon}
 \| u \|^{2^*}_{A}\,.
\end{align*}
If we take $\| u \|_{A}= \rho$ and choose $\varepsilon > 0$ and $\rho > 0$
small enough, we obtain $(ii)$ is satisfied for  suitable $\rho_0, \delta_0 > 0$.

(iii) Let $v_0 \in H^1_{A}(\mathbb{R}^N, \mathbb{C}) \setminus \{0 \}$ arbitrary  and
$t > 0$. By Remark \ref{rmkfFq}  it is
$$
I_A(t v_0) \leq \frac 12 \, t^2 \| v_0 \|^2_{A} -  \frac 1q  \lambda  t^q
\int_{\mathbb{R}^N} |u|^{q} \, dx.
$$
Since $q > 2$, we obtain  $I_A(t \, v_0) \to - \infty$ as $t \to + \infty$ thus,
taken $u_0 = t v_0$ for $t$ sufficiently large (iii) is proved.
\end{proof}

 Therefore, we can  define the following minimax value or mountain pass level
of $I_A$, i.e.
$$
b_A = \inf_{\gamma \in \Gamma_A} \max_{t \in [0, 1]} I_A(\gamma(t)) > 0,
$$
where
$$
\Gamma_A = \{\gamma \in C\left([0, 1], H^1_{A}(\mathbb{R}^N, \mathbb{C}) \right):
\gamma(0) =0 \ \text{and} \ I_A(\gamma(1)) < 0  \}.
$$


 At this point, it is useful to consider the real scalar problem
\begin{equation}\label{LE2}
\begin{gathered}
- \Delta  u + u  =f(|u|^2) u \quad \text{in } \mathbb{R}^N \\
u \in H^1(\mathbb{R}^N, \mathbb{R})
\end{gathered}
\end{equation}
with $f {:} \mathbb{R}^+ \to \mathbb{R}$ satisfying  assumptions
 (A4), (A5)-(A7) when $N \geq 3$. Here below we  give a brief description
 of the results which we will exploit in next sections and which have been
established  in Alves, Souto and Montenegro \cite{AlvesSM}
by using  the ideas in Berestycki and Lions \cite{BL}, Coleman,
Glazer and Martin \cite{Colglazmart} and in  Jeanjean and Tanaka \cite{JJTan}.
The functional  $I_{0} \in C^1(H^1(\mathbb{R}^N, \mathbb{R}), \mathbb{R})    $
 associated with \eqref{LE2}    is
$$
I_{0}(u) = \frac 12 \| u \|^2 - \frac 12 \int_{\mathbb{R}^N} F(|u|^2) \, dx \quad
\text{for }u \in H^1(\mathbb{R}^N, \mathbb{R}),
$$
where
$$
\| u \|  =  \Big(\int_{\mathbb{R}^N} ( |\nabla u|^2 + |u|^2 ) \, dx   \Big)^{1/2}
\quad \text{for  } u \in H^1(\mathbb{R}^N, \mathbb{R}).
$$
The authors  investigated the existence of a ground state solution
to \eqref{LE2},  which  means a solution to
$u \in H^1(\mathbb{R}^N, \mathbb{R})$ such that $I_{0}(u) \leq I_{0}(v)$
for every nontrivial solution $v \in H^1(\mathbb{R}^N, \mathbb{R})$ of \eqref{LE2}.
Denoting
$$
m_{0} = \inf \{I_{0}(u): u \text{is a nontrivial solution to \eqref{LE2}} \}
 $$
and taking into  consideration  the set of non-zero critical point of
$I_{0}$, namely
$$
Sigma_{0} = \{  u \in H^1(\mathbb{R}^N, \mathbb{R}) \setminus \{0 \}: I'_{0}(u) =0 \}, $$
it follows that
$$
m_{0} = \inf_{u \in \Sigma_{0}} I_{0}(u).
$$
Let
$$
D_{0} = \inf_{u \in \mathcal{M}_{0}} T_{0}(u)  \quad\text{with}\quad 
  T_{0}(u) = \frac 12 \int_{\mathbb{R}^N}  |\nabla u |^2 dx \quad
 \text{for  }u \in H^1(\mathbb{R}^N, \mathbb{R})
$$
and the   $C^1$ manifold
$$
\mathcal{M}_{0}= \big\{ u \in H^1(\mathbb{R}^N, \mathbb{R}) \setminus \{0 \}:
\int_{\mathbb{R}^N} G(u) \, dx =1 \big\} \quad\text{with } N \geq 3.
$$
It has been taken into account  the Pohozaev identity manifold
\begin{align*}
\mathcal{P}_{0}
&= \big\{u \in H^1(\mathbb{R}^N, \mathbb{R}) \setminus \{ 0 \}:
 \frac{N-2}{2} \int_{\mathbb{R}^N}  |\nabla u|^2 dx
= N  \int_{\mathbb{R}^N} G(u)  \, dx \big\} \\
&= \big\{u \in H^1(\mathbb{R}^N, \mathbb{R}) \setminus \{ 0 \}:
  (N-2) \int_{\mathbb{R}^N}  |\nabla u|^2 dx
=  N \int_{\mathbb{R}^N} \left(F(|u|^2)- |u|^2 \right) \, dx \big\}.
\end{align*}
and
$p_{0} = \inf_{u \in \mathcal{P}_{0}} I_{0}(u)$.
Since $I_{0}$ has a mountain pass geometry,   we  define
$$
_{0} = \inf_{\overline{\gamma} \in \Gamma}
\max_{t \in [0, 1]} I_{0}(\overline{\gamma}(t)),
$$
where
$$
\Gamma_{0}= \{\overline{\gamma} \in C\left([0, 1], H^1(\mathbb{R}^N, \mathbb{R}) \right):
 \overline{\gamma}(0) =0 \text{ and }  I_{0}(\overline{\gamma}(1)) < 0  \}.
$$
In  \cite{AlvesSM} the authors  showed that   $D_{0}$ is attained in
$H^1_{rad}(\mathbb{R}^N, \mathbb{R})$ and the following least energy characterizations holds
$$
m_{0} = b_{0} = p_{0} =  \frac 1N \big(\frac{N-2}{2N} \big)^{(N-2)/2} (2D_{0})^{N/2}
  \quad\text{if } N \geq 3\,;
 $$
so that \eqref{LE2} has a nontrivial  ground state solution.

These existence results have been established without assuming two
widely used conditions, that is the monotonicity  condition
$$
\frac{f(s^2)s}{s} \quad \text{is increasing in } (0, + \infty)
 $$
and the Ambrosetti-Rabinowitz condition:
there exists a constant $\theta > 2$ such that
$$
0 < \theta F(s^2) \leq f(s^2)s^2\quad \text{for any }s \in \mathbb{R}\setminus \{0 \}.
$$
Therefore, the paper   \cite{AlvesSM} complements the results obtained
in the subcritical case by Jeanjean and Tanaka \cite{JJTan}  and improves
previous results established under the two previous  conditions.
Furthermore, we stress that assumption (A7) ensures
that there exists $s > 0$ such that $G(s) > 0$,
which is a necessary condition for the existence of a solution to \eqref{LE2}
since it allows to exploit Pohozaev's identity then Pohozaev identity
manifold $\mathcal{P}_{0}$ as done by  Berestycki and Lions in
\cite[Proposition 1]{BL}.

As  in \cite{AlvesSM} which we follow, from the moment we extend the existence
of a solution to \eqref{LE2} to the magnetic case ($A(x) \neq 0$),
we also improve previous results obtained  under the two previous
conditions in absence of a magnetic field.
Moreover, we complement all the papers in literature treating equation
\eqref{LE1} in the subcritical case.

\section{Preliminary results}\label{prelimN3}


Here we establish some preliminary  results which will be used
for proving Theorem \ref{groundN3} in Section \ref{sectmainres}.

\begin{lemma}\label{Mmanif}
Under hypothesis  {\rm (A1),  (A4), (A5)--(A7)},  the following assertions hold:
\begin{itemize}
\item[(a)] $\mathcal{M}_A$ is not empty;
\item[(b)] $\mathcal{M}_A$ is a $C^1$ manifold.
\end{itemize}
\end{lemma}

\begin{proof}
(a) By remarks \ref{rmkfF} and \ref {rmkfFq},
for $u\neq 0$ one has $J(tu)<0$ if $t>0$ is small and $J(tu)\to +\infty$
if $t \to +\infty$. Hence, $J(t_0)$=1, for some $t_0>0$.
\smallskip

\noindent (b) By the definitions of  $J$ and $J'$ given in Section \ref{tools},
for every $u \in \mathcal{M}_A$ we have
\begin{align*}
J'(u) u &=   \int_{\mathbb{R}^N} \left( f(|u|^2) |u|^2   -  |u|^2 \right) \, dx \\
&=   \int_{\mathbb{R}^N} \left( f(|u|^2) |u|^2  -  F(|u|^2) \right) \, dx
+ \int_{\mathbb{R}^N} \left( F(|u|^2)  -  |u|^2 \right) \, dx.
\end{align*}
By (A6) and $J(u)=1$, it follows that
\begin{align*}
J'(u) u  \geq   \int_{\mathbb{R}^N} \left( F(|u|^2)   -  |u|^2 \right) \, dx
= 2 \int_{\mathbb{R}^N} G(u) \, dx = 2 > 0
\end{align*}
for all $u \in \mathcal{M}_A$. Then, $J'(u) u \neq 0$ for any $u \in \mathcal{M}_A$.
\end{proof}

\begin{lemma}\label{bdedminseq}
Let assumptions  {\rm (A1)--(A6)}  be satisfied. Then,  any minimizing
sequence $\{ u_n \}$ for $D_A$ is bounded in $H^1_{A}(\mathbb{R}^N, \mathbb{C})$.
The same assertion holds  in particular in $H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$.
\end{lemma}

\begin{proof}
Taken  $\{ u_n \}  $  a minimizing sequence for $D_A$ in
$H^1_{A}(\mathbb{R}^N, \mathbb{C})$, we obtain
$$
\frac 12 \int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx \to D_A \quad \text{as }
n \to + \infty
 $$
and
$$
\int_{\mathbb{R}^N} G(u_n) \, dx = 1, \quad \text{that is }
 \frac 12 \int_{\mathbb{R}^N} \left( F(|u_n|^2)   -  |u_n|^2 \right) \, dx =1.
$$
It follows that
\begin{equation}\label{boundgradA}
\int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx
\leq C \quad \text{for all $n \in \mathbb{N}$ and for some constant $C > 0$ }
\end{equation}
and
$$
\int_{\mathbb{R}^N}  F(|u_n|^2)  \, dx  = 2 + \int_{\mathbb{R}^N}   |u_n|^2 dx.
$$
By  Remark \ref{rmkfF} with $\varepsilon = \frac 12$ we obtain
$$
2 + \int_{\mathbb{R}^N}   |u_n|^2 dx
\leq  \frac 12 \int_{\mathbb{R}^N}   |u_n|^2 dx
+ \overline{C}_{1/2} \int_{\mathbb{R}^N}   |u_n|^{2^*}  \, dx.
$$
Then, by \eqref{boundgradA}, for every $n \in \mathbb{N}$,
$$
\frac 12 \int_{\mathbb{R}^N}   |u_n|^2 dx
\leq  \overline{C}_{1/2} \int_{\mathbb{R}^N}   |u_n|^{2^*}  \, dx
\leq \frac{\overline{C}_{1/2}}{S^{2^* /2}}
 \Big(\int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx \Big)^{2^* /2}\leq \overline C.
$$
Consequently, $\{ u_n \}  $ is bounded in $L^2(\mathbb{R}^N, \mathbb{C})$ and this implies
that  $\{ u_n \}  $ is bounded in $H^1_{A}(\mathbb{R}^N, \mathbb{C})$.
Without any difficulty, these arguments work also in
$H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$.
\end{proof}

Clearly, by Sobolev imbeddings we obtain any minimizing sequence
$\{ u_n \}$ for $D_A$ is bounded also in   $L^m(\mathbb{R}^N, \mathbb{C})$
for every $m \in [2, 2^*]$.
It is useful to establish some lemmas involving the level $D_{0}$ associated
with \eqref{LE2} and  the min-max level $b_A$ of $I_A$.

\begin{lemma}\label{blowbound}
Under assumptions {\rm (A1), (A4)--(A7)}, it holds
$$
\frac 1N \Big(\frac{N-2}{2N} \Big)^{(N-2)/2} (2D_{0})^{N/2} \leq b_A.
$$
\end{lemma}

\begin{proof}
First let us prove  that, by diamagnetic inequality \eqref{diam}, we have
$b_{0} \leq b_A$,
where $b_{0}$ and $b_A$  are respectively  the min-max levels for the
functionals $I_{0}$ and $I_A$, that is
$$
b_{0} = \inf_{\overline{\gamma} \in \Gamma_{0}}
\max_{t \in [0, 1]} I_{0}(\overline{\gamma}(t)) \quad \text{and} \quad
b_{A} = \inf_{\gamma \in \Gamma_A} \max_{t \in [0, 1]} I_A(\gamma(t)).
$$
To do this, we take an arbitrary $\gamma \in \Gamma_A$.
 Then, since $\gamma(t) \in H^1_{A}(\mathbb{R}^N, \mathbb{C})$  for every
$t \in [0, 1]$, by diamagnetic inequality \eqref{diam} we obtain
$| \gamma(t) | \in H^1(\mathbb{R}^N, \mathbb{R})$ and
$$
\int_{\mathbb{R}^N}  |\nabla (|\gamma(t)|) |^2 dx
 \leq  \int_{\mathbb{R}^N}  |\nabla_A (\gamma(t)) |^2 dx
$$
which implies
\begin{equation}\label{maxx1}
 I_{0}(|\gamma(t)|) \leq I_{A}(\gamma(t))
 \end{equation}
for any $t \in [0, 1]$. Therefore,
\begin{equation}\label{maxx2}
\max_{t \in [0, 1]} I_{0}(|\gamma(t)|) \leq \max_{t \in [0, 1]} I_{A}(\gamma(t)).
\end{equation}
Now,  since $\gamma \in C\left([0, 1], H^1_{A}(\mathbb{R}^N, \mathbb{C})\right)$
 we obtain  $|\gamma| \in C\left([0, 1], H^1(\mathbb{R}^N, \mathbb{R})\right)$; moreover,
$\gamma(0)=0$ implies $|\gamma(0)| =0$ and by   $I_A(\gamma(1)) < 0$
and \eqref{maxx1} we obtain $I_{0}(|\gamma(1)|) < 0$. Consequently,
$| \gamma | \in \Gamma_{0}$ which easily gives
$$
b_{0} \leq \max_{t \in [0, 1]} I_{0}(|\gamma(t)|).
$$
Now, taking  the infimum over  all  $\gamma \in \Gamma_A$ by \eqref{maxx2}
we conclude $b_{0} \leq b_A$.
Therefore,  to obtain the claim, it is sufficient to prove that
$$
\frac 1N \Big(\frac{N-2}{2N} \Big)^{(N-2)/2} (2D_{0})^{N/2}  \leq b_{0}.
$$
By exploiting the results obtained in \cite{AlvesSM} (see Section \ref{tools})
 we know in particular that
$$
\frac 1N \Big(\frac{N-2}{2N} \Big)^{(N-2)/2} (2D_{0})^{N/2} = p_{0} =   b_{0}.
$$
For the reader's convenience, we sketch here the proof of
$p_{0} \leq   b_{0}$.  Indeed, from \cite{JJTan}, for each
$\overline \gamma \in \Gamma_{0}$ with
$$
\Gamma_{0} = \{\overline{\gamma} \in C\left([0, 1], H^1(\mathbb{R}^N, \mathbb{R}) \right):
\overline{\gamma}(0) =0 \ \text{and} \ I_{0}(\overline{\gamma}(1)) < 0  \}
$$
it results   $\overline{\gamma}([0, 1]) \cap \mathcal{P}_{0} \neq \emptyset$.
Then, there exists $t_0 \in [0, 1]$ such that
$\overline{\gamma}(t_0) \in \mathcal{P}_0$. So
$$
p_{0} \leq I(\overline{\gamma}(t_0)) \leq \max_{t \in [0, 1]}
I(\overline{\gamma{0}}(t))
$$
implies $p_{0} \leq  b_{0}$. Consequently, since $b_{0} \leq b_A$ we obtain
the result.
\end{proof}

\begin{lemma}\label{Dposit}
Suppose {\rm (A1), (A4)--(A6)} hold. Then the number $D_A$ given
in \eqref{DefD} is positive, namely, $ D_A > 0$.
\end{lemma}

\begin{proof}
By definition  $D_A \geq 0$. Suppose, by contradiction, that $D_A =0$.
If  $\{ u_n \}  $  is a minimizing sequence for $D_A =0$ in
$H^1_{A}(\mathbb{R}^N, \mathbb{C})$, from (A3),  without loss of generality,
 we can suppose that $\{ u_n \}  $  is a minimizing sequence for
$D_A =0$ in $H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$. Then
$$
\frac 12 \int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx \to 0
\quad \text{as } n \to + \infty $$
and
$$
1 = \int_{\mathbb{R}^N} G(u_n) \, dx
= \frac 12 \int_{\mathbb{R}^N} \left( F(|u_n|^2)   -  |u_n|^2 \right) \, dx.
$$
Then, by Remark \ref{rmkfF} it follows that
$$
2 +  \int_{\mathbb{R}^N}  |u_n|^2 dx =  \int_{\mathbb{R}^N} F(|u_n|^2)  \, dx
\leq \varepsilon \int_{\mathbb{R}^N}  |u_n|^2 dx
+ \overline{C}_{\varepsilon} \int_{\mathbb{R}^N}  |u_n|^{2^*}  \, dx
$$
so  we obtain
$$
2 + (1 - \varepsilon) \int_{\mathbb{R}^N}  |u_n|^2 dx
\leq \overline{C}_{\varepsilon} \int_{\mathbb{R}^N}  |u_n|^{2^*}  \, dx
\leq  \frac{\overline{C}_{\varepsilon}}{S^{2^* /2}}  \,
\Big(\int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx \Big)^{2^* /2}
$$
By choosing $\varepsilon = \frac 12$ we obtain
$$
2 \leq      \frac{\overline{C}_{1/2}}{S^{2^* /2}}
\Big(\int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx \Big)^{2^* /2}
   \to 0 \ \text{as } n \to + \infty
$$
which is an absurd.
\end{proof}

\begin{remark}\label{rmkekeland} \rm
Assume (A1), (A4)--(A6) hold. By Ekeland Variational Principle stated
in \cite[Theorem 8.5]{Willem}, we can suppose that the minimizing sequence
$\{ u_n \}  \subset \mathcal{M}_A $ to $D_A$ is a Palais-Smale sequence,
namely, there exists a Lagrange multipliers sequence
$\{ \lambda_n \} \subset \mathbb{R} $  such that
\begin{gather*}
\frac 12  \int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx
\to  D_A \quad \text{as }n \to + \infty, \\
T'_A(u_n) - \lambda_n J'(u_n) \to 0   \text{ in
$H^{-1}_{A}(\mathbb{R}^N, \mathbb{C})$  as $n \to + \infty$},
\end{gather*}
where we recall that, for simplicity of notation,
in Section \ref{tools} we set
$T_A(u) = \frac 12 \int_{\mathbb{R}^N}  |\nabla_A u|^2 dx$
for any $u \in H^1_{A}(\mathbb{R}^N, \mathbb{C})$.
\end{remark}

\begin{lemma}\label{lambdanD}
Let  {\rm (A1), (A4)--(A6)} be satisfied. Then,  the sequence
$\{ \lambda_n \}$ of Lagrange multipliers (see  Remark \ref{rmkekeland})
 is bounded from above and
$$
0 < \liminf_{n \to + \infty} \lambda_n
\leq \limsup_{n \to + \infty} \lambda_n \leq D_A.
$$
\end{lemma}

\begin{proof}
Since $T'_A(u_n) - \lambda_n J'(u_n) \to 0$ in   $H^{-1}_{A}(\mathbb{R}^N, \mathbb{C})$
as $n \to + \infty$, we obtain
\begin{align*}
&T'_A(u_n)u_n - \lambda_n J'(u_n)u_n  \\
&=  \int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx
  - \lambda_n \operatorname{Re} \int_{\mathbb{R}^N} g(u_n) \overline{u_n} \, dx \\
&= \int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx
 - \lambda_n \int_{\mathbb{R}^N} \left(f(|u_n|^2)| u_n|^2 - |u_n|^2 \right) \, dx
 = o_n(1).
\end{align*}
This is equivalent to
\begin{align*}
  \int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx
&=   \lambda_n \int_{\mathbb{R}^N} \left(f(|u_n|^2)| u_n|^2 - F(|u_n|^2) \right)
 \, dx  \\
&  \quad  +  \lambda_n \int_{\mathbb{R}^N} \left(F(|u_n|^2) - |u_n|^2 \right) \, dx
+ o_n(1).
\end{align*}
By (A6) we obtain
\begin{align*}
\int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx
&\geq       \lambda_n \int_{\mathbb{R}^N} \left(F(|u_n|^2) - |u_n|^2 \right) \, dx
  + o_n(1) \\
 &=   2 \lambda_n \int_{\mathbb{R}^N} G(u_n) \, dx + o_n(1) = 2 \lambda_n + o_n(1).
\end{align*}
Then, from Remark \ref{rmkekeland}
$$
 \limsup_{n \to + \infty} \lambda_n
\leq \frac 12 \limsup_{n \to + \infty} \int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx =D_A
 $$
and the right-hand inequality in the thesis  follows. \\
Now, we prove that $ \liminf_{n \to + \infty} \lambda_n > 0$.
First observe that,  by Remark \ref{rmkfF}, Sobolev imbeddings and
the boundedness of the  minimizing sequence $\{u_n \}$ stated
in Lemma \ref{bdedminseq}, we obtain
\begin{align*}
|J'(u_n) u_n |
& \leq  \int_{\mathbb{R}^N} \Big( |f(|u_n|^2|) \,  |u_n|^2 + |u_n|^2   \Big) \, dx  \\
& \leq  \int_{\mathbb{R}^N} \left( (\varepsilon +1) |u_n|^2
+ C_{\varepsilon} |u_n|^{2^*} \right) \, dx \leq C.
\end{align*}
Then, we  conclude that
$$
2 T_A(u_n) = T'_A(u_n) u_n = \lambda_n J'(u_n)u_n  +  o_n(1)
$$
and $2 T_A(u_n) \to 2 D_A > 0$ where   the positivity of $D_A$
has been established in Lemma \ref{Dposit}.
\end{proof}

As observed in the following remark,  to study the compactness of a
minimizing sequence $\{u_n \}$ to $D_A$, it is sufficient to consider
 the sequence $\{|u_n| \}$ of the  modula of $u_n $.
So we can exploit the diamagnetic inequality \eqref{diam} and apply
to $\{|u_n| \}$ the Concentration-Compactness technique of
Lions \cite{Lionsibero} which is based on a measure representation.

\begin{remark}\label{rmklions} \rm
If a sequence $\{ u_n \}$ is bounded in $H^1_{A}(\mathbb{R}^N, \mathbb{C})$,
then by diamagnetic inequality \eqref{diam} the sequence
 $\{| u_n| \}$ of its modula is bounded in $H^1(\mathbb{R}^N, \mathbb{R})$.
Then, by Lions Concentration Compactness principle \cite[Lemma 1.2]{lionsI},
there are a countable index set $\Lambda$, nonnegative  finite measures
$\mu$ and $\nu$ and families $\{ \mu_i \}, \{ \nu_i \} \subset  (0, + \infty)$
and $\{ x_i \} \subset \mathbb{R}^N$ such that
\begin{itemize}
\item[(j)] $|\nabla |u_n||^2 \rightharpoonup \mu \geq |\nabla |u||^2
 +  \sum_{i \in \Lambda} \delta_{x_i} \mu_i $ (weak$^*$ sense of measures);

\item[(jj)] $ |u_n|^{2^*} \rightharpoonup \nu =  | u|^{2^*}
 +  \sum_{i \in \Lambda} \delta_{x_i} \nu_i $ (weak$^*$ sense of measures)

\item[(jjj)] $\mu_i \geq S \, \nu_i^{2/2^*}$ for every $i \in \Lambda$.
\end{itemize}
 where $|u|$ is the  weak limit of $|u_n|$ and $\delta_{x_i}$ are Dirac
measures at $x_i$.  This remark can be employed in the next lemma.
\end{remark}

\begin{lemma}\label{nuSD}
Suppose that {\rm  (A1), (A4)--(A6)} are satisfied.
If $\nu_i > 0$ for some index $i$, then
$$
\nu_i \geq \Big(\frac{S}{D_A}\Big)^{N/2}.
$$
\end{lemma}

\begin{proof}
Let $\{u_n \}$ be a minimizing sequence to $D_A$ which we know it is
bounded by  Lemma \ref{bdedminseq}. Let
 $\varphi \in C_0^{\infty}(\mathbb{R}^N, [0, 1])$ be such that
\[
\varphi(x)= \begin{cases}
1 &  \text{if } x \in B_1(0) \\
0 &   \text{if } x \in \mathbb{R}^N \setminus B_2(0)
\end{cases}
\]
and $| \nabla \varphi (x)| \leq 2$ for every $x \in \mathbb{R}^N$.
 Then, we can consider
$$
\varphi_{\varepsilon, x_i}(x) = \varphi \Big( \frac{x - x_i}{\varepsilon} \big)
$$
for $\varepsilon > 0$ and $ x_i $ a  singular point of the measures
 $ \sum_{i \in \Lambda} \delta_{x_i} \mu_i $   and
$ \sum_{i \in \Lambda} \delta_{x_i} \nu_i $   whose existence is ensured
in Remark \ref{rmklions} by Lions Concentration Compactness Principle.
Clearly, $\varphi_{\varepsilon, x_i} \in C_0^{\infty}(\mathbb{R}^N, [0, 1])$
satisfies
\[
\varphi_{\varepsilon, x_i}(x)= \begin{cases}
1 &  \text{if } x \in B_{\varepsilon}(x_i) \\
0 &  \text{if } x \in \mathbb{R}^N \setminus B_{2 \varepsilon}(x_i)
\end{cases}
\]
and $| \nabla \varphi_{\varepsilon, x_i} (x)| \leq 2/ \varepsilon$
for every $x \in \mathbb{R}^N$. Since $\{u_n \varphi_{\varepsilon, x_i}\}$
is bounded in  $ H^1_{A}(\mathbb{R}^N, \mathbb{C})$, by Remark \ref{rmkekeland} we obtain
 $$
T'_A(u_n)(u_n \varphi_{\varepsilon, x_i})
= \lambda_n J'(u_n)(u_n \varphi_{\varepsilon, x_i}) + o_n(1);
$$
that is,
\begin{equation}\label{reunphie}
\begin{aligned}
&\operatorname{Re}  \int_{\mathbb{R}^N}
 \nabla_A u_n \cdot \overline{\nabla_A (u_n \varphi_{\varepsilon, x_i}) }  \, dx  \\
& = \lambda_n \, \operatorname{Re}   \int_{\mathbb{R}^N}
\left(f(|u_n|^2) |u_n|^2 - |u_n|^2 \right)
 \overline{ \varphi_{\varepsilon, x_i} }  \, dx +  o_n(1).
\end{aligned}
\end{equation}
Since $\varphi_{\varepsilon, x_i}$ takes real values and by direct
calculations,
$$
\overline{\nabla_A (u_n\varphi_{\varepsilon, x_i}) }
= \overline{\nabla_A u_n} \varphi_{\varepsilon, x_i}
+ \overline{u_n} \nabla \varphi_{\varepsilon, x_i}\,.
$$
Then we can write the term on the left-hand side in \eqref{reunphie} as
\begin{equation} \label{rerere}
\begin{aligned}
&\operatorname{Re}  \int_{\mathbb{R}^N}  \nabla_A u_n
\cdot \overline{\nabla_A (u_n \varphi_{\varepsilon, x_i}) }  \, dx \\
& = \operatorname{Re}  \int_{\mathbb{R}^N}  \nabla_A u_n \cdot
\left(  \overline{\nabla_A u_n} \varphi_{\varepsilon, x_i}
 + \overline{u_n}  \nabla \varphi_{\varepsilon, x_i} \right)   \, dx \\
&=   \int_{\mathbb{R}^N}  |\nabla_A u_n|^2  \varphi_{\varepsilon, x_i} \, dx
+  \operatorname{Re}  \int_{\mathbb{R}^N} \overline{u_n}    \nabla_A u_n \cdot
 \nabla \varphi_{\varepsilon, x_i}  \, dx.
\end{aligned}
\end{equation}
Now, observe that
\begin{equation}\label{reconj}
\operatorname{Re} (\overline u_n   \nabla_A u_n) = |u_n| \nabla |u_n|.
\end{equation}
Indeed,
\begin{align*}
\operatorname{Re} (\overline u_n \nabla_A u_n)
&=  \operatorname{Re} \left(\overline u_n
  \left(\nabla u_n + i A(x) u_n \right) \right) \\
&=  \operatorname{Re} \left(\overline u_n \nabla u_n + i A(x) |u_n|^2 \right)
  = \operatorname{Re} (\overline u_n \nabla u_n) \\
&=  |u_n| \, \operatorname{Re} \Big(\frac{\overline u_n}{|u_n|} \nabla u_n \Big)
= |u_n| \nabla |u_n|.
 \end{align*}
By substituting \eqref{reconj} in \eqref{rerere} and by replacing in
turn \eqref{reconj} in \eqref{reunphie} we obtain by  diamagnetic inequality
\begin{equation}\label{diamrere}
\begin{aligned}
&\int_{\mathbb{R}^N}  |\nabla |u_n||^2 \varphi_{\varepsilon, x_i} \, dx
 + \int_{\mathbb{R}^N}  |u_n| \nabla |u_n| \cdot \nabla \varphi_{\varepsilon, x_i}
  \, dx  \\
&\leq     \lambda_n   \int_{\mathbb{R}^N}   f(|u_n|^2) |u_n|^2
 \varphi_{\varepsilon, x_i}   \, dx + \lambda_n \int_{\mathbb{R}^N}
  |u_n|^2 \varphi_{\varepsilon, x_i}  \, dx   + o_n(1).
\end{aligned}
\end{equation}
Note that by (A4), there exists $\delta>0$ such that
$$
|f(s)|\leq 1, \quad \text{for all } 0\leq  s\leq \delta.
$$
From (A5), there exists $K>0$ such that
$$
|f(s)|\leq |s|^{(2^* -2)/2}, \quad \text{for all }  s\geq K.
$$
Moreover, for $s\in [\delta,k]$, we obtain
$$
|f(s)|\leq K_{1}, \quad \text{for some }  \ K_1 >0.
$$
Then
$$
|f(s)| \leq   (1+K_1) +  |s|^{(2^* -2)/2} ,   \quad\text{for every } s \geq 0.
$$
In particular,
$$
f(|u_n|^2) |u_n|^2 \leq  (1+K_1) |u_n|^2  +   |u_n|^{2^*},
 \quad\text{for every } n \in \mathbb{N}.
$$
Therefore, by \eqref{diamrere} it follows that
\begin{align*}
&\int_{\mathbb{R}^N}  |\nabla |u_n||^2 \varphi_{\varepsilon, x_i} \, dx
 + \int_{\mathbb{R}^N}  |u_n| \nabla |u_n| \cdot \nabla \varphi_{\varepsilon, x_i}
 \, dx \\
& \leq  \lambda_n \Big( (2+K_1) \int_{\mathbb{R}^N}
  |u_n|^2 \varphi_{\varepsilon, x_i}  \, dx
 +       \int_{\mathbb{R}^N}   |u_n|^{2^*} \varphi_{\varepsilon, x_i}  \, dx
 \Big) +  o_n(1).
\end{align*}
Now, we prove that
\begin{equation}\label{limsupen}
\limsup_{\varepsilon \to 0} \limsup_{n \to + \infty} \int_{\mathbb{R}^N}
|u_n| \nabla |u_n| \cdot \nabla \varphi_{\varepsilon, x_i}  \, dx =0.
\end{equation}
Indeed, by H\"older's inequality,
\begin{align*}
&\limsup_{n \to + \infty} \Big(\int_{\mathbb{R}^N}  |u_n| \nabla |u_n|
 \cdot \nabla \varphi_{\varepsilon, x_i}  \, dx \Big)  \\
&\leq \limsup_{n \to + \infty} \Big(\int_{\mathbb{R}^N}   |\nabla |u_n||\,
 |u_n| \, |\nabla \varphi_{\varepsilon, x_i}|  \, dx \Big)\\
&\leq  \limsup_{n \to + \infty}
 \Big(\int_{\mathbb{R}^N}   |\nabla |u_n||^2 dx \Big)^{1/2}
 \Big(\int_{\mathbb{R}^N} |u_n|^2  |\nabla \varphi_{\varepsilon, x_i}|^2 dx
 \Big)^{1/2}  \\
&\leq   C_1  \limsup_{n \to + \infty}
 \Big(\int_{\mathbb{R}^N} |u_n|^2  |\nabla \varphi_{\varepsilon, x_i}|^2 dx
 \Big)^{1/2} \\
&=  C_1 \Big(\int_{\mathbb{R}^N} |u|^2  |\nabla \varphi_{\varepsilon, x_i}|^2 dx \Big)^{1/2}
\end{align*}
where $ C_1 = \sup_{n} \big(\int_{\mathbb{R}^N}   |\nabla |u_n||^2 dx \big)^{1/2}$.
 By using H\"older's inequality again it follows that
\begin{align*}
\int_{\mathbb{R}^N} |u|^2  |\nabla \varphi_{\varepsilon, x_i}|^2 dx
&\leq \Big(\int_{B_{2 \varepsilon}(x_i)} |u|^{2^*} \, dx \Big)^{2/2^*}
 \Big(\int_{\mathbb{R}^N}   |\nabla \varphi_{\varepsilon, x_i}|^N \, dx
 \Big)^{2/N} \\
&\leq C_N  \Big(\int_{B_{2 \varepsilon}(x_i)} |u|^{2^*} \, dx \Big)^{2/2^*}
\end{align*}
with $C_N > 0$ a suitable constant depending on $N$.
Letting $\varepsilon \to 0$ we obtain \eqref{limsupen} is satisfied.
Now, since
$u_n \to u$ in $L^2_{\rm loc}(\mathbb{R}^N)$ it easily follows
$$
\limsup_{\varepsilon \to 0} \limsup_{n \to + \infty}
\int_{\mathbb{R}^N}  |u_n|^2 \varphi_{\varepsilon, x_i}  \, dx =0
$$
which together with  the assertion
$$
\limsup_{n \to + \infty} \lambda_n \leq D_A,
$$
in Lemma \ref{lambdanD},
allow us to have
$ \mu_i \leq  D_A \,  \nu_i$.
By Remark \ref{rmklions} (jjj) we have $  S \nu_i^{2/2^*} \leq \mu_i$.
Then
$$
S \nu_i^{2/2^*} \leq D_A  \nu_i
$$
and $\nu_i \geq \big(\frac{S}{D_A}\big)^{N/2}$.
This completes the proof.
\end{proof}

\begin{lemma}\label{D2S}
Suppose that {\rm (A1),  (A4)--(A6)} are satisfied. 
If $\nu_i > 0$ for some index $i$, then
$D_A \geq 2^{- 2/N} S$.
\end{lemma}

\begin{proof}
From  $ 0 \leq  \varphi_{\varepsilon, x_i} \leq 1$, it follows that
$$ 
\int_{\mathbb{R}^N}   |u_n|^{2^*} \varphi_{\varepsilon, x_i} \, dx 
\leq \int_{\mathbb{R}^N}  |u_n|^{2^*}  \, dx 
\leq S^{- 2^{*}/2}  \Big( \int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx \Big)^{2^*/2}.
$$
Passing to the limit as $n \to + \infty$, we have
$$ 
\nu_i \leq  S^{- 2^{*}/2} \left(2D_A \right)^{2^{*}/2}.
$$
From Lemma \ref{nuSD}, since $\nu_i \geq \left(\frac{S}{D_A}\right)^{N/2}$, 
it follows that $D_A \geq 2^{- 2/N} S$;
thus completing the proof.
\end{proof}

Now, let us recall the next result relating the constants $D_A$ and $D_{0}$.

\begin{lemma}\label{DAD0r}
Under assumptions {\rm (A1), (A4)--(A6)}, we have
$D_A = D_{0}$.
\end{lemma}


\begin{proof}
We follows the arguments used in \cite{ariolisz, EstLions}, and 
for the sake of completeness, we give here the details of the proof. 
By diamagnetic inequality \eqref{diam} we obtain
$$
D_{0} \leq  \int_{\mathbb{R}^N} \left( |\nabla |u||^2 + |u|^2 \right)   \, dx  
\leq \int_{\mathbb{R}^N} \left( |\nabla_A u|^2 + |u|^2 \right)   \, dx , 
$$
which implies $D_{0} \leq D_A $. 

Now, we show the reversed inequality $D_A  \leq D_{0}$ also holds. 
Taking  $\varepsilon > 0$ an arbitrarily small constant,  we consider  
$\varphi_{\varepsilon} \in C_c^{\infty}(\mathbb{R}^N, \mathbb{R})$  
whose $\operatorname{supp}(\varphi_{\varepsilon})$  is a (compact) neighborhood 
of $x_0 \in \mathbb{R}^N$ (for simplicity, we can assume $x_0 =0$) satisfying
$$
\int_{\mathbb{R}^N} G(\varphi_{\varepsilon}) \, dx
 = 1 \quad \text{and} \quad 
 \int_{\mathbb{R}^N} |\nabla \varphi_{\varepsilon}|^2 dx
\leq D_{0} + \varepsilon. 
$$
We can define a function $u_{\varepsilon} = e^{i \chi_{x_0}(x)} \varphi_{\varepsilon}$ 
with $\chi_{x_0}(x)= A(x_0) \cdot x$. 
Since $\varphi_{\varepsilon} \in H^1(\mathbb{R}^N, \mathbb{R})$, for a direct calculation,  
we obtain $u_{\varepsilon} \in H^1_A(\mathbb{R}^N, \mathbb{C})$ and
\begin{align*}
\int_{\mathbb{R}^N} G(u_{\varepsilon}) \, dx 
&= \int_{\mathbb{R}^N} G( e^{i \chi_{x_0}(x)} \varphi_{\varepsilon}) \, dx \\
&= \frac 12 \int_{\mathbb{R}^N} \Big( F( | e^{i \chi_{x_0}(x)} 
 \varphi_{\varepsilon} |^2) - | e^{i \chi_{x_0}(x)} \varphi_{\varepsilon} |^2 \Big)
  \, dx \\
&= \frac 12 \int_{\mathbb{R}^N} \Big( F( |  \varphi_{\varepsilon} |^2) 
 - | \varphi_{\varepsilon} |^2 \Big) \, dx \\
&= \int_{\mathbb{R}^N} G(\varphi_{\varepsilon}) \, dx =1.
\end{align*}
Now,  by the continuity assumption in  (A1) we obtain 
$|A(x) + A(x_0)|^2 \leq c$ in the $\operatorname{supp}(\varphi_{\varepsilon})$, 
if we suppose $\|\varphi_{\varepsilon} \|_2 = o(\varepsilon)$,  we deduce
\begin{align*}
D_A \leq \int_{\mathbb{R}^N} |\nabla_A u_{\varepsilon}|^2 dx
&=  \int_{\mathbb{R}^N} |\nabla \varphi_{\varepsilon}|^2 dx
 + \int_{\mathbb{R}^N} |A(x) + A(x_0)|^2 | \varphi_{\varepsilon}|^2 dx \\
& \leq  D_{0} + o(\varepsilon)\,
\end{align*}
which completes the proof.
\end{proof}

From Lemmas \ref{D2S} and \ref{DAD0r} we obtain easily the next result.

\begin{lemma}\label{DAD0r1}
Let  {\rm  (A1),  (A4)--(A6)} be satisfied. If $\nu_i > 0$ for some index $i$, then
$D_{0} \geq   D_A \geq 2^{- 2/N} S$.
\end{lemma}

\begin{lemma}\label{lambdabbound}
Suppose {\rm (A1),  (A4)--(A7)} be satisfied. If
\begin{equation}\label{lamcq}
 \lambda > \Big[ 2^{(2-N)/2} S^{-N/2} N \Big(\frac{2N}{N-2} \Big)^{(N-2)/2}
\Big]^{(q-2)/2}
 c_{A}^{(q-2)/2},
\end{equation}
then
$$
b_A <  \frac 1N \Big(\frac{N-2}{2N} \Big)^{(N-2)/2} 2^{(N-2)/2} S^{N/2}.
$$
\end{lemma}

\begin{proof}
Take $\varphi \in H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$ a solution of
 \eqref{Subcritical}. From the definition 
$  b_A = \inf_{\gamma \in \Gamma_A} \max_{t \in [0, 1]} I_A(\gamma(t))$,
 \eqref{Sara1}, \eqref{Sara2} and (A7), it follows that
\[
b_A \leq \max_{ t \geq 0} I_A(t \varphi) 
 \leq  \max_{t \geq 0} \big\{ \frac{t^2}{2} - \lambda \frac{t^q}{q} \big\}
c_{A}\frac{2q}{q-2} 
=   \frac{c_{A}}{\lambda^{2/(q-2)}}.
\]
By using the lower bound on  $\lambda$ assumed in  hypothesis \eqref{lamcq}
we obtain
\begin{align*}
b_A <  \frac 1N \Big(\frac{N-2}{2N} \Big)^{(N-2)/2} 2^{(N-2)/2} S^{N/2}\,.
\end{align*}
This completes the.
\end{proof}


\begin{lemma}\label{weaklim}
Assume  {\rm (A1)--(A7)} are satisfied. If  \eqref{lamcq} holds, namely 
$$ 
\lambda >  \Big[ 2^{(2-N)/2} S^{-N/2} N \Big(\frac{2N}{N-2} \Big)^{(N-2)/2}
\Big]^{(q-2)/2}  c_{A}^{(q-2)/2},
$$
then the weak limit $u$ of any  minimizing sequence $\{ u_n \}$ to $D_A$ 
is nontrivial.
\end{lemma}

\begin{proof}
Let $\{ u_n \}$ be  a minimizing sequence to $D_A$. 
Then Lemma \ref{bdedminseq} states $\{ u_n \}$ is bounded in
$H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$. Then
$\{ |u_n| \}$ is bounded in $H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{R})$ and 
there exists $u \in H^1_{A}(\mathbb{R}^N, \mathbb{C})$ such that 
$|u_n| \rightharpoonup |u|$ in $H^1_{A}(\mathbb{R}^N, \mathbb{R})$.  
Suppose by contradiction that  $|u| =0$.  By  Remark \ref{rmklions} (jj)
based on  Lions Concentration Compactness principle we obtain
\begin{equation}\label{sumdeltai}
|u_n|^{2^*} \rightharpoonup d\nu =  \sum_i \delta_{x_i} \nu_i   
\quad \text{(in the weak$^*$ sense of measures)}.
\end{equation}
Since $\{|u_n|\} \subset  H^1_{rad}(\mathbb{R}^N, \mathbb{R})$, 
 by  \cite[Radial Lemma A.II]{BL} there exist a radius $R = R(N) > 0$ 
and a constant $C = C(N) > 0$ both independent of $n$  such that
$$
|u_n(x)| \leq C \, |x|^{-(N-1)/2} \quad \text{for }|x| \geq R, 
$$
or equivalently
$$
|u_n(r)| \leq C \, r^{-(N-1)/2} \quad \text{for } r \geq R. 
$$
Then the  sequence  $ \{ |u_n| \}$ is bounded in $L^{\infty}(B_R^c(0), \mathbb{R})$  
for every $R >0$ or equivalently there exists a constant $M > 0$ such that
\begin{equation}\label{boundinfty}
 \|u_n \|_{L^{\infty}(B_R^c(0), \mathbb{R})} \leq M \quad\text{for every $n \in \mathbb{N}$}.
 \end{equation}
This implies  $ \{ |u_n| \}$ converges strongly to $0$ in 
$L^{m}(B_R^c(0), \mathbb{R})$ for all $m > 2$ and for any $R > 0$.
We prove that also $\nu_{i_0} =0$. If on the contrary $\nu_{i_0}  > 0$,  
by Lemma \ref{DAD0r1}  we obtain $D_{0} \geq D_A \geq 2^{- 2/N} S$. 
Since by Lemma \ref{blowbound},
$$
\frac 1N \Big(\frac{N-2}{2N} \Big)^{(N-2)/2} (2D_{0})^{N/2} \leq b_A,
$$
we obtain
$$ 
\frac 1N \Big(\frac{N-2}{2N} \Big)^{(N-2)/2} 2^{(N-2)/2} S^{N/2} \leq b_A.  
$$
But this last inequality  contradicts Lemma \ref{lambdabbound}. 
Then there is no $\nu_i > 0$ for every $i \in \Lambda$. 
Consequently, by \eqref{sumdeltai},
$$  
|u_n|^{2^*} \rightharpoonup 0   \quad\text{(in the weak$^*$ sense of  measures)},
$$
hence
$$ 
\int_{\mathbb{R}^N} |u_n|^{2^*} \varphi \, dx \to  0,   \quad\text{for every }
\varphi \in C_0^{\infty}(\mathbb{R}^N, \mathbb{C}). 
$$
This implies
\begin{equation}\label{2*loc}
 u_n \to 0    \quad\text{in $ L_{\rm loc}^{2^*}(\mathbb{R}^N, \mathbb{C})$}.
 \end{equation}
Using the same argument,  we have
$$  
u_n \to 0    \quad\text{in $ L^{2^*}(B_R^c(0), \mathbb{C})$} 
\quad\text{for any }R > 0, 
 $$
which together with  \eqref{2*loc} implies
$u_n \to 0$ in $ L^{2^*}(\mathbb{R}^N, \mathbb{C})$. 

Now, we can follow again  the arguments used in Lemma \ref{Dposit}.
 Indeed, since $\{u_n\} \subset \mathcal{M}_A$ implies
$$ 
1 = \int_{\mathbb{R}^N} G(u_n) \, dx
 = \frac 12 \int_{\mathbb{R}^N} \left( F(|u_n|^2)   -  |u_n|^2 \right) \, dx, 
$$
by Remark \ref{rmkfF} it follows that
$$
2 +  \int_{\mathbb{R}^N}  |u_n|^2 dx
=  \int_{\mathbb{R}^N} F(|u_n|^2)  \, dx 
\leq \varepsilon \int_{\mathbb{R}^N}  |u_n|^2 dx
 + C_{\varepsilon} \int_{\mathbb{R}^N}  |u_n|^{2^*}  \, dx  
$$
so
$$
2 + (1 - \varepsilon) \int_{\mathbb{R}^N}  |u_n|^2 dx
\leq C_{\varepsilon}  \int_{\mathbb{R}^N} |u_n|^{2^*}  \, dx.
$$
By choosing $\varepsilon = 1/2$ we obtain
$$
2 \leq C_{1/2}  \int_{\mathbb{R}^N}  | u_n|^{2^*}  \, dx \to 0 \quad
 \text{as $n \to + \infty$} $$
which is an absurd. Thus, we can conclude that $u \neq 0$.
\end{proof}

\section{Proof of Theorem \ref{groundN3}}\label{sectmainres}


Under the  assumptions in Theorem \ref{groundN3} we show that $D_A$ 
is attained by $u$, where $u$ is the non trivial weak limit of the 
minimizing sequence $\{ u_n \}$ to $D_A$. 
Indeed, since $\{u_n \}$ is bounded  by  Lemma \ref{bdedminseq}, 
we have   $u_n \rightharpoonup u$ in $H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$ 
and being  the weak limit $u$  not trivial thanks to  Lemma \ref{weaklim}, we deduce
that
\begin{equation}\label{TD}
T_A(u) = \frac 12\int_{\mathbb{R}^N}  |\nabla_A u|^2 dx
\leq \liminf_{n \to + \infty} \frac 12 \int_{\mathbb{R}^N}  |\nabla_A u_n|^2 dx =D_A.
\end{equation}
Now, by  Lemmas \ref{blowbound}, \ref{DAD0r1} and  \ref{lambdabbound}
 we obtain that $\nu_i =0$ for every $i$. 
It remains to prove that $u \in \mathcal{M}$. 
To do this, first observe that  the uniform decay at infinity of  
$\{u_n\} \subset  H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$  
together with (A4) imply the existence of a radius $R > 0$ such that
$$
|u_n|^2 - F(|u_n|^2) \geq 0 \quad 
\text{for any $n \in \mathbb{N}$ and in $\mathbb{R}^N \setminus B_R$} 
$$
where $B_R$ denotes the ball of radius $R$ centered in $0$. 
Since $u_n \to u \ \text{in $L^{2^*}(B_R, \mathbb{C})$}$,
from \cite[Theorem 4.9, Section 4]{Brezis}, up to a subsequence,  
$u_{n} \to u $ a.e. in $B_{R}$ and there exists $v\in L^{2^{*}}(B_{R})$ 
such that $|u_{n}(x)|\leq v(x)$ a.e.\ in $B_{R}$.
Moreover we have $F(|u_{n}(x)|^2) \to F(|u(x)|^2)$  a.e.\ and, by 
 Remark \ref{rmkfF}, in correspondence of  any $\varepsilon > 0$ 
we obtain the existence of $\overline{C}_{\varepsilon} > 0$ such that
$$ 
|F(|u_n|^2) | \leq \varepsilon |u_n|^2 
+ \overline{C}_{\varepsilon} |u_n|^{2^*} 
\leq \varepsilon |v|^2 + \overline{C}_{\varepsilon} |v|^{2^*} 
$$
for all $n \in \mathbb{N}$. By the arbitrariness of $\varepsilon$ and the
Dominated Convergence Theorem, we obtain
$$
\int_{B_{R}} F(|u_{n}|^2)\, dx \to \int_{B_{R}} F(|u|^2) \, dx. 
$$
Now, since
$$
\frac 12 \int_{B_R}F(|u_n|^2)\,dx
= \frac{1}{2}\int_{B_R} |u_n|^{2} \, dx
+ \frac 12 \int_{\mathbb{R}^N \setminus B_R}\left(|u_n|^{2}-F(|u_n|^2)\right) \, dx +1,
$$
taking into account the above considerations, the  properties 
of limit inferior with respect to the sum of sequences and  Fatou's Lemma
 we infer that
$$
\int_{B_R} F(|u|^2)\, dx 
\geq \frac{1}{2}\int_{B_R} |u^{2}| \,dx
+   \frac 12 \int_{\mathbb{R}^N \setminus B_R}\left(|u|^{2}-F(|u|^2)\right)\,dx +1;
$$
that is,
$ \int_{\mathbb{R}^N}  G(u) \, dx  \geq 1$.
If we prove that also that
\begin{equation}\label{minug}
 \int_{\mathbb{R}^N}  G(u) \, dx  \leq  1,
 \end{equation}
 then we obtain
$u \in \mathcal{M}_A$ and $ T_A(u) =D_A $,
or equivalently,
\[
T_A(u) =D_A 
=  \min \Big\{\frac 12 \int_{\mathbb{R}^N}  |\nabla_A u|^2 dx:
 u \in H^1_{A}(\mathbb{R}^N, \mathbb{C})\setminus \{0 \},\; 
  \int_{\mathbb{R}^N} G(u) \, dx =1 \Big\}.
\]
To show \eqref{minug}, suppose by contradiction that
$$ 
\int_{\mathbb{R}^N}  G(u) \, dx  > 1.
$$
We define $h{:} [0, 1] \to \mathbb{R}$ by
$$
h(t) = \int_{\mathbb{R}^N}  G(tu) \, dx 
=  \frac 12 \int_{\mathbb{R}^N}  \left(F(|tu|^2) - |tu|^2 \right)  \, dx 
\quad \text{for every }t \in [0, 1]. 
 $$
Now we show that  $h(t) < 0$ for $t$ close to $0$. 
Indeed, by Remark \ref{rmkfF} we obtain
\begin{align*}
h(t) &=  \frac 12 \int_{\mathbb{R}^N}  \left(F(t^2  |u|^2)
-  t^2   |u|^2 \right)  \, dx \\
& \leq  \frac 12   \int_{\mathbb{R}^N} \left( \varepsilon 
  t^2  |u|^2 + \overline{C}_{\varepsilon}   t^{2^*}   |u|^{2^*}  \right) \, dx
  - \frac 12   \int_{\mathbb{R}^N} t^2  |u|^2 dx \\
&=  \frac 12   \overline{C}_{\varepsilon}  t^{2^*}\int_{\mathbb{R}^N}  |u|^{2^*}   
\, dx - \frac 12 \,   t^2  (1- \varepsilon)  \int_{\mathbb{R}^N}  |u|^2 dx.
\end{align*}
Choosing $\varepsilon > 0$ sufficiently small, e.g.\
  $\varepsilon < 1$, we obtain    $h(t) < 0$ for $ t > 0$ small enough. 
 Clearly, $h(1) = \int_{\mathbb{R}^N}  G(u) \, dx > 1$. 
Then, by the continuity of $h$, there exists $t_0 \in (0, 1)$ such that 
$h(t_0)=1$ which gives
$$
\int_{\mathbb{R}^N}  G(t_0 u) \, dx = 1      \Longleftrightarrow 
t_0 u \in \mathcal{M}_A.
$$
Consequently, by \eqref{TD},
$$
D_A \leq T_A(t_0 u) = \frac{t_0^2}{2} \int_{\mathbb{R}^N} |\nabla_A u|^2 dx 
= t_0^2 T_A(u) \leq t_0^2 D_A < D_A $$
which is absurd. Thus,
$T_A(u) =D_A$ and $ u \in \mathcal{M}_A$;
that is,
$$
D_A = \frac{1}{2} \int_{\mathbb{R}^N} |\nabla_A  u |^2 dx, \quad 
J( u)= \int_{\mathbb{R}^N} G( u) \, dx = \frac 12 \int_{\mathbb{R}^N}  
\left(F(|u|^2) - | u|^2 \right) \, dx =1.
$$
By Lagrange Multipliers Theorem, there exists a multiplier $\theta \in \mathbb{R}$ such that
$$ 
T'_A( u) = \theta  J'( u); 
$$
in particular, for every $v \in H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$ we obtain
$ T'_A( u)v = \theta  J'( u)v$,
namely
\begin{equation}\label{rethetaequat}
\operatorname{Re} \int_{\mathbb{R}^N} \nabla_A  u  \cdot \overline{\nabla_A v }  
\, dx  = \theta  \operatorname{Re} \int_{\mathbb{R}^N} \left(f(| u|^2) u -  u \right) 
\overline v \, dx.
\end{equation}
By adapting  the arguments in Berestycki and Lions \cite{BL}, we are able 
to prove that $\theta > 0$. Indeed, first remark that $\theta \neq 0$; 
if not, namely if $\theta = 0$ we would have  $ T'_A( u) =0$ and in
 particular $\int_{\mathbb{R}^N} |\nabla_A  u |^2 dx = 0$. 
Therefore, $u =0$ which is impossible. 

Specifically, it results that $\theta > 0$. Indeed, suppose  by 
contradiction that $\theta < 0$. Moreover, observe that $J'(u) \neq 0$; otherwise,
$$
J'(u) v =  \operatorname{Re} \int_{\mathbb{R}^N}  
\left(f(|u|^2)u - u \right) \overline v \, dx =0
$$
would imply $f(|u|^2)u -u=0 $ then $F(|u|^2) -|u|^2=0 $ which leads to 
a  contradiction with $J(u) =1$. 

Now let us consider a test function $w$ such that
$$
J'(u) w =  \operatorname{Re} \int_{\mathbb{R}^N}  
\left(f(|u|^2)u - u \right) \overline w \, dx > 0.
 $$
Since $J(u + \varepsilon w) \cong J(u) + \varepsilon J'(u) w$ and
$$
T_A(u + \varepsilon w) \cong T_A(u) +  \varepsilon \theta J'(u) w 
\quad \text{for $\varepsilon \to 0$  and $\theta < 0$}, 
$$
it is possible to choose $\varepsilon > 0$ small enough so that 
$v = u + \varepsilon w$ satisfies $J(v) > J(u) =1$ and
$T_A(v) < T_A(u) = D_A$. Now, by a scale change 
$v_{\sigma}(x) = v(x/ \sigma)$, there exists $0 < \sigma < 1$ such that
$$
J(v_{\sigma}) =  \frac 12 \int_{\mathbb{R}^N}  
\left(F(|v_{\sigma}|^2) - | v_{\sigma}|^2 \right) \, dx  = \sigma^N \, J(v) = 1
$$
and, thanks to assumption (A2) we obtain
$$
T_A(v_{\sigma}) = \frac{1}{2} \int_{\mathbb{R}^N} |\nabla_A v_{\sigma}|^2 dx 
= \sigma^{N-2} T_A(v) < D_A
$$
which is absurd. Then, $\theta > 0$.  
Then   $u$  in $H^1_{A, {\rm rad}}(\mathbb{R}^N, \mathbb{C})$
satisfies (in the weak sense) 
\begin{equation}\label{thetaequat}
 - \Delta_A   u    = \theta \left( f(| u|^2) u - u \right).
 \end{equation}
Now, we aim to prove that by exploiting a suitable change of 
variable the re-scaled $ u$, say ${u}_{\theta}$, satisfies
\begin{equation}\label{reuthetaequat}
\operatorname{Re} \int_{\mathbb{R}^N} \nabla_A  u_{\theta}  
\cdot \overline{\nabla_A v }  \, dx  
=   \operatorname{Re} \int_{\mathbb{R}^N} \left(f(| u_{\theta}|^2) u_{\theta} 
-  u_{\theta} \right) \overline v \, dx
\end{equation}
namely, ${u}_{\theta}$ satisfies (in the weak sense)
\begin{equation}\label{uthetaequat}
 - \Delta_A  { u}_{\theta}+  { u}_{\theta} = f(|{ u}_{\theta}|^2) { u}_{\theta}
 \end{equation}
so that  ${ u}_{\theta}$ is a solution to \eqref{LE1}.
Indeed, since
\begin{align*}
\operatorname{Re} \int_{\mathbb{R}^N} \nabla_A u \cdot \overline{\nabla_A v} \, dx 
& =  \operatorname{Re} \int_{\mathbb{R}^N} \nabla u \cdot \overline{\nabla v} \, dx
+ \operatorname{Re} \int_{\mathbb{R}^N} \nabla u \cdot \overline{i  A(x) v} \, dx \\
&\quad + \operatorname{Re} \int_{\mathbb{R}^N} i  A(x) u \cdot \overline{\nabla v} 
\, dx + \operatorname{Re} \int_{\mathbb{R}^N} i  A(x) u \cdot \overline{i  A(x) v} 
\, dx,
\end{align*}
by substituting $u(x) = u_{\theta}(\sqrt{\theta} x)$ 
(that is, $u_{\theta}(x) = u\left(\frac{x}{\sqrt{\theta}} \right)$) 
in  \eqref{rethetaequat} and by exploiting the change of variable 
$y= \sqrt{\theta} x$,  we obtain
\begin{equation} \label{Afvariable}
\begin{aligned}
&\operatorname{Re} \int_{\mathbb{R}^N}  \sqrt{\theta} \nabla_y u_{\theta}(y)
 \cdot       \overline{ \sqrt{\theta} \nabla_y v\big(\frac{y}{\sqrt{\theta}}\big)}
 \frac{1}{(\sqrt{\theta})^N} \, dy \\
&+  \operatorname{Re} \int_{\mathbb{R}^N}  \sqrt{\theta} \nabla_y u_{\theta}(y)
  \cdot \overline{i  A \big(\frac{y}{\sqrt{\theta}}\big) v
\big(\frac{y}{\sqrt{\theta}}\big)} \frac{1}{(\sqrt{\theta})^N} \, dy \\
&+  \sqrt{\theta} \operatorname{Re} \int_{\mathbb{R}^N}
 i  A \big(\frac{y}{\sqrt{\theta}}\big) u_{\theta}(y)
 \cdot \overline{\sqrt{\theta} \nabla_y v \big(\frac{y}{\sqrt{\theta}}\big)} 
 \frac{1}{(\sqrt{\theta})^N} \, dy  \\
&+  \operatorname{Re} \int_{\mathbb{R}^N}
 i  A \big(\frac{y}{\sqrt{\theta}}\big) u_{\theta}(y)
 \cdot   \overline{i  A \big(\frac{y}{\sqrt{\theta}}\big) v
 \big(\frac{y}{\sqrt{\theta}}\big)} \frac{1}{(\sqrt{\theta})^N} \, dy    \\
& = \theta \operatorname{Re} \int_{\mathbb{R}^N}
\big( f(| u_{\theta}(y)|^2) u_{\theta}(y) - u_{\theta}(y) \big)
\overline{v \big(\frac{y}{\sqrt{\theta}}\big)} \frac{1}{(\sqrt{\theta})^N}  \, dy.
\end{aligned}
\end{equation}
If we simplify $1/(\sqrt{\theta})^N$ and  put in evidence $\theta$, it follows that
\begin{align*}
&\operatorname{Re} \int_{\mathbb{R}^N}
\Big(\nabla_y + i \frac{1}{\sqrt{\theta}} A \big(\frac{y}{\sqrt{\theta}}\big) \Big)
 u_{\theta}(y) \cdot  \overline{\Big( \nabla_y + i \frac{1}{\sqrt{\theta}}
  A \big(\frac{y}{\sqrt{\theta}}\big) \Big) v
 \big(\frac{y}{\sqrt{\theta}}\big) } \, dy  \\
&= \operatorname{Re} \int_{\mathbb{R}^N}
\left( f(| u_{\theta}(y)|^2) u_{\theta}(y) - u_{\theta}(y) \right)
\overline{v \big(\frac{y}{\sqrt{\theta}}\big)}  \, dy.
\end{align*}
By assumption (A2), and since
$\frac{1}{\sqrt{\theta}} A \big(\frac{y}{\sqrt{\theta}}\big) = A(y)$
for every $y \in \mathbb{R}^N$, we have
\[
\operatorname{Re} \int_{\mathbb{R}^N} \nabla_{A} u_{\theta}(y) \cdot
 \overline{\nabla_A v \big(\frac{y}{\sqrt{\theta}}\big) } \, dy
 = \operatorname{Re} \int_{\mathbb{R}^N}
\left( f(| u_{\theta}(y)|^2) u_{\theta}(y) - u_{\theta}(y) \right)
 \overline{v \big(\frac{y}{\sqrt{\theta}}\big)}  \, dy;
\]
thus we can conclude $u_{\theta}$ satisfies \eqref{uthetaequat}.

\subsection*{Acknowledgments}
S. Barile was partially supported by the INdAM-GNAMPA Project
2017 ``Metodi variazionali per fenomeni non-locali''.
G. M. Figueiredo was supported by CNPQ, CAPES, FAP-DF.

\begin{thebibliography}{00}

\bibitem{alvesfig} C. O. Alves, G. M. Figueiredo; 
\emph{Multiple Solutions for a Semilinear Elliptic Equation with Critical 
Growth and Magnetic Field}, Milan J. Math., \textbf{82} (2) (2014), 389-405.

\bibitem{AlvesSM} C. O. Alves,  M. Montenegro, M. A. S. Souto; 
\emph{Existence of a ground state solution for a nonlinear scalar field 
equation with critical growth}, Calc. Var., \textbf{43} (2012), 537-554.

\bibitem{ariolisz} G. Arioli, A. Szulkin;
\emph{A semilinear     Schr\"odinger equations in the presence of a magnetic field},
  Arch. Ration. Mech. Anal., \textbf{170} (2003), 277--295.

\bibitem{Barcs} S. Barile, S. Cingolani, S. Secchi;
\emph{Single-peaks for a magnetic Schr\"odinger equation with critical growth}, 
Adv. Differential Equations, \textbf{11} (10) (2006), 1135-1166.

\bibitem{N2magnetic} S. Barile, G. M.  Figueiredo;
\emph{An existence result for Schr\"odinger equations with magnetic 
fields and exponential critical growth},
 J. Elliptic and Parabolic Equations, 3(2017), pp 105-125.

\bibitem{BL} H. Berestycki, P. L. Lions;
\emph{Nonlinear scalar field equations I and II},
 Arch. Rat. Mech. Anal. \textbf{82} (1983), 313-375.

\bibitem{Brezis} H. Brezis;
\emph{Functional Analysis, Sobolev Spaces and Partial Differential Equations},
 Springer, New York (2010).

\bibitem{chsz} J. Chabrowski, A. Szulkin;
\emph{On the Schr\"odinger equation involving a critical Sobolev exponent and 
magnetic field}, Top. Meth. Nonlinear Anal.  \textbf{25} (2005), 3--21.

\bibitem{Colglazmart} S. Coleman, V. Glazer, A. Martin; 
\emph{Action minima among solutions to a class of Euclidean scalar
 field equations}, Comm. Math. Phys., \textbf{58} (1978), 211-221.

\bibitem{dingliu} Y. Ding, X. Liu; 
\emph{Semiclassical solutions of Schr\"odinger equations with magnetic 
fields and critical nonlinearities},  Manuscripta Math.,
 \textbf{140} (1-2) (2013), 51-82.

\bibitem{EstLions}  M. Esteban, P. L. Lions;
\emph{Stationary solutions of nonlinear Schr\"odinger equations with an 
external magnetic  field}, PDE and Calculus of Variations, Vol. 
I, 401--449, Progr. Nonlinear Differential Equations Appl. 1, 
Birkh\"{a}user Boston, MA, 1989.

\bibitem{pigong} P. Han; 
\emph{Solutions for singular critical growth Schr\"odinger equations with 
magnetic field},  Port. Math. (N.S.), \textbf{63} (2006), 37-45.

\bibitem{JJTan} A. Jeanjean, K. Tanaka; 
\emph{A remark on least energy solutions in $\mathbb{R}^N$}, 
Proc. Amer. Math. Soc., \textbf{131} (2002),  2399-2408.

\bibitem{liangsongBVP}  S. Liang, Y. Song; 
\emph{Multiplicity of solutions of perturbed Schr\"odinger equation with
 electromagnetic fields and critical nonlinearity in $\mathbb{R}^N$}, 
Bound. Value Probl. \textbf{240} (2014), 1-14.

\bibitem{liangzhang}  S. Liang, J.  Zhang; 
\emph{Solutions of perturbed Schr\"odinger equations with electromagnetic 
fields and critical nonlinearity}, Proc. Edinb. Math. Soc.,
  \textbf{54} (1) (2011), 131-147.

\bibitem{LL} E. H. Lieb, M. Loss; 
\emph{Analysis}, Graduate Studies in Mathematics, Vol. 14, Amer. Math. Soc., 1997.

\bibitem{lionsI}  P. L. Lions; 
\emph{The concentration-compactness principle in the calculus of variations. 
The locally compact case. Part I.},  Ann. Inst. Henri Poincar\'e, Anal. 
Non Lin\'eaire, \textbf{1} (1984),  109-145.

\bibitem{Lionsibero}  P. L. Lions; 
\emph{The concentration-compactness principle in the calculus of variations. 
The limit case, part 1}, Rev. Mat. Iberoamericana \textbf{1} (1985), 145-201.

\bibitem{Lionsibero2}  P. L. Lions; 
\emph{The concentration-compactness principle in the calculus of variations. 
The limit case, part 2}, Rev. Mat. Iberoamericana, \textbf{1} (1985), 45-121.

\bibitem{Schectzou} M. Schechter, W. Zou; 
\emph{Critical point theory and its applications}, Springer, New York 2006.

\bibitem{wangF} F. Wang; 
\emph{On an electromagnetic Schr\"odinger equation with critical growth}, 
 Nonlinear Anal., \textbf{69} (2008), 40884098.

\bibitem{Willem} M. Willem; 
\emph{Minimax Theorems}, Birkh\"{a}user, Boston (1996).

\bibitem{ZZou} J. Zhang, W. Zou; 
\emph{A Berestycki-Lions theorem revisited}, 
Comm. Contemp. Math., \textbf{14} (5) (2012), 1-14.

\bibitem{tangwang} Z. M. Tang, Y. L. Wang; 
\emph{Least energy solutions for semilinear Schr\"odinger equation with 
electromagnetic fields and  critical growth},
Sci. China Math. \textbf{58} (2015), 2317-2328.

\end{thebibliography}

\end{document}
