Comparative studies of young and old exoplanet populations offer a glimpse into how planets may form and evolve with time. We present an occurrence rate study of short-period (<12 days) planets between 1.8 and 10 R ⊕ around 1374 FGK stars in nearby (200 pc) young clusters (<1 Gyr), utilizing data from the Transiting Exoplanet Survey Satellite mission. These planets represent a population closer to their primordial state. We find that the occurrence rate of young planets is higher ( 6 4 − 22 + 32 %) compared to the Gyr-old population observed by Kepler ( 7.9 8 − 0.35 + 0.37 %). Dividing our sample into bins of young (10–100 Myr) and intermediate (100 Myr–1 Gyr) ages, we also find that the occurrence distribution in orbital period remains unchanged, while the distribution in planet radius changes with time. Specifically, the radius distribution steepens with age, indicative of a larger planet population shrinking due to the atmospheric thermal cooling and mass loss. We also find evidence for an increase (1.9 σ ) in occurrence after 100 Myr, possibly due to tidal migration driving planets inside of 12 days. While evidence suggests that postdisk migration and atmospheric mass loss shape the population of short-period planets, more detections of young planets are needed to improve statistical comparisons with older planets. Detecting long-period young planets and planets <1.8 R ⊕ will help us understand these processes better. Additionally, studying young planetary atmospheres provides insights into planet formation and the efficiency of atmospheric mass-loss mechanisms on the evolution of planetary systems.
ABSTRACT Determining the precise ages of young (tens to a few hundred Myr) kinematic (‘moving’) groups is important for placing star, protoplanetary disc, and planet observations on an evolutionary timeline. The nearby ∼25 Myr-old β Pictoris Moving Group (BPMG) is an important benchmark for studying stars and planetary systems at the end of the primordial disc phase. Gaia DR3 astrometry and photometry, combined with ground-based observations and more sophisticated stellar models, permit a systematic re-evaluation of BPMG membership and age. We combined Gaia astrometry with previously published radial velocities to evaluate moving group membership in a Bayesian framework. To minimize the effect of unresolved stellar multiplicity on age estimates, we identified and excluded multistar systems using Gaia astrometry, ground-based adaptive optics imaging, and multi-epoch radial velocities, as well as literature identifications. We estimated age using isochrone and lithium-depletion-boundary fitting with models that account for the effect of magnetic activity and spots on young, rapidly rotating stars. We find that age estimates are highly model-dependent; Dartmouth magnetic models with ages of 23 ± 8 and 33$^{+9}_{-11}$ Myr provide best fits to the lithium depletion boundary and Gaia MG versus BP–RP colour–magnitude diagram, respectively, whereas a Dartmouth standard model with an age of 11$^{+4}_{-3}$ Myr provides a best fit to the 2-Micron All-Sky Survey-Gaia$M_{K_S}$ versus BP–RP colour–magnitude diagram.
The binary star NY Hya is a bright, detached, double-lined eclipsing system with an orbital period of just under five days with two components each nearly identical to the Sun and located in the solar neighbourhood. The objective of this study is to test and confront various stellar evolution models for solar-type stars based on accurate measurements of stellar mass and radius. We present new ground-based spectroscopic and photometric as well as high-precision space-based photometric and astrometric data from which we derive orbital as well as physical properties of the components via the method of least-squares minimisation based on a standard binary model valid for two detached components. Classic statistical techniques were invoked to test the significance of model parameters. Additional empirical evidence was compiled from the public domain; the derived system properties were compared with archival broad-band photometry data enabling a measurement of the system's spectral energy distribution that allowed an independent estimate of stellar properties. We also utilised semi-empirical calibration methods to derive atmospheric properties from Strömgren photometry and related colour indices. Data was used to confront the observed physical properties with classic and magnetic stellar evolution models.
We report the discovery and characterization of a nearby (∼85 pc), older (27 ± 3 Myr), distributed stellar population near Lower Centaurus Crux (LCC), initially identified by searching for stars comoving with a candidate transiting planet from TESS (HD 109833; TOI 1097). We determine the association membership using Gaia kinematics, color–magnitude information, and rotation periods of candidate members. We measure its age using isochrones, gyrochronology, and Li depletion. While the association is near known populations of LCC, we find that it is older than any previously found LCC subgroup (10–16 Myr), and distinct in both position and velocity. In addition to the candidate planets around HD 109833, the association contains four directly imaged planetary-mass companions around three stars, YSES-1, YSES-2, and HD 95086, all of which were previously assigned membership in the younger LCC. Using the Notch pipeline, we identify a second candidate transiting planet around HD 109833. We use a suite of ground-based follow-up observations to validate the two transit signals as planetary in nature. HD 109833 b and c join the small but growing population of <100 Myr transiting planets from TESS. HD 109833 has a rotation period and Li abundance indicative of a young age (≲100 Myr), but a position and velocity on the outskirts of the new population, lower Li levels than similar members, and a color–magnitude diagram position below model predictions for 27 Myr. So, we cannot reject the possibility that HD 109833 is a young field star coincidentally nearby the population.
Low-mass stars spin-down over time due to angular momentum loss from magnetized stellar winds. This spin-down is interrupted by an epoch of stalling where stars halt their spin-down for 0.1–2 Gyr. Duration of the stalling epoch is mass-dependent, with higher mass stars spending less time stalled than the lowest mass stars. We demonstrate that the mass-dependent age at which stars resume spinning down correlates with the age where low-mass stars are predicted to achieve their minimum core temperature on the main sequence. This suggests that the stalling epoch is commensurate with a phase in low-mass main-sequence star evolution where stars develop a small convective core and their convective envelopes are growing in mass.
The demographics of young exoplanets can shed light on their formation and evolution processes. Exoplanet properties are derived from the properties of their host stars. As such, it is important to accurately characterize the host stars since any systematic biases in their derivation can negatively impact the derivation of planetary properties. Here we present a uniform catalog of photometrically derived stellar effective temperatures, luminosities, radii, and masses for 4865 young (<1 Gyr) stars in 31 nearby clusters and moving groups within 200 pc. We compared our photometrically derived properties to a subset of those derived from spectra and found them to be in good agreement. We also investigated the effect of stellar properties on the detection efficiency of transiting short-period young planets with TESS as calculated in Fernandes et al. (2022) and found an overall increase in the detection efficiency when the new photometrically derived properties were taken into account. Most notably, there is a 1.5 × increase in the detection efficiencies for sub-Neptunes/Neptunes (1.8–6 R ⊕ ) implying that, for our sample of young stars, better characterization of host star properties can lead to the recovery of more small transiting planets. Our homogeneously derived catalog of updated stellar properties, along with a larger unbiased stellar sample and more detections of young planets, will be a crucial input to the accurate estimation of the occurrence rates of young short-period planets.
Age is a stellar parameter that is both fundamental and difficult to determine. Among middle-aged M dwarfs, the most prolific hosts of close-in and detectable exoplanets, gyrochronology is the most promising method to assign ages, but requires calibration by rotation-temperature sequences (gyrochrones) in clusters of known ages. We curated a catalog of 249 late K- and M-type ((T_eff=3200-4200K) exoplanet host stars with established rotation periods, and applied empirical, temperature-dependent rotation-age relations based on relevant published gyrochrones, including one derived from observations of the 4 Gyr-old open cluster M67. We estimated ages for 227 of these stars, and upper limits for 8 others, excluding 14 which have too rapidly rotating or are otherwise outside the valid parameter range of our gyrochronology. We estimated uncertainties based on observed scatter in rotation periods in young clusters, error in the gyrochrones, and uncertainties in temperature and non-solar metallicity. For those stars with measured metallicities, we provide but do not incorporate a correction for the effects of deviation from solar metallicity. The age distribution of our sample declines to near zero at 10 Gyr, the age of the Galactic disk, with the handful of outliers explainable by large uncertainties. Continued addition or extension of cluster rotation sequences to more thoroughly calibrate the gyrochronology in time and temperature space, more precise and robust measurement of rotation periods, and more accurate stellar parameter measurements will enable continued improvements in the age estimates of these important exoplanet host stars.
Studies of T Tauri disks inform planet formation theory; observations of variability due to occultation by circumstellar dust are a useful probe of unresolved, planet-forming inner disks, especially around faint M dwarf stars. We report observations of 2M0632, an M dwarf member of the Carina young moving group that was observed by TESS over two one-year intervals. The combined light curve contains >300 dimming events, each lasting a few hours, and as deep as 40% (0.55 magnitudes). These stochastic events are correlated with a distinct, stable 1.86-day periodic signal that could be stellar rotation. Concurrent ground-based, multi-band photometry show reddening consistent with ISM-like dust. The star's excess emission in the infrared and emission lines in optical and infrared spectra, reveal a T Tauri-like accretion disk around the star. We confirm membership of 2M0632 in the Carina group by a Bayesian analysis of its Galactic space motion and position. We combine stellar evolution models with Gaia photometry and constraints on Teff, luminosity, and the absence of detectable lithium in the photosphere to constrain the age of the group and 2M0632 to 40-60 Myr, consistent with earlier estimates. 2M0632 joins a handful of long-lived disks which challenge the canon that disk lifetimes are <10 Myr. All known examples surround M dwarfs, suggesting that lower X-ray/UV irradiation and slower photoevaporation by these stars can dramatically affect disk evolution. The multi-planet systems spawned by long-lived disks probably experienced significant orbital damping and migration into close-in, resonant orbits, and perhaps represented by the TRAPPIST-1 system.
Main-sequence stars with convective envelopes often appear larger and cooler than predicted by standard models of stellar evolution for their measured masses. This is believed to be caused by stellar activity. In a recent study, accurate measurements were published for the K-type components of the 1.62-day detached eclipsing binary EPIC 219511354, showing the radii and temperatures for both stars to be affected by these discrepancies. This is a rare example of a system in which the age and chemical composition are known, by virtue of being a member of the well-studied open cluster Ruprecht 147 (age~3 Gyr, [Fe/H] = +0.10). Here, we report a detailed study of this system with nonstandard models incorporating magnetic inhibition of convection. We show that these calculations are able to reproduce the observations largely within their uncertainties, providing robust estimates of the strength of the magnetic fields on both stars: 1600 ± 130 G and 1830 ± 150 G for the primary and secondary, respectively. Empirical estimates of the magnetic field strengths based on the measured X-ray luminosity of the system are roughly consistent with these predictions, supporting this mechanism as a possible explanation for the radius and temperature discrepancies.
The Gaia M-dwarf gap is a significant underdensity of stars observed near M-G = 10.2 in a color-magnitude diagram for stars within 200 pc of the Sun. It has been proposed that the gap is the manifestation of structural instabilities within stellar interiors due to nonequilibrium He-3 fusion prior to some stars becoming fully convective. To test this hypothesis, we use Dartmouth stellar evolution models, MARCS model atmospheres, and simple stellar population synthesis to create synthetic M-G-(G(BP) - G(RP)) color-magnitude diagrams. We confirm that the proposed He-3 instability is responsible for the appearance of the M-dwarf gap. Our synthetic gap shows qualitatively similar features to the observed gap, including its vertical extent in M-G, its slope in the color-magnitude diagram, and its relative prominence at bluer colors as compared to redder colors. Furthermore, corresponding overdensities of stars above the gap are reproduced by the models. While qualitatively similar, the synthetic gap is approximately 0.2 magnitudes bluer, and when this color offset is accounted for, it is 0.16 magnitudes brighter than the observed gap. Our results reveal that the Gaia M-dwarf gap is sensitive to conditions within cores of M-dwarf stars, making the gap a powerful tool for testing the physics of M-dwarf stars and potentially using M dwarfs to understand the local star formation history.
ABSTRACT M dwarf stars are high-priority targets for searches for Earth-size and potentially Earth-like planets, but their planetary systems may form and evolve in very different circumstellar environments than those of solar-type stars. To explore the evolution of these systems, we obtained transit spectroscopy and photometry of the Neptune-size planet orbiting the ≈650-Myr-old Hyades M dwarf K2-25. An analysis of the variation in spectral line shape induced by the Doppler ‘shadow’ of the planet indicates that the planet’s orbit is closely aligned with the stellar equator ($\lambda =-1.7_{-3.7}^{+5.8}$ deg), and that an eccentric orbit found by previous work could arise from perturbations by another planet on a coplanar orbit. We detect no significant variation in the depth of the He i line at 1083 nm during transit. A model of atmospheric escape as an isothermal Parker wind with a solar composition shows that this non-detection is not constraining compared to escape rate predictions of ∼0.1 M⊕ Gyr−1; at such rates, at least several Gyr are required for a Neptune-like planet to evolve into a rocky super-Earth.
The Hertzsprung-Russell diagram (HRD) is fully examined using the Fourier analysis. This work shows more stars are above the gap in the lower main sequence than below it, and this implies that stars spend more time above the gap while they undergo variability associated with the He-3 instability. The enhanced HRD also shows the width of the gap is not linear and depends on the G(BP) -G(RP) color up until G(BP) - G(RP)=2.7. Beyond this color limit, the gap is hardly seen. Besides, a new low density region is revealed for the first time centered at M-G approximate to 10.7 and G(BP) - G(RP) approximate to 2.8, which is below the lower right corner of the gap. This work also shows that the main sequence appears to have fine stripes where stellar densities are relatively low or high compared to their adjacent regions on the main sequence. These stripes can be seen throughout the main sequence of stars redder than G(BP) - G(RP) = 0.8 and are not limited to any specific color or spectral type. Slopes of these features are different from the main sequence, but are pretty consistent throughout the main sequence, with a few exceptions. We are perplexed by these new features, but the complexities of stellar atmospheric features and opacities of dwarfs may have caused these patterns.
The mass–luminosity relation for late-type stars has long been a critical tool for estimating stellar masses. However, there is growing need for both a higher-precision relation and a better understanding of systematic effects (e.g., metallicity). Here we present an empirical relationship between and M* spanning 0.075 M⊙ < M* < 0.70 M⊙. The relation is derived from 62 nearby binaries, whose orbits we determine using a combination of Keck/NIRC2 imaging, archival adaptive optics data, and literature astrometry. From their orbital parameters, we determine the total mass of each system, with a precision better than 1% in the best cases. We use these total masses, in combination with resolved KS magnitudes and system parallaxes, to calibrate the –M* relation. The resulting posteriors can be used to determine masses of single stars with a precision of 2%–3%, which we confirm by testing the relation on stars with individual dynamical masses from the literature. The precision is limited by scatter around the best-fit relation beyond measured M* uncertainties, perhaps driven by intrinsic variation in the –M* relation or underestimated uncertainties in the input parallaxes. We find that the effect of [Fe/H] on the –M* relation is likely negligible for metallicities in the solar neighborhood (0.0% ± 2.2% change in mass per dex change in [Fe/H]). This weak effect is consistent with predictions from the Dartmouth Stellar Evolution Database, but inconsistent with those from MESA Isochrones and Stellar Tracks (at 5σ). A sample of binaries with a wider range of abundances will be required to discern the importance of metallicity in extreme populations (e.g., in the Galactic halo or thick disk).
The accuracy of masses of pre-main-sequence stars derived from their locations on the Hertzsprung–Russell diagram (HRD) can be tested by comparison with accurate and precise masses determined independently. We present 29 single stars in the Taurus star-forming region (SFR) and 3 in the Ophiuchus SFR with masses measured dynamically to a precision of at least 10%. Our results include 9 updated mass determinations and 3 that have not had their dynamical masses published before. This list of stars with fundamental, dynamical masses, Mdyn, is drawn from a larger list of 39 targets in the Taurus SFR and 6 in the Ophiuchus SFR. Placing the stars with accurate and precise dynamical masses on HRDs that do not include internal magnetic fields underestimates the mass compared to Mdyn by about 30%. Placing them on an HRD that does include magnetic fields yields mass estimates in much better agreement with Mdyn, with an average difference between Mdyn and the estimated track mass of 0.01 ± 0.02 M⊙. The ages of the stars, 3–10 MY on tracks that include magnetic fields, is older than the 1–3 MY indicated by the nonmagnetic models. The older ages of T Tauri stars predicted by the magnetic models increase the time available for evolution of their disks and formation of the giant gas exoplanets. The agreement between our Mdyn values and the masses on the magnetic field tracks provides indirect support for these older ages.
Young exoplanets are snapshots of the planetary evolution process. Planets that orbit stars in young associations are particularly important because the age of the planetary system is well constrained. We present the discovery of a transiting planet larger than Neptune but smaller than Saturn in the 45 Myr Tucana–Horologium young moving group. The host star is a visual binary, and our follow-up observations demonstrate that the planet orbits the G6V primary component, DS Tuc A (HD 222259A, TIC 410214986). We first identified transits using photometry from the Transiting Exoplanet Survey Satellite ( TESS ; alerted as TOI 200.01). We validated the planet and improved the stellar parameters using a suite of new and archival data, including spectra from Southern Astrophysical Research/Goodman, South African Extremely Large Telescope/High Resolution Spectrograph and Las Cumbres Observatories/Network of Robotic Echelle Spectrographs; transit photometry from Spitzer ; and deep adaptive optics imaging from Gemini/Gemini Planet Imager. No additional stellar or planetary signals are seen in the data. We measured the planetary parameters by simultaneously modeling the photometry with a transit model and a Gaussian process to account for stellar variability. We determined that the planetary radius is 5.70 ± 0.17 R ⊕ and that the orbital period is 8.1 days. The inclination angles of the host star’s spin axis, the planet’s orbital axis, and the visual binary’s orbital axis are aligned within 15° to within the uncertainties of the relevant data. DS Tuc Ab is bright enough ( V = 8.5) for detailed characterization using radial velocities and transmission spectroscopy.