{"id":"b58ef89e-8e05-4808-b905-baeca5eac888","arxiv_id":"1908.03104","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"HST/STIS resolves the sub-kiloparsec warm outflow in ULIRG F14394+5332E, revealing a steep velocity gradient and large turbulent line widths consistent with an AGN-driven wind shredding dense clouds.","lead":"Using Hubble's STIS spectrograph, astronomers resolved the compact, fast-moving gas outflow in the merging galaxy F14394+5332E for the first time, showing gas accelerating from about 500 to 1800 km/s across a region less than 1 kpc across. The result tests how much energy supermassive black holes deposit into their host galaxies during major mergers.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The headline radii and velocity gradient are measured from an assumed, unquantified AGN position hidden in a dust lane; a plausible 0.2 arcsec offset could move r_max out of the sub-kpc regime.","rationale":"The reader correctly identified the assumed AGN position as the weakest assumption. I agree that this is the single most load-bearing choice in the paper: every radial extent, the velocity gradient in Figure 3, and the geometric discussion in Section 4.1 are defined relative to the cross in Figure 1. The authors are transparent about the uncertainty, and the qualitative picture—a bright, broad-line, blueshifted region to the north embedded in a fainter, narrower, less ionized halo—is likely robust to small shifts. However, the abstract's quantitative claims ('r_max ~ 0.9 kpc', '~500 to 1800 km/s', 'sub-kpc') could change by tens of percent for a 0.1–0.2 arcsec offset, and the paper does not provide a sensitivity analysis or uncertainty on the nuclear position. This does not invalidate the work, but it means acceptance should be conditional on adding such an analysis or on obtaining an independent AGN location from mid-IR or radio observations. The data reduction, fitting procedures, and modeling are otherwise careful and internally consistent, and the paper's self-flagged limitations are appropriate. Therefore I recommend CONDITIONAL rather than outright rejection: the central claim is plausible and well-supported, but its headline numbers need a robustness check against the assumed origin.","tokens_in":117,"tokens_out":9330,"duration_ms":223405,"concrete_test":"Using the existing STIS data and the line-free HST continuum image, recompute the spatial offsets and the ΔV versus radius relations for a grid of nuclear positions shifted by ±0.1, ±0.2, and ±0.3 arcsec in both the along-slit and perpendicular directions, placing the nucleus at plausible locations within the dust lane. If for any offset within ±0.2 arcsec the [OIII] r_max becomes >1 kpc or the monotonic inner-to-outer blueshift gradient disappears, then the headline sub-kpc and gradient claims need to be revised or qualified; if all offsets preserve the qualitative result within uncertainties, the central claim is robust.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central result—a compact high-ionization outflow with r_max ~ 0.9 kpc and a velocity gradient increasing from ~500 to 1800 km/s—is derived from distances measured from the cross in Figure 1, which the authors adopt as the AGN position. Section 2 explicitly states: 'lacking the high resolution mid-IR imaging that might locate the nucleus more accurately, it is clear that there remains some uncertainty about the true nuclear position.' Because the AGN is hidden in a dust lane and the outflow is detected only to the north, a shift of the assumed nucleus by 0.1–0.2 arcsec (0.19–0.37 kpc) along the slit changes every radial distance used in Figures 3 and 4. For example, a 0.2 arcsec southward shift changes the [OIII] detection limit from 0.87 kpc to ~1.24 kpc, no longer sub-kpc, while a northward shift steepens the gradient; a perpendicular shift changes which rows are assigned to which radius. The monotonicity of the velocity gradient may survive moderate offsets, but its slope and the sub-kpc characterization are not robust until this sensitivity is quantified. This is the most load-bearing assumption because it enters the abstract's headline numbers and the geometric interpretation in Section 4.1.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents HST/STIS long-slit spectroscopy of the ULIRG F14394+5332E, resolving the warm AGN-driven outflow at sub-kpc scales for the first time in a ULIRG. The observations show a compact, high-ionization [OIII] outflow region (r_max ~ 0.9 kpc) embedded in a more extended Halpha-emitting halo (r_max ~ 1.4 kpc) that is kinematically quiescent. The [OIII] lines are broad and blueshifted, with a steep velocity gradient in which the blueshift increases from ~750 to ~1800 km/s along the central slit. The paper interprets the kinematics as evidence for local acceleration and hydrodynamic destruction of clouds in a hot AGN wind, and uses the resolved measurements to estimate mass outflow rates and kinetic powers, finding values within a factor ~2-3 of earlier spatially integrated estimates.","tokens_in":14185,"tokens_out":9252,"duration_ms":99542,"significance":"The result is significant because it provides the first spatially resolved warm outflow in a ULIRG, directly constraining the radial extent, geometry, and velocity structure that earlier ground-based studies could only model. The new STIS data are independent of the model assumptions, and the analysis is careful: standard CALSTIS reduction, explicit treatment of doublets, an instrumental-width correction estimated from narrow-line regions, and candid discussion of Halpha+[NII] fitting degeneracies. The paper also makes falsifiable predictions for future molecular and X-ray observations. If the nuclear-position issue raised below is addressed, the conclusions about compact outflows and modest energetics would support the emerging picture that warm AGN outflows in local ULIRGs are compact and energetically modest.","major_comments":[{"comment":"The zero point for all radial distances and for the velocity gradient is the assumed AGN position at the cross in Figure 1. The text acknowledges that this position is uncertain, but the uncertainty is never quantified. A shift of 0.1-0.2 arcsec (0.19-0.37 kpc) along the slit, which is comparable to the 0.1 arcsec slit separation and several 0.05078 arcsec pixels, changes the [OIII] detection limit from 0.88 kpc to 1.06-1.24 kpc if the true nucleus lies further south, directly affecting the abstract's 'sub-kpc' characterization; it also changes the slope and amplitude of the velocity gradient shown in Figure 3 and the geometric interpretation in Section 4.1. Please add a quantitative sensitivity analysis (for example, shift the zero point by +/-0.1 and +/-0.2 arcsec and recompute the r_max values and gradient fits), and either propagate the resulting systematic uncertainties into the headline claims or temper the 'sub-kpc' and specific radius statements accordingly.","section":"2 (Observations and Reductions), final paragraph"},{"comment":"The open-symbol estimates in Figures 7 and 8 are described as upper limits, but they assume v_out = 1800 km/s and FWHM = 1200 km/s for all spatial regions, even though the measured [OIII] FWHM exceeds 2000 km/s in the compact NE region of the E slit (Section 3.1) and the true velocity could exceed the maximum projected value if projection effects are important. These estimates should be described as 'maximal estimates under a stated set of assumptions' rather than strict upper limits, or the assumptions should be tightened so that the label 'upper limit' is valid. This affects the strength of the conclusion that the kinetic powers remain modest (<0.5% L_bol).","section":"4.2, Figures 7 and 8"}],"minor_comments":[{"comment":"The abstract states '600 < FWHM < 1500 km/s' throughout the outflow region, but the E-slit [OIII] measurement in Section 3.1 is FWHM > 2000 km/s (with large uncertainty). Please qualify the abstract or explicitly exclude that low-S/N region from the quoted range.","section":"Abstract and Section 3.1"},{"comment":"The sentence describing the spatial registration ('fitting Gaussians to the peak emission in spatial continuum slices') is unclear; please specify how the continuum centroid was measured as a function of wavelength and how the offsets to the assumed nucleus position were determined.","section":"Section 2"},{"comment":"The name 'Baldwin, Philips & Terlevich' should be 'Baldwin, Phillips & Terlevich'.","section":"Section 3.2"},{"comment":"The in-text citation 'Mellema et al. 2003' is listed as 'Mellema, Kurk & Rottgering 2002' in the reference list, and 'Revalaski et al. 2018' should be 'Revalski et al. 2018'. Please verify these citations.","section":"References"},{"comment":"No error bars are visible on the plotted kinematic measurements; if uncertainties are not shown, please state this in the captions and give representative uncertainties in the text.","section":"Figures 3 and 4"},{"comment":"The sentence 'and the Delta-V = -990 km/s assumed by Spence et al.' appears to be missing a comparison word; presumably the broad Halpha velocity shifts are lower in amplitude than -990 km/s.","section":"Section 4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is fundamentally sound and the new data are valuable. The main substantive obstacle is the unquantified assumed AGN position, which enters the headline sub-kpc radius and the velocity gradient. This is addressable with a sensitivity analysis and adjusted wording, so I recommend major revision rather than rejection. The energetics 'upper limit' wording should also be corrected. If the authors provide the requested sensitivity analysis, I would support publication."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Solid paper, worth a serious referee. It gives the first spatially resolved STIS spectroscopy of a warm AGN-driven outflow in a ULIRG, and the new numbers — line widths 600–1500 km/s, blueshifts up to 1800 km/s, steep radial gradient — are real measurements, not derived from models. The data reduction is transparent, the Gaussian fitting is constrained by physics, and the authors flag their main degeneracies (Hα+[NII] blending, assumed instrumental width). The comparison with their earlier ground-based and imaging work is honest and the mass/energy estimates are presented as a range with clear assumptions. That is the good part.\n\nThe soft spot is the nuclear position. The AGN is hidden in a dust lane and the authors simply adopt the cross in Figure 1. They admit the uncertainty but never quantify it. That matters because every radius and the velocity gradient are measured from that point. A 0.1–0.2 arcsec offset along the slit (0.19–0.37 kpc) would change the [OIII] detection limit from ~0.9 kpc to ~1.2 kpc — still compact, but no longer 'sub-kpc' in the headline sense. A perpendicular shift changes which rows fall at which radius. The velocity gradient probably survives, but its slope and the compactness claim are not fully robust until this sensitivity is tested. This is a minor-to-moderate issue, not a fatal one: the qualitative picture of a compact, turbulent, AGN-driven outflow is well supported, and the paper already contains most of what you'd need to do the check.\n\nOther assumptions — uniform extinction E(B-V)=0.63, uniform n_e=3500 cm^-3 — are reasonable and explicitly tied to prior work; the authors know they are simplifications. The main conclusions about mass outflow rates and kinetic powers are appropriately hedged, and they conclude that the spatial resolution changes the integrated estimates by only a factor ~2–3. That is a useful, quantitative statement.\n\nWho should read it: anyone working on AGN feedback in ULIRGs or on the geometry of warm outflows. It deserves peer review. My recommendation: accept after a minor revision that quantifies the effect of the assumed nuclear position on the radial extents and the velocity gradient — a simple shift test would settle it.","headline":"First resolved sub-kpc warm outflow in a ULIRG, but the assumed AGN position is the one real soft spot — a 0.2 arcsec shift could push the headline radius over 1 kpc; still a strong, careful paper.","tokens_in":14727,"tokens_out":2385,"would_cite":true,"duration_ms":25908,"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":"HST/STIS spectroscopy resolves the sub-kpc warm outflow in the ULIRG F14394+5332E, revealing a compact, turbulent flow accelerated by the AGN wind.","keywords":["ultraluminous infrared galaxies","AGN-driven outflows","warm ionized gas","emission-line kinematics","galaxy mergers","active galactic nuclei","HST/STIS spectroscopy","F14394+5332"],"falsifier":"High-resolution mid-infrared imaging that pins down the true AGN position would test the central claim: if the nucleus is significantly offset from the assumed cross, the measured radial extents and the velocity gradient change. Alternatively, a molecular-gas map or X-ray observation could test the wind-shock interpretation, since the paper predicts any molecular outflow has an extent similar to or smaller than the warm outflow and predicts X-ray emission from hot post-shock gas.","tokens_in":13755,"feed_emoji":"🔭","tokens_out":6484,"duration_ms":60992,"temperature":0.7,"pith_summary":"This paper uses HST/STIS long-slit spectroscopy to resolve, for the first time, the warm ionized outflow in a ULIRG (the merging galaxy F14394+5332) at sub-kiloparsec scales. It finds a compact, high-ionization outflow region roughly 0.9 kpc across, with line widths between 600 and 1500 km/s and blueshifts that grow from about 500 km/s near the nucleus to 1800 km/s at the outer edge. The outflow sits inside a more extended, kinematically quiet halo of lower-ionization gas out to about 1.4 kpc. These observations matter because ground-based studies cannot resolve such outflows, leaving their radial extents, geometries, and true energetics uncertain; the resolved picture supports an AGN-driven wind that accelerates and shreds dense clouds rather than a galaxy-wide outflow.","feed_headline":"First sub-kpc warm outflow resolved in a ULIRG","feed_subtitle":"HST/STIS maps the outflow's speed and size, sharpening estimates of how much AGN power feeds the gas.","key_machinery":"The key object is the spatially resolved emission-line map built from three parallel STIS slits, which lets the authors measure [OIII] and H-alpha+[NII] velocity shifts, line widths, and ionization ratios as functions of projected radius. The load-bearing diagnostic is the steep outward increase in blueshift along with the uniformly large line widths, which distinguishes between true radial acceleration, projection geometry, and local turbulent acceleration. The paper also uses BPT line-ratio diagrams to separate the high-ionization outflow from the lower-ionization quiescent halo, and a mass-conservation equation for the resolved outflow to convert H-alpha luminosities, assumed density, and velocity into local mass outflow rates and kinetic powers.","core_discovery":"The central claim is that the warm, AGN-driven outflow in F14394+5332E is compact and turbulent rather than galaxy-scale. At the assumed nucleus, the [OIII] emission is broad and strongly blueshifted throughout a region extending roughly 0.9 kpc to the north, with the projected blueshift increasing steeply outward from about 500 to 1800 km/s while the line width stays high (600–1500 km/s); a fainter H-alpha-emitting halo reaches about 1.4 kpc and is kinematically quiescent with lower ionization. The paper interprets the kinematics as local acceleration of dense clouds by the forward shock of a hot, low-density AGN wind, with the large widths signalling hydrodynamic destruction. A corollary is that the sharp edge of the disturbed region marks the current extent of the shock, so the outflow has not yet engulfed the host galaxy's warm and cool gas.","pith_inferences":["If the projection interpretation is right, studies that use lower mean velocities from unresolved spectra will systematically underestimate the warm outflow's kinetic power, while the compact size limits how much of the galaxy the outflow directly affects.","The same slit-mapping method applied to more ULIRGs could turn compactness from a single-object finding into a census, and the predicted similarity between warm and molecular outflow extents is directly testable with CO mapping.","The apparent 'hole' in narrow H-alpha emission near the nucleus suggests a region where all warm gas has been accelerated; comparing such holes across merger stages may trace how the outflow progressively sweeps up the ISM."],"forward_implications":["Warm AGN-driven outflows in ULIRGs can be compact (sub-kpc) even in rapidly evolving merger remnants, so they do not necessarily fill the host galaxy's warm and cool gas.","Mass outflow rates from resolved regions fall around 0.3–0.7 solar masses per year, with upper-limit estimates around 1–3 solar masses per year, only a factor of 2–3 different from spatially integrated values.","Kinetic coupling efficiencies lie between about 0.0025% and 0.5% of the AGN bolometric luminosity, leaving little room for the warm phase to be the dominant feedback channel.","The large line widths throughout the outflow favour local acceleration and destruction of clouds by the hot wind over a simple radially accelerating outflow.","If the outer blueshifts reflect projection geometry, the 1800 km/s velocity measured at the outflow edge is the closest available estimate of the true deprojected outflow speed."],"supporting_citations":[{"why":"Previous ground-based spectroscopy and trans-auroral density measurements that set the galaxy redshift, electron density, and spatially integrated outflow estimates the STIS results are compared against.","marker":"Spence et al. 2018"},{"why":"HST narrow-band [OIII] imaging that identified the compact extended emission and selected the STIS slit positions.","marker":"Tadhunter et al. 2018"},{"why":"Paper I of the series, providing the sample context and the electron-density method used to convert fluxes to masses.","marker":"Rose et al. 2018"},{"why":"Mapping of the two nuclei and the merger structure used to associate the AGN with F14394+5332E.","marker":"Kim, Veilleux & Sanders 2002"},{"why":"Theoretical wind-driven outflow model with a hot shock that the paper invokes to explain cloud acceleration and destruction.","marker":"King & Pounds 2015"},{"why":"Cloud-crushing and Kelvin-Helmholtz instability calculations used to argue that clouds are destroyed while being accelerated.","marker":"Klein, McKee & Colella 1994"},{"why":"Comparable STIS observations of nearby quasar outflows that support the compact-outflow picture.","marker":"Fischer et al. 2018"},{"why":"Another STIS outflow study used as a consistency check on resolved outflow kinematics in AGN.","marker":"Revalski et al. 2018"}],"fun_headline_variants":["Resolved AGN outflow in ULIRG shows clouds at 1800 km/s","Hubble resolves sub-kpc outflow, sees 1800 km/s blueshift","Compact AGN outflow mapped: shock accelerates, destroys clouds","First sub-kpc warm outflow imaged in ULIRG's turbulent core"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The paper assumes the AGN lies at a specific point marked on its image, and all radial distances, the outward velocity gradient, and the inferred outflow geometry would shift if the true nucleus were anywhere else.","fun_headline_variants_meta":{"raw":{"variants":["Resolved AGN outflow in ULIRG shows clouds at 1800 km/s","Hubble resolves sub-kpc outflow, sees 1800 km/s blueshift","Compact AGN outflow mapped: shock accelerates, destroys clouds","First sub-kpc warm outflow imaged in ULIRG's turbulent core"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000645,"raw_usage":{"total_tokens":2997,"prompt_tokens":1011,"completion_tokens":1986,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":627,"completion_tokens_details":{"reasoning_tokens":1905}},"tokens_in":627,"tokens_out":1986,"duration_ms":18032,"temperature":1.0,"reasoning_tokens":1905,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T14:23:32.524273+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"High-resolution mid-infrared imaging that pins down the true AGN position would test the central claim: if the nucleus is significantly offset from the assumed cross, the measured radial extents and the velocity gradient change. Alternatively, a molecular-gas map or X-ray observation could test the wind-shock interpretation, since the paper predicts any molecular outflow has an extent similar to or smaller than the warm outflow and predicts X-ray emission from hot post-shock gas.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Previous ground-based spectroscopy and trans-auroral density measurements that set the galaxy redshift, electron density, and spatially integrated outflow estimates the STIS results are compared against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"HST narrow-band [OIII] imaging that identified the compact extended emission and selected the STIS slit positions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Paper I of the series, providing the sample context and the electron-density method used to convert fluxes to masses."},{"cited_title":"C., Veilleux, S., Sanders D","cited_arxiv_id":null,"evidence_quote":"Mapping of the two nuclei and the merger structure used to associate the AGN with F14394+5332E."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Theoretical wind-driven outflow model with a hot shock that the paper invokes to explain cloud acceleration and destruction."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Cloud-crushing and Kelvin-Helmholtz instability calculations used to argue that clouds are destroyed while being accelerated."},{"cited_title":"et al., 2018, ApJ, 856, 102","cited_arxiv_id":null,"evidence_quote":"Comparable STIS observations of nearby quasar outflows that support the compact-outflow picture."},{"cited_title":"et al., 2018, ApJ, 867, 88","cited_arxiv_id":null,"evidence_quote":"Another STIS outflow study used as a consistency check on resolved outflow kinematics in AGN."}],"review_version":1}