pith. sign in

arxiv: 2604.11765 · v1 · submitted 2026-04-13 · 🌌 astro-ph.SR · astro-ph.GA

NGC 1647: A young open cluster with a broad main sequence observed with LAMOST

A. Frasca , M. Qin , J. ALonso-Santiago , G. Catanzaro , J. N. Fu , J. Y. Zhang , A. Bragaglia This is my paper

Pith reviewed 2026-05-10 15:11 UTC · model grok-4.3

classification 🌌 astro-ph.SR astro-ph.GA
keywords open clustersmain-sequence turn-offdifferential reddeninglithium abundancestellar rotationspectroscopic parametersNGC 1647
0
0 comments X

The pith

Differential reddening across NGC 1647 produces its broad main sequence, not fast rotation.

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

The paper shows that the apparent spread in the main-sequence turn-off of this young cluster arises from varying amounts of interstellar dust along different lines of sight. Spectra and photometry together demonstrate that the color offset of stars from the lower main-sequence boundary tracks the local reddening derived from spectral-energy-distribution fits, while showing no relation to measured rotation speeds. This distinction matters because many young clusters display similar extended turn-offs, and attributing them to reddening rather than intrinsic stellar differences simplifies models of cluster formation and age determination.

Core claim

Medium-resolution LAMOST spectra of 158 members yield effective temperatures, gravities, metallicities, radial velocities, and projected rotation velocities. TESS light curves provide rotation periods for 160 stars. Spectral-energy-distribution fitting reveals an average visual extinction of 1.1 magnitudes that varies across the field. The local extinction values correlate tightly with the color offset from the lower main-sequence boundary, whereas the same offset shows no correlation with vsini. Lithium equivalent widths give an age of 203 plus or minus 27 million years, consistent with both gyrochronological and isochrone ages.

What carries the argument

Differential reddening measured from spectral-energy-distribution fits, shown to correlate with color offset while rotation velocity does not.

If this is right

  • The extended main-sequence turn-off in NGC 1647 is explained without invoking multiple stellar populations or rotation-driven effects.
  • Lithium depletion, gyrochronology, and isochrone fitting all converge on an age near 200 million years for this cluster.
  • Similar broad main sequences in other young open clusters should first be tested for spatially varying reddening before other explanations are adopted.
  • The cluster radial-velocity dispersion of 1.6 km/s and mean metallicity of -0.08 dex align with expectations for a typical solar-neighborhood open cluster of this age.

Where Pith is reading between the lines

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

  • Surveys of other young clusters could use the same reddening-versus-color-offset test to decide whether their extended turn-offs require exotic interpretations.
  • If differential reddening proves common, it would reduce the need to invoke rapid early mass loss or binary interactions to explain apparent age spreads.
  • Gyrochronology calibrated on clusters like this one can be applied with greater to field stars once reddening effects are isolated.

Load-bearing premise

The observed color offset from the lower main-sequence boundary is produced entirely by differential reddening and is not significantly affected by unresolved binaries, stellar activity, or inaccuracies in the inactive photospheric templates.

What would settle it

High-resolution imaging or spectroscopy that finds a substantial population of unresolved binaries concentrated at the turn-off, or a re-analysis after reddening correction that still shows a strong vsini-color correlation.

Figures

Figures reproduced from arXiv: 2604.11765 by A. Bragaglia, A. Frasca, G. Catanzaro, J. ALonso-Santiago, J. N. Fu, J. Y. Zhang, M. Qin.

Figure 1
Figure 1. Figure 1: Left panel: Spatial distribution of the NGC 1647 members according to Hunt & Reffert (2024, magenta dots), Cantat-Gaudin et al. (2018, blue dots), and Qin et al. (2026, green dots). Symbol size scales with the G magnitude. The red circles highlight the data points that correspond to objects observed with LAMOST MRS. The black squares enclose the three stars with UVES archive observations. The meaning of th… view at source ↗
Figure 2
Figure 2. Figure 2 [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of the atmospheric parameters derived from the blue- and red-arm LAMOST MRS spectra using ROTFIT. Panels, from left to right, show the effective temperature (Teff ), surface gravity (log g), metallicity ([Fe/H]), and projected rotation velocity (v sin i). The symbols are consistent with those presented in [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Radial velocity distribution obtained from all ana￾lyzed LAMOST MRS spectra of the cluster members (green his￾togram). Separate distributions are shown for the slowly rotating and fast rotating stars, as indicated in the legend. The best￾fit Gaussian function applied specifically to the “slowly-rotating stars” histogram is overplotted as a solid cyan line; its center (µ) and dispersion (σ) are also marked.… view at source ↗
Figure 6
Figure 6. Figure 6: Extinction distribution for all the NGC 1647 members with LAMOST data (green histogram) and for the stars ana￾lyzed by Zdanavičius et al. (2005) (Zda05, orange histogram). The average AV values from both datasets are marked with ver￾tical solid lines and reported with the corresponding color in the upper left corner of the box. The AV values determined by us via the SED fitting of the 160 sources with know… view at source ↗
Figure 8
Figure 8. Figure 8: HR diagram for NGC 1647. The MS stars are denoted with dots (red for LAMOST MRS and cyan for UVES spec￾tra). The orange squares represent the two giants studied by Carrera et al. (2022). PARSEC isochrones (Nguyen et al. 2022) at 100, 150, 200, and 300 Myr for models with the cluster metal￾licity ([Fe/H] = −0.08 dex) are overlaid. Green squares enclose the lower probability members indicated by arrows in [… view at source ↗
Figure 9
Figure 9. Figure 9: Hα luminosity ratio (R ′ Hα) versus effective tempera￾ture (Teff ). The plot distinguishes between the ‘golden sample’ of likely members and the lower-probability targets using dif￾ferent colors, as specified in the legend. Upper limits on R ′ Hα are indicated by downward gray arrows. The straight dashed line represents the empirical boundary derived by Frasca et al. (2015) separating purely chromospheric … view at source ↗
Figure 11
Figure 11. Figure 11: Rotation periods versus the dereddened color index (GBP − GRP)0 for the members of NGC 1647 (black crosses). The periods of Pleiades (τ ≃ 125 Myr) members based on K2 (Rebull et al. 2016) and TESS (Frasca et al. 2025) are overplot￾ted with orange dots, while those of NGC 3532 (τ ≃ 300 Myr) members (Fritzewski et al. 2021) are also shown by blue dots for comparison. After obtaining the de-reddened colors, … view at source ↗
Figure 12
Figure 12. Figure 12: Projected rotational velocity (v sin i) versus the Gaia color shift ∆(GBP − GRP). The color shift is measured relative to the lower envelope of the MS strip (Fig B.6). The Pearson’s correlation coefficients (ρ) is marked in the upper right corner of the panel. cently, eMSTOs and broadened main sequences have been identified within Galactic OCs (e.g., Cordoni et al. 2018, 2024, and references therein). How… view at source ↗
read the original abstract

In this work we present the results of our analysis of medium-resolution LAMOST spectra of candidate members of the cluster NGC 1647 with the aim of determining the stellar parameters, activity level, lithium abundance, and to study the cluster properties. We used the code ROTFIT to determine the atmospheric parameters (Teff, logg, and [Fe/H]), radial velocity (Vr), and projected rotation velocity (vsini) for 158 cluster members. Moreover, for solar-type and cooler stars (Teff< 6500 K), we calculated the H-alpha and LiI-6708 net equivalent width by means of the subtraction of inactive photospheric templates. We determined the rotation periods for 160 stars by analyzing the available TESS photometry. We found four double-lined spectroscopic systems for which we provide the radial velocities of the two components. The Vr distribution of the cluster members peaks at -5.3 km/s with a dispersion of 1.6 km/s, while the average metallicity is [Fe/H]=-0.08$\pm$0.08 dex, in line with previous determinations. From the fitting of the spectral energy distribution of 160 likely members we infer the existence of a differential reddening across the cluster field with an average value of $A_V$=1.1 mag. The $A_V$ values show a distinct correlation with the color offset from the lower boundary of the main sequence; conversely, this offset appears to be uncorrelated with vsini. These two findings confirm that differential reddening is the primary driver behind the observed extended Main-Sequence Turn-Off (eMSTO) in this cluster. The age of NGC 1647, obtained from the lithium abundance, is 203$\pm$27 Myr, which is compatible with the values inferred from a gyrochronological approach and the isochrone fitting.

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.

Referee Report

2 major / 3 minor

Summary. The manuscript analyzes medium-resolution LAMOST spectra of 158 candidate members of the young open cluster NGC 1647. Using the ROTFIT code, the authors derive Teff, log g, [Fe/H], Vr, and vsini; measure net equivalent widths of H-alpha and Li I 6708 via template subtraction for stars cooler than 6500 K; identify four SB2 systems; extract rotation periods from TESS photometry for 160 stars; and perform SED fitting to obtain individual AV values. They report a correlation between AV and the photometric color offset from the lower main-sequence boundary, but no correlation between the offset and vsini. These results are used to conclude that differential reddening (average AV = 1.1 mag) is the primary driver of the observed extended main-sequence turn-off. The lithium-based age is given as 203 ± 27 Myr, stated to be consistent with gyrochronological and isochrone ages.

Significance. If the central claim holds after addressing alternatives, the work would demonstrate that spatially variable extinction can produce an apparent eMSTO in a ~200 Myr cluster without requiring rotation or age spreads. The multi-method age consistency and the direct AV-color correlation provide a useful case study for interpreting broadened main sequences in young open clusters.

major comments (2)
  1. [Abstract and CMD/reddening analysis] Abstract and the section presenting the AV-color offset correlation: the statement that the AV correlation and lack of vsini correlation 'confirm that differential reddening is the primary driver' of the eMSTO is not supported by the analysis. The manuscript identifies only four SB2 systems among 158 members but provides no binary-fraction estimate, no binary-cleaned subsample, and no population synthesis to show whether unresolved binaries could produce a comparable redward color spread or mimic the observed AV correlation. Binaries are known to shift stars redward in the CMD independently of reddening; without this test the primacy claim cannot be established.
  2. [Lithium analysis and age determination] Section on lithium abundance and age: the age of 203 ± 27 Myr is derived from Li I equivalent widths after template subtraction, yet the manuscript does not specify the exact inactive templates employed, the error budget on the net EWs, or the precise Li-depletion calibration (including reference models) used to convert abundance to age. This prevents quantitative assessment of the quoted uncertainty and of the stated compatibility with the gyrochronological and isochrone ages.
minor comments (3)
  1. [Abstract and membership section] The abstract states that AV values were obtained from SED fitting of 160 likely members, yet the main text reports 158 spectroscopic members; a clear statement of the final membership list and any overlap or selection cuts is needed.
  2. [Reddening and CMD analysis] The correlation between AV and color offset is described qualitatively; reporting the Spearman or Pearson coefficient together with its significance would strengthen the result.
  3. [Photometric analysis] The definition of the 'lower boundary of the main sequence' used to compute the color offset is not stated; a brief description or reference to the fitting procedure would improve reproducibility.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The comments highlight important points that we address below. We will revise the manuscript to strengthen the analysis and provide additional details where needed.

read point-by-point responses

  1. Referee: [Abstract and CMD/reddening analysis] Abstract and the section presenting the AV-color offset correlation: the statement that the AV correlation and lack of vsini correlation 'confirm that differential reddening is the primary driver' of the eMSTO is not supported by the analysis. The manuscript identifies only four SB2 systems among 158 members but provides no binary-fraction estimate, no binary-cleaned subsample, and no population synthesis to show whether unresolved binaries could produce a comparable redward color spread or mimic the observed AV correlation. Binaries are known to shift stars redward in the CMD independently of reddening; without this test the primacy claim cannot be established.

    Authors: We agree that unresolved binaries can contribute to the color spread in the CMD and that a dedicated test would strengthen the conclusion. However, the observed correlation is specifically between the color offset and AV values derived independently from SED fitting; binaries are not expected to induce a systematic correlation with extinction. The absence of correlation with vsini further supports that rotation is not the driver. We identified only four SB2 systems in the sample of 158 members. In the revised manuscript we will add a dedicated discussion of the possible role of binaries, provide a rough binary-fraction estimate based on the spectroscopic detections and available TESS photometry, and clarify that differential reddening is the primary driver while binaries remain a possible secondary contributor. revision: partial

  2. Referee: [Lithium analysis and age determination] Section on lithium abundance and age: the age of 203 ± 27 Myr is derived from Li I equivalent widths after template subtraction, yet the manuscript does not specify the exact inactive templates employed, the error budget on the net EWs, or the precise Li-depletion calibration (including reference models) used to convert abundance to age. This prevents quantitative assessment of the quoted uncertainty and of the stated compatibility with the gyrochronological and isochrone ages.

    Authors: We acknowledge that additional methodological details are required for full reproducibility and assessment of the uncertainty. In the revised manuscript we will specify the inactive templates used for the subtraction (selected from the ROTFIT library of slowly rotating stars with no detectable Li absorption or Hα emission), describe how the error budget on the net EWs was evaluated (including residuals from the subtraction and photon-noise contributions), and provide the exact Li-depletion calibration together with the reference models employed to convert the measured EWs into age. revision: yes

Circularity Check

0 steps flagged

No significant circularity; claims rest on independent observational correlations and measurements

full rationale

The paper derives stellar parameters via ROTFIT on LAMOST spectra, computes net equivalent widths by template subtraction, extracts rotation periods from TESS light curves, and performs SED fitting to obtain A_V values. The reported correlation between A_V and photometric color offset from the lower main-sequence boundary, together with the absence of correlation with vsini, constitutes a direct empirical test rather than a quantity defined by the fit itself. The lithium-based age is obtained from abundance measurements and stated to be compatible with independent gyrochronological and isochrone results. No step reduces by construction to its inputs, no fitted parameter is relabeled as a prediction, and no load-bearing self-citation or uniqueness theorem is invoked. The analysis chain is therefore self-contained against external photometric and spectroscopic benchmarks.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The analysis assumes standard stellar-atmosphere models and template spectra are adequate for the cluster's metallicity and activity levels; no new physical entities are introduced.

free parameters (1)
  • average AV
    Inferred from SED fitting of 160 members and reported as 1.1 mag.
axioms (1)
  • domain assumption Color offset from the lower main-sequence boundary is produced by differential reddening
    Invoked to interpret the correlation between AV and color offset.

pith-pipeline@v0.9.0 · 5677 in / 1206 out tokens · 23439 ms · 2026-05-10T15:11:19.698815+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

2 extracted references · 2 canonical work pages

  1. [1]

    2020, arXiv e-prints, arXiv:2005.07210, doi: 10.48550/arXiv.2005.07210

    Abdurro’uf, Accetta, K., Aerts, C., et al. 2022, ApJS, 259, 35 Allard, F. 2014, in Exploring the Formation and Evolution of Plan- etary Systems, ed. M. Booth, B. C. Matthews, & J. R. Graham, Vol. 299, 271–272 Alonso-Santiago, J., Frasca, A., Catanzaro, G., et al. 2021, A&A, 656, A149 Bai, Y., Liu, J., Wang, Y., & Wang, S. 2020, AJ, 159, 84 Barnes, S. A. 2...

    work page arXiv 2022

  2. [2]

    The differences∆Vr=V LAMOST r − V Gaia r are displayed in the bottom panel along with their aver- age,<∆V r>, and standard deviation, rms

    The one-to-one relation is shown by the solid blue line. The differences∆Vr=V LAMOST r − V Gaia r are displayed in the bottom panel along with their aver- age,<∆V r>, and standard deviation, rms. Fig. B.3.Top:Example of an SED fitting for the hot star J044636.9+190649 (= HD285997).Bottom:χ 2-contour map of the fitting. The red contour corresponds to the 1...

    work page arXiv 2022