pith. machine review for the scientific record. sign in

arxiv: 2605.02138 · v1 · submitted 2026-05-04 · 🌌 astro-ph.GA

Recognition: 2 theorem links

· Lean Theorem

Fragmentation in the Serpens/Aquila Star-forming Region

Authors on Pith no claims yet

Pith reviewed 2026-05-08 17:51 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords dense coresstarless coresALMA observationsfragmentationturbulent core collapseAquila regionstar formationmagnetohydrodynamical simulations
0
0 comments X

The pith

ALMA observations detect two completely starless dense cores in Aquila, matching the roughly one detection predicted by turbulent collapse simulations.

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

The paper studies fragmentation in the 100 most unstable dense cores in the Aquila region through ALMA 12m continuum observations at 106 GHz. It identifies 66 sources, two of which show no signs of nearby protostars and thus appear fully starless, plus nine additional starless substructures inside protostellar cores. Synthetic observations drawn from turbulent magnetohydrodynamical simulations of collapsing starless cores forecast 1.19 such detections once central density and observational selection are taken into account. The observed count of two aligns with this forecast, supporting the turbulent core collapse picture. The work further shows that cores with more fragmented parent structures on larger scales develop higher multiplicity on the smallest scales.

Core claim

We present a population study of ALMA Cycle 6 observations of the 100 most gravitationally unstable dense cores in Aquila using a simple mass versus size analysis. We identify 66 continuum sources from ALMA 12m observations at 106 GHz and, through comparisons with known protostellar catalogs, find that two of these detected dense cores appear to be completely starless. The turbulent core collapse model is tested by conducting synthetic observations of turbulent magnetohydrodynamical simulations of collapsing starless cores, which predict at least one (1.19) detection given realistic central densities and density profiles. This prediction is consistent with the two detections of ALMA 12m emis

What carries the argument

Synthetic observations of turbulent magnetohydrodynamical simulations of collapsing starless cores, used to predict the number of detectable starless cores on the basis of central density and density profile.

If this is right

  • The turbulent core collapse model successfully reproduces the observed incidence of detectable starless cores in Aquila.
  • Dense cores frequently fragment into mixed populations of starless and protostellar substructures.
  • The degree of fragmentation on larger spatial scales directly increases the multiplicity observed on the smallest scales.
  • Turbulence shapes substructure development as dense cores collapse to form new star systems.

Where Pith is reading between the lines

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

  • If the match holds in other regions, the fraction of starless cores detectable at millimeter wavelengths may be used to calibrate the time spent in the prestellar phase.
  • Higher-resolution follow-up could test whether the simulated density profiles match the observed cores at radii smaller than those probed here.
  • The scaling of multiplicity with parent fragmentation suggests that initial conditions on cloud scales influence the final stellar multiplicity function.

Load-bearing premise

The turbulent magnetohydrodynamical simulations accurately capture the central densities, density profiles, and observational selection effects of the real Aquila cores without post-hoc adjustments.

What would settle it

A comparable survey of unstable dense cores that yields either zero starless ALMA detections or a number well above two would contradict the simulation prediction of 1.19.

Figures

Figures reproduced from arXiv: 2605.02138 by Helen Kirk, Michael Dunham, Samuel Fielder, Stella Offner.

Figure 1
Figure 1. Figure 1: Column density map of the Serpens South portion of the Aquila Region at 18. ′′2 resolution, from V. K¨onyves et al. (2015). The red circles represent a portion of the 100 dense cores observed by ALMA, with the diameters of the circles representative of the 12 m primary beam. We overplot a scale bar, computed with an adopted distance to Aquila of 436 pc. We use the line-free spectral channels of this last s… view at source ↗
Figure 2
Figure 2. Figure 2: ALMA 12 m mosaic field 183135.8-020349Mosaic containing candidate starless core A65 ext. The color bar ranges linearly from −0.4 − 1.0 mJy beam−1 . The white contours correspond to the associated ACA 7 m continuum emission at the corresponding levels of 3σ (dashed lines), 5σ, 7σ, and 9σ (solid lines), where the 1σ rms is 1.0 mJy beam−1 . All detections are labeled with their source number in white, and pro… view at source ↗
Figure 3
Figure 3. Figure 3: ALMA 12 m mosaic field 183109.7-020622Mosaic, with starless candidate sources A62 ext and A63 ext. See view at source ↗
Figure 4
Figure 4. Figure 4: shows a histogram of the peak flux for the candidate starless and protostellar detections in our ob￾servations. In general, the distribution of protostellar sources extend to higher central densities, and, there￾fore, should have higher peak fluxes compared to their starless counterparts. The majority of our candidate starless core detections lie at lower peak flux bins agree￾ing well with both F24 and H. … view at source ↗
Figure 5
Figure 5. Figure 5: Distribution of ALMA 12 m detection total inte￾grated fluxes. Detections with an integrated flux > 30 mJy are plotted in the final bin shown. Due to the lack of Gaus￾sian fits of the extended sources ( ext), only five of the star￾less detections have reported data. core mass in Orion B North’s was 1.4 M⊙ according to F24). The mean dense core mass in the Aquila HGBS catalog is 0.66 M⊙, however, our observe… view at source ↗
Figure 6
Figure 6. Figure 6: Synthetic ALMA 106 GHz observations of the 4.0 M⊙ simulation, given at six time steps. Each panel indicates the simulation time, along with the central number density. The synthesized beam and scale bar are given in the first panel. Dashed contours represent the 3σ level of emission, while solid contours represent the 5σ level and increase by 2σ, where 1σ rms ∼ 0.08 mJy beam−1 . (see view at source ↗
Figure 7
Figure 7. Figure 7: Top panel: cumulative completeness (red) and associated cumulative expected number of detections (blue; see Equation 4) as a function of number density of the observed prestellar core population (N = 48 cores). Bottom panel: the distribution of peak number density (peak flux divided by observed core area) of all prestellar cores in the HGBS dataset (line) and the ALMA observed prestellar cores (solid). A v… view at source ↗
Figure 8
Figure 8. Figure 8: Mass vs. observed radius for the Group D cores. Blue arrows show the pathways between the HGBS cores and their associated 7 m substructures, while the red arrows show the pathways between the 7 m substructures and the 12 m substructures. The Group D cores across the three spatial scales are highlighted in a red border, while the other empty markers show the rest of the detections. In the HGBS dataset, circ… view at source ↗
Figure 9
Figure 9. Figure 9: Mass vs. observed radius for the Group E cores. See view at source ↗
Figure 10
Figure 10. Figure 10: The ratio of the substructure mass(es) to the parent mass, as a function of the total substructure mass. Circular markers show the prestellar core population, while the star-shaped markers show the protostellar population. Left panel: ACA 7 m detection(s) compared to their parent HGBS cores. Right panel: ALMA 12 m detections compared to their parent ACA 7 m cores. For robust 12 m substructures that are ac… view at source ↗
Figure 11
Figure 11. Figure 11: This figure presents examples of adjacent candidate starless core substructures discussed in Section B.1, illustrating groupings of starless and protostellar components that may represent small clusters in the early stages of formation. Top panel: ALMA 12 m mosaic field 182933.8-015055Mosaic, with starless candidate source A14. Bottom panel: ALMA 12 m mosaic field 182958.2-020115Mosaic, containing the sta… view at source ↗
Figure 12
Figure 12. Figure 12: This figure presents examples of adjacent candidate starless core substructures discussed in Appendix B.1, illustrating groupings of starless and protostellar components that may represent small clusters in the early stages of formation. Top panel: ALMA 12 m mosaic field 183143.8-020440Mosaic, with starless candidate source A66 ext. Bottom panel: ALMA 12 m field 182908.3-013046, containing the starless ca… view at source ↗
Figure 13
Figure 13. Figure 13: HGBS dense cores detection by the ALMA 12 m array (right panels), with no comparable detection with the ALMA-ACA 7 m array (left panels). Plotting conventions generally follow those found in view at source ↗
Figure 14
Figure 14. Figure 14: ALMA 12 m emission showing potential streamer-like morphologies. Top panel: ALMA 12 m field 183139.4-021655 containing the starless candidate source A46. Bottom panel: ALMA 12 m field 183004.1-020305 containing the starless candidate sources A26, A29, A60 ext, and A61 ext. See view at source ↗
Figure 15
Figure 15. Figure 15: HGBS protostellar dense cores only detected in ACA 7 m emission. Panel titles indicate the HGBS dense core and associated ALMA dataset. Plotting conventions generally follow those found in view at source ↗
Figure 16
Figure 16. Figure 16: HGBS protostellar dense cores only detected in ACA 7 m emission. Panel titles indicate the HGBS dense core and associated ALMA dataset. Plotting conventions generally follow those found in view at source ↗
Figure 17
Figure 17. Figure 17: CO outflow signatures for ALMA 12 m detections without direct protostellar associations. Left panels: ALMA 12 m continuum emission is shown in grayscale ranging linearly from −0.4 − 1.0 mJy beam−1 . The synthesized beam is given in the bottom left corner and a scalebar in the bottom right corner. The blue and red contours represent the velocity-shifted components of the 12CO(2 − 1) data, with integrated v… view at source ↗
Figure 18
Figure 18. Figure 18: Left panel: zoomed-in view of the ALMA field 183116.7-020709Mosaic containing source A64 ext. The integrated velocity ranges of the blue- and red-shifted components are −30−5.0 km s−1 and 10.0−29.5 km s−1 , respectively. The contours are drawn at 3σ, 5σ, 7σ, where the noise is 1σ ∼ 0.135 mJy beam−1 . See view at source ↗
Figure 19
Figure 19. Figure 19: Synthetic ALMA 106 GHz observations of the constructed BE-sphere models, given at six central number densities. Each panel indicates the central number density. The synthesized beam and scale bar are given in the first panel. Dashed contours represent the 3σ level of emission, while solid contours represent the 5σ level and increase by 2σ, where 1σ rms is ∼ 0.09 mJy beam−1 view at source ↗
Figure 20
Figure 20. Figure 20: All ALMA 12 m continuum detections not previously shown in the text. ALMA 12 m field 183116.7-020709Mosaic containing continuum detection(s). See view at source ↗
Figure 21
Figure 21. Figure 21: A four-panel view of the observed dense core(s), centered on the ACA 7 m field 182908.3-013046. The top-left panel contains the Herschel H2 column density map, the top-right panel contains the SCUBA-2 450 µm emission map, the bottom-left panel contains ACA 7 m continuum emission map, and the bottom-right panel contains the ALMA 12 m continuum emission map. All panels are linearly normalized to their respe… view at source ↗
read the original abstract

We present a population study of Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 6 observations of the 100 most gravitationally unstable dense cores in Aquila using a simple mass versus size analysis. We identify 66 continuum sources from ALMA 12m observations at 106GHz and through comparisons with known protostellar catalogs; two of these detected dense cores appear to be completely starless, without any accompanying/nearby protostar detections. Additionally, we find nine other starless ALMA 12m detections within protostellar cores that have fragmented into a mixture of starless and protostellar substructures. We test the turbulent core collapse model by conducting synthetic observations of turbulent magnetohydrodynamical simulations of collapsing starless cores in order to predict how many starless cores should be detected given their central density and density profile. The simulations predict at least one (1.19) detection, consistent with our two detections of ALMA 12m emission within completely starless cores. We also use a combination of ALMA Compact Array Cycle 4 observations and the Herschel Gould Belt Survey data to analyze how mass is distributed on three distinct spatial scales, in order to understand how turbulence shapes the evolution of substructure development as dense cores collapse to form new star systems. We find an increase in multiplicity at the smallest scales when the parent larger-scale structure also has a higher degree of fragmentation.

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 paper reports ALMA Cycle 6 12m observations at 106 GHz targeting the 100 most gravitationally unstable dense cores in Aquila, identifying 66 continuum sources including two completely starless cores and nine additional starless substructures within protostellar cores. It tests the turbulent core collapse model by performing synthetic observations on turbulent MHD simulations of collapsing starless cores, which predict 1.19 ALMA detections; this is stated to be consistent with the two observed starless cores. The work also combines ALMA Compact Array Cycle 4 data with Herschel Gould Belt Survey maps to examine mass distribution across three spatial scales and reports an increase in multiplicity at the smallest scales when the parent structure shows higher fragmentation.

Significance. If the simulation-observation agreement is robust, the work supplies a quantitative consistency test of the turbulent core collapse model through predicted versus observed detection rates of starless cores. The multiplicity trend across scales provides observational constraints on how turbulence drives substructure evolution during core collapse. The use of synthetic observations from independent MHD runs is a positive feature that makes the comparison more falsifiable than purely empirical claims.

major comments (2)
  1. [Section describing the turbulent MHD simulations and synthetic observations] The central consistency claim (two observed starless detections versus 1.19 predicted) depends on the simulated cores reproducing the central densities, radial density profiles, and observational selection effects (sensitivity, beam, spatial filtering) of the real Aquila sample. No quantitative metrics comparing simulated and observed core properties (e.g., central n_H2 or power-law index of the density profile) are provided to establish this match.
  2. [Methods section on core selection and simulation setup] The mapping from the observed sample of the 100 most unstable cores (selected via mass-size analysis) to the initial conditions of the MHD simulations is not specified. Without this, it is unclear whether the synthetic observations properly incorporate the same selection biases and core property distribution as the ALMA data.
minor comments (2)
  1. [Abstract] The abstract reports the numerical agreement (2 vs. 1.19) without accompanying uncertainties or details on source extraction criteria; these should be stated explicitly in the main text and abstract if space permits.
  2. [Title and introduction] The title references the Serpens/Aquila region while the abstract and analysis focus exclusively on Aquila; clarify the scope in the title or introduction.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive review and positive assessment of the significance of our work. We address each of the major comments point by point below. Where the comments identify areas that would strengthen the manuscript, we have made revisions to incorporate the requested details and clarifications.

read point-by-point responses
  1. Referee: [Section describing the turbulent MHD simulations and synthetic observations] The central consistency claim (two observed starless detections versus 1.19 predicted) depends on the simulated cores reproducing the central densities, radial density profiles, and observational selection effects (sensitivity, beam, spatial filtering) of the real Aquila sample. No quantitative metrics comparing simulated and observed core properties (e.g., central n_H2 or power-law index of the density profile) are provided to establish this match.

    Authors: We agree that explicit quantitative metrics would make the comparison more robust. The MHD simulations were initialized with central densities and radial density profiles chosen to be representative of the gravitationally unstable cores in the Aquila sample (as derived from the Herschel Gould Belt Survey data used in our core selection). In the revised manuscript, we have added a new subsection with direct comparisons: a table reporting the mean central n_H2 from the simulations (1.2 x 10^6 cm^-3) versus the observed range for the 100 cores, and a figure overlaying the average simulated density profile (power-law index ~1.8 in the inner regions) against the observed profiles. We also detail the synthetic observation pipeline, which applies the exact ALMA 12m Cycle 6 uv-coverage, primary beam, and noise levels from our data to ensure the same sensitivity and spatial filtering effects are modeled. revision: yes

  2. Referee: [Methods section on core selection and simulation setup] The mapping from the observed sample of the 100 most unstable cores (selected via mass-size analysis) to the initial conditions of the MHD simulations is not specified. Without this, it is unclear whether the synthetic observations properly incorporate the same selection biases and core property distribution as the ALMA data.

    Authors: The simulations are not one-to-one mappings to each of the 100 cores but use initial conditions drawn from the statistical distribution of the observed unstable cores. Specifically, the initial core masses, sizes, and central densities are sampled from the mass-size relation and density range of the 100 most unstable cores identified in our Herschel-based selection. In the revised Methods section, we have expanded the description to explicitly state this: the turbulent Mach number and magnetic field strengths are set to values consistent with Aquila observations, and the synthetic observations are performed on an ensemble of 50 simulation runs to capture the population variance. This ensures the predicted detection rate of 1.19 incorporates the same selection biases as the real sample. revision: yes

Circularity Check

0 steps flagged

No circularity: simulation prediction is independent of the new observations

full rationale

The paper's central prediction of 1.19 detections is obtained by running synthetic observations on turbulent MHD simulations of collapsing starless cores (chosen to represent the turbulent core collapse model) and counting how many would be detected at ALMA 12m sensitivity given their central densities and profiles. This number is then compared to the two starless-core detections found in the new Cycle 6 data on the 100 most unstable Aquila cores. No equation, selection criterion, or parameter in the provided text shows the simulation initial conditions, density profiles, or detection threshold being fitted or adjusted to reproduce the observed count of two; the 1.19 figure is therefore not forced by construction. The mass-size analysis, protostar catalog cross-matching, and multi-scale fragmentation statistics are likewise derived directly from the ALMA and Herschel data without reduction to prior self-citations or renamed empirical patterns. The consistency statement is a genuine external test rather than a tautology.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, new entities, or ad-hoc axioms are stated. The work relies on standard assumptions of the turbulent core collapse model and public survey data.

axioms (1)
  • domain assumption Turbulent core collapse model with given central density and density profile accurately represents real cores
    Invoked to generate synthetic observations and predict detection numbers.

pith-pipeline@v0.9.0 · 5560 in / 1278 out tokens · 120523 ms · 2026-05-08T17:51:26.294488+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

74 extracted references · 71 canonical work pages · 1 internal anchor

  1. [1]

    2024, in Comets III, ed

    Andr´ e, P., Di Francesco, J., Ward-Thompson, D., et al. 2014, University of Arizona Press, 27, doi: 10.2458/azu uapress 9780816531240-ch002 Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33, doi: 10.1051/0004-6361/201322068 Astropy Collaboration, Price-Whelan, A. M., Sip˝ ocz, B. M., et al. 2018, AJ, 156, 123, doi: 10....

  2. [2]

    2011, A&A, 527, A145, doi: 10.1051/0004-6361/201015733

    Belloche, A., Schuller, F., Parise, B., et al. 2011, A&A, 527, A145, doi: 10.1051/0004-6361/201015733

  3. [3]

    E., Zhao, B., et al

    Caselli, P., Pineda, J. E., Zhao, B., et al. 2019, ApJ, 874, 89, doi: 10.3847/1538-4357/ab0700

  4. [4]

    534, Protostars and Planets VII, ed

    Chevance, M., Krumholz, M. R., McLeod, A. F., et al. 2023, in ASP Conference Series, Vol. 534 (Kyoto, Japan: Astronomical Society of the Pacific), 1, doi: 10.48550/arXiv.2203.09570

  5. [5]

    2009, MNRAS, 396, 1383, doi:10.1111/j.1365-2966.2009.14843.x

    Dapp, W. B., & Basu, S. 2009, MNRAS, 395, 1092, doi: 10.1111/j.1365-2966.2009.14616.x di Francesco, J., Evans, II, N. J., Caselli, P., et al. 2007, University of Arizona Press, 951, 17, doi: 10.48550/arXiv.astro-ph/0602379

  6. [6]

    M., Crapsi, A., Evans, II, N

    Dunham, M. M., Crapsi, A., Evans, II, N. J., et al. 2008, ApJS, 179, 249, doi: 10.1086/591085

  7. [7]

    M., et al., 2015, @doi [ ] 10.1088/0067-0049/220/1/11 , https://ui.adsabs.harvard.edu/abs/2015ApJS..220...11D 220, 11

    Dunham, M. M., Allen, L. E., Evans, II, N. J., et al. 2015, ApJS, 220, 11, doi: 10.1088/0067-0049/220/1/11

  8. [8]

    M., Offner, S

    Dunham, M. M., Offner, S. S. R., Pineda, J. E., et al. 2016, ApJ, 823, 160, doi: 10.3847/0004-637X/823/2/160

  9. [9]

    2020, ApJS, 251, 20, doi: 10.3847/1538-4365/abba26

    Dutta, S., Lee, C.-F., Liu, T., et al. 2020, ApJS, 251, 20, doi: 10.3847/1538-4365/abba26

  10. [10]

    , keywords =

    Feigelson, E. D., Townsley, L. K., Broos, P. S., et al. 2013, ApJS, 209, 26, doi: 10.1088/0067-0049/209/2/26

  11. [11]

    2026, CADC, doi: 10.11570/25.0103

    Fielder, S., & Kirk, H. 2026, CADC, doi: 10.11570/25.0103

  12. [12]

    D., Kirk, H., Dunham, M

    Fielder, S. D., Kirk, H., Dunham, M. M., & Offner, S. S. R. 2024, ApJ, 968, 10, doi: 10.3847/1538-4357/ad3d56

  13. [13]

    Garrett, J., Luis, E., H. -H. Peng, et al. 2023,, Zenodo doi: 10.5281/ZENODO.4106649

  14. [14]

    Generozov, A., Offner, S. S. R., Kratter, K. M., et al. 2025, Nature, doi: 10.1038/s41550-025-02686-5

  15. [15]

    N., Offner, S

    Guszejnov, D., Raju, A. N., Offner, S. S. R., et al. 2023, MNRAS, 518, 4693, doi: 10.1093/mnras/stac3268

  16. [16]

    A., Bourke, T

    Gutermuth, R. A., Bourke, T. L., Allen, L. E., et al. 2008, ApJ, 673, L151, doi: 10.1086/528710

  17. [17]

    2011, A&A, 533, A34, doi: 10.1051/0004-6361/201117039

    Hacar, A., & Tafalla, M. 2011, A&A, 533, A34, doi: 10.1051/0004-6361/201117039

  18. [18]

    2013, A&A, 554, A55, doi: 10.1051/0004-6361/201220090

    Hacar, A., Tafalla, M., Kauffmann, J., & Kov´ acs, A. 2013, A&A, 554, A55, doi: 10.1051/0004-6361/201220090

  19. [19]

    1999, Astronomische Nachrichten, 320, 362, doi: 10.1002/1521-3994(199908)320:4/5⟨362::AID- ASNA362⟩3.0.CO;2-W

    Hamaguchi, K., Terada, H., Bamba, A., & Koyama, K. 1999, Astronomische Nachrichten, 320, 362, doi: 10.1002/1521-3994(199908)320:4/5⟨362::AID- ASNA362⟩3.0.CO;2-W

  20. [20]

    Hartmann, L., Ballesteros-Paredes, J., & Bergin, E. A. 2001, ApJ, 562, 852, doi: 10.1086/323863

  21. [21]

    2024, ApJ, 961, 123, doi: 10.3847/1538-4357/ad09e2

    Hirano, N., Sahu, D., Liu, S.-Y., et al. 2024, ApJ, 961, 123, doi: 10.3847/1538-4357/ad09e2

  22. [22]

    , keywords =

    Ishihara, K., Nakamura, F., Sanhueza, P., & Saito, M. 2025, A&A, 695, L25, doi: 10.1051/0004-6361/202452427

  23. [23]

    2024, AJ, 974, 95, doi: 10.3847/1538-4357/ad630f

    Ishihara, K., Sanhueza, P., Nakamura, F., et al. 2024, AJ, 974, 95, doi: 10.3847/1538-4357/ad630f

  24. [24]

    L., Evans, N

    Kauffmann, J., Bertoldi, F., Bourke, T. L., Evans, II, N. J., & Lee, C. W. 2008, A&A, 487, 993, doi: 10.1051/0004-6361:200809481

  25. [25]

    M., Di Francesco, J., et al

    Kirk, H., Dunham, M. M., Di Francesco, J., et al. 2017, ApJ, 838, 114, doi: 10.3847/1538-4357/aa63f8

  26. [26]

    2018, ApJS, 238, 8, doi: 10.3847/1538-4365/aada7f

    Kirk, H., Hatchell, J., Johnstone, D., et al. 2018, ApJS, 238, 8, doi: 10.3847/1538-4365/aada7f

  27. [27]

    G., & Iurevich, L

    Kolesnik, I. G., & Iurevich, L. V. 1983, Astrofizika, 19, 761

  28. [28]

    C., Caselli, P., et al

    Kong, S., Tan, J. C., Caselli, P., et al. 2017, ApJ, 834, 193, doi: 10.3847/1538-4357/834/2/193 K¨ onyves, V., Andr´ e, Ph., Men’shchikov, A., et al. 2015, A&A, 584, A91, doi: 10.1051/0004-6361/201525861

  29. [29]

    A., de Souza, R

    Kuhn, M. A., de Souza, R. S., Krone-Martins, A., et al. 2021, ApJS, 254, 33, doi: 10.3847/1538-4365/abe465

  30. [30]

    Larson, R. B. 1969, MNRAS, 145, 271, doi: 10.1093/mnras/145.3.271 Le Gouellec, V. J. M., Hull, C. L. H., Maury, A. J., et al. 2019, ApJ, 885, 106, doi: 10.3847/1538-4357/ab43c2

  31. [31]

    T., Offner, S

    Lee, A. T., Offner, S. S. R., Kratter, K. M., Smullen, R. A., & Li, P. S. 2019, ApJ, 887, 232, doi: 10.3847/1538-4357/ab584b

  32. [32]

    2021, The Journal of Open Source Software, 6, 3771, doi: 10.21105/joss.03771 Mac Low, M.-M

    Li, P., Cunningham, A., Gaches, B., et al. 2021, The Journal of Open Source Software, 6, 3771, doi: 10.21105/joss.03771 Mac Low, M.-M. 1999, ApJ, 524, 169, doi: 10.1086/307784

  33. [33]

    , keywords =

    Mallick, K. K., Kumar, M. S. N., Ojha, D. K., et al. 2013, ApJ, 779, 113, doi: 10.1088/0004-637X/779/2/113

  34. [34]

    V., Paladini, R., et al

    Marton, G., T´ oth, L. V., Paladini, R., et al. 2016, MNRAS, 458, 3479, doi: 10.1093/mnras/stw398

  35. [35]

    2011, A&A, 535, A77, doi: 10.1051/0004-6361/201117132

    Bontemps, S. 2011, A&A, 535, A77, doi: 10.1051/0004-6361/201117132

  36. [36]

    and Ostriker, Eve C

    McKee, C. F., & Ostriker, E. C. 2007, ARA&A, 45, 565, doi: 10.1146/annurev.astro.45.051806.110602 Men’shchikov, A., Andr´ e, Ph., Didelon, P., et al. 2012, A&A, 542, A81, doi: 10.1051/0004-6361/201218797 Aquila Clusters41

  37. [37]

    Moon, S., & Ostriker, E. C. 2025, ApJ, 988, 82, doi: 10.3847/1538-4357/ade239

  38. [38]

    M., van Dishoeck, E

    Murillo, N. M., van Dishoeck, E. F., Hacar, A., Harsono, D., & Jørgensen, J. K. 2022, A&A, 658, A53, doi: 10.1051/0004-6361/202141250

  39. [39]

    M., van Dishoeck, E

    Murillo, N. M., van Dishoeck, E. F., Tobin, J. J., & Fedele, D. 2016, A&A, 592, A56, doi: 10.1051/0004-6361/201628247

  40. [40]

    Offner, S. S. R., Capodilupo, J., Schnee, S., & Goodman, A. A. 2012, MNRAS, 420, L53, doi: 10.1111/j.1745-3933.2011.01194.x

  41. [41]

    Offner, S. S. R., & Chaban, J. 2017, ApJ, 847, 104, doi: 10.3847/1538-4357/aa8996

  42. [42]

    Offner, S. S. R., Moe, M., Kratter, K. M., et al. 2023, in ASP Conference Series, Vol. 534 (Kyoto, Japan: Astronomical Society of the Pacific), 275, doi: 10.48550/arXiv.2203.10066

  43. [43]

    Offner, S. S. R., Taylor, J., & Grud´ ıc, M. Y. 2025, ApJ, 982, 138, doi: 10.3847/1538-4357/adb71d

  44. [44]

    Offner, S. S. R., Taylor, J., Markey, C., et al. 2022, MNRAS, 517, 885, doi: 10.1093/mnras/stac2734 Ortiz-Le´ on, G. N., Dzib, S. A., Loinard, L., et al. 2023, A&A, 673, L1, doi: 10.1051/0004-6361/202346369 Ortiz-Le´ on, G. N., Loinard, L., Dzib, S. A., et al. 2018, ApJ, 869, L33, doi: 10.3847/2041-8213/aaf6ad

  45. [45]

    1994, A&A, 291, 943

    Ossenkopf, V., & Henning, Th. 1994, A&A, 291, 943

  46. [46]

    , keywords =

    Palmeirim, P., Andr´ e, Ph., Kirk, J., et al. 2013, A&A, 550, A38, doi: 10.1051/0004-6361/201220500

  47. [47]

    2024, in SPIE Astronomical Telescopes + Instrumentation, Vol

    Petry, D., D´ ıaz Trigo, Mar´ ıa., Kneissl, R., et al. 2024, in SPIE Astronomical Telescopes + Instrumentation, Vol. 13098 (pan: SPIE), 130980P, doi: 10.1117/12.3012662

  48. [48]

    Nature Astronomy , keywords =

    Pineda, J. E., Segura-Cox, D., Caselli, P., et al. 2020, Nature, 4, 1158, doi: 10.1038/s41550-020-1150-z

  49. [49]

    and Arzoumanian, Doris and Andr

    Pineda, J. E., Arzoumanian, D., Andre, P., et al. 2023, in ASP Conference Series, Vol. 534 (Kyoto, Japan: Astronomical Society of the Pacific), 233, doi: 10.48550/arXiv.2205.03935

  50. [50]

    , keywords =

    Plunkett, A. L., Fern´ andez-L´ opez, M., Arce, H. G., et al. 2018, A&A, 615, A9, doi: 10.1051/0004-6361/201732372

  51. [51]

    C., Dunham, M

    Pokhrel, R., Myers, P. C., Dunham, M. M., et al. 2018, ApJ, 853, 5, doi: 10.3847/1538-4357/aaa240

  52. [52]

    A., Betti , S

    Pokhrel, R., Gutermuth, R. A., Betti, S. K., et al. 2020, ApJ, 896, 60, doi: 10.3847/1538-4357/ab92a2

  53. [53]

    T., Gutermuth, R

    Pokhrel, R., Megeath, S. T., Gutermuth, R. A., et al. 2023, ApJS, 266, 32, doi: 10.3847/1538-4365/acbfac

  54. [54]

    S., Kuhn, M

    Povich, M. S., Kuhn, M. A., Getman, K. V., et al. 2013, ApJS, 209, 31, doi: 10.1088/0067-0049/209/2/31

  55. [55]

    P., Whitney, B

    Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJS, 169, 328, doi: 10.1086/512039

  56. [56]

    P., Whitney, B

    Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 167, 256, doi: 10.1086/508424

  57. [57]

    D., Getman, K

    Romine, G., Feigelson, E. D., Getman, K. V., Kuhn, M. A., & Povich, M. S. 2016, AJ, 833, 193, doi: 10.3847/1538-4357/833/2/193

  58. [58]

    2016, MNRAS, 460, 4150, doi: 10.1093/mnras/stw1100

    Rumble, D., Hatchell, J., Pattle, K., et al. 2016, MNRAS, 460, 4150, doi: 10.1093/mnras/stw1100

  59. [59]

    2021, ApJ, 907, L15, doi: 10.3847/2041-8213/abd3aa

    Sahu, D., Liu, S.-Y., Liu, T., et al. 2021, ApJ, 907, L15, doi: 10.3847/2041-8213/abd3aa

  60. [60]

    2023, ApJ, 945, 156, doi: 10.3847/1538-4357/acbc26

    Sahu, D., Liu, S.-Y., Johnstone, D., et al. 2023, ApJ, 945, 156, doi: 10.3847/1538-4357/acbc26

  61. [61]

    J., Xue, R., & Fourkas, A

    Sheehan, P., Tobin, J. J., Xue, R., & Fourkas, A. 2025,, Zenodo doi: 10.5281/ZENODO.17871742

  62. [62]

    Shu, F. H. 1977, ApJ, 214, 488, doi: 10.1086/155274

  63. [63]

    , keywords =

    Shuping, R. Y., Vacca, W. D., Kassis, M., & Yu, K. C. 2012, ApJ, 144, 116, doi: 10.1088/0004-6256/144/4/116

  64. [64]

    F., Cutri, R

    Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163, doi: 10.1086/498708 Straiˇ zys, V.,ˇCernis, K., & Bartaˇ si¯ ut˙ e, S. 1996, Baltic Astronomy, 5, 125, doi: 10.1515/astro-1996-0106 Straiˇ zys, V., Cernis, K., & Bartasiute, S. 2003, A&A, 405, 585, doi: 10.1051/0004-6361:20030599

  65. [65]

    E., Caselli, P., et al

    Taniguchi, K., Pineda, J. E., Caselli, P., et al. 2024, ApJ, 965, 162, doi: 10.3847/1538-4357/ad2fa1 The CASA Team, Bean, B., Bhatnagar, S., et al. 2022, PASP, 134, 114501, doi: 10.1088/1538-3873/ac9642

  66. [66]

    2014, ApJ, 789, L4, doi: 10.1088/2041-8205/789/1/L4

    Tokuda, K., Onishi, T., Saigo, K., et al. 2014, ApJ, 789, L4, doi: 10.1088/2041-8205/789/1/L4

  67. [67]

    2016, ApJ, 826, 26, doi: 10.3847/0004-637X/826/1/26

    Tokuda, K., Onishi, T., Matsumoto, T., et al. 2016, ApJ, 826, 26, doi: 10.3847/0004-637X/826/1/26

  68. [68]

    Menten, K. M. 2015, A&A, 575, A27, doi: 10.1051/0004-6361/201322955

  69. [69]

    2022 , note =

    Valdivia-Mena, M. T., Pineda, J. E., Segura-Cox, D. M., et al. 2022, A&A, 667, A12, doi: 10.1051/0004-6361/202243310

  70. [70]

    2007, , 119, 855, 10.1086/521277

    Ward-Thompson, D., Di Francesco, J., Hatchell, J., et al. 2007, PASP, 119, 855, doi: 10.1086/521277

  71. [71]

    2000, , 143, 9, 10.1051/aas:2000332

    Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9, doi: 10.1051/aas:2000332

  72. [72]

    L., Eisenhardt, P

    Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868, doi: 10.1088/0004-6256/140/6/1868

  73. [73]

    Y., Laureijs, R

    Zhang, C. Y., Laureijs, R. J., Clark, F. O., & Wesselius, P. R. 1988, A&A, 199, 170

  74. [74]

    2015, ApJ, 804, 141, doi: 10.1088/0004-637X/804/2/141

    Zhang, Q., Wang, K., Lu, X., & Jim´ enez-Serra, I. 2015, ApJ, 804, 141, doi: 10.1088/0004-637X/804/2/141