Detection of an Atmospheric Outflow from the Young Hot Saturn TOI-1268b

Photoevaporative mass-loss rates are expected to be highest when planets are young and the host star is more active, but to date there have been relatively few measurements of mass-loss rates for young gas giant exoplanets. In this study we measure the present-day atmospheric mass-loss rate of TOI-1268b, a young (110–380 Myr) and low density (0.71 −0.13+0.17 g cm−3) hot Saturn located near the upper edge of the Neptune desert. We use Palomar/WIRC to search for excess absorption in the 1083 nm helium triplet during two transits of TOI-1268b. We find that it has a larger transit depth ( 0.285−0.050+0.048% excess) in the helium bandpass than in the TESS observations, and convert this excess absorption into a mass-loss rate by modeling the outflow as a Parker wind. Our results indicate that this planet is losing mass at a rate of logṀ=10.2±0.3 g s−1 and has a thermosphere temperature of 6900−1200+1800 K. This corresponds to a predicted atmospheric lifetime much larger than 10 Gyr. Our result suggests that photoevaporation is weak in gas giant exoplanets even at early ages.


INTRODUCTION
Planets on close-in orbits receive extreme amounts of high-energy radiation from their host star, which can cause their atmospheres to undergo hydrodynamic escape (Owen 2019).Young stars typically have high activity levels (Johnstone et al. 2021;King & Wheatley 2021), and young planets that are still radiating away residual heat from their formation have low surface gravities.As a result, most atmospheric mass loss is thought to take place at relatively early times (< Gyr; e.g., Kubyshkina & Fossati 2022;Ketzer & Poppenhaeger 2023).However, radiative transfer modeling of photoevaporative outflows also predicts that the stellar highenergy flux does work powering the outflow, thus modifying the mass loss efficiency (e.g., Murray-Clay et al. 2009;Salz et al. 2016;Caldiroli et al. 2022).For young planets on close-in orbits with higher XUV fluxes the corresponding mass loss efficiency will be reduced (see equation 1).
Where Ṁ is the mass loss rate, ε is the mass loss efficiency, R p is the planet's radius, F XUV is the incident XUV flux from the host star, G is the gravitational constant and M p is the planet's mass.Young stars also have fast and dense winds, which can strongly sculpt the planetary outflow geometry and may alter the corresponding mass-loss rates (e.g., Kubyshkina & Fossati 2021;Kubyshkina & Vidotto 2021;MacLeod & Oklopčić 2022;Wang & Dai 2021;Carolan et al. 2021).In order to understand the net effect of these competing processes, it is therefore critical to measure the mass-loss rates of planets orbiting young stars.
There are currently very few observational constraints on mass-loss rates for young transiting exoplanets (e.g., Dos Santos 2022).We can quantify the atmospheric mass-loss rates of transiting planets by measuring the excess absorption in either the Lyman α or the 1083 nm helium absorption triplet (Dos Santos 2022; Owen 2019).Both of these features are strongly absorbing at low pressures, resulting in a significantly deeper transit when there is an atmospheric outflow (for a more detailed explanation see e.g.Vidal-Madjar et al. 2003;Salz et al. 2015;Spake et al. 2018;Oklopčić & Hirata 2018).In a series of recent papers, Zhang et al. (2022bZhang et al. ( ,a, 2023) used a combination of Lyman α and metastable helium absorption to measure the mass-loss rates for a sample of five sub-Neptune-sized (2−3 R ⊕ ) planets orbiting young stars.However, three of these planets do not have measured masses, while the remaining two (HD 63433c and TOI-560b) have relative uncertainties greater than 20% on their measured masses (Barragán et al. 2021;Mallorquín et al. 2023).This complicates the interpretation of their atmospheric absorption signatures during the transit.Stellar activity poses an additional challenge when measuring mass-loss rates for sub-Neptune-sized planets orbiting young stars, as the stellar Lyα emission and 1083 nm metastable helium absorption features can vary by amounts comparable to the white-light transit depth (Zhang et al. 2022b(Zhang et al. , 2023)).Young transiting gas giant planets provide a promising alternative pathway to study mass loss in young planetary systems.Their masses are easier to measure than those of young transiting sub-Neptunes, while their larger planet-to-star radius ratios allow for more precise constraints on their atmospheric mass-loss rates.Of the sample of gas giant exoplanets with published mass-loss rates, only two (WASP-52b and WASP-80b) have -ages less than a Gyr (400  Triaud et al. 2013 respectively).While these estimates should be interpreted with caution as angular momentum transfer between close-in Jupiters and their host stars can obfuscate gyrochronology (Lanza 2010;Poppenhaeger & Wolk 2014;Mancini et al. 2017), intriguingly, both of these planets appear to have weaker-than-predicted outflows (Vissapragada et al. 2022;Kirk et al. 2022;Fossati et al. 2022).
In this work, we utilize the 1083 nm metastable helium absorption triplet (Spake et al. 2018;Nortmann et al. 2018) to measure the present-day mass-loss rate of the Jupiter-sized (0.81 ± 0.05R J ; Šubjak et al. 2022), Saturn-mass (0.303±0.026MJ ; Šubjak et al. 2022) planet TOI-1268b.This is the youngest (110 − 380 Myr) hot Saturn-mass transiting exoplanet currently known, and one of the only young transiting planets with a well-measured mass.It orbits an early K star with a mass of 0.90 ± 0.13M ⊙ and a radius of 0.86 ± 0.02R ⊙ (Dong et al. 2022).With a planetary density of just 0.71 +0.17 −0.13 g cm −3 ( Šubjak et al. 2022) and a J-band magnitude of 9.40 ± 0.02 (Skrutskie et al. 2006), TOI-1268b is an ideal target for helium mass loss studies.Rossiter-McLaughlin measurements also indicate that the planet's orbit is aligned with its host star (Dong et al. 2022), which along with the system's relatively young age appears to disfavor at least some forms of high-eccentricity migration (Albrecht et al. 2022).
This planet is also notable because it resides at the edge of the 'Neptune desert' (see Figure 1), which is a deficit of Neptune-sized planets on close-in orbits (Szabó & Kiss 2011;Beauge & Nesvornỳ 2012;Lundkvist et al. 2016).The mechanisms that create this desert are linked to the origins of these planets, which may have been formed in situ or migrated inward to their current orbits (Bailey & Batygin 2018;Fortney et al. 2021).Several studies have suggested that the upper edge of the desert might be explained by high eccentricity migration followed by tidal disruption of planets with periastron distances that bring them too close to the star (Owen & Lai 2018;Matsakos & Königl 2016).Other studies have suggested that photoevaporation may clear out the desert (Kurokawa & Nakamoto 2014;Thorngren et al. 2023), but we previously surveyed photoevaporation rates for planets near the upper edge of the desert (see Figure 1), and found that atmospheric mass loss for planets in this region are too low to sculpt this gap (Vissapragada et al. 2020(Vissapragada et al. , 2022;;Paragas et al. 2021).One limitation of this study was the relatively old ages of the planets in our sample, which forced us to extrapolate back in time in order to calculate the cumulative mass-loss rates for the planets in our sample.TOI-1268b is located near the upper edge of the desert, and therefore provides us with an opportunity to circumvent this limitation by probing mass loss directly and quantitatively during the period when it may be most important for planetary evolution.
In this study, we characterize the mass-loss rate of TOI-1268b for the first time.We observed one full transit and one partial transit with the the Wide-Field Infrared Camera (WIRC) on the Hale 200" Telescope on Palomar Observatory using a narrow-band filter centered on the 1083 nm helium line.In §2 we describe our observations.In §3 we perform a joint fit with the phased TESS light curve to better constrain the transit shape and provide a baseline white-light transit depth for comparison.In §4, we discuss the implications of our measurement for this planet's mass loss history.

Palomar/WIRC
We observed two transits of TOI-1268b on the nights of UT 11 March 2022 and 29 April 2022.During the first night conditions were mostly good, with thin clouds that modestly reduced the total number of counts during the second half of the night.At the start of the night the airmass value was 2.41, at the end it was 1.37 and the target passed through the zenith.Conditions during the second night were more strongly affected by clouds, which prevented us from taking data during the transit ingress.At the start of the night the airmass value was 1.32, at the end it was 1.87 and the target passed through the zenith.We additionally discarded the first nine images from the second night, as the cloud cover during these exposures was thick enough that the telescope tracking algorithm struggled to calculate the centroid of the the point spread functions (PSFs).
We utilized a beam-shaping diffuser to reshape the PSFs of stars into a 3 ′′ top hat, which helps control timecorrelated systematics (Stefansson et al. 2017;Vissapragada et al. 2020).Following the procedure described in Vissapragada et al. (2020) and Paragas et al. (2021), we first performed the standard astronomical image calibration steps to create dark-subtracted, flat-fielded images.We then stacked a set of nine calibrated dithered images to create a combined background frame, which we used to remove the OH telluric emission lines from each image; these lines form bright radial arcs on the detector as a result of our narrow bandpass.We extracted photometric light curves for our target and a set of comparison stars using aperture photometry with radii ranging from 1 to 25 pixels (0. ′′ 25 to 6. ′′ 25) using the photutils package (Bradley et al. 2022).We selected the optimal aperture for each night by choosing the aperture which minimized the rms in our fit residuals (0.39% and 0.45% for the first and second night, respectively; Vissapragada et al. 2020).We found optimal aperture sizes of 12 and 14 pixels (3 ′′ and 3. ′′ 5) on the first and second nights, respectively and used the same 3 comparison stars for each night.
To choose the comparison stars, we considered all detected sources which had signal-to-noise ratio values larger than 50 in the first image.We then selected three stars whose extracted photometry most closely matched the target star.

Transiting Exoplanet Survey Satellite
The Transiting Exoplanet Survey Satellite (TESS ) is a space observatory that was designed to carry out a dedicated all-sky survey for transiting extrasolar planets (Ricker et al. 2015).We downloaded data for TOI-1268 from TESS sectors 15, 21, 22, 41, 48 and 49 using the Mikulski Archive for Space Telescopes (MAST).We used the Presearch Data Conditioning Simple Aperture Photometry (PDCSAP) fluxes obtained from MAST.These fluxes were extracted using the Science Processing Operations Center (SPOC) pipeline (Jenkins et al. 2016).We then removed outliers using a 5σ clip of the data using a moving median filter of a width of 31 data points for 5 iterations.
As discussed by Šubjak et al. (2022) and Dong et al. (2022), this relatively young star exhibits a large rotational modulation (P = 10.9 ± 0.5 days, 1 − 2% variability) in the TESS data due to star spots.In order to fit the transit shape, we must first remove this variability from the light curve.We calculated mid transit times for all transits using the linear ephemeris from Šubjak et al. (2022).We masked out all of the transits using a 1.04 transit duration (T 14 = 4.001 +0.025  −0.024 hr) wide mask (± 0.52 transit duration from the transit center) and then divided the data into segments, where we considered any gap longer than 1 day as a division between segments.We mirrored the data at the beginning and end of each segment to eliminate edge effects on the smoothing filter (Greklek-McKeon et al. 2023), and used a Gaussian kernel convolution with a kernel width of 1 transit duration to generate a smoothed trend representing the rotational modulation.We divided this trend out of our TESS data to normalize and flatten the light curve.We also found that transit numbers 13 and 15 in the TESS photometry had star spot crossings, which we identified by first fitting the stacked profile of all transits and then flagging transits with 5 or more consecutive points that are more than 3σ from the best-fit stacked transit model.We removed these transits from the stacked profile so as to not bias the joint fit.et al. 2021).Our model included the quadratic limb darkening coefficients (u 1 , u 2 ), impact parameter b, orbital period P , mid-transit time T 0 , transit depth R p /R star , stellar radius R ⋆ and semi-major axis a/R star as free parameters.We also included a set of additional free parameters to account for time-varying systematic noise sources, where we selected the optimal set of de-trending model parameters for each night of WIRC data by minimizing the Bayesian Information Criterion (BIC, Schwarz 1978).Our final set of detrending parameters for both WIRC nights included a set of linear weights for the three comparison star light curves and the airmass of the target star.During the second night, we also added a linear function of the water proxy described in Paragas et al. (2021), which tracks the time-varying telluric water absorption.We considered models with linear weights on the PSF width and the target centroid position, but found that they did not improve the quality of the fit as tracked by the BIC.
To inform our priors for the joint WIRC+TESS fit, we first fit the TESS light curve independently, using a model that includes the quadratic limb darkening coefficients (u 1 , u 2 ), impact parameter b, orbital period P , mid-transit time T 0 , transit depth R p /R star , stellar radius R ⋆ and semi major axis a/R star as free parameters.We list the priors and posteriors for each fit parameter in Table 1.
We then fit the light curves for each WIRC night individually using the NUTS sampler in PyMC3 (Salvatier et al. 2016).We ran 4 chains with 1500 tuning and 1000 sampling steps per chain and we found the Gelman-Rubin parameter R < 1.01 for all sampled parameters, indicating convergence.We obtained a transit depth (R p /R star ) 2 of 0.98 +0.05 −0.04 % for the first night and 1.0 +0.2 −0.3 % for the second night.These values are consistent within 1σ, indicating that the amount of helium absorption over a timescale of two months is constant at the level of our measurement errors.
Finally, we fit the two WIRC light curves simultaneously along with the phased TESS light curve.We fit for global values of the impact parameter b, the orbital period P , the predicted mid-transit time T 0 , the stellar radius R ⋆ and the semi-major axis a/R star parameters.We allowed the transit depths R p /R star and the quadratic limb darkening coefficients (u 1 , u 2 ) to vary independently in each bandpass.We included a TESS error scaling parameter describing the deviation from photon noise, and we included a jitter term for each night of WIRC data log(σ extra ), accounting for the discrepancy between the photon noise and the true variance in the data (i.e., σ 2 = σ 2 photon + σ 2 extra ).Using the NUTS sampler (running 4 chains with 1000 tuning and 1500 sampling steps per chain) we obtain the posterior probability distributions for the joint fit with R < 1.01.The detrended light curve and residuals are displayed in Figure 2. The Allan deviation plot for the joint fit of the WIRC and TESS data is displayed in Figure 3.The priors used and the posteriors obtained are given in Table 1.

RESULTS AND MASS-LOSS MODELLING
Our results indicate that TOI 1268b does indeed have an escaping atmosphere; in this section we explore the implications of our measurement for the planet's longterm evolution.Following the approach in Linssen et al. (2022), we convert the measured excess absorption in the helium bandpass to an atmospheric mass-loss rate by modelling the outflow with a 1-dimensional Parker wind (Oklopčić & Hirata 2018;Lampón et al. 2020) using Cloudy (Ferland et al. 1998(Ferland et al. , 2017)).For a given mass-loss rate Ṁ , a constant thermosphere temperature T , and a hydrogen fraction for the outflow f H that we keep fixed at 0.9, we can compute the corresponding density and velocity profiles for the atmospheric wind.Using the MUSCLES spectrum from ϵ Eridani (France et al. 2016;Loyd et al. 2016;Youngblood et al. 2016), which has a similar spectral type and age to that of TOI-1268, we can then run Cloudy to calculate a non constant temperature structure that we compare to the assumed constant value to constrain a sensible thermo-spheric temperature.Cloudy also calculates the radially varying metastable helium density, from which we calculate the wavelength-dependent absorption signal.We average this signal over our 0.635 nm bandpass, and compare the predicted excess absorption to our measured value.We compute a grid of models spanning mass-loss rates between 10 9 < Ṁ < 10 11 g• s −1 and temperatures between 5000 < T < 11000 K.We then compare the predictions of this model grid to our measured excess absorption signal (see Figure 4).We refer to sections 3.2 and 4.2 in Linssen et al. (2022) for additional details on the model framework.We find that our observations are best-matched by a mass-loss rate of log( Ṁ ) = 10.2 ± 0.3 g s −1 and a corresponding thermosphere temperature of 6900 +1800 −1200 K. Following the approach described in section 5.3 of Vissapragada et al. ( 2022), we can convert this mass-loss rate into the corresponding mass loss efficiency ε (see Eq. 2 of Vissapragada et al. 2022).However, this is dependent on the incident XUV flux from the host star (see Eq. 1).If we assume that the mass loss efficiency has remained constant over time, we can then estimate the planet's cumulative mass loss from formation to the present day.This calculation is likely an overestimate of the true mass loss, as the mass loss efficiency should be lower at early times when the star is more active and the outflow is more highly ionized (Owen & Jackson 2012).Note-BTJD = BJDTDB -2457000.We omitted the detrending weights for each night.
Integrating over the typical expression for the energylimited mass-loss rate (see Eq. 2 of Vissapragada et al. 2022) we can calculate the mass of the planet as a function of the fraction of mass lost to date, f (by defining the planet's initial mass as M 0 = M p /f ): where ε is the mass loss efficiency, R p is the radius of the planet, E XUV is the integrated stellar XUV luminosity, a is the semi-major axis, K is the Erkaev et al. (2007) Roche lobe correction factor, and G is the gravitational constant.
We estimated the time-integrated high energy flux E XUV for the star using the fast rotator track for K dwarfs from Johnstone et al. (2021) as a guide, and adopt a value of 10 46 erg.Using this value, we find that TOI-1268b has only lost 0.45% of its initial mass.This is also likely an overestimate, as the E XUV values reported in this study are integrated up to 1 Gyr, much older than TOI-1268b's estimated system age of 110-380 Myr ( Šubjak et al. 2022).Furthermore, in our calculation we have adopted a fixed mass-radius relation as do (Owen & Lai 2018).Thorngren et al. (2023) argue that for hot planets (> 1000K), radius inflation is also time-dependent and should be integrated over alongside the luminosity in Eq. 2. However, at the T eq ∼ 900K of this planet ( Šubjak et al. 2022), we are able to circumvent this issue, which adds some confidence to our interpretation.
We can also look forward in time to determine whether or not TOI-1268b is expected to survive over the rest of the star's main sequence lifetime.If we assume a constant mass-loss rate going forward (also an overestimate, as the star's activity level should decrease as it ages), we find that the time required for this planet to lose all of its mass to atmospheric escape is 584 Gyr.This is comparable to the predicted atmospheric lifetimes of other close-in gas giant planets presented in Vissapragada et al. (2022) and Kirk et al. (2022).We conclude that TOI-1268b is quite stable against photoevaporation over the main-sequence lifetime of its host star.

DISCUSSION AND CONCLUSIONS
In this study we presented the first detection of atmospheric mass loss from TOI-1268b, a young Saturnmass planet located near the upper edge of the Neptune desert ( Šubjak et al. 2022).We detect excess helium absorption during the transit with a significance of 5.7σ, and translate this excess absorption measurement into a corresponding mass-loss rate of log( Ṁ ) = 10.2 ± 0.3 g s −1 and a thermosphere temperature of 6900 +1800 −1200 K.We estimate that TOI-1268b has lost less than 1% of its total primordial envelope mass, and if it continues losing mass at the current rate, it will lose less than 1% of its current envelope mass over the main-sequence lifetime of its host star.This is consistent with the results of Vissapragada et al. (2022) and other literature mass loss studies (Owen & Lai 2018;Ionov et al. 2018;Dos Santos 2022), which find that photoevaporation is negligible for gas giant planets near the upper edge of the Neptune desert.With a density of 0.71 +0.17 −0.13 g cm −3 ( Šubjak et al. 2022), our results are also consistent with predictions from Thorngren et al. ( 2023) that planets with densities higher than 0.1 g cm −3 should be stable against photoevaporation.
Our results indicate that TOI-1268b could have either formed in situ or arrived at its current location very early on via disk migration and still easily survive to the present day.However, it is also possible that it has only recently migrated inward from a more distant formation location; this would further reduce its total atmospheric mass loss.This planet appears to have a modest orbital eccentricity (0.105 +0.040 −0.039 , Šubjak et al. 2022), potentially indicating that it arrived via high eccentricity migration and has only partially circularized.However, its orbit is also well-aligned with that of its host star (Dong et al. 2022), which would tend to disfavor Kozai-type migration channels.In order to better constrain its past migration history, it would be helpful to obtain additional radial velocity measurements.By extending the current radial velocity baseline, these observations could also be used to search for the presence of outer companions in the system that might have perturbed TOI-1268b's orbit.Furthermore, we calculate the tidal realignment timescale for TOI-1268b, τ R (Albrecht et al. 2012).We obtain a timescale of 6.35 • 10 5 Gyr which is two orders of magnitude larger than the estimated age of the planet (110-380 Myr).Therefore realignment is highly unlikely.
Our observations of TOI-1268b represent a meaningful addition to the small sample of young transiting gas giant planets with measured atmospheric outflows.In the future, spectroscopically-resolved observations of TOI-1268b's helium absorption signal could be used to explore how the strong winds, magnetic fields and high activity levels of young stars act to sculpt the 3D outflow geometry (Schreyer et al. 2023).Similar observations of outflows from young sub-Neptune-sized planets suggest that they may have complex structures that are not readily matched by standard 3D models (Zhang et al. 2022a(Zhang et al. ,b, 2023)); it would be interesting to know if young gas giant planet outflows possess similar features.Although we did not find any evidence of an extended comet-like tail in our data, spectroscopic observations should be more sensitive to these kinds of features (Nortmann et al. 2018;Spake et al. 2018Spake et al. , 2021;;Kirk et al. 2020).As the ongoing TESS survey continues to identify new young transiting planet systems, it is also likely that there will be additional targets like TOI-1268b that are amenable to mass loss studies using metastable helium.By obtaining a larger sample of mass loss measurements for young transiting gas giant planets, we can obtain a clearer view of early atmospheric mass loss processes on close-in exoplanets.

Figure 1 .
Figure 1.Mass versus orbital period for the population of transiting exoplanets with fractional mass uncertainties smaller 30%.The Neptune desert boundaries from Mazeh et al. (2016) are overplotted as dashed lines.Planets orbiting stars with 4000 K< T eff < 5400 K are the most favorable targets for 1083 nm metastable helium observations; these systems are outlined in red.The seven planets with previously published Palomar/WIRC helium measurements tabulated in Vissapragada et al. (2022) and TOI-1268b are outlined in blue.Grey shading indicates the magnitude of each star, and the point size is proportional to the logarithm of the predicted signal-to-noise ratio in the WIRC helium bandpass.Figure adapted from Vissapragada et al. (2022).
3. LIGHT CURVE MODELINGWe modeled each individual night of Palomar/WIRC data using the exoplanet package (Foreman-Mackey

Figure 2 .
Figure 2. Left: Combined WIRC light curve for TOI-1268b including data from both Palomar nights, with circle markers for night 1 and triangle markers for night 2. The blue curve is the best-fit transit model with the radius ratio from the TESS bandpass, and the red curve is the best-fit transit model with the radius ratio in the WIRC helium bandpass, with the 1σ uncertainty in the transit depth shown as light red shading.Right: TESS data with the joint WIRC+TESS model showed in blue.Unbinned photometry is shown in grey, while black points are binned to five minutes.

Figure 3 .
Figure 3. Allan deviation plot for the joint WIRC+TESS fit.The black points indicate the root mean square error of the binned residuals for each light-curve, solid red line indicates the expected noise from photon noise statistics and the dashed red line indicates the photon noise scaled up the root mean square error of the binned residuals.

Figure 4 .
Figure 4. Combined constraint on the mass-loss rate and temperature.The orange dark and light areas are 1σ and 2σ contours from the temperature constrained profile.The blue dark and light areas are the best-fit models, 1σ and 2σ contours of the mass loss-rate profiles.The 1σ contour of joint constraint is shown in black.The 1σ constraints of simulations ran at three times and one third of the EUV flux are shown in gray

Table 1 .
Priors and posteriors for Palomar/WIRC joint and TESS -only fits.