In this paper, we implement the adaptive moving mesh method (AMM) to the solution of initial value problems involving the Schrödinger equation, and more specifically the Schrödinger-Poisson system of equations. This method is based on the solution of the problem on a discrete domain, whose resolution is coordinate and time-dependent, and allows us to dynamically assign numerical resolution in terms of desired refinement criteria. We apply the method to solve various test problems involving stationary solutions of the SP system, and toy scenarios related to the disruption of subhalo s made of ultralight bosonic dark matter traveling on top of host galaxies.

• Received 14 December 2021
• Accepted 14 March 2022

DOI:https://doi.org/10.1103/PhysRevD.105.083521