Expansion in higher harmonics of boson stars using a generalized Ruffini-Bonazzola approach. Part 1. Bound states

Abstract

The method pioneered by Ruffini and Bonazzola (RB) to describe boson stars involves an expansion of the boson field which is linear in creation and annihilation operators. In the nonrelativistic limit, the equation of motion of RB is equivalent to the nonlinear Schrödinger equation. Further, the RB expansion constitutes an exact solution to a non-interacting field theory, and has been used as a reasonable ansatz for an interacting one. In this work, we show how one can go beyond the RB ansatz towards an exact solution of the interacting operator Klein-Gordon equation, which can be solved iteratively to ever higher precision. Our Generalized Ruffini-Bonazzola approach takes into account contributions from nontrivial harmonic dependence of the wavefunction, using a sum of terms with energy k E0, where k⩾1 and E0 is the chemical potential of a single bound axion. The method critically depends on an expansion in a parameter Δ ≡ √1 − E02/m2 < 1, where m is the mass of the boson. In the case of the axion potential, we calculate corrections which are relevant for axion stars in the transition or dense branches of solutions. We find with high precision the local minimum of the mass, Mmin≈ 463 f2/m, at Δ≈0.27, where f is the axion decay constant. This point marks the crossover from the transition branch to the dense branch of solutions, and a corresponding crossover from structural instability to stability.