{"id":"16e3ef5b-7433-4c56-9e32-a8bdffe2ec77","arxiv_id":"1908.05504","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A new grid of photoionization models calibrates UV, optical, and IR line ratios to ISM pressure and electron density, revealing strong metallicity and ionization-zone effects on traditional density diagnostics.","lead":"This paper uses new photoionization model grids to turn emission-line ratios into measurements of interstellar medium pressure and electron density in star-forming galaxies. It shows that the standard [S II] and [O II] density indicators depend on metallicity and trace different zones, so using them interchangeably can bias results.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Abundance scaling and unpublished Fe-depletion anchor are the load-bearing assumption; no sensitivity analysis is given for high-z grids.","rationale":"The reader's weakest assumption—the Nicholls et al. (2017) abundance scaling and fixed Fe depletion—is exactly the load-bearing external assumption. The paper's internal logic is self-consistent: given that abundance pattern, MAPPINGS v5.1 produces the grids, and the grids do show strong metallicity dependence and distinct ionization zones. The concern is not an internal inconsistency; it is that the advertised applicability to high-redshift galaxies assumes a local abundance pattern and an unpublished depletion anchor. The paper itself flags the sensitivity to metallicity calibration in Section 5.4, so the authors are aware that the choice matters. I agree with the CONDITIONAL verdict: the central claim is plausible and useful within the modeling framework, but a quantitative sensitivity test to the abundance scaling and depletion would settle whether the stated 0.2–0.4 dex systematic effects are robust. The placeholder zero values in Tables 1–2 are a secondary reproducibility issue and should be replaced with non-detection symbols, but they do not change the scientific conclusion. No change to the reader's verdict is needed.","tokens_in":96177,"tokens_out":6048,"duration_ms":68658,"concrete_test":"Recompute the [S II] λ6717/λ6731 and [O II] λ3729/λ3726 columns of Table 1 at log(P/k)=5.0, log(q)=7.5, Z=8.23 and 8.93 with MAPPINGS v5.1, changing only the gas-phase Fe depletion from -1.5 dex to -1.0 and -2.0 dex, and in a second arm using an α-enhanced abundance pattern ([α/Fe]=+0.3 at fixed O/H). If the resulting line ratios change by more than the 0.02–0.05 dex level implied by the figures, the published grids need abundance-pattern sensitivity bounds; if the change exceeds 0.4 dex in derived ne, the central caution about [S II]/[O II] interchangeability must be re-quantified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Every calibration curve in Tables 1–2 and Figures 5–17 inherits the nebular abundance pattern of Nicholls et al. (2017): the local B-star abundance scale, linear α-element scaling, and a single gas-phase Fe depletion of -1.5 dex quoted from an unpublished analysis ('Dopita et al., in prep', Section 2.2). The central claim is that these grids hold from local to high-redshift galaxies over log(O/H)+12 = 7.63–9.23. But the line ratios that define the grids are set by the electron temperature structure, which is set by the coolants; if the real abundance pattern at fixed O/H (S/O, N/O, C/O, α-enhancement) or the dust depletion differs from the assumed local pattern, the same observed ratio maps to a different P or ne. The paper's own Section 5.4 states that different metallicity calibrations disagree by up to ~1 dex, and Section 4.7.5 reports a 0.4 dex metallicity effect on [S II]-derived density; Sections 4.3 and 4.8 report ~1 dex pressure effects. No grid variants, sensitivity runs, or propagated uncertainties are provided, and the Fe-depletion anchor is an unpublished reference. The internal model sequence is coherent, but the external claim—that the diagnostics are valid at high redshift—is only as strong as the assumed abundance scaling.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents theoretical grids of photoionization models (MAPPINGS v5.1 with Starburst99 ionizing spectra) that calibrate emission-line ratios as diagnostics of ISM pressure (P/k) and electron density (ne) in H II regions and star-forming galaxies. Three model families are computed: constant-pressure models (log P/k = 4 to 9), constant-density models (log ne = 0 to 5), and single-atom fixed-temperature models, each spanning metallicities log(O/H)+12 = 7.63 to 9.23 and ionization parameters log q = 6.5 to 8.5. The calibrations cover UV ratios ([Si III]/Si III], [C III]/C III], [Al II]/Al II], [Ne IV], [N IV]/N IV]), optical ratios ([O II], [S II], [Ar IV], [N I], [Cl III]), and infrared fine-structure ratios ([S III], [O III], [N II]), with tabulated flux ratios in Tables 1 and 2 and calibration curves in Figures 5-17. The main conclusions are that the diagnostics depend strongly on gas-phase metallicity (up to ~1 dex in log P/k) with a weaker ionization-parameter dependence; that different ratios probe different ionization zones, with the [S II] and [O II] doublets tracing different regions so that they should not be used interchangeably unless the electron temperature is known to be constant; and that ISM pressure is a more robust quantity to measure than a single electron density in nebulae with gradients or clumps.","tokens_in":96454,"tokens_out":20872,"duration_ms":190016,"significance":"If the calibrations stand, this is a timely community resource: a single, internally consistent set of pressure and density grids from the UV to the far-IR for JWST- and ALMA-era spectroscopy. The construction is genuinely forward-modeling — the grids are computed from photoionization theory, not fitted to the pressures or densities they calibrate — so no circularity arises. Notable strengths include the [S II] calibration agreeing with the classical Osterbrock (1989) relation near Z ~ 8.9-9.2 (Section 4.7.5); the plane-parallel diagnostics recovering the density structure of a spherical model to within a few percent (Section 5.1, Figure 19); and the zone-separation claim directly supported by the model emission profiles in Figure 21, making the [S II]/[O II] warning a falsifiable prediction. The quantified metallicity sensitivity — 0.4 dex on [S II]-derived densities, ~1 dex on pressures — is itself an important caution for the high-redshift literature.","major_comments":[{"comment":"The calibration grids are computed with a single nebular abundance pattern — the local B-star scale of Nieva & Przybilla (2012) with the scaling of Nicholls et al. (2017) — and a single gas-phase Fe depletion of -1.5 dex anchored to an unpublished analysis cited as 'Dopita et al., in prep'. This choice is load-bearing for the central claim that Tables 1-2 and Figures 5-17 map line ratios to P and ne from local to high-redshift galaxies over Z = 7.63-9.23. The paper itself quantifies metallicity effects of up to ~1 dex in log(P/k) for the pressure diagnostics (Sections 4.3 and 4.8.2) and 0.4 dex for the [S II] density diagnostic (Section 4.7.5), and Section 5.4 concedes that competing metallicity calibrations disagree by up to 1 dex; an error in the abundance pattern or depletion at fixed O/H therefore propagates almost directly into the derived P or ne, with the Si- and Al-based UV diagnostics (Sections 4.3-4.5) affected directly because those elements are depleted. No grid variants, sensitivity runs, or propagated depletion uncertainties are presented. I request bracketing models spanning, for example, Fe depletions of -1.0 and -2.0 dex and at least one alternative alpha-element or S/O scaling, with the resulting shifts in the calibration curves quantified, so that the high-redshift and super-solar validity claims are supported rather than assumed.","section":"Section 2.2: abundance scaling and Fe depletion anchor."},{"comment":"No uncertainty budget is given for the calibrations. Tables 1 and 2 list flux ratios to five significant figures and Figures 5-17 plot single calibration curves, yet Section 3 shows that the Case-2 ratios (e.g., C III], Si III]) are 'extremely sensitive' to electron temperature within their usable density ranges, that the IR Case-3 ratios have strong temperature dependence in the low-density limit, and that Case-1 ratios shift by up to ~0.15 dex between log Te = 3.5 and 4.5; Section 5.3 acknowledges that differing atomic datasets can shift derived densities (Copetti & Writzl 2002; Proxauf et al. 2014); and Section 4.7.1 attributes the shape difference between its [O II] calibration and Osterbrock (1989) to changes in oxygen atomic data. A user therefore cannot tell whether a disagreement between two diagnostics is a physical zone effect, as argued in Section 5.2, or an atomic-data offset. I request that the paper propagate collision-strength and A-value uncertainties on at least the [O II], [S II], [S III], [O III], and [N II] calibrations, or state an explicit accuracy floor for the tables.","section":"Tables 1-2 and Section 5.3: uncertainty budget."}],"minor_comments":[{"comment":"Equation (3) gives the radiative depopulation rate as R_ji^rad = ne Nj Aji, but spontaneous radiative decay does not depend on the electron density and the expression is dimensionally inconsistent (cm^-6 s^-1 on the right-hand side). It should read R_ji^rad = Nj Aji; Equation (4) already contains the required ne factor for collisional de-excitation. The subsequent critical-density definition is the standard one, so the model grids are unaffected, but the equation as printed is incorrect.","section":"Section 3.1, Eq. (3)."},{"comment":"The density-recovery percentages are reported inconsistently: the text states that '[S II] and [N II]' match the true density to within 3% and that '[O II] and [N II]' agree to within 5% and 7%, while the Figure 19 caption says '[O II] and [N II]' provide the best fit, within 3%. One of the [N II] mentions is likely a typo for [N I]; please reconcile the text and caption so that each of [S II], [O II], [N II], and [N I] is assigned its measured recovery accuracy.","section":"Section 5.1 and Figure 19 caption."},{"comment":"Unit labels are garbled in the electron-density ranges: Section 4.8.1 gives '2.5 < log(ne cm3) < 4' (presumably log(ne/cm^-3)) and Section 4.8.3 quotes the [N II] density range as '1 < log(P/k) < 2.5' and 'log(P/k) < 1', where the quantity must be log(ne), not log(P/k).","section":"Sections 4.8.1 and 4.8.3."},{"comment":"In Table 1 the [O III] column header lists 'lambda52um lambda88um', whereas the text and Figure 16 use the ratio [O III] 88um/[O III] 52um; the tabulated values (greater than unity at low pressure) are consistent with 88/52, so the header ordering appears inconsistent. Please state the ratio convention explicitly and apply it uniformly to the table columns and figures.","section":"Table 1, [O III] column."},{"comment":"The tables contain many entries '...' and exact '0.0000' values, but the meaning is not documented. If '...' flags lines below the stated detectability threshold (flux < 1e-5 Hbeta) and '0.0000' flags a zero or undefined ratio, a footnote is needed; as printed, a reader cannot tell whether a quoted 0.0000 is an undetected numerator line, a saturated ratio, or a missing model.","section":"Tables 1 and 2."},{"comment":"Several notation slips should be cleaned up: in Section 1 the pressure-density relation appears as 'n = P Tek' and should be n = P/(Te k); in Section 2.2 the units of P/k are given as 'cm^-3k' and should be K cm^-3; and the same section should specify that the relative metallicity scaling is applied to the nebular abundances, not the stellar tracks, to avoid ambiguity.","section":"Sections 1 and 2.2, notation."},{"comment":"Several load-bearing inputs are cited as unpublished: the Fe-depletion anchor ('Dopita et al., in prep', Section 2.2), the matching metallicity surface ('Kewley et al., ARAA, submitted', Section 5.4), and the temperature diagnostics ('Nicholls et al., in prep', Section 3.2). For a calibration paper whose tables will be used at face value, at least the depletion anchor should be backed by a citable public analysis by the time of publication.","section":"Reference list."}],"recommendation":"major_revision","confidential_remarks":"The paper's main data product (Tables 1 and 2) will be used as published calibrations, and at least one load-bearing input — the -1.5 dex Fe depletion in Section 2.2 — is cited to an in-preparation paper by the same author group. I would ask the editor to require that this anchor be publicly documented or bracketed by sensitivity runs before acceptance. The manuscript otherwise fits the journal's scope well, and both major comments are addressable with additional calculations rather than conceptual reworking."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"David,\n\nThis is a method paper from Kewley et al. that builds a unified grid of UV-optical-IR pressure and density diagnostics from MAPPINGS v5.1 + Starburst99. The genuinely new piece is the breadth and internal consistency: one grid, five metallicities, nine ionization parameters, and pressure and density tracks, all tabulated for a dozen line ratios. The old individual line diagnostics (Si III, C III, etc.) are properly cited, and the paper adds a coherent treatment of how each ratio depends on temperature, pressure, and metallicity. The claim that [S II] and [O II] are produced in different zones - and should not be used interchangeably - is well supported by the model emission profiles in Figure 21. That is a caution worth taking seriously.\n\nThe internal model logic is sound. The [S II] relation broadly matches Osterbrock (1989) at the right metallicities. The paper is honest about the known contamination sources - DIG, shocks, atomic-data limitations - and it explicitly warns, in Section 5.4, that metallicity calibrations can disagree by up to 1 dex and that the diagnostics should only be used with consistent metallicity calibrations.\n\nWhere it gets softer: the high-redshift reach is the load-bearing extrapolation. The abundance patterns come from the local B-star scale of Nieva & Przybilla (2012), with a linear alpha-element scaling and a single Fe depletion of -1.5 dex that is anchored to an unpublished analysis (Dopita et al., in prep). No sensitivity runs show what happens if S/O, N/O, C/O or the depletion differ at fixed O/H, and the paper itself reports 0.4-1 dex pressure/density shifts from metallicity. The warning to use 'consistent' calibrations is reasonable but not a substitute for a quantitative error budget. A referee should push for at least one or two grid variants on the abundance scaling and for a statement of atomic-data uncertainty propagation, or at least a clear acknowledgment that the quoted calibrations are not uncertainty-quantified.\n\nOne practical issue: the preprint tables contain placeholder zero entries in the machine-readable data that could seriously mislead anyone interpolating. Those need to be flagged as non-detections or blanked before release.\n\nBottom line: this is a useful paper for the JWST/ELT planning era. It deserves peer review. I would accept it if the authors add a sensitivity analysis for abundance scaling, fix the table placeholders, and are more careful about separating what is robust (the relative behavior of the ratios) from what is assumed (the absolute abundance pattern at high-z). I'd be comfortable citing the grids once those fixes are in.\n\n- [Your name]","headline":"A solid, internally consistent calibration-grid paper that is worth refereeing; the high-redshift claims rest on an unquantified abundance-scaling assumption, but the inner logic holds.","tokens_in":96998,"tokens_out":3370,"would_cite":true,"duration_ms":34546,"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 derives self-consistent theoretical grids that turn UV, optical, and infrared emission-line ratios into ISM pressure and electron density for star-forming galaxies, and demonstrates that the [S II] and [O II] density…","keywords":["galaxies: starburst","galaxies: abundances","galaxies: fundamental parameters","ISM pressure","electron density","emission-line diagnostics","photoionization models","H II regions"],"falsifier":"Measure [S II], [O II], and the electron temperature in a set of spatially resolved H II regions spanning metallicities $\\log(\\mathrm{O/H})+12 = 7.6$ to $9.2$. The claim predicts that at fixed electron temperature the two ratios give the same density and shift with metallicity by up to 0.4 dex; observing that they disagree at constant temperature by more than that would falsify the zoning and metallicity calibration.","tokens_in":1688,"feed_emoji":"🔭","tokens_out":3887,"duration_ms":123043,"temperature":0.7,"pith_summary":"The paper is trying to establish a theoretically grounded way to read the physical state of ionized gas in star-forming galaxies from its emission lines, usable from local H II regions to high redshift. It produces grids that convert UV, optical, and IR line ratios into ISM pressure and electron density, and argues that pressure is often a more meaningful quantity than a single density because real nebulae contain temperature and density structure. It also demonstrates a specific hazard: the widely used [S II] and [O II] density diagnostics trace different ionization zones, so treating them interchangeably can mislead unless the electron temperature is constant.","feed_headline":"New grids map emission lines to ISM pressure and density","feed_subtitle":"Calibrations cover the full pressure range of H II regions and show why [S II] and [O II] density probes disagree.","key_machinery":"The load-bearing machinery is a grid of constant-pressure photoionization models in which electron temperature and density are solved step by step through the nebula, together with companion constant-density grids and simple fixed-temperature atomic models that isolate temperature effects. Each diagnostic line ratio considered—[Si III], [C III], [Al II], [Ne IV], [N IV], [O II], [Ar IV], [N I], [Cl III], [S II], [S III], [O III], and [N II]—is mapped over five metallicities and nine ionization parameters to give pressure and density calibration curves, including the density range over which each ratio is useful and its temperature sensitivity.","core_discovery":"The central claim is that a new grid of self-consistent photoionization models, run at constant pressure and separately at constant density, provides reliable calibrations that turn emission-line ratios across the ultraviolet, optical, and infrared into ISM pressure and electron density for metallicities $\\log(\\mathrm{O/H})+12 = 7.63$ to $9.23$, pressures $4 < \\log(P/k) < 9$, and densities $0 < \\log(n_e/\\mathrm{cm^{-3}}) < 5$. The grids, presented in Tables 1 and 2 and Figures 5 through 17, show a large metallicity dependence: for a fixed line ratio the inferred pressure can shift by up to about an order of magnitude across the model metallicity range, and the traditional [S II] density estimate can shift by up to 0.4 dex. The paper further establishes, from the computed ionization structure, that [S II] is produced in the outer partially ionized zone while [O II] is produced closer to the ionizing source, so the two diagnostics should not be used interchangeably unless the electron temperature is known to be constant.","pith_inferences":["If the zoning claim is right, published high-redshift electron densities that mix [S II] and [O II] without metallicity corrections may contain systematic offsets that mimic or hide evolutionary trends; re-analysis with zone- and metallicity-aware calibrations could shift those measurements.","The same theoretical grids could be inverted to diagnose gas-phase metallicity from pressure-sensitive ratios when pressure is constrained independently, a use the paper does not explore.","The fixed iron depletion of -1.5 dex is a lever: if future observations show depletion varies with metallicity or redshift, the pressure calibrations can be rescaled, but the qualitative zone-separation result would likely survive.","The ratio combinations can be applied to Monte Carlo photoionization simulations of clumpy media to derive corrections for unresolved clumping, connecting the diagnostics to the paper's discussion of density structure."],"forward_implications":["ISM pressure can be measured directly from calibrated line ratios rather than inferred from density under a fixed-temperature assumption.","Electron-density estimates from [S II] should be corrected for gas-phase metallicity; differences up to 0.4 dex in density can arise solely from abundance changes.","[S II] and [O II] densities should not be averaged or used interchangeably because the two ratios probe different ionized zones within a nebula.","Combining ratios that trace different ionization zones lets observers map pressure stratification in resolved H II regions, from outer [N I] and [S II] zones to inner high-ionization zones.","The calibrations extend to UV and IR lines observable at high redshift, so future resolved observations can measure ISM pressure in early-universe galaxies."],"supporting_citations":[{"why":"Supplies the relative element abundance scaling used to set nebular abundances at every sub-solar metallicity in the grids.","marker":"Nicholls et al. (2017)"},{"why":"Provides the local B-star photospheric abundance measurements that anchor the adopted present-day metallicity scale.","marker":"Nieva & Przybilla (2012)"},{"why":"Provides the parametric dust depletion model used for the fixed -1.5 dex iron depletion in the nebulae.","marker":"Jenkins (2009)"},{"why":"One of the references describing the photoionization code used to compute the model grids.","marker":"Sutherland & Dopita (1993)"},{"why":"Describes the current version of the photoionization code and its dust and radiation-pressure treatment used for the models.","marker":"Dopita et al. (2013)"},{"why":"Supplies the atomic data, collision strengths, and transition probabilities used for the 30 elements in the models.","marker":"Del Zanna et al. (2015)"},{"why":"Provides the standard single-temperature [S II] and [O II] density relations that the new calibrations are compared against.","marker":"Osterbrock (1989)"},{"why":"Documents updated photoionization-based density calibrations and atomic data differences that motivate the new models.","marker":"Proxauf et al. (2014)"}],"fun_headline_variants":["Emission-line ratios now trace ISM pressure across cosmic time","Metallicity shifts pressure and density probes in H II regions","Why [S II] and [O II] give different electron densities","Pressure diagnostics for H II regions from UV to IR"],"cache_read_input_tokens":99200,"weakest_assumption_plain":"The load-bearing premise is that the adopted nebular abundance scaling—based on local B-star photospheric abundances and a fixed amount of iron depletion onto dust—is valid at all grid metallicities and at high redshift.","fun_headline_variants_meta":{"raw":{"variants":["Emission-line ratios now trace ISM pressure across cosmic time","Metallicity shifts pressure and density probes in H II regions","Why [S II] and [O II] give different electron densities","Pressure diagnostics for H II regions from UV to IR"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00056,"raw_usage":{"total_tokens":2731,"prompt_tokens":1086,"completion_tokens":1645,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":702,"completion_tokens_details":{"reasoning_tokens":1574}},"tokens_in":702,"tokens_out":1645,"duration_ms":11540,"temperature":1.0,"reasoning_tokens":1574,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:12:18.812939+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure [S II], [O II], and the electron temperature in a set of spatially resolved H II regions spanning metallicities $\\log(\\mathrm{O/H})+12 = 7.6$ to $9.2$. The claim predicts that at fixed electron temperature the two ratios give the same density and shift with metallicity by up to 0.4 dex; observing that they disagree at constant temperature by more than that would falsify the zoning and metallicity calibration.","supporting_citations":[{"cited_title":"Upgrading electron temperature and electron density diagnostic diagrams of forbidden line emission","cited_arxiv_id":"1311.5041","evidence_quote":"Documents updated photoionization-based density calibrations and atomic data differences that motivate the new models."}],"review_version":1}