K2-138 g: Spitzer Spots a Sixth Planet for the Citizen Science System

$K2$ greatly extended $Kepler$'s ability to find new planets, but it was typically limited to identifying transiting planets with orbital periods below 40 days. While analyzing $K2$ data through the Exoplanet Explorers project, citizen scientists helped discover one super-Earth and four sub-Neptune sized planets in the relatively bright ($V=12.21$, $K=10.3$) K2-138 system, all which orbit near 3:2 mean motion resonances. The $K2$ light curve showed two additional transit events consistent with a sixth planet. Using $Spitzer$ photometry, we validate the sixth planet's orbital period of $41.966\pm0.006$ days and measure a radius of $3.44^{+0.32}_{-0.31}\,R_{\oplus}$, solidifying K2-138 as the $K2$ system with the most currently known planets. There is a sizeable gap between the outer two planets, since the fifth planet in the system, K2-138 f, orbits at 12.76 days. We explore the possibility of additional non-transiting planets in the gap between f and g. Due to the relative brightness of the K2-138 host star, and the near resonance of the inner planets, K2-138 could be a key benchmark system for both radial velocity and transit timing variation mass measurements, and indeed radial velocity masses for the inner four planets have already been obtained. With its five sub-Neptunes and one super-Earth, the K2-138 system provides a unique test bed for comparative atmospheric studies of warm to temperate planets of similar size, dynamical studies of near resonant planets, and models of planet formation and migration.


INTRODUCTION
The NASA K2 mission searched for exoplanets in different fields spanning the ecliptic plane, subsequent to the loss of two reaction wheels, which inhibited the Kepler spacecraft's ability to precisely point at the original Kepler field for extended durations. Using solar pres-sure and thrusters, Kepler was able to point to fields along the ecliptic plane for a period of ∼83 days each before the spacecraft was rotated to prevent sunlight from entering the telescope (Putnam & Wiemer 2014;Howell et al. 2014). K2 has so far enabled the discovery of 425 new planets and an additional 889 planet candidates. 1 K2-138 was the first K2 planet system discovered by citizen scientists through the Exoplanet Explorers 2 program on the Zooniverse 3 (Christiansen et al. 2018). The citizen scientists were able to identify four sub-Neptune sized planets by visual inspection of the light curve. A closer inspection of the diagnostic plots from the TERRA algorithm 4 (Petigura et al. 2013b,a) elucidated a super-Earth interior to the orbits of the other four planets. Using LcTools (Kipping et al. 2015;Schmitt et al. 2019), Christiansen et al. (2018) also identified two additional transits 41.97 days apart, indicating a possible sixth planet for the system. Lopez et al. (2019) obtained radial velocity (RV) measurements of K2-138 with HARPS, yielding mass measurements of 3.1 ± 1.1, 6.3 +1.1 −1.2 , 7.9 +1.4 −1.3 , and 13.0 ± 2.0 M ⊕ for planets b, c, d, and e, respectively. Precise masses for K2-138 f and the putative planet K2-138 g were not measured. K2-138 f has an orbital period of 12.76 days, about half of the 24.7 ± 2.2 day stellar rotation period, and its signal was likely absorbed by the Gaussian process regression used to remove stellar activity. This process also likely muted the signal of K2-138 g. Lopez et al. (2019) placed upper limits at 99% confidence of 8.7 and 25.5 M ⊕ on K2-138 f and g, respectively. Due to the near 3:2 orbital resonances, K2-138 is amenable to transit timing variation (TTV) measurements to constrain planet masses. Using their measured masses and assuming zero eccentricity, Lopez et al. (2019) computed TTV amplitudes between 2.0 and 7.3 minutes for the inner five planets, similar to the amplitudes computed by Christiansen et al. (2018). Though Christiansen et al. (2018) were not able to detect significant TTVs in the 30 minute cadence K2 data, higher cadence observations with instruments such as CHEOPS, which were scheduled for late 2020 (Program ID 017 (EP); PI: T. Lopez), should allow TTV mass measurements of planets c, d, and e, making K2-138 an important benchmark system for comparing TTV and RV masses. Since RV mass measurements are currently limited to host stars brighter than V 13, TTVs enable mass measurements for a much wider pool of planets (Holczer et al. 2016). However, fewer than 10 systems have both RV and TTV mass measurements, and detection sensitivity may bias RV measurements for planets with orbital periods larger than 11 days (Mills & Mazeh 2017;Petigura et al. 2018). These reasons highlight the importance of adding new TTV/RV benchmark systems 2 https://www.zooniverse.org/projects/ianc2/exoplanet-explorers 3 https://www.zooniverse.org/ 4 https://github.com/petigura/terra in order to cross-check masses between measurement techniques.
In this paper we verify the outermost planet K2-138 g with an orbital period of 41.96645 +0.00603 −0.00665 days. This adds to the nine systems with six or more planets currently known, makes K2-138 the K2 discovered system with the most planets, 5 and yields one of the longest period K2 planets. Using the Spitzer Space Telescope, we observed a third transit of K2-138 g within one hour of the time predicted from the K2 ephemeris. We present our observations and data reduction in Section 2 and discuss our results in Section 3.

Stellar Classification
We obtained a 0.38 to 0.7 µm spectrum of K2-138 using the Goodman spectrograph (Clemens et al. 2004) on the Southern Astrophysical Research Telescope (Program ID 2019A-0364; PI: K. Hardegree-Ullman), and a 0.7 to 2.4 µm spectrum using the SpeX spectrograph (Rayner et al. 2003) on the NASA Infrared Telescope Facility (Program ID 2017A-106; PI: K. Hardegree-Ullman). We followed the procedures outlined in § 2.1-2.3 of Hardegree-Ullman et al. (2019) for observing the target and reducing the data. We compared the combined optical and infrared spectra between 0.38 and 1 µm to optical SDSS spectral templates from Kesseli et al. (2017) following the procedures outlined in § 4.1 in Hardegree-Ullman et al. (2020), which yielded a spectral type of G8 V, consistent with the spectral type found by Lopez et al. (2019). Figure 1 shows our 0.38 to 1 µm spectrum compared to G7 V, G8 V, and G9 V template spectra.
In Table 1 Lopez et al. (2019), andHardegree-Ullman et al. (2020). The measured parameters are all consistent to within 1σ, except for one measurement of log(g) which is within 2σ. It is reassuring that different data sets and pipelines yield similar results, but when it comes to calculating planet parameters, small differences in stellar parameters can have a large impact. For large surveys and population studies of exoplanets, it is crucial to have a uniformly   derived set of stellar parameters (e.g., Fulton et al. 2017;Berger et al. 2020;Hardegree-Ullman et al. 2020). For individual systems, however, it is typical to choose a single set of stellar parameters, which may be susceptible to systematic bias. Rather than cherry picking measurements from different references, we combined all the available measurements for T eff , log(g), [Fe/H], and M ⋆ . Instead of using a weighted mean, which would produce uncharacteristically small uncertainties, 6 we instead employed the following Monte Carlo method. For each measurement with symmetric uncertainties, we randomly drew 10 4 values from a Gaussian distribution, and for asymmetric uncertainties we drew 10 4 values from a split normal distribution. The posterior distributions were concatenated and we took the median, 16th, and 84th percentiles of the resultant distribution as our measurement and errors.
For our analysis, we used the instrumental systematics-corrected light curve ( Figure 2) produced by the k2phot 7 pipeline (Petigura et al. 2015;Aigrain et al. 2016). We compared the k2phot light curve to those produced by EVEREST  and K2SFF (Vanderburg & Johnson 2014), and found that the k2phot light curves had the lowest overall RMS scatter and the fewest outliers. We first masked out data that was flagged in the k2phot pipeline as a thruster fire event or an outlier in background flux. Periodic transit signals were initially found by flattening the light curve with a Savitsky-Golay filter over a window of 101 points (∼50 hours), then running a box least squares periodogram, iteratively masking out the higher signal-to-noise transits until there were no more convincing planet signals in the data. This search gave estimates of planet periods, transit times, and transit depths. Next, we made use of the exoplanet 8 toolkit to model stellar variability using a Gaussian process with a simple harmonic oscillator kernel, while simultaneously fitting planet transits as described in Foreman-Mackey et al. (2017). In order to simultaneously fit the K2 and Spitzer data ( § 2.4), we used the flattened the light curve by subtracting the Gaussian process stellar model without fitting out the planet transits, which is shown in the middle panel of Figure 2.

IRAC Photometry
We observed K2-138 with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope (DDT 13253; PI J. Christiansen) at the predicted transit time of the putative sixth planet. We used Channel 2 (4.5 µm) since it is less affected by intrapixel sensitivity variations than Channel 1 (3.6 µm). The observation began with a 30 minute pre-observation stare which was discarded in the analysis, but included in the Astronomical Observation Request to allow the telescope and instruments to settle after slewing (Grillmair et al. 2012). To minimize the pixel-phase effect and achieve a pointing accuracy to within ∼0.1 pixel, we conducted pre-observations in peak-up mode using the Pointing Calibration and Reference Sensor (Ingalls et al. 2012).
Observations of K2-138 were conducted between 2018 March 15 and 2018 March 16 for a total duration of 11 7 https://github.com/petigura/k2phot 8 https://exoplanet.dfm.io/en/stable/ hours centered near the predicted time of mid-transit from the K2 ephemeris of the sixth planetary signal found by Christiansen et al. (2018). Individual frame exposure times were set to two seconds to stay in the linear regime of the detector for this bright target. The subarray mode was used to minimize readout times and data volume. In total, 19,840 individual frames were taken.
We performed centroiding and aperture photometry using photutils (Bradley et al. 2019), fitting a 2D Gaussian to each image. To select the optimal aperture radius, we computed photometry from the centroid positions using fixed radii between 1.5 and 3.0 pixels in 0.1 pixel increments. Background levels were found by the method described by Knutson et al. (2011). This process entails masking out the regions within a radius of 12 pixels from the centroid along with the central two rows and columns, then finding the median background value of the pixels after clipping 3σ outliers.
We modeled systematics in the Spitzer light curves using pixel-level decorrelation (PLD; Deming et al. 2015).
PLD has become a premier technique for correcting Spitzer systematics in planet transit analyses (e.g., Beichman et al. 2016;Benneke et al. 2017;Dressing et al. 2018;Feinstein et al. 2019;Livingston et al. 2019;Berardo et al. 2019), and was developed to account for intra-pixel sensitivity variations which produce intensity fluctuations in the photometry. In our analysis, we used PLD to model the Spitzer systematics simultaneously with the exoplanet system parameters. The full model is described by: where w i are individual time-independent pixel weights in the n selected pixels in the region centered on the star, D i (t) is the observed flux (or counts) in the individual pixels of the selected region for each time step t, m is the slope of a linear temporal ramp, and M tr (θ, t) is the transit model with model parameters θ. The first part of this equation normalizes the individual pixel intensities so their sum at each time step is unity. In our analysis we used a 3×3 pixel region centered around the brightest pixel.

Transit Fitting
Using emcee (Foreman-Mackey et al. 2013), we simultaneously model the K2 and Spitzer data, computing the posterior probability distributions for six transiting planets and the Spitzer systematics. To model the transits we used batman , which solves the analytic equations for an exoplanet transit as derived in Mandel & Agol (2002). We computed posterior probability distributions for the mid-transit times T 0 , orbital periods P , the ratios of planet to star radii R p /R ⋆ , the scaled semi-major axes a/R ⋆ , impact parameters b, two sets of quadratic limb darkening coefficients q 1 and q 2 (one set for K2 and another for Spitzer ), and the nine pixel weights w i and linear slope m from Equation 1. We performed an autocorrelation analysis 9 to ensure chain convergence. Due to our large set of parameters, we used 500 walkers and 250,000 steps.
The resultant K2 light curve fit using the median values of the posteriors is shown in Figure 2, with the phase folded light curves shown in Figure 3. Figure 4 shows a 9 https://emcee.readthedocs.io/en/latest/tutorials/autocorr/ clear transit event in the Spitzer data for K2-138 g, confirming the existence of a sixth planet in the K2-138 system. We note that Christiansen et al. (2018) obtained high resolution AO imaging of the K2-138 system, ruling out nearby stellar companions that could contaminate or mimic a planet signal. Further, since K2-138 is a multi-planet system, it is more likely that additional transit-like signals come from another planet (validation by multiplicity, e.g., Lissauer et al. 2014;Sinukoff et al. 2016). Table 2 lists all the derived planet parameters for the K2 and Spitzer data. We compare the K2 and Spitzer light curves for K2-138 g in Figure 5. The transit durations for the light curves are nearly identical, but the transit depth posterior distributions show a slightly larger radius in the Spitzer data, although the difference is < 1σ.   The Spitzer radius is larger, but it is still consistent with the K2 radius within 1σ. It is also possible that systematics could bias the radius measurements. For example, having only two transits in the 30 minute cadence K2 data means that any outliers could skew the measured transit depth. Additional transits at these and other wavelengths will be necessary to constrain atmospheric properties of this planet. · · · · · · · · · · · · · · · · · · · · · † Equilibrium temperatures were computed assuming a Bond albedo of 0.3. ‡ Period, a, b, and i were computed jointly with the K2 data.

Near Resonances and Gap Planets
The ratio of orbital periods between successive K2-138 planets are: c:b = 1.513, d:c = 1.518, e:d = 1.529, f:e = 1.544, g:f = 3.290. In order to determine how close to 3:2 resonance these planets are, we estimated the mean motion resonance widths using the program from Volk & Malhotra (2020) 10 , which is based on the analytical derivations of resonance widths in the singleplanet limit from Murray & Dermott (2000). For this calculation, we used the masses of planets K2-138 b, c, d, and e from Lopez et al. (2019), and estimated masses of K2-138 f and g (6.72 +8.04 −3.86 and 8.94 +12.89 −5.91 M ⊕ ) from mass-radius relationships (Ning et al. 2018). The results of this calculation, out to fourth order mean motion resonances, are shown in Figure 6. Within the upper and lower planet mass limits, K2-138 b, c, d, and e are near (within a few half-widths) their mutual 3:2 resonances at low eccentricity, but the outer pair of planets are not near any low-order resonances.
The sizeable gap between K2-138 f and g leads to speculation that there could be additional non-transiting planets in the system. Indeed, Gilbert & Fabrycky (2020) suggest ∼20% of high multiplicity planet systems host additional planets in the gaps between detected planets. Each consecutive planet pair of K2-138 has period ratios that slip further away from 3:2, and assuming the orbital period ratios continued at 1.544 (f:e), planets could be expected with orbital periods near 19.70, 30.42, and 46.98 days. However, if the K2-138 planets were all in perfect 3:2 resonance with planet b, there would orbits at 17.87, 26.80, and 40.21 days. Without additional data, we are unable to conclude whether or not K2-138 g would be near a 3:2 resonance with a planet in the gap.
Multi-planet systems have been found to be highly coplanar (e.g., Fabrycky et al. 2014;Zhu et al. 2018;Gilbert & Fabrycky 2020), however, the more distant a planet orbits, the closer to 90 • inclination it must be to be in a transiting geometry. Assuming orbital periods of 19.70 and 30.42 days, planets around K2-138 would need to be at inclinations above 88.9 • and 89.2 • , respectively, for us to observe them in transit. Even within the Solar System, the planets are nearly coplanar, yet they still have mutual inclinations between 0. • 33 and 6.3 • (Winn & Fabrycky 2015).
We further explore the possibility of planets within the gap between planets f and g using DYNAMITE 11 , which uses population statistics to predict previously Figure 6. The location and analytically estimated widths of mean-motion resonances for the K2-138 system. Each planet is plotted in relative size to the other planets along the discontinuous y-axis indicating orbital period. Eccentricity is given along the x-axis, and extending from each planet is a line out to the eccentricity at which the planet would cross another planet's orbit. Horizontal dashed lines indicate the locations of interior (e.g., 3-b:2) and exterior (e.g., 3:2c) resonances up to fourth order, color coded to match the label. The shaded regions surrounding each resonance line are the resonance widths corresponding to the lower (darker) and upper (lighter) planet mass limits. Planets b, c, d, and e are sufficiently near their mutual 3:2 resonances at low eccentricity for their dynamics to be affected, likely inducing TTVs (see Section 3.2). undetected planets (Dietrich & Apai 2020). This model takes inputs of stellar parameters (radius, mass, temperature) and known planet parameters (inclination, radius, period), and yields probability distributions where the population models predict a planet or planets might exist. We considered four different scenarios as inputs to DYNAMITE, which are shown in Figure 7:

Masses and TTVs
Due to its distance from first order resonance, TTV measurements for K2-138 g would be difficult.  Lopez et al. (2019) and the aforementioned estimated masses of K2-138 f and g. We also assumed zero eccentricity. Our average six minute (1σ) K2 timing precision was insufficient to measure TTVs for this system, however, higher cadence (one minute) observations with CHEOPS should improve the timing precision enough to allow detection of TTVs of the inner five planets.
In measuring the masses of the inner four K2-138 planets, Lopez et al. (2019) did not identify additional planets, though additional signals might have been absorbed by their Gaussian process to fit out stellar activity at the 5.6 m s −1 level. The mass measurement of K2-138 f was hindered by its orbital period of 12.8 days, near half the 24.7 day stellar rotation period. The stellar rotation period might also hinder detection of a planet in orbit near the next 3:2 resonance beyond planet f around 20 days. Future planet searches and mass measurements for this system would likely benefit from simultaneous photometric and RV observations, as Kosiarek & Crossfield (2020) suggest this could enhance the precision of RV measurements. Lopez et al. (2019) were unable to reliably measure a mass of K2-138 g either, but assuming a mass of 8.94 +12.89 −5.91 M ⊕ , we predict an RV semi-amplitude of 1.79 +2.56 −1.18 m s −1 . If there were planets between f and g of similar masses to the other planets in the system, we would expect them to have RV semi-amplitudes between 1.5 and 2.5 m s −1 , which would make them similarly difficult to detect due to stellar activity levels.
We note that the outer five planets of K2-138 are all sub-Neptunes similar in size, and planet b is likely a rocky super-Earth with a density of 5.01 +2.73 −2.00 g cm −3 . Common sizing of multi-planet systems has previously been found for Kepler systems (e.g. Millholland et al. 2017;Wang 2017;Weiss et al. 2018;Gilbert & Fabrycky 2020). From a planet formation standpoint, Adams et al. (2020) found that energy optimization occurs when planets are nearly equal in mass for low-mass (super-Earth/sub-Neptune) planet systems, which is consistent with what we see with K2-138. Though, we note that the outer planets of K2-138 have larger radii than the inner planets, a trend consistent with the findings of Ciardi et al. (2013), Millholland et al. (2017), Kipping (2018), and Weiss et al. (2018), possibly the result of enhanced photoevaporation closer to the star. We plot the planet radii with respect to incident stellar flux for the K2-138 planets, compared to the population of K2 planets (shown as density contours) from Hardegree-Ullman et al. (2020) in Figure 8. K2-138 b has incident flux (F ⊕ ) over 400 times higher than Earth and is the only planet in the system below the planet radius valley. The other planets in the system receive less than 250 F ⊕ , apparently low enough to retain an atmosphere.
From Figure 8, it appears that many of the K2-138 planets are inflated relative to their counterparts with similar incident stellar flux. If the system was relatively young, we would expect the planets to still be undergo- ing mass loss. Lopez et al. (2019) computed an age of 2.3 +0.44 −0.36 Gyr for K2-138 based on chromospheric emission, and 2.8 +3.8 −1.7 Gyr from their joint radial velocity, light curve, and spectral energy distribution analysis. Similarly, we input photometry, stellar parameters, and a rotation period of 24.7 days into the isochrone fitting with gyrochronology package stardate 12 and compute an age of 2.8 ± 0.3 Gyr, consistent with Lopez et al. (2019). However, a visual assessment of the raw flux in Figure 2 and a Lomb-Scargle periodogram yields a significant peak corresponding to a period of ∼12.5 days. This rotation period corresponds to a younger age closer to ∼0.9 Gyr. We note that even at this younger age, it is unlikely the planets are still undergoing significant mass loss since this process occurs within the first few hundred Myrs (Lopez et al. 2012).
Another possibility for these relatively large planets is tidally induced radius inflation (Millholland 2019). The Kepler mission unveiled a statistical overabundance of planet pairs just outside of first-order mean motion resonances, specifically 2:1 and 3:2 (Lissauer et al. 2011;Millholland & Laughlin 2019). As noted in Section 3.1, most of the planet pairs of K2-138 fall just outside a 3:2 resonance. Tidal forces from the host star can push planets into near-resonant configurations, but host-star tides alone cannot explain how all the energy from this process is dissipated to keep planets in this configuration. Millholland & Laughlin (2019) showed that obliq- Figure 9. The orbits and planets of K2-138 and HD 158259, highlighting the similarities of the two systems. The orbits and the star sizes are to scale, but the planets are enlarged by 50× to show detail. Orbital distances of 0. 034, 0.046, 0.060, 0.080, 0.105, and 0.135 AU for HD 158259 b, c, d, e, f, and tentative planet (g), respectively, were computed from Kepler's third law using stellar mass and planet orbital periods from Hara et al. (2020). uity tides may be the source of energy dissipation that helps sculpt these near-resonant systems. Consequently, these tidal forces heat the planet interiors, leading to atmospheric inflation (Millholland 2019). We posit that K2-138 is a strong candidate for planet radius inflation due to obliquity tides.

Comparison to Other Multi-planet Systems
To date, there have only been nine other exoplanet systems with six or more confirmed planets 13 , including radial velocity discovered systems HD 10180 (6 planets), HD 219134 (6 planets), and HD 34445 (6 planets), and transiting systems Kepler-11 (6 planets), Kepler-20 (6 planets), Kepler-80 (6 planets), Kepler-90/KOI-351 (8 planets), TRAPPIST-1 (7 planets), and TOI-178 (6 planets). Perhaps most similar to K2-138, however, is the HD 158259 system, with five confirmed planets and a sixth candidate outer planet (Hara et al. 2020), all near 3:2 orbital mean motion resonances. Four of the confirmed planets were detected in radial velocity data with the SOPHIE spectrograph, and the innermost planet was found to be transiting in TESS data. The outermost candidate planet orbits every 17.4 days, close to the stellar rotation period, complicating confirmation of this planet. The five innermost planets of K2-138 and HD 158259 are each located at nearly identical distances to their host stars, as shown in Figure 9. We estimated HD 158259 planet radii for the five non-transiting planets using the planet masses from Hara et al. (2020) and the mass-radius relationships of Chen & Kipping (2017). These non-transiting planets are also all similar-sized sub-Neptunes, with estimated radii larger than 2 R ⊕ -again consistent with the aforementioned common sizing of multi-planet systems. Each respective planet in HD 158259 is slightly smaller than its counterpart K2-138 planet, which could be the result of HD 158259 being a larger host star (1.08 ± 0.10 M ⊙ ) that is more efficient at stripping away planetary atmospheres by intense irradiation (Ehrenreich et al. 2015). We note, however, that there are significant uncertainties in planetary mass-radius relationships. Without transit data, it is difficult to test whether or not this system undergoes tidal radius inflation as mentioned in Section 3.2.
We qualitatively compared the orbital spacing (a/R ⋆ ) of these high-multiplicity systems with transiting planets ( Figure 10). In addition to the K2-138 system, there is a sizeable gap between the outermost detected transiting planets of the Kepler-11, Kepler-20, and Kepler-80 Figure 10. Orbital spacing of systems with six or more planets, and at least one transiting planet. Systems are arranged from largest (top) to smallest (bottom) stellar host, and the regions are colored according to the host temperatures (Harre & Heller 2021). The width of the colored regions are scaled to the stellar radii. Planets are to scale with the stellar radii but enlarged by 10× for clarity, and are placed at their respective transiting inclination angles (randomly distributed above and below the stellar mid-point; gray line). Non-transiting planet locations are labeled in blue. For clarity, we did not plot HD 219134 h on this scale, but note that it is located at a/R⋆ = 857. In addition to K2-138, the systems Kepler-11, Kepler-20, HD 219134, and Kepler-80 have notable gaps between their outermost planets.
systems. HD 219134 has two transiting planets and four non-transiting planets detected via RV measurements, again with a large gap between the two outermost planets. This large outermost planet gap is also present in the RV system HD 34445. Notably, a non-transiting planet was identified in the gap between outer planets Kepler-20 f and d with RV data (Buchhave et al. 2016). As noted in Section 3.1, planets orbiting further out must be closer to i = 90 • to be in a transiting geometry, but RV and TTV data may uncover unseen planets. We encourage further investigations of this outer planet gap feature in high-multiplicity planet systems in order to disambiguate whether it is caused by observational biases or planet formation processes.

JWST, ARIEL, and Future Prospects
We computed the transmission spectroscopy metric (TSM) for the K2-138 planets as defined by Equations 1, 2, and 3 of Kempton et al. (2018). The TSM is the expected signal-to-noise for a 10 hour observing program with JWST/NIRISS. For planets b, c, d, and e, we used the planet masses measured by Lopez et al. (2019), and for planets f and g, we used our estimated masses. The equilibrium temperature was calculated assuming zero albedo and full day-night heat redistribution. The resultant TSM values are listed in Table 2, and are ∼20 for the outer planets, falling to 2.62 for the innermost planet b. These values are well below the recommended threshold of TSM > 90 for high quality atmospheric characterization of sub-Neptune sized planets. For now, the K2-138 planets are unlikely to be selected as highpriority targets for JWST observations. The European Space Agency Atmospheric Remotesensing Infrared Exoplanet Large-survey (ARIEL) space mission aims to gather transmission spectra of 1000 exoplanets during its four year mission (expected to launch in 2028) in order to study their composition, formation, and evolution. Edwards et al. (2019) compiled a list of potential targets for ARIEL, taking into account currently known stellar and planet parameters. K2-138 falls very near the average star system considered for this target list. The inner five planets of K2-138 also fall within the range of planets considered for the potential target list, but very few planets with orbital periods beyond ∼20 days will likely be considered, all but ruling out observations of K2-138 g. However, since K2-138 contains five similarly-sized sub-Neptunes with a ∼500 K range of equilibrium temperatures from warm to temperate, these planets might provide a unique test bed for comparative sub-Neptune atmosphere studies.
We have confirmed the existence of K2-138 g, solidifying K2-138 as the largest K2 multi-planet system. K2-138 g breaks the continuous near 3:2 mean motion resonance of the inner five planets, but the sizeable gap between K2-138 f and g hints at the possibility there could be additional non-transiting planets in this system. We encourage future observations of this potential key benchmark system to (1) constrain TTVs of the inner planets, (2) enable more precise masses and potential discovery of additional planets with simultaneous photometric and RV measurements, and (3) facilitate comparative atmospheric studies of warm to temperate sub-Neptune planets.