{"id":"49ec8b20-1385-4308-b3d7-0b2dafa5c4aa","arxiv_id":"1908.05622","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A stacked analysis of 315 dense cores around 53 HII regions finds an r^-3 volume density profile of cores, a boundary excess, an interior deficit, and no significant heating of clumps by the HII regions or OB stars.","lead":"Stacking 315 dense cores from 38 SCUBA-2 fields shows that HII regions sit at the centers of extended core clusters, with surface density falling as r^-2 and an extra ring of cores just outside each ionized bubble. The paper also finds that mature HII regions and their OB stars do not heat nearby dust clumps, and that most measured cores are warmer than their surrounding clouds.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The r^-3 volume-density claim is underdetermined: the projected stacked counts are Abel-inverted under spherical symmetry, while the paper's own §4 admits filamentary geometries can fit the same data. Kinematic or isotropy checks are needed before the central claim is accepted.","rationale":"The Reader's weakest-assumption analysis and my stress-test converge: the Abel inversion to n(r)∝r^-3 assumes spherical symmetry, and the paper itself flags the filamentary alternative in Section 4. I examined other candidate weaknesses—map-edge incompleteness at large Θ, circle-fitting to non-spherical HII regions, and the SFE-count inconsistency in the abstract—but none is as load-bearing as the 3D-geometry degeneracy. The empirical histogram and boundary excess can survive, yet the headline r^-3 and 'at the center of a cluster' statement would be wrong if the cores trace a network of filaments. The paper's own suggestion of radial-velocity follow-up is the natural test. Consequently the conditional-accept verdict stands: the work is honest and valuable, but the central 3D interpretation is not uniquely established.","tokens_in":16359,"tokens_out":9528,"duration_ms":99540,"concrete_test":"Obtain C18O or N2H+ radial velocities for the cores in 3–5 representative fields (e.g., Sh-2 104, Sh-2 305, Sh-2 168). If the velocities reveal coherent filamentary gradients, multiple velocity components, or a flat, non-central dispersion profile instead of a smoothly decreasing, roughly isotropic velocity dispersion expected for a relaxed spherical cluster, then the Abel-inverted r^-3 volume-density interpretation fails. If the velocity field is centrally symmetric and centrally concentrated, the spherical interpretation is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that each HII region lies at or near the center of a cluster of dense cores whose volume density falls as r^-3. This is obtained in two steps: the stacked projected histogram gives N(Θ) = (31.6±7.3) Θ^{-1.1±0.2} for Θ≥1; converting ring counts to surface density gives Σ∝R^{-2.1}; and Abel inversion under spherical symmetry gives n(r)∝r^{-3.1}. The weakest link is the last step. A radial histogram is a one-dimensional projection of the actual two-dimensional core distribution, and many 3D geometries produce exactly the same projected profile. The authors concede this in Section 4: filamentary structures with a distribution of lengths and non-uniform core positions can reproduce the observed counts, and with that freedom one could fit almost any number-count distribution. Thus the data do not uniquely determine that the cores form a spherical cluster; the r^-3 volume-density profile is an interpretive assumption rather than a measurement. The same assumption also underlies the paper's statement that the HII region is 'at the center' of the core cluster. No kinematic information is used, and the paper itself recommends radial-velocity observations to settle the geometry. Without such a test, the strongest claim is plausible but not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a stacking analysis of 315 dense cores detected in SCUBA-2 450/850 μm images around 53 galactic HII regions. The authors scale each core's separation by the radio radius of its associated HII region and construct a stacked radial count histogram, which they fit with a power law N = (31.6±7.3) Θ^(−1.1±0.2) cores per bin for Θ≥1. Interpreting the projected counts as a spherical distribution, they argue that the surface density of cores falls as r^(−2.1) and the volume density as r^(−3.1), that there is an excess of cores near the HII region boundary, and that the interior of the HII regions is largely devoid of cores. They also report no significant heating of clouds and cores by the HII regions or their OB stars, a cloud temperature–column density anti-correlation, and star formation efficiencies of 1–9% for the most reliable systems.","tokens_in":16538,"tokens_out":7813,"duration_ms":67609,"significance":"If the interpretation holds, the result that mature HII regions sit at the centers of extended, r^(−3) core clusters would be a novel and important constraint on feedback and star formation triggering, and the boundary excess would support collect-and-collapse as a general mechanism rather than a rare phenomenon. The paper's main strength is its large, homogeneous sample of 53 HII regions and 315 cores, which allows a stacking approach that case-by-case analyses cannot provide. The power-law count distribution is shown to be robust to binning choices, and the authors are transparent about many of the assumptions underlying their conversions, including the spherical-symmetry assumption and the alternative filamentary interpretation.","major_comments":[{"comment":"The conversion of the projected count histogram to a volume density profile via Abel inversion assumes that the three-dimensional distribution of cores is spherically symmetric about each HII region. The paper's own Section 4 acknowledges that \"by using more complex filamentary structures, such as a distribution of the lengths of the filaments and/or a non-uniform distribution of the dense cores along the filaments, it is possible to fit these number counts\" and that \"with this level of freedom to choose parameters one could fit almost any number count distribution.\" Therefore the projected data alone do not uniquely determine a spherical geometry; the r^(−3) volume-density profile (abstract and Section 3.2) is an interpretive assumption rather than a measured result. The authors should either provide an isotropy or kinematic test (e.g., radial velocity observations, as they themselves suggest) or explicitly restate the central claim as conditional on spherical symmetry, with the filamentary alternative given equal weight in the abstract and conclusions.","section":"§3.2 and §4"},{"comment":"The stacking procedure divides each core's separation by the radio radius of its associated HII region, which assumes that all HII regions have self-similar core distributions. The sample, however, spans physical radii from 0.35 pc to 20.9 pc (Table 1), and the paper does not test whether the scaled profile is robust to this normalization. Without such a test—for example, splitting the sample at the median radius and comparing the stacked profiles—the fitted power-law index could partly reflect the mixing of different physical scales. The authors should demonstrate that the stacked result is insensitive to the radius distribution of the sample.","section":"§3.2"},{"comment":"The background/foreground contamination is estimated by fitting a toy model N = N0 π Θ^2 to the outer tail (20 ≤ Θ ≤ 25) of the same dataset used in the main analysis, under the assumption that these bins contain only background/foreground cores. This is a self-referential estimate, and although the resulting contamination of 2–6 cores is small, the paper does not demonstrate that the outer bins are uncontaminated by the extended core population described by the power-law fit. An independent estimate using control fields away from HII regions would make the claim that \"the issue of background/foreground contamination [is] insignificant\" more convincing.","section":"§3.2"},{"comment":"The significance of the boundary excess is assessed by comparing the counts in the bins 1 ≤ Θ ≤ 2 with the integral of a power law fitted to all bins Θ ≥ 1, including those same boundary bins. The paper itself notes that fitting only Θ ≥ 2 gives a lower curve and a less well-constrained index, so the quoted excess significance is optimistic because the fit absorbs some of the excess. The authors should quote the significance from a fit that excludes the boundary bins, or else correct for this circularity, when claiming \"significantly more than the expected amount\" in Section 3.2.","section":"§3.2"}],"minor_comments":[{"comment":"The citation \"Sandford et al. (1982)\" does not match the reference list entry \"Sandford M.T., Whitaker R.W., & Klein R.I. 1984, ApJ, 282, 178\"; please correct the year or the reference.","section":"§1 and References"},{"comment":"The abstract states that the star formation efficiency of \"the 7 HII region systems with the most reliable mass budgets ranged between 1% and 9%,\" while Section 3.4 and Table 2 list nine systems with SFE below 10%; please reconcile this inconsistency.","section":"§3.4 and Abstract"},{"comment":"The fitted relation N_H2 = (7.14×10^25±7.70×10^25) T^(−3.94±0.43) is given without units or a definition of the variable N_H2; the text should state that this is the average H2 column density in cm^-2 and specify the units of the prefactor.","section":"§3.3"},{"comment":"The text cites \"Rumble et al. (2014)\" for the MCW 297 result, but the reference list contains only Rumble et al. (2015) and Rumble et al. (2016); please add the missing reference or correct the citation.","section":"§4 and References"},{"comment":"The entry \"HII Region ne (cm−3)\" should read \"n_e\" with a subscript for clarity.","section":"Table 1"},{"comment":"The caption says \"Two Histograms\" but only one of the two panels is described in the text as excluding the shell-like HII regions; please clarify the caption to state what each panel shows.","section":"Figure 3 caption"},{"comment":"In the reference list, the Kirk et al. entry incorrectly contains both \"2016\" and \"2014\" before the journal reference; please correct.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is a solid observational study with a valuable dataset, and the projected core distribution measurement is robust. The main obstacle to acceptance is the over-interpretation of the projected profile as a spherical r^-3 volume-density distribution, which the authors themselves know is not unique. I would encourage the editor to request a major revision in which the spherical-symmetry assumption is either tested or prominently caveated, and the self-similarity of the stacking is validated. The inconsistency in the number of SFE systems (seven vs nine) should also be fixed."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bobotsis and Fich stack 53 HII regions from SCUBA-2 archives and find that projected core counts fall as N proportional to Theta^-1.1 beyond the ionized boundary, with an excess near the boundary and a deficit inside. This is a genuinely new population-level result; prior work focused on individual regions or small samples. The stacking method is straightforward, the power-law fit is robust to binning, and removing the two obvious shell-like regions does not kill the boundary excess. The null heating result is also worth attention, though it is limited by the lack of very early O stars and by projection effects.\n\nThe soft spots are in the interpretation. The paper converts the projected counts into a volume density n proportional to r^-3 by Abel inversion, which assumes spherical symmetry about each HII region. The authors acknowledge in Section 4 that a distribution of filament lengths and non-uniform core positions could also fit the counts; with that freedom, the spherical model is not unique. Without kinematic data, the r^-3 law is a plausible interpretation, not a measurement. The abstract states it as a result. The background correction is fitted to the outer tail of the same histogram, but it yields a negligible contamination (2-6 cores), so this is a minor concern. There is also a small internal inconsistency: the abstract says the SFE of the 7 most reliable systems is 1-9%, while the body reports 9 systems with SFE below 10%. And the phrase 'consistent with no cores' in the interior is stronger than the data warrant given possible incompleteness in projection.\n\nNone of these issues undermine the main counting result. The paper is honest about its limitations and suggests the right follow-up (radial velocities). I would send it to a good referee: the stacked radial distribution and boundary excess are worth publishing, and a referee can help clarify the distance between the measured counts and the deprojected profile. The paper will be most useful to people working on HII region feedback and dense core statistics, who can treat the projected counts as solid and the r^-3 law as a working hypothesis. I would cite the paper for the stacking result, but not for the r^-3 law without a caveat. For a reading group, it is a useful case study in how a large archival sample can yield a statistical result and where geometric assumptions enter.","headline":"A novel stacking of 53 HII regions shows a projected core profile and boundary excess that are robust, but the r^-3 volume-density interpretation rests on an untested spherical-symmetry assumption and should be presented as a model, not a measurement.","tokens_in":17162,"tokens_out":3499,"would_cite":true,"duration_ms":33528,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Mature HII regions sit at the centers of dense-core clusters whose volume density falls as distance cubed.","keywords":["HII regions","dense cores","SCUBA-2","submillimeter dust emission","collect and collapse","star formation efficiency","stellar feedback","radial distribution"],"falsifier":"Measure the radial velocities of the dense cores around a sample of these HII regions: a spherical $r^{-3}$ cluster predicts a velocity-dispersion profile that falls with radius and no preferred streaming direction, while a filamentary arrangement predicts coherent velocity gradients along the filaments.","tokens_in":16055,"feed_emoji":"🌟","tokens_out":6786,"duration_ms":61387,"temperature":0.7,"pith_summary":"The paper argues that each mature HII region is at or near the center of an extended cluster of dense cores, and that the cluster extends far beyond the ionized gas. Stacking SCUBA-2 450 and 850 micron images of 53 HII regions yields a projected core count that falls as a power law, which the paper converts into a volume density falling as $r^{-3}$. It also finds an excess of cores just outside the HII region boundary and a near-empty interior. A sympathetic reader would care because this would make the boundary excess a general signature of collect-and-collapse triggering rather than a curiosity of a few shell-like regions, and because the same data argue against significant heating of cores by the HII regions or their OB stars.","feed_headline":"Stacked images place HII regions at centers of core clusters","feed_subtitle":"Core density falls as distance cubed, with an extra shell at the ionized edge, across 53 HII regions","key_machinery":"The load-bearing object is the stacked, radius-scaled separation histogram. Each core's distance from its nearest HII region is divided by that region's radio radius to give $\\Theta_{\\rm scaled}$; individual systems have too few cores for a fit, but when stacked together the bins form a power law. The paper then uses Abel inversion, a standard projection-to-volume conversion, to turn the projected count per ring into a volume number density under spherical symmetry. The same histogram's departure from the power law near $\\Theta=1$ supplies the boundary-excess and interior-deficit results. The scaling by radio radius is what lets heterogeneous HII regions be combined into one statistically meaningful profile.","core_discovery":"On the paper's own terms, the central discovery is that a mature HII region is not at the edge of its parent cloud but at the center of a cluster of dense cores: stacking 315 cores from 53 HII regions and scaling separations by the radio radius gives a projected count $N = (31.6\\pm 7.3)\\Theta^{-1.1\\pm0.2}$ for $\\Theta\\ge1$, which after Abel inversion corresponds to a volume density falling as $r^{-3}$. The histogram also shows an excess of cores in the annulus $1\\le\\Theta\\le2$, with 90 counted versus $70\\pm3$ expected, and this excess survives when the two obvious shell-like regions are removed. Interior to $\\Theta=1$ the counts drop to near zero, consistent with dense cores being destroyed inside the ionized region. A secondary set of results concerns temperature: cloud and core temperatures do not correlate with distance to the HII region or its OB stars, clouds are hotter at low column density, and 74 percent of cores are warmer than their surrounding cloud; the paper interprets this as strong shielding from external heating and as evidence that most cores near HII regions have begun collapse. Star formation efficiencies calculated from dust-based gas masses and a Kroupa IMF are 1 to 9 percent for the seven most reliable systems.","pith_inferences":["A testable consequence the paper does not develop: if the extended core cluster is gravitationally bound, its velocity dispersion should fall with radius; molecular-line radial velocities of cores in a few fields would distinguish this from a filamentary arrangement.","The self-similar scaling assumes all HII regions share the same dimensionless environment; if larger regions have already swept up more of the cluster, the stacked $\\Theta$-profile may mix different evolutionary stages, and splitting by radius or age would reveal how the profile sharpens with time.","The filamentary alternative noted in the paper implies the $r^{-3}$ law is not uniquely determined by the projected counts; combining the SCUBA-2 counts with core velocity or shape data would settle which geometry is real.","If the boundary excess is a general collect-and-collapse signature, then star formation efficiency estimates that ignore the extended core population will overstate the gas reservoir's depletion; including all cores and clumps in the mass budget lowers the efficiency."],"forward_implications":["If the $r^{-3}$ profile is real, mature HII regions trace the centers of dense-core clusters, so massive-star formation in these regions need not be edge-triggered; the OB star sits inside a pre-existing or surviving cluster.","The excess of cores just outside the ionized boundary becomes a general observational signature of collect-and-collapse triggering, applying to most HII regions and not only the few prominent shell-like objects.","The near-empty interior implies dense neutral cores inside mature HII regions should be rare, because cores that start there are destroyed by ionizing radiation.","The lack of temperature-distance correlations implies external heating by the HII region and OB stars is largely shielded by cloud outer layers, so dust-based temperatures trace internal processes more than feedback.","Dust-based total gas mass budgets yield star formation efficiencies of order a few percent, lower than earlier estimates that counted only the most massive clump, so feedback-boosted star formation may be less efficient than previously claimed."],"supporting_citations":[{"why":"Sharpless catalogue defines the sample of mature, optically visible HII regions.","marker":"Sharpless 1959"},{"why":"BFS catalogue supplies additional mature HII regions for the same sample.","marker":"Blitz, Fich, & Stark 1982"},{"why":"VLA radio observations provide positions, sizes, and electron densities used to define each HII region's radius for scaling.","marker":"Fich 1986, 1993"},{"why":"Distances convert angular separations to physical radii, needed for the scaled coordinate and physical densities.","marker":"Foster & Brunt 2015"},{"why":"Dust opacity model converts 450 and 850 micron fluxes into masses and column densities.","marker":"Ossenkopf & Henning 1994"},{"why":"Gould Belt survey recipes supply the photometry and property-calculation methods applied to clouds and cores.","marker":"Buckle et al. 2015"},{"why":"Collect-and-collapse mechanism that the boundary excess is interpreted as evidence for.","marker":"Elmegreen and Lada 1977"},{"why":"Observation that massive dense cores sit near the center of core clusters in Orion B, supporting the central-cluster interpretation.","marker":"Kirk et al. 2016"},{"why":"Independent mass and star formation efficiency budget for the Sh-2 254 complex used to validate the IMF-based method.","marker":"Chavarria et al. 2008"},{"why":"IMF used to extrapolate the stellar mass contributed by low- and intermediate-mass stars in efficiency estimates.","marker":"Kroupa 2001"}],"fun_headline_variants":["HII regions sit at centers of dense core clusters","Dense core density falls as r^-3 around HII regions","Cores cluster around HII regions, but not inside them","HII regions shielded cores heat up, spark star formation","HII regions at heart of core clusters, density drops r^-3"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The conversion of the observed projected counts into an $r^{-3}$ volume density assumes the cores are distributed spherically around each HII region; the paper acknowledges that filamentary structures with adjustable lengths and core positions could also reproduce the counts.","fun_headline_variants_meta":{"raw":{"variants":["HII regions sit at centers of dense core clusters","Dense core density falls as r^-3 around HII regions","Cores cluster around HII regions, but not inside them","HII regions shielded cores heat up, spark star formation","HII regions at heart of core clusters, density drops r^-3"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000333,"raw_usage":{"total_tokens":1886,"prompt_tokens":1019,"completion_tokens":867,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":635,"completion_tokens_details":{"reasoning_tokens":781}},"tokens_in":635,"tokens_out":867,"duration_ms":7712,"temperature":1.0,"reasoning_tokens":781,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:07:54.893510+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the radial velocities of the dense cores around a sample of these HII regions: a spherical $r^{-3}$ cluster predicts a velocity-dispersion profile that falls with radius and no preferred streaming direction, while a filamentary arrangement predicts coherent velocity gradients along the filaments.","supporting_citations":[{"cited_title":"1959, ApJS, vol","cited_arxiv_id":null,"evidence_quote":"Sharpless catalogue defines the sample of mature, optically visible HII regions."},{"cited_title":"1982, A&AS, 49, 183","cited_arxiv_id":null,"evidence_quote":"BFS catalogue supplies additional mature HII regions for the same sample."},{"cited_title":"1986, AJ, 92, 787","cited_arxiv_id":null,"evidence_quote":"VLA radio observations provide positions, sizes, and electron densities used to define each HII region's radius for scaling."},{"cited_title":"2015, AJ, 150, 13","cited_arxiv_id":null,"evidence_quote":"Distances convert angular separations to physical radii, needed for the scaled coordinate and physical densities."},{"cited_title":"1994, A&A, 291, 943","cited_arxiv_id":null,"evidence_quote":"Dust opacity model converts 450 and 850 micron fluxes into masses and column densities."},{"cited_title":"2015, MNRAS, 449, 2472","cited_arxiv_id":null,"evidence_quote":"Gould Belt survey recipes supply the photometry and property-calculation methods applied to clouds and cores."},{"cited_title":"1977, ApJ, 214, 725","cited_arxiv_id":null,"evidence_quote":"Collect-and-collapse mechanism that the boundary excess is interpreted as evidence for."},{"cited_title":"2016, 2014, ApJ, 821, 18","cited_arxiv_id":null,"evidence_quote":"Observation that massive dense cores sit near the center of core clusters in Orion B, supporting the central-cluster interpretation."},{"cited_title":"2008, ApJ, 682, 445","cited_arxiv_id":null,"evidence_quote":"Independent mass and star formation efficiency budget for the Sh-2 254 complex used to validate the IMF-based method."}],"review_version":1}