Self-consistent construction of virialized wave dark matter halos
Using C++ and MPI, we solve the Schrödinger-Poisson equation numerically by treating the wavefunction as composed of noninteracting eigenstates and assigning their amplitudes based on a distribution function, assuming spherically symmetric. This method provides a way to construct Milky Way-sized wave dark matter halos.