Scaling K2. VII. Evidence For a High Occurrence Rate of Hot Sub-Neptunes at Intermediate Ages

The NASA K2 mission obtained high-precision time-series photometry for four young clusters, including the near-twin 600–800 Myr old Praesepe and Hyades clusters. Hot sub-Neptunes are highly prone to mass-loss mechanisms, given their proximity to the host star and the weakly bound gaseous envelopes, and analyzing this population at young ages can provide strong constraints on planetary evolution models. Using our automated transit detection pipeline, we recover 15 planet candidates across the two clusters, including 10 previously confirmed planets. We find a hot sub-Neptune occurrence rate of 79%–107% for GKM stars in the Praesepe cluster. This is 2.5–3.5σ higher than the occurrence rate of 16.54−0.98+1.00 % for the same planets orbiting the ∼3–9 Gyr old GKM field stars observed by K2, even after accounting for the slightly supersolar metallicity ([Fe/H] ∼ 0.2 dex) of the Praesepe cluster. We examine the effect of adding ∼100 targets from the Hyades cluster and extending the planet parameter space under examination, and we find similarly high occurrence rates in both cases. The high occurrence rate of young, hot sub-Neptunes could indicate either that these planets are undergoing atmospheric evolution as they age, or that planetary systems that formed when the Galaxy was much younger are substantially different than from today. Under the assumption of the atmospheric mass-loss scenario, a significantly higher occurrence rate of these planets at the intermediate ages of Praesepe and Hyades appears more consistent with the core-powered mass-loss scenario for the origin of the planet radius valley, compared to the photoevaporation scenario.


INTRODUCTION
After successive hardware failures rendered the NASA Kepler spacecraft (Borucki et al. 2010) incapable of continuing to point with high precision and stability at its original field of view, the K2 mission repurposed the spacecraft to perform a four year survey of 18 fields, called 'campaigns', around the ecliptic plane (Howell et al. 2014).This significant expansion in sky coverage compared to the original Kepler mission led to a corresponding increase in both the number and variety of targets for which high precision time-series photometry could be obtained.Where the goal of Kepler was to examine the planet occurrence rates for solar-type stars, the target list for K2 was provided by the community to meet a variety of science goals.The goal of the Scaling K2 project is the exploitation of this expanded target set to search for planet candidates in new areas of stellar parameter space beyond that probed by Kepler.As the K2 campaigns were limited to 70-80 days in duration, and the photometric noise achieved with K2 was higher than that achieved with Kepler, our sensitivity is necessarily limited to short-period (<40 day) planets down to super-Earth-size.
In Papers I-V we presented the uniformly derived set of stellar parameters for the K2 targets comprising our search sample (Hardegree-Ullman et al. 2020), the automated detection pipeline for transiting candidates, with wellcharacterized completeness and reliability (Zink et al. 2020), a pilot study on a single campaign (Zink et al. 2020), the full, homogeneous planet candidate catalog from campaigns 1-8 and 10-18 (Zink et al. 2021), and 60 newly validated planets from the catalog (Christiansen et al. 2022).In Paper VI (Zink et al. 2023) we compared the FGK occurrence rates for Kepler and K2 in their regions of stellar and planet parameter overlap and found good agreement between the two samples, validating our approach to date, and providing independent confirmation of the trends observed in Kepler data.In addition, we found that the occurrence rate of small (sub-Neptunes and super-Earths), short (1-40 d) period planets is inversely correlated with the galactic oscillation amplitude of the host star-the higher a star reaches above the plane in its galactic orbit, the less likely it is to host small, short-period planets.
Another axis along which K2 can probe is that of stellar age-in its four-year mission it observed four young clusters: Upper Scorpius (5-10 Myr; ∼1100 members observed), Pleiades (120 Myr; ∼850 members), Hyades (600-800 Myr; ∼300 members) and Praesepe (600-800 Myr; ∼1000 members).As different planet formation and evolution scenarios are associated with different timescales, measuring occurrence rates of planets as a function of age provides particularly useful constraints on the viability of these scenarios.One particular current area of interest is the origin of the planet radius valley and sub-Neptune desert for small, short-period planets (Fulton et al. 2017;Van Eylen et al. 2018;Hardegree-Ullman et al. 2020).The leading candidates for the explanation of this feature are planetary atmosphere evolution in the form of either photoevaporation (Lopez et al. 2012;Lopez & Fortney 2013;Owen & Wu 2013, 2017;Rogers & Owen 2021) or core-powered mass loss (Ginzburg et al. 2016;Gupta & Schlichting 2019, 2020, 2021).These mechanisms are expected to act over very different timescales; Ginzburg et al. (2016) suggested that analyzing occurrence rates of young planetary systems would be a useful way to discriminate between the two.In this paper, we measure the occurrence rate for small, short-period planets (hot sub-Neptunes) in the Praesepe cluster, subsequently augmenting the analysis with the significantly smaller number of Hyades cluster members.The paper is organized as follows: in Section 2 we describe the cluster stellar sample and corresponding planet candidates retrieved by our pipeline; in Section 3 we describe the occurrence rate calculations; and in Section 4 we discuss the implications of the results.

Praesepe
Praesepe (NGC 2632, M44) is a 600-800 Myr-old open cluster located about 180 pc away.The age, distance, and metallicity of the cluster have been studied extensively (see, e.g., Brandt & Huang 2015;Gossage et al. 2018, and references therein).The cluster has slightly super-solar metallicity, with estimates ranging from [Fe/H] = 0.12±0.04(Boesgaard et al. 2013) to 0.21±0.01(D'Orazi et al. 2020), depending on the calibration and methods used.Fortuitously lying in the ecliptic, it was observed by K2 in campaigns 5, 16, and 18 (the latter of which had the same pointing as campaign 5).We use the member list compiled by Rebull et al. (2017) in their analysis of the rotation periods of Praesepe members in K2 data, which originally included the targets observed in campaign 5 (and 18), and was subsequently augmented by additional members observed in campaign 16. Figure 1 shows the spatial distribution of the 1030 Praesepe targets observed by K2 across the sky.
As part of our search of the full K2 dataset, presented in Zink et al. (2021), the Praesepe targets were processed through our standard K2 data reduction and transit detection pipeline detailed in Zink et al. (2020); in brief, light curves are extracted using EVEREST (Luger et al. 2016;Luger et al. 2018), detrended using Gaussian process regression and harmonic removal, searched for periodic transit signals using TERRA (Petigura et al. 2013), and the resulting signals uniformly vetted using the EDI-Vetter tool (Zink et al. 2020).Although there is some overlap between the campaigns, for this analysis we treat each campaign separately, rather than combine light curves from multiple campaigns; this limits our period sensitivity to below 40 days but greatly simplifies the completeness and reliability corrections in the  2020), and were derived using a random forest regression on photometry, trained on a large set of targets with LAMOST spectra; where the machine learning parameters were unavailable, they were supplemented with parameters from isochrone fitting.subsequent analysis.We compared the noise properties of the Praesepe targets to the field K2 field star dataset (after removing the known young stars in Upper Scorpius, Pleiades, Hyades and Praesepe)-young stars rotate more rapidly and therefore, although the Praesepe members studied here are largely on the main sequence, their light curves could potentially have increased correlated noise.This could make transit signals more difficult to detect than would be reflected by the detection efficiency previously measured for the pipeline across all targets, and potentially necessitate an updated detection efficiency measurement specifically for the Praesepe targets.We examine the Combined Differential Photometric Precision (CDPP) values (Christiansen et al. 2012) established for Kepler light curves.The left panel of Figure 2 shows the 8-hour root-mean-square (rms) CDPP values as a function of the Kepler magnitude for the Praesepe targets in blue, overlaid on an orange heatmap showing the 8-hour rms CDPP values for our full K2 dataset.The Praesepe targets largely follow the broad features of the overall sample, except at the very faintest magnitudes.Zink et al. (2021) noted that stars with an 8-hour rms CDPP above 1200 ppm do not contribute meaningfully to the derived occurrence rates, as they have no sensitivity to the planet parameter range under study 1 .Applying the same CDPP cut as Zink et al. (2021) removes 429 targets, largely very faint M dwarfs smaller than 0.33 R ⊙ .The 601 remaining Praesepe targets are shown in the right hand panel of Figure 1.
In addition to the overall scatter in a light curve, we can also measure the extent of the correlated noise.Burke & Catanzarite (2017) showed that calculating the CDPP "slope," by measuring how the rms CDPP values change as a function of the timescale over which they are calculated, could capture this information.For white noise, the rms CDPP values should decrease with increasing duration as N short /N long , where N short and N long are the shorter and longer durations over which the CDPP is calculated, respectively.The right panel of Figure 2 shows the normalized distributions of the CDPP slopes for the 2-hour and 8-hour rms CDPP values of the 601 Praesepe targets in blue and the remainder of the K2 targets in orange.White noise would produce a CDPP slope of 2/8 = 0.5, and both populations peak at just above that value (and drop off rapidly below it), indicating largely well-behaved light curves.The Praesepe stars have a slightly broader distribution up to a CDPP slope of ∼1, indicating a marginally higher proportion of light curves with correlated noise than the full dataset, although above ∼1.1 (where binning over a longer timescale actually increases the noise, indicating highly correlated noise) the fractional proportion of the two populations is very similar.These comparisons indicate that the noise properties of the Praesepe targets are similar to For white noise, the CDPP slope is expected to be 0.5-both distributions peak just above 0.5 and lie primarily between 0.5 and 1.0, indicating non-pathological light curves.The Hyades targets we investigated follow comparable noise patterns.
the full K2 dataset and that the detection pipeline completeness and reliability measured on the full dataset in Zink et al. (2020) are applicable here.
The pipeline detected 10 planet candidates orbiting nine of the 601 Praesepe targets, originally published in the full K2 catalog provided in Zink et al. (2021).The candidates are listed in Table 1-they range from 1.2-3.9R ⊕ , with periods from 1.67-21.2days, typical of the large population of warm sub-Neptunes and super-Earths revealed by Kepler.Given the interest in detecting and characterizing planets earlier in their formation and evolution histories, the K2 Praesepe data have been extensively searched for individual planets, and seven of the 10 candidates are already confirmed in the literature, as indicated in Table 1.However, this paper presents the first statistical sample of planet candidates from Praesepe for population analysis.We note that one of the unconfirmed planet candidates (EPIC 211885995.01)has a consistency score falling below our cut-off of 0.75 (see Section 3.13 of Zink et al. 2020), and indeed a closer examination of the light curve indicates that the transit signal is likely related to stellar variability; we remove this target from further analysis, leaving nine remaining candidates shown in Figure 3.

Hyades
In campaigns 4 and 13, K2 also observed the Hyades open cluster, a smaller, closer, near-twin of Praesepe.Age estimates range from 625±25 Myr (Perryman et al. 1998), to 680 Myr (Gossage et al. 2018), to 800 Myr (David & Hillenbrand 2015;Brandt & Huang 2015).It also has a slightly super-solar metallicity of [Fe/H] = 0.14±0.05dex (Perryman et al. 1998), and has also been meticulously combed for planets.A list of potential Hyades members was assembled from membership lists in the literature (e.g.Douglas et al. 2016;Lodieu et al. 2019;Röser et al. 2011Röser et al. , 2019;;Goldman et al. 2013, among others), the list of Hyades members submitted to the original K2 proposal call, and a run of BANYAN Σ (Gagné et al. 2018), then compared to the list of objects with K2 light curves (Rebull et al. in prep), for a total of 300 targets.This is significantly smaller than our Praesepe target list, in part because Hyades is closer and therefore subtends a larger area on the sky, a smaller fraction of which falls in the K2 field of view; many Hyades cluster members were not observed by K2.Applying the same noise and stellar parameter cuts as before leaves 101 Hyades targets to add to our sample.We examined the noise properties of these targets, as above for Praesepe, and found them to be comparable.Similarly to our Praesepe sample, the Hyades targets are on average cooler than the full K2 sample.
Our pipeline recovered five planet candidates from these 101 Hyades targets, listed in Table 1.Three are known planets, K2-155 b and c (Hirano et al. 2018) and K2-136 c (Mann et al. 2018), although in both cases our pipeline missed known smaller or longer period additional known planets in the system.The final two are a pair of small  3. OCCURRENCE RATE ANALYSIS

Parametric Model
In order to calculate occurrence rates, we use the ExoMult forward modeling software (Zink et al. 2019(Zink et al. , 2020)), which, in brief, produces synthetic K2 planet samples drawn from an underlying parametric planet distribution function, and simulates their detection throughput using the completeness measured in Zink et al. (2021).The synthetic population of planets and the empirical sample were compared using the two-sample Andersen-Darling test statistic and then optimized for each relevant observable.Since the aim of the analysis is to compare the intermediate-age planet sample to the field K2 sample and to isolate any age correlation, careful attention must be paid to the form of the parametric model to ensure it captures the differences between the samples to the extent possible.
Given its significantly larger dataset, we begin by examining the Praesepe sample.The Praesepe candidates orbit dwarf stars from G1-M3, so we restrict our population analysis to GKM stars spanning the range 0.33 < R * < 1.5 R ⊙ , 3500 < T eff < 6500 K, −0.5 < [Fe/H] < 0.5 dex, and logg > 4.0.Our stellar parameters were uniformly computed for the K2 sample as described in Hardegree-Ullman et al. (2020), using a random forest regression on photometric colors trained on stars with LAMOST spectra, ultimately providing spectral types, effective temperatures, surface gravities, metallicities, radii, and masses.Where the machine learning parameters were unavailable, they were supplemented with parameters from isochrone grid matching, as per Zink et al. (2023).For the well-behaved Praesepe targets, this reduces the dataset to nine planet candidates and 487 stars.For the well-behaved (8hr CDPP <1200 ppm) targets in the full K2 field star dataset (with the stars in the four identified young clusters removed), this leaves 520 planet candidates and 108,502 stars, refered to as the 'field K2 sample'.The rest of the analysis will build on these two samples.The systematic differences in the stellar samples necessitate an appropriate parametric model that can account for any known correlations.
The two stellar samples have significantly different distributions in stellar temperature and metallicity, shown in Figure 4.There are well established correlations between the frequencies of planets over these stellar parameters.The median iron abundance of our Praesepe sample is 0.146 dex, consistent with the slightly super-solar metallicity of the cluster.Although sub-Neptunes do not have as steep a metallicity correlation as giant planets (Petigura et al. 2018;Zink et al. 2023), there is still a positive correlation that could lead to an enhancement in the frequency of small planets in Praesepe compared to the field (slightly sub-solar) K2 sample if not taken into account.Similarly, the field K2 sample includes a much larger fraction of higher mass stars than the Praesepe sample, and stellar mass has been shown to be anti-correlated with the frequency of small, short-period planets (e.g.Mullally et al. 2015).To account for these differences we include host star metallicity and stellar effective temperature terms in our parametric model.
For the parameter space probed by the nine Praesepe candidates (1-4 R ⊕ , 1-30 days), previous analyses using the older and much larger Kepler sample have found that the population is well represented by a broken power law in orbital period with a break at ∼10 days (Mullally et al. 2015), and a single power law in planet radius (e.g.Petigura et al. 2018;Mulders et al. 2018).With only nine Praesepe planet candidates, our statistical power to simultaneously fit many parameters is low, and indeed an attempt to fit the full parametric model with a broken power law in orbital period and single power law planet radius was unable to meaningfully constrain the overall occurrence rate, let alone the exponents of the various power laws.However, by concentrating on hot sub-Neptunes, we can simplify the parametric fit in order to make more robust measurements.For the following analysis, we examine candidates with orbital periods 1-10 days, inwards of the expected period break and therefore able to be modeled with a single power law, and planet radii 1.7-4 R ⊕ , capturing the sub-Neptune population above the planet radius valley (Fulton et al. 2017;Van Eylen et al. 2018;Hardegree-Ullman et al. 2020).These planets are particularly susceptible to mass-loss mechanisms, due to their proximity to their host stars, and their weakly-bound gaseous upper atmospheres.In addition, the same parameter range is investigated in Zink et al. (2023), which allows for a straightforward comparison to their results.
In this space, we fit a parametric model n of the form where f captures the number of planets per star within the range of our sample, and α, β, γ, and λ are all free parameters, as in Zink et al. (2023).We start by fitting the field K2 sample, allowing all tunable parameters to float.
In our targeted stellar and planet parameter space, this leaves 149 planet candidates from Zink et al. (2021).Examining the corner plot shows convergence of the MCMC parameter chains and highlights their independence.The fitted values are given in Table 2.For the field K2 sample, the occurrence rate is 16.54 +1.00 −0.98 %, and we recover the strong negative radius scaling (α = −2.6 ± 0.2) and positive period power-law scaling (β = 2.1 ± 0.1) seen in previous studies (e.g.α = −2.7 and β = 1.97 from Zink et al. (2023).We also see the expected weak but positive dependence on stellar metallicity (0.19 ± 0.09), and the negative dependence on stellar effective temperature (−0.072 +0.018  −0.017 ) 3.2.Praesepe In the Praesepe sample, we have three candidates in this parameter space (K2-100 b, K2-104 b, and K2-264 b), shown in the red dashed box in Figure 3.Even our simplified model is poorly constrained with only three candidates.We attempt to fit all five parameters (f , α, β, γ, and λ), which gives an occurrence rate of 95%, however, not unexpectedly, the fit does not converge, specifically for β (the period power-law coefficient).This is perhaps driven by the fact that the three candidates occupy only the middle ∼half of the 1-10 day range.We proceed by first fixing the radius power law coefficient, α, while fitting the remaining parameters, then fixing β while fitting the remaining parameters, and finally fixing both α and β, using the values for α and β from the full K2 field sample.Under the assumption that the radius and period power law coefficients are independent, which is likely appropriate over this limited region of parameter space, this allows us to explore the extent to which we can test for differences between these values for Praesepe compared to the full K2 data set, and in the overall occurrence rates.
The results for all fits are shown in Table 2.The resulting occurrence rates vary from 79 +31 −25 % to 107 +32 −26 % of GKM stars hosting a hot sub-Neptune 2 .The uncertainties are large, due to the small size of the sample, but the rates are still 2.5-3.5σhigher than the 17% occurrence rate recovered from the full K2 sample.In each case, the parameters describing the correlations with period, planet radius, stellar effective temperature and metallicity (α, β, λ, and γ) lie within 1.3-1.6σ of the values from the full K2 sample, indicating that our (small) Praesepe sample cannot demonstrate any significant differences in the strength of these correlations at younger ages.The only significant difference between the two samples is in the overall occurrence rate, which is significantly higher for the young small (1.7-4 R ⊕ ), short-period (1-10 day) hot sub-Neptunes in Praesepe than for field stars.To test this, we perform a final fit where we fix all but the occurrence rate, f , assuming the null hypothesis that there is no difference in the correlations between the intermediate age and field stars, and find an occurrence rate of 100%.

Praesepe + Hyades
Given the high occurrence rate of hot sub-Neptunes found in Praesepe, we combine the Praesepe and Hyades targets to create a larger target list and test the result.The stellar effective temperature and metallicity distributions of the Hyades targets are shown in Figure 5.The median metallicity of our 101 Hyades targets is [Fe/H] = -0.002,which is consistent with the cluster being slightly super-solar on average, given the typical uncertainty (0.235 dex) on our metallicity values and the small size of the sample.In the planet parameter space defined in Section 3.1, we add one additional Hyades planet (K2-155 b) to our sample.We repeat the same sequence of fits as above, using the parametric model described in Section 3.1.We first allow all five tunable parameters to float, which again does not converge for β (the period power-law coefficient), the additional 6.3-day Hyades planet being close in period to the 5.8-day Praesepe planet and likely not providing additional constraints over the 1-10 day period range.We then fix the α parameter, the β parameter, both parameters, and finally all four non-f parameters.The fitted parameters α, λ and γ lie within 0.4-1.8σ of the values from the field K2 sample; the β value is poorly constrained.As before, we cannot exclude the null hypothesis that the shape of the underlying planet population is the same for intermediate age stars as for field stars.The occurrence rates range from 90 +31 −26 % to 121 +30 −27 %, again indicating a significantly (2.8-3.8σ)higher rate of hot sub-Neptunes around these intermediate-age stars.

DISCUSSION
We investigated the occurrence rate of hot sub-Neptunes in the intermediate-age Praesepe and Hyades clusters observed by K2.We derived significantly higher occurrence rates (79-121%) for the small (1.7-4 R ⊕ ), short-period (1-10 day) planets in our Praesepe and combined Praesepe+Hyades datasets than for our full field K2 sample (17%).Although less work has been done to age-date K2 field stars than Kepler stars, examining the ages of the field K2 planet host stars reveals of median of 5.6 Gyr, and a 68 percentile range of 2.8-9.0Gyr (NASA Exoplanet Archive 2023b), which is similar to the range of ages found for Kepler targets (see, e.g.David et al. 2022).We therefore identify a significant correlation between stellar age and the occurrence rate of hot sub-Neptunes.2023) also found a hint of a higher occurrence rate for small, short-period planets around young stars.They analyzed 5 nearby young clusters observed by TESS, with ages spanning from 15-450 Myr, and a median age of ∼45-50 Myr.Over a slightly larger planet parameter space (planet radii from 1.8-6 R ⊕ and orbital periods from 1-12.5 days), and a slightly higher mass sample of stars than the GKM sample studied here (FGK stars from 0.55-1.63M ⊙ ), they find an occurrence rate of 93 ± 38%3 , with the caution that they began by analyzing the clusters in which planets had already been discovered, which may have biased their derived occurrence rate.We repeat our analysis described above in this larger parameter space, as shown in Table 2 (albeit for our slightly different stellar mass range).This has the effect of removing the smallest planet from our previous sample (K2-104b; 1.77 R ⊕ ), and adding a longer period Praesepe planet (K2-95 b, 10.13 days; Obermeier et al. 2016).First, we fit the full K2 sample, finding an occurrence rate slightly increased from 17% to 21%.With the addition of the longer period planet, we are able to fit all five tunable parameters, and do so for the Praesepe and Praesepe+Hyades samples, finding occurrence rates of 87 +33 −25 % and 89 +31 −25 %, respectively.These are 2.7σ higher than the full K2 occurrence rate, and are in excellent agreement with the value from the younger TESS sample, shown in Figure 6, lending credence to the result that the occurrence rate of hot sub-Neptune planets is significantly higher around younger stars (<1000 Myr) than around 3-9 Gyr-old field stars.For consistency with our previous analysis, we also perform fits where the α and β parameters are fixed to the values found for the full K2 sample, shown in Table 2.For both the Praesepe and Praesepe+Hyades samples this has the effect of increasing the overall occurrence rate and increasing the strength of the negative trend with radius.
Identifying a correlation between stellar age and planet occurrence rates could imply one of several things.Potentially, there is some ongoing process after the birth of a planetary system, independent of when that system is born, that acts to sculpt and evolve planets over time, such that the planetary systems we observe in young clusters evolve to become the planetary systems we see in the field K2 sample.Or, that stars born at different ages of our Galaxy are born in  2023), to the intermediate age (600-800 Myr) K2 Praesepe GKM sample studied here, to the field K2 GKM stars (3-9 Gyr).The age uncertainties capture the width of the age distributions in each case.The shaded regions indicate the putative timescales for the photoevaporation (<100 Myr) and core-powered mass loss processes (0.5-2 Gyr).The persistent high rate of hot sub-Neptunes at ∼700 Myr is more consistent with predictions of core-powered mass loss sculpting this population than with photoevaporation.
sufficiently different conditions as to create different types and architectures of planetary systems in the first place.Here, we discuss these possibilities.

Planet atmosphere evolution
One intriguing property of small, short-period planets is their bimodal distribution in planet radius (Fulton et al. 2017;Van Eylen et al. 2018;Hardegree-Ullman et al. 2020), or perhaps more accurately trimodal in planet density (Luque & Pallé 2022).Many of the ideas proposed to explain this distribution involve mass loss from planets early in their histories, including through photoevaporation (Lopez et al. 2012;Lopez & Fortney 2013;Owen & Wu 2013, 2017;Rogers & Owen 2021) or core-powered mass-loss (Ginzburg et al. 2018;Gupta & Schlichting 2019, 2020, 2021).Perhaps the very high rate of hot sub-Neptunes found in Praesepe and Hyades is in the process of reducing to the smaller rate we see in the K2 field stars.
In the photoevaporation scenario, the gaseous atmosphere is stripped away by high energy (XUV) radiation from the host star.The extent of the mass loss, and the timescale over which it occurs, are driven by the initial planetary core mass, the initial envelope mass fraction, the orbital separation of the planet, and the stellar spectral energy distribution.Lopez et al. (2012) investigated the role of photoevaporation mass-loss for a specific system, finding that the 1.8 R ⊕ planet Kepler-11 b could have lost a factor of 2-3 in radius and 1.5-5 in mass between ∼100 Myr and 10 Gyr.Lopez & Fortney (2013) extended this study, predicting that planets in the range 1.8-4.0R ⊕ should become significantly less common on orbits <10 days.Subsequently, Rogers & Owen (2021) concluded that although some super-Earths are born rocky, approximately four times as many are born with large H/He atmospheres that are subsequently stripped by photoevaporation, a number that is tantalizingly in agreement with our finding of 4-5 times as many young sub-Neptunes as old sub-Neptunes, if a significant fraction of the deficit evolve to become super-Earths.
In the core-powered mass loss scenario, the gaseous atmosphere can be blown away by both the planet's primordial energy from formation and bolometric heating from the star.If the formation energy is similar to or greater than the gravitational binding energy of the atmosphere, the cooling luminosity of the planet may blow off a H/He atmosphere entirely.For planets with masses between Earth and Neptune, the core temperatures can be 10,000-100,000 K, which can take up to several Gyr to cool (Ginzburg et al. 2016).Gupta & Schlichting (2020) found that the average size of sub-Neptunes decreases significantly with age, and commensurately that the relative occurrence rate of sub-Neptunes decreased by a factor of 2-2.5 over the timescale of the mass loss (see their Fig.10), which is consistent at 1-2σ level with our results.
One obvious question is whether our finding, that the occurrence rate of hot sub-Neptunes is still very high at ∼700 Myr, before dropping by a factor of up to 5 at several Gyr, agrees with the putative timescales of these mass-loss mechanisms.Owen & Wu (2013, 2017) investigated the timescales of photoevaporation, finding that most atmospheric erosion happens within the first 100 Myr (see, e.g., Fig. 2 of Owen & Wu 2017).This would imply that planets shrink in size somewhat sooner than our Praesepe+Hyades result, but is not at odds with the younger TESS result.Core-powered mass loss is expected to occur over much longer timescales than photoevaporation, typically 0.5-2 Gyr (Ginzburg et al. 2016;Gupta & Schlichting 2019, 2020), which is much more consistent with our finding that at 0.7 Gyr these planets are still large.Gupta & Schlichting (2020) state "we expect a drastic change in the planet size distribution between stars that are younger and older than the typical core-powered mass-loss time-scale, which is of the order of a Gyr.Due to these long mass-loss time-scales, we predict the transformation of sub-Neptunes into super-Earths to continue over Gyr time-scales".Indeed, early examination of individual planets being detected in the Praesepe and Hyades K2 observations noted that they appeared to be larger, on average, than planets orbiting older stars of the same mass observed with Kepler (see, e.g.Fig. 12 of Rizzuto et al. 2018), which could have been evidence for ongoing thermal contraction of their atmospheres, and/or could have hinted at a preference for core-powered mass loss (Gupta & Schlichting 2020).However, without understanding the completeness of the survey detecting those planets, this could just as easily have been a selection bias, where the noisier K2 light curves limited the detection to larger planets when compared with the more well-behaved Kepler light curves.Berger et al. (2020) and Sandoval et al. (2021) used isochronal ages to divide the Kepler planet sample into coarse age bins, and found that the relative frequency of sub-Neptunes to super-Earths decreases with increasing age, i.e. that either the number of sub-Neptunes was decreasing, the number of super-Earths was increasing, or both.Here we extend the Kepler ages to the young clusters observed by K2 and find that there is indeed a "drastic change" (a factor of ∼4-5) in the absolute frequency of sub-Neptunes between stars that are younger and older than ∼1 Gyr.This result, specific to the Praesepe and Hyades clusters studied here, is more consistent with the published predictions of the core-powered mass loss origin for the planet radius valley rather than a photoevaporation origin, or perhaps points to some missing process in the photoevaporation model that would act to slow down the expected mass loss rate.Indeed, recent work by Owen & Schlichting (2023) conversely theorizes that planets that have undergone core-powered mass loss can then transition to photoevaporation, depending on the initial envelope mass fraction; whichever the order, the bulk of the mass loss is indicated to happen after the first 0.5 Gyrs.
Another factor that must be considered is the natal disk opacity, which plays an important role in both photoevaporation and core-powered mass loss processes.Planets that have formed from material that is richer in metals have the ability to maintain their envelopes for a longer period of time due to their increased heat capacity, enabling efficient cooling of their outer envelopes.Given that our cluster planet sample is more metal-rich than the [Fe/H] ∼ −0.2 of Kepler (Dong et al. 2014) and the [Fe/H] ∼ 0.0 of the field K2 sample, it is important to consider how this may affect our conclusions.According to Owen & Murray-Clay (2018), the envelope mass loss rate ( Ṁ ) dictated by the photoevaporation model scales with the disk metallicity (Z) as Ṁ ∝ Z −0.77 , while the mass loss originating from the residual heat in the planet's core scales as Ṁ ∝ Z −1 (Gupta & Schlichting 2020).This implies that the removal of sub-Neptune atmospheres will be delayed in a metal-rich environment, increasing the expected timescales of photoevaporation and core-power mass loss processes by 1.4× and 1.6×, respectively.4However, the surplus of hot sub-Neptunes at 600-800 Myr is more than 6× the fiducial timescale for photoevaporative mass loss (∼ 100 Myr).Therefore this potential delay, arising from the cluster's metallicity enhancement, cannot account for the abundance of hot sub-Neptunes identified in this study.
Finally, our sample covers a wide range of stellar masses, from 0.3 to 1.3 solar masses.The evolutionary timescales discussed above have been derived for an FGK sample of stars.The inclusion of low mass stars in our sample may modify the time frame over which hot sub-Neptunes are converted to super-Earths.Under the action of the photoevaporation framework, small M dwarfs experience a heightened integrated XUV flux, accelerating and/or increasing the magnitude of the mass loss process (Rogers et al. 2021).Qualitatively, this would reduce the expected timescale needed to dissociate the planetary envelopes of sub-Neptunes, culminating in an overall occurrence reduction in the 600-800 Myr clusters.However, our results show that the opposite is true.In contrast, the core-powered mass loss evolutionary process relies on the star's total and therefore expects decelerated and/or reduced mass loss around less massive stars (Gupta & Schlichting 2020).The inclusion of M dwarfs should increase the expected mass loss timescale and leave a larger reservoir of short-period sub-Neptunes at 600-800 Myr, matching the observation of this survey.Overall, the complexities of stellar mass and composition unique to our sample are directionally consistent with a core-powered mass loss mechanism, providing further evidence in favor of this formation pathway.

Primordial planet formation
An alternative explanation for the overabundance of hot sub-Neptunes in Praesepe and Hyades compared to the field K2 sample is that there could be systematic evolution of planetary system formation and architectures with different environmental/Galactic conditions over cosmic time.There are a number of known correlations between small, short-period planet occurrence rates and their host star properties-including properties that have evolved over the lifetime of the Galaxy, as the natal conditions for each generation of stars is impacted by previous generations.One such property is the host star metallicity-stars that are born today are typically formed in giant molecular clouds that have been enriched by the supernova remnants of earlier stars, in a more heavy-element-rich environment than in early Galactic times.This is reflected in the slightly super-solar average metallicities of Praesepe and Hyades, but was also explicitly fit as part of our parametric model, so presumably what we are seeing here is an additional effect.Both observations and models of galaxy formation in Milky Way-like galaxies suggest that stars born at earlier times are typically born in more strongly dynamically clustered environments (Lee & Hopkins 2020;Grudić et al. 2023), are much more strongly irradiated externally by elevated galactic star formation rates (Tacconi et al. 2020;Guszejnov et al. 2020), with likely larger ionizing cosmic ray fluxes (Hopkins et al. 2021) and systematically warmer natal gas clouds (Guszejnov et al. 2020).There is a long history of work on how this may or may not influence the low-mass end of the stellar initial mass function, but until recently little work on its possible effects on planet formation.Possibly, early and frequent interactions with nearby stars destabilize planetary systems to the extent that the overall occurrence rate, even for planets on the close-in orbits probed by our analysis, is reduced.Zink et al. (2023) recently discovered a strong negative correlation between the occurrence rate of 1.7-4 R ⊕ , 1-10 day planets (the same planet parameter space probed here) and the maximum height above the Galactic plane that the host stars reach as they oscillate through the plane.This height can be used as a proxy for thick disk (older, more metal poor, more likely to have formed in more clustered stellar environments) versus thin disk (younger, more metal rich, less likely to have formed in more clustered environment) membership, and since Zink et al. (2023) also accounted for metallicity, points to an additional line of evidence that older stars have fewer small, short-period planets when compared to younger stars, although as here there was no systematic difference identified in the shape of the planet population.Although there has been less work on quantitative predictions for these primordial planet formation models compared to the atmosphere evolution models, nothing in this work disfavors such explanations as a category.

CONCLUSIONS
We analyze K2 observations of the Praesepe and Hyades clusters, and measure a significantly higher occurrence rate for hot sub-Neptunes orbiting intermediate (600-800 Myr) GKM stars than for the older (3-9 Gyr) GKM field stars.Notably, our result persists after the inclusion of additional data, when varying the boundaries of the planet parameter space under inspection, and regardless of whether the shape of the population is assumed to be the same as that of the field stars or allowed to vary.If the decline in the number of hot sub-Neptunes with age is a result of planet evolution, the high rate of these planets at intermediate ages is more consistent with the predicted timescales of the core-powered mass loss mechanism than the photoevaporation mechanism.If it is instead a reflection of the primordial planet populations in these clusters compared to the field stars, it may point to some other trend with planet formation as a function of Galactic age, dynamical or radiation history, or cluster environment.One interesting test would be to search Praesepe and Hyades for longer period planets, such as with radial velocity, astrometry, or direct imaging.The putative mechanisms for planet evolution happen predominately at close-in distances to the host star-if the outer regions of the intermediate age planetary systems look similar to the older planetary systems, that could point to evolution being the explanation, rather than primordial differences in planetary system architectures.Facilities: Kepler Software: astropy (Astropy Collaboration et al. 2013, 2018), ExoMult (Zink et al. 2019), BANYAN Σ (Gagné et al. 2018)

Figure 1 .
Figure 1.Left: The spatial distribution of the Praesepe cluster members on the sky observed by K2 ; targets observed in campaigns 5 and 18 are shown in blue, and targets observed in campaign 16 are shown in green.Right: The temperature, radii, and metallicity of the Praesepe cluster members with 8-hour rms CDPP < 1200 ppm.The stellar parameters are from the uniform catalog of K2 stellar parameters published in Hardegree-Ullman et al. (2020), and were derived using a random forest regression on photometry, trained on a large set of targets with LAMOST spectra; where the machine learning parameters were unavailable, they were supplemented with parameters from isochrone fitting.

Figure 2 .
Figure 2. Left: The distribution of 8-hour rms CDPP values for the Praesepe targets in blue, and for the full K2 dataset in orange.The 1200-ppm cut-off is shown as the dashed line-it primarily removes faint, late-type M dwarfs.Right: The distribution of 8-hour/2-hour CDPP slope values for the targets below the 1200-ppm cut-off, with Praesepe shown in blue and the full K2 dataset in orange.For white noise, the CDPP slope is expected to be 0.5-both distributions peak just above 0.5 and lie primarily between 0.5 and 1.0, indicating non-pathological light curves.The Hyades targets we investigated follow comparable noise patterns.

Figure 3 .
Figure 3. Praesepe planet candidates are shown as red stars, and Hyades planet candidates are shown as yellow plus symbols.The full K2 planet candidate catalog of Zink et al. (2020) is shown in blue, and all confirmed planets at the NASA Exoplanet Archive (2023a) are shown in grey.The red box shows the hot sub-Neptune parameter space (1-10 days, 1.7-4 R⊕) analysed in Section 3.

Figure 4 .
Figure 4.The normalized distributions of the Praesepe GKM sample in orange and the field K2 GKM sample in blue, for stellar effective temperature (left) and stellar metallicity (right).The vertical dashed lines in the right panel are the median values of the distributions.The known super-solar metallicity of the Praesepe cluster is clear, although the distribution of our photometrically-derived values is wider than what is seen spectroscopically(Boesgaard et al. 2013;D'Orazi et al. 2020).The systematic differences in the stellar samples necessitate an appropriate parametric model that can account for any known correlations.

Figure 5 .
Figure 5.As for Figure 4, but for the significantly smaller Hyades GKM sample.

Figure 6 .
Figure 6.The occurrence rates of hot sub-Neptunes with time, from the young TESS FGK (15-450 Myr) from Fernandes & Hardegree-Ullman et al. (2023), to the intermediate age (600-800 Myr) K2 Praesepe GKM sample studied here, to the field K2 GKM stars (3-9 Gyr).The age uncertainties capture the width of the age distributions in each case.The shaded regions indicate the putative timescales for the photoevaporation (<100 Myr) and core-powered mass loss processes (0.5-2 Gyr).The persistent high rate of hot sub-Neptunes at ∼700 Myr is more consistent with predictions of core-powered mass loss sculpting this population than with photoevaporation.
This research has made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under Exoplanet Exploration Program.This paper includes data collected by the Kepler mission and obtained from the MAST data archive at the Space Telescope Science Institute (STScI).Funding for the Kepler mission is provided by the NASA Science Mission Directorate.STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555.This material is based upon work supported by the National Aeronautics and Space Administration under Agreement No. 80NSSC21K0593 for the program "Alien Earths".The results reported herein benefited from collaborations and/or information exchange within NASA's Nexus for Exoplanet System Science (NExSS) research coordination network sponsored by NASA's Science Mission Directorate.JZ acknowledges the support provided by NASA through the Hubble Fellowship grant HST-HF2-51497.001 awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, In., for NASA, under the contract NAS 5-26555.Support for PFH was provided byNSF Research Grants 1911233 & 20009234, NSF CAREER grant 1455342, NASA grants 80NSSC18K0562, HST-AR-15800.001-A.

Table 1 .
Summary parameters of the 15 planet candidates and their host stars identified in Praesepe and Hyades.C# = campaign in which the signal was detected by our pipeline.CP name = confirmed planet name in the literature.a Mann et al. (2017); b Obermeier et al. (2016); c Rizzuto et al. (2018); d Hirano et al. (2018); e Mann et al. (2018); * removed from further analysis due to low consistency score.

Table 2 .
The measured parameters of the parametric model for the field K2 the Praesepe sample, and the combined Praesepe+Hyades sample, over the two planet parameter ranges defined in the analysis.Bold values indicate parameters that were fixed in each analysis.†: A repeat of the analysis on the line immediately above without the CDPP<1200 ppm cut, to confirm that it does not significantly affect the result.