Abstract We consider normalised solutions of the stationary Gross–Pitaevskii–Poisson (GPP) equation with a defocusing local nonlinear term, aligned - u+ u+|u|²u = (I_ *|u|²) u in R³, ₑ³u²dx= ², aligned - Δ u + λ u + | u | 2 u = (I α ∗ | u | 2) u in R 3, ∫ R 3 u 2 d x = ρ 2, where ²>0 ρ 2 > 0 is the prescribed mass of the solutions, R λ ∈ R is an a-priori unknown Lagrange multiplier, and I_ (x) =A_ |x|^3- I α (x) = A α | x | 3 - α is the Riesz potential of order (0, 3) α ∈ (0, 3). When =2 α = 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
Molle et al. (Sat,) studied this question.