2025/07/03 by Jiao, Cheng-Liang, Zhao, Er-gang, Zhu, Liying +1
#Astrophysics of Galaxies (astro-ph.GA) #FOS: Physical sciences #High Energy Astrophysical Phenomena (astro-ph.HE) #Solar and Stellar Astrophysics (astro-ph.SR)
paper · doi:10.48550/arxiv.2507.02621
Spherically symmetric accretion incorporating self-gravity constitutes a three-point boundary value problem (TPBVP) governed by constraints at the outer boundary, sonic point, and accretor surface. Previous studies have two limitations: either employing an incorrect formula for self-gravity potential in analytical treatments, or introducing additional input parameters in numerical implementations to circumvent solving the full TPBVP. To address these issues, we present a self-consistent TPBVP formulation, solved using the relaxation method. We also derive approximate analytical formulae that enable rapid estimates of self-gravity effects. Our analysis identifies a dimensionless parameter β≡ 2G ρ rout2/aout2 that characterizes the strength of self-gravity, where ρ and rout are the mean density and outer radius of the flow, respectively, and aout is the adiabatic sound speed of the external medium. For practical estimation, ρ may be approximated by the external medium density ρout. We identify an upper limit for β, beyond which steady accretion becomes unsustainable -- a behavior consistent with classical gravitational instability that previous studies failed to capture. The accretion rate enhancement decreases monotonically as the adiabatic index γ increases. For γ=5/3, self-gravity ceases to augment the accretion rate. These theoretical predictions are validated by our numerical solutions. We further apply our results to two astrophysical scenarios: hyper-Eddington accretion onto supermassive black hole seeds in the early Universe, where self-gravity is significant; and accretion onto stellar-mass objects embedded in active galactic nuclei (AGN) disks, where self-gravity is non-negligible under certain conditions and should be evaluated using β.