{"id":"15f358dc-c909-43d1-a26c-a531384cfb9c","arxiv_id":"2502.08619","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"JWST observations reveal extended [Ne VI] emission in the Phoenix cluster, tracing a recent burst of cooling at 10^5.5 K with an inferred rate of 5,000-23,000 solar masses per year.","lead":"Using JWST's MIRI instrument, the authors map the [Ne VI] emission line in the Phoenix cluster and find a large cloud of 10^5.5 K gas that coincides with the cooling peak, cold gas, and star formation. This is the first large-scale map of this temperature phase in a cluster core, and it implies a recent, short-lived cooling spike of 5,000-23,000 solar masses per year.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Cooling-rate inference depends on the unobscured AGN SED assumption; a spatially resolved test of the [Ne VI]/[Ne V] ratio against the AGN-illumination model would settle whether the rate is robust.","rationale":"The detection of extended [Ne VI] is well supported by carefully described PSF subtraction, alternate-pipeline sanity checks, ionization-diagram comparisons against shock/stellar models, and the observed bump in the [Ne VI]/[Ne V] ratio that is insensitive to density. Those results support the qualitative conclusion that gas is cooling through ~10^5.5 K in the cluster core. The weak link is the translation from line luminosity to a mass cooling rate, which depends on the assumed AGN illumination and the neon abundance. The paper's own systematic-error analysis (Supplementary Information §3.2) varies L_bol and Z_Ne by factors of two, but not the SED shape or a heavily obscured line of sight; the text admits obscuration is 'potentially different' from ours. Because the [Ne VI] transition is excited by hard photons and the correction is nonlinear, a large attenuation could reduce the inferred rate by an order of magnitude, which would move the result from a 'spike' to a rate closer to the X-ray-inferred long-term average. This is not an internal inconsistency, but it is the single most load-bearing uncertainty on the number. The reader's weakest_assumption points to the same issue, so agreement is 'agree.'","tokens_in":23393,"tokens_out":1547,"duration_ms":15062,"concrete_test":"Recompute CLOUDY cooling rates for the northern cloud with (a) the nominal unobscured hyperluminous SED, (b) the same SED attenuated by a column of ~1e23 cm^-2, and (c) the minimum possible AGN luminosity consistent with the observed [O III] or narrow-line luminosity after de-reddening, at the same distance and gas density. If the inferred [Ne VI]-derived cooling rate decreases by more than a factor of 3, the headline rate is not robust to the illumination assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim — a short-lived cooling spike of 5,000–23,000 Msun/yr — rests on photoionization corrections computed for gas that sees the full hyperluminous quasar SED. The paper itself flags this as a key systematic in Supplementary Information §3.2, noting that obscuration along the cooling gas line of sight is unknown. If the gas sees a heavily obscured or factor-fewer fainter SED, the inferred cooling rate changes substantially. The single-component fit quoted for [Ne VI] is 25,000 ± 5,000 Msun/yr, and the quoted systematic range is derived from scaling L_bol and Z_Ne by factors of two; it does not include SED-shape/obscuration variations by orders of magnitude, which are plausible for lines of sight near a torus. The cospatiality and detection claims are well supported independently; the cooling-rate scale is not.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents JWST MIRI/MRS integral-field spectroscopy of the Phoenix cluster core, mapping the [Ne VI] λ7.652 μm coronal line over tens of kiloparsecs. The authors develop a spectral-spatial decomposition to subtract the bright QSO PSF and measure extended [Ne VI] emission, finding it cospatial with the lowest-entropy ICM, the cool ionized and molecular gas, and star-forming filaments. A radial surface-brightness profile shows a localized (≲5 kpc) enhancement above an r^-2 AGN-photoionized background, coincident with an increase in the [Ne VI]/[Ne V] ratio. Using CLOUDY cooling models with AGN photoionization, they convert the observed line luminosities into mass cooling rates, obtaining ~25,000 ± 5,000 M⊙ yr−1 for [Ne VI], a multi-line average of ~10,000 ± 5,000 M⊙ yr−1, and a 'conservative' range of 5,000–23,000 M⊙ yr−1; mixing-layer models give 15,000 ± 2,000 M⊙ yr−1 with a range of 7,000–36,000 M⊙ yr−1. They also present the line-of-sight velocity profile, which is consistent with gas in the wake of a buoyantly rising X-ray bubble, and discuss alternative ionization sources (shocks, stellar photoionization, cosmic rays, turbulent mixing).","tokens_in":23576,"tokens_out":13336,"duration_ms":131630,"significance":"If the detection and interpretation hold, this is the first large-scale map of 10^5–10^6 K gas in a cluster core, directly imaging an intermediate-temperature phase of a cooling flow and strengthening the case that the Phoenix cluster is undergoing a rapid, short-lived cooling episode. The paper's strengths include careful data reduction (custom background subtraction, striping removal, PSF subtraction from a stellar template), public code (LOKI on GitHub) and data availability, and a systematic consideration of alternative ionization mechanisms with an ionization diagram. The central detection and the spatial correlations are well supported. The quantitative cooling rate, however, rests on several model choices—most notably the unobscured quasar SED, the fitted isochoric fraction, and the assumed initial conditions—and the quoted systematic range does not capture all of these. The paper is an important advance regardless, but the headline number should be treated as model-dependent until the photoionization systematics are addressed.","major_comments":[{"comment":"The quantitative cooling-rate scale rests on the assumption that the cooling gas is illuminated by the unobscured, full hyperluminous quasar SED. The paper states this assumption and notes in SI §3.2 that the obscuration along the line of sight to the cooling gas is unknown. The systematic range quoted (5,000–23,000 M⊙ yr−1) is derived only from factor-of-two variations in L_bol and Z_Ne; it does not include SED-shape or obscuration changes, which can plausibly alter the photoionization correction by orders of magnitude for lines of sight near the torus. Please test this assumption with the spatially resolved [Ne VI]/[Ne V] ratio map by fitting the incident SED or the ionizing photon flux as a free parameter at each radius, or run a grid of torus-obscured SEDs to bracket the cooling rate. Without such a test, the abstract's numerical range is not robust, even though the detection and cospatiality claims are.","section":"Main text ('We adopt a generic SED template...') and SI §3.2"},{"comment":"The isochoric fraction x is optimized by fitting the observed line ratios under the assumption that Mdot1/Mdot2 ≈ 1 (Eq. 6), and the same x is then used in Eq. (7) to convert each observed line luminosity into a cooling rate. This procedure forces the line-by-line cooling rates toward mutual agreement; the 'median absolute deviation' quoted in the main text is therefore not an independent validation of a single cooling rate. Please fit x and the absolute normalization simultaneously to the observed line luminosities, or propagate the uncertainty in x into each Mdot, and report the best-fit line-ratio residuals. The current treatment makes the multi-line average partly a product of the fitting ansatz.","section":"SI §3.1, Eqs (6)–(7)"},{"comment":"The cooling rates in Figure 3 are computed from the observed [Ne VI] surface-brightness profile, 'assuming that all of the luminosity in the northern cloud is contained within a small Δr at each radius.' However, Figure 2 shows that this profile is decomposed into an r^-2 power-law component, attributed to AGN photoionization of a constant-density medium, and a localized ≤5 kpc box component, attributed to cooling. If the r^-2 component is not actually cooling gas but gas in photoionization equilibrium, then including it in the luminosity entering L = Mdot Γ will overestimate the cooling rate. Please compute Mdot using only the box component, or justify explicitly that the diffuse photoionized gas is also part of the cooling flow.","section":"Main text ('The total [Ne VI] flux...') and Figure 3 caption"},{"comment":"The CLOUDY cooling models assume a single initial density (n_e = 0.42 cm−3) and temperature (kT = 2.3 keV) for the cooling parcel, taken from X-ray measurements. The systematic uncertainty analysis in SI §3.2 varies only L_bol and Z_Ne; the sensitivity of the inferred cooling rate to the adopted initial density and temperature is not explored. Since the ionization parameter and hence the photoionization boost depend on density, please quantify how Mdot changes over the plausible range of initial conditions, or justify explicitly that the result is insensitive to these choices.","section":"Main text ('The cooling gas starts with an initial density...') and SI §3.2"}],"minor_comments":[{"comment":"Typo: 'A raidal surface brightness profile' should read 'A radial surface brightness profile'.","section":"Figure 2 caption"},{"comment":"The sentence 'This is a few times higher than the observed 2.1 ± 0.3 × 10^10 M⊙ yr−1[21]' has incorrect units for a molecular gas mass; it should be M⊙, not M⊙ yr−1.","section":"Main text, 'forming 0.5–2 ×10^11 M⊙ of molecular gas'"},{"comment":"The text includes O VI in the list of lines used for the average cooling rate of 10,000 ± 5,000 M⊙ yr−1, but then states that O VI is an outlier 'most likely due to uncertainties in the UV extinction correction.' Please clarify whether the quoted average includes O VI, and if so, how the outlier affects the median absolute deviation; consider quoting the average with and without O VI.","section":"Main text, 'When we consider all high-ionization lines...'"},{"comment":"The F-test threshold of 0.3% (3σ) is applied spaxel-by-spaxel without correction for multiple testing; given the large number of spaxels, a small number of spurious two-component fits may be expected. Please state whether the two-component kinematic results are robust to a more conservative threshold or a spatial-coherence requirement.","section":"SI §1.4"}],"recommendation":"major_revision","confidential_remarks":"The paper is a strong observational advance with a careful reduction and a convincing detection of extended coronal emission. The main risk is that the abstract's headline cooling rate is presented as a robust range when it depends on an untested unobscured-SED assumption and on a fitting procedure that partly enforces consistency across lines. These issues are addressable within the manuscript's scope, so major revision rather than rejection seems appropriate. I would encourage the editor to request the spatially resolved SED test and a clearer separation of the diffuse and cloud components in the cooling-rate calculation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The headline: this paper gives the first resolved, kiloparsec-scale map of 10^5.5 K gas in a cluster core, using JWST/MIRI to map [Ne VI] in Phoenix. That map is real, and the spatial coincidence with the lowest-entropy ICM, the cold molecular gas, and the starburst is compelling. The cooling-rate number attached to it, 5,000-23,000 M_sun/yr, is plausible but much softer than the map.\n\nWhat's genuinely new: previous coronal-line detections in clusters were unresolved or confined to a few spaxels. Here they have hundreds of spaxels with careful PSF subtraction and a full mid-IR spectral decomposition. The checks against shocks, stellar photoionization, and mixing are honest and mostly convincing. Shocks can't produce the flux, stellar UV can't reach the ionization potential, and the [Ne VI]/[Ne V] bump rules out a pure density effect in the northern cloud. The line ratios are consistent with AGN photoionization plus a cooling component. I buy the qualitative result: gas is cooling through 10^5.5 K in that cloud, and the cospatiality with the cooling peak and star formation is the strongest such case to date.\n\nThe soft spot is the conversion from line luminosity to cooling rate. That conversion goes through CLOUDY, with an assumed AGN SED that the gas sees unobscured along the jet axis, a fitted isochoric fraction x=0.13, and an X-ray-measured Ne abundance. The paper is transparent about the L_bol and Z_Ne scaling, but the quoted 5,000-23,000 range brackets only factor-of-two variations on those two parameters. It does not bracket the possibility that the cooling gas sees an obscured or factor-of-several fainter SED, which would change the photoionization correction and shift the rate much more than the quoted range. The paper itself flags exactly this in Supplementary Section 3.2. So the map and the qualitative cooling claim should be taken as solid; the absolute rate should be treated as model-dependent until there is an independent constraint on the radiation field seen by the gas. The mixing-layer model fits the line ratios better but adds parameters, and the authors are appropriately cautious about it.\n\nThis is for cluster cooling-flow people, AGN feedback modelers, and anyone who wants a new JWST-based window on the 10^5.5 K phase. It deserves a serious referee: the detection will stand, and the rate question is precisely what referees should push on. I would like to see the SED/obscuration issue addressed with a spatially resolved test, for example comparing the [Ne VI]/[Ne V] ratio against the AGN-illumination model at different radii. Send it to review, and ask for hard numbers on the illumination geometry.","headline":"A genuinely new resolved map of 10^5.5 K gas in a cluster core, with a cooling-rate estimate that is plausible but more model-dependent than the quoted range suggests.","tokens_in":24213,"tokens_out":2130,"would_cite":true,"duration_ms":21629,"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 JWST/MIRI observations of [Ne VI] 7.652 µm, the paper maps gas at $\\sim 10^{5.5}$ K across tens of kiloparsecs in the Phoenix cluster core and claims this coronal-phase gas is the cooling flow: it sits with the lowest-entropy X-ray…","keywords":["cooling flows","galaxy clusters","intracluster medium","coronal line emission","Ne VI emission","AGN feedback","Phoenix cluster","JWST MIRI/MRS"],"falsifier":"Re-fit the observed [Ne VI]/[Ne V] and [Mg VII]/[Ne VI] ratios in the northern cloud with CLOUDY models that replace the unobscured hyperluminous-quasar SED by one attenuated by the torus; if those obscured models match the observed ratios equally well, the photoionization correction shrinks and the cooling rate could fall from the claimed 5,000–23,000 $M_\\odot$ yr$^{-1}$ toward the X-ray long-term average.","tokens_in":23181,"feed_emoji":"🌌","tokens_out":7861,"duration_ms":80466,"temperature":0.7,"pith_summary":"The paper claims to have directly imaged a cooling flow in the core of the Phoenix cluster, the first time such a flow has been mapped through the intermediate-temperature gas phase in any galaxy cluster. Using the mid-infrared [Ne VI] line at 7.652 µm, which traces gas near $10^{5.5}$ K, the authors map extended coronal emission over tens of kiloparsecs and show it is cospatial with the lowest-entropy X-ray gas, the coolest optical and molecular gas phases, and sites of active star formation. From this they infer a recent, short-lived spike in the cooling rate of 5,000–23,000 $M_\\odot$ yr$^{-1}$ through $10^{5.5}$ K, far above the long-term X-ray average of a few thousand solar masses per year. If correct, this shows that the \"cooling flow problem\" is not an absolute absence of cooling but a matter of timing: cooling happens in brief bursts, apparently promoted rather than only suppressed by the AGN's radio bubble.","feed_headline":"JWST maps the Phoenix cluster's cooling flow directly","feed_subtitle":"The mapped gas sits amid star formation, implying a short, rapid cooling spike in the cluster core.","key_machinery":"The load-bearing tool is the mid-infrared coronal line [Ne VI] at 7.652 µm, observed with JWST's MIRI/MRS integral-field spectrograph; this line emits near $10^{5.5}$ K with essentially no extinction and no stellar continuum contamination, and a custom spectral decomposition separates the quasar point-spread function from host-galaxy emission across hundreds of spaxels. The conversion from line flux to mass cooling rate is carried out with CLOUDY simulations of a cooling gas parcel exposed to the quasar's radiation field, using a composite of isobaric and isochoric cooling and, in an alternative model, mixing between the hot ICM and warm ISM. The radial surface-brightness profile's bump, together with the bump in the [Ne VI]/[Ne V] ratio, identifies the cooling region as a cloud of size $\\lesssim 5$ kpc located roughly 10 kpc from the nucleus.","core_discovery":"The central discovery is that the $[\\mathrm{Ne\\,VI}]$ $\\lambda 7.652\\,\\mu$m emission in the Phoenix cluster is extended and spatially coherent, forming two clouds to the north of the nucleus that coincide with the minimum entropy of the hot atmosphere, the optical line-emitting gas, the cold molecular gas, and a region of enhanced star formation. Because [Ne VI] forms near $10^{5.5}$ K, this is the first large-scale map of gas between $10^5$ and $10^6$ K in a cluster core. The paper interprets the emission as gas actually cooling out of the hot intracluster medium while being illuminated by the central quasar, with a local bump in the surface brightness and in the density-insensitive $[\\mathrm{Ne\\,VI}]/[\\mathrm{Ne\\,V}]$ ratio marking the switch from AGN photoionization to cooling as the dominant process. The resulting mass cooling rate through $10^{5.5}$ K is 5,000–23,000 $M_\\odot$ yr$^{-1}$, with a mixing-layer variant of the model giving 7,000–36,000 $M_\\odot$ yr$^{-1}$; the authors argue that this spike is short-lived, forming 0.5–$2\\times10^{11}$ $M_\\odot$ of molecular gas and sustained by the turbulent wake of the buoyant X-ray bubble.","pith_inferences":["The same MIRI/MRS technique applied to other cool-core clusters should reveal whether extended [Ne VI] emission is common; clusters with weak or absent coronal-phase nebulae would test the claim that Phoenix is an unusual, short-lived state.","If the cooling gas truly sees an unobscured quasar, then variations in the quasar's ionizing luminosity should imprint on the photoionization-dominated parts of the [Ne VI] map while leaving the cooling component stable; multi-epoch observations could separate these contributions.","The mixing-layer models imply that neon abundances and dust depletion in the $10^{5.5}$ K phase dominate the systematic uncertainty; future mid-infrared neon-line observations in other systems could measure these quantities directly instead of relying on X-ray abundances averaged over 300 kpc.","If short-lived cooling spikes with rates of order $10^4$ $M_\\odot$ yr$^{-1}$ are generic in massive clusters, the molecular gas reservoirs they build should be detectable as extended CO or dust emission in other starbursting brightest cluster galaxies."],"forward_implications":["The Phoenix cluster would become the first system in which a cooling flow has been mapped continuously from the hot intracluster medium down to star-forming temperatures, closing a long-standing observational gap.","The inferred cooling rate through $10^{5.5}$ K exceeds the classical X-ray-inferred, gigayear-averaged rate by a factor of a few, implying that cooling in cluster cores is episodic rather than steady.","The model predicts that this cooling spike forms $0.5$–$2\\times10^{11}$ $M_\\odot$ of molecular gas, several times the observed $2.1\\times10^{10}$ $M_\\odot$, suggesting that much of the cold gas is destroyed by feedback or re-mixed before reaching the molecular phase.","The coincidence of the cooling gas with the wake of the buoyant X-ray bubble indicates that AGN feedback can promote cooling as well as regulate it, providing a concrete mechanism for fueling future feedback cycles."],"supporting_citations":[{"why":"Establishes Phoenix as a massive cooling-flow-induced starburst and supplies the quasar bolometric luminosity used in the models.","marker":"[17]"},{"why":"Supplies the X-ray entropy map, bubble locations, and the identification of the low-entropy gas cospatial with the [Ne VI] clouds.","marker":"[20]"},{"why":"Provides the cold molecular gas mass and filaments used to compare with the predicted cooling products.","marker":"[21]"},{"why":"Gives the previous O VI-based cooling rate and the X-ray data used for comparison and for the ICM pressure/density profiles.","marker":"[23]"},{"why":"Provides the XMM-RGS cooling rate and the hot-phase metal abundances adopted in the CLOUDY models.","marker":"[24]"},{"why":"Maps the optical [O II] and warm ionized gas used to establish cospatiality with active star formation.","marker":"[32]"},{"why":"Cloudy code used to simulate the cooling parcel with AGN illumination and compute line emissivities.","marker":"[37]"},{"why":"Supplies the mean SED template for hyperluminous quasars adopted for the AGN radiation field.","marker":"[38]"},{"why":"Gives the composite isobaric/isochoric cooling prescription used to set the cooling rate.","marker":"[40]"}],"fun_headline_variants":["JWST maps Phoenix cluster's cooling flow directly","First large-scale map of Phoenix cluster's warm gas","JWST catches Phoenix cluster gas mid-cooling","Cooling flow revealed in Phoenix cluster by JWST","Phoenix cluster's cooling spike seen in new light"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The quantitative cooling rate assumes that the cooling gas sees the quasar's full unobscured radiation field; if the gas is shielded by the dusty torus or a fainter source, the photoionization correction shrinks and the inferred cooling rate would drop substantially.","fun_headline_variants_meta":{"raw":{"variants":["JWST maps Phoenix cluster's cooling flow directly","First large-scale map of Phoenix cluster's warm gas","JWST catches Phoenix cluster gas mid-cooling","Cooling flow revealed in Phoenix cluster by JWST","Phoenix cluster's cooling spike seen in new light"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000345,"raw_usage":{"total_tokens":1977,"prompt_tokens":1111,"completion_tokens":866,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":727,"completion_tokens_details":{"reasoning_tokens":792}},"tokens_in":727,"tokens_out":866,"duration_ms":9870,"temperature":1.0,"reasoning_tokens":792,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-08T04:27:03.217315+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Re-fit the observed [Ne VI]/[Ne V] and [Mg VII]/[Ne VI] ratios in the northern cloud with CLOUDY models that replace the unobscured hyperluminous-quasar SED by one attenuated by the torus; if those obscured models match the observed ratios equally well, the photoionization correction shrinks and the cooling rate could fall from the claimed 5,000–23,000 $M_\\odot$ yr$^{-1}$ toward the X-ray long-term average.","supporting_citations":[{"cited_title":"A Massive, Cooling-Flow-Induced Starburst in the Core of a Highly Luminous Galaxy Cluster","cited_arxiv_id":"1208.2962","evidence_quote":"Establishes Phoenix as a massive cooling-flow-induced starburst and supplies the quasar bolometric luminosity used in the models."},{"cited_title":"Anatomy of a Cooling Flow: The Feedback Response to Pure Cooling in the Core of the Phoenix Cluster","cited_arxiv_id":"1904.08942","evidence_quote":"Supplies the X-ray entropy map, bubble locations, and the identification of the low-entropy gas cospatial with the [Ne VI] clouds."},{"cited_title":"Deep Chandra, HST-COS, and Megacam Observations of the Phoenix Cluster: Extreme Star Formation and AGN Feedback on Hundred Kiloparsec Scales","cited_arxiv_id":"1508.05941","evidence_quote":"Gives the previous O VI-based cooling rate and the X-ray data used for comparison and for the ICM pressure/density profiles."},{"cited_title":"AGN feedback in the Phoenix cluster","cited_arxiv_id":"1808.02872","evidence_quote":"Provides the XMM-RGS cooling rate and the hot-phase metal abundances adopted in the CLOUDY models."},{"cited_title":"The State of the Warm and Cold Gas in the Extreme Starburst at the Core of the Phoenix Galaxy Cluster (SPT-CLJ2344-4243)","cited_arxiv_id":"1311.0867","evidence_quote":"Maps the optical [O II] and warm ionized gas used to establish cospatiality with active star formation."},{"cited_title":"Shaping the X-ray spectrum of galaxy clusters with AGN feedback and turbulence","cited_arxiv_id":"1410.7769","evidence_quote":"Gives the composite isobaric/isochoric cooling prescription used to set the cooling rate."}],"review_version":1}