{"id":"4cdd7973-2dc0-453e-ab18-11f5e09a1932","arxiv_id":"2607.27382","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Spectral-line stacking can recover aggregate 13C-isotopologue emission without false positives, and existing TMC-1 observations could distinguish a local from an enriched carbon-isotope ratio in PAHs if the needed lab spectra existed.","lead":"This paper tests whether radio astronomers can add together many faint spectral lines from molecules differing only in containing the rare carbon-13 isotope, and whether this can create false detections. It finds the method is safe and that existing telescope data on the cold cloud TMC-1 could distinguish locally formed from inherited aromatic molecules once laboratory spectra of the isotopic variants exist.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Unvalidated 13C spectral predictions make the claimed 12C/13C≈30 detection threshold, and hence the discrimination forecast, an open question.","rationale":"The paper is a clean proof-of-concept: the false-positive/false-negative tests in §3.1 are convincing, and the HC9N demonstration uses real laboratory spectra, showing that the stacking/matched-filter method works when line lists are known. The unresolved step is the extrapolation to cyanonaphthalene/cyanopyrene 13C isotopologues, whose laboratory spectra do not yet exist. The simulations in §3.3 measure the detectability of a signal generated from the calculated line lists. They do not measure the effect of any error in those lists. The paper asserts the scaled constants are 'sufficiently accurate' but provides no quantitative error analysis and no test against measured 13C spectra for a related molecule. Frequency errors are the most dangerous because the matched-filter stack relies on precise alignment; a systematic offset of even a few tenths of a km/s would reduce the coherent SNR and could move the 5σ threshold upward. The central claim is explicitly conditional on laboratory spectra, so a skeptical reader might say the claim will be settled only when those spectra are measured. But the paper's purpose is to motivate those laboratory efforts; for that, an unvalidated proxy is weak support. The proposed concrete test using HC9N (or benzonitrile, where 13C data may also exist) would directly quantify the scaling method's accuracy and allow a corrected sensitivity forecast. This is the single check that would decide whether the discrimination claim is robust. The reader's weakest_assumption flagged the same frequency-accuracy issue, and also mentioned the hidden-line power-law extrapolation; we agree that the frequency accuracy is the more load-bearing of the two because it affects the very signal being searched for. Therefore our verdict remains CONDITIONAL, consistent with the reader's.","tokens_in":17761,"tokens_out":8982,"duration_ms":114245,"concrete_test":"Take a molecule with known laboratory spectra for all singly substituted 13C isotopologues (e.g., HC9N from McCarthy et al. 2000). Apply the exact §2.3.1 scaling method (M06-2X/6-31+G(d) parent-constant scaling, unchanged distortion/hyperfine) to predict its 13C isotopologue spectra. Measure the rms and maximum velocity offsets between predicted and measured line frequencies. Then re-run the §3.3 signal-injection pipeline, offsetting each injected PAH line by an error drawn from this distribution, and recompute the 12C/13C ratio at which the median matched-filter response first exceeds 5σ. If the threshold moves above ~69, the discrimination claim fails. If the offsets are negligible relative to the 0.4 km/s FWHM and the threshold remains below ~35, the forecast is supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The Abstract's central forecast is that, with laboratory rotational spectra in hand, current GOTHAM data could discriminate a bulk 12C/13C~69 from an enriched ratio below ~30. To support this, §3.3 injects aggregate 13C-isotopologue signals into synthetic noise and recovers them with the same calculated line lists. This is a self-consistency test, not a test of the fidelity of those line lists. The predicted spectra rest on parent-constant scaling of M06-2X/6-31+G(d) rotational constants and the assumption that centrifugal distortion and 14N hyperfine constants are unchanged from the parent (§2.3.1). No validation against experimental 13C spectra of any test molecule is provided. Residual errors in scaled rotational constants for a ~170 amu molecule could be ~1 MHz; at 10 GHz this is a ~30 km/s velocity offset, orders of magnitude larger than the 0.4 km/s line width used in the stacks. Even if errors are an order of magnitude smaller, systematic offsets change the co-added line pattern and the SNR of the matched-filter response. The threshold ratio at which 5σ is reached could shift upward, potentially above the local bulk ratio of 69, invalidating the claimed discrimination. Because the same line list is used for injection and extraction, the simulations cannot reveal this failure mode. The paper explicitly defers to future laboratory data, but the conditional claim is meaningful only if the proxy line lists are representative; that representativeness is currently unestablished.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper develops and stress-tests a spectral line stacking / matched-filtering method for recovering aggregate emission from many singly substituted 13C isotopologues in the GOTHAM TMC-1 survey. Synthetic injection-recovery experiments with n=100 trials per line density show that the matched-filter response is centered near 0σ for absent signals and degrades to ~1σ under complete line confusion, implying that buried lines cannot inflate false positives above the adopted 5σ threshold (§3.1, Eqs. 1–2). The method is demonstrated on HC9N isotopologues in real GOTHAM data, giving a 14.5σ aggregate stack, with a rough 12C/13C ratio of 130 (§3.2). The main forecast is for cyanonaphthalene and cyanopyrene: using quantum-chemically predicted 13C spectra, the authors simulate that at 12C/13C≲30 current observations would yield >5σ aggregate detections, whereas at the local bulk ratio ~69 they would not, potentially distinguishing in-situ from inherited PAH formation (§3.3).","tokens_in":18113,"tokens_out":5135,"duration_ms":49026,"significance":"If the forecast is correct, the paper provides a clear observational route to discriminating in-situ from inherited PAH formation and justifies costly laboratory measurements. Its strengths are the controlled injection-recovery design with known ground truth, the large number of Monte Carlo trials, the explicit false-positive/false-negative analysis, and the candid discussion of caveats (power-law extrapolation, need for laboratory spectra). However, the headline discrimination claim is only as strong as the unvalidated quantum-chemical 13C line lists and the circular HC9N calibration; these are load-bearing and need further work or qualification.","major_comments":[{"comment":"The central PAH forecast is a self-consistency test: the same quantum-chemically predicted 13C line lists (M06-2X/6-31+G(d), parent-constant scaling, unchanged centrifugal/hyperfine constants; §2.3.1) are used both to inject and to recover signal in §3.3. No validation against experimental 13C spectra is provided. The text itself concedes that 'laboratory accuracy will still be required.' How large a frequency error is tolerable? A 1 MHz offset at 10 GHz is ~30 km/s, two orders of magnitude larger than the 0.4 km/s line width; even 0.1 MHz offsets would materially degrade the matched-filter response. Because the same line list appears on both sides, the simulated >5σ thresholds at 12C/13C≈30 cannot establish that real laboratory spectra would produce a detection. Please add a sensitivity analysis over plausible scaling errors, or a validation against measured 13C spectra of a known molec","section":"§3.3, §2.3.1"},{"comment":"The HC9N application is partly circular. The aggregate simulation's 12C/13C ratio is 'empirically scaled ... to a value of 130 to match the spectral stack intensity by eye' (§3.2), and the Conclusions then state 'the lower limit of the carbon isotope ratio for HC9N estimated from our presented analysis is consistent with previous estimates.' A value adjusted to fit the stack cannot be independently validated with that same stack, and no uncertainty is quoted. As a proof of aggregate recovery the section is convincing; as an isotopic constraint it is not. Either remove the interpreting sentence and present the 130 as a calibration parameter, or incorporate the forthcoming MCMC fit (Burkhardt 2026) with quoted uncertainties.","section":"§3.2, Conclusions"},{"comment":"The paper is admirably explicit that the hidden-line density estimate is 'highly speculative' and that the power-law extrapolation has 'no compelling physical reason.' This caveat matters because vertical line A in Fig. 2 and the claim that GOTHAM's SNR is 'slightly depressed' depend on that extrapolation. The controlled injection-recovery experiments (n=100) stand on their own, but the application to GOTHAM's actual line-confusion regime is model-dependent. The text should mark the GOTHAM-specific depression claim as an illustration under an assumed line density, not a measurement.","section":"§3.1.3, Appendix A"}],"minor_comments":[{"comment":"The text refers to 'all four isomers' of cyanopyrene but enumerates three (1-, 2-, and 4-cyanopyrene) and the counts 17+17+11=45. Please correct 'four' to 'three'.","section":"§3.3.2"},{"comment":"The symbols S, N_c, N_w are defined in the sentence after the equation; moving these definitions before the equation would improve readability.","section":"Eq. (1)"},{"comment":"The expression 'n σσ local,i' appears to contain a typographical duplication of σ; it should likely be nσ × σ_local,i.","section":"Appendix A"},{"comment":"The caption describes the pink dataset as 'showing false negative detection rates' but does not explicitly say that the orange dataset shows false positives; consider making this explicit for clarity.","section":"Figure 2 caption"}],"recommendation":"major_revision","confidential_remarks":"This is a well-executed proof-of-concept whose headline claim is gated by the unvalidated quantum-chemical 13C line lists. The stress-test concern is real: the same line lists are used for injection and recovery, so the simulated 5σ thresholds cannot speak to the fidelity of the predicted spectra. The HC9N 'by eye' scaling in §3.2 is not a rigorous isotopic constraint and should not be interpreted as one. Both issues are fixable within the manuscript's scope by adding a sensitivity analysis or explicit validation, and by softening the abstract/conclusions to match what the simulations actually establish."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth knowing: this paper doesn't claim a detection. It claims that if you had lab spectra of all singly substituted 13C isotopologues of the cyanonaphthalenes and cyanopyrenes, the GOTHAM data could discriminate bulk vs enriched 12C/13C. That's a conditional forecast, and the load-bearing assumption is that quantum-chemically scaled line lists are representative. That assumption is never tested. So read the '12C/13C~30' threshold as a placeholder until lab spectra arrive.\n\nWhat's genuinely new: the systematic characterization of how line density pushes stacking and matched filtering toward false negatives (and why it can't manufacture false positives), and the aggregate multi-isotopologue recovery idea. The injection-recovery tests are clean: 100 trials, known ground truth, and the analytic SNR→1 in the fully confused limit is a nice sanity check. The HC9N demonstration is the right way to show the aggregate stacking idea works when lab line lists actually exist. The paper is also honest about the biggest weaknesses: the power-law line density extrapolation is flagged as having no physical basis, and the by-eye HC9N scaling is disclosed.\n\nSoft spots, in order. First, the HC9N 'lower limit of 130' is not independent. They scaled the aggregate simulation to match the stack by eye, then the conclusions use that same number as evidence consistent with in-situ formation. They say a proper MCMC fit is coming (Burkhardt 2026), but as written the conclusion outruns the evidence. The fix is easy: present the by-eye value as a tuning parameter, not a measurement.\n\nSecond, and more important, the PAH sensitivity forecast is a self-consistency test. The same scaled line lists are used for injection and extraction, so any systematic frequency error cancels out. Errors in the M06-2X scaled constants could plausibly be ~1 MHz, which at 10 GHz is a ~30 km/s velocity offset, orders of magnitude above the 0.4 km/s line width. Even if the errors are ten times smaller, the co-added response would smear. The paper defers to future lab spectra, which is reasonable, but the abstract's discrimination claim should be framed more carefully as 'assuming the line lists are accurate.' The stress-test note lands: unvalidated spectral predictions make the threshold an open question, not a result.\n\nMinor: no code or data release. For a methods paper, that would help.\n\nBottom line: this is a solid proof-of-concept for the false-positive/negative characterization, and the HC9N demonstration is credible. The PAH forecast is a well-posed thesis that needs lab validation. It deserves a serious referee, with the expectation of major revision: MCMC for HC9N, a quantitative error budget for the scaled spectra, and an abstract that matches the actual conditionality.","headline":"Clean stress-test of stacking methods, but the PAH sensitivity forecast rests on unvalidated line lists — conditional proof-of-concept, not a result.","tokens_in":18665,"tokens_out":2787,"would_cite":true,"duration_ms":26543,"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":"This proof-of-concept shows that stacking the predicted lines of all singly 13C-substituted isotopologues of a PAH can recover a single aggregate signal from TMC-1 data, and that existing observations would separate a locally formed aromati","keywords":["spectral line stacking","matched filtering","13C isotopologues","polycyclic aromatic hydrocarbons","TMC-1","interstellar chemistry","radio astronomy","carbon isotope ratios"],"falsifier":"Measure the laboratory rotational spectrum of a single singly 13C-substituted isotopologue of 1-cyanonaphthalene, shift it to TMC-1's velocity, and run it through the paper's stacking-and-filtering pipeline in place of the predicted spectrum; if the recovered impulse response drops below the simulated 5σ threshold at 12C/13C ≈ 30, the central forecast fails.","tokens_in":17636,"feed_emoji":"📡","tokens_out":5582,"duration_ms":48309,"temperature":0.7,"pith_summary":"This paper asks whether spectral line stacking and matched filtering can be pushed to recover, in aggregate, the emission of all singly 13C-substituted isotopologues of a polycyclic aromatic hydrocarbon. The authors argue that the answer is yes: the combined signal of many weak isotopologue lines can be retrieved as one stacked spectrum, buried interfering lines depress rather than inflate the recovered significance, and false positives are effectively impossible at the adopted 5σ threshold. They demonstrate the approach on real GOTHAM observations of HC9N, recovering a strong aggregate signal, and then use synthetic data to show that if laboratory rotational spectra of the 13C isotopologues of cyanonaphthalene and cyanopyrene existed, current TMC-1 observations would distinguish a locally formed aromatic population (12C/13C near the bulk value of about 69) from one enriched in 13C by circumstellar inheritance (ratio below about 30). If right, this gives a direct isotopic test of where interstellar PAHs come from.","feed_headline":"Stacking faint lines can unmask 13C-rich PAHs in TMC-1","feed_subtitle":"Aggregate isotopologue stacking recovers a 13C signal and would flag circumstellar inheritance at 12C/13C below ~30.","key_machinery":"Spectral line stacking and matched filtering: relevant rotational transitions of a molecule (or of all isotopologues summed into one aggregate spectrum) are extracted from a broadband survey, re-gridded into velocity space with line centers aligned, and SNR-weighted; the stacked simulated spectrum then serves as a matched filter cross-correlated with the stacked data. The load-bearing statistic is the apparent SNR, SNR_app = (S + N_c)/(N_w + N_c), which in the limit of strong line confusion approaches unity, so buried lines can suppress but never manufacture a 5σ detection. The 5σ impulse-response threshold is the adopted detection criterion.","core_discovery":"The central claim is a conditional sensitivity forecast: if laboratory rotational spectra of all singly 13C-substituted isotopologues of the cyanonaphthalene and cyanopyrene isomers were available, the GOTHAM observations of TMC-1 would already be sensitive enough to discriminate, at a 5σ matched-filter threshold, between a local, bulk 12C/13C ratio of about 69 and an enhanced 13C abundance with a ratio at or below about 30. The paper also establishes, through synthetic line-confusion stress tests, that stacking many predicted lines into one matched filter cannot be fooled into a false positive by hidden sub-noise lines—the recovered significance instead collapses toward 1σ as line density a","pith_inferences":["The test is asymmetric in a way the authors note but do not emphasize: an enriched-13C detection cleanly supports inheritance, but a bulk-ratio detection does not cleanly rule it out, because circumstellar PAHs need not always be 13C-enriched; the useful null result is the non-detection at the enrichment threshold.","The forecast's weakest link—the accuracy of the scaled quantum-chemical spectra—can be tested incrementally: measuring even one 13C isotopologue in the laboratory and checking that its strongest lines fall inside the matched-filter passband would validate or invalidate the entire sensitivity argument before large observational effort is spent.","The same aggregate-stacking logic should transfer to other isotopic substitutions (15N, D) and to families of related molecules, potentially turning isotopic stacking into a general provenance tool for interstellar complex molecules.","The paper's own caution about the hidden-line power-law extrapolation (Appendix A: 'no compelling physical reason') implies that the 5σ thresholds carry an unquantified systematic uncertainty; a direct line-count comparison with the higher-frequency QUIJOTE survey at frequencies where cyanopyrene lines cluster would test the hidden-line model."],"forward_implications":["If the forecast holds, a >5σ aggregate detection of 13C-cyanonaphthalene or 13C-cyanopyrene at the local bulk ratio would argue for (though not prove) a circumstellar, top-down inheritance origin, since in-situ formation in TMC-1 would be expected to dilute rather than enrich 13C.","A non-detection would weigh against an enhanced-13C inheritance scenario, because simulations show that at 12C/13C ≤ ~30 a robust detection should be obtained regardless of the noise realization.","The demonstrated recovery of aggregate 13C signal from HC9N shows the technique already works on real observations and can yield preliminary, loose isotope-ratio constraints in advance of detailed MCMC fitting.","The line-confusion stress tests establish that stacking results cannot be artificially inflated by hidden lines, so previously reported stacked detections are not vulnerable to this particular false-positive pathway.","The sensitivity curves can be used to prioritize which isotopologue laboratory spectra to measure: molecules with strong predicted aggregate signal are the ones where laboratory effort is most likely to result in a detection or a meaningful limit."],"fun_headline_variants":["Stacking 13C lines could reveal PAH origins in TMC-1","If lab spectra existed, stacked lines would flag 13C-rich PAHs","Line stacking proves able to unmask 13C enrichment in TMC-1","Synthetic test: stacked 13C lines won't fake a detection"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The forecast depends on the assumption that the rotationally scaled quantum-chemical spectra predict the true 13C isotopologue line frequencies and intensities closely enough that real lines fall inside the matched filter's passband—if the predicted positions or relative strengths are off, the simulated 5σ discrimination thresholds do not transfer to actual observations.","fun_headline_variants_meta":{"raw":{"variants":["Stacking 13C lines could reveal PAH origins in TMC-1","If lab spectra existed, stacked lines would flag 13C-rich PAHs","Line stacking proves able to unmask 13C enrichment in TMC-1","Synthetic test: stacked 13C lines won't fake a detection"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000198,"raw_usage":{"total_tokens":1269,"prompt_tokens":871,"completion_tokens":398,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":615,"completion_tokens_details":{"reasoning_tokens":314}},"tokens_in":615,"tokens_out":398,"duration_ms":4187,"temperature":1.0,"reasoning_tokens":314,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T08:21:49.430685+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the laboratory rotational spectrum of a single singly 13C-substituted isotopologue of 1-cyanonaphthalene, shift it to TMC-1's velocity, and run it through the paper's stacking-and-filtering pipeline in place of the predicted spectrum; if the recovered impulse response drops below the simulated 5σ threshold at 12C/13C ≈ 30, the central forecast fails.","supporting_citations":[],"review_version":1}