Twinkle Twinkle Little Star, Roman Sees Where You Are: Predicting Exoplanet Transit Yields in the Rosette Nebula with the Nancy Grace Roman Space Telescope

Andrew Vanderburg; Johanna M. Vos; Mary Anne Limbach; Melinda Soares-Furtado; Nishanth Ramanujam; Ritvik Sai Narayan

arxiv: 2510.23708 · v2 · submitted 2025-10-27 · 🌌 astro-ph.EP · astro-ph.GA· astro-ph.SR

Twinkle Twinkle Little Star, Roman Sees Where You Are: Predicting Exoplanet Transit Yields in the Rosette Nebula with the Nancy Grace Roman Space Telescope

Ritvik Sai Narayan , Melinda Soares-Furtado , Mary Anne Limbach , Nishanth Ramanujam , Andrew Vanderburg , Johanna M. Vos This is my paper

Pith reviewed 2026-05-18 02:51 UTC · model grok-4.3

classification 🌌 astro-ph.EP astro-ph.GAastro-ph.SR

keywords exoplanetstransitsRosette NebulaRoman Space Telescopeyoung planetsplanet migrationstar forming regions

0 comments

The pith

The Roman Space Telescope could detect 33 young transiting exoplanets in a month-long survey of the Rosette Nebula.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper evaluates how many young exoplanets transiting their host stars could be found by pointing the Roman Space Telescope at the Rosette Nebula for a few weeks. Young stars are variable and often in crowded, dusty areas, making planet detection hard, so the authors run detailed simulations to predict yields while including those effects. A sympathetic reader would care because these planets are at an age when they are still changing size and position due to migration, and we have almost no examples to study. The result is that 30 or so detections are plausible, mostly small planets on short orbits, opening a window on early planet evolution.

Core claim

Using Monte Carlo simulations of transit injections and recoveries with Roman's sensitivity, the authors predict 33 plus or minus 9 transiting young exoplanets in a 30-day survey of the Rosette Nebula and 29 plus or minus 8 in 14 days. Detections would mostly be 1 to 2 Earth-radius planets with periods of 8 days or less. This sample would greatly increase the number of known transiting planets less than 20 million years old, in the phase where radii are inflated and migration is likely still happening.

What carries the argument

Monte Carlo injection-recovery simulations that incorporate the Roman Exposure Time Calculator, nebular extinction, and models of youth-driven stellar variability to determine detectable planets.

If this is right

The main gain from a longer baseline is better detection of longer-period planets around FGK-type stars.
Planets around M dwarfs are detected nearly as well in two weeks as in a month.
Detections are expected to be dominated by super-Earths and sub-Neptunes smaller than 2 Earth radii on orbits shorter than 8 days.
This would provide targets for further study with JWST and other observatories to understand early planetary evolution.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

If the survey succeeds it would test whether close-in planets form and migrate quickly in the first 10 million years.
The approach could be extended to other star-forming regions with Roman or future telescopes.
Improved measurements of young star variability from other surveys would make these yield predictions more accurate.

Load-bearing premise

Planet occurrence rates around young stars and the detailed properties of their variability are assumed based on sparse existing observations or scaled from older populations.

What would settle it

Performing the survey and recovering a number of transiting planets much different from the predicted 33 would show that the assumed occurrence rates or variability amplitudes are inaccurate.

read the original abstract

Young stars host only a small fraction of the known exoplanet population because their photometric variability, magnetic activity, and frequent placement in dense, poorly-resolved regions hamper exoplanet detections. Yet, measuring planets at these ages is crucial since these phases are when dynamical processes that drive planetary migration are most active. We assess the expected yield of a hypothetical Nancy Grace Roman Space Telescope transit survey of the Rosette Nebula, a ${\sim}10$ Myr star-forming region with a dense and diverse stellar population. Using the Roman Exposure Time Calculator to quantify sensitivity to Rosette members, we establish detection thresholds for companions and evaluate yields via Monte Carlo injection-recovery simulations, accounting for nebular extinction and youth-driven stellar variability. We predict the detection of $33 \pm 9$ young transiting exoplanets orbiting stellar hosts in a month-long survey, and $29 \pm 8$ in a two-week survey. The extended baseline primarily improves sensitivity to longer-period planets orbiting FGK stars, while most M dwarf detections are well-sampled within two weeks. Irrespective of the temporal baseline, transit detections are dominated by of 1-2 $R_\oplus$ super-Earths and sub-Neptunes with $P\lesssim8$ days. Such a sample would substantially expand the census of only three detected transiting planets younger than 20 Myr around stars less massive than the Sun, probing an age regime in which planetary radii remain inflated, the stability of close-in orbits is uncertain, and planetary migration may still be ongoing. This survey offers a path to constrain early planetary evolution and establish prime follow-up targets for JWST, Rubin, and the Habitable Worlds Observatory.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper forecasts 33 young transiting planets from a month-long Roman survey of the Rosette but the quoted errors probably understate the effect of uncertain occurrence rates and youth variability.

read the letter

The main thing to know is that this work predicts a Roman transit survey of the Rosette Nebula could detect 33 young planets in a month or 29 in two weeks. That would be a real increase over the three known transiting cases under 20 Myr around solar-mass or smaller stars, and it targets the age range where migration and radius inflation are still active. The calculation uses the published Roman Exposure Time Calculator, adds nebular extinction, and runs Monte Carlo injection-recovery that includes stellar variability from youth. Most detections come out as 1-2 Earth-radius planets with periods under eight days, and the extra baseline mainly helps longer-period cases around FGK stars while M-dwarf detections are already well covered in two weeks. The approach is a straightforward forward model applied to a new target rather than a new method, and the breakdown by host type and period is clear enough to be useful for planning. The soft spot sits in the error bars. The ±9 and ±8 appear to come from the simulation scatter, but the planet occurrence rates and the amplitude or timescale of variability are drawn from limited samples or extrapolations from older stars. At 10 Myr, ongoing dynamical evolution could alter occurrence, and spot or flare activity might exceed the modeled levels. A factor-of-two shift in either input would move the central yields outside the reported uncertainties, yet the abstract does not detail the exact priors or number of realizations. This leaves the robustness harder to judge without the full methods. The paper is aimed at people planning Roman observations or studying early exoplanet evolution. It gives quantitative numbers for one specific young region that were not in the literature before, so readers who need yield estimates for survey design or who want to compare against other clusters will get direct value. I would send it for peer review. The core simulation setup is solid and the topic fills a practical gap, even if the input assumptions need closer scrutiny on the assumptions and sensitivity tests.

Referee Report

2 major / 2 minor

Summary. The manuscript uses the Roman Exposure Time Calculator and Monte Carlo injection-recovery simulations to forecast the number of young transiting exoplanets detectable in a hypothetical survey of the Rosette Nebula (~10 Myr star-forming region). Accounting for nebular extinction and youth-driven stellar variability, the authors predict yields of 33 ± 9 planets in a one-month survey and 29 ± 8 in a two-week survey, dominated by 1–2 R_⊕ super-Earths and sub-Neptunes with P ≲ 8 days. The work emphasizes the scientific value for studying early planetary migration and radii inflation, and for identifying follow-up targets for JWST and other facilities.

Significance. If the central yield predictions hold under reasonable input assumptions, the result would be significant for planning Roman observations of young clusters and for expanding the tiny existing sample of transiting planets younger than 20 Myr. The quantitative forecasts, grounded in the published Roman ETC, provide a concrete basis for assessing survey feasibility and highlight the dominance of short-period small planets even in a dense, variable young population.

major comments (2)

[§4] §4 (Monte Carlo simulation description): The reported uncertainties (±9 and ±8) appear to reflect only Poisson sampling from the injection-recovery runs. No sensitivity analysis is shown for systematic variations in the two key free parameters (planet occurrence rate per star and stellar variability amplitude/timescale). A factor-of-two change in occurrence rate or 50% increase in variability amplitude—both plausible for ~10 Myr stars—would scale the headline yields proportionally and likely move them outside the quoted error bars.
[§3.1] §3.1 (occurrence rate and variability priors): The manuscript adopts occurrence rates and variability prescriptions from limited empirical samples or extrapolations calibrated on older stars. For the Rosette's ~10 Myr population, ongoing migration and enhanced spot/flare activity could alter these inputs; the paper should quantify the resulting range in predicted yields rather than treating the inputs as fixed.

minor comments (2)

[Figure 3] Figure 3 (or equivalent yield histogram): The caption should explicitly state the number of Monte Carlo realizations performed and whether the plotted error bars include only counting statistics.
[Abstract] Abstract and §5: The phrase 'dominated by of 1-2 R_⊕' contains a grammatical error; rephrase for clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript's significance and for the constructive comments on the treatment of uncertainties and input parameters. We have revised the paper to incorporate sensitivity analyses and expanded discussion of the priors, as detailed below.

read point-by-point responses

Referee: [§4] §4 (Monte Carlo simulation description): The reported uncertainties (±9 and ±8) appear to reflect only Poisson sampling from the injection-recovery runs. No sensitivity analysis is shown for systematic variations in the two key free parameters (planet occurrence rate per star and stellar variability amplitude/timescale). A factor-of-two change in occurrence rate or 50% increase in variability amplitude—both plausible for ~10 Myr stars—would scale the headline yields proportionally and likely move them outside the quoted error bars.

Authors: We agree that the quoted uncertainties capture primarily the statistical variation across our Monte Carlo realizations (driven by the finite number of injected planets and the random sampling of detection outcomes). To address this limitation, we have added a dedicated sensitivity analysis in the revised §4. We reran the injection-recovery suite with occurrence rates scaled by factors of 0.5 and 2.0 relative to the fiducial value, and with stellar variability amplitudes increased by 50%. The results show that yields scale nearly linearly with occurrence rate (as expected from the simulation design) while a 50% variability increase reduces overall detection efficiency by ~20–30%, primarily affecting the smallest planets. We have inserted a new table and accompanying text summarizing the resulting yield ranges (e.g., 17–66 planets for the one-month survey under these variations) and have updated the abstract and conclusions to present these as plausible systematic bounds rather than relying solely on the Poisson errors. revision: yes
Referee: [§3.1] §3.1 (occurrence rate and variability priors): The manuscript adopts occurrence rates and variability prescriptions from limited empirical samples or extrapolations calibrated on older stars. For the Rosette's ~10 Myr population, ongoing migration and enhanced spot/flare activity could alter these inputs; the paper should quantify the resulting range in predicted yields rather than treating the inputs as fixed.

Authors: The fiducial occurrence rates are taken from the most recent young-cluster transit surveys (primarily <20 Myr samples), and the variability model is calibrated on TESS light curves of pre-main-sequence stars in comparable regions. We acknowledge that these are extrapolations and that ongoing migration or heightened magnetic activity at exactly ~10 Myr could shift the inputs. In the revision we have expanded §3.1 to explicitly discuss the sample-size limitations and possible youth-specific effects. We then link directly to the new sensitivity results added in §4, which now quantify how plausible changes in these parameters propagate to the final yield predictions. This provides the requested range without overstating the precision of the current empirical constraints. revision: yes

Circularity Check

0 steps flagged

No significant circularity in yield predictions

full rationale

The paper derives its headline yield numbers (33 ± 9 and 29 ± 8) exclusively via forward Monte Carlo injection-recovery simulations that take as inputs the Roman Exposure Time Calculator sensitivity curves, literature-based planet occurrence rates, empirical or extrapolated youth variability amplitudes/timescales, and an adopted extinction map. None of these inputs are fitted to Rosette data within the paper, nor do any equations or self-citations reduce the final counts to a tautological re-expression of the same quantities. The derivation chain remains open to external benchmarks and is therefore self-contained.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The yield prediction depends on several inputs whose values are not derived within the paper: planet occurrence rates around young stars, the statistical properties of photometric variability for pre-main-sequence stars, and the spatial distribution of extinction within the Rosette. These are treated as known or extrapolated quantities.

free parameters (2)

Planet occurrence rate per star
Used to set the number of injected planets; value not stated in abstract but must be chosen to produce the reported yields.
Stellar variability amplitude and timescale
Youth-driven variability is explicitly included; its functional form and parameters determine detection completeness for short-period planets.

axioms (2)

domain assumption The Roman Exposure Time Calculator accurately predicts photometric precision for Rosette members after nebular extinction correction.
Invoked when setting detection thresholds; no independent verification of the ETC for this crowded, extincted field is described.
domain assumption Current empirical occurrence rates and variability statistics for young stars can be extrapolated to the Rosette population without large systematic error.
Central to scaling the Monte Carlo results to absolute yields.

pith-pipeline@v0.9.0 · 5885 in / 1589 out tokens · 32264 ms · 2026-05-18T02:51:22.237343+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We develop a Monte Carlo injection-recovery framework that creates a sample population and models our detectability of transiting exoplanets... using spectral-type-dependent occurrence rates... age-dependent radius inflation... 7σ detection threshold.
IndisputableMonolith/Foundation/DimensionForcing.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

We predict the detection of 33 ± 9 young transiting exoplanets... 29 ± 8 in a two-week survey.

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.