Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Statistical 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.

At least 1,135 records · Page 63

Power and Limitations of Linear Programming Decoder for Quantum LDPC Codes

Decoding quantum error-correcting codes is a key challenge in enabling fault-tolerant quantum computation. In the classical setting, linear programming (LP) decoders offer provable performance guarantees and can leverage fast practical optimization algorithms. Although LP decoders have been proposed for quantum codes, their performance and limitations remain relatively underexplored. In this work, we uncover a key limitation of LP decoding for quantum low-density parity-check (LDPC) codes: certain constant-weight error patterns lead to ambiguous fractional solutions that cannot be resolved through independent rounding. To address this issue, we incorporate a post-processing technique known as ordered statistics decoding (OSD), which significantly enhances LP decoding performance in practice. Our results show that LP decoding, when augmented with OSD, can outperform belief propagation with the same post-processing for intermediate code sizes of up to hundreds of qubits. These findings suggest that LP-based decoders, equipped with effective post-processing, offer a promising approach for decoding near-term quantum LDPC codes.

Gu, Shouzhen [Yale U.]↗

DESI DR2 results. I. Baryon acoustic oscillations from the Lyman alpha forest

We present the baryon acoustic oscillation (BAO) measurements with the Lyman-𝛼 (Ly⁢𝛼) forest from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our BAO measurements include both the autocorrelation of the Ly⁢𝛼 forest absorption observed in the spectra of high-redshift quasars and the cross-correlation of the absorption with the quasar positions. The total sample size is approximately a factor of 2 larger than the DR1 dataset, with forest measurements in over 820,000 quasar spectra and the positions of over 1.2 million quasars. We describe several significant improvements to our analysis in this paper, and two supporting papers describe improvements to the synthetic datasets that we use for validation and how we identify damped Ly⁢𝛼 absorbers. Our main result is that we have measured the BAO scale with a statistical precision of 1.1% along and 1.3% transverse to the line of sight, for a combined precision of 0.65% on the isotropic BAO scale at 𝑧 eff =2.33. This excellent precision, combined with recent theoretical studies of the BAO shift due to nonlinear growth, motivated us to include a systematic error term in Ly⁢𝛼 BAO analysis for the first time. We measure the ratios 𝐷 𝐻 ⁡(𝑧 eff )/𝑟 𝑑 = 8.632 ± 0.098 ± 0.026 and 𝐷 𝑀 ⁡(𝑧 eff )/𝑟 𝑑 = 38.99 ± 0.52 ± 0.12, where 𝐷 𝐻 = 𝑐/𝐻⁡(𝑧) is the Hubble distance, 𝐷 𝑀 is the transverse comoving distance, 𝑟 𝑑 is the sound horizon at the drag epoch, and we quote both the statistical and the theoretical systematic uncertainty. The companion paper presents the BAO measurements at lower redshifts from the same dataset and the cosmological interpretation.

baryon acoustic oscillations↗

The coyote universe extended: Precision emulation of the matter power spectrum

Modern sky surveys are returning precision measurements of cosmological statistics such as weak lensing shear correlations, the distribution of galaxies, and cluster abundance. To fully exploit these observations, theorists must provide predictions that are at least as accurate as the measurements, as well as robust estimates of systematic errors that are inherent to the modeling process. In the nonlinear regime of structure formation, this challenge can only be overcome by developing a large-scale, multi-physics simulation capability covering a range of cosmological models and astrophysical processes. As a first step to achieving this goal, we have recently developed a prediction scheme for the matter power spectrum (a so-called emulator), accurate at the 1% level out to k ∼ 1 Mpc{sup –1} and z = 1 for wCDM cosmologies based on a set of high-accuracy N-body simulations. It is highly desirable to increase the range in both redshift and wavenumber and to extend the reach in cosmological parameter space. To make progress in this direction, while minimizing computational cost, we present a strategy that maximally reuses the original simulations. We demonstrate improvement over the original spatial dynamic range by an order of magnitude, reaching k ∼ 10 h Mpc{sup –1}, a four-fold increase in redshift coverage, to z = 4, and now include the Hubble parameter as a new independent variable. To further the range in k and z, a new set of nested simulations run at modest cost is added to the original set. The extension in h is performed by including perturbation theory results within a multi-scale procedure for building the emulator. This economical methodology still gives excellent error control, ∼5% near the edges of the domain of applicability of the emulator. A public domain code for the new emulator is released as part of the work presented in this paper.

79 ASTRONOMY AND ASTROPHYSICS↗

DESI DR2 Results I: Baryon Acoustic Oscillations from the Lyman Alpha Forest

We present the Baryon Acoustic Oscillation (BAO) measurements with the Lyman-alpha (LyA) forest from the second data release (DR2) of the Dark Energy Spectroscopic Instrument (DESI) survey. Our BAO measurements include both the auto-correlation of the LyA forest absorption observed in the spectra of high-redshift quasars and the cross-correlation of the absorption with the quasar positions. The total sample size is approximately a factor of two larger than the DR1 dataset, with forest measurements in over 820,000 quasar spectra and the positions of over 1.2 million quasars. We describe several significant improvements to our analysis in this paper, and two supporting papers describe improvements to the synthetic datasets that we use for validation and how we identify damped LyA absorbers. Our main result is that we have measured the BAO scale with a statistical precision of 1.1% along and 1.3% transverse to the line of sight, for a combined precision of 0.65% on the isotropic BAO scale at $z_{eff} = 2.33$. This excellent precision, combined with recent theoretical studies of the BAO shift due to nonlinear growth, motivated us to include a systematic error term in LyA BAO analysis for the first time. We measure the ratios $D_H(z_{eff})/r_d = 8.632 \pm 0.098 \pm 0.026$ and $D_M(z_{eff})/r_d = 38.99 \pm 0.52 \pm 0.12$, where $D_H = c/H(z)$ is the Hubble distance, $D_M$ is the transverse comoving distance, $r_d$ is the sound horizon at the drag epoch, and we quote both the statistical and the theoretical systematic uncertainty. The companion paper presents the BAO measurements at lower redshifts from the same dataset and the cosmological interpretation.

79 ASTRONOMY AND ASTROPHYSICS↗

THE LARGE-SCALE THREE-POINT CORRELATION FUNCTION OF SLOAN DIGITAL SKY SURVEY LUMINOUS RED GALAXIES

We present new measurements of the redshift-space three-point correlation function (3PCF) of Luminous Red Galaxies (LRGs) from the Sloan Digital Sky Survey (SDSS). Using the largest data set to date, the Data Release 7 LRGs, and an improved binning scheme compared to previous measurements, we measure the LRG 3PCF on large scales up to {approx}90 h{sup -1} Mpc, from the mildly nonlinear to quasi-linear regimes. Comparing the LRG correlations to the dark matter two- and three-point correlation functions, obtained from N-body simulations we infer linear and nonlinear bias parameters. As expected, LRGs are highly biased tracers of large-scale structure, with a linear bias b{sub 1} {approx} 2; the LRGs also have a large positive nonlinear bias parameter, in agreement with predictions of galaxy population models. The use of the 3PCF to estimate biasing helps to also make estimates of the cosmological parameter {sigma}{sub 8}, as well as to infer best-fit parameters of the halo occupation distribution parameters for LRGs. We also use a large suite of public mock catalogs to characterize the error covariance matrix for the 3PCF and compare the variance among simulation results with jackknife error estimates.

79 ASTRONOMY AND ASTROPHYSICS↗

Addressing Interdependency in a Multimodel Ensemble by Interpolation of Model Properties

The diverse set of Earth system models used to conduct the CMIP5 ensemble can partly sample the uncertainties in future climate projections. However, combining those projections is complicated by the fact that models developed by different groups share ideas and code and therefore biases. The authors propose a method for combining model results into single or multivariate distributions that are more robust to the inclusion of models with a large degree of interdependency. This study uses a multivariate metric of present-day climatology to assess both model performance and similarity in two recent model intercomparisons, CMIP3 and CMIP5. Model characteristics can be interpolated and then resampled in a space defined by independent climate properties. A form of weighting can be applied by sampling more densely in the region of the space close to the projected observations, thus taking into account both model performance and interdependence. Furthermore, the choice of the sampling distribution’s parameters is a subjective decision that should reflect the researcher’s prior assumptions as to the acceptability of different model errors.

58 GEOSCIENCES↗

SAMBA: Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos

A new arbitrary Polynomial Chaos (aPC) method is presented for moderately high-dimensional problems characterised by limited input data availability. The proposed methodology improves the algorithm of aPC and extends the method, that was previously only introduced as tensor product expansion, to moderately high-dimensional stochastic problems. The fundamental idea of aPC is to use the statistical moments of the input random variables to develop the polynomial chaos expansion. This approach provides the possibility to propagate continuous or discrete probability density functions and also histograms (data sets) as long as their moments exist, are finite and the determinant of the moment matrix is strictly positive. For cases with limited data availability, this approach avoids bias and fitting errors caused by wrong assumptions. In this work, an alternative way to calculate the aPC is suggested, which provides the optimal polynomials, Gaussian quadrature collocation points and weights from the moments using only a handful of matrix operations on the Hankel matrix of moments. It can therefore be implemented without requiring prior knowledge about statistical data analysis or a detailed understanding of the mathematics of polynomial chaos expansions. The extension to more input variables suggested in this work, is an anisotropic and adaptive version of Smolyak's algorithm that is solely based on the moments of the input probability distributions. It is referred to as SAMBA (PC), which is short for Sparse Approximation of Moment-Based Arbitrary Polynomial Chaos. It is illustrated that for moderately high-dimensional problems (up to 20 different input variables or histograms) SAMBA can significantly simplify the calculation of sparse Gaussian quadrature rules. SAMBA's efficiency for multivariate functions with regard to data availability is further demonstrated by analysing higher order convergence and accuracy for a set of nonlinear test functions with 2, 5 and 10 different input distributions or histograms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Rapid neutron and gamma-ray source localization using machine learning

Rapid localization of radiation sources is critical for applications including nuclear emergency response, safeguards, and security. However, conventional imaging systems such as neutron scatter cameras and Compton cameras depend on rare coincidence events, which often result in long acquisition times. In this work, we address the challenge of rapid source localization by developing a machine learning approach to predict the direction of a single radiation source using only count rates from an array of neutron and gamma-ray detectors. The proposed model is a fully connected neural network (FCNN) trained using Monte Carlo simulation data from a 252 Cf source. The model hyperparameters are optimized with a small set of routine 252 Cf measurements. We benchmarked the performance of the trained and optimized machine learning model using additional 252 Cf , 137 Cs , and PuBe measurements under laboratory conditions with varying source-detector configurations. For these measurements, the machine learning model achieved a mean localization error smaller than 30° with 3 x 10 3 system counts, corresponding to 8 s measurement time for the imaging system used in this work. In this low-statistics regime, the method outperformed traditional scatter-based imaging by more than 75% in localization accuracy for the evaluated measurement configurations. These results demonstrate that a machine learning-based approach can significantly reduce the time required for accurate single-source localization, providing a robust and computationally efficient alternative to traditional imaging systems in time-critical nuclear security and emergency response scenarios.

Gamma-ray imaging↗

Detector-Response Correction of Two-Dimensional γ -Ray Spectra from Neutron Capture

The neutron-capture reaction produces a large variety of γ-ray cascades with different γ-ray multiplicities. A measured spectral distribution of these cascades for each γ-ray multiplicity is of importance to applications and studies of γ-ray statistical properties. The DANCE array, a 4π ball of 160 BaF2 detectors, is an ideal tool for measurement of neutron-capture γ-rays. The high granularity of DANCE enables measurements of high-multiplicity γ-ray cascades. The measured two-dimensional spectra (γ-ray energy, γ-ray multiplicity) have to be corrected for the DANCE detector response in order to compare them with predictions of the statistical model or use them in applications. The detector-response correction problem becomes more difficult for a 4π detection system than for a single detector. A trial and error approach and an iterative decomposition of γ-ray multiplets, have been successfully applied to the detector-response correction. Applications of the decomposition methods are discussed for two-dimensional γ-ray spectra measured at DANCE from γ-ray sources and from the 10B(n, γ) and 113Cd(n, γ) reactions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

High-Fidelity Velocity and Concentration Measurements of Turbulent Buoyant Jets

Accurate models of turbulent buoyant flows are essential for the design of the cooling circuit of nuclear reactors and passive safety systems. However, available models fail to fully capture the physics of turbulent mixing when buoyancy becomes predominant with respect to momentum. Therefore, high-fidelity experiments of well-controlled fundamental flows are needed to develop and validate more accurate models. We analyze experiments of positive and negative turbulent buoyant jets, both in uniform and stratified environments, with the aim of understanding the thermal hydraulics of turbulent mixing with variable density and providing high-fidelity data for the development and validation of turbulence models. Non-intrusive, simultaneous particle image velocimetry and laser-induced fluorescence measurements were carried out to acquire instantaneous velocity and concentration fields on a vertical section parallel to the axis of a jet in the self-similar region. The refractive index matching method was applied to measure high-resolution buoyant jets with up to 8.6% density difference. These data are free of the typical errors that characterize optical measurements of buoyancy-driven flows (e.g. natural and mixed convection) where the refractive index of the fluid is inhomogeneous throughout the measurement domain. Turbulent statistics and entrainment of buoyant jets in uniform and stratified environments are presented. These data are compared with non-buoyant jets in a uniform environment, as a reference to investigate the effects of buoyancy and stratification on turbulent mixing. The results will be used for the assessment of current turbulence models and as a basis for the development of a new model that captures turbulent mixing.

laser induced fluorescence↗

Composite-dimensional topological codes with boundaries and defects

We introduce new algorithms and provide example constructions of stabilizer models for the gapped boundaries, domain walls, and 0D defects of Abelian composite-dimensional twisted quantum doubles. Using the physically intuitive concept of condensation, our algorithm explicitly describes how to construct the boundary and domain-wall stabilizers starting from the bulk model. This extends the utility of Pauli stabilizer models in describing nontranslationally invariant topological orders with gapped boundaries. To highlight this utility, we provide a series of examples, including a new family of quantum error-correcting codes where the double of ℤ4 is coupled to instances of the double semion (DS) phase. We discuss the codes' utility in the burgeoning area of quantum error correction with an emphasis on the interplay between deconfined anyons, logical operators, error rates, and decoding. We also augment our construction, built using algorithmic tools to describe the properties of explicit stabilizer layouts at the microscopic lattice level, with dimensional counting arguments and macroscopic-level constructions building on pants decompositions. The latter outlines how such codes' representation and design can be automated. Our results are validated by a series of error-correcting threshold calculations comparing our codes' performance with that of standard surface codes. To do so, we introduce a composite-dimensional belief-propagation decoder with ordered statistics that utilizes combination sweeps. Going beyond our worked-out examples, we expect our explicit step-by-step algorithms to pave the path for higher-dimensional codes to be discovered and implemented in near-future architectures that take advantage of various hardware platforms.

Mousa, Mohamad [Purdue University]↗

Geoneutrinos and Hydridic Earth (or primordially hydrogen-rich planet)

Geoneutrino is a new channel of information about geochemical composition of the Earth. We analyzed here the following problem. What statistics do we need to distinguish between predictions of Bulk Silicate Earth model and Hydridic Earth model for Th/U signal ratio? We obtained the simple formula for estimation of error of Th/U signal ratio. Our calculations show that we need more than 22 kt · year exposition for Gran-Sasso underground laboratory and Sudbury Neutrino Observatory. We need more than 27 kt · year exposition for Kamioka site in the case of stopping of all Japanese nuclear power plants.

54 ENVIRONMENTAL SCIENCES↗

Application of Linear Additive Conditions for Near-Infrared Diffuse Reflectance Absorption Spectroscopy

Determining the homogeneity of material mixing in real time during product processing is critical for quality control. According to the Kubelka–Munk (K-M) function of diffuse reflectance absorption spectrum, absorbance (A) is approximately linear with the content of the components when the sample scattering coefficient (S) is in a certain range. The S is determined by the particle size of powder samples. Therefore, this study determined particle size ranges that satisfy linear additivity in near-infrared diffuse reflectance spectroscopy (NIRDRS). Thus, the proposed NIRDRS analysis technique can be used to determine the homogeneity of material mixes or analyze the percentages of the components in the mixture. In this study, vitamin B3 and vitamin C were used for preparing mixed samples with varying percentages. The experimental results revealed that linear additivity is satisfied when the powder particle size is in the range of less than 280, 280–450, and 450–900 μm. When the confidence level is 0.01, the actual mixed spectra are not significantly different from the “simulated mixed spectra” constructed by linear addition, with their relative deviations less than 1.08%. The absolute errors of the actual and analytic percentages were within 2.98% for each component in the mixtures. The above conclusions also hold for sorghum, which has a complex material composition. Statistical models cannot analyze the percentages of components in the mixture. In contrast, linear addition and direct calibration approach avoids the use of a large number of samples for statistical modeling and analyze the percentages of mixed samples. Meanwhile, it can be used to discriminate and analyze the material mixing uniformity by building a mechanistic model.

Feng, Zhiyue↗

Offshore wind energy forecasting sensitivity to sea surface temperature input in the Mid-Atlantic

Abstract. As offshore wind farm development expands, accurate wind resource forecasting over the ocean is needed. One important yet relatively unexplored aspect of offshore wind resource assessment is the role of sea surface temperature (SST). Models are generally forced with reanalysis data sets, which employ daily SST products. Compared with observations, significant variations in SSTs that occur on finer timescales are often not captured. Consequently, shorter-lived events such as sea breezes and low-level jets (among others), which are influenced by SSTs, may not be correctly represented in model results. The use of hourly SST products may improve the forecasting of these events. In this study, we examine the sensitivity of model output from the Weather Research and Forecasting model (WRF) 4.2.1 to different SST products. We first evaluate three different data sets: the Multiscale Ultrahigh Resolution (MUR25) SST analysis, a daily, 0.25∘ × 0.25∘ resolution product; the Operational Sea Surface Temperature and Ice Analysis (OSTIA), a daily, 0.054∘ × 0.054∘ resolution product; and SSTs from the Geostationary Operational Environmental Satellite 16 (GOES-16), an hourly, 0.02∘ × 0.02∘ resolution product. GOES-16 is not processed at the same level as OSTIA and MUR25; therefore, the product requires gap-filling using an interpolation method to create a complete map with no missing data points. OSTIA and GOES-16 SSTs validate markedly better against buoy observations than MUR25, so these two products are selected for use with model simulations, while MUR25 is at this point removed from consideration. We run the model for June and July of 2020 and find that for this time period, in the Mid-Atlantic, although OSTIA SSTs overall validate better against in situ observations taken via a buoy array in the area, the two products result in comparable hub-height (140 m) wind characterization performance on monthly timescales. Additionally, during hours-long flagged events (< 30 h each) that show statistically significant wind speed deviations between the two simulations, both simulations once again demonstrate similar validation performance (differences in bias, earth mover's distance, correlation, and root mean square error on the order of 10−1 or less), with GOES-16 winds validating nominally better than OSTIA winds. With a more refined GOES-16 product, which has been not only gap-filled but also assimilated with in situ SST measurements in the region, it is likely that hub-height winds characterized by GOES-16-informed simulations would definitively validate better than those informed by OSTIA SSTs.

17 WIND ENERGY↗

p-Wave Cold Collisions in an Optical Lattice Clock

We study ultracold collisions in fermionic ytterbium by precisely measuring the energy shifts they impart on the atoms' internal clock states. Exploiting Fermi statistics, we uncover p-wave collisions, in both weakly and strongly interacting regimes. With the higher density afforded by two-dimensional lattice confinement, we demonstrate that strong interactions can lead to a novel suppression of this collision shift. In addition to reducing the systematic errors of lattice clocks, this work has application to quantum information and quantum simulation with alkaline-earth atoms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Understanding boundary effects in quantum state tomography – One qubit case

For classical and quantum estimation with finite data sets, the estimation error can deviate significantly from its asymptotic (large data set) behavior. In quantum state tomography, a major reason for this is the existence of a boundary in the parameter space imposed by constraints, such as the positive semidefiniteness of density matrices. Intuitively, we should be able to reduce the estimation error by using our knowledge of these constraints. This intuition is correct for maximumlikelihood estimators, but the size of the reduction has not been evaluated quantitatively. In this proceeding, we evaluate the improvement in one qubit state tomography by using mathematical tools in classical statistical estimation theory. In particular, we show that the effect of the reduction decreases exponentially with respect to the number of data sets when the true state is mixed, and it remains at arbitrarily large data set when the true state is pure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Advances in statistical methods for cancer surveillance research: an age-period-cohort perspective

Background: Analysis of Lexis diagrams (population-based cancer incidence and mortality rates indexed by age group and calendar period) requires specialized statistical methods. However, existing methods have limitations that can now be overcome using new approaches. Methods: We assembled a “toolbox” of novel methods to identify trends and patterns by age group, calendar period, and birth cohort. We evaluated operating characteristics across 152 cancer incidence Lexis diagrams compiled from United States (US) Surveillance, Epidemiology and End Results Program data for 21 leading cancers in men and women in four race and ethnicity groups (the “cancer incidence panel”). Results: Nonparametric singular values adaptive kernel filtration (SIFT) decreased the estimated root mean squared error by 90% across the cancer incidence panel. A novel method for semi-parametric age-period-cohort analysis (SAGE) provided optimally smoothed estimates of age-period-cohort (APC) estimable functions and stabilized estimates of lack-of-fit (LOF). SAGE identified statistically significant birth cohort effects across the entire cancer panel; LOF had little impact. As illustrated for colon cancer, newly developed methods for comparative age-period-cohort analysis can elucidate cancer heterogeneity that would otherwise be difficult or impossible to discern using standard methods. Conclusions: Cancer surveillance researchers can now identify fine-scale temporal signals with unprecedented accuracy and elucidate cancer heterogeneity with unprecedented specificity. Birth cohort effects are ubiquitous modulators of cancer incidence in the US. The novel methods described here can advance cancer surveillance research.

60 APPLIED LIFE SCIENCES↗

BLENDING OF CEPHEIDS IN M33

A precise and accurate determination of the Hubble constant based on Cepheid variables requires proper characterization of many sources of systematic error. One of these is stellar blending, which biases the measured fluxes of Cepheids and the resulting distance estimates. We study the blending of 149 Cepheid variables in M33 by matching archival Hubble Space Telescope data with images obtained at the Wisconsin-Indiana-Yale-NOAO (WIYN) 3.5 m telescope, which differ by a factor of 10 in angular resolution. We find that 55% {+-} 4% of the Cepheids have no detectable nearby companions that could bias the WIYN V-band photometry, while the fraction of Cepheids affected below the 10% level is 73% {+-} 4%. The corresponding values for the I band are 60% {+-} 4% and 72% {+-} 4%, respectively. We find no statistically significant difference in blending statistics as a function of period or surface brightness. Additionally, we report all the detected companions within 2'' of the Cepheids (equivalent to 9 pc at the distance of M33) which may be used to derive empirical blending corrections for Cepheids at larger distances.

79 ASTRONOMY AND ASTROPHYSICS↗