pith. sign in

arxiv: 2604.22573 · v1 · submitted 2026-04-24 · 🌌 astro-ph.GA · astro-ph.IM· astro-ph.SR

A homogeneous three-dimensional view of Molecular Cloud kinematics out to 2.5 kpc. Using Young Stellar Objects and Open Clusters as complementary tracers

Xabier P\'erez-Couto , Santiago Torres , Nuria Miret-Roig , Friedrich Anders , Alexander J. Mustill , Eva Villaver , Minia Manteiga This is my paper

Pith reviewed 2026-05-08 10:49 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.IMastro-ph.SR
keywords molecular cloudskinematicsGaia DR3young stellar objectsopen clustersMilky Waystar formation3D velocities
0
0 comments X

The pith

Young stellar objects and open clusters trace the same bulk motions of 15 molecular cloud complexes out to 2.5 kpc

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

The paper assembles Gaia DR3 astrometry and radial velocities for over 24,000 young stars and clusters to test whether open clusters retain the same large-scale velocities as young stellar objects in their parent molecular clouds. The two populations agree closely, with a typical offset of only 2 km/s, which justifies treating both as reliable tracers of cloud motion. This agreement produces a uniform three-dimensional velocity field for 15 cloud complexes within 2.5 kpc. The clouds move at a median 8.7 km/s relative to galactic rotation, with clear internal expansion in Orion and Ophiuchus and rotation in at least seven complexes. The resulting map also reconstructs the Sun's past passage through Orion and links several clouds to a common origin in the Sco-Cen region.

Core claim

Using Gaia DR3 astrometry together with complementary spectroscopic surveys for radial velocities, we compiled a unified sample of 24,732 stellar tracers. We applied robust clustering in proper motion space to identify co-moving YSOs and derived cloud-averaged motions via Monte Carlo sampling. These were compared with the kinematics of OCs younger than 30 Myr. We derive homogeneous 3D kinematics for 15 MC complexes within 2.5 kpc. YSOs and OCs exhibit strongly consistent kinematics, with a median spatial velocity offset of ≃ 2 km s^{-1}, confirming that both populations trace the bulk motion of their parent clouds. The resulting cloud kinematics show a median peculiar velocity of ≃ 8.7 km s⁻

What carries the argument

Clustering of YSOs in proper motion space combined with Monte Carlo averaging of 3D velocities for each cloud complex and orbital integration in a realistic galactic potential

If this is right

  • The 3D velocities allow back-integration to show the Sun's passage through the Orion complex and the shared origin of Lupus, Ophiuchus and Corona Australis in Sco-Cen.
  • Orion and Ophiuchus exhibit statistically significant internal expansion at the 5-sigma level.
  • At least seven cloud complexes display coherent internal rotation.
  • The method supplies a homogeneous kinematic reference for studying how molecular clouds interact with galactic rotation and spiral structure.

Where Pith is reading between the lines

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

  • The same stellar-tracer approach could be applied to more distant clouds once deeper surveys provide comparable young-star samples.
  • The measured peculiar velocities of the clouds offer a direct test for models that link cloud motion to galactic spiral arms or supernova feedback.
  • Future Gaia releases may permit extension to somewhat older clusters, allowing the method to track how cloud kinematics evolve after star formation ends.

Load-bearing premise

That identified groups of young stars and clusters younger than 30 million years still carry the overall velocity of the gas they formed from without having drifted apart or been mixed with unrelated stars.

What would settle it

Independent radio maps of the gas velocities in the same clouds showing a systematic offset larger than 4 km/s from the stellar velocities would show that the tracers do not share the bulk motion.

Figures

Figures reproduced from arXiv: 2604.22573 by Alexander J. Mustill, Eva Villaver, Friedrich Anders, Minia Manteiga, Nuria Miret-Roig, Santiago Torres, Xabier P\'erez-Couto.

Figure 1
Figure 1. Figure 1: Galactic map of the molecular clouds from view at source ↗
Figure 2
Figure 2. Figure 2: Resulting clusters of comoving YSOs (blue) in the 2D (µ view at source ↗
Figure 2
Figure 2. Figure 2: (continued) trieved the mean 3D motions of 15 molecular cloud complexes, as shown in view at source ↗
Figure 3
Figure 3. Figure 3: Evolution of the spatial compactness of the Orion (left) and Ophiuchus (right) complexes, showing the median pairwise view at source ↗
Figure 4
Figure 4. Figure 4: Detection probability of a 30 Myr, 100 M⊙ OC located be￾hind the four molecular clouds without any associated OC (see view at source ↗
read the original abstract

Understanding the large-scale dynamics of molecular clouds (MCs) is crucial for constraining the processes that govern star formation and the structure and evolution of the Galaxy. While gas tracers have traditionally been used to map MC kinematics, stellar tracers such as young stellar objects (YSOs) and open clusters (OCs) provide a complementary approach that enables direct comparisons between the stellar and gaseous components. We aim to validate OCs as complementary tracers by testing whether they retain the same bulk kinematic imprint as YSOs, and to reconstruct the three-dimensional (3D) motions of the main MC complexes within 2.5 kpc of the Sun using YSOs and young OCs as tracers. Using Gaia DR3 astrometry together with complementary spectroscopic surveys for radial velocities, we compiled a unified sample of 24,732 stellar tracers. We applied robust clustering in proper motion space to identify co-moving YSOs and derived cloud-averaged motions via Monte Carlo sampling. These were compared with the kinematics of OCs younger than 30 Myr. Finally, we performed orbital integrations in a realistic Galactic potential to trace the past evolution of the clouds and quantify their expansion and rotation. We derive homogeneous 3D kinematics for 15 MC complexes within 2.5 kpc. YSOs and OCs exhibit strongly consistent kinematics, with a median spatial velocity offset of $\simeq 2$ km s$^{-1}$, confirming that both populations trace the bulk motion of their parent clouds. The resulting cloud kinematics show a median peculiar velocity of $\simeq 8.7$ km s$^{-1}$ with respect to Galactic rotation. We trace back the Solar System's voyage through the Orion cloud and the common origin of Lupus, Ophiuchus, and Corona Australis in Sco-Cen. Internally, we detect significant expansion in Orion and Ophiuchus ($5\sigma$) and coherent rotation in at least seven complexes.

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 / 2 minor

Summary. The manuscript presents a homogeneous three-dimensional kinematic analysis of 15 molecular cloud complexes within 2.5 kpc of the Sun, using a compiled sample of 24,732 YSOs and young open clusters (<30 Myr) as tracers from Gaia DR3 astrometry supplemented by spectroscopic radial velocities. Co-moving YSO groups are identified via clustering in proper-motion space, cloud-averaged 3D motions are derived through Monte Carlo sampling, and these are compared to OC kinematics; orbital integrations in a Galactic potential are then used to trace past evolution, expansion, and rotation, yielding a median peculiar velocity of ≃8.7 km s^{-1} and specific findings such as 5σ expansion in Orion and Ophiuchus.

Significance. If the results hold, the work offers a valuable homogeneous dataset of nearby molecular cloud kinematics that complements gas-tracer studies and validates young open clusters as reliable complementary tracers to YSOs. The multi-tracer approach, large sample, and orbital tracing provide useful constraints on Galactic dynamics, star formation processes, and the dynamical history of regions such as Orion and Sco-Cen.

major comments (2)
  1. [Methods section on YSO sample compilation and proper-motion clustering] The central claim that YSOs and OCs exhibit strongly consistent kinematics (median spatial velocity offset ≃2 km s^{-1}) rests on YSO groups identified solely by proper-motion clustering plus parallax cuts. No quantitative interloper fraction is reported, nor is there a 6D membership probability re-analysis or radial-velocity consistency check to bound field-star contamination. Such contamination could bias the Monte Carlo-averaged 3D velocities toward the local standard of rest, artifacting the reported offset and the conclusion that both populations trace bulk cloud motion.
  2. [Results section on cloud kinematics comparison] Although the abstract positions the stellar tracers as complementary to traditional gas tracers, the manuscript provides no direct validation of the derived cloud velocities against independent molecular-gas measurements (e.g., CO line velocities) for the 15 complexes. This cross-check is load-bearing for confirming that the stellar kinematics reflect parent-cloud bulk motions without significant dynamical separation.
minor comments (2)
  1. [Methods] The description of the clustering algorithm would benefit from explicit parameters (e.g., DBSCAN epsilon and min_samples) and a brief statement on how the Monte Carlo sampling propagates Gaia uncertainties and radial-velocity errors.
  2. [Figures] Figure captions for velocity-vector or orbital-path plots should include explicit reference-frame definitions and error-bar conventions to improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed comments on our manuscript. We address each major comment point by point below, indicating where revisions will be incorporated to strengthen the work.

read point-by-point responses

  1. Referee: [Methods section on YSO sample compilation and proper-motion clustering] The central claim that YSOs and OCs exhibit strongly consistent kinematics (median spatial velocity offset ≃2 km s^{-1}) rests on YSO groups identified solely by proper-motion clustering plus parallax cuts. No quantitative interloper fraction is reported, nor is there a 6D membership probability re-analysis or radial-velocity consistency check to bound field-star contamination. Such contamination could bias the Monte Carlo-averaged 3D velocities toward the local standard of rest, artifacting the reported offset and the conclusion that both populations trace bulk cloud motion.

    Authors: We agree that a quantitative interloper fraction is not reported in the current version. In the revised manuscript we will add a dedicated paragraph in the Methods section estimating the expected contamination level from the local field-star density in proper-motion space and the performance of the clustering algorithm. The close agreement between YSO-derived velocities and those of young open clusters (median offset ≃2 km s^{-1}) provides supporting evidence that any residual bias is small, as OCs are subject to different selection effects. A full 6D membership re-analysis is not possible with the present data because radial velocities are unavailable for the majority of the YSO sample; we will explicitly note this limitation and state that the Monte Carlo sampling already propagates the available astrometric uncertainties. We do not expect a systematic pull toward the LSR, given that the derived peculiar velocities remain consistent with independent expectations for the local Galactic environment. revision: partial

  2. Referee: [Results section on cloud kinematics comparison] Although the abstract positions the stellar tracers as complementary to traditional gas tracers, the manuscript provides no direct validation of the derived cloud velocities against independent molecular-gas measurements (e.g., CO line velocities) for the 15 complexes. This cross-check is load-bearing for confirming that the stellar kinematics reflect parent-cloud bulk motions without significant dynamical separation.

    Authors: We concur that a direct comparison to molecular-gas velocities is an important validation step. In the revised manuscript we will add a new table (and accompanying text in the Results section) that lists, for each of the 15 complexes, our stellar-derived 3D velocity together with the corresponding literature CO line velocity (or range) when available. This will allow readers to assess the level of agreement and confirm that the stellar tracers trace the bulk cloud motion. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is observational and data-driven

full rationale

The paper compiles Gaia DR3 astrometry plus spectroscopic radial velocities for 24,732 tracers, applies clustering in proper-motion space to identify co-moving YSO groups, computes Monte-Carlo-averaged 3D velocities, and compares them empirically to young OCs. No load-bearing step reduces the reported median 2 km s^{-1} offset, the 8.7 km s^{-1} peculiar velocities, or the expansion/rotation detections to a fitted parameter, self-citation chain, or ansatz that is defined in terms of the target result. Orbital integrations employ a standard Galactic potential independent of the present sample. The analysis therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The work rests on standard domain assumptions about data accuracy and tracer fidelity rather than new free parameters or invented entities.

axioms (2)
  • domain assumption Gaia DR3 astrometry and complementary radial velocities provide sufficiently accurate 3D velocities for YSOs and OCs to trace molecular cloud bulk motions
    Invoked when compiling the unified sample of 24,732 tracers and deriving cloud-averaged motions.
  • domain assumption Open clusters younger than 30 Myr retain the same bulk kinematic imprint as their parent molecular clouds
    Used to validate OCs as complementary tracers and to compare with YSO kinematics.

pith-pipeline@v0.9.0 · 5702 in / 1330 out tokens · 53174 ms · 2026-05-08T10:49:04.947270+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

6 extracted references · 1 canonical work pages · 1 internal anchor

  1. [1]

    The evolution of velocity dispersion in the Sco-Cen OB association

    Abdurro’uf, Accetta, K., Aerts, C., et al. 2022, Astrophys. J. Suppl. Ser., 259, 35 Akhmetov, V . S., Bucciarelli, B., Crosta, M., et al. 2024, MNRAS, 530, 710 Alves, J., Zucker, C., Goodman, A. A., et al. 2020, Nature, 578, 237 Andrae, R., Fouesneau, M., Sordo, R., et al. 2023, Astron. Astrophys., 674, A27 André, P., Di Francesco, J., Ward-Thompson, D., ...

    work page internal anchor Pith review arXiv 2022

  2. [2]

    These values are mainly aligned with those obtained with the YSOs, while the OCs show a slightly larger difference

    authors obtained values ofµ α∗ =−2.71±0.02 mas yr −1 and µδ =−4.24±0.02 mas yr −1. These values are mainly aligned with those obtained with the YSOs, while the OCs show a slightly larger difference. The RV obtained in this work shows large dispersions for both OCs and YSOs, and within one standard deviation of the value of−11±3 km s −1 reported by Lim et ...

    2019

  3. [3]

    Similarly, the meanRVreported by Galli et al

    As a result, our proper motions are shown to be in very good agreement with those calculated in these previous works. Similarly, the meanRVreported by Galli et al. (2020) for the global cloud is 1.4±0.5 km s −1, but ranging from−0.7± 2.2 in Lupus 1 to 2.1±0.7 km s −1 in Lupus 4 (with RV data for Lupus 5 and 6 not available). These values are within one st...

    2020

  4. [4]

    (2023) to figure out that the results of the YSOs and the OC restricted to these subregions align very well with those of Ortiz-León et al

    Therefore, we used the limits (l,b) of each subregion provided by Dharmawar- dena et al. (2023) to figure out that the results of the YSOs and the OC restricted to these subregions align very well with those of Ortiz-León et al. (2018). The same study published mean Galactic (U,V,W) veloc- ities for these subregions. Using the centroid coordinates (l,b) o...

    2023

  5. [5]

    (2021) used theGaiaastrometry of 403 OB stars, disk-bearing YSOs, and X-ray sources to trace the proper motion of NGC

    Lim et al. (2021) used theGaiaastrometry of 403 OB stars, disk-bearing YSOs, and X-ray sources to trace the proper motion of NGC

    2021

  6. [6]

    On the other hand, Muži´c et al

    The values provided in their work are µα∗ =−1.731 mas yr −1 and µδ =0.312 mas yr −1. On the other hand, Muži´c et al. (2022) used 688 X-ray sources and infrared excess candidates to provide a mean µα∗ =−1.75±0.41 mas yr−1, and µδ =0.25±0.44 mas yr −1. As can be seen, in both cases, the proper motions of those authors match those provided here. We note, ho...

    2022