{"id":"7d47dd0b-090a-4d61-a476-bb397e7a8244","arxiv_id":"2608.08275","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Electron densities from [Ar IV] lines in 16 type 2 QSOs imply outflow mass rates and kinetic powers roughly 1 dex lower than assuming n_e=1000 cm^-3, plus a positive density-velocity correlation above 630 km/s.","lead":"Using VLT/MUSE spectra of 18 nearby obscured quasars, this paper measures gas densities with the [Ar IV] doublet and finds that ionized outflow rates and kinetic powers are about ten times lower than usually assumed. The study also reports a correlation between gas density and outflow velocity that could inform simulations of AGN feedback.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 1-dex result assumes the integrated r=1.5'' [Ar IV] density represents outflowing gas across tens of kpc; since the ratio is emission-weighted and the paper says it traces the galaxy overall, this is the weakest load-bearing point.","rationale":"The reader's weakest assumption is exactly the one I find most load-bearing: a single integrated-aperture [Ar IV] density is applied to all outflow spaxels, while the paper itself cautions that this density traces the galaxy overall rather than the outflow specifically. My concern adds a sharper technical point: the [Ar IV] doublet ratio is emission-measure-weighted, so in a multiphase medium it will be biased high relative to the mass-weighted density required by Eq. 1. Because all derived masses, rates, and powers scale as n_e^-1, this bias directly threatens the headline '1 dex lower' result. The paper is otherwise careful: the arithmetic is transparent, the [Ar IV] measurements are cross-checked against [Ne V] results in overlapping sources, and the authors explicitly flag the density caveat in Sec. 5.3. The correlation claim in Fig. 7 has additional weaknesses (post-hoc velocity cut, non-independent global/peak points, no fit uncertainties), but the density representativeness issue is more fundamental because it also controls the energetics. Since the reader already identified this as the weakest assumption and chose CONDITIONAL, my stress-test does not change the verdict; it does, however, sharpen the condition: the authors should demonstrate, rather than assume, that the measured density is representative of the outflowing gas.","tokens_in":32275,"tokens_out":4909,"duration_ms":50590,"concrete_test":"Stack the spectra of outflow-dominated spaxels (W80 above the Table 2 threshold, excluding the central r=1.5'' aperture) and measure the [Ar IV] 4711/4740 ratio; where S/N is insufficient, use [S II] 6717/6731 on the same stacks. Insert the resulting outflow-specific n_e into Eq. 1 and recompute Table 3. If the stacked density is about 1000 cm^-3 or below for a majority of the sample, the claimed ~1 dex reduction is an artifact of applying a nuclear, emission-weighted density to diffuse outflow gas; if it remains above 2000 cm^-3, the conclusion is robust.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline quantitative claim—that [Ar IV]-based n_e lowers outflow rates and kinetic powers by ~1 dex relative to n_e=1000 cm^-3—is only as good as the assumption that the density measured in the r=1.5'' integrated spectrum is the density of the W80-selected outflow gas. Section 4.2 states the lines are 'only observable in the integrated r=1.5'' spectra' and 'the derived electron density is likely tracing the overall density of the galaxy, not specifically the outflow density.' Eq. 1 scales all masses as n_e^-1, so every outflow spaxel, including those at tens of kpc, inherits this single nuclear value. In addition, a line-ratio density in a multiphase medium is emission-measure-weighted and therefore biased toward the densest clumps, not the mass-weighted density that Eq. 1 requires. The paper itself concedes in Sec. 5.3 that 'high density gas (which we are tracing here) might be a small fraction of the total mass that is outflowing.' If the true outflowing gas is closer to 1000 cm^-3, the [Ar IV]-based rates are underestimates by up to 0.8 dex, and the paper's '1 dex lower' message would be inverted rather than confirmed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents VLT/MUSE integral-field observations of 18 type 2 quasars from the Quasar Feedback Survey (QFeedS), characterising ionised gas kinematics via non-parametric [O III]5007 measurements and identifying outflow-dominated spaxels using a galaxy-dependent W80 threshold. The new quantitative contribution is the measurement of electron densities from the [Ar IV] 4711,4740 doublet in r=1.5 arcsec nuclear spectra, yielding 500–27700 cm^-3 (median ~5900 cm^-3), and the resulting ionised outflow masses, mass-loss rates, and kinetic powers (0.03–4.17 Msun/yr and 10^39.25–10^41.91 erg/s for the global method), which are ~1 dex lower than the commonly assumed n_e=10^3 cm^-3. A second claim is a positive correlation between n_e and V_out for V_out > 630 km/s, interpreted as gas compression at higher outflow velocities. The paper also compares radio and W80 morphologies and places the outflow energetics in the context of literature samples and coupling efficiencies.","tokens_in":32510,"tokens_out":4283,"duration_ms":40607,"significance":"If the density assumption holds, the paper provides a useful demonstration that using a high-density tracer such as [Ar IV] rather than an assumed n_e=10^3 cm^-3 changes ionised outflow energetics by about an order of magnitude, with direct implications for feedback coupling efficiencies and for comparisons with simulations. The homogeneous 18-object sample, consistent non-parametric analysis, and the global/peak decomposition are strengths, as is the careful comparison with earlier work on overlapping objects. However, the statistical support for the headline n_e–V_out correlation is currently weak, and the representativeness of the single nuclear density for the spatially extended outflow is not established, so the quantitative conclusions are conditional on an assumption the paper itself questions.","major_comments":[{"comment":"The single electron density measured from the r=1.5 arcsec integrated nuclear spectrum is used in Eq. 1 for every W80-selected outflow spaxel, including gas at projected radii of tens of kpc (e.g., the 45–50 kpc features in J1356+1026). The paper itself states that the [Ar IV] density is 'likely tracing the overall density of the galaxy, not specifically the outflow density.' Since M_out scales as n_e^-1, the headline result that [Ar IV]-based rates and powers are ~1 dex lower than for n_e=10^3 cm^-3 (Table 3 vs. Table A.1) is contingent on the outflow gas having the same density as the nuclear-averaged gas. If the true outflowing phase is closer to 10^3 cm^-3, the reported values would be underestimated by up to ~0.8 dex and the '1 dex lower' message would be inverted. Please either provide a quantitative justification (e.g., spatially resolved density constraints or photoionisation-model checks) or explicitly frame the numbers as conditional on this assumption.","section":"Sec. 4.2 and Eq. 1"},{"comment":"The claimed positive n_e–V_out correlation is fit only to points with V_out > 630 km/s, a threshold selected after inspecting the same plot, and the fit combines global and peak values from the same galaxies as independent measurements. The text notes that the correlation is weaker when global or peak values are considered alone. To support the 'clear positive correlation' statement, the paper should report the correlation coefficient and its significance (e.g., Spearman rank with p-value) for the combined sample and for global-only and peak-only subsamples, and should account for the non-independence of the two points per galaxy (e.g., by bootstrapping over galaxies or fitting one point per galaxy). As presented, the correlation is not statistically established.","section":"Sec. 5.4 and Fig. 7"},{"comment":"The most extreme outflow-velocity galaxy, J1347+1217 (V_out ~2466 km/s global), is excluded from the n_e–V_out fit because its [Ar IV] doublet is not detected at 3 sigma. If its non-detection is not simply a signal-to-noise effect, the missing high-velocity point could weaken or alter the fitted correlation. Please discuss how this exclusion affects the correlation and whether the non-detection is consistent with the fitted trend.","section":"Sec. 5.4 and Table 2"}],"minor_comments":[{"comment":"The sentence 'the amplitude of unnecessary components components tends to zero' contains a duplicated word.","section":"Sec. 3.1"},{"comment":"The text 'clouds of cool emitting gas are are eventually destroyed' has a duplicated 'are'.","section":"Sec. 5.2"},{"comment":"The captions use 'colobars' instead of 'colorbars'.","section":"Fig. A.2–A.16 captions"},{"comment":"The text says the grey region contains points with V_out < 630 km/s, while the body text says 'We exclude the 4 galaxies at the lower V_out end'; please make the description of the threshold and the excluded points consistent.","section":"Sec. 5.4 / Fig. 7 caption"},{"comment":"The entries for J1114+1939 and J1347+1217 would be clearer with explicit 'not detected' notes rather than a dash, especially since the notes below the table already explain the reason.","section":"Table 2"}],"recommendation":"major_revision","confidential_remarks":"The paper's own caveats in Sec. 4.2 and Sec. 5.3 substantially weaken the headline '1 dex lower' claim, and the n_e–V_out correlation is presented more strongly in the abstract and conclusions than the body supports. The revisions needed are within scope: reframing the density claim as conditional and adding a proper statistical treatment of the correlation. The manuscript is otherwise a solid observational contribution appropriate for A&A."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the 1-dex reduction is real arithmetic and the paper deserves a serious referee; the n_e-V_out correlation is a promising hint, not a locked result.\n\nWhat's new: first systematic [Ar IV] densities for this sample of type 2 QSOs, with a clean MUSE/VLA dataset and non-parametric W80 kinematics. The finding that densities (median ~5900 cm^-3) push outflow mass rates and kinetic powers down by ~1 dex relative to the common n_e=1000 assumption is well supported by the tables and follows directly from Eq. 1. The paper is honest about the caveats: densities come from integrated r=1.5'' spectra, no outflow-specific component, and they even flag that high-density gas may be a small fraction of the outflowing mass. That is more than many papers do.\n\nWhere it gets shaky: the correlation in Fig. 7. It is fitted after cutting V_out < 630 km/s, combining global and peak measurements from the same galaxies as if independent, and no fit uncertainties are reported. The abstract says 600, the text says 630. With 15 galaxies and two methods per galaxy, the effective number of independent points is small. The correlation may well be real, but the case is not yet made; it needs a proper covariance treatment or a split-sample check. This is the main thing I would ask for in revision.\n\nThe bigger conceptual caveat is the one the stress-test note raises: an emission-line ratio in a multiphase medium is emission-measure-weighted, biased to dense clumps, while Eq. 1 needs a mass-weighted density. If the bulk of the outflowing gas is at lower density, the measured rates are not 1 dex too high relative to the truth, they are 1 dex too low relative to what n_e=1000 would give. The paper's wording is careful, but the '1 dex lower' message is easier to over-read than the caveat.\n\nBottom line: this is a solid, useful paper for anyone working on ionised outflows in AGN. The density set is a real addition, and the comparison with previous samples is valuable. The correlation claim needs more work, but the paper's central result holds as far as the measurements reach. I'd send it to a good referee, ask them to focus on the correlation and the density-representativeness caveat, and cite it myself for the density values. Yes, I'd bring it to reading group, mostly for the discussion of density tracer systematics.","headline":"The 1-dex density correction is solid arithmetic; the n_e-V_out correlation is a promising hint that needs a proper fit.","tokens_in":33201,"tokens_out":2471,"would_cite":true,"duration_ms":23070,"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":"Measuring electron densities from the [Ar IV] doublet in 18 type 2 quasars lowers ionised outflow mass rates and kinetic powers by about an order of magnitude and reveals a positive density–velocity correlation for fast outflows.","keywords":["AGN feedback","ionised outflows","electron density","type 2 quasars","[Ar IV] doublet","integral field spectroscopy","kinetic coupling efficiency","quasar-driven winds"],"falsifier":"A spatially resolved density measurement in outflow spaxels of one of these quasars (for example, with deeper MUSE or JWST/NIRSpec observations that detect [Ar IV] or transauroral [O II]/[S II] lines outside the nuclear aperture) would settle the claim: if densities at kiloparsec radii approach $10^3\\,\\mathrm{cm}^{-3}$, the roughly one-dex reduction and the $n_e$–$V_{\\rm out}$ correlation would not survive.","tokens_in":32034,"feed_emoji":"🔭","tokens_out":15339,"duration_ms":118932,"temperature":0.7,"pith_summary":"This paper studies ionised gas outflows in 18 luminous, obscured quasars (type 2 QSOs) at redshifts 0.085–0.2 using MUSE integral-field spectroscopy. Its central quantitative claim is that when electron densities are measured from the [Ar IV] 4711,4740 doublet rather than assumed to be 1000 $\\mathrm{cm}^{-3}$, the derived outflow mass rates (0.03–4.17 $M_\\odot\\,\\mathrm{yr}^{-1}$) and kinetic powers ($10^{39.25}$–$10^{41.91}$ erg s$^{-1}$) fall by roughly one order of magnitude. The paper also reports a positive correlation between electron density and outflow velocity for velocities above about 630 km s$^{-1}$, which it interprets as compression of the ionised gas clouds by faster outflows. A sympathetic reader would care because outflow mass and energy are the quantities used to judge whether AGN feedback can affect galaxy evolution, and the electron density has been one of the largest uncertainties in those estimates.","feed_headline":"Measured gas densities cut quasar outflow rates tenfold","feed_subtitle":"With densities from the [Ar IV] doublet, outflow powers fall about 1 dex and density tracks velocity.","key_machinery":"The load-bearing measurement is the [Ar IV] 4711,4740 Å doublet ratio, an electron-density diagnostic for fully ionised gas that is sensitive from roughly $10^2$ to $6\\times10^4\\,\\mathrm{cm}^{-3}$ at an assumed electron temperature of 15,000 K; the ratio is converted to $n_e$ with the PyNeb routine. Because the doublet is detected only in integrated $r=1.5$ arcsecond nuclear spectra, the single resulting density is inserted into the [O III] mass formula $M = 8\\times10^7 M_\\odot\\,(L_{\\rm [OIII]}/10^{44}\\,\\mathrm{erg\\,s}^{-1})(n_e/500\\,\\mathrm{cm}^{-3})^{-1}$ for every outflow spaxel. Outflow spaxels are selected with a per-galaxy $W_{80}$ threshold set from the unresolved nuclear velocity dispersion, and outflow velocities and radii are computed with flux-weighted global and radial peak methods. This chain carries the argument from line ratios to outflow energetics.","core_discovery":"The paper's central discovery is that [Ar IV]-based electron densities in the nuclear regions of these quasars are high (500 to 27,700 $\\mathrm{cm}^{-3}$, median about 5,900 $\\mathrm{cm}^{-3}$), and that adopting these densities instead of the customary $10^3\\,\\mathrm{cm}^{-3}$ lowers outflow masses, mass-loss rates, and kinetic powers by about one dex: global mass outflow rates fall to $0.03$–$4.17\\,M_\\odot\\,\\mathrm{yr}^{-1}$ and kinetic powers to $10^{39.25}$–$10^{41.91}$ erg s$^{-1}$. It further identifies a positive correlation between $n_e$ and $V_{\\rm out}$ for $V_{\\rm out}>630$ km s$^{-1}$, suggesting that faster gas is denser, consistent with gas compression in quasar-driven outflows. The paper also reports that the radio emission peaks coincide with the largest [O III] line widths ($W_{80}$) in most galaxies, and that mass outflow rates peak in the innermost $\\sim$2 kpc and decline outward, consistent with cloud destruction or incomplete ionisation at large radii.","pith_inferences":["If the $n_e$–$V_{\\rm out}$ correlation holds beyond this sample, unresolved single-fibre surveys could use $W_{80}$ to statistically correct outflow densities and energetics in large quasar samples where only integrated spectra exist.","The roughly one-dex drop in kinetic power brings observed outflow coupling closer to the low coupling efficiencies used in some galaxy-formation models, although the comparison is indirect because only the kinetic energy in one gas phase is measured.","A natural test is to look for the same density–velocity trend in molecular outflows of these galaxies: if the dense optical gas is a small fraction of the outflowing mass, the mass rates could be higher while the compression link between speed and density remains.","Re-observing a few of these quasars with higher spatial resolution and longer integrations could map $n_e$ directly inside the outflow cones, separating nuclear-averaged density from true outflow density — the assumption on which the current result rests."],"forward_implications":["Outflow masses and kinetic powers published for type 2 quasars with an assumed $n_e=1000\\,\\mathrm{cm}^{-3}$ are likely overestimated by about an order of magnitude.","The kinetic coupling efficiencies of these outflows ($\\log \\varepsilon_f \\approx -6.8$ to $-4.2$ of $L_{\\rm bol}$) are low, and are lower limits because only one gas phase at optical wavelengths and no turbulent term are included.","The positive $n_e$–$V_{\\rm out}$ correlation at $V_{\\rm out}>630$ km s$^{-1}$, if confirmed, implies that faster outflows are denser, favouring models in which the outflow compresses the ionised clouds.","The decreasing radial mass-outflow profiles with a peak inside 2 kpc indicate that most of the outflowing ionised mass is released in the nuclear region, with material farther out either destroyed or not excited by the AGN.","Radio peaks spatially matching the $W_{80}$ peaks support a jet or radio-driven disturbance of the ionised gas even in radio-quiet quasars."],"supporting_citations":[{"why":"Supplies the [O III]-based ionised gas mass formula adopted to convert luminosity and density into outflow mass.","marker":"Carniani et al. 2015"},{"why":"Provides the non-parametric $W_{80}$ and velocity-percentile technique used for all kinematic measurements and the $W_{80}=2.563\\sigma$ conversion.","marker":"Liu et al. 2013"},{"why":"Establishes the [Ar IV] 4711,4740 doublet as an electron-density tracer for AGN narrow-line regions with its probed density range.","marker":"Binette et al. 2024"},{"why":"Provides the PyNeb routine used to convert the [Ar IV] flux ratio into electron density.","marker":"Luridiana et al. 2015"},{"why":"Supplies the per-galaxy $W_{80}$ outflow threshold method and the global/peak outflow property definitions used here.","marker":"Riffel et al. 2023"},{"why":"Documents the limited density range and biases of [S II]-based densities, motivating the choice of [Ar IV].","marker":"Holden et al. 2026"},{"why":"Simulations used to interpret the $n_e$–$V_{\\rm out}$ correlation as compression and to predict density–luminosity relations.","marker":"Almeida et al. 2026"},{"why":"Simulations showing that high-density gas may carry a small fraction of outflow mass, bounding the energetics interpretation.","marker":"Ward et al. 2024"}],"fun_headline_variants":["Quasar outflow rates drop tenfold with measured gas densities","Electron density data cut quasar outflow power tenfold","Quasar outflows: faster gas is denser, power drops 10x","Quasar winds: density tracks velocity, power cut tenfold","Denser gas in quasar outflows cuts energy by 10x"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The single electron density measured from the integrated nuclear spectrum is assumed to apply to all gas classified as outflowing, even tens of kiloparsecs away; if those clouds are less dense, the derived masses, rates, and the reported density–velocity correlation would all be biased.","fun_headline_variants_meta":{"raw":{"variants":["Quasar outflow rates drop tenfold with measured gas densities","Electron density data cut quasar outflow power tenfold","Quasar outflows: faster gas is denser, power drops 10x","Quasar winds: density tracks velocity, power cut tenfold","Denser gas in quasar outflows cuts energy by 10x"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001168,"raw_usage":{"total_tokens":4975,"prompt_tokens":1230,"completion_tokens":3745,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":846,"completion_tokens_details":{"reasoning_tokens":3655}},"tokens_in":846,"tokens_out":3745,"duration_ms":22895,"temperature":1.0,"reasoning_tokens":3655,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T00:11:33.860554+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A spatially resolved density measurement in outflow spaxels of one of these quasars (for example, with deeper MUSE or JWST/NIRSpec observations that detect [Ar IV] or transauroral [O II]/[S II] lines outside the nuclear aperture) would settle the claim: if densities at kiloparsec radii approach $10^3\\,\\mathrm{cm}^{-3}$, the roughly one-dex reduction and the $n_e$–$V_{\\rm out}$ correlation would not survive.","supporting_citations":[],"review_version":1}