A quadratic function is a better model with a BIC of 56. The period derivative obtained from a quadratic function is −6.2 × 10−9 ± 2.9 ×10−10 per orbit, indicating an orbital decay timescale 1.4 Myr.