{"id":"497525ba-d42a-4ea5-b1d9-1e423b3c9316","arxiv_id":"2504.14001","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"In simulated high-redshift galaxies, the [CII]-to-gas conversion factor varies with local density, metallicity, and starburstiness, so a single global value biases gas surface densities and depletion times.","lead":"Using the SERRA simulations of early galaxies, this paper maps how the [CII]-to-gas conversion factor changes across a galaxy, from dense star-forming centers to faint outskirts. It provides new recipes for turning [CII] maps into gas masses and shows that a single conversion factor can misestimate gas depletion times by up to four times.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Sub-grid density PDF assumption is the load-bearing link for Eq. 10; the paper's power-law-tail caveat is untested, and the Fig. 7 recovery is circular.","rationale":"The reader's weakest assumption (sub-grid density PDF) is also the most load-bearing in my read. The abstract-level statement that a constant alpha_CII flattens the KS relation and overestimates t_dep by up to 4x is a consequence of W_CII decreasing with Sigma_CII (Eq. 10). That decrease is computed from Cloudy models integrated over a log-normal density distribution whose dispersion is tied to the Mach number. [CII] emission is not a linear tracer: it saturates near n_crit ~ 10^3 cm^-3, so the relative weight of densities around and above n_crit in the sub-grid PDF directly sets how much a pixel emits per unit gas mass, which is exactly W_CII. An untested change in that PDF can change the sign or magnitude of the fitted anti-correlation, and with it the correction to the KS slope and t_dep. The paper explicitly flags the power-law tail possibility but dismisses it in one sentence; this is an unsupported assertion, not a robustness test. I also note two reinforcing issues: (i) the Fig. 7 recovery of the KS relation using Eq. 10 is circular because the same pixels were used to fit Eq. 10, so it cannot validate the relation; (ii) the bias is computed at 30 pc resolution, far finer than current ALMA high-z observations, so the observational amplitude is not yet demonstrated. Both are downstream of the PDF question: if the PDF sensitivity test passes, Eq. 10 is anchored in the simulation and these become application/validation caveats; if it fails, the whole bias claim needs re-quantification. Given the missing sensitivity study and the absence of released code/data, the reader's CONDITIONAL verdict is appropriate; I would not upgrade to ACCEPT or move to REJECT because the F19 analytical model and global alpha_CII comparison provide partial support, and the requested check is well-defined and feasible.","tokens_in":21313,"tokens_out":11973,"duration_ms":109038,"concrete_test":"Re-run the Cloudy [CII] post-processing on a representative set of SERRA cells (or 2-3 full galaxies) replacing the log-normal sub-grid PDF with (a) a log-normal plus a power-law tail with slope from self-gravity (e.g., Federrath & Klessen 2013), and (b) a log-normal with dispersion inflated/deflated by 50% at fixed Mach number; refit Eq. 10 on the resulting maps and recompute the KS-slope and t_dep bias of Fig. 7. If A remains below -0.25 and the t_dep overestimate stays above roughly 2x in all variants, the PDF assumption is not load-bearing; if A or the factor changes materially, the central claim must be reported as conditional on the assumed density distribution.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central bias claim follows from the resolved conversion-factor relation (Eq. 10), whose negative slope in Sigma_CII is produced by the Cloudy post-processing. In Sec. 2.2/4.1, each 30 pc cell is assigned a log-normal density PDF with dispersion set by the local Mach number; [CII] is density-sensitive near n_crit ~ 10^3 cm^-3, so the fitted W_CII(Sigma_CII,Z) and hence the KS flattening and ~4x t_dep bias in Fig. 7 depend on the assumed high-density tail of that PDF. The paper states (Sec. 4.1) that a power-law tail has 'minimal' impact because the log-normal samples n_crit for M>10, but no calculation or test is shown; for a typical cell with n~10^2 cm^-3 and M=15, n_crit is roughly 2 sigma above the log-normal median, and [CII] weighting of the tail is not obviously negligible. The Sec. 5.1 'recovery' of the KS relation by applying Eq. 10 is not independent evidence, since Eq. 10 was fitted to the same SERRA pixels; this makes the PDF sensitivity test the key missing check.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper uses the SERRA cosmological zoom-in simulation, which features on-the-fly radiative transfer and ~30 pc spatial resolution, to compute [CII] emission with Cloudy post-processing for 98 galaxies at 4<z<9. The authors derive the global conversion factor alpha_CII = M_gas/L_CII and the spatially resolved W_CII = Sigma_gas/Sigma_CII, finding that alpha_CII decreases with metallicity and that W_CII anticorrelates with [CII] surface brightness and metallicity. Their central quantitative result is the fitted relation in Eq. (10), log W_CII = -0.355 log Sigma_CII - 0.324 log Z + 3.37, and they argue that applying a constant alpha_CII overestimates Sigma_gas in bright [CII] patches, thereby flattening the resolved Kennicutt-Schmidt relation and overestimating the depletion time by up to a factor of about 4. They further provide relations connecting global alpha_CII to galaxy compactness and effective [CII] surface brightness.","tokens_in":21656,"tokens_out":4789,"duration_ms":45994,"significance":"If the central relation holds, the paper provides a practical, physically motivated correction for spatially resolved [CII]-based gas mass measurements in the Epoch of Reionization, with direct relevance for ALMA and JWST observing programs. The analysis builds on a high-resolution simulated sample with on-the-fly radiative transfer, an explicit Cloudy treatment of PDRs, and a physically transparent comparison against the F19 analytic model; the fitting formulae in Table 1 are ready to use. The paper also makes a falsifiable prediction: the slope of the resolved KS relation derived from [CII] should steepen when a spatially varying W_CII is adopted. However, the central quantitative claim currently depends on the assumed sub-grid density PDF and on sign-constrained regression, so the significance of the paper is conditional on the robustness checks requested below.","major_comments":[{"comment":"","section":"Sec. 4.1 and Eq. (10)"},{"comment":"","section":"Sec. 4.2, Eq. (10)"},{"comment":"","section":"Sec. 5.1, Fig. 7"}],"minor_comments":[{"comment":"","section":"Sec. 2.2"},{"comment":"","section":"Sec. 3"},{"comment":"","section":"Table 1"},{"comment":"","section":"Fig. 5"},{"comment":"","section":"Sec. 5.1, Fig. 7"}],"recommendation":"major_revision","confidential_remarks":"This is a solid simulation-based study with a clearly stated central claim and useful fitting relations. The main risk is that the sign of the key resolved relation is imposed by priors and that the sub-grid density PDF, which controls the [CII] weighting, is untested. If the authors can provide the density-PDF robustness test and unconstrained or out-of-sample checks, the paper would be suitable for publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague — this is a solid paper on an important practical problem: converting resolved [CII] maps into gas surface densities at z~4-9. The genuinely new pieces are Eq. 10, a resolved W[CII](Sigma[CII], Z) relation from the SERRA simulation, and the demonstration in Sec. 5.1 that applying a fixed alpha[CII] flattens the Kennicutt-Schmidt slope and overestimates depletion times by up to 4x. That warning is useful for ALMA/JWST users. The paper also gives a clean physical decomposition using the F19 analytic model, and the global alpha-Z trend is consistent with prior work, so the citation context is fair.\n\nWhat is good: the analysis is transparent, fits are reported with scatter, and the figures support the stated trends. The two-regime behavior (low-Z vs high-Z) is physically intuitive, and the compactness relations for global alpha could be practically useful. I found no sign of sloppy reasoning or overclaimed novelty.\n\nThe soft spots, in order of importance. First, Eq. 10 is fit to SERRA pixels and then \"validated\" by recovering the KS relation from the same SERRA pixels. That is a consistency check, not independent evidence, and the paper should say so more plainly. Second, the sub-grid density distribution is a lognormal scaled by Mach number, and [CII] emission is density-sensitive near n_crit ~ 1e3 cm^-3. The stress-test note is right: the claim that a power-law tail has minimal impact is asserted but not tested. For a typical cell, the critical density sits roughly 2 sigma above the lognormal median, so the tail is not obviously negligible. Since no maps or code are released, a reader cannot easily probe this. Third, the sign-constrained priors (A<0, B<0) mean the fitted slopes are not a blind measurement, though the posteriors do look consistent with the plotted scatter. Minor: the Mach number definition is terse.\n\nIs the central claim wrong? Probably not. The direction of the bias is robust in any model where brighter [CII] regions are denser and more metal-rich, and hence have lower W. But the factor-of-4 number and the exact slope of Eq. 10 should be treated as simulation-dependent until an external observational check or a density-PDF sensitivity study is done.\n\nWho is this for: anyone using [CII] to infer gas masses or resolved KS relations at high redshift. I would send it to a serious referee. It deserves publication after a revision that either releases the data or adds a density-PDF sensitivity test, and makes clear that the KS \"recovery\" is not independent validation. My own verdict: conditional accept, not reject.","headline":"Useful, internally consistent simulation-based paper giving a resolved [CII]-to-gas recipe and a clear warning about constant-alpha bias; treat the fitted relation as simulation-dependent until the sub-grid density PDF is stress-tested.","tokens_in":22255,"tokens_out":2452,"would_cite":true,"duration_ms":24332,"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":"A constant [CII]-to-gas conversion factor flattens the resolved star-formation law and inflates depletion times up to 4x in high-redshift galaxies.","keywords":["galaxies: evolution","galaxies: high-redshift","galaxies: ISM","[CII] 158 micron emission","gas conversion factor","Kennicutt-Schmidt relation","depletion time","Epoch of Reionization"],"falsifier":"Compare spatially resolved [CII] maps with independent gas surface density maps (from rest-frame dust continuum, CO, or dynamical modelling) across bright and faint patches of a $z\\approx5-7$ galaxy: if $W_{[\\rm CII]}$ does not decrease with increasing $\\Sigma_{[\\rm CII]}$ at fixed metallicity, or if the Kennicutt-Schmidt slope measured with a constant $\\alpha_{[\\rm CII]}$ is not flatter than the slope from the independent gas tracer, the central claim fails.","tokens_in":21139,"feed_emoji":"🌌","tokens_out":10258,"duration_ms":82092,"temperature":0.7,"pith_summary":"The paper asks whether a single, constant [CII]-to-gas conversion factor can recover gas surface densities in galaxies at redshift 4–9, and answers no. Using a cosmological zoom-in simulation with roughly 30-parsec resolution, it shows that the resolved conversion factor $W_{[\\rm CII]} = \\Sigma_{\\rm gas}/\\Sigma_{[\\rm CII]}$ falls in bright, dense, metal-rich, and bursty patches of the interstellar medium, following the fitted relation $\\log W_{[\\rm CII]} = -0.355\\log\\Sigma_{[\\rm CII]} - 0.324\\log Z + 3.37$. Applying a constant global factor in resolved maps therefore overestimates gas surface density exactly where [CII] is brightest, flattening the star-formation law (the Kennicutt-Schmidt relation) and overestimating the gas depletion time by up to a factor of about four. The paper provides ready-to-use brightness- and metallicity-dependent conversion relations so that ALMA and JWST observations of early galaxies can avoid this bias.","feed_headline":"Fixed [CII]-to-gas factor overestimates depletion time 4x","feed_subtitle":"The true conversion factor falls in bright, metal-rich patches; a metallicity-aware formula recovers the KS slope.","key_machinery":"The object that carries the argument is the spatially resolved conversion factor $W_{[\\rm CII]}\\equiv \\Sigma_{\\rm gas}/\\Sigma_{[\\rm CII]}$, together with the fitted empirical relation $\\log W_{[\\rm CII]} = -0.355\\log\\Sigma_{[\\rm CII]} - 0.324\\log Z + 3.37$ (Eq. 10). The physical machinery behind it is an analytic model for [CII] surface brightness that expresses $\\Sigma_{[\\rm CII]}$ as a function of gas density, metallicity, and burstiness, yielding $W_{[\\rm CII]}\\propto 1/(nZ)$ at $Z<0.2\\,Z_\\odot$ and a $\\kappa_s$-dependent scaling above that threshold. The relation is calibrated on [CII] emission computed with Cloudy PDR grids, interpolated over a sub-grid log-normal density distribution whose dispersion is set by the local Mach number.","core_discovery":"On the paper's own terms, the central discovery is that the [CII]-to-gas conversion factor is not constant across the interstellar medium of early galaxies: $W_{[\\rm CII]}=\\Sigma_{\\rm gas}/\\Sigma_{[\\rm CII]}$ anticorrelates with local gas density, metallicity, and, above $Z\\approx 0.2\\,Z_\\odot$, with the star-formation burstiness $\\kappa_s$, because the same conditions that make a patch bright in [CII] also raise its emissivity per unit gas. The consequence is $W_{[\\rm CII]} \\propto \\Sigma_{[\\rm CII]}^{-0.5}$ over the simulated sample. Replacing the standard constant $\\alpha_{[\\rm CII]}=31\\,M_\\odot/L_\\odot$ with the resolved relation restores the intrinsic ~1.4 Kennicutt-Schmidt slope and lowers the inferred mean depletion time from about 0.4 Gyr to 0.1 Gyr.","pith_inferences":["The same bias likely affects resolved [CII] studies at lower redshift: any galaxy whose bright centre dominates the [CII] emission will show a flattened KS slope when a constant conversion factor is applied, so the prediction can be checked against local samples with independent gas tracers.","A direct observational falsifier is to measure $W_{[\\rm CII]}$ in bright and faint patches of a single $z\\sim5-7$ galaxy using independent gas surface densities from dust continuum or dynamical modelling; the slope of $\\log W_{[\\rm CII]}$ versus $\\log\\Sigma_{[\\rm CII]}$ should be about $-0.36$.","If the relation holds, the inferred cold-gas masses of compact $z>10$ [CII]-undetected galaxies become even lower than a constant-$\\alpha$ analysis suggests, sharpening the claimed tension between their stellar masses and depletion times.","The fitted relation could be ported to other density-sensitive fine-structure lines (e.g. [O III] 88 μm), predicting that their conversion factors also drop in bright, metal-rich patches and similarly bias resolved KS studies."],"forward_implications":["Resolved ALMA observations of [CII] at $z\\approx4-9$ should convert $\\Sigma_{[\\rm CII]}$ to gas surface density using the fitted $W_{[\\rm CII]}(\\Sigma_{[\\rm CII]},Z)$ rather than a constant $\\alpha_{[\\rm CII]}$.","Kennicutt-Schmidt slopes measured from [CII] with a single conversion factor will come out artificially flat; the true, steeper slope is recovered once the brightness-dependence of $W_{[\\rm CII]}$ is included.","Gas depletion times inferred from unresolved [CII] in bursty galaxies can be overestimated by up to roughly $4\\times$, so the census of how quickly high-redshift galaxies consume their gas may need revising.","More compact galaxies, with smaller $r_{e,[\\rm CII]}$ and higher effective [CII] surface brightness, have systematically lower global $\\alpha_{[\\rm CII]}$; unresolved gas masses should therefore be derived with a compactness-dependent conversion factor.","The non-detection of [CII] in compact, metal-poor JWST-selected galaxies at $z>10$ sets tight upper limits on their cold gas content, pointing to very short depletion times."],"supporting_citations":[{"why":"Supplies the analytic [CII] surface-brightness model used to derive the regime-dependent scalings of W[CII] with density, metallicity, and burstiness.","marker":"Ferrara et al. (2019)"},{"why":"Provides the canonical constant α[CII]=31 M_sun/L_sun conversion factor used as the biased baseline in the KS and depletion-time comparisons.","marker":"Zanella et al. (2018)"},{"why":"Defines the sub-grid log-normal density distribution parameterized by Mach number that sets the density interpolation for the Cloudy [CII] models.","marker":"Vallini et al. (2018)"},{"why":"Presents the SERRA zoom-in simulation suite from which the 98-galaxy sample and ~30 pc resolved ISM properties are taken.","marker":"Pallottini et al. (2022)"},{"why":"Provides the Cloudy photoionization/PDR grid used to compute [CII] emission cell by cell.","marker":"Ferland et al. (2017)"},{"why":"Sets the Hα-based SFR calibration and the reference Kennicutt-Schmidt relation against which the bias is quantified.","marker":"Kennicutt (1998)"}],"fun_headline_variants":["Variable [CII] factor fixes early galaxy gas estimates","Metallicity-aware [CII] factor cuts depletion time 4x","Bright [CII] patches mislead gas mass if factor fixed","Variable [CII] conversion restores Kennicutt-Schmidt slope"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the unresolved gas inside each 30-parsec cell has a log-normal density distribution whose width is set by the local Mach number; because [CII] emission is most sensitive near gas densities of about a thousand particles per cubic centimetre, a different small-scale density structure would change the fitted $W_{[\\rm CII]}$ relation and the size of the predicted biases.","fun_headline_variants_meta":{"raw":{"variants":["Variable [CII] factor fixes early galaxy gas estimates","Metallicity-aware [CII] factor cuts depletion time 4x","Bright [CII] patches mislead gas mass if factor fixed","Variable [CII] conversion restores Kennicutt-Schmidt slope"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000568,"raw_usage":{"total_tokens":2799,"prompt_tokens":1161,"completion_tokens":1638,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":777,"completion_tokens_details":{"reasoning_tokens":1562}},"tokens_in":777,"tokens_out":1638,"duration_ms":10653,"temperature":1.0,"reasoning_tokens":1562,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:58:31.610535+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compare spatially resolved [CII] maps with independent gas surface density maps (from rest-frame dust continuum, CO, or dynamical modelling) across bright and faint patches of a $z\\approx5-7$ galaxy: if $W_{[\\rm CII]}$ does not decrease with increasing $\\Sigma_{[\\rm CII]}$ at fixed metallicity, or if the Kennicutt-Schmidt slope measured with a constant $\\alpha_{[\\rm CII]}$ is not flatter than the slope from the independent gas tracer, the central claim fails.","supporting_citations":[],"review_version":1}