# Towards an Analytic Construction of the Wavefunction of Boson Stars

## TL;DR

This paper introduces a semi-analytic, highly accurate method for constructing the wavefunction profile of non-interacting boson stars using the Gross-Pitaevskii-Poisson equations in the Newtonian limit.

## Contribution

It provides a novel semi-analytic solution to the equations governing boson stars, improving stability and accuracy over previous approaches.

## Key findings

- Solution is stable and numerically precise with better than 10^{-6} accuracy.
- Applicable to non-interacting scalar fields forming boson stars.
- Enhances understanding of boson star profiles in the Newtonian regime.

## Abstract

Light scalar fields can form gravitationally bound compact objects called boson stars. The profile of boson stars in the Newtonian limit is described by the Gross-Pitaevskii-Poisson equations. We present a semi-analytic solution to these equations and construct the profile of boson stars formed by a non-interacting scalar field. Our solution is stable with respect to numerical errors and has accuracy better than $10^{-6}$ over the entire range.