{"id":"a50221dd-0120-4c0a-850b-200db76cb51c","arxiv_id":"2411.14253","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The BE99 diagnostic can be extended with more line ratios and higher ionisation fractions, but for dense jets like Par Lup 3-4 a multi-line fit gives different and more consistent gas parameters.","lead":"This paper extends the classic BE99 method for reading gas conditions in protostellar jets by adding more emission lines and allowing higher ionisation states. Applying it to the dense outflow Par Lup 3-4 shows the old method gives inconsistent answers, so the authors recommend a multi-line model fit instead.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Par Lup 3-4 high-density solution rests on UVB/NIR [S II] flux ratios that the paper itself finds inconsistent with any extinction law; until this calibration problem and the neglected H-atom collisions are tested, the factor-10 ne claim is not secure.","rationale":"The reader's conditional verdict is appropriate, and the weakest-assumption label points in the right direction: the five-level, hydrogen-collision-free excitation model is the main source of model dependence. I would sharpen the concern to a more specific and more actionable weakness. The paper's own Section 6.2.2 flags an unresolved inconsistency in the observed UVB/NIR [S II] ratios, and those are precisely the ratios that push the multi-line fit to ne ~ 5e4 cm^-3. Since the central claim is that the classical BE99 method underestimates ne by a factor of ten in dense gas, a calibration error or missing H-atom collisions in the very lines that establish the high-density solution would directly invalidate that claim. The proposed test would settle this by showing whether the high-ne solution survives removal or rescaling of the suspect ratios, and by checking the relative flux calibration with an independent standard. The paper's methodological contributions (time-domain integration, extended BE99e stripes, and the demonstration that BE99 converges faster than the hydrogen recombination time) are not affected by this concern and remain valuable. Therefore the reader's CONDITIONAL verdict stands; the analysis should be accepted only after the calibration and H-collision sensitivity checks are addressed.","tokens_in":33390,"tokens_out":7296,"duration_ms":70465,"concrete_test":"Re-run the multi-line fit (Eq. 51) on the Table A.8 fluxes under two perturbations: (i) exclude the UVB/NIR [S II] ratios involving [S II] lambda4068/4076 and [S II] lambda10287/10320/10336/10370 from the fit, and (ii) multiply the NIR [S II] fluxes by gain factors of 1.3 and 0.7 to simulate an unresolved calibration offset. If the best-fit ne moves below 2e4 cm^-3, or if Te and xe change by more than the quoted ranges, the claimed factor-10 discrepancy in ne is not robust. A complementary observational check is to re-reduce the archival X-Shooter UVB and NIR arms (programme 085.C-0238(A)) using a same-night telluric standard to independently verify the UVB-to-NIR relative flux calibration, which would settle whether the Miller-ratio inconsistency is instrumental.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline result for Par Lup 3-4 is the multi-line fit (ne = 45,000-53,000 cm^-3, Te = 7,600-8,000 K, xe = 0.027-0.036), which is also the evidence that the classical BE99 method underestimates ne by roughly a factor of ten. That result is not anchored by an independent diagnostic: the fit minimizes Eq. 51 over all line ratios, and the high-density branch is driven by the [S II] lambda4068/lambda6731 and [S II] lambda4068/lambda10320 ratios from Table A.8. Section 6.2.2 reports that the same UVB/NIR [S II] ratios are inconsistent with any reddening law ('no extinction value is consistent with the observed Miller line ratios'), leaving the cause unresolved and suggesting either a flux-calibration offset between the UVB and NIR arms or an unmodeled excitation effect. The grid fit can absorb this inconsistency into a high-density model, since [S II] lambda4068/6731 rises with ne while the UVB/NIR [S II] ratios are also density-sensitive. A second load-bearing assumption is the excitation model itself: at the fitted xe = 0.03 and ne = 5e4 cm^-3, the neutral hydrogen density is n(H0) ~ 1.6e6 cm^-3, yet Appendix A's five-level model includes only electron collisions; Section 5.2 concedes that H-atom collisions are missing and affect line emissivities 'to an unknown amount'. Either the flux-calibration problem or the missing H collisions could shift the fitted ne by more than the quoted 45,000-53,000 range, and neither is tested in the paper.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript revisits the BE99 diagnostic method (Bacciotti & Eislöffel 1999), which derives electron density ne, electron temperature Te, and hydrogen ionisation fraction xe in protostellar jets from three forbidden-line ratios. The authors (i) integrate the BE99 ionisation network in the time domain and conclude that the BE99 method converges on timescales shorter than the hydrogen recombination time even before the reaction equilibrium is reached; (ii) propose extensions ('BE99e', 'BE99e+') that add blue and near-infrared line ratios and higher ionisation states; and (iii) apply the scheme to two objects. For the low-excitation outflow Par Lup 3-4, the classical method yields ne ≈ 4.3×10^3 cm^-3, Te ≈ 3.1×10^4 K, xe ≈ 3×10^-3, while a multi-line excitation-model fit yields ne ≈ 4.5–5.3×10^4 cm^-3, Te ≈ 7.6–8.0×10^3 K, xe ≈ 0.027–0.036; the paper concludes that the classical BE99 method underestimates ne by roughly an order of magnitude because the optical [S II] lines are quenched at high density. For the 244-440 Proplyd, the extended diagram gives xe = 0.58 ± 0.05 at knot E3 without using [N II] lines.","tokens_in":33697,"tokens_out":47212,"duration_ms":399744,"significance":"If the factor-of-ten claim holds, it is a genuinely important result: the overlapping-stripe solution of the BE99 method is widely used as the adopted gas state, and Par Lup 3-4 would be a documented counterexample in dense gas near the driving source. The paper has real strengths that should be credited. The reaction network and complete atomic data are tabulated in the appendices, making the calculations reproducible; the 192-model convergence study (Section 2.6) is systematic; the argument against the hot classical solution based on the [O I]λ5577 non-detection is sharp and falsifiable; and the Par Lup 3-4 non-overlap of the stripes is an observational result, not an artifact of synthetic spectra, so the core of the method critique is not circular. The BE99e extension is practical for X-Shooter/MUSE-era data. The weakness is that the quantitative Par Lup 3-4 parameters rest on a five-level, electron-collision-only excitation model and on a flux set that the authors themselves find internally inconsistent; as argued in the major comments, the factor-of-ten magnitude is plausible but not yet secured.","major_comments":[{"comment":"The multi-line fit that produces the headline parameters is performed on a flux set that the authors themselves flag as internally inconsistent, and the reported fit quality is contradicted by the paper's own atomic data. In the five-level model, [S II]λ4068.6 and [S II]λ10320.5 share the upper level 2Po3/2, so with the A-values of Table A.5 their emissivity ratio is the fixed branching ratio A(4068.6)E(4068.6)/A(10320.5)E(10320.5) ≈ 3.1 at every grid point. The observed ratio in Table A.8 is 39.48/6.06 ≈ 6.5, a residual of ≈0.32 dex, i.e., about 5σ at the quoted errors, which no (ne, Te, xe) can remove. Consequently, the statement in Section 6.2.6 that 'all observed lines are consistent with the best excitation model within 10%' cannot be correct, and the quoted parameters (ne ≈ 45,000–53,000 cm^-3, Te = 7,600–8,000 K, xe ≈ 0.027–0.036) are at least partly a least-squares compromise that absorbs the UVB/NIR inconsistency described in Section 6.2.2. Testing non-zero extinction cannot fix this either, since the observed ratio exceeds the intrinsic branching ratio and would require negative AV. The density result is plausibly anchored by the UVB/VIS ratios ([S II]λ6716/λ6731 ≈ 0.50, [S II]λ4068/λ4076 ≈ 3.6, [S II]λ4068/(λ6716+λ6731) ≈ 1.5), which are mutually consistent and agree with W14; the problem is the global fit's quality claim and the derived Te and xe, which absorb the unmodelable NIR ratios. The authors should either resolve the calibration problem (e.g., telluric correction near 1.03 μm, independent UVB–NIR flux verification) or show explicitly that the fitted parameters are unchanged when the offending ratios are excluded.","section":"Sections 6.2.2 and 6.2.6; Tables A.5 and A.8; Eq. (51)"},{"comment":"The headline parameters are derived from an excitation model that includes only electron-impact collisions, but at the fitted conditions (ne ≈ 5×10^4 cm^-3, xe ≈ 0.03), Eq. (22) gives n(H0) ≈ 1.6×10^6 cm^-3. With neutral hydrogen this abundant, H0-impact excitation is likely comparable to or dominant over electron-impact excitation for several fitted lines (notably [O I] and [N I], whose upper states have known H-collision channels), even for rate coefficients near 10^-10 cm^3 s^-1. Section 5.2 concedes that the missing H-atom collisions affect line emissivities 'to an unknown amount', and no such collisions appear in the five-level model of Appendix A. Because the same model sets the positions of the BE99e stripes and defines the multi-line fit, the derived Te and xe, and hence the normalisation of the density comparison with the classical BE99 result, rest on an acknowledged but unquantified approximation. A sensitivity test (e.g., adding H-collision rates of plausible magnitude and re-fitting, or at least estimating the effect on the [O I] and [N I] lines) is required before the BE99 bias can be attributed specifically to [S II] quenching.","section":"Section 5.2; Appendix A (Eqs. A.40–A.41); Eq. (22)"},{"comment":"The fitted multi-line solution is in tension with the paper's treatment of equilibrium. Eq. (25) and Fig. 2 give xeq_e ≈ 3–4×10^-4 at Te = 7,600–8,000 K, whereas the best fit has xe ≈ 0.027–0.036 at the same temperature, about two orders of magnitude above the reaction-equilibrium value. Section 6.2.6 nevertheless rules out non-equilibrium as the cause of the stripe mismatch ('neither extinction nor non-equilibrium effects are likely responsible'), and Section 6.2.4 argues that the gas 'should have had enough time to stay close enough to the equilibrium'. If the gas were close to the BE99 equilibrium, the fitted xe at Te ≈ 7,800 K should be near 4×10^-4; if the observed ratios genuinely require xe ≈ 0.03 at low Te, then the gas is far from that equilibrium, and non-equilibrium is a live alternative (or additional) explanation for the classical BE99 bias. The dynamical-age estimate of ≈2 yr in Section 6.2.4 is also shorter than the fiducial convergence time τBE ≈ 3 yr of Section 2.6 (with the 192-model spread extending to ≈32 yr), so the timescale argument does not clearly exclude this possibility. The attribution of the BE99 failure to quenching alone is therefore not uniquely supported and should be discussed explicitly.","section":"Sections 6.2.4–6.2.6; Eq. (25); Fig. 2"}],"minor_comments":[{"comment":"'We aim to extent the BE99 method' should read 'extend'.","section":"Abstract; Section 1"},{"comment":"The abstract and conclusions state that 'the BE99 equilibrium is reached faster than the hydrogen recombination time', which contradicts the quantitative example in Section 2.5 (τeq ≈ 8×10^4 yr against τrec ≈ 200 yr for HH34 knot J). This should be reworded to refer explicitly to the BE99-method convergence time τBE introduced in Section 2.6.","section":"Abstract; Section 7"},{"comment":"'a visual extinction ABE V in Eq. 3' should refer to Eq. (30), the reddening relation; Eq. (3) is a charge-exchange reaction.","section":"Section 3.1"},{"comment":"The line labels '[OvII]λ7320' and '[SvIII]λ9530' are typos for [O II] and [S III].","section":"Table A.8"},{"comment":"The quoted parameter ranges (ne ≈ 45,000–53,000 cm^-3, Te = 7,600–8,000 K, xe ≈ 0.027–0.036) are not derived from any stated statistical criterion: the grid spacings are Δne = 100 cm^-3, ΔTe = 100 K, Δxe = 0.001, and no ΔC threshold for Eq. (51), confidence region, or error propagation from Table A.8 is given. The ranges should be defined or stated as approximate.","section":"Section 6.2.6"},{"comment":"The claim that the Par Lup 3-4 gas 'should have had enough time to stay close enough to the equilibrium' is not quantified; given the ≈2 yr dynamical age computed in the same paragraph and the fiducial τBE ≈ 3 yr of Section 2.6, the authors should provide a τBE estimate at the fitted ne ≈ 5×10^4 cm^-3.","section":"Section 6.2.4"},{"comment":"The synthetic demonstrations (Figs. A.3–A.7) generate the spectra with the same model that is then used to interpret them, so the overlap of stripes in the equilibrium snapshots is in part a self-consistency check; the manuscript should state this explicitly for the extended diagrams, as it does implicitly for the classic stripes.","section":"Sections 3.2 and Appendix C"},{"comment":"'may not be accessible able to any BE99 method' contains a typo ('accessible able').","section":"Section 5.2"}],"recommendation":"major_revision","confidential_remarks":"The most serious issue is Major Comment 1: the claimed fit quality ('within 10%') is verifiably inconsistent with the paper's own atomic data (Table A.5) and measured fluxes (Table A.8), because the [S II]λ4068/λ10320 ratio is a fixed branching ratio. This should be checked carefully before the numbers propagate, since either the A-values, the flux calibration, or the reported fit quality is wrong. The remaining concerns (H-collision omission, equilibrium interpretation) are requests for tests that are within the authors' reach. The appendices and the reproducible model description are strengths of the manuscript."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear X,\n\nThis paper is a genuinely useful extension of the BE99 diagnostic, not just a repackaging. The authors integrate the BE99 reaction network in time and show that the method converges (tau_BE ~ 10^8 s) well before hydrogen recombination equilibrium – that's new, and it actually strengthens the case for using BE99 in outflows. The extension to extra line-ratio stripes (BE99e/BE99e+) is straightforward and sensible, and the test on the 244-440 proplyd shows you can get xe even when [N II] is missing. The re-reduction of the Par Lup 3-4 X-Shooter data also catches a real flux-calibration inconsistency in Whelan et al. (2014); that kind of archival housekeeping is worth having.\n\nThe soft spot is the headline Par Lup 3-4 result. The multi-line fit gives ne ~ 5e4 cm^-3, Te ~ 7800 K, xe ~ 0.03, a factor-ten change in density from classic BE99. But the fit is driven by [S II] 4068/6731 and the UVB/NIR [S II] ratios, and Section 6.2.2 reports that those same ratios are inconsistent with any extinction law. The authors are upfront about this – they say no extinction value works – but they don't resolve it. An unmodeled flux offset between UVB and NIR arms would be absorbed by the grid into a high-density solution. On top of that, the excitation model includes only electron collisions; Section 5.2 concedes H-atom collisions are missing and could affect emissivities \"to an unknown amount\". At the fitted xe and ne, n(H0) is large enough that H collisions plausibly matter. So the factor-ten claim is not secure.\n\nThe synthetic-diagram demonstrations are partly circular, as the authors acknowledge: the stripes are built from the same model used to interpret them. That's fine for illustration; the Par Lup non-overlap is a real observational falsification.\n\nBottom line: as a methods paper, it deserves a serious referee. The time-domain result and the extended stripe library are solid contributions. But the Par Lup application should be presented as model-dependent until the calibration inconsistency and H-collision physics are addressed. I'd recommend acceptance after major revision.","headline":"Useful methodological extension of BE99 with a genuinely new timescale result, but the headline Par Lup 3-4 density claim is not yet anchored because the driving [S II] ratios are internally inconsistent and the excitation model omits H-atom collisions.","tokens_in":34280,"tokens_out":2589,"would_cite":true,"duration_ms":22912,"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":"The paper argues that the standard six-line BE99 diagnostic underestimates electron density in dense jet gas by roughly a factor of ten, and that a multi-line extension recovers the true conditions.","keywords":["protostellar jets","forbidden emission lines","BE99 diagnostic","ionisation fraction","electron density","line quenching","shock-excited gas","X-Shooter spectroscopy"],"falsifier":"Measure the [O II] 3726/3729 doublet and the [O I] 5577 line in Par Lup 3-4 with deeper spectroscopy: densities near 50,000 cm$^{-3}$ and an absent 5577 line would support the quenching story, while densities near the classical 4,000 cm$^{-3}$ or a detected 5577 line would refute it.","tokens_in":33134,"feed_emoji":"🔭","tokens_out":9794,"duration_ms":83928,"temperature":0.7,"pith_summary":"The paper sets out to test and extend the BE99 method, a standard diagnostic that uses six bright forbidden lines of [S II], [N II], and [O I] to read electron density, temperature, and hydrogen ionisation fraction in shocked protostellar jets. The authors find that in dense gas the classical method can be badly wrong: for the outflow Par Lup 3-4, the classic six lines give $n_e \\approx 4\\,300$ cm$^{-3}$ and $T_e \\approx 31\\,000$ K, while a multi-line excitation model fit gives $n_e = 45\\,000$\\,--\\,$53\\,000$ cm$^{-3}$, $T_e = 7\\,600$\\,--\\,$8\\,000$ K, and $x_e = 0.027$\\,--\\,$0.036$. They trace the discrepancy to quenching of the optical [S II] lines at densities near their critical density, which biases the standard ratios. They also show that the BE99 equilibrium is reached faster than the hydrogen recombination time, and they extend the stripe-diagram method to additional blue and near-infrared lines and to higher ionisation fractions ($x_e > 0.3$). The payoff would be a more reliable way to measure gas conditions in jets, especially close to the driving source where densities are high.","feed_headline":"Classic jet-gas diagnostic underreads density tenfold","feed_subtitle":"Multi-line fits put Par Lup 3-4 near 50,000 cm⁻³, not ~4,300; optical sulfur lines are quenched.","key_machinery":"The central object is the BE99 diagnostic diagram: stripes in the ($x_e$, $T_e$) plane drawn from ratios of forbidden lines, generated by a five-level collisional-excitation model that shares the ionisation network of the original method. The paper's machinery is a time-dependent reaction network of hydrogen, oxygen, nitrogen, and sulphur, including charge exchange, collisional ionisation, and radiative or dielectronic recombination, integrated until equilibrium, with synthetic spectra reddened and fed back into the BE99 stripes. This same network is extended with O$^{2+}$, N$^{2+}$, S$^{2+}$, and S$^{3+}$ for higher ionisation fractions, and the five-level model supplies emissivities for lines from 3500 to 11000 Å. The load-bearing piece is the consistency check this enables: if all stripes cross in one place, the assumed equilibrium and extinction are consistent; if they do not, as for Par Lup 3-4, the classical result is suspect and a multi-line grid fit is used instead.","core_discovery":"In the classical BE99 diagram, three observed line ratios OI/NII, OI/SII, and NII/SII each trace a stripe in the ($x_e$, $T_e$) plane, and their common overlap is taken as the gas state. The paper's central discovery is that this overlap is not a reliable check: the three classical stripes always meet even when the gas is out of equilibrium or extinction is misjudged, and adding one extra stripe from [N I] 5198+5200 to the Par Lup 3-4 data makes the stripes fail to intersect. Using a grid of five-level excitation models and fitting all observed lines at once, the authors obtain $n_e = 45\\,000$\\,--\\,$53\\,000$ cm$^{-3}$, $T_e = 7\\,600$\\,--\\,$8\\,000$ K, and $x_e = 0.027$\\,--\\,$0.036$, which they state substantially differ from the classical BE99 values. The explanation offered is quenching of the optical [S II] lines in the high-density gas, which makes the classical method underestimate the electron density by roughly an order of magnitude. For the higher-excitation proplyd 244-440, the extended method recovers the ionisation fraction $x_e = 0.58 \\pm 0.05$ at knot E3 from [S II], [O I], and [O II] ratios even though [N II] is unavailable.","pith_inferences":["The same stripe-overlap test could serve as a cheap validation step for archival X-Shooter and MUSE spectra: if the classical three stripes overlap but an extra stripe does not, no parameter set from the six-line method should be quoted without a caveat.","The factor-of-ten density bias implies that published jet surveys that used BE99 in dense knots may have systematically underestimated electron densities, and with them the ionisation fractions and mass-flux estimates that depend on those densities.","Because the blue [O II] 3726/3729 doublet has a critical density comparable to the optical [S II] lines, the paper's own caveat suggests that in the very densest gas near the source none of the optical BE99-type diagnostics are usable; near-infrared [Fe II] or [S II] 1.03 μm ratios would be the natural test.","The unresolved [S II] blue-versus-near-infrared inconsistency in Par Lup 3-4 points to a calibration or excitation effect beyond extinction; a cross-check with a different instrument would settle whether the blue [S II] lines are over-luminous or the near-infrared ones under-luminous."],"forward_implications":["Any BE99 result obtained in a jet knot with densities near or above the critical density of the optical [S II] lines should be checked against another density diagnostic before being trusted.","Observers can now add stripes from [N I] 5198+5200, [S II] 4068+4076, [O II] 3726/3729, and near-infrared [S II] and [N I] lines to the same ($x_e$, $T_e$) diagram; failure of the stripes to overlap flags bad extinction, non-equilibrium gas, or an inconsistent parameter set.","The equilibrium assumption is safer than previously thought: in the tested models the BE99 method converges within about 3 years, well before the hydrogen recombination time, even when the full chemical equilibrium is far off.","The extension to higher ionisation states lets the method work on high-excitation objects where [N II] is missing, as demonstrated by the ionisation fraction derived for Proplyd 244-440.","For Par Lup 3-4 the multi-line fit implies an electron density about six times higher and a temperature about three times lower than previously reported, so outflow models built on the older BE99 values would need to be re-evaluated."],"supporting_citations":[{"why":"Defines the BE99 stripe method whose assumptions, time scales, and failure modes this paper tests and extends.","marker":"Bacciotti & Eislöffel 1999"},{"why":"Supplies the earlier Par Lup 3-4 BE99-derived gas parameters and X-Shooter spectra that this paper re-reduces and questions.","marker":"Whelan et al. 2014"},{"why":"Provides the multi-line excitation-model fitting approach used to obtain the new Par Lup 3-4 parameters.","marker":"Hartigan & Morse 2007"},{"why":"Supplies critical densities and diagnostic diagrams used to identify quenching of the optical sulphur and oxygen lines.","marker":"Giannini et al. 2019"},{"why":"Provides the MUSE line fluxes of the 244-440 proplyd used in the high-excitation application.","marker":"Kirwan et al. 2023"},{"why":"Supplies charge-transfer rate coefficients for nitrogen and sulphur in the ionisation network.","marker":"Kingdon & Ferland 1996"},{"why":"Supplies the oxygen charge-transfer rates used to keep the oxygen-hydrogen coupling accurate in the network.","marker":"Stancil et al. 1999"},{"why":"Provides the collisional ionisation and recombination rate expressions for the time-dependent network.","marker":"Landini & Monsignori Fossi 1990"},{"why":"Provides the S II collisional strengths that determine the optical and near-infrared sulphur line emissivities, which are central to the quenching argument.","marker":"Tayal & Zatsarinny 2010"}],"fun_headline_variants":["BE99 stripes always cross, hiding gas density error","Classic jet diagnostic quenched: density 10x higher","Multi-line fits reveal jets denser than BE99 says","Old method's overlap is a trap: jets denser","Adding one line breaks BE99 overlap, fixes density"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The analysis stands on the atomic model's predicted line ratios being accurate, because both the stripe test and the multi-line fit would shift if the rate coefficients, critical densities, or missing hydrogen-atom collisions were wrong.","fun_headline_variants_meta":{"raw":{"variants":["BE99 stripes always cross, hiding gas density error","Classic jet diagnostic quenched: density 10x higher","Multi-line fits reveal jets denser than BE99 says","Old method's overlap is a trap: jets denser","Adding one line breaks BE99 overlap, fixes density"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000252,"raw_usage":{"total_tokens":1745,"prompt_tokens":1317,"completion_tokens":428,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":933,"completion_tokens_details":{"reasoning_tokens":348}},"tokens_in":933,"tokens_out":428,"duration_ms":5135,"temperature":1.0,"reasoning_tokens":348,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:21:48.612626+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the [O II] 3726/3729 doublet and the [O I] 5577 line in Par Lup 3-4 with deeper spectroscopy: densities near 50,000 cm$^{-3}$ and an absent 5577 line would support the quenching story, while densities near the classical 4,000 cm$^{-3}$ or a detected 5577 line would refute it.","supporting_citations":[{"cited_title":"& Eisl \\\"o ffel , J","cited_arxiv_id":null,"evidence_quote":"Defines the BE99 stripe method whose assumptions, time scales, and failure modes this paper tests and extends."},{"cited_title":"T., Bonito , R., Antoniucci , S., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the earlier Par Lup 3-4 BE99-derived gas parameters and X-Shooter spectra that this paper re-reduces and questions."},{"cited_title":"& Morse , J","cited_arxiv_id":null,"evidence_quote":"Provides the multi-line excitation-model fitting approach used to obtain the new Par Lup 3-4 parameters."},{"cited_title":"2019, , 631, A44","cited_arxiv_id":null,"evidence_quote":"Supplies critical densities and diagnostic diagrams used to identify quenching of the optical sulphur and oxygen lines."},{"cited_title":"F., Whelan , E","cited_arxiv_id":null,"evidence_quote":"Provides the MUSE line fluxes of the 244-440 proplyd used in the high-excitation application."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies charge-transfer rate coefficients for nitrogen and sulphur in the ionisation network."},{"cited_title":"& Monsignori Fossi , B","cited_arxiv_id":null,"evidence_quote":"Provides the collisional ionisation and recombination rate expressions for the time-dependent network."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the S II collisional strengths that determine the optical and near-infrared sulphur line emissivities, which are central to the quenching argument."}],"review_version":1}