Stellar Obliquities in Long-period Exoplanet Systems (SOLES) I: The Spin-Orbit Alignment of K2-140 b

Obliquity measurements for stars hosting relatively long-period giant planets with weak star-planet tidal interactions may play a key role in distinguishing between formation theories for shorter-period hot Jupiters. Few such obliquity measurements have been made to date due to the relatively small sample of known wide-orbiting, transiting Jovian-mass planets and the challenging nature of these targets, which tend to have long transit durations and orbit faint stars. We report a measurement of the Rossiter-McLaughlin effect across the transit of K2-140 b, a Jupiter-mass planet with period $P=6.57$ days orbiting a $V=12.6$ star. We find that K2-140 is an aligned system with projected spin-orbit angle $\lambda=0.5\pm9.7$ degrees, suggesting a dynamically cool formation history. This observation builds towards a population of tidally detached giant planet spin-orbit angles that will enable a direct comparison with the distribution of close-orbiting hot Jupiter orbital configurations, elucidating the prevalent formation mechanisms of each group.


INTRODUCTION
The obliquity of a star, or the degree of alignment between the star's spin axis and its companion planets' net orbital angular momentum axis, provides crucial insights into the dynamical evolution of the surrounding system. Trends in stellar obliquity can constrain the prevalence of various processes crafting the observed distribution of extrasolar planets (e.g. Winn et al. 2010;Albrecht et al. 2012;Winn & Fabrycky 2015). * NSF Graduate Research Fellow While the obliquities of many hot-Jupiter-hosting stars have been determined through measurements of the Rossiter-McLaughlin effect (Rossiter 1924;McLaughlin 1924), only a handful of these measurements have been made for systems with planets orbiting at larger distances from their host star, where starplanet tidal interactions are too weak to significantly influence the planets' orbital alignment. We refer to these wide-orbiting planets as "tidally detached." Because the probability p of transit for a given planet with semimajor axis a falls as p ∝ 1/a, fewer tidally detached transiting giant planets have been discovered to date as compared with closer-in hot Jupiters. Due to the volume-limited sample of bright stars, those tidally detached Jupiters that have been discovered are often challenging Rossiter-McLaughlin targets, orbiting stars too faint to obtain a sufficient number of high-resolution spectroscopic observations across a single transit with any but the largest existing ground-based telescopes.
The vast majority of spin-orbit angle measurements have, consequently, been made for planets on tight orbits with a/R * 12. However, the obliquities of stars hosting wider-separation giant planets may provide the evidence necessary to distinguish between the various proposed formation mechanisms for hot Jupiters (Winn & Fabrycky 2015;Dawson & Johnson 2018). Furthermore, in the classical quiescent formation framework for hot and warm Jupiters, the long tidal realignment timescales of these systems enable a direct measurement of the primordial dispersion in protoplanetary disk misalignments.
We present a Rossiter-McLaughlin measurement across the transit of K2-140 b with the High Resolution Echelle Spectrometer (HIRES; Vogt et al. 1994) on the 10-meter Keck I telescope. K2-140 (EPIC 228735255) is a V = 12.6 G5 star hosting a 1.019±0.070M J planet, K2-140 b, with a P = 6.57−day orbital period (Giles et al. 2018). With a/R * = 12.88, K2-140 b is one of the widest-separation tidally detached planets to date with a measured spin-orbit angle. This is the first measurement in our Stellar Obliquities in Long-period Exoplanet Systems (SOLES) survey designed to expand the sample of Rossiter-McLaughlin measurements for wide-separation exoplanets.
K2-140 b was first characterized using photometry from the K2 mission (Howell et al. 2014), a re-purposed extension of the Kepler mission (Borucki et al. 2010), after two of its reaction wheels failed. Giles et al. (2018) used radial velocity (RV) observations from CORALIE (Queloz et al. 2000) and the High Accuracy Radial velocity Planet Searcher (HARPS; Pepe et al. 2000) to confirm the planetary nature of K2-140 b and to constrain its physical and orbital properties.
The planet was independently characterized in Korth et al. (2019), which incorporated RV data from the FIbre-fedÉchelle Spectrograph (FIES) (Telting et al. 2014) within their analysis. While Giles et al. (2018) found a low but nonzero eccentricity for the planet (e = 0.120 +0.056 −0.046 ), Korth et al. (2019) found that the planet was consistent with e = 0. Data from both of these past studies is incorporated within our joint analysis.

OBSERVATIONS
We obtained 18 radial velocity measurements of K2-140 with the Keck/HIRES instrument from 9:25-15:40 UT on Feb 24, spanning a full transit of K2-140 b. Conditions were favorable throughout most of the observing period, with typical seeing ranging from 1.0 − 1.3 . A spike in humidity during the post-transit baseline observations led to a telescope dome closure at 14:40 UT, resulting in a ∼40-minute gap in data before the last radial velocity measurement.
All RV observations were obtained using the C2 decker (14 × 0.861 , R = 60, 000) and an iodine absorption cell, which imprints a dense forest of molecular iodine features onto each spectrum to enable high Doppler precision (Butler et al. 1996). The 14 length of the C2 decker allows for direct sky subtraction, improving RV precision for faint stars such as K2-140. The median exposure time was 1, 119 seconds, with ∼34k exposure meter counts per spectrum.
Our dataset was reduced using the California Planet Search pipeline outlined in Howard et al. (2010), and we obtained a typical signal-to-noise ratio of 71 per pixel from the reduced spectra. The HIRES radial velocity results and uncertainties can be found in Table 1 and are shown in the rightmost panel of Figure 1. We include the S-index and associated uncertainty at each observation in Table 1 for reference.
We also obtained a 45-minute iodine-free HIRES exposure of K2-140 using the B3 decker (14.0 × 0.574 , R = 72, 000) two nights after the measurement of the Rossiter-McLaughlin effect, during UT Feb 26. This template observation was used to calibrate our RVs and to precisely determine stellar parameters (see Section 4). Conditions were favorable during this measurement, and seeing was 1.2 . Our reduced template had a signal-to-noise ratio of 123 per pixel (106k exposure meter counts).

OBLIQUITY MODELING
To determine the sky-projected spin-orbit angle λ for K2-140 b, we used the allesfitter Python package (Günther & Daylan 2020) to jointly model the in-transit HIRES radial velocity data together with photometry from K2 and archival radial velocity datasets available from the FIES, CORALIE, and HARPS spectrographs.
All fitted parameters listed in Table 2 were allowed to vary, and each parameter was initialized with uniform priors. Initial guesses for P , T 0 , cos i, R p /R , (R + R p )/a, K, √ e cos ω, and √ e sin ω were obtained using values from Giles et al. (2018). The two limb darkening coefficients q 1 and q 2 were each initialized with values of 0.5. We accounted for potential radial velocity offsets between each separate spectrograph, with priors bounded by ±1000 m/s. Jitter terms were modeled separately for each instrument and added in quadrature to  the instrumental uncertainties. λ was allowed to vary between −180 • and +180 • . We ran an affine-invariant Markov Chain Monte Carlo (MCMC) analysis with 100 walkers to sample the posterior distributions of all model parameters. The best-fit model parameters and their associated 1σ uncertainties were extracted after obtaining 500,000 accepted steps per walker. Our results are listed in Table 2 and are in good agreement with the associated values obtained by Giles et al. (2018) and Korth et al. (2019).
The best-fit joint model is shown in Figure 1 together with each dataset included in the analysis, as well as the residuals of each fit. The fitted and derived parameters corresponding to this model are provided in Table 2. We obtain a low but nonzero eccentricity e = 0.069 +0.042 −0.028 for K2-140 b, in agreement with the value derived by Giles et al. (2018). K2-140 is consistent with alignment, with λ = 0.5 • ± 9.7 • and v sin i * = 2.51 ± 0.38 km/s.

STELLAR PARAMETERS
An understanding of host star properties can help to contextualize the evolutionary pathways through which a system may have reached its current state. We extracted stellar parameters from our Keck/HIRES template spectrum of K2-140 using the data-driven spectroscopic modeling program The Cannon (Ness et al. 2015;Casey et al. 2016), following the methods of Rice & Brewer (2020).
Given a set of uniformly processed input training spectra and associated stellar "labels" -that is, stellar parameters and elemental abundances -The Cannon constructs a generative model describing the probability density function of flux at each wavelength as a function of the labels. The model can then be applied to a new set of spectra, uniformly processed in the same manner as the training set, to obtain the associated stellar labels.
We trained The Cannon using the uniformly analyzed Spectral Properties of Cool Stars (SPOCS) catalogue (Brewer et al. 2016) of 18 stellar labels, including 3 global stellar parameters (T eff , log g, v sin i * ), and 15 elemental abundances: C, N, O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Fe, Ni, and Y. The sample of 1202 Keck/HIRES spectra vetted in Rice & Brewer (2020) was applied as our training/test set. The continuum baseline of each spectrum was uniformly fit and divided out using the iterative polynomial fitting procedure outlined in . Then, we split the SPOCS sample into an 80%/20% training/test split, applied the telluric mask from Rice & Brewer (2020) to all spectra, and trained the model, using the scatter of the test set results to determine the uncertainties of each extracted parameter. Finally, the trained model was applied to the newly acquired K2-140 Keck/HIRES template spectrum. A segment of the obtained model spectrum is shown in comparison with the HIRES template data in Figure 2.
Our results are provided in Table 3 together with estimates from previous works for reference. Our values for T eff , log g, and [Fe/H] are in agreement with those acquired in previous studies. The v sin i * value obtained using The Cannon is lower than previous estimates but is in agreement with the RM joint fit results and Korth et al. (2019) within 1σ. We find that K2-140 is metalenriched relative to solar abundances, consistent with past evidence showing that short-period giant planets are more common around metal-rich stars . At T eff = 5585 K, K2-140 should have a convective envelope. In the framework of equilibrium tides, the planet's timescale for realignment from turbulent friction would thus follow where τ CE is the realignment timescale for host stars with convective envelopes, and M p /M * is the planet-tostar mass ratio (Zahn 1977;Albrecht et al. 2012). For K2-140 b, τ CE = 1.2 × 10 13 yr, longer than the age of the Universe. As a result, we conclude that K2-140 b was likely aligned at the time of protoplanetary disk dispersal.
The theory of equilibrium tides presented in Equation 1 is employed as a simplified heuristic for a broader theoretical framework that has substantially advanced over recent years (Ogilvie 2014). A key problem in the equilibrium tides framework is that, under standard assumptions, hot Jupiters should experience rapid orbital decay. Lai (2012) demonstrated that one component of the tidal potential (the "obliquity tide") can excite inertial waves in the convective envelopes of cool stars, which, when damped, can enhance dissipation of the stellar obliquity without shrinking the companion's orbit. Obliquity tides have been further explored in additional work (Ogilvie 2013;Lin & Ogilvie 2017;Anderson et al. 2021). While obliquity tides are not included in our analysis, previous work has demonstrated that Equation 1 is a useful heuristic revealing that lowobliquity systems tend to have shorter tidal timescales than high-obliquity systems (Albrecht et al. 2012).
K2-140 b is one of only 13 giant planets with a measured spin-orbit angle at a/R * > 12, as shown in the left panel of Figure 3. Of these planets, it is one of only a few with an aligned orbit despite its long tidal realignment timescale. The alignment of K2-140 b suggests that at least some hot Jupiters form through quiescent pathways, such as in-situ formation or disk migration in an initially aligned disk.
The highly misaligned planets at large a/R * each have high eccentricities (right panel of Figure 3), indicative of strong dynamical interactions that may have produced both elevated obliquities and eccentricities. The coexistence of this population with the dynamically quiescent K2-140 system at large a/R * suggests that there are multiple hot Jupiter formation channels.
No neighboring planets have yet been found in the K2-140 system, despite previous observations showing that inner, coplanar companions are common for slightly longer-period warm Jupiters (10 < P < 200 days; Huang et al. 2016) and predictions that outer, mutuallyinclined companions with P 100 days should be regularly produced by in-situ hot Jupiter formation via core accretion (Batygin et al. 2016). However, additional low-mass or distant planets in the system cannot be ruled out by the existing observations, with radial velocity residuals of tens of m/s and an observing baseline of only ∼4 years.

Motivation for Additional Obliquity Measurements in Tidally Detached Systems through SOLES
Existing Rossiter-McLaughlin measurements of wideorbiting planets suggest that this population may have an intrinsically different spin-orbit angle distribution from the shorter-period hot Jupiter population. In particular, the trend of systematically lower obliquities observed for stars at temperatures below the Kraft break ( T eff ≈ 6100 K; Kraft 1967;Winn et al. 2010;Schlaufman 2010), which may result from tidal damping (e.g. Wang et al. 2021), is not immediately evident at longer orbital periods, as shown in Figure 4.
Conversely, Figure 4 reveals tentative (2.79σ) evidence that relatively long-period planets (P > 5 days) around hot stars are preferentially less misaligned than their shorter-period counterparts. Furthermore, the most misaligned P > 5 day planet in the middle panel, KELT-6 b (P = 7.8 days, |λ| = 36 ± 11 • ; Damasso et al. 2015), orbits a star directly bordering the nominal Kraft break (T eff = 6102 K; Collins et al. 2014), making it an ambiguous member of its group. We note that hot stars typically have larger radii and thus stronger tidal dissipation at a given orbital period; however, a similar trend has also been previously suggested for a smaller population of planets at large a/R * (Yu et al. 2018).
Longer-period planets have comparatively long tidal alignment timescales, especially around hot stars. The small spin-orbit angles observed for exoplanets around hot stars, therefore, may suggest that protoplanetary disks tend to be aligned at the time of gas dispersal (in contrast with a primordial spin-orbit misalignment from a tilted disk ;Batygin 2012;Spalding & Batygin 2015). In this case, hot Jupiters would need to obtain their misalignments after the protoplanetary disk has dispersed,  favoring high-eccentricity migration as an important formation mechanism (e.g. Wu & Murray 2003;Fabrycky & Tremaine 2007). Most of the misaligned planets at P > 5 days are on eccentric orbits, suggesting that interactions with stellar or planetary companions in the systems may have induced both elevated obliquities and eccentricities.
If a growing sample of spin-orbit angles at P > 5 days reveals that cool stars below the Kraft break are preferentially more misaligned than hotter stars above the Kraft break, this may also be suggestive of interactions with additional companions. In the presence of an external, inclined Jovian-mass companion, spin-orbit misalignments can be excited preferentially in cool star systems due to a secular resonance between the host star's spin axis precession frequency and nodal precession induced by interactions with the companion (Anderson & Lai 2018  Measurements for companion planets with M < 0.3MJ are shown at lower opacity. The cumulative sums in the bottom panel include only planets above the Kraft break, comparing the P > 5 day population (purple) with 5000 P < 5 day samples (gray), randomly sampled without replacement. The histogram on the bottom right shows a vertical cut through the final cumulative sum distribution. The P > 5 day cumulative sum is a 2.79σ outlier from the random draw distribution.
continue to rapidly rotate over their lifetimes. As a result, irrespective of the external giant planet companion rate around hot stars, this mechanism should occur only in cool star systems.
The substantially inclined companions necessary to induce this mechanism would also be capable of exciting the inner planet's orbital eccentricity. This may be reflected by the apparent increase in both misalignments and eccentricities for long-period planets around cool stars, shown in the top panel of Figure 4. Conversely, the alignment of K2-140 b and similar planets on loweccentricity orbits indicates that they should not have nearby giant planet companions with large (> 10 • ) mutual inclinations within ∼ 2 au.
If the trend of large misalignments for cool stars hosting tidally detached planets persists, while misalignments remain small for their hot star counterparts, further monitoring would be warranted to constrain the long-period giant planet companion rate for misaligned cool star systems. Additional observations are needed to parse the emerging relationship between stellar temperature, eccentricity, and obliquity for wide-separation planets.