Radial oscillations of boson stars made of ultralight repulsive dark matter

We compute the lowest frequency radial oscillation modes of boson stars. It is assumed that the object is made of pseudo-Goldstone bosons subjected to a scalar potential that leads to a repulsive self-interaction force, and which is characterized by two unknown mass scales m (mass of the particle) and F (decay constant). First we integrate the Tolman-Oppenheimer-Volkoff equations for the hydrostatic equilibrium of the star, and then we solve the Sturm-Liouville boundary value problem for the perturbations using the shooting method. The effective potential that enters into the Schrödinger-like equation as well as several associated eigenfunctions are shown as well. Moreover, we found that the large frequency separation, i.e. the difference between consecutive modes, is proportional to the square root of the mass of the star and the cube of the mass scale defined by Λ ≡√{ mF }.