Do pulsar radio fluxes violate the inverse-square law?

Singleton et al. (arXiv:0912.0350, 2009) have argued that the flux of pulsars measured at 1400 MHz shows an apparent violation of the inverse-square law with distance (r), and instead the flux scales as 1/r. They deduced this from the fact that the convergence error obtained in reconstructing the luminosity function of pulsars using an iterative maximum likelihood procedure is about 105 times larger for a distance exponent of two (corresponding to the inverse-square law) compared to an exponent of one. When we applied the same technique to this pulsar dataset with two different values for the trial luminosity function in the zeroth iteration, we find that neither of them can reproduce a value of 105 for the ratio of the convergence error between these distance exponents. We then reconstruct the differential pulsar luminosity function using Lynden-Bell's C-method after positing both inverse-linear and inverse-square scalings with distance. We show that this method cannot help in discerning between the two exponents. Finally, when we tried to estimate the power-law exponent with a Bayesian regression procedure, we do not get a best-fit value of one for the distance exponent. The model residuals obtained from our fitting procedure are larger for the inverse-linear law compared to the inverse-square law. Moreover, the observed pulsar flux cannot be parameterized only by power-law functions of distance, period, and period derivative. Therefore, we conclude from our analysis using multiple methods that there is no evidence that the pulsar radio flux at 1400 MHz violates the inverse-square law or that the flux scales inversely with distance.