{"id":"38e948f4-c7e5-4cbc-a4ed-adcdaa86c690","arxiv_id":"2509.02698","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Red main-sequence stars in the young LMC clusters NGC 1866 and NGC 1856 rotate at mean equatorial velocities near 300 km/s, about 70-80% of critical, linking rotation to the split main sequences and high Be star fractions.","lead":"Astronomers measured how fast stars spin in two young clusters in the Large Magellanic Cloud, NGC 1866 and NGC 1856, using MUSE spectroscopy. They found that the redder main-sequence stars rotate at about 300 km/s, close to 70-80% of their break-up speed, which can explain the large number of Be stars seen there.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline Veq~300 km/s rests on a newly implemented direct V sin i fit that is never validated end-to-end; an unquantified LSF/template bias would shift the central claim directly.","rationale":"The reader's weakest-assumption analysis focused on isotropy of spin axes and the single-Gaussian Veq assumption. Those are real limitations, but they are acknowledged in the paper and mostly affect the tail-based Be-star statistics, not the bulk Veq measurement. The more upstream, least-secure link is the V sin i calibration itself: the paper introduces a direct rotational-broadening fit in Spexxy and gives no independent evidence that this fit is accurate at MUSE resolution. Because the entire quantitative result is a velocity scale, a systematic error in V sin i would propagate directly into the headline Veq~300 km/s and into the SYCLIST-derived V/Vcrit and Be-star threshold arguments. I am not claiming the measurements are wrong; rather, the supporting validation is missing. An injection test would settle the question quickly. If that test passes, the current ACCEPT verdict is justified; until then, I would treat the acceptance as conditional on this calibration check.","tokens_in":28740,"tokens_out":18402,"duration_ms":194057,"concrete_test":"Run an end-to-end recovery test: embed rotationally broadened synthetic spectra of 2.7 Msun, Z=0.006 B-type stars with known V sin i = 0, 50, 100, 200, 300, and 400 km/s into the actual MUSE cubes (or into noise-matched synthetic cubes with the same LSF, sampling, and S/N), then extract and fit them with the same PampelMuse/Spexxy pipeline used in Secs. 2.4-2.5. Compare recovered versus injected V sin i per star. If the recovered values are biased by more than ~5% or ~15 km/s, the reported ⟨Veq⟩ and V/Vcrit values need revision; if unbiased within these thresholds, the concern is resolved. Report the resulting systematic error from LSF and template choice alongside Table 1.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim (Secs. 4.2 and 5.1) depends entirely on the V sin i values produced by Spexxy using a newly implemented rotational-broadening recipe (Sec. 2.5). The paper reports no recovery test, no comparison with independent V sin i measurements, and no systematic error budget for this step. At MUSE resolution (LSF FWHM ~2.5 Å), the rotational signal at 300 km/s is only ~4-5 Å at the Balmer lines, so small errors in the adopted LSF or in the assumed limb-darkening/template physics translate directly into a shifted V sin i scale. The template grid is spherical and non-rotating, and log g is fixed from isochrones; gravity darkening is not modeled in the fitting, although the stars are fast rotators near 70-80% critical. A systematic offset of ~10% in V sin i would move the derived ⟨Veq⟩ away from ~300 km/s and change the SYCLIST-based V/Vcrit ~0.7-0.8 and the 85%-critical Be-star threshold. The mock checks in Fig. 9 validate the inclination deprojection conditional on the V sin i values; they do not validate the V sin i values themselves.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper analyzes MUSE integral-field spectroscopy of two young massive LMC clusters, NGC 1866 (~200 Myr) and NGC 1856 (~300 Myr), to measure projected rotational velocities (V sin i) along their main sequences. The authors find a strong correlation between photometric color and V sin i, with red main-sequence stars having a mean equatorial velocity of ~300 km/s after statistical inclination correction, corresponding to 70-80% of the SYCLIST-model critical velocity. They also identify large populations of Be stars, with maximum fractions of 85% in NGC 1866 and 78% in NGC 1856, and argue that if ~85% of critical rotation triggers disk formation, the single-star evolutionary channel can explain the majority of these Be stars. The paper discusses implications for the origin of the split main sequences, the Be-star phenomenon, and the physical scenarios proposed for the stellar spin bimodality in young clusters.","tokens_in":28990,"tokens_out":2955,"duration_ms":29523,"significance":"If correct, the paper delivers the first direct spectroscopic measurement of the equatorial velocities of the red main-sequence population in these two clusters, providing strong evidence that stellar rotation, rather than age spreads, drives the split main sequence and extended main-sequence turn-off. The large samples, careful proper-motion membership selection, and the explicit likelihood model with MCMC uncertainties and mock-data checks (Table 1, Fig. 9) are strengths. The comparison of Be-star fractions across four clusters is valuable, and the argument that photometric searches underestimate Be fractions in older clusters is important. However, the central measurement rests on a newly implemented V sin i fitting recipe that is not validated end-to-end, making the headline Veq values vulnerable to systematic errors; this is the main factor limiting the current assessment.","major_comments":[{"comment":"The central claim of mean equatorial velocities Veq ~300 km/s depends entirely on the V sin i values produced by the newly implemented rotational-broadening recipe in Spexxy. The paper reports no recovery test, no comparison with independent V sin i measurements, and no systematic error budget for the adopted MUSE LSF, the spherical non-rotating template grid, the fixed log g, or the neglect of gravity darkening. At MUSE resolution the rotational signal at ~300 km/s is only a few Angstroms, so a small LSF or template mismatch would shift the V sin i scale and directly change the derived Veq, the SYCLIST-based V/Vcrit, and the inferred Be-star threshold. The mock checks in Fig. 9 validate only the inclination deprojection conditional on the V sin i values, not the V sin i measurements themselves. I request an end-to-end validation: inject artificial spectra with known V sin i through the same extraction and fitting pipeline, compare measured V sin i against high-resolution literature values for stars in common, and quantify the sensitivity of the results to the assumed LSF and to the choice of limb-darkening/template treatment.","section":"Sec. 2.5 and Secs. 4.2/5.1"},{"comment":"The inference of Veq assumes (i) isotropically distributed spin axes and (ii) a single Gaussian distribution of intrinsic equatorial velocities. The authors themselves note that the model slightly underestimates the number of NGC 1866 stars with V sin i below ~100 km/s and that the NGC 1856 peak is broader than predicted, which they attribute to a possible second population and to a genuinely broader intrinsic distribution. Because the fitted mean Veq is the load-bearing quantity for the paper's conclusions, the possible presence of a second population (e.g., photometric binaries migrating from the blue main sequence) should be quantitatively tested. I recommend fitting a two-component model, excluding the slow tail and refitting, or otherwise demonstrating that the derived Veq is robust to these assumptions. The current single-Gaussian fit is a simplified model whose systematic uncertainty is not captured in the quoted errors of Table 1.","section":"Sec. 4.2, Eq. (2)"},{"comment":"The conversion of the measured Veq into fractions of critical rotation is entirely model-dependent. The SYCLIST models are initialized at Ω/Ωcrit = 0.95 specifically to reproduce the measured Veq, and the paper itself notes that different model suites (MESA, PARSEC) give different fractional critical velocities, as illustrated by the NGC 1866 comparison with Gossage et al. (2019) and Wang et al. (2023a). The abstract's statement that the stars rotate at \"70-80% of the critical values predicted for such stars by current stellar models\" is therefore only one model's estimate, and this caveat should be carried through the title-level summary. I do not see this as circularity, since Veq is measured directly, but the fractional values and the inferred ~85% Be-onset threshold should be presented as SYCLIST-specific rather than as a robust empirical result.","section":"Sec. 4.3 and Sec. 5.4.4"}],"minor_comments":[{"comment":"The description of the new rotational-broadening implementation would benefit from specifying the functional form used (e.g., Gray 2005 formulation), the adopted limb-darkening coefficient, and whether any tests were performed on the numerical accuracy of the convolution.","section":"Sec. 2.5"},{"comment":"The caption states that the grey lines show the expected V sin i distributions from the model, but the text clarifies that these are the median values from 1000 mock samples. Please also show the scatter (e.g., shaded regions) in the figure or caption, since the current visual impression of model-data agreement is based only on the median.","section":"Fig. 9"},{"comment":"The shell-star fractions and derived disk half-opening angles are interesting, but the statement that the fractions agree with the 22.8% field value (Hanuschik 1996) while also finding a possible age/mass trend should be briefly rephrased to avoid an apparent contradiction; the small-number statistics may not support a trend.","section":"Sec. 3.3"},{"comment":"The text notes that the V sin i values of Be stars are not considered reliable due to emission contamination, yet Fig. 8 shows these values and the narrative in Sec. 4.1 draws attention to them exceeding critical velocities. Please either include a clear caveat in the figure caption or move the Be-star V sin i discussion fully to the section where the unreliability is stated.","section":"Sec. 5.4.3"},{"comment":"The abstract quotes a Be-star fraction of '>~50%' at the turn-off while the text reports maxima of 85% and 78%. Using the same values in the abstract would avoid a perceived inconsistency.","section":"Abstract and Sec. 6"},{"comment":"Typographical issues: 'aknowledges' should be 'acknowledges', and 'Caii' should be 'Ca ii' in two places.","section":"Acknowledgements"}],"recommendation":"major_revision","confidential_remarks":"The paper presents a compelling and well-written observational study, and I believe the central measurement is likely correct. However, the lack of end-to-end validation of the newly implemented V sin i fitting is a genuine load-bearing gap that should be addressed before publication. The authors may already have internal tests that were not included in the manuscript; if so, adding them would substantially strengthen the paper. I also encourage the editor to ensure that the model-dependence of the V/Vcrit values is clearly flagged in the abstract and conclusions, as these numbers are being quoted as empirical results."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nYou should know two things about this paper. First, it delivers the first spectroscopic rotation maps of NGC 1866 and NGC 1856, and the headline is that the red main-sequence stars in both clusters have mean equatorial velocities around 300 km/s, about 70-80% of the SYCLIST critical velocity. Second, the measurement pipeline has an untested step: the authors implemented a new direct V sin i fitting routine in Spexxy and never validate it against independent measurements or synthetic spectra.\n\nWhat's genuinely good: the sample selection is careful (proper-motion membership, quality cuts), the likelihood model for deprojecting V sin i to Veq is transparent and checked with mock drawings, and the Be star census is a real step up from narrowband photometry—it suggests photometric studies missed a large fraction of Be stars in older clusters. They are also appropriately cautious about model dependence when converting to V/Vcrit, and they explicitly refuse to interpret V sin i for Be stars because of emission contamination. The qualitative story—slow rotators on the blue main sequence, fast rotators on the red—is consistent with earlier work, so the new quantitative layer is what matters.\n\nThe soft spot is real: the direct V sin i fit is load-bearing, and nothing calibrates it. The templates are spherical non-rotating models, log g is fixed from isochrones, gravity darkening is not included in the fitting, and the MUSE LSF is adopted rather than measured. At 2.5 Å resolution, a 300 km/s rotational signal is only a few Å at the Balmer lines. A systematic 10% offset in the LSF or in the template physics moves the headline Veq away from 300 km/s and changes the 85%-critical threshold for the Be argument. The mock checks in Fig. 9 validate the inclination model, not the V sin i scale.\n\nThe single-Gaussian Veq distribution is a further simplification, acknowledged by the authors when they see a slow-rotator tail in NGC 1866. That isn't fatal, but it means the quoted mean and sigma are model-dependent.\n\nVerdict: this deserves a serious referee. The measurement is relevant and the analysis is mostly transparent, but the authors need to either validate the V sin i scale (comparison with literature V sin i for common stars, or a synthetic recovery test) or quantify the systematic uncertainty. As it stands, I'd treat the qualitative picture as solid and the 300 km/s number as provisional. The paper is for the cluster-rotation and Be-star community. I'd support peer review and would push for the validation before publication.","headline":"Careful, useful rotation measurements in two more LMC clusters, but the new V sin i fitting lacks end-to-end validation and that uncertainty sits right on the headline number.","tokens_in":29598,"tokens_out":3809,"would_cite":true,"duration_ms":36011,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In two young LMC clusters, red main-sequence stars spin at about 300 km/s, 70–80% of their break-up speed.","keywords":["stellar rotation","split main sequence","extended main-sequence turn-off","young massive clusters","Large Magellanic Cloud","Be stars","MUSE integral field spectroscopy","projected rotation velocity"],"falsifier":"Take high-resolution spectra of the red main-sequence stars in NGC 1866 in the blue, where rotation broadening is cleanly separable and no disk emission contaminates the lines, and reconstruct the distribution of true equatorial velocities without assuming a single Gaussian; if a noticeable slow-rotator population appears, the reported mean of about 300 km/s and the 70–80% critical fraction would be overestimated.","tokens_in":28490,"feed_emoji":"🔄","tokens_out":8361,"duration_ms":75207,"temperature":0.7,"pith_summary":"This paper sets out to test whether stellar rotation, rather than star-formation age spreads, produces the split main sequences and extended turn-offs seen in young massive clusters. Using MUSE spectroscopy of NGC 1866 and NGC 1856, it finds a tight colour–rotation link: stars on the red side of the split rotate fast, blue-side stars rotate slowly, and the red main-sequence population has a mean equatorial velocity of about $300\\,\\mathrm{km\\,s^{-1}}$ in both clusters, corresponding to 70–80% of the break-up speed predicted by current models. The paper also reports very high Be-star fractions at the turn-off (up to 85% in NGC 1866) and argues that single-star evolution can explain them if rotation at about 85% of critical is enough to form a decretion disk. If correct, rotation becomes the main driver of these photometric anomalies and the leading explanation for Be stars in young massive clusters.","feed_headline":"Red main-sequence stars spin at 300 km/s in two LMC clusters","feed_subtitle":"Rotation, not age spread, explains the split main sequences and may create most Be stars in these clusters.","key_machinery":"The measurement chain rests on MUSE integral-field spectroscopy combined with full-spectrum fitting that now fits rotational line broadening directly, yielding projected velocities $V\\sin i$ for thousands of cluster members. The statistical deprojection uses a model in which the intrinsic equatorial velocities form a single Gaussian and spin axes are isotropic, marginalising over inclination. The SYCLIST rotating stellar models translate the recovered $V_{\\rm eq}$ into a fraction of the critical (break-up) velocity and predict its evolution along the main sequence, which is how the paper converts $300\\,\\mathrm{km\\,s^{-1}}$ into the 70–80% claim and links it to the onset of Be disks.","core_discovery":"Using MUSE spectroscopy of more than 1,200 member stars in NGC 1866 and 1,300 in NGC 1856, the paper establishes a direct, monotonic link between photometric colour along the split main sequence and projected rotation velocity $V\\sin i$: blue-side stars are slow rotators (mostly below $100\\,\\mathrm{km\\,s^{-1}}$), while red-side stars cluster around $200\\,\\mathrm{km\\,s^{-1}}$ in projection. After statistically correcting for the unknown inclination of each spin axis, the intrinsic equatorial-velocity distribution of the red main-sequence population peaks at $\\langle V_{\\rm eq}\\rangle = 302.5^{+4.5}_{-4.7}\\,\\mathrm{km\\,s^{-1}}$ in NGC 1866 and $291.3^{+6.0}_{-6.0}\\,\\mathrm{km\\,s^{-1}}$ in NGC 1856, i.e. roughly $300\\,\\mathrm{km\\,s^{-1}}$ in both clusters. Comparing with SYCLIST rotating models, these values correspond to 70–80% of the critical (break-up) velocity for stars of about $2.7\\,M_\\odot$. The paper further reports that Be stars constitute up to 85% of the turn-off population in NGC 1866 and 78% in NGC 1856, and argues that this is the expected outcome if the fast rotators approach $\\sim$85% of critical rotation before leaving the main sequence, so that most Be stars in these clusters can form through single-star evolution rather than binary mass transfer.","pith_inferences":["Reading the four spectroscopically studied young LMC clusters together, the same fast-rotator population appears at roughly 300 km/s regardless of age; a natural extension would be to test whether clusters of different metallicities share the same equatorial velocity rather than the same fraction of critical rotation.","If the approximately 85%-critical threshold for disk formation holds, the Be duty cycle in these clusters must be near unity; a testable consequence is that few fast-rotating B stars should be found without disks near the turn-off.","The gradual, rather than sharply bimodal, colour–$V\\sin i$ relation favours spin distributions set by a continuum of initial conditions, such as disk-locking, over merger or binary scenarios that predict a clean split into two velocity groups.","The apparent age trend in spatial segregation of fast rotators, strongest in the youngest cluster and absent in the oldest, could be a dynamical mixing signature worth testing with cluster simulations.","A direct discriminator between the single-star and binary channels would be unbiased rotation measurements of the Be stars at blue wavelengths, since the emission lines that contaminate the current fits preclude comparing their spins with those of ordinary B stars."],"forward_implications":["If the rotation scenario is right, the split main sequences and extended turn-offs of young massive clusters are primarily spin artefacts, not evidence for tens of millions of years of star formation, so inferred cluster ages and formation histories need revision.","Fast rotators near the turn-off should routinely approach the critical velocity, turning a large fraction of upper main-sequence stars into Be stars; the observed 50–85% fractions become the expected outcome rather than an anomaly.","Photometric searches for Be stars miss a large share of the population in clusters older than about 100 Myr, because the H-alpha emission is weaker and partly shifted out of narrow-band filters; spectroscopic surveys are required for accurate Be fractions.","Stellar evolution model fits that ignore direct rotation measurements cannot uniquely determine spin distributions, so future isochrone fitting should incorporate $V\\sin i$ measurements.","Because the older cluster NGC 1856 has a broader equatorial-velocity distribution and sits closer to critical rotation, it shows Be stars at magnitudes where NGC 1866 does not, tying the Be phenomenon to the age-dependent drop in critical velocity."],"supporting_citations":[{"why":"Prior MUSE study of NGC 1850 that established the colour–$V\\sin i$ link and the analysis methods adapted here.","marker":"Kamann et al. (2023)"},{"why":"First spectroscopic detection of faster rotation among red main-sequence stars in NGC 1818, the baseline result this paper extends.","marker":"Marino et al. (2018)"},{"why":"Provides the SYCLIST rotating models used to convert measured equatorial velocities into fractions of critical rotation.","marker":"Georgy et al. (2013)"},{"why":"Disk-locking scenario for slow and fast rotators whose spatial segregation prediction is tested in this paper.","marker":"Bastian et al. (2020)"},{"why":"Model of Be-star production from single-star evolution used to compare predicted and observed Be fractions.","marker":"Hastings et al. (2020)"},{"why":"Upper limits on binary-channel Be fractions used to argue that binaries alone cannot explain the observed high fractions.","marker":"Hastings et al. (2021)"},{"why":"Rotating MESA model comparison showing model dependence of critical-velocity fractions and supporting moderate rotation rates.","marker":"Wang et al. (2023a)"},{"why":"Photometric Be fractions and split main-sequence catalogs that the spectroscopic results are compared against.","marker":"Milone et al. (2018)"},{"why":"Provides the HST photometry and proper motions used to define cluster members and verticalised main sequences.","marker":"Niederhofer et al. (2024)"}],"fun_headline_variants":["Fast rotators explain split main sequence in LMC clusters","300 km/s spin links colour to rotation in young clusters","Red main-sequence stars in LMC clusters rotate near break-up","Single-star fast rotation makes most Be stars in LMC clusters","Rotation, not age, drives split main sequence in two LMC clusters"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 300 km/s number assumes that the stars' spin axes point in random directions and that all red main-sequence stars belong to one single group with a simple bell-shaped spread of true speeds, so a hidden second, slower group would change the answer.","fun_headline_variants_meta":{"raw":{"variants":["Fast rotators explain split main sequence in LMC clusters","300 km/s spin links colour to rotation in young clusters","Red main-sequence stars in LMC clusters rotate near break-up","Single-star fast rotation makes most Be stars in LMC clusters","Rotation, not age, drives split main sequence in two LMC clusters"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000198,"raw_usage":{"total_tokens":1457,"prompt_tokens":1123,"completion_tokens":334,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":739,"completion_tokens_details":{"reasoning_tokens":248}},"tokens_in":739,"tokens_out":334,"duration_ms":3827,"temperature":1.0,"reasoning_tokens":248,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T16:36:31.224082+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take high-resolution spectra of the red main-sequence stars in NGC 1866 in the blue, where rotation broadening is cleanly separable and no disk emission contaminates the lines, and reconstruct the distribution of true equatorial velocities without assuming a single Gaussian; if a noticeable slow-rotator population appears, the reported mean of about 300 km/s and the 70–80% critical fraction would be overestimated.","supporting_citations":[],"review_version":2}