Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model Counting”

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 19 records

A simple electron counting model for half-Heusler surfaces

Heusler compounds are a ripe platform for discovery and manipulation of emergent properties in topological and magnetic heterostructures. In these applications, the surfaces and interfaces are critical to performance; however, little is known about the atomic-scale structure of Heusler surfaces and interfaces or why they reconstruct. Using a combination of molecular beam epitaxy, core-level and angle-resolved photoemission, scanning tunneling microscopy, and density functional theory, we map the phase diagram and determine the atomic and electronic structures for several surface reconstructions of CoTiSb (001), a prototypical semiconducting half-Heusler. At low Sb coverage, the surface is characterized by Sb-Sb dimers and Ti vacancies, while, at high Sb coverage, an adlayer of Sb forms. The driving forces for reconstruction are charge neutrality and minimizing the number of Sb dangling bonds, which form metallic surface states within the bulk bandgap. We develop a simple electron counting model that explains the atomic and electronic structure, as benchmarked against experiments and first-principles calculations. We then apply the model to explain previous experimental observations at other half-Heusler surfaces, including the topological semimetal PtLuSb and the half-metallic ferromagnet NiMnSb. The model provides a simple framework for understanding and predicting the surface structure and properties of these novel quantum materials.

36 MATERIALS SCIENCE↗

Poisson Log-Normal Process for Count Data Prediction

Modeling count data is important in physics and other scientific disciplines, where measurements often involve discrete, non-negative quantities such as photon or neutrino detection events. Traditional parametric approaches can be trained to generate integer-count predictions but may struggle with capturing complex, non-linear dependencies often observed in the data. Gaussian process (GP) regression provides a robust non-parametric alternative to modeling continuous data; however, it cannot generate integer outputs. We propose the Poisson Log-Normal (PoLoN) process, a framework that employs GP to model Poisson log-rates. As in GP regression, our approach relies on the correlations between data points captured via GP kernel structure rather than explicit functional parameterizations. We demonstrate that the PoLoN predictive distribution is Poisson-LogNormal and provide an algorithm for optimizing kernel hyperparameters. Furthermore, we adapt the PoLoN approach to the problem of detecting weak localized signals superimposed on a smoothly varying background - a task of considerable interest in many areas of science and engineering. Our framework allows us to predict the strength, location and width of the detected signals. We evaluate PoLoN's performance using both synthetic and real-world datasets, including the open dataset from CERN which was used to detect the Higgs boson at the Large Hadron Collider. Our results indicate that the PoLoN process can be used as a non-parametric alternative for analyzing, predicting, and extracting signals from integer-valued data.

Saha, Anushka [Rutgers U., Piscataway]↗

A galaxy model from two micron all sky survey star counts in the whole sky, including the plane

We use the star count model of Ortiz and Lépine to perform an unprecedented exploration of the most important Galactic parameters comparing the predicted counts with the Two Micron All Sky Survey observed star counts in the J, H, and K{sub S} bands for a grid of positions covering the whole sky. The comparison is made using a grid of lines of sight given by the HEALPix pixelization scheme. The resulting best-fit values for the parameters are: 2120 ± 200 pc for the radial scale length and 205 ± 40 pc for the scale height of the thin disk, with a central hole of 2070{sub −800}{sup +2000} pc for the same disk, 3050 ± 500 pc for the radial scale length and 640 ± 70 pc for the scale height of the thick disk, 400 ± 100 pc for the central dimension of the spheroid, 0.0082 ± 0.0030 for the spheroid to disk density ratio, and 0.57 ± 0.05 for the oblate spheroid parameter.

79 ASTRONOMY AND ASTROPHYSICS↗

State-to-State Molecular Reactions in the Ultracold Regime (Final Scientific Report)

We have achieved the proposed goal to experimentally probe the AB + CD and AB + C types of reactions with state-to-state resolution, which we also compared to advanced theoretical calculations to help elucidate the role of quantum mechanics in the processes of bond breakage and formation. Our approach uses reactants that are prepared at ultracold temperatures (< 1µK) such that the quantum effects of translational motion are an important factor. Specific example reactions, including the potassium-rubidium metathesis reaction KRb + KRb → K 2 + Rb 2 as well as the atom exchange reaction Rb + KRb → Rb2 + K, are chosen because the technology of quantum internal and motional state control of these types of molecules is particularly advanced. The results for the entire funding period are fruitful. For the majority of this grant, we have constructed a one-of-the-kind quantum degenerate gas apparatus that integrates ion detection and velocity map imaging capabilities, allowing us to explore the KRb + KRb → K 2 + Rb 2 bimolecular reaction in detail. Specifically, we first verified such a reaction indeed proceed at ultracold temperatures by direct detection of reaction products. We then mapped out the complete product state distribution, which was compared to a state-counting model based on statistical theory. Our results show an overall agreement with the statistical state counting model, but also reveal several deviating state-pairs. An exact quantum calculation for molecule-molecule collisions, that is needed to understand these deviations, is however beyond the current state-of-the-art. Beside scrutinizing the reaction products, we also directly observe the reaction intermediate complex, which was quite a surprise to us. The intermediate complexes are long-lived and can interact with the inferred light that we use to trap the ultracold gas. After molecule-molecule collisions, we then explored the more theoretically tractable Rb + KRb reaction, which is endothermic. Surprisingly, we observed an exceedingly long-lived KRb$^*_2$ collisional complexes, with our experimentally measured complex lifetime deviating from conventional theoretical calculations by five orders of magnitude. This discrepancy has motivated many explorations of possible underlying causes, though no model yet captures this phenomenon completely. In the final year and the work that continues today, we extend upon these atom-molecule collision experiments to explore the origin of the long-lived KRb$^*_2$ complex lifetime and develop means to control the outcome of the reaction complex. The 5-year funded work advanced our understanding of chemical reactions at the lowest possible temperatures and at the same time opened up many new questions that are beyond our initial imaginations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Building an Efficient Cluster Cosmology Software Package for Modeling Cluster Counts and Lensing

We introduce a software suite developed for galaxy cluster cosmological analysis with the Dark Energy Survey Data. Cosmological analyses based on galaxy cluster number counts and weak-lensing measurements need efficient software infrastructure to explore an increasingly large parameter space, and account for various cosmological and astrophysical effects. Our software package is designed to model the cluster observables in a wide-field optical survey, including galaxy cluster counts, their averaged weak-lensing masses, or the cluster's averaged weak-lensing radial signals. To ensure maximum efficiency, this software package is developed in C++ in the CosmoSIS software framework, making use of the CUBA integration library. We also implement a testing and validation scheme to ensure the quality of the package. We demonstrate the effectiveness of this development by applying the software to the Dark Energy Survey Year 1 galaxy cluster cosmological data sets, and acquired cosmological constraints that are consistent with the fiducial Dark Energy Survey analysis.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Density split statistics: Joint model of counts and lensing in cells

We present density split statistics, a framework that studies lensing and counts-in-cells as a function of foreground galaxy density, thereby providing a large-scale measurement of both 2-point and 3-point statistics. Our method extends our earlier work on trough lensing and is summarized as follows: given a foreground (low redshift) population of galaxies, we divide the sky into subareas of equal size but distinct galaxy density. We then measure lensing around uniformly spaced points separately in each of these subareas, as well as counts-in-cells statistics (CiC). The lensing signals trace the matter density contrast around regions of fixed galaxy density. Through the CiC measurements this can be related to the density profile around regions of fixed matter density. Together, these measurements constitute a powerful probe of cosmology, the skewness of the density field and the connection of galaxies and matter. In this paper we show how to model both the density split lensing signal and CiC from basic ingredients: a non-linear power spectrum, clustering hierarchy coefficients from perturbation theory and a parametric model for galaxy bias and shot-noise. Furthermore, using N-body simulations, we demonstrate that this model is sufficiently accurate for a cosmological analysis on year 1 data from the Dark Energy Survey.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Time bomb or hidden treasure? Characteristics of junk TVs and of the US households who store them

Highlights: ► US households are storing 84.1 million broken or obsolete (junk) TVs. ► They represent 2.12 million metric tons of scrap. ► The value of these materials is approximately $21 per TV. ► Our count models characterize US households who store junk TVs. ► Our results are useful for designing more effective TV recycling programs. - Abstract: Within the growing stockpile of electronic waste (e-waste), TVs are especially of concern in the US because of their number (which is known imprecisely), their low recycling rate, and their material content: cathode ray tube televisions contain lead, and both rear projection and flat panel displays contain mercury, in addition to other potentially toxic materials. Based on a unique dataset from a 2010 survey, our count models show that pro-environmental behavior, age, education, household size, marital status, gender of the head of household, dwelling type, and geographic location are statistically significant variables for explaining the number of broken or obsolete (junk) TVs stored by US households. We also estimate that they are storing approximately 84.1 million junk TVs, which represents 40 pounds of scrap per household. Materials in each of these junk TVs are worth $21 on average at January 2012 materials prices, which sets an upper bound on collecting and recycling costs. This information should be helpful for developing more effective recycling strategies for TVs in the e-waste stream.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The local hole: a galaxy underdensity covering 90 per cent of sky to ≈200 Mpc

We investigate the ‘Local Hole’, an anomalous underdensity in the local galaxy environment, by extending our previous galaxy K-band number-redshift and number-magnitude counts to ≈90 per cent of the sky. Our redshift samples are taken from the 2MASS Redshift Survey (2MRS) and the 2M++ catalogues, limited to K < 11.5. We find that both surveys are in good agreement, showing an ≈21–22 per cent underdensity at z < 0.075 when compared to our homogeneous counts model that assumes the same luminosity function (LF) and other parameters as in our earlier papers. Using the Two Micron All Sky Survey (2MASS) for n(K) galaxy counts, we measure an underdensity relative to this model of 20±2 per cent at K < 11.5, which is consistent in both form and scale with the observed n(z) underdensity. To examine further the accuracy of the counts model, we compare its prediction for the fainter n(K) counts of the Galaxy and Mass Assembly (GAMA) survey. We further compare these data with a model assuming the parameters of a previous study where little evidence for the Local Hole was found. At 13 < K < 16, we find a significantly better fit for our galaxy counts model, arguing for our higher LF normalization. Although our implied underdensity of ≈20 per cent means local measurements of the Hubble Constant have been overestimated by ≈3 per cent, such a scale of underdensity is in tension with a global ΛCDM cosmology at an ≈3σ level.

79 ASTRONOMY AND ASTROPHYSICS↗

FOREGROUND PREDICTIONS FOR THE COSMIC MICROWAVE BACKGROUND POWER SPECTRUM FROM MEASUREMENTS OF FAINT INVERTED RADIO SOURCES AT 5 GHz

We present measurements of a population of matched radio sources at 1.4 and 5 GHz down to a flux limit of 1.5 mJy in 7 deg{sup 2} of the NOAO Deep Field South. We find a significant fraction of sources with inverted spectral indices that all have 1.4 GHz fluxes less than 10 mJy and are therefore too faint to have been detected and included in previous radio source count models that are matched at multiple frequencies. Combined with the matched source population at 1.4 and 5 GHz in 1 deg{sup -2} in the ATESP survey, we update models for the 5 GHz differential number counts and distributions of spectral indices in 5 GHz flux bins that can be used to estimate the unresolved point source contribution to the cosmic microwave background temperature anisotropies. We find a shallower logarithmic slope in the 5 GHz differential counts than in previously published models for fluxes {approx}< 100 mJy as well as larger fractions of inverted spectral indices at these fluxes. Because the Planck flux limit for resolved sources is larger than 100 mJy in all channels, our modified number counts yield at most a 10% change in the predicted Poisson contribution to the Planck temperature power spectrum. For a flux cut of 5 mJy with the South Pole Telescope and a flux cut of 20 mJy with the Atacama Cosmology Telescope, we predict a {approx}30% and {approx}10% increase, respectively, in the radio source Poisson power in the lowest frequency channels of each experiment relative to that predicted by previous models.

79 ASTRONOMY AND ASTROPHYSICS↗

Misclassification in Workers’ Telecommuting Frequency Choices Using a Generalized Extreme Value Model

Telecommuting frequency is a response variable collected in travel surveys and is, therefore, prone to errors leading to mismeasurements or misclassification. Misclassification of explanatory variables is a common risk when using statistical modeling techniques. We define “misclassification” as a response reported or recorded in the wrong category; for example, a variable is recorded as a 1 when it should be 0. Here, in this context, this study aims to develop a statistical model to analyze telecommuting data which accounts for potential misclassification errors by building on existing literature in econometrics. The empirical analysis was undertaken using the 2017 National Household Travel Survey (NHTS) and the general extreme value (GEV) models available in the literature. Specifically, the frequency of telecommuting days was analyzed using the negative binomial (NB) model recast as the multinomial logit (MNL) model. By nature—and consistent with other studies—NHTS data are prone to errors that can be classified as intentional or unintentional misinformation provided by the person being interviewed. Ignoring these errors while modeling telecommuting frequencies using standard discrete count models can result in biased parameter estimates. The misclassification parameter was calculated for both over-reporting and under-reporting scenarios. The misclassification errors can be as high as 14% over-reported and 10% under-reported, particularly for the neighboring values. Statistical fit comparison between the models shows that models that ignore misclassification have worse data fit and biased parameter estimates with significant policy implications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Spatial Multiplication Model as an alternative to the Point Model in Neutron Multiplicity Counting

The point model is commonly used in neutron multiplicity counting to relate the correlated neutron detection rates (singles, doubles, triples) to item properties (mass, (α,n) reaction rate and neutron multiplication). The point model assumes that the probability that a neutron will induce fission is a constant across the physical extent of the item. However, in reality, neutrons near the center of an item have a greater probability of inducing fission then items near the edges. As a result, the neutron multiplication has a spatial distribution.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scatter and crosstalk corrections for {sup 99m}Tc/{sup 123}I dual-radionuclide imaging using a CZT SPECT system with pinhole collimators

Purpose: The energy spectrum for a cadmium zinc telluride (CZT) detector has a low energy tail due to incomplete charge collection and intercrystal scattering. Due to these solid-state detector effects, scatter would be overestimated if the conventional triple-energy window (TEW) method is used for scatter and crosstalk corrections in CZT-based imaging systems. The objective of this work is to develop a scatter and crosstalk correction method for {sup 99m}Tc/{sup 123}I dual-radionuclide imaging for a CZT-based dedicated cardiac SPECT system with pinhole collimators (GE Discovery NM 530c/570c). Methods: A tailing model was developed to account for the low energy tail effects of the CZT detector. The parameters of the model were obtained using {sup 99m}Tc and {sup 123}I point source measurements. A scatter model was defined to characterize the relationship between down-scatter and self-scatter projections. The parameters for this model were obtained from Monte Carlo simulation using SIMIND. The tailing and scatter models were further incorporated into a projection count model, and the primary and self-scatter projections of each radionuclide were determined with a maximum likelihood expectation maximization (MLEM) iterative estimation approach. The extracted scatter and crosstalk projections were then incorporated into MLEM image reconstruction as an additive term in forward projection to obtain scatter- and crosstalk-corrected images. The proposed method was validated using Monte Carlo simulation, line source experiment, anthropomorphic torso phantom studies, and patient studies. The performance of the proposed method was also compared to that obtained with the conventional TEW method. Results: Monte Carlo simulations and line source experiment demonstrated that the TEW method overestimated scatter while their proposed method provided more accurate scatter estimation by considering the low energy tail effect. In the phantom study, improved defect contrasts were observed with both correction methods compared to no correction, especially for the images of {sup 99m}Tc in dual-radionuclide imaging where there is heavy contamination from {sup 123}I. In this case, the nontransmural defect contrast was improved from 0.39 to 0.47 with the TEW method and to 0.51 with their proposed method and the transmural defect contrast was improved from 0.62 to 0.74 with the TEW method and to 0.73 with their proposed method. In the patient study, the proposed method provided higher myocardium-to-blood pool contrast than that of the TEW method. Similar to the phantom experiment, the improvement was the most substantial for the images of {sup 99m}Tc in dual-radionuclide imaging. In this case, the myocardium-to-blood pool ratio was improved from 7.0 to 38.3 with the TEW method and to 63.6 with their proposed method. Compared to the TEW method, the proposed method also provided higher count levels in the reconstructed images in both phantom and patient studies, indicating reduced overestimation of scatter. Using the proposed method, consistent reconstruction results were obtained for both single-radionuclide data with scatter correction and dual-radionuclide data with scatter and crosstalk corrections, in both phantom and human studies. Conclusions: The authors demonstrate that the TEW method leads to overestimation in scatter and crosstalk for the CZT-based imaging system while the proposed scatter and crosstalk correction method can provide more accurate self-scatter and down-scatter estimations for quantitative single-radionuclide and dual-radionuclide imaging.

60 APPLIED LIFE SCIENCES↗

Meso-modeling of Carbon Fiber Composite for Crash Safety Analysis

In the conventional approach, the material properties for crash safety simulations are typically obtained from standard coupon tests, where the test results only provide single layer material properties used in crash simulations. However, the lay-up effects for the failure behaviors of the real structure were not considered in numerical simulations. Hence, there was discrepancy between the crash simulations and experimental tests. Consequently, an intermediate stage is required for accurate predictions. Some component tests are required to correlate the material models in the intermediate stage. In this paper, a Mazda Tube under high-impact velocity is chosen as an example for the crash safety analysis. The tube consists of 24 layers of uni-directional (UD) carbon fiber composite materials, in which 4 layers are perpendicular to, while the other layers are parallel to the impact direction. An LS-DYNA meso-model was constructed with orthotropic material models counting for the single-layer material behaviors. Between layers, a node-based tie-break contact was used for modeling the delamination of the composite material. Since fiber directions are not single-oriented, the lay-up effects could be an important effect. From the first numerical trial, premature material failure occurred due to the use of material parameters obtained directly from the coupon tests. Some parametric studies were conducted to identify the cause of the numerical instability. The finding is that the material failure strength used in the numerical model needs to be enlarged to stabilize the numerical model. Some hypothesis was made to provide the foundation for enlarging the failure strength and the corresponding experiments will be conducted to validate the hypothesis.

Lin, Shih-Po↗

Extracting the gamma-ray source-count distribution below the Fermi-LAT detection limit with deep learning

We reconstruct the extra-galactic gamma-ray source-count distribution, or dN/dS, of resolved and unresolved sources by adopting machine learning techniques. Specifically, we train a convolutional neural network on synthetic 2-dimensional sky-maps, which are built by varying parameters of underlying source-counts models and incorporate the Fermi-LAT instrumental response functions. The trained neural network is then applied to the Fermi-LAT data, from which we estimate the source count distribution down to flux levels a factor of 50 below the Fermi-LAT threshold. We perform our analysis using 14 years of data collected in the (1,10) GeV energy range. The results we obtain show a source count distribution which, in the resolved regime, is in excellent agreement with the one derived from cataloged sources, and then extends as dN/dS ~ S -2 in the unresolved regime, down to fluxes of 5 · 10 -12 cm -2 s -1 . The neural network architecture and the devised methodology have the flexibility to enable future analyses to study the energy dependence of the source-count distribution.

79 ASTRONOMY AND ASTROPHYSICS↗

Practical Probabilistic Programming

Recent advances in probabilistic programming languages (PPLs) have provided the capability for exact inference: computing a closed-form probability distribution for a given probabilistic program. In particular, the new language Roulette uses a language oriented programming (LOP) approach, wherein analysts build new programming languages on top of a set of primitives provided by Roulette, which then translates these structures into a weighted model counting problem which can be solved by automated reasoning tools. However, because Roulette provides few convenience features, developing these new languages is challenging even for expert users. We developed a standard library of common probability functions for Roulette with the goal of improved usability. This included approximation of continuous probability density functions using discrete probability mass functions. We demonstrated this approach by modeling a cosmic ray striking a RAM controller. We found that Roulette provides a powerful interface for highly expressive probabilistic programs to be generated. In collaboration with the NNSA Advanced Simulation and Computing program, which resulted in development of a tool called Circulette, we were able to model complex circuits expressed in Verilog using probabilistic programs with an expressivity not previously possible. Our research question that motivated the development of a Roulette standard library was to determine whether non-experts could use a PPL to model relevant problems regarding radiation effects on microelectronics. This standard library improved the expressivity of Roulette by implementing common probability density functions, mathematical operators on distributions, and support for empirical distributions. While Roulette is a powerful modeling language, the untyped, LOP approach makes error messages difficult to understand and requires expert aid. We recommend further research on Roulette, especially with its error messages, to enable improved usability. At the same time, this project demonstrated that for users familiar with Roulette and the LOP approach, Roulette provides powerful new capabilities that can be integrated with other Sandia modeling capabilities.

97 MATHEMATICS AND COMPUTING↗

Precision test of statistical dynamics with state-to-state ultracold chemistry

Chemical reactions represent a class of quantum problems that challenge both the current theoretical understanding and computational capabilities1. Reactions that occur at ultralow temperatures provide an ideal testing ground for quantum chemistry and scattering theories, because they can be experimentally studied with unprecedented control, yet display dynamics that are highly complex. Here we report the full product state distribution for the reaction 2KRb → K 2 + Rb 2 . Ultracold preparation of the reactants allows us complete control over their initial quantum degrees of freedom, whereas state-resolved, coincident detection of both products enables the probability of scattering into each of the 57 allowed rotational state-pairs to be measured. Our results show an overall agreement with a state-counting model based on statistical theory, but also reveal several deviating state-pairs. In particular, we observe a strong suppression of population in the state-pair closest to the exoergicity limit as a result of the long-range potential inhibiting the escape of products. Finally, the completeness of our measurements provides a benchmark for quantum dynamics calculations beyond the current state of the art.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Switchable electronic and enhanced magnetic properties of CrI 3 edges

Owing to its novel electronic and magnetic properties, two-dimensional CrI 3 has great potential in the application of spintronic devices. However, as an inevitable line defect, the properties of the edges of CrI 3 remain elusive. Here, via first-principles calculations with spin–orbit coupling, we investigated the thermodynamic stabilities, electronic and magnetic properties of thirteen CrI 3 edges with different structures. We showed that zigzag edges are more stable than armchair edges, and a CrI 3 nanoribbon can be either metallic or insulating depending on its chemical growth conditions. Here, the edge stability and associated electronic properties can be understood in terms of the octahedron ligand field and electron counting model. In most cases, both the magnetic moment and Curie temperature can be enhanced by edges, which are in startle contrast to the surfaces of three-dimensional ferromagnetic materials, where a magnetic dead layer is often observed.

36 MATERIALS SCIENCE↗