Thomas Donlon II, Sukanya Chakrabarti, Lawrence M. Widrow, Sophia Vanderwaal, Scott Ransom, Enrico Ramirez-Ruiz

We present a novel method that enables us to estimate the acceleration of individual millisecond pulsars (MSPs) using only their spin period and its time derivative. For our binary MSP sample, we show that one can obtain an empirical calibration of the magnetic braking term that relies only on observed quantities. We find that such a model for magnetic braking is only valid for MSPs with small surface magnetic field strengths ($<3\times10^8$ G) and large characteristic ages ($>$ 5 Gyr). With this method we are able to effectively double the number of pulsars with line-of-sight acceleration measurements, from 28 to 54 sources. This expanded dataset leads to an updated measurement of the total density in the midplane, which we find to be $\rho_0$ = 0.086 $\pm$ 0.001 stat. $\pm$ 0.006 sys. M$_\odot$/pc$^3$, and an updated measurement of the local dark matter density, which we calculate to be $\rho_{0,\mathrm{DM}}$ = 0.014 $\pm$ 0.005 stat. $\pm$ 0.006 sys. M$_\odot$/pc$^3$ (0.53 $\pm$ 0.30 GeV/cm$^3$). This updated value for $\rho_{0,\mathrm{DM}}$ is in good agreement with literature values derived from kinematic estimates. We show that each new acceleration measurement improves the precision of $\rho_{0,\mathrm{DM}}$ by the same amount as roughly 10$^5$ stars. The pulsar accelerations are very asymmetric above and below the disk; we show that the shape and size of this asymmetry can be largely explained by the north-south asymmetry of disk star counts and the offset in the Milky Way disk and halo centers of mass due to the Large Magellanic Cloud.