Merger of a white dwarf-neutron star binary to 10²⁹ carat diamonds: Origin of the pulsar planets

We show that the merger and tidal disruption of a carbon/oxygen (C/O) white dwarf (WD) by a neutron star (NS) binary companion provides a natural formation scenario for the PSR B1257+12 planetary system. Starting with initial conditions for the debris disc produced of the disrupted WD, we model its long-term viscous evolution, including for the first time the effects of mass and angular momentum loss during the early radiatively inefficient accretion flow (RIAF) phase and accounting for the unusual C/O composition on the disc opacity. For plausible values of the disc viscosity α ~ 10^-3-10^-2 and the RIAF mass-loss efficiency, we find that the disc mass remaining near the planet formation radius at the time of solid condensation is sufficient to explain the pulsar planets. Rapid rocky planet formation via gravitational instability of the solid carbon dominated disc is facilitated by the suppression of vertical shear instabilities due to the high solid-to-gas ratio. Additional evidence supporting a WD-NS merger scenario includes (1) the low observed occurrence rate of pulsar planets (≲1 per cent of NS birth), comparable to the expected WD-NS merger rate; (2) accretion by the NS during the RIAF phase is sufficient to spin PSR B1257+12 up to its observed 6 ms period; (3) similar models of 'low angular momentum' discs, such as those produced from supernova fallback, find insufficient mass reaching the planet formation radius. The unusually high space velocity of PSR B1257+12 of ≳326 km s^-1 suggests a possible connection to the calcium-rich transients, dim supernovae which occur in the outskirts of their host galaxies and were proposed to result from mergers of WD-NS binaries receiving supernova kicks. The C/O disc composition implied by our model likely results in carbon-rich planets with diamond interiors.