## TESS discovery of two super-Earths orbiting the M-dwarf stars TOI-6002 and TOI-5713 near the radius valley

M. Ghachoui 1 2 , , B.V. Rackham 3 4 , M. Dévora-Pajares , 5 6 , , J. Chouqar 1 , M. Timmermans 1 , L. Kaltenegger 7 , D. Sebastian , F.J. Pozuelos 8 9 , J.D. Eastman 10 , A.J. Burgasser 11 , F. Murgas 12 13 , , K.G. Stassun 14 , M. Gillon 1 , Z. Benkhaldoun 2 , E. Palle 12 13 , , L. Delrez 1 15 , , J.M. Jenkins 16 , K. Barkaoui 1 3 12 , , , N. Narita 17 18 12 , , , J. P. de Leon 19 , M. Mori 18 20 , , A. Shporer 21 , P. Rowden 22 , V. Kostov 23 , G. F˝ urész 4 , K.A. Collins 10 , R.P. Schwarz 10 , D. Charbonneau 10 , N.M. Guerrero 24 25 , , G. Ricker 4 26 , , E. Jehin 15 , A. Fukui 17 12 , , Y. Kawai 19 , Y. Hayashi 19 , E. Esparza-Borges 12 13 , , H. Parviainen 13 12 , , C.A. Clark 27 28 , , D.R. Ciardi 28 , A.S. Polanski 29 , J. Schleider 23 , E.A. Gilbert 27 , I. J.M. Crossfield 29 , T. Barclay 23 , C.D. Dressing 30 , P.R. Karpoor 11 , E. Softich 11 , R. Gerasimov 31 32 , , and F. Davoudi 1

(A ffi liations can be found after the references)

September 17, 2024

## ABSTRACT

We present the validation of two TESS super-Earth candidates transiting the mid-M dwarfs TOI-6002 and TOI-5713 every 10.90 and 10.44 days, respectively. The first star (TOI-6002) is located 32 038 . ± 0 019 pc away, with a radius of 0 2409 . . + 0 0066 . -0 0065 . R , a mass of 0 2105 ⊙ . + 0 0049 . -0 0048 . M , and an ⊙ e ff ective temperature of 3229 + 77 -57 K. The second star (TOI-5713) is located 40 946 . ± 0 032 pc away, with a radius of 0 2985 . . + 0 0073 . -0 0072 . R , a mass ⊙ of 0 2653 . ± 0 0061 . M , and an e ⊙ ff ective temperature of 3225 + 41 -40 K. We validated the planets using TESS data, ground-based multi-wavelength photometry from many ground-based facilities, as well as high-resolution AO observations from Keck NIRC2. TOI-6002 b has a radius of 1 65 / . + 0 22 . -0 19 . R ⊕ and receives 1 77 . + 0 16 . -0 11 . S ⊕ . TOI-5713 b has a radius of 1 77 . + 0 13 . -0 11 . R ⊕ and receives 2 42 . ± 0 11 . S ⊕ . Both planets are located near the radius valley and near the inner edge of the habitable zone of their host stars, which makes them intriguing targets for future studies to understand the formation and evolution of small planets around M-dwarf stars.

Key words. Planets and satellites - Techniques: photometric - Methods: numerical

## 1. Introduction

The Transiting Exoplanets Survey Satellite (TESS; Ricker et al. 2015) has significantly increased the sample of exoplanets smaller than Neptune and larger than Earth (broadly called "super-Earths" and "sub-Neptunes") known to orbit nearby stars. The absence of planets in this size range in our own Solar System limits us to understand how these planets form and evolve, making them compelling targets to investigate. Although it is arguably important to focus on a wide range of stellar types to study exoplanets, M dwarfs provide appealing targets for our study of super-Earths due to their abundance in our galaxy (Henry et al. 2019; Reylé et al. 2021) and the significantly more favorable planet-to-star size ratio compared to Sun-analog stars. Especially planets at the inner edge of the Habitable Zone (HZ) provide compelling targets because of the transition from hot, potentially habitable rocky worlds to Venus-like worlds (Kaltenegger et al. 2023a). Especially planets close to the inner edge of the HZ of M dwarfs, which could arguably show extended secondary atmospheres resulting from a runaway greenhouse e ff ect, thus provide very interesting targets for observations (Kaltenegger et al. 2023b; Delrez et al. 2022; Turbet et al. 2019)

The abundance and the distribution of planets with radii below 4 R ⊕ for orbits shorter than 100 days allows the identification of specific features that can link their abundance to their formation and evolution like the radius valley (see e.g. Fulton et al. 2017; Cloutier &amp; Menou 2020) and the density valley (Luque &amp; Pallé 2022), based on their radii and densities distributions.

The radius valley features a lack in planets ranging from 1.5 to 2 R . In this regard, some studies argue that rocky super-Earths ⊕ are a byproduct of atmospheric escape experienced by gaseous sub-Neptunes driven by the energy radiated by the host stars (see e.g. Lopez &amp; Fortney 2013; Owen &amp; Wu 2013; Jin et al. 2014; Chen &amp; Rogers 2016) and (or) the energy released by their cooling planetary cores (see e.g. Owen &amp; Wu 2013; Ginzburg et al. 2016, 2018). Other studies argue that some planets are formed as rocky super-Earths in a gas-depleted (see e.g. Lopez &amp; Rice 2018; Cloutier &amp; Menou 2020; Cherubim et al. 2023) or gaspoor environment with thin gaseous envelopes (see e.g. Lee &amp; Connors 2021; Cherubim et al. 2023). Additional theoretical studies predict that the sub-population of sub-Neptunes consists typically of water-rich worlds formed beyond the water ice-line and migrate toward their host stars (se e.g. Venturini et al. 2020; Burn et al. 2024; Venturini et al. 2024). The density valley features three peaks in the density distribution defining rocky superEarths, water-rich worlds, and pu ff y sub-Neptunes. Luque &amp; Pallé (2022) suggest that the pebble accretion is the main mechanism forming these planets with super-Earths formed within the ice line and water-rich world formed outside it (se e.g. Venturini et al. 2020; Burn et al. 2024; Venturini et al. 2024). For now the mechanisms forming super-Earths remain unclear, making observations of planets in this mass-radius regime a critical part to understand it.

TESS is increasing the number of small exoplanets amenable to mass and atmospheric characterization (Ricker et al. 2015; Sullivan et al. 2015). Detailed studies of such exoplanets will allow insights into the origin of the radius and density valleys in

the near future. Such analyses should shed further light on the mechanisms dominating the formation and evolution of superEarths. In addition, the analyses of exoplanets with incident irradiation between Earth and Venus will allow key insights into the limits of the HZ (e.g. Kaltenegger et al. 2023a).

This paper reports the validation of two super-Earths, TOI6002 b (1 65 . + 0 22 . -0 19 . R ) and TOI-5713 b (1 77 ⊕ . + 0 13 . -0 11 . R ) orbiting ⊕ two mid- type M-dwarf stars, TOI-6002 and TOI-5713 . They were first detected by TESS, and validated in this study using various ground-based observations (see Sect. 2). The masses of both planets are left unconstrained by our observations, but their radii suggest rocky or water-rich composition. Note that precise mass measurements should be possible with state-of-the art Doppler spectographs to derive their masses, as discussed in more detail in Sect. 4.

The paper is structured as follows. Sect. 2 introduces the data from TESS and all ground-based observations. Stellar characterization, validation of the planetary nature of the two TESS candidates, and global data analyses are covered in Sect. 3. Results, potential for mass measurements, and atmospheric characterization are discussed for both planets in Sect. 4. Sect. 5 summarizes our findings.

## 2. Observations

This section presents all observations of TOI-6002 and TOI5713 obtained from space and ground-based instruments. In Table 1, we summarize all the ground-based, photometric timeseries observations, as well as the observation parameters. We note that the di ff erential photometry was performed using uncontaminated circular apertures for all ground-based observations.

## 2.1. TESS photometry

TOI-6002 was observed by TESS in sectors 14 (18 July to 14 August 2019), 41 (24 July to 20 August 2021), and 54-55 (9 July to 1 September 2022). TOI-5713 was observed by TESS in sectors 16 (12 September to 06 October 2019), 22-23 (19 February to 15 April 2020), and 49 (26 February to 25 March 2022). Note that this target was also observed in sectors 74 and 75 from 03 January to 26 February 2024. These two sectors were not included in our analyses here (see Sect. 3.6), but the data appears consistent with our results and will be used as part of future research on these planets.

The TESS time-series observations were processed by the TESS Science Processing Operations Center (SPOC) pipeline (Jenkins et al. 2016, 2020). SPOC calibrates the image data, carries out quality control, labels bad data, and performs simple aperture photometry (SAP) for each star in the field of view. Then, it corrects for instrumental systematic o ff sets using a Presearch Data Conditioning module (PDC: Smith et al. 2012; Stumpe et al. 2012, 2014) and produces pre-search data conditioned light curves. Appendix A.1 shows the target pixel files (TPFs) and the SAP apertures used for photometry, along with the overlaid locations of nearby Gaia DR2 (Gaia Collaboration et al. 2018) sources.

The SPOC pipeline executed a transit search of the combined light curve from sectors 14-55 on 11 September 2022 for TOI6002 and from sectors 16-49 on 05 May 2020 for TOI-5713 (Jenkins 2002; Jenkins et al. 2020), yielding a threshold crossing event (TCE, a statistically significant periodic, transit-like signal on the target star) with a 10.9 day period for the first target and a TCE with a 10.44 day period for the second. Initial limb-darkened transit models were fitted (Li et al. 2019) and a suite of diagnostic tests were performed to prove or refute the planetary origin of the signals (Twicken et al. 2018) for both TCEs. The di ff erence image centroiding tests conducted for the sectors 14-55 search of TOI-6002 and the sectors 16-49 search of TOI-5713 constrained the location of the target star to within 4.6 ± 4.1 ′′ and 3.4 ± 3.2 ′′ of the transit source, respectively. The two transit signatures passed all diagnostic tests presented in the SPOC data validation reports to discard di ff erent possible false positives, mainly eclipsing binaries and blended eclipsing binaries. The TESS Science O ffi ce (TSO) evaluated the vetting information and issued an alert for TOI-6002 b on 01 December 2022 and for TOI-5713 b on 8 August 2022.

For our analyses, we downloaded the 2-minute PDC-SAP light curves from the Mikulski Archive for Space Telescopes (MAST 1 ). Using the "QUALITY" flag, we identified all the badquality data points that were removed from the data files. The year two data of TESS were subject to sky background overestimation in the SPOC pipeline for dim and (or) crowded targets, a ff ecting the transit depth and, thus, the planet radius. TOI-6002 is a crowded-field target and was subjected to this e ff ect in sector 14. We corrected the PDC-SAP light curves of sector 14 using equation 2 described in section 4.2 of the data release notes of sector 27 2 (see section 4.2, equation 2). This e ff ect was negligible for TOI-5713 due to the low level of crowding. The light curve of TOI-5713 showed some clear flare-like signals, which we removed using an upper 3 σ clipping. We used wot n ¯ a (Hippke et al. 2019) to detrend the light curves for stellar variability, utilizing a bi-weight filter with a window 2.5 times longer than each planets' transit duration. The obtained lightcurves are presented in Appendix A.2, with the locations of the transit signals highlighted with di ff erent colors.

## 2.2. Ground-based photometry

## 2.2.1. LCOGT-1m0 photometry

We used the 4096 × 4096 SINISTRO CCD camera on the 1.0-m robotic telescope of the Las Cumbres Observatory Global Telescope (LCOGT; Brown et al. 2013) network node at Teide Observatory (Tenerife, Spain) to observe two full transits of TOI6002 b on UT 2023 May 13 and 24. Table 1 provides the observation parameters. We used the standard LCOGT BANZAI pipeline (McCully et al. 2018) to calibrate the science images, and we used AstroImageJ (AIJ: Collins et al. 2017) for the photometric extraction.

## 2.2.2. LCOGT MuSCAT3 photometry

We used the MuSCAT3 multi-band imager (Narita et al. 2020) on the LCOGT 2m Faulkes Telescope North at Haleakala Observatory on Maui (Hawaii) to observe one full transit of TOI6002 b on UTC August 19, 2023 (see Table 1 for more details). The multi-band observation covered g r , , i , and z s bands simultaneously. The data reduction and aperture photometry were performed in a similar way as described in the previous section.

1 https://archive.stsci.edu/

2

https://tasoc.dk/docs/release\_notes/tess\_sector\_27\_

drn38\_v02.pdf

Table 1. Ground-based facilities and their corresponding parameters used in the time-series photometric observations of TOI-6002 and TOI-5713.

.

| Date (UT)      | Filter              | Facility               | Pixel Scale ( ′′ / pixel)   | Exposure Time [s]   | Number of Exposures   | Airmass Range   | Aperture Size ( ′′ )   | Detrending Parameters               |
|----------------|---------------------|------------------------|-----------------------------|---------------------|-----------------------|-----------------|------------------------|-------------------------------------|
| TOI-6002       |                     |                        |                             |                     |                       |                 |                        |                                     |
| 13 May 2023    | i ′                 | LCO-SINISTRO (1.0 m)   | 0.389                       | 86                  | 93                    | [1.00,1.40]     | 2.33                   | Airmass + Background                |
| 24 May 2023    | i ′                 | LCO-SINISTRO (1.0 m)   | 0.389                       | 86                  | 83                    | [1.04,1.89]     | 2.33                   | Airmass + FWHM                      |
| 19 August 2023 | r ′ , z s g ′ , i ′ | LCOGT MuSCAT3 (2.0 m)  | 0.265 ...                   | 46, 32 207, 28      | 251, 339 35, 375      | [1.03,1.30] ... | 2.38 ...               | Airmass + FWHM Airmass + Background |
| TOI-5713       |                     |                        |                             |                     |                       |                 |                        |                                     |
| 17 June 2023   | I + z               | TRAPPIST-North (0.6 m) | 0.640                       | 50                  | 209                   | [1.25,2.88]     | 5.40                   | Airmass + FWHM                      |
| 08 July 2023   | r ′ , z s i ′       | TCS-MuSCAT2 (1.52 m)   | 0.440 ...                   | 50, 25 15           | 195, 457 234          | [1.16,1.94] ... | 6.96 ...               | Airmass + FWHM Airmass + Background |

Aperture size column presents the optimum uncontaminated apertures used for extracting the light curves. Detrending parameters were chosen on the basis of maximizing the log-likelihood of each transit.

## 2.2.3. TRAPPIST-North photometry

We used the CCD camera on TRAPPIST-North, a 0.6-m robotic telescope at Oukaimeden Observatory in Morocco (Jehin et al. 2011; Gillon et al. 2013; Barkaoui et al. 2019; Ghachoui et al. 2023), to observe a full transit of TOI-5713 b on UTC June 17, 2023 (see Table 1). We observed with the non-standard I + z ′ filter. The data reduction and photometry were performed using AstroImageJ .

## 2.2.4. MuSCAT2 photometry

For TOI-5713 b, we observed one full transit on UTC July 8, 2023 (see Table 1), using the multi-color imager ( g , r , i , and z s bands) MuSCAT2 (Narita et al. 2019), mounted on the 1.52-m Telescopio Carlos Sánchez (TCS) at Teide Observatory, Spain. We used the g band images for guiding, as the star is faint in the blue optical ( B ∼ 17 mag). We calibrated and reduced the raw data using the code and procedure of Parviainen et al. (2019). Due to the short exposure time and the faintness of the target, the MuSCAT2 g band photometry presented a much larger scatter than in the other bands, and it was not used in the analysis presented below.

We also observed a full transit of TOI-6002 b with MuSCAT2 on UTC October 01, 2023. However, this observation was a ff ected by clouds around the mid-transit time, and thus was omitted in the global analyses.

## 2.3. Spectroscopic observations

## 2.3.1. IRTF/SpeX

We used the SpeX spectrograph (Rayner et al. 2003) on the 3.2m NASA Infrared Telescope Facility (IRTF) to collect a nearinfrared spectrum of TOI-5713 on 2023 May 16 (UT) and TOI6002 on 2023 Oct 31 (UT). On both nights, conditions were clear with seeing of 0 ′′ . 6-0 7. We used the short-wavelength cross-′′ . dispersed (SXD) mode, and aligned the 0 ′′ . 3 × 15 ′′ slit to the parallactic angle. The resulting spectra span 0.80-2.42 µ mand have a resolving power of R ∼ 2000. We collected six integrations of 139.9 s on TOI-5713 and 299.8 s on TOI-6002, using an ABBA nod pattern. Immediately after the science observations, we gathered a standard set of flat-field and arc-lamp calibrations and observed A0 V standard stars. We reduced the data with Spextool v4.1 (Cushing et al. 2004), following the standard instructions from the Spextool User's Manual 3 . The final spectrum of TOI5713 has a median per-pixel S N of 92 with peak values in the / J , H , and K bands of of 116, 131, and 121, respectively. For TOI-6002, the final spectrum has a median S N per pixel of 156 / with peak values of 195, 219, and 205 in the J , H , and K bands, respectively.

## 2.3.2. Shane telescope/kast double spectrograph observations

We observed TOI-5713 and TOI-6002 with the Kast Double Spectrograph (Miller &amp; Stone 1994) on the 3-m Shane Telescope at Lick Observatory. TOI-5713 was observed on 2023 July 7 (UT) in clear conditions with 1.1 ′′ seeing; TOI-6002 was observed on 2023 November 5 (UT) in partly cloudy conditions with 1 ′′ . 2 seeing. We used both the red and blue channels of Kast split by the D57 dichroic at 5700 Å and the 1.5 ′′ -wide slit for both observations. In the red channel, we used the 600 7500 / grating providing 5900-9000 Å wavelength coverage at an average resolution of λ/ ∆ λ = 1450. In the blue channel, we used the 600 4310 grism providing 3600-5600 Å wavelength cover-/ age at an average resolution of λ/ ∆ λ = 1100. For TOI-5713, we obtained one blue exposure of 1200 s and two red exposures of 600 s each at an average airmass of 1.07; for TOI-6002, we obtained one blue exposure of 900 s and two red exposures of 450 s each at an average airmass of 1.03. We observed a nearby G2 V star immediately after each source for calibration of telluric absorption, and spectrophotomeric flux standards PG 1708 + 60 and BD + 28 4211 (Hamuy et al. 1992, 1994) during each night for flux calibration. Flat-field lamps, HeNeAr arc lamps, and HgCd arc lamps were also observed at the start of each night for pixel response and wavelength calibration.

All data were reduced using the kastredux package , which 4 included image calibration (pixel response correction, conversion to electrons second, bad pixel masking), boxcar extraction / with median background subtraction, wavelength calibration, flux calibration, and correction of telluric absorption in the red channel. The final spectra have median signals-to-noise of 140 in both red channels (near 7450 Å) and 45 in the blue channel of TOI-6002 (near 5100 Å). Due to calibration issues, we discarded the blue observations of TOI-5713.

3

Available at

http://irtfweb.ifa.hawaii.edu/~spex/

observer/

.

4 kastredux: https://github.com/aburgasser/kastredux .

To constrain possible contamination by bound or line-of-sight companions (Ciardi et al. 2015), We observed TOI-5713 and TOI-6002 with near-infrared adaptive optics (AO) imaging at the Keck (2023-Jun-10 UT) and at the Palomar (2023-Jun-08 UT) Observatory.

TOI-5713 was observed with the NIRC2 instrument on Keck-II (10 m) using the natural guide star AO system (Wizinowich et al. 2000) utilizing a standard 3-point dither pattern. Each position in the dither pattern is o ff set from the previous position by 3 ′′ . The 3-point dither pattern was performed 3 times, each pattern o ff set from the previous pattern by 0 5 . ′′ , resulting in a total of 9 on-target exposures. NIRC2 was used in the narrowangle mode with the full field of view ( ∼ 10 ′′ ) and a pixel scale of roughly 0 0099442 . ′′ per pixel. Observations were made in the Kcont filter ( λ o = 2 2706; . ∆ λ = 0 0296 . µ m) with an integration time of 15 seconds per dither position, yielding a total of 135 seconds on-source time.

We observed TOI-6002 with the PHARO instrument (Hayward et al. 2001) on the Palomar Hale (5m) using the P3K natural guide star AO system (Dekany et al. 2013) in the narrowband Brγ filter ( λ o = 2 1686; . ∆ λ = 0 0326 . µ m). The PHARO pixel scale is 0 025 . ′′ per pixel. We twice repeated a standard 5-point quincunx dither pattern with steps of 5 ′′ , separating each repeat by 0.5 ′′ . We used a single frame integration time of 15 sec, yielding a total integration of 225 seconds.

Flat fields were taken on-sky, dark-subtracted, and median averaged. We generated sky frames using the median average of the dithered science frames. We then sky-subtracted and flatfielded each science frame and combined the calibrated science frames into a single mosaiced image. The final resolutions of the combined dithers, determined from the full-width halfmaximum (FWHM) of the point spread functions, were 0.057 ′′ and 0.108 ′′ for TOI-5713 and TOI-6002, respectively.

We assessed the sensitivities of the final combined AO image by injecting simulated sources azimuthally around the primary target, following the approach of Furlan et al. (2017). Fig. 1 shows the resulting sensitivity curves. We do not detect close-in ( ≲ 1 ′′ ) stellar companions for either star. The Keck and Palomar data have sensitivities of δ mag ≈ 7 5 mag and 6 9 mag at . . 0.5 ′′ , respectively. This corresponds to a spatial limit ≲ 32 au for TOI-6002 and ≲ 40 94 au for TOI-5713. .

## 2.5. Gaia assessment

In addition to high-resolution imaging, we used Gaia to pinpoint any widely separated, bound stellar companions. The Gaia catalog does not identify any nearby sources that have the same distance and proper motion as TOI-5713 or TOI-6002 (see also Mugrauer &amp; Michel 2020, 2021). There is a 4 ′′ neighbor source to TOI-6002, but it is a line-of-sight source at &gt; 2000 pc and unrelated to TOI-6002.

We also used the Gaia renormalized unit weight error (RUWE) to assess the possibility of any close-in companions. RUWE values ≳ 1 4 indicate excess astrometric noise, which . can potentially be caused by a bound companion (Ziegler et al. 2020). The Gaia DR3 RUWE values for TOI-5713 and TOI6002 are 1.3 and 1.1, respectively, indicating astrometric solutions that are consistent with single-star models.

## 3. Analyses

## 3.1. Stellar characterization

## 3.1.1. Spectroscopic analysis

Fig. 2 shows the SpeX SXD spectra of TOI-5713 and TOI-6002. To assign spectral types, we compared the spectra to single-star spectral standards in the IRTF Spectral Library (Cushing et al. 2005; Rayner et al. 2009) using the SpeX Prism Library Analysis Toolkit (SPLAT, Burgasser &amp; Splat Development Team 2017). For both stars, we find the single best match to the M3.5 dwarf GJ273, with M3 and M4 dwarf standards providing the next closest matches. Accordingly, we adopt an infrared spectral type of M3.5 ± 0.5 for both targets. Following Delrez et al. (2022), we used the Mann et al. (2013a) relation between the equivalent widths of the K -band Na i and Ca i doublets, the H2O-K2 index (Rojas-Ayala et al. 2012), and stellar metallicity to estimate an iron abundance of [Fe / H] = -0 03 . ± 0 13 for TOI-5713 and . [Fe H] / = -0 18 . ± 0 12 for TOI-6002. .

The Shane Kast optical spectra of TOI-5713 and TOI-6002 / are shown in Fig. 3. Comparison to SDSS M dwarf templates from Bochanski et al. (2007) shows a best-match to the M4 template in both cases, with index-based classifications from Reid et al. (1995); Gizis (1997); and Lépine et al. (2003) indicating equivalent classifications of M3.5-M4, also consistent with the near-infrared classifications. We detect H α emission in TOI-5713 with an equivalent width of 1.94 ± 0.06 Å, corresponding to log 10 LH α/ Lbol = -4.16 ± 0.07 based on the χ factor relation of Douglas et al. (2014). For TOI-6002, we see no evidence of H I emission, with an H α equivalent width limit of 0.6 Å. These emission signatures indicate ages less than and greater than ∼ 4.5 Gyr, respectively (West et al. 2008). We measured metallicity-sensitive ζ indices (Lépine et al. 2007; Mann et al. 2013b) of 1.07 ± 0.01 and 1.02 ± 0.01 for TOI-5713 and TOI-6002, respectively, consistent with near-solar metallicities of [Fe H] / = + 0.08 ± 0.20 and [Fe H] / = + 0.04 ± 0.20 (Mann et al. 2013b), again formally consistent with the metallicities inferred from the near-infrared spectra.

## 3.1.2. SED and empirical relations within the global modeling

We performed a Spectral Energy Distribution (SED) fit jointly with transit modeling (see Sect. 3.6) using the latest version of EXOFASTv2 5 (Eastman et al. 2019) software package. The latest version of EXOFASTv2 comes with the possibility of simultaneously utilizing the empirical relations of Mann et al. (2015) to determine the stellar radius R ∗ , and of Mann et al. (2019) to determine the stellar mass M ∗ . We activated this feature using the two corresponding flags, MANRAD and MANMASS . The R ∗ and M ∗ values obtained by EXOFASTv2 are optimized and self-consistent compromises between the constraints from the SED, empirical relations, and transits (via stellar density). The stellar parameters obtained in this case are more precise than when using the MIST stellar evolution models (Dotter 2016), which was disabled by setting the NOMIST flag as recommended in Eastman et al. (2019). For the SED fit, we used the JHKS magnitudes from 2MASS , the W1-W4 magnitudes from WISE (Cutri et al. 2013), the G BP G RP magnitudes from Gaia DR3 (Gaia Collaboration 2020), and the Pan-STARRS (Chambers et al. 2016) yz magnitudes (see Table 2). We applied a Gaussian prior on the Gaia EDR3parallax, which we corrected for systematics by subtracting -0.044991895 mas from the nominal value according to

5 https://github.com/jdeast/EXOFASTv2/branches

Fig. 1. NIR AO imaging and sensitivity curve for TOI 5713 (Keck) and TOI 6002 (Palomar). Inset: Images of the central portion of the data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-4b10bd92-684f-4fe1-a835-f932ee8c684c.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-1bf3b5e0-4d2e-4df2-a6d5-af8380587098.jpg)

Fig. 2. SpeX SXD spectra of TOI-5713 (top) and TOI-6002 (bottom). Spectra of the targets (blue) and the M3.5V spectral standard GJ 273 (grey) are shown. Annotations indicate strong spectral features of M dwarfs, and shaded regions cover wavelengths with strong telluric absorption.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-b555438a-2fc1-4fb4-b356-51547b5dae06.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-974a5d8b-ea1a-4d19-98ad-2e60158b57f1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-0447894f-ef93-48a8-a343-c7671253f658.jpg)

Lindegren et al. (2021). We used an upper limit of AV = 8 1537 . for TOI-6002 and AV = 0 0477 for TOI-5713 on the extinctions . from the dust maps of Schlafly &amp; Finkbeiner (2011) and a Gaus-

Fig. 3. Shane Kast blue and red optical spectra (black lines) of TOI-/ 5713 (top) and TOI-6002 (bottom) compared to their best-fit M4 SDSS standard templates (magenta line; Bochanski et al. 2007). Data are normalized at 7400 Å, and the gap between 5600 Å and 5900 Å in the spectrum of TOI-6002 corresponds to the dichroic split between the Kast blue and red channels. Key spectral features across the 4200-8900 Å region are labeled. Inset boxes highlight the regions around the 6563 Å H α and 6708 Å Li I lines.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-4e824bf5-9ce1-48ca-927b-0b7c227e8864.jpg)

sian prior on the stellar metallicity from IRTF SpeX (see Table / 3). Also, a Gaussian prior was applied on the 2MASS mK band to be used by Mann et al. (2015) and Mann et al. (2019) relations. We run the fit with R ∗ , M ∗ , and the e ff ective temperature Tef f as free parameters. The SED fits are presented in Fig. 4, and results are reported in Table 3.

Fig. 4. Spectral energy distribution (SED) of TOI-5713 (left) and TOI-6002 (right). Black lines represent the SED best-fit models. Blue circles represent the average broadband flux measurements based on the SED model. Horizontal error bars show the filters bandwidth, and the vertical error bars show the uncertainty on the flux measurements.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-3c4988aa-cd0d-4b0d-bd3f-4d777395e423.jpg)

## 3.1.3. Galactic kinematics

We used the Gaia DR3 data for each target, including coordinates, proper motions, parallaxes, and radial velocities (Gaia Collaboration 2020), to infer their UVW galactic velocities in the local standard of rest (LSR) frame. We used astropy.coordinates (Astropy Collaboration et al. 2013, 2018, 2022) to convert the Gaia data to UVW velocities, utilizing the correction of Schönrich et al. (2010) for the barycentric velocity relative to the LSR. The resulting velocities are summarized in Table 2. Following Equations 6 and 7 of Li &amp; Zhao (2017), we find these UVW velocities correspond to thindisk, thick-disk, and halo membership probabilities of { 63 9 . ± 1 3% 36 1 . , . ± 1 3% . , &lt; 0 1% . } and { 97 2 . ± 0 1% 2 9 . , . ± 0 1% . , &lt; 0 1% . } for TOI-6002 and TOI-5713, respectively. The likely thin-disk membership of both stars implies ages of ≲ 8 2 Gyr (Kilic et al. . 2017).

## 3.2. Discarding possible false-positive scenarios

Two primary cases of astrophysical false positives can masquerade as exoplanet transits in TESS data. The first case is background foreground eclipsing binaries (EBs), especially those in / grazing configurations. For both targets, we searched for archival imaging, retrieving images from POSS I DSS (Minkowski &amp; / Abell 1963) and POSS II DSS Lasker et al. (1996). These im-/ ages, in addition to those taken from our recent observations from MusCAT3 (for TOI-6002) and MusCAT2 (for TOI-5713) in 2023, span more than 70 years, as shown in Fig. 5. TOI6002 has a high proper motion of ∼ 344.42 mas yr / and has moved by ∼ 24.4 ′′ from 1955 to 2023. This allowed us to confirm that there was no bright star in the current position of TOI6002. On the other hand, TOI-5713 has a relatively low proper motion ( ∼ 110.24 mas yr) and has moved by / ∼ 7.7 ′′ from 1953 to 2023, making it challenging to confirm the absence of any unresolved companion in its current position. However, highresolution adaptive optics imaging from Keck-II (for TOI-5713) and Palomar (for TOI-6002) confirm that no stellar companions were detected close to ( ≲ 1 ′′ ) the targets as discussed in Sect. 2.4. Additionally, the renormalized unit weight error (RUWE) from Gaia DR3 is below 1.4 for both stars, indicating that the stars are singles (see Sect. 2.5).

The second astrophysical false positive case is that of bound eclipsing binaries. Many findings argue against this scenario for both stars: (1) the spectral energy distribution is fit well by single-star models (see Fig. 4); (2) the spectra collected by SpeX SXD and Shane Kast, for both stars, show no secondary spec-/ trum (see Sect. 2), and Fig. 2 and 3; (3) the achromatic nature of the transits, where their depths in all bands agree within 1 σ for each planet (see Table 5); (4) the stellar densities determined from the transits alone are in good agreement with the those determined from the the basic parameters of each star (see Table 3). Considering these arguments, we discard these false-positive scenarios.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-da32ff2b-69ce-4093-986d-444d110ad157.jpg)

1953-07-11 (possl red)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-166ac34b-8097-49c2-97a0-f4f5ee22070b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-07190c7f-69fc-489d-92ff-48b49a56e421.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-1aca2da6-57f8-407d-9e41-6328b136618d.jpg)

Fig. 5. Archival images cropped with a field of view of 1' × 1' around TOI-5713 (top row) and TOI-6002 (bottom row). Circles in red show the current position of the stars, and the circles in yellow, in the bottom row, show the positions of TOI-6002 at the dates of exposures.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-ee41688f-22bb-4614-835a-e743e6ff47ec.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-621e7d72-2a65-4bc5-8da5-86c5ca00bec7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-2f332ffa-3508-4022-8aa7-75818b101c83.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-54c210f4-2480-4ee4-a422-74c44fa99979.jpg)

## 3.3. Statistical validation

In addition to ruling out false-positives via the arguments in the previous section, we used the Tool for Rating Interesting Candidate Exoplanets and Reliability Analysis of Transits Originating from Proximate Stars ( TRICERATOPS: Giacalone et al. 2021) to statistically validate TOI-6002 b and TOI-5713 b. TRICERATOPS simulates potential astrophysical false positives due to bound stellar companions, nearby stars, and line-of-sight stars. To assess the planetary nature of the transit signal, it calculates two quantities: the false positive probability (FPP), which quantifies the probability that the signal originates from other

Table 2. Stellar astrometric and photometric properties of TOI-6002 and TOI-5713.

| Parameter                    | TOI-6002            | TOI-5713            | Source   |
|------------------------------|---------------------|---------------------|----------|
| TIC                          | 102734241           | 219041246           | 1        |
| 2MASS                        | J19582747 + 3217197 | J13534988 + 5236102 | 2        |
| UCAC 4                       | 612-092217          | 714-051779          | 3        |
| Gaia DR3                     | 2034047656176989952 | 1560676770454324224 | 4        |
|                              | Photometry          |                     |          |
| TESS                         | 12.6396 ± 0.0074    | 12.7022 ± 0.0073    | 1        |
| V                            | 14.6 ± 0.2          | 15.355 ± 0.018      | 3        |
| B                            | ...                 | 17.035 ± 0.011      | 3        |
| BP                           | 15.510 ± 0.030      | 15.58 ± 0.02        | 4        |
| Gaia                         | 13.930 ± 0.030      | 13.66 ± 0.02        | 4        |
| RP                           | 12.710 ± 0.030      | 12.76 ± 0.02        | 4        |
| J                            | 10.997 ± 0.023      | 11.069 ± 0.02       | 2        |
| H                            | 10.425 ± 0.022      | 10.467 ± 0.02       | 2        |
| K                            | 10.175 ± 0.019      | 10.217 ± 0.02       | 2        |
| WISE 3.4 µ m                 | 9.940 ± 0.030       | 10.041 ± 0.03       | 5        |
| WISE 4.6 µ m                 | 9.786 ± 0.030       | 9.877 ± 0.030       | 5        |
| WISE 12 µ m                  | 9.459 ± 0.30        | 9.723 ± 0.038       | 5        |
| WISE 22 µ m                  | 8.650 ± 0.407       | 9.019 ± 0.378       | 5        |
| yPS                          | ...                 | 12.307 ± 0.030      | 6        |
| zPS                          | ...                 | 12.634 ± 0.030      | 6        |
|                              | Astrometry          |                     |          |
| RA (J2000)                   | 19 58 27.46         | 13 53 50.10         | 4        |
| DEC (J2000)                  | + 32 17 19.18       | + 52 36 10.24       | 4        |
| RA PM (mas / yr)             | -37.926 ± 0.015     | 108.546 ± 0.015     | 4        |
| DEC PM (mas / yr)            | -346.509 ± 0.016    | -19.281 ± 0.017     | 4        |
| Radial velocity (km / s)     | -10.20 ± 2.21       | -14.69 ± 1.64       | 4        |
| Parallax (mas)               | 31.1664 ± 0.0149    | 24.3773 ± 0.0162    | 4        |
| U galactic velocity (km / s) | 3.78 ± 0.91         | 14.43 ± 0.88        | 6        |
| V galactic velocity (km / s) | -36.4 ± 1.6         | -32.3 ± 0.48        | 6        |
| W galactic velocity (km / s) | -42.9 ± 1.2         | -6.9 ± 1.3          | 6        |

1. Stassun et al. (2018), 2. Cutri et al. (2003), 3. Zacharias et al. (2012), 4. Gaia Collaboration (2020), 5. Cutri et al. (2013). 6. This work.

astrophysical sources; and the nearby false positive probability (NFPP), which quantifies the probability that the transit signal originates from a resolved nearby star rather than the target star. The thresholds for statistically validating a planet candidate are FPP &lt; 0 01 and NFPP . &lt; 0 001. .

## 3.4. Planet search and detection limits

We ran TRICERATOPS on TESS data supplied with contrast curves from NIRC2 speckle imaging for TOI-5713 and from Palomar for TOI-6002. As our ground-based photometry confirmed the transit events were on target, removing the possibility of a nearby false positive, we set the NFPP to 0. The inferred FPP is 0.0024 ± 0.001 for TOI-5713 b and 0.00783 ± 0.00016 for TOI-6002 b. Independently, we also ran TRICERATOPS on the data obtained with MuSCAT3 (for TOI-6002 b) and MuSCAT2 (for TOI-5713 b) because they present a higher photometric precision. Also supplied with contrast curves, this led to FPP and NFPP lower than 6 × 10 -4 for both candidates. Given these results using both TESS and ground-based data independently, we consider the two candidates as statistically validated planets.

As mentioned in Sect. 2.1, TOI-5713 and TOI-6002 have both been observed in multiple sectors of TESS - respectively, sectors 16, 22, 23, 46, and sectors 14, 41, 54, 55. Since both planets were detected by the TESS pipelines with a low signal-to-noise ratio, it is possible additional planets remained unnoticed below the threshold set to trigger an alert (Multiple Event Statistic = 7.1 σ ) (Tenenbaum et al. 2012). In this context, we performed a search for additional planet candidates in the TESS data with the SHERLOCK 6 pipeline (Pozuelos et al. 2020; Demory et al. 2020; Dévora-Pajares et al. 2024) and assessed the detection limits with injection and recovery tests with the MATRIX 7 ToolKit (Pozuelos et al. 2020; Dévora-Pajares &amp; Pozuelos 2022).

The SHERLOCK pipeline is an open-source package that provides six di ff erent modules to find and validate transit signals

6 https://github.com/franpoz/SHERLOCK

7 https://github.com/PlanetHunters/tkmatrix

Table 3. Stellar parameters. Parameters in bold are the adopted stellar values in the analyses.

| Parameter          | TOI-6002                                    | TOI-5713                               | Source                                                            |
|--------------------|---------------------------------------------|----------------------------------------|-------------------------------------------------------------------|
| Sp. Type           | M3.5 ± 0.5 M4 ± 0.5                         | M3.5 ± 0.5 M4 ± 0.5                    | IRTF / SpeX a Shane / Kast a                                      |
| T e ff / K         | 3241 + 80 - 60 3477 ± 100 3314 ± 80         | 3228 ± 41 3277 ± 100 3300 ± 80         | Joint fit b Mann et al. (2015) Stefan-Boltzmann                   |
| [Fe / H]           | -0.18 ± 0.12 + 0.04 ± 0.20                  | -0.03 ± 0.13 + 0.08 ± 0.20             | IRTF / SpeX c Shane / Kast d                                      |
| M ⋆ / M ⊙          | 0.2409 + 0 . 0066 - 0 . 0065                | 0.2654 ± 0.0061                        | Joint fit b                                                       |
| R ⋆ / R ⊙          | 0.2105 + 0 . 0049 - 0 . 0048                | 0.2985 + 0 . 0074 - 0 . 0073           | Joint fit b                                                       |
| L ⋆ . 10 - 3 / L ⊙ | 5.80 + 0 . 51 - 0 . 38 6.39126 ± 0.5621     | 8.73 + 0 . 27 - 0 . 29 9.1564 ± 0.8340 | Joint fit b BC K (Mann et al. 2015)                               |
| log g ⋆ / dex      | 4.998 ± 0.025                               | 4.912 ± 0.023                          | Joint fit b                                                       |
| ρ⋆ / g cm - 3      | 21.2 + 1 . 9 - 1 . 7 23 . 4 + 1 . 4 - 3 . 2 | 14.1 + 1 . 1 - 1 . 0 11 . 4 + 7 - 6    | Joint fit b Transits fit alone                                    |
| Age / Gyr          | ≳ 4 . 5 ≲ 8 . 2 > 2 . 19 ± 0 . 10           | ≲ 4 . 5 ≲ 8 . 2 < 1 . 81 ± 0 . 10      | H α emission Kinematics Mass-spindown relation (Pass et al. 2024) |

a Classification based on spectral templates from Bochanski et al. (2007) and spectral type index relations from Reid et al. (1995), / Martín et al. (1999), Lépine et al. (2003) and Riddick et al. (2007).

b Ajoint fit of the spectral energy distribution (SED) and transits with the use of Mann et al. (2015) and Mann et al. (2019) empirical relations (see Sect. 3.1.2 and 3.6).

c Metallicity based on measurement of K -band Na i and Ca i doublets and the H2O-K2 index (Rojas-Ayala et al. 2012), and the calibrations of Mann et al. (2014).

d Metallicity based on measurement of the ζ index Lépine et al. (2007); Mann et al. (2013b) and calibrations from Mann et al. (2013b).

in photometric time-series observations: (1) acquiring automatically the data from an online database, such as MAST in the case of TESS data; (2) searching for planetary transit signals; (3) performing a vetting of the detected signals; (4) conducting a statistical validation of the vetted signals; (5) modeling them to retrieve their ephemerides; and (6) finding the upcoming transit windows observable from ground-based observatories. See Delrez et al. (2022) and Pozuelos et al. (2023) for recent applications and further details.

number of detrends to 7, set the maximum number of runs to 5, and also kept the S N threshold to 5. /

SHERLOCK obtains the PDC-SAP fluxes from the MAST archive, and applies a bi-weight filter using the W¯tan o package (Hippke et al. 2019) to detrend the light curve with various window sizes. The transit search is then performed over the nominal light curve as well as the detrended ones using the Transit Least Square (TLS) package (Hippke &amp; Heller 2019). This strategy maximizes the signal-to-noise ratio (S N) and the sig-/ nal detection e ffi ciency (SDE) of the transit signals found. In the case of TOI-6002 and TOI-5713, we downloaded the shortcadence data from TESS and applied seven di ff erent window sizes between 0.2 and 1.0 days to detrend the data. As the transit search is performed in an iterative way, if a signal is detected with a S N / ⩾ 5, it is masked, and a new search is implemented. The operation is repeated until the algorithm no longer finds signals above the threshold or reaches the number of loops (socalled 'runs') set by the user. For all our analyses, we kept the

We performed a first test to see if the first planet is recovered well by SHERLOCK for both systems with a period grid between 1 and 15 days. For TOI-6002, the signal was not recovered in the first instance. We then applied a Savitzky-Golay digital filter (Savitzky &amp; Golay 1964) and repeated the test, where we indeed recovered the planet. TOI-6002 b was recovered in the second run and in seven of the detrended light curves with a maximum S N of 12.8 and an SDE of 13.0. Similarly, TOI-5713 b was re-/ covered in the first run and in six of the detrended light curves and, at best, with an S N of 11.8 and an SDE of 23.1. In order / to assess the presence of additional candidates, we performed a search for the period grid between 0.5 and 25 days on 7 light curves detrended by a bi-weight filter with window sizes from 0.2 and 1.2 days, allowing a maximum of 5 runs. Neither of these analyses yielded positive planetary-like signals; all signals could be attributed to noise, variability, or systematics.

The lack of additional transiting planet candidates in the TESS data for both systems does not necessarily exclude the existence of extra planets. Several scenarios could explain the lack of interesting signals in our transit search with SHERLOCK : (1) the systems only host one planet; (2) if other planets exist, they are not transiting; (3) the data we examined do not have the photometric precision required to detect the additional tran-

siting planets, or they have a longer period than the period grid searched in our analyses (Wells et al. 2021). In the case of scenario (2), radial velocity measurements may be able to detect them, as discussed in Sect. 4.2. In the case of scenario (3), a ground-based high-precision photometric follow-up could possibly detect more planets, such as was the case for LP 890-9 c (Delrez et al. 2022). However, to test scenario (3) against the TESS data, we performed injection and recovery tests over the available short-cadence data for both systems.

MATRIX is an injection and recovery software specifically built to establish the retrieval thresholds of hypothetical transiting exoplanet signals in space-based data for a given star. It is built on top of the same core features as SHERLOCK so that it can reproduce the same algorithms and techniques implemented in the pipeline, and their results can be compared. That is, MATRIX sets detection limits to SHERLOCK , and by extension, to any similar pipeline using preliminary detrending and phase-folded least squares search techniques. To do so, three-dimensional grids (orbital period, planet size, orbital phase) are built to create synthetic scenarios. For each of these scenarios, MATRIX runs the pre-processing of the curve (smooth, noise-masks, detrending) and searches for the injected planet either using TLS or Box Least Squares (BLS; Kovács et al. 2002).

Fig. 6. Top panel: TOI-5713 injection and recovery test, where TOI5713 b is rendered with a red marker. Bottom panel: TOI-6002 injection and recovery test, where TOI-6002 b is rendered with a red marker.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-04a39e40-3c84-4f1e-aa64-045bfc9d14b6.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-c2d35503-747c-4c74-b5d8-c80e0fcadb40.jpg)

As both systems' planets are very similar in orbital periods and sizes, we planned the same scenarios for both of them. We wanted to explore the recoverability of planets with sizes between 0.5 and 3 R ⊕ and orbital periods from 1 up to 25 days. Therefore, the final injection and recovery test cases included a grid of 25 periods (from 1 to 25 days), 15 planet radii (from 0.5 to 3 R ⊕ ), and 3 orbital phases, summing a total of 1125 independent injected synthetic transiting planets for each target.

Our tests show that TOI-5713 b and TOI-6002 b lie at the detectability limits found for each of their systems (see Fig. 6). This is consistent with the low S N obtained for both planets, / either in SHERLOCK ( S / N TOI -5713b = 11 8, . S / N TOI -6002b = 12 8) . or SPOC ( S / N TOI -5713b = 7 7, . S / N TOI -6002b = 8 1), and the low . number of SHERLOCK detrends where they were found. It is es- sential to recognize that the S N reported by the SHERLOCK / algorithm is derived from the TLS package, which employs a straightforward estimation approach that is not directly comparable to the metrics produced by alternative pipelines, such as SPOC. Instead, SHERLOCK utilizes these values internally to evaluate the detected signals and identify the most significant ones among them.

From these results, we rule out the existence of transiting planets larger than 2 3 . R ⊕ with orbital periods up to 25 days around TOI-6002, whilst the same statement applies for planets larger than 3 R ⊕ in TOI-5713. If additional transiting exoplanets were present in the analyzed period range, they would probably be within the limits or below the super-Earth radius regime in TOI-6002. On the other side, there could also be undetected sub-Neptune transiting exoplanets in TOI-5713, especially for outermost orbits above 14 days, where planets with a radius between 2 R ⊕ to 2 75 . R ⊕ could remain undetected with detection rates between 33-66%.

## 3.5. Stellar activity and rotation period

If the host star is magnetically active and has an inhomogeneous photosphere, its rotation can pose a challenge for accurately measuring an exoplanet's mass via radial velocity (RV) observations, especially if the rotation time scale is similar to the orbital period of the planet. Therefore, we examined the e ff ect of stellar activity on RV observations.

We first visually inspected the PDC-SAP data light curves of the two targets. We observed no rotational modulations nor flares for TOI-6002, while we found clear flares in the four observed sectors of TOI-5713 in addition to hints of rotational modulation.

We then used the Systematics-Insensitive Periodogram ( TESS-SIP 8 ) (Hedges et al. 2020) to build a periodogram using the whole TESS photometry (4 sectors) for each of the two targets. This tool simultaneously creates a Lomb-Scargle periodogram and detrends TESS systematics as described in Angus et al. (2016). The SPOC pipeline is known to attenuate rotational periods, thus TESS-SIP reproduces simple aperture photometry (SAP) light curves from the TPFs and the SPOC apertures. We applied this procedure on all four sectors observed by TESS. Searching for periods between 2 and 15 days, we applied TESS-SIP on the target and on all the background pixels around the SPOC apertures. We found no significant rotational signal for TOI-6002; however, we found a significant rotational period at 9.32 ± 0.20 days for TOI-5713. Thus, we report the SIP power spectrum for only TOI-5713 in Fig. 7. We also found that this signal remains on the PDC-SAP light curves of each individual sector, which underscores the robustness of the signal. Beyond 15 days, we did not find any rotational signal.

Given these rotation-period constraints, we can estimate the ages of the stars using the mass-dependent spindown relation from Pass et al. (2024). The 'jump' times at which fully convective M dwarfs transition from rapidly to slowly rotating are 2 19 . ± 0 10 Gyr and 1 81 . . ± 0 10 Gyr for stars with the masses . of TOI-6002 and TOI-5713, respectively (see Pass et al. 2024, Equation 2). Assuming that our non-detection of rotational modulations for TOI-6002 owes to a long rotational period ( &gt; 10 d), we can constrain its age to &gt; 2 19 . ± 10 Gyr. By contrast, the relatively rapid rotation of TOI-5713 implies an age of &lt; 1 813 . ± 0 10 Gyr. These estimates are consistent with the age estimates . we derive from their H α emission (for TOI-5713) or lack of it

8 https://github.com/christinahedges/TESS-SIP

(for TOI-6002) as well as the likely thin-disk membership of both stars.

Fig. 7. TESS-SIP power spectrum of TOI-5713 based on the four TESS sectors. We calculate the periodogram for both the corrected light curve (top panel) and the background (BKG) pixels (bottom panel). The star's periodogram shows some power excess compared to the BKG.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-ef613ca7-0039-4e30-8816-d4558ba19352.jpg)

## 3.6. Global modeling

For each planet, we used the latest version of EXOFASTv2 (Eastman et al. 2019) software package to globally analyze the TESS and ground-based data described in Sect. 2, jointly with the star's SED fit and the empirical relations Mann et al. (2015) and Mann et al. (2019) to characterize the host stars (see description in Sect. 3.1.2).

In the priors file, we set starting values of M ∗ , R ∗ and Tef f . We set lower and upper limits priors on the orbital period ( P ± 0 1 . P d) and the transit epoch ( Tc ± 0 33 . P ) from the values determined by SPOC and reported in ExoFOP. For the LD parameters, we disabled the priors based on the Claret tables (Claret 2017; Claret &amp; Bloemen 2011) by setting the NOCLARET flag given the low mass of our stars where these tables are not reliable (see Eastman et al. 2019, and references therin). Instead, we fitted for quadratic LD parameters assuming Gaussian priors of variance 0.1, computed with the PyLDTk code (Parviainen 2015) for each passband (see Table 4). Regarding the TESS data, we used the detrended PDC-SAP light curves (see Sect. 2.1) which are corrected for the contribution of neighboring stars and subsequently corrected for sky background bias for TOI-6002 in sector 14 (see Sect. 2.1), but we still fit for a possible dilution by applying a Gaussian prior of 0 ± 10% of the contamination ratio reported in ExoFOP. This is to propagate any uncertainty in the correction. Also, Gaussian priors were applied on the metallicities from IRTF SpeX (see Table 3). For each planet, we fitted / for circular and eccentric orbits to check for orbital eccentricity evidences using only photometric data. The two fits were run until convergence. We found the di ff erence in the loglikelihood is ∆ ln Z = ln Z eccentric -ln Z circular ⩾ 2 for both TOI-6002 b and TOI-5713 b, which suggests that eccentric orbits are marginally

Table 4. Priors assumed for the quadratic limb-darkening coe ffi cients for each passband for both stars TOI-6002 and TOI-5713.

| Filter   | u 1          | u 2          |
|----------|--------------|--------------|
| TOI-6002 |              |              |
| TESS     | 0.2534 ± 0.1 | 0.3798 ± 0.1 |
| g ′      | 0.4895 ± 0.1 | 0.4039 ± 0.1 |
| r ′      | 0.4992 ± 0.1 | 0.3384 ± 0.1 |
| i ′      | 0.2856 ± 0.1 | 0.4168 ± 0.1 |
| z s      | 0.1907 ± 0.1 | 0.4718 ± 0.1 |
| TOI-5713 |              |              |
| TESS     | 0.2806 ± 0.1 | 0.3069 ± 0.1 |
| r ′      | 0.5388 ± 0.1 | 0.2842 ± 0.1 |
| i ′      | 0.3186 ± 0.1 | 0.3422 ± 0.1 |
| z s      | 0.2383 ± 0.1 | 0.3334 ± 0.1 |
| I + z    | 0.2222 ± 0.1 | 0.2880 ± 0.1 |

favored for both planets. Table 5 presents the fits' results for both the circular and eccentric cases. Fig. 8 shows the detrended and modeled gound-based light curves for both TOI-6002 b and TOI5713 b, and Fig. 9 shows the TESS phase-folded light curves.

We also re-fitted the data for a circular orbit and with the same priors as explained above but without any Gaussian prior on the stellar mass. This is to check for consistency between the stellar densities determined from R ∗ and M ∗ (from SED + empirical relations) and those determined from the transits via the a R / ∗ (the semi-major axis over the stellar radius) ratio. We found ρ ∗ = 23 4 . + 1 4 . -3 2 . g cm -3 for TOI-6002 b and 11 4 . + 7 -6 g cm -3 for TOI-5713 b, which are consistent with the previously determined stellar densities (see Table 3). This provides additional assurance that the transit signals for both candidates originate from transiting planets and not astrophysical false positives.

## 4. Results and discussion

## 4.1. Radius valley

The radius valley, discovered by Fulton et al. (2017) from 1.5 to 2 R ⊕ for solar-type stars and by Cloutier &amp; Menou (2020) from 1.4 to 1.7 R ⊕ for low-mass stars (i.e., M ≲ 0 65 . M ), has ⊙ attracted attention because it is important to understand planet formation and evolution. There are several theories to explain the origin of the radius valley: (1) thermally-driven mass loss, where planets lose their atmospheric envelopes because of stellar high energetic X-ray and ultra-violet (UV) radiation (Lopez &amp; Fortney 2013; Owen &amp; Wu 2013; Jin et al. 2014; Chen &amp; Rogers 2016) and of the energy emerging from the cooling planetary cores (core-powered mass loss, Owen &amp; Wu 2013; Ginzburg et al. 2016, 2018); (2) gas-poor formation, where some low-mass planetary cores accrete a thin gaseous envelope (see e.g. Lee &amp; Connors 2021; Cherubim et al. 2023). These planets can still experience thermally-driven mass loss; and (3) gas-depleted formation, where the population of super-Earths is formed late in a gas-depleted environment (see e.g. Lopez &amp; Rice 2018; Cloutier &amp; Menou 2020; Cherubim et al. 2023); (4) Other studies based on coupled formation and evolution models predict that the subNeptunes are rich in water formed outside the ice-line and migrate toward the host stars (Venturini et al. 2020; Burn et al. 2024; Venturini et al. 2024).

Fig. 8. Phase-folded detrended ground-based transit light curves with time in hours from mid-transit of TOI-6002 b (left) and TOI-5713 b (right). The unbinned light curves are shown with gray points. The binned points (bin size = 6 min) are shown in black with the corresponding error bars. The solid-colored lines corresponds to the best-fit transit model from the final joint fit.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-55eb03ff-cfb4-4fc1-99f9-7d71f8c9ae91.jpg)

Fig. 9. TESS phase-folded detrended light curve with the time in hours from mid-transit of TOI-6002 b (left) and TOI-5713 b (right). The red solid lines represent the best-fit transit model from the global fit. The gray dots are the 2-min TESS data and the black dots with error bars show the data binned every 5 minutes. The residuals are plotted in the bottom panel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-2afbcfc4-cb5e-4394-9777-a1275a6ab233.jpg)

Table 5. Median values and 68% confidence intervals for the parameters of TOI-6002 b and TOI-5713 b obtained using EXOFASTv2 for circular and eccentric orbits.

| Parameter             | Description                                                  | TOI-6002 b                            | TOI-6002 b                            | TOI-5713 b                            | TOI-5713 b                            |
|-----------------------|--------------------------------------------------------------|---------------------------------------|---------------------------------------|---------------------------------------|---------------------------------------|
|                       |                                                              | Circular                              | Eccentric                             | Circular                              | Eccentric                             |
| Stellar Parameters:   | Stellar Parameters:                                          | Stellar Parameters:                   | Stellar Parameters:                   | Stellar Parameters:                   | Stellar Parameters:                   |
| M ∗ . . . .           | Mass ( M ⊙ ) . . . . . . . . . . . . . . . . . . . . . . .   | 0 . 2105 + 0 . 0049 - 0 0048          | 0 . 2101 + 0 . 0049 - 0 0048          | 0 . 2653 + 0 . 0062 - 0 0061          | 0 . 2653 ± 0 . 0061                   |
| R ∗ . . . .           | Radius ( R ⊙ ) . . . . . . . . . . . . . . . . . . . . . .   | 0 . 2409 + 0 . 0066 - 0 0065          | 0 . 2432 + 0 . 0063 - 0 0065          | 0 . 2984 ± 0 . 0073                   | 0 . 2992 + 0 . 0073 - 0 0072          |
| L ∗ . . . .           | Luminosity (10 - 3 L ⊙ ) . . . . . . . . . . . . .           | . 5 . 80 + 0 . 51 - 0 38              | . 5 . 82 + 0 . 48 - 0 36              | 8 . 74 + 0 . 27 - 0 28                | . 8 . 74 + 0 . 28 - 0 29              |
| ρ ∗ . . . .           | Density (cgs) . . . . . . . . . . . . . . . . . . . . .      | 21 . 2 + 1 . 9 - 1 7                  | 20 . 6 + 1 . 8 - 1 5                  | 14 . 1 + 1 . 1 - 1 0                  | 13 . 96 + 1 . 1 - 0 99                |
| log g . .             | Surface gravity (cgs) . . . . . . . . . . . . . .            | 4 . 998 ± 0 . 025                     | 4 . 989 + 0 . 025 - 0 024             | 4 . 912 ± 0 . 023                     | 4 . 910 + 0 . 023 - 0 022             |
| T e ff . . .          | E ff ective Temperature (K) . . . . . . . . .                | 3241 + 82 - 60                        | 3229 + 77 - 57                        | 3229 + 41 - 40                        | 3225 + 41 - 40                        |
| Planetary Parameters: | Planetary Parameters:                                        | Planetary Parameters:                 | Planetary Parameters:                 | Planetary Parameters:                 | Planetary Parameters:                 |
| P . . . . .           | Orbital period (days) . . . . . . . . . . . . . .            | 10 . 904821 + 0 . 000021 - 0 . 000018 | 10 . 904821 + 0 . 000021 - 0 . 000018 | 10 . 441988 + 0 . 000014 - 0 . 000013 | 10 . 441989 + 0 . 000015 - 0 . 000014 |
| R P . . . .           | Planet radius ( R ⊕ ) . . . . . . . . . . . . . . . .        | 1 . 63 + 0 . 21 - 0 19                | 1 . 65 + 0 . 22 - 0 19                | 1 . 75 + 0 . 16 - 0 14                | 1 . 77 + 0 . 17 - 0 14                |
| T C . . . .           | Transit time (BJD TDB ) . . . . . . . . . . . . .            | 2458692 . 7636 + 0 . 0024 - 0 0028    | 2458692 . 7639 + 0 . 0024 - 0 0028    | 2458745 . 6773 + 0 . 0014 - 0 0016    | 2458745 . 6776 + 0 . 0015 - 0 0017    |
| a . . . . .           | Semi-major axis (AU). . . . . . . . . . . . .                | 0 . 05725 ± 0 . 00044                 | 0 . 05722 ± 0 . 00044                 | 0 . 06008 ± 0 . 00046                 | 0 . 06008 ± 0 . 00046                 |
| i . . . . . .         | Inclination (Degrees). . . . . . . . . . . . . .             | 89 . 657 + 0 . 14 - 0 099             | 89 . 68 ± 0 . 21                      | 89 . 385 + 0 . 10 - 0 088             | 89 . 37 + 0 . 36 - 0 21               |
| e . . . . .           | Eccentricity . . . . . . . . . . . . . . . . . . . . .       | 0 (Fixed)                             | 0 . 15 + 0 . 42 - 0 11                | 0 (Fixed)                             | 0 . 24 + 0 . 40 - 0 17                |
| τ circ                | Tidal circularization timescale (Gyr)                        | 6200 + 4200 - 2400                    | 3500 + 3900 - 3300                    | 4600 + 2300 - 1500                    | 2000 + 2600 - 2000                    |
| ω ∗                   | Argument of Periastron (Degrees)                             | 90 (fixed)                            | - 160 + 140 - 110                     | 90 (fixed)                            | - 139 + 120 - 99                      |
| T eq . . .            | Equilibrium temperature a (K) . . . . . .                    | 321 . 1 + 7 . 0 - 5 6                 | 321 . 4 + 6 . 6 - 5 2                 | 347 . 2 + 3 . 0 - 3 1                 | 347 . 2 + 3 . 0 - 3 2                 |
| R P / R ⋆             | Planet-to-star radius ratio . . . . . . . . . .              | . 0 . 0622 + 0 . 0080 - 0 0070        | . 0 . 0622 + 0 . 0080 - 0 0070        | . 0 . 0539 + 0 . 0048 - 0 0039        | . 0 . 0541 + 0 . 0050 - 0 0041        |
| a / R ⋆ . .           | Semi-major axis in stellar radii . . . . .                   | 51 . 1 + 1 . 5 - 1 4                  | 50 . 6 + 1 . 4 - 1 3                  | 43 . 3 ± 1 . 1                        | 43 . 2 + 1 . 1 - 1 0                  |
| δ . . . . .           | ( R P / R ⋆ ) 2 (ppt). . . . . . . . . . . . . . . . . . . . | . 3 . 87 + 1 . 1 - 0 82               | . 3 . 87 + 1 . 1 - 0 83               | 2 . 90 + 0 . 54 - 0 40                | . 2 . 92 + 0 . 56 - 0 42              |
| δ g ′ . . . .         | Transit depth in g ′ (ppt) . . . . . . . . . . .             | 4 . 7 + 1 . 3 - 1 0                   | 4 . 7 + 1 . 3 - 1 0                   | ...                                   | ...                                   |
| δ i ′ . . . .         | Transit depth in i ′ (ppt) . . . . . . . . . . . .           | 4 . 48 + 1 . 2 - 0 95                 | 4 . 47 + 1 . 2 - 0 97                 | 3 . 32 + 0 . 63 - 0 48                | 3 . 30 + 0 . 63 - 0 47                |
| δ r ′ . . . .         | Transit depth in r ′ (ppt). . . . . . . . . . . .            | 5 . 0 + 1 . 4 - 1 1                   | 5 . 0 + 1 . 4 - 1 1                   | 3 . 76 + 0 . 72 - 0 54                | 3 . 70 + 0 . 74 - 0 56                |
| δ z ′ . . . .         | Transit depth in z ′ (ppt). . . . . . . . . . . .            | 4 . 16 + 1 . 1 - 0 89                 | 4 . 16 + 1 . 1 - 0 90                 | 3 . 16 + 0 . 60 - 0 45                | 3 . 16 + 0 . 60 - 0 44                |
| δ TESS . .            | Transit depth in TESS (ppt) . . . . . . . .                  | 4 . 22 + 1 . 5 - 0 63                 | 4 . 39 + 1 . 2 - 0 95                 | 3 . 28 + 0 . 61 - 0 47                | 3 . 27 + 0 . 61 - 0 46                |
| τ . . . . .           | Ingress / egress transit duration (days)                     | 0 . 00444 + 0 . 00062 - 0 00054       | 0 . 00442 + 0 . 00077 - 0 00058       | 0 . 00469 + 0 . 00052 - 0 00045       | 0 . 00457 + 0 . 0021 - 0 00087        |
| T 14 . . .            | Total transit duration (days) . . . . . . . .                | 0 . 06913 + 0 . 00081 - 0 00080       | 0 . 06904 + 0 . 00095 - 0 00088       | 0 . 0728 + 0 . 0016 - 0 0018          | 0 . 0728 ± 0 . 0023                   |
| b . . . . .           | Transit impact parameter . . . . . . . . . .                 | 0 . 306 + 0 . 078 - 0 12              | 0 . 26 + 0 . 18 - 0 17                | 0 . 464 + 0 . 057 - 0 071             | 0 . 44 + 0 . 21 - 0 27                |
| ⟨ F ⟩ . . .           | Incident Flux ( S ⊕ ) . . . . . . . . . . . . . . . .        | 1 . 7705 + 0 . 1616 - 0 1175          | 1 . 6823 + 0 . 1689 - 0 3673          | 2 . 4170 ± 0 . 0881                   | 2 . 2554 + 0 . 1689 - 0 6024          |
| d / R ∗ . .           | Separation at mid-transit . . . . . . . . . .                | 51 . 1 + 1 . 5 - 1 4                  | 48 . 2 + 4 . 0 - 7 7                  | 43 . 3 ± 1 . 1                        | 40 . 5 + 6 . 0 - 7 8                  |
| Predicted Parameters: | Predicted Parameters:                                        | Predicted Parameters:                 | Predicted Parameters:                 | Predicted Parameters:                 | Predicted Parameters:                 |
| M P . . .             | Planet mass b ( M ⊕ ) . . . . . . . . . . . . . . . .        | 3 . 8 + 1 . 6 - 1 0                   | 3 . 8 + 1 . 6 - 1 0                   | 4 . 3 + 1 . 7 - 1 1                   | 4 . 3 + 1 . 7 - 1 1                   |
| K . . . . .           | RV semi-amplitude (m / s) . . . . . . . . . .                | 3 . 10 + 1 . 3 - 0 85                 | 3 . 4 + 1 . 8 - 1 0                   | 3 . 02 + 1 . 2 - 0 76                 | 3 . 40 + 1 . 8 - 0 96                 |
| TSM..                 | Transmission spectroscopy metric . .                         | 60.0                                  | 50.0                                  | 38.0                                  | 40.0                                  |
| M P / M ∗             | Planet-to-star mass ratio ( × 10 - 5 ) . . .                 | 5 . 4 + 2 . 3 - 1 5                   | 5 . 5 + 2 . 4 - 1 5                   | 4 . 8 + 1 . 9 - 1 2                   | 4 . 9 + 1 . 9 - 1 2                   |
| ρ P . . . .           | Planet mean density (cgs). . . . . . . . . .                 | 4 . 8 + 1 . 9 - 1 . 3                 | 4 . 7 + 1 . 9 - 1 . 2                 | 4 . 3 + 1 . 6 - 1 . 0                 | 4 . 3 + 1 . 6 - 1 . 0                 |

a The equilibrium temperature is estimated supposing a null-albedo and a 100% e ffi cient heat recirculation. b The planetary mass is estimated using the Chen &amp; Kipping (2017) mass-radius relation. The radial velocity semi-amplitude ( K ) is predicted using the estimated mass.

Observational and theoretical studies showed that the position of the radius valley is period-dependent (see Fig. 10). Its location has been found trending to smaller planet radii for long periods in both thermally-driven mass loss and gas-poor formation models for solar FGK stars (see e.g. Fulton et al. 2017; Van Eylen et al. 2018; Lopez &amp; Rice 2018; Lee &amp; Connors 2021). For Mdwarfs, its location has been found trending to larger radii for long periods in the gas-depleted theory (see e.g Lopez &amp; Rice 2018; Cloutier &amp; Menou 2020). Conversely, other studies on small planets around M dwarfs found a negative slope (see e.g. Van Eylen et al. 2021; Bonfanti et al. 2024), something similar to that of FGK stars.

As shown in Fig. 10, a wedge of opposing predictions emerged because of the di ff erence in the slopes of the radius valley. Planets in this region are known as "keystone planets." They are predicted to be purely rocky by thermally-driven mass loss and gas-poor formation, and potentially rocky with small gaseous envelopes by the gas-depleted formation model (see e.g. Cherubim et al. 2023). Only a few planets with measured densities are in this wedge. Having a large sample of keystone planets is important to understand the models dominating the formation and evolution of small planets. TOI-5713 b ( Rp = 1 75 . + 0 14 . -0 16 . R ⊕ and P = 10 441989 . + 0 000015 . -0 000014 . d) is located within the wedge of keystone planets with the possibility of being at its upper edge, given the uncertainty on the planetary radius. TOI-6002 b ( Rp = 1 63 . + 0 21 . -0 19 . R ⊕ and P = 10 904821 . + 0 000021 . -0 000018 . d) is located at the lower edge of this wedge with the possibility of being inside it, given the uncertainty on the planetary radius. When their masses are determined with RV observations (see Sect.4.2), these planets will join the small sample of keystone planets with density measurements. This sample will, in the future, help disentangle which model is dominating the formation of small (i.e., Rp ≲ 4 R ⊕ ) planets.

In addition, for M dwarfs, a later study of Luque &amp; Pallé (2022) on small planets (with mass and radius precision better than 25% and 8%, respectively) identified instead a density valley separating rocky planets, with density peaked at ρ = 0 94 . ± 0 13 . ρ ⊕ , and water-rich worlds, with density peaked at ρ = 0 47 . ± 0 05 . ρ ⊕ . A third peak was identified at ρ = 0 24 . ± 0 04 . ρ ⊕ , attributed to gas-rich planets. This study favored the pebble accretion model as the main mechanism that forms these populations of planets. In this model, water-rich planets are formed outside the ice-line and migrate toward the host stars. This finding was theoretically predicted by Venturini et al. (2020), and newly supported by Burn et al. (2024) and Venturini et al. (2024) using coupled formation and evolution models. Conversely, Parviainen et al. (2024) stand against the existence of distinct population of water-rich planets around M dwarf stars. However, in all previous studies on radius and density valleys, strong conclusions require larger sample of precisely characterized exoplanets.

Given their predicted masses estimated from their radii using the relations of Chen &amp; Kipping (2017), these planets are predicted to have densities of 0 855 . + 0 345 . -0 218 . ρ ⊕ for TOI-6002 b, and 0 782 . + 0 2911 . -0 182 . ρ ⊕ for TOI-5713 b. We also used the recent SPRIGHT 9 python package of Parviainen et al. (2024), where we found median values of the masses (and then densities) consistent with those obtained using Chen &amp; Kipping (2017) relations. To further guess their compositions, we referred to the radius space of the study of Luque &amp; Pallé (2022) and the corresponding density regimes in the density space. Given their radii, TOI5713 b falls slightly towards the range of water-rich planets, and TOI-6002 b falls in the region where we find both rocky and

9 https://github.com/hpparvi/spright

water-rich planets. Therefore, it is uncertain whether the two planets could be rocky or water-rich worlds. However, even if they are not of intermediate composition, they might then be either some of the largest rocky planets or some of the smallest water-rich worlds known today. As for the radius valley, both planets straddle near the opposing side of the region of "keystone planets," where it is di ffi cult to guess whether they lost their gaseous envelopes or not. When their masses are determined, these two planets can be added to the sample of "keystone planets," which will be used to conduct in-depth statistical inferences on the radius and density valleys. Such studies will help to refine the relative dominance of the various mechanisms proposed to be responsible for the formation and evolution of super-Earths and sub-Neptunes around M-dwarf stars.

Fig. 10. Distribution of planet radii and orbital periods for all confirmed small planets hosted by low-mass stars (M ∗ ≲ 0 65 M . ⊙ ). The solid line represents the predicted location of the radius valley based on the gasdepleted formation model (Cloutier &amp; Menou 2020). The dashed line shows the predicted location of the valley for the thermally-driven photoevaporation and gas-poor formation models (Van Eylen et al. 2018). The red and blue dots represents the planets TOI-6002b and TOI-5713b, respectively. The 1D radius distribution, with the location of the two planets, is shown on the right panel.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-c0e65082-ec02-4c3e-94bb-a9c225963124.jpg)

## 4.2. Prospects for a radial velocity follow-up

In Table 5, we show the predicted semi-amplitudes for TOI6002 b and TOI-5713 b for both circular and eccentric orbits. With both planets inducing reflex motions to their host-stars with semi-amplitudes on the order of 3 ms -1 , high-resolution, wellstabilized spectrographs are needed for their measurement. Due to the relatively faint apparent magnitudes of both stars, reaching such precision will be challenging for instruments on 4-m class telescopes, such as HARPS-N (Cosentino et al. 2012) in the optical. This is also true for NIR instruments such as CARMENES (Quirrenbach et al. 2020) or SPIRou (Donati et al. 2018) which would allow to reach a RV precision in the order of 10 ms -1 or less. The MAROON-X instrument (Seifahrt et al. 2020) at the 8.1-m Gemini North telescope o ff ers the necessary stability and throughput to successfully measure the reflex motions of both stars. Assuming exposure times of 900 s, the instrument's red arm will achieve a maximum signal-to-noise-ratio larger than 95 and 65 for TOI-6002 and TOI-5713, respectively.

Following Engle &amp; Guinan (2023), we can expect negligible line broadening and a low activity level for TOI-6002 due to the

age of this star. Assuming an additional noise term of 1 ms -1 accounting for stellar activity, we expect that the reflex motion can be detected at a 5 σ level with less than 11h telescope time using MAROON-X. This would not only allow us to measure the mass of TOI-6002 b, but also to derive the eccentricity of the planet's orbit.

In Sect. 3.5, we determine the stellar rotation period of TOI5713 from photometry. This translates to a maximum rotational line broadening ( v sin i ∗ ) on the order of 1 6 km s . -1 , which is below the resolution limit for most high-resolution spectrographs, and thus does not a ff ect the radial-velocity precision. Nevertheless, we observe moderate activity and a photometric rotational variability close to the orbital period of TOI-5713 b, making activity mitigation challenging. Assuming an optimistic noise term of 3 ms -1 accounting for this stellar activity after successful mitigation, a 5 σ detection of the radial velocity signal due to TOI5713 b would need at least 16 hr telescope time with MAROONX.

## 4.3. Potential for atmospheric characterization

TOI-6002b and TOI-5713b, with radii of 1.63 R ⊕ and 1.75 R ⊕ and incident stellar fluxes of 1.77 and 2.42 times that of Earth respectively, could be classified as 'super-Venus' planets if they have rocky compositions. Venus, in comparison, receives about 1.9 times the flux of Earth and is 95% the size of Earth. According to Kane et al. (2014), Venus analogs are defined as predominantly rocky planets within the "Venus-Zone (VZ)," receiving insolation fluxes between approximately 0.95 and 25 times that of Earth. This classification suggests that both TOI-6002b and TOI-5713b fall within this category. With incident flux slightly lower (TOI-6002b) and slightly higher (TOI-5713b) than modern Venus, these two exoplanets are intriguing candidates to expand our understanding of the evolution of rocky planets in the state from a temperate Earth-like to a hot Venus-like planet. The incident flux levels indicate that both planets might experience or might have experienced strong runaway greenhouse e ff ects, depending on factors like the initial water inventory, planetary evolution and atmospheric escape rates (Kaltenegger et al. 2023a). Combined with the estimated ages of the stars, this could indicate atmospheric conditions similar to modern Venus. However, studying planets with incident irradiation higher than modern Earth to modern Venus levels and beyond, that may be experiencing strong greenhouse e ff ects or may have evolved into a post-runaway greenhouse state is crucial, as it can provide insights into the factors that cause some rocky planets to resemble Venus-analogs while others may maintain Earth-analog characteristics (Ehrenreich et al. 2012; Chouqar et al. 2023; Jordan et al. 2021; Kaltenegger et al. 2023a).

To investigate two cases of potential atmospheric characteristics of these planets, we constructed models simulating a) cloudless and 10 mbar cloudy, 92-bar Venus-like, and b) 1-bar Earth-like atmospheres. We link the one-dimensional VULCAN chemical kinetic code (Tsai et al. 2017) to the one-dimensional radiative transfer model PETITRADTRANS to predict the transmission spectra for these clear atmospheres (Mollière, P. et al. 2019). The radiative-convective temperature-pressure profiles were computed with the HELIOS code (Malik et al. 2018). The species considered for our atmospheric models include N2, O2, O3 H2O, CO2, CO, HCN, CH4, NO, NH3, and C2H2 for the Earth-like atmosphere (rich in N2). Additionally, Rayleigh scattering for N2 and O2, as well as collision-induced absorption for the N2-N2 and O2-O2 pairs, are considered. For the Venus-like atmosphere (rich in CO2), we include H2S, SO2,

OCS, H2O, CO2, CO, CH4, HCN, and C2H2 as well as CO2CO2 collision-induced absorption, as opacity sources in the atmosphere. Rayleigh scattering is included for CO2.

Ehrenreich et al. (2012) showed that during Venus's 2012 transit, its transmission spectra indicated that the cloud and haze layers obstruct probing below an altitude of 80 km. Pidhorodetska et al. (2021) modeled the transmission spectra of the L 9859 planets, a benchmark system for studying the VZ planets (Ostberg et al. 2023). They showed that if the L 98-59 planets possess a clear Venus-like atmosphere, NIRSpec could detect CO2 within 26 transits for each planet. However, the presence of H2SO4 clouds would significantly suppress CO2 absorption. A new study by Kaltenegger et al. (2023a) shows that there is a wide range of possible atmospheres for a rocky planet experiencing a runaway greenhouse e ff ect, which changes the atmosphere and observable features significantly. Model transmission spectra of the planet LP 890-9c (Delrez et al. 2022), another benchmark system for studying the VZ planets (Kaltenegger et al. 2023a) show that depending on the state of the planet, James Webb Space Telescope (JWST) observations could infer evidence of H2O with 3 transit (at 3 σ confidence) for a full runaway greenhouse scenario for LP 890-9c, and CO2-dominated atmospheres resembling Venus without high-altitude terminator clouds with 8 transits (Gomez Barrientos et al. 2023). A 3D Venus model for the same planet with higher, thick cloud coverage increases the amount of transits needed to be able to explore the atmospheric composition. Note that all predictions could be complicated by the impact of clouds and (or) unocculted starspots, especially for cool host stars.

These results show that observing these exoplanets is of critical importance to understanding their evolutionary stage and the e ff ect of increasing solar irradiation on rocky worlds. Additionally, planetary systems like the TRAPPIST-1 system could also allow insight into the e ff ect of the increase of irradiation on the environment on a planet (e.g. Payne &amp; Kaltenegger (2024)). However, even for a specific orbital distance and stellar irradiation, the observable spectral features in transit for an exoplanet depend on their specific atmospheric evolution and cloud coverage (e.g., Kaltenegger (2017); Kaltenegger et al. (2020)), requiring a wider range of observations to explore the characteristic changes of rocky worlds with stellar irradiation. TOI-6002b and TOI-5713b provide two important datapoints in that exploration.

For TOI-6002b and TOI-5713b, the synthetic transmission spectra generated from our models (Earth-like surface temperatures are 524 K and 580 K, respectively; Venus-like surface temperatures are 1053 K and 1110 K, respectively) show key di ff erences between Earth-like and Venus-like atmospheres (see Fig. 11).

Models of TOI-6002b and TOI-5713b for clear hot Earthlike atmospheric compositions are dominated by water vapor. Water vapor is a strong absorber across the infrared for Earthlike atmospheric compositions, with amplitudes of about 300 ppm in the mid-infrared (see Fig. 11, top panel). Note that Fig. 11 (top panel) also shows a notable CO2 feature at 4-5 µ with a depth of 400 ppm. However, we would assume the atmospheric CO2 concentration to decrease during a runaway greenhouse stage, reducing that feature strength. Because we are not modeling the possible evolution of the planet's atmosphere, we have not changed the CO2 concentration here, but note that the features could be much smaller than shown.

Models of TOI-6002b and TOI-5713b using clear Venus-like atmospheric composition are primarily characterized by CO2 opacity, with a depth of about 250 ppm in the mid-infrared, and a prominent SO2 feature at 4 µ and between 7 and 10 µ , with

depths above ≈ 100 ppm. With a cloud deck at 10 mbar, the CO2 feature at approximately 4.2 µ m remains discernible, exhibiting a depth greater than 200 ppm.

These models show that TOI-6002b and TOI-5713b represent promising targets for atmospheric characterization with JWST, as they may provide valuable insights into the diverse evolutionary pathways and climate states of hot rocky or waterrich exoplanets. Mass measurements could help provide better constraints on the viability of atmospheric studies and would be possible with the MAROON-X instrument at the 8.1-m Gemini North telescope.

## 5. Conclusion

In this paper, we reported the TESS discovery and initial characterization of TOI-6002 b and TOI-5713 b, two super-Earth planets orbiting the mid-M-dwarf stars TOI-6002 (V = 14.6) and TOI-5713 (V = 15.355). We used the 2-min-cadence TESS observations from four sectors for each target, ground-based photometry, high-angular resolution imaging, and spectroscopic observations to validate the planetary nature of the detected transit events. We jointly analyzed the transit light curves observed by TESS and ground-based facilities to characterize the planets.

TOI-6002 b has a radius of 1 65 . + 0 22 . -0 19 . R , ⊕ an orbital period of 10 904821 . + 0 000021 . -0 000018 . days, and receives 1 77 . + 0 16 . -0 11 . S ⊕ . TOI5713 b has a radius of 1 77 . + 0 13 . -0 11 . R , ⊕ an orbital period of 10 441989 . + 0 000015 . -0 000014 . days, and receives 2 42 . ± 0 11 . S ⊕ . Both planets are located just outside but near the inner edge of the habitable zone around their host stars, making them interesting targets for future studies to explore the evolution of exoplanets from hot, potentially habitable to Venus-like worlds. These two planets provide compelling additions to the sample of planets critical for understanding the formation and evolution of small worlds on the basis of the radius and density valleys.

Acknowledgements. This research was carried out in part at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D0004). This article is based on observations made with the MuSCAT2 instrument, developed by ABC, at Telescopio Carlos Sánchez operated on the island of Tenerife by the IAC in the Spanish Observatorio del Teide. TRAPPIST-North is a project funded by the University of Liege (Belgium), in collaboration with Cadi Ayyad University of Marrakesh (Morocco). We acknowledge the use of public TESS data from pipelines at the TESS Science O ffi ce and at the TESS Science Processing Operations Center. This work makes use of observations from the LCOGT network. Part of the LCOGT telescope time was granted by NOIRLab through the Mid-Scale Innovations Program (MSIP). MSIP is funded by NSF. This research has made use of the Exoplanet Follow-up Observation Program (ExoFOP; DOI: 10.26134 ExoFOP5) website, which is operated by the California Institute of / Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. We acknowledge the use of public TESS data from pipelines at the TESS Science O ffi ce and at the TESS Science Processing Operations Center. This paper includes data collected by the TESS mission that are publicly available from the Mikulski Archive for Space Telescopes (MAST). Funding for the TESS mission is provided by NASA's Science Mission Directorate. KAC acknowledges support from the TESS mission via subaward s3449 from MIT. Resources supporting this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center for the production of the SPOC data products. This paper is based on observations made with the MuSCAT3instrument, developed by the Astrobiology Center and under financial supports by JSPS KAKENHI (JP18H05439) and JST PRESTO (JPMJPR1775), at Faulkes Telescope North on Maui, HI, operated by the Las Cumbres Observatory. This work was partially supported by a grant from the Erasmus + International Credit Mobility program (M. Ghachoui). This publication benefits from the support of the French Community of Belgium in the context of the FRIA Doctoral Grant awarded to MT. MG is F.R.S-FNRS Research Director. The postdoctoral fellowship of KB is funded by F.R.S.-FNRS grant T.0109.20 and by the Francqui Foundation. Author F.J.P acknowledges financial support from the Severo Ochoa grant CEX2021-001131-S funded by MCIN AEI 10.13039 501100011033 and / / /

Ministerio de Ciencia e Innovación through the project PID2022-137241NBC43. H.P. acknowledges support by the Spanish Ministry of Science and Innovation with the Ramon y Cajal fellowship number RYC2021-031798-I. This work is partly supported by JSPS KAKENHI Grant Number JP18H05439 and JST CREST Grant Number JPMJCR1761. E.E-B. acknowledges financial support from the European Union and the State Agency of Investigation of the Spanish Ministry of Science and Innovation (MICINN) under the grant PRE2020-093107 of the Pre-Doc Program for the Training of Doctors (FPI-SO) through FSE funds We acknowledge financial support from the Agencia Estatal de Investigación of the Ministerio de Ciencia e Innovación MCIN AEI 10.13039 501100011033 and / / / the ERDF 'A way of making Europe' through project PID2021-125627OB-C32, and from the Centre of Excellence 'Severo Ochoa' award to the Instituto de Astrofisica de Canarias. This material is based upon work supported by the National Aeronautics and Space Administration under Agreement No. 80NSSC21K0593 for the program 'Alien Earths'. The results reported herein benefited from collaborations and (or) information exchange within NASA's Nexus for Exoplanet System Science (NExSS) research coordination network sponsored by NASA's Science Mission Directorate. This research was supported by Wallonia-Brussels International.

## References

- Aller, A., Lillo-Box, J., Jones, D., Miranda, L. F., &amp; Barceló Forteza, S. 2020, A&amp;A, 635, A128 Angus, R., Foreman-Mackey, D., &amp; Johnson, J. A. 2016, apj, 818, 109 Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167 Astropy Collaboration, Price-Whelan, A. M., Sip˝cz, B. M., et al. 2018, AJ, 156, o 123 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&amp;A, 558, A33 Barkaoui, K., Burdanov, A., Hellier, C., et al. 2019, aj, 157, 43 Bochanski, J. J., West, A. A., Hawley, S. L., &amp; Covey, K. R. 2007, AJ, 133, 531 Bonfanti, A., Brady, M., Wilson, T. G., et al. 2024, A&amp;A, 682, A66 Brown, T. M., Baliber, N., Bianco, F. B., et al. 2013, Publications of the Astronomical Society of the Pacific, 125, 1031 Burgasser, A. J. &amp; Splat Development Team. 2017, in Astronomical Society of India Conference Series, Vol. 14, Astronomical Society of India Conference Series, 7-12 Burn, R., Mordasini, C., Mishra, L., et al. 2024, Nature Astronomy, 8, 463 Chambers, K. C., Magnier, E. A., Metcalfe, N., et al. 2016, arXiv e-prints, arXiv:1612.05560 Chen, H. &amp; Rogers, L. A. 2016, The Astrophysical Journal, 831, 180 Chen, J. &amp; Kipping, D. 2017, apj, 834, 17 Cherubim, C., Cloutier, R., Charbonneau, D., et al. 2023, AJ, 165, 167 Chouqar, J., Lustig-Yaeger, J., Benkhaldoun, Z., et al. 2023, Monthly Notices of the Royal Astronomical Society, 522, 648 Ciardi, D. R., Beichman, C. A., Horch, E. P., &amp; Howell, S. B. 2015, ApJ, 805, 16 Claret, A. 2017, A&amp;A, 600, A30 Claret, A. &amp; Bloemen, S. 2011, A&amp;A, 529, A75 Cloutier, R. &amp; Menou, K. 2020, The Astronomical Journal, 159, 211 Collins, K. A., Kielkopf, J. F., Stassun, K. G., &amp; Hessman, F. V. 2017, aj, 153, 77 Cosentino, R., Lovis, C., Pepe, F., et al. 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 8446, Ground-based and Airborne Instrumentation for Astronomy IV, ed. I. S. McLean, S. K. Ramsay, &amp;H. Takami, 84461V Cushing, M. C., Rayner, J. T., &amp; Vacca, W. D. 2005, ApJ, 623, 1115 Cushing, M. C., Vacca, W. D., &amp; Rayner, J. T. 2004, PASP, 116, 362 Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog, II 246 / Cutri, R. M., Wright, E. L., Conrow, T., et al. 2013, Explanatory Supplement to the AllWISE Data Release Products, Explanatory Supplement to the AllWISE Data Release Products Dekany, R., Roberts, J., Burruss, R., et al. 2013, ApJ, 776, 130 Delrez, L., Murray, C. A., Pozuelos, F. J., et al. 2022, A&amp;A, 667, A59 Demory, B. O., Pozuelos, F. J., Gómez Maqueo Chew, Y., et al. 2020, A&amp;A, 642, A49 Dévora-Pajares, M. &amp; Pozuelos, F. J. 2022, MATRIX: Multi-phAse Transits Recovery from Injected eXoplanets, Zenodo Dévora-Pajares, M., Pozuelos, F. J., Thuillier, A., et al. 2024, MNRAS[ arXiv:2407.14602 ] Donati, J.-F., Kouach, D., Lacombe, M., et al. 2018, in Handbook of Exoplanets, ed. H. J. Deeg &amp; J. A. Belmonte, 107 Dotter, A. 2016, ApJS, 222, 8 Douglas, S. T., Agüeros, M. A., Covey, K. R., et al. 2014, ApJ, 795, 161

TOI-6002b

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-016ba8b1-0bbb-40d7-8188-bd27ec68de13.jpg)

TOI-5713b

Fig. 11. Simulated JWST transmission spectra for TOI-6002b and TOI-5713b. The top panel shows spectra for a cloudless, 1-bar Earth-like atmosphere. The bottom panel displays spectra for a 92-bar Venus-like atmosphere, both cloudless (black) and with a cloud deck at 10 mbar (green). Error bars are shown for a transit observation simulated at R = 30.

Eastman, J. D., Rodriguez, J. E., Agol, E., et al. 2019, arXiv e-prints, arXiv:1907.09480

Jordan, S., Rimmer, P. B., Shorttle, O., &amp; Constantinou, T. 2021, The Astrophysical Journal, 922, 44

- Ehrenreich, D., Vidal-Madjar, A., Widemann, T., et al. 2012, Astronomy &amp; Astrophysics, 537, L2
- Engle, S. G. &amp; Guinan, E. F. 2023, ApJ, 954, L50
- Fulton, B. J., Petigura, E. A., Howard, A. W., et al. 2017, aj, 154, 109
- Furlan, E., Ciardi, D. R., Everett, M. E., et al. 2017, AJ, 153, 71
- Gaia Collaboration. 2020, VizieR Online Data Catalog, I 350 /

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&amp;A, 616, A1

Ghachoui, M., Soubkiou, A., Wells, R. D., et al. 2023, A&amp;A, 677, A31

Giacalone, S., Dressing, C. D., Jensen, E. L. N., et al. 2021, AJ, 161, 24

- Gillon, M., Anderson, D. R., Collier-Cameron, A., et al. 2013, aap, 552, A82 Ginzburg, S., Schlichting, H. E., &amp; Sari, R. 2016, The Astrophysical Journal, 825, 29

Ginzburg, S., Schlichting, H. E., &amp; Sari, R. 2018, Monthly Notices of the Royal Astronomical Society, 476, 759

Gizis, J. E. 1997, AJ, 113, 806

Gomez Barrientos, J., Kaltenegger, L., &amp; MacDonald, R. J. 2023, MNRAS, 524, L5

- Hamuy, M., Suntze ff , N. B., Heathcote, S. R., et al. 1994, PASP, 106, 566
- Hamuy, M., Walker, A. R., Suntze ff , N. B., et al. 1992, PASP, 104, 533
- Hayward, T. L., Brandl, B., Pirger, B., et al. 2001, PASP, 113, 105
- Hedges, C., Angus, R., Barentsen, G., et al. 2020, Research Notes of the American Astronomical Society, 4, 220
- Henry, T. J., Jao, W. C., Winters, J. G., et al. 2019, VizieR Online Data Catalog, J AJ 155 265 / / /
- Hippke, M., David, T. J., Mulders, G. D., &amp; Heller, R. 2019, The Astronomical Journal, 158, 143

Hippke, M., David, T. J., Mulders, G. D., &amp; Heller, R. 2019, AJ, 158, 143 Hippke, M. &amp; Heller, R. 2019, A&amp;A, 623, A39

Jehin, E., Gillon, M., Queloz, D., et al. 2011, The Messenger, 145, 2

Jenkins, J. M. 2002, ApJ, 575, 493

- Jenkins, J. M., Tenenbaum, P., Seader, S., et al. 2020, Kepler Data Processing Handbook: Transiting Planet Search, Kepler Science Document KSCI-19081-
- 003

Jenkins, J. M., Twicken, J. D., McCauli ff , S., et al. 2016, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 9913, Software and Cyberinfrastructure for Astronomy IV, ed. G. Chiozzi &amp; J. C. Guzman, 99133E

Jin, S., Mordasini, C., Parmentier, V., et al. 2014, The Astrophysical Journal, 795, 65

- Kaltenegger, L. 2017, ARA&amp;A, 55, 433
- Kaltenegger, L., Lin, Z., &amp; Rugheimer, S. 2020, ApJ, 904, 10

Kaltenegger, L., Payne, R., Lin, Z., Kasting, J., &amp; Delrez, L. 2023a, Monthly

Notices of the Royal Astronomical Society: Letters, 524, L10

- Kaltenegger, L., Payne, R. C., Lin, Z., Kasting, J., &amp; Delrez, L. 2023b, Monthly Notices of the Royal Astronomical Society: Letters, 524, L10
- Kane, S. R., Kopparapu, R. K., &amp; Domagal-Goldman, S. D. 2014, The Astrophysical Journal Letters, 794, L5

Kilic, M., Munn, J. A., Harris, H. C., et al. 2017, ApJ, 837, 162

Kovács, G., Zucker, S., &amp; Mazeh, T. 2002, A&amp;A, 391, 369

Lasker, B. M., Doggett, J., McLean, B., et al. 1996, in Astronomical Society of the Pacific Conference Series, Vol. 101, Astronomical Data Analysis Software and Systems V, ed. G. H. Jacoby &amp; J. Barnes, 88

Lee, E. J. &amp; Connors, N. J. 2021, The Astrophysical Journal, 908, 32

Lépine, S., Rich, R. M., &amp; Shara, M. M. 2003, AJ, 125, 1598

Lépine, S., Rich, R. M., &amp; Shara, M. M. 2007, ApJ, 669, 1235

Li, C. &amp; Zhao, G. 2017, ApJ, 850, 25

Li, J., Tenenbaum, P., Twicken, J. D., et al. 2019, PASP, 131, 024506

Lindegren, L., Bastian, U., Biermann, M., et al. 2021, A&amp;A, 649, A4

Lopez, E. D. &amp; Fortney, J. J. 2013, The Astrophysical Journal, 776, 2

- Lopez, E. D. &amp; Rice, K. 2018, Monthly Notices of the Royal Astronomical Society, 479, 5303
- Luque, R. &amp; Pallé, E. 2022, Science, 377, 1211
- Malik, M., Grosheintz, L., Mendonça, J. M., et al. 2018, Astrophysics Source Code Library, ascl
- Mann, A. W., Brewer, J. M., Gaidos, E., Lépine, S., &amp; Hilton, E. J. 2013a, AJ, 145, 52
- Mann, A. W., Brewer, J. M., Gaidos, E., Lépine, S., &amp; Hilton, E. J. 2013b, AJ, 145, 52
- Mann, A. W., Deacon, N. R., Gaidos, E., et al. 2014, AJ, 147, 160
- Mann, A. W., Dupuy, T., Kraus, A. L., et al. 2019, apj, 871, 63
- Mann, A. W., Feiden, G. A., Gaidos, E., Boyajian, T., &amp; von Braun, K. 2015, apj, 804, 64
- Martín, E. L., Delfosse, X., Basri, G., et al. 1999, AJ, 118, 2466
- McCully, C., Volgenau, N. H., Harbeck, D.-R., et al. 2018, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 10707, Software and Cyberinfrastructure for Astronomy V, ed. J. C. Guzman &amp; J. Ibsen, 107070K
- Miller, J. S. &amp; Stone, R. P. S. 1994, The Kast Double Spectrograph, Tech. Rep. 66, University of California Lick Observatory Technical Reports

Minkowski, R. L. &amp; Abell, G. O. 1963, in Basic Astronomical Data: Stars and Stellar Systems, ed. K. A. Strand, 481

- Mollière, P., Wardenier, J. P., van Boekel, R., et al. 2019, A&amp;A, 627, A67
- Mugrauer, M. &amp; Michel, K.-U. 2020, Astronomische Nachrichten, 341, 996

Mugrauer, M. &amp; Michel, K.-U. 2021, Astronomische Nachrichten, 342, 840

- Narita, N., Fukui, A., Kusakabe, N., et al. 2019, Journal of Astronomical Telescopes, Instruments, and Systems, 5, 015001
- Narita, N., Fukui, A., Yamamuro, T., et al. 2020, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 11447, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 114475K

Ostberg, C., Kane, S. R., Li, Z., et al. 2023, The Astronomical Journal, 165, 168 Owen, J. E. &amp; Wu, Y. 2013, The Astrophysical Journal, 775, 105

- Parviainen, H. 2015, MNRAS, 450, 3233
- Parviainen, H., Luque, R., &amp; Palle, E. 2024, MNRAS, 527, 5693
- Parviainen, H., Tingley, B., Deeg, H. J., et al. 2019, A&amp;A, 630, A89
- Pass, E. K., Charbonneau, D., Latham, D. W., et al. 2024, arXiv e-prints, arXiv:2401.10167
- Payne, R. C. &amp; Kaltenegger, L. 2024, MNRAS, 530, L13
- Pidhorodetska, D., Moran, S. E., Schwieterman, E. W., et al. 2021, The Astronomical Journal, 162, 169
- Pozuelos, F. J., Suárez, J. C., de Elía, G. C., et al. 2020, A&amp;A, 641, A23
- Pozuelos, F. J., Timmermans, M., Rackham, B. V., et al. 2023, A&amp;A, 672, A70 Quirrenbach, A., CARMENES Consortium, Amado, P. J., et al. 2020, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 11447, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 114473C
- Rayner, J. T., Cushing, M. C., &amp; Vacca, W. D. 2009, ApJS, 185, 289
- Rayner, J. T., Toomey, D. W., Onaka, P. M., et al. 2003, PASP, 115, 362
- Reid, I. N., Hawley, S. L., &amp; Gizis, J. E. 1995, AJ, 110, 1838
- Reylé, C., Jardine, K., Fouqué, P., et al. 2021, A&amp;A, 650, A201
- Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, Journal of Astronomical Telescopes, Instruments, and Systems, 1, 014003

Riddick, F. C., Roche, P. F., &amp; Lucas, P. W. 2007, MNRAS, 381, 1067

Rojas-Ayala, B., Covey, K. R., Muirhead, P. S., &amp; Lloyd, J. P. 2012, ApJ, 748,

93

- Savitzky, A. &amp; Golay, M. J. E. 1964, Analytical Chemistry, 36, 1627 Schlafly, E. F. &amp; Finkbeiner, D. P. 2011, ApJ, 737, 103
- Schönrich, R., Binney, J., &amp; Dehnen, W. 2010, MNRAS, 403, 1829
- Seifahrt, A., Bean, J. L., Stürmer, J., et al. 2020, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 11447, Society of

Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 114471F

Smith, J. C., Stumpe, M. C., Van Cleve, J. E., et al. 2012, PASP, 124, 1000

- Stassun, K. G., Oelkers, R. J., Pepper, J., et al. 2018, AJ, 156, 183
- Stumpe, M. C., Smith, J. C., Catanzarite, J. H., et al. 2014, PASP, 126, 100
- Stumpe, M. C., Smith, J. C., Van Cleve, J. E., et al. 2012, pasp, 124, 985
- Sullivan, P. W., Winn, J. N., Berta-Thompson, Z. K., et al. 2015, ApJ, 809, 77
- Tenenbaum, P., Christiansen, J. L., Jenkins, J. M., et al. 2012, The Astrophysical Journal Supplement Series, 199, 24
- Tsai, S.-M., Lyons, J. R., Grosheintz, L., et al. 2017, ApJS, 228, 20 Turbet, M., Ehrenreich, D., Lovis, C., Bolmont, E., &amp; Fauchez, T. 2019, A&amp;A, 628, A12
- Twicken, J. D., Catanzarite, J. H., Clarke, B. D., et al. 2018, PASP, 130, 064502 Van Eylen, V., Agentoft, C., Lundkvist, M. S., et al. 2018, MNRAS, 479, 4786 Van Eylen, V., Astudillo-Defru, N., Bonfils, X., et al. 2021, MNRAS, 507, 2154 Venturini, J., Guilera, O. M., Haldemann, J., Ronco, M. P., &amp; Mordasini, C. 2020, A&amp;A, 643, L1

Venturini, J., Ronco, M. P., Guilera, O. M., et al. 2024, A&amp;A, 686, L9 Wells, R. D., Rackham, B. V., Schanche, N., et al. 2021, A&amp;A, 653, A97 West, A. A., Hawley, S. L., Bochanski, J. J., et al. 2008, AJ, 135, 785 Wizinowich, P., Acton, D. S., Shelton, C., et al. 2000, PASP, 112, 315 Zacharias, N., Finch, C. T., Girard, T. M., et al. 2012, VizieR Online Data Catalog, I 322A /

Ziegler, C., Tokovinin, A., Briceño, C., et al. 2020, AJ, 159, 19

- 1 Astrobiology Research Unit, Université de Liège, Allée du 6 Août 19C, B-4000 Liège, Belgium
- 2 Oukaimeden Observatory, High Energy Physics and Astrophysics Laboratory, Cadi Ayyad University, Marrakech, Morocco
- 3 Department of Earth, Atmospheric and Planetary Science, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA
- 4 Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- 5 Dpto. Física Teórica y del Cosmos, Universidad de Granada, 18071, Granada, Spain
- 6 Avature Machine Learning, Spain
- 7 Carl Sagan Institute at Cornell University, Space Science Building, Ithaca, NY 14850, USA
- 8 School of Physics &amp; Astronomy, University of Birmingham, Edgbaston, Birmimgham B15 2TT, UK
- 9 Instituto de Astrofísica de Andalucía (IAA-CSIC), Glorieta de la Astronomía s n, 18008 Granada, Spain /
- 10 Center for Astrophysics | Harvard &amp; Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA
- 11 Department of Astronomy &amp; Astrophysics, University of California San Diego, La Jolla, CA 92093, USA
- 12 Instituto de Astrofísica de Canarias (IAC), Calle Vía Láctea s n, / 38200, La Laguna, Tenerife, Spain
- 13 Departamento de Astrofísica, Universidad de La Laguna (ULL), 38206 La Laguna, Tenerife, Spain
- 14 Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA
- 15 Space Sciences, Technologies and Astrophysics Research (STAR) Institute, Université de Liège, Allé du 6 Août 19C, B-4000 Liège, Belgium
- 16 NASA Ames Research Center, Mo ff ett Field, CA 94035, USA
- 17 Komaba Institute for Science, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan
- 18 Astrobiology Center, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan
- 19 Department of Multi-Disciplinary Sciences, Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo 153-8902, Japan
- 20 National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan
- 21 Department of Physics and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA02139, USA
- 22 Royal Astronomical Society, Burlington House, Piccadilly, London W1J 0BQ, UK
- 23 NASAGoddard Space Flight Center, 8800 Greenbelt Rd, Greenbelt, MD20771, USA
- 24 Crow Observatory, Portalegre, Portugal
- 25 Bryant Space Science Center, Department of Astronomy, University of Florida, Gainesville, FL 32611, USA
- 26 Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
- 27 Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109 USA
- 28 NASA Exoplanet Science Institute, IPAC, California Institute of Technology, Pasadena, CA 91125 USA
- 29 Department of Physics and Astronomy, University of Kansas, Lawrence, KS 66045, USA
- 30 Department of Astronomy, University of California Berkeley, Berkeley, CA 94720, USA
- 31 Department of Physics, Center for Astrophysics &amp; Space Sciences, University of California San Diego, La Jolla, CA 92093, USA
- 32 Department of Physics and Astronomy, University of Notre Dame, South Bend, IN 46556, USA
- 33 AIM, CEA, CNRS, Université Paris-Saclay, Université de Paris, F91191 Gif-sur-Yvette, France
- 34 Paris Region Fellow, Marie Sklodowska-Curie Action
- 35 Cavendish Laboratory, JJ Thomson Avenue, Cambridge, CB3 0HE, UK

clearpage

## Appendix A: Additional figures

/s)

Flux ×10 3  (e

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-ecf1b831-e90f-411a-9c7d-376845fa615d.jpg)

TIC 102734241 - Sector 54

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-e4aafd9b-40d0-4471-a303-1cff098a8a8b.jpg)

TIC 219041246 - Sector 16

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-b8c54c5b-df74-48f3-9088-8431b9d78577.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-d8c05a40-4e54-4a65-bd9e-2a01db8fd5f6.jpg)

TIC 102734241 - Sector 55

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-da347add-432f-4a9b-8c06-21f13baca24e.jpg)

TIC 219041246 - Sector 22

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-3dc680e0-ed3c-46af-96e7-d9b2e38d2d37.jpg)

TIC 219041246 - Sector 23

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-7f166c3f-4972-4ed1-9ad3-41bcd46a0b53.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-4db90473-52bd-472b-9501-11a72145bea1.jpg)

TIC 219041246 - Sector 49

Fig. A.1. Target pixel files (TPFs) of TOI-6002 (TIC 102734241) in TESS sectors 14, 41, 54, 55, and TOI-5713 (TIC 219041246) in sectors 16, 22, 23, and 49 created with tpfplotter (Aller et al. 2020). Regions in orange highlight the aperture used to extract the SPOC photometry. Red circles represent the Gaia DR2 sources, with sizes representing the magnitude contrasts with respect to the targets.

Fig. A.2. Detrended TESS photometry of TOI-6002 b and TOI-5713 b. The flares in the TOI-573 light curves were removed using an upper 3 σ clipping. Gray points show the detrended PDC-SAP fluxes from the SPOC pipeline. Orange and red points highlight the transits of TOI-6002 b and TOI-5713 b, respectively

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-4f651976-5765-40ec-90fd-c0bb0209d480.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00709v4-figure-339a199e-e137-4cc1-9a87-8907fc786727.jpg)

.
M. Ghachoui B. V. Rackham M. Dévora-Pajares J. Chouqar M. Timmermans L. Kaltenegger D. Sebastian F. J. Pozuelos J. D. Eastman A. J. Burgasser F. Murgas K. G. Stassun M. Gillon Z. Benkhaldoun E. Palle L. Delrez J. M. Jenkins K. Barkaoui N. Narita J. P. de Leon M. Mori A. Shporer P. Rowden V. Kostov G. Fűrész K. A. Collins R. P. Schwarz D. Charbonneau N. M. Guerrero G. Ricker E. Jehin A. Fukui Y. Kawai Y. Hayashi E. Esparza-Borges H. Parviainen C. A. Clark D. R. Ciardi A. S. Polanski J. Schleider E. A. Gilbert I. J. M. Crossfield T. Barclay C. D. Dressing P. R. Karpoor E. Softich R. Gerasimov F. Davoudi
TESS discovery of two super-Earths orbiting the M-dwarf stars TOI-6002 and TOI-5713 near the radius valley
We present the validation of two TESS super-Earth candidates transiting the
mid-M dwarfs TOI-6002 and TOI-5713 every 10.90 and 10.44 days, respectively.
The first star (TOI-6002) is located $32.038\pm0.019$ pc away, with a radius of
$0.2409^{+0.0066}_{-0.0065}$ \rsun, a mass of $0.2105^{+0.0049}_{-0.0048}$
\msun, and an effective temperature of $3229^{+77}_{-57}$ K. The second star
(TOI-5713) is located $40.946\pm0.032$ pc away, with a radius of
$0.2985^{+0.0073}_{-0.0072}$ \rsun, a mass of $0.2653\pm0.0061$ \msun, and an
effective temperature of $3225^{+41}_{-40}$ K. We validated the planets using
TESS data, ground-based multi-wavelength photometry from many ground-based
facilities, as well as high-resolution AO observations from Keck/NIRC2.
TOI-6002 b has a radius of $1.65^{+0.22}_{-0.19}$ \re\ and receives
$1.77^{+0.16}_{-0.11} S_\oplus$. TOI-5713 b has a radius of
$1.77_{-0.11}^{+0.13} \re$ and receives $2.42\pm{0.11} S_\oplus$. Both planets
are located near the radius valley and near the inner edge of the habitable
zone of their host stars, which makes them intriguing targets for future
studies to understand the formation and evolution of small planets around
M-dwarf stars.