Open Access
Issue
A&A
Volume 666, October 2022
Article Number A4
Number of page(s) 19
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202243381
Published online 27 September 2022

© K. Kawauchi et al. 2022

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.

1 Introduction

Small planets with radii between 1 and 4 R are extremely common in the Milky Way, but do not exist in our solar system. NASA’s Kepler space telescope revealed that this population has a bimodal distribution that separates the planets into super-Earth (1–1.5R) and sub-Neptune (2–4R) populations (Fulton et al. 2017) for planets with P < 100 days. This “radius valley” is consistent with the transition from rocky to non-rocky planets and is predicted by photoevaporation models (e.g., Owen & Wu 2013, 2017; Jin et al. 2014; Lopez & Rice 2018; Mordasini 2020) and/or core-powered mass-loss models (Ginzburg et al. 2016, 2018; Gupta & Schlichting 2019, 2021).

Furthermore, the location of the radius valley is revealed to be dependent on the planetary period, insolation, and stellar type (e.g., Van Eylen et al. 2018; Martinez et al. 2019; Wu 2019). The position of the radius gap in the radius-period distribution plane around Sun-like stars decreases with the orbital period following a power law RgapP−0.11, which is consistent with photoevaporation and core-powered mass-loss model predictions (Martinez et al. 2019). On the other hand, it has been claimed that the position of the radius valley around late-type stars increases with the orbital period following a power law RgapP+0.058, which is consistent with gas-poor formation models (Cloutier & Menou 2020). However, a recent study found that the correlation between the position of the radius gap and orbital period is negative around M-type stars whose masses are above 0.5 M (RgapP−0.12), and this correlation does not seem to vary with the stellar type (Petigura et al. 2022). To further explore the nature and origins of the radius valley around M dwarfs, it is important to discover new small planets and to measure their masses as well as radii, in order to have a more complete sample.

Sub-Neptune-sized planets are likely to have H2O-dominated ices and/or fluids (water worlds) in addition to a rocky core enveloped in a H2-He gas (e.g., Zeng et al. 2019). Recent studies indicated the possibility that sub-Neptune-sized planets could have large oceans with habitable conditions under H2-rich atmospheres (Hycean worlds; Madhusudhan et al. 2021). The direct observation of their atmospheres helps our understanding of the composition of the planets.

The infrared helium triplet is a good indicator of a primary atmosphere accreted from the protoplanetary nebula (Elkins-Tanton & Seager 2008). Previous studies have detected the He i 10830 Å absorption lines in the atmosphere of about a half-dozen Jupiter- and Neptune-sized planets (e.g., Spake et al. 2018; Nortmann et al. 2018; Allart et al. 2018), but only in three sub-Neptunes: GJ3470b (Palle et al. 2020), GJ 1214b (Orell-Miquel et al. 2022), and TOI 560.01 (Zhang et al. 2022). Other studies of sub-Neptunes have been able to place only upper limits on He i absorption, but they still help to constrain the mass-loss rate (e.g., Krishnamurthy et al. 2021).

Here we report the discovery of TOI-2136b, a sub-Neptune-sized planet around an M dwarf, using photometry from NASA’s Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015) and ground-based follow-up observations. We also search for helium in its atmosphere with high-resolution transmission spectroscopy. The paper is organized as follows. In Sects. 2 and 3, we introduce the observation of TESS and ground-based telescope. In Sects. 4, 5, 6, and 7, we describe the analysis and results for stellar parameters, planet parameters, and transmission spectroscopy. We discuss the physical properties in Sect. 8, and we summarize this work in Sect. 9.

thumbnail Fig. 1

TESS target pixel file image of TOI-2136 observed in Sector 26 (left) and Sector 40 (right). The pixels in red show the photometric aperture used by TESS to create the PDCSAP light curves. The position of nearby stars and their magnitudes are represented by the red circles. This image was produced using tpfplotter (Aller et al. 2020).

2 TESS photometry

TIC 336128819 (TOI-2136) was observed with a two-minute cadence by TESS in Sector 26 from UT 9 June 2020 to 4 July 2020, and in Sector 40 from UT 24 June 2021 to 23 July 2021. For each sector there are 0.919 days downtime to download data and 23.95 days of science data in total. For both TESS sectors, the target star was positioned on charge-coupled device (CCD) 1 Camera 1. Figure 1 shows the TESS images around TOI-2136 for Sectors 26 and 40.

TESS raw images were processed by the Science Processing Operations Center (SPOC) at NASA Ames Research Center. After the data were reduced by the SPOC pipeline (Jenkins et al. 2016), a preliminary TESS light curve was produced using simple aperture photometry (SAP; Morris et al. 2020) and instrumental systematic effects were removed using the presearch data conditioning (PDC) pipeline module (Smith et al. 2012; Stumpe et al. 2014). A transit search (Jenkins 2002; Jenkins et al. 2020) was carried out on the TESS light curves and a signal of ~7.8 days was detected. After the transit signal was validated by the TESS team (Twicken et al. 2018), it was assigned a TESS object of interest number (TOI-2136.01) and announced to the community.

3 Follow-up observations

3.1 MuSCATs

We observed a total of five full transit events with the Multicolor Simultaneous Camera for studying Atmospheres of Transiting exoplanets (MuSCAT; Narita et al. 2015), MuSCAT2 (Narita et al. 2019) and MuSCAT3 (Narita et al. 2020). MuSCAT is mounted on the 1.88 m telescope of the National Astronomical Observatory of Japan (NAOJ) in Okayama, Japan. MuSCAT2 is mounted on the 1.52 m Telescópio Carlos Sanchez (TCS) at the Teide Observatory, in the Canary Islands, Spain. MuSCAT3 is mounted on the 2 m Faulkes Telescope North at the Haleakala Observatory, in Maui, Hawai’i. These instruments have g, r, i and zs bands (only g, r, and zs bands for MuSCAT). MuSCAT and MuSCAT2 are equipped with 1024 × 1024 pixel CCDs, and MuSCAT3 is equipped with 2048 × 2048 pixel CCDs. MuSCAT provides a field of view (FOV) of 6.1 × 6.1 arcmin2 with a pixel scale of 0.36 arcsec per pixel. MuSCAT2 provides a field of view of 7.4 × 7.4 arcmin2 with a pixel scale of 0.44 arcsec per pixel. MuSCAT3 provides a field of view of 9.1 × 9.1 arcmin2 with a pixel scale of 0.27 arcsec per pixel.

The first transit was observed by MuSCAT3 on UT 12 June 2021, and the exposure times were set to 60, 26, 23, and 20 s for g, r, i, and zs, respectively. The fourth transit was observed by MuSCAT on UT 23 September 2021, and the exposure times were set to 30 and 30 s for r and zs; for that night the g-band CCD of MuSCAT did not work. The rest of the transits were observed by MuSCAT2 on UT 22 and 30 August 2021, and 24 October 2021. For the observations of 22 August 2021, the exposure times were set to 30, 15, and 15 s for g, i, and zs, respectively; for that night the r-band CCD of MuSCAT2 did not work and we also discarded the g-band data due to the large photometric scatter present in the light curve. For the MuSCAT2 observations of 30 August 2021, the exposure times were the same for each band as the night of 22 August 2021; the r-band CCD did not work that night also, but the g-band data were much better and we decided to include them in the analysis. For the night of 24 October 2021, the exposure times were set to 30, 20, 15, and 15 s for the g, r, i, and zs bands, respectively. We used a fixed aperture radius of 4.3 arcsec for each band for the MuSCAT photometry, and 10.9 arcsec for each band for the MuSCAT2 photometry. For the MuSCAT3 photometry, we used fixed aperture radii of 3.7, 4.3, 4.3, and 4.8 arcsec for g, r, i, and zs, respectively.

Table 1

Ground-based photometric observations.

3.2 LCO 1 m/Sinistro

We observed one almost-full transit of TOI-2136b on UT 9 October 2021 simultaneously with two 1 m telescopes of Las Cumbres Observatory (LCO) at the McDonald Observatory in the USA. Each telescope is equipped with a single-band imager Sinistro, which has a 4k × 4k CCD with a pixel scale of 0″.389 pixel−1 and an FOV of 26′.5 × 26′.5. We selected the i-band filter for one telescope and the Z-band filter for the other. We set the exposure time at 30 s with the full-frame readout mode for both telescopes. The telescopes were slightly defocused (by 0.2 mm) to avoid detector saturation. The obtained raw images were processed by the BANZAI pipeline (McCully et al. 2018) to perform dark-image and flat-field corrections. Aperture photometry was then performed using a custom pipeline (Fukui et al. 2011). A summary of the ground-based photometric observations, facilities, instruments, and filters used in this work is presented in Table 1.

3.3 Subaru/IRD

We obtained a total of 38 high-resolution spectra of TOI-2136 from UT 27 September 2020 to 27 October 2021 with the InfraRed Doppler (IRD) instrument on the Subaru 8.2 m telescope under the Subaru-IRD TESS intensive follow-up program (ID: S20B-088I, S21B-118I). The IRD is a fiber-fed, near-infrared (NIR) spectrometer covering from 970 to 1750 nm with a spectral resolution of 70 000 (Tamura et al. 2012; Kotani et al. 2018). The exposure times were set to 750-1800 s depending on the sky condition. We also observed a telluric standard star (A0 or A1 star) to correct for telluric lines in the template spectrum, for the radial velocity (RV) analysis.

Raw IRD data were reduced with IRAF (Tody 1993) and our custom code to correct for the bias pattern of the detectors, and also to perform wavelength calibration with the Th-Ar lamp and laser-frequency comb (Kuzuhara et al. 2018; Hirano et al. 2020). The reduced one-dimensional spectra have signal-to-noise ratios (S/Ns) of 40-94 per pixel around 1000 nm. We measured the precision RV for each frame with the RV analysis pipeline for the IRD described in Hirano et al. (2020). The typical RV internal errors are 3–5 m s−1. The RV individual measurements are given in Table A.1.

The IRD data also include the observation of a full transit of TOI-2136b on 27 September 2020, with the goal of searching for an excess helium absorption by the planetary atmosphere. We obtained 11 consecutive spectra with a 900 s exposure time from UT 05:01 to 07:51, out of which five frames are out-of-transit spectra and six frames are in-transit spectra. The S/Ns of the extracted spectra are 50 – 63 per pixel around 1000 nm, and the airmass varied from 1.05 to 1.33 during the observations.

3.4 Gemini North/‘Alopeke

Exoplanet parameters, such as the radius determination, depend on the transit depth as well as the assumption that the star is single. A spatially close companion (bound or line-of-sight) can create a false-positive transit signal, or cause “third-light” flux, leading to an underestimation of the planet radius (Ciardi et al. 2015). Such companion stars can even cause nondetections of small planets residing within the same exoplanetary system (Lester et al. 2021). Nearly 50% of FGK-type stars are binary or multiple star systems (Matson et al. 2018), providing crucial information toward our understanding of exoplanetary formation, dynamics, and evolution (Howell et al. 2021). Thus, to search for close-in (<1″) companions unresolved in TESS images or other ground-based follow-up observations, we obtained high-resolution imaging speckle observations of TOI-2136.

TOI-2136 was observed on UT 17 October 2021 using the ‘Alopeke speckle instrument on the Gemini North 8-m telescope (Scott et al. 2021). ‘Alopeke provides simultaneous speckle imaging in two bands (562 nm and 832 nm) with output data products including a reconstructed image with robust contrast limits on companion detections (e.g., Howell et al. 2016). Five sets of 1000 × 0.06 s exposures were collected for TOI-2136 and subjected to Fourier analysis in our standard reduction pipeline (see Howell et al. 2011). Figure 2 shows our final contrast curves and the 832 nm reconstructed speckle image. We find that TOI-2136 is a single star with no companion dimmer than four to seven magnitudes above that of the target star from the diffraction limit (20 mas) out to 1.2. At the distance of TOI-2136 (d = 33 pc), these angular limits correspond to spatial limits of 0.7–40 au.

4 Stellar properties

4.1 Stellar parameters

We determined the effective temperature Teff and metallicity of the host star from the template spectrum, which was constructed from the IRD spectra. The analysis is based on a line-by-line comparison of equivalent widths (EWs) between the observed spectra and synthetic spectra, following the procedure described in Sect. 3.1.2 of Hirano et al. (2021). The synthetic spectra were calculated by a one-dimensional LTE spectral synthesis code, which was based on the same assumptions as the model atmosphere program of Tsuji (1978), with the interpolated grid of MARCS model atmosphere (Gustafsson et al. 2008). Through spectral analysis, a surface gravity log g value of 4.88 ± 0.0008 was adopted from the TESS Input Catalog (TIC) version 8 (Stassun et al. 2019), and the microturbulent velocity was fixed at 0.5 ± 0.5 km s−1 for simplicity.

The selected 47 FeH molecular lines in the Wing-Ford band at 990–1020 nm were used for the Teff estimation. Ishikawa et al. (2022) provides the selection criteria and the list of the 47 lines. The original data of these line transitions are available on the MARCS1 web page. We measured the EW of each FeH line by fitting the Gaussian profile to find the Teff at which the synthetic spectra best reproduce the EW by an iterative search.

For the stellar metallicity, we used the analysis tool developed by Ishikawa et al. (2020). We used the atomic absorption lines of Na, Mg, Ca, Ti, Cr, Mn, Fe, and Sr based on the selection criteria of: (1) not suffering from blending of other absorption lines, (2) sensitive to elemental abundances, and (3) continuum level can be reasonably determined. The original data of the atomic lines are taken from the Vienna Atomic Line Database (VALD; Kupka et al. 1999; Ryabchikova et al. 2015). The EWs were measured by fitting the synthetic spectra on a line-by-line basis. We searched for an elemental abundance until the synthetic EW matched the observed one for each line, and took the average for all the lines to estimate [X/H] for an element X. The [M/H] is calculated as an error-weighted average of abundances of all the eight elements.

As a first step, we estimated Teff assuming the solar metallicity. The average of the Teff estimates for the 47 lines was adopted as the best estimate. We temporarily determined the elemental abundances using the Teff estimate and, as a second step, we redetermined Teff adopting the iron abundance [Fe/H] as the metallicity of the atmospheric model grid. Its uncertainty was given as the sum of the statistical error (8 K), which was calculated by dividing the standard deviation of the 47 lines by the square root of the number of lines (σ/N$\sigma {\rm{/}}\sqrt N $), and the possible systematic error (100 K). We then adopted the resulting Teff to finally determined the elemental abundances. The procedure up to this point allowed the results of Teff and abundances to converge well within the measurement errors. Consequently, we obtained Teff = 3373 ± 108 K, [Fe/H] = 0.02 ± 0.14 dex, and [M/H] = 0.06 ± 0.06 dex.

Based on these parameters, as well as the literature values in Table 2, we estimated the other physical parameters, including the stellar radius, which is required to infer the planet radius. Making use of the empirical relations for M dwarfs by Mann et al. (2015, Mann et al. 2019), we derived the stellar mass M*, radius R*, and density ρ*. with their uncertainties estimated by the Monte Carlo approach using Gaussian distributions for the Gaia parallax, 2MASS magnitudes, and [Fe/H]. The stellar luminosity L* was also derived from the Stefan-Boltzmann law using Teff in Table 2.

thumbnail Fig. 2

Gemini/‘Alopeke 5-sigma contrast curves in 562 nm (blue) and 832 nm (red). The inset figure is the 832 nm reconstructed speckle image.

Table 2

Stellar parameters.

4.2 Stellar rotation

We obtained publicly available photometry of TOI-2136 from the MEarth project (Berta et al. 2012) and the Zwicky Transient Facility (ZTF, Masci et al. 2019). The MEarth project uses eight 40 cm telescopes equipped with RG715 filters and it mostly focuses on observing M stars to search for transits (Nutzman & Charbonneau 2008). The project has two observing sites in both hemispheres. One site is located at the Fred Lawrence Whipple Observatory, on Mount Hopkins in Arizona, in the USA. The other is located at the Cerro Tololo Inter-American Observatory (CTIO) in Chile. From the MEarth project, we obtained TOI-2136 (LSPMJ1844+3633) photometry from the 2011 to 2020 north target light curves.

The ZTF uses a 47 squared degree camera mounted on the Palomar 48-inch Schmidt Telescope to search and study transient and variable objects in the sky. Using its large observing area, ZTF is capable of scanning the northern hemisphere sky every two nights. From ZTF we obtained g– and r–band photometry of TOI-2136 from the ZTF data release 8 (DR8).

Visual inspection of the MEarth and ZTF data shows that TOI-2136 presented periodic photometric variations (see Fig. 3). These variations can be produced by inhomogeneities in the stellar surface (e.g., spots, plages) that appear and disappear out of sight as the star rotates. We performed a joint fit to the MEarth and ZTF photometry to establish the rotation period of the star. We used the package for Gaussian Processes Celerite (Foreman-Mackey et al. 2017) and fit the photometry with the following kernel: kij Phot=B2+Ce| titj |/L[ cos(2π| titj |Prot)+(1+C) ],${k_{i\,j\,{\rm{Phot}}}} = {B \over {2 + C}}{{\rm{e}}^{ - \left| {{t_i} - {t_j}} \right|/L}}\left[ {\cos \left( {{{2\pi \left| {{t_i} - {t_j}} \right|} \over {{P_{{\rm{rot}}}}}}} \right) + \left( {1 + C} \right)} \right],$(1)

where |titj| is the difference between two epochs or observations, B, C, and L are positive constants, and Prot is the stellar rotational period (see Foreman-Mackey et al. 2017 for details). For the joint fit, each data set had the constants B, C, and L as free parameters, but shared a common rotational period. The prior functions and parameter limits used in the joint fit are in Table B.1. After global optimization of a log likelihood function, we sampled the posterior distribution of the kernel used to fit the photometry with emcee (Foreman-Mackey et al. 2013), using 80 chains and 20 000 iterations. The final parameter values (median and 1σ uncertainties) were estimated from the posterior distribution.

Figure 3 shows the MEarth and ZTF photometry, and the best fitted model using the kernel described in Eq. (1). The generalized Lomb–Scargle (GLS, Zechmeister & Kürster 2009) periodograms for each photometric data set are shown in the right panel. The three data sets present peaks in their respective periodograms for periods between 70 and 90 days. The stellar rotational period found by the best fitted model is Prot = 82.56 ± 0.45 days. This rotation period is in agreement with Newton et al. (2016), which found a rotational period of Prot = 82.97 days for this star, using MEarth data.

Engle & Guinan (2018) investigated the relationship between the stellar age and rotation of M dwarfs. Using their equation, the stellar age of TOI-2136 is calculated as ~5 Gyr. We also estimated the stellar age from its correlation with the Galactic space velocities (U, V, W). We computed the UVW velocities using the Gaia Early Data Release (EDR) 3 information and the absolute RV measured from the IRD spectra using the procedure described in Johnson & Soderblom (1987; Table 2). Following Hirano et al. (2021), we obtained an age of 6.23.1+3.7$6.2_{ - 3.1}^{ + 3.7}$ Gyr when the age prior based on the Geneva-Copenhagen Survey catalog (Casagrande et al. 2011) is imposed, and 6.34.1+5.0$6.3_{ - 4.1}^{ + 5.0}\quad $ Gyr for the uniform age prior (0 < age < l4 Gyr). These results indicate that the stellar rotation and UVW velocities do not constrain the age of the star.

5 Methods

5.1 Transit light curve modeling

For the transit and RV modeling, we followed the procedure described in Murgas et al. (2021). The light curves were modeled using the python package PyTransit2 (Parviainen 2015). To fit the transit events analyzed here, we used PyTransit’s Mandel & Agol (2002) analytic models, assuming a quadratic limb-darkening law.

Before performing the joint fit of all the data available, the ground-based transit observations were fitted simultaneously using a common transit model and for each night the systematic effects were accounted for independently. To model the instrumental noise, we used a linear model with one free parameter dependent on the airmass of the target, two free parameters dependent on the star position in the detector (X- and Y-axis), and a term dependent on the full width at half maximum (FWHM) of the point spread function (PSF) as a proxy for seeing variations. For the final joint fit we used the lCo, MuSCAT, MuSCAT2 and MuSCAT3 detrended light curves (i.e., without systematic effects).

For the joint fit, the common transit free parameters used in the transit modeling were the planet-to-star radius ratio Rp/R*, the central time of the transit Tc, the stellar density ρ*, and the transit impact parameter b. A constant orbital period (i.e., assuming no transit timing variations) was used as a common parameter for the transit light curves and RV measurements. The passband dependent quadratic limb-darkening (LD) coefficients u1 and u2 were set free. During the fitting process, we converted the LD coefficients (u1, u2) to the parametrization proposed by Kipping (2013), (q1,q2) with qi ∈ [0,1]. The fitted LD coefficients were weighted against the predicted coefficients delivered by the Python Limb Darkening Toolkit3 (ldtk, Parviainen & Aigrain 2015), computed using the TOI-2136 stellar parameters from Table 2.

To model the stellar variability and residual systematic noise present in the TESS data, we used Gaussian processes (GPs; e.g., Rasmussen & Williams 2010, Gibson et al. 2012, Ambikasaran et al. 2015). For the GPs used in the TESS data, we chose the Matern 3/2 kernel implementation of Celerite (Foreman-Mackey et al. 2017): kij TESS=c12(1+3| titj |τ1)exp(3| titj |τ1),${k_{i\,j\,{\rm{TESS}}}} = c_1^2\left( {1 + {{\sqrt 3 \left| {{t_i} - {t_j}} \right|} \over {{\tau _1}}}} \right)\exp \left( { - {{\sqrt 3 \left| {{t_i} - {t_j}} \right|} \over {{\tau _1}}}} \right),$(2)

where |titj| is the difference between two epochs or observations, c1 is the amplitude of the flux variation, and τ1 is a characteristic timescale. For each TESS data set available, the constants c1 and τ1 were set as free parameters in the fit.

5.2 Radial velocity modeling

The RV data were modeled using the Radial Velocity Modeling Toolkit RadVel4 (Fulton et al. 2018). The free parameters used to model the RVs were: the planet-induced RV semiamplitude (KRV), the host star systemic velocity (γ0), the orbital eccentricity and argument of the periastron (using the parametrization of e$\sqrt e $ cos ω, e$\sqrt e $ sin ω), and the instrumental RV jitter (σRV jitter). The orbital period and central transit time were also set free, but were taken to be global parameters in common with the light curves. To account for systematic noise present in the RV time series, we used GPs with an exponential squared kernel (i.e., a Gaussian kernel): kij RV=c22exp((titj)2τ22),${k_{i\,j\,{\rm{RV}}}} = c_2^2\exp \left( { - {{{{\left( {{t_i} - {t_j}} \right)}^2}} \over {\tau _2^2}}} \right),$(3)

where titj is the difference between two epochs or observations, c2 is the amplitude of the exponential squared kernel, and τ2 is a characteristic timescale. The constants c2 and τ2 were set as free parameters.

We tested how the fitted parameter values changed by fitting a single-planet Keplerian to the data without taking into account the red noise in the RV measurements, and by modeling the sys-tematics using GPs. We then computed the model comparison metrics Bayesian information criterion (BIC, Schwarz 1978) and Akaike information criterion (AIC, Akaike 1974), defined by BIC=klnn2lnmax,${\rm{BIC}} = k\ln n - 2\ln {{\cal L}_{\max }},$(4)

and AIC=2k2lnmax.$AIC = 2k - 2ln{{\cal L}_{\max }}.$(5)

Where n is the number of observed data points, k is the number of fitted parameters, and max is the maximized likelihood. We found that the RV model that includes GPs is slightly preferred, although the difference between the criterion values for the RV model including GPs and without GPs (ΔBIC and ΔAIC) is not large enough to select the single Keplerian model plus the GPs model according to the criteria of Raftery (1995) (ΔBIC < 2.0 and ΔAIC < 2.5).

thumbnail Fig. 3

Ground-based long-term photometric observations of TOI-2136 using MEarth and ZTF, and their respective GLS periodograms (Zechmeister & Kürster 2009) for each data set. The stellar photometric variability is detected with a stronger signal in MEarth data. Fitting a kernel with a periodic term to all the data sets using Gaussian processes, we find a stellar rotation period of 82.56 ± 0.45 days.

5.3 Joint fit

For the joint fit of all the data available, we employed a Bayesian approach. We performed an uninformative transit search in the TES S time series (from Sectors 26 and 40) using Transit Least Squares (TLS, Hippke & Heller 2019). From TLS we obtained an estimate of the period and epoch of the central time of the transit with their respective uncertainties; we used these values as priors for the joint fit. Then, we implemented a Markov chain Monte Carlo (MCMC) to explore the posterior distribution of the parameters using emcee (Foreman-Mackey et al. 2013). With emcee we evaluated a likelihood plus a prior function for several iterations. The likelihood function was the sum of the log likelihood for each transit time series plus a log likelihood for the star LD coefficients and the log likelihood of the RV observations. A total of 26 free parameters were sampled in the joint fit. The fitting process started with the optimization of the global likelihood with PyDE5, and once the optimization process ended we used the optimal parameter distribution to start the MCMC. The MCMC consisted of 200 chains and ran for 2000 iterations as a burn-in, and the main MCMC ran for 8000 iterations.

The final parameter values and their respective uncertainties were computed from the posterior distributions, adopting the median as parameter values and using the 1σ levels of the distribution as uncertainties. The prior functions and parameter limits used in the joint fit are in Table C.1. Figure C.1 shows the correlation plots of the fitted orbital parameters, excluding the limb-darkening coefficients and parameters related to the modeling of red noise present in the data.

thumbnail Fig. 4

TOI-2136 light curves of TESS Sector 26 (top panel) and Sector 40 (bottom panel). The black points represent TESS binned photometry. The gray points are the TESS PDCSAP data. The best fitting transit model model including systematic effects is shown in dark blue. Light-blue points are the TESS data after correcting the systematic effects. The best fitting transit model is shown in red. Transit events of TOI-2136b are marked with black vertical lines.

6 Results

We present the fitted and derived parameters for the joint fit of the data with and without the use of GPs to model the systematics present in the RV measurements (Table 3). For the case of the fit without GPs to model the red noise, the RV jitter parameter increased by a factor of two when compared to the GP modeling results; an expected result since there was no component to model the red noise in this approach. For the non-GP approach, we find an upper limit for the planetary mass of 8.3 M while for the fit including GPs in the RV modeling, we find Mp < 9.9 M. As discussed in Sect. 5.2, the fits performed with and without using GPs for the RV measurements are not statistically different according to the difference between the model selection metrics estimated using the BIC and AIC definitions. We choose to adopt the results of the fit using GPs, since this approach takes into account the red noise in the data and the results provide a more conservative upper mass limit for the mass of the planet.

Figure 4 presents the photometric light curves from TESS Sectors 26 and 40, and the best fitting model (including systematic effects) found by the joint fit. Figure 5 shows the detrended and phase-folded TESS light curves for both TESS sectors analyzed in this work. Figure 6 shows the multicolor ground-based transit follow-up observations of TOI-2136b used in the joint fit. Finally, Fig. 7 presents the IRD RV measurements and the best fitting model.

Our joint fit finds that TOI-2136b has a period of P = 7.851925 ± 0.000016 days and an orbital eccentricity consistent with 0. We also determined that TOI-2136b has a radius of Rp = 2.20 ± 0.07 R and an equilibrium temperature of Teq = 378 ± 13 K. We note TOI-2136 light curves of TESS Sector 26 have a nontrivial bias by sky background noise correction. As this bias is about 1.54% in the data of Sector 26, it is only about 0.19% in our results, which is below the size of our error bars.

From Fig. 7, it is clear that the RV measurements present some variability, the origin of which it is not clear. Looking at the long-term photometric follow-up (see Sect. 4.2) and the TESS light curves, the star seems to present a low level of activity in the form of flares that could affect the precision of the RV measurements. To be conservative, we decided to fit a simple single-planet model to the RV data and assume that the rest of the signal is attributable to systematic noise modeled by the GPs. Figure 8 presents the posterior distribution for the planetary mass. The upper-limit (95% confidence level) for the mass of TOI-2136b is 9.9 M. The planetary mass estimate from the joint fit is of Mp=4.72.6+3.1M${M_{\rm{p}}} = 4.7_{ - 2.6}^{ + 3.1}{M_ \oplus }$.

The left panel in Fig. 9 shows the median mass and radius values of TOI-2136b derived in this work compared to known transiting exoplanets with precisely determined bulk properties. The parameters for the known exoplanets were taken from TEPcat database (Southworth 2011). We only show planets having mass measurements with a precision better than 30%. In the panel, we show theoretical mass-radius relationships for planets with five different bulk compositions, including bare rocky planets with an Earth-like composition (32.5% Fe plus 67.5% MgSiO3), 100% water worlds, and Earth-like rocky cores with 0.1%, 1%, and 5% H + He envelopes. For bare rocky planets and water worlds, we take the data provided by Zeng et al. (2019) and use the intermediate value of the radii for 300 K and 500 K at a given mass. For the envelopes, we numerically integrate convective structure with Saumon et al. (1995)’s adiabats for a 75% H + 25% He mixture up to 1 kbar. Above this we use the analytic solutions from Matsui & Abe (1986) for radiative atmospheres with equilibrium temperature of 400 K (see Kurosaki & Ikoma 2017, for details). Since the mass measurement of TOI-2136b presents a large uncertainty, when compared to models, TOI-2136b could be a low-density planet (ρp ~ 1.0 g cm−3) with a significant gaseous envelope in its low mass limit, or a rocky planet with an Earth-like composition and a 1% H2 gaseous envelope in its upper mass limit.

To evaluate the potential for follow-up observations with the aim of detecting the presence of an atmosphere in TOI-2136b, we computed the transmission spectroscopy metric (TSM) of this planet following Kempton et al. (2018). This metric is an estimate of the S/N achievable with the Near Infrared Imager and Slitless Spectrograph (NIRISS) instrument on the James Webb Space Telescope (JWST, Gardner et al. 2006), using ten hours of observation time. For the median mass value of TOI-2136b, we find a TSM of ~93, and using the low and upper 1er mass limits, we find TSM values of 212 (for Mp low = 2.04 M) and 56 (for Mp up = 7.75 M). We compared the TSM obtained with the median mass value with those of previously discovered exoplanets that have a radius smaller than 4 R (see Fig. 9, right panel). We find that TOI-2136b has a relatively large TSM when compared with other planets with radii around 2 R and a similar J-band magnitude, and it also has an almost comparable TSM to hycean planets (labeled in Fig. 9, right panel), making it a very suitable target for JWST atmospheric exploration.

Table 3

Planet parameters.

thumbnail Fig. 5

TESS sector 26 and 40 folded light curves, and the best fitting model. The photometric variability present in both TESS sectors have been removed.

thumbnail Fig. 6

Ground-based transit observations of TOI-2136b. The figure shows the photometry obtained with MuSCAT (M1), MuSCAT2 (M2), MuS-CAT3 (M3), and LCO facilities, and the best fitting model (black line). The systematic effects presented in each data set were fitted and subtracted before performing the joint fit. The dates for each observation are in Table 1. The oldest observations are in the top row, and the most recent are in bottom row.

thumbnail Fig. 7

Radial velocity measurements of TOI-2136 taken with the IRD. The top panel presents the time series and best fitting model with the use of GPs to model the red noise (the fit without GPs is presented in Fig. C.2). The middle panel presents the residuals of the fit. The bottom panel shows the RV measurements in phase.

7 Transmission spectroscopy

TOI-2136b likely has a H/He gaseous envelope from the mass-radius relationship. Using the IRD high-resolution spectra taken during transit, it is possible to search for an extended atmosphere of H/He. In this particular case focusing on the helium infrared triplet at 10829.08, 10830.25, and 10830.34 Å in air.

Around the He absorption lines, there are OH emission lines and telluric H2O absorption lines that sometimes contaminate a possible planetary absorption feature (Nortmann et al. 2018; Salz et al. 2018; Alonso-Floriano et al. 2019; Palle et al. 2020). To obtain the transmission spectrum, we corrected for these telluric absorption lines, as well as for the stellar RV shift and the planet RV. To obtain an OH emission model for our data, we fitted the strongest OH emission line with a Gaussian to estimate the central wavelength, and we compared it with previous works (Palle et al. 2020; Orell-Miquel et al. 2022) to estimate the central wavelength of other OH emission lines. The region of the strongest emission was then masked.

The telluric H2O absorption lines were removed from the observed spectrum using a telluric standard star, and the correlation between line depth and airmass, following Kawauchi et al. (2018). We used a telluric standard star observed on 29 September 2020 as there were no observations of a telluric standard star on the same day. We created a telluric template spectrum by fitting the spectrum of the telluric standard star with the first spectrum of our transit time series. Figure 10 shows that, with this correction, the telluric absorption lines are eliminated from the spectra down to the noise level.

The stellar system is shifted due to the Earth's rotation and revolution. These stellar RVs in the barycentric reference frame were estimated from RA, Dec and observational time (UT) with IRAF task rvcorrect. We shifted the stellar spectra and corrected them (see the left panel of Fig. 11). However, these spectra were not aligned with the model stellar spectra at air and vacuum wavelength.

To correct the wavelength, we estimated the additional stellar RV shift from the model spectra with least-squares deconvolution (LSD) method (Donati et al. 1997). Using this method, we created a line profile of the observed spectra in velocity space, convolved with the delta functions obtained from the center wavelength and the strength of atomic and molecular lines created in the model stellar atmosphere (Watanabe et al. 2020; Collier Cameron et al. 2010). We used the observed spectra in three orders around He lines (10681.55–10979.75 Å), and the 426 absorption lines, assuming the stellar temperature and log g from VALD.

After we created a line profile of each frame, we fitted it with a flat line and a Gaussian to estimate the shifted velocity from the model. We obtained a value of −52.6 km s−1 of the shifted velocity, and we shifted the observed spectrum by that value (see right panel of Fig. 11). As a result, the observed spectra are consistent with the model spectra in air wavelength.

After the above correction, we created the template stellar spectrum to combine all out-of-transit spectra (five frames) and divided each frame by it. The planet signals were shifted in wavelength by planet orbital motion. We calculated the RV of the planet at each frame with Eq. (8) of Khalafinejad et al. (2017). The RV of TOI-2136b changes from +1.7 km s−1 to −1.5 km s−1 during transit. We shifted each frame by the calculated velocity and combined all frames during transit to create the final transmission spectrum.

Figure 12 shows the final transmission spectrum of TOI-2136 over the He i triplet spectral region. A possible He i signal at the triplet wavelengths (10829.081, 10830.25 and, 10830.34 Å) is not apparent in the data, although some spectral features fall near the expected wavelengths. In the presence of an atmospheric outflow caused by stellar wind, the He i absorption lines could be blue-shifted, but red-shifted helium absorption has been detected in HAT-P-32b (Czesla et al. 2022) and TOI 560.01 (Zhang et al. 2022). Therefore, in order to estimate the detectability of He I, we searched for the helium absorption line in a large range around the theoretical center of wavelength, and assumed that the largest absorption feature in this region of the spectra could be associated with the He triplet.

To further explore this spectral feature at ~ 10830.8 Å (which is red-shifted by ~ 15 km s−1 from the expected wavelength of the strongest triplet line), we fitted a Gaussian with the central wavelength, depth, sigma, and baseline as free parameters, using the MCMC algorithm emcee. We used a normal prior for the central wavelength between 10829.21 and 10831.38 Å, which correspond to about ±30 km s−1 from the center of the stronger He i lines, and a uniform prior for sigma between 0 and 0.344, which corresponds to the thermal broadening at a temperature of 30 000 K. We obtained a statistically nonsignificant signal with a depth of 2.2 ± 1.1%, a sigma of 0.08 ± 0.06 and an EW of 4.3 ± 2.0 mÅ (Fig. 12). However, this spectral feature is only detected at the 2 σ level, and cannot be identified as significant or planetary in origin. Thus, our data only allow us to place upper limits to the presence of He i in TOI-2136b’s atmosphere. Any possible He i signal would have an EW <7.8 mÅ, with 95% confidence, and an absorption signal <1.44%, with 95% confidence.

thumbnail Fig. 8

Planetary mass posterior distribution. We can place an upper limit on the mass of TOI-2136b of 9.9 M (95% confidence level, dashed blue line).

8 Discussion

In previous sections, we validated the planetary nature of TOI-2136.01. Figure 13 shows its location in a period–radius diagram, together with all other known planets around M-type stars whose radius measurements were better than 8% precise. In the figure, the black line indicates the radius valley measured around low-mass stars (Cloutier & Menou 2020), which is consistent with the gas-poor formation model (Lopez & Rice 2018). The dashed black line indicates the radius valley measured around Sun-like stars (Martinez et al. 2019), which is consistent with the photoevaporation (Lopez & Rice 2018) and core powered mass-loss models (Gupta & Schlichting 2019). TOI-2136b lies above both radius valleys, so it would be expected to be non-rocky in composition and possess some faction of volatile elements in its atmosphere. This prediction is consistent with our results in the mass-radius diagram (Fig. 9).

Additionally, thanks to the multiple ground-based photometric observations, we obtained the radius of TOI-2136b with an uncertainty of ~3%. Presently, there are only a few planets measured with such high precision around M-type stars. Among them, LTT3780c (2.42 ± 0.10R, 6.29 ± 0.63 M, 2.88 ± 0.28 S; Nowak et al. 2020), L321-32d (2.13 ± 0.06 R, 4.78 ± 0.43 M, 3.74 ± 0.61 S; Van Eylen et al. 2021), and TOI-776C (2.02 ± 0.14R, 5.30 ± 1.80M, 4.9 ± 0.2S; Luque et al. 2021) have similar parameters to TOI-2136b. These planets are candidates to be hycean worlds, which are composed of water-rich interiors with massive oceans underlying H2-rich atmospheres (Madhusudhan et al. 2021). Hycean planets have received attention since a possible biosignature might be detectable in their atmospheres with a relatively modest amount of JWST observation time.

We also obtained a transmission spectrum of the planet around the He i 10 830 Å absorption lines. In other works, helium absorption is detected on the atmosphere of three sub-Neptunes (TOI-560.01; Zhang et al. 2022, GJ3470b; Palie et al. 2020, and GJ1214b; Orell-Miquel et al. 2022). To compare our results with previous He i detections in the literature, we calculated the equivalent height of the atmosphere normalized by the scale height. This value for TOI-2136b is < 114, adopting a mass of Mp = 4.7 M and is comparable with the value of TOI 560.01 (Fig. 4 in Orell-Miquel et al. 2022). TOI 560.01 has a detection of He i (Zhang et al. 2022), meaning that TOI-2136b likely has a less extended or a less rich H/He atmosphere than TOI 560.01. For GJ1214b and GJ3470b, the ratio between the equivalent height and atmospheric scale height is ~57 for both planets, about half of our upper limit for TOI-2136b. Despite our non-detection of He i on this planet, the atmospheric metrics and its place in the mass-radius space suggest that TOI-2136b is an ideal candidate to search for atmospheric signals using ground-based high-resolution spectroscopic observations taken during a transit event and space-based observations with JWST.

During the submission and revision stage of this manuscript, two studies on the TOI-2136 system were announced: Gan et al. (2022); Beard et al. (2022). Gan et al. (2022) used TESS data, ground-based follow-up observations, and the Spirou instrument on the Canada-France-Hawaii Telescope (CFHT) to obtain RV measurements; the study presents a planetary radius of 2.19 ± 0.17 R and a mass of 6.4 ± 2.4 M for the transiting planet. Beard et al. (2022) used TESS data, ground-based follow-up observations, and the Hobby-Eberly Telescope (HET) instrument Habitable-Zone Planet Finder (HPF) instrument to obtain RV measurements. Beard et al. (2022) presents a planetary radius of 2.09 ± 0.08 R and a mass upper limit of < 15 M, with a median planet mass of 4.64 M. Our results for the planetary mass and radius of TOI-2136b are consistent with both studies (taking into account the uncertainties).

thumbnail Fig. 9

Mass-radius diagram and TSM for TOI-2136b. Left panel: mass and radius values for known exoplanets (data taken from TEPcat database. Southworth 2011). The position of TOI-2136b is marked by the red star. In the figure we only show planets having mass measurements with a precision better than 30%. The lines represent theoretical mass-radius relationships for five different bulk compositions, including bare rocky planets with an Earth-like composition, pure water worlds, and Earth-like rocky cores with 0.1%, 1%, and 5% H + He envelopes with radiative equilibrium temperature of 400 K (see text). The orange points represent planets orbiting M-type stars, and the gray squares are planets around other types of stars. Right panel: apparent magnitude in J-band vs. the TSM by Kempton et al. (2018) for TOI-2136b and known transiting planets with radii smaller than 4R. The other labeled planets are the hycean candidates.

thumbnail Fig. 10

Example of a single IRD spectrum of TOI-2136 around the helium triplet before (red) and after (black) the telluric correction. Blue and orange lines indicate the telluric absorption lines and OH emission lines. This figure indicates the position of the emission lines and the correction of the telluric absorption lines down to the noise level.

thumbnail Fig. 11

Example of a single IRD spectrum of TOI-2136b (black) compared with the BT-Settl model (Allard 2014) spectrum in vacuum (blue) and air (red). The left figure indicates the raw IRD spectrum and the right figure indicates the IRD-spectrum corrected shifted velocity. Vertical green dash lines in the left and right panels indicate the center of the predicted He triplet in air.

thumbnail Fig. 12

Transmission spectrum of TOI-2136 around the He triplet. The blank region is where we mask the strongest OH emission. The red line is contaminated by the other shallower OH emission lines. The orange line is the best-fitted model with a Gaussian function. Vertical green dashed lines indicate the center of the predicted He triplet in air, and the horizontal blue lines are standard deviations between 10820 and 10840 Å.

9 Summary

In this work, we reported the follow-up efforts on the planet candidate TOI-2136b, a sub-Neptune planet orbiting around a nearby M3 dwarf. The initial observations made by TESS led to the detection of the transiting planet candidate. Multicolor, ground-based, follow-up photometric observations and RV measurements made it possible to confirm the planetary nature of the candidate and rule out false positives.

Using the IRD instrument on the Subaru 8.2 m telescope, we obtained high-resolution spectra and derived stellar parameters for the planet host star. We find that TOI-2136 is a nearby (d = 33.361 ± 0.019 pc) M dwarf with Teff = 3373 ± 108 K, [Fe/H] = 0.02 ± 0.14 dex, and an estimated stellar mass of M* = 0.3272 ± 0.0082 M and a radius of R* = 0.3440 ± 0.0099 R. We also measured the rotation period of the star using MEarth and ZTF long-term photometric observations, and found a value of Prot = 82.56 ± 0.45 days.

We analyzed the available transit and RV observations of the planet candidate fitting the data jointly. The transit observations and RV measurements from IRD were fitted simultaneously using an MCMC procedure that included Gaussian processes to model the systematic effects present in the TESS and RV measurements. We find that TOI-2136b has a planetary radius of Rp = 2.20 ± 0.07 R and a mass of Mp=4.72.6+3.1M${M_{\rm{p}}} = 4.7_{ - 2.6}^{ + 3.1}{M_ \oplus }$ (with an upper limit on the mass value of < 9.9 M, with 95% confidence). Our model supports a circular orbit (e=0.070.05+0.09$e = 0.07_{ - 0.05}^{ + 0.09}$) with an orbital period of 7.851925 ± 0.000016 days.

Finally, we observed a transit event using high-resolution spectroscopy in an effort to detect the existence of an extended atmosphere around this planet. In particular we focused on the He i 10 830 Å absorption lines. We found no statistical significant extra absorption around the He i lines that could be attributable to the planet, but we can place an upper limit on the presence of a He-rich extended atmosphere: an absorption depth of < 1.44% and an EW of <7.8 mÅ, with 95% confidence.

TOI-2136b is a small sub-Neptune planet orbiting a relatively bright (J = 10.1 mag) M star. Based on the TSM, this planet presents a very valuable opportunity for further atmospheric characterization studies and the search for additional planets in the system.

thumbnail Fig. 13

Period-radius diagram for known exoplanets around M-type stars with radius measurements with a precision better than 8% (orange points). The red star represents TOI-2136b. The black line and the dashed line indicate the radius valley around low-mass stars (Cloutier & Menou 2020) and Sun-like stars (Martinez et al. 2019), respectively.

Acknowledgements

This research is based on data collected at Subaru Telescope, which is operated by the National Astronomical Observatory of Japan. We are honored and grateful for the opportunity of observing the Universe from Maunakea, which has the cultural, historical and natural significance in Hawaii. The part of our data analysis was carried out on common use data analysis computer system at the Astronomy Data Center, ADC, of the National Astronomical Observatory of Japan. This paper makes use of data from the MEarth Project, which is a collaboration between Harvard University and the Smithsonian Astro-physical Observatory. The MEarth Project acknowledges funding from the David and Lucile Packard Fellowship for Science and Engineering, the National Science Foundation under grants AST-0807690, AST-1109468, AST-1616624 and AST-1004488 (Alan T. Waterman Award), the National Aeronautics and Space Administration under Grant No. 80NSSC18K0476 issued through the XRP Program, and the John Templeton Foundation. Based on observations obtained with the Samuel Oschin 48-inch Telescope at the Palomar Observatory as part of the Zwicky Transient Facility project. Z.T.F. is supported by the National Science Foundation under Grant No. AST-1440341 and a collaboration including Caltech, IPAC, the Weizmann Institute for Science, the Oskar Klein Center at Stockholm University, the University of Maryland, the University of Washington, Deutsches Elektronen-Synchrotron and Humboldt University, Los Alamos National Laboratories, the TANGO Consortium of Taiwan, the University of Wisconsin at Milwaukee, and Lawrence Berkeley National Laboratories. Operations are conducted by COO, IPAC, and UW. Our data reductions benefited from PyRAF and PyFITS that are the products of the Space Telescope Science Institute, which is operated by AURA for NASA. G.M. has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No. 895525. R.L. acknowledges financial support from the Spanish Ministerio de Ciencia e Innovación, through project PID2019-109522GB-C52, and the Centre of Excellence “Severo Ochoa” award to the Instituto de Astrofísica de Andalucía (SEV-2017-0709). J.K. gratefully acknowledge the support of the Swedish National Space Agency (SNSA; DNR 2020-00104). This work is partly supported by JSPS KAKENHI Grant Numbers JP21K20376, JP21K13975, JP21H00035, JP21K20388, JP20K14518, JP20K14521, JP19K14783, JP18H05439, JP17H04574, by Grant-in-Aid for JSPS Fellows Grant Number JP20J21872, by JST CREST Grant Number JPMJCR1761, by Astrobiology Center PROJECT Research AB031014, by Astrobiology Center SATELLITE Research project AB022006, and the Astrobiology Center of National Institutes of Natural Sciences (NINS) (Grant Number AB031010). M.T. is supported by MEXT/JSPS KAKENHI grant Nos. 18H05442, 15H02063, and 22000005. Some of the observations in the paper made use of the High-Resolution Imaging instrument Alopeke obtained under Gemini LLP Proposal Number: GN/S-2021A-LP-105. Alopeke was funded by the NASA Exoplanet Exploration Program and built at the NASA Ames Research Center by Steve B. Howell, Nic Scott, Elliott P. Horch, and Emmett Quigley. Alopeke was mounted on the Gemini North (and/or South) telescope of the international Gemini Observatory, a program of NSF’s OIR Lab, which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation, on behalf of the Gemini partnership: the National Science Foundation (United States), National Research Council (Canada), Agencia Nacional de Investigación y Desarrollo (Chile), Mnisterio de Ciencia, Tecnología e Innovación (Argentina), Ministério da Ciência, Tecnologia, Inovações e Comunicações (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea).

Appendix A Radial velocities

Table A.1

Subaru IRD RV measurements

Appendix B Stellar rotation priors

Table B.1

Stellar rotation period fitted parameters, prior functions, and limits.

Appendix C Light curve and radial velocity joint fit

Table C.1

Global fit parameters, prior functions, and limits.

thumbnail Fig. C.1

Correlation plot for the fitted orbital parameters. Limb-darkening coefficients and systematic effect parameters were intentionally left out for easy viewing. The blue lines give the median values of the parameters.

thumbnail Fig. C.2

RV measurements of TOI-2136 taken with IRD. The top panel presents the time series and best fitting model without GPs. The middle panel presents the residuals of the fit. The bottom panel shows the RV measurements in phase after subtracting the red noise.

thumbnail Fig. C.3

Correlation plot for the fitted orbital parameters without GPs. Limb-darkening coefficients and systematic effect parameters were intentionally left out for easy viewing. The blue lines give the median values of the parameters.

References

  1. Akaike, H. 1974, IEEE Trans. Automatic Control, 19, 716 [CrossRef] [Google Scholar]
  2. Allard, F. 2014, in Exploring the Formation and Evolution of Planetary Systems, eds. M. Booth, B. C. Matthews, & J. R. Graham, 299, 271 [Google Scholar]
  3. Allart, R., Bourrier, V., Lovis, C., et al. 2018, Science, 362, 1384 [Google Scholar]
  4. Aller, A., Lillo-Box, J., Jones, D., Miranda, L. F., & Barceló Forteza, S. 2020, A&A, 635, A128 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Alonso-Floriano, F. J., Snellen, I. A. G., Czesla, S., et al. 2019, A&A, 629, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Ambikasaran, S., Foreman-Mackey, D., Greengard, L., Hogg, D. W., & O’Neil, M. 2015, IEEE Trans. Pattern Anal. Mach. Intell., 38, 252 [Google Scholar]
  7. Beard, C., Robertson, P., Kanodia, S., et al. 2022, AJ, 163, 286 [NASA ADS] [CrossRef] [Google Scholar]
  8. Berta, Z. K., Irwin, J., Charbonneau, D., Burke, C. J., & Falco, E. E. 2012, AJ, 144, 145 [Google Scholar]
  9. Casagrande, L., Schönrich, R., Asplund, M., et al. 2011, A&A, 530, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Ciardi, D. R., Beichman, C. A., Horch, E. P., & Howell, S. B. 2015, ApJ, 805, 16 [NASA ADS] [CrossRef] [Google Scholar]
  11. Cloutier, R., & Menou, K. 2020, AJ, 159, 211 [NASA ADS] [CrossRef] [Google Scholar]
  12. Collier Cameron, A., Guenther, E., Smalley, B., et al. 2010, MNRAS, 407, 507 [NASA ADS] [CrossRef] [Google Scholar]
  13. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
  14. Czesla, S., Lampón, M., Sanz-Forcada, J., et al. 2022, A&A, 657, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Donati, J. F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, MNRAS, 291, 658 [Google Scholar]
  16. Elkins-Tanton, L. T., & Seager, S. 2008, ApJ, 685, 1237 [NASA ADS] [CrossRef] [Google Scholar]
  17. Engle, S. G., & Guinan, E. F. 2018, RNAAS, 2, 34 [NASA ADS] [Google Scholar]
  18. Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
  19. Foreman-Mackey, D., Agol, E., Ambikasaran, S., & Angus, R. 2017, AJ, 154, 220 [Google Scholar]
  20. Fukui, A., Narita, N., Tristram, P. J., et al. 2011, PASJ, 63, 287 [NASA ADS] [Google Scholar]
  21. Fulton, B. J., Petigura, E. A., Howard, A. W., et al. 2017, AJ, 154, 109 [Google Scholar]
  22. Fulton, B. J., Petigura, E. A., Blunt, S., & Sinukoff, E. 2018, PASP, 130, 044504 [Google Scholar]
  23. Gan, T., Soubkiou, A., Wang, S. X., et al. 2022, MNRAS, 514, 4120 [CrossRef] [Google Scholar]
  24. Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, Space Sci. Rev., 123, 485 [Google Scholar]
  25. Gibson, N. P., Aigrain, S., Roberts, S., et al. 2012, MNRAS, 419, 2683 [NASA ADS] [CrossRef] [Google Scholar]
  26. Ginzburg, S., Schlichting, H. E., & Sari, R. 2016, ApJ, 825, 29 [Google Scholar]
  27. Ginzburg, S., Schlichting, H. E., & Sari, R. 2018, MNRAS, 476, 759 [Google Scholar]
  28. Gupta, A., & Schlichting, H. E. 2019, MNRAS, 487, 24 [Google Scholar]
  29. Gupta, A., & Schlichting, H. E. 2021, MNRAS, 504, 4634 [NASA ADS] [CrossRef] [Google Scholar]
  30. Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Hippke, M., & Heller, R. 2019, A&A, 623, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Hirano, T., Kuzuhara, M., Kotani, T., et al. 2020, PASJ, 72, 93 [Google Scholar]
  33. Hirano, T., Livingston, J. H., Fukui, A., et al. 2021, AJ, 162, 161 [NASA ADS] [CrossRef] [Google Scholar]
  34. Howell, S. B., Everett, M. E., Sherry, W., Horch, E., & Ciardi, D. R. 2011, AJ, 142, 19 [Google Scholar]
  35. Howell, S. B., Everett, M. E., Horch, E. P., et al. 2016, ApJ, 829, L2 [NASA ADS] [CrossRef] [Google Scholar]
  36. Howell, S. B., Matson, R. A., Ciardi, D. R., et al. 2021, VizieR Online Data Catalog: J/AJ/161/164 [Google Scholar]
  37. Ishikawa, H. T., Aoki, W., Kotani, T., et al. 2020, PASJ, 72, 102 [CrossRef] [Google Scholar]
  38. Ishikawa, H. T., Aoki, W., Hirano, T., et al. 2022, AJ, 163, 72 [NASA ADS] [CrossRef] [Google Scholar]
  39. Jenkins, J. M. 2002, ApJ, 575, 493 [Google Scholar]
  40. Jenkins, J. M., Twicken, J. D., McCauliff, S., et al. 2016, SPIE Conf. Ser., 9913, 99133E [Google Scholar]
  41. Jenkins, J. M., Tenenbaum, P., Seader, S., et al. 2020, Kepler Data Processing Handbook: Transiting Planet Search, Kepler Science Document KSCI-19081-003 [Google Scholar]
  42. Jin, S., Mordasini, C., Parmentier, V., et al. 2014, ApJ, 795, 65 [Google Scholar]
  43. Johnson, D. R. H., & Soderblom, D. R. 1987, AJ, 93, 864 [Google Scholar]
  44. Kawauchi, K., Narita, N., Sato, B., et al. 2018, PASJ, 70, 84 [NASA ADS] [CrossRef] [Google Scholar]
  45. Kempton, E. M. R., Bean, J. L., Louie, D. R., et al. 2018, PASP, 130, 114401 [Google Scholar]
  46. Khalafinejad, S., von Essen, C., Hoeijmakers, H. J., et al. 2017, A&A, 598, A131 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Kipping, D. M. 2013, MNRAS, 435, 2152 [Google Scholar]
  48. Kotani, T., Tamura, M., Nishikawa, J., et al. 2018, SPIE Conf. Ser., 10702, 1070211 [Google Scholar]
  49. Krishnamurthy, V., Hirano, T., Stefánsson, G., et al. 2021, AJ, 162, 82 [CrossRef] [Google Scholar]
  50. Kupka, F., Piskunov, N., Ryabchikova, T. A., Stempels, H. C., & Weiss, W. W. 1999, A&AS, 138, 119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Kurosaki, K., & Ikoma, M. 2017, AJ, 153, 260 [NASA ADS] [CrossRef] [Google Scholar]
  52. Kuzuhara, M., Hirano, T., Kotani, T., et al. 2018, SPIE Conf. Ser., 10702, 1070260 [NASA ADS] [Google Scholar]
  53. Lépine, S., & Shara, M. M. 2005, AJ, 129, 1483 [Google Scholar]
  54. Lester, K. V., Matson, R. A., Howell, S. B., et al. 2021, AJ, 162, 75 [NASA ADS] [CrossRef] [Google Scholar]
  55. Lopez, E. D., & Rice, K. 2018, MNRAS, 479, 5303 [NASA ADS] [CrossRef] [Google Scholar]
  56. Luque, R., Serrano, L. M., Molaverdikhani, K., et al. 2021, A&A, 645, A41 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Madhusudhan, N., Piette, A. A. A., & Constantinou, S. 2021, ApJ, 918, 1 [NASA ADS] [CrossRef] [Google Scholar]
  58. Mandel, K., & Agol, E. 2002, ApJ, 580, L171 [Google Scholar]
  59. Mann, A. W., Feiden, G. A., Gaidos, E., Boyajian, T., & von Braun, K. 2015, ApJ, 804, 64 [Google Scholar]
  60. Mann, A. W., Dupuy, T., Kraus, A. L., et al. 2019, ApJ, 871, 63 [Google Scholar]
  61. Martinez, C. F., Cunha, K., Ghezzi, L., & Smith, V. V. 2019, ApJ, 875, 29 [Google Scholar]
  62. Masci, F. J., Laher, R. R., Rusholme, B., et al. 2019, PASP, 131, 018003 [Google Scholar]
  63. Matson, R. A., Howell, S. B., Horch, E. P., & Everett, M. E. 2018, AJ, 156, 31 [NASA ADS] [CrossRef] [Google Scholar]
  64. Matsui, T., & Abe, Y. 1986, Nature, 322, 526 [CrossRef] [Google Scholar]
  65. McCully, C., Volgenau, N. H., Harbeck, D.-R., et al. 2018, SPIE Conf. Ser., 10707, 107070K [Google Scholar]
  66. Mordasini, C. 2020, A&A, 638, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. Morris, R. L., Twicken, J. D., Smith, J. C., et al. 2020, Kepler Data Processing Handbook: Photometric Analysis, Kepler Science Document KSCI-19081-003 [Google Scholar]
  68. Murgas, F., Astudillo-Defru, N., Bonfils, X., et al. 2021, A&A, 653, A60 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Narita, N., Fukui, A., Kusakabe, N., et al. 2015, J. Astron. Telescopes Instrum. Syst., 1, 045001 [CrossRef] [Google Scholar]
  70. Narita, N., Fukui, A., Kusakabe, N., et al. 2019, J. Astron. Telescopes Instrum. Syst., 5, 015001 [Google Scholar]
  71. Narita, N., Fukui, A., Yamamuro, T., et al. 2020, SPIE Conf. Ser., 11447, 114475K [NASA ADS] [Google Scholar]
  72. Newton, E. R., Irwin, J., Charbonneau, D., et al. 2016, ApJ, 821, 93 [Google Scholar]
  73. Nortmann, L., Pallé, E., Salz, M., et al. 2018, Science, 362, 1388 [Google Scholar]
  74. Nowak, G., Luque, R., Parviainen, H., et al. 2020, A&A, 642, A173 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. Nutzman, P., & Charbonneau, D. 2008, PASP, 120, 317 [Google Scholar]
  76. Orell-Miquel, J., Murgas, F., Pallé, E., et al. 2022, A&A, 659, A55 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  77. Owen, J. E., & Wu, Y. 2013, ApJ, 775, 105 [Google Scholar]
  78. Owen, J. E., & Wu, Y. 2017, ApJ, 847, 29 [Google Scholar]
  79. Palle, E., Nortmann, L., Casasayas-Barris, N., et al. 2020, A&A, 638, A61 [EDP Sciences] [Google Scholar]
  80. Parviainen, H. 2015, MNRAS, 450, 3233 [Google Scholar]
  81. Parviainen, H., & Aigrain, S. 2015, MNRAS, 453, 3821 [Google Scholar]
  82. Petigura, E. A., Rogers, J. G., Isaacson, H., et al. 2022, AJ, 163, 179 [NASA ADS] [CrossRef] [Google Scholar]
  83. Raftery, A. E. 1995, Sociol. Methodol., 25, 111 [CrossRef] [Google Scholar]
  84. Rasmussen, C., & Williams, C. 2010, The MIT Press, 122, 935 [Google Scholar]
  85. Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telescopes Instrum. Syst., 1, 014003 [Google Scholar]
  86. Riello, M., De Angeli, F., Evans, D. W., et al. 2021, A&A, 649, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Ryabchikova, T., Piskunov, N., Kurucz, R. L., et al. 2015, Phys. Scr, 90, 054005 [Google Scholar]
  88. Salz, M., Czesla, S., Schneider, P. C., et al. 2018, A&A, 620, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Saumon, D., Chabrier, G., & van Horn, H. M. 1995, ApJS, 99, 713 [NASA ADS] [CrossRef] [Google Scholar]
  90. Schwarz, G. 1978, Ann. Stat., 6, 461 [Google Scholar]
  91. Scott, N. J., Howell, S. B., Gnilka, C. L., et al. 2021, Front. Astron. Space Sci., 8, 138 [NASA ADS] [CrossRef] [Google Scholar]
  92. Smith, J. C., Stumpe, M. C., Van Cleve, J. E., et al. 2012, PASP, 124, 1000 [Google Scholar]
  93. Southworth, J. 2011, MNRAS, 417, 2166 [Google Scholar]
  94. Spake, J. J., Sing, D. K., Evans, T. M., et al. 2018, Nature, 557, 68 [Google Scholar]
  95. Stassun, K. G., Oelkers, R. J., Paegert, M., et al. 2019, AJ, 158, 138 [Google Scholar]
  96. Stumpe, M. C., Smith, J. C., Catanzarite, J. H., et al. 2014, PASP, 126, 100 [Google Scholar]
  97. Tamura, M., Suto, H., Nishikawa, J., et al. 2012, SPIE Conf. Ser., 8446, 84461T [NASA ADS] [Google Scholar]
  98. Tody, D. 1993, in Astronomical Society of the Pacific Conference Series, Astronomical Data Analysis Software and Systems II, eds. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes, 52, 173 [NASA ADS] [Google Scholar]
  99. Tsuji, T. 1978, A&A, 62, 29 [NASA ADS] [Google Scholar]
  100. Twicken, J. D., Catanzarite, J. H., Clarke, B. D., et al. 2018, PASP, 130, 064502 [Google Scholar]
  101. Van Eylen, V., Agentoft, C., Lundkvist, M. S., et al. 2018, MNRAS, 479, 4786 [Google Scholar]
  102. Van Eylen, V., Astudillo-Defru, N., Bonfils, X., et al. 2021, MNRAS, 507, 2154 [NASA ADS] [Google Scholar]
  103. Watanabe, N., Narita, N. & Johnson, M. C. 2020, PASJ, 72, 19 [CrossRef] [Google Scholar]
  104. Wu, Y. 2019, ApJ, 874, 91 [NASA ADS] [CrossRef] [Google Scholar]
  105. Zechmeister, M. & Kürster, M. 2009, A&A, 496, 577 [CrossRef] [EDP Sciences] [Google Scholar]
  106. Zeng, L., Jacobsen, S. B., Sasselov, D. D., et al. 2019, Proc. Natl. Acad. Sci. U.S.A., 116, 9723 [Google Scholar]
  107. Zhang, M., Knutson, H. A., Wang, L., Dai, F. & Barragán, O. 2022, AJ, 163, 67 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Ground-based photometric observations.

Table 2

Stellar parameters.

Table 3

Planet parameters.

Table A.1

Subaru IRD RV measurements

Table B.1

Stellar rotation period fitted parameters, prior functions, and limits.

Table C.1

Global fit parameters, prior functions, and limits.

All Figures

thumbnail Fig. 1

TESS target pixel file image of TOI-2136 observed in Sector 26 (left) and Sector 40 (right). The pixels in red show the photometric aperture used by TESS to create the PDCSAP light curves. The position of nearby stars and their magnitudes are represented by the red circles. This image was produced using tpfplotter (Aller et al. 2020).

In the text
thumbnail Fig. 2

Gemini/‘Alopeke 5-sigma contrast curves in 562 nm (blue) and 832 nm (red). The inset figure is the 832 nm reconstructed speckle image.

In the text
thumbnail Fig. 3

Ground-based long-term photometric observations of TOI-2136 using MEarth and ZTF, and their respective GLS periodograms (Zechmeister & Kürster 2009) for each data set. The stellar photometric variability is detected with a stronger signal in MEarth data. Fitting a kernel with a periodic term to all the data sets using Gaussian processes, we find a stellar rotation period of 82.56 ± 0.45 days.

In the text
thumbnail Fig. 4

TOI-2136 light curves of TESS Sector 26 (top panel) and Sector 40 (bottom panel). The black points represent TESS binned photometry. The gray points are the TESS PDCSAP data. The best fitting transit model model including systematic effects is shown in dark blue. Light-blue points are the TESS data after correcting the systematic effects. The best fitting transit model is shown in red. Transit events of TOI-2136b are marked with black vertical lines.

In the text
thumbnail Fig. 5

TESS sector 26 and 40 folded light curves, and the best fitting model. The photometric variability present in both TESS sectors have been removed.

In the text
thumbnail Fig. 6

Ground-based transit observations of TOI-2136b. The figure shows the photometry obtained with MuSCAT (M1), MuSCAT2 (M2), MuS-CAT3 (M3), and LCO facilities, and the best fitting model (black line). The systematic effects presented in each data set were fitted and subtracted before performing the joint fit. The dates for each observation are in Table 1. The oldest observations are in the top row, and the most recent are in bottom row.

In the text
thumbnail Fig. 7

Radial velocity measurements of TOI-2136 taken with the IRD. The top panel presents the time series and best fitting model with the use of GPs to model the red noise (the fit without GPs is presented in Fig. C.2). The middle panel presents the residuals of the fit. The bottom panel shows the RV measurements in phase.

In the text
thumbnail Fig. 8

Planetary mass posterior distribution. We can place an upper limit on the mass of TOI-2136b of 9.9 M (95% confidence level, dashed blue line).

In the text
thumbnail Fig. 9

Mass-radius diagram and TSM for TOI-2136b. Left panel: mass and radius values for known exoplanets (data taken from TEPcat database. Southworth 2011). The position of TOI-2136b is marked by the red star. In the figure we only show planets having mass measurements with a precision better than 30%. The lines represent theoretical mass-radius relationships for five different bulk compositions, including bare rocky planets with an Earth-like composition, pure water worlds, and Earth-like rocky cores with 0.1%, 1%, and 5% H + He envelopes with radiative equilibrium temperature of 400 K (see text). The orange points represent planets orbiting M-type stars, and the gray squares are planets around other types of stars. Right panel: apparent magnitude in J-band vs. the TSM by Kempton et al. (2018) for TOI-2136b and known transiting planets with radii smaller than 4R. The other labeled planets are the hycean candidates.

In the text
thumbnail Fig. 10

Example of a single IRD spectrum of TOI-2136 around the helium triplet before (red) and after (black) the telluric correction. Blue and orange lines indicate the telluric absorption lines and OH emission lines. This figure indicates the position of the emission lines and the correction of the telluric absorption lines down to the noise level.

In the text
thumbnail Fig. 11

Example of a single IRD spectrum of TOI-2136b (black) compared with the BT-Settl model (Allard 2014) spectrum in vacuum (blue) and air (red). The left figure indicates the raw IRD spectrum and the right figure indicates the IRD-spectrum corrected shifted velocity. Vertical green dash lines in the left and right panels indicate the center of the predicted He triplet in air.

In the text
thumbnail Fig. 12

Transmission spectrum of TOI-2136 around the He triplet. The blank region is where we mask the strongest OH emission. The red line is contaminated by the other shallower OH emission lines. The orange line is the best-fitted model with a Gaussian function. Vertical green dashed lines indicate the center of the predicted He triplet in air, and the horizontal blue lines are standard deviations between 10820 and 10840 Å.

In the text
thumbnail Fig. 13

Period-radius diagram for known exoplanets around M-type stars with radius measurements with a precision better than 8% (orange points). The red star represents TOI-2136b. The black line and the dashed line indicate the radius valley around low-mass stars (Cloutier & Menou 2020) and Sun-like stars (Martinez et al. 2019), respectively.

In the text
thumbnail Fig. C.1

Correlation plot for the fitted orbital parameters. Limb-darkening coefficients and systematic effect parameters were intentionally left out for easy viewing. The blue lines give the median values of the parameters.

In the text
thumbnail Fig. C.2

RV measurements of TOI-2136 taken with IRD. The top panel presents the time series and best fitting model without GPs. The middle panel presents the residuals of the fit. The bottom panel shows the RV measurements in phase after subtracting the red noise.

In the text
thumbnail Fig. C.3

Correlation plot for the fitted orbital parameters without GPs. Limb-darkening coefficients and systematic effect parameters were intentionally left out for easy viewing. The blue lines give the median values of the parameters.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.