Phase-mixed comet surface density (${\Sigma }_c(r)$)

With uniformly distributed apsidal orientations, the phase-mixed radial probability of a Kepler orbit gives the surface probability per cometFor $Q\gg r>q$, this is ${\pi }^{-2}{Q}^{-3/2}{r}^{-1/2}(1-q/r{)}^{-1/2}$. Multiply by the population size for a number surface density of a disk.