{"id":"fe47ae41-0871-4892-885c-1af82289e85a","arxiv_id":"2501.06156","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In models, moderate external ultraviolet radiation reshapes the outer disk chemistry of a protoplanetary disk while leaving the innermost 25 au largely unchanged.","lead":"Protoplanetary disks around young stars are bathed in ultraviolet light from nearby stars, and this paper models what happens when that light is only moderately stronger than average, 10 to 100 times the galactic background. The models predict that the outer disk chemistry changes noticeably, especially for C+, N2H+, C, and CS, while the inner 25 au stays almost identical to an isolated disk.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing point is the Section 2.3 treatment of external UV penetration: the isotropic assumption and the scalar 3D-averaged optical depth determine where the outer ionization front sits, and a direction-resolved rerun of the 200 au/100 G0 model would test whether the species rankings are…","rationale":"The reader's weakest-assumption identification is correct: Section 2.3 is the sole channel by which the external UV field enters the chemistry, and all of the abstract's affected species are direct readouts of how far that field penetrates. My attack sharpens the same weakness into a checkable radiation-transfer issue. If the exact mean transmission is similar to the scalar approximation in this geometry, the concern does not land and the central claims stand under the stated assumption; if not, the quantitative ordering and emission ratios in Section 3.3 change. I do not view the missing H2(v) channel or radial drift as equally load-bearing because the paper explicitly flags them and they do not enter the abstract's primary mechanism, whereas the UV penetration treatment is used without sensitivity analysis. The inner-disk invariance claim is less vulnerable because the inner disk is optically thick to all plausible external geometries. Thus the reader's CONDITIONAL verdict is appropriate; I would attach this specific numerical check as an additional condition rather than change the verdict.","tokens_in":19516,"tokens_out":7519,"duration_ms":77271,"concrete_test":"Modify the external UV module in Section 2.3 for the 200 au disk at 100 G0: replace the scalar 3D-averaged optical depth with the exact angular integral I_ext = (1/4pi) ∫ I0(omega) exp[-tau(omega)] dOmega over the same 150 rays, and also run an anisotropic case with I0 nonzero only in the hemisphere facing the cluster center. Recompute the C+/C boundary, the CO photodesorption region, and the N2H+ and CS column densities between 50 and 200 au. If the front location or any column density shifts by more than a factor of ~2 relative to the published 100 G0 curve, the species rankings and proposed line ratios are geometry-dependent and must be restated as isotropic-field results or as upper limits.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that moderate external UV reconfigures outer disk chemistry while leaving the inner disk unchanged. The species most affected (C+, C, N2H+, CS, HCN) are set by the depth to which the external field penetrates, which is controlled entirely by the Section 2.3 treatment. Two approximations are load-bearing. First, the field is treated as isotropic. In the intended application (edge of a cluster such as the ONC), radiation arrives preferentially from the cluster center; isotropic averaging dilutes the low-column direction through which photons can actually reach the outer midplane, so the outer ionization front will move. Second, even for an isotropic field, the attenuated intensity is the angular average of the transmission, <exp(-tau)>, not generally exp(-<tau>). The text says a 3D-averaged optical depth is computed over 150 rays and applied to compute the extincted ISRF spectrum; if the code uses exp(-<tau>) (or an effective optical depth defined from <tau>), the result can differ by orders of magnitude in a flared disk where column density is strongly direction dependent. Because C+ production, CO photodissociation/photodesorption, N2H+ destruction, and CS/HCN responses all trace the UV penetration depth, the predicted species rankings and emission ratios inherit this uncertainty. This is not a stability question about the chemical network; it is the boundary condition that creates the phenomenon. No code or model outputs are released and no sensitivity run addresses this, so the central quantitative predictions are currently unverified against the radiation-transfer approximation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a systematic grid of 2D gas-grain chemical models of a 0.01 M_sun T Tauri disk (outer radii 50, 100, and 200 au) exposed to external interstellar radiation fields of G0 = 1, 10, 50, and 100. The external UV is treated as isotropic and attenuated by a 3D-averaged optical depth following Cleeves (2016). The central result is that moderate external UV produces an outer ionization front that significantly changes the abundances, column densities, and predicted line emission of C+, C, N2H+, CS, and HCN, while leaving the inner ~25 au essentially unchanged. The authors explicitly list the main simplifications (isotropic constant ISRF, no radial drift, no vibrationally excited H2, no accretion heating) and post-process a subset of species with LIME to suggest observational diagnostics, including C and the C/N2H+ and C/CS ratios.","tokens_in":19834,"tokens_out":6249,"duration_ms":60663,"significance":"The paper targets a genuinely understudied regime: G0 ~ 10-100, which is argued to be the most common disk-harboring environment. The model construction uses established codes and reaction networks, the G0 values and disk radii are inputs rather than fitted to the computed chemistry, and the predictions are falsifiable through disk-integrated line ratios and resolved line profiles. The explicit statement of limitations is a strength. The main risk is the Section 2.3 treatment of external UV penetration, which determines where the ionization front sits and therefore which species respond; if the averaging operation or the isotropy assumption changes the effective penetration depth, the quantitative rankings and diagnostic ratios would shift, even though the qualitative inner/outer dichotomy is likely robust.","major_comments":[{"comment":"The external UV penetration treatment is the load-bearing boundary condition for the entire study. The text says a 3D-averaged optical depth is estimated from 150 rays and that 'this optical depth is then applied to compute the extincted ISRF spectrum,' but it does not state whether the code evaluates exp(-<tau>) or the angular average <exp(-tau)>. In a flared disk these differ substantially, and for the intended cluster application the field is anisotropic rather than isotropic. Because C+ production, CO photodissociation/photodesorption, N2H+ destruction, and the CS/HCN responses all trace the UV penetration depth, the species rankings and the proposed diagnostic ratios inherit this uncertainty. Please report the exact averaging used and provide a sensitivity test (for example, a direction-resolved rerun of the 200 au/100 G0 model) to show that the qualitative and quantitative conclusions are robust to the choice.","section":"2.3"},{"comment":"The integrated fluxes in Table 3 do not track the reported total particle numbers in several rows. For CS at 50 au, the total particle number drops by roughly a factor of 5 from 1 G0 to 100 G0 while the integrated flux remains essentially constant (4.3 to 4.5 mJy km/s). For C2H at 200 au, the particle number rises from 4.1e38 to 5.5e38 while the flux is non-monotonic (2.9, 5.3, 3.3, 3.7 mJy km/s). Please check the LIME outputs or explain the optical-depth/excitation effects that create these plateaus; as presented, the table weakens the observational claims derived from the same spectra.","section":"3.3, Table 3"}],"minor_comments":[{"comment":"The text says the CS column at 175 au shows 'almost an order of magnitude decrease' while giving a ratio of ~2.8, and the C column at the edge is 'an order of magnitude higher' with a ratio of ~3; these descriptions should be reconciled with the quoted numbers.","section":"3.1.2, 3.1.3"},{"comment":"The caption states that C2H and HCN spectra include hyperfine components; please clarify whether CN also includes hyperfine structure, and if not, why not.","section":"Figure 10"},{"comment":"The phrase 'gas columm density' contains a typo.","section":"3.1 (near Figure 3)"},{"comment":"Releasing the chemical model outputs or input files would improve reproducibility; the current text refers to prior codes for details but does not provide enough information to rerun the exact 150-ray averaging.","section":"2.3, reproducibility"}],"recommendation":"major_revision","confidential_remarks":"The isotropy/averaging issue in Section 2.3 is the main scientific risk; I would be willing to accept after a sensitivity analysis. The self-citations to Cleeves (2016) and Anderson et al. (2021) are appropriate because the methods are from those papers, and I do not see a citation-pattern concern. The manuscript fits an astrochemistry/disk evolution journal well."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read of Gross & Cleeves. The genuinely new thing is systematic modeling of disk chemistry in the G0=10–100 regime, between isolated disks and the extreme cluster center. They keep the disk mass fixed and vary the outer radius, which cleanly separates density-structure effects from UV strength. The main qualitative finding—outer disk chemistry is reworked by an external UV ionization front, while the inner ~25 au is nearly unchanged—is credible and likely robust. The synthetic spectra and proposed line ratios (C/N2H+, C/CS) are sensible, observationally actionable diagnostics.\n\nThe paper is honest about its limits: no radial drift, no vibrationally excited H2, isotropic ISRF, no accretion heating. It uses established codes (TORUS, Bethell & Bergin, the Fogel/Cleeves/Anderson network) and proper self-shielding for CO, H2, and N2. The comparisons to Walsh et al. at high G0 and to ONC observations are useful.\n\nThe real soft spot is the external UV penetration treatment in Section 2.3. They compute a 3D-averaged optical depth over 150 lines of sight and apply that scalar to attenuate the ISRF. The paper does not say whether they average the optical depth and exponentiate, i.e. exp(-<tau>), or average the transmission, <exp(-tau)>. In a flared disk the two can differ by orders of magnitude along low-column sightlines. Since the positions of the C+/C/N2H+/CS/HCN response features all trace the penetration depth, the quantitative column ratios inherit that uncertainty. The isotropic assumption is also acknowledged, but at a cluster edge the radiation comes preferentially from one direction; a direction-resolved run would likely shift the outer ionization front. These are not fatal to the qualitative story, but they mean the specific numbers should not be read as precise predictions. The stress-test note is on target here.\n\nA second issue is reproducibility: no code or model outputs are released. For a modeling paper this is a real hinderance, and it amplifies the penetration uncertainty because no one can rerun the 100 G0/200 au case with a different angular treatment. A sensitivity run on the UV averaging and on photodesorption yields would go a long way.\n\nOverall, this is a solid parameter-space exploration that deserves peer review. The authors should be asked to clarify or patch the UV transmission issue and, ideally, tabulate or release the output columns. I would not desk-reject it. For a reading group, it is worth a slot if anyone cares about cluster disk chemistry, though it is not a landmark.","headline":"A useful first systematic look at disk chemistry under moderate external UV; the outer-disk effects are probably real, but the UV penetration treatment needs a sensitivity test before the quantitative predictions are trusted.","tokens_in":20389,"tokens_out":3492,"would_cite":true,"duration_ms":35087,"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":"Moderate external ultraviolet, 10-100 times the galactic average, redraws the outer chemistry of planet-forming disks while leaving the inner 25 au effectively unchanged.","keywords":["protoplanetary disks","astrochemistry","external UV irradiation","interstellar radiation field","chemical modeling","disk structure","ionization front","molecular line diagnostics"],"falsifier":"Map the radial and vertical distribution of C I 1-0, N2H+ 3-2, and CS 6-5 emission in a moderately irradiated (G0 ~ 10-100) disk at a cluster edge using a submillimeter interferometer. The model predicts C brightens while N2H+ and CS fade beyond the outer tens of au in an extended disk and stay flat in a compact disk; an observation showing no such structural dependence, or the same trends in all disks, would contradict the ionization-front picture.","tokens_in":19303,"feed_emoji":"🪐","tokens_out":8111,"duration_ms":74710,"temperature":0.7,"pith_summary":"The paper asks what a modestly elevated ultraviolet background, 10 to 100 times the galactic average and typical of the outskirts of a young cluster like Orion, does to the chemistry of a planet-forming disk. Using time-dependent gas-grain models of T Tauri disks with three outer radii and a fixed mass, the authors find that the external UV creates an outer ionization front that raises the abundances and column densities of C and C+ while suppressing N2H+, CS, and HCN, with CO and HCO+ barely changing. The effect depends strongly on disk structure: a compact, dense 50 au disk shields its chemistry, whereas a more extended 200 au disk lets the UV penetrate and reprocess the cold outer reservoir. The inner 25 au remains chemically indistinguishable from an isolated disk, suggesting that even in clusters the planet-forming zone can be studied using nearby quiescent systems. Synthetic spectra suggest that C and the C/N2H+ and C/CS ratios are usable diagnostics of external UV exposure.","feed_headline":"Moderate UV rewires outer disk chemistry, spares inner 25 au","feed_subtitle":"Model shows cluster-edge starlight reshapes outer C, N, S chemistry while the inner planet-forming zone stays calm.","key_machinery":"The outer ionization front is the central mechanism: the boundary in the outer disk where external FUV photons, attenuated to each point by an optical depth averaged over 150 lines of sight across 4π steradians, become strong enough to drive photodissociation and photodesorption faster than freeze-out and recombination. This averaged extinction treatment converts the isotropic ISRF into local photo-rates within the gas-grain chemical network, and it is what couples disk structure (outer radius and density) to the species response curves.","core_discovery":"On its own terms, the paper establishes that moderate external UV irradiation (G0 = 10-100) is not a scaled-down version of the extreme Orion-center case: it produces a distinct outer ionization front that reconfigures the cold disk chemistry rather than simply eroding it. Photodissociation and photodesorption of CO in the outer, low-density regions release carbon, boosting C and C+, while N2H+ is destroyed both by the released CO and by increased electron recombination, and CS and HCN are depleted; the CN/HCN column density ratio rises as HCN is photodissociated. CO itself keeps a nearly constant column density despite a changing midplane abundance, and a secondary CO snowline appears at the outer edge of the 200 au model from UV-driven photodesorption. The same physical models show that the inner disk within roughly 25 au is optically thick to the external field, so its ionization is governed by the star's own UV and X-rays; consequently, inner-disk chemistry should resemble that of isolated disks. The magnitude of every effect is set by the disk's density structure, so compact disks are nearly immune while extended disks are strongly altered.","pith_inferences":["If the isotropic assumption is relaxed to a single illuminating source, the outer ionization front would likely sit deeper on the side facing the massive stars and shallower on the far side, producing azimuthally asymmetric species distributions that the current models cannot capture.","The result that the inner 25 au is buffered suggests the volatile inventory feeding planet formation is set by the host star rather than the birth environment; extending the models to include radial drift of photodesorbed material could link the altered outer reservoir to inner-disk composition.","A quick observational test would be to compare the C/N2H+ ratio in a sample of cluster-edge disks against the predicted monotonic rise with G0.","The prediction that HCO+ is insensitive to moderate UV means its brightness in externally irradiated disks cannot by itself indicate environmental irradiation, so surveys relying on HCO+ detections may be selecting disks whose chemistry is set internally."],"forward_implications":["Observers can use the C/N2H+ and C/CS line flux ratios as a gauge of external UV exposure, because C brightens while N2H+ and CS fade as G0 rises.","N2H+ loses reliability as a CO-snowline tracer under moderate external UV, since its column-density peak shifts inward by 10-20 au as electron recombination suppresses the outer edge.","Inner-disk (<25 au) chemistry in cluster disks should match that of isolated disks like those in Taurus and Lupus, so nearby quiescent systems can serve as proxies for planet-forming chemistry in clusters.","CO column density is a poor UV diagnostic; only its emitting height and midplane abundance respond, so CO intensity alone will not reveal external irradiation.","Compact, dense disks (like those seen near the ONC center) are the most resistant to external UV, meaning observed chemical differences in those disks likely come from other environmental factors."],"supporting_citations":[{"why":"It supplies the 3D-averaged optical depth method used to attenuate the external ISRF at every disk point.","marker":"Cleeves (2016)"},{"why":"It supplies the three constant-mass disk models and the updated photodesorption yields used in the chemistry.","marker":"Anderson et al. (2021)"},{"why":"It is the base 2D gas-grain chemical network from which the abundance models are evolved.","marker":"Fogel et al. (2011)"},{"why":"It adds the secondary CO snowline and photodesorption behavior, and benchmarks the isolated-disk G0 levels.","marker":"Cleeves et al. (2016b)"},{"why":"It is the extreme-UV comparison model that motivates the CN/HCN diagnostic and the moderate-field extension.","marker":"Walsh et al. (2013)"},{"why":"It is the non-LTE excitation and radiative transfer code used to produce the synthetic line spectra.","marker":"Brinch & Hogerheijde (2010)"},{"why":"It provides the molecular collision rate coefficients used in the radiative transfer calculations.","marker":"Schöier et al. (2005)"},{"why":"It supplies the wavelength-dependent photoionization and photodissociation rates driving the UV chemistry.","marker":"van Dishoeck et al. (2006)"}],"fun_headline_variants":["Mild UV from cluster edges reshapes outer disk chemistry","Moderate starlight triggers outer ionization front in disks","Outer disk chemistry shifts under moderate UV, inner zone calm","Cluster-edge UV rewires disk chemistry beyond inner 25 au","Moderate external UV reorders disk chemistry, inner disk safe"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The models assume the external UV arrives isotropically and is attenuated by a direction-averaged optical depth, whereas real cluster radiation comes from specific massive stars, so the true depth and shape of the ionization front could differ.","fun_headline_variants_meta":{"raw":{"variants":["Mild UV from cluster edges reshapes outer disk chemistry","Moderate starlight triggers outer ionization front in disks","Outer disk chemistry shifts under moderate UV, inner zone calm","Cluster-edge UV rewires disk chemistry beyond inner 25 au","Moderate external UV reorders disk chemistry, inner disk safe"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00024,"raw_usage":{"total_tokens":1575,"prompt_tokens":1062,"completion_tokens":513,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":678,"completion_tokens_details":{"reasoning_tokens":429}},"tokens_in":678,"tokens_out":513,"duration_ms":5747,"temperature":1.0,"reasoning_tokens":429,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T21:05:30.065468+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Map the radial and vertical distribution of C I 1-0, N2H+ 3-2, and CS 6-5 emission in a moderately irradiated (G0 ~ 10-100) disk at a cluster edge using a submillimeter interferometer. The model predicts C brightens while N2H+ and CS fade beyond the outer tens of au in an extended disk and stay flat in a compact disk; an observation showing no such structural dependence, or the same trends in all disks, would contradict the ionization-front picture.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It supplies the 3D-averaged optical depth method used to attenuate the external ISRF at every disk point."}],"review_version":1}