Long-period Jovian Tilts the Orbits of Two sub-Neptunes Relative to Stellar Spin Axis in Kepler-129

We present the discovery of Kepler-129 d ($P_d={7.2}_{-0.3}^{+0.4}$ yr, $m\sin i_d={8.3}_{-0.7}^{+1.1}{M}_{Jup}$, $e_d={0.15}_{-0.05}^{+0.07}$) based on six years of radial-velocity observations from Keck/HIRES. Kepler-129 also hosts two transiting sub-Neptunes: Kepler-129 b (P _b = 15.79 days, r _b = 2.40 ± 0.04 R _⊕) and Kepler-129 c (P _c = 82.20 days, r _c = 2.52 ± 0.07 R _⊕) for which we measure masses of m _b < 20 M _⊕ and $m_c={43}_{-12}^{+13}{M}_{\oplus }$. Kepler-129 is a hierarchical system consisting of two tightly packed inner planets and a massive external companion. In such a system, two inner planets precess around the orbital normal of the outer companion, causing their inclinations to oscillate with time. Based on an asteroseismic analysis of Kepler data, we find tentative evidence that Kepler-129 b and c are misaligned with stellar spin axis by ≳38°, which could be torqued by Kepler-129 d if it is inclined by ≳19° relative to inner planets. Using N-body simulations, we provide additional constraints on the mutual inclination between Kepler-129 d and inner planets by estimating the fraction of time during which two inner planets both transit. The probability that two planets both transit decreases as their misalignment with Kepler-129 d increases. We also find a more massive Kepler-129 c enables the two inner planets to become strongly coupled and more resistant to perturbations from Kepler-129 d. The unusually high mass of Kepler-129 c provides a valuable benchmark for both planetary dynamics and interior structure, since the best-fit mass is consistent with this 2.5 R _⊕ planet having a rocky surface.