Distant Companions and Planets around Millisecond Pulsars
  We present a general method for determining the masses and orbital parameters of binary millisecond pulsars with long orbital periods (P_orb >> 1 yr), using timing data in the form of pulse frequency derivatives. We apply our method to analyze the properties of the second companion in the PSR B1620-26 triple system. We use the latest timing data for this system to constrain the mass and orbital parameters of the second companion. We find that all possible solutions have a mass m_2 in the range 2.4 10^-4 M_sun <= m_2 sin i_2 <= 1.2 10^-2 M_sun, i.e., almost certainly excluding a second companion of stellar mass and suggesting instead that the system contains a planet or a brown dwarf. Using Monte-Carlo realizations of the triple configuration in three dimensions we find the most probable value of m_2 to be 0.010(5) M_sun, corresponding to a distance of 38(6) AU from the center of mass of the inner binary (the errors indicate 80% confidence intervals). We also apply our method to analyze the planetary system around PSR B1257+12, where a distant, giant planet may be present in addition to the three well-established Earth-mass planets. We find that the simplest interpretation of the frequency derivatives implies the presence of a fourth planet with a mass of ~100 M_earth in a circular orbit of radius ~40 AU. 
K. J. Joshi and F. A. Rasio