{"id":"8f73b22a-209f-44a7-9358-ef39407ea1fd","arxiv_id":"2507.22482","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"Density-based clustering of molecular emission features reproduces known c-C3H2/CH3OH segregation and claims a new c-C3H2/CH3CCH segregation, but the new claim is weakened by biased input sampling and missing significance tests.","lead":"This paper applies DBSCAN and HDBSCAN clustering to molecular emission maps of three starless cores, using physical features such as intensity and velocity instead of spatial coordinates. It claims to reveal a chemical segregation between c-C3H2 and CH3CCH that is not visible by eye.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Case 2 clusters may be an artifact of duplicated CH3CCH transitions per pixel; an equal-per-pixel control is needed before the claimed segregation is accepted.","rationale":"The reader's weakest-assumption analysis identifies the same load-bearing flaw: the Case 2 input data duplicate CH3CCH transitions at each pixel in B68 and L1521E, inflating local density and biasing the density-based clustering. This is not a peripheral methodological quibble; it directly controls whether the paper's headline discovery, the c-C3H2/CH3CCH segregation that is 'not apparent' from the maps, is real. The paper's own Section 5.1 admits the duplication raises the likelihood of clusters at CH3CCH-rich locations, but it does not run the decisive control. Because a REJECT verdict was already issued on this basis, my read does not move the verdict. The concrete test proposed here would settle the matter: if the segregation vanishes under equal per-pixel sampling, the manuscript should be rejected or substantially revised; if it survives, the central claim would need to be reassessed in light of the control.","tokens_in":38470,"tokens_out":2714,"duration_ms":35817,"concrete_test":"Re-run Case 2 exactly as in Section 3.2 and 3.3, but collapse the CH3CCH data to one point per spatial pixel, for example by selecting the highest-S/N transition or taking the mean of the normalized transition intensities, while keeping c-C3H2 at one point per pixel in B68 and L1521E. Re-optimize DBSCAN and HDBSCAN hyperparameters using the same DBCV and coverage criteria, then compare cluster molecular ratios and spatial maps. If the imbalanced c-C3H2/CH3CCH clusters disappear or lose their spatial separation, the claimed segregation is an artifact of duplicated transitions; if the clusters persist, the central claim survives this specific concern.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central new claim, that clustering reveals a segregation between c-C3H2 and CH3CCH in all three cores, rests entirely on Case 2. In B68 and L1521E, that dataset includes two or three CH3CCH transitions per spatial pixel while c-C3H2 contributes only one point per pixel (Table 4). Because DBSCAN and HDBSCAN are density-based, these duplicated CH3CCH points artificially raise the local point density in feature space, preferentially creating clusters in CH3CCH-bright regions. The paper itself acknowledges this mechanism in Section 5.1: it increases the likelihood of forming clusters at locations with enhanced CH3CCH data density. The load-bearing assumption is that these duplicate line measurements are independent physical samples. No test is provided that controls for per-pixel sampling, so the reported segregation in B68 and L1521E, and the interpretation that CH3CCH traces an inner layer while c-C3H2 traces an outer shell, are not established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper applies the density-based clustering algorithms DBSCAN and HDBSCAN to IRAM 30 m observations of c-C3H2, CH3OH, and CH3CCH toward the starless cores B68 and L1521E and the prestellar core L1544. The spatial coordinates of each emission pixel are discarded and replaced by six physical features (integrated intensity, velocity offset, linewidth, H2 column density, distance to the dust peak, and H2 column density gradient). Four case studies are considered: c-C3H2 versus CH3OH (Case 1), c-C3H2 versus CH3CCH (Case 2), CH3OH versus CH3CCH (Case 3), and all three molecules together (Case 4). The analysis reproduces the known c-C3H2/CH3OH segregation in all three cores and, as the main new result, claims to identify a segregation between c-C3H2 and CH3CCH that is not apparent from the emission maps. The authors also measure CH3CCH abundances at the dust peaks of several cores and use a gas-grain chemical model to argue that the CH3CCH peak in L1544 traces the landing point of freshly accreted gas.","tokens_in":38568,"tokens_out":6170,"duration_ms":69624,"significance":"If the central new claim is robust, the paper demonstrates that density-based clustering on small molecular datasets can uncover subtle chemical segregation, and the c-C3H2/CH3CCH segregation in the less-evolved cores B68 and L1521E would be an observationally new result with implications for chemical layering and accretion flows. The successful reproduction of the previously known c-C3H2/CH3OH segregation serves as a useful positive control, and the measured evolutionary trend of CH3CCH abundances from starless to prestellar cores is an interesting addition. The authors also make their detailed clustering figures publicly available on Zenodo, which improves reproducibility. However, the principal new claim rests on a sampling assumption that is acknowledged but not tested, and the imbalance criterion used throughout lacks a statistical significance test; these issues need to be resolved before the result can be accepted.","major_comments":[{"comment":"The central new claim of a c-C3H2/CH3CCH segregation in B68 and L1521E is not established because Case 2 includes two or three CH3CCH transitions per spatial pixel while c-C3H2 contributes only one point per pixel. Since DBSCAN and HDBSCAN are density-based, these duplicated CH3CCH points inflate the local density in feature space at CH3CCH-bright positions, biasing cluster formation toward those regions. The paper acknowledges this effect in §5.1 but does not provide an equal-per-pixel control (e.g., using a single representative transition per molecule, or averaging transitions, or down-weighting duplicated points). Until such a control is performed, the reported segregation and the interpretation that CH3CCH traces an inner layer while c-C3H2 traces an outer shell are not robust against this sampling artifact.","section":"§3.3, Table 4; §5.1"},{"comment":"The definition of an 'imbalanced' cluster as a deviation of at least 10% from the initial molecular ratio is not accompanied by any significance test. Several clusters contain very few points (N of order 6-10; e.g., Table C.1, combi 2, L1521E, cluster 2 and combi 8, B68, cluster 4), where a 10% shift can be consistent with Poisson counting noise. The paper should either restrict the imbalance designation to clusters whose ratios differ significantly from the input ratio (e.g., via a chi-square or permutation test) or state the expected binomial scatter around the input ratio. This matters because the qualitative summaries in §4.3.2 and the final conclusions on molecular segregation rely on these labels.","section":"§4.3, Table C.1"},{"comment":"The clustering treats each emission-pixel sample as an independent data point even though adjacent pixels are strongly correlated given the 8 arcsec pixel size and 32 arcsec beam, and it treats multiple transitions of the same molecule as independent feature vectors. This pseudo-replication can affect the local density estimates that DBSCAN and HDBSCAN rely on, particularly when transitions of one molecule are duplicated in Case 2 and Case 3. The manuscript should discuss this limitation explicitly and, ideally, include a test in which the data are smoothed or thinned so that each independent beam is represented by a single sample per molecule.","section":"§3.1, §3.3"}],"minor_comments":[{"comment":"The abstract and conclusions list integrated intensity, velocity offset, H2 column density, and H2 column density gradient as the key features driving the clustering, but the paper does not present a quantitative feature-relevance analysis; this claim appears to be based on visual inspection of two-dimensional feature projections. Please clarify the basis for this statement or add a feature-importance measure.","section":"Abstract, §5.1"},{"comment":"The excitation temperature is fixed at 8 K for all cores, and the statement that a lower or higher excitation temperature only shifts the abundances without changing the overall trend is not quantified. Please provide the magnitude of the shift for a reasonable range (e.g., 5-10 K) to support the robustness claim.","section":"§4.4"},{"comment":"The transition labels are inconsistent: Table 2 lists CH3OH 20,2-10,1 (E2) and 21,2-11,1 (A+), while Fig. A.1 labels the maps '202-101E0' and '212-111E2'. Please harmonize the notation.","section":"Table 2, Fig. A.1"},{"comment":"The caption of Fig. 2 states that dashed contours represent 30%, 50%, and 90% of the H2 column density peak, but the citation to Spezzano et al. (2020) appears after the period in an awkward way. Please fix the sentence structure.","section":"Fig. 2 caption"},{"comment":"Table 4 gives the molecular ratios as percentages but does not list the total number of data points per dataset. Adding N would help the reader evaluate the statistical weight of the clusters shown in Table C.1.","section":"Table 4"},{"comment":"For Case 4, the text says the dataset combines Case 1 and Case 2, but it is not stated explicitly which transitions are used for each molecule in that combined dataset. Please specify this to avoid ambiguity, especially because the number of CH3CCH transitions differs between cores.","section":"§3.3"}],"recommendation":"major_revision","confidential_remarks":"The reader's report is largely fair: the duplicated CH3CCH transitions in Case 2 are a genuine, acknowledged, and load-bearing issue for the paper's main claim, and the absence of a significance test for cluster-ratio imbalances is also a real weakness. However, I do not see these as irreparable. Both problems can be addressed with additional control experiments and statistical tests that fall within the scope of the manuscript. The positive control of the known c-C3H2/CH3OH segregation and the abundance trend for CH3CCH suggest the underlying data and methods are suitable for A&A. I therefore recommend major revision rather than rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the method demonstration is worth a look, but the paper's headline result—a new segregation between c-C3H2 and CH3CCH—doesn't hold up under the current input design. The authors themselves acknowledge the mechanism in Section 5.1: in B68 and L1521E they feed two or three CH3CCH transitions per pixel against one c-C3H2 point, so CH3CCH gets a density bonus in DBSCAN/HDBSCAN. No control is shown that the clustering survives when each pixel contributes one point per molecule. For a density-based algorithm, this is a load-bearing flaw, not a cosmetic one.\n\nWhat's good: the reproduction of the known c-C3H2/CH3OH segregation acts as a useful positive control; the feature engineering (physical environment rather than position) is sensible and genuinely helpful for cross-core comparison; the preprocessing is careful; the maps and data are available. The chemical discussion about CH3CCH in L1544 is plausible but secondary—it rests on an assumed Tex and a qualitative model comparison, and it is not needed to assess the clustering claim.\n\nWhere the soft spots are: (1) the duplicated transitions issue above; (2) the 10% imbalance threshold for cluster significance has no statistical justification—there's no test that the molecular ratios differ from the input ratio by more than sampling noise; (3) the claim that the c-C3H2/CH3CCH segregation is 'not apparent' from the maps is not credible. Figures A.1/A.2 show CH3CCH is clearly less extended and peaks at different positions in B68 and L1521E. The clustering may be useful as a way to quantify that morphology, but calling it a new, hidden segregation overstates it.\n\nThe reader's verdict of reject is harsh if read as 'the method has no value'; as a statement about the central discovery, it is fair. The paper deserves a serious referee, because the method—with an equal-per-pixel control and proper significance tests—could support a modest version of the claim. My recommendation: send to review, but the referee should ask for the control and a softening of the conclusions.","headline":"The method demonstration is solid as a positive control, but the claimed new c-C3H2/CH3CCH segregation is likely an artifact of duplicated transitions per pixel and does not survive scrutiny.","tokens_in":39244,"tokens_out":1992,"would_cite":false,"duration_ms":23336,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Density-based clustering of molecular-line pixels, with spatial coordinates removed, reveals a chemical segregation between c-C3H2 and CH3CCH in all three cores that is not visible in the emission maps.","keywords":["chemical segregation","dense cores","DBSCAN","HDBSCAN","c-C3H2","CH3CCH","molecular line emission","H2 column density gradient"],"falsifier":"The decisive test is a clustering rerun on B68 and L1521E with one data point per molecule per pixel, for example averaging the CH3CCH transitions before building the feature space; if the previously imbalanced clusters persist, the segregation is chemical, and if they vanish, it is an artifact of duplicated points.","tokens_in":38177,"feed_emoji":"🌌","tokens_out":9881,"duration_ms":98781,"temperature":0.7,"pith_summary":"The paper sets out to show that unsupervised density-based clustering (DBSCAN and HDBSCAN) can uncover chemical structure in dense cores from very small molecular datasets. Each emission pixel is described by its physical environment—integrated intensity, velocity offset, linewidth, H2 column density, distance to the dust peak, and column-density gradient—rather than by its position on the sky. The analysis reproduces the known segregation between c-C3H2 and CH3OH in B68, L1521E, and L1544, and it further claims a new segregation between c-C3H2 and CH3CCH in all three cores, a separation that is not apparent from the emission maps themselves. The authors interpret this as evidence that c-C3H2 traces lower-density outer gas while CH3OH and CH3CCH trace inner or accreting gas, and they use the clustering plus chemical modelling to argue that the CH3CCH peak in L1544 marks the landing point of chemically fresh accreted material.","feed_headline":"Clustering reveals hidden chemical split in starless cores","feed_subtitle":"Density-based grouping separates c-C3H2 from CH3CCH, pointing to layered chemistry invisible in emission maps.","key_machinery":"The machinery is density-based clustering—DBSCAN, which groups points whose mutual distance is below a threshold epsilon and labels smaller groupings as noise, and HDBSCAN, its hierarchical extension that handles clusters of varying density. The input is constructed by removing spatial coordinates and representing each observed pixel by six physical features: integrated intensity, velocity offset from the source's systemic velocity, linewidth, H2 column density, projected distance to the dust peak, and the magnitude of the H2 column-density gradient. This feature-space representation lets the same analysis run across cores of different sizes and locations, and the clusters are then read by looking at which molecule dominates each imbalanced cluster. Chemical simulations of the CH3CCH formation and destruction network (formation via dissociative recombination of C3H5+ and destruction by atomic carbon) provide the mechanism invoked to explain why the molecule survives at the L1544 accretion landing point.","core_discovery":"At the paper's center is a claim about observable chemistry: in the starless cores B68 and L1521E and the prestellar core L1544, the carbon-chain molecules c-C3H2 and CH3CCH occupy different physical layers even where their projected emission maps overlap. The argument is made by combining both molecules without labels in a feature space built from each pixel's intensity, velocity offset, linewidth, H2 column density, distance to the dust peak, and column-density gradient; clusters that come out imbalanced in one molecule correspond to regions of chemical segregation. This approach reproduces the already known c-C3H2/CH3OH segregation, validates the method, and then reveals a c-C3H2/CH3CCH segregation in all three cores. The paper also reports that CH3OH and CH3CCH cluster similarly relative to c-C3H2, that the most informative features are intensity, velocity offset, and column-density-related quantities, and that in L1544 the CH3CCH peak sits at a shielded northwest location where fresh accreted gas can form CH3CCH before atomic carbon destroys it. Abundance measurements at dust peaks add an evolutionary reading: CH3CCH is about one order of magnitude more abundant in starless than in prestellar cores, with L1544 as the exception.","pith_inferences":["A natural robustness check: repeat the clustering with one data point per molecule per pixel in B68 and L1521E to test whether the new c-C3H2/CH3CCH segregation survives removal of duplicated transitions; the paper acknowledges the duplication but does not perform this test.","The same feature-space recipe could be applied to archival multi-molecule maps of other cores to hunt for hidden segregation wherever projected emission overlaps, giving a census of chemical layers across many clouds.","If the accretion interpretation holds, c-C3H2-to-CH3CCH or c-C3H2-to-CH3OH ratio maps might serve as a practical tracer of inflowing, chemically fresh gas, complementing kinematic methods.","The prominence of velocity offset suggests a testable dynamical prediction: synthetic observations from collapsing-core models with time-dependent chemistry should reproduce the same cluster separations only when accretion and photochemistry are included."],"forward_implications":["Dense cores are chemically layered in a way that single-molecule maps obscure: c-C3H2 traces a lower-density outer shell while CH3CCH and CH3OH concentrate in inner or freshly accreted gas.","Chemical-differentiation studies need not wait for large line surveys; a handful of targeted molecular transitions can expose segregation with density-based clustering.","Because velocity offset dominates the cluster splits, static chemical models are inadequate for anisotropic chemical structures, and dynamical models including accretion must be used.","CH3CCH abundance, about one order of magnitude lower in prestellar than starless cores except at L1544, can serve as a combined probe of evolutionary stage and local accretion environment.","Using the H2 column density gradient as a feature in place of sky coordinates lets cores in different environments be compared directly, and the varying relevance of this feature across cores reflects different illumination conditions."],"supporting_citations":[{"why":"Supplies the DBSCAN algorithm used for all clusterings.","marker":"Ester et al. (1996)"},{"why":"Supplies the HDBSCAN algorithm and its density-based validation approach.","marker":"Campello et al. (2013)"},{"why":"Provides the L1544 c-C3H2 and CH3OH data and the known segregation that the method must reproduce.","marker":"Spezzano et al. (2016)"},{"why":"Identifies the molecular families in L1544, including the anomalous CH3CCH peak that the clustering analysis explains.","marker":"Spezzano et al. (2017)"},{"why":"Provides the L1521E data and reports the c-C3H2/CH3OH differentiation used as a baseline.","marker":"Nagy et al. (2019)"},{"why":"Supplies B68 and L1521E data and the H2 column density maps used to build the clustering features.","marker":"Spezzano et al. (2020)"},{"why":"Shows that c-C3H2 traces lower-density gas, the interpretation applied to the new segregation.","marker":"Lin et al. (2022)"},{"why":"Provides the chemical network used to model CH3CCH formation and destruction in L1544.","marker":"Wakelam et al. (2015)"},{"why":"Provides the Gaussian-derivative method used to compute the H2 column-density gradient feature.","marker":"Soler et al. (2013)"},{"why":"Provides the column-density formula used for the abundance comparison at dust peaks.","marker":"Mangum & Shirley (2015)"}],"fun_headline_variants":["Clustering finds new chemical split in starless cores","Unsupervised clustering reveals hidden core chemistry","Starless cores: clustering spots unseen molecular layers","Algorithm exposes distinct chemical zones in starless cores"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result depends on treating each spectral transition of a molecule as an independent sample, which in B68 and L1521E duplicates CH3CCH at the same pixels and can create clusters through point density rather than physical chemistry—a choice the paper flags but does not correct.","fun_headline_variants_meta":{"raw":{"variants":["Clustering finds new chemical split in starless cores","Unsupervised clustering reveals hidden core chemistry","Starless cores: clustering spots unseen molecular layers","Algorithm exposes distinct chemical zones in starless cores"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0004,"raw_usage":{"total_tokens":2206,"prompt_tokens":1175,"completion_tokens":1031,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":791,"completion_tokens_details":{"reasoning_tokens":972}},"tokens_in":791,"tokens_out":1031,"duration_ms":11406,"temperature":1.0,"reasoning_tokens":972,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T11:37:38.224911+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"The decisive test is a clustering rerun on B68 and L1521E with one data point per molecule per pixel, for example averaging the CH3CCH transitions before building the feature space; if the previously imbalanced clusters persist, the segregation is chemical, and if they vanish, it is an artifact of duplicated points.","supporting_citations":[{"cited_title":"1996, Proceedings of the Second International Conference on Knowledge Discovery and Data Mining, KDD-96 (AAAI Press), 226","cited_arxiv_id":null,"evidence_quote":"Supplies the DBSCAN algorithm used for all clusterings."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the HDBSCAN algorithm and its density-based validation approach."},{"cited_title":"E., et al","cited_arxiv_id":null,"evidence_quote":"Supplies B68 and L1521E data and the H2 column density maps used to build the clustering features."},{"cited_title":"C., Herbst , E., et al","cited_arxiv_id":null,"evidence_quote":"Provides the chemical network used to model CH3CCH formation and destruction in L1544."},{"cited_title":"D., Hennebelle , P., Martin , P","cited_arxiv_id":null,"evidence_quote":"Provides the Gaussian-derivative method used to compute the H2 column-density gradient feature."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the column-density formula used for the abundance comparison at dust peaks."}],"review_version":1}