{"id":"9f0aff37-f165-40c5-8494-d9ac95b44990","arxiv_id":"2411.17173","paper_version":3,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":9,"one_line_summary":"Modeling the cool CGM as ~10^3 misty cloud complexes of radius ~10 kpc matches MgII absorption data, and the product of average column density and covering fraction yields a robust cool gas mass estimate.","lead":"This paper proposes that the cool gas around galaxies is gathered into thousands of cloud complexes, each made of a mist of tiny cloudlets, rather than spread smoothly. The authors show that a simple model with a total cool gas mass near 10^10 solar masses can match observed magnesium absorption and offers a way to estimate that mass from two observed quantities.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fiducial turbulent-broadening ansatz (hot-phase Mach 0.5, §3.4.2) is unvalidated and directly sets both the MgII EW distribution and the EW-threshold covering fraction used in the mass proxy; adopting Mach 0.2 or 1.0 could shift the inferred cool CGM mass substantially.","rationale":"The reader's weakest assumption identifies the same load-bearing point: the turbulent broadening ansatz in §3.4.2 is an unvalidated input that directly controls the EW distribution and the threshold-based covering fraction, and therefore also affects the mass-proxy calibration. My independent reading of the manuscript confirms that this is the most fragile link in the central claim. The paper has real strengths: the mist-limit treatment is analytically tractable, the Monte Carlo machinery is described in enough detail to be reproduced, the public code is cited, and Section 5 provides an independent but computationally limited check that the fiducial parameters survive when cloudlets are modeled explicitly. The circularity of fitting the same COS-Halos data is a secondary concern, partially mitigated by the comparison with Huang et al. (2021) covering fractions, but the mass estimate itself is still tied to the same sample. The decisive issue remains that a factor ~2.5 change in b_turb,CC, well within observational uncertainty, would change saturated EWs and threshold covering fractions by order-unity amounts; without a test at Mach ~0.2 or ~1.0, the claimed robustness of the 10^10 M_sun mass is not established. The recommended concrete test would settle whether this concern actually lands, so the reader's conditional verdict should stand unchanged.","tokens_in":36446,"tokens_out":8077,"duration_ms":83470,"concrete_test":"Run the fiducial Monte Carlo (§3.4) twice more with the hot-phase turbulent Mach number set to 0.2 and 1.0 (σ_3D,turb,CGM ≈ 43 and 214 km/s), keeping all other parameters fixed, and regenerate Fig. 6 (EW distribution and EW>0.3 Å covering fraction) and Fig. 7 (mass-proxy contours). If the COS-Halos/Huang data then prefer a different M_cool by more than a factor of 2, or if the product ⟨N_MgII⟩×f_A no longer collapses onto the 10^10 M_sun contour, the turbulent-broadening ansatz is load-bearing and the quoted mass is not robust. An even more direct check is to repeat Section 5's explicit cloudlet modeling with the same two Mach numbers; if the fiducial parameters (N_CC=10^3, R_CC=10 kpc, M_cool=10^10 M_sun) survive, the concern is mitigated.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central comparison with observed MgII equivalent widths and the calibration of the covering-fraction mass proxy both run through the assumed turbulent broadening across a cloud complex (§3.4.2). The paper sets σ_3D,turb,CGM ≈ 107 km/s by assuming a hot-phase Mach number of 0.5 at ~2×10^6 K, then scales to CC scale with a Kolmogorov factor (R_CC/R_CGM)^{1/3}, giving b_turb,CC ≈ 20 km/s. No direct observation of hot CGM turbulence in L* halos is cited; the only empirical anchor (Hitomi) is ICM Mach ~0.2. Because MgII column densities near 10^13 cm^-2 lie on the flat part of the curve of growth (Fig. 3), the predicted EW is sensitive to b_turb, not just N_MgII. The EW-threshold covering fraction (Fig. 6, bottom; EW>0.3 Å) therefore moves with b_turb, and §6.3.1 explicitly recommends using such an EW threshold to define the covering fraction entering the mass proxy. Thus the proxy advertised as 'independent of other model parameters' is not independent of b_turb. If the true hot-phase Mach number is 0.2 instead of 0.5, b_turb,CC falls to ~8 km/s; the EW of saturated components shrinks and the threshold covering fraction drops, so matching the observed covering fraction would require a different M_cool. The claimed ~10^10 M_sun cool CGM mass and the 'reproduces' validation are therefore contingent on an unvalidated factor ~2.5 in b_turb.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper introduces the 'mCC' (misty cloud complex) model, in which the cool (~10^4 K) circumgalactic medium around L* galaxies is not a uniform mist but is concentrated in ~10^3 clumpy complexes of radius ~10 kpc, each internally in the mist limit with a large number of tiny cloudlets. Analytic expressions are derived for the mean MgII column density and its variance along lines of sight for uniform and power-law radial distributions of complexes, and Monte Carlo realizations are used to generate predicted distributions of column density, equivalent width, and covering fraction. Comparing with COS-Halos and M3/Huang et al. data, the authors conclude that M_cool ~ 10^10 M_sun, N_CC ~ 10^3, R_CC ~ 10 kpc, and alpha = 1 reproduce the observed trends; they further propose that the product of area-averaged MgII column density and area covering fraction is a robust mass proxy. The paper also presents an 'advanced' model with distributed complex sizes and masses, direct cloudlet-level modeling over the full CGM, and an analogous OVI prediction for warm gas.","tokens_in":36827,"tokens_out":7251,"duration_ms":67628,"significance":"If validated, the mCC framework would be a useful and computationally cheap bridge between analytic misty-CGM models and cloudlet-resolved models such as CloudFlex. The paper's strengths include publicly released code and data, closed-form analytic expressions (Eqs. 6 and 9, and the variance estimate in Appendix C), a careful Monte Carlo description, a concrete mist-limit criterion (Appendix A), and an explicit demonstration that line blending of small cloudlets converges to the misty-CC absorption profile (Section 4). The main physical claims, however, rest on two load-bearing choices that the manuscript itself flags as uncertain: the turbulent broadening ansatz of Section 3.4.2, which has no direct observational calibration for L* halos, and the calibration of the fiducial model on the same COS-Halos data that are later used to 'infer' the cool CGM mass. These issues do not invalidate the framework as a phenomenological tool, but they currently prevent the paper from delivering a robust quantitative constraint on the cool CGM mass.","major_comments":[{"comment":"The turbulent broadening across a cloud complex is the main link between the model's column densities and the observed equivalent widths, yet it is set by an unvalidated ansatz: a hot-phase Mach number of 0.5 at ~2x10^6 K, scaled to the CC scale with a Kolmogorov factor. Footnote 6 concedes that no direct observations of hot CGM turbulence in Milky Way-mass halos exist, and the only empirical anchor cited (Hitomi) gives ICM Mach ~0.2. Because MgII column densities near 10^13 cm^-2 lie on the flat part of the curve of growth (Fig. 3), the predicted EW and the EW-threshold covering fraction (Fig. 6, bottom panel; also the recommendation in Section 6.3.1) depend directly on the adopted broadening. If the true hot-phase Mach number were 0.2 instead of 0.5, the 1D turbulent dispersion across a 10 kpc CC would drop from 20 to ~8 km/s and b_turb,CC from ~28 to ~11 km/s; the saturated-component EWs would shrink and the EW>0.3 A covering fraction would fall, so matching the observed covering fraction would require a different M_cool and/or different CC parameters. The claimed ~10^10 M_sun mass and the statement that the model 'reproduces' the observations are therefore contingent on a factor-of-2.5 uncertainty in the broadening parameter that is not propagated into the quoted constraint.","section":"Section 3.4.2 and Section 3.6"},{"comment":"The inference of M_cool ~ 10^10 M_sun is circular as presented. The fiducial mass is selected in the top panels of Fig. 5 by matching the COS-Halos column density distribution, and the same COS-Halos data are then placed on the mass-proxy diagram of Fig. 7 to conclude that the cool CGM mass is ~10^10 M_sun. The green point in Fig. 7 thus confirms a value that was already built into the model through Fig. 5, not an independent measurement. The mass-proxy relation itself is a model prediction, but its application to COS-Halos is not a test. To support the title's 'constraining' claim, the model should be calibrated on one sample (or on a subset of sightlines) and validated on an independent one (for example, the M3/Huang et al. EW data, or a different ion such as OVI with an independent mass estimate), or the manuscript should be reframed as a demonstration of the framework rather than a measurement of the cool CGM mass.","section":"Section 3.5 and Section 3.7"},{"comment":"The claim that the fiducial parameter set 'best matches' or 'most consistent with' the COS-Halos data is based on qualitative, by-eye comparison of scatter plots with highly heterogeneous data (detections, upper limits, and lower limits). No quantitative goodness-of-fit statistic, likelihood, or sensitivity of the visual ranking to the adopted binning is provided. Given the large intrinsic scatter and the presence of censored measurements, a quantitative comparison (for instance, a two-dimensional Anderson-Darling or KS test with proper handling of limits, or a likelihood over the impact-parameter distribution) is needed before 'reproduces' can be taken as a rigorous claim. This is fixable and would materially strengthen the paper's central conclusion.","section":"Section 3.5 and Fig. 5"},{"comment":"The advertised robustness of the mass proxy is overstated. Fig. 7 shows that the product of the detection-averaged column density and the covering fraction is approximately constant for different N_CC and R_CC at fixed M_cool, which is a useful degeneracy-breaking result. However, the product is not independent of all other model parameters: it depends on the assumed cool-gas density normalization (Eq. 13), metallicity (0.3 Z_sun), and the MgII ion fraction model (Appendix B), all of which are fixed assumptions rather than fitted quantities. Moreover, Section 6.3.1 recommends defining the covering fraction with an EW threshold of 0.3 A, which reintroduces the dependence on the turbulent broadening ansatz discussed above. The text should state explicitly that the proxy is independent only of the CC configuration parameters (N_CC, R_CC, alpha) conditional on the assumed gas physics, and that the EW-threshold version is sensitive to b_turb.","section":"Section 3.7 and Section 6.3.1"}],"minor_comments":[{"comment":"The notation for the turbulent broadening parameter is confusing: the text gives sigma_3D,turb,CC = 35 km/s, then sigma_1D = 20 km/s, and then defines b_turb,CC = sqrt(2) sigma_1D, which gives b_turb,CC ~ 28 km/s, not 20 km/s. Please define a single symbol for the Doppler broadening parameter and use it consistently in Figs. 3 and 4 and in the text.","section":"Section 3.4.2"},{"comment":"The dotted lines representing the expected standard deviation around the mean column density are not labeled in the legend; please add an explicit legend entry and state in the caption whether the spread is computed from Eq. C2 or from the Monte Carlo realizations.","section":"Fig. 2, bottom panel"},{"comment":"The cloudlet-generation cylinder in Section 4 has a height of 20 kpc while the CC radius is 10 kpc; please justify why the cylinder extends beyond the CC, or state explicitly that this is a numerical convenience for the line-blending test.","section":"Section 4"},{"comment":"The volume-fraction estimate in Eq. (10) uses n_cool = 10^-2 cm^-3, while the adopted density profile in Eq. (13) is radius-dependent and gives n_cool = 10^-3 cm^-3 at R_CGM = 280 kpc. Please clarify which density is used for the fiducial numerical values in Eq. (10).","section":"Eq. (10)"},{"comment":"The description of the M3 sample states a stellar mass range of 2x10^8 - 4x10^11 M_sun with median 4x10^10 M_sun, but the comparison in Figs. 6 and 12 treats the sample as comparable to COS-Halos; a brief statement on the stellar-mass and redshift matching between the two samples would help the reader assess the validity of the joint comparison.","section":"Section 2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is well within the scope of a general astrophysics journal and has a genuinely useful, low-cost modeling framework with released code. My main concerns are the circular calibration on COS-Halos, the unvalidated turbulent broadening ansatz and its direct effect on the EW and covering-fraction comparisons, and the absence of a quantitative goodness-of-fit test. These are fixable with reframing and additional analysis, but they are central to the paper's strongest claims, so I recommend major revision rather than rejection. If the authors can validate the mass proxy on an independent data set and propagate the uncertainty in b_turb, the resulting paper would be a solid contribution."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The mCC model is a genuinely practical addition to the CGM toolkit. Treating the cool gas as ~10^3 misty cloud complexes and averaging over cloudlets analytically makes Monte Carlo generation of MgII column densities, EWs, and covering fractions cheap and reproducible. The code and data are public, the analytic estimates in Eqs. 6 and 9 are straightforward, and the line-blending demonstration in Section 4 is the clearest part of the paper: it shows how unresolved cloudlets produce effective turbulent broadening at CC scale. The proposed mass proxy, <N_MgII> times the area covering fraction, is also a useful idea, and Fig. 7 shows it is fairly insensitive to how the CCs are arranged for a fixed cool mass. Those pieces deserve credit, and the paper is worth reading for anyone working on quasar absorption statistics.\n\nThe soft spots are exactly where the paper oversells itself. The fiducial parameters (M_cool = 10^10 Msun, N_CC = 10^3, R_CC = 10 kpc, alpha = 1) are chosen to match the COS-Halos and M3 data, and then the same data are used to say the inferred cool CGM mass is ~10^10 Msun. That is circular for the absolute value. The proxy itself may still be useful for comparing samples or tracking relative mass, but the paper's abstract and Section 6.3.1 claim it is “independent of other model parameters.” That is not supported. The covering fraction in Fig. 6 is defined with an EW threshold, and the EW depends on b_turb through the flat part of the curve of growth. The b_turb value comes from an unvalidated ansatz: a hot-phase Mach number of 0.5, scaled to CC radius with Kolmogorov scaling. The only direct empirical anchor, Hitomi in the ICM, gives Mach ~0.2. Dropping to Mach 0.2 reduces b_turb by roughly a factor of 2.5, which shifts the predicted EW distribution and the threshold covering fraction; matching the observed covering fraction would then require a different N_CC or M_cool. So the claimed 10^10 Msun mass and the “reproduces observations” statement are conditional on that turbulence assumption. The paper should tone down the robustness claim and either calibrate b_turb against simulations or treat it as a free parameter in the inference.\n\nMinor but worth noting: metallicity and ion fraction are fixed rather than marginalized, which also affects the column density and hence the mass estimate. The advanced model partially addresses CC size/mass distributions, but it does not break the circularity.\n\nWho is this for? CGM observers and modelers who want a quick phenomenological model for MgII and OVI absorption statistics. It is not a definitive measurement of the cool CGM mass, but it is a solid framework that can be improved and applied to independent data (MAGG, CUBS, etc.). I would send it to peer review, and I would ask for a sensitivity analysis on b_turb and metallicity, plus an application to a survey not used for calibration. The paper is honest about many idealizations, and the central framework is sound; the circularity and the turbulence dependence are the main things standing between this and the strong claims made in the abstract.","headline":"A useful, computationally cheap toy model for the clumpy cool CGM, but the headline mass is calibrated on the same data it claims to infer.","tokens_in":37417,"tokens_out":2756,"would_cite":true,"duration_ms":30872,"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":"The cool circumgalactic medium around Milky Way-like galaxies can be described as roughly 1,000 ten-kiloparsec cloud complexes, each a mist of tiny cloudlets, holding about 10^10 solar masses of gas.","keywords":["circumgalactic medium","cloud complexes","MgII absorption","quasar absorption lines","covering fraction","multiphase gas","turbulent broadening","cool gas mass"],"falsifier":"Measure the line-of-sight velocity dispersion of the hot CGM around a Milky Way-mass galaxy at ~10 kpc scales with high-resolution X-ray spectroscopy; if the one-dimensional dispersion is much larger or smaller than the assumed ~20 km/s, the predicted MgII equivalent-width distribution from the fiducial mCC model will not match the observed spectra, falsifying the turbulent-broadening ansatz.","tokens_in":36203,"feed_emoji":"☁️","tokens_out":6571,"duration_ms":57991,"temperature":0.7,"pith_summary":"The paper argues that the ~$10^{4}$ K gas in the circumgalactic medium of Milky Way-like galaxies is not spread uniformly but is gathered into 'cloud complexes'—clumps a few to tens of kiloparsec across, each containing a mist of tiny cloudlets. Treating each complex in the mist limit (so a sightline through it always sees cool gas) and placing ~$10^{3}$ such complexes by Monte Carlo, the authors find they can reproduce the MgII column densities, equivalent widths, and covering fractions measured by the COS-Halos and M3 surveys. Their fiducial solution is a power-law radial distribution dN_CC/dR ∝ $R^{{-1}}$ with ~$10^{3}$ complexes of radius ~10 kpc and a total cool-gas mass of ~$10^{10}$ M_sun. A second result is that the product of the area-averaged MgII column density and the area covering fraction is an observationally accessible proxy for the cool CGM mass, nearly independent of how the mass is split into complexes. If correct, the model would turn large sightline-to-sightline scatter in absorption from a nuisance into a diagnostic of cloud structure.","feed_headline":"Cool halo gas sits in ~1,000 giant cloud complexes","feed_subtitle":"A new model matches MgII absorption in galaxy halos and yields a simple way to weigh the cool gas.","key_machinery":"The load-bearing object is the misty cloud complex: a spherical region of radius R_CC that contains cool cloudlets dense enough that any sightline through the complex intersects cool gas (unit area covering fraction), allowing the complex to be smoothed into an average density ⟨n_gas⟩ = 3M_CC/($4πR_CC^{3}$ μ m_p). The model's second mechanism is a Kolmogorov turbulent broadening ansatz, σ_turb,CC = σ_turb,CGM (R_CC/R_CGM)^{1/3}, anchored to a hot-phase Mach number ~0.5, which sets the Doppler parameter b_tot and therefore the equivalent width from the curve of growth. The primary observable identity is ⟨N_MgII⟩ × f_CC^A ≈ const for fixed cool mass, which connects the data products to M_cool. The advanced model replaces identical complexes with power-law distributions of complex size, mass, and radius, but preserves the same mass-proxy relation.","core_discovery":"The central discovery is that the observed clumpiness of the cool CGM can be captured by a two-level 'misty cloud complex' (mCC) model: within each complex, tiny cloudlets are so numerous that the area covering fraction is unity, while the complexes themselves cover only a small fraction of the CGM volume. Monte Carlo realizations with $10^{3}$ complexes of radius 10 kpc, a power-law radial distribution with index α=1, and total cool mass $10^{10}$ M_sun match the MgII column density and equivalent-width trends with impact parameter in the COS-Halos sample and the covering fraction data of Huang et al. 2021. The authors further show that ⟨N_MgII⟩ × f_CC^A, the average column density times the area covering fraction, is nearly constant for a fixed cool-gas mass across wide variations in complex number, size, and mass, so the pair of observables can be used to weigh the cool CGM. Directly placing parsec-scale cloudlets inside complexes shows that absorption from many cloudlets blends to reproduce the same turbulent broadening that the analytic mist prescription assumes.","pith_inferences":["Future surveys with many quasar–galaxy pairs could map cool CGM mass as a function of stellar mass and environment simply by binning average column density and covering fraction without velocity-resolved fits; this extension is not developed in the paper.","The turbulent broadening ansatz is the least protected link: a direct measurement of hot CGM turbulence around Milky Way-mass halos would either validate the ~20 km/s CC-scale dispersion or require re-fitting the CC size and number, since EW and column density are degenerate there.","The two-level mist picture suggests that unresolved multiphase gas in cosmological simulations could be represented statistically by CCs rather than resolved cloudlets, a subgrid scheme the paper mentions but does not implement.","The same ⟨N⟩ × f_cov proxy could in principle be applied to FRB dispersion measures to cross-check the cool-mass estimate against the total electron column, since dispersion measure is insensitive to metallicity and ionization corrections."],"forward_implications":["For Milky Way-like galaxies, the cool CGM mass is constrained to ~10^10 M_sun, with about 10^3 cloud complexes of radius ~10 kpc distributed with a dN/dR ∝ R^{-1} profile.","Large intrinsic scatter in MgII column density at fixed impact parameter is a natural prediction of the patchy CC distribution, not an observational artifact.","The product of average MgII column density and covering fraction can be used as a mass estimator for cool CGM without needing to know the detailed cloud geometry.","Spectra toward a single complex should show blended absorption that mimics turbulent broadening even if individual cloudlets are tiny, implying that resolved components do not directly map to 3D structures.","The same machinery applied to OVI reproduces the COS-Halos OVI column spread, suggesting the CC picture extends to warm gas."],"supporting_citations":[{"why":"Supplies the COS-Halos MgII column density and equivalent width data against which the fiducial model is compared.","marker":"Werk et al. 2013"},{"why":"Provides the best-fit MgII equivalent width–impact parameter relation and the covering-fraction measurements used to validate the model.","marker":"Huang et al. 2021"},{"why":"The mCGM model this paper extends; its uniform-mist baseline column density motivates the need for patchy cloud complexes.","marker":"Dutta et al. 2024b"},{"why":"Introduces the CloudFlex cloudlet-within-complex approach that the mCC model reproduces in the mist limit.","marker":"Hummels et al. 2024"},{"why":"Defines the KS18 ultraviolet background used in CLOUDY computations of the MgII ion fraction.","marker":"Khaire & Srianand 2019"},{"why":"Provides the CLOUDY photoionization code used to compute MgII ion fractions at z≈0.2.","marker":"Ferland et al. 2017"},{"why":"Motivates the n_cool ∝ R^{-1} cool-gas density profile adopted for the radial variation of ion fractions.","marker":"Stern et al. 2016"}],"fun_headline_variants":["Cool halo gas: 1,000 cloudy blobs match MgII","Weigh cool halo gas with one simple observable pair","Misty cloud complexes explain patchy halo absorption","New model for cool CGM: tiny cloudlets, huge complexes","Patchy cool halo gas: model gives mass estimate"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The equivalent widths and covering fractions that validate the model depend on an assumed internal velocity spread for each cloud complex, estimated from hot-gas turbulence whose strength has not been directly measured; if that spread is wrong, the inferred complex size and number change.","fun_headline_variants_meta":{"raw":{"variants":["Cool halo gas: 1,000 cloudy blobs match MgII","Weigh cool halo gas with one simple observable pair","Misty cloud complexes explain patchy halo absorption","New model for cool CGM: tiny cloudlets, huge complexes","Patchy cool halo gas: model gives mass estimate"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000203,"raw_usage":{"total_tokens":1476,"prompt_tokens":1125,"completion_tokens":351,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":741,"completion_tokens_details":{"reasoning_tokens":269}},"tokens_in":741,"tokens_out":351,"duration_ms":4713,"temperature":1.0,"reasoning_tokens":269,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:25:15.545724+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the line-of-sight velocity dispersion of the hot CGM around a Milky Way-mass galaxy at ~10 kpc scales with high-resolution X-ray spectroscopy; if the one-dimensional dispersion is much larger or smaller than the assumed ~20 km/s, the predicted MgII equivalent-width distribution from the fiducial mCC model will not match the observed spectra, falsifying the turbulent-broadening ansatz.","supporting_citations":[],"review_version":1}