{"id":"6c2ff777-4659-485c-9a45-9a3a475d5e85","arxiv_id":"2506.14931","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Spatially resolved BPT maps of NGC 5005 reveal a Seyfert-like nucleus, a thin LINER-like cocoon, and jet-associated star-forming regions.","lead":"Using Hubble Space Telescope narrowband images, astronomers classified the ionized gas in the galaxy NGC 5005 by its emission-line ratios and found a compact Seyfert-like nucleus surrounded by a thin LINER-like cocoon. The maps also show star-forming regions near a radio jet, raising the possibility that the jet triggers star formation in this low-luminosity active galaxy.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The thin LINER cocoon and the compact Seyfert 'cones' rest on a post-hoc [OIII]/Hbeta split applied to a handful of pixels with no error propagation; a forward-model test can determine whether they are PSF/threshold artifacts.","rationale":"The reader's weakest-assumption analysis identifies exactly the load-bearing issue: the thin cocoon and the Seyfert 'cones' may be artifacts of an eye-chosen threshold and PSF smearing. My independent reading agrees, with the additional emphasis that no uncertainty propagation is performed on the line ratios, making it impossible to know whether the handful of Seyfert and cocoon pixels are statistically distinct from the surrounding LINER population. The paper's own text acknowledges reliance on a forthcoming Chandra study for AGN confirmation and explicitly calls for high-resolution IFU kinematics to disentangle ionization mechanisms, which supports a conditional verdict rather than a rejection. The large-scale results (the r~4 kpc HII ring, the central r~1 kpc LINER region, and the extended LINER zone) are based on many more pixels and are less vulnerable to PSF/threshold effects, so they should stand. The proposed forward-model test is straightforward with public HST data and would settle whether the cocoon is a real physical structure or a numerical artifact. I therefore retain the reader's CONDITIONAL verdict, pending this targeted check.","tokens_in":17351,"tokens_out":4324,"duration_ms":45372,"concrete_test":"Forward-model the nuclear line emission: place an unresolved Seyfert-like source at the optical centroid with [OIII]/Hbeta and [SII]/Halpha ratios matching the observed Seyfert pixels, embed it in a smooth extended LINER-like background with the median green-LINER ratios, convolve all line images with the WFC3 PSF, add Poisson noise matching the exposure times, apply the same 3-sigma mask and the same log([OIII]/Hbeta) >= -0.28 split, and compare the resulting 'cocoon' pixels with Map 3 of Figure 4. If the model reproduces a ~10-20 pc blue annulus around a 7-pixel Seyfert cone, the cocoon is a PSF/threshold artifact rather than a physical shell; if it does not, the cocoon survives this specific challenge.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's most distinctive claims are the ~20 pc LINER cocoon and the Seyfert-like biconical nucleus (Section 4.1.1). Both depend on per-pixel S-BPT classifications made from narrowband images with no propagated uncertainties. In Section 3.2, LINER pixels are split into 'blue' and 'green' groups using a fixed threshold, log([OIII]/Hbeta) >= -0.28, chosen because a diagonal branch is visible in the S-BPT diagram. This threshold is post-hoc and has no physical or statistical justification. The cocoon and 'cones' are tiny: the primary Seyfert 'cone' is 7 pixels and the secondary is 3 pixels; cocoon-like pixels are 0.54% and 2.69% of Maps 2 and 3. With WFC3 PSF FWHM ~0.1\" (~10 pc) and pixel scale 0.04\" (~4 pc), a claimed cocoon thickness of 10-20 pc is only 1-2 PSF widths. Because the nucleus is bright in [OIII], PSF smearing plus a fixed [OIII]/Hbeta threshold can produce an annular 'cocoon' at the radius where the PSF wings cross the threshold even if the underlying gas is a smoothly stratified LINER region. The Seyfert-like 'biconical' morphology is therefore not established, and the positive-feedback HII regions are supported only by spatial coincidence with radio features, not by kinematics or stellar population ages. The large-scale maps and the r~4 kpc HII ring are likely robust, but the nuclear cocoon and jet-triggered star-formation claims are not yet secured.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents HST/WFC3 narrowband imaging of the LINER I galaxy NGC 5005, derives per-pixel S-BPT classifications using [O III]/Hβ versus [S II]/Hα, and constructs spatially resolved excitation maps from 4 pc pixels over an 8 kpc field. It reports a compact Seyfert-like nucleus, a ~20 pc thick LINER cocoon, a central ~1 kpc LINER region, an extended >2 kpc LINER zone, and H II regions including a ~4 kpc star-forming ring. The authors interpret the inner H II regions as possibly triggered by jet–ISM interactions and the cocoon as shock-excited gas.","tokens_in":17700,"tokens_out":5612,"duration_ms":47411,"significance":"If the small-scale claims hold, this is the first spatially resolved BPT study of a LINER-dominated galaxy at ~10 pc resolution and would be an important step in linking low-luminosity AGN to spatially resolved ionization structure. The data reduction is standard and the large-scale maps (Map 1) are a useful resource. However, the central new claims—the Seyfert biconical nucleus, the LINER cocoon, and the positive-feedback interpretation—rest on per-pixel classifications that are sensitive to the assumed Hα fraction and a post-hoc threshold, and on spatial scales comparable to the PSF. The paper provides no propagated uncertainties and no forward-model test, so these claims are not yet secure.","major_comments":[{"comment":"The classification depends on the assumed Hα fraction in the F658N bandpass. The paper adopts 45% but reports that using 30% shifts the fractions substantially (Map 3 Seyfert-like from 0.41% to 0.11%, cocoon-like from 2.69% to 3.92%) and, per §4.1.2, moves the secondary Seyfert-like cone into the cocoon class. Since the Seyfert-like cone and the cocoon are the core new results, the paper must show that the central conclusions are invariant under this systematic, or provide a full uncertainty map.","section":"§2.3, §3.2"},{"comment":"The blue/green LINER split at log([O III]/Hβ) ≥ −0.28 is chosen after visually identifying a diagonal branch in the S-BPT diagram. This makes the cocoon a post-hoc selection rather than a prediction. A statistical or physical justification (e.g., a two-component fit to the LINER distribution, a bimodality test, or an independent diagnostic) is needed before the cocoon can be described as a distinct physical structure.","section":"§3.2, Figure 4"},{"comment":"The claimed cocoon thickness (~10–20 pc) is only 1–2 times the WFC3/UVIS PSF FWHM (~10 pc), and the \"biconical\" Seyfert-like region comprises just 10 pixels. With a bright nuclear [O III] source and a hard threshold in [O III]/Hβ, PSF wings can produce an annular cocoon-like structure even if the underlying excitation profile is smooth. The authors should forward-model the PSF (e.g., smooth a smooth radial profile and apply the same classification pipeline) or show that the cocoon persists with different thresholds and after deconvolution.","section":"§4.1.1"},{"comment":"The \"positive feedback\" interpretation is based only on spatial coincidence between radio jet contours and H II-like pixels; no stellar population ages, extinction-corrected equivalent widths, or kinematic data are used to demonstrate recent star formation triggered by the jet. The abstract and conclusions present this as a finding despite the speculative wording in §4.1.2. Either the authors should add independent evidence or clearly label this as a working hypothesis to be tested with future IFU/age-sensitive observations.","section":"§4.1.2, §5"}],"minor_comments":[{"comment":"The continuum filters F547M and F814W have 695 s exposures, much shorter than the line-filter exposures; the S/N threshold is applied to the line images but not to the continuum-subtracted products. The noise contribution from continuum subtraction should be included in the uncertainty budget.","section":"§2.2"},{"comment":"The pixel fractions quoted in the text and in Figure 4 would be easier to interpret if they included the number of pixels in each class for Maps 1–3, given the large differences in map area.","section":"§3.2"},{"comment":"The large-scale LINER fraction increases from 2.68% to ~15% when the lower Hα limit is adopted; this should be stated alongside the first mention of the large-scale LINER emission to make the systematic uncertainty clear.","section":"§4.2.2"},{"comment":"Equation (A7) defines E(B−V) using a reference color C, but the text does not state how the reference region (the green circle in Figure 1) was selected or how the results depend on that choice. A brief discussion of this sensitivity would help.","section":"Appendix A.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is potentially publishable after a major revision. The large-scale maps and the approach are sound, but the small-scale claims (Seyfert biconical structure, LINER cocoon, positive feedback) need robustness tests against the Hα fraction, the post-hoc threshold, and PSF smearing. I would not reject, because these issues are addressable with additional analysis and a more cautious presentation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth a look, with the cocoon claims taken as provisional. The genuinely new thing here is the first spatially resolved BPT map of a LINER-dominated galaxy: HST narrowband images resolved a compact Seyfert-like core, embedded LINER-like emission, and an outer HII ring out to ~8 kpc. That map is valuable, the data are public, and the methods are described well enough to reproduce. The large-scale picture – central LINER region, r~4 kpc star-forming ring, and overall excitation structure – looks robust.\n\nThe problems start at the inner ~200 pc. The Seyfert-like 'cones' are 7 and 3 pixels. The thin LINER cocoon is defined by splitting LINER pixels at log([OIII]/Hbeta) = -0.28, a threshold chosen after the authors saw a diagonal branch in the S-BPT diagram. With WFC3 PSF FWHM ~10 pc and the claimed cocoon thickness 10–20 pc, PSF smearing plus this fixed threshold can plausibly produce an apparent shell around a bright nucleus even if the underlying excitation is smooth. No uncertainties are propagated into the per-pixel classifications, and the constant 45% Halpha fraction in F658N is a known knob: the authors show a 30% lower limit shifts many pixels between classes. These issues do not kill the paper, but they mean the cocoon and the positive-feedback HII regions are not established. The HII regions are supported only by spatial coincidence with radio features, not by kinematics or stellar ages.\n\nThe authors are honest about these limitations – they flag the shock alternative and the need for high-resolution IFU data. The work is a serious observational study, not a fitted-prediction paper, and the circularity burden is mild. I would send it to referees, and I'd probably cite the resolved map in any LINER/feedback discussion, but I would not yet hang a physical story on the cocoon. It merits follow-up with PSF forward modeling and matched-resolution IFU spectroscopy.","headline":"First resolved S-BPT map of a LINER I nucleus is a solid addition, but the thin cocoon and jet-triggered star formation are threshold- and PSF-sensitive.","tokens_in":18313,"tokens_out":2449,"would_cite":true,"duration_ms":25751,"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":"Using HST narrowband imaging, this paper claims that NGC 5005's nucleus is a faint Seyfert-like AGN wrapped in a thin shock-excited LINER cocoon, with nearby star-forming clumps possibly triggered by the radio jet.","keywords":["LINER galaxies","AGN photoionization","shock excitation","BPT diagram","narrow-line region","jet-ISM interaction","positive feedback","NGC 5005"],"falsifier":"Re-map the central 200 pc with an instrument whose PSF is at least a factor of two smaller than the claimed cocoon thickness (or after PSF deconvolution), and with matched IFU kinematics: if the ring of pixels with $\\log([\\mathrm{O\\,III}]/\\mathrm{H}\\beta) \\ge -0.28$ does not remain as a distinct ~10–20 pc shell, and if no shock-broadened line profiles appear at that radius, the cocoon and its positive-feedback interpretation lose their primary support.","tokens_in":17151,"feed_emoji":"🔭","tokens_out":7119,"duration_ms":65928,"temperature":0.7,"pith_summary":"With ~4 pc pixels from HST narrowband images, this paper maps which ionization mechanism dominates each part of NGC 5005's narrow-line region, from the nucleus out to 8 kpc. It finds a compact Seyfert-like nucleus consistent with a low-luminosity AGN, surrounded by a thin (~20 pc) higher-excitation LINER cocoon that it attributes to shock heating of the interstellar medium. Around the cocoon sit two larger LINER zones, one photoionized by the AGN and one possibly powered by post-AGB stars plus inflow shocks, and an outer H II ring at ~4 kpc. Inside 500 pc, clumpy H II regions overlap the radio jet and a bubble-like Hα feature, suggesting jet-ISM interactions that may locally compress gas and trigger star formation. If right, the paper is the first to resolve the transition from AGN photoionization to shock-excited cocoon to surrounding LINER and star-forming gas in a LINER-dominated galaxy, and it adds a concrete case of AGN positive feedback.","feed_headline":"Thin shock cocoon found around NGC 5005's faint AGN","feed_subtitle":"HST line-ratio maps trace Seyfert, LINER, and star-forming zones, linking jet-ISM shocks to local star formation.","key_machinery":"The spatially resolved S-BPT (Baldwin–Phillips–Terlevich) diagram built per WFC3 pixel from narrowband [O III], Hβ, Hα, and [S II] images. Each pixel's $\\log([\\mathrm{O\\,III}]/\\mathrm{H}\\beta)$ vs $\\log([\\mathrm{S\\,II}]/\\mathrm{H}\\alpha)$ position assigns it to Seyfert, LINER, or H II classes using Kewley et al. division lines. The new step is splitting LINER pixels at $\\log([\\mathrm{O\\,III}]/\\mathrm{H}\\beta) = -0.28$ into green and blue groups, which spatially separates the inner LINER body from the surrounding high-excitation cocoon. This per-pixel mapping is what allows the paper to locate the cocoon and connect excitation classes to radio jet and CO morphologies.","core_discovery":"The paper's central claim is that NGC 5005's inner few hundred parsecs contain four distinct ionization regimes arranged in a nested structure. At the center, Seyfert-like pixels form a biconical (two opposed cone-like) structure colocated with the optical nucleus and a hard X-ray source; this is photoionization by a low-luminosity AGN, with an estimated bolometric luminosity between $10^{41}$ and $2\\times10^{42}$ erg s$^{-1}$. Enclosing those cones is a thin LINER-like cocoon, 10–20 pc thick, defined by pixels with $\\log([\\mathrm{O\\,III}]/\\mathrm{H}\\beta) \\ge -0.28$ within the LINER region of the S-BPT diagram; the paper argues this cocoon is shock-excited gas, possibly from a $\\sim200$ km s$^{-1}$ outflow or from filtered AGN radiation outside the ionization cones. Beyond it, a $\\sim1$ kpc LINER zone is attributed to AGN photoionization, and a $\\gtrsim2$ kpc LINER zone to post-AGB stars and inflow shocks. Clumpy H II regions at 90–500 pc overlap the 1.5 GHz jet and an Hα bubble, which the paper interprets as jet- or outflow-driven compression that triggers localized star formation, while a $\\sim4$ kpc H II ring is ordinary star formation in a bar-driven molecular ring.","pith_inferences":["A concrete testable extension: if the cocoon is shock-excited, its line ratios should vary with position along the shell and its kinematics should show broadening or velocity gradients of tens to ~200 km/s; this can be checked with IFU spectroscopy at HST resolution.","An implicit consequence the paper does not pursue: the same cocoon-splitting threshold applied to other LINER galaxies could reveal whether thin high-excitation shells are generic around low-luminosity AGN or rare, which would discriminate between outflow-shock and filtered-radiation models.","The positive-feedback interpretation predicts that the young stellar populations in the H II clumps near the jet are younger than stars elsewhere in the ring; resolved color-magnitude or SED fitting of those clumps would test the triggering claim."],"forward_implications":["LINER galaxies can host a faint AGN whose photoionization is visible only within ~20 pc, so nuclear Seyfert-like signatures in BPT maps do not require a powerful quasar.","The thin cocoon provides a resolved example of shock excitation in the ISM, so similar cocoons should be sought around other low-luminosity AGN with HST-class resolution.","Jet-ISM interaction is a viable positive-feedback channel: the H II clumps and the Hα bubble inside 500 pc may be star formation triggered by the radio jet or outflow, not by the AGN radiation field alone.","The kiloparsec LINER zones need not be AGN-powered; in NGC 5005, post-AGB stars and inflow shocks can account for the extended LINER-like emission outside ~2 kpc.","The H II ring at ~4 kpc is consistent with star formation along bar-driven molecular gas streams, implying that large-scale morphology and small-scale AGN feedback coexist in the same galaxy."],"supporting_citations":[{"why":"Supplies the pilot method of spatially resolved BPT mapping and the cocoon morphology that this work extends to a LINER galaxy.","marker":"W. P. Maksym et al. 2016"},{"why":"Extends cocoon identification across Seyfert 2 galaxies and provides methodology comparisons for classifying cocoon-like pixels.","marker":"J. Ma et al. 2021"},{"why":"Supplies the S-BPT classification lines used to separate Seyfert, LINER, and H II pixels.","marker":"L. J. Kewley et al. 2006"},{"why":"Provides the 1.5 GHz eMERLIN radio jet maps used to compare jet-ISM interactions with excitation regions.","marker":"R. D. Baldi et al. 2018"},{"why":"Provides CO(1-0) molecular gas maps and the bar/inflow model used to interpret large-scale LINER and ring emission.","marker":"K. Sakamoto et al. 2000"},{"why":"Provides the cLIER context and post-AGB star ionization scenario for central and extended LINER emission.","marker":"F. Belfiore et al. 2016"},{"why":"Supplies IFU Hα/[N II] measurements used to correct contamination in the F658N filter.","marker":"E. E. Richards et al. 2015"},{"why":"Hydrodynamic simulations of inclined jets interacting with galactic disks are used to interpret the Hα bubble and triggered star formation.","marker":"D. Mukherjee et al. 2018"}],"fun_headline_variants":["Thin shock cocoon encircles NGC 5005's AGN","Jet shocks trigger star formation near NGC 5005's AGN","NGC 5005's AGN wrapped in thin LINER cocoon","Four ionization zones revealed around NGC 5005's nucleus","Shock cocoon and star-forming ring map NGC 5005's heart"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"That the thin LINER cocoon is a genuine physical shell of shock-excited gas rather than an artifact of the eye-chosen $\\log([\\mathrm{O\\,III}]/\\mathrm{H}\\beta) \\ge -0.28$ split acting on gas blurred by the ~10 pc point-spread function, since the claimed shell thickness is only 10–20 pc.","fun_headline_variants_meta":{"raw":{"variants":["Thin shock cocoon encircles NGC 5005's AGN","Jet shocks trigger star formation near NGC 5005's AGN","NGC 5005's AGN wrapped in thin LINER cocoon","Four ionization zones revealed around NGC 5005's nucleus","Shock cocoon and star-forming ring map NGC 5005's heart"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001716,"raw_usage":{"total_tokens":6922,"prompt_tokens":1212,"completion_tokens":5710,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":828,"completion_tokens_details":{"reasoning_tokens":5617}},"tokens_in":828,"tokens_out":5710,"duration_ms":32724,"temperature":1.0,"reasoning_tokens":5617,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T00:10:11.953982+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-map the central 200 pc with an instrument whose PSF is at least a factor of two smaller than the claimed cocoon thickness (or after PSF deconvolution), and with matched IFU kinematics: if the ring of pixels with $\\log([\\mathrm{O\\,III}]/\\mathrm{H}\\beta) \\ge -0.28$ does not remain as a distinct ~10–20 pc shell, and if no shock-broadened line profiles appear at that radius, the cocoon and its positive-feedback interpretation lose their primary support.","supporting_citations":[{"cited_title":"P., Fabbiano, G., Elvis, M., et al","cited_arxiv_id":null,"evidence_quote":"Supplies the pilot method of spatially resolved BPT mapping and the cocoon morphology that this work extends to a LINER galaxy."},{"cited_title":"J., Groves, B., Kauffmann, G., & Heckman, T","cited_arxiv_id":null,"evidence_quote":"Supplies the S-BPT classification lines used to separate Seyfert, LINER, and H II pixels."},{"cited_title":"D., Williams, D","cited_arxiv_id":null,"evidence_quote":"Provides the 1.5 GHz eMERLIN radio jet maps used to compare jet-ISM interactions with excitation regions."},{"cited_title":"J., & Scoville, N","cited_arxiv_id":null,"evidence_quote":"Provides CO(1-0) molecular gas maps and the bar/inflow model used to interpret large-scale LINER and ring emission."},{"cited_title":"E., van Zee, L., Barnes, K","cited_arxiv_id":null,"evidence_quote":"Supplies IFU Hα/[N II] measurements used to correct contamination in the F658N filter."},{"cited_title":"J., Bicknell, G","cited_arxiv_id":null,"evidence_quote":"Hydrodynamic simulations of inclined jets interacting with galactic disks are used to interpret the Hα bubble and triggered star formation."}],"review_version":1}