arxiv: v1 [astro-ph.ep] 1 Apr 2019

Transcription

1 Astronomy & Astrophysis aepte ESO 219 April 2, 219 The Transit Least Squares Survey I. Disovery an valiation of an Earth-size planet in the four-planet system K2-32 near the 1:2:5:7 resonane René Heller 1, Kai Roenbek 1, 2, an Mihael Hippke 3 arxiv: v1 [astro-ph.ep] 1 Apr Max Plank Institute for Solar System Researh, Justus-von-Liebig-Weg 3, 3777 Göttingen, Germany, heller@mps.mpg.e 2 Inst. for Astrophysis, Georg-August-Univ. Göttingen, Frierih-Hun-Platz 1, 3777 Göttingen, Germany, roenbek@mps.mpg.e 3 Sonneberg Observatory, Sternwartestraße 32, Sonneberg, Germany, mihael@hippke.org Reeive 15 February 219 / Aepte 3 Marh 219 ABSTRACT We apply, for the first time, the Transit Least Squares (TLS) algorithm to searh for new transiting exoplanets. TLS has been evelope as a suessor to the Box Least Squares (BLS) algorithm, whih has serve as a stanar tool for the etetion of perioi transits. In this proof-of-onept paper, we emonstrate how TLS fins small planets that have previously been misse. We showase TLS s apabilities using the K2 EVEREST-etrene light urve of the star K2-32 (EPIC ), whih has been known to have three transiting planets. TLS etets these known Neptune-size planets K2-32 b,, an in an iterative searh an fins an aitional transit signal with a high signal etetion effiieny (SDE TLS ) of 26.1 at a perio of We show that this aitional signal remains etetable (SDE TLS = 13.2) with TLS in the K2SFF light urve of K2-32, whih inlues a less optimal etrening of the systemati trens. The signal is below ommon etetion threshols, however, if searhe with BLS in the K2SFF light urve (SDE BLS = 8.9) as in previous searhes. Markov Chain Monte Carlo sampling with the emee software shows that the raius of this aniate is R. We analyze its phase-fole transit light urve using the vespa software an alulate a false positive probability FPP = Taking into aount the multipliity boost of the system, we estimate an FPP < , formally valiating K2-32 e as a planet. K2-32 now hosts at least four planets that are very lose to a 1:2:5:7 mean motion resonane hain. The offset of the orbital perios of K2-32 e an b from a 1:2 mean motion resonane is in very goo agreement with the sample of transiting multi-planet systems from Kepler, lening further reene to the planetary nature of K2-32 e. We expet that TLS an fin many more transits of Earth-size an even smaller planets in the Kepler an K2 ata that have hitherto remaine unetete with algorithms that searh for box-like signals. Key wors. elipses methos: ata analysis planets an satellites: etetion planets an satellites: iniviual: K2-32 stars: planetary systems tehniques: photometri 1. Introution Both the ata from the Kepler primary mission (K1; Boruki et al. 21), whih operate from 29 to 213, an from the repurpose K2 mission (Howell et al. 214) that worke from 214 to 218 have been subjet to extensive transit searhes. Most of their onfirme or valiate planets (2338 from K1 an 359 from K2) 1 an of their aniates yet to be onfirme (2423 from K1 an 536 from K2) have been foun using the Box Least Squares (BLS) transit searh algorithm (Kovás et al. 22) or similar algorithms searhing for box-like flux ereases in stellar light urves (Batalha et al. 213; Vanerburg et al. 216; Crossfiel et al. 216; Mayo et al. 218; Livingston et al. 218b,a; Yu et al. 218; van Sluijs & Van Eylen 218; Crossfiel et al. 218). We evelope the Transit Least Squares (TLS) algorithm (Hippke & Heller 219) as the suessor of BLS in orer to be even more sensitive to smaller, possibly sub-earth-size planets. Rather than searhing for box-like flux ereases in the light urve, TLS is base on an analytial transit moel with stellar limb arkening (Manua et al. 1977; Manel & Agol 22). The signal etetion effiieny of TLS is onsequently improve 1 on 28 Marh 219 ompare to BLS, while the false positive rate is also suppresse (Hippke & Heller 219). Here we use TLS to searh for hitherto unknown planets in the K2 ata of K2-32 (EPIC ) an we present our first isovery from our new ata analysis ampaign, the TLS Survey. Three planets have previously been reporte aroun K2-32 by Vanerburg et al. (216), who formally esignate them as aniates, an by Crossfiel et al. (216), who valiate them as planets K2-32 b,, an using aitional high-resolution spetrosopy an inepenent stellar photometry to fee the statistial vetting software vespa (Morton 212, 215). Mayo et al. (218), also using vespa an aitional aaptive optis observations, again etete the three transiting objets an etermine their probabilities of being an elipsing binary, bakgroun elipsing binary, or hierarhial elipsing binary eah to be < 1 4. Most important, all these previous etetions of K2-32 b,, an were ahieve in searhes for box-like transit signals. Vanerburg et al. (216) use BLS, Crossfiel et al. (216) use the TERRA software (Petigura et al. 213a), whih inlues a searh for box-like transit signals very muh like BLS, an Mayo et al. (218) also use BLS. K2-32 b,, an have also been onfirme via stellar raial veloity measurements by Dai et al. (216) an Petigura et al. (217). Artile number, page 1 of 8

2 A&A aepte Raw flux (normalize) Detrene flux (ppm) (a) (b) Time (ays) Fig. 1. K2 light urve of K2-32. (a) After orretion for systemati effets with EVEREST. The re line shows our running meian filter. (b) Detrene light urve obtaine by iviing the EVEREST light urve by the running meian. Transits etete with TLS are highlighte with green (K2-32 b), magenta (K2-32 ), blue (K2-32 ), an strethe re (K2-32 e) vertial bars. 2. Methos 2.1. Target seletion In this paper, we move on from the testing phase of TLS (Hippke & Heller 219) an apply TLS to real light urves. In this initial phase of the TLS survey, we restrit ourselves to multi-planet systems from the Kepler mission sine they have been shown to exhibit extremely low false positive probabilities (FPP) (Lissauer et al. 212). Among the systems that we stuie, K2-32 with its three previously known planets stoo out with our etetion of a highly signifiant fourth transit signal, the methoology of whih we esribe in the following Transit searh The light urve of K2-32 from ampaign 2 ontains 77.4 of almost uninterrupte observations with 3527 useful exposures at a aene of 3 min. We ignore the first 64 exposures in the light urve as they are affete by strong systematis. We remove outliers, efine as ata points > 3 σ above the running mean, from the publily available 2 K2 light urve of K2-32, whih has been orrete for instrumental effets with EVEREST (Luger et al. 216, 218) (Fig. 1a). We then remove stellar variability an other trens using a meian filter with a winow size of 1. The resulting etrene light urve is shown in Fig. 1(b). We teste other winow sizes an foun that a with of 1 offers the best ompromise between both suffiient removal of unwante variability an preservation of transit signals with urations of up to several hours. Kepler s thir law of motion (Kepler 1619) preits that planets with orbital perios < 8 aroun sun-like stars have transit urations shorter than about 8 hr (see Fig. 5 in Hippke & Heller 219). For K2-32 in partiular, the maximum uration of a planet with a single transit in the 8 K2 light urve on a irular orbit is 5.7 hr, an for a planet to exhibit two transits (hene P < 4 ) the maximum transit u- 2 hlsp_everest_k2_ll_ _kepler_v2._l.fits ration is 4.8 hr, suggesting that a 1 with for our meian filter harly affets physially plausible transit signals. We applie the publily available 3 python implementation of TLS (Hippke & Heller 219) using the stellar limb arkening, mass, an raius estimates available from the EPIC atalog (Huber et al. 216). We use TLS version in its efault parameterization but set the maximum trial perio equal to the length of the light urve in orer to searh for possible single transits Markov Chain Monte Carlo analysis of the transits To refine the planetary parameters we utilize the Markov Chain Monte Carlo (MCMC) sampler emee (Foreman-Makey et al. 213). As inputs to emee, we provie the times of the mipoints of the first transit (T ) an the orbital perio (P) obtaine with TLS for eah planet. P, T, the planet-to-star raius ratio (R p /R s ), an the transit impat parameter (b) serve as moel parameters for eah planet, while the stellar ensity (ρ s ) an two limb arkening oeffiients for a quarati limb arkening law (Kipping 213) were global parameters for all transits. We ran the MCMC analysis with 1 walkers with eah walker performing 2, steps. The first half of eah walk was isare to ensure that we preserve only burne-in MCMC hains. 3. Results 3.1. Transit etetion an haraterization In Fig. 2 we show the signal etetion effiieny (SDE TLS ) perioograms of our iterative transit searh with TLS in the K2 light urve that has been extrate with EVEREST. Eah iteration ignores the in-transit ata orresponing to transits etete in previous iterations. TLS first foun planet b (top panel) sine it is both the largest planet an it exhibits more transits than planets an in the light urve. Next, TLS foun planet (seon panel), whih is the seon largest planet in the system. Planet 3 Artile number, page 2 of 8

3 René Heller et al.: The Transit Least Squares Survey I. Disovery an valiation of K2-32 e Table 1. Charaterization of the new planet K2-32 e from MCMC moel fitting to the full set of transits in the EVEREST light urve of K2-32. Planet e Planet b Planet Planet ML a meian +err. err. ML a meian +err. err. ML a meian +err. err. ML a meian err. +err. P [] T [] b R p /R s b Notes. (a) ML = value with maximum likelihoo. (b) T = BKJD 265 with the Baryentri Kepler Julian Day BKJD = BJD - 2,454,833.. Table 2. System parameterization of K2-32 an its planets an omparison to literature values. Star Planet e Planet b Planet Planet Ref. P [] C D P M H19 R 5.62 R 3.32 R 3.77 R C16 a R R R R D R R R R P R R R R M R R R R H19 b M M M < 8.1 M (95 % onf.) < 35. M (95 % onf.) D M M < 12.1 M (95 % onf.) M P17 a [AU] C16 a P H19 FPP C16 a < < S16 < < < M H19 Notes. C16: Crossfiel et al. (216); D16: Dai et al. (216); S16: Sinukoff et al. (216); P17: Petigura et al. (217); M18: Mayo et al. (218); H19: this stuy. (a) Base on stellar lassifiation (.92 ±.7 R,.88 ±.3 M ) by Huber et al. (216). (b) Planetary raii an 68.3 % onfiene intervals base on our measure planet-to-star raius ratio an stellar raius estimates from P17. () Semimajor axes an 68.3 % onfiene intervals base on our measurements of the orbital perios an stellar mass estimates from P17. () FPP = false positive probability. Values o not ontain the multipliity boost (reution in FPP) but are etaile in S16. was etete in our thir iteration with TLS (thir panel). In a fourth run then, we foun another strong peak (SDE TLS = 26.1) at a perio of , whih we preliminary referre to as a aniate ubbe K2-32 e. We also exeute a fifth searh after masking out all transits from the four planets, but we i not fin any other signifiant signals. Figure 3 illustrates the phase-fole light urves of K2-32 e (top panel) an of the previously isovere planets b,, an from top to bottom. The shallow transit epth of about 2 parts per million (ppm) is immeiately suggestive of an Earth-size planet, given that K2-32 is a somewhat subsolar-size K warf star. Note that the orinate in the top panel, whih shows the transit ata of the new objet, is an 18-fol zoom ompare to the bottom panel ontaining the ata of planets b,, an. To further haraterize the new planet aniate, we applie the MCMC sampler to the entire light urve. The stellar ensity was foun to be times the solar ensity an the best-fit limb arkening parameters are q 1 = an q 2 = Our results for K2-32 e are shown in Table 1 an our results from the MCMC sampling for the previously known planets K2-32 b,, an are ompare to the literature values in Table 2. Figure 4 is a visualization of the physial an orbital harateristis. K2-32 is iniate as a roppe orange shae irle an its planets are enote by empty irles with raii to sale Artile number, page 3 of 8

4 A&A aepte K2-32 b 2 e 4 2 Flux (ppm) Flux (ppm) Time from mi-transit (hrs) b e Fig. 3. Phase-fole transit light urves of all planets. Dots show the K2 ata, lines represent our best-fit MCMC moels. Note the ifferent sales of the orinates in the top an bottom panels. In the top panel, the transit ip of K2-32 e is shown together with a sliing mean (ashe line) of 11 aenes in with. K2-32 e b Perio (ays) Fig. 2. Signal etetion effiienies of our suessive transit searh after masking out previously etete transits. with the star. The vertial bars enote the unertainties in b, an the transit urves are to-sale representations of the best-fit transit moels shown in Fig. 3. Orbital resonanes are iniate at the bottom of this hart, with the 1:1 resonane referring to the innermost planet, K2-32 e. Interestingly, we foun that this fourplanet system is lose to being in a 1:2:5:7 mean motion resonane (MMR) hain. At the same time, however, the system is learly out-of this resonane ompare to our measurement unertainties False positive vetting an valiation We use the publily available vespa software (Morton 212, 215) to evaluate the FPP of K2-32 e. In brief, vespa takes the phase-fole transit light urve together with the elestial oorinates an stellar parameters to alulate the probabilities of the ata being ause by non-assoiate blene elipsing binaries, hierarhial triples, genuine elipsing binaries, an non-assoiate stars with transiting planets. We supplie vespa with the elestial oorinates of K2-32, the orbital perio an planet-to-star raius ratio of our aniate signal (Table 1), the phase-fole transit light urve (Fig 3, top panel), the stellar parameters as etermine spetrosopially by Petigura et al. (217) (Table 2), an the 2MASS broaban photometry (J = 1.44 ±.24, H = ±.25, K = ±.19; Cutri et al. 23) following the Lessons Learne setion in 1:1 1:2 1:5 1:7 Fig. 4. System arhiteture. Stellar an planetary raii are to sale. Planetary istanes to the star are mutually to sale but not with respet to the raii. The shapes of the transit light urves are to sale as well. Orbital resonanes are iniate with respet to the innermost planet K2-32 e. The error bars enote our unertainties in the transit impat parameter (Table 1). Petigura et al. (217). As a limiting aperture within whih the transits are observe, we referre to aaptive optis an K2 light urve analyses of Sinukoff et al. (216), who foun that the transits are loalize within 8 arseons of K2-32. vespa returne an FPP of for K2-32 e. Although this number formally valiates K2-32 e as a planet alreay (the ommonly use FFP threshol for valiation is 1 %), it still oes not onsier the multipliity boost (reution in FPP) inherent to planet aniates in multi-planet systems. Planet aniates in systems known to alreay harbor planets have a muh lower FPP than single aniates (Lissauer et al. 212). We use the values of the homogeneous K2 exoplanet survey by Vanerburg et al. (216) with n t = 59,174 as the number of target stars (K2 ampaigns to 3), n = 234 as the number of K2 planet aniates, an n m = 26 as the number of aniate multi-planet systems an plugge them into Eq. (6) of Lissauer et al. (212). We foun that the total number (the expetation value) of K2 systems like K2-32 with two or more known plan- Artile number, page 4 of 8

5 René Heller et al.: The Transit Least Squares Survey I. Disovery an valiation of K2-32 e BLS K2-32 b TLS K2-32 b e 6 e Perio (ays) Perio (ays) Fig. 5. Iterative transit searh using K2SFF ata with BLS (left) an TLS (right). Planets K2-32 b,, an are suessively etete with strong signals (see labels) with both BLS an TLS. Using a ommon etetion threshol value of 9, however, planet e is not etete with BLS (SDE BLS = 8.9) but it is etete with TLS (SDE TLS = 13.2). The major peak near P =.4 with (SDE BLS = 1.5) in the fourth iteration of BLS is a false positive (blue marker). It overshaows the true transit signal near P = 4.3 with (SDE BLS = 8.9). ets an one aitional false positive ranges between.1, assuming a true positive rate (or planet fielity) of 9 %, an.5, assuming a true positive rate of 5 %. These estimates inrease our onfiene in K2-32 e being a true planet beause both values are muh smaller than 1. Sinukoff et al. (216) showe that, base on the statistial framework of Lissauer et al. (212), the multipliity boost is at least an orer of magnitue in FPP. As a onsequene, a onservative estimate of the FPP for K2-32 e inluing the multipliity boost is < With the epth of the seonary transit (δ 2 ) being of the orer of (R p /a) 2, or.56 ppm in the ase of K2-32 e, we expet this phenomenon to be invisible to K2. Thus, as an aitional test for a false positive, for example, ause by an elipsing binary, we measure δ 2 in that part of the light urve where the seonary transit an be expete on a irular orbit. We etete a 8 ± 8 ppm ip, that is to say, a 1 σ signal that is statistially ompatible with noise Comparison of the BLS an TLS performanes One of the most pressing questions is: why has K2-32 e not been etete in previous transit searhes, e.g. the ones by Vanerburg et al. (216) an Crossfiel et al. (216)? For our omparison of the transit searh results obtaine with BLS an TLS we have to keep in min that the efinitions of the signal etetion effiieny introue by Kovás et al. (22) for BLS (SDE BLS ) an by Hippke & Heller (219) for TLS (SDE TLS ) are ifferent. 4 While SDE BLS is more losely relate to the mean transit epth, SDE TLS is erive from a χ 2 statisti an therefore naturally tens to proue higher values for a given transit-like signal. This prinipal ifferene between SDE BLS an SDE TLS must also be note in the ontext of the heuristially hosen etetion threshols of typially between 8 an 9 in other stuies. We have shown in Hippke & Heller (219) using injete transits of Earth-like planets in the light urves of Sun-like stars with white noise that a false positive rate of 1 % is ahieve at virtually the same signal etetion effiieny value of 7 for both BLS (here SDE BLS = 7) an TLS (here SDE TLS = 7). Although the statis of BLS an TLS are ifferent, they proue the same SDE values for this partiular set of simulate ata. Most important, however, the true positive rates were 76 % for BLS an 93 % for TLS. The TERRA algorithm uses yet another searh statisti, the signal-to-noise ratio (S/N) with a etetion threshol set to a nominal value of 12 by trial an error (Petigura et al. 213b). 4 See Eqs. (5) an (6) in Kovás et al. 22 an Eqs. (3) an (5) in Hippke & Heller 219, respetively. Artile number, page 5 of 8

6 A&A aepte Both Vanerburg et al. (216) an Crossfiel et al. (216) use the K2 light urve proue with the K2 self-flat-fieling (K2SFF) pipeline, whih removes most of the instrumental jitter (Vanerburg & Johnson 214). This light urve ontains a somewhat stronger noise ompare to the light urve extrate with the EVEREST pipeline an so we nee to examine whether our isovery of K2-32 e is owing to improvements of the K2 light urve extration an systemati etrening or ue to the use of TLS instea of BLS. We thus repeate our analysis using the K2SFF ata an teste both TLS an BLS, the results of whih are shown in Fig. 5. In the left olumn of Fig. 5, we present an iterative transit searh in the K2SFF light of K2-32 with BLS. We suessively etete planets b,, an with very strong SDE BLS signals (see labels) whereas planet e proue an SDE BLS of just 8.9 in the fourth iteration, whih is below the wiely use etetion threshol value of 9 (Vanerburg et al. 216; Crossfiel et al. 216). We verifie that the aition of a smoothing filter to the SDE BLS spetrum, akin to what is one in TLS to the SDE TLS spetrum, oes not lift this signal over the etetion threshol. In fat, none of the previous K2 transit surveys by Vanerburg et al. (216), Crossfiel et al. (216), Petigura et al. (218), an Mayo et al. (218) atually applie suh a smoothing filter to their SDE BLS or S/N spetra. Moreover, the signature of K2-32 e is only the seon strongest signal, while the strongest signal near P =.4 with an SDE BLS of 1.5 is a false positive, as we verifie. We foun that this false positive signal is proue by a ombination of two harateristis of BLS; the quantifiation in transit urations in units of aenes an the linear perio gri typially use in BLS appliations. For example, Petigura et al. (213b) use transit urations that are integer numbers of long aene measurements (...) with a gri of e.g., T = [3, 5, 7, 1, 14, 18] long aene measurements. Aliasing effets of uration an perio integer multiples ause aitional jitter in the BLS spetrum, an effet that is absent in TLS sine it oes not apply any binning. In aition, the perio gri in TLS is not linear (Ofir 214) beause the frequeny gri is linear, whih intrinsially suppresses the possibility of aliasing. Hene, we are able to reproue the results of previous searhes with the traitional BLS algorithm an onfirm that K2-32 e remains unetete in the K2SFF light urve. In the right olumn of Fig. 5, we show the SDE TLS results for an iterative transit searh with TLS in the K2SFF ata. We foun planets b,, an with SDE TLS values omparable to the ones obtaine with BLS, but most interestingly planet e also passe the etetion limit with a signifiant SDE TLS of This is a fasinating real-worl example of the simulation-base finings by Hippke & Heller (219) that the strongest improvement in using TLS rather than BLS ours for shallow transits or, in this ase, for Earth-size planets. We also foun that planet e an be reovere in the EVEREST ata with both BLS (SDE BLS = 21.3) an TLS (SDE TLS = 26.1). We onlue that if either the TLS searh algorithm or the EVEREST ata reution pipeline woul have been available to Vanerburg et al. (216) an Crossfiel et al. (216), they oul have foun K2-32 e. Most important though, the ombination of EVEREST an TLS leverages the optimal etrening an signal etetion proeures. We also examine the ifferent methos that were use to remove longer-term stellar trens. Both Vanerburg et al. (216) an Crossfiel et al. (216) use spline fits, while we use a running meian by efault. When we use a spline fit instea of a running meian the SDE TLS values obtaine with TLS were virtually iential with variations of about Disussion 4.1. The K2-32 system For now, the mass of K2-32 e remains unonstraine. With an Earth raius in size, a first guess for the mass of K2-32 e is about an Earth mass. Although planets b,, an have ensities ranging between those of Saturn an Neptune, suggesting large an massive gaseous envelopes, it an safely be assume that K2-32 e oes not arry large amounts of gas. If K2-32 e onsiste mostly of gas, its low ensity woul imply a sub- Earth-mass planet. Suh a very-low-mass (an very-low-gravity) planet, however, oul harly hol on to a signifiant gas envelope, in partiular uner the effets of extreme stellar irraiation. Assuming an Earth-mass (M ) for K2-32 e, a nominal stellar mass of.856 M, an a semimajor axis of.4951 AU (Table 2), we estimate a stellar RV amplitue of.44 m s 1, whih seems just out of reah of moern spetrographs (Dai et al. 216; Petigura et al. 217). Even if K2-32 e were mae up fully of iron, its mass of then roughly 1.4 M woul imply an RV amplitue of just about.6 m s 1 an harly surpass the RV noise (or RV jitter ) of several m s 1 that is typial for solar type main-sequene stars (Wright 25) We also stuie if transit timing variations (TTVs) oul be an alternative avenue to etermine the planetary masses aroun K2-32 (Holman & Murray 25; Agol et al. 25). Using our MCMC haraterization of the system (Table 1), a mass of 16.5 M for planet K2-32 b (Petigura et al. 217), an a nominal Earth mass for K2-32 e, we fin that Eq. (8) from Lithwik et al. (212) preits a TTV amplitue of 14 s or.4 hr for K2-32 e. This is probably unetetable in the available Kepler ata (see top panel of Fig. 3) but oul be teste in a etaile follow-up stuy. The near 1:2:5:7 MMR hain of K2-32 e that we report is somewhat reminisent of other MMR hains, suh as the 3:4:6:8 MMR hain observe in the Kepler-223 system (Mills et al. 216). While the Kepler-223 system exhibits a very preise MMR tuning that has been interprete as a footprint of planet migration, the fat that K2-32 is substantially non-synhronize with the 1:2:5:7 MMR hain suggests that aitional proesses have been at work. The origin of these eviations from preise ommensurability might be in long-term star-planet tial interation that was only overrule by the MMR-reating effets of the protoplanetary isk in the very early stages of the system (Papaloizou 211; Heller 218). K2-32 joins the family of planetary systems from the Kepler mission that are just wie of exat ommensurability (Lee et al. 213), although K2-32 represents a partiular ase with four rather than just two near-resonant planets. Lissauer et al. (211) efine a new variable ( 1 ζ 1 = 3 P 1 Roun( 1 ) ) (1) P 1 as a measure of the ifferene between an observe first-orer MMR perio ratio, where P = P o /P i is the orbital perio ratio between the outer an the inner planet. These authors foun statistially signifiant eviations of the observe ζ 1 istribution from a ranom perio ratio istribution for about 1 transiting exoplanets from Kepler known at the time. In fat, their results showe that for a ranomly rawn pair of Kepler planets Artile number, page 6 of 8

7 René Heller et al.: The Transit Least Squares Survey I. Disovery an valiation of K2-32 e Perio (ays) Perio (ays) Fig. 6. Left: Zoom into the bottom left panel of Fig. 5. The major peak near P =.4 (blue marker) in the SDE BLS spetrum is a false positive that is ause by an aliasing effet of the isrete aene-base uration gri employe by most BLS implementations, here using 1, trial perios on a linear gri between.3 an 5. Right: Same as in the left panel but now using more than ten times as many trial perios an a non-linear perio gri with the spaing inreasing as P 1/3. Although the false positive peak near P =.4 has isappeare, the signal of K2-32 e near P = 4.3 is still formally a false negative with SDE BLS = 8.9. ζ is most likely to range between -.1 an -.2. Taking the meian orbital perios from our MCMC sampling for K2-32 e an b (Table 1), whih are very lose to the 2:1 MMR resonane (P o /P i = 2.676), we fin ζ 1 =.19. Although this result annot be use as a means of verifiation for K2-32 e it is in very goo agreement with a large sample of exoplanets from Kepler (Lissauer et al. 211) an, hene, serves as an iniret piee of eviene in favor of the planetary nature of K2-32 e Aliasing effets in BLS Figure 5 shows our results obtaine with the BLS referene implementation in astropy. Nevertheless, the question arises if the ourrene of a false positive near P =.4 oul be avoie with ifferent gris for the trial perios an trial urations. Hene, we varie the resolution an spaing of both the uration an the perio gris for BLS to etermine whether we an remove this alias peak. The results are illustrate in Fig. 6 an isusse in the following. First, the gri of trial transit urations use in BLS must be in multiples of a onstant aene. As an intrinsi property of BLS, any given aene is teste for one of two test flux values that relate to either the out-of-transit flux or to the in-transit flux (Kovás et al. 22). As a onsequene, when BLS is applie with a perio gri of onstant spaing, short trial perios may happen to be integer multiples of the trial uration. For example, a trial perio of.416 is suffiiently lose to the 2-fol multiple of a 3 min trial transit uration (a single long aene exposure of a K2 light urve), an a tenfol multiple of a 6 min transit uration (worth two aenes), et. As the transit uration grows with the orbital perio to the power of one thir (Eq. (1) in Hippke & Heller 219) everything else being equal, aliasing beomes less of an issue for longer trial perios. Nevertheless, aliasing annot be avoie ompletely by using an infinitely fine perio gri. Instea, it an be avoie by using a non-linear perio gri as employe by TLS. Seon, the gri of trial urations teste with BLS is usually the same for all trial perios, see for example the BLS implementation in astropy. Although users an efine a maximum trial uration of, for example,.21 ays (or 1 aenes), whih is suffiiently long for most long-perio planets in K2, this trial uration is as muh as half of an orbital perio for the most shortperio planets that are physially plausible. Again, lustering of ata an generate aliases suh as those visible near P =.3 in Fig. 6 (left panel) an eventually lea to false positive etetions. In this example, the most reliable way for us to remove the alias peaks was to use a hyper fine non-linear gri of > 1, trial perios in ombination with a uration gri that uses a maximum trial uration less than half of the shortest perio. Nevertheless, an maybe most important, the signal of K2-32 e near P = 4.3 remaine a false negative for the BLS searh. In pratie, it oul be possible to evelop more appropriate trial perio an trial uration gris for BLS similar to the one implemente in TLS, but this is beyon the sope of this paper an the isavantage of the suboptimal etetor shape (box versus limb-arkene transit) woul remain. 5. Conlusions We etermine with high signifiane (SDE TLS = 26.1) a fourth sequene of perioi transits in the light urve of K2-32. Our MCMC simulations suggest that this signal has a perio of an the transiting planet has a raius of R, making it one of the smallest planets foun with K2 so far. Our false positive vetting of K2-32 e as an iniviual transiting objet yiels FPP = Fatoring in the planetary multipliity aroun K2-32, we fin FPP < , formally valiating K2-32 e as a planet. This new planet reveals a near 1:2:5:7 MMR hain of now four planets aroun K2-32. While being very lose to this MMR hain, however, the planets are in fat just wie of exat ommensurability, thereby joining a growing family of this type of systems from the Kepler mission. We fin that the offset from exat ommensurability between the new planet an the previously know K2-32 b is in very goo agreement with the offsets foun in other multi-planet transiting systems, aing more eviene in favor of the planetary nature of K2-32 e. K2-32 also joins the list of K2 systems with four or more transiting planet aniates, about a ozen of whih are known as of toay. For now, K2-138 is the only system with five (Christiansen et al. 218) or potentially even six (Haregree-Ullman & Christiansen Artile number, page 7 of 8

8 219) transiting planets, all of whih are in the super-earth to sub-neptune raius regime. Our isovery onfirms that TLS an fin sub-earth-size planets that have previously been misse with searh algorithms looking for box-like transit signals suh as BLS (Kovás et al. 22). We valiate that there are two reasons for why previous searhes have misse K2-32 e. First, we use the K2 light urve of K2-32 that was subjet to the highly effiient removal of systemati effets with EVEREST (Luger et al. 216), while previous searhes use the light urve etrene with K2SFF that has slightly less favorable noise properties. Seon, TLS is intrinsially more effiient in fining shallow transits, a fat that is simply roote in the searh funtion being a transit rather than a box (Hippke & Heller 219). Interestingly, we fin that K2-32 e oul have been etete in the K2SFF light urve with TLS. Aknowlegements. The authors thank the referee for a helpful report. This researh has mae use of the NASA Exoplanet Arhive, whih is operate by the California Institute of Tehnology, uner ontrat with the National Aeronautis an Spae Aministration uner the Exoplanet Exploration Program. This work mae use of NASA s ADS Bibliographi Servies. RH is supporte by the German spae ageny (Deutshes Zentrum für Luft- un Raumfahrt) uner PLATO Data Center grant 5OO151. KR is a member of the International Max Plank Researh Shool for Solar System Siene at the University of Göttingen. KR performe the MCMC analysis of the light urve. A&A aepte Referenes Agol, E., Steffen, J., Sari, R., & Clarkson, W. 25, MNRAS, 359, 567 Batalha, N. M., Rowe, J. F., Bryson, S. T., et al. 213, ApJS, 24, 24 Boruki, W. J., Koh, D., Basri, G., et al. 21, Siene, 327, 977 Christiansen, J. L., Crossfiel, I. J. M., Barentsen, G., et al. 218, AJ, 155, 57 Crossfiel, I. J. M., Ciari, D. R., Petigura, E. A., et al. 216, ApJS, 226, 7 Crossfiel, I. J. M., Guerrero, N., Davi, T., et al. 218, ApJS, 239, 5 Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 23, VizieR Online Data Catalog, 2246 Dai, F., Winn, J. N., Albreht, S., et al. 216, ApJ, 823, 115 Foreman-Makey, D., Hogg, D., Lang, D., & Gooman, J. 213, PASP, 125, 36 Haregree-Ullman, K. & Christiansen, J. 219, in Amerian Astronomial Soiety Meeting Abstrats, Vol. 233, Amerian Astronomial Soiety Meeting Abstrats #233, #164.7 Heller, R. 218, arxiv e-prints [arxiv: ] Hippke, M. & Heller, R. 219, A&A, 623, A39 Holman, M. J. & Murray, N. W. 25, Siene, 37, 1288 Howell, S. B., Sobek, C., Haas, M., et al. 214, PASP, 126, 398 Huber, D., Bryson, S. T., Haas, M. R., et al. 216, ApJS, 224, 2 Kepler, J. 1619, Harmonies muni libri V Kipping, D. M. 213, MNRAS, 435, 2152 Kovás, G., Zuker, S., & Mazeh, T. 22, A&A, 391, 369 Lee, M. H., Fabryky, D., & Lin, D. N. C. 213, ApJ, 774, 52 Lissauer, J. J., Mary, G. W., Rowe, J. F., et al. 212, ApJ, 75, 112 Lissauer, J. J., Ragozzine, D., Fabryky, D. C., et al. 211, ApJS, 197, 8 Lithwik, Y., Xie, J., & Wu, Y. 212, ApJ, 761, 122 Livingston, J. H., Crossfiel, I. J. M., Petigura, E. A., et al. 218a, AJ, 156, 277 Livingston, J. H., Enl, M., Dai, F., et al. 218b, AJ, 156, 78 Luger, R., Agol, E., Kruse, E., et al. 216, AJ, 152, 1 Luger, R., Kruse, E., Foreman-Makey, D., Agol, E., & Sauners, N. 218, AJ, 156, 99 Manel, K. & Agol, E. 22, ApJ, 58, L171 Manua, A., Bell, R. A., & Gustafsson, B. 1977, A&A, 61, 89 Mayo, A. W., Vanerburg, A., Latham, D. W., et al. 218, AJ, 155, 136 Mills, S. M., Fabryky, D. C., Migaszewski, C., et al. 216, Nature, 533, 59 Morton, T. D. 212, ApJ, 761, 6 Morton, T. D. 215, VESPA: False positive probabilities alulator, Astrophysis Soure Coe Library Ofir, A. 214, A&A, 561, A138 Papaloizou, J. 211, CeMDA, 111, 83 Petigura, E. A., Crossfiel, I. J. M., Isaason, H., et al. 218, AJ, 155, 21 Petigura, E. A., Howar, A. W., & Mary, G. W. 213a, PNAS, 11, Petigura, E. A., Mary, G. W., & Howar, A. W. 213b, ApJ, 77, 69 Petigura, E. A., Sinukoff, E., Lopez, E. D., et al. 217, AJ, 153, 142 Sinukoff, E., Howar, A. W., Petigura, E. A., et al. 216, ApJ, 827, 78 van Sluijs, L. & Van Eylen, V. 218, MNRAS, 474, 463 Vanerburg, A. & Johnson, J. A. 214, PASP, 126, 948 Vanerburg, A., Latham, D. W., Buhhave, L. A., et al. 216, ApJS, 222, 14 Wright, J. T. 25, PASP, 117, 657 Yu, L., Crossfiel, I. J. M., Shlieer, J. E., et al. 218, AJ, 156, 22 Artile number, page 8 of 8