pith. sign in

arxiv: 2604.02717 · v1 · submitted 2026-04-03 · 🌌 astro-ph.GA · astro-ph.SR

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

Jia-Peng Li , Hai-Jun Tian , Chen Wang , Xiang-Ming Yang , Fan Wang This is my paper

Pith reviewed 2026-05-13 18:41 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.SR
keywords snake-like stellar structuresgiant molecular cloudshierarchical star formationearly stellar feedbackopen clustersGaia DR3co-evolution
0
0 comments X

The pith

Snake-like stellar chains are relics of step-by-step star formation inside giant molecular clouds.

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

The study maps a winding chain of young stars spanning hundreds of parsecs using positions, motions, and distances from Gaia together with gas maps from CO surveys and spectra from LAMOST. It finds that the stars and their parent gas clouds are physically linked, with star formation progressing from denser to less dense regions as time passes. Early stars push and compress the gas, triggering more stars while eventually clearing out the cloud. This provides a live example of how giant clouds build stars hierarchically over time, leaving behind filamentary star patterns.

Core claim

Combining multi-band data from Gaia DR3, MWISP CO, and LAMOST DR11, the analysis identifies 5683 member stars with a median age of 7.6 Myr and 12 embedded open clusters in the Snake III structure. Molecular cloud density increases along the structure, older clusters sit in cavities near denser areas, and young stars form in current high-density zones. CO temperature, velocity, and dispersion plus H-alpha emission indicate early feedback compresses edges then disperses clouds, with the youngest cluster ASCC 125 near the densest gas showing shell-like perturbations from stellar winds and a possible supernova.

What carries the argument

5-D phase-space selection of stars combined with BEEP distances and CO velocities to associate the stellar complex with its parent molecular clouds.

Load-bearing premise

The 5-D phase-space selection and BEEP distances with CO velocities correctly identify physically associated stars and clouds without major contamination from unrelated objects along the line of sight.

What would settle it

If the spatial and velocity overlaps between the selected stars and CO clouds are no stronger than expected from random alignments in the galactic plane.

Figures

Figures reproduced from arXiv: 2604.02717 by Chen Wang, Fan Wang, Hai-Jun Tian, Jia-Peng Li, Xiang-Ming Yang.

Figure 1
Figure 1. Figure 1: Spatial distribution of the member stars of Snake III. The color bar encodes distance. Black circles mark the centers and angular sizes of the open clusters. The black arrow in the top-right corner indicates the median tangential velocity of the entire sample, while the remaining arrows show the median tangential velocities of the individual cluster samples, with their lengths proportional to the velocity … view at source ↗
Figure 2
Figure 2. Figure 2: Spatial distribution of Snake III member stars (gray dots) and cluster member stars (shown as dots in var￾ious colors) in Cartesian coordinates (X, Y, Z). The coordi￾nates of the Sun are (X, Y, Z) = (0.0, 0.0, 0.0) kpc than 100 member stars (Alessi Teutsch 5, ASCC 125, IC 1396 and UBC 178), their ample membership num￾bers enable relatively reliable age estimates, which are also the key targets of our follo… view at source ↗
Figure 3
Figure 3. Figure 3: Tangential velocity in (l, b) distribution of Snake III member stars (gray dots) and cluster member stars (shown as dots in various colors). All velocities are given with respect to the LSR. sumption of LTE (J. E. Pineda et al. 2008; J. L. Pineda et al. 2010), these lines yield a host of physical parame￾ters. 12CO is generally optically thick (P. F. Goldsmith & W. D. Langer 1999); its excitation temperatur… view at source ↗
Figure 4
Figure 4. Figure 4: Distribution of the observed Vrs for all member candidates along Galactic longitude l. The Vrs are obtained from the multiple surveys (color-coded as specified in the legend). The histograms of Vrs are displayed in the right sub-panel. The median values are marked with dashed lines. The color-coding is the same as the main panel. 6.2 6.4 6.6 6.8 7.0 7.2 7.4 7.6 7.8 0 500 1000 1500 2000 2500 3000 N 6.850±0.… view at source ↗
Figure 6
Figure 6. Figure 6: The background shows the Green 3D dust-extinc￾tion map integrated over 0.7–0.9 kpc (upper) and 0.9–1.1 kpc (lower). Overlaid are the MWISP 12CO integrated-intensity contours for integration range of [-20,5] km s−1 . G106.1+00.5, G106.5+04.0, G106.5+01.6, G106.6+01.0, G107.7+02.9, G108.4+00.3, G109.0- 00.1, G110.9+03.5, G111.7+00.0, G114.5-00.1 and G117.0+03.7. In the integrated intensity map of 12CO shown … view at source ↗
Figure 5
Figure 5. Figure 5: log Age distribution histograms of the Snake I member stars (upper panel) and Snake III member stars (lower panel) fitting by Sagitta. Contour lines of different colors delineate intervals of PMS probability (green for pms < 0.01, red for pms > 0.1, and blue for pms > 0.8), and the median log Age and number of stars calculated for each interval is marked in the corresponding color. backdrop to interpret th… view at source ↗
Figure 7
Figure 7. Figure 7: Galactic longitude-latitude map of 12CO integrated intensity over velocity (upper panel) and longitude–velocity map of 12CO integrated intensity over Galactic latitude (lower panel) overlaid with Snake III member stars (gray dots). Some molecular cloud distance and velocity released by the BEEP method marked as a, b, c, d, e, f, g, h, i, j, k, and l; and the IC 1396 H II region marked as m. The four open c… view at source ↗
Figure 8
Figure 8. Figure 8: Upper panel: H2 column density background map with stellar age l − b distribution plot; the color bar shows the distribution of single-star age estimates fitted using Sagitta (pms > 0.1). Lower panel: H2 column density background map with local stellar age l − b distribution plot; the color bar shows the distribution of median local stellar age estimates fitted using Sagitta (pms > 0.1), and the rectangula… view at source ↗
Figure 9
Figure 9. Figure 9: Contour plot (upper panel) and box plot distri￾bution (lower panel) of Sagitta estimated stellar ages (pms > 0.1) versus local H2 column density values at their positions; both age and H2 column density are taken in logarithmic scale. mostly already merged into the four identified clusters. This reveals the statistical result that as the gas density at the stellar location increases, the stellar age struct… view at source ↗
Figure 10
Figure 10. Figure 10: The 12CO integrated intensity, excitation temperature (Tex) and centroid velocity (Vcen) of clouds in Region A and the scatter plot of IC 1396 in Galactic coordinates (l, b). All molecular cloud-related maps use a velocity integration range of [−15, −5] km s−1 . The dark-red solid box marks the Region A molecular clouds. The black solid circle marks IC 1396 (circle size scaled to the cluster radius). The … view at source ↗
Figure 11
Figure 11. Figure 11: The 12CO integrated intensity, excitation temperature (Tex) and centroid velocity (Vcen) of clouds in Regions B, C, and D, the scatter plot of stars and LAMOST MRA-N Hα emission in Galactic coordinates (l, b) among sub-panels (a), (b), and (c). All molecular cloud-related maps use a velocity integration range of [−15, −5] km s−1 . The dark-red solid boxes mark the Regions B, C, and D molecular clouds. The… view at source ↗
Figure 12
Figure 12. Figure 12: Histogram comparing the velocity dispersion (σv) distributions selected in Regions B, C, D. Area I (green, unper￾turbed) and Area II (red, perturbed) are plotted in different colors, with their median σv values and uncertainties calculated respectively. The DKS and p values shown in the label represent the K-S tests performed for each comparison region, all demonstrating the significance of the difference… view at source ↗
Figure 13
Figure 13. Figure 13: The cumulative distribution function (CDF) plots of Sagitta ages (pms > 0.1) for the four open clusters, with logarithmic age on the x-axis and cumulative probabil￾ity (0–1) on the y-axis. The p-values from K-S tests between other clusters and ASCC 125 are labeled in the figure to in￾dicate its significance (p ≪ 0.001). Alessi Teutsch 5, and IC 1396) formed in progressively less dense portions of the orig… view at source ↗
Figure 14
Figure 14. Figure 14: The 12CO integrated intensity map of the bub￾ble in Region E (left panel) and the position-velocity (P-V) diagram of the selected PV-belt (right panel). In the left panel, the rectangular dashed box marks the extent of the PV-belt, with an arrow inside the box indicating the start￾ing direction; the solid circle marks the physical scale of the E-bubble and its estimated radius (R ≈ 0.75◦ ). In the right p… view at source ↗
read the original 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$\alpha$ 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.

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

1 major / 1 minor

Summary. The manuscript investigates the co-evolution of stars and molecular clouds in the Snake III stellar structure by combining Gaia DR3, MWISP CO, and LAMOST DR11 data. It identifies 5683 member stars (median age 7.6 Myr) and 12 open clusters using 5D phase-space selection in a volume of approximately 300 × 500 × 175 pc³. Molecular clouds are associated using BEEP distances and 12CO velocities, revealing that cloud density regulates the star-formation sequence, with early feedback compressing cloud edges and dispersing parent clouds. The results suggest that snake-like structures are filamentary relics of hierarchical star formation modulated by cloud density and feedback.

Significance. If the physical associations hold, this provides direct observational evidence for hierarchical star formation within giant molecular clouds and the joint modulation by cloud density and early feedback. The large multi-tracer sample offers a clear young laboratory for star-cloud co-evolution studies, with potential to inform models of density-regulated progression and feedback effects.

major comments (1)
  1. [Abstract] Abstract: The central claim that snake-like structures are filamentary relics of hierarchical star formation modulated by cloud density and early feedback rests on the 5-D phase-space selection plus BEEP distances and 12CO velocities accurately identifying physically associated members without significant projection effects. Over the 300×500×175 pc³ volume, no quantitative contamination fraction, completeness estimate, or control-field test is referenced, leaving open the possibility that field stars or unrelated clouds with overlapping velocities contaminate the sample and affect the reported age–density–feedback correlations.
minor comments (1)
  1. [Abstract] Abstract: The reported numbers (5683 stars, median age 7.6 Myr) would benefit from accompanying error bars, membership probability thresholds, or uncertainty ranges to allow readers to assess robustness directly from the summary.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive and detailed review of our manuscript. We address the major comment below and will incorporate the suggested improvements in the revised version to better substantiate the robustness of our member selection and physical associations.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The central claim that snake-like structures are filamentary relics of hierarchical star formation modulated by cloud density and early feedback rests on the 5-D phase-space selection plus BEEP distances and 12CO velocities accurately identifying physically associated members without significant projection effects. Over the 300×500×175 pc³ volume, no quantitative contamination fraction, completeness estimate, or control-field test is referenced, leaving open the possibility that field stars or unrelated clouds with overlapping velocities contaminate the sample and affect the reported age–density–feedback correlations.

    Authors: We agree that explicit quantitative validation of the sample is necessary to support the central claims. The current manuscript describes the 5D phase-space selection (proper motions, parallaxes from Gaia DR3, radial velocities from LAMOST) in Section 2 but does not include the requested contamination or completeness metrics. In the revised manuscript, we will add a new subsection (Section 2.3) that: (1) applies identical selection criteria to a control field at similar Galactic latitude but offset by ~5° in longitude to quantify the contamination fraction; (2) estimates completeness using the Gaia DR3 selection function and synthetic populations matched to our volume; and (3) demonstrates the kinematic coherence by comparing the 12CO velocity distribution of associated clouds against the control field and field stars. These results will be summarized in the abstract and used to qualify the age–density–feedback correlations. The combination of precise BEEP distances and tight velocity matching already reduces projection effects, but the new tests will provide the quantitative support requested. revision: yes

Circularity Check

0 steps flagged

No significant circularity; claims rest on external survey data

full rationale

The paper identifies 5683 member stars and 12 clusters via 5-D phase-space selection from Gaia DR3, then associates them with molecular clouds using BEEP distances and 12CO velocities from MWISP. These steps rely on direct observational cross-matching over the stated volume rather than any fitted parameter, self-defined quantity, or prediction that reduces to the input selection by construction. No equations appear that rename or force the reported density-age-feedback sequence; the co-evolution interpretation follows from the spatial-kinematic associations themselves. The analysis is therefore self-contained against the cited external catalogs with no load-bearing self-citation or ansatz smuggling required for the central claim.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claims rest on standard domain assumptions about kinematic membership and cloud association rather than new free parameters or invented entities.

axioms (2)
  • domain assumption 5-D phase-space selection accurately isolates true member stars of the complex
    Used to identify the 5683 member stars and 12 clusters
  • domain assumption BEEP distances plus CO velocities correctly link molecular clouds to the stellar structure in 3D
    Required to map density variations and feedback signatures

pith-pipeline@v0.9.0 · 5633 in / 1307 out tokens · 25707 ms · 2026-05-13T18:41:44.550388+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

3 extracted references · 3 canonical work pages

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

    work page

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

    work page

  3. [3]

    Yw P ' W;8+! ݣ W _r+L > bl

    thebibliography [1] 20pt to REFERENCES 6pt =0pt \@twocolumntrue 12pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key o...

    work page doi:10.3847/1538-4357/abbf4b 2017