SEARCH · Engineering Papers
Results for “Systematic error”
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.
Universal framework for simultaneous tomography of quantum states and SPAM noise
We present a general denoising algorithm for performing simultaneous tomography of quantum states and measurement noise. This algorithm allows us to fully characterize state preparation and measurement (SPAM) errors present in any quantum system. Our method is based on the analysis of the properties of the linear operator space induced by unitary operations. Given any quantum system with a noisy measurement apparatus, our method can output the quantum state and the noise matrix of the detector up to a single gauge degree of freedom. We show that this gauge freedom is unavoidable in the general case, but this degeneracy can be generally broken using prior knowledge on the state or noise properties, thus fixing the gauge for several types of state-noise combinations with no assumptions about noise strength. Such combinations include pure quantum states with arbitrarily correlated errors, and arbitrary states with block independent errors. This framework can further use available prior information about the setting to systematically reduce the number of observations and measurements required for state and noise detection. Our method effectively generalizes existing approaches to the problem, and includes as special cases common settings considered in the literature requiring an uncorrelated or invertible noise matrix, or specific probe states.
Validation of the DESI-DR1 3x2-pt analysis: scale cut and shear ratio tests
Combined survey analyses of galaxy clustering and weak gravitational lensing (3x2-pt studies) will allow new and accurate tests of the standard cosmological model. However, careful validation is necessary to ensure that these cosmological constraints are not biased by uncertainties associated with the modelling of astrophysical or systematic effects. In this study we validate the combined 3x2-pt analysis of the Dark Energy Spectroscopic Instrument Data Release 1 (DESI-DR1) spectroscopic galaxy clustering and overlapping weak lensing datasets from the Kilo-Degree Survey (KiDS), the Dark Energy Survey (DES), and the Hyper-Suprime-Cam Survey (HSC). By propagating the modelling uncertainties associated with the non-linear matter power spectrum, non-linear galaxy bias and baryon feedback, we design scale cuts to ensure that measurements of the matter density and the amplitude of the matter power spectrum are biased by less than 30% of the statistical error. We also test the internal consistency of the data and weak lensing systematics by performing new measurements of the lensing shear ratio. We demonstrate that the DESI-DR1 shear ratios can be successfully fit by the same model used to describe cosmic shear correlations, and analyse the additional information that can be extracted about the source redshift distributions and intrinsic alignment parameters. This study serves as crucial preparation for the upcoming cosmological parameter analysis of these datasets.
SU-E-P-36: Evaluation of MLC Positioning Errors in Dynamic IMRT Treatments by Analyzing Dynalog Files
Purpose: To assess the accuracy of MLC positioning in Varian linear accelerator, in dynamic IMRT technique, from the analysis of dynalog files generated by the MLC controller. Methods: In Clinac accelerators (pre-TrueBeam technology), control system has an approximately 50ms delay (one control cycle time). Then, the system compares the measured position to the planned position corresponding to the next control cycle. As it has been confirmed by Varian technical support, this effect causes that measured positions appear in dynalogs one cycle out of phase with respect to the planned positions. Around 9000 dynalogs have been analyzed, coming from the three linear accelerators of our center (one Trilogy and two Clinac 21EX) equipped with a Millennium 120 MLC. In order to compare our results to recent publications, leaf positioning errors (RMS and 95th percentile) are calculated with and without delay effect. Dynalogs have been analyzed using a in-house Matlab software. Results: The RMS errors were 0.341, 0.339 and 0.348mm for each Linac; being the average error 0.343 mm. The 95th percentiles of the error were 0.617, 0.607 and 0.625; with an average of 0.617mm. A recent multi-institution study carried out by Kerns et al. found a mean leaf RMS error of 0.32mm and a 95th percentile error value of 0.64mm.Without delay effect, mean leaf RMS errors obtained were 0.040, 0.042 and 0.038mm for each treatment machine; being the average 0.040mm. The 95th percentile error values obtained were 0.057, 0.058 and 0.054 mm, with an average of 0.056mm. Conclusion: Results obtained for the mean leaf RMS error and the mean 95th percentile were consistent with the multi-institution study. Calculated error statistics with delay effect are significantly larger due to the speed proportional and systematic leaf offset. Consequently it is proposed to correct this effect in dynalogs analysis to determine the MLC performance.
Evaluation of KDP Estimation Algorithm Performance in Rain Using a Known-Truth Framework
Accurate estimation of specific differential phase ( K DP ) is necessary for rain rate estimation, attenuation correction, and hydrometeor classification algorithms. There are numerous published methods to process polarimetric radar observations of propagation differential phase shift (Φ DP ) and estimate K DP , but the corresponding K DP estimate uncertainty is unquantified. This study provides guidance on how commonly used K DP estimation algorithms perform in various environments. Here, we create numerous synthetic (“true”) K DP profiles, integrate over them to obtain “smoothed” Φ DP , and then add noise typical of S-band operational weather radar measurements. Each algorithm is applied to our noisy Φ DP profiles and compared to the true K DP profile such that the errors and uncertainty are quantified. The synthetic K DP profiles are Gaussian in shape, which allows systematic variations in their magnitude and width to determine how each algorithm performs in smooth, slowly changing K DP profiles, as well as steep profiles. Results demonstrate that algorithm performance is dependent on the Φ DP field received. These results are further supported by an error analysis of each algorithm for two more complicated synthetic K DP profiles. Some K DP algorithms allow users to change various tuning parameters; a subset of these tuning parameters is tested to provide guidance on how changing these parameters impacts algorithm performance. We then provide evidence that our known-truth framework provides insight into algorithm performance in observed data through two case studies.
Dependence of Cosmological Constraints on Gray Photometric Zero-point Uncertainties of Supernova Surveys
Abstract Type Ia supernova (SN) measurements of the Hubble constant, H 0 ; cosmic mass density, Ω M ; and dark energy equation-of-state parameter, w , rely on heterogeneous SN surveys across three decades of observation. These distinct surveys may have undiagnosed, relative photometric zero-point errors. We determine the sensitivities of the SH0ES+Pantheon+ cosmological constraints to unknown gray systematics in the photometric zero-point calibration between the 19 surveys that comprise the Pantheon+ SN compendium. Varying the surveys’ gray zero-points simultaneously with cosmological parameters, we determine that the SH0ES+Pantheon+ measurement of H 0 is robust against gray intersurvey photometric miscalibration. Specifically, uncalibrated intersurvey systematics could represent a source of uncertainty no larger than 0.39 km s −1 Mpc −1 for H 0 . This modest increase in H 0 uncertainty could not account for the 7 km s −1 Mpc −1 “Hubble tension” between the SH0ES measurement of H 0 and the Planck ΛCDM-based inference of H 0 . However, the SH0ES+Pantheon+ best-fit values of Ω M and w are not robust against gray zero-point error, slipping by up to 0.16 and −0.63. Because measurements of Ω M and w depend on intrasurvey cross-band calibration, the hypothetical gray miscalibration underestimates the sensitivity of these measurements to zero-point miscalibration. Because the Pantheon+ compendium contains many surveys that share low- z Hubble flow (HF) and Cepheid-paired SNe, intersurvey photometric calibration errors do not significantly impede the joint use of SH0ES and Pantheon+ to measure H 0 to 1% accuracy. However, H 0 constraints that rely on one HF survey but numerous galactic distance calibration surveys are susceptible to intersurvey photometric miscalibration.
Quantitative 14 N NMR with Monte Carlo Uncertainty Analysis of Nitrate/Nitrite in Alkaline Nuclear Waste
While monitoring of nitrate and nitrite concentrations is important for managing corrosion in nuclear waste systems, existing analytical methods are hindered by turbidity, spectral interference, and delays from sample handling. Here, we demonstrate quantitative 14 N nuclear magnetic resonance (qNMR) spectroscopy as a direct, matrix-tolerant approach for nitrate and nitrite detection at natural abundance. Monte Carlo resampling was integrated into the workflow to quantify random error, establish precision–time tradeoffs, and separate noise-limited uncertainty from systematic bias arising from shimming, transmitter offset, or excitation pulse conditions. Quantification of nitrate and nitrite were validated in controlled alkaline matrix challenges and in 18-component Hanford-type simulants. These results establish 14 N qNMR as a practical, uncertainty-bounded tool for monitoring redox-active nitrogen species in chemically complex environments and provide a generalizable framework for quantitative analysis of quadrupolar nuclei.
STEPs-SOL, a Peptoid Force Field Parameterization to Include Solvent Effects
As peptoids (N-substituted glycines) continue to gain popularity as a class of biomimetic polymers, the importance and demand for accurate force fields in molecular simulations also grow. Building on the vacuum-optimized Systematic and Extensible Force Field for Peptoids (STEPs) force field, here we present STEPs-SOL, a novel peptoid force field parametrization that effectively incorporates solvent effects to enhance the accuracy of peptoid simulations. The development of STEPs-SOL is based on the need for precise electrostatic modeling achieved through solvent-specific partial charge optimization. Here, our systematic approach significantly improves agreement with experimental measurements, reducing the mean absolute error in cis/trans ratio predictions (ΔG c/t ) by an average of 38% across multiple peptoid residues and solvent environments. This improved parametrization addresses computational challenges associated with nonbonded energies while maintaining a workflow that relies on high-level quantum mechanical data rather than depending solely on limited experimental equilibrium properties. By evaluating the effects of conformational bias in restrained electrostatic potential (RESP) charge generation and examining their impact on peptoid conformations in various solvents, we enhance our understanding of peptoid structural dynamics while providing a more accurate modeling framework.
Toward Accurate Spin–Orbit Splittings from Relativistic Multireference Electronic Structure Theory
Most nonrelativistic electron correlation methods can be adapted to account for relativistic effects, as long as the relativistic molecular spinor integrals are available, from either a four-, two-, or one-component mean-field calculation. Furthermore, relativistic multireference correlation methods remain a relatively unexplored area, with mixed evidence regarding the improvements brought by perturbative treatments. We report, for the first time, the implementation of state-averaged four-component relativistic multireference perturbation theories to second and third order based on the driven similarity renormalization group (DSRG). With our methods, named 4c-SA-DSRG-MRPT2 and 3, we find that the dynamical correlation included on top of 4c-CASSCF references can significantly improve the spin-orbit splittings in p-block elements and potential energy surfaces when compared to 4c-CASSCF and 4c-CASPT2 results. We further show that 4c-DSRG-MRPT2 and 3 are applicable to these systems over a wide range of the flow parameter, with systematic improvement from second to third order in terms of both improved error statistics and reduced sensitivity with respect to the flow parameter.
Nonlinear optical higher-order effects in an optical lattice clock
The development of optical frequency standards with a relative uncertainty of reproducing the time and frequency units at a level as low as 10{sup −17}–10{sup −18} calls for an unprecedented accuracy in estimating the role of higher orders of optical nonlinearity, caused by the influence of the optical lattice on the frequency shift of the ‘clock transition’. This paper presents a systematic calculation of the contributions of multipole nonlinear anharmonic effects to the error of clocks based on optical lattices for alkaline-earth-like Sr, Yb, and Hg atoms. (paper)
COSMIC SHEAR RESULTS FROM THE DEEP LENS SURVEY. I. JOINT CONSTRAINTS ON {Omega} {sub M} AND {sigma}{sub 8} WITH A TWO-DIMENSIONAL ANALYSIS
We present a cosmic shear study from the Deep Lens Survey (DLS), a deep BVRz multi-band imaging survey of five 4 deg{sup 2} fields with two National Optical Astronomy Observatory (NOAO) 4 m telescopes at Kitt Peak and Cerro Tololo. For both telescopes, the change of the point-spread-function (PSF) shape across the focal plane is complicated, and the exposure-to-exposure variation of this position-dependent PSF change is significant. We overcome this challenge by modeling the PSF separately for individual exposures and CCDs with principal component analysis (PCA). We find that stacking these PSFs reproduces the final PSF pattern on the mosaic image with high fidelity, and the method successfully separates PSF-induced systematics from gravitational lensing effects. We calibrate our shears and estimate the errors, utilizing an image simulator, which generates sheared ground-based galaxy images from deep Hubble Space Telescope archival data with a realistic atmospheric turbulence model. For cosmological parameter constraints, we marginalize over shear calibration error, photometric redshift uncertainty, and the Hubble constant. We use cosmology-dependent covariances for the Markov Chain Monte Carlo analysis and find that the role of this varying covariance is critical in our parameter estimation. Our current non-tomographic analysis alone constrains the {Omega} {sub M}-{sigma}{sub 8} likelihood contour tightly, providing a joint constraint of {Omega} {sub M} = 0.262 {+-} 0.051 and {sigma}{sub 8} = 0.868 {+-} 0.071. We expect that a future DLS weak-lensing tomographic study will further tighten these constraints because explicit treatment of the redshift dependence of cosmic shear more efficiently breaks the {Omega} {sub M}-{sigma}{sub 8} degeneracy. Combining the current results with the Wilkinson Microwave Anisotropy Probe 7 year (WMAP7) likelihood data, we obtain {Omega} {sub M} = 0.278 {+-} 0.018 and {sigma}{sub 8} = 0.815 {+-} 0.020.
Accurate shear estimation with fourth-order moments
ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.
A measurement of CMB cluster lensing with SPT and DES year 1 data
Clusters of galaxies gravitationally lens the cosmic microwave background (CMB) radiation, resulting in a distinct imprint in the CMB on arcminute scales. Measurement of this effect offers a promising way to constrain the masses of galaxy clusters, particularly those at high redshift. We use CMB maps from the South Pole Telescope Sunyaev-Zel'dovich (SZ) survey to measure the CMB lensing signal around galaxy clusters identified in optical imaging from first year observations of the Dark Energy Survey. The cluster catalog used in this analysis contains 3697 members with mean redshift of $\bar{z} = 0.45$. We detect lensing of the CMB by the galaxy clusters at $8.1\sigma$ significance. Using the measured lensing signal, we constrain the amplitude of the relation between cluster mass and optical richness to roughly $17\%$ precision, finding good agreement with recent constraints obtained with galaxy lensing. The error budget is dominated by statistical noise but includes significant contributions from systematic biases due to the thermal SZ effect and cluster miscentering.
Moments of nucleon isovector structure functions in 2+1+1-flavor QCD
We present results on the isovector momentum fraction, $\langle x \rangle_{u–d}$ , helicity moment, $\langle x \rangle_{Δ u–Δd}$ , and the transversity moment, $\langle x \rangle_{δu–δd}$, of the nucleon obtained using nine ensembles of gauge configurations generated by the MILC Collaboration using 2 + 1 + 1 -flavors of dynamical highly improved staggered quarks. The correlation functions are calculated using the Wilson-Clover action, and the renormalization of the three operators is carried out nonperturbatively on the lattice in the RI'–MOM scheme. The data have been collected at lattice spacings a ≈ 0.15 , 0.12, 0.09, and 0.06 fm and $M_π$ ≈ 310 , 220, and 135 MeV, which are used to obtain the physical values using a simultaneous chiral-continuum-finite-volume fit. The final results, in the $\overline{\text{MS}}$ scheme at 2 GeV, are $\langle$ x $\rangle$ u – d = 0.173 ( 14 ) ( 07 ) , $\langle x \rangle_{Δu–Δd}$ = 0.213 ( 15 ) ( 22 ) , and $\langle x \rangle_{δu–δd}$ = 0.208 ( 19 ) ( 24 ) , where the first error is the overall analysis uncertainty and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with other recent lattice calculations and phenomenological global fit values.
Cosmology from cosmic shear with Dark Energy Survey Science Verification data
We present the first constraints on cosmology from the Dark Energy Survey (DES), using weak lensing measurements from the preliminary Science Verification (SV) data. We use 139 square degrees of SV data, which is less than 3% of the full DES survey area. Using cosmic shear 2-point measurements over three redshift bins we find σ 8 (m=0.3) 0.5 = 0:81 ± 0:06 (68% confidence), after marginalising over 7 systematics parameters and 3 other cosmological parameters. Furthermore, we examine the robustness of our results to the choice of data vector and systematics assumed, and find them to be stable. About 20% of our error bar comes from marginalising over shear and photometric redshift calibration uncertainties. The current state-of-the-art cosmic shear measurements from CFHTLenS are mildly discrepant with the cosmological constraints from Planck CMB data. Our results are consistent with both datasets. Our uncertainties are ~30% larger than those from CFHTLenS when we carry out a comparable analysis of the two datasets, which we attribute largely to the lower number density of our shear catalogue. We investigate constraints on dark energy and find that, with this small fraction of the full survey, the DES SV constraints make negligible impact on the Planck constraints. The moderate disagreement between the CFHTLenS and Planck values of σ 8 (Ω m =0.3) 0.5 is present regardless of the value of w.
Algorithm-Based Fault Tolerance for Convolutional Neural Networks
Convolutional neural networks (CNNs) are becoming more and more important for solving challenging and critical problems in many fields. CNN inference applications have been deployed in safety-critical systems, which may suffer from soft errors caused by high-energy particles, high temperature, or abnormal voltage. Of critical importance is ensuring the stability of the CNN inference process against soft errors. Traditional fault tolerance methods are not suitable for CNN inference because error-correcting code is unable to protect computational components, instruction duplication techniques incur high overhead, and existing algorithm-based fault tolerance (ABFT) techniques cannot protect all convolution implementations. In this paper, we focus on how to protect the CNN inference process against soft errors as efficiently as possible, with the following three contributions. (1) We propose several systematic ABFT schemes based on checksum techniques and analyze their fault protection ability and runtime thoroughly. Unlike traditional ABFT based on matrix-matrix multiplication, our schemes support any convolution implementations. (2) We design a novel workflow integrating all the proposed schemes to obtain a high detection/correction ability with limited total runtime overhead. (3) We perform our evaluation using ImageNet with well-known CNN models including AlexNet, VGG-19, ResNet-18, and YOLOv2. Here, experimental results demonstrate that our implementation can handle soft errors with very limited runtime overhead (4%~8% in both error-free and error-injected situations).
Metals and Quantum Materials with Spin-orbit Interactions by Quantum Monte Carlo methods
The key goals of this project were as follows: 1) Analysis and benchmarks of electron correlation effects recovered in the fixed-node approximation that is inherent to quantum Monte Carlo (QMC) method as applied to metallic states; 2) development of new algorithms for electron spin-degrees of freedom to be treated as explicit quantum variables; 3) designing electronic structure QMC algorithm for efficient evaluation of spin-orbit effects in systems with heavy atoms; 4) adapting the algorithm to complex wave functions and developing corresponding fixed-phase approximation; 5) design and testing of algorithm for valence-only non-local spin-orbit operators; 6) analysis of fixed-node vs fixed-phase errors and their comparisons. The key accomplishments: i) We carried out a systematic study of Li systems by the fixed-node diffusion Monte Carlo method. This involved Li atom, molecule, cluster and solid calculated by the full range of QMC methods including fixed-node QMC.
Thermonuclear Reaction Rates and Primordial Nucleosynthesis
Assuming the best numerical value for the cosmic baryonic density and the existence of three neutrino flavors, standard Big Bang nucleosynthesis is a parameter-free model. It is important to assess if the observed primordial abundances can be reproduced by simulations. Numerous studies have shown that the simulations overpredict the primordial 7 Li abundance by a factor of ≈3 compared to the observations. The discrepancy may be caused by unknown systematics in 7 Li observations, poorly understood depletion of lithium in stars, errors in thermonuclear rates that take part in the lithium and beryllium synthesis, or physics beyond the standard model. Here, we focus on the likelihood of a nuclear physics solution. The status of the key nuclear reaction rates is summarized. Big Bang nucleosynthesis simulations are performed with the most recent reaction rates, and the uncertainties of the predicted abundances are established using a Monte Carlo technique. Correlations between abundances and reaction rates are investigated based on the metric of mutual information. The rates of four reactions impact the primordial 7 Li abundance: 3 He(α,γ) 7 Be, d(p,γ) 3 He, 7 Be(d,p)2α, and 7 Be(n,p) 7 Li. We employ a genetic algorithm to search for simultaneous rate changes in these four reactions that may account for all observed primordial abundances. Furthermore, when the search is performed for reaction rate ranges that are much wider than recently reported uncertainties, no acceptable solutions are found. Based on the currently available evidence, we conclude that it is highly unlikely for the cosmological lithium problem to have a nuclear physics solution.