Line 309 of particle_info_ops.f90 in the BondiAccretionRate routine has a bug where it calculates the relative velocity of the sink particle and gas in the sink's host cell. The array particle%imom(:) was divided by the particle mass, however imom is not the particle momentum (bad naming convention here), it is actually the particle's velocity vector. So I removed the division by the particle mass.