We consider normalised solutions of the stationary Gross–Pitaevskii–Poisson (GPP) equation with a defocusing local nonlinear term, \(\begin{aligned} -\Delta u+\lambda u+|u|^2u =(I_\alpha *|u|^2)u\quad \text { in } \mathbb {R}^3,\qquad \int _{\mathbb {R}^3}u^2dx=\rho ^2, \end{aligned}\) where \(\rho ^2>0\) is the prescribed mass of the solutions, \(\lambda \in \mathbb {R}\) is an a-priori unknown Lagrange multiplier, and \(I_\alpha (x)=A_\alpha |x|^{3-\alpha }\) is the Riesz potential of order \(\alpha \in (0,3)\) . When \(\alpha =2\) this problem appears in the models of self–gravitating Bose–Einstein condensates, which were proposed in cosmology and astrophysics to describe Cold Dark Matter and Boson Stars.
We establish the existence of branches of normalised solutions to the GPP equation, paying special attention to the shape of the associated mass–energy relation curves and to the limit profiles of solutions at the endpoints of these curves. The behaviour of normalised solutions depends sensitively on whether \(\alpha \) is greater than, equal to, or less than one. The main novelty in this work is in the derivation of sharp asymptotic estimates on the mass–energy curves. These estimates allow us to show that after appropriate rescalings, the constructed normalised solutions converge either to a ground state of the Choquard equation or to a compactly supported radial ground state of the integral Thomas–Fermi equation. A major difficulty is the case \(\alpha <1\) when analysis of limit profiles involves estimates along the branch of unstable mountain pass solutions.
