We adapted the Multinest code to execute the dynamical fits as well (e.g., Nesvorný et al. 2013).

[**From David K’s email on 1/15/2012:**]{}

Cool, with such a complex and high SNR signal, I wonder if this could be the first example of a perturber being uniquely determined from TTV alone? Exciting stuff!

I am running refined TTV fits, they will take 1-2 days minimum.

The next day David K sent me KOI-872’s TTVs from the [*Kepler*]{} quarters 1-6, and I performed a more complete dynamical analysis immediately. Each dynamical fit was executed by the analytic codes in minutes (Figure 2). All the effort I have put into the code over the past three years finally paid off. Doing this with N-body is possible but takes days on a supercomputer.

[**From David N’s email on 1/16/2012:**]{}

It is looking good!

The analytic search routine in 5D (planet’s mass,a,e,varpi,capm) found promising solutions that fit TTVs very nicely (TDVs were ignored at this step). Most of the solutions could not have been fine tuned with my numerical (and more) precise code in 7D (mass,a,e,inc,capom,varpi,capm), so I discarded them. Two remained.

Solution 1: mass$=$0.000934 Msun, sema$=$0.2998 AU, ecc$=$0.0560

This solution is related to the 