In this article, we study the existence and multiplicity of high energy solutions to the problem proposed as a model for the dynamics of galaxies: −Δu+V(x)u=|u|2∗−2u|y|,x=(y,z)∈Rm×Rn−m,where n>4, 2≤m<n, 2∗≔2(n−1)n−2 and potential function V(x):Rn→R. Benefiting from a global compactness result, we show that there exist at least two positive high energy solutions. Our proofs are based on barycenter function, quantitative deformation lemma and Brouwer degree theory.