The subJovian desert is a region in the massperiod and radiusperiod parameter space that typically encompasses shortperiod ranges between superEarths and hot Jupiters and exhibits an intrinsic dearth of planets. This scarcity is likely shaped by photoevaporation caused by the stellar irradiation received by giant planets that have migrated inward. We report the detection and characterization of TOI3568 b a transiting superNeptune with a mass of 26.4  1.0 MSUBSUB a radius of 5.30  0.27 RSUBSUB a bulk density of 0.98  0.15 g cmSUP3SUP and an orbital period of 4.417965 5 d situated in the vicinity of the subJovian desert. This planet orbiting a K dwarf star with solar metallicity was identified photometrically by the Transiting Exoplanet Survey Satellite TESS. It was characterized as a planet by our highprecision radialvelocity RV monitoring program using MAROONX at Gemini North supplemented with additional observations from the SPICE large program with SPIRou at CFHT. We performed a Bayesian MCMC joint analysis of the TESS and groundbased photometry and MAROONX and SPIRou RVs to measure the orbit radius and mass of the planet as well as a detailed analysis of the highresolution flux and polarimetric spectra to determine the physical parameters and elemental abundances of the host star. Our results reveal TOI3568 b to be a hot superNeptune rich in hydrogen and helium with a core of heavier elements of between 10 and 25 MSUBSUB in mass. We analyzed the photoevaporation status of TOI3568 b and find that it experiences one of the highest extremeultraviolet EUV luminosities among planets with a mass of MSUBpSUB lt 2 MSUBNepSUB yet it has an evaporation lifetime exceeding 5 Gyr. Positioned in the transition between two significant populations of exoplanets on the massperiod and energy diagrams this planet presents an opportunity to test theories concerning the origin of the subJovian desert.
['techniques: photometric', 'techniques: radial velocities', 'planetary systems', 'stars: individual: TOI-3568', 'Astrophysics - Earth and Planetary Astrophysics', 'Astrophysics - Solar and Stellar Astrophysics']
TOI3568 b A superNeptune in the subJovian desert
{'TOI-3568 b: a super-Neptune in the sub-Jovian desert': 'E. Martioli 1 , 2 , R. P. Petrucci 3 , 4 , E. Jofré 3 , 4 , G. Hébrard 2 , 5 , L. Ghezzi 6 , Y. Gómez Maqueo Chew 7 , R. F. Díaz 8 , H. D. Perottoni 1 , L. H. Garcia 3 , D. Rapetti 9 , 10 , A. Lecavelier des Etangs 2 , L. de Almeida 1 , L. Arnold 11 , É. Artigau 12 , R. Basant 13 , J. L. Bean 13 , A. Bieryla 14 , I. Boisse 15 , X. Bonfils 16 , M. Brady 13 , C. Cadieux 12 , A. Carmona 16 , N. J. Cook 12 , X. Delfosse 16 , J.-F. Donati 17 , R. Doyon 12 , E. Furlan 18 , S. B. Howell 9 , J. M. Jenkins 9 , D. Kasper 13 , F. Kiefer 19 , 2 , D. W. Latham 14 , A. M. Levine 20 , D. Lorenzo-Oliveira 1 , R. Luque 13 , K. K. McLeod 21 , J. Melendez 22 , C. Moutou 17 , Y. Netto 22 , T. A. Pritchard 23 , P. Rowden 24 , A. Seifahrt 13 , G. Stefánsson 25 , J. Stürmer 26 , and J D. Twicken 27 , 9  \n(A ffi liations can be found after the references)  \nReceived xxxx ; accepted xxxx', 'ABSTRACT': "The sub-Jovian desert is a region in the mass-period and radius-period parameter space, typically encompassing short-period ranges between super-Earths and hot Jupiters, that exhibits an intrinsic dearth of planets. This scarcity is likely shaped by photoevaporation caused by the stellar irradiation received by giant planets that have migrated inward. We report the detection and characterization of TOI-3568 b, a transiting superNeptune with a mass of 26 . 4 ˘ 1 . 0 M C , a radius of 5 . 30 ˘ 0 . 27 R C , a bulk density of 0 . 98 ˘ 0 . 15 g cm ' 3 , and an orbital period of 4.417965 (5) d situated in the vicinity of the sub-Jovian desert. This planet orbiting a K dwarf star with solar metallicity, was identified photometrically by the Transiting Exoplanet Survey Satellite (TESS). It was characterized as a planet by our high-precision radial velocity monitoring program using MAROON-X at Gemini North, supplemented by additional observations from the SPICE large program with SPIRou at CFHT. We performed a Bayesian MCMC joint analysis of the TESS and ground-based photometry, MAROON-X and SPIRou radial velocities, to measure the orbit, radius, and mass of the planet, as well as a detailed analysis of the high-resolution flux and polarimetric spectra to determine the physical parameters and elemental abundances of the host star. Our results reveal TOI-3568 b as a hot super-Neptune, rich in hydrogen and helium with a core of heavier elements with a mass between 10 and 25 M C . We analyzed the photoevaporation status of TOI-3568 b and found that it experiences one of the highest EUV luminosities among planets with a mass M p ă 2 MNep, yet it has an evaporation lifetime exceeding 5 Gyr. Positioned in the transition between two significant populations of exoplanets on the mass-period and energy diagrams, this planet presents an opportunity to test theories concerning the origin of the sub-Jovian desert.  \nKey words. stars: planetary systems - stars: individual: TOI-3568 - techniques: photometric, radial velocity", '1. Introduction': 'Exoplanets with sizes ranging between Jupiters and super-Earths and orbiting their host stars at short periods ă 5 days are notably rare, giving rise to what is termed the sub-Jovian or Neptunian desert (Szabó & Kiss 2011; Mazeh et al. 2016). The prevailing hypotheses to explain this intrinsic deficit in the planet population involve high-eccentricity migration and stellar irradiation (Owen & Lai 2018). This combination is believed to cause planetary inflation or the erosion of primary atmospheres of exoplanets through photoevaporation, imposing constraints on planet mass and radius depending on their distance from the star.  \nRecent confirmations of several exoplanets residing within the sub-Jovian desert (see Szabó et al. 2023; Kálmán et al. 2023; Frame et al. 2023) showed that this region is not completely barren. Instead, it represents a less probable condition for planets to exist. The sparse presence of planets within the sub-Jovian desert suggests that the formation and evolutionary trajectories of these exoplanets may have taken a unique path, diverging from those observed in more densely populated regions of parameter space. Ongoing surveys optimized for redder stars, such as the Transiting Exoplanet Survey Satellite (TESS, Ricker et al. 2015), are detecting an increasing number of planets in the subJovian desert, particularly around smaller stars (e.g., Eigmüller et al. 2017; West et al. 2019; Eigmüller et al. 2019; Murgas et al. 2021; Mori et al. 2022). Increasing this sample is important to test physical mechanisms involved in planet formation across  \na range of stellar parameters, including spectral types, masses, metallicities, and galactic populations.  \nHere, we present the discovery and characterization of TOI3568 b, a hot super-Neptune orbiting a K-dwarf star located in the sub-Jovian desert, within a transitional region between the populations of hot Jupiters and short-period super-Earths, where planets are notably scarce. This discovery emerged from a program aimed at identifying and characterizing planetary systems within the thick-disk galactic population. Stars from distinct galactic populations exhibit di ff erences in kinematics and chemical composition, potentially leading to variations in the frequency of giant planets as suggested by the planet-metallicity correlation (Fischer & Valenti 2005; Johnson et al. 2010; Buchhave et al. 2014; Wang & Fischer 2015; Ghezzi et al. 2018). We selected TOI-3568.01 as a planet candidate from the TESS Object of Interest (TOI) catalog (Guerrero et al. 2021), based on its high thick disk to thin disk membership probability (TD / D = 3.58) determined through kinematic classification (Carrillo et al. 2020). However, our analysis showed that the nature of this system is more consistent with what is typically associated with the thin disk.  \nThis paper is organized as follows. In Section 2 we present the observations and data reduction; in Sections 3 and 4 we present the characterizations of the star and the planet, respectively; in Section 5 we discuss the characteristics of this new planetary system in the context of the population of exoplanets; and we conclude in Section 6.  \nTable 1. Log of TESS observations of TOI-3568.', '2.1. TESS photometry': "TOI-3568 (TIC 160390955) was first observed by the TESS in sector 15 with a cadence of 30 minutes. A planet candidate with a 4.42 d period was identified in the MIT Quick Look Pipeline (QLP) faint star transit search (Kunimoto et al. 2022). An alert for TOI-3568.01 was issued by the TESS Science O ffi ce on 23 June 2021. We performed a pre-analysis of the 30-minute cadence TESS Full-Frame Image (FFI) data using the community Python package L ightkurve 1 (Lightkurve Collaboration et al. 2018) to obtain the photometric time series. We inspected the light curve data and analyzed the transits, where we employed the methods described in Martioli et al. (2020) and obtained a well-constrained model for the planetary parameters. Thus, we concluded that the events observed by TESS were probably planetary. We submitted TOI-3568 to the TESS Director's Discretionary Targets (DDT 062, PI: E. Martioli), where we were able to include it for observations in sectors 55 and 56 in the 2-min cadence mode. A subsequent search of the 2-min data from sectors 55 and 56 by the TESS Science Processing Operations Center (SPOC) pipeline detected the planetary signature of TOI-3568 b. The di ff erence image centroiding test (Twicken et al. 2018) located the host star to lie within 1 . 1 ˘ 2 . 9 arcsec of the source of the transits. Table 1 shows the log of TESS observations of TOI-3568.  \nFor the observations of TOI-3568 obtained in sectors 55 and 56, we first used the Presearch Data Conditioning (PDC) flux time series (Smith et al. 2012; Stumpe et al. 2012, 2014) processed by the SPOC pipeline at NASA Ames Research Center (Jenkins et al. 2016) obtained from the TESS data products available in the Mikulski Archive for Space Telescopes (MAST) 2 . Figure 1 shows the TOI-3568's target pixel files (TPF) for sectors 15, 55, and 56. It highlights the pixels used in the aperture for photometry and marks the positions of sources in the field from Gaia DR3's catalog (Gaia Collaboration et al. 2016, 2023; Katz et al. 2023).  \nWecalculated our final TESS photometry using an optimized systematic error correction algorithm, following the methodology of Rapetti et al. (in prep.). We use an adaptation of the Pixel Level Decorrelation (hereafter PLD; Deming et al. 2015; Luger et al. 2016, 2018) technique implemented in the PLDC orrector class of L ightkurve . This method employs (i) a spline polynomial fit to describe stellar variability; (ii) Principal Component Analysis (PCA) eigenmodes to model the background light; and (iii) the PLD technique to account for pointing and mechanical e ff ects.  \nTo account for the background as described in (ii), we wish to begin with calibrated pixels that are not corrected for the background. We thus add the background flux estimated by SPOC into the calibrated and background-removed pixel values in the original TPF before the correction. Since PLD might preserve the mean of the uncorrected light curve after the regression, to re-", 'TIC 160390955 Sector 15': "TIC 160390955 Sector 55  \n<!-- image -->  \nFig. 1. Target pixel files for TOI-3568 observed in TESS sectors 15 (top panel), 55 (middle panel), and 56 (bottom panel). The red circles represent the Gaia DR3 sources, with TOI-3568 marked by a cross and labeled as index 1. Symbol sizes correlate with the Gaia DR3 source's brightness compared to TOI-3568's brightness. The red-shaded squares indicate the pixels used in aperture photometry. The three images were produced with tpfplotter (Aller et al. 2020).  \n<!-- image -->  \n<!-- image -->  \nPixel Column Number  \ncover the true mean flux level of the corrected light curve we apply a flux level adjustment. We adjust the corrected flux by subtracting a constant level calculated as the third dimmest median pixel flux value times the number of pixels. This is inspired by the scalar background bias method applied in the SPOC pipeline. We then adjust the flux for crowding by non-target stars and for the fraction of the target star flux captured in the photometric aperture using the methods of and the crowding and flux fraction values provided by the SPOC pipeline. We also use the flux fraction to scale the flux errors, but not the crowding since its e ff ect on the flux errors is negligible.  \nBefore applying the PLD corrector, we add the background flux and errors estimated by the SPOC pipeline back onto the Simple Aperture Photometry (SAP) light curve. Flux level, fraction and crowding adjustments are applied to the corrected light curves. To automatically optimize the selection of parameter values for the corrector, we evaluate the resulting light curve using the Savitsky-Golay Combined Di ff erential Photometric Precision (sgCDPP) proxy algorithm (Gilliland et al. 2011; Van Cleve et al. 2016) implemented in L ightkurve , for durations of 30, 60, 120, 160, and 200 minutes. For a grid of corrector parameter values (for further details on the parameters and the grid, see Rapetti et al. (in prep.)), we calculate the harmonic mean (HM) of these quantities and select the corrected light curve that minimizes the HM. Using this analysis, we were able to recover segments of the data that were initially excluded by the SPOC pipeline due to scattered light. The TESS light curves are illustrated in Figure 2 along with the results of our model, as outlined in Section 4.1.", '2.2. Ground-based photometry': "Ground-based time-series photometry was collected as part of the TESS Follow-up Observing Program Sub Group 1 3 (Collins 2019) which uses a customized TESS version of the TAPIR software package (Jensen 2013) for observation planning. The data reduction and aperture photometry were performed using AstroImageJ (Collins et al. 2017) and the light curves are available on ExoFOP 4 .  \nWe observed a full transit event of TOI-3568 b on 2022-1213 using the 0.7 m telescope at the Wellesley College Whitin Observatory 5 (WCWO) in MA, USA. Images were taken in an Sloan Digital Sky Survey (SDSS) ' r 1 filter using 30 s exposures, and photometry was extracted using a circular aperture with a radius of 2.7 arcsec. The aperture was small enough to exclude the light from the two nearest GAIA DR3 stars with projected separations of 3.7 arcsec and 5.7 arcsec, the latter of which is bright enough to have caused the TESS event if it had been an eclipsing binary. These uncontaminated r 1 data are used along with the TESS light curves in the joint analysis in Section 4.", '2.3. High contrast imaging': "High-angular resolution observations of candidate systems hosting transiting exoplanets can aid in identifying blended sources within sub-arcsecond scales. These sources might produce a false positive transit signal, particularly if the source is an eclipsing binary. The EXOFOP-TESS website 6 reports five high contrast imaging observations, spanning from the optical range at 562 nm to the near-infrared (NIR) at 2.2 µ m, as summarized in  \nFigure 3 shows the 5σ sensitivity curves derived from observations made with the 'Alopeke dual-channel speckle imaging instrument on Gemini-N (PI: Howell). These observations were obtained with a pixel scale of 0.01 arcsec per pixel and a full width at half maximum (FWHM) resolution of 0.02 arcsec. The data were processed with the speckle pipeline (Howell et al. 2011). 'Alopeke performed simultaneous speckle imaging at 562 (54) nm and 832 (40) nm. The results from these observations e ff ectively eliminate the possibility of any companion with a contrast of ∆ mag ă 6 . 42 at 0.5 arcsec separation at 832 nm. This wavelength range is particularly pertinent to the photometry and spectroscopic observations presented in this paper, as it overlaps with the spectral sensitivity of TESS and MAROON-X.  \n<!-- image -->  \n<!-- image -->  \nFig. 2. TESS light curve of TOI-3568. The gray points show the TESS photometry data normalized by the median flux of each sector. The black points show the same data binned by the weighted average with bin sizes of 0.1 d. The red line shows the baseline GP model that was fitted to the binned data multiplied by the best-fit transit model for TOI3568 b. The red triangles show the central time of TOI-3568 b transits.  \n<!-- image -->  \nTable 2. These observations indicate the absence of a close-in companion with su ffi cient brightness to generate a false positive signal.  \nTable 2. Summary of high-contrast imaging observations of TOI-3568.  \nFig. 3. Contrast sensitivity for TOI-3568 as a function of angular separation at 562 nm (blue line) and at 832 nm (red line) obtained from the 'Alopeke / Gemini-N speckle imaging observations. The inset panel displays the reconstructed speckle image at 832 nm.  \n<!-- image -->", '2.4. MAROON-X spectroscopy': "MAROON-X is a high-resolution spectrograph ( λ { ∆ λ ' 85000) operating in the optical range (500-920 nm) and installed on the 8.1-m Gemini North telescope atop Maunakea, Hawaii (Seifahrt et al. 2018). The spectrograph is fiber-fed, highly stabilized, and bench-mounted. It was designed to achieve sub-m s ' 1 radial velocity (RV) precision.  \nWe obtained 34 spectra of TOI-3568 under program ID GN2022A-Q-207 / -Q-113 (PI: R. Petrucci) with an exposure time of 720 s and using MAROON-X in its single mode of operation. These spectra were collected over ten di ff erent nights spanning from 2022-04-08 to 2022-07-26. The strategy of obtaining ' 3 spectra per visit was precautionary, to mitigate the impact of potential outliers. On average, these observations yielded a peak signal-to-noise ratio (S / N) per spectral element of 51 ˘ 11 in the blue arm and 70 ˘ 16 in the red arm.  \nThe MAROON-X raw data have been reduced using the standard procedure implemented in the instrument Python3 pipeline (Seifahrt et al. 2020). This procedure involved bias and background subtraction, order tracing and the extraction of one-dimensional wavelength-calibrated spectra. Wavelength solutions and instrumental drift corrections were based on the simultaneous calibration data of a stabilized Fabry-Pérot etalon (Stürmer et al. 2017), which allows order-by-order drift corrections at the sub-m s ' 1 level. The flux-weighted midpoint of each observation was used to calculate the barycentric corrections.  \nWe analyzed the spectra using the SpEctrum Radial Velocity AnaLyser (SERVAL) pipeline (Zechmeister et al. 2018), which employs the template-matching algorithm to extract precise relative RVs. The blue and red channels of MAROON-X are reduced separately, producing independent RVs, as presented in Table 3. The blue channel RVs show a root mean square (RMS) disper-  \nTable 3. MAROON-X relative radial velocities for the red and blue channels.  \nsion of 9 . 7ms ' 1 and a median error of 1 . 9ms ' 1 , whereas the red channel shows an RMS dispersion of 9 . 5ms ' 1 and a median error of 3 . 5ms ' 1 .  \nWe co-added all individual MAROON-X observations shifted to the same stellar reference frame to obtain a master spectrum with high S / N. Our spectroscopic analysis in Section 3 uses this master spectrum.  \nIn addition, we employ a reference solar spectrum obtained from observations of sunlight reflected by the asteroid Vesta on the night of 2022-04-27 under program ID GN-2022A-Q-227 (PI: Y. Netto). These observations were carried out adopting the same MAROON-X setup (S / N ' 400 at 600 nm) to ensure precision in determining the stellar parameters and chemical abundances of TOI-3568 through a di ff erential analysis.", '2.5. SPIRou spectro-polarimetry': "TOI-3568 was observed by the SpectroPolarimètre InfraRouge (SPIRou) 7 under the large program SPIRou Legacy Survey - Consolidation & Enhancement (SPICE 8 ; PI: JeanFrançois Donati) on nights spanning from 2022-11-14 to 2022-  \nℓ  \nTable 4. SPIRou relative radial velocities and the B data.  \n11-21. SPIRou is a stabilized high-resolution near infrared (NIR) spectropolarimeter (Donati et al. 2020) mounted on the 3.6 m Canada-France-Hawaii Telescope (CFHT) atop Maunakea, Hawaii. It is designed for high-precision velocimetry to detect and characterize exoplanets and it provides a full coverage of the NIR spectrum from 950 nm to 2500 nm at a spectral resolving power of λ { ∆ λ ' 70000.  \nWe observed TOI-3568 with SPIRou / CFHT at five di ff erent epochs, where we obtained a total of 20 spectra with an individual exposure time of 900 s. These observations were carried out in the circular polarization mode (Stokes V), where each set of four exposures provides a polarimetric spectrum. The peak S / N per spectral element varied between 40 and 65, with a median of 59. The air mass of the observations ranged from 1.1 to 1.4 and the Barycentric Earth Radial Velocity (BERV) ranged from -20.0 to 20.7 km s ' 1 .  \nThe SPIRou data have been reduced using the APERO pipeline v.0.7.284 (Cook et al. 2022), which produces 1D optimally extracted fluxes that underwent detector gain and artifact corrections, wavelength calibration, blaze correction, and correction for telluric atmospheric absorption. Additionally, the pipeline computed polarimetric Stokes V and null spectra.  \nThe flux spectra have been analyzed using the line-by-line (LBL) method of Artigau et al. (2022), wherein a high-S / N template spectrum of HD 189733 observed by SPIRou was used as a reference to obtain the RVs. The Table 4 shows the SPIRou RVs. These RVs show an RMS of 12.8 m s ' 1 , and a median error of 12.4 m s ' 1 .", '3. Stellar characterization': 'We carried out a study to derive the host star properties as will be detailed in the next sections. Table 5 presents a summary of the stellar parameters of TOI-3568.', '3.1. Atmospheric parameters': "The fundamental atmospheric parameters (Te ff , log g , [Fe / H], and vt ) of TOI-3568 were determined by imposing a strictly line-by-line di ff erential spectroscopic equilibrium of neutral and singly-ionized iron lines relative to the Sun (e.g., Jofré et al. 2021). To perform this process automatically, we employed the q2 program 9 (Ramírez et al. 2014). The iron line list, as well as the atomic parameters, namely the excitation potential (EP) and oscillator strengths (log gf ), are the same as in Jofré et al. (2021), and the equivalent widths (EWs) of the MAROON-X spectra of TOI-3568 and the Sun (reflected sunlight from Vesta) were manually measured by fitting Gaussian profiles using the splot task in IRAF 10 . The solar values were kept fixed at (Te ff ,  \nlog g , [Fe / H], vt ) = (5777 K, 4.44 dex, 0.0 dex, 1.0 km s ' 1 ). The resulting stellar parameters are Te ff = 4969 ˘ 45 K, log g ' 4 . 63 ˘ 0 . 08 dex, r Fe { H s ' ' 0 . 01 ˘ 0 . 02 dex, and vt ' 0 . 77 ˘ 0 . 13 km s ' 1 .  \nAs an independent check, we derived fundamental parameters also using the excitation and ionization equilibria of Fe I and Fe II lines technique but without using a strictly line-byline di ff erential analysis (e.g., Ghezzi et al. 2018, 2021; Jofré et al. 2020). We measured the EWs automatically with the ARES code (e.g., Sousa et al. 2015) and the pipeline from Ghezzi et al. (2018) that uses the 2017 version of the MOOG code 11 (Sneden 1973) and a line list with solar log gf values derived using the same setup (see Ghezzi et al. 2018 for further details). In excellent agreement with the parameters obtained with the di ff erential technique we obtained: Te ff ' 4897 ˘ 47 K, log g ' 4 . 52 ˘ 0 . 11 dex, r Fe { H s ' 0 . 02 ˘ 0 . 02 dex, vt ' 0 . 52 ˘ 0 . 16 km s ' 1 .  \nWedetermined the projected star rotation velocity ( v rot sin i ‹ ) based on the spectral synthesis of relatively isolated iron lines using the code i S pec (Blanco-Cuaresma 2019) and following the procedure of Carlberg et al. (2012). Adopting the calibration of Takeda & UeNo (2017) to determine a macroturbulence velocity of 1.06 km s ' 1 , we find v rot sin i ‹ ' 1 . 4 ˘ 0 . 4 kms ' 1 . However, considering the resolving power of MAROON-X (R = 85,000), we adopt an upper limit of 2 km s ' 1 .", '3.2. Mass, radius, and age': "We derived stellar mass, radius, and age using Yonsei-Yale (YY) stellar isochrones (Yi et al. 2001; Demarque et al. 2004), as described in Jofré et al. (2021). This was accomplished via the q2 pipeline, using as input the spectroscopic Te ff and [Fe / H] obtained from the di ff erential method, Gaia DR3 parallax, and V -mag (corrected for extinction 12 ). We obtained an age of 6 . 1 ˘ 3 . 7 Gyr, a mass of 0 . 780 ˘ 0 . 021 M d , and a radius of 0 . 720 ˘ 0 . 013 R d . This code also provides the trigonometric gravity, which allows us to perform a consistency check on the spectroscopic log g values. Here, q2 yields log g ' 4 . 62 ˘ 0 . 02 dex, which is in excellent agreement with the estimates found from the spectroscopic equilibrium.  \nAs a consistency check, we also employed the 1.3 version of the PARAM web interface 13 that performs a Bayesian estimation of stellar parameters (da Silva et al. 2006; Miglio et al. 2013) based on PARSEC isochrones (Bressan et al. 2012). As input, we employed the same parameters as those used above in q2 . In good agreement, within the errors, with our estimations from q2 , PARAM returned an age of 4 . 2 ˘ 3 . 7 Gyr, a mass of 0 . 75 ˘ 0 . 02 M d , a radius of 0 . 69 ˘ 0 . 01 R d , and log g ' 4 . 60 ˘ 0 . 02 dex.", '3.3. Longitudinal magnetic field': "We performed a least squares deconvolution (LSD) analysis and computed the longitudinal magnetic field ( B ℓ ) on individual SPIRou spectra for Stokes I, Stokes V, and null polarizations,  \nTable 5. Summary of the stellar parameters for TOI-3568.References. (1) Gaia Collaboration (2020); (2) Stassun et al. (2018, 2019), sourced from the EXOFOP-TESS website a ; (3) Cutri et al. (2003); (4) Wright et al. (2010).  \nusing the methodologies introduced by Donati et al. (1997) and implemented by Martioli et al. (2020). The resulting mean LSD profiles are depicted in the top panels of Figure 4, while the time series of B ℓ is illustrated in the bottom panel. Table 4 shows the values of B ℓ . The Stokes V profile is featureless, indicating the absence of a Zeeman signature for this star. The average longitudinal magnetic field, B ℓ ' ' 2 . 5 ˘ 14 . 6 G, is consistent with a null detection of a magnetic field in TOI-3568, suggesting its magnetic inactivity.", '3.4. Activity': "To investigate activity in TOI-3568, we measured three spectral index proxies for chromospheric activity in the MAROON-X spectra: the Ca infrared triplet (Ca IRT), the Hα (available in both the blue and red channels), and the Na D1 and D2 doublet. We found no significant correlations between these indices and  \nRVs, nor with the FWHM in either the blue or red channels, indicating that the RV variations are likely not caused by stellar activity. We measured the median values and standard deviations of the CaIRT and NaD indices as follows: CaIRT1 ' 0 . 446 ˘ 0 . 004, CaIRT2 ' 0 . 332 ˘ 0 . 004, CaIRT3 ' 0 . 447 ˘ 0 . 005, NaD1 ' 0 . 175 ˘ 0 . 004, and NaD2 ' 0 . 206 ˘ 0 . 002. For Hα , we measured 0 . 281 ˘ 0 . 005 in the blue channel and 0 . 283 ˘ 0 . 007 in the red channel, resulting in a mean value of 0 . 282 ˘ 0 . 004. This represents a variation in Hα of 1.4% over a baseline of 100 days, indicating that TOI-3568 appears to have low levels of chromospheric activity.  \nWe also looked for signs of stellar variability in the 2-min TESS PLD data from sectors 55 and 56. To do so, we carried out the analysis of the light curve with the tools in the lightkurve package. As a first step, based on the transit parameters obtained in our analysis (see Section 4), we removed all the points falling within the transits of TOI-3568 b. Then, we ran two algorithms on the resulting light curve: the Lomb-Scargle periodogram (LS;  \nFig. 4. LSD analysis and longitudinal magnetic field from SPIRou spectropolarimetric data. The top panel shows the Stokes I (red points) with a Voigt profile model fit (green line), Stokes V (blue points), and null (dark cyan points) LSD mean profiles. The absence of the characteristic 'S'-shape in the Stokes V profile indicates a non-detection of a Zeeman signature. In the bottom panel, the longitudinal magnetic field (green points) measured across the five epochs observed by SPIRou displays an average value of B ℓ ' ' 2 . 5 ˘ 14 . 6 G, which is also consistent with a null value.  \n<!-- image -->  \nScargle 1982) and, an additional independent method, the AutoCorrelation Function (ACF; McQuillan et al. 2013). After applying the criteria described in Petrucci et al. (2024) to assess if a detected signal is real, we determined that no significant peak indicating periodic variability was found and, hence, there is no evidence for rotational modulation in the data. This may be a consequence of a small spot coverage on the stellar surface, or due to the existence of a rotation period longer than ' 14 d, which is likely undetected in the TESS data. Additionally, no flare candidate was detected on the 2-min cadence light curve by the A ltaipony code (Davenport 2016; Ilin et al. 2021), optimized to search for sporadic events. A careful by-eye inspection of the TESS photometry confirms these results.  \nThe absence of signs of variability suggests that TOI-3568 is an inactive star. This aligns with the mature age value obtained in this study (see Table 5) and with the low levels of chromospheric and magnetic activity measured from our spectropolarimetric data.", '3.5. Chemical composition': "We measured line-by-line di ff erential abundances relative to solar ([X / H]) abundances of 20 elements other than iron, including C, O, Na, Mg, Al, Si, S, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, and Ba. This was achieved through EW measurements and by employing the curve-of-growth approach with the MOOG program ( abfind driver) using the q2 code. The EWs were manually measured using the splot task in IRAF and the adopted line list and atomic parameters were taken from Jofré et al. (2021). Hyperfine splitting was taken into account for Sc, V, Mn, Co, Cu Y, and Ba. The O abundance was computed from the 7771-5 Å IR triplet, adopting the non-LTE corrections by Ramírez et al. (2007).  \nSimilar to the previous section, as an independent check, we also performed a non-strictly di ff erential analysis to determine elemental abundances for TOI-3568 following Schuler et al. (2015) and Teske et al. (2015). The EWs were automatically measured for elements with more than three spectral lines using the code ARES (Sousa et al. 2015) and an updated line list based on solar log gf values that will be described in a future paper 14 that investigates possible correlations between planetary properties and the chemical abundances of their stellar hosts. For elements with three or less lines, we manually measured the EWs with the splot task in IRAF. The abundances were calculated using the MOOG code and abfind and blends drivers for the elements without and with hyperfine splitting, respectively. The blends driver was also used to determine the abundance from the [O I] 6300 Å line in order to take the contamination from Ni into account. The non-LTE corrections for the O I 7771-5 Å IR triplet were taken from Amarsi et al. (2019).  \nTable 6 lists both line-by-line di ff erential abundances and non-di ff erential values. The uncertainties were determined considering the standard deviation of the mean abundances (for elements with three or more lines) as well as the contributions of the uncertainties on the atmospheric parameters. In order to calculate these contributions, we vary each atmospheric parameter in turn by ˘ 1 σ and calculate new abundances for all elements. Wethen subtract these new values from the original ones and obtain two di ff erences caused by the variation of each atmospheric parameter. The contribution of each atmospheric parameter for the final uncertainty is taken as the maximum of these two differences. Finally, we add in quadrature the contributions of each of the four atmospheric parameters as well as the standard deviation of the mean abundances. For elements with only one or two lines, the latter is not considered. Figure 5 shows the abundances as a function of atomic number. There is a very good agreement between both sets of results, which are consistent within 1 σ . We also notice the good agreement between the O abundances obtained from di ff erent indicators as well as between neutral and ionized species for Sc, Ti and Cr, which further shows that the ionization equilibrium achieved during the determination of the atmospheric parameters is robust.  \nWe also determined the lithium abundance by performing a spectral synthesis of the Li I feature at 6707.8 Å. We adopted the line list from Ghezzi et al. (2010) and a similar methodology, but with two di ff erences: (1) we determined the Gaussian broadening (that considers the combined e ff ects of the instrumental profile, stellar rotation and macroturbulence) using the Fe I line at 6703.567 Å and (2) we kept the abundances of Fe, C and Si as free parameters (using the previously determined values as ini-  \nTable 6. Strictly di ff erential abundances (relative to the Sun) and non-di ff erential abundances of TOI-3568  \nFig. 5. Elemental abundances as a function of atomic number for TOI3568 from the analysis of MAROON-X spectra. The methods are described in the legend, and each element's symbol is depicted in the plot. Dashed line marks the weighted average metallicity obtained from the di ff erential abundances and the shaded area indicates the standard deviation.  \n<!-- image -->  \nal guesses) for the fit. As we can see in Figure 6, the best fit shows no distinguishable Li feature and the residuals are smaller than 0.5%. We are only able to determine an upper limit of A(Li) ď 0.45 and this value is consistent with those of stars with similar e ff ective temperatures (e.g., Ghezzi et al. 2010; Ramírez et al. 2012).", '3.6. Stellar population membership': "As we mention on the introduction, TOI-3568 was classified as a thick disk star by Carrillo et al. (2020). To independently check this previous classification, we determined the galactic population membership in two ways. First, we performed a kinematic classification by computing the Toomre diagram (e.g., Bensby  \nFig. 6. Spectral synthesis for the Li I feature at 6707.8 Å for TOI3568. Top panel: Best fit between the observed (red circles) and synthetic (black line) spectra. The lithium line is not detected and we determined an upper limit of A(Li) ď 0.45. Bottom panel: Residuals of the fit, which are within 0.5%.  \n<!-- image -->  \net al. 2003). To identify regions dominated by each galactic population in the v ϕ vs. b v 2 R ' v 2 z plane, we utilized the G alaxia model (Sharma et al. 2011) with updated velocity distributions and local fractions of the thick disk and stellar halo (Amarante et al. 2020). This is illustrated in Figure 7, where the thin disk, thick disk, and halo are represented by red, grey, and blue shaded areas, respectively. The kinematic characteristics of TOI-3568 are consistent with those typically associated with nearby thindisk stars. We also applied a formalism similar to that of Bensby  \nFig. 7. Toomre diagram for the galactic components based on the G alaxia model. TOI-3568 is represented by the yellow star. The thin disk, thick disk, and halo are represented by the red, gray and blue colours, respectively. The varying darkness of shading indicates the fraction of each population at specific locations on the diagram (75%, 85% or 95%).  \n<!-- image -->  \net al. (2003) (see Perottoni et al. 2021 for more details) to kinematically classify the galactic component of TOI-3568. Stars are classified as thick disk in this method when they show thickdisk-to-thin-disk (TD / D) membership ratios higher than 10, a condition not met by TOI-3568 ( TD { D ' 1 . 9).  \nFurthermore, we estimated the orbital properties of TOI3568, following the description provided by Perottoni et al. (2021). TOI-3568 is on prograde orbit with a high angular momentum Lz ' 1 . 2 ˆ 10 3 kpc km s ' 1 , mild eccentricity e ' 0 . 37, and it reaches Z max ' 0 . 08 kpc from the galactic plane. Additionally, TOI-3568 is situated in a region of E -LZ space predominantly occupied by disk stars. Although TOI-3568 exhibits an uncommon eccentricity, the other orbital parameters are similar to those of thin disk stars. Understanding the mechanism that led to this star's eccentricity is beyond the scope of this paper.  \nOnthe other hand, we also performed a classification of TOI3568 based on its chemical composition. Generally, thick disk stars are metal-poor and enhanced in α elements (e.g., Fuhrmann 1998; Reddy et al. 2006; Haywood 2008; Adibekyan et al. 2012). Figure 8 shows the [ α { Fe] 15 versus [Fe / H] for TOI-3568 in comparison with 1111 FGK dwarfs observed within the context of the HARPS GTO planet search program for which precise abundances of α -elements are available (Adibekyan et al. 2012). In this figure the objects are chemically separated, by the dashed line, between thick disk stars, overabundant in α elements and thin disk stars with less content of α -elements for a given [Fe / H].  \nIn line with the kinematic classification, the α content of TOI-3568, with a [ α { Fe] value of 0 . 061 ˘ 0 . 029 dex from the di ff erential method (and r α { Fe s ' 0 . 05 ˘ 0 . 06 dex for the nondi ff erential approach), does not exhibit significant enhancement for its metallicity. Therefore, TOI-3568 falls within the thin-disk region, albeit near the transition zone. For comparison, in Figure 8 we include bona fide thick disk stars with transiting planets for which there are alpha-element abundances available derived from high-resolution optical spectra.  \nAdditional support to our kinematic and chemical Galaxy population classification of TOI-3568 is based on its age. Gen-  \nFig. 8. Iron abundance vs. [ α { Fe] for stars in the HARPS GTO program. The position of TOI-3568 is marked by a blue star (di ff erential abundances) and by a green star (non-di ff erential abundances). The black dashed line separates the galactic thin-disk (bottom) and thick-disk populations (top). For comparison, we also include the bona fide thick disk stars with transiting planets Kepler-444 (Campante et al. 2015), TOI561 (Weiss et al. 2021), K2-111 (Mortier et al. 2020), and WASP-21 (Bouchy et al. 2010).  \n<!-- image -->  \nerally, thick-disk field stars of the Milky Way are older than about 10-11 Gyr (e.g., Bensby et al. 2005; Reddy et al. 2006; Adibekyan et al. 2011). Hence, the age of TOI-3568, estimated from the isochrones analysis (approximately in the range of 19 Gyr considering the errors in both methods, Sec. 3.2) would be more consistent with those of thin-disk stars ( À 9 Gyr, e.g., Bernkopf et al. 2001; Nissen 2004; Kilic et al. 2017). Moreover, as a further check to the ages from isochrones estimated above, we utilize the age-[Y / Mg] relation from Tucci Maia et al. (2016) to obtain the age of TOI-3568 employing our measured abundances. This yields an age of 7 . 6 ˘ 1 . 2 Gyr, consistent with the values derived from the isochrones analysis. Therefore, our kinematic and chemical analysis indicate that TOI-3568 is a thin disk (or thin / thick disk transition object) rather than a thick disk star as indicated previously by Carrillo et al. (2020).", '4.1. Analysis of TESS photometry data': "We analyzed the TESS data from sectors 15, 55, and 56 using methods outlined in Martioli et al. (2021). This involved fitting a transit model along with a baseline polynomial within selected windows around each transit of TOI-3568 b. We selected six transits with a 30-minute cadence from sector 15 and thirteen transits with a 2-minute cadence from sectors 55 and 56. On average, each transit window contained 32 data points in sector 15 and 468 data points in sectors 55 and 56. This window size includes approximately twice as many out-of-transit data points on each side as in-transit ones. This was found to be a reasonable balance for a non-active star, providing enough baseline data to constrain transit parameters accurately and to allow for precise modeling of the baseline using a low-order polynomial. In this case, we use a first-order polynomial. As illustrated in Figure 2, we have also modeled the TESS photometry data with a Gaussian Process (GP) regression using a quasi-periodic kernel as in Martioli et al. (2022, 2023). The primary purpose of employing the GP in this analysis is for detrending, as it does not constrain any significant periodicity. Since the GP approach doesn't show  \nFig. 9. Generalized Lomb-Scargle periodogram for TOI-3568 MAROON-X and SPIRou RVs data. The highest power is detected at 4.4178 d (dashed red vertical line) with FAP ă 10 ' 14 .  \n<!-- image -->  \nsignificant improvements compared to the window approach, we opt for the latter for the rest of our analysis for the sake of simplicity.  \nOur transit model is calculated using the code BATMAN (Kreidberg 2015), adopting a linear limb-darkening law. The quality of our ground-based WCWO photometry is not su ffi -cient to constrain the limb-darkening on its own. To simplify, we adopted a single coe ffi cient for both TESS and WCWO photometry, which is a reasonable approximation given the significant overlap between the TESS and r-band passes. The posterior distribution of transit parameters is sampled using a Bayesian Monte Carlo Markov Chain (MCMC) framework with the package emcee (Foreman-Mackey et al. 2013). We use uninformative priors for the transit parameters, shown in Table 7.", "4.2. Joint analysis of RVs and photometry to obtain the system's parameters": "We calculate the Generalized Lomb-Scargle (GLS, Zechmeister &Kürster 2009) periodogram for the MAROON-X and SPIRou RVs (see Tables 3 and 4), as illustrated in Figure 9. The highest power is detected at 4.4178 d with a false alarm probability (FAP) of FAP ă 10 ' 14 , coinciding precisely with the periodicity of transits observed in the TESS photometry data alone. This suggests that our RV data detects the signal of the star's reflex motion induced by the orbit of the planet TOI-3568 b.  \nWe therefore perform a Bayesian MCMC joint analysis of RVs and photometry data to determine the system's parameters employing the same approach as in Martioli et al. (2022, 2023). We calculate the log-likelihood for a model that includes both the transits and the RV orbit using a Keplerian model as in Martioli et al. (2010), allowing for the simultaneous fitting of the TESS photometry data within the transit windows, the groundbased single-transit photometry data, and the RV data. The RV model includes an independent systemic velocity for each RV data set. In the case of MAROON-X blue and red data, this velocity represents only a systematic o ff set with respect to the template spectrum and therefore should be close to zero. On the other hand, the SPIRou systemic RV should be close to the real radial velocity of the system, although the value obtained in our analysis is not absolutely calibrated. Initially, the parameters are fitted using an optimization least-squares (OLS) code with initial parameters obtained from an iterative preliminary analysis. We include a white noise jitter term for each RV dataset, which is fitted only in the OLS analysis. This is followed by sampling the  \nposteriors using the Bayesian MCMC framework implemented with the emcee package. We run 20,000 iterations with 50 random walkers, discarding the first 5,000 samples as burn-in. The priors and posteriors for each parameter are presented in Table 7, where the best-fit values are considered to be the mode of the distribution, and the errors are the 34th percentile on each side of the median. Using these fit parameters and the stellar parameters from Table 5, we derived other planetary quantities, which are also listed in Table 7. Note that we employ bounded priors for the time of conjunction and orbital period, which have su ffi -ciently wide bounds in the uniform distribution to ensure minimal impact on the posterior distribution. The best-fit orbit model for TOI-3568 b shows an eccentricity of 0 . 035 ˘ 0 . 021, consistent within 2 . 5 σ with a circular orbit, which is not uncommon for close-in exoplanets.  \nFigure 10 illustrates both the TESS and ground-based photometry data, and the best-fit model for the selected windows around the transits. Figures 11 and 12 illustrate the MAROONX and SPIRou RV data and the orbit model. The final RMS of photometry residuals is 3.6 ppt, while for RV residuals, it is 3.5 m s ' 1 . The MCMC samples, their correlations, and the posterior distributions for each parameter are illustrated in Figure A.1 in Appendix A.", '4.3. Limits on additional planets from TESS photometry': "We explored the detection limits of additional transiting planets by performing an injection-recovery test in the TESS light curve of TOI-3568. To do so, we employed the 2-min TESS PLD photometry residuals from sectors 55 and 56 obtained after removing the baseline GP model multiplied by the best-fit transit model for TOI-3568 b determined in Section 4.1.  \nWe used the BATMAN code (Kreidberg 2015) to generate synthetic transit signals that were injected into the TESS PLD photometry residuals. For all of these simulated planets, a limbdarkening linear law, equatorial transits ( b ' 0) and circular orbits ( e ' 0) were assumed. The value adopted for the limbdarkening coe ffi cient, u 0 ' 0 . 85, was extracted from Table 7. We surveyed the planetary radius-orbital period parameter space, RP -P , in the ranges of 0-11 R C with steps of 1 R C and 1-34 d with a 3 d step, respectively, adopting a multi-phase approach that allows five di ff erent values of Tc for each RP -P combination. As in previous works (see Jofré et al. 2021), to detect the injected signals, we ran the Transit Least Squares code (TLS; Hippke & Heller 2019), an optimized algorithm to search for periodic transits from time-series photometry. A positive planet detection was considered when the recovered orbital period is within 5% of any half multiple of the injected period. In Figure 13, we present the detectability map resulting from our injectionrecovery test.  \nFrom this plot it can be seen that planets with a radius larger than ' 4.0 R C and periods À 25 d have recovery rates of more than 80%, hence, we can exclude the presence of such additional objects in the system. For similar-sized planets but orbital periods longer than 25 d, the chances of detection are between 0 and 100%. As we consider longer periods, fewer transits of a given planet are expected. Here, the chances of detection can drop abruptly, from high to low percentages, if one or more transits are not detected (for example, if the event falls in the gap between orbits). In the space parameter corresponding to planet size smaller than ' 4.0 R C , most of the recovery rates are lower than 20%. This indicates that still might exist small planetary objects that would remain undetected in the present data.  \nTable 7. Priors and posteriors of the parameters of TOI-3568 b obtained from a joint analysis of TESS and WCWO photometry, MAROON-X and SPIRou RVs.  \nNotes. ( a ) semi-major axis derived from the fit period and Kepler's law. ( b ) assuming a uniform heat redistribution and an arbitrary geometric albedo of 0.1. ( c ) the MAROON-X RVs are relative to the template and, therefore, the systemic RVs represent a systematic o ff set relative to our orbital fit.", '5.1. Characterization of TOI-3568 b': "Considering the planet-to-star radius ratio of Rp { R ‹ ' 0 . 067 ˘ 0 . 003 and the stellar radius obtained from our stellar analysis, we derive the true physical radius of TOI-3568 b as 5 . 30 ˘ 0 . 27 R C , which is approximately 1 . 37 times the size of Neptune. The RV semi-amplitude of 12 . 1 ˘ 0 . 4 ms ' 1 implies a planet mass of 26 . 4 ˘ 1 . 0 M C , approximately 50% larger than Neptune's mass. We estimate a bulk density of 0 . 98 ˘ 0 . 15 g cm ' 3 . In Figure 14 we show the mass-radius diagram for the known exoplanets with masses in the range 0.5-500 M C and radii in the range 0.725 R C . Comparison of the evolutionary models by Fortney et al. (2007) for Hydrogen-Helium (H / He) rich planets at 0.045 au and ages between 1 and 10 Gyr reveals that TOI-3568 b is likely a H / He-dominated planet with a core of heavier elements, with a mass between 10 and 25 M C .  \nWith an orbital period of 4.4 days and at an orbital distance of 0 . 0485 ˘ 0 . 0004 au, we estimated the equilibrium temperature for TOI-3568 b as in Heng & Demory (2013), assuming a  \nuniform heat redistribution and an arbitrary geometric albedo of 0.1, which gives T eq ' 899 ˘ 12 K. Therefore this planet belongs to the hot super-Neptune class, a rare type of planet.", '5.2. TOI-3568 b in the sub-Jovian desert': 'In Figure 15, we explore the mass-period diagram, illustrating the population of exoplanets. The color contrasts in this plot represent the detection rate rather than the actual occurrence rate, which ideally should reflect the fundamental physics of planet formation. However, selection e ff ects, primarily due to the limited sensitivity of radial velocity and transit surveys, are significant only in the regime of long periods and small masses. Therefore, the scarcity of planets at short periods and masses above a few hundredths of a Jupiter mass should reflect the intrinsic planetary population in this region. We highlight the boundaries of the sub-Jovian desert as defined by Mazeh et al. (2016). TOI3568 b lies within the boundaries of the desert, albeit in a region that has a low detection rate, yet is not entirely devoid of planets. TOI-3568 b falls near the lower boundary of the subJovian desert, which is thought to be caused by photoevapora-  \nFig. 10. Photometric observations of the transits of TOI-3568 b. The top panel shows the nineteen transits observed by TESS, while the bottom panel exhibits the single transit observed by our ground-based follow-up using the 0.7 m telescope at the Wellesley College Whitin Observatory. The light blue points show the photometry data around the transits of TOI-3568 b, with the times in the top panel being relative to the central time of each transit. The red lines show the best-fit transit model and the green points show the residuals plus an arbitrary o ff set for better visualization. The dark blue and green points in the top panel represent the weighted averages of bins with a size of 0.01 d.  \n<!-- image -->  \ntion (Owen & Lai 2018). As pointed out by Frame et al. (2023), this lower boundary is becoming more blurred as we detect more planets around that region, raising the question whether this need to be reevaluated. Regardless of the existence of a real desert in this region, TOI-3568 b is notably positioned right within the transition between the populations of hot-Jupiters and superEarths, thus having potential importance in investigating the origin of this natural segmentation of planet populations.', '5.3. Evaporation status of TOI-3568 b': "Photoevaporation is thought to play an important role in the distribution of planets in the mass-period diagram (Owen & Lai 2018; Mordasini 2020; Vissapragada et al. 2022). However, this diagram does not take into account di ff erent stellar types. To assess the evaporation status of TOI-3568 b, we plot it on the energy diagram proposed by Lecavelier Des Etangs (2007). This diagram explores the relationship between the potential energy  \nof the planet and the extreme ultraviolet (EUV) luminosity it receives from its parent star. These two competing sources of energy control the evaporation rate of the planet throughout its lifetime.  \nAs in Eq. 10 of Lecavelier Des Etangs (2007), we consider the planet's potential energy Ep 1 , including tidal forces from the gravitational interaction with the parent star, and calculate the EUV luminosity following his recipe. For simplicity, we assume a constant mean EUV luminosity and adopt a correction factor of γ ' 6 to account for the time variation in the energy flux received by the planet. This approximation becomes more problematic for very young systems and very hot stars, but both types of objects account for a small fraction of the exoplanets that we present in our analysis.  \nFigure 16 reproduces the energy plot from Lecavelier Des Etangs (2007), showing up-to-date exoplanet data from the exoplanet.eu catalog for transiting planets with measured masses and highlighting planets in three mass regimes: Jupiterlike (M p ą 2MNep), Neptune-like (0 . 25MNep ă M p ă 2MNep), and Earth-like (M p ă 0 . 25MNep). TOI-3568 b stands out as one of the super-Neptunes (M p ' 1 . 54 ˘ 0 . 06 MNep) with the highest levels of EUV luminosity, receiving about dE EUV { dt ą 10 40 . 6 erg Gyr ' 1 . The dashed lines in Figure 16 represent planet lifetimes of 0.1, 5, and 10 Gyr. The regions below these lines are considered evaporation-forbidden regions, where planets receive more EUV energy than is needed to fill their potential well, and thus would evaporate in less than 0.1, 5, or 10 Gyr, respectively.  \nAs in the mass-period diagram, this energy diagram also reveals two distinct populations of exoplanets: one population consists of more massive planets with high potential energy and high EUVluminosity. These planets accumulate in the top-left part of the diagram, above the forbidden region at 5 Gyr, with the lower boundary clearly sculpted by photoevaporation. A second population of less massive planets (M p ă 2MNep) is accumulated in a lower EUV luminosity regime in the bottom-right part of the diagram. TOI-3568 b lies at the transition between these two populations, making it di ffi cult to clearly classify it into one category or the other. Although being in a region that is not significantly populated, TOI-3568 b is above a lifetime of 5 Gyr, which is consistent with the mature age of the system. However, the question remains: why is this region less populated? As pointed out by Owen & Lai (2018), there are likely other formation mechanisms sculpting the lower mass end, causing this natural gap in the population. TOI-3568 b appears to be an important planet for probing this gap. However, a more detailed analysis of this topic is beyond the scope of this paper.", '6. Conclusions': "We report the discovery of the transiting exoplanet TOI-3568 b, an inflated hot super-Neptune situated in the sub-Jovian desert. Using the TESS and ground-based photometry, MAROON-X optical spectra, and SPIRou NIR spectropolarimetry, we determined the orbit of the planet and the physical parameters of the system. The star is a quiet and mature K dwarf with an e ff ective temperature of 4969 ˘ 45 K, with nearly solar metallicity. Our analysis identifies this star as part of the transitional population between the galactic thin and thick disk, exhibiting characteristics more consistent with the thin disk population. Although our observations did not confirm this candidate as a member of the galactic thick disk, it emerged as an interesting discovery of a rare super-Neptune situated in a region of the mass-period diagram with a low occurrence rate of planets. TOI-3568 b is likely an H / He-dominated planet with a core of heavier elements  \nFig. 11. TOI-3568 MAROON-X and SPIRou RVs. In the top panel, the red and blue points show the MAROON-X RV data for the red and blue channels, and the orange triangles show the SPIRou RV data. The green line shows the best-fit orbit model for TOI-3568 b. The bottom panel shows the residuals.  \n<!-- image -->  \nFig. 12. TOI-3568 RVs in phase with the orbit of the planet. In the top panel, the red and blue points show the MAROON-X RV data for the red and blue channels, and the orange points show the SPIRou RV data. The green line shows the best-fit orbit model for TOI-3568 b. The bottom panel shows the residuals.  \n<!-- image -->  \nwith a mass between 10 and 25 M C , and it lies at the lowest point between the two significant populations of hot-Jupiters and super-Earths. We analyzed the photoevaporation status of TOI3568 b, finding that this planet experiences a high regime of EUV luminosity for its mass range, making it one of the planets with the highest EUV luminosities among those with a mass M p ă 2 MNep. However, this planet is not in an evaporationforbidden region, as its status is still consistent with a planet having an evaporation lifetime exceeding 5 Gyr. Perhaps the most interesting aspect of this planet is that it is not a common type of planet. It lies in a transition region in both the mass-period and energy diagrams, an area with a dearth of planets that cannot be explained by photoevaporation.  \nAcknowledgements. This work is based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the US National Science Foundation (NSF), the Cana-  \nFig. 13. Injection-recovery test performed on the TESS PLD light curve residuals from sectors 55 and 56 of TOI-3568 to check the detectability of additional planets in the system. High and low recovery rates are presented in light yellow-green and dark blue colors, respectively. From this detectability map, we can conclude that the existence of additional planets with radii Á 4.0 R C and periods À 25 d should be excluded, given that the recovery rates range from 80 to 100%. However, small planets with sizes À 4.0 R C would remain undetected for almost the entire set of periods investigated.  \n<!-- image -->  \ndian National Research Council (NRC), the Chilean Agencia Nacional de Investigación y Desarrollo (ANID), the Brazilian Ministério da Ciência, Tecnologia e Inovação, the Argentinean Ministerio de Ciencia, Tecnología e Innovación, and the Korea Astronomy and Space Institute (KASI). This work is based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated from the summit of Maunakea by the National Research Council of Canada, the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii. Based on observations obtained with SPIRou, an international project led by Institut de Recherche en Astrophysique et Planétologie, Toulouse, France. E.M. acknowledges funding from Fundação de Amparo à Pesquisa do Estado de Mi-  \nFig. 14. Mass-radius diagram. The grey circles show the exoplanet data from exoplanets.eu and the red stars show the solar system planets. The black point corresponds to the measured value for TOI-3568 b. Internal structure models for planets with pure rock and pure ice (H2O) compositions from Zeng et al. (2019) are depicted by the solid red and orange lines, respectively. Additionally, evolutionary models for H / He dominated planets from Fortney et al. (2007) are shown for a pure H / He without a core (green line), and for a core of 10 M C (blue line) and 25 M C (purple line) of heavier elements. Iso-density curves are also represented by the black lines.  \n<!-- image -->  \nFig. 15. Mass-period diagram. The grey circles show the exoplanet data from exoplanet.eu , where the size of the circles is proportional to the planetary radius. The color map depicts the planet detection rate in percentage normalized to the sum of all known exoplanets in the mass range between 0.003 MJup and 20 MJup and orbital periods between 0.3 d and 200 d. The red dashed lines show the innermost limits of the sub-Jovian desert as defined by Mazeh et al. (2016). The orange star shows TOI-3568 b, which lies within the desert right in the transition between the two populations of the most frequent types of exoplanets, hot-Jupiters (top) and warm super-Earths (bottom).  \n<!-- image -->  \n(FAPEMIG) under project number APQ-02493-22 and research productivity grant (PQ) number 309829 / 2022-4 awarded by the National Council for Scientific and Technological Development (CNPq), Brazil. The University of Chicago group acknowledges funding for the MAROON-X project from the David and Lucile Packard Foundation, the Heising-Simons Foundation, the Gordon and Betty Moore Foundation, the Gemini Observatory, the NSF (award  \nFig. 16. Diagram of planet potential energy versus mean EUV luminosity for transiting exoplanets from the exoplanet.eu catalog with measured masses. Exoplanets are represented by di ff erent symbols for three mass ranges: M p ą 2MNep (blue hollow circles), 0 . 25MNep ă M p ă 2MNep (blue filled circles), and M p ă 0 . 25MNep (orange circles). TOI3568 b is highlighted by the red circle. Dashed lines represent planet lifetimes of 0.1 Gyr (blue), 5 Gyr (orange), and 10 Gyr (green).  \n<!-- image -->  \nnumber 2108465), and NASA (grant number 80NSSC22K0117). The Gemini observations are associated with programs ID GN-2022A-Q-207 / -Q-113. Funding for the TESS mission is provided by NASA's Science Mission Directorate. We acknowledge the use of public TESS data from pipelines at the TESS Science O ffi ce and at the TESS Science Processing Operations Center. Resources supporting this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center for the production of the SPOC data products. DR was supported by NASA under award number NNA16BD14C for NASA Academic Mission Services. TESS data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST) at the Space Telescope Science Institute. This work has made use of data from the European Space Agency (ESA) mission Gaia ( https://www.cosmos.esa.int/gaia ), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www. cosmos.esa.int/web/gaia/dpac/consortium ). Funding for the DPAC has been provided by national institutions, in particular, the institutions participating in the Gaia Multilateral Agreement. This work made use of tpfplotter by J. Lillo-Box (publicly available in www.github.com / jlillo / tpfplotter), which also made use of the python packages astropy , lightkurve , matplotlib and numpy . KKMacknowledges support from the New York Community Trust Fund for Astrophysical Research.", 'References': "Adibekyan, V. Z., Santos, N. C., Sousa, S. G., & Israelian, G. 2011, A&A, 535, L11  \nAdibekyan, V. Z., Sousa, S. G., Santos, N. C., et al. 2012, A&A, 545, A32 Aller, A., Lillo-Box, J., Jones, D., Miranda, L. F., & Barceló Forteza, S. 2020, A&A, 635, A128  \n- Amarante, J. A. S., Smith, M. C., & Boeche, C. 2020, MNRAS, 492, 3816\n- Amarsi, A. M., Nissen, P. E., & Skúladóttir, Á. 2019, A&A, 630, A104\n- Arenou, F., Grenon, M., & Gomez, A. 1992, A&A, 258, 104\n- Artigau, É., Cadieux, C., Cook, N. J., et al. 2022, AJ, 164, 84\n- Bensby, T., Feltzing, S., & Lundström, I. 2003, A&A, 410, 527\n- Bensby, T., Feltzing, S., Lundström, I., & Ilyin, I. 2005, A&A, 433, 185  \nBernkopf, J., Fidler, A., & Fuhrmann, K. 2001, in Astronomical Society of the  \nPacific Conference Series, Vol. 245, Astrophysical Ages and Times Scales,  \ned. T. von Hippel, C. Simpson, & N. Manset, 207  \nBlanco-Cuaresma, S. 2019, MNRAS, 486, 2075  \nBouchy, F., Hebb, L., Skillen, I., et al. 2010, A&A, 519, A98  \nBressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127  \n- Buchhave, L. A., Bizzarro, M., Latham, D. W., et al. 2014, Nature, 509, 593\n- Campante, T. L., Barclay, T., Swift, J. J., et al. 2015, ApJ, 799, 170  \nCarlberg, J. K., Cunha, K., Smith, V. V., & Majewski, S. R. 2012, ApJ, 757, 109  \n- Carrillo, A., Hawkins, K., Bowler, B. P., Cochran, W., & Vanderburg, A. 2020, MNRAS, 491, 4365\n- Collins, K. 2019, in American Astronomical Society Meeting Abstracts, Vol. 233, American Astronomical Society Meeting Abstracts #233, 140.05 Collins, K. A., Kielkopf, J. F., Stassun, K. G., & Hessman, F. V. 2017, AJ, 153, 77\n- Cook, N. J., Artigau, É., Doyon, R., et al. 2022, PASP, 134, 114509 Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data\n- Catalog, II / 246\n- da Silva, L., Girardi, L., Pasquini, L., et al. 2006, A&A, 458, 609\n- Davenport, J. R. A. 2016, ApJ, 829, 23\n- Demarque, P., Woo, J.-H., Kim, Y.-C., & Yi, S. K. 2004, ApJS, 155, 667\n- Deming, D., Knutson, H., Kammer, J., et al. 2015, ApJ, 805, 132\n- Donati, J. F., Kouach, D., Moutou, C., et al. 2020, MNRAS, 498, 5684\n- Donati, J. F., Semel, M., Carter, B. D., Rees, D. E., & Collier Cameron, A. 1997, MNRAS, 291, 658\n- Eigmüller, P., Chaushev, A., Gillen, E., et al. 2019, A&A, 625, A142\n- Eigmüller, P., Gandolfi, D., Persson, C. M., et al. 2017, AJ, 153, 130\n- Fischer, D. A. & Valenti, J. 2005, ApJ, 622, 1102\n- Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306\n- Fortney, J. J., Marley, M. S., & Barnes, J. W. 2007, ApJ, 659, 1661  \nFrame, G., Armstrong, D. J., Cegla, H. M., et al. 2023, MNRAS, 523, 1163  \nFuhrmann, K. 1998, A&A, 338, 161  \n- Gaia Collaboration. 2020, VizieR Online Data Catalog, I / 350  \nGaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1  \n- Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1  \nGhezzi, L., Cunha, K., Smith, V. V., & de la Reza, R. 2010, ApJ, 724, 154  \nGhezzi, L., Martinez, C. F., Wilson, R. F., et al. 2021, ApJ, 920, 19  \nGhezzi, L., Montet, B. T., & Johnson, J. A. 2018, ApJ, 860, 109  \nGilliland, R. L., Chaplin, W. J., Dunham, E. W., et al. 2011, ApJS, 197, 6  \nGuerrero, N. M., Seager, S., Huang, C. X., et al. 2021, ApJS, 254, 39  \nHaywood, M. 2008, MNRAS, 388, 1175  \n- Heng, K. & Demory, B.-O. 2013, ApJ, 777, 100\n- Hippke, M. & Heller, R. 2019, Astronomy and Astrophysics, 623, A39\n- Howell, S. B., Everett, M. E., Sherry, W., Horch, E., & Ciardi, D. R. 2011, AJ, 142, 19\n- Ilin, E., Schmidt, S. J., Poppenhäger, K., et al. 2021, A&A, 645, A42  \nJenkins, J. M., Twicken, J. D., McCauli  \nff  \n, S., et al. 2016, in Proc. SPIE, Vol.  \n- 9913, Software and Cyberinfrastructure for Astronomy IV, 99133E\n- Jensen, E. 2013, Tapir: A web interface for transit / eclipse observability, Astrophysics Source Code Library, record ascl:1306.007\n- Jofré, E., Almenara, J. M., Petrucci, R., et al. 2020, A&A, 634, A29  \nJofré, E., Petrucci, R., Maqueo Chew, Y. G., et al. 2021, AJ, 162, 291  \n- Johnson, J. A., Aller, K. M., Howard, A. W., & Crepp, J. R. 2010, PASP, 122, 905  \nKálmán, S., Szabó, G. M., Borsato, L., et al. 2023, MNRAS, 522, 488  \nKatz, D., Sartoretti, P., Guerrier, A., et al. 2023, A&A, 674, A5  \nKilic, M., Munn, J. A., Harris, H. C., et al. 2017, ApJ, 837, 162  \nKreidberg, L. 2015, PASP, 127, 1161  \nKunimoto, M., Daylan, T., Guerrero, N., et al. 2022, ApJS, 259, 33  \nLecavelier Des Etangs, A. 2007, A&A, 461, 1185  \n- Lightkurve Collaboration, Cardoso, J. V. d. M., Hedges, C., et al. 2018, Lightkurve: Kepler and TESS time series analysis in Python, Astrophysics Source Code Library\n- Luger, R., Agol, E., Kruse, E., et al. 2016, AJ, 152, 100\n- Luger, R., Kruse, E., Foreman-Mackey, D., Agol, E., & Saunders, N. 2018, AJ, 156, 99\n- Martioli, E., Hébrard, G., Correia, A. C. M., Laskar, J., & Lecavelier des Etangs, A. 2021, A&A, 649, A177\n- Martioli, E., Hébrard, G., de Almeida, L., et al. 2023, arXiv e-prints, arXiv:2311.07011\n- Martioli, E., Hébrard, G., Fouqué, P., et al. 2022, A&A, 660, A86\n- Martioli, E., Hébrard, G., Moutou, C., et al. 2020, A&A, 641, L1\n- Martioli, E., McArthur, B. E., Benedict, G. F., et al. 2010, ApJ, 708, 625  \nMazeh, T., Holczer, T., & Faigler, S. 2016, A&A, 589, A75  \n- McQuillan, A., Aigrain, S., & Mazeh, T. 2013, MNRAS, 432, 1203\n- Miglio, A., Chiappini, C., Morel, T., et al. 2013, MNRAS, 429, 423 Mordasini, C. 2020, A&A, 638, A52\n- Mori, M., Livingston, J. H., Leon, J. d., et al. 2022, AJ, 163, 298\n- Mortier, A., Zapatero Osorio, M. R., Malavolta, L., et al. 2020, MNRAS, 499, 5004  \nMurgas, F., Astudillo-Defru, N., Bonfils, X., et al. 2021, A&A, 653, A60  \nNissen, P. E. 2004, in Origin and Evolution of the Elements, ed. A. McWilliam  \n&M. Rauch, 154  \n- Owen, J. E. & Lai, D. 2018, MNRAS, 479, 5012  \nPerottoni, H. D., Amarante, J. A. S., Limberg, G., et al. 2021, ApJ, 913, L3  \n- Petrucci, R. P., Gómez Maqueo Chew, Y., Jofré, E., Segura, A., & Ferrero, L. V. 2024, MNRAS, 527, 8290\n- Ramírez, I., Allende Prieto, C., & Lambert, D. L. 2007, A&A, 465, 271\n- Ramírez, I., Fish, J. R., Lambert, D. L., & Allende Prieto, C. 2012, ApJ, 756, 46 Ramírez, I., Meléndez, J., Bean, J., et al. 2014, A&A, 572, A48\n- Reddy, B. E., Lambert, D. L., & Allende Prieto, C. 2006, MNRAS, 367, 1329 Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, Journal of Astronomical Telescopes, Instruments, and Systems, 1, 014003\n- Scargle, J. D. 1982, ApJ, 263, 835\n- Schuler, S. C., Vaz, Z. A., Katime Santrich, O. J., et al. 2015, ApJ, 815, 5 Seifahrt, A., Bean, J. L., Stürmer, J., et al. 2020, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 11447, Groundbased and Airborne Instrumentation for Astronomy VIII, ed. C. J. Evans, J. J. Bryant, & K. Motohara, 114471F\n- Seifahrt, A., Stürmer, J., Bean, J. L., & Schwab, C. 2018, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 10702, Ground-based and Airborne Instrumentation for Astronomy VII, ed. C. J. Evans, L. Simard, & H. Takami, 107026D\n- Sharma, S., Bland-Hawthorn, J., Johnston, K. V., & Binney, J. 2011, ApJ, 730, 3 Smith, J. C., Stumpe, M. C., Van Cleve, J. E., et al. 2012, PASP, 124, 1000 Sneden, C. A. 1973, PhD thesis, University of Texas, Austin\n- Sousa, S. G., Santos, N. C., Adibekyan, V., Delgado-Mena, E., & Israelian, G. 2015, A&A, 577, A67\n- Stassun, K. G., Oelkers, R. J., Paegert, M., et al. 2019, AJ, 158, 138\n- Stassun, K. G., Oelkers, R. J., Pepper, J., et al. 2018, AJ, 156, 102\n- Stumpe, M. C., Smith, J. C., Catanzarite, J. H., et al. 2014, PASP, 126, 100  \nStumpe, M. C., Smith, J. C., Van Cleve, J. E., et al. 2012, PASP, 124, 985  \n- Stürmer, J., Seifahrt, A., Schwab, C., & Bean, J. L. 2017, Journal of Astronomical Telescopes, Instruments, and Systems, 3, 025003\n- Szabó, G. M., Kálmán, S., Borsato, L., et al. 2023, A&A, 671, A132\n- Szabó, G. M. & Kiss, L. L. 2011, ApJ, 727, L44\n- Takeda, Y. & UeNo, S. 2017, PASJ, 69, 46\n- Teske, J. K., Ghezzi, L., Cunha, K., et al. 2015, ApJ, 801, L10  \nTucci Maia, M., Ramírez, I., Meléndez, J., et al. 2016, A&A, 590, A32 Twicken, J. D., Catanzarite, J. H., Clarke, B. D., et al. 2018, PASP, 130, 064502 Van Cleve, J. E., Howell, S. B., Smith, J. C., et al. 2016, PASP, 128, 075002 Vissapragada, S., Knutson, H. A., Greklek-McKeon, M., et al. 2022, AJ, 164, 234  \nWang, J. & Fischer, D. A. 2015, AJ, 149, 14 Weiss, L. M., Dai, F., Huber, D., et al. 2021, AJ, 161, 56 West, R. G., Gillen, E., Bayliss, D., et al. 2019, MNRAS, 486, 5094 Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 Yi, S., Demarque, P., Kim, Y.-C., et al. 2001, ApJS, 136, 417 Zechmeister, M. & Kürster, M. 2009, A&A, 496, 577  \n- Zechmeister, M., Reiners, A., Amado, P. J., et al. 2018, A&A, 609, A12\n- Zeng, L., Jacobsen, S. B., Sasselov, D. D., et al. 2019, Proceedings of the National Academy of Science, 116, 9723\n- 1 Laboratório Nacional de Astrofísica, Rua Estados Unidos 154, 37504-364, Itajubá - MG, Brazil, e-mail: martioli@iap.fr\n- 2 Institut d'Astrophysique de Paris, CNRS, UMR 7095, Sorbonne Université, 98 bis bd Arago, 75014 Paris, France\n- 3 Universidad Nacional de Córdoba - Observatorio Astronómico de Córdoba, Laprida 854, X5000BGR, Córdoba, Argentina\n- 4 Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Godoy Cruz 2290, CABA, CPC 1425FQB, Argentina\n- 5 Observatoire de Haute Provence, St Michel l'Observatoire, France\n- 6 Universidade Federal do Rio de Janeiro, Observatório do Valongo, Ladeira do Pedro Antônio, 43, Rio de Janeiro, RJ 20080-090, Brazil\n- 7 Instituto de Astronomía, Universidad Nacional Autónoma de México, Ciudad Universitaria, Ciudad de México, 04510, México\n- 8 International Center for Advanced Studies (ICAS) and ICIFI (CONICET), ECyT-UNSAM, Campus Miguelete, 25 de Mayo y Francia, (1650) Buenos Aires, Argentina.\n- 9 NASA Ames Research Center, Mo ff ett Field, CA 94035, USA\n- 10 Research Institute for Advanced Computer Science, Universities Space Research Association, Washington, DC 20024, USA\n- 11 Canada-France-Hawaii Telescope, CNRS, 96743 Kamuela, Hawaii, USA\n- 12 Université de Montréal, Département de Physique, IREX, Montréal, QC H3C 3J7, Canada\n- 13 Department of Astronomy & Astrophysics, University of Chicago, Chicago, IL, USA\n- 14 Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA  \n- 15 Aix Marseille Univ, CNRS, CNES, LAM, 38 rue Frédéric JoliotCurie, 13388 Marseille, France\n- 16 Université Grenoble Alpes, CNRS, IPAG, 414 rue de la Piscine, 38400 St-Martin d'Hères, France\n- 17 Université de Toulouse, CNRS, IRAP, 14 avenue Belin, 31400 Toulouse, France\n- 18 NASA Exoplanet Science Institute, Caltech IPAC, 1200 E. California Blvd., Pasadena, CA 91125, USA\n- 19 LESIA, Observatoire de Paris, Université PSL, CNRS, Sorbonne Université, Université Paris Cité, 5 place Jules Janssen, 92195 Meudon, France\n- 20 Department of Physics and Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, Cambridge, MA02139, USA\n- 21 Department of Astronomy, Wellesley College, Wellesley, MA 02481, USA\n- 22 Universidade de São Paulo, Instituto de Astronomia, Geofísica e Ciências Atmosféricas (IAG), Departamento de Astronomia, Rua do Matão 1226, Cidade Universitária, 05508-900, SP, Brazil\n- 23 NASA Goddard Space Flight Center, 8800 Greenbelt Road, Greenbelt, MD 20771, USA\n- 24 Royal Astronomical Society, Burlington House, Piccadilly, London W1J 0BQ, United Kingdom\n- 25 Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands\n- 26 Landessternwarte, Zentrum für Astronomie der Universität Heidelberg, Königstuhl 12, D-69117 Heidelberg, Germany\n- 27 SETI Institute, Mountain View, CA 94043, USA", 'Appendix A: Posterior distributions of model parameters': 'This appendix presents, in Figure A.1, the MCMC samples and the final posterior distributions of the free parameters used in our joint analysis of the TESS and ground-based photometry, as well as the MAROON-X and SPIRou RVs of TOI-3568.  \nFig. A.1. Pairs plot showing the MCMC samples and posterior distributions of the free parameters in our joint analysis of the TESS and groundbased photometry and the MAROON-X and SPIRou RVs of TOI-3568. The contours mark the 1 σ , 2 σ , and 3 σ regions of the distribution. The blue crosses indicate the best-fit values for each parameter obtained by the mode, and the dashed vertical lines in the projected distributions show the median values and the 1 σ uncertainty (34% on each side of the median).  \n<!-- image -->'}