Scaling relations and tidal disruption in spin $s$ ultralight dark matter models

Abstract: We explore the impact of spin 0, spin 1 and spin 2 Ultra-Light Dark Matter (ULDM) on small scales by numerically solving the Schr\"odinger-Poisson system using the time-split method. We perform simulations of ULDM for each spin, starting with different numbers of identical initial solitons and analyse the properties of the resulting halos after they merge and relax in a steady-state. Our findings reveal that higher spin values lead to broader, less dense final halo with more prominent Navarro-Frenk-White (NFW) tails, a characteristic that persists regardless of the number of initial solitons involved. We identify scaling relations that describe the density profile, core and NFW tail of spin $s$ ULDM halos as a function of the number of initial solitons $N_{sol}$. These relations allow us to construct equivalent halos based on average density or total mass, for arbitrarily large $N_{sol}$, without having to simulate those systems. We simulate the orbit of a ULDM satellite in a constructed halo treated as an external potential, and find that for host halos having the same average density, the orbital decay time of the satellite is as predicted for uniform sphere host halo regardless of the spin. However, satellites orbiting haloes having the same mass for each spin, result in faster disruption in the case of spin 0, while satellites orbiting haloes having the same core size result in faster disruption in the case of spin 2.