TTVs are derived from Q1-Q17 Kepler data. x-axes: “Observed Tc” (Mid-Transit Time): EXOFAST’s best-fits from Normalized PDCSAP_FLUX Kepler light flux vs. time (BJD_tdb - 2454900) data. y-axes: “(O – C)”: difference between Observed Tc and the Calculated Tc from the graphically obtained linear ephemeris.
Kepler KOI-4571 (KIC-6761777) 4-(or more?)-Planet System
Discussion: In the example of 4571.01, after a sinusoidal curve-fit of the (O-C) vs. Time data showed (Figure 1. below) a periodicity (P_ttv) of 198.93 ± 3.77 days (198.46 days was observed in the Lomb-Scargle Periodogram (LSP) of the same data), a plot (Figure 2.) of the Residuals also gave a good sinusoidal curve-fit with a periodicity of 494.54 ± 26.45 days (LSP: 492.43 days). Similarly, the Residuals-of-the-Residuals also gave a good sinusoidal curve-fit (Figure 3.) with a periodicity of 1309.05 ± 318.95 days (LSP: 1360.26 days). The added combination of all three sinusoidal curves, arrayed in Figure 4. produces a complex overall curvature reasonably consistent with the initial data. While it is certainly possible (see recent work of Lithwick and others) that some of this unusual curvature obtains from eccentric orbits (and maybe even precessing eccentric orbits) of planetary objects in this system, it is also possible that at least 4 planets (with only 4571.01 transiting) in near-circular orbits are mutually-interacting to give the TTV distribution observed.