Search NASA⌕ Search

SEARCH · Search NASA

Results for “Subtraction Technique”

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.

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

At least 253 records · Page 14

Exospheric temperatures deduced from 7320- to 7330-A /O/+//2P/ - O/+//2D// twilight observations

A technique developed to deduce exospheric temperatures from the 7320- to 7330-A emission measured by the visible airglow experiment on board the AE-E satellite is considered. An excess emission in the measured 7320- to 7330-A brightness is noticed as a result of the interaction between the spacecraft and the atmosphere. The observed brightnesses are corrected for this effect. The galactic background emission is also carefully subtracted. The deduced temperatures exhibit a positive correlation with solar activity. It varies from approximately 700 K in late 1976 to approximately 1700 K at the peak of this solar cycle. The presence of a nonthermal oxygen corona is considered inconclusive.

Yee, J. H.↗

Decomposition Technique for Remaining Useful Life Prediction

The prognostic tool disclosed here decomposes the problem of estimating the remaining useful life (RUL) of a component or sub-system into two separate regression problems: the feature-to-damage mapping and the operational conditions-to-damage-rate mapping. These maps are initially generated in off-line mode. One or more regression algorithms are used to generate each of these maps from measurements (and features derived from these), operational conditions, and ground truth information. This decomposition technique allows for the explicit quantification and management of different sources of uncertainty present in the process. Next, the maps are used in an on-line mode where run-time data (sensor measurements and operational conditions) are used in conjunction with the maps generated in off-line mode to estimate both current damage state as well as future damage accumulation. Remaining life is computed by subtracting the instance when the extrapolated damage reaches the failure threshold from the instance when the prediction is made.

Saha, Bhaskar↗

Preshaping clear glass at low temperatures

Advances in available glass chemistries and glass processing methods have accompanied and enabled some of the biggest technology revolutions, from the development and mass production of light bulbs to low-loss fiber optics and durable smartphone touchscreens. An emerging generation of low-temperature processing technologies aims to continue this important trend and make a broader array of glass components mass producible. In the issue, Mader et al. (1) describe one such innovation in glass processing—the use of low-temperature injection molding to preshape silica particle–filled composites that can later be transformed into transparent fused silica glass objects. Traditionally, transparent glass objects are manufactured in high volume from molten or softened glass, which is floated, drawn, blown, cast, or blow-molded to a desired shape (see the figure, top). The glass composition and processing technique dictate the working temperature, which is usually quite high (near 1000°C) and often restricts the choice of compatible equipment or limits the choice of glass composition. Because geometry-specific capital investment is required for production, drastic or frequent component design changes or small batches may be cost prohibitive. Alternatively, transparent glass components can also be shaped at ambient temperature from solid glass by a series of subtractive processes, including cutting or multiple stages of grinding, followed by slower processing steps, such as polishing or etching. This approach is somewhat less amenable to mass production, and certain geometries containing tool-inaccessible regions cannot be fabricated in this way. Several emerging glass-shaping technologies aim to reduce the required manufacturing temperatures and still provide access to a broader range of glass compositions and component geometries (see the figure). These approaches use a three-step process. First, a desired shape is preformed at low temperature from a glass-forming, organic-inorganic composite. Next, the preform is dried, and organic materials used to bind particles are removed. Finally, the preform is heated (sintered) well below the glass-melting temperature to densify to transparent glass. Although the second and third steps do occur at increased temperatures, only standard, geometry-agnostic driers and furnaces are required. This strategy builds on the well-studied sol-gel approach to forming monolithic glass, where silica network–forming chemical solutions are poured into molds, slowly dried, and condensed into dense glass without melting (2). In a departure from the sol-gel process, these new technologies use solvents, cross-linkers, and polymers to formulate organic-inorganic composites tuned for compatibility with a particular shaping process, with formats ranging from photocurable liquids to shear-thickening pastes to solids. The composite inorganic loadings are also typically higher than those in the pure sol-gel approach, which drastically reduces shrinkage in comparison.

36 MATERIALS SCIENCE↗

21 cm Power Spectrum Analysis of North Celestial Pole Observations with the Tianlai Dish Pathfinder Array

The Tianlai Dish Pathfinder Array (TDPA) is a radio interferometer designed to test techniques for 21 cm intensity mapping in the post-reionization Universe as a means of measuring large-scale cosmic structure. Using nine nights of observations targeting the North Celestial Pole field, totaling approximately 107 hr of integration time, we analyze data in the frequency range 700–800 MHz (corresponding to redshift z ∼ 0.9). We do the data format conversion, radio frequency interference flagging, calibration, imaging and point source subtraction, and foreground removal via Singular Value Decomposition. The spherically averaged power spectrum Δ 2 (k) is obtained. Furthermore, this work successfully establishes and validates a comprehensive data analysis framework for the TDPA. We identify key improvements including sky model refinement, increased integration time, and pipeline optimization that will enable future detection of the 21 cm signal through auto-correlation and cross-correlation with optical galaxy surveys.

cosmology: large-scale structure of universe↗

Transverse momentum dependent PDFs at N3LO

We compute the quark and gluon transverse momentum dependent parton distribution functions at next-to-next-to-next-to-leading order (N 3 LO) in perturbative QCD. Our calculation is based on an expansion of the differential Drell-Yan and gluon fusion Higgs production cross sections about their collinear limit. This method allows us to employ cutting edge multiloop techniques for the computation of cross sections to extract these universal building blocks of the collinear limit of QCD. The corresponding perturbative matching kernels for all channels are expressed in terms of simple harmonic polylogarithms up to weight five. As a byproduct, we confirm a previous computation of the soft function for transverse momentum factorization at N 3 LO. Our results are the last missing ingredient to extend the q T subtraction methods to N 3 LO and to obtain resummed q T spectra at N 3 LL' accuracy both for gluon as well as for quark initiated processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computation of potential flows with embedded vortex rings and applications to helicopter rotor wakes

A finite difference scheme for solving the motion of a number of vortex rings is developed. The method is an adaptation of the 'cloud-in-cell' technique to axisymmetric flows, and is thus a combined Eulerian-Lagrangian technique. A straightforward adaptation of the cloud-in-cell scheme to an axisymmetric flow field is shown to introduce a grid dependent self-induced velocity to each vortex ring. To correct this behavior the potential is considered to consist of two parts, a local and a global field. An improved difference formula is derived, allowing the accurate calculation of the potential at points near vortex locations. The local potential is then subtracted before calculating the velocity, leaving only the influences of the remaining vortices. The correct self-induced velocity is then explicitly added to the vortex velocity. Calculations of the motion of one and two vortex rings are performed, demonstrating the ability of the new method to eliminate the grid dependence of the self-induced velocity. The application of the method to the calculation of helicopter rotor flows in hover is attempted.

Roberts, T. W.↗

Analysis and interpretation of MAGSAT anomalies over north Africa

Crustal anomaly detection with MAGSAT data is frustrated by inherent resolving power of the data and by contamination from external and core fields. Quality of the data might be tested by modeling specific tectonic features which produce anomalies that fall within proposed resolution and crustal amplitude capabilities of MAGSAT fields. To test this hypothesis, north African hotspots associated with Ahaggar, Tibesti and Darfur were modeled as magnetic induction anomalies. MAGSAT data were reduced by subtracting external and core fields to isolate scalar and vertical component crustal signals. Of the three volcanic areas, only the Ahaggar region had an associated anomaly of magnitude above error limits of the data. Hotspot hypothesis was tested for Ahaggar by seeing if predicted magnetic signal matched MAGSAT anomaly. Predicted model magnetic signal arising from surface topography of the uplift and the Curie isothermal surface was calculated at MAGSAT altitudes by Fourier transform technique modified to allow for variable magnetization. Curie isotherm surface was calculated using a method for temperature distribution in a moving plate above a fixed hotspot. Magnetic signal was calculated for a fixed plate as well as a number of plate velocities and directions.

Phillips, R. J.↗

The inverse gravimetric problem in gravity modelling

One of the main purposes of geodesy is to determine the gravity field of the Earth in the space outside its physical surface. This purpose can be pursued without any particular knowledge of the internal density even if the exact shape of the physical surface of the Earth is not known, though this seems to entangle the two domains, as it was in the old Stoke's theory before the appearance of Molodensky's approach. Nevertheless, even when large, dense and homogeneous data sets are available, it was always recognized that subtracting from the gravity field the effect of the outer layer of the masses (topographic effect) yields a much smoother field. This is obviously more important when a sparse data set is bad so that any smoothing of the gravity field helps in interpolating between the data without raising the modeling error, this approach is generally followed because it has become very cheap in terms of computing time since the appearance of spectral techniques. The mathematical description of the Inverse Gravimetric Problem (IGP) is dominated mainly by two principles, which in loose terms can be formulated as follows: the knowledge of the external gravity field determines mainly the lateral variations of the density; and the deeper the density anomaly giving rise to a gravity anomaly, the more improperly posed is the problem of recovering the former from the latter. The statistical relation between rho and n (and its inverse) is also investigated in its general form, proving that degree cross-covariances have to be introduced to describe the behavior of rho. The problem of the simultaneous estimate of a spherical anomalous potential and of the external, topographic masses is addressed criticizing the choice of the mixed collection approach.

Sanso, F.↗

In-Situ Reaction Monitoring for Lunar Applications Utilizing a Single Quadrupole Residual Gas Analyzer

With a revived focus to create a permeant presence on the moon and in preparation for future Mars exploration, it is imperative that all re-sources are utilized to their fullest potential. In-Situ Resource Utilization (ISRU) will be critical for future mission success as it would enable independent operation of missions while minimally relying on the complex supply chain created between the Earth, moon and Mars. One of the most critical resources that have been identified for ISRU is the creation of liquid oxygen (O2) for not only breathable air, but also for rocket fuel. On the moon and potentially Mars, one of the most O2 rich resources is in the presence of regolith. Although there are numerous minerals within regolith and various processes to extract O2, this is not the scope of this paper. Since the O2 extraction method can vary, a critical standardized analytical method is needed to verify O2 (or oxygenated precursor compounds) for extraction efficiency and purity. Our group utilizes a modified consumer off-the-shelf (COTS), Residual Gas Analyzer (RGA) Single Quadrupole Mass Spectrometer (QMS). Due to the current field mission requirements, the gas analysis module will need to rely only on mass-to-charge ratios (m/z) and peak intensities to differentiate and quantitate targeted gas-phase reactants or products, i.e. carbon monoxide (CO), carbon dioxide (CO2), methane (CH4) and hydrogen (H2). However, due to solar wind implanted volatiles and various other compounds present in lunar regolith, such as nitrogen (N2), helium (He), and deuterium (D2), deconvolution and quantification of isobaric compounds becomes quite difficult. Specifically, the presence of N2 makes quantification of CO, an oxygenated precursor, difficult due to both compounds having a peak at m/z 28. Approaches to differentiate isobaric components: Initially, the solution to overcome isobaric overlap was to use isotopic ratios, particularly be-tween carbon-13 and nitrogen-15. Although limited literature discuss isotopic composition and attempt to quantify the ratios, variation is particularly high and limited to a few lunar regolith samples. The following were identified as potential approaches to quantify isobaric species, N2 and CO, using our RGA. Threshold Ionization. CO and N2 both have peaks at m/z 28, but only N2 is capable of producing a peak at m/z 7. While a N2 peak at m/z 7 is generally rare, instrument parameters can be optimized in order to maximize the peak intensity. Typically a standard ionization energy of 70 eV is used to monitor a batch process and would produce reliable results for a majority of the reactants or products that we are interested in quantifying. However, to analyze the rare molecular fragment of N2 (N2+ at m/z 7), a different ionization energy may be beneficial. Threshold Ionization Mass Spectrometry (TIMS) will be explored to optimize the N2+ peak at m/z 7 by varying parameters such as ionization energy and pressure. Therefore, a quantitative value at m/z 7 will be used to determine the product gas concentration of N2, which will be subtracted from the total concentration determined at m/z 28 enabling us to approximate the amount of CO. Secondary Detector (infrared gas detector, IR). Another approach to quantitate and differentiate be-tween isobaric components is to leverage alternative modes of detection. Unlike mass spectrometry, which generally uses a separation step prior to analysis (i.e. gas chromatography), IR produces energy that is specifically absorbed depending on the bond and configuration of the molecule. This technique could be leveraged since it is generally only active on non-symmetric molecules, i.e. CO. Utilizing concentrations of CO determined from another detector could help correct concentrations of CO obtained from the QMS. However, the inclusion of another detector adds additional mission requirements, such as power, software changes and avionics updates that add additional cost and schedule requirements.

mass spectrometry↗

Impacts of Land Carbon Flux Variability on Variations in Atmospheric CO2

To better understand the role of atmospheric CO2 in the global carbon cycle, it is important to understand the factors that control its variability in time and space. Here we investigate these factors using a coupled land-atmosphere modeling system (the NASA GEOS-5 model) fitted with land carbon physics and atmospheric CO2 transport. To separate the influences of land carbon flux variability and atmospheric transport variability on the variability of CO2 concentrations, we conducted two carbon-fitted GEOS-5 AGCM simulations run in replay mode (a technique that guides the model's weather to match that of the MERRA-2 reanalysis): (i) a control simulation of 15-year carbon cycle dynamics and climate, and (ii) a simulation in which the climatological seasonal cycles of net biosphere production (NBP), as determined from diagnostics produced in the control experiment, are applied at the surface to the atmosphere instead. The impact of land flux variability on atmospheric CO2 variability is then isolated by subtracting the variability inherent in the second simulation (as induced, e.g., by strong El Nino conditions) from that in the control. These results are also presented in the context of supplemental simulations that have examined the impact of imposed drought in a free running AGCM simulation on carbon fluxes and atmospheric CO2 variability.

Lee, Eunjee↗

Long term frequency stability analysis of the GPS NAVSTAR 6 Cesium clock

Time domain measurements, taken between the NAVSTAR 6 Spacecraft Vehicle (SV) and the Vandenberg Global Positioning System (GPS) Monitor Site, by a pseudo random noise receiver, were collected over an extended period of time and analyzed to estimate the long term frequency stability of the NAVSTAR 6 onboard frequency standard, referenced to the Vandenberg MS frequency standard. The technique employed separates the clock offset from the composite signal by first applying corrections for equipment delays, ionospheric delay, tropospheric delay, Earth rotation and the relativistic effect. The data are edited and smoothed using the predicted SV ephemeris to calculate the geometric delay. Then all available passes from each of the four GPS monitor stations, are collected at 1-week intervals and used to calculate the NAVSTAR orbital elements. The procedure is then completed by subtracting the corrections and the geometric delay, using the final orbital elements, from the composite signal, thus leaving the clock offset and random error.

Mccaskill, T. B.↗

Analysis prcedures for ROSAT XRT/PSPC observations of extended objects and the diffuse background

Proper reduction of ROSAT observations of extended sources or the diffuse background requires techniques that are quite different from the standard point-source analysis. We describe in detail an appropriate set of procedures for this purpose for the X-ray Telescope (XRT), Position Sensitive Proportional Counter (PSPC). We define a standard set of pulse-height bands recommended for analyzing extended structure and describe the use of detector efficiency maps for proper flat-fielding of images. Sources of noncosmic background contamination and their modeling and subtraction are discussed and demonstrated. IDL routines for implementing these data reduction procedures are being made available to the community through the ROSAT IDL library maintained at Goddard Space Flight Center (GSFC).

Snowden, S. L.↗

Additively Manufactured Interdigital Filters for Ultra- Wideband Radar

Abstract— In this paper, we present the development of ultrawideband (UWB) inter-digital cavity filters manufactured with an additive technique. We employed a selective laser melting process with stainless steel via a “flipping and fixturing” approach to achieve dimensional accuracy and frequency selectivity (no postfabrication tuning) comparable to that of conventional micromachining. The filters operate at a center frequency of 3.75 GHz and 6 GHz, respectively, with an instantaneous bandwidth of 2 GHz. This paper presents design and fabrication details, a discussion of measured characteristics for both filters, and a performance comparison with respect to subtractively manufactured parts.

Keywords— Interdigital filters, UWB radar, additiv↗

Direct Detection Doppler Lidar for Spaceborne Wind Measurement

Aerosol and molecular based versions of the double-edge technique can be used for direct detection Doppler lidar spaceborne wind measurement. The edge technique utilizes the edge of a high spectral resolution filter for high accuracy wind measurement using direct detection lidar. The signal is split between an edge filter channel and a broadband energy monitor channel. The energy monitor channel is used for signal normalization. The edge measurement is made as a differential frequency measurement between the outgoing laser signal and the atmospheric backscattered return for each pulse. As a result the measurement is insensitive to laser and edge filter frequency jitter and drift at a level less than a few parts in 10(exp 10). We have developed double edge versions of the edge technique for aerosol and molecular-based lidar measurement of the wind. Aerosol-based wind measurements have been made at Goddard Space Flight Center and molecular-based wind measurements at the University of Geneva. We have demonstrated atmospheric measurements using these techniques for altitudes from 1 to more than 10 km. Measurement accuracies of better than 1.25 m/s have been obtained with integration times from 5 to 30 seconds. The measurements can be scaled to space and agree, within a factor of two, with satellite-based simulations of performance based on Poisson statistics. The theory of the double edge aerosol technique is described by a generalized formulation which substantially extends the capabilities of the edge technique. It uses two edges with opposite slopes located about the laser frequency at approximately the half-width of each edge filter. This doubles the signal change for a given Doppler shift and yields a factor of 1.6 improvement in the measurement accuracy compared to the single edge technique. The use of two high resolution edge filters substantially reduces the effects of Rayleigh scattering on the measurement, as much as order of magnitude, and allows the signal to noise ratio to be substantially improved in areas of low aerosol backscatter. We describe a method that allows the Rayleigh and aerosol components of the signal to be independently determined using the two edge channels and an energy monitor channel. The effects of Rayleigh scattering may then subtracted from the measurement and we show that the correction process does not significantly increase the measurement noise for Rayleigh to aerosol ratios up to 10. We show that for small Doppler shifts a measurement accuracy of 0.4 m/s can be obtained for 5000 detected photon, 1.2 m/s for 1000 detected photons, and 3.7 m/s for 50 detected photons for a Rayleigh to aerosol ratio of 5. Methods for increasing the dynamic range of the aerosol-based system to more than +/- 100 m/s are given.

Korb, C. Laurence↗

Identifying Transient Candidates in the Dark Energy Survey Using Convolutional Neural Networks

The ability to discover new transient candidates via image differencing without direct human intervention is an important task in observational astronomy. For these kind of image classification problems, machine learning techniques such as Convolutional Neural Networks (CNNs) have shown remarkable success. In this work, we present the results of an automated transient candidate identification on images with CNNs for an extant data set from the Dark Energy Survey Supernova program, whose main focus was on using Type Ia supernovae for cosmology. By performing an architecture search of CNNs, we identify networks that efficiently select non-artifacts (e.g., supernovae, variable stars, AGN, etc.) from artifacts (image defects, mis-subtractions, etc.), achieving the efficiency of previous work performed with random Forests, without the need to expend any effort in feature identification. The CNNs also help us identify a subset of mislabeled images. Performing a relabeling of the images in this subset, the resulting classification with CNNs is significantly better than previous results, lowering the false positive rate by 27% at a fixed missed detection rate of 0.05.

79 ASTRONOMY AND ASTROPHYSICS↗

Waves in Radial Gravity Using Magnetic Fluid

Terrestrial laboratory experiments studying various fluid dynamical processes are constrained, by being in an Earth laboratory, to have a gravitational body force which is uniform and unidirectional. Therefore fluid free-surfaces are horizontal and flat. Such free surfaces must have a vertical solid boundary to keep the fluid from spreading horizontally along a gravitational potential surface. In atmospheric, oceanic, or stellar fluid flows that have a horizontal scale of about one-tenth the body radius or larger, sphericity is important in the dynamics. Further, fluids in spherical geometry can cover an entire domain without any sidewall effects, i.e. have truly periodic boundary conditions. We describe spherical body-force laboratory experiments using ferrofluid. Ferrofluids are dilute suspensions of magnetic dipoles, for example magnetite particles of order 10 nm diameter, suspended in a carrier fluid. Ferrofluids are subject to an additional body force in the presence of an applied magnetic field gradient. We use this body force to conduct laboratory experiments in spherical geometry. The present study is a laboratory technique improvement. The apparatus is cylindrically axisymmetric. A cylindrical ceramic magnet is embedded in a smooth, solid, spherical PVC ball. The geopotential field and its gradient, the body force, were made nearly spherical by careful choice of magnet height-to-diameter ratio and magnet size relative to the PVC ball size. Terrestrial gravity is eliminated from the dynamics by immersing the "planet" and its ferrofluid "ocean" in an immiscible silicone oil/freon mixture of the same density. Thus the earth gravity is removed from the dynamics of the ferrofluid/oil interface and the only dynamically active force there is the radial magnetic gravity. The entire apparatus can rotate, and waves are forced on the ferrofluid surface by exterior magnets. The biggest improvement in technique is in the wave visualization. Fluorescing dye is added to the oil/freon mixture and an argon ion laser generates a horizontal light that can be scanned vertically. Viewed from above, the experiment is a black circle with wave deformations surrounded by a light background. A contour of the image intensity at any light sheet position gives the surface of the ferrofluid "ocean" at that "latitude". Radial displacements of the waves as a function of longitude are obtained by subtracting the contour line positions from a no-motion contour at that laser sheet latitude. The experiments are run by traversing the forcing magnet with the laser sheet height fixed and images are frame grabbed to obtain a time-series at one latitude. The experiment is then re-run with another laser-sheet height to generate a full picture of the three-dimensional wave structure in the upper hemisphere of the ball as a function of time. We concentrate here on results of laboratory studies of waves that are important in Earth's atmosphere and especially the ocean. To get oceanic scaling in the laboratory, the experiment must rotate rapidly (4-second rotation period) so that the wave speed is slow compared to the planetary rotation speed as in the ocean. In the Pacific Ocean, eastward propagating Kelvin waves eventually run into the South American coast. Theory predicts that some of the wave energy should scatter into coastal-trapped Kelvin waves that propagate north and south along the coast. Some of this coastal wave energy might then scatter into mid-latitude Rossby waves that propagate back westward. Satellite observations of the Pacific Ocean sea-surface temperature and height seem to show signatures of westward propagating mid-latitude Rossby waves, 5 to 10 years after the 1982-83 El Nino. The observational data is difficult to interpret unambiguously owing to the large range of motions that fill the ocean at shorter timescales. This series of reflections giving eastward, north- ward, and then westward traveling waves is observed cleanly in the laboratory experiments, confirming the theoretical expectations.

Ohlsen, D. R.↗

General bounds on holographic complexity

We prove a positive volume theorem for asymptotically AdS spacetimes: the maximal volume slice has nonnegative vacuum-subtracted volume, and the vacuum-subtracted volume vanishes if and only if the spacetime is identically pure AdS. Under the Complexity=Volume proposal, this constitutes a positive holographic complexity theorem. The result features a number of parallels with the positive energy theorem, including the assumption of an energy condition that excludes false vacuum decay (the AdS weak energy condition). Our proof is rigorously established in broad generality in four bulk dimensions, and we provide strong evidence in favor of a generalization to arbitrary dimensions. Our techniques also yield a holographic proof of Lloyd’s bound for a class of bulk spacetimes. We further establish a partial rigidity result for wormholes: wormholes with a given throat size are more complex than AdS-Schwarzschild with the same throat size.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Optimization of foreground moment deprojection for semi-blind CMB polarization reconstruction

Abstract Upcoming Cosmic Microwave Background (CMB) experiments, aimed at measuring primordial CMB polarization B-modes, require exquisite control of instrumental systematics and Galactic foreground contamination. Blind minimum-variance techniques, like the Needlet Internal Linear Combination (NILC), have proven effective in reconstructing the CMB polarization signal and mitigating foregrounds and systematics across diverse sky models without suffering from foreground mismodelling errors. Still, residual foreground contamination from NILC may bias the recovered CMB polarization at large angular scales when confronted with the most complex foreground scenarios.By adding constraints to NILC to deproject statistical moments of the Galactic emission, the Constrained Moment ILC (cMILC) method has been demonstrated to further enhance foreground subtraction, albeit with an associated increase in overall noise variance. Faced with this trade-off between foreground bias reduction and overall variance minimization, there is still no recipe on which moments to deproject and which are better suited for blind variance minimization. To address this, we introduce the optimized cMILC (ocMILC) pipeline, which performs full automated optimization of the required number and set of foreground moments to deproject, pivot parameter values, and deprojection coefficients across the sky and angular scales, depending on the actual sky complexity, available frequency coverage, and experiment sensitivity. The optimal number of moments for deprojection, before paying significant noise penalty, is determined through a data diagnosis inspired by the Generalized NILC (GNILC) method.Validated on B-mode simulations of thePICOspace mission concept with four challenging foreground models, ocMILC exhibits lower Galactic foreground contamination compared to NILC and cMILC at all angular scales, with limited noise penalty. This multi-layer optimization enables the ocMILC pipeline to achieve unbiased posteriors of the tensor-to-scalar ratio, regardless of foreground complexity.

Astronomy & Astrophysics↗