We present the construction of the ground state of the Gross–Pitaevskii–Poisson equations using genetic algorithms. By employing numerical solutions, we develop an empirical formula for the density that works within the considered parameter space. Through the analysis of both numerical and empirical solutions, we investigate the stability of these ground-state solutions. Our findings reveal that while the numerical solution outperforms the empirical formula, both solutions lead to similar oscillation modes. We observe that the stability of the solutions depends on specific values of the central density and the nonlinear self-interaction term and establish an empirical criterion delineating the conditions under which the solutions exhibit stability or instability.
