Open Access
Issue
A&A
Volume 691, November 2024
Article Number A12
Number of page(s) 14
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202450934
Published online 25 October 2024

© The Authors 2024

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

The sub-Saturn classification of exoplanets refers to planets larger than Neptune but smaller than Saturn, typically falling within the range of 4–8 R (Petigura et al. 2016). Also known as Super-Neptunes, they are roughly ten times rarer than planets between the sizes of Earth and Neptune (Brady et al. 2018). Sub-Saturns are often described as failed gas giants, possessing equally massive cores but significantly smaller total masses due to their much smaller accreted envelopes (Millholland et al. 2020). These planets serve as valuable laboratories for studying the physics of envelope accretion, as Neptunes usually have envelope mass fractions (fenv) of less than 10%, while gas giants are almost entirely composed of envelopes (fenv ~ 100%) (Petigura et al. 2017). The absence of these planets in our Solar System highlights the variety of planetary systems possible, and studying them around other stars offers useful insights into this diversity (Petigura et al. 2017; Fairnington et al. 2024). Sub-Saturn planets show mass diversity regardless of size, spanning the mass range of 10–100 M, strongly correlated to the host star metallicity (Petigura et al. 2017).

The mass distribution between the core and envelope of sub-Saturns provides clues regarding their formation and evolution. A few theories have been proposed for the formation of sub-Saturns. One of them, proposed by Petigura et al. (2017), suggests the failure of runaway accretion, where the surrounding disk around a young star simply lacks sufficient material for the planet to become a full gas giant like Jupiter or Saturn. Alternatively, Lee & Chiang (2016) and Lee et al. (2018) proposed formation in a gas-depleted disk. In this scenario, the disk itself is depleted of gas, possibly due to stellar radiation or winds, leaving a sub-Saturn with limited building materials to reach the size of Jupiter or Saturn. Another proposed mechanism is formation through atmospheric mass loss, as suggested by Hallatt & Lee (2022). This theory proposes that sub-Saturns initially form as giant planets but lose their atmosphere over time, causing them to shrink down (Owen & Lai 2018). Two main mechanisms are proposed for atmospheric mass loss: photoevaporation (e.g., Owen & Lai 2018; Owen 2019; Owen & Jackson 2012; Jin et al. 2014; Lopez 2017) and high-eccentricity orbital migration followed by tidal interactions with the star, also known as tidal stripping (e.g., Matsakos & Königl 2016; Owen & Lai 2018). The formation pathways followed by a planet in its life among these two mechanisms result in diverse densities and masses.

So far1, only 74 sub-Saturns are known with mass and radius constraints better than 50 and 20%, respectively, and 41 of them have orbital periods of less than 10 days (hot sub-Saturns). The noticeable scarcity of sub-Jovian planets on short orbits (P < 10 days) is known as the “sub-Jovian desert” or “hot Neptunian desert” and that has been extensively discussed in the literature (e.g., Lecavelier Des Etangs 2007; Davis & Wheatley 2009; Szabó & Kiss 2011; Beaugé & Nesvorný 2013; Mazeh et al. 2016). The dearth is not due to observational biases, as close-in super-Neptunes or sub-Saturns are readily detectable through both transit and Radial Velocity (RV) methods. The origin of the Neptunian desert is also associated with the formation and evolution mechanism. For example, the upper edge of the Neptunian desert (see Section 4.3) is linked with the tidal heating mechanism as it is argued that denser and more massive sub-Saturns tend to lose their atmospheres primarily through tidal heating rather than photoevaporation. Photoevaporation remains the dominant mechanism for atmospheric loss in less massive and less dense sub-Saturns (Vissapragada et al. 2022; Petigura et al. 2017). Thorngren et al. (2023) offer a different viewpoint suggesting that photoevaporation might not be entirely negligible even for massive sub-Saturns, implying a possible interplay between photoevaporation and tidal heating. Alternatively, a deeper understanding can also be gained by detailed estimation of the planet’s internal structure and composition. This approach is particularly valuable because the formation pathway of these planets is closely linked to their core mass and the fenv (Fortney et al. 2007; Petigura et al. 2016). Despite our current knowledge, gaps regarding the formation and evolution of these planets and their absence in our own solar system make them even more intriguing planets for further study.

Here, we present the discovery and characterization of TOI-6651b, a dense hot sub-Saturn transiting a metal-rich G-type star, utilizing data obtained from PARAS-2 and TESS. PARAS-2 is a high-resolution spectrograph (R ≈ 107 000) attached to the recently installed Physical Research Laboratory (PRL) 2.5 m telescope at Mount Abu Observatory, located at Gurushikhar, Mount Abu, India (Chakraborty et al. 2018, 2024). In addition to TESS and PARAS-2, we have utilized speckle imaging data in this work, which is explained in Section 2. In Section 3, our focus is directed toward the host star and the comprehensive characterization of the star-planet system, achieved through the combined analysis of transit photometry and RV data. Section 4 is dedicated to exploring the implications of our findings, while Section 5 provides a summary of the paper.

2 Observations

2.1 TESS photometry

TOI-6651 (TIC 174328045) was observed with TESS at a 30-minute cadence in sector 17 from October 8 to November 2, 2019, and at a 200-second cadence in sector 57 from September 30 to October 29, 2022. Additionally, TOI-6651 is slated for observation in sector 84 from October 1 to October 26, 2024, as indicated by the TESS-point Web Tool2. TOI-6651 was initially identified as a community TESS object of interest (CTOI; Guerrero et al. 2021) by Montalto (2022) using the DIAmante3 pipeline, which detected transit signals in light curve (LC) derived from TESS Full Frame Images (FFIs) of sector 17. Its status was subsequently elevated to TOI in September 2023 after LC generated from FFIs of sector 57 were scrutinized for transit signals by the Science Processing Operations Center (SPOC; Jenkins et al. 2016) and Quick Look Pipeline (QLP; Huang et al. 2020a,b).

For photometric analysis, we used the Pre-search Data Conditioning Simple Aperture Photometry (PDCSAP) flux time series processed by the SPOC pipeline (Jenkins et al. 2016), publicly accessible through the Mikulski Archive for Space Telescopes (MAST4). For detrending and normalizing the LCs, we employed a Python package citlalicue (Barragán et al. 2021), which involves Gaussian Process regression to remove systematics and stellar variability. The resulting LCs, post-procedure, are illustrated in Figure 1. We detected a total of 9 transits: 4 in sector 17 and 5 in sector 57. However, due to telescope reorientation during data transmission, one transit in sector 17 was omitted.

We also use tpfplotter (Aller et al. 2020) to overplot the Gaia DR3 (Gaia Collaboration 2023) sources onto the TESS target pixel file (TPF). This enabled us to identify potential sources of dilution in the TESS photometry, with a magnitude difference (Δm) limit of up to 6. The tpfplotter generated images are shown in Figure 2. We found one additional source at the edge of the aperture mask of the target in sector 17. However, given its approximately six-magnitude difference compared to TOI-6651, this source is too faint to yield any significant dilution. No additional sources contributing to dilution were found in the observations of sector 57. Additionally, the pipeline analysis revealed that 99.63% of the light in the optimal aperture originates from the target itself rather than other stellar sources, indicating negligible dilution from faint background stars.

2.2 Speckle observations

Labeyrie (1970) introduced the concept of speckle interferometric imaging as a technique to minimize the effects of atmospheric turbulence that primarily limits the resolution of ground-based telescopes. This method involves acquiring multiple short-exposure images during which the atmosphere is partially frozen, which are then processed to reconstruct a single diffraction-limited image. Speckle imaging serves as a valuable follow-up observation to exoplanet detection studies by providing such high angular resolution images as it aids in disentangling signals from close companions to the host star, thus reducing the likelihood of false detections or incorrect estimation of stellar and planetary parameters (e.g., Ciardi et al. 2015; Hirsch et al. 2017; Furlan & Howell 2017). Next, we present a discussion on the speckle imaging analysis of TOI-6651 using two different speckle cameras.

2.2.1 Speckle imaging with PRL 2.5 m telescope

Speckle imaging of TOI-6651 was conducted using the speckle imager on the PRL 2.5m telescope at Mount Abu Observatory in January 2024. The speckle imager is attached to one of the side ports of the telescope. The imager is equipped with a TRIUS PRO-814 CCD detector (model ICX814AL)5 that has a physical dimension of 3388 × 2712 pixels, each pixel measuring 3.69 μm × 3.69 μm. It offers a field of view (FOV) of 2.15′ × 1.7′, with a pixel scale of 38 mas pixel−1, and acquires images in the Bessel-V filter. A sequence of ~2500 speckle frames of TOI-6651, each of 10 ms exposure time, were acquired. The sky conditions during the observation were favorable, with an average seeing close to 1″. A custom data reduction pipeline written in PYTHON was employed to analyze the speckle frames, following a procedure akin to methodologies outlined in Ziegler et al. (2020) and Tokovinin (2018). A brief outline of the steps is discussed below.

The speckle dataset is first dark-subtracted using a master dark derived from a series of dark frames obtained on the same night. Then, each of these dark-subtracted speckle images is cropped to a 128 × 128 pixel region centered on the image centroid. An average background value determined from the area outside the cropped region is subtracted from each frame. Following this, a normalized mean Power Spectrum Density (PSD) is estimated from the squared modulus of the Fourier Transform computed for each frame. The photon noise is determined from the PSD by averaging its values for spatial frequencies f > fc. The cutoff frequency is fcDλ$\[f_c \approx \frac{D}{\lambda}\]$, where D is the telescope diameter, and λ is the central wavelength of the filter. This photon noise value is then subtracted from the observed PSD.

The pipeline next repeats the above steps for a set of speckle frames acquired for a standard star observed close to the target. The PSD of the science target is then divided by the PSD of the standard star to deconvolve the speckle transfer function (STF). However, due to the unavailability of speckle imaging data for a standard star on the night of TOI-6651 observation, we employed a synthetic model of the STF for the deconvolution process, as discussed in Tokovinin et al. (2010). The filtered observed PSD is then fitted with a model to extract both astrometric and photometric parameters (Tokovinin 2018).

An autocorrelation function (ACF), also referred to as an autocorellogram, is computed from the PSD. A 5σ contrast curve is then generated using this ACF. Concentric annuli centered on the primary star are defined within the ACF, and the detection limit is set as 5σ brighter than the local mean within each annulus. The estimated contrast for TOI-6651 is ΔV = 2.76 and ΔV = 4.98 at 0.25″ and 1″, respectively. No stellar companions were detected for TOI-6651. Figure 3 shows the 5σ contrast curve along with the ACF for TOI-6651.

thumbnail Fig. 1

Detrended TESS light curve (LC) from sector 17 (1800-second cadence) and 57 (200-second cadence) shown in green and blue points, respectively. The upper panel displays the full TESS LC plotted against time, while the lower panel presents the phase-folded LC with pink dots representing 25-minute binned data points. The black solid line in both panels represents the best-fit transit model for TOI-6651b from our joint fit analysis (see Section 3.5).

thumbnail Fig. 2

Target pixel file for TOI-6651 in sector 17 (on the left) and sector 57 (on the right) generated with tpfplotter (Aller et al. 2020). The squared region is the aperture mask used in the photometry, whereas the size of the individual red dot is the magnitude contrast (Δm) from TOI-6651. The position of TOI-6651 is marked with “1”.

thumbnail Fig. 3

Contrast curve in V band for TOI-6651 obtained from the speckle imager for PRL 2.5 m telescope (on the left). The contrast curve from SPP speckle analysis in the Ic band (on the right). The speckle ACF is displayed as an inset. No stellar companions are detected.

2.2.2 Speckle interferometry with SPP

TOI-6651 was observed with the SPeckle Polarimeter (SPP; Strakhov et al. 2023) in the Ic band on September 29, 2023. SPP is a facility instrument of the 2.5m telescope at the Caucasian Observatory of Sternberg Astronomical Institute (SAI) of Lomonosov Moscow State University. A fast low-noise CMOS Hamamatsu ORCA-quest was used as a detector, with a pixel scale of 20.6 mas pixel−1, an angular resolution of 89 mas, and a field of view of 5″ × 5″ centered on the star. The power spectrum was estimated from 4000 frames with 30 ms exposure and the atmospheric dispersion compensator was employed. We do not detect stellar companions, with detection limits of ΔI = 4.7 and ΔI = 6.4 at stellocentric distances of 0.25″ and 1.0″, respectively. The final contrast curve for TOI-6651 is shown in Figure 3.

2.3 Precise radial velocities

The RV observations of TOI-6651 were made using PARAS-2 (Chakraborty et al. 2018, 2024) spectrograph attached with the PRL 2.5 m telescope at Mount Abu Observatory. PARAS-2 is a fiber-fed high-resolution (R ≈ 107 000) cross-dispersed echelle spectrograph maintained under highly stable conditions of temperature and pressure and works in the wavelength range of 3800–6900 Å. The spectrograph uses the simultaneous referencing method using the Uranium Argon (UAr) hollow cathode lamp for wavelength calibration and precise RV calculations. Consistently, it has shown off-sky instrumental precision of 30–50 cm s−1 during the course of one night (~12 hours). One can refer to Chakraborty et al. (2018, 2024) for extensive details about the spectrograph. The methodology for the data reduction and RV calculations for the PARAS-2 spectra are described in Section 2.4. We took a total of 27 spectroscopic observations of TOI-6651 spanning 64 days from November 17, 2023, to January 20, 2024, all with an exposure time of 3600 seconds, with S/N at 5500 Å ranging from 15 to 41 and mean S/N of 26. The corresponding photon noise in the RVs is in the range from 2.62 to 13.50 m s−1 with a mean of 6.05 m s−1, calculated using the techniques mentioned in Chaturvedi et al. (2016, 2018). All the details of RVs and their corresponding errors without RV jitter are given in Table A.1. Along with TOI-6651, we have also observed an RV standard star HD 109358 (Bouchy et al. 2013) for a span of ~35 days from December 13, 2023, to January 16, 2024, acquiring a total of 37 RVs. The measured daily RV dispersion is 2.65 m s−1 (see Figure 4).

2.4 PARAS-2 data reduction and analysis pipeline

The PARAS-2 data reduction pipeline is written in IDL and is based on the PARAS-1 pipeline (Chakraborty et al. 2014), which uses the REDUCE package (Piskunov & Valenti 2002) for extraction of cross-dispersed echelle spectra. The package has been adapted and evolved for the PARAS-2 in the same manner it is being done for the PARAS-1 pipeline.

Using the various calibration and science frames, the pipeline is capable of optimally extracting the spectra for science targets. In brief, a master bias frame is generated by a median combining all the bias frames, which is then subtracted from all non-biased frames. A master flat for each fiber is similarly created by median combining all the corresponding flat frames. These master images are further used for echelle order tracing, by finding the order curvature and its locations. The order tracing algorithm associated with REDUCE mainly uses a two-dimensional clustering process, in which it selects the probable order pixels and probes the level of clustering in each of them before merging the most probable orders selected. Since the instrument is very stable, we need not repeat this exercise on a regular basis, instead we use a master order trace. In case of large shifts in the cross-dispersion direction in the flat frames, we need to create the master trace again using the clustering method from REDUCE.

Before doing the optimal reduction, the pipeline does the cosmic ray rejection of the science frames using the Laplacian Edge detection method as described in van Dokkum (2001). It is also important to estimate the scattered light beneath the orders extracted by the pipeline. For this purpose, the pipeline interpolates the background between orders after isolating the inter-order noise. A bleeding map of bright Argon lines is also created by acquiring the single fiber illumination with the UAr lamp. Finally, the spectra are extracted using the swath-based decomposition method.

The spectra need to be wavelength-calibrated for any further analysis. A coarse wavelength mapping of the spectrograph is done using the uranium linelist given in Sharma & Chakraborty (2021) by identifying their pixel positions on the detector. Then, for each order, a fourth-order polynomial is fitted between the central wavelength and the pixel position. This gives the wavelength solution for every order (Sharma & Chakraborty 2021), which is the blueprint of the wavelength solution, available at every instant. The relative wavelength calibration is done using this wavelength solution as the spectrograph is very stable. The wavelength-calibrated spectra of TOI-6651 are cross-correlated with a numerically weighted mask of G2 type star, similarly as described in Baranne et al. (1996) and Pepe et al. (2002).

thumbnail Fig. 4

RVs of the standard star HD 109358 measured over a period of ~35 days. The star has a V-band apparent magnitude of 4.25 and a spectral type of G0V.

3 Result and analysis

3.1 Stellar parameters from spectral synthesis

We estimated the stellar parameters of TOI-6651 using a combined spectrum with S/N ≈ 80, which was generated by co-adding 8 individual high S/N spectra acquired with PARAS-2 (see Section 2.3). The Zonal Atmospheric Stellar Parameters Estimator (ZASPE; Brahm et al. 2017) was used for spectral synthesis. ZASPE compares the co-added spectra with a library of synthetic spectra generated from ATLAS9 model atmospheres (Castelli & Kurucz 2003), finding the best match through least square minimization. Additionally, a high S/N spectrum with S/N ≈ 50 was obtained with the TRES instrument (R ≈ 44 000) on the FLWO 1.5m telescope in September 2023, covering wavelengths from 3850–9096 Å. The stellar parameters were derived using the Stellar Parameter Classification tool (SPC; Buchhave et al. 2010; Buchhave et al. 2012; Buchhave et al. 2014) through cross-correlation with a grid of synthetic spectra based on Kurucz atmospheric models (Kurucz 1992). Both spectra were analyzed individually, and the weighted average of stellar parameter resulted in effective temperature (Teff) of 5933 ± 57 K, surface gravity (log g) of 4.04 ± 0.32, metallicity ([Fe/H]) of 0.213 ± 0.030 dex, and rotational velocity (v sin i*) of 4.88 ± 0.14 km s−1. These parameters indicate that the host star is a sub-giant, metal-rich G-type star. We further use these parameters, along with the broadband spectral energy distribution (SED) and the evolutionary tracks, for a complete stellar characterization, which is elaborated in Section 3.5.

3.2 Galactic kinematics

We calculate the Galactic space velocity components UVW in the barycentric frame utilizing the gal_uvw6 function. These values are listed in Table 1, with UVW being positive in the directions of the Galactic center, Galactic rotation, and the north Galactic pole, respectively. We also provide these velocities with respect to the local standard of rest (LSR) using the solar velocities from Schönrich et al. (2010) (see Table 1). Based on our analysis, we found that TOI-6651 is associated with the thin disk (Leggett 1992; Bensby et al. 2014). The BANYAN ∑ algorithm estimates cluster membership probabilities based on sky positions, proper motions, parallaxes, and RV, classifying TOI-6651 as field stars (>99% chance of membership) without associations to known young clusters (Gagné et al. 2018).

3.3 Rotational period of the star

Using the v sin i* = 4.88 ± 0.14 km s−1, estimated with the spectral synthesis (see Section 3.1) and the stellar radius R* = 1.721 ± 0.065 R, the rotational period of the star is Prot = 17.9 ± 0.8 days, calculated by assuming that star is observed equator-on. We used the archival long-term photometry data from Wide Angle Search for Planets (SuperWASP; Butters et al. 2010) to find the periodically modulated signals attributed to the rotation of the host star. SuperWASP observations were conducted using a broadband filter with a passband ranging from 4000–7000 Å. The dataset utilized for TOI-6651 in this study encompasses 13560 measurements spanning a four-year baseline from June 11, 2004, to July 23, 2008. Since the high cadence TESS data is available for only one sector, we did not account for it for Lomb-Scargle (LS). Due to significant noise in the SuperWASP data, the resulting periodogram is highly uncertain and difficult to interpret. We identified multiple peaks in the LS periodogram exceeding the 0.1% FAP threshold. The most prominent peaks are clustered around a period of 16.5 days, consistent with the rotational period derived from the star’s rotational velocity and radius.

Table 1

Astrometry, photometry, and kinematics properties of TOI-6651.

3.4 Periodogram analysis

We compute the Generalized Lomb-Scargle (GLS) periodogram (Zechmeister & Kürster 2009) on PARAS-2 RV data (panel 1 of Figure 5) and calculate the theoretical 0.1, 1.0, and 10% false alarm probability (FAP) levels, indicated by the horizontal dotted lines. The baseline of our observations is ~64 days, corresponding to a frequency resolution of 1/64 ≈ 0.015 d−1. The most significant period is found to be 5.070 days, which is consistent with the period (≈5.057 days) found in the TESS transit signal. The other two less significant peaks can be seen around ≈0.83 and ≈1.24 days, above the 0.1% FAP level. After fitting for the 5.070 days signal, we do not observe any additional significant peak in the periodogram of residuals (panel 2 of Figure 5), suggesting that these two peaks are aliases of the most significant periodic signal. Furthermore, we compute the periodograms for the Window function, CCF FWHM, and bisector inverse slope (BIS) of PARAS-2 data, and these are shown in panels 3, 4, and 5 of Figure 5. These diagnostics serve as indicators of stellar activity, as they measure line asymmetries resembling Doppler shifts. However, our analysis did not reveal any statistically significant signals of stellar activity in the periodograms.

thumbnail Fig. 5

GLS periodogram for the RVs, residual RVs, window function, FWHM, and bisector span of TOI-6651 are shown in panels 1–5 (top to bottom), respectively. The most significant period in RVs is found to be 5.070 days, denoted by a pink dashed vertical line, which is consistent with the period 5.057 days found in the TESS transit signal. The FAP levels corresponding to 0.1, 1.0, and 10.0% are shown in horizontal blue, orange, and green dashed lines, respectively.

3.5 Global modeling

We determined the stellar and planetary parameters of the TOI-6651 system with simultaneous fitting of the TESS LC and the PARAS-2 RV data using EXOFASTv2 (Eastman et al. 2019), a global modeling software tailored for exoplanet fitting and developed in IDL. This software uses differential Evolution Markov Chain Monte Carlo (dMCMC) analysis to fit multiple datasets simultaneously. The fitting process continues until it meets the standard convergence criteria of EXOFASTv2: number of independent draws, Tz > 1000 and Gelman-Rubin statistic (GR) < 1.01 (Gelman & Rubin 1992; Ford 2006).

3.5.1 Stellar parameters

To model the host star within the EXOFASTv2 framework, a combination of SED (Stassun & Torres 2016) fitting using the Kurucz stellar atmosphere model (Kurucz 1979) and MESA Isochrones and Stellar Tracks (MIST) stellar evolutionary models (Choi et al. 2016; Dotter 2016) were employed. This combination, along with transit photometry, allows for precise estimation of the star’s mass, radius, logg, and age (Torres et al. 2008; Eastman et al. 2023). Therefore, we used the TESS LC of TOI-6651 to model the host star. The broadband photometric magnitudes G, GBP, GRP from Gaia, J, H, K from 2MASS, and W1, W2, W3, W4 from ALLWISE were incorporated into the SED fitting process, as given in Table 1. The uncertainties for these magnitudes are inflated to compensate for systematic floors in measuring accurate absolute photometry, as described in Lindegren et al. (2021). Priors were carefully selected to constrain the parameters effectively. Gaussian priors were set on [Fe/H] based on results from spectral synthesis (see Section 3.1) and on parallax measurements from Gaia DR3 (Gaia Collaboration 2023) with systematic correction from Lindegren et al. (2021). We did not impose a prior constraint on log g and Teff but provided the initial values derived from the spectral analysis (see Section 3.1). Additionally, a uniform prior was imposed to enforce an upper limit on V-band extinction from Schlafly & Finkbeiner (2011) dust maps at the location of TOI-6651.

The resulting posterior distribution of stellar parameters like M* and age show bimodality after all the chains converged. This bimodality arises due to the degeneracy between the MIST isochrones in the Teff − log g plane (e.g., Carmichael et al. 2021; Rodriguez et al. 2021; Khandelwal et al. 2022). The bimodal peaks for TOI-6651 occur at stellar masses of 1.323 M and 1.168 M. We use the splitpdf tool within EXOFASTv2 to separate the posterior distribution into two solutions (referred to as high-mass and low-mass solutions) and report the corresponding parameters in Table 2. The relative probabilities for the high-mass and low-mass solutions are 58.67 and 41.33%, respectively. We adopt the high-mass solution for our further analysis and discussion. The Teff versus log g corresponding to the high-mass solution is marked in Figure 6, where the most likely MIST evolutionary track for TOI-6651 is illustrated as a solid line. The resulting best-fit SED model with broadband photometry fluxes is shown in Figure 7. The stellar parameters derived using EXOFASTv2, with median values and a 68% confidence interval for each posterior (1σ uncertainty), are summarized in Table 2. From our analysis, we find the most probable host star parameters to be: M=1.3230.045+0.050 M,R=1.7210.062+0.065 R$\[M_*=1.323_{-0.045}^{+0.050} ~M_{\odot}, R_*=1.721_{-0.062}^{+0.065} ~R_{\odot}\]$, Teff = 5940 ± 110 K, log g = 4.0870.032+0.035$\[4.087_{-0.032}^{+0.035}\]$ dex, [Fe/H]=0.2250.045+0.044$\[[\mathrm{Fe} / \mathrm{H}]=0.225_{-0.045}^{+0.044}\]$ at an age of 3.710.78+0.68$\[3.71_{-0.78}^{+0.68}\]$ Gyr. These parameters are consistent with those derived from spectral synthesis (see Section 3.1).

thumbnail Fig. 6

MIST evolutionary track for TOI-6651 shown as a solid black line. The black point indicates the Teff and log g, while the red asterisk denotes the current age of TOI-6651.

thumbnail Fig. 7

Spectral energy distribution of TOI-6651. The red symbols show the photometric measurements observed, while the horizontal bars indicate the effective width of the pass band. The blue points represent the model fluxes. The residuals are shown in the lower panel.

3.5.2 Planetary parameters

The joint modeling of TESS photometry with PARAS-2 RVs is done by keeping all orbital fitting parameters (such as b, i, a, Rp, K, e, and ω) as free variables and only supplying initial values for P and Tc, derived from the TESS SPOC pipeline. The transit model is generated using the method outlined by Mandel & Agol (2002), while the RV data are modeled with a standard non-circular Keplerian orbit. We use the default quadratic limb-darkening law for the TESS passband, and its coefficients (u1 and u2) are computed based on tables reported in Claret & Bloemen (2011) and Claret (2017). TESS observations of TOI-6651 are from two sectors with different exposure times, so we resampled the transit model for sector 17 LC, which is 1800-second exposures over 10 steps to account for smearing. The joint modeling reveals a sub-Saturn planet with a most probable radius of RP=5.090.26+0.27 R$\[R_P=5.09_{-0.26}^{+0.27} ~R_{\oplus}\]$ and a mass of MP=61.07.9+7.6 M$\[M_P=61.0_{-7.9}^{+7.6} ~M_{\oplus}\]$, corresponding to the RV semi-amplitude of K = 18.9 ± 2.3 m s−1 in an eccentric orbit (e=0.0910.062+0.096)$\[\left(e=0.091_{-0.062}^{+0.096}\right)\]$. The best-fitted orbital and planetary parameters for both high and low mass solutions are listed in Table 2 with their median values and 68% confidence interval of the posterior (1-σ uncertainty). The TESS LCs and PARAS-2 RVs with their best-fitted models are shown in Figures 1 and 8, respectively.

3.6 Tidal circularization

We estimated the timescales for tidal circularization using equation 2 from Adams & Laughlin (2006), and considering the tidal quality factor Qp = 105, which is within the range of values estimated for the Neptune to Saturn-like planets (Lainey et al. 2009; Šubjak et al. 2022). The stellar tidal quality factor influences the circularization timescale only if the value of Qp is larger than 105.5 (Jackson et al. 2008). For TOI-6651b, the estimated circularization timescale is around 4.4 ± 1.6 Gyr. Compared to the age of the host star 3.710.78+0.68$\[3.71_{-0.78}^{+0.68}\]$ Gyr, these estimates suggest that TOI-6651b could maintain its eccentric orbit to the present day.

Table 2

Summary of EXOFASTv2 fitted and derived parameters with 68% confidence interval for TOI-6651 system.

thumbnail Fig. 8

TOI-6651 RVs (black dots) from 27 spectroscopic observations with PARAS-2 are plotted with time (on the left). TOI-6651 RVs plotted with orbital phase (on the right). In both plots, the best-fit RV model with EXOFASTv2 is represented by the red line, and residuals between the best-fit model and the data are shown in their respective lower panel.

4 Discussion

4.1 Mass radius diagram and internal structure

We plot the mass-radius (M-R) diagram for known exoplanets and display the position of TOI-6651b in Figure 9. The sample of exoplanets is taken from the TepCat database (Southworth 2011), and we only select planets with equilibrium temperature (Teq) between 300 and 2500 K, and with mass and radius constraints better than 50 and 20%, respectively. We have also plotted the M-R curves for pure iron, pure rock, pure high-pressure ices (Zeng et al. 2021), and pure H2 core composition (Becker et al. 2014).

It can be seen in the Figure that TOI-6651b occupies a position between pure high-pressure ices and pure H2 core compositions, suggesting the presence of a significant envelope mass fraction (fenv). Additionally, its high density (ρP=2.520.44+0.52 g cm3)$\[(\rho_P=2.52_{-0.44}^{+0.52} \mathrm{~g} \mathrm{~cm}^{-3})\]$ positions it among the densest known sub-Saturns, indicating a massive core and a smaller fenv. We used the photoevolver7 python code (Fernández et al. 2023) to model the internal structure of TOI-6651b. To first order, sub-Saturns can be approximated as two-component planets consisting of a rocky heavy-element core surrounded by an H/He envelope, as the derived fenv is not sensitive to the detailed composition of the planet core (Petigura et al. 2016). The photoevolver code provides the models (functions) that define the planet’s internal structure as well as its evolution. For internal structure, it defines four crucial parameters: the radius (Rcore) and mass (Mcore) of the core, and the radius (Renv) and fenv of the surrounding envelope. These parameters are related to the total radius RP and mass MP of the planet as RP=Rcore+Renv and fenv=MpMcoreMp$\[R_P=R_{core}+R_{env} \text { and } f_{env}=\frac{M_p-{M}_{core}}{M_p}\]$.

We utilized the empirical M-R relationships provided by Otegi et al. (2020) to infer the core M-R relationship. For the H/He envelope, we adopted the envelope structure model of Chen & Rogers (2016) to relate the Renv and fenv of the surrounding envelope. Given the total mass and radius of the planet and the irradiation flux, the photoevolver employs the above equations and model relations to solve for the internal structure parameters at the given system age. The resulting internal structure of TOI-6651b is summarized in Table 3. These results indicate that ≈87% of the planet’s mass consists of dense materials such as rock and iron in the core, with a core radius of 3.26 ± 0.19 R, which occupies only about 26% of the planet’s volume. The remaining mass comprises a low-density envelope of hydrogen and helium, which occupies most of the planet’s volume.

thumbnail Fig. 9

Planetary radius vs. mass of known exoplanets, with histograms in the right and top panels depicting the distribution of radii and masses, respectively. The sample is from the TepCat database (Southworth 2011) and includes exoplanets with mass and radius constraints better than 50 and 20%, respectively, and equilibrium temperatures of planets (Teq) between 300 and 2500 K. The light pink highlighted region represents the sub-Saturn classification, while the black dashed line represents the iso-density curve for ρP = 2.52. The red, orange, green, and blue solid curves represent the M-R diagram for pure-iron, pure-rock, pure high-pressure ices (Zeng et al. 2021) and pure H2 core composition (Becker et al. 2014), respectively. The Solar system planets are represented as black dots. TOI-6651b is positioned above the pure ice composition, indicating the presence of a significant fraction of a H/He envelope.

Table 3

Internal structure of TOI-6651b.

4.2 The formation scenario

TOI-6651b has a core mass of ≈53 M (87% of the total mass), making it massive enough to trigger runaway gas accretion and become a gas giant planet (Pollack et al. 1996). If this were the case, then TOI-6651b might have undergone significant atmospheric mass loss (see Section 1) earlier in its history. However, due to its high core mass, it is expected to be immune to significant photoevaporation (Petigura et al. 2017), wherein a star’s radiation (X-ray and EUV) strips away a planet’s atmosphere (Owen & Wu 2017). The other possibility is that the planet’s orbit became highly elliptical, and it underwent high-eccentricity migration, during which the orbit was circularized by tidal damping at periastron, leading to the loss of its atmosphere through tidal heating (Tidal Stripping; Owen & Lai 2018; Matsakos & Königl 2016; Vick et al. 2019). Dense and massive sub-Saturns like TOI-6651b are more likely to lose their atmosphere due to tidal heating compared to photoevaporation, which is typically the dominant mechanism for less massive and less dense sub-Saturns (Vissapragada et al. 2022).

It has also been proposed that the formation of massive dense sub-Saturns is likely the result of dynamical instabilities (Petigura et al. 2017; Brady et al. 2018), and TOI-6651b may never have undergone runaway accretion of gas. Initially, TOI-6651b might have had neighboring planets but was situated in a dynamically unstable architecture, eventually resulting in close encounters and planet-planet scattering. Given its close proximity to the star (P ≈ 5 days), TOI-6651b is deeply embedded in the gravitational field of its host star, making scattering events more likely to result in mergers rather than ejections from the system. The maximum velocity a planet can impart to its neighbor is the escape velocity (Goldreich et al. 2004). For TOI-6651b, we calculated the escape velocity and the orbital velocity to be 38.8 ± 2.7 km/s and 137 ± 3 km/s, respectively. Since the escape velocity is smaller than the orbital velocity (ratio ≈ 0.28), single scattering events cannot lead to ejections. Therefore, any previous dynamical instabilities would likely have caused planet mergers, increasing the total mass of TOI-6651b, and the slightly eccentric orbit today might be a relic of that process.

Another possible scenario involves the merger of several sub-critical cores during the disk dispersal stage, and the resulting core accreted the envelope from the gas-depleted disk. Even if the merger took place after the gas disk was fully dispersed and sub-critical cores had their own gas envelopes, then under certain conditions, mergers can form a massive core by breaking apart the lighter bound gaseous envelopes of sub-critical cores (Liu et al. 2015). The formation of super-massive Neptune-sized planet TOI-1853b, through mergers and scattering, was also presented by Naponiello et al. (2023). Although it is difficult to ascertain the precise detailed formation history of the planet, further detection and in-depth studies of such planets can provide more insight into their origins.

thumbnail Fig. 10

Planetary radii vs. orbital period (upper panel) and planetary masses vs. orbital period (lower panel) of known exoplanets from the NASA Exoplanet Archive (Akeson et al. 2013) with mass and radius constraints better than 50 and 20%, respectively. The light pink highlighted region represents the hot Neptunian desert with boundaries defined by Mazeh et al. (2016). The position of TOI-6651b is shown as a pink star.

4.3 Position in the Neptunian desert

The hot Neptunian desert (also known as the sub-Jovian desert) is a region of radius-period and mass-period parameter space where sub-Jovian planets with P < 10 days have been rarely found, and this scarcity cannot be attributed to observational biases (Szabó & Kiss 2011). The probable mechanism for the origin of the Neptunian desert is atmospheric mass loss via photoevaporation and tidal heating following high-eccentricity migration. TOI-6651b (P ≈ 5.06 days) resides at the edge of the desert, considering the boundary up to an orbital period of ~5 days. The latter part of the desert, spanning orbital periods from 5 to 10 days, lacks a comprehensive explanation (Mazeh et al. 2016; Hartman et al. 2019; Esposito et al. 2019; Šubjak et al. 2022; Osborn et al. 2023), which makes TOI-6651b a compelling candidate for studying the origin of the Neptunian desert.

Figure 10 illustrates the position of TOI-6651b in both the mass-period and radius-period parameter spaces along with other known exoplanets, as well as the Neptunian desert with boundaries defined by Mazeh et al. (2016). We find that TOI-6651b is close to the upper boundary of the desert in the mass-period parameter space. Various theories are proposed to elucidate the boundaries of the desert (e.g., Owen & Lai 2018; Vissapragada et al. 2022; Thorngren et al. 2023). Owen & Lai (2018) argue that photoevaporation may account for the lower boundary, while tidal effects following high-eccentricity migration could be responsible for the upper boundary. However, Thorngren et al. (2023) suggests a slightly different perspective, indicating that photoevaporation may partially contribute to carving out the upper boundary. While the detailed discussion of these theories falls beyond the scope of our current study, we contribute by presenting a new hot sub-Saturn located near the upper boundary of the mass-period desert, with a precise mass determination of ≈13%. This finding will be useful for future studies on the functional dependence of the upper boundary on stellar and planetary properties.

thumbnail Fig. 11

Bulk density of known sub-Saturns plotted against their radius (on the left), with TOI-6651b being the third-highest density among them (after Kepler-111c, Dalba et al. 2024, and Kepler-413b, Kostov et al. 2014). Mass of the sub-Saturns as a function of their host star metallicity (on the right), showing a positive correlation between them. The color of individual points is according to the equilibrium temperature of the planet. The position of TOI-6651b in both plots is shown as a pink edge-colored star.

4.4 Position of TOI-6651b in sub-Saturns

In this section, we compare TOI-6651b in the context of other known sub-Saturns. We used data from the NASA Exoplanet Archive (Akeson et al. 2013), and if a planet is studied multiple times, we selected the parameters from the study that provided the most robust constraints on basic planetary parameters. So far8, 74 sub-Saturns have been discovered, with radius and mass constraints better than 20 and 50 percent, respectively. Notably, 22 of these discoveries were made with TESS. Among these sub-Saturns, 41 are classified as “hot sub-Saturns” (P ≤ 10 days). Furthermore, a select few sub-Saturns have densities exceeding 2.0 g cm 3, a rarity in this class of exoplanets. Among these dense sub-Saturns, TOI-6651b emerges as the third-highest density sub-Saturn (after Kepler-111c, Dalba et al. 2024 and Kepler-413b, Kostov et al. 2014) and the most dense among those detected with TESS. The density of known sub-Saturns as a function of their radius is shown in the left panel of Figure 11.

Massive sub-Saturns tend to be found around metal-rich stars, as there were more solids present in the protoplanetary disk to form a massive core. Alternatively, slightly more massive cores can lead to dynamical instabilities, which may result in mergers that form even more massive cores capable of accreting larger gaseous envelopes (Lee & Chiang 2015; Petigura et al. 2017). The correlation between the mass of sub-Saturns and the metallicity of their host stars has been previously observed and discussed by Petigura et al. (2017). In the right panel of Figure 11, we plot the mass of sub-Saturns as a function of their host star’s metallicity. We found that the MP − [Fe/H] positive correlation is consistent for the sub-Saturns detected so far, including the recent discoveries from TESS, in addition to those from Kepler and K2.

4.5 Future follow-up opportunities

Given the equilibrium temperature of planet Teq = 1493 ± 15 K and the surface gravity gP = 22.9 ± 3.6 m s−2, along with the assumption of a H/He atmosphere with a mean molecular mass of μ = 2.3 amu, the atmospheric scale height is calculated to be Hb = 235 ± 37 km using the formula Hb=kTeqμgp$\[H_b=\frac{k T_{e q}}{\mu g p}\]$ (Madhusudhan et al. 2014). TOI-6651b has a low Transmission Spectroscopic Metric (TSM, Kempton et al. 2018) of 19.6 ± 4.3. With a low transit depth of 0.737 ppt and a slowly rotating host with v sin i* = 4.88 km s−1, the Rossiter-McLaughlin (RM) effect (Rossiter 1924, McLaughlin 1924) semi-amplitude for TOI-6651b is estimated to be ≈2 m s−1. Further RV observations and transit-timing variations (TTV) from upcoming TESS data can be used to check for the presence of other planets in the system.

5 Summary

We presented the discovery and characterization of a dense sub-Saturn exoplanet, TOI-6651b, transiting a metal-rich G-type sub-giant star. Our analysis combined RV data from the PARAS-2, transit photometry data from TESS, and speckle observations from the PRL and SAI telescopes. This enabled us to establish the planetary nature of TOI-6651b and determine its mass and radius with ≈13 and ≈5% precision, respectively. Notably, TOI-6651b is the densest sub-Saturn detected with TESS and the third densest among all known sub-Saturns to date. The discovery supports the theory of a positive correlation between planet mass and host star metallicity. TOI-6651b is estimated to have a core mass of ≈53 M, composed predominantly of dense materials like rock and iron, which make around 87% of the planet’s total mass. The remaining mass fraction, with fenv ≈ 0.13, consists of a low-density H/He envelope. The formation history of TOI-6651b is unusual, possibly involving processes such as dynamical instabilities leading to mergers, tidal heating following high-eccentricity migration, and envelope accretion in a gas-depleted disk. Located at the edge of the Neptunian desert, TOI-6651b will be helpful in understanding the factors that define desert boundaries. Further observational and theoretical studies of dense sub-Saturns like TOI-6651b can provide deeper insights into their formation and evolution.

Acknowledgements

We extend our gratitude for the generous support provided by PRL-DOS (Department of Space, Government of India), as well as the Director, PRL. Their support has been instrumental in funding the PARAS-2 spectrograph for our exoplanet discovery project, as well as facilitating research grants for S.B., A.K., and S.N.D. B.S.S. and I.A.S. acknowledge the support of M.V. Lomonosov Moscow State University Program of Development. We express our gratitude to all the Mount Abu Observatory staff for their invaluable assistance throughout the observations. S.B. also acknowledges Mr. Trinesh Sana for his assistance in executing Mathematica code for the M-R diagram. This research has made use of the SIMBAD database and the VizieR catalog access tool, operated at CDS, Strasbourg, France. Additionally, this study utilized the Exoplanet Follow-up Observation Program (ExoFOP; DOI: 10.26134/ExoFOP5) website, managed by the California Institute of Technology under contract with the National Aeronautics and Space Administration as part of the Exoplanet Exploration Program. This paper includes data collected through the TESS mission, obtained from the MAST data archive at the Space Telescope Science Institute (STScI), and utilizes the TepCat and NASA Exoplanet Archive catalog. The authors also express gratitude to the anonymous referee for their numerous helpful suggestions, which significantly improved the quality of the paper.

Appendix A Additional table

Table A.1

Radial velocity measurements of TOI-6651 with PARAS-2(*).

Appendix B Additional figure

thumbnail Fig. B.1

Corner plot showing the covariances for the fitted parameters from EXOFASTv2, corresponding to the high-mass solution (most probable) after splitting the bimodal distribution. The inner and outer contours around the median value represent the 68% and 95% confidence intervals, respectively.

References

  1. Adams, F. C., & Laughlin, G. 2006, ApJ, 649, 1004 [Google Scholar]
  2. Akeson, R. L., Chen, X., Ciardi, D., et al. 2013, PASP, 125, 989 [Google Scholar]
  3. 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]
  4. Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Barragán, O., Aigrain, S., Rajpaul, V. M., & Zicher, N. 2021, MNRAS, 509, 866 [CrossRef] [Google Scholar]
  6. Beaugé, C., & Nesvorný, D. 2013, ApJ, 763, 12 [Google Scholar]
  7. Becker, A., Lorenzen, W., Fortney, J. J., et al. 2014, ApJS, 215, 21 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bensby, T., Feltzing, S., & Oey, M. S. 2014, A&A, 562, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bouchy, F., Díaz, R. F., Hébrard, G., et al. 2013, A&A, 549, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Brady, M. T., Petigura, E. A., Knutson, H. A., et al. 2018, AJ, 156, 147 [NASA ADS] [CrossRef] [Google Scholar]
  11. Brahm, R., Jordán, A., Hartman, J., & Bakos, G. 2017, MNRAS, 467, 971 [NASA ADS] [Google Scholar]
  12. Buchhave, L. A., Bakos, G. Á., Hartman, J. D., et al. 2010, ApJ, 720, 1118 [NASA ADS] [CrossRef] [Google Scholar]
  13. Buchhave, L. A., Latham, D., Johansen, A., et al. 2012, Nature, 486, 375 [NASA ADS] [CrossRef] [Google Scholar]
  14. Buchhave, L. A., Bizzarro, M., Latham, D. W., et al. 2014, Nature, 509, 593 [Google Scholar]
  15. Butters, O. W., West, R. G., Anderson, D. R., et al. 2010, A&A, 520, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Carmichael, T. W., Quinn, S. N., Zhou, G., et al. 2021, AJ, 161, 97 [Google Scholar]
  17. Castelli, F., & Kurucz, R. L. 2003, IAU Symp., 210, A20 [Google Scholar]
  18. Chakraborty, A., Mahadevan, S., Roy, A., et al. 2014, PASP, 126, 133 [CrossRef] [Google Scholar]
  19. Chakraborty, A., Thapa, N., Kumar, K., et al. 2018, SPIE, 10702, 107026G [NASA ADS] [Google Scholar]
  20. Chakraborty, A., Bharadwaj, K. K., Prasad, N. J. S. S. V., et al. 2024, Bull. Liege R. Soc. Sci., 93, 12 [Google Scholar]
  21. Chaturvedi, P., Chakraborty, A., Anandarao, B. G., Roy, A., & Mahadevan, S. 2016, MNRAS, 462, 554 [Google Scholar]
  22. Chaturvedi, P., Sharma, R., Chakraborty, A., Anandarao, B. G., & Prasad, N. J. S. S. V. 2018, AJ, 156, 27 [NASA ADS] [CrossRef] [Google Scholar]
  23. Chen, H., & Rogers, L. A. 2016, ApJ, 831, 180 [Google Scholar]
  24. Choi, J., Dotter, A., Conroy, C., et al. 2016, ApJ, 823, 102 [Google Scholar]
  25. Ciardi, D. R., Beichman, C. A., Horch, E. P., & Howell, S. B. 2015, ApJ, 805, 16 [NASA ADS] [CrossRef] [Google Scholar]
  26. Claret, A. 2017, A&A, 600, A30 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Claret, A., & Bloemen, S. 2011, A&A, 529, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
  29. Cutri, R. M., Wright, E. L., Conrow, T., et al. 2021, VizieR Online Data Catalog: II/328 [Google Scholar]
  30. Dalba, P. A., Kane, S. R., Isaacson, H., et al. 2024, ApJS, 271, 16 [NASA ADS] [CrossRef] [Google Scholar]
  31. Davis, T. A., & Wheatley, P. J. 2009, MNRAS, 396, 1012 [CrossRef] [Google Scholar]
  32. Dotter, A. 2016, ApJS, 222, 8 [Google Scholar]
  33. Eastman, J. D., Rodriguez, J. E., Agol, E., et al. 2019, arXiv e-prints [arXiv:1907.09480] [Google Scholar]
  34. Eastman, J. D., Diamond-Lowe, H., & Tayar, J. 2023, AJ, 166, 132 [NASA ADS] [CrossRef] [Google Scholar]
  35. Esposito, M., Armstrong, D. J., Gandolfi, D., et al. 2019, A&A, 623, A165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Fairnington, T. R., Nabbie, E., Huang, C. X., et al. 2024, MNRAS, 527, 8768 [Google Scholar]
  37. Fernández, J. F., Wheatley, P. J., & King, G. W. 2023, MNRAS, 522, 4251 [CrossRef] [Google Scholar]
  38. Ford, E. B. 2006, ApJ, 642, 505 [Google Scholar]
  39. Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661 [Google Scholar]
  40. Furlan, E., & Howell, S. B. 2017, AJ, 154, 66 [NASA ADS] [CrossRef] [Google Scholar]
  41. Gagné, J., Mamajek, E. E., Malo, L., et al. 2018, ApJ, 856, 23 [Google Scholar]
  42. Gaia Collaboration (Brown, A. G. A., et al.) 2021, A&A, 649, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Gaia Collaboration (Vallenari, A., et al.) 2023, A&A, 674, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Gelman, A., & Rubin, D. B. 1992, Stat. Sci., 7, 457 [Google Scholar]
  45. Goldreich, P., Lithwick, Y., & Sari, R. 2004, ARA&A, 42, 549 [NASA ADS] [CrossRef] [Google Scholar]
  46. Guerrero, N. M., Seager, S., Huang, C. X., et al. 2021, ApJS, 254, 39 [NASA ADS] [CrossRef] [Google Scholar]
  47. Hallatt, T., & Lee, E. J. 2022, ApJ, 924, 9 [NASA ADS] [CrossRef] [Google Scholar]
  48. Hartman, J. D., Bakos, G. Á., Bayliss, D., et al. 2019, AJ, 157, 55 [Google Scholar]
  49. Hirsch, L. A., Ciardi, D. R., Howard, A. W., et al. 2017, AJ, 153, 117 [NASA ADS] [CrossRef] [Google Scholar]
  50. Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 [Google Scholar]
  51. Huang, C. X., Vanderburg, A., Pál, A., et al. 2020a, Res. Notes AAS, 4, 204 [NASA ADS] [CrossRef] [Google Scholar]
  52. Huang, C. X., Vanderburg, A., Pál, A., et al. 2020b, Res. Notes AAS, 4, 206 [NASA ADS] [CrossRef] [Google Scholar]
  53. Jackson, B., Greenberg, R., & Barnes, R. 2008, ApJ, 678, 1396 [Google Scholar]
  54. Jenkins, J., Twicken, J., McCauliff, S., et al. 2016, SPIE, 9913, 99133E [NASA ADS] [Google Scholar]
  55. Jin, S., Mordasini, C., Parmentier, V., et al. 2014, ApJ, 795, 65 [Google Scholar]
  56. Kempton, E. M.-R., Bean, J. L., Louie, D. R., et al. 2018, PASP, 130, 114401 [CrossRef] [Google Scholar]
  57. Khandelwal, A., Chaturvedi, P., Chakraborty, A., et al. 2022, MNRAS, 509, 3339 [NASA ADS] [Google Scholar]
  58. Kostov, V. B., McCullough, P. R., Carter, J. A., et al. 2014, ApJ, 784, 14 [NASA ADS] [CrossRef] [Google Scholar]
  59. Kurucz, R. L. 1979, ApJS, 40, 1 [NASA ADS] [CrossRef] [Google Scholar]
  60. Kurucz, R. L. 1992, IAU Symp., 149, 225 [NASA ADS] [Google Scholar]
  61. Labeyrie, A. 1970, A&A, 6, 85 [Google Scholar]
  62. Lainey, V., Arlot, Jean-Eudes, Karatekin, Özgür, & Van Hoolst, Tim. 2009, Nature, 459, 957 [NASA ADS] [CrossRef] [Google Scholar]
  63. Lecavelier Des Etangs, A. 2007, A&A, 461, 1185 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Lee, E. J., & Chiang, E. 2015, ApJ, 811, 41 [Google Scholar]
  65. Lee, E. J., & Chiang, E. 2016, ApJ, 817, 90 [Google Scholar]
  66. Lee, E. J., Chiang, E., & Ferguson, J. W. 2018, MNRAS, 476, 2199 [NASA ADS] [CrossRef] [Google Scholar]
  67. Leggett, S. K. 1992, ApJS, 82, 351 [Google Scholar]
  68. Lindegren, Klioner, S. A., Hernández, J., et al. 2021, A&A, 649, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  69. Liu, S.-F., Hori, Y., Lin, D. N. C., & Asphaug, E. 2015, ApJ, 812, 164 [CrossRef] [Google Scholar]
  70. Lopez, E. D. 2017, MNRAS, 472, 245 [NASA ADS] [CrossRef] [Google Scholar]
  71. Madhusudhan, N., Knutson, H., Fortney, J. J., & Barman, T. 2014, Exoplanetary Atmospheres (Tucson: University of Arizona Press) [Google Scholar]
  72. Mandel, K., & Agol, E. 2002, ApJ, 580, L171 [Google Scholar]
  73. Matsakos, T., & Königl, A. 2016, ApJ, 820, L8 [Google Scholar]
  74. Mazeh, T., Holczer, T., & Faigler, S. 2016, A&A, 589, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  75. McLaughlin, D. B. 1924, ApJ, 60, 22 [Google Scholar]
  76. Millholland, S., Petigura, E., & Batygin, K. 2020, ApJ, 897, 7 [NASA ADS] [CrossRef] [Google Scholar]
  77. Montalto, M. 2022, MNRAS, 518, L31 [NASA ADS] [CrossRef] [Google Scholar]
  78. Naponiello, L., Mancini, L., Sozzetti, A., Bonomo, A. S., Morbidelli, A., & et al. 2023, Nature, 622, 799 [Google Scholar]
  79. Osborn, A., Armstrong, D. J., Fernández Fernández, J., et al. 2023, MNRAS, 526, 548 [NASA ADS] [CrossRef] [Google Scholar]
  80. Otegi, J. F., Bouchy, F., & Helled, R. 2020, A&A, 634, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Owen, J. E. 2019, Ann. Rev. Earth Planet. Sci., 47, 67 [Google Scholar]
  82. Owen, J. E., & Jackson, A. P. 2012, MNRAS, 425, 2931 [Google Scholar]
  83. Owen, J. E., & Lai, D. 2018, MNRAS, 479, 5012 [Google Scholar]
  84. Owen, J. E., & Wu, Y. 2017, ApJ, 847, 29 [Google Scholar]
  85. Pepe, F., Mayor, M., Galland, F., et al. 2002, A&A, 388, 632 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  86. Petigura, E. A., Howard, A. W., Lopez, E. D., et al. 2016, ApJ, 818, 36 [Google Scholar]
  87. Petigura, E. A., Sinukoff, E., Lopez, E. D., et al. 2017, AJ, 153, 142 [NASA ADS] [CrossRef] [Google Scholar]
  88. Piskunov, N. E., & Valenti, J. A. 2002, A&A, 385, 1095 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  89. Pollack, J. B., Hubickyj, O., Bodenheimer, P., et al. 1996, Icarus, 124, 62 [NASA ADS] [CrossRef] [Google Scholar]
  90. Rodriguez, J. E., Quinn, S. N., Zhou, G., et al. 2021, AJ, 161, 194 [NASA ADS] [CrossRef] [Google Scholar]
  91. Rossiter, R. A. 1924, ApJ, 60, 15 [Google Scholar]
  92. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [Google Scholar]
  93. Schönrich, R., Binney, J., & Dehnen, W. 2010, MNRAS, 403, 1829 [CrossRef] [Google Scholar]
  94. Sharma, R., & Chakraborty, A. 2021, J. Astron. Telesc. Instrum. Syst., 7, 038005 [NASA ADS] [CrossRef] [Google Scholar]
  95. Southworth, J. 2011, MNRAS, 417, 2166 [Google Scholar]
  96. Stassun, K. G., & Torres, G. 2016, ApJ, 831, L6 [NASA ADS] [CrossRef] [Google Scholar]
  97. Stassun, K. G., Oelkers, R. J., Pepper, J., et al. 2018, AJ, 156, 102 [Google Scholar]
  98. Strakhov, I. A., Safonov, B. S., & Cheryasov, D. V. 2023, Astrophys. Bull., 78, 234 [NASA ADS] [CrossRef] [Google Scholar]
  99. Šubjak, J., Endl, M., Chaturvedi, P., et al. 2022, A&A, 662, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  100. Szabó, G. M., & Kiss, L. L. 2011, ApJ, 727, L44 [Google Scholar]
  101. Thorngren, D. P., Lee, E. J., & Lopez, E. D. 2023, ApJ, 945, L36 [NASA ADS] [CrossRef] [Google Scholar]
  102. Tokovinin, A. 2018, PASP, 130, 035002 [Google Scholar]
  103. Tokovinin, A., Mason, B. D., & Hartkopf, W. I. 2010, AJ, 139, 743 [NASA ADS] [CrossRef] [Google Scholar]
  104. Torres, G., Winn, J. N., & Holman, M. J. 2008, ApJ, 677, 1324 [Google Scholar]
  105. van Dokkum, P. G. 2001, PASP, 113, 1420 [Google Scholar]
  106. Vick, M., Lai, D., & Anderson, K. R. 2019, MNRAS, 484, 5645 [NASA ADS] [Google Scholar]
  107. Vissapragada, S., Knutson, H. A., Greklek-McKeon, M., et al. 2022, AJ, 164, 234 [NASA ADS] [CrossRef] [Google Scholar]
  108. Zechmeister, M., & Kürster, M. 2009, A&A, 496, 577 [CrossRef] [EDP Sciences] [Google Scholar]
  109. Zeng, L., Jacobsen, S. B., Hyung, E., et al. 2021, ApJ, 923, 247 [NASA ADS] [CrossRef] [Google Scholar]
  110. Ziegler, C., Tokovinin, A., Briceño, C., et al. 2020, AJ, 159, 19 [Google Scholar]

1

According to the NASA Exoplanet Archive (Akeson et al. 2013) as of July 17, 2024.

8

According to the NASA Exoplanet Archive (Akeson et al. 2013) as of July 17, 2024.

All Tables

Table 1

Astrometry, photometry, and kinematics properties of TOI-6651.

Table 2

Summary of EXOFASTv2 fitted and derived parameters with 68% confidence interval for TOI-6651 system.

Table 3

Internal structure of TOI-6651b.

Table A.1

Radial velocity measurements of TOI-6651 with PARAS-2(*).

All Figures

thumbnail Fig. 1

Detrended TESS light curve (LC) from sector 17 (1800-second cadence) and 57 (200-second cadence) shown in green and blue points, respectively. The upper panel displays the full TESS LC plotted against time, while the lower panel presents the phase-folded LC with pink dots representing 25-minute binned data points. The black solid line in both panels represents the best-fit transit model for TOI-6651b from our joint fit analysis (see Section 3.5).

In the text
thumbnail Fig. 2

Target pixel file for TOI-6651 in sector 17 (on the left) and sector 57 (on the right) generated with tpfplotter (Aller et al. 2020). The squared region is the aperture mask used in the photometry, whereas the size of the individual red dot is the magnitude contrast (Δm) from TOI-6651. The position of TOI-6651 is marked with “1”.

In the text
thumbnail Fig. 3

Contrast curve in V band for TOI-6651 obtained from the speckle imager for PRL 2.5 m telescope (on the left). The contrast curve from SPP speckle analysis in the Ic band (on the right). The speckle ACF is displayed as an inset. No stellar companions are detected.

In the text
thumbnail Fig. 4

RVs of the standard star HD 109358 measured over a period of ~35 days. The star has a V-band apparent magnitude of 4.25 and a spectral type of G0V.

In the text
thumbnail Fig. 5

GLS periodogram for the RVs, residual RVs, window function, FWHM, and bisector span of TOI-6651 are shown in panels 1–5 (top to bottom), respectively. The most significant period in RVs is found to be 5.070 days, denoted by a pink dashed vertical line, which is consistent with the period 5.057 days found in the TESS transit signal. The FAP levels corresponding to 0.1, 1.0, and 10.0% are shown in horizontal blue, orange, and green dashed lines, respectively.

In the text
thumbnail Fig. 6

MIST evolutionary track for TOI-6651 shown as a solid black line. The black point indicates the Teff and log g, while the red asterisk denotes the current age of TOI-6651.

In the text
thumbnail Fig. 7

Spectral energy distribution of TOI-6651. The red symbols show the photometric measurements observed, while the horizontal bars indicate the effective width of the pass band. The blue points represent the model fluxes. The residuals are shown in the lower panel.

In the text
thumbnail Fig. 8

TOI-6651 RVs (black dots) from 27 spectroscopic observations with PARAS-2 are plotted with time (on the left). TOI-6651 RVs plotted with orbital phase (on the right). In both plots, the best-fit RV model with EXOFASTv2 is represented by the red line, and residuals between the best-fit model and the data are shown in their respective lower panel.

In the text
thumbnail Fig. 9

Planetary radius vs. mass of known exoplanets, with histograms in the right and top panels depicting the distribution of radii and masses, respectively. The sample is from the TepCat database (Southworth 2011) and includes exoplanets with mass and radius constraints better than 50 and 20%, respectively, and equilibrium temperatures of planets (Teq) between 300 and 2500 K. The light pink highlighted region represents the sub-Saturn classification, while the black dashed line represents the iso-density curve for ρP = 2.52. The red, orange, green, and blue solid curves represent the M-R diagram for pure-iron, pure-rock, pure high-pressure ices (Zeng et al. 2021) and pure H2 core composition (Becker et al. 2014), respectively. The Solar system planets are represented as black dots. TOI-6651b is positioned above the pure ice composition, indicating the presence of a significant fraction of a H/He envelope.

In the text
thumbnail Fig. 10

Planetary radii vs. orbital period (upper panel) and planetary masses vs. orbital period (lower panel) of known exoplanets from the NASA Exoplanet Archive (Akeson et al. 2013) with mass and radius constraints better than 50 and 20%, respectively. The light pink highlighted region represents the hot Neptunian desert with boundaries defined by Mazeh et al. (2016). The position of TOI-6651b is shown as a pink star.

In the text
thumbnail Fig. 11

Bulk density of known sub-Saturns plotted against their radius (on the left), with TOI-6651b being the third-highest density among them (after Kepler-111c, Dalba et al. 2024, and Kepler-413b, Kostov et al. 2014). Mass of the sub-Saturns as a function of their host star metallicity (on the right), showing a positive correlation between them. The color of individual points is according to the equilibrium temperature of the planet. The position of TOI-6651b in both plots is shown as a pink edge-colored star.

In the text
thumbnail Fig. B.1

Corner plot showing the covariances for the fitted parameters from EXOFASTv2, corresponding to the high-mass solution (most probable) after splitting the bimodal distribution. The inner and outer contours around the median value represent the 68% and 95% confidence intervals, respectively.

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.