{"id":"054be137-7003-4171-8685-246d4e36608e","arxiv_id":"1909.00849","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"A two-component model of the Milky Way hot halo, fit to X-ray emission measures, predicts a Galactic halo dispersion measure of 30 to 245 pc cm^-3, with a sky average of 43 pc cm^-3.","lead":"This paper models the hot gas halo around the Milky Way using diffuse X-ray observations and estimates the halo's contribution to the dispersion measure, a radio-signal delay, of distant fast radio bursts. The model gives a Galactic halo correction of about 43 pc cm^-3 on average, a quantity observers need to subtract before using fast radio bursts to measure cosmic distances.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The full-sky DM map and the 245 pc cm^-3 maximum rest on an untested low-latitude extrapolation of the disk-like halo fit; the practical 4% formula is therefore not yet supported in the plane.","rationale":"The paper is a useful, clearly written two-component model; the mean halo DM of 43 pc cm^-3 is a reasonable benchmark, and the authors are explicit in Section 5.2 that systematics dominate. The most load-bearing weakness is the extrapolation from the Suzaku footprint to the low-latitude sky. This is not an accusation of error: the model may be correct, but the data used do not constrain the region that generates the largest claimed corrections. The reader's weakest_assumption identified exactly this extrapolation, and I agree. A direct low-latitude X-ray comparison, or a sensitivity check that truncates or renormalizes the disk at |b| < 15, would settle whether the 245 pc cm^-3 end of the range is physical or an artifact of the assumed exponential profile. Because this concern affects the precision claim of Eq. (8) more than the headline mean, it supports the reader's CONDITIONAL verdict without moving it to rejection: the model is a valid benchmark but not yet a validated direction-dependent correction.","tokens_in":16610,"tokens_out":14187,"duration_ms":156780,"concrete_test":"Compare the model's predicted EM and DM against an independent all-sky X-ray map that includes low latitudes, e.g. the 3/4 keV band EM map of Henley & Shelton (2013) or eROSITA/eFEDS data, in longitude bins at |b| < 15 deg, especially l near 0 deg. Compute the ratio of predicted to observed EM in each bin after masking known foreground and point sources; if the mean ratio deviates from unity by more than the ~0.4 dex scatter of the Suzaku sample, the exponential disk extrapolation underlying the 245 pc cm^-3 peak and the low-latitude coefficients of Eq. (8) is not supported and should be replaced by a flagged extrapolation or a bounded upper limit.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing step is the extrapolation in Section 4 of the disk-like halo component fitted in Section 3.1 to N18 Suzaku EM data that cover only 75 deg < l < 285 deg and |b| > 15 deg. The headline range 30-245 pc cm^-3 is dominated by low-latitude sightlines, in particular l ~ 0 where the line of sight crosses the dense central part of the exponential disk twice (Figure 5). Those directions lie entirely outside the fitting region, so neither the in-plane normalization n0 nor the R0, z0 shapes in Eq. (2) are tested there. The paper's uncertainty discussion in Section 5.2 varies the spherical mass and integration limit, but it does not bound the low-latitude disk extrapolation; and the 0.2 dex DM scatter inferred from EM scatter at |b| > 15 cannot certify directions where the model predicts its largest corrections. This is a calibration gap rather than an internal inconsistency, so the mean 43 pc cm^-3 remains a plausible benchmark, but the claim that Eq. (8) estimates the halo DM along any line of sight to 4% is not yet supported where the correction matters most.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs a two-component model of the Milky Way's hot halo gas: a spherical, isothermal, hydrostatic component normalized to the cosmological baryon budget of a 10^12 solar mass halo, plus an exponential disk-like component whose central density and scale lengths are fitted by MCMC to 107 Suzaku X-ray emission-measure sightlines from Nakashima et al. (2018) covering 75 deg < l < 285 deg and |b| > 15 deg. The combined model is used to integrate electron density along lines of sight to the virial radius, yielding a mean halo DM of 43 pc cm^-3, a claimed full-sky range of 30-245 pc cm^-3, and a 0.2 dex rms fluctuation inherited from EM scatter. The authors provide a seventh-order polynomial formula (Eq. 8, Table 1) for DMhalo(l,b) and claim that it reproduces the model to better than 4% (better than 1% over 98% of the sky). The model is checked against LMC pulsar DMs and O VII column measurements, and applied to FRB 180924 and FRB 190523 to derive constraints on f_IGM.","tokens_in":16838,"tokens_out":6070,"duration_ms":63934,"significance":"If the calibration holds, this is a useful and needed ingredient for FRB cosmology: it converts a commonly ignored or spherically approximated foreground into a direction-dependent correction, with a simple fitting formula. The MCMC procedure is standard, priors and errors are stated, and the independent checks (LMC pulsars, O VII) are appropriate and go in the right direction; these are genuine strengths. The manuscript is also honest about systematic uncertainties in f_b and the integration limit. However, the quantitative headline claims, especially the 245 pc cm^-3 maximum and the any-line-of-sight validity of Eq. (8), rest on an extrapolation of the disk-like component into a region of parameter space (low |b|, l near 0) not sampled by the fitted EM data. As it stands, the mean DM benchmark is plausible, but the practical formula needs additional support before it can be used with the claimed accuracy.","major_comments":[{"comment":"The headline full-sky range DMhalo = 30-245 pc cm^-3 and the practical validity of Eq. (8) 'along any line of sight' depend on the disk-like component at low Galactic latitudes, where the model produces its largest values (e.g., the innermost curves in Figure 5). However, the disk parameters in Eq. (2) were fitted exclusively to N18 Suzaku EM sightlines with 75 deg < l < 285 deg and |b| > 15 deg (Section 3.1). The directions that dominate the upper end of the claimed range lie entirely outside that fitting region, and no independent low-latitude measurement (e.g., in-plane X-ray absorption, pulsar DMs, or other tracers) is used to validate the exponential extrapolation in R and z. Section 5.2 varies f_b and the integration limit but does not bound this low-latitude extrapolation. The result is a calibration gap, not an internal inconsistency: the mean 43 pc cm^-3 may remain a reasonable benchmark, but the claim that Eq. (8) estimates the halo DM along any line of sight, and in particular the 245 pc cm^-3 maximum, is not yet supported in the directions where the correction matters most.","section":"Section 4, Eq. (8), Figure 5"},{"comment":"The statement that Eq. (8) reproduces the theoretical prediction within 4% accuracy refers only to the accuracy of the polynomial approximation to the model, not to the accuracy of the model as a description of the Milky Way. The physical uncertainties quoted in Section 5.2 are much larger: 21-50 pc cm^-3 for the mean from f_b in [0,1], 14-26% from the integration limit, and a 0.2 dex rms scatter from EM fluctuations. Because Table 1 lists only coefficients and no uncertainty map, a user computing DMhalo with Eq. (8) has no way to propagate model uncertainty into the derived DMIGM or source redshift. The paper should either provide an uncertainty map for DMhalo(l,b) or explicitly restrict the '4% accuracy' claim to the polynomial representation of the fiducial model.","section":"Section 4, Table 1; Section 5.2"},{"comment":"The sentence stating that 'statistical uncertainties in the best-fit parameters of the disk-like halo component is negligible' is not demonstrated. The MCMC errors quoted in Section 3.1 are roughly 20-30% in n_disk0 and 10-20% in R0 and z0; these propagate into the disk DM, particularly along low-latitude sightlines where the integral passes through the high-density central region. No calculation is shown that the resulting DM uncertainty is negligible compared with the systematic effects discussed in Section 5.2, and in the extrapolated low-latitude region the parameter errors could be larger. This should be quantified if the claim of negligible statistical uncertainty is retained.","section":"Section 5.2"}],"minor_comments":[{"comment":"There are repeated typos: 'descirbed' should be 'described' and 'siteline' should be 'sightline' (also in Section 6).","section":"Section 3.1 and Conclusions"},{"comment":"The domain of l should be stated explicitly, since the polynomial in |l| and |b| is not manifestly periodic; in particular, define how l near 360 deg maps to |l|, as done in the Figure 1 caption.","section":"Equation (8)"},{"comment":"The notation EMN18,☉ is introduced somewhat abruptly after Eq. (5); a consistent subscript such as EM_N18,☉ would improve readability, and Eq. (5) would benefit from an explicit reminder that nH = chi_H n_e.","section":"Section 2, Eq. (5)"},{"comment":"The statement that larger integration limits r = 1.5rvir-2.0rvir increase the mean DMhalo 'only by 14%-26%' should specify whether this is relative to the fiducial mean of 43 pc cm^-3, to avoid ambiguity.","section":"Section 5.2"},{"comment":"The reference 'Pietrzyski' should be spelled 'Pietrzyński'.","section":"References"}],"recommendation":"major_revision","confidential_remarks":"The paper is potentially very useful, but the extrapolation gap is the key blocker. The authors can address it by restricting their claims to the fitted sky region, by adding a bounded uncertainty for the low-latitude extrapolation, or by validating with independent low-latitude data (e.g., inner-Galaxy pulsar DMs or X-ray absorption near the plane). I do not see a circularity problem: the halo DM is not used to derive the model parameters, and the MCMC fit to EM data with a baryon-budget normalization is a reasonable procedure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a genuinely useful paper for the FRB community, giving a direction-dependent Milky Way halo DM correction with a simple polynomial formula. The mean value of about 43 pc cm^-3 is a reasonable benchmark, but the 4% accuracy claim applies only to the polynomial's fit to the model, not to the model's reliability in the low-latitude directions where the correction is largest.\n\nWhat's actually new: the paper combines a disk-like hot halo component, fitted to 107 Suzaku emission-measure sightlines from Nakashima et al., with a spherical isothermal component normalized to the Milky Way's baryon budget. Neither piece is original, but the explicit two-component combination and the all-sky analytic formula in Eq. (8) are not in the cited literature. The fit itself is standard MCMC, with stated priors and errors, and the model is sensibly checked against LMC pulsar DMs and O VII absorption. The discussion of why the hot disk-like halo is unlikely to be double-counted in NE2001/YMW16 is plausible and worth having.\n\nThe soft spots are in the extrapolation, not the machinery. The Suzaku EM data cover only |b| > 15° and longitudes 75°–285°. The model then predicts its highest halo DMs (up to 245 pc cm^-3) in exactly the low-latitude, uncovered-longitude directions where the exponential disk is least constrained. That's a calibration gap, not an internal inconsistency, but it means the full-sky range is speculative in the regime that matters most for near-plane FRBs. There is no uncertainty map to accompany Eq. (8); the 4% is the polynomial's fitting error, not the model's. Also, the spherical component's consistency with the cosmic baryon budget is built in by construction, so it's not an independent check. The authors acknowledge some of this in Section 5.2, but they stop short of quantifying how wrong the disk extrapolation could be.\n\nWho this is for: FRB observers and people interpreting FRB DMs, especially those working on nearby low-DM bursts. The paper deserves a serious referee. I'd send it out, with the expectation that the authors add an uncertainty map, discuss the low-latitude gap more quantitatively, and soften the claim that Eq. (8) gives better than 4% accuracy for any line of sight. The central benchmark is probably fine; the precision claim is not yet supported.","headline":"Useful direction-dependent halo DM formula for FRB analyses, but the low-latitude extrapolation makes the 4% accuracy claim premature.","tokens_in":17417,"tokens_out":3017,"would_cite":true,"duration_ms":27524,"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":"The Milky Way's hot gas halo adds 30-245 pc/cm3 to FRB dispersion measures, and a new formula makes the subtraction practical.","keywords":["fast radio bursts","dispersion measure","Milky Way halo","hot gas halo","X-ray emission measure","dark matter halo","intergalactic medium","Galactic coordinates"],"falsifier":"Take a set of FRBs or pulsars with independent distance estimates toward low-latitude sightlines ($|b|<15^\\circ$), subtract the warm-ISM contribution using NE2001 and YMW16, and compare the residual with Eq. (8). If the residuals do not rise toward the plane as predicted, or systematically exceed the $30$–$245\\,{\\rm pc\\,cm^{-3}}$ envelope, the disk-like halo extrapolation is falsified. A complementary check is an all-sky X-ray emission-measure map covering the longitude gap and low latitudes, which should show the same exponential disk if the model is right.","tokens_in":16350,"feed_emoji":"🌌","tokens_out":9833,"duration_ms":99460,"temperature":0.7,"pith_summary":"This paper argues that the hot gas halo of the Milky Way is not a negligible or spherically symmetric correction for fast radio bursts. Its two-component model — a spherical, isothermal hot halo plus a compact disk-like hot component — predicts a mean halo dispersion measure of $43\\,{\\rm pc\\,cm^{-3}}$ and a full-sky range of $30$–$245\\,{\\rm pc\\,cm^{-3}}$, with the disk-like component dominating in most directions. Because the correction depends on where on the sky the burst is seen, the paper supplies a polynomial formula in Galactic coordinates that reproduces the model's halo dispersion measure to better than 4 percent. If the model holds, observers can subtract the Milky Way's halo contribution direction by direction and obtain cleaner estimates of the intergalactic and host-galaxy components of FRB dispersion.","feed_headline":"Milky Way halo gas adds up to 245 pc/cm3 to FRB dispersion","feed_subtitle":"A direction-dependent model lets observers subtract the Galaxy's halo when converting burst signals into distances.","key_machinery":"The load-bearing object is the two-component electron density profile $n_e = n_{\\mathrm{disk}} + n_{\\mathrm{sphe}}$. The disk-like component is an exponential disk, $n_{\\mathrm{disk}}(R,z) = n_0 \\exp[-(R/R_0 + |z|/z_0)]$, with best-fit $n_0 = 7.4\\times10^{-3}\\,(Z_{\\mathrm{halo}}/Z_\\odot)^{-1}\\,{\\rm cm^{-3}}$, $R_0 = 4.9\\,{\\rm kpc}$, and $z_0 = 2.4\\,{\\rm kpc}$, determined by fitting X-ray emission measures. The spherical component is isothermal gas at $kT = 0.3\\,{\\rm keV}$ in hydrostatic equilibrium with an NFW dark matter potential, normalized to a total baryon mass of $1.2\\times10^{11}\\,M_\\odot$. Integrating this density along any line of sight gives the halo DM sky map, which is then compressed into the analytic formula $\\mathrm{DM_{halo}} = \\sum_{i,j} c_{ij} |l|^i |b|^j$ with the coefficients listed in the paper's Table 1.","core_discovery":"The paper's central claim is that the directional variation seen in diffuse X-ray emission from the Milky Way's hot gas demands a disk-like halo component in addition to the extended spherical halo, and that the resulting electron density model predicts halo dispersion measures of $30$–$245\\,{\\rm pc\\,cm^{-3}}$ across the sky with a mean of $43\\,{\\rm pc\\,cm^{-3}}$. The disk-like component, fitted to X-ray emission measures, contributes between $0.4$ and $9$ times as much dispersion as the spherical component, so the halo DM is strongly non-isotropic. The model stays consistent with LMC pulsar dispersion measures after subtracting warm ISM models, and with O VII absorption column estimates, and the paper packages the result as a seventh-order polynomial in $|l|$ and $|b|$ with tabulated coefficients.","pith_inferences":["A statistical prediction worth testing: if the disk-like halo is real, FRB sightlines at low $|b|$ should show systematically larger residual DM after warm-ISM subtraction, and stacking FRBs by Galactic latitude could reveal the halo's signature without needing individual host redshifts.","The $0.2$ dex density fluctuation inferred from X-ray scatter suggests that a single polynomial value underweights sightline-to-sightline variance; future FRB samples should treat the halo DM as a distribution, not a point prediction.","The same two-component reasoning may apply to other galaxies, since intervening galaxy halos are omitted from the standard DM budget; if their halos resemble the Milky Way's, some of what is attributed to the IGM could actually be accumulated halo gas.","A decisive check is low-latitude X-ray spectroscopy: mapping emission measures at $|b|<15^\\circ$ would test whether the exponential extrapolation that produces the largest corrections is physically present."],"forward_implications":["FRB distance estimates can now include a direction-dependent Milky Way halo subtraction of typically $30$–$50\\,{\\rm pc\\,cm^{-3}}$, rising to $245\\,{\\rm pc\\,cm^{-3}}$ near the plane.","Adding Eq. (8) to the NE2001 or YMW16 warm-ISM models gives the total Milky Way electron contribution, isolating the intergalactic plus host-galaxy remainder.","For the host-identified FRB 180924 and FRB 190523, the lower halo DM weakens the upper bound on the ionized IGM fraction ($f_{\\rm IGM} < 0.79$–$0.96$ and $<0.99$–$1$ in the paper's estimates).","The scatter in X-ray emission measures implies a roughly $0.2$ dex rms fluctuation in halo DM, so the smooth formula should be read as a mean correction rather than an exact value for any single sightline.","Nearby, low-dispersion FRBs are where the halo term matters most, because the Milky Way's electrons can dominate their total observed DM."],"supporting_citations":[{"why":"supplies the 107-sightline diffuse X-ray emission measures fitted by the disk-like halo component.","marker":"N18"},{"why":"provides the Milky Way dark matter halo virial mass and radius used to normalize the spherical component.","marker":"Klypin et al. (2002)"},{"why":"defines the NFW profile that generates the gravitational potential in which the isothermal gas is in hydrostatic equilibrium.","marker":"Navarro et al. (1997)"},{"why":"supports the concentration value $c_{\\rm vir}=12$ adopted for the NFW halo.","marker":"Bullock et al. (2001)"},{"why":"supplies the fiducial Galactic baryon fraction $f_b = 0.75$ and serves as the high-halo comparison model.","marker":"Prochaska & Zheng (2019)"},{"why":"provides the NE2001 warm-ISM model used to subtract disk electron contributions from pulsar and FRB dispersion measures.","marker":"Cordes & Lazio (2002)"},{"why":"provides the YMW16 warm-ISM model used as the alternative subtraction in the same comparisons.","marker":"Yao et al. (2017)"},{"why":"supplies the adiabatic spherical-halo comparison model (F13) that the two-component model is contrasted with.","marker":"Fang et al. (2013)"},{"why":"supplies the multi-phase isothermal comparison model (F17) used in density and mass profile comparisons.","marker":"Faerman et al. (2017)"},{"why":"provides the cosmological simulation-based halo DM estimate (about $30\\,{\\rm pc\\,cm^{-3}}$) that the paper compares its mean to.","marker":"Dolag et al. (2015)"}],"fun_headline_variants":["Halo gas contributes up to 245 pc/cm3 to FRB dispersion","Halo DM varies 30–245 pc/cm3 by line of sight","New model: MW halo DM 30–245 pc/cm3 for FRBs","X-ray data fit: halo gas adds 30–245 pc/cm3 to FRB DM","Milky Way halo DM 30–245 pc/cm3 for FRB distances"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The disk-like exponential profile is fitted to X-ray data at $|b|>15^\\circ$ over a limited longitude range and then extrapolated to the entire sky, including the low-latitude directions where the model predicts its largest dispersion measures; if the hot gas distribution there differs from this extrapolation, the full-sky range and the fitting formula fail exactly in the directions that matter most.","fun_headline_variants_meta":{"raw":{"variants":["Halo gas contributes up to 245 pc/cm3 to FRB dispersion","Halo DM varies 30–245 pc/cm3 by line of sight","New model: MW halo DM 30–245 pc/cm3 for FRBs","X-ray data fit: halo gas adds 30–245 pc/cm3 to FRB DM","Milky Way halo DM 30–245 pc/cm3 for FRB distances"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000919,"raw_usage":{"total_tokens":3965,"prompt_tokens":991,"completion_tokens":2974,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":607,"completion_tokens_details":{"reasoning_tokens":2867}},"tokens_in":607,"tokens_out":2974,"duration_ms":21486,"temperature":1.0,"reasoning_tokens":2867,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:34:27.466029+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take a set of FRBs or pulsars with independent distance estimates toward low-latitude sightlines ($|b|<15^\\circ$), subtract the warm-ISM contribution using NE2001 and YMW16, and compare the residual with Eq. (8). If the residuals do not rise toward the plane as predicted, or systematically exceed the $30$–$245\\,{\\rm pc\\,cm^{-3}}$ envelope, the disk-like halo extrapolation is falsified. A complementary check is an all-sky X-ray emission-measure map covering the longitude gap and low latitudes, which should show the same exponential disk if the model is right.","supporting_citations":[],"review_version":1}