Search NASA⌕ Search

SEARCH · Search NASA

Results for “Function Estimation”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Random insights into the complexity of two-dimensional tensor network calculations

Projected entangled pair states (PEPS) offer memory-efficient representations of some quantum many-body states that obey an entanglement area law and are the basis for classical simulations of ground states in two-dimensional (2d) condensed matter systems. However, rigorous results show that exactly computing observables from a 2d PEPS state is generically a computationally hard problem. Yet approximation schemes for computing properties of 2d PEPS are regularly used, and empirically seen to succeed, for a large subclass of (“not too entangled”) condensed matter ground states. Adopting the philosophy of random matrix theory, in this work, we analyze the complexity of approximately contracting a 2d random PEPS by exploiting an analytic mapping to an effective replicated statistical mechanics model that permits a controlled analysis at a large bond dimension. Through this statistical-mechanics lens, we argue that (i) although approximately sampling wave-function amplitudes of random PEPS faces a computational-complexity phase transition above a critical bond dimension, and (ii) one can generically efficiently estimate the norm and correlation functions for any finite bond dimension. Furthermore, these results are supported numerically for various bond-dimension regimes. It is an important open question whether the above results for random PEPS apply more generally also to PEPS representing physically relevant ground states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Computational Optimization of 133m Xe Production via Neutron Irradiation in a TRIGA Reactor

Here, the Comprehensive Nuclear-Test-Ban Treaty bans all nuclear tests worldwide. As part of treaty compliance, the concentration of radioactive nuclides in the atmosphere is monitored to detect nuclear weapons tests. Radioactive noble gas fission products, specifically radioxenon, can vent into the atmosphere after a nuclear weapons test, even if the test is well contained underground or underwater. Radioxenon thus serves as a signal for nuclear weapons tests. All atmospheric monitoring systems require samples of radioxenon isotopes for detector calibration, quality control, and certification. Here, we present a novel, improved method for creating samples of 133m Xe via neutron irradiation of 132 Xe in the Washington State University TRIGA reactor. 132 Xe neutron absorption results in either 133 Xe or 133m Xe—thermal neutron absorption results in 133m Xe 12% of the time, while fast neutron absorption (above ~1 MeV) results in 133m Xe ~50% of the time. To optimize the production of 133m Xe via neutron absorption in 132 Xe in the thermal TRIGA reactor, spectral tuning using an irradiation chamber is required to maximize the fraction of fast neutrons being absorbed and minimize the number of thermal neutrons interacting with the 132 Xe. We used MCNP to tally 132 Xe absorptions with the isotopic tally function, flux tallies and neutron attenuation to estimate the number of neutrons reaching the 132 Xe through the irradiation chamber, and the adjoint importance function to improve the source strength estimate. Additionally, we performed a heat transfer analysis for safety considerations. It was determined that the use of a 96% enriched 10 B boron carbide chamber, placed next to the fuel elements in reactor position D8, increases the 133m Xe/ 133 Xe activity ratio from a baseline value of 0.3 to 1.0, a 233% increase. Additionally, it was determined that the alpha heating produced in the boron does not become an unmanageable problem in the Washington State University reactor.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2022

Update of the Nuclear Criticality Slide Rule: Review of the Estimation of the Number of Fissions

IRSN (France), LLNL (USA) and ORNL (USA) began a long-term collaboration effort in 2015 to update the nuclear criticality Slide Rule (https://ncsp.llnl.gov/analytical-methods/criticality-sliderule) for the emergency response to a nuclear criticality accident. The Slide Rule permits the estimation of neutron and gamma dose rates and integrated doses based upon estimated fission yields, as a function of distance from the fission source, and time after criticality accidents for different critical systems. This paper presents the review of the section “estimation of the number of fissions”, named “fission yield” in the Slide Rule document. This estimation is important because the other information provided by the Slide Rule is directly proportional to the number of fissions. The previous models, based on Hansen and Barbry formulae, will be presented with their associated assumptions and limitations. Then, new proposals will be presented with a comparison to past criticality accidents and experimental data. The new formulae, based on the heat energy equation, cover many types of criticality accidents (solution, dry or low-moderated powder, rods/assemblies in water and dry metal) and require limited information in an emergency situation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Inverse regression for ridge recovery: a data-driven approach for parameter reduction in computer experiments

Parameter reduction can enable otherwise infeasible design and uncertainty studies with modern computational science models that contain several input parameters. In statistical regression, techniques for sufficient dimension reduction (SDR) use data to reduce the predictor dimension of a regression problem. A computational scientist hoping to use SDR for parameter reduction encounters a problem: a computer prediction is best represented by a deterministic function of the inputs, so data comprised of computer simulation queries fail to satisfy the SDR assumptions. To address this problem, here we interpret SDR methods sliced inverse regression (SIR) and sliced average variance estimation (SAVE) as estimating the directions of a ridge function, which is a composition of a low-dimensional linear transformation with a nonlinear function. Within this interpretation, SIR and SAVE estimate matrices of integrals whose column spaces are contained in the ridge directions’ span; we analyze and numerically verify convergence of these column spaces as the number of computer model queries increases. Moreover, we show example functions that are not ridge functions but whose inverse conditional moment matrices are low-rank. Consequently, the computational scientist should beware when using SIR and SAVE for parameter reduction, since SIR and SAVE may mistakenly suggest that truly important directions are unimportant.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Maximum Likelihood Estimation: Some Basics

The maximum likelihood estimation is a general estimation procedure. It is often compared to estimation procedures like the ordinary least squares regression or generalized method of moments, to name a few. We discuss some basics about the maximum likelihood estimation, its advantages and disadvantages, and provide an example application to a gamma distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Maximum Likelihood Estimation: Some Basics

The maximum likelihood estimation is a general estimation procedure. It is often compared to estimation procedures like the ordinary least squares regression or generalized method of moments, to name a few. We discuss some basics about the maximum likelihood estimation, its advantages and disadvantages, and provide an example application to a gamma distribution function.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Analysis of the Behavior of Individual Emission Sites on the Surface of a Multi-Tip Field Cathode

A method of analysis of the surface of multi-tip field cathodes using data obtained with the aid of a computer-controlled field emission projector and a system of fast high-voltage scanning is described. An investigation of a polymer/carbon nanotube composite field electron emitter by this method was used to determine the distribution of individual emission sites over the surface of a multi-tip field cathode and to estimate variation of the work function of separate emission sites depending on the applied voltage. The contributions of separate groups of emission sites to a change of the shape of current–voltage characteristic of a multi-tip cathode have been numerically estimated.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Visualizing Orbital Content of Electronic Bands in Anisotropic 2D Semiconducting ReSe 2

Many properties of layered materials change as they are thinned from their bulk forms down to single layers, with examples including indirect-to-direct band gap transition in 2H semiconducting transition metal dichalcogenides as well as thickness-dependent changes in the valence band structure in post-transition-metal monochalcogenides and black phosphorus. Here, we use angle-resolved photoemission spectroscopy to study the electronic band structure of monolayer ReSe 2 , a semiconductor with a distorted 1T structure and in-plane anisotropy. By changing the polarization of incoming photons, we demonstrate that for ReSe 2 , in contrast to the 2H materials, the out-of-plane transition metal dz2 and chalcogen pz orbitals do not contribute significantly to the top of the valence band, which explains the reported weak changes in the electronic structure of this compound as a function of layer number. We estimate a band gap of 1.7 eV in pristine ReSe 2 using scanning tunneling spectroscopy and explore the implications on the gap following surface doping with potassium. Additionally, a lower bound of 1.4 eV is estimated for the gap in the fully doped case, suggesting that doping-dependent many-body effects significantly affect the electronic properties of ReSe 2 . Our results, supported by density functional theory calculations, provide insight into the mechanisms behind polarization-dependent optical properties of rhenium dichalcogenides and highlight their place among two-dimensional crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

How accurate are static polarizability predictions from density functional theory? An assessment over 132 species at equilibrium geometry

Static polarizabilities are the first response of the electron density to electric fields, and are therefore important for predicting intermolecular and molecule-field interactions. They also offer a global measure of the accuracy of the treatment of excited states by density functionals in a formally exact manner. We have developed a database of benchmark static polarizabilities for 132 small species at equilibrium geometry, using coupled cluster theory through triple excitations (extrapolated to the complete basis set limit), for the purpose of developing and assessing density functionals. The performance of 60 popular and recent functionals are also assessed, which indicates that double hybrid functionals perform the best, having RMS relative errors in the range of 2.5-3.8%. Many hybrid functionals also give quite reasonable estimates with 4-5% RMS relative error. A few meta-GGAs like mBEEF and MVS yield performance comparable to hybrids, indicating potential for improved excited state predictions relative to typical local functionals. Some recent functionals however are found to be prone to catastrophic failure (possibly as a consequence of overparameterization), indicating a need for caution in applying these.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

The two-point correlation function covariance with fewer mocks

We present FITCOV an approach for accurate estimation of the covariance of two-point correlation functions that requires fewer mocks than the standard mock-based covariance. This can be achieved by dividing a set of mocks into jackknife regions and fitting the correction term first introduced in Mohammad & Percival (2022), such that the mean of the jackknife covariances corresponds to the one from the mocks. This extends the model beyond the shot-noise limited regime, allowing it to be used for denser samples of galaxies. We test the performance of our fitted jackknife approach, both in terms of accuracy and precision, using lognormal mocks with varying densities and approximate EZmocks mimicking the Dark Energy Spectroscopic Instrument LRG and ELG samples in the redshift range of z = [0.8, 1.1]. We find that the Mohammad–Percival correction produces a bias in the two-point correlation function covariance matrix that grows with number density and that our fitted jackknife approach does not. We also study the effect of the covariance on the uncertainty of cosmological parameters by performing a full-shape analysis. We demonstrate that our fitted jackknife approach based on 25 mocks can recover unbiased and as precise cosmological parameters as the ones obtained from a covariance matrix based on 1000 or 1500 mocks, while the Mohammad–Percival correction produces uncertainties that are twice as large.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Thin-Film Thermal Conductivity Measurements Using Superconducting Nanowires

We present a simple experimental scheme for estimating the cryogenic thermal transport properties of thin films using superconducting nanowires. In a parallel array of nanowires, the heat from one nanowire in the normal state changes the local temperature around adjacent nanowires, reducing their switching current. Calibration of this change in switching current as a function of bath temperature provides an estimate of the temperature as a function of displacement from the heater. This provides a method of determining the contribution of substrate heat transport to the cooling time of superconducting nanowire single-photon detectors. Understanding this process is necessary for successful electrothermal modeling of superconducting nanowire systems.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Extraction of the specific shear viscosity of quark-gluon plasma from two-particle transverse momentum correlations

The specific shear viscosity, η/s, of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions at RHIC and LHC is estimated based on the progressive longitudinal broadening of transverse momentum two-particle correlators, G 2 , reported as a function of collision centrality by the STAR and ALICE experiments. Estimates are computed as a function of collision centrality using the Gavin ansatz which relates the G 2 longitudinal broadening to the specific shear viscosity. Freeze out times required for the use of the ansatz are computed using a linear fit of freeze out times reported as a function of the cubic root of the charged particle pseudorapidity density (dN ch /dη) 1/3 . Estimates of η/s based on ALICE data exhibit little to no dependence on collision centrality at LHC energy, while estimates obtained from STAR data hint that η/s might be a function of collision centrality at top RHIC energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

High-dimensional multivariate autoregressive model estimation of human electrophysiological data using fMRI priors

Multivariate autoregressive (MVAR) model estimation enables assessment of causal interactions in brain networks. However, accurately estimating MVAR models for high-dimensional electrophysiological recordings is challenging due to the extensive data requirements. Hence, the applicability of MVAR models for study of brain behavior over hundreds of recording sites has been very limited. Prior work has focused on different strategies for selecting a subset of important MVAR coefficients in the model to reduce the data requirements of conventional least-squares estimation algorithms. Here we propose incorporating prior information, such as resting state functional connectivity derived from functional magnetic resonance imaging, into MVAR model estimation using a weighted group least absolute shrinkage and selection operator (LASSO) regularization strategy. The proposed approach is shown to reduce data requirements by a factor of two relative to the recently proposed group LASSO method of Endemann et al (Neuroimage 254:119057, 2022) while resulting in models that are both more parsimonious and more accurate. The effectiveness of the method is demonstrated using simulation studies of physiologically realistic MVAR models derived from intracranial electroencephalography (iEEG) data. The robustness of the approach to deviations between the conditions under which the prior information and iEEG data is obtained is illustrated using models from data collected in different sleep stages. This approach allows accurate effective connectivity analyses over short time scales, facilitating investigations of causal interactions in the brain underlying perception and cognition during rapid transitions in behavioral state.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Generalized Multiscale Finite Element Methods with energy minimizing oversampling

In this work, we propose a general concept for constructing multiscale basis functions within Generalized Multiscale Finite Element Method, which uses oversampling and stable decomposition. The oversampling refers to using larger regions in constructing multiscale basis functions and stable decomposition allows estimating the local errors. The analysis of multiscale methods involves decomposing the error by coarse regions, where each error contribution is estimated. In this estimate, we often use oversampling techniques to achieve a fast convergence. We demonstrate our concepts in the mixed, the Interior Penalty Discontinuous Galerkin, and Hybridized Discontinuous Galerkin discretizations. One of the important features of the proposed basis functions is that they can be used in online Generalized Multiscale Finite Element Method, where one constructs multiscale basis functions using residuals. In these problems, it is important to achieve a fast convergence, which can be guaranteed if we have a stable decomposition. In our numerical results, we present examples for both offline and online multiscale basis functions. Our numerical results show that one can achieve a fast convergence when using online basis functions. Moreover, we observe that coupling using Hybridized Discontinuous Galerkin provides a better accuracy compared with Interior Penalty Discontinuous Galerkin, which is due to using multiscale glueing functions.

42 ENGINEERING↗
DOE OSTI2023

Automating methods for estimating metabolite volatility

The volatility of metabolites can influence their biological roles and inform optimal methods for their detection. Yet, volatility information is not readily available for the large number of described metabolites, limiting the exploration of volatility as a fundamental trait of metabolites. Here, we adapted methods to estimate vapor pressure from the functional group composition of individual molecules (SIMPOL.1) to predict the gas-phase partitioning of compounds in different environments. We implemented these methods in a new open pipeline called volcalc that uses chemoinformatic tools to automate these volatility estimates for all metabolites in an extensive and continuously updated pathway database: the Kyoto Encyclopedia of Genes and Genomes (KEGG) that connects metabolites, organisms, and reactions. We first benchmark the automated pipeline against a manually curated data set and show that the same category of volatility (e.g., nonvolatile, low, moderate, high) is predicted for 93% of compounds. We then demonstrate how volcalc might be used to generate and test hypotheses about the role of volatility in biological systems and organisms. Specifically, we estimate that 3.4 and 26.6% of compounds in KEGG have high volatility depending on the environment (soil vs. clean atmosphere, respectively) and that a core set of volatiles is shared among all domains of life (30%) with the largest proportion of kingdom-specific volatiles identified in bacteria. With volcalc , we lay a foundation for uncovering the role of the volatilome using an approach that is easily integrated with other bioinformatic pipelines and can be continually refined to consider additional dimensions to volatility. The volcalc package is an accessible tool to help design and test hypotheses on volatile metabolites and their unique roles in biological systems.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Scale-invariant Cosmology and CMB Temperatures as a Function of Redshifts

Cosmological models assuming the scale invariance of the macroscopic empty space show an accelerated expansion, without calling for some unknown particles. Several comparisons between models and observations (tests on distances, m−z diagram, Ω{sub Λ} versus Ω{sub m} plot, age versus H {sub 0}, H(z) versus z, transition braking-acceleration) have indicated an impressive agreement. We pursue the tests with the CMB temperatures T{sub CMB} as a function of redshifts z. CO molecules in DLA systems provide the most accurate excitation temperatures T{sub exc} up to z≈2.7. Such data need corrections for local effects like particle collisions, optical depths, UV radiation, etc., We estimate these corrections as a function of the (CO/H{sub 2}) ratios from far-UV observations of CO molecules in the Galaxy. The results show that it is not sufficient to apply theoretical collisional corrections to get the proper values of T{sub CMB} versus z. Thus, the agreement often found with the standard model may be questioned. The T{sub CMB}(z) relation needs further careful attention and the same for the scale-invariant cosmology in view of its positive tests.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Mortality attributable to PM 2.5 from wildland fires in California from 2008 to 2018

In California, wildfire risk and severity have grown substantially in the last several decades. Research has characterized extensive adverse health impacts from exposure to wildfire-attributable fine particulate matter (PM 2.5 ), but few studies have quantified long-term outcomes, and none have used a wildfire-specific chronic dose-response mortality coefficient. Here, we quantified the mortality burden for PM 2.5 exposure from California fires from 2008 to 2018 using Community Multiscale Air Quality modeling system wildland fire PM 2.5 estimates. We used a concentration-response function for PM 2.5 , applying ZIP code–level mortality data and an estimated wildfire-specific dose-response coefficient accounting for the likely toxicity of wildfire smoke. We estimate a total of 52,480 to 55,710 premature deaths are attributable to wildland fire PM 2.5 over the 11-year period with respect to two exposure scenarios, equating to an economic impact of $\$432$ to $\$456$ billion. These findings extend evidence on climate-related health impacts, suggesting that wildfires account for a greater mortality and economic burden than indicated by earlier studies.

60 APPLIED LIFE SCIENCES↗