SEARCH · Engineering Papers

Results for “likelihood ratio”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2024

High-dimensional and permutation invariant anomaly detection

Methods for anomaly detection of new physics processes are often limited to low-dimensional spaces due to the difficulty of learning high-dimensional probability densities. Particularly at the constituent level, incorporating desirable properties such as permutation invariance and variable-length inputs becomes difficult within popular density estimation methods. In this work, we introduce a permutation-invariant density estimator for particle physics data based on diffusion models, specifically designed to handle variable-length inputs. We demonstrate the efficacy of our methodology by utilizing the learned density as a permutation-invariant anomaly detection score, effectively identifying jets with low likelihood under the background-only hypothesis. To validate our density estimation method, we investigate the ratio of learned densities and compare to those obtained by a supervised classification algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Galactic Chemical Evolution of Short-lived Radioactive Isotopes Produced by Explosive Nucleosynthesis: 60 Fe and 53 Mn

Several short-lived radionuclides (SLRs) are know to have existed in the early solar system (ESS). These species, which typically decay with half-lives of the order of a few million years, can be used to probe the timescales of events preceding the birth of the Sun. We investigate the ESS origin of 53 Mn, produced by core-collapse supernovae (CCSNe) and Type Ia supernovae (SNe Ia), and 60 Fe, produced exclusively by CCSNe. We model the evolution of the radioactive-to-stable abundance ratios of these SLRs with a Galactic chemical evolution (GCE) framework accounting for different SN yields, SN Ia delay times, and other Galactic features (K). A further set of models is calculated assuming that SNe Ia did not contribute any 53 Mn to the ESS. The predicted ratios are compared to meteoritic ratios to derive a distribution of solar isolation times that includes uncertainties due to stochastic chemical enrichment and measurements of the ESS values. The isolation times are then compared to those of 107 Pd and 182 Hf calculated in previous work. A self-consistent solution can be found within the current uncertainties, especially when using the GCE setups with K = 1.6 and 2.3, although the maximum likelihood for the 60 Fe distribution is typically ∼4–5 Myr shorter than for 53 Mn. The predicted 60 Fe/ 53 Mn ratio, instead, is completely inconsistent with the ESS value; this could be resolved using a larger fraction of faint CCSNe than usually considered in GCE models.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Almost medium-free measurement of the Hoyle state direct-decay component with a TPC

The structure of the Hoyle state, a highly α -clustered state at 7.65 MeV in 12 C, has long been the subject of debate. Understanding if the system comprises of three weakly interacting α particles in the 0s orbital, known as an α-condensate state, is possible by studying the decay branches of the Hoyle state. The direct decay of the Hoyle state into three α particles, rather than through the 8 Be ground state, can be identified by studying the energy partition of the three α particles arising from the decay. This paper provides details on the breakup mechanism of the Hoyle stating using a new experimental technique. Method: By using β-delayed charged-particle spectroscopy of 12 N using the Texas active target time-projection chamber, a high-sensitivity measurement of the direct 3α decay ratio can be performed without contributions from pileup events. A Bayesian approach to understanding the contribution of the direct components via a likelihood function shows that the direct component is <0.043% at the 95% confidence level. This value is in agreement with several other studies, and, here, we can demonstrate that a small nonsequential component with a decay fraction of about 10 –4 is most likely. Here, the measurement of the non-sequential component of the Hoyle state decay is performed in an almost medium-free reaction for the first time. The derived upper limit is in agreement with previous studies and demonstrates sensitivity to the absolute branching ratio. Further experimental studies would need to be combined with robust microscopic theoretical understanding of the decay dynamics to provide additional insight into the idea of the Hoyle state as an α condensate.

6 ≤ A ≤ 19↗

DOE OSTI2020

Prospects for Characterizing the Haziest Sub-Neptune Exoplanets with High-resolution Spectroscopy

Observations to characterize planets larger than Earth but smaller than Neptune have led to largely inconclusive interpretations at low spectral resolution due to hazes or clouds that obscure molecular features in their spectra. However, here we show that high-resolution spectroscopy (R ∼ 25,000–100,000) enables one to probe the regions in these atmospheres above the clouds where the cores of the strongest spectral lines are formed. We present models of transmission spectra for a suite of GJ 1214b–like planets with thick photochemical hazes covering 1–5 μm at a range of resolutions relevant to current and future ground-based spectrographs. Furthermore, we compare the utility of the cross-correlation function that is typically used with a more formal likelihood-based approach, finding that only the likelihood-based method is sensitive to the presence of haze opacity. We calculate the signal-to-noise ratio (S/N) of these spectra, including telluric contamination, Required to robustly detect a host of molecules such as CO, CO{sub 2}, H{sub 2}O, and CH{sub 4} and photochemical products like HCN as a function of wavelength range and spectral resolution. Spectra in the M band require the lowest S/N{sub res} to detect multiple molecules simultaneously. CH{sub 4} is only observable for the coolest models (T {sub eff} = 412 K) and only in the L band. We quantitatively assess how these requirements compare to what is achievable with current and future instruments, demonstrating that characterization of small cool worlds with ground-based high-resolution spectroscopy is well within reach.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Constraints on directionality effect of nuclear recoils in a liquid argon time projection chamber

The direct search for dark matter in the form of weakly interacting massive particles (WIMP) is performed by detecting nuclear recoils produced in a target material from the WIMP elastic scattering. The experimental identification of the direction of the WIMP-induced nuclear recoils is a crucial asset in this field, as it enables unmistakable modulation signatures for dark matter. The Recoil Directionality (ReD) experiment was designed to probe for such directional sensitivity in argon dual-phase time projection chambers (TPC), that are widely considered for current and future direct dark matter searches. The TPC of ReD was irradiated with neutrons at the INFN Laboratori Nazionali del Sud. Data were taken with nuclear recoils of known directions and kinetic energy of 72 keV, which is within the range of interest for WIMP-induced signals in argon. The direction-dependent liquid argon charge recombination model by Cataudella et al. was adopted and a likelihood statistical analysis was performed, which gave no indications of significant dependence of the detector response to the recoil direction. The aspect ratio R of the initial ionization cloud is $R < 1.072$ with 90 % confidence level.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

A Constraint on Primordial B-modes from the First Flight of the Spider Balloon-borne Telescope

We present the first linear polarization measurements from the 2015 long-duration balloon flight of Spider, which is an experiment that is designed to map the polarization of the cosmic microwave background (CMB) on degree angular scales. The results from these measurements include maps and angular power spectra from observations of 4.8% of the sky at 95 and 150 GHz, along with the results of internal consistency tests on these data. While the polarized CMB anisotropy from primordial density perturbations is the dominant signal in this region of sky, Galactic dust emission is also detected with high significance. Galactic synchrotron emission is found to be negligible in the Spider bands. We employ two independent foreground-removal techniques to explore the sensitivity of the cosmological result to the assumptions made by each. The primary method uses a dust template derived from Planck data to subtract the Galactic dust signal. A second approach, which constitutes a joint analysis of Spider and Planck data in the harmonic domain, assumes a modified-blackbody model for the spectral energy distribution of the dust with no constraint on its spatial morphology. Using a likelihood that jointly samples the template amplitude and r parameter space, we derive 95% upper limits on the primordial tensor-to-scalar ratio from Feldman–Cousins and Bayesian constructions, finding r < 0.11 and r < 0.19, respectively. Roughly half the uncertainty in r derives from noise associated with the template subtraction. New data at 280 GHz from Spider's second flight will complement the Planck polarization maps, providing powerful measurements of the polarized Galactic dust emission.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

Photometry of Outer Solar System Objects from the Dark Energy Survey. I. Photometric Methods, Light-curve Distributions, and Trans-Neptunian Binaries

We report the methods of and initial scientific inferences from the extraction of precision photometric information for the >800 trans-Neptunian objects (TNOs) discovered in the images of the Dark Energy Survey (DES). Scene-modeling photometry is used to obtain shot-noise-limited flux measures for each exposure of each TNO, with background sources subtracted. Comparison of double-source fits to the pixel data with single-source fits are used to identify and characterize two binary TNO systems. A Markov Chain Monte Carlo method samples the joint likelihood of the intrinsic colors of each source as well as the amplitude of its flux variation, given the time series of multiband flux measurements and their uncertainties. A catalog of these colors and light-curve amplitudes A is included with this publication. We show how to assign a likelihood to the distribution q(A) of light-curve amplitudes in any subpopulation. Using this method, we find decisive evidence (i.e., evidence ratio <0.01) that cold classical (CC) TNOs with absolute magnitude 6 < H$_{r}$ < 8.2 are more variable than the hot classical (HC) population of the same H$_{r}$, reinforcing theories that the former form in situ and the latter arise from a different physical population. Resonant and scattering TNOs in this H$_{r}$ range have variability consistent with either the HCs or CCs. DES TNOs with H$_{r}$ < 6 are seen to be decisively less variable than higher-H$_{r}$ members of any dynamical group, as expected. More surprising is that detached TNOs are decisively less variable than scattering TNOs, which requires them to have distinct source regions or some subsequent differential processing.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2025

Systematic Bayesian evaluation of resonance parameters in 19 Ne for the 15 O ⁡(𝛼, 𝛾)⁢ 19 Ne and 18 F ⁡(𝑝, 𝛼)⁢ 15 O reactions

Here, we present a comprehensive evaluation of the nuclear structure properties of 19 Ne using a novel and rigorous Bayesian statistical framework. Precise characterization of 19 Ne resonance parameters is critical for accurately determining reaction rates of the astrophysically significant 15 O ⁡(𝛼, 𝛾)⁢ 19 Ne and 18 F ⁡(𝑝, 𝛼)⁢ 15 O reactions, which govern breakout from the hot CNO cycle in x-ray bursts and influence 𝛾-ray emission in novae, respectively. By reconstructing likelihood functions from published experimental data—including asymmetric uncertainties and upper or lower limits—we derive posterior distributions for resonance energies, decay widths, and branching ratios. Our Bayesian approach systematically incorporates previously reported discrepancies among measurements, providing a statistically robust and consistent treatment of these uncertainties. The evaluated resonance parameters and associated uncertainties provide crucial input for stellar nucleosynthesis modeling, contributing to a refined understanding of explosive astrophysical phenomena.

Kim, Sohyun H. [Sungkyunkwan Univ., Suwon (Republi↗

DOE OSTI2021

Detection and Bulk Properties of the HR 8799 Planets with High-resolution Spectroscopy

Using the Keck Planet Imager and Characterizer, we obtained high-resolution (R ~ 35,000) K-band spectra of the four planets orbiting HR 8799. We clearly detected H{sub 2}O and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and H{sub 2}O in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral resolution to date when considering both their angular separations and flux ratios. We developed a forward-modeling framework that allows us to jointly fit the spectra of the planets and the diffracted starlight simultaneously in a likelihood-based approach and obtained posterior probabilities on their effective temperatures, surface gravities, radial velocities, and spins. We measured vsin (i) values of 10.1{sub -2.7}{sup +2.8} km s{sup -1} for HR 8799 d and 15.0{sub -2.6}{sup +2.3} km s{sup -1} for HR 8799 e, and placed an upper limit of <14 km s{sup -1} of HR 8799 c. Under two different assumptions of their obliquities, we found tentative evidence that rotation velocity is anticorrelated with companion mass, which could indicate that magnetic braking with a circumplanetary disk at early times is less efficient at spinning down lower-mass planets.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2025

Multimodality in the Search for New Physics in Pulsar Timing Data and the Case of Kination-amplified Gravitational-wave Background from Inflation

We investigate the kination-amplified inflationary gravitational-wave background (GWB) interpretation of the signal recently reported by various pulsar timing array (PTA) experiments. Kination is a post-inflationary phase in the expansion history dominated by the kinetic energy of some scalar field, characterized by a stiff equation of state w = 1. Within the inflationary GWB model, we identify two modes that can fit the current data sets (NANOGrav and EPTA) with equal likelihood: the kination-amplification (KA) mode and the ordinary, no-kination-amplification (no-KA) mode. The multimodality of the likelihood motivates a Bayesian analysis with nested sampling. We analyze the free spectra of current PTA data and mock free spectra constructed with higher signal-to-noise ratios using nested sampling. The analysis of the mock spectrum designed to be consistent with the best fit to the NANOGrav 15 yr (NG15) data successfully reveals the expected bimodal posterior for the first time while excluding the reheating mode that appears in the fit to the current NG15 data, making a case for our correct and comprehensive treatment of potential multimodal posteriors arising from future PTA data sets. The resultant Bayes factor is $\mathcal{B}$ $\equiv$ Z no–KA /Z KA = 2.9 ± 1.9, indicating comparable statistical significance between the two modes. Given the theoretical model-building challenges of producing highly blue-tilted primordial tensor spectra, the KA mode has the advantage of requiring less blue primordial spectra, compared with the no-KA mode. The synergy between future cosmic microwave background polarization, pulsar timing, and laser interferometer measurements of gravitational waves will help resolve the ambiguity implied by the multimodal posterior in PTA-only searches.

Cosmology↗

DOE OSTI2021

Measurement of the cross section of top quark pairs produced in association with a Photon in Lepton +Jets events at $\sqrt{s}$ = 13 TeV with full Run II CMS data

The inclusive production cross section of top quark pairs in association with a photon is measured in proton-proton collisions at the LHC with 13 TeV energy using the full RunII data collected by CMS in 2016, 2017, and 2018 with a total corresponding integrated luminosity of 137 fb$^{-1}$. The relative fraction of t$ \bar{t} \gamma$ events normalized to inclusive t$ \bar{t}$ production is measured. The cross section measurement provides important information about the electromagnetic coupling of the standard model top quark and is sensitive to physics beyond the standard model. The analysis is carried out in the in semileptonic decay channel with a well isolated high P$ \_{T}$ lepton (electron and muon), at least four jets out of which at least one must be b-tagged, and an isolated photon. Photons may be emitted from initial state radiation, top quarks, and decay products of top quarks. The fiducial phase space with a photon of transverse momentum $ \ge$ 20 GeV and j\x11$ \eta$ <1.4441 is used. A simultaneous likelihood fit of control regions with the signal region is done to constraint the backgrounds and to extract the t$ \bar{t} \gamma$ cross section. The measurement of the ratio of t$ \bar{t} \gamma$ to t$ \bar{t}$ is $0.02055 \pm 0.00099$ (syst.) $ \pm 0.00099$ (stat.) in $e$ channel , 0.02156 $ \pm$ 0.00068 (syst.) $ \pm$ 0.00068 (stat.) in $ \mu$ channel , and 0.02203 $ \pm$ 0.00064 (syst.) $ \pm$ 0.00064 (stat.) in e + $ \mu$ channel , respectively. The measured inclusive cross section is 3.81 $ \pm$ 0.15 (syst.) $ \pm$ 0.10 (stat.) pb in e channel , 3.87 $ \pm$ 0:11 (syst.) $ \pm$ 0:07 (stat.) pb in $ \mu$ channel , and 3.96 $ \pm$ 0.10 (syst.) $ \pm$ 0.06 (stat.) pb in e + $ \mu$ channel , respectively for full RunII data. The total inclusive cross section with a photon of transverse momentum $ \ge$ 20 GeV is found to be 3.96 $ \pm$ 0.10 (syst.) $ \pm$ 0.06 (stat.) pb with full RunII CMS data. The results are in agreement with the standard model next to leading order prediction.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Test of light-lepton universality in τ decays with the Belle II experiment

We present a measurement of the ratio R µ = $\mathcal{B}$($τ^- → µ^-\overline{ν}_µν_τ$)/$\mathcal{B}(τ^- → e^-\overline{ν}_eν_τ$) of branching fractions B of the τ lepton decaying to muons or electrons using data collected with the Belle II detector at the SuperKEKB e + e - collider. The sample has an integrated luminosity of 362 ± 2 fb -1 at a centre-of-mass energy of 10.58 GeV. Using an optimised event selection, a binned maximum likelihood ft is performed using the momentum spectra of the electron and muon candidates. The result, R µ = 0.9675 ± 0.0007 ± 0.0036, where the first uncertainty is statistical and the second is systematic, is the most precise to date. It provides a stringent test of the light-lepton universality, translating to a ratio of the couplings of the muon and electron to the W boson in τ decays of 0.9974 ± 0.0019, in agreement with the standard model expectation of unity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Dark Energy Survey Year 3 results: redshift calibration of the weak lensing source galaxies

ABSTRACT Determining the distribution of redshifts of galaxies observed by wide-field photometric experiments like the Dark Energy Survey (DES) is an essential component to mapping the matter density field with gravitational lensing. In this work we describe the methods used to assign individual weak lensing source galaxies from the DES Year 3 Weak Lensing Source Catalogue to four tomographic bins and to estimate the redshift distributions in these bins. As the first application of these methods to data, we validate that the assumptions made apply to the DES Y3 weak lensing source galaxies and develop a full treatment of systematic uncertainties. Our method consists of combining information from three independent likelihood functions: self-organizing map p(z) (sompz), a method for constraining redshifts from galaxy photometry; clustering redshifts (WZ), constraints on redshifts from cross-correlations of galaxy density functions; and shear ratios (SRs), which provide constraints on redshifts from the ratios of the galaxy-shear correlation functions at small scales. Finally, we describe how these independent probes are combined to yield an ensemble of redshift distributions encapsulating our full uncertainty. We calibrate redshifts with combined effective uncertainties of σ〈z〉 ∼ 0.01 on the mean redshift in each tomographic bin.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

A Full Accounting of the Visible Mass in SDSS MaNGA Disk Galaxies

We present a study of the ratio of visible mass to total mass in spiral galaxies to better understand the relative amount of dark matter present in galaxies of different masses and evolutionary stages. Using the velocities of the Hα emission line measured in spectroscopic observations from the Sloan Digital Sky Survey (SDSS) MaNGA Data Release 17 (DR 17), we evaluate the rotational velocity of over 5500 disk galaxies at their 90% elliptical Petrosian radii, R 90 . We compare this to the velocity expected from the total visible mass, which we compute from the stellar, H i, H 2 , and heavy metals and dust masses. H 2 mass measurements are available for only a small subset of galaxies observed in SDSS MaNGA DR17, so we derive a parameterization of the H 2 mass as a function of absolute magnitude in the r band using galaxies observed as part of SDSS DR7. With these parameterizations, we calculate the fraction of visible mass within R 90 that corresponds to the observed velocity. Based on statistically analyzing the likelihood of this fraction, we conclude that the null hypothesis (no dark matter) cannot be excluded at a confidence level better than 95% within the visible extent of the disk galaxies. We also find that when all mass components are included, the ratio of visible to total mass within the visible extent of star-forming disk galaxies increases with galaxy luminosity.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Improved Constraints on Primordial Gravitational Waves using Planck , WMAP, and BICEP/ Keck Observations through the 2018 Observing Season

Here, we present results from an analysis of all data taken by the BICEP2, Keck Array, and BICEP3 CMB polarization experiments up to and including the 2018 observing season. We add additional Keck Array observations at 220 GHz and BICEP3 observations at 95 GHz to the previous 95/150/220 GHz dataset. The Q/U maps now reach depths of 2.8, 2.8, and 8.8 μK CMB arcmin at 95,150, and 220 GHz, respectively, over an effective area of ≈ 600 square degrees at 95 GHz and ≈ 400 square degrees at 150 and 220 GHz. The 220 GHz maps now achieve a signal-to-noise ratio on polarized dust emission exceeding that of Planck at 353 GHz. We take auto- and cross-spectra between these maps and publicly available WMAP and Planck maps at frequencies from 23 to 353 GHz and evaluate the joint likelihood of the spectra versus a multicomponent model of lensed ΛCDM + r + dust + synchrotron + noise . The foreground model has seven parameters, and no longer requires a prior on the frequency spectral index of the dust emission taken from measurements on other regions of the sky. This model is an adequate description of the data at the current noise levels. The likelihood analysis yields the constraint r 0.05 < 0.036 at 95% confidence. Running maximum likelihood search on simulations we obtain unbiased results and find that σ(r) = 0.009 . These are the strongest constraints to date on primordial gravitational waves.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Measurement of the $c$-jet mistagging efficiency in $t\bar{t}$ events using $pp$ collision data at $\sqrt{s}=13$ $\text {TeV}$ collected with the ATLAS detector

A technique is presented to measure the efficiency with which c-jets are mistagged as b-jets (mistagging efficiency) using tt¯ events, where one of the W bosons decays into an electron or muon and a neutrino and the other decays into a quark–antiquark pair. The measurement utilises the relatively large and known W → cs branching ratio, which allows a measurement to be made in an inclusive c-jet sample. The data sample used was collected by the ATLAS detector at √s = 13 TeV and corresponds to an integrated luminosity of 139 fb –1 . Events are reconstructed using a kinematic likelihood technique which selects the mapping between jets and tt¯ decay products that yields the highest likelihood value. The distribution of the b-tagging discriminant for jets from the hadronic W decays in data is compared with that in simulation to extract the mistagging efficiency as a function of jet transverse momentum. The total uncertainties are in the range 3–17%. The measurements generally agree with those in simulation but there are some differences in the region corresponding to the most stringent b-jet tagging requirement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Estimate of the Mass and Radial Profile of the Orphan–Chenab Stream's Dwarf-galaxy Progenitor Using MilkyWay@home

We fit the mass and radial profile of the Orphan–Chenab Stream's (OCS) dwarf-galaxy progenitor by using turnoff stars in the Sloan Digital Sky Survey and the Dark Energy Camera to constrain N-body simulations of the OCS progenitor falling into the Milky Way on the 1.5 PetaFLOPS MilkyWay@home distributed supercomputer. We infer the internal structure of the OCS's progenitor under the assumption that it was a spherically symmetric dwarf galaxy composed of a stellar system embedded in an extended dark matter halo. We optimize the evolution time, the baryonic and dark matter scale radii, and the baryonic and dark matter masses of the progenitor using a differential evolution algorithm. The likelihood score for each set of parameters is determined by comparing the simulated tidal stream to the angular distribution of OCS stars observed in the sky. We fit the total mass of the OCS's progenitor to (2.0 ± 0.3) × 10 7 M ⊙ with a mass-to-light ratio of γ = 73.5 ± 10.6 and (1.1 ± 0.2) × 10 6 M ⊙ within 300 pc of its center. Within the progenitor's half-light radius, we estimate a total mass of (4.0 ± 1.0) × 10 5 M ⊙ . We also fit the current sky position of the progenitor's remnant to be (α, δ) = ((166.0 ± 0.9)°, (–11.1 ± 2.5)°) and show that it is gravitationally unbound at the present time. The measured progenitor mass is on the low end of previous measurements and, if confirmed, lowers the mass range of ultrafaint dwarf galaxies. Our optimization assumes a fixed Milky Way potential, OCS orbit, and radial profile for the progenitor, ignoring the impact of the Large Magellanic Cloud.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

CMB-S4: Foreground-cleaning Pipeline Comparison for Measuring Primordial Gravitational Waves

We compare multiple foreground-cleaning pipelines for estimating the tensor-to-scalar ratio, r, using simulated maps of the planned CMB-S4 experiment within the context of the South Pole Deep Patch. To evaluate robustness, we analyze bias and uncertainty on r across various foreground suites using map-based simulations. The foreground-cleaning methods include: a parametric maximum likelihood approach applied to auto- and cross-power spectra between frequency maps; a map-based parametric maximum-likelihood method; and a harmonic-space internal linear combination using frequency maps. We summarize the conceptual basis of each method to highlight their similarities and differences. To better probe the impact of foreground residuals, we implement an iterative internal delensing step, leveraging a map-based pipeline to generate a lensing B-mode template from the large aperture telescope frequency maps. Our results show that the performance of the three approaches is comparable for simple and intermediate-complexity foregrounds, with σ(r) ranging from 3–5 ×10 −4 . However, biases at the 1σ–2σ level appear when analyzing more complex forms of foreground emission. By extending the baseline pipelines to marginalize over foreground residuals, we demonstrate that contamination can be reduced to within statistical uncertainties, albeit with a pipeline-dependent impact on σ(r), which translates to a detection significance between 2σ and 4σ for an input value of r = 0.003. These findings suggest varying levels of maturity among the tested pipelines, with the auto- and cross-spectra-based approach demonstrating the best stability and overall performance. Moreover, given the extremely low noise levels, mutual validation of independent foreground-cleaning pipelines is essential to ensure the robustness of any potential detection.

astronomy data analysis↗