{"id":"4cf450b3-0d79-456a-ad5d-0467fc955dee","arxiv_id":"2412.10506","paper_version":3,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"Across 31 nearby galaxies, HCN/CO rises and SFR/HCN falls with increasing stellar mass surface density, molecular gas surface density, and pressure, implying the Gao-Solomon scatter is physically driven.","lead":"This paper maps how much dense gas galaxies contain and how efficiently that gas forms stars, using the largest sample of resolved galaxy observations to date. It shows that dense gas fractions rise and star formation efficiencies fall in high-pressure, high-density galactic environments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'demonstrate a physical origin' claim is undercut by §3.3's own concession that α_HCN and HCN excitation may vary with environment; line-ratio trends alone cannot yet distinguish physical SFE variations from tracer artifacts.","rationale":"The reader's conditional verdict is appropriate, but the most load-bearing concern is not the CO-to-H2 conversion factor used for the Σmol and PDE axes. That factor affects the quantitative slopes and would shift them if wrong, yet the qualitative trends (and the independent Σ⋆ axis) would likely survive. The deeper problem is the 'demonstrate physical origin' claim, which requires excluding the possibility that HCN excitation or α_HCN variations produce the SFR/HCN trends. The paper itself flags this in §3.3, saying that a similarly likely explanation is that HCN is not a robust tracer of dense gas in centres, and that a lower α_HCN in centres would make inferred efficiencies comparable to discs. Without a formal test (e.g., alternative dense gas tracers or an HCN excitation diagnostic), the abstract's 'demonstrate' is not supported; 'reinforce' or 'are consistent with' would be more accurate. This agrees with the reader's rationale but not with the reader's nominated weakest assumption, hence 'partial'. The empirical dataset, public code, and homogeneous SFR calibration are strong, and the trends are likely real as line-ratio trends; the issue is specifically the interpretive claim. A conditional acceptance requiring language softening or an additional test would be reasonable, so the verdict remains unchanged from the reader's CONDITIONAL.","tokens_in":22833,"tokens_out":5391,"duration_ms":53088,"concrete_test":"Re-run the same spectral-stacking analysis measuring SFR/HCO+ and SFR/CS (or SFR/HCN after applying a physically motivated α_HCN(Σ⋆) correction) in the same bins of Σ⋆, Σmol, and PDE. If the anti-correlation with environment persists for HCO+ and CS, the physical-origin interpretation is supported; if it weakens or vanishes, the HCN-specific excitation/conversion artifact is the more likely explanation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract and §3.2 interpret the decreasing SFR/HCN with increasing Σ⋆, Σmol, and PDE as demonstrating that the scatter in the Gao–Solomon relation is physically caused. But §3.3 explicitly acknowledges that HCN may not robustly trace dense gas in galaxy centres because of increased optical depth, IR pumping, or electron excitation, and that α_HCN could plausibly be lower in centres, which would raise the inferred SFR/M_dense and make centre values comparable to discs. If α_HCN varies with environment, the observed SFR/HCN anti-correlation could be produced without any change in true dense-gas star formation efficiency. The paper presents only line ratios and fixed conversion factors; it includes no HCN excitation diagnostic (e.g., HCN(3-2)/(1-0)) or independent dense-gas mass calibration to break this degeneracy. Moreover, no formal statistical test shows that the environmental trends actually reduce or explain the 0.4 dex scatter in the Gao–Solomon relation; a significant correlation with an environmental variable does not, by itself, demonstrate that the scatter is physical in origin. Since the paper's own limitations section raises this alternative, the central 'demonstrate' wording goes beyond what the data establish.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This Letter presents resolved kiloparsec-scale measurements of HCN(1-0)/CO(1-0) and SFR/HCN(1-0) as functions of stellar mass surface density, molecular gas surface density, and dynamical equilibrium pressure, combining the ALMA ALMOND survey with the IRAM 30m EMPIRE survey for a total of 31 nearby galaxies. The authors spectrally stack HCN and CO in bins of each environmental quantity, fit linear relations with LinMix while treating non-detections as limits, and report that HCN/CO rises (slope ~0.5) and SFR/HCN falls (slope ~ -0.6) with increasing surface density and pressure. They also place ALMOND/EMPIRE in the context of a literature compilation of the Gao-Solomon relation. The central claim is that the environmental trends demonstrate that the scatter in the Gao-Solomon relation has a physical origin rather than being an artifact of the tracer.","tokens_in":23044,"tokens_out":4189,"duration_ms":40763,"significance":"The paper is a valuable empirical contribution: it more than triples the number of galaxies with resolved HCN maps, homogenizes SFR, CO, and HCN calibrations across two surveys, and provides public data products and analysis scripts. The stacking and limit-aware fitting procedures are standard and clearly described, and the cross-survey consistency checks in Appendix A strengthen confidence in the measurements. If the trends are robust, they constrain models of how the dense gas fraction and dense-gas star formation efficiency depend on galactic environment. However, the headline interpretation goes beyond what the data alone can establish, given the paper's own admission in §3.3 that HCN excitation and the HCN-to-dense-gas conversion factor may vary with environment.","major_comments":[{"comment":"The abstract states that the results 'demonstrate that the scatter in the Gao–Solomon relation (SFR/HCN) has a physical origin,' and §4 says this 'reinforces the notion' of physical origin. However, §3.3 explicitly concludes that HCN may not be a robust dense-gas tracer in galaxy centres because of increased optical depth, IR pumping, electron excitation, and a plausibly varying alpha_HCN, which would raise SFE_dense in centres and make them comparable to discs. The paper presents only line ratios and fixed conversion factors; it includes no HCN excitation diagnostic (e.g., HCN(3-2)/(1-0)) and no test of whether the environmental trends actually reduce the 0.4 dex scatter of the Gao–Solomon relation. A significant correlation with an environmental variable does not by itself demonstrate that the scatter is physical in origin, because the alternative of a tracer-ratio artifact is not excluded. I recommend softening the abstract and conclusions to 'consistent with a physical origin' and, ideally, adding a quantitative comparison of the scatter before and after removing the environmental trend.","section":"Abstract, §3.3, §4"},{"comment":"The relations of HCN/CO and SFR/HCN against Sigma_mol and P_DE share W_CO (or quantities derived from it, via Eq. (2) and Eq. (1)) on both axes, so correlated noise and calibration errors in CO can bias the fitted slopes and produce artificial correlations. The paper acknowledges this issue only in §3.3, where it chooses Sigma_star for the centre–disc comparison because it 'have[s] uncorrelated axes,' yet the abstract and §3.2 present all three x-axis variables as equivalent measurements of the same phenomenon. The Sigma_star relation avoids the covariance, but the quantitative slopes for Sigma_mol and P_DE should either be accompanied by an estimate of the covariance bias (e.g., a Monte Carlo error analysis that propagates CO uncertainties on both axes) or be de-emphasized in the headline claims.","section":"§3.2, Eq. (2), Table 2"},{"comment":"The Pearson correlation coefficients and the statement that 'all p-values are much less than 0.01' are computed from stacked bins that are not independent: each galaxy contributes multiple bins, and neighbouring bins within a galaxy share the underlying physical conditions and measurement systematics. The effective sample size is therefore substantially smaller than the raw number of stacked points, so the quoted p-values overstate the significance of the correlations. A per-galaxy bootstrap, a hierarchical regression, or a randomization test that preserves the within-galaxy structure is needed to support the significance claims.","section":"Table 2"}],"minor_comments":[{"comment":"The name 'Frank Bigiel' appears twice in the author list; one occurrence appears to be a duplication.","section":"Author list"},{"comment":"The abstract quotes a scatter of ~0.4 dex for SFR/HCN, while Table 2 reports scatter values sigma = 0.28–0.34 dex for the three relations; please harmonize these numbers or clarify what the 0.4 dex refers to.","section":"Abstract vs. Table 2"},{"comment":"The variable alpha_CO prescription is adopted without an explicit treatment of its uncertainty; given that this factor enters the Sigma_mol and P_DE axes, a brief statement on the adopted uncertainty or at least a reference to its origin would help readers gauge the systematic error in the slopes.","section":"Appendix C, Eq. (C.2)"},{"comment":"The caption's phrase 'median (all S/N)' is ambiguous; the text in §3.2 clarifies that this includes non-detections, so please make the caption explicit that the cyan line includes S/N < 3 data.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The empirical content of the paper is solid and the data products are a valuable community resource. The main issue is that the abstract and conclusions overclaim what can be concluded from line ratios alone, and the paper itself contains the limitation needed to see the overclaim. I recommend major revision to recalibrate the language and add the missing statistical support, but I do not see a fatal flaw that would require rejection."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"This is the field's largest resolved HCN(1-0) scaling-relation sample: 31 galaxies with a homogeneous SFR calibration, stacking in environmental bins, and fits that include upper/lower limits. The paper deserves credit for the careful homogenization of ALMOND and EMPIRE and for making data and scripts public. The HCN/CO versus Sigma_star relation is the cleanest result because it avoids sharing CO on both axes; that trend is likely robust.\n\nThe main soft spot is the abstract's claim that the results 'demonstrate' a physical origin for the scatter in the Gao-Solomon relation. The paper's own §3.3 concedes that alpha_HCN and HCN excitation may vary with environment, which could lower SFR/HCN in centers without any change in true dense-gas SFE. The paper presents no HCN excitation diagnostic or independent dense-gas mass calibration to break that degeneracy, and no formal test that the environmental trends actually reduce the 0.4 dex scatter. The discussion correctly uses 'reinforces,' but the abstract and §3.2 go further than the data allow. The center-disc SFR/HCN comparison also overlaps at roughly 1-sigma, so 'factor of a few lower' in the abstract should carry that uncertainty.\n\nThe Sigma_mol and PDE relations depend on a variable alpha_CO prescription and share CO on both axes, so those quantitative slopes are more model-dependent; the paper flags this. The Sigma_star result stands on its own. I do not see a load-bearing flaw: the measurements, stacking, and fits are described clearly, and the limitations are honestly acknowledged.\n\nFor a reviewer, I would require softening 'demonstrate' to 'support' or 'are consistent with,' and ask for a sentence acknowledging that alpha_HCN variations could explain part of the trend without changing true SFE. A test of whether including the environment actually reduces residual scatter would be a plus, but not a blocker. This paper deserves a regular referee, and I'd bring it to the reading group.","headline":"Largest resolved HCN scaling-relation sample to date, but the 'physical origin' claim outruns the line-ratio evidence.","tokens_in":23732,"tokens_out":2997,"would_cite":true,"duration_ms":27065,"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":"This paper establishes that HCN/CO increases and SFR/HCN decreases with local surface density and pressure across 31 galaxies, giving the Gao-Solomon scatter a physical origin.","keywords":["dense gas tracers","HCN(1-0)","Gao-Solomon relation","star formation efficiency","molecular gas surface density","dynamical equilibrium pressure","spectral stacking","nearby galaxies"],"falsifier":"A dust-based molecular gas surface-density map of the same 31 galaxies, independent of CO, that disagrees systematically with the adopted variable conversion factor would recalibrate the molecular gas and pressure axes; if the HCN/CO and SFR/HCN trends disappeared under that recalibration, the claimed environmental relations for those two axes would be refuted, while the stellar-mass-based relation would survive.","tokens_in":22581,"feed_emoji":"🌌","tokens_out":6152,"duration_ms":49361,"temperature":0.7,"pith_summary":"This paper argues that the ratio of dense gas (HCN) to bulk molecular gas (CO) and the star formation rate per unit HCN both vary systematically with the local galactic environment, not randomly. Using the ALMA ALMOND and IRAM 30m EMPIRE surveys, it stacks HCN(1-0) and CO(1-0) spectra across 31 nearby spiral galaxies and measures how these ratios change with molecular gas surface density, stellar mass surface density, and dynamical equilibrium pressure. The result is that HCN/CO rises with a slope of about 0.5 and SFR/HCN falls with a slope of about -0.6 as these environmental quantities increase, with scatter of roughly 0.2 and 0.4 dex. A sympathetic reading of the paper takes this as evidence that the scatter in the Gao-Solomon relation, SFR against HCN luminosity, has a physical origin in the varying state of molecular gas. This matters because it links observable line ratios to the conditions that regulate star formation across galaxy discs.","feed_headline":"HCN/CO rises, SFR/HCN falls with surface density in 31 galaxies","feed_subtitle":"The Gao-Solomon scatter traces real environment-driven changes in dense gas fraction and star formation efficiency.","key_machinery":"Spectral stacking of HCN(1-0) and CO(1-0) in bins of environmental conditions. For each galaxy, spectra are co-added in logarithmic bins of stellar mass surface density, molecular gas surface density, and dynamical equilibrium pressure, using the high-signal-to-noise CO line to set the velocity window; non-detections are kept as upper and lower limits and the relations are fitted with a linear regression that incorporates those limits. A variable CO-to-H2 conversion factor that depends on metallicity and stellar mass surface density calibrates the molecular gas surface density and pressure axes.","core_discovery":"The central discovery is that the two ratios that connect dense gas to star formation are not universal constants but continuous functions of kiloparsec-scale environment. Across 31 local galaxies, HCN/CO increases and SFR/HCN decreases with increasing stellar mass surface density, molecular gas surface density, and dynamical equilibrium pressure, with galaxy centres lying at the high-density, high-pressure end of the same trend. The slopes are significant, about 0.5 dex per dex for HCN/CO and about -0.6 for SFR/HCN, and the residual scatter is modest for HCN/CO (0.2 dex) and larger for SFR/HCN (0.4 dex). The authors interpret this as showing that deeper gravitational potentials and more abundant gas produce denser molecular clouds, that dense gas in these environments is less efficient at forming stars per unit mass, and that the well-known scatter in the Gao-Solomon relation therefore reflects real environmental physics rather than measurement noise.","pith_inferences":["If the environmental trends hold at cloud scale, resolved HCN observations at sub-kiloparsec resolution should find the same HCN/CO and SFR/HCN gradients within single galaxies, offering a direct test of the claimed continuity.","A dust-based molecular gas map of the same galaxies, independent of CO, would test the variable conversion-factor calibration; the stellar-mass relation would remain robust because it does not use that calibration.","The lower SFR/HCN in centres could also reflect HCN tracing more bulk gas there through optical depth or IR pumping; if so, a different dense gas tracer such as HCO+ might show a weaker anti-correlation, distinguishing environmental efficiency from tracer bias.","Applying the same stacking method to the HCO+ and CS data already observed in ALMOND could test whether the trends are specific to HCN or general to dense gas tracers."],"forward_implications":["The Gao-Solomon relation's ~0.5 dex scatter is not random; it correlates with measurable local conditions, so SFR/HCN predictions can be improved by including surface density or pressure.","Galaxy centres naturally host a higher fraction of HCN-bright gas than discs, following the same continuous environmental trend without needing a distinct centre mechanism.","SFR/HCN decreases in high-surface-density, high-pressure regions, implying the efficiency of star formation per unit dense gas varies by a factor of a few between centres and discs.","The similar trends seen with cloud-scale CO imaging and with kiloparsec-scale HCN/CO spectroscopy support a picture where molecular cloud properties vary with galactic environment.","Because ALMOND and EMPIRE agree on overlapping galaxies, the combined 31-galaxy sample provides a homogeneous benchmark for future resolved dense gas surveys."],"supporting_citations":[{"why":"Supplies the EMPIRE HCN and CO maps and the previous nine-galaxy scaling relations that this work extends.","marker":"Jiménez-Donaire et al. (2019)"},{"why":"Provides the ALMOND HCN maps and the original ALMOND methodology for SFR estimation and stacking.","marker":"Neumann et al. (2023b)"},{"why":"Supplies the ALMA CO(2-1) maps used for ALMOND galaxies to define velocity windows and compute CO intensities.","marker":"Leroy et al. (2021b)"},{"why":"Defines the Gao-Solomon SFR-HCN relation whose scatter is the paper's target of explanation.","marker":"Gao & Solomon (2004)"},{"why":"Previous resolved study showing environmental trends in HCN/CO and SFR/HCN in M51, providing the baseline these relations extend.","marker":"Usero et al. (2015)"},{"why":"Provides the variable CO-to-H2 conversion factor used to calibrate molecular gas surface density and pressure axes.","marker":"Schinnerer & Leroy (2024)"},{"why":"Supplies the gas velocity dispersion and the dynamical equilibrium pressure estimation approach.","marker":"Sun et al. (2018)"},{"why":"Provides the CO(2-1)/CO(1-0) line ratio calibration used to convert ALMOND CO intensities to CO(1-0).","marker":"den Brok et al. (2021)"}],"fun_headline_variants":["Dense gas ratios shift with environment across 31 galaxies","Star formation efficiency drops as galaxy density rises","Gao-Solomon scatter is physical, not noise, in 31 galaxies","HCN/CO slope 0.5, SFR/HCN slope -0.6: environment rules","Galactic centers follow same dense gas trend as disks"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The adopted variable CO-to-H2 conversion factor, which sets how molecular gas surface density and pressure are inferred from CO emission, correctly captures how those quantities vary with metallicity and stellar mass across the 31 galaxies.","fun_headline_variants_meta":{"raw":{"variants":["Dense gas ratios shift with environment across 31 galaxies","Star formation efficiency drops as galaxy density rises","Gao-Solomon scatter is physical, not noise, in 31 galaxies","HCN/CO slope 0.5, SFR/HCN slope -0.6: environment rules","Galactic centers follow same dense gas trend as disks"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00026,"raw_usage":{"total_tokens":1645,"prompt_tokens":1054,"completion_tokens":591,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":670,"completion_tokens_details":{"reasoning_tokens":498}},"tokens_in":670,"tokens_out":591,"duration_ms":5506,"temperature":1.0,"reasoning_tokens":498,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T15:53:28.383932+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A dust-based molecular gas surface-density map of the same 31 galaxies, independent of CO, that disagrees systematically with the adopted variable conversion factor would recalibrate the molecular gas and pressure axes; if the HCN/CO and SFR/HCN trends disappeared under that recalibration, the claimed environmental relations for those two axes would be refuted, while the stellar-mass-based relation would survive.","supporting_citations":[{"cited_title":"J., Bigiel , F., Leroy , A","cited_arxiv_id":null,"evidence_quote":"Supplies the EMPIRE HCN and CO maps and the previous nine-galaxy scaling relations that this work extends."},{"cited_title":"K., Walter , F., et al","cited_arxiv_id":null,"evidence_quote":"Previous resolved study showing environmental trends in HCN/CO and SFR/HCN in M51, providing the baseline these relations extend."},{"cited_title":"New Constraints on the $^{12}$CO(2-1)/(1-0) Line Ratio Across Nearby Disc Galaxies","cited_arxiv_id":"2103.10442","evidence_quote":"Provides the CO(2-1)/CO(1-0) line ratio calibration used to convert ALMOND CO intensities to CO(1-0)."}],"review_version":1}