{"id":"ee57c812-ee4e-435b-b52b-648a08989bd0","arxiv_id":"2501.09083","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Cloud-crushing simulations of CGM clouds predict multiphase absorption but fail to produce O VI alongside low ions, implying observed O VI systems need physics beyond these models.","lead":"This paper runs hydro simulations of cold gas clouds flying through the hot gas around Milky-Way-like galaxies and computes the ions those clouds would produce. It finds that different ions can appear in the same absorption feature, but the simulated clouds cannot produce the O VI absorption seen together with cooler gas in real quasar spectra.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equilibrium ionization is baked into both the cooling function (Eq. 3, §2.1) and the post-processing (§4.1), so the O VI deficit may be an artifact of the very assumption the paper concedes could resolve it.","rationale":"The reader correctly identified the equilibrium, optically thin, solar-metallicity post-processing in Section 4.1 as the weakest assumption. My stress-test agrees but sharpens the concern: the same equilibrium assumption is already embedded in the hydrodynamics through the Wiersma et al. (2009) cooling tables used in Eq. (3), which determine the thermal structure of the mixing layers. Because O VI is collisionally ionized in these simulations, the amount of O VI is set by the mass of gas near 10^5.5 K, and that mass is controlled by the assumed cooling. Thus the negative result is conditional on an assumption that affects both the dynamics and the post-processing. This does not change the reader's CONDITIONAL verdict, but it makes the condition more specific and more difficult to dismiss as a mere post-processing caveat. The paper is otherwise careful: the methodology is clearly described, the synthetic spectra are not tuned to observations, and the authors explicitly acknowledge the missing physics. The absence of public code and data is a secondary concern, not the load-bearing one.","tokens_in":26572,"tokens_out":6375,"duration_ms":71030,"concrete_test":"Run one matched cloud-crushing simulation with the same initial conditions and resolution as T0.3_v1700_chi300, but with self-consistent non-equilibrium cooling and time-dependent ionization (e.g., using the MAIHEM network or an equivalent), then repeat the Section 4.2–5.6 absorber fitting and ion-matching analysis at t50, z = 0.540, and θ = 60°. If the O VI column-density distribution shifts upward by ≳0.5 dex, or the H I–O VI matched fraction rises above a few percent, the central negative claim does not survive the equilibrium assumption.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central negative result (Section 6) is that the simulated clouds, after equilibrium, optically thin post-processing with HM12 and solar abundances, underproduce O VI and its coexistence with low ions. This conclusion rests on the equilibrium assumption in two places, not one. First, in the hydrodynamics, the cooling function entering Eq. (3) is taken from Wiersma et al. (2009) equilibrium Cloudy tables, so the gas temperatures in the mixing layers—where O VI forms—are themselves determined assuming ionization equilibrium and solar abundances. Second, in Section 4.1, ion fractions are again computed under ionization equilibrium, an optically thin HM12 background, and solar metallicity. Since Section 5.5 shows O VI is essentially collisionally ionized, the O VI column density is controlled by the mass of gas near T ~ 10^5.5 K. If non-equilibrium cooling modifies that thermal structure—as found in other mixing-layer studies—then the deficit of O VI and its low matched fraction with low ions could be an artifact of the equilibrium assumption applied at both stages. The authors list non-equilibrium effects as a possible explanation, but because the dynamics already assumed equilibrium cooling, this is not merely a post-processing caveat; it is a missing process in the physical model itself. A dedicated non-equilibrium simulation is needed to determine whether the negative claim survives.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"Yang et al. present a parameter study of cloud-crushing simulations in which cool, dense clouds are embedded in a hot supersonic flow, with parameters chosen to approximate clouds in the CGM of Milky Way-like galaxies. They vary the Mach number, initial density contrast, and thermal-conduction level; post-process the snapshots under the assumptions of ionization equilibrium, solar metallicity, an optically thin HM12 photoionizing background, and no self-shielding; generate mock absorption spectra; fit individual absorbers; and study the distributions of column densities, line widths, and ion-ion matching fractions. The main positive result is that a single cloud can produce multiphase absorption, with low- and mid-ionization ions sharing similar velocity structure. The main negative result, stated in the abstract and in Section 6, is that the simulations cannot reproduce O VI and its coexistence with low ions as seen in many observed QSO absorption systems.","tokens_in":26823,"tokens_out":7404,"duration_ms":76005,"significance":"If the negative result is robust, the paper provides a useful constraint: single cold clouds in a uniform hot flow, treated with equilibrium cooling and photoionization, cannot by themselves explain O VI-bearing multiphase absorbers. The methods are described in enough detail to be reproduced, the parameter space is explored systematically, and the statistical predictions (covering areas, line-width distributions, matching fractions) are falsifiable against QSO absorber samples. Section 6 is also commendably explicit in listing non-equilibrium ionization, metallicity, magnetic fields, and collections of clouds as possible ways to produce O VI. The significance is nevertheless bounded by the idealized single-cloud setup and by the equilibrium assumption, which enters both the hydrodynamics and the post-processing.","major_comments":[{"comment":"The cooling function in Eq. (3) is taken from Wiersma et al. (2009) equilibrium Cloudy tables, and the ion fractions in §4.1 are computed under the same ionization-equilibrium assumption. Section 5.5 shows that O VI and Ne VIII are almost entirely collisionally ionized, so the O VI column density is set by the mass of gas near the O VI peak temperature in the mixing layer. If non-equilibrium cooling changes the thermal structure of that layer, as found in other mixing-layer studies, the paper's central negative result—the lack of O VI and of O VI–low-ion coexistence—could be an artifact of applying the equilibrium assumption at both stages. The paper lists non-equilibrium effects as a possible explanation in Section 6, but because the dynamics already assume equilibrium cooling, this is not merely a post-processing caveat. I would like to see a dedicated non-equilibrium simulation, or at least a quantitative estimate of the expected change in O VI column densities, before the negative conclusion is stated as robust.","section":"§2.1, Eq. (3); §4.1; §5.5; §6"},{"comment":"The central negative claim is not tied to a specific observed benchmark. The paper does not state the observed O VI column densities, matching fractions, or line-width thresholds that the simulations fail to reproduce; Figure 14 shows internal matching fractions, but no observed values are overplotted. Without such a benchmark, 'unable to explain' is difficult to evaluate or falsify. Please add a quantitative comparison with a defined observational sample (e.g., column-density thresholds and a matching tolerance) or explicitly limit the claim to 'cannot reproduce the O VI column densities and coexistence fractions found in [specific sample].'","section":"§6 and abstract"},{"comment":"The ion fractions are computed for solar metallicity, an optically thin HM12 background, and no self-shielding, and these choices are not varied. Section 4.2's rescaling of column densities is a uniform multiplicative shift and therefore cannot alter relative ion ratios; it cannot mimic super-solar metallicity, a different UV background, or self-shielding. Since Section 6 itself identifies super-solar metallicity and self-shielding as possible ways to produce more O VI, the abstract's unqualified statement that the simulations are 'unable to explain high ions like O VI' overstates the scope of the study. Either include a small variation of metallicity/radiation field (or self-shielding) or add a qualifier that the conclusion applies to the specific assumptions used.","section":"§4.1, §4.2, §6"}],"minor_comments":[{"comment":"Table 2 lists d25 values up to ~200 kpc for a simulation box of only 8.7 kpc; this is presumably the distance the ambient wind has swept past the cloud, not the displacement of the cloud within the box. Please state this explicitly to avoid confusion.","section":"Table 2"},{"comment":"The affiliation line contains 'Universist¨ at' and several author names with umlauts rendered via LaTeX escapes; these should be corrected in the final typeset version.","section":"Author affiliations"},{"comment":"The statement that H I and Mg II are 'destroyed by photoionization' is a shorthand; more precisely, photoionization lowers their neutral or singly-ionized fractions. Please rephrase for clarity.","section":"§5.5"},{"comment":"The color bar in Figure 14 is labeled 'P', but the text defines the quantity as a matching probability. Use 'matched fraction' or 'P_match' to avoid ambiguity.","section":"Figure 14"},{"comment":"The velocity tolerance Δv = 10 km/s used for ion matching is first introduced in §5.6; it would be clearer to define it in §4.2 where the absorber-fitting method is described.","section":"§4.2"}],"recommendation":"major_revision","confidential_remarks":"The paper is honest and well structured, but the central negative claim about O VI is more fragile than the abstract suggests. The equilibrium assumption enters both the cooling function and the post-processing, so the O VI deficit may be an artifact of the physical model rather than a robust prediction. A non-equilibrium simulation or a clearly scoped, quantitative benchmark would address this. I think major revision is appropriate; if the authors add such a test or meaningfully narrow the claim, the paper could become publishable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this for the most systematic mock-absorber statistics so far from low-density cloud-crushing simulations. The headline negative result—underproduction of O VI and O VI+low-ion coexistence—holds for their model, but the equilibrium assumption is load-bearing in both the hydrodynamics and the post-processing, so the result is more fragile than the paper's framing suggests.\n\nWhat's new: they move cloud crushing into a CGM-like regime (ρ_c ~ 1e-26 g/cm^3, R_c ~ 545 pc, pressure matched to a MW-like halo), select cloud material by tracer fraction rather than density, use off-axis sightlines, and fit per-component Gaussians instead of mean projection column densities. That combination is a real step beyond Cottle et al. 2018, de la Cruz et al. 2021, and Casavecchia et al. 2024. The analysis is careful throughout: mass-loss curves, phase-space maps, parameter dependence across Mach number, density contrast, conduction, projection angle, and redshift. The matched-fraction analysis of ion coexistence is a genuinely useful way to quantify multiphase structure.\n\nThe soft spot is the one the stress-test flags, and I think it lands. The cooling function in Eq. (3) uses Wiersma et al. (2009) equilibrium Cloudy tables, and the ion fractions in §4.1 are again computed assuming equilibrium, optically thin gas, HM12 background, solar metallicity. Since §5.5 shows O VI and Ne VIII are almost entirely collisionally ionized, their columns trace the mass of gas near T ~ 1e5.5 K. Non-equilibrium cooling can change that thermal structure—mixing-layer studies have shown this—so the O VI deficit could be an artifact of the dynamics, not just the ionization post-processing. The authors do list non-equilibrium in §6, but as one possible explanation for the mismatch with observations. That undersells it: it's a missing process in the physical model, and a dedicated non-equilibrium run is needed before the negative claim is solid. Minor issues: no public data or code, and the optical-depth normalization before fitting is sensible for separating physical components but makes direct observational comparison indirect, though they handle this with scaling fairly.\n\nWho gets value: people calibrating subgrid wind models like PhEW, and observers wanting a theory-side statistical baseline for CGM absorbers. This deserves a serious referee. I'd send it, with a major-revision request for either a non-equilibrium simulation or a substantially harder caveat in the abstract.","headline":"The O VI deficit is real for their model, but equilibrium ionization is baked into both the dynamics and the post-processing, so the negative result is more fragile than the paper lets on.","tokens_in":27394,"tokens_out":3041,"would_cite":true,"duration_ms":31126,"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":"Cloud-crushing simulations reproduce multiphase QSO absorbers, but fail to produce enough O VI to match observations.","keywords":["circumgalactic medium","cloud-crushing simulations","QSO absorption lines","multiphase gas","ionization equilibrium","O VI","thermal conduction","galactic outflows"],"falsifier":"Take the same simulation snapshots and recompute ion densities under non-equilibrium ionization or with super-solar metallicity; if the O VI column densities and their alignment with low ions then reach the observed values, the paper's central negative conclusion would be overturned, whereas if they remain low, the conclusion that single CGM-like clouds cannot produce the observed O VI systems would be confirmed.","tokens_in":26345,"feed_emoji":"🔭","tokens_out":5909,"duration_ms":51711,"temperature":0.7,"pith_summary":"This paper asks whether a single cold cloud launched into the hot circumgalactic medium can explain the multiphase quasar absorption systems observed along many sightlines. Using cloud-crushing simulations with densities typical of a Milky Way-like CGM, the authors show that the disrupted cloud produces absorbers in which low-ionization species and middle-ionization species coexist at similar velocities and line widths, so a single outflow cloud can mimic a multiphase absorber. The same simulations, however, produce very little O VI and almost no O VI absorbers that align with low ions, so the model cannot account for the observed O VI systems. The outcome matters because it defines what a single-cloud outflow model can and cannot do, and it points to which physics must be added to match observations.","feed_headline":"Crushed clouds make multiphase absorbers, miss O VI","feed_subtitle":"A single cold outflow cloud aligns low- and mid-ion absorption, but cannot explain the observed O VI systems","key_machinery":"The central machinery is the cloud-crushing simulation: a uniform spherical cloud of mass $10^5\\,\\mathrm{M}_\\odot$, density $10^{-26}\\,\\mathrm{g\\,cm^{-3}}$, and temperature $10^4\\,\\mathrm{K}$ moving through a hot ambient flow in thermal pressure equilibrium, characterized by the Mach number $M_{\\rm hot}$ and the initial density contrast $\\chi_0$, with radiative cooling and either no, weak, or full thermal conduction. Ion densities are post-processed from the simulated gas using ionization-equilibrium tables under the HM12 photoionizing background, and mock absorption spectra are generated and decomposed into individual Gaussian absorbers. The mixing layers created by Kelvin-Helmholtz instabilities, and the suppression of those instabilities by conduction, are what make different ions occupy different spatial regions while still appearing aligned in velocity space.","core_discovery":"Working in the framework of cloud-crushing simulations, the paper's central discovery is that a single cool cloud with CGM-like density, crushed by a hot supersonic flow, generates absorption systems with ions of widely different ionization potentials appearing in the same velocity component and with comparable line widths, yet the same simulations fail to produce significant O VI, and in particular fail to produce O VI aligned with low-ionization species. The authors present this as a positive and a negative result: the multiphase absorber phenomenon can emerge from the mixing layers and fragments of one cloud, but the observed coexistence of O VI with H I and low ions requires physics beyond their default assumptions of ionization equilibrium, an optically thin HM12 photoionizing background, and solar metallicity.","pith_inferences":["If non-equilibrium ionization or super-solar metallicity were included in the same snapshots, O VI production could plausibly rise enough to match observations; this is a direct, testable extension of the paper's own post-processing pipeline.","The paper's negative O VI result indirectly suggests that much of the observed O VI may trace the hot ambient CGM itself rather than the cool outflow clouds, with velocity alignment arising from shared outflow kinematics rather than co-location.","The absorber-matching method used here could be turned into an observational diagnostic: the predicted distribution of velocity offsets and line-width ratios between low ions and O VI differs noticeably between conductive and non-conductive models, providing a way to infer the thermal conduction state of the CGM.","Because the simulations use a uniform ambient medium, a next step would be to embed a crushed cloud in a realistic, density-stratified CGM; a lower ambient density encountered as the cloud travels outward would likely prolong cloud survival and change the ion ratios."],"forward_implications":["A single disrupted cloud can yield the aligned low- and mid-ion absorbers seen in QSO spectra, so multiphase absorption does not by itself require multiple physically separate gas phases along the sightline.","The orientation of the sightline matters: head-on projections show much higher column densities than all other angles, so down-the-barrel and transverse QSO observations of the same cloud should differ systematically.","Thermal conduction changes the absorber population: non-conductive clouds fragment and produce broader, multi-component absorption, while conductive clouds stay coherent and produce narrower, higher-column-density low-ion components.","Photoionization from the HM12 background destroys H I and Mg II, boosts O II, Si III, and C III, and leaves O VI and Ne VIII essentially collisionally ionized, so the high-ion content of these absorbers is set by temperature rather than radiation.","The failure to produce O VI alongside low ions implies that observed O VI systems must come from larger or more massive clouds, super-solar metallicity, non-equilibrium ionization, magnetic-field effects, or sightlines through multiple clouds."],"supporting_citations":[{"why":"Supplies the cloud-crushing simulation setup and dimensionless parameters ($M_{\\rm hot}$, $\\chi_0$) that this work adapts to lower CGM-like densities.","marker":"Scannapieco & Brüggen (2015)"},{"why":"Provides the equations of motion, the thermal-conduction treatment, and the mass-loss phenomenology that the simulations are built on.","marker":"Brüggen & Scannapieco (2016)"},{"why":"Defines the cloud crushing time $t_{\\rm cc}$ used to set the simulation snapshots.","marker":"Klein et al. (1994)"},{"why":"Establishes the prior approach for computing ion densities and generating absorption lines from cloud-crushing simulations that this paper extends.","marker":"Cottle et al. (2018)"},{"why":"Supplies the HM12 photoionizing background used in the ionization-equilibrium tables.","marker":"Haardt & Madau (2012)"},{"why":"Provides the cooling tables with solar abundances used in the hydrodynamical evolution.","marker":"Wiersma et al. (2009)"},{"why":"Provides the tool used to generate mock absorption spectra and fit individual absorbers.","marker":"Hummels et al. (2017)"},{"why":"Provides the photoionization code version used to compute ionization fraction tables under ionization equilibrium.","marker":"Chatzikos et al. (2023)"}],"fun_headline_variants":["One cloud crash yields multiphase absorbers, no O VI","Cloud-crushing simulations align ions, miss O VI","Single cold cloud mimics CGM absorption, lacks O VI","Crushed cloud shows mixed ions, fails to explain O VI"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The whole negative result about O VI rests on the post-processing assumption that the gas is in ionization equilibrium, optically thin, and solar-metallicity, with the HM12 photoionizing background; if any of these is wrong, O VI could be produced in the same simulated clouds.","fun_headline_variants_meta":{"raw":{"variants":["One cloud crash yields multiphase absorbers, no O VI","Cloud-crushing simulations align ions, miss O VI","Single cold cloud mimics CGM absorption, lacks O VI","Crushed cloud shows mixed ions, fails to explain O VI"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0003,"raw_usage":{"total_tokens":1667,"prompt_tokens":816,"completion_tokens":851,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":432,"completion_tokens_details":{"reasoning_tokens":781}},"tokens_in":432,"tokens_out":851,"duration_ms":9410,"temperature":1.0,"reasoning_tokens":781,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T20:10:47.596037+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take the same simulation snapshots and recompute ion densities under non-equilibrium ionization or with super-solar metallicity; if the O VI column densities and their alignment with low ions then reach the observed values, the paper's central negative conclusion would be overturned, whereas if they remain low, the conclusion that single CGM-like clouds cannot produce the observed O VI systems would be confirmed.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the cloud-crushing simulation setup and dimensionless parameters ($M_{\\rm hot}$, $\\chi_0$) that this work adapts to lower CGM-like densities."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the HM12 photoionizing background used in the ionization-equilibrium tables."},{"cited_title":"B., Smith B","cited_arxiv_id":null,"evidence_quote":"Provides the tool used to generate mock absorption spectra and fit individual absorbers."}],"review_version":1}