{"id":"f65e01fc-2517-41dd-9781-9e2591f0e5f8","arxiv_id":"2607.21432","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":6,"one_line_summary":"Time-domain template subtraction using dedicated 183 GHz-line-wing monitor detectors suppresses simulated atmospheric CMB contamination about 4,000 times better than filter-bin map-making at ℓ=30–300.","lead":"Ground-based CMB telescopes lose the largest angular scales to water-vapor noise in the atmosphere. This paper proposes dedicating focal-plane detectors to monitoring the atmosphere and subtracting fitted templates in the time domain — in simulations the approach suppresses atmospheric power about 4,000 times better than standard polynomial filtering, but only for an atmosphere model the method assumes by construction.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. 7 as written has a sign inconsistency: the GLS weights from Eq. 6 applied to zero-padded templates produce the negative of the optimal atmosphere predictor, so Eq. 9 would add rather than subtract atmosphere. The manuscript needs a sign correction or a revised derivation.","rationale":"The reader's weakest assumption (simulated atmosphere exactly matches the pipeline's model) is a legitimate scope limitation and the authors concede it. However, the more immediately load-bearing issue is internal: Eq. 7 as written is inconsistent with the GLS derivation. Since the reported simulation results show good cleaning, the code almost certainly uses the opposite sign, but the manuscript must be corrected. This does not overturn the central claim if the code is right, so the verdict remains CONDITIONAL with an additional mandatory correction. I would also encourage the model-mismatch tests the reader requested, but the sign inconsistency should be fixed first.","tokens_in":12574,"tokens_out":13449,"duration_ms":135712,"concrete_test":"Re-derive Eq. 7 from Eq. 6 for a single monitor template: let C = [[Var(d), Cov(d,z)], [Cov(d,z), Var(z)]]. Compute w^T and w^T [0; z]. If the result is the negative of the conditional mean Cov(d,z)/Var(z) * z, the sign error is confirmed. Alternatively, run the published pipeline as written (with no sign flip) on a simulation; if the cleaned PSD is worse than the raw PSD, the equation is inconsistent with the claimed results.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section 3.2.2 defines the interpolation weights in Eq. 6 as the GLS estimator of the sky signal in Eq. 4: w^T = (M^T C^{-1} M)^{-1} M^T C^{-1}. For a block covariance C with first entry the detector data and remaining entries the seven monitor templates, the row vector is [1, -C_{12}C_{22}^{-1}]. Applying these weights to the zero-padded template vector in Eq. 7 gives \\hat{s}^atm_i = -C_{12}C_{22}^{-1} [\\hat{s}^atm_1 ...]^T, the negative of the minimum-variance prediction of the detector atmosphere from the monitor templates. The correct interpolated template should be +C_{12}C_{22}^{-1} z (or equivalently d_i - w^T y_i). As written, subtracting this in Eq. 9 would inject correlated atmosphere into the cleaned TOD. The reported residual power spectra are therefore inconsistent with the equations unless the code uses the opposite sign, which would mean Eq. 7 is a typo. This is load-bearing because the central claim depends on the template-subtraction step; the manuscript cannot be reproduced as written.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a time-domain, multi-frequency atmosphere mitigation pipeline for ground-based CMB intensity observations. A 37-pixel focal plane includes seven 'tetrachroic' atmosphere-monitor pixels observing two extra bands in the wings of the 183 GHz water line (179 and 188 GHz). The pipeline (i) performs generalized least-squares component separation on the monitor timelines to form atmosphere templates, (ii) interpolates these templates onto each science detector using data–template covariance weights, (iii) Wiener-filters the interpolated template to avoid injecting monitor white noise, and (iv) subtracts the filtered template from the science TOD. Using TOAST simulations with MERRA-2-calibrated weather and a frozen-flow single-component atmosphere, plus a simplified CMB+dust sky, the authors compare the method with order-20 filter-bin polynomial cleaning. They define an atmospheric decontamination factor A^atm_ell and report an average of 1.3e-5 over 30<=ell<=300 for template subtraction versus 5.6e-2 for filter-bin, a factor of about 4000, and state that the residual scales as the inverse survey time.","tokens_in":12876,"tokens_out":7700,"duration_ms":83774,"significance":"If correct, the method offers a genuine alternative to filter-bin and maximum-likelihood map-making for large-scale CMB intensity: it avoids the transfer-function loss of polynomial filtering and recovers sky modes that filter-bin removes. The simulation setup is concrete and reproducible in its use of TOAST, MERRA-2 weather statistics, and the ATM atmospheric model, and the new metric A^atm_ell is a useful way to quantify residual atmospheric contamination. The genuinely novel content is the time-domain architecture: dedicated line-wing monitor bands, GLS component separation, covariance-weighted template interpolation, and Wiener filtering before subtraction. However, the demonstration is a matched test: the simulated atmosphere is exactly the single-component, constant-scaling model the pipeline assumes, and the sky mixing matrix is taken as known. The quantitative 4000x factor should therefore be read as a proof-of-concept result for idealized conditions, not as a demonstrated robustness property for real atmospheric complexity.","major_comments":[{"comment":"There is a sign inconsistency in the derivation of the interpolated template. For the block covariance C partitioned with first entry the detector data and remaining entries the seven monitor templates, the GLS weights in Eq. (6) are w_i^T = [1, -C_{12} C_{22}^{-1}]. Applying these weights to the zero-padded template vector [0, s^atm_1, ..., s^atm_7]^T in Eq. (7) gives s^atm_i = -C_{12} C_{22}^{-1} [s^atm_1 ...]^T, the negative of the minimum-variance prediction of the detector atmosphere from the monitors. The correct interpolated template is +C_{12} C_{22}^{-1} z, equivalently d_i - w_i^T y_i. As written, Eq. (9) subtracts the negative template and therefore adds correlated atmosphere to the cleaned TOD. This is load-bearing because the central decontamination claim depends on this subtraction step; the manuscript cannot be reproduced as written unless the code uses the opposite sign,","section":"Section 3.2.2, Eqs. (6)-(7)"},{"comment":"The headline result is a matched-test demonstration. The atmosphere in Section 2.3 is a single frozen-flow 3D field converted to detector signal by a per-band absorption scaling; the pipeline estimates that same constant mixing vector G by regressing slopes between TODs and projects it out, with the sky mixing matrix F assumed known. The manuscript itself concedes in Section 5 that the mixing vector is 'exact by construction' in the simulation. The reported A^atm = 1.3e-5 therefore measures how well the pipeline inverts its own generative model, not how it behaves under the violations most relevant for real data: multiple atmosphere layers at different velocities, PWV-dependent emission near the 183 GHz line, or mis-estimated F. Because the factor-4000 claim is the central result, the authors should add a robustness test, e.g. two atmosphere components with different wind velocities, a P","section":"Section 3.2.1 and Section 5"}],"minor_comments":[{"comment":"The definition of A^atm_ell contains a typographical 'B' where '≡' or '=' is intended. The formula is garbled as printed.","section":"Equation (13)"},{"comment":"The symbol N in Eq. (3) is introduced earlier as a per-sample diagonal cross-frequency noise covariance, but this should be stated explicitly at Eq. (3). Also clarify whether the covariance C_i in Eq. (5) is computed over the full one-hour scan or separately for each scan direction, since the text says the scan directions are treated separately.","section":"Eq. (3) and Eq. (5)"},{"comment":"The claim that A^atm scales as the inverse of the survey time is presented without error bars or a fitted scaling law. It is an expected consequence of coadding uncorrelated weather realizations, but it would strengthen the paper to state explicitly whether this is an analytic prediction and to show the A^atm values for the 10/20/40-hour cases.","section":"Section 5 and Fig. 10"},{"comment":"The captions use '10 hours, 20 hours, 40 hours' to describe integration time, but the survey consists of four 10-hour seasons combined with inverse-variance weighting in Eq. (12). Please define the integration-time convention and clarify how the 20-hour and 40-hour cases are formed from the four seasons.","section":"Figs. 8-10"},{"comment":"The manuscript contains no code or data availability statement. Since the paper is a methods proof of concept built on TOAST and public sky models, releasing the pipeline scripts and TOAST configuration would substantially improve reproducibility.","section":"Code availability"}],"recommendation":"major_revision","confidential_remarks":"The sign error in Eq. (7) is the main blocker: either the equations are wrong or the code uses a different sign, and the reported residual spectra cannot be reproduced from the text as written. The matched-simulation limitation is also important but is acknowledged by the authors; a robustness test with a more realistic atmosphere model would materially raise confidence. If the sign is corrected and a robustness test added, the paper would be a solid proof-of-concept for a new time-domain atmosphere mitigation architecture."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Julien and colleagues have an idea worth taking seriously: dedicate a handful of tetrachroic pixels to atmosphere monitoring, do GLS component separation on those time streams, interpolate the atmosphere template across the focal plane, and Wiener-filter before subtraction. The proof of concept on TOAST/MERRA-2 simulations is credible, the statistics are standard, and the authors are explicit about their assumptions. The A^atm metric is useful, and the claimed ~4000× improvement over order-20 filter-bin in the matched simulation is plausible if the pipeline works.\n\nBut there is a load-bearing problem in the written equations. Eq. 6 gives the GLS weights for estimating the sky signal from the stacked vector y_i = (d_i, templates). With M = (1,0,...,0), those weights come out to (1, -C12 C22^{-1}). Applying them to the zero-padded vector (0, templates) in Eq. 7 yields the negative of the minimum-variance atmosphere prediction, not the prediction itself. Eq. 9 then subtracts that from the data, which would inject correlated atmosphere into the cleaned TOD. The reported results cannot be reproduced from the equations as written. If the code uses the opposite sign, Eq. 7 needs a sign correction or a revised derivation. This is not cosmetic; it sits at the center of the paper's claim.\n\nThe other weaknesses are real but secondary. The simulation is a matched filter in two senses: the atmosphere is generated as a single component with a constant per-band scaling (exactly what the pipeline assumes), and the sky mixing matrix F is taken as known. The authors concede this in Section 5, but the headline residual of 1.3e-5 is conditional on it. Only four realizations are run, with no mismatch tests (multi-layer wind, PWV-dependent line shape, mis-estimated foreground SEDs), and no polarization. The 5 Hz decorrelation excess is a hint that the real atmosphere will be messier.\n\nStill, the idea deserves a serious referee. If the sign error is a typo, and if the authors add robustness tests and release the pipeline code, this could be a useful contribution to ground-based CMB analysis. Send it to review, with the expectation of heavy revision.","headline":"Genuinely new idea for atmosphere subtraction, but the written equations have a sign error that breaks the central claim; worth referee time only if the correction and robustness tests are forthcoming.","tokens_in":13474,"tokens_out":5586,"would_cite":false,"duration_ms":51527,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"By treating the atmosphere as a time-domain component, this paper recovers CMB maps with 4000 times less atmospheric contamination than filtering.","keywords":["cosmic microwave background","atmospheric contamination","time-domain component separation","multichroic focal plane","water-vapor line","atmospheric decontamination factor","map-making","large-scale CMB modes"],"falsifier":"Run the same pipeline on a simulation with two independent atmospheric layers moving at different velocities (or with a water-vapor-dependent emission law across the 183 GHz line) and recompute A^atm averaged over 30<=ell<=300; if it rises much above 10^-3 or the residual excess near 5 Hz grows, the single-component template assumption is the weak point. A complementary empirical check is to apply the method to real multi-frequency focal-plane data and compare the cleaned large-scale map spectrum against an independent instrument.","tokens_in":12317,"feed_emoji":"🔭","tokens_out":11389,"duration_ms":96906,"temperature":0.7,"pith_summary":"The paper sets out to replace low-frequency filtering in ground-based CMB analysis with a time-domain cleaning step. The idea is to give a few focal-plane pixels extra frequency channels that straddle the 183 GHz water-vapor line, use generalized least-squares component separation on those monitor pixels to build atmosphere-only templates, interpolate the templates to every science detector using the measured covariance between detectors and monitors, and Wiener-filter before subtracting. On simulations with realistic weather and a 40-hour survey, the method leaves only 1.3e-5 of the atmospheric angular power over multipoles 30-300, against 5.6e-2 for a filter-bin pipeline, and the residual falls as the inverse of survey time. If this carries to real data, ground-based experiments could recover large-scale intensity modes that filtering currently discards, with no transfer-function correction.","feed_headline":"Atmosphere templates beat CMB filtering by 4000x in simulations","feed_subtitle":"Dedicated water-line monitors build templates that remove atmosphere without the transfer-function loss of filtering.","key_machinery":"The load-bearing object is the time-domain atmosphere template. Seven tetrachroic monitor pixels observe the science bands plus two narrow channels at 179 and 188 GHz that straddle the 183 GHz water-vapor emission line. Generalized least-squares component separation with the sky mixing matrix taken as known and the atmosphere mixing vector fitted as the slope between frequency time streams projects the atmosphere out of the monitor data. The resulting templates are interpolated to each science detector with minimum-variance weights built from the data-template covariance matrix, then Wiener-filtered so monitor white noise is not injected during subtraction. The atmospheric decontamination fa","core_discovery":"The central claim, demonstrated in simulation, is that atmosphere contamination can be treated as an additional time-domain component and separated from sky signal before map-making, rather than removed by filtering that also removes sky signal. The authors define the atmospheric decontamination factor A^atm_ell as the residual atmosphere power in the output map (after white-noise subtraction and transfer-function correction) divided by the binned atmosphere power. They report a factor of about 4000 improvement over filter-bin, with A^atm averaged over 30<=ell<=300 equal to 1.3e-5 for template subtraction and 5.6e-2 for filter-bin. They also find that the residual scales as the inverse of su","pith_inferences":["Editorial extension: the factor-of-4000 is computed for a single frozen-flow atmosphere with one constant mixing vector; a natural stress test is a two-layer atmosphere with different wind velocities, which would reveal how much of the advantage survives real atmospheric structure.","Editorial extension: the 179/188 GHz monitor channels could double as a continuous precipitable-water-vapor diagnostic, potentially replacing part of the weather-cut bookkeeping with in-situ calibration.","Editorial extension: because the sky mixing matrix is assumed known in the paper, a full pipeline on real data would need an outer iteration that re-estimates the sky components from first-pass maps; the paper gestures at this but does not demonstrate convergence.","Editorial extension: the method's transfer-function-free property suggests it could be combined with maximum-likelihood map-making rather than inverse-variance binning to propagate the correlated residuals properly, a step the paper lists as future work."],"forward_implications":["Large-scale CMB intensity modes can in principle be recovered from the ground without the transfer-function loss that polynomial filtering imposes; the cleaned maps have T_ell = 1 by construction.","Longer surveys directly improve large-scale sensitivity: residual atmosphere in the cleaned maps scales as the inverse of total survey time.","Only a small set of dedicated monitor pixels (7 out of 37 in the toy focal plane) is sufficient to clean all science detectors in a 1-degree field of view in this simulation.","Because the pipeline is linear and each contribution can be propagated separately, the same simulation framework can isolate how much of the residual comes from atmosphere, noise, or sky-model error.","The method applies to intensity; reducing unpolarized atmosphere in the timeline also lowers the intensity-to-polarization leakage in polarization measurements, though quantifying that requires an instrument model."],"fun_headline_variants":["CMB atmosphere removal 4000x better with template separation","Dedicated atmosphere detectors boost CMB sensitivity 4000x","Time-domain component separation cuts CMB atmosphere 4000x","Atmosphere templates beat filter-bin in CMB by 4000x"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The result rests on the assumption that the atmosphere seen by the monitor and science detectors is exactly one component with a fixed frequency scaling that the pipeline knows and inverts perfectly, plus a known sky mixing matrix; if the real atmosphere has multiple layers, a frequency-dependent emission law, or a mis-estimated sky signal, the template projection leaks and the 4000-fold gain shrinks.","fun_headline_variants_meta":{"raw":{"variants":["CMB atmosphere removal 4000x better with template separation","Dedicated atmosphere detectors boost CMB sensitivity 4000x","Time-domain component separation cuts CMB atmosphere 4000x","Atmosphere templates beat filter-bin in CMB by 4000x"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000202,"raw_usage":{"total_tokens":1224,"prompt_tokens":752,"completion_tokens":472,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":496,"completion_tokens_details":{"reasoning_tokens":398}},"tokens_in":496,"tokens_out":472,"duration_ms":5239,"temperature":1.0,"reasoning_tokens":398,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T07:29:01.406713+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same pipeline on a simulation with two independent atmospheric layers moving at different velocities (or with a water-vapor-dependent emission law across the 183 GHz line) and recompute A^atm averaged over 30<=ell<=300; if it rises much above 10^-3 or the residual excess near 5 Hz grows, the single-component template assumption is the weak point. A complementary empirical check is to apply the method to real multi-frequency focal-plane data and compare the cleaned large-scale map spectrum against an independent instrument.","supporting_citations":[],"review_version":1}