{"id":"5a4db533-19dc-4041-9ff1-091816009a1f","arxiv_id":"2608.10227","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Using ALMA's ACES survey, the authors produce two-dimensional maps of electron temperature and emission measure for HII regions across the Milky Way's central molecular zone, finding a mean electron temperature near 5900 K.","lead":"This paper maps the temperature and density of ionized gas around newborn massive stars across the Milky Way's central molecular zone using ALMA observations. It provides the first complete two-dimensional view of electron temperature and emission measure in this region, which is key to understanding star formation near the Galactic center.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Moment-0 clipping bias in the integrated H40α map can inflate Te and the CMZ mean; the single-Gaussian concern is secondary because Eq. 4 uses integrated line intensity, not peak width.","rationale":"The paper does a transparent job of applying a standard RRL method and its mean Te agrees with prior measurements, so ACCEPT-with-revision (CONDITIONAL) is the right posture. However, the Reader's stated weakest assumption (single Gaussian vs. multiple components affecting Δv) is not the most load-bearing: Eq. 3/4 depend on the integrated line intensity, so the line width cancels; multiple components would only matter through EM-weighted temperature blending, not through a Δv overestimate. The real soft spot is the threshold-clipped moment-0 map. The manuscript explicitly notes the clipping edge effect but asserts it does not affect Te/EM, which is an empirical claim that should have been tested. A quantitative comparison against full-cube integrated intensities would settle it. If the bias is small, the headline mean stands and the paper is acceptable with minor revision; if large, the mean and gradients are threshold artifacts. Since this is testable and the authors have the cube in hand, the reader's CONDITIONAL verdict remains appropriate.","tokens_in":23588,"tokens_out":26508,"duration_ms":242256,"concrete_test":"Take the H40α data cube and recompute the integrated line intensity without threshold clipping: after continuum subtraction, integrate each spectrum over the full velocity range that contains line emission (or over the fitted total of one/two Gaussian components), for a representative sample spanning low, medium, and high peak S/N, including the Sgr B2 Main gradient region. Re-run the TeEM equations with these corrected IL maps and compare the new <Te>_CMZ, the radial profile in Table 1, and the Sgr B2 Main cross-section gradient to the published values. If <Te> shifts by more than ~500 K or the gradient by more than ~1000 K pc^-1, the clipping bias is a dominant systematic; also tabulate ΔIL/IL versus peak S/N to show where the bias matters.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (maps of Te and EM, mean <Te>=5872 K) rests on the fidelity of the integrated H40α line intensity IL inserted into Eq. 4. The released moment-0 map is produced by clipping the cube above a threshold (Section 2.3.6). For faint emission or broad line wings, real flux is removed, so IL is underestimated and Te=(1.845e5)(IC/IL)^(1/1.15) is biased high at those pixels. The paper acknowledges that clipping produces sharp edges in Te but asserts without quantification that it 'does not affect' the Te/EM calculation. This is consequential: the mean is an average over many low-S/N pixels (Sgr C is excluded for exactly this reason), so a systematic IL deficit would raise the headline mean and could manufacture the claimed steep gradients (e.g., Sgr B2 Main ~4400 K pc^-1) if the threshold clips a spatially varying amount of line flux. The reader's identified mechanism (multiple velocity components inflating Δv) is less directly load-bearing because Eq. 3/4 use integrated line intensity, not the peak-to-integrated width; the vulnerability is missing flux in the clipped integral, not the Gaussian line-width assumption. A second-order effect is line-of-sight temperature mixing, where the IC/IL ratio is an EM-weighted blend and can be biased cool, but the clipping bias is the more concrete systematic.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the 'TeEM' method applied to ALMA ACES Band 3 data (99.6 GHz continuum, H40α peak intensity, and H40α integrated-intensity maps) to derive two-dimensional electron temperature and emission measure maps of HII regions in the CMZ. From standard radio recombination line and free-free continuum relations (Eqs. 3 and 8), it obtains a CMZ-wide mean <Te>_CMZ = 5872 ± 78 (SE) ± 3682 (SD) K, EM values from ~10^5 to ~3×10^8 pc cm^-6, a radial Te profile, and region-by-region measurements for Sgr B2, Sgr B1, the Sickle, the Pistol, the Bridges, Sgr A HII regions, and the Minispiral. Sgr C is excluded because of insufficient line signal-to-noise. The central claim is that TeEM applied directly to 2D moment maps yields reliable Te and EM maps and a relatively uniform CMZ electron temperature near 5900 K, with steep local gradients in Sgr B2 Main and the Minispiral.","tokens_in":23811,"tokens_out":7484,"duration_ms":78548,"significance":"The strength of the paper is its simplicity and transparency: it applies well-established RRL/free-free formulas (Quireza et al. 2006; Oster 1961) to a public, high-resolution survey, and the resulting mean Te is consistent with independent VLA and ALMA measurements of individual regions. If the maps are validated, they would provide the first systematic ~0.1 pc resolution Te and EM maps across the CMZ, a useful benchmark for star-formation-rate and metallicity-gradient studies. The paper also makes falsifiable predictions, notably the steep Te gradient in Sgr B2 Main and the very high Te in the Minispiral's north-south arm. The main weakness is that several systematic effects are acknowledged but not quantified, so the headline mean and gradients must be treated as conditional on the adopted masking and integration thresholds.","major_comments":[{"comment":"The statement that clipping of the moment-0 map 'does not affect the calculation of Te and EM' is not supported and is likely incorrect. Because Te = 1.845×10^5 (IC/IL)^(1/1.15), any H40α flux removed by the intensity threshold lowers IL and therefore raises Te. The sharp edges visible in the Te maps are direct evidence that the clipped moment-0 map is entering Eq. (4). The pixels most affected are low-S/N pixels, which also dominate the sample and are the stated reason for excluding Sgr C. I request a quantitative test: recompute IL by directly integrating the line cube with a range of thresholds (or down to the noise floor) and report the resulting variation in <Te>_CMZ and in the Sgr B2 Main gradient. Without this test, the headline mean may be systematically biased high.","section":"§2.3.6, Eq. (4)"},{"comment":"All quoted uncertainties are statistical only (SE or SD). The headline mean and the per-parsec gradient claims need a systematic error budget that includes at least the moment-0 clipping threshold (§2.3.6), the dust-mask threshold (§2.3.3), the continuum zero-level and beam mismatch between the 2.45'' line beam and 2.14'' continuum beam, molecular-line contamination (§2.3.5), and the fixed assumptions [He]/[H]=0.07 and a=0.9 in Eq. (3). As written, the reported ±78 K standard error gives a false impression of the reliability of <Te>_CMZ; the systematic uncertainties are likely much larger than the statistical ones.","section":"§3.3, Tables 1 and 2"},{"comment":"The dust-excess mask is applied at a spectral index threshold of α > +2, but no sensitivity test is shown. Since dust contamination raises IC, it biases Te high wherever the mask is imperfect. The comparison of panels F and G in Fig. 5 shows non-negligible changes in Sgr B2, and the text says 'some differences are found in peaky strong continuum sources.' Please report the number of masked pixels and repeat the <Te>_CMZ calculation for a range of thresholds (e.g., α = 1, 2, 3), or replace the ad hoc mask with a direct free-free/dust decomposition using the 86.6 and 99.6 GHz maps.","section":"§2.3.3 and Fig. 5"},{"comment":"The single-Gaussian argument in §2.3.1 is not the right justification for Eq. (4), which uses the integrated line intensity and does not require a Gaussian line shape. The more relevant concern is line-of-sight superposition: if two HII components with different Te fall within the beam, IC/IL is an EM-weighted mixture that is biased toward dense, cool gas. The few profiles shown in Fig. 3 do not establish that this effect is negligible over the whole field. I recommend a validation in which a sample of spectra (including Sgr B2 Main and the Minispiral) are fit with single- and multi-component models, and the TeEM results are compared with the fit-based Te values; the expected bias direction and magnitude should be stated.","section":"§2.3.1 and Eq. (3)"}],"minor_comments":[{"comment":"The panel labels are inconsistent: the Fig. 1 caption and §2.1 say Panels B, C, D are continuum, integrated intensity, and peak intensity, while §3.1 describes Panel A as continuum and Panel B as peak intensity. Please unify the numbering.","section":"Fig. 1, §2.1, §3.1"},{"comment":"The citation 'Rohlsfs and Wilson 2000' is a typo; it should be 'Rohlfs and Wilson 2000.'","section":"Eq. (3)"},{"comment":"The text refers to both the 'NS arm' and the 'SN arm' of the Minispiral; the latter is the correct expression for the north-south arm. Please use one nomenclature consistently.","section":"§4.6"},{"comment":"The stated reason for excluding Sgr C—that the integrated intensity is too weak to calculate the velocity width by dividing by peak intensity—is confusing, because Eq. (4) uses the integrated intensity directly. Please rephrase to say that the H40α integrated-intensity map has no significant detection above the clipping threshold in the Sgr C field.","section":"§4.8"},{"comment":"The Sgr B2 Main gradient is quoted as ~3000 K pc^-1 in §4.1 and as ~4400 K pc^-1 in §5.3. Make the values consistent or state explicitly that they are measured along different cuts or with different averaging.","section":"§4.1 and §5.3"},{"comment":"The electron densities in Table 2 use L ~ 2r with no uncertainty or discussion of the geometry; please state explicitly that n_e is an order-of-magnitude estimate, or propagate a plausible range of L.","section":"§2.6 and Table 2"}],"recommendation":"major_revision","confidential_remarks":"The paper is a useful application of a standard method to a public survey, and the core derivation is sound, but the systematic biases—especially the moment-0 clipping effect—need to be quantified before the headline mean can be taken at face value. The authors should also check that all companion ACES papers are cited in their final published versions, as several are listed as submitted or arXiv preprints."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague, this is the first paper to push integrated-intensity-based Te mapping across the whole ACES CMZ footprint, and the resulting Te and EM maps are genuinely new data products. The method is standard radio recombination line physics (Quireza et al. 2006) packaged as a simple algorithmic shortcut — no per-pixel fitting — and applied to good, publicly available ALMA data. The headline mean of ~5900 K is consistent with prior measurements, and the paper is transparent about its assumptions, even listing the ones it is uneasy about. That is real value, especially as a guide to where HII gas sits in the CMZ. I would not want to lose the maps.\n\nThe soft spots are concentrated in one place: the moment-0 clipping. Section 2.3.6 states that the threshold used to make the integrated H40α map creates sharp edges in Te but 'does not affect the calculation.' That is not argued, and it matters. If the clipping removes real line flux — broad wings, faint extended emission — then IL is too small and Eq. (4) pushes Te up. Almost all of the pixels entering the global mean are low-S/N pixels; Sgr C is excluded for exactly that reason. The stress-test note has it right: the single-Gaussian concern is secondary because Eqs. (3)-(4) use integrated line intensity, not the peak-to-width relation. The real vulnerability is missing flux in the clipped integral. I would like a comparison between TeEM and per-pixel Gaussian fits on a sample of sight lines, plus a recomputation with the clipping threshold varied, before trusting the headline number or the quoted steep gradients. The SE of 78 K on the mean is also not the right error bar; the pixel-to-pixel SD of 3682 K shows how much structure there is, but neither captures the systematic. The dust mask threshold α>+2 is ad hoc, though the authors do show the unmasked map for Sgr B2, which is the right instinct. Sgr C exclusion is clearly stated and defensible, but it means the map is not the whole CMZ.\n\nNone of this kills the paper. The central derivation is standard, the comparison to earlier measurements is honest, and the maps have obvious survey utility. Citation practice looks appropriate: the heavy ACES citations are to the survey papers whose data products are being used. I would send this to a referee. The right outcome is probably publication after a revision that quantifies the clipping bias and adds a spot-check against spectral fitting. Who benefits: anyone working on Galactic Center star formation, radio recombination line methods, or the ACES survey. I would bring it to a reading group.","headline":"Useful first CMZ-wide Te/EM maps from a simple method, but the headline mean is hostage to an unquantified moment-0 clipping bias; referee with a request for systematics.","tokens_in":24650,"tokens_out":3045,"would_cite":true,"duration_ms":27694,"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 paper maps electron temperature and emission measure across the Galactic center's HII regions, finding a nearly uniform mean temperature of 5,872 K and emission measures spanning three orders of magnitude.","keywords":["Galaxy: center","HII regions","ISM: lines and bands","ISM: abundance","radio continuum: ISM","radio lines: ISM","stars: formation","electron temperature"],"falsifier":"Fit Gaussian or multi-Gaussian profiles to the H40α data cubes in Sgr B2 Main, the Minispiral, and Sgr B1 and compare the per-pixel line widths and Te values with the TeEM maps; if multi-component fits yield systematically higher Te, the single-Gaussian assumption would be falsified.","tokens_in":23350,"feed_emoji":"🔭","tokens_out":9192,"duration_ms":73347,"temperature":0.7,"pith_summary":"By applying the TeEM method to ALMA's ACES survey of the Central Molecular Zone, the paper derives two-dimensional maps of electron temperature and emission measure for the entire population of HII regions in the Galactic center at about 0.1 parsec resolution. It claims a remarkably uniform mean electron temperature of 5,872 ± 78 (standard error) K across the CMZ, with most regions between 4,000 and 8,000 K. Emission measure, by contrast, varies by orders of magnitude, from roughly $10^{5}$ to 3×$10^{8}$ pc $cm^{-6}$. The maps reveal steep temperature gradients of several thousand kelvin per parsec inside sources like Sgr B2 Main and the Minispiral, and an east-west asymmetry of star-forming regions about Sgr A*. If correct, this gives an extinction-free census of ionized gas in the Milky Way's nuclear region and a direct measurement of the physical conditions of star formation there.","feed_headline":"ALMA maps: Galactic center HII gas sits near 5,900 K","feed_subtitle":"Full survey of the CMZ yields first 0.1-pc maps of electron temperature and emission measure","key_machinery":"The machinery is the TeEM ('Te-EM mapping') algorithm, which converts three 2D ACES maps — continuum intensity IC, H40α peak intensity, and integrated intensity IL — into maps of electron temperature Te and emission measure EM. The line width is approximated as Δv = IL / IL_peak, and the electron temperature follows from the standard LTE recombination-line relation, Te/K ≈ 1.845×$10^{5}$ (IC(νL)/IL)^(1/1.15) (Equation 4), with EM then derived from the free-free optical depth (Equation 8). The method's practical power is that it needs only moment maps, skipping spectral cube fitting; its practical safeguard is a dust-contamination mask that drops pixels whose 99.6/86.6 GHz continuum spectral index exceeds +2.","core_discovery":"The paper's central claim is that the electron temperature and emission measure of HII regions throughout the CMZ can be read directly from three two-dimensional ACES maps — 99.6 GHz continuum intensity, H40α peak intensity, and H40α integrated intensity — using the TeEM algorithm, without fitting individual spectra. From these, the paper derives a grand mean temperature for the CMZ of 5872 ± 78 (SE) ± 3682 (SD) K, finds that the temperature is nearly flat across the region at the ~6000 K level, and finds emission measures spanning $10^{5}$ to 3×$10^{8}$ pc $cm^{-6}$. It further reports local temperature gradients as steep as ~4400 K $pc^{-1}$ in Sgr B2 Main and ~15,000 K in the north-south arm of the Minispiral nearest Sgr A*, interpreting these as internal structure and possible excitation by the central source rather than metallicity variations. The paper also reports that HII regions are strongly concentrated at positive longitudes relative to Sgr A*, indicating asymmetric current star formation in the CMZ.","pith_inferences":["If the ~5900 K mean holds up, any CMZ HII region with a Te measurement far above or below this value becomes a candidate for exotic heating or cooling (e.g., X-ray or shock ionization), and one could run that filter automatically over the maps.","A natural stress test of the method is to compare TeEM temperatures with per-pixel Gaussian fitting of the same cube in Sgr B2 Main and the Minispiral; if the fitted widths are systematically narrower, the claimed internal gradients would shrink, pointing to line blending rather than physical temperature structure.","The same two-frequency dust separation used here could be combined with lower-frequency data to construct a synchrotron-corrected EM map, refining densities in the Radio Arc region."],"forward_implications":["The CMZ's HII regions are essentially isothermal at ~5900 K, so average electron temperature is not a strong function of environment within the central ~100 pc.","Emission measure maps, combined with assumed line-of-sight depths, give electron densities ranging from ~10^3 to ~4×10^4 cm^-3, tracing the densest ionized clumps.","Combined with Galactic-disk Te data, the CMZ value anchors a temperature gradient of ~300 K kpc^-1 from the center outward, which can be read as a metallicity gradient.","The TeEM method can be applied to any interferometric recombination-line plus continuum survey to produce Te and EM maps without per-spectrum fitting, so it may become a standard processing step for such data."],"supporting_citations":[{"why":"Provides the ACES survey mosaics from which the continuum and line maps are drawn.","marker":"Longmore et al. 2026"},{"why":"Supplies the 99.6 GHz continuum map used as the free-free tracer.","marker":"Ginsburg et al. 2025"},{"why":"Supplies the H40α line peak and integrated intensity maps.","marker":"Hsieh et al. 2026"},{"why":"Gives the electron-temperature formula that converts the continuum-to-line ratio into Te.","marker":"Quireza et al. 2006"},{"why":"Establishes the LTE recombination-line/continuum relation underlying the method.","marker":"Mezger & Hoglund 1967"},{"why":"Demonstrates the Te measurement technique at cm/mm wavelengths in the CMZ that the paper extends to full 2D maps.","marker":"Tsuboi et al. 2017"},{"why":"Quantifies dust emission at mm wavelengths, motivating the spectral-index mask.","marker":"Schmiedeke et al. 2016"},{"why":"Provides earlier Sgr B2 electron-temperature measurements used as consistency checks.","marker":"Mehringer et al. 1993"},{"why":"Supplies the 1.3 GHz MeerKAT continuum map used to show synchrotron contamination is small at 99.6 GHz.","marker":"Heywood et al. 2022"}],"fun_headline_variants":["CMZ ionized gas: ACES maps reveal ~5,900 K temperature","Galactic center HII regions mapped: mean Te ~5,900 K, EM up to 3e8","ACES: Te-EM maps of CMZ HII regions, mean 5,872 K","First 0.1-pc Te-EM maps of Galactic center HII gas"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result hinges on the assumption that every sight line's H40α line is a single simple Gaussian, so that the integrated-to-peak intensity ratio reliably measures the line width; if multiple velocity components blend together, the width is overestimated and the temperature underestimated.","fun_headline_variants_meta":{"raw":{"variants":["CMZ ionized gas: ACES maps reveal ~5,900 K temperature","Galactic center HII regions mapped: mean Te ~5,900 K, EM up to 3e8","ACES: Te-EM maps of CMZ HII regions, mean 5,872 K","First 0.1-pc Te-EM maps of Galactic center HII gas"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000624,"raw_usage":{"total_tokens":3038,"prompt_tokens":1243,"completion_tokens":1795,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":859,"completion_tokens_details":{"reasoning_tokens":1696}},"tokens_in":859,"tokens_out":1795,"duration_ms":11132,"temperature":1.0,"reasoning_tokens":1696,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T04:10:20.623062+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit Gaussian or multi-Gaussian profiles to the H40α data cubes in Sgr B2 Main, the Minispiral, and Sgr B1 and compare the per-pixel line widths and Te values with the TeEM maps; if multi-component fits yield systematically higher Te, the single-Gaussian assumption would be falsified.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the LTE recombination-line/continuum relation underlying the method."},{"cited_title":"T., Bania, T","cited_arxiv_id":null,"evidence_quote":"Gives the electron-temperature formula that converts the continuum-to-line ratio into Te."},{"cited_title":"2017, ApJ, 842, 2,","cited_arxiv_id":null,"evidence_quote":"Demonstrates the Te measurement technique at cm/mm wavelengths in the CMZ that the paper extends to full 2D maps."},{"cited_title":"2016, A&A, 588, A143","cited_arxiv_id":null,"evidence_quote":"Quantifies dust emission at mm wavelengths, motivating the spectral-index mask."},{"cited_title":"M., Palmer, P., Goss, W","cited_arxiv_id":null,"evidence_quote":"Provides earlier Sgr B2 electron-temperature measurements used as consistency checks."},{"cited_title":"2022, ApJ, 925,","cited_arxiv_id":null,"evidence_quote":"Supplies the 1.3 GHz MeerKAT continuum map used to show synchrotron contamination is small at 99.6 GHz."}],"review_version":1}