Following Ruffini and Bonazzola, we use a quantized boson field to describe condensates of axions forming compact objects. Without substantial modifications, the method can only be applied to axions with decay constant, f a , satisfying δ = (f a / M P ) 2 1, where M P is the Planck mass. Similarly, the applicability of the Ruffini-Bonazzola method to axion stars also requires that the relative binding energy of axions satisfies1, where E a and m a are the energy and mass of the axion. The simultaneous expansion of the equations of motion in δ and ∆ leads to a simplified set of equations, depending only on the parameter, λ = √ δ / ∆ in leading order of the expansions. Keeping leading order in ∆ is equivalent to the infrared limit, in which only relevant and marginal terms contribute to the equations of motion. The number of axions in the star is uniquely determined by λ. Numerical solutions are found in a wide range of λ. At small λ the mass and radius of the axion star rise linearly with λ. While at larger λ the radius of the star continues to rise, the mass of the star, M , attains a maximum at λ max 0.58. All stars are unstable for λ > λ max . We discuss the relationship of our results to current observational constraints on dark matter and the phenomenology of Fast Radio Bursts.