Precision X-ray Timing of RX J0806.3+1527 with CHANDRA: Evidence for Gravitational Radiation from an Ultracompact Binary

RX J0806.3+1527 is a candidate double degenerate binary with possibly the shortest known orbital period. The source shows an $\approx 100 \%$ X-ray intensity modulation at the putative orbital frequency of 3.11 mHz (321.5 s). If the system is a detached, ultracompact binary gravitational radiation should drive spin-up with a magnitude of $\dot\nu \sim 10(exp -16)$ Hi s$(exp -l}$. Efforts to constrain the X-ray frequency evolution to date have met with mixed success, principally due to the sparseness of earlier observations. Here we describe the results of the first phase coherent X-ray monitoring campaign on RX J0806.3+1527 with {\it Chandra) . We obtained a total of 70 ksec of exposure in 6 epochs logarithmically spaced over 320 days. With these data we conclusively show that the X-ray frequency is increasing at a rate of $3.77 \pm\ 0.8 \times l0(exp -16)$ Hz s$(exp -l)$. Using the ephemeris derived from the new data we are able to phase up all the earlier {\it Chandra) and ROSAT data and show they are consistent with a constant $\dot\nu = 3.63 \pm 0.06 \times 10(exp -16)$ Hz s$(exp -1)s over the past decade. This value appears consistent with that recently derived by Israel et al. largely from monitoring of the optical modulation, and is in rough agreement with the solutions reported initially by Hakala et al., based on ground-based optical observations. The large and stable $\dot\nu$ over a decade is consistent with gravitational radiation losses driving the evolution. An intermediate polar (IP) scenario where the observed X-ray period is the spin period of an accreting white dwarf appears less tenable because the observed $\dot\nu$ requires an $\dot m \approx 2 \times 10(exp -8)$ $M-{\odot)$ yr$(exp -l)$, that is much larger than that inferred from the observed X-ray luminosity (although this depends on the uncertain distance and bolometric corrections), and it is difficult to drive such a high $\dot m$ in a binary system with parameters consistent with all the multiwavelength data. If the ultracompact scenario is