Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bias correction”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The most luminous, merger-free AGNs show only marginal correlation with bar presence

The role of large-scale bars in the fuelling of active galactic nuclei (AGNs) is still debated, even as evidence mounts that black hole growth in the absence of galaxy mergers cumulatively dominates and may substantially influence disc (i.e. merger-free) galaxy evolution. We investigate whether large-scale galactic bars are a good candidate for merger-free AGN fuelling. Specifically, we combine slit spectroscopy and Hubble Space Telescope imagery to characterize star formation rates (SFRs) and stellar masses of the unambiguously disc-dominated host galaxies of a sample of luminous, Type-1 AGN with 0.02 < z < 0.24. After carefully correcting for AGN signal, we find no clear difference in SFRs between AGN hosts and a stellar mass-matched sample of galaxies lacking an AGN (0.013 < z < 0.19), although this could be due to small sample size (nAGN = 34). We correct for SFR and stellar mass to minimize selection biases, and compare the bar fraction in the two samples. We find that AGNs are marginally (~1.7σ) more likely to host a bar than inactive galaxies, with AGN hosts having a bar fraction, $f_{\mathrm{bar}}=0.59^{+0.08}_{-0.09}$ and inactive galaxies having a bar fraction, $f_{\mathrm{bar}}=0.44^{+0.08}_{-0.09}$. However, we find no further differences between SFR- and mass-matched AGNs and inactive samples. While bars could potentially trigger AGN activity, they appear to have no further, unique effect on a galaxy’s stellar mass or SFR.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Planet Occurrence Rate Correlated to Stellar Dynamical History: Evidence from Kepler and Gaia

The dynamical history of stars influences the formation and evolution of planets significantly. To explore the influence of dynamical history on the planet formation and evolution using observations, we assume stars that experienced significantly different dynamical histories tend to have different relative velocities. Utilizing the accurate Gaia–Kepler Stellar Properties Catalog, we select single main-sequence stars and divide these stars into three groups according to their relative velocities, i.e., high-V, medium-V, and low-V stars. After considering the known biases from Kepler data and adopting prior and posterior correction to minimize the influence of stellar properties on planet occurrence rate, we find that high-V stars have a lower occurrence rate of super-Earths and sub-Neptunes (1–4 R {sub ⊕}, P < 100 days) and a higher occurrence rate of sub-Earth (0.5–1 R {sub ⊕}, P < 30 days) than low-V stars. Additionally, high-V stars have a lower occurrence rate of hot Jupiter-sized planets (4–20 R {sub ⊕}, P < 10 days) and a slightly higher occurrence rate of warm or cold Jupiter-sized planets (4–20 R {sub ⊕}, 10 < P < 400 days). After investigating multiplicity and eccentricity, we find that high-V planet hosts prefer a higher fraction of multiplanet systems and lower average eccentricity, which are consistent with the eccentricity–multiplicity dichotomy of Kepler planetary systems. All of these statistical results favor the scenario that high-V stars with large relative velocity may experience fewer gravitational events, while low-V stars may be influenced by stellar clustering significantly.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Type Ia Supernova Growth-rate Measurement with LSST Simulations: Intrinsic Scatter Systematics

Measurement of the growth rate of structures (fσ 8 ) with Type Ia supernovae (SNe Ia) will improve our understanding of the nature of dark energy and enable tests of general relativity. In this paper, we generate simulations of the 10 yr SN Ia data set of the Rubin-LSST survey, including a correlated velocity field from an N-body simulation and realistic models of SNe Ia properties and their correlations with host-galaxy properties. We find, similar to SN Ia analyses that constrain the dark energy equation-of-state parameters w 0 w a , that constraints on fσ 8 can be biased depending on the intrinsic scatter of SNe Ia. While for the majority of intrinsic scatter models we recover fσ 8 with a precision of ∼13%–14%, for the most realistic dust-based model, we find that the presence of non-Gaussianities in Hubble diagram residuals leads to a bias on fσ 8 of ∼ −20%. When trying to correct for the dust-based intrinsic scatter, we find that the propagation of the uncertainty on the model parameters does not significantly increase the error on fσ 8 . We also find that while the main component of the error budget of fσ 8 is the statistical uncertainty (>75% of the total error budget), the systematic error budget is dominated by the uncertainty on the damping parameter, σ u , that gives an empirical description of the effect of redshift space distortions on the velocity power spectrum. Our results motivate a search for new methods to correct for the non-Gaussian distribution of the Hubble diagram residuals, as well as an improved modeling of the damping parameter.

Carreres, Bastien [Duke Univ., Durham, NC (United ↗
DOE OSTI2021

Simulation of radon-222 with the GEOS-Chem global model: Emissions, seasonality, and convective transport

Radon-222 (222Rn) is a short-lived radioactive gas naturally emitted from land surface, and has long been used to assess convective transport in atmospheric models. In this study, we simulate 222Rn using the GEOS-Chem chemical transport model with aims to improve our understanding of 222Rn emissions and surface concentration seasonality, and to characterize convective transport associated with two Goddard Earth Observing System (GEOS) meteorological products, MERRA and GEOS-FP. We evaluate four available 222Rn emission scenarios by comparing model results with surface observations at 51 global surface sites. The default emission scenario in GEOS-Chem yields a moderate agreement with global surface observations (< 70% data within a factor of 2) and a large underestimate of wintertime surface 222Rn concentrations at Northern Hemisphere mid- and high-latitudes due to an oversimplified formulation of 222Rn emission fluxes (1 atom cm-2 s-1 over land with a reduction by a factor of 3 under freezing conditions). We compose a new global 222Rn emission scenario based on Zhang et al. (2011) and show its potential to improve simulated surface 222Rn concentrations and seasonality. The regional components of this emission scenario include spatially and temporally varying emission fluxes derived from previous measurements of soil radium content and soil exhalation models, which are key factors to 222Rn emission flux rates. However, large model underestimates of surface 222Rn concentrations still exist in Asia, suggesting unusually high regional 222Rn emissions. We propose a conservative up-scaling factor of 1.2 for 222Rn emission fluxes in China, as also constrained by the observed deposition fluxes of 210Pb (decay daughter of 222Rn). With this modification, the model shows better agreement with the observations in Europe and North America (>80% data within a factor of 2), and reasonable agreement in Asia (close to 70%). Further constraints on 222Rn emissions would require additional observations of surface 222Rn concentrations and emission fluxes in central U.S., Canada, Africa, and Asia. We also compare and assess convective transport in model simulations driven by MERRA and GEOS-FP using observed 222Rn vertical profiles during northern mid-latitude summertime and from two short-term airborne campaigns. While the simulations with both GEOS products are able to capture the observed vertical gradient of 222Rn concentrations in the lower troposphere (0-4 km), neither correctly represents the level of convective detrainment, resulting in biases in the middle and upper troposphere. Compared to GEOS-FP, MERRA leads to stronger convective transport of 222Rn, which is partially compensated by its weaker large-scale vertical advection, resulting in similar global vertical distributions of 222Rn concentrations between the two simulations.

Zhang, Bo↗
DOE OSTI2020

Biases in genome reconstruction from metagenomic data

Background Advances in sequencing, assembly, and assortment of contigs into species-specific bins has enabled the reconstruction of genomes from metagenomic data (MAGs). Though a powerful technique, it is difficult to determine whether assembly and binning techniques are accurate when applied to environmental metagenomes due to a lack of complete reference genome sequences against which to check the resulting MAGs. Methods We compared MAGs derived from an enrichment culture containing ~20 organisms to complete genome sequences of 10 organisms isolated from the enrichment culture. Factors commonly considered in binning software—nucleotide composition and sequence repetitiveness—were calculated for both the correctly binned and not-binned regions. This direct comparison revealed biases in sequence characteristics and gene content in the not-binned regions. Additionally, the composition of three public data sets representing MAGs reconstructed from the Tara Oceans metagenomic data was compared to a set of representative genomes available through NCBI RefSeq to verify that the biases identified were observable in more complex data sets and using three contemporary binning software packages. Results Repeat sequences were frequently not binned in the genome reconstruction processes, as were sequence regions with variant nucleotide composition. Genes encoded on the not-binned regions were strongly biased towards ribosomal RNAs, transfer RNAs, mobile element functions and genes of unknown function. Our results support genome reconstruction as a robust process and suggest that reconstructions determined to be >90% complete are likely to effectively represent organismal function; however, population-level genotypic heterogeneity in natural populations, such as uneven distribution of plasmids, can lead to incorrect inferences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Joint inference of multiplicative and additive systematics in galaxy density fluctuations and clustering measurements

Galaxy clustering measurements are a key probe of the matter density field in the Universe. With the era of precision cosmology upon us, surveys rely on precise measurements of the clustering signal for meaningful cosmological analysis. However, the presence of systematic contaminants can bias the observed galaxy number density, and thereby bias the galaxy two-point statistics. As the statistical uncertainties get smaller, correcting for these systematic contaminants becomes increasingly important for unbiased cosmological analysis. We present and validate a new method for understanding and mitigating both additive and multiplicative systematics in galaxy clustering measurements (two-point function) by joint inference of contaminants in the galaxy overdensity field (one-point function) using a maximum-likelihood estimator (MLE). We test this methodology with Kilo-Degree Survey-like mock galaxy catalogues and synthetic systematic template maps. We estimate the cosmological impact of such mitigation by quantifying uncertainties and possible biases in the inferred relationship between the observed and the true galaxy clustering signal. Our method robustly corrects the clustering signal to the sub-percent level and reduces numerous additive and multiplicative systematics from 1.5σ to less than 0.1σ for the scenarios we tested. In addition, we provide an empirical approach to identifying the functional form (additive, multiplicative, or other) by which specific systematics contaminate the galaxy number density. Even though this approach is tested and geared towards systematics contaminating the galaxy number density, the methods can be extended to systematics mitigation for other two-point correlation measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Version 4 CALIPSO Imaging Infrared Radiometer ice and liquid water cloud microphysical properties – Part I: The retrieval algorithms

Following the release of the version 4 Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) data products from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, a new version (version 4; V4) of the CALIPSO Imaging Infrared Radiometer (IIR) Level 2 data products has been developed. The IIR Level 2 data products include cloud effective emissivities and cloud microphysical properties such as effective diameter and ice or liquid water path estimates. Dedicated retrievals for water clouds were added in V4, taking advantage of the high sensitivity of the IIR retrieval technique to small particle sizes. This paper (Part I) describes the improvements in the V4 algorithms compared to those used in the version 3 (V3) release, while results will be presented in a companion (Part II) paper. The IIR Level 2 algorithm has been modified in the V4 data release to improve the accuracy of the retrievals in clouds of very small (close to 0) and very large (close to 1) effective emissivities. To reduce biases at very small emissivities that were made evident in V3, the radiative transfer model used to compute clear-sky brightness temperatures over oceans has been updated and tuned for the simulations using Modern-Era Retrospective analysis for Research and Applications version 2 (MERRA-2) data to match IIR observations in clear-sky conditions. Furthermore, the clear-sky mask has been refined compared to V3 by taking advantage of additional information now available in the V4 CALIOP 5 km layer products used as an input to the IIR algorithm. After sea surface emissivity adjustments, observed and computed brightness temperatures differ by less than ±0.2 K at night for the three IIR channels centered at 08.65, 10.6, and 12.05 µm, and inter-channel biases are reduced from several tens of Kelvin in V3 to less than 0.1 K in V4. We have also improved retrievals in ice clouds having large emissivity by refining the determination of the radiative temperature needed for emissivity computation. The initial V3 estimate, namely the cloud centroid temperature derived from CALIOP, is corrected using a parameterized function of temperature difference between cloud base and top altitudes, cloud absorption optical depth, and CALIOP multiple scattering correction factor. As shown in Part II, this improvement reduces the low biases at large optical depths that were seen in V3 and increases the number of retrievals. As in V3, the IIR microphysical retrievals use the concept of microphysical indices applied to the pairs of IIR channels at 12.05 and 10.6 µm and at 12.05 and 08.65 µm. The V4 algorithm uses ice look-up tables (LUTs) built using two ice habit models from the recent “TAMUice2016” database, namely the single-hexagonal-column model and the eight-element column aggregate model, from which bulk properties are synthesized using a gamma size distribution. Four sets of effective diameters derived from a second approach are also reported in V4. Here, the LUTs are analytical functions relating microphysical index applied to IIR channels 12.05 and 10.6 µm and effective diameter as derived from in situ measurements at tropical and midlatitudes during the Tropical Composition, Cloud, and Climate Coupling (TC4) and Small Particles in Cirrus Science and Operations Plan (SPARTICUS) field experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Universal scaling of weak localization in graphene due to bias-induced dispersion decoherence

The differential conductance of graphene is shown to exhibit a zero-bias anomaly at low temperatures, arising from a suppression of the quantum corrections due to weak localization and electron interactions. A simple rescaling of these data, free of any adjustable parameters, shows that this anomaly exhibits a universal, temperature- (T) independent form. According to this, the differential conductance is approximately constant at small voltages (V < kBT/e), while at larger voltages it increases logarithmically with the applied bias. For theoretical insight into the origins of this behaviour, which is inconsistent with electron heating, we formulate a model for weak-localization in the presence of nonequilibrium transport. According to this model, the applied voltage causes unavoidable dispersion decoherence, which arises as diffusing electron partial waves, with a spread of energies defined by the value of the applied voltage, gradually decohere with one another as they diffuse through the system. The decoherence yields a universal scaling of the conductance as a function of eV/kBT, with a logarithmic variation for eV/kBT > 1, variations in accordance with the results of experiment. Our theoretical description of nonequilibrium transport in the presence of this source of decoherence exhibits strong similarities with the results of experiment, including the aforementioned rescaling of the conductance and its logarithmic variation as a function of the applied voltage.

42 ENGINEERING↗
DOE OSTI2024

Absorptive corrections to the electromagnetic form factor in high-energy elastic proton-proton scattering

Recently, it was noted that absorptive corrections to the electromagnetic form factor in high-energy proton-proton scattering are important for the theoretical interpretation of the p ↑ p and p ↑ A analyzing power A N ( t ) measurements with the Hydrogen Jet Target polarimeter (HJET) at RHIC. Here, a concise expression for the absorptive correction was derived within the eikonal approach. The resulting analysis reveals a systematic bias, nearly independent of the beam energy, in the experimental determination of the real-to-imaginary ratio ρ when absorption effects are overlooked in the data analysis. Quantification of this bias, as ρ meas = ρ + ( 0.036 ± 0.016 ) bias , was achieved using a Regge fit applied to available proton-proton measurements of ρ meas ( s ) and σ tot meas ( s ) . Considering the potential impact of such an effect on the experimentally determined A N ( t ) , one may enhance consistency between the HJET and STAR measurements of the hadronic spin-flip amplitude. While the sign of the bias in the value of ρ aligns with the anticipated effective increase in the proton charge radius in p p scattering due to absorption, it amplifies the observed discrepancy between σ tot meas and ρ meas values at s = 13 TeV as measured in the TOTEM experiment. Evaluation (using published TOTEM data) of the measured proton-proton d σ / d t dependence on the absorptive corrections indicated that possible soft photon corrections to the hadronic amplitude slope may be essential for such data analysis. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Dark Energy Survey Year 6 Results: improved mitigation of spatially varying observational systematics with masking

As photometric surveys reach unprecedented statistical precision, systematic uncertainties increasingly dominate large-scale structure probes relying on galaxy number density. Defining the final survey footprint is critical, as it excludes regions affected by artefacts or suboptimal observing conditions. For galaxy clustering, spatially varying observational systematics, such as seeing, are a leading source of bias. Template maps of contaminants are used to derive spatially dependent corrections, but extreme values may fall outside the applicability range of mitigation methods, compromising correction reliability. The complexity and accuracy of systematics modelling depend on footprint conservativeness, with aggressive masking enabling simpler, robust mitigation. We present a unified approach to define the DES Year 6 joint footprint, integrating observational systematics templates and artefact indicators that degrade mitigation performance. This removes extreme values from an initial seed footprint, leading to the final joint footprint. By evaluating the DES Year 6 lens sample MagLim++ plus plus on this footprint, we enhance the Iterative Systematics Decontamination (ISD) method, detecting non-linear systematic contamination and improving correction accuracy. While the mask's impact on clustering is less significant than systematics decontamination, it remains non-negligible, comparable to statistical uncertainties in certain w(theta) scales and redshift bins. Supporting coherent analyses of galaxy clustering and cosmic shear, the final footprint spans 4031.04 deg2, setting the basis for DES Year 6 1x2pt, 2x2pt, and 3x2pt analyses. This work highlights how targeted masking strategies optimise the balance between statistical power and systematic control in Stage-III and -IV surveys.

Rodríguez-Monroy, M. [Madrid, IFT; IJCLab, Orsay]↗
DOE OSTI2022

In-flight polarization angle calibration for LiteBIRD: blind challenge and cosmological implications

We present a demonstration of the in-flight polarization angle calibration for the JAXA/ISAS second strategic large class mission, LiteBIRD, and estimate its impact on the measurement of the tensor-to-scalar ratio parameter, r, using simulated data. We generate a set of simulated sky maps with CMB and polarized foreground emission, and inject instrumental noise and polarization angle offsets to the 22 (partially overlapping) LiteBIRD frequency channels. Our in-flight angle calibration relies on nulling the EB cross correlation of the polarized signal in each channel. This calibration step has been carried out by two independent groups with a blind analysis, allowing an accuracy of the order of a few arc-minutes to be reached on the estimate of the angle offsets. Both the corrected and uncorrected multi-frequency maps are propagated through the foreground cleaning step, with the goal of computing clean CMB maps. We employ two component separation algorithms, the Bayesian-Separation of Components and Residuals Estimate Tool (B-SeCRET), and the Needlet Internal Linear Combination (NILC). We find that the recovered CMB maps obtained with algorithms that do not make any assumptions about the foreground properties, such as NILC, are only mildly affected by the angle miscalibration. However, polarization angle offsets strongly bias results obtained with the parametric fitting method. Once the miscalibration angles are corrected by EB nulling prior to the component separation, both component separation algorithms result in an unbiased estimation of the r parameter. While this work is motivated by the conceptual design study for LiteBIRD, its framework can be broadly applied to any CMB polarization experiment. In particular, the combination of simulation plus blind analysis provides a robust forecast by taking into account not only detector sensitivity but also systematic effects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Path-Integrated X-Ray Digital Image Correlation using Synthetic Reference Images

X-rays can provide images when an object is visibly obstructed, allowing for motion measurements via x-ray digital image correlation (DIC). However, x-ray images are path-integrated and contain data for all objects between the source and detector. If multiple objects are present in the x-ray path, conventional DIC algorithms may fail to correlate the x-ray images. A new DIC algorithm called path-integrated (PI)-DIC addresses this issue by reformulating the matching criterion for DIC to account for multiple, independently-moving objects. PI-DIC requires a set of reference x-ray images of each independent object. However, due to experimental constraints, such reference images might not be obtainable from the experiment. Here, this work focuses on the reliability of synthetically-generated reference images, in such cases. A simplified exemplar is used for demonstration purposes, consisting of two aluminum plates with tantalum x-ray DIC patterns undergoing independent rigid translations. Synthetic reference images based on the “as-designed” DIC patterns were generated. However, PI-DIC with the synthetic images suffered some biases due to manufacturing defects of the patterns. A systematic study of seven identified defect types found that an incorrect feature diameter was the most influential defect. Synthetic images were re-generated with the corrected feature diameter, and PI-DIC errors were improved by a factor of 3-4. Final biases ranged from 0.00-0.04 px, and standard uncertainties ranged from 0.06-0.11 px. In conclusion, PI-DIC accurately measured the independent displacement of two plates from a single series of path-integrated x-ray images using synthetically-generated reference images, and the methods and conclusions derived here can be extended to more generalized cases involving stereo PI-DIC for arbitrary specimen geometry and motion. This work thus extends the application space of x-ray imaging for full-field DIC measurements of multiple surfaces or objects in extreme environments where optical DIC is not possible.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A Non‐Intrusive Machine Learning Framework for Debiasing Long‐Time Coarse Resolution Climate Simulations and Quantifying Rare Events Statistics

Abstract Due to the rapidly changing climate, the frequency and severity of extreme weather is expected to increase over the coming decades. As fully‐resolved climate simulations remain computationally intractable, policy makers must rely on coarse‐models to quantify risk for extremes. However, coarse models suffer from inherent bias due to the ignored “sub‐grid” scales. We propose a framework to non‐intrusively debias coarse‐resolution climate predictions using neural‐network (NN) correction operators. Previous efforts have attempted to train such operators using loss functions that match statistics. However, this approach falls short with events that have longer return period than that of the training data, since the reference statistics have not converged. Here, the scope is to formulate a learning method that allows for correction of dynamics and quantification of extreme events with longer return period than the training data. The key obstacle is the chaotic nature of the underlying dynamics. To overcome this challenge, we introduce a dynamical systems approach where the correction operator is trained using reference data and a coarse model simulation nudged toward that reference. The method is demonstrated on debiasing an under‐resolved quasi‐geostrophic model and the Energy Exascale Earth System Model (E3SM). For the former, our method enables the quantification of events that have return period two orders longer than the training data. For the latter, when trained on 8 years of ERA5 data, our approach is able to correct the coarse E3SM output to closely reflect the 36‐year ERA5 statistics for all prognostic variables and significantly reduce their spatial biases.

Barthel Sorensen, B.↗
DOE OSTI2023

A Consistent Representation of Cloud Overlap and Cloud Subgrid Vertical Heterogeneity

Many global climate models underestimate the cloud cover and overestimate the cloud albedo, especially for low-level clouds. We determine how a correct representation of the vertical structure of clouds can fix part of this bias. We use the 1D McICA framework and focus on low-level clouds. Using Large Eddy Simulations results as reference, we propose a method based on exponential-random overlap that represents the cloud overlap between layers and the subgrid cloud properties over several vertical scales, with a single value of the overlap parameter. Starting from a coarse vertical grid, representative of atmospheric models, this algorithm is used to generate the vertical profile of the cloud fraction with a finer vertical resolution, or to generate it on the coarse grid but with subgrid heterogeneity and cloud overlap that ensures a correct cloud cover. Doing so we find decorrelation lengths are dependent on the vertical resolution, except if the vertical subgrid heterogeneity and interlayer overlap are taken into account coherently. We confirm that the frequently used maximum-random overlap leads to a significant error by underestimating the low-level cloud cover with a relative error of about 50%, that can lead to an error of SW cloud albedo as big as 70%. Not taking into account the subgrid vertical heterogeneity of clouds can cause a relative error of 20% in brightness, assuming the cloud cover is correct.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Weak lensing mass-richness relation of redMaPPer clusters in LSST DESC DC2 simulations

Cluster scaling relations are key ingredients in cluster abundance-based cosmological studies. In optical cluster cosmology, where clusters are detected through their richness, cluster-weak gravitational lensing has proven to be a powerful tool to constrain the cluster mass-richness relation. This work is conducted as part of the Dark Energy Science Collaboration (DESC), which aims to analyze the Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory, starting in 2026. Cluster properties inferred from weak lensing, such as mass, suffer from several sources of bias. In this paper, we aim to test the impact of modeling choices and observational systematics in cluster lensing on the inference of the mass-richness relation. We constrained the mass-richness relation of 3600 clusters detected by the redMaPPer algorithm in the cosmoDC2 extragalactic mock catalog of the LSST DESC DC2 simulation, covering 440 deg 2 , using number count measurements and either stacked weak lensing profiles or mean cluster masses in several intervals of richness (20 ≤ λ ≤ 200) and redshift (0.2 ≤ z ≤ 1). We provide the first constraints on the redMaPPer cluster mass-richness relation detected in cosmoDC2. We find that for an LSST-like source galaxy density, our constraints are robust to changes in the concentration-mass relation, as well as the dark matter density profile modeling choices, when source redshifts and shapes are perfectly known. We find that photometric redshift uncertainties can introduce bias at the 1σ level, which could be mitigated by an overall correction factor fitted jointly with the scaling parameters. We find that including positive shear-richness covariance in the fit shifts the results by up to 0.5σ. Our constraints also offer a fair comparison to a fiducial mass-richness relation, obtained from matching cosmoDC2 halo masses to redMaPPer-detected cluster richness results.

galaxy clusters↗
DOE OSTI2023

Real-world CMB lensing quadratic estimator power spectrum response

I describe a method to estimate response matrices of Cosmic Microwave Background (CMB) lensing power spectra estimators to the true sky power under realistic conditions. Applicable to all lensing reconstruction pipelines based on quadratic estimators (QE), it uses a small number of Gaussian CMB Monte-Carlos and specially designed QE's in order to obtain sufficiently accurate matrices with little computational effort. This method may be used to improve the modelling of CMB lensing band-powers by incorporating at least some of the non-idealities encountered in CMB lensing reconstruction. These non-idealities always include masking, and often inhomogeneous filtering, either in the harmonic domain or pixel space. I obtain these matrices for Planck latest lensing reconstructions, and then show that the residual couplings induced by masking explain very well the residual multiplicative bias seen on the Planck simulations, removing the need for an empirical correction.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Electronic band gaps from quantum Monte Carlo methods

We develop a method for calculating the fundamental electronic gap of semiconductors and insulators using grand canonical quantum Monte Carlo simulations. Further, we discuss the origin of the bias introduced by supercell calculations of finite size and show how to correct the leading and subleading finite size errors either based on observables accessible in the finite-sized simulations or from density-functional theory calculations. Our procedure is applied to carbon, silicon, and molecular hydrogen crystals, and compared to experiment for carbon and silicon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Calibration and Validation of the SAIL Radar

Calibration of a weather radar is not a simple process, as is well described and documented in several articles (Chandrasekar, V, L Baldini, N Bharadwaj, and PL Smith. 2015. "Calibration procedures for global precipitation-measurement ground-validation radars." URSI Radio Science Bulletin 2015(355): 45–73, https://10.23919/URSIRSB.2015.7909473.). Direct hardware-based characterization is called calibration whereas indirect check with other sources such as disdrometer, or comparison with other radars or indirect means, is called validation. The validation process has its own measurement error and the comparison should fall within its limit. In the following we describe calibration and validation processes for the SAIL radar. Based on the comparison of calibration and validation, an addition of 2 dB is suggested for the SAIL radar observations, with a nuanced Z dr calibration correction, as presented. If a fixed Z dr calibration is preferred, a bias addition of 0.5 dB is suggested.

54 ENVIRONMENTAL SCIENCES↗