KOI-134 c: Resonant, Inclined Exoplanet

• KOI-134 c is a non-transiting, sub-Saturn exoplanet in a 2:1 resonance with KOI-134 b, noted for its significant mutual inclination and dynamic orbital interaction.
• Its orbital parameters—such as a ~33.95-day period, moderate eccentricity, and precise semimajor axis measurements—are robustly determined using advanced TTV and TDV modeling.
• The complex architecture, combining resonant capture with sizeable mutual inclination, challenges standard formation theories and offers a stringent test for planetary evolution models.

KOI-134 c is a dynamically significant, non-transiting exoplanet in the KOI-134 system, interior to and in 2:1 mean-motion resonance with the transiting warm Jupiter KOI-134 b. Its presence and detailed orbital architecture are inferred exclusively through the pronounced transit timing variations (TTVs) and transit duration variations (TDVs) recorded in KOI-134 b's Kepler photometry. Unlike the majority of Kepler multi-planet systems, which are tightly coplanar, KOI-134 b and c exhibit a substantial mutual inclination, presenting a rare instance of a resonant, significantly non-coplanar planetary configuration. The architecture, mass, and orbital parameters of KOI-134 c provide a stringent test for dynamical theories of planet formation and evolution, as no single canonical mechanism can simultaneously account for the observed resonance, mutual inclination, moderate eccentricities, and long-term stability (Nabbie et al., 3 Jul 2025).

1. Orbital Parameters and Dynamical State

The best-fit $N$-body solution places KOI-134 c with a period $P_c = 33.950^{+0.013}_{-0.020}$ days, interior to KOI-134 b's $P_b = 67.1277^{+0.0045}_{-0.0057}$ days, yielding a period ratio $P_b/P_c \approx 1.9772$, narrowly inside exact 2:1 commensurability. The normalized distance to resonance is $\Delta = -0.011$, indicating the pair is "well within the resonant domain" as defined by first-order resonance conventions. The orbital semimajor axis for c is $a_c = 0.2305^{+0.0069}_{-0.0081}$ AU (non-photodynamic solution) and $a_c = 0.22534^{+0.00052}_{-0.00050}$ AU (photodynamical check), both closely consistent. The measured eccentricity is $e_c = 0.24^{+0.12}_{-0.03}$, moderate but robustly nonzero, with high-precision integrations adopting $e_c = 0.285108$.

Table 1: Adopted KOI-134 c Orbital Elements | Parameter | Value (Best-fit) | Unit | |------------------------|-----------------------------------|---------| | Period ($P_c$) | $P_c = 33.950^{+0.013}_{-0.020}$0 | days | | Semimajor axis ($P_c = 33.950^{+0.013}_{-0.020}$1) | $P_c = 33.950^{+0.013}_{-0.020}$2 | AU | | Eccentricity ($P_c = 33.950^{+0.013}_{-0.020}$3) | $P_c = 33.950^{+0.013}_{-0.020}$4 | – | | Argument of periastron ($P_c = 33.950^{+0.013}_{-0.020}$5) | $P_c = 33.950^{+0.013}_{-0.020}$6 | deg | | Longitude of node ($P_c = 33.950^{+0.013}_{-0.020}$7) | $P_c = 33.950^{+0.013}_{-0.020}$8 | deg | | Inclination ($P_c = 33.950^{+0.013}_{-0.020}$9) | $P_b = 67.1277^{+0.0045}_{-0.0057}$0 | deg | | Mean anomaly ($P_b = 67.1277^{+0.0045}_{-0.0057}$1) | $P_b = 67.1277^{+0.0045}_{-0.0057}$2 | deg |

KOI-134 c’s orbital parameters are inferred from the dynamical imprint on KOI-134 b’s transits, as c does not itself transit.

2. Dynamical Mass Determination via TTV–TDV Modeling

KOI-134 c’s mass is robustly measured as $P_b = 67.1277^{+0.0045}_{-0.0057}$3 from a suite of $P_b = 67.1277^{+0.0045}_{-0.0057}$4-body fits matching both the timing and duration variations observed in KOI-134 b. An independent photodynamical analysis obtained $P_b = 67.1277^{+0.0045}_{-0.0057}$5, confirming a sub-Saturn mass. Radial velocity data (TRES) are dominated by stellar jitter and only exclude stellar-mass companions.

KOI-134 b exhibits TTVs of $P_b = 67.1277^{+0.0045}_{-0.0057}$6 hours—among the largest in Kepler multiplanet systems—and TDVs at the $P_b = 67.1277^{+0.0045}_{-0.0057}$7 minute level, inconsistent with a constant model at $P_b = 67.1277^{+0.0045}_{-0.0057}$8. The mass inference for c derives from forward-model $P_b = 67.1277^{+0.0045}_{-0.0057}$9-body integrations (REBOUND/IAS15), optimized to match both the observed mid-transit times and durations. The best-fit solution is found only when KOI-134 c and b are locked in a 2:1 resonance with a modest mutual inclination. An alternate, higher mass/higher inclination solution is rejected on dynamical stability grounds.

The central mathematical relationship for mass inference is operational rather than analytic: the full $P_b/P_c \approx 1.9772$0-body integration simulates the coupled TTV-TDV signatures induced by c’s gravitational perturbations. The architecture is especially sensitive to the values of the resonant angles: $P_b/P_c \approx 1.9772$1 where $P_b/P_c \approx 1.9772$2 and $P_b/P_c \approx 1.9772$3. The mutual inclination, reported as $P_b/P_c \approx 1.9772$4, is calculated using the standard formula derived from sky-plane inclinations and nodal angles.

3. Mutual Inclination, TTV/TDV Degeneracy Breaking, and Inclination Variability

KOI-134 c is not coplanar with KOI-134 b. The joint TTV–TDV fit requires a mutual inclination of $P_b/P_c \approx 1.9772$5, placing KOI-134 well outside the canonical $P_b/P_c \approx 1.9772$6–$P_b/P_c \approx 1.9772$7 mutual inclinations typical of compact Kepler multiplanet systems. The necessity for mutual inclination arises from the strong TDVs in KOI-134 b, reflecting secular changes in the observed impact parameter associated with nodal precession induced by the inclined companion.

The modeling pipeline entailed extraction of individual transit times and durations (using the batman light-curve model), subtraction of a linear ephemeris to yield TTVs, exploratory $P_b/P_c \approx 1.9772$8-body integrations of parameterized two-planet configurations, numerical transit search and duration computation at each simulated transit, and joint likelihood evaluation against the TTV and TDV data. Only models with substantial mutual inclination—combined with resonance—match both TTV and TDV datasets.

Over long timescales, the mutual inclination varies significantly: 10 Myr integrations show $P_b/P_c \approx 1.9772$9 oscillating between $\Delta = -0.011$0 and $\Delta = -0.011$1, with associated secular precession periods of $\Delta = -0.011$2 years.

4. Resonant Dynamics and Long-Term Stability

KOI-134 c and b are in a strong 2:1 mean-motion resonance as established by both the proximity in period ratio and the behavior of the resonant angles. Specifically, the eccentricity-type resonant angle $\Delta = -0.011$3 librates about $\Delta = -0.011$4 with an amplitude of $\Delta = -0.011$5 and a period of $\Delta = -0.011$6 years, while $\Delta = -0.011$7 does not librate. This marks an asymmetric resonant configuration distinct from many classic cases.

The system is also notable for significant "free" and "forced" eccentricity components. For c: $\Delta = -0.011$8 and $\Delta = -0.011$9, with the corresponding values for b being $a_c = 0.2305^{+0.0069}_{-0.0081}$0 and $a_c = 0.2305^{+0.0069}_{-0.0081}$1. The non-negligible free eccentricity further distinguishes the dynamical state from pure damped-migration outcomes.

Notably, secular inclination evolution is so strong that KOI-134 b's transit chord will move off the stellar disk, making b unobservable as a transiting planet in approximately 100 years (predicted disappearance of transits around the year 2059). The fraction of time KOI-134 b is observable in transit from Earth is only $a_c = 0.2305^{+0.0069}_{-0.0081}$2, while c is predicted never to transit due to its sky-plane inclination.

5. Constraints on Formation and Evolutionary Pathways

The combined properties of KOI-134 c—a resonance, substantial mutual inclination, and moderate eccentricity—pose a challenge to single-process formation scenarios. Conventional smooth disk migration can readily produce resonance locking but almost invariably yields mutual inclinations $a_c = 0.2305^{+0.0069}_{-0.0081}$3, unless augmented by passage through higher-order inclination resonances (e.g., the 4:2 resonance). However, such passages require much higher eccentricities (e.g., $a_c = 0.2305^{+0.0069}_{-0.0081}$4 for the inner planet), inconsistent with the measured $a_c = 0.2305^{+0.0069}_{-0.0081}$5 and $a_c = 0.2305^{+0.0069}_{-0.0081}$6 in KOI-134. Simulations with parametrized migration, eccentricity damping, and inclination damping show that only under unrealistically low damping or high eccentricity can inclinations reach the observed $a_c = 0.2305^{+0.0069}_{-0.0081}$7 level.

Alternative channels, such as stellar $a_c = 0.2305^{+0.0069}_{-0.0081}$8-driven secular resonance sweeping or disk-induced nodal precession during disk dispersal, are also disfavored. KOI-134 c's semimajor axis is too large for stellar $a_c = 0.2305^{+0.0069}_{-0.0081}$9 effects and the required disk mass for nodal precession is implausibly high. As such, the observed architecture likely reflects a composite evolutionary history: initial resonance capture during disk presence followed by a later episode—possibly involving planet-planet scattering, external companion-induced tilting, or other dynamical mechanisms—that excited mutual inclination without breaking the resonance. The persistence of resonance amidst significant inclination argues for a dynamically cold epoch followed by impulsive reorientation (Nabbie et al., 3 Jul 2025).

6. Summary and Broader Significance

KOI-134 c is a non-transiting, sub-Saturn-mass planet $a_c = 0.22534^{+0.00052}_{-0.00050}$0 interior to the warm Jupiter KOI-134 b, with orbital period $a_c = 0.22534^{+0.00052}_{-0.00050}$1 days, semimajor axis $a_c = 0.22534^{+0.00052}_{-0.00050}$2 AU, and eccentricity $a_c = 0.22534^{+0.00052}_{-0.00050}$3. The system is in or very near exact 2:1 mean-motion resonance, with c’s eccentricity-type resonant angle librating, but host to a large mutual inclination ($a_c = 0.22534^{+0.00052}_{-0.00050}$4). C’s gravitational perturbations generate large TTVs and TDVs in KOI-134 b, robustly constraining c’s mass and orbital parameters through dynamical modeling.

The KOI-134 system exemplifies a rare but dynamically important planetary architecture: a resonant, non-coplanar giant pair. This challenges current formation models, as no single canonical process can reproduce all aspects of the present architecture. The system provides a benchmark for testing the interplay of resonant capture, inclination-excitation mechanisms, and long-term stability in multiplanet dynamics. Its distinctive properties also lead to pronounced observational effects, including the predicted future disappearance of KOI-134 b’s transits, and the persistent non-transiting status of c. As such, KOI-134 c is a key system for both theoretical and observational studies of planetary system evolution (Nabbie et al., 3 Jul 2025).