NGC 1647: Young Open Cluster Properties

• NGC 1647 is a moderately young open cluster (≈200 Myr) in Taurus–Auriga, noted for its broad main sequence and extended MSTO primarily driven by differential reddening.
• Comprehensive Gaia astrometry and LAMOST spectroscopy have established secure membership, accurate kinematics, and near-solar metallicity for the cluster.
• Joint asteroseismic modeling of δ Sct pulsators yields a seismic age of 178+11/–9 Myr, complementing gyrochronology and lithium depletion analyses.

NGC 1647 is a moderately young open cluster in Taurus–Auriga, projected just behind the Taurus dark clouds and centered at $\mathrm{RA} = 71.48^\circ$, $\mathrm{Dec} = 19.08^\circ$ (J2000.0). Gaia-based analyses place it at about $587$ pc, consistent with a mean parallax of about $1.68$–$1.70$ mas. Its Gaia color–magnitude diagram (CMD) exhibits a conspicuously broad main sequence and an extended main-sequence turn-off (eMSTO), features that have made its age and reddening historically uncertain. Recent work has substantially clarified the cluster through LAMOST spectroscopy, Gaia astrometry and photometry, TESS and K2 time series, and joint asteroseismic modeling of $\delta$ Sct p-mode pulsators (Frasca et al., 13 Apr 2026, Qin et al., 20 Mar 2026).

1. Global properties and membership

NGC 1647 has long had discrepant age estimates in the literature. Gaia-CMD-based studies cited in recent work span values from about $117$ to $363$ Myr, while other compilations gave ages around $260$ Myr or $360$–$\mathrm{Dec} = 19.08^\circ$0 Myr; both recent studies identify reddening, membership, and model assumptions as the main sources of this spread (Frasca et al., 13 Apr 2026, Qin et al., 20 Mar 2026).

A concise summary of the currently reported cluster-scale parameters is given below.

Quantity	Reported value	Basis
Distance	$\mathrm{Dec} = 19.08^\circ$1 pc	Cantat-Gaudin et al. value used in (Qin et al., 20 Mar 2026)
Mean parallax	$\mathrm{Dec} = 19.08^\circ$2–$\mathrm{Dec} = 19.08^\circ$3 mas	Gaia astrometry in (Qin et al., 20 Mar 2026)
Mean metallicity	$\mathrm{Dec} = 19.08^\circ$4 dex	LAMOST cool-star weighted mean in (Frasca et al., 13 Apr 2026)
Spectroscopic metallicity	$\mathrm{Dec} = 19.08^\circ$5 dex	19-star weighted mean in (Qin et al., 20 Mar 2026)
Seismic metallicity	$\mathrm{Dec} = 19.08^\circ$6 dex	Joint $\mathrm{Dec} = 19.08^\circ$7 Sct modeling in (Qin et al., 20 Mar 2026)
Cluster radial velocity	$\mathrm{Dec} = 19.08^\circ$8, $\mathrm{Dec} = 19.08^\circ$9	Slow-rotator subsample in (Frasca et al., 13 Apr 2026)
Mean extinction	$587$0 mag	SED fitting in (Frasca et al., 13 Apr 2026)
Extinction range	$587$1 to $587$2 mag	SED fitting in (Frasca et al., 13 Apr 2026)
Lithium age	$587$3 Myr	EAGLES fit in (Frasca et al., 13 Apr 2026)
Seismic age	$587$4 Myr	Joint p-mode fit in (Qin et al., 20 Mar 2026)

Membership has been constructed in different but complementary ways. One study defined a “golden sample” of $587$5 high-confidence candidates common to the Gaia-based catalogs of Cantat-Gaudin et al. (2018), Hunt & Reffert (2024), and Qin et al. (2026), with low expected contamination among LAMOST targets (Frasca et al., 13 Apr 2026). Another began from all Gaia DR3 sources within a projected radius of $587$6 pc around the adopted center, applied astrometric and photometric quality cuts to obtain Sample 1 of $587$7 stars, used HDBSCAN in $587$8 space to isolate Sample 2 of $587$9 stars, and then adopted $1.68$0 to obtain $1.68$1 candidates; radial-velocity clipping reduced this to $1.68$2 high-confidence members (Qin et al., 20 Mar 2026).

These membership constructions are not numerically identical because they serve different purposes and start from different parent samples. Taken together, they indicate that NGC 1647 is now characterized by a relatively secure central membership set and a larger surrounding candidate population.

2. Spectroscopic characterization and cluster kinematics

The most extensive spectroscopic characterization to date is based on LAMOST medium-resolution spectra with nominal $1.68$3, measured here as $1.68$4 and $1.68$5, with wavelength coverage $1.68$6–$1.68$7 Å in the blue arm and $1.68$8–$1.68$9 Å in the red arm. Cross-matching Gaia members with LAMOST DR12 yielded $1.70$0 targets with $1.70$1 MRS spectra; after $1.70$2 cuts, $1.70$3 stars with $1.70$4 co-added spectra were analyzed, and three bright upper-main-sequence stars were added through archival VLT/UVES spectra (Frasca et al., 13 Apr 2026).

Atmospheric parameters were derived with ROTFIT. For stars cooler than about $1.70$5 K, the template grid consisted of ELODIE spectra of $1.70$6 slowly rotating, low-activity FGKM stars degraded to LAMOST resolution. For hotter stars, BT-Settl synthetic spectra with solar metallicity and $1.70$7–$1.70$8 K, $1.70$9–$\delta$0 dex were used, with $\delta$1 fixed at $\delta$2 dex. The final products for individual stars were $\delta$3, $\delta$4, $\delta$5, $\delta$6, and $\delta$7, with a floor at $\delta$8 set by resolution and sampling (Frasca et al., 13 Apr 2026).

The analyzed members span $\delta$9 K to $117$0 K, and $117$1 ranges from below the detection limit to $117$2 for the hottest B-type stars. This places NGC 1647 in the regime where precise spectroscopy is possible for cool dwarfs, while the upper main sequence is dominated by broad-lined rapid rotators (Frasca et al., 13 Apr 2026).

For the high-quality slow-rotator subsample with $117$3 ($117$4 stars), the cluster radial-velocity distribution is described by

$117$5

Including fast rotators gives $117$6 and $117$7. The intrinsic dispersion is stated to be narrower than the median formal uncertainty from ROTFIT ($117$8); this suggests that the internal velocity dispersion is probably below a few $117$9, consistent with a dynamically cold open cluster (Frasca et al., 13 Apr 2026).

Metallicity determinations are driven by the cool-star sample. A Gaussian fit to the $363$0 distribution gives $363$1 dex, while the weighted mean is

$363$2

A separate compilation based on $363$3 LAMOST-MRS members, one APOGEE star, and one Nordic Optical Telescope/FIES spectrum yields $363$4 dex, and the two results are explicitly described as consistent within uncertainties (Frasca et al., 13 Apr 2026, Qin et al., 20 Mar 2026).

The spectroscopic survey also revealed $363$5 stars with $363$6 in multi-epoch radial-velocity variability, flagged as “RVvar”, and identified four double-lined spectroscopic binaries through double peaks in the cross-correlation function. For these SB2s, component radial velocities were measured, while atmospheric parameters were excluded from the main parameter table because they are not reliable for composite spectra (Frasca et al., 13 Apr 2026).

3. Rotation, chromospheric activity, and lithium chronology

Rotation periods were derived from TESS photometry. One study reports that $363$7 candidate members have TESS light curves, with $363$8 s cadence in sectors 43–44 and, for some targets, $363$9 s cadence in sectors 70–71. Lomb–Scargle periodograms were computed with attention to frequencies $260$0, and visual inspection of light curves, periodograms, and phase-folded curves led to $260$1 stars being classified as having reliable rotational periods (Frasca et al., 13 Apr 2026).

In the color–period diagram, NGC 1647 lies between the $260$2 Myr Pleiades and the $260$3 Myr cluster NGC 3532. The slow-rotator I-sequence overlaps the Pleiades upper sequence at bluer colors and approaches the NGC 3532 sequence at redder colors, while the fast-rotator C-sequence resembles the Pleiades more than NGC 3532. No explicit analytic gyrochronology relation was fitted, but the morphology was reported as consistent with an intermediate age compatible with $260$4 Myr (Frasca et al., 13 Apr 2026).

For stars with $260$5 K, chromospheric activity was measured through H$260$6 after subtraction of non-active photospheric templates matched in $260$7, $260$8, and $260$9, shifted to the stellar radial velocity and broadened to the measured $360$0. Net H$360$1 equivalent widths were then converted into $360$2 and

$360$3

All analyzed NGC 1647 stars lie below the accretion boundary derived in Frasca et al. (2015) and cluster around the same locus as Pleiades members in the $360$4–$360$5 plane, indicating purely chromospheric activity rather than accretion (Frasca et al., 13 Apr 2026).

Lithium provides the sharpest non-seismic age diagnostic in the spectroscopic study. For stars with $360$6 K, Li I $360$7 Å equivalent widths were measured from template-subtracted red-arm spectra, with uncertainties estimated from

$360$8

where $360$9, $\mathrm{Dec} = 19.08^\circ$00 is the integration window, and $\mathrm{Dec} = 19.08^\circ$01. Lithium abundances $\mathrm{Dec} = 19.08^\circ$02 were derived from the non-LTE curves of growth of Lind et al. (2009), and the age was fitted in the $\mathrm{Dec} = 19.08^\circ$03–$\mathrm{Dec} = 19.08^\circ$04 plane with EAGLES using $\mathrm{Dec} = 19.08^\circ$05 stars, including both detections and upper limits. The resulting lithium age is

$\mathrm{Dec} = 19.08^\circ$06

The same study states that this is compatible with both gyrochronology and isochrone fitting (Frasca et al., 13 Apr 2026).

The activity, lithium, and rotation diagnostics therefore converge on a cluster that is no longer pre-main-sequence, but has not yet reached complete rotational convergence across the full mass range.

4. Differential reddening and the broadened main sequence

The defining morphological feature of NGC 1647 is its broad main sequence and eMSTO in the Gaia CMD. In the spectroscopic analysis, the lower envelope of the main sequence was defined empirically by fitting a fourth-order polynomial to the bluest edge of the observed sequence over the relevant magnitude range. For each star, the color offset relative to that lower envelope was then defined as

$\mathrm{Dec} = 19.08^\circ$07

The extinction vector for $\mathrm{Dec} = 19.08^\circ$08 mag was shown to be roughly aligned with the direction of the observed broadening (Frasca et al., 13 Apr 2026).

SED fitting of $\mathrm{Dec} = 19.08^\circ$09 likely members yielded individual extinctions from $\mathrm{Dec} = 19.08^\circ$10 to $\mathrm{Dec} = 19.08^\circ$11 mag and a weighted mean

$\mathrm{Dec} = 19.08^\circ$12

Independent photometric extinctions from Zdanavičius et al. (2005) gave $\mathrm{Dec} = 19.08^\circ$13 mag, and the two sets are well correlated for the $\mathrm{Dec} = 19.08^\circ$14 stars in common, with Pearson coefficient $\mathrm{Dec} = 19.08^\circ$15. Assuming $\mathrm{Dec} = 19.08^\circ$16, the corresponding mean color excess is

$\mathrm{Dec} = 19.08^\circ$17

Spatially, $\mathrm{Dec} = 19.08^\circ$18 correlates with IRAS $\mathrm{Dec} = 19.08^\circ$19 dust emission across the field (Frasca et al., 13 Apr 2026).

The key CMD result is the correlation between extinction and color offset. Using SED-based $\mathrm{Dec} = 19.08^\circ$20, the Pearson coefficient between $\mathrm{Dec} = 19.08^\circ$21 and $\mathrm{Dec} = 19.08^\circ$22 is $\mathrm{Dec} = 19.08^\circ$23; using the Zdanavičius et al. (2005) extinctions, it is $\mathrm{Dec} = 19.08^\circ$24. By contrast, the correlation between $\mathrm{Dec} = 19.08^\circ$25 and color offset is essentially null, with $\mathrm{Dec} = 19.08^\circ$26. The reported interpretation is therefore unambiguous: differential reddening is the primary driver of the broad main sequence and the apparent eMSTO in NGC 1647, whereas rotation is not the dominant cause in this cluster (Frasca et al., 13 Apr 2026).

This conclusion is reinforced by the asteroseismic study, which initially adopted a single global reddening $\mathrm{Dec} = 19.08^\circ$27 from Guerrero et al. (2011) for CMD fitting, but later noted that differential extinction of $\mathrm{Dec} = 19.08^\circ$28–$\mathrm{Dec} = 19.08^\circ$29 mag can shift $\mathrm{Dec} = 19.08^\circ$30 by about $\mathrm{Dec} = 19.08^\circ$31 mag and $\mathrm{Dec} = 19.08^\circ$32 by about $\mathrm{Dec} = 19.08^\circ$33 mag. In that framework, the CMD broadening was explicitly treated as a major source of age degeneracy in photometric isochrone fitting (Qin et al., 20 Mar 2026).

A common misconception in the interpretation of eMSTO morphologies is that such structures necessarily imply intrinsic age spreads or rotationally modified turn-offs. In NGC 1647, the available evidence does not require either explanation: a coeval population with spatially variable line-of-sight extinction is sufficient to reproduce the observed CMD spread.

5. Variable-star population and asteroseismic age determination

The time-domain survey of NGC 1647 identified a rich pulsating and rotationally modulated population. Using $\mathrm{Dec} = 19.08^\circ$34 K2 light curves and $\mathrm{Dec} = 19.08^\circ$35 TESS light curves, processed with local $\mathrm{Dec} = 19.08^\circ$36 clipping, detrending with wotan, concatenation, normalization, and Lomb–Scargle periodograms, one study found $\mathrm{Dec} = 19.08^\circ$37 periodic variables among $\mathrm{Dec} = 19.08^\circ$38 members. The census comprises $\mathrm{Dec} = 19.08^\circ$39 rotational variables, $\mathrm{Dec} = 19.08^\circ$40 slowly pulsating B-type stars, $\mathrm{Dec} = 19.08^\circ$41 $\mathrm{Dec} = 19.08^\circ$42 Dor stars, $\mathrm{Dec} = 19.08^\circ$43 $\mathrm{Dec} = 19.08^\circ$44 Sct stars, $\mathrm{Dec} = 19.08^\circ$45 $\mathrm{Dec} = 19.08^\circ$46 Sct–$\mathrm{Dec} = 19.08^\circ$47 Dor hybrids, $\mathrm{Dec} = 19.08^\circ$48 eclipsing binaries, and $\mathrm{Dec} = 19.08^\circ$49 variables of uncertain type (Qin et al., 20 Mar 2026).

Among the nine p-mode pulsators, five stars showed sufficiently regular high-frequency patterns to permit measurement of a large separation $\mathrm{Dec} = 19.08^\circ$50: TIC 18438742, TIC 18442420, TIC 18480935, TIC 18547844, and TIC 385502143. Frequency extraction used TESS short-cadence data and the MultiModes package for iterative pre-whitening, while mode pattern recognition used échelle diagrams generated with the echelle package. For these stars, measured $\mathrm{Dec} = 19.08^\circ$51 values are typically $\mathrm{Dec} = 19.08^\circ$52–$\mathrm{Dec} = 19.08^\circ$53, and identified ridges include multiple $\mathrm{Dec} = 19.08^\circ$54 and $\mathrm{Dec} = 19.08^\circ$55 modes, with a few $\mathrm{Dec} = 19.08^\circ$56 or $\mathrm{Dec} = 19.08^\circ$57 assignments in two objects (Qin et al., 20 Mar 2026).

The asteroseismic analysis treated age and $\mathrm{Dec} = 19.08^\circ$58 as common cluster parameters and stellar masses as star-by-star nuisance parameters. Evolutionary tracks were computed with MESA and oscillation frequencies with GYRE over

$\mathrm{Dec} = 19.08^\circ$59

and ages from the zero-age main sequence to $\mathrm{Dec} = 19.08^\circ$60 Myr, sampled every $\mathrm{Dec} = 19.08^\circ$61 Myr. The models did not include rotation. The per-star objective function combined the large-separation term and the frequency-matching term, and the total cluster statistic was the average of the five per-star minima over mass (Qin et al., 20 Mar 2026).

The final inference used a likelihood framework based on $\mathrm{Dec} = 19.08^\circ$62, a kernel density estimate in the $\mathrm{Dec} = 19.08^\circ$63 plane, and a likelihood-ratio $\mathrm{Dec} = 19.08^\circ$64 contour defined by

$\mathrm{Dec} = 19.08^\circ$65

Within the resulting high-probability region, the adopted seismic parameters are

$\mathrm{Dec} = 19.08^\circ$66

The study reports mean $\mathrm{Dec} = 19.08^\circ$67 and a mean difference between observed and modeled $\mathrm{Dec} = 19.08^\circ$68 of about $\mathrm{Dec} = 19.08^\circ$69, together with good ridge alignment in the échelle diagrams (Qin et al., 20 Mar 2026).

The asteroseismic age lies inside the isochrone-based range $\mathrm{Dec} = 19.08^\circ$70–$\mathrm{Dec} = 19.08^\circ$71, corresponding to $\mathrm{Dec} = 19.08^\circ$72–$\mathrm{Dec} = 19.08^\circ$73 Myr, but is much more precise than classical CMD fitting. The same study explicitly attributes the residual photometric age degeneracy to the eMSTO, binaries, rotation, and especially differential reddening (Qin et al., 20 Mar 2026).

6. Astrophysical interpretation and significance

The combined spectroscopic and asteroseismic picture places NGC 1647 between the Pleiades and NGC 3532 in evolutionary state. Its rotation distribution retains both fast and slow branches, its H$\mathrm{Dec} = 19.08^\circ$74 activity is chromospheric and Pleiades-like over the overlapping $\mathrm{Dec} = 19.08^\circ$75 range, and its lithium depletion pattern is consistent with a near-solar-metallicity cluster near $\mathrm{Dec} = 19.08^\circ$76 Myr (Frasca et al., 13 Apr 2026).

The most important structural conclusion is that NGC 1647 is a counter-example to the assumption that eMSTO-like CMD morphologies are primarily rotational in origin. In this cluster, the observed CMD spread correlates with $\mathrm{Dec} = 19.08^\circ$77, spatially follows the dust distribution seen in IRAS $\mathrm{Dec} = 19.08^\circ$78 emission, and does not correlate with $\mathrm{Dec} = 19.08^\circ$79. This establishes differential reddening as the dominant effect shaping the broad main sequence and apparent turn-off extension (Frasca et al., 13 Apr 2026).

At the same time, NGC 1647 is an important demonstration case for ensemble $\mathrm{Dec} = 19.08^\circ$80 Sct asteroseismology. Joint modeling of five cluster p-mode pulsators yields a seismic age of $\mathrm{Dec} = 19.08^\circ$81 Myr and a seismic metallicity consistent with spectroscopy, showing that cluster asteroseismology can bypass some of the reddening and CMD-degeneracy limitations that affect purely photometric determinations (Qin et al., 20 Mar 2026).

Methodological limitations remain explicit. The asteroseismic models are non-rotating, despite measured $\mathrm{Dec} = 19.08^\circ$82 values such as $\mathrm{Dec} = 19.08^\circ$83 for TIC 18547844, and core overshooting was not separately tuned. The quoted uncertainties therefore capture statistical errors at fixed physics more directly than full model-systematic errors. Even with those caveats, the convergence of spectroscopic metallicity, lithium depletion, rotational morphology, isochrone placement, and p-mode seismology marks NGC 1647 as a well-calibrated open cluster with mildly subsolar composition, strong and spatially variable extinction, and an age securely in the young-to-intermediate regime (Qin et al., 20 Mar 2026, Frasca et al., 13 Apr 2026).