Multiple normalized solutions for a Sobolev critical Schrödinger-Poisson-Slater equation

We look for solutions to the Schr\"{o}dinger-Poisson-Slater equation $$- \Delta u + \lambda u - \gamma (|x|^{-1} * |u|^2) u - a |u|^{p-2}u = 0 \quad \text{in} \quad \mathbb{R}^3, $$ which satisfy \begin{equation*} \int_{\mathbb{R}^3}|u|^2 \, dx = c \end{equation*} for some prescribed $c>0$. Here $ u \in H^1(\mathbb{R}^3)$, $\gamma \in \mathbb{R},$ $ a \in \mathbb{R}$ and $p \in (\frac{10}{3}, 6]$. When $\gamma >0$ and $a > 0$, both in the Sobolev subcritical case $p \in (\frac{10}{3}, 6)$ and in the Sobolev critical case $p=6$, we show that there exists a $c_1>0$ such that, for any $c \in (0,c_1)$, the equation admits two solutions $u_c^+$ and $u_c^-$ which can be characterized respectively…