We study a nonlinear Schrödinger–Poisson system which reduces to the nonlinear and nonlocal PDE -Δu+u+λ21ω|x|N-2⋆ρu2ρ(x)u=|u|q-1ux∈RN,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} - \Delta u+ u + \lambda ^2 \left( \frac{1}{\omega |x|^{N-2}}\star \rho u^2\right) \rho (x) u = |u|^{q-1} u \quad x \in {{\mathbb {R}}}^N, \end{aligned}$$\end{document}where ω=(N-2)|SN-1|,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\omega = (N-2)|{\mathbb {S}}^{N-1} |,$$\end{document}λ>0,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\lambda >0,$$\end{document}q∈(1,2∗-1),\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$q\in (1,2^{*} -1),$$\end{document}ρ:RN→R\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\rho :{{\mathbb {R}}}^N \rightarrow {{\mathbb {R}}}$$\end{document} is nonnegative, locally bounded, and possibly non-radial, N=3,4,5\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$N=3,4,5$$\end{document} and 2∗=2N/(N-2)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2^*=2N/(N-2)$$\end{document} is the critical Sobolev exponent. In our setting ρ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\rho $$\end{document} is allowed as particular scenarios, to either (1) vanish on a region and be finite at infinity, or (2) be large at infinity. We find least energy solutions in both cases, studying the vanishing case by means of a priori integral bounds on the Palais–Smale sequences and highlighting the role of certain positive universal constants for these bounds to hold. Within the Ljusternik–Schnirelman theory we show the existence of infinitely many distinct pairs of high energy solutions, having a min–max characterisation given by means of the Krasnoselskii genus. Our results cover a range of cases where major loss of compactness phenomena may occur, due to the possible unboundedness of the Palais–Smale sequences, and to the action of the group of translations.
