{"id":"d0e0f5e3-fef7-4c5c-b4bc-861bd6c7b1d7","arxiv_id":"2610.31175","paper_version":1,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Masking radio sources above about 0.1 Jy and adding one power-law residual component lets ground-based CMB experiments recover unbiased B-mode parameters.","lead":"This paper simulates how polarized radio galaxies and dusty galaxies contaminate CMB B-mode searches, then tests masking and modeling strategies to remove the contamination. It finds that masking radio sources above about 0.1 Jy and adding one power-law component cleans the bias at negligible cost to parameter precision.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"The 0.1 Jy recommendation is validated only against foregrounds that exactly match the separation model; because the point-source marginalizer is spectrally similar to synchrotron, real SED complexity could change the bias balance before it is tested.","rationale":"The reader identified the matched-foreground assumption as the weakest structural premise, and I agree. The paper's own sentences b144, b228, b309, and b312 acknowledge the main idealizations, so the issue is not a hidden flaw; it is an untested boundary of the central claim. What makes it load-bearing rather than merely an acknowledged simplification is the spectral similarity between the marginalizing component and synchrotron: both are power laws in Rayleigh-Jeans temperature units, and the separation between them in the fiducial setup relies on the synthetic synchrotron having exactly the model's fixed SED. If real synchrotron SED curvature or spatial variation decorrelates the synchrotron signal across the six SO bands, the point-source component can absorb part of it, changing the reported biases in A_lens, alpha_s, beta_s, and A_s. A PySM s1/d1 test would settle this directly. I also note a secondary wording issue: Table 3 shows that adding the point-source component roughly doubles the A_lens uncertainty (0.005 to 0.010 extended; 0.017 to 0.042 baseline) and increases the alpha_s uncertainty, so the abstract's 'negligible cost in parameter uncertainties' overstates the cost for parameters other than r. This does not undermine the main scientific conclusion about r, but the recommendation should be stated as conditional on realistic foreground complexity.","tokens_in":33629,"tokens_out":10145,"duration_ms":78455,"concrete_test":"Generate synthetic cross-frequency spectra from PySM s1/d1 maps (spatially varying, curved SEDs) over the SO SAT footprint, inject the same ERS residual as in Section 4.2.3 at S_lim=0.1 Jy with f_FWHM=1, and run the paper's BBPower pipeline including the power-law point-source component. If the recovered A_lens, alpha_s, beta_s, or A_s shift by more than their 1-sigma uncertainties, the matched-foreground assumption is load-bearing and the 0.1 Jy recommendation should be restated as conditional on foreground complexity.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"Section 3.1 simulates the Galactic foregrounds with the same parametric form assumed by the component-separation model: pure power-law D_l spectra with fixed SEDs, Eqs. (30)-(31), and the likelihood uses exactly this form, Eq. (36). The central recommendation (abstract and Section 4.2.4) is therefore a matched-model validation. This is the right first test, but it is load-bearing because the additional component used to marginalize over residual ERSs is also a power-law SED with a white angular spectrum (A_p, beta_p; Section 4.2.4), the same functional family as synchrotron in frequency. On a real sky, synchrotron and dust SEDs vary along the line of sight and across the sky, so the effective synchrotron SED is not the single power law of Eq. (31). Under that mismatch, the power-law point-source component can partially absorb unmodelled synchrotron signal, and the reported near-zero biases in alpha_s, A_s, beta_s, and A_lens at S_lim=0.1 Jy may not survive. The authors explicitly acknowledge the optimistic foreground assumption (sentence b144), but the quantitative headline—the 0.1 Jy threshold and the claim that adding the point-source component restores unbiased foreground parameters—has not been tested against this failure mode. The DSFG and polarization-fraction robustness checks do not address it.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"This paper assesses the impact of polarised extragalactic point sources on CMB B-mode parameter inference for ground-based experiments such as the Simons Observatory Small Aperture Telescopes, and proposes masking and modelling strategies that keep the recovered parameters unbiased. The authors forecast source counts and residual polarised cross-frequency power spectra for extragalactic radio sources (ERSs) using the Toffolatti et al. (1998) and Tucci et al. (2011) number-count and SED models with a constant polarisation fraction of 2.8%, and for dusty star-forming galaxies (DSFGs) using the two-population model of Zagatti et al. (2024) with a 1.4% polarisation fraction. Masking uses a common mask across frequencies, defined at 27 GHz, with the mask radius parametrised by f_FWHM and residual flux fractions computed analytically. The contamination is propagated through a power-spectrum component-separation pipeline (BBPower) applied to noise-free synthetic data containing CMB (r=0, A_lens=0.3) and power-law Galactic dust and synchrotron with fixed SEDs (Eqs. 30-31); the covariance includes instrument noise, narrow-kernel-approximation mode coupling, and sky-fraction rescaling due to the point-source mask.","tokens_in":34134,"tokens_out":29640,"duration_ms":182234,"significance":"If the results hold, the paper provides directly actionable guidance for ongoing B-mode analyses: a concrete ERS flux threshold, a recommended mask radius at 27 GHz, and a two-parameter marginalising component that absorbs residual contamination without inflating sigma(r). The asymmetry the authors find—r robust while A_lens and the synchrotron tilt are strongly biased—is non-obvious and worth communicating, as is the physically motivated DSFG null result. The analysis has genuine strengths: the synthetic data are noise-free, so the reported offsets are genuine model bias; the input source-count and SED models are published and are propagated with careful treatment of stochastic spectral parameters; the covariance uses the public NaMaster NKA implementation and a realistic SO SAT footprint; the parameter studies span three noise configurations, a grid of mask radii, flux cuts from 0.01 to 1 Jy, and a polarisation-fraction scatter scenario; and the forecasting code is publicly available. The paper also states its main structural assumption explicitly (the matched foreground model of Eqs. 30-31), which makes the scope of the forecast checkable.","major_comments":[{"comment":"The reported sky-fraction loss is not consistent with the stated mask geometry. For the fiducial case (S_lim=0.1 Jy, f_FWHM=1, nu_m=27 GHz), Eq. (7) gives R_mask=45.5 arcmin and Omega_mask=pi R^2 about 5.5×10^-4 sr; combining this with N_PS=1699 sr^-1 in Eq. (34) yields f_sky=0.097 exp(-0.935)=3.8%, whereas Table 3 reports 7.8%. The reported value instead corresponds to R about 22.8 arcmin, half of the Eq. (7) radius. The same factor of two reproduces essentially every row with a quoted f_sky (for example, the high-resolution rows match R about 6.75 arcmin rather than the 13.5 arcmin from Eq. (7) with f_FWHM=3 and theta_FWHM=9 arcmin), and the S_lim=0.01 Jy rows are not mutually consistent because N_PS is listed as 25144 in the baseline configuration and 21544 in the extended configuration although the number counts depend only on the flux cut. Since Eq. (33) rescales the covariance by the f_sky ratio, all quoted uncertainties and sigma-level bias significances depend on this quantity; if the true remaining sky fraction at S_lim=0.1 Jy is about 3.8% rather than 7.8%, parameter errors would be roughly 40% larger and the reported biases correspondingly smaller in sigma units. Please state the mask radius actually used in the public code and correct either Eq. (7) or the f_sky computation.","section":"§4.2.2, Eqs. (7), (33)–(34), Table 3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The central recommendation is validated only against synthetic Galactic foregrounds that are exactly of the form assumed by the separation model (pure power-law D_l, fixed SEDs, Eqs. 30-31); the authors acknowledge this optimistic assumption in §3.1. The residual point-source marginaliser (A_p, beta_p; power-law SED in Rayleigh-Jeans units, white angular spectrum) belongs to the same frequency-scaling family as the synchrotron model, so under real foreground SED complexity—curvature, spatial variation, frequency decorrelation—the marginaliser could partially absorb unmodelled synchrotron signal, and the near-zero biases that support the 0.1 Jy recommendation might not persist. The robustness checks in §4.2.6 (polarisation-fraction scatter) and §4.3 (DSFGs) do not exercise this failure mode. In addition, §3.3 states that frequency decorrelation was tested “checking that such a term is not detected”, but no quantitative result for that test is presented anywhere in the paper. I recommend either adding a robustness run with a more complex synchrotron model (for example, a running spectral index or a PySM-style d1/s1 realisation with spatially varying SED) propagated through the same pipeline and marginaliser, or explicitly scoping the abstract's claim as conditional on the assumed foreground model; and reporting or removing the decorrelation test.","section":"§3.1, §3.3, §4.2.4, §4.2.6, §4.3","authors_acknowledge":"We note that the true Galactic sky is likely more complex than this, with each foreground component not exactly described by a power lawDℓ model and more complex spectral dependence, with foreground SEDs varying across the Galaxy, both across the sky and along any given line-of-sight (e.g. Krachmalnicoff et al. 2018b; PlanckCollaboration2020b;Martireetal.2022;Weilandetal.2022).","authors_acknowledge_p":0.77}],"minor_comments":[{"comment":"The abstract's claim of a “negligible cost in parameter uncertainties” is stronger than what Table 3 and §4.2.4 show for the recommended configuration: adding the marginalising component at S_lim=0.1 Jy roughly doubles the error on A_lens (extended configuration: ±0.005 to ±0.010) and increases the error on alpha_s by about 30%, while sigma(r) is indeed essentially unchanged. I suggest qualifying the sentence so that the negligible cost refers to r.","section":"Abstract, §4.2.4, Table 3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The manuscript contains several writing and typesetting artifacts: the sentence beginning “And find that, depending on the observing frequencies...” in §1 is a grammatical fragment; “we revised such forecasts” should read “reviewed”; and identifiers such as “Tuccietal.2011” appear with missing spaces in §2.3. These should be corrected in the published version.","section":"§1, §2.3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The adopted DSFG polarisation fraction <Pi_IR>=1.4% (De Zotti et al. 2018) sits below the stacked detections of Bonavera et al. (2017) that the paper cites (3.10±0.75% at 217 GHz and 3.65±0.66% at 353 GHz) and below the Trombetti et al. (2018) upper limits. Since the negligibility of DSFGs is one of the headline results, I suggest either justifying the choice against these values or adding a rescaling check (for example, a factor of five in DSFG power) to confirm that the null result is robust.","section":"§2.4.2, §4.3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The high-resolution row for S_lim=0.1 Jy with the marginalising component reads “0.0093±0.00071”, with the chi-squared entry concatenated into the same cell; the table should be reformatted and the precision of the uncertainties made uniform across rows.","section":"Table 3","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision"},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"Worth engaging. This paper does something prior forecasts (Puglisi, Lagache) did not: it pushes ERS and DSFG polarized emission through a multi-frequency power-spectrum pipeline and reports biases on r, A_lens, alpha_s, and the foreground parameters, plus concrete masking recipes. The design is clean: noiseless synthetic data ensure reported offsets are genuine model bias, and the covariance handles instrument noise, mode coupling, and sky-area loss. The code is public. The DSFG conclusion is credible within the model, and the nuance that r stays largely unbiased while A_lens and alpha_s are the vulnerable parameters is useful and non-obvious.\n\nThe main soft spot is exactly the stress-test concern: the simulated Galactic foregrounds have the same functional form as the component-separation model—power-law D_l spectra with fixed SEDs—and the residual point-source component used for marginalization is also a power-law SED with a white spectrum, the same frequency family as synchrotron. On a real sky with spatially varying SEDs, that component could absorb unmodelled synchrotron complexity, and the reported near-zero biases at S_lim = 0.1 Jy might not survive. The authors acknowledge the optimistic foreground assumption (b144), but the quantitative 0.1 Jy recommendation and the claim that the power-law component restores unbiased foreground parameters have not been tested against that failure mode. I don't think this is fatal—a matched-model forecast is the right first test—but the conclusions should present it as a matched-model result, and a robustness check against more realistic PySM-style foregrounds would materially strengthen the paper before publication.\n\nOther soft spots are minor and mostly disclosed: the fixed polarization fraction is tested at a higher value; E-to-B leakage from complex masks is explicitly deferred; and scaling the DSFG radio component by the same polarization fraction as the dust is acknowledged as an approximation. The pivot-frequency degeneracy in the marginalizer is handled sensibly with nu_p0 = 93 GHz. Citation practice looks standard: published source-count and SED models, the BBPower pipeline, and prior forecasts are all credited. No circularity.\n\nAudience: CMB B-mode forecasters and Simons Observatory analysis people. I'd send this to a serious referee; the referee should ask for a foreground-complexity robustness test before acceptance, but this paper earns the referee time.","headline":"Solid, honest forecast that fills a real gap—propagating polarized point sources through full component separation—but its headline 0.1 Jy recipe is only tested against foregrounds that match the separation model.","tokens_in":34608,"tokens_out":1613,"duration_ms":15662},"pith_extraction":{"msc":[],"pacs":["98.70.Vc","98.80.Es"],"model":"deepseek/deepseek-v4-flash","headline":"Unmasked polarised radio galaxies would bias the synchrotron and lensing B-mode parameters in a ground-based CMB experiment; masking above 0.1 Jy plus modelling the residue removes the bias at negligible cost.","keywords":["cosmic microwave background","primordial B-modes","tensor-to-scalar ratio","extragalactic radio sources","dusty star-forming galaxies","point source masking","component separation","foreground contamination"],"falsifier":"Measure the effective polarisation fraction $\\langle\\Pi^2\\rangle^{1/2}$ of radio sources near $0.1\\,$Jy in the 90–150 GHz bands: if it is closer to 4.94% than to 2.8%, the paper's own scaling (Section 4.2.6) implies the $0.1\\,$Jy residual biases on $A_{\\rm lens}$ and $\\alpha_s$ grow to roughly $3\\sigma$, overturning the recommendation; a complementary check is to re-run the pipeline on simulations with spatially varying synchrotron spectral indices and verify that the power-law residual component still separates cleanly.","tokens_in":35153,"feed_emoji":"📡","tokens_out":14939,"duration_ms":85563,"temperature":0.7,"pith_summary":"This paper asks whether polarised extragalactic point sources — radio galaxies and dusty star-forming galaxies — can corrupt the search for primordial $B$ modes in the cosmic microwave background, and what a real analysis must do about them. The authors inject state-of-the-art models of both source populations into synthetic cross-frequency $B$-mode spectra for a small-aperture ground-based experiment targeting noise levels near $\\sigma_P \\simeq 1\\,\\mu\\mathrm{K}\\,\\mathrm{arcmin}$, then run a full component-separation and parameter-estimation pipeline. They find that dusty galaxies are too faint to matter, while unmasked radio sources bias the recovered synchrotron amplitude and tilt and the lensing $B$-mode amplitude by up to several $\\sigma$, even though the tensor-to-scalar ratio $r$ itself stays unbiased. The proposed remedy is concrete: mask all radio sources with fluxes above roughly $0.1\\,$Jy at the lowest channel and model whatever remains as a Poisson component with a power-law spectrum. That recipe recovers unbiased cosmological and foreground parameters at a cost of roughly ten to forty percent in uncertainty.","feed_headline":"Mask radio sources above 0.1 Jy to keep B-modes unbiased","feed_subtitle":"Unmasked point sources bias synchrotron and lensing parameters by up to 12 sigma; the fix costs 10-40 percent in errors.","key_machinery":"The forecasting machinery is a pair of integrals (Eqs. 2 and 6) that turn source number counts, SED distributions, and a masking criterion into the number of masked sources and the residual cross-frequency polarised power spectra; the masking fraction $f_\\nu$ (Eq. 9) is the part that encodes both the flux cut and how much of each source's flux survives a finite masking radius. The analysis machinery is a power-spectrum-based Bayesian component separation that models all observed cross-spectra as a linear combination of CMB, dust, and synchrotron components, each with a fixed frequency SED and angular shape, and samples the parameters with MCMC. The load-bearing addition is the effective point-source component: a flat (Poisson-like) angular spectrum with a power-law SED, which is deliberately simple yet, with its pivot set at the most sensitive channel, absorbs the residual contamination well enough that all parameter biases drop below one $\\sigma$.","core_discovery":"The central claim is that residual polarised extragalactic radio sources (ERSs) are a genuine contaminant for $B$-mode analyses at these noise levels, not primarily because they threaten $r$, but because they masquerade as other components: their power-law spectra in Rayleigh-Jeans units resemble Galactic synchrotron, and their Poisson (white) angular spectrum resembles the lensing $B$-mode template on the scales studied ($32<\\ell\\le600$). Left unmasked at the flux limits reachable by small-aperture telescopes, they bias the synchrotron tilt $\\alpha_s$ by up to $12\\sigma$ and the lensing amplitude $A_{\\rm lens}$ by up to $7\\sigma$ in the most sensitive configuration considered. The paper shows the biases are avoidable: masking all ERSs above roughly $0.1\\,$Jy at 27 GHz, with a mask radius of one beam FWHM, reduces them to about one $\\sigma$ while inflating uncertainties by only 10–40 percent, and adding one extra component — a white angular spectrum $C_\\ell = A_p$ with a power-law SED $f_p(\\nu)=(\\nu/\\nu_{p0})^{\\beta_p}$, pivoted at the most sensitive channel — removes the residual bias entirely, even in the laxest masking case, at the price of roughly doubling the errors on $A_{\\rm lens}$ and $\\alpha_s$ while leaving $r$ essentially unchanged. By contrast, polarised dusty star-forming galaxies, modelled with a two-population halo model, produce no significant bias at these noise levels and require no masking.","pith_inferences":["Editorial extension: the margin of safety on the $0.1\\,$Jy recommendation is set by the assumed polarisation fraction — the paper's own test shows that raising the effective fraction from 2.8% to 4.94% triples the residual power and pushes the $0.1\\,$Jy biases to about $3\\sigma$; direct measurements of faint-source polarisation near 90–150 GHz would settle how much margin remains.","Editorial extension: because the residual component and synchrotron are both power laws in Rayleigh-Jeans units, their separation is a spectral-index game; with real, spatially varying synchrotron SEDs the fitted $\\beta_p$ could absorb genuine Galactic complexity, so the unbiased-$\\alpha_s$ result is best read as proven for the matched foreground model.","Editorial extension: the analysis implies a division of labour for future $B$-mode pipelines — protect $r$ with moderate masking, protect $A_{\\rm lens}$ and synchrotron science with the residual-component fit — which suggests that delensing performance, not the tensor-to-scalar ratio, is the most contamination-sensitive step in the search."],"forward_implications":["A ground-based $B$-mode experiment at roughly $1\\,\\mu$K-arcmin noise can protect its results with one achievable choice: mask radio sources above $\\sim0.1\\,$Jy at the lowest channel and fit a single extra Poisson-like component.","Unmasked radio sources would be absorbed by the synchrotron and lensing parameters, biasing $\\alpha_s$ by up to $12\\sigma$ and $A_{\\rm lens}$ by up to $7\\sigma$ in the most sensitive configuration, while $r$ remains unbiased; the vulnerable science is therefore delensing and synchrotron foreground characterisation, not $r$ itself.","The $0.1\\,$Jy mask costs only 10–40 percent in parameter uncertainties, and the residual-component fit roughly doubles the errors on $A_{\\rm lens}$ and $\\alpha_s$ while barely touching $r$.","Dusty star-forming galaxies need no masking at the noise levels considered: even completely unmasked they cause no significant bias in any model parameter.","High-resolution instruments are a different regime: they require a stricter $0.01\\,$Jy flux cut together with the same residual-component fit, because small beams make the residual point-source signal far more damaging near 90–150 GHz."],"supporting_citations":[{"why":"Supplies the flat+inverted ERS number counts at 5 GHz, the spectral-index distributions, and the high-frequency SED models used for steep, flat, and inverted populations.","marker":"Tucci et al. (2011)"},{"why":"Supplies the total radio-source number counts at 5 GHz from which the steep population is derived by difference.","marker":"Toffolatti et al. (1998)"},{"why":"Supplies the 2.8% polarisation fraction, the per-band flux limits used as the initial masking benchmarks, and the method for converting intensity counts to polarised point-source spectra.","marker":"Lagache et al. (2020)"},{"why":"Supplies the scatter-enhanced polarisation fraction (4.94%) used in the robustness test of Section 4.2.6.","marker":"Puglisi et al. (2018)"},{"why":"Supplies the two-population DSFG model — luminosity functions, SED templates, and halo occupation distributions — from which DSFG counts and polarised spectra are computed.","marker":"Zagatti et al. (2024)"},{"why":"Supplies the 1.4% mean polarisation degree adopted for dusty star-forming galaxies.","marker":"De Zotti et al. (2018)"},{"why":"Supplies the power-spectrum-based component separation and likelihood methodology on which the forecasting pipeline is built.","marker":"Wolz et al. (2024)"},{"why":"Supplies the frequency bands, noise model, beams, and sky footprint used to define the synthetic experiment.","marker":"The Simons Observatory Collaboration (2019)"}],"fun_headline_variants":["Mask radio sources above 0.1 Jy to keep B-mode parameters unbiased","Point-source bias hits synchrotron and lensing; a 0.1 Jy cut fixes it","Dusty galaxies are harmless; radio sources drive the B-mode bias risk","Radio-source masking above 0.1 Jy removes B-mode contamination bias","Unmasked radio sources bias synchrotron by 12 sigma; 0.1 Jy mask helps"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The load-bearing premise is that the simulated Galactic dust and synchrotron have exactly the functional form the fitting model assumes — power-law angular spectra with fixed, spatially constant SEDs — so the unbiased recovery at $S_{\\rm lim}=0.1\\,$Jy is demonstrated against a perfectly matched foreground, and the authors themselves flag this as an optimistic scenario.","fun_headline_variants_meta":{"raw":{"variants":["Mask radio sources above 0.1 Jy to keep B-mode parameters unbiased","Point-source bias hits synchrotron and lensing; a 0.1 Jy cut fixes it","Dusty galaxies are harmless; radio sources drive the B-mode bias risk","Radio-source masking above 0.1 Jy removes B-mode contamination bias","Unmasked radio sources bias synchrotron by 12 sigma; 0.1 Jy mask helps"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":0.000409,"raw_usage":{"cost":0.0004085718,"is_byok":false,"cost_details":{"upstream_inference_cost":0.0004085718,"upstream_inference_prompt_cost":3.12858e-05,"upstream_inference_completions_cost":0.000377286},"total_tokens":3300,"prompt_tokens":1227,"completion_tokens":2073,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":1104,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":1973}},"tokens_in":1227,"tokens_out":2073,"cost_source":"openrouter_usage","duration_ms":14704,"image_pages":[],"temperature":1.0,"billed_attempts":1,"reasoning_tokens":1973,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-23T10:43:06.219561+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"pipeline_version":"pith-pipeline@v1.4.1-daily-deepseek-or-medium","input_kind":null,"input_chars":null,"input_total_chars":null,"passes_missing":null,"published_standing":null,"published_confidence":null,"objection_audit":null,"run_record":null,"basis":{"kind":"not_recorded","source":null,"chars_given":null,"chars_total":null},"receipt":null,"falsifier":"Measure the effective polarisation fraction $\\langle\\Pi^2\\rangle^{1/2}$ of radio sources near $0.1\\,$Jy in the 90–150 GHz bands: if it is closer to 4.94% than to 2.8%, the paper's own scaling (Section 4.2.6) implies the $0.1\\,$Jy residual biases on $A_{\\rm lens}$ and $\\alpha_s$ grow to roughly $3\\sigma$, overturning the recommendation; a complementary check is to re-run the pipeline on simulations with spatially varying synchrotron spectral indices and verify that the power-law residual component still separates cleanly.","supporting_citations":[{"cited_title":"The Simons Observatory: pipeline comparison and validation for large-scale B-modes","cited_arxiv_id":"2302.04276","evidence_quote":"Supplies the power-spectrum-based component separation and likelihood methodology on which the forecasting pipeline is built."}],"lane_accuracy":{"status":"not yet measured","lane":{"model":"deepseek/deepseek-v4-flash","input_kind":null},"url":"https://pith.science/accuracy"},"review_version":1}