arxiv: 2605.12222 · v1 · submitted 2026-05-12 · 🌌 astro-ph.GA

Recognition: no theorem link

A tool of Hierarchical cOre ideNtification and Kinematic property AssIgnment (HONKAI) for Dense Cores

Jiawei Liu , Zhiyuan Ren , Di Li , Jinjin Xie , Gary A. Fuller , Yuchen Xing , Xin Lyu , Fengwei Xu

show 3 more authors

Chen Wang Fanyi Meng Sihan Jiao

Authors on Pith no claims yet

Pith reviewed 2026-05-13 04:08 UTC · model grok-4.3

classification 🌌 astro-ph.GA

keywords infrared dark cloudsdense coreshigh-mass star formationvirial ratiocore mass functionhierarchical structure

0 comments

The pith

Dense cores in infrared dark clouds are mostly self-gravitating yet rarely meet the conditions to form high-mass stars without later mass growth.

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

The paper presents HONKAI, an automated procedure that identifies hierarchical cores and clumps in IRDCs, disentangles velocity components from spectral data, measures physical properties, and produces catalogs. Applied to three IRDCs using 850 micron dust continuum from JCMT and 13CO data from the PMO 14-m telescope, the method locates 193 dense cores within 16 clumps. Most cores (136 of 193) exhibit virial ratios above 1, showing they are self-gravitating, but their mass-size relations fall below the threshold associated with high-mass star formation potential, while the core mass function displays a steeper slope at the high-mass end than seen in more evolved samples. These findings together indicate that only a small fraction of the cores currently appear capable of forming high-mass stars.

Core claim

Although many IRDC cores are self-gravitating, only a small fraction are seemingly possible to form high-mass stars. In subsequent core evolution, some further mass assembly trend may be involved to facilitate the high-mass star formation.

What carries the argument

HONKAI, the hierarchical core identification and kinematic property assignment procedure that resolves elemental components, disentangles velocities, measures properties, and generates catalogs from multi-band data.

If this is right

Only a small fraction of IRDC cores currently satisfy conditions for high-mass star formation.
Core mass functions in the earliest stages are steeper at high masses than in evolved regions.
Additional mass assembly is required after the initial core stage to enable high-mass star formation.
Many self-gravitating cores in IRDCs will need external or internal processes to increase mass before forming massive stars.

Where Pith is reading between the lines

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

Star formation models for IRDCs must include continued accretion or merging after the dense core identification stage.
Applying similar automated tools to larger IRDC samples could determine whether the small fraction of high-mass candidates is a general feature.
Comparisons with hydrodynamic simulations of core evolution could check whether the observed mass-size deficit persists without added mass input.

Load-bearing premise

The chosen virial ratio threshold above 1 and mass-size relation cutoff accurately flag limited high-mass star formation potential, and core identification plus velocity disentanglement introduce no major biases in the statistics.

What would settle it

Direct observation of mass growth onto these cores in later evolutionary stages, or a match between the current core properties and the final stellar mass distribution without added accretion, would test the claim.

Figures

Figures reproduced from arXiv: 2605.12222 by Chen Wang, Di Li, Fanyi Meng, Fengwei Xu, Gary A. Fuller, Jiawei Liu, Jinjin Xie, Sihan Jiao, Xin Lyu, Yuchen Xing, Zhiyuan Ren.

**Figure 1.** Figure 1: The maps display the intensity map of three representative star-forming regions from ALOHA observations: sdc19p2 (left), sdc35p0 (center), and sdc38p7 (right). Bright cyan contours outline specific peak-intensity levels. These panels clearly reveal both dense cores and clumps. between tint = 3 and 12 hours. We used the standard procedure in Starlink package10 to calibrate data and produce images. The noise… view at source ↗

**Figure 2.** Figure 2: Visualization of the dust structure dendrogram hierarchy and its spatial footprint. (a) The ACORNS dendrogram built from the dust map, where branches are colored by hierarchical level. (b) The dust intensity map as the background, with multi-scale structure outlines overlaid in the same level-based color corresponding to the dendrogram. The colorbars indicate dust intensity and dendrogram level, respective… view at source ↗

**Figure 3.** Figure 3: The maps display the peak intensity map of the same regions (as [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

**Figure 4.** Figure 4: Identification results of ALOHA 2D dust structures and the MWISP 3D 13CO structures in sdc35p0. The red and cyan outlines display the boundaries of the identified small scale and large scale 2D dust structures. Left: The background image displays the ALOHA dust-continuum emission. The translucent colored patches represent the footprints of different 13CO velocity components identified in PPV space. Right: … view at source ↗

**Figure 5.** Figure 5: The cross-matching between 2D dust map and 13CO in PPV space. From left to right: sdc38p7, sdc19p2, sdc35p0. The background is the 2D dust emission map, dense cores are outlined in red. Different 13CO structures are marked with colors In PPV space, and the component matched to dense cores are highlighted with red. (1) The column density is derived from the dust continuum emission under the assumption of op… view at source ↗

**Figure 6.** Figure 6: The core mass-size relationship. The colors denote the virial mass ratio Rvir, with green for Rvir < 1, red for 1 < Rvir < 3, and black for Rvir > 3. The solid curve represents the empirical threshold for massive star formation (J. Kauffmann & T. Pillai 2010). the threshold for massive star formation. These results suggest that although most cores in our IRDC sample are gravitationally bound and likely uns… view at source ↗

**Figure 7.** Figure 7: Left: Cumulative Core Mass Function with the two power laws fitted to different parts of the data. The power-law index for each fitted line is labeled in the plot. Right: Differential Core Mass Function. The dashed line shows the −2.35 slope of a Salpeter stellar IMF, and the solid line is the best fit to the binned differential CMF with a slope of −1.1. confinement at a very early evolutionary stage. Howe… view at source ↗

read the original abstract

Infrared dark clouds (IRDCs) contains cold dense gas at the earliest stage of massive star and cluster formation. In studying the IRDCs, a universal and fundamental task is to resolve their internal hierarchical structures. Various packages and algorithms were developed for this purpose, but with most of them mainly focused on certain individual steps in data processing. In this work, we build a more automatic procedure for multi-band structure measurement HONKAI (Hierarchical cOre ideNtification and Kinematic property AssIgnment), which can resolve the elemental components including cores and clumps, disentangle the velocity components in spectral data, measure their physical properties, and generate a catalogue for all the measured properties. We use {\sc honkai} for a joint study towards three IRDCs observed in 850 $\mu$m dust continuum with James Clerk Maxwell Telescope (JCMT) and the $^{13}CO$ $(1-0)$ data cube with the Purple Mount Observatory 14-m telescope. 193 dense cores in 16 clumps are identified. As major dynamical properties, a large amount of the cores (136 out of 193) are measured to have large virial ratio of $R_{\rm vir}>1$, but their mass-size relation is bellow the threshold for massive star formation. Meanwhile, core mass function (CMF) also exhibits a steeper slope towards high-mass end compared to more evolved core samples. These three properties in accordance suggest that although many IRDC cores are self-gravitating, only a small fraction are seemingly possible to form high-mass stars. In subsequent core evolution, some further mass assembly trend may be involved to facilitate the high-mass star formation.

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.

Desk Editor's Note private letter to a colleague

HONKAI is a practical automated pipeline for core finding and property measurement in IRDCs that delivers a usable catalog, but the virial ratio claims rest on an unclear definition that creates tension with the self-gravitating interpretation.

read the letter

HONKAI is a new automated pipeline that combines core identification, velocity component separation, and property measurement for IRDC studies, and the application to three clouds gives a catalog of 193 cores with some interesting dynamical trends. The work does a good job integrating existing methods into a single procedure that reduces manual steps. Applying it to JCMT continuum and PMO spectral data produces specific results: 193 dense cores, most with R_vir >1, mass-size relations below high-mass thresholds, and a steeper CMF at the high end. This supports their point that while cores may be self-gravitating, few seem set up for high-mass stars right away, possibly needing more mass assembly later. The main soft spot is the virial ratio part. The abstract reports 136 cores with R_vir >1 but still calls many self-gravitating without spelling out the definition or formula for R_vir. If R_vir follows the usual convention where values over 1 indicate unbound gas, that creates a direct tension with the self-gravitating claim. Even if they mean something else, the lack of explicit equations or error estimates on the ratios makes the dynamical conclusions harder to assess. Completeness tests for the core finding would also strengthen the statistics. This is a paper for people in the dense core and IRDC community who might adopt the tool or use the catalog for comparisons. It does not change the big picture of star formation but offers a practical advance in data processing. I think it deserves peer review. The tool is new and applied usefully, and referees can help clarify the virial definitions and add validation details.

Referee Report

2 major / 2 minor

Summary. The paper introduces HONKAI, an automated pipeline for hierarchical identification of cores and clumps, velocity disentanglement from spectral cubes, and measurement of physical properties in infrared dark clouds. Applied to JCMT 850 μm continuum and PMO 13CO(1-0) observations of three IRDCs, it reports 193 dense cores in 16 clumps. Key results are that 136 cores have R_vir >1, the mass-size relation lies below the high-mass star formation threshold, and the high-mass end of the core mass function is steeper than in more evolved samples. These are interpreted as evidence that many IRDC cores are self-gravitating yet only a small fraction can form high-mass stars, with further mass assembly possibly required in later evolution.

Significance. If the processing steps and dynamical interpretations are robust, the integrated tool could reduce subjectivity in core extraction and kinematic assignment for IRDC studies. The reported statistics on virial ratios, mass-size relations, and CMF slopes would supply useful initial-condition constraints for models of massive star formation.

major comments (2)

[Abstract] Abstract: The central claim that 'although many IRDC cores are self-gravitating' is not supported by the reported statistic. The manuscript states that 136 of 193 cores have R_vir >1, yet the standard virial parameter α >1 indicates unbound cores (M < M_vir). No definition of R_vir, no equation for its calculation (e.g., involving velocity dispersion, radius, and mass), and no explicit self-gravitation threshold (R_vir <1 or <2) are provided. This internal tension must be resolved by stating the exact definition and classification criterion used.
[Abstract] Abstract: The statements that the mass-size relation lies 'below the threshold for massive star formation' and that the CMF 'exhibits a steeper slope towards high-mass end' lack quantitative detail. The specific threshold value and its literature reference, the measured slope value with uncertainty, and the comparison sample(s) are not given. These omissions prevent assessment of whether the three properties are mutually consistent with the headline conclusion.

minor comments (2)

[Abstract] Typo: 'bellow' should read 'below'.
[Abstract] Grammar: 'Infrared dark clouds (IRDCs) contains' should be 'contain'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review. The two major comments both concern the abstract, and we will revise it to add the requested definitions, equations, quantitative values, references, and to resolve the noted internal tension while preserving the scientific content of the manuscript.

read point-by-point responses

Referee: [Abstract] Abstract: The central claim that 'although many IRDC cores are self-gravitating' is not supported by the reported statistic. The manuscript states that 136 of 193 cores have R_vir >1, yet the standard virial parameter α >1 indicates unbound cores (M < M_vir). No definition of R_vir, no equation for its calculation (e.g., involving velocity dispersion, radius, and mass), and no explicit self-gravitation threshold (R_vir <1 or <2) are provided. This internal tension must be resolved by stating the exact definition and classification criterion used.

Authors: We agree that the abstract as currently worded contains an internal inconsistency and omits necessary technical details. In the revised manuscript we will (1) state the exact definition of R_vir, (2) give the equation used to compute it from velocity dispersion, radius and mass, and (3) specify the numerical threshold adopted for self-gravitation. We will then rephrase the concluding sentence so that it is unambiguously supported by the reported statistic of 136 cores with R_vir >1. If our convention for R_vir differs from the standard α (for example if R_vir is defined as the inverse), this will be made explicit; otherwise the wording will be adjusted to reflect that most cores are unbound according to the classical criterion. These changes will be implemented in the abstract and cross-referenced to the methods section. revision: yes
Referee: [Abstract] Abstract: The statements that the mass-size relation lies 'below the threshold for massive star formation' and that the CMF 'exhibits a steeper slope towards high-mass end' lack quantitative detail. The specific threshold value and its literature reference, the measured slope value with uncertainty, and the comparison sample(s) are not given. These omissions prevent assessment of whether the three properties are mutually consistent with the headline conclusion.

Authors: We accept that the abstract would be clearer with explicit numbers and citations. In the revised version we will insert (1) the precise mass-size threshold value together with its literature reference, (2) the measured high-mass-end CMF slope and its uncertainty, and (3) the specific comparison samples against which the slope is judged steeper. These additions will allow readers to evaluate the mutual consistency of the three reported properties with the headline interpretation. revision: yes

Circularity Check

0 steps flagged

No circularity; all reported statistics derive directly from data processing without reduction to self-defined or fitted quantities.

full rationale

The paper introduces the HONKAI pipeline to extract cores, clumps, velocities, masses, sizes, and virial ratios from JCMT 850 μm continuum maps and PMO 13CO(1-0) cubes. The central statistics—193 cores identified, 136 with R_vir > 1, mass-size relations below the high-mass threshold, and a steeper high-mass CMF—are computed outputs of this pipeline applied to the input observations. No equation in the derivation chain defines a target result in terms of a fitted parameter that is then re-labeled as a prediction, nor does any load-bearing step rely on a self-citation whose justification loops back to the present work. The interpretive claim that many cores are self-gravitating yet few can form high-mass stars follows from comparing these independently measured quantities against external literature thresholds; it does not collapse to a tautology or to quantities defined by the authors' own ansatz.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The work relies on standard radio-astronomy definitions of cores and clumps plus common virial-parameter calculations; no new physical entities or ad-hoc constants are introduced beyond routine choices in structure-finding algorithms.

axioms (1)

domain assumption Standard definitions and extraction methods for dense cores from dust continuum and spectral-line data cubes
Invoked throughout the tool description and application section; these are community conventions rather than paper-specific inventions.

pith-pipeline@v0.9.0 · 5650 in / 1253 out tokens · 88107 ms · 2026-05-13T04:08:48.639658+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

25 extracted references · 25 canonical work pages

[1]

A new look at the statistical model identification.IEEE Transactions on Automatic Control19, 716–723 (1974)

Akaike, H. 1974, IEEE Transactions on Automatic Control, 19, 716, doi: 10.1109/TAC.1974.1100705 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.3847/1538-3881/aabc4f Astropy Collaboration, ...

work page doi:10.1109/tac.1974.1100705 1974
[2]

M., et al

Beuther, H., Schilke, P., Menten, K. M., et al. 2002, ApJ, 566, 945, doi: 10.1086/338334

work page doi:10.1086/338334 2002
[3]

Bonnor, W. B. 1956, MNRAS, 116, 351, doi: 10.1093/mnras/116.3.351

work page doi:10.1093/mnras/116.3.351 1956
[4]

2009, ApJS, 181, 360, doi: 10.1088/0067-0049/181/2/360

Simon, R. 2009, ApJS, 181, 360, doi: 10.1088/0067-0049/181/2/360

work page doi:10.1088/0067-0049/181/2/360 2009
[5]

2013, A&A, 552, A40, doi: 10.1051/0004-6361/201219567

Chira, R.-A., Beuther, H., Linz, H., et al. 2013, A&A, 552, A40, doi: 10.1051/0004-6361/201219567

work page doi:10.1051/0004-6361/201219567 2013
[6]

J., Moore, T

Eden, D. J., Moore, T. J. T., Plume, R., et al. 2017, MNRAS, 469, 2163, doi: 10.1093/mnras/stx874

work page doi:10.1093/mnras/stx874 2017
[7]

J., Liu, T., Kim, K.-T., et al

Eden, D. J., Liu, T., Kim, K.-T., et al. 2019, MNRAS, 485, 2895, doi: 10.1093/mnras/stz574

work page doi:10.1093/mnras/stz574 2019
[8]

S., White, S

Frenk, C. S., White, S. D. M., & Davis, M. 1983, ApJ, 271, 417, doi: 10.1086/161209

work page doi:10.1086/161209 1983
[9]

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

work page doi:10.1051/0004-6361/201220090 2013
[10]

D., Longmore, S

Henshaw, J. D., Longmore, S. N., Kruijssen, J. M. D., et al. 2016, MNRAS, 457, 2675, doi: 10.1093/mnras/stw121

work page doi:10.1093/mnras/stw121 2016
[11]

D., Ginsburg, A., Haworth, T

Henshaw, J. D., Ginsburg, A., Haworth, T. J., et al. 2019, MNRAS, 485, 2457, doi: 10.1093/mnras/stz471

work page doi:10.1093/mnras/stz471 2019
[12]

2010, ApJL, 723, L7, doi: 10.1088/2041-8205/723/1/L7

Kauffmann, J., & Pillai, T. 2010, ApJL, 723, L7, doi: 10.1088/2041-8205/723/1/L7

work page doi:10.1088/2041-8205/723/1/l7 2010
[13]

1990, Stellar Structure and Evolution

Kippenhahn, R., & Weigert, A. 1990, Stellar Structure and Evolution

work page 1990
[14]

F., & Langer, W

Li, D., Velusamy, T., Goldsmith, P. F., & Langer, W. D. 2007, ApJ, 655, 351, doi: 10.1086/509736

work page doi:10.1086/509736 2007
[15]

2018, ARA&A, 56, 41, doi: 10.1146/annurev-astro-091916-055235

Motte, F., Bontemps, S., & Louvet, F. 2018, ARA&A, 56, 41, doi: 10.1146/annurev-astro-091916-055235

work page doi:10.1146/annurev-astro-091916-055235 2018
[16]

Peretto, N., & Fuller, G. A. 2009, A&A, 505, 405, doi: 10.1051/0004-6361/200912127

work page doi:10.1051/0004-6361/200912127 2009
[17]

M., & Kr¨ ugel, E

Pillai, T., Wyrowski, F., Menten, K. M., & Kr¨ ugel, E. 2006, A&A, 447, 929, doi: 10.1051/0004-6361:20042145

work page doi:10.1051/0004-6361:20042145 2006
[18]

M., Jackson, J

Rathborne, J. M., Jackson, J. M., & Simon, R. 2006, ApJ, 641, 389, doi: 10.1086/500423

work page doi:10.1086/500423 2006
[19]

A., & Wilson, C

Reid, M. A., & Wilson, C. D. 2005, ApJ, 625, 891, doi: 10.1086/429790

work page doi:10.1086/429790 2005
[20]

J., Menten, K

Reid, M. J., Menten, K. M., Brunthaler, A., et al. 2019, ApJ, 885, 131, doi: 10.3847/1538-4357/ab4a11

work page doi:10.3847/1538-4357/ab4a11 2019
[21]

Goodman, A. A. 2008, ApJ, 679, 1338, doi: 10.1086/587685

work page doi:10.1086/587685 2008
[22]

2019, ApJ, 886, 102, doi: 10.3847/1538-4357/ab45e9

Sanhueza, P., Contreras, Y., Wu, B., et al. 2019, ApJ, 886, 102, doi: 10.3847/1538-4357/ab45e9

work page doi:10.3847/1538-4357/ab45e9 2019
[23]

D., et al

Schinnerer, E., Emsellem, E., Henshaw, J. D., et al. 2023, ApJL, 944, L15, doi: 10.3847/2041-8213/acac9e

work page doi:10.3847/2041-8213/acac9e 2023
[24]

2019, ApJS, 240, 9, doi: 10.3847/1538-4365/aaf1c8

Su, Y., Yang, J., Zhang, S., et al. 2019, ApJS, 240, 9, doi: 10.3847/1538-4365/aaf1c8

work page doi:10.3847/1538-4365/aaf1c8 2019
[25]

2014, MNRAS, 439, 3275, doi: 10.1093/mnras/stu127

Wang, K., Zhang, Q., Testi, L., et al. 2014, MNRAS, 439, 3275, doi: 10.1093/mnras/stu127

work page doi:10.1093/mnras/stu127 2014