The initial data is obtained by solving an eigenvalue problem very similar to the one for boson stars, with similar assumptions and boundary conditions. The solutions also share similar features and display a similar behavior. A nice property of the Newtonian limit is that all the solutions can be obtained by rescaling from one known solution ,
The possibility of including self-interaction terms in the potential was considered in , studying also the gravitational cooling (i.e., the relaxation and virialization through the emission of scalar field bursts) of spherical perturbations. Non-spherical perturbations were further studied in , showing that the final state is a spherically symmetric configuration. Single Newtonian boson stars were studied in , either when they are boosted with/without an external central potential. Rotating stars were first successfully constructed in  within the Newtonian approach. Numerical evolutions of binary boson stars in Newtonian gravity are discussed in Section 4.2.
Recent work by Chavanis with Newtonian gravity solves the Gross–Pitaevskii equation, a variant of Eq. (58) which involves a pseudo-potential for a Bose–Einstein condensate, to model either dark matter or compact alternatives to neutron stars [46, 45, 47].
Much recent work considers boson stars from a quantum perspective as a Bose-Einstein condensate involving some number, , of scalar fields. Ref.  studies the collapse of boson stars mathematically in the mean field limit in which . Ref.  argues for the existence of bosonic atoms instead of stars. Ref.  uses numerical methods to study the mean field dynamics of BSs.
Living Rev. Relativity 15, (2012), 6
This work is licensed under a Creative Commons License.