The Stellar "Snake"-III: Co-evolution of Stars and Molecular Clouds Unveiled by Gaia, MWISP, and LAMOST

Abstract: By combining multi-band data from Gaia DR3, MWISP CO, and LAMOST DR11 LSR/MSR, we investigate the co-evolution of stars and their parent molecular cloud in a snake-like stellar structure, named Snake III. Based on 5-D phase-space selection, we identified 5683 member stars (median age 7.6 Myr) across approximately $300 \times 500 \times 175$ pc$^3$ volume, along with 12 embedded open clusters. Then we use BEEP distances combined with $^{12}$CO velocities to clearly identify the molecular clouds associated with the stellar complex in spatial and kinematics. The molecular cloud density increases with Galactic longitude, with older open clusters forming in cavities near higher-density regions (except ASCC 125), while young field stars currently form preferentially in present-day high-density environments, indicating that cloud density regulates the star-formation sequence. $^{12}$CO excitation temperature, centroid velocity, velocity dispersion and H$α$ emission reveal that early feedback first compresses cloud edges to trigger new stars, then sweeps and disperses the parent clouds. The extremely young cluster (ASCC 125, 4.4 Myr) lies near the densest region yet is surrounded by a shell with bidirectional density-velocity perturbations, consistent with a delayed-triggering scenario under the combined influence of UBC 178 stellar-wind feedback and a suspected supernova blast. Our results naturally demonstrate that snake-like stellar structures are filamentary relics of hierarchical star formation within giant molecular clouds. They provide direct observational evidence that cloud density and early feedback jointly modulate the progression of star formation, offering a clear and young laboratory for studying star-cloud co-evolution.

• The paper identifies the co-evolution of a young stellar population and molecular clouds within the Snake III complex across a 300×500×175 pc³ volume.
• It employs multi-wavelength data and advanced 5-D phase-space clustering to reveal spatial, kinematic, and age correlations among 5683 stars and 12 open clusters.
• The study highlights feedback-driven processes, with quantitative velocity dispersion analysis confirming the impact of stellar feedback on shaping star formation.

Co-Evolution of Stars and Molecular Clouds in the Stellar “Snake”–III Complex

Introduction

The paper "The Stellar 'Snake'-III: Co-evolution of Stars and Molecular Clouds Unveiled by Gaia, MWISP, and LAMOST" (2604.02717) offers a comprehensive, multi-wavelength analysis of a hierarchical, filamentary stellar structure (Snake III) in the Galactic plane, integrating astrometric data from Gaia DR3, CO molecular line maps from the MWISP survey, and spectroscopic measurements from LAMOST and APOGEE. This study targets the co-evolutionary dynamics between a young, spatially extensive population of stars (5683 members, median age 7.6 Myr, identified via 5-D phase-space clustering) and their natal molecular cloud environment. The authors dissect star-cloud interplay across $\sim$300×500×175 pc$^3$ volume, tracing how spatial and kinematic coherence, local gas density, and stellar feedback mechanisms jointly modulate star formation pathways.

Structural and Kinematic Characterization

Snake III manifests as a prominent stellar filament embedded in the Galactic disk ($90^\circ < l < 115^\circ$, $0^\circ < b < 5^\circ$). The FoF algorithm (ROCKSTAR) isolates member stars using five-dimensional phase space, supplemented with stringent parallax and proper-motion quality cuts. The spatial distribution reveals 12 embedded open clusters, each showing distinct coherence both in Cartesian space and tangential velocity planes.

Figure 1: Spatial distribution of Snake III member stars colored by distance, with marked open clusters and median tangential velocities.

Figure 2: Cartesian $(X, Y, Z)$ distribution of all Snake III stars and clusters relative to the Sun.

Figure 3: Tangential velocity parameters ($V_l$, $V_b$) elegantly separate cluster populations from the complex.

The kinematic profile, including radial velocity histograms, confirms a global median $V_r = -20.57 \pm 0.51$ km/s, underscoring physical homogeneity and common origin.

Figure 4: Radial velocity distribution of member stars along Galactic longitude, consolidated by multiple spectroscopic datasets.

Age Estimation and Population Segregation

Robust age determination hinges on pre-main sequence (PMS) identification via the Sagitta neural model, cross-validated against isochrone fitting for open clusters (ASteCA, PARSEC). PMS probability thresholds ($pms > 0.1$) ensure reliable log(Age) estimates. The four richest clusters (Alessi Teutsch 5, ASCC 125, IC 1396, UBC 178) exhibit ages spanning $4.4-10$ Myr, providing key temporal anchors across the spatial density gradient.

Figure 5: Age histograms for Snake I and III, delineated by PMS probability contours and cluster median log(Age).

Molecular Cloud Mapping and Physical Correlation

MWISP $^3$0CO, $^3$1CO, and C$^3$2O observations, coupled with 3D dust extinction maps (Green et al.), define the molecular architecture coincident with Snake III. High-precision distance anchoring (BEEP method) correlates molecular clouds (G105.2+05.0 to G117.0+03.7) spatially and kinematically with Snake III, forming the Cep OB3 complex. Cloud density, excitation temperature, centroid velocity, and velocity dispersion metrics are derived via LTE formalism and CO-$^3$3 factor conversion.

Figure 6: 3D extinction integrated dust maps overlaid with MWISP $^3$4CO intensity contours at key distance intervals.

Figure 7: Galactic $^3$5 and $^3$6 maps of $^3$7CO intensity, open clusters, and identified clouds; cluster distances annotated.

Cloud column density and mass mapping identifies a monotonic density increase towards higher Galactic longitude, with young clusters forming in higher-density regions (exceptions discussed below).

Figure 8: H$^3$8 column density background (upper: single-star ages, lower: local median ages) with key molecular regions and cluster positions highlighted.

Quantitative Star–Cloud Correlation

Statistical analysis of cluster and field star ages against local gas density reveals clear trends:

• Older clusters (Alessi Teutsch 5, UBC 178) reside in gas-depleted zones, younger (ASCC 125, IC 1396) in high-density regions.
• Field star population age and density contours, box plots, and density bins show young stars preferentially form in present-day dense gas, while older stars dominate lower-density regions.

  Figure 9: Contour and box plots showing log(Age) versus log(H$^3$9 column density) for PMS stars, supporting density-regulated star formation.

Feedback, Dynamical Coupling, and Velocity Structure

Cloud properties (integrated intensity, $90^\circ < l < 115^\circ$0, centroid velocity) and LAMOST H$90^\circ < l < 115^\circ$1 emission maps are synthesized across cluster–cloud regions (A–D), revealing distinct feedback signatures (shells, bubbles, BRCs). Cluster-driven feedback compresses local gas, modulates velocity dispersions, and triggers second-generation star formation. K-S tests on velocity dispersion distributions between perturbed and unperturbed regions yield $90^\circ < l < 115^\circ$2 values significantly above CO sound speeds, indicating supersonic turbulence induced by feedback.

Figure 10: Region A molecular cloud properties with IC 1396, CO intensity, excitation temperature, and centroid velocity.

Figure 11: Regions B–D molecular maps with open clusters, feedback directions, perturbed/unperturbed areas, and H$90^\circ < l < 115^\circ$3 emission overlays.

Figure 12: Velocity dispersion ($90^\circ < l < 115^\circ$4) histograms for each region, showing significant K-S test separation between feedback-perturbed and unperturbed areas.

Second-Generation Cluster Formation and the Case of ASCC 125

ASCC 125 is notably young ($90^\circ < l < 115^\circ$5 Myr) and centrally located in the densest gas, but its formation is delayed relative to neighboring clusters. Detailed K-S testing of age distributions confirms statistically significant divergence. The cluster appears second-generation, triggered by combined feedback from UBC 178 and a possible supernova-induced bubble (“E-bubble”), with detailed P–V analysis supporting a dynamical age for the shell of $90^\circ < l < 115^\circ$6 Myr and measurable expansion velocities.

Figure 13: CDFs for Sagitta ages of four clusters, with K-S $90^\circ < l < 115^\circ$7-values confirming ASCC 125’s uniquely young profile.

Figure 14: $90^\circ < l < 115^\circ$8CO integrated intensity map and P–V diagram of the E-bubble structure, evidencing dynamical impact and shell expansion.

Implications, Future Directions, and Observational Biases

This work validates hierarchical filamentary star formation within GMCs, with direct evidence for density-regulated star formation and feedback-driven propagation of stellar generations. The findings highlight the necessity of high-resolution, multi-wavelength observations to resolve star–gas interactions and to recalibrate mass function and binary fraction estimates in the presence of strong extinction and selection biases. Enhanced coverage via multi-band platforms (e.g., JWST) and interferometric imaging will sharpen constraints on SFE, IMF slope variations, and feedback efficiency.

Conclusion

Snake III exemplifies a physically and kinematically coherent young stellar complex spatially embedded within residual molecular gas, demonstrating:

• Strong numerical evidence for physical homogeneity and common origin across both stellar and molecular phases.
• Statistical correlation between stellar age and gas density, modulated by local feedback, with quantifiable velocity dispersion enhancements in feedback zones.
• The formation of anomalously young, high-density clusters (ASCC 125) through second-generation triggering by cluster winds and supernova remnants.

The study provides a template for dissecting star–cloud co-evolution in hierarchical Galactic structures. It underlines the interplay of initial cloud density gradients and stellar feedback as primary regulators of star formation, offering a uniquely young laboratory for quantifying co-evolutionary histories in star-forming regions.