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 217 records · Page 12

IDENTIFYING STAR STREAMS IN THE MILKY WAY HALO

We develop statistical methods for identifying star streams in the halo of the Milky Way that exploit observed spatial and radial velocity distributions. Within a great circle, departures of the observed spatial distribution from random provide a measure of the likelihood of a potential star stream. Comparisons between the radial velocity distribution within a great circle and the radial velocity distribution of the entire sample also measure the statistical significance of potential streams. The radial velocities enable construction of a more powerful joint statistical test for identifying star streams in the Milky Way halo. Applying our method to halo stars in the Hypervelocity Star (HVS) survey, we detect the Sagittarius stream at high significance. Great circle counts and comparisons with theoretical models suggest that the Sagittarius stream comprises 10%-17% of the halo stars in the HVS sample. The population of blue stragglers and blue horizontal branch stars varies along the stream and is a potential probe of the distribution of stellar populations in the Sagittarius dwarf galaxy prior to disruption.

79 ASTRONOMY AND ASTROPHYSICS↗

Quantifying mass segregation and new core radii for 54 Milky Way globular clusters

We present core radii for 54 Milky Way globular clusters determined by fitting King-Michie models to cumulative projected star count distributions. We find that fitting star counts rather than surface brightness profiles produces results that differ significantly due to the presence of mass segregation. The sample in each cluster is further broken down into various mass groups, each of which is fit independently, allowing us to determine how the concentration of each cluster varies with mass. The majority of the clusters in our sample show general agreement with the standard picture that more massive stars will be more centrally concentrated. We find that core radius versus stellar mass can be fit with a two-parameter power law. The slope of this power law is a value that describes the amount of mass segregation present in the cluster, and is measured independently of our distance from the cluster. This value correlates strongly with the core relaxation time and physical size of each cluster. Supplementary figures are also included showing the best fits and likelihood contours of fit parameters for all 54 clusters.

79 ASTRONOMY AND ASTROPHYSICS↗

An empirical method for improving the quality of RXTE PCA spectra

We fitted all of the several hundred RXTE PCA spectra of the Crab individually to a simple power-law model; the total number of counts in the composite spectrum is >10{sup 9}. We then used the spectrum of residuals to derive a calibration tool, called pcacorr, that we apply to large samples of spectra for GX 339-4, H1743-322, and XTE J1550-564. Application of the tool improved the quality of all the fits, and the improvement is dramatic for spectra with ≳ 10{sup 7} counts. The Crab residual spectrum is somewhat different for each of the five Proportional Counter Array (PCA) detectors, but it was relatively stable over the course of the mission. We recommend that pcacorr be routinely applied to spectra with ≳ 10{sup 6} counts and that one include a systematic error of 0.1%, rather than the 0.5-1% value that has customarily been used. We expect that application of the tool will result in an increase in sensitivity of the PCA to faint spectral features by up to an order of magnitude.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-D-BRF-01: Applications of Photonuclear Activation of Biological Tissues in Clinical High-Energy X-Ray Beams

Purpose: The commercial availability of high-energy accelerators opens new therapeutic opportunities for X-rays with energies in excess of 15MV. Three clinical beams (Varian 18MV, Elekta SL25, and Top Grade LA45) were compared by the production of photo-activated positron emitters in five types of biological tissues. Methods: The activation studies were performed using FLUKA2011 Monte Carlo suite with beam models designed for the three accelerators. Absolute activity density (Bq/ml) distribution in space was obtained from this study. Additionally, the temporal evolution of all activated species was monitored for cooling times up to 30 minutes. Results: The relative activation contrast of tissue pairs was evaluated for 2Gy of dose at 10cm depth for five tissues (normal, hypoxic, adipose, bone, and lung). In bone the sort-lived isotopes O-15 and P-30 dominated the activity at the early cooling stages. In all tissues 15 minutes post-irradiation the C-11 activity became dominant. Tissues with higher carbon-to-oxygen ratio in their elemental composition became clearly visible in a PET scan at longer cooling times. Radiation treatment of a lung tumor with a hypoxic core was simulated in an anthropomorphic phantom using the LA45 beam model. Increase in the PET counts by more than 50% was measured in the hypoxic volume 15 minutes post-irradiation, demonstrating the benefits of the relative activation contrast concept. Conclusion: LA45 beam is found to produce measurable activation distribution for 2Gy of dose, suitable for tissue type discrimination studies. Maximum activation contrast of 1.7 between hypoxic and normal tissues was measured in a simulated treatment. The absolute activity density values obtained from these Monte Carlo studies suggest that the window of opportunity for a PET scan exists up to 60 minutes after 2Gy of dose is deposited by 45MV X-ray beam. Lung-to-tissue activation contrast can be explored for treatment QA purposes as well.

60 APPLIED LIFE SCIENCES↗

Pseudogap behavior in the emery model for the electron-doped superconductor Nd{sub 2-x}Ce{sub x}CuO{sub 4}: Multiband LDA + DMFT + {Sigma}{sub k} approach

We propose a generalization of the LDA + DMFT + {Sigma}{sub k} approach to the multiband case, in which correlated and uncorrelated states are present in the model simultaneously. Using the multiband version of the LDA + DMFT + {Sigma}{sub k} approach, we calculate the density of states and spectral functions for the Emery model in a wide energy interval around the Fermi level. We also obtain the Fermi surfaces for the electron-doped high-temperature superconductor Nd{sub 2-x}Ce{sub x}CuO{sub 4} in the pseudogap phase. The self-energy part {Sigma}{sub k} introduced additionally to take into account pseudogap fluctuations describes the nonlocal interaction of correlated electrons with collective Heisenberg short-range spin fluctuations. To solve the effective impurity model, the numerical renorm-group (NRG) method is used for the DMFT equations. Good qualitative agreement of the Fermi surfaces calculated using the LDA + DMFT + {Sigma}{sub k} approach and experimental angle-resolved photoemission spectroscopic data is attained. The stability of the dielectric solution with charge transfer in the Emery model with correction for double counting is analyzed in the Appendix.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

INTERPRETING FLUX FROM BROADBAND PHOTOMETRY

We discuss the transformation of observed photometry into flux for the creation of spectral energy distributions (SED) and the computation of bolometric luminosities. We do this in the context of supernova studies, particularly as observed with the Swift spacecraft, but the concepts and techniques should be applicable to many other types of sources and wavelength regimes. Traditional methods of converting observed magnitudes to flux densities are not very accurate when applied to UV photometry. Common methods for extinction and the integration of pseudo-bolometric fluxes can also lead to inaccurate results. The sources of inaccuracy, though, also apply to other wavelengths. Because of the complicated nature of translating broadband photometry into monochromatic flux densities, comparison between observed photometry and a spectroscopic model is best done by forward modeling the spectrum into the count rates or magnitudes of the observations. We recommend that integrated flux measurements be made using a spectrum or SED which is consistent with the multi-band photometry rather than converting individual photometric measurements to flux densities, linearly interpolating between the points, and integrating. We also highlight some specific areas where the UV flux can be mischaracterized.

79 ASTRONOMY AND ASTROPHYSICS↗

An Integral-based Technique to Accelerate the Monte Carlo Radiative Transfer Computation for Supernovae

We present an integral-based technique (IBT) algorithm to accelerate supernova (SN) radiative transfer calculations. The algorithm utilizes “integral packets,” which are calculated by the path integral of the Monte Carlo (MC) energy packets, to synthesize the observed spectropolarimetric signal at a given viewing direction in a 3D time-dependent radiative transfer program. Compared to the event-based technique (EBT) proposed by M. Bulla et al., our algorithm significantly reduces the computation time and increases the MC signal-to-noise ratio (S/N). Using a 1D spherical symmetric Type Ia SN ejecta model DDC10 and its derived 3D model, the IBT algorithm has successfully passed the verification of spherical symmetry and cross comparison on a 3D SN model with the direct-counting technique and EBT. Notably, with our algorithm implemented in the 3D MC radiative transfer code SEDONA, the computation time is faster than EBT by a factor of 10−30, and the S/N is better by a factor of 1.5−3, with the same number of MC quanta.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimation of the Performance of Multiple Active Neutron Interrogation Signatures for Detecting Shielded HEU

A comprehensive modeling study has been carried out to evaluate the utility of multiple active neutron interrogation signatures for detecting shielded highly enriched uranium (HEU). The modeling effort focused on varying HEU masses from 1 kg to 20 kg; varying types of shields including wood, steel, cement, polyethylene, and borated polyethylene; varying depths of the HEU in the shields, and varying engineered shields immediately surrounding the HEU including steel, tungsten, and cadmium. Neutron and gamma-ray signatures were the focus of the study and false negative detection probabilities versus measurement time were used as a performance metric. To facilitate comparisons among different approaches an automated method was developed to generate receiver operating characteristic (ROC) curves for different sets of model variables for multiple background count rate conditions. This paper summarizes results or the analysis, including laboratory benchmark comparisons between simulations and experiments. The important impact engineered shields can play towards degrading detectability and methods for mitigating this will be discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Locating Equitable Solar Opportunities by Census Tract: A Guide to the Screening Tool for Equitable Adoption and Deployment of Solar (STEADy Solar)

The Screening Tool for Equitable Adoption and DeploYment of Solar (STEADy Solar) is a database and mapping tool that indicates locations that may be eligible for the Investment Tax Credit bonus adders defined in the 2022 Inflation Reduction Act (IRA). The tool combines publicly available information on demographics, solar technical potential, solar economics (modeled net present value), building counts by use-type, and eligibility for tax credit adders. It can be used by states, municipalities, community-based organizations, developers, and researchers to identify sites where solar projects may be economical and where federal incentives may be available to support equitable adoption of solar. This report describes the STEADy dataset and presents high level insights from the data.

census tract↗

Pulse pileup statistics for energy discriminating photon counting x-ray detectors

Purpose: Energy discriminating photon counting x-ray detectors can be subject to a wide range of flux rates if applied in clinical settings. Even when the incident rate is a small fraction of the detector's maximum periodic rate N{sub 0}, pulse pileup leads to count rate losses and spectral distortion. Although the deterministic effects can be corrected, the detrimental effect of pileup on image noise is not well understood and may limit the performance of photon counting systems. Therefore, the authors devise a method to determine the detector count statistics and imaging performance. Methods: The detector count statistics are derived analytically for an idealized pileup model with delta pulses of a nonparalyzable detector. These statistics are then used to compute the performance (e.g., contrast-to-noise ratio) for both single material and material decomposition contrast detection tasks via the Cramer-Rao lower bound (CRLB) as a function of the detector input count rate. With more realistic unipolar and bipolar pulse pileup models of a nonparalyzable detector, the imaging task performance is determined by Monte Carlo simulations and also approximated by a multinomial method based solely on the mean detected output spectrum. Photon counting performance at different count rates is compared with ideal energy integration, which is unaffected by count rate. Results: The authors found that an ideal photon counting detector with perfect energy resolution outperforms energy integration for our contrast detection tasks, but when the input count rate exceeds 20%N{sub 0}, many of these benefits disappear. The benefit with iodine contrast falls rapidly with increased count rate while water contrast is not as sensitive to count rates. The performance with a delta pulse model is overoptimistic when compared to the more realistic bipolar pulse model. The multinomial approximation predicts imaging performance very close to the prediction from Monte Carlo simulations. The monoenergetic image with maximum contrast-to-noise ratio from dual energy imaging with ideal photon counting is only slightly better than with dual kVp energy integration, and with a bipolar pulse model, energy integration outperforms photon counting for this particular metric because of the count rate losses. However, the material resolving capability of photon counting can be superior to energy integration with dual kVp even in the presence of pileup because of the energy information available to photon counting. Conclusions: A computationally efficient multinomial approximation of the count statistics that is based on the mean output spectrum can accurately predict imaging performance. This enables photon counting system designers to directly relate the effect of pileup to its impact on imaging statistics and how to best take advantage of the benefits of energy discriminating photon counting detectors, such as material separation with spectral imaging.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Downfolding from ab initio to interacting model Hamiltonians: comprehensive analysis and benchmarking of the DFT+cRPA approach

Abstract Model Hamiltonians are regularly derived from first principles to describe correlated matter. However, the standard methods for this contain a number of largely unexplored approximations. For a strongly correlated impurity model system, here we carefully compare a standard downfolding technique with the best possible ground-truth estimates for charge-neutral excited-state energies and wave functions using state-of-the-art first-principles many-body wave function approaches. To this end, we use the vanadocene molecule and analyze all downfolding aspects, including the Hamiltonian form, target basis, double-counting correction, and Coulomb interaction screening models. We find that the choice of target-space basis functions emerges as a key factor for the quality of the downfolded results, while orbital-dependent double-counting corrections diminish the quality. Background screening of the Coulomb interaction matrix elements primarily affects crystal-field excitations. Our benchmark uncovers the relative importance of each downfolding step and offers insights into the potential accuracy of minimal downfolded model Hamiltonians.

Chemistry↗

Using the D-Wave 2X Quantum Computer to Explore the Formation of Global Terrorist Networks

Social networks with signed edges (+/-) play an important role in an area of social network theory called structural balance. In these networks, edges represent relationships that are labeled as either friendly (+) or hostile (-). A signed social network is balanced only if all cycles of three or more nodes in the graph have an odd number of hostile edges. A fundamental property of a balanced network is that it can be cleanly divided into 2 factions, where all relationships within each faction are friendly, and all relationships between members of different factions are hostile. The more unbalanced a network is, the more edges will fail to adhere to this rule, making factions more ambiguous. Social theory suggests unbalanced networks should be unstable, a finding that has been supported by research on gangs, which shows that unbalanced relationships are associated with greater violence, possibly due to this increased ambiguity about factional allegiances (Nakamura et al). One way to estimate the imbalance in a network, if only edge relationships are known, is to assign nodes to factions that minimize the number of violations of the edge rule described above. This problem is known to be computationally NP-hard. However, Facchetti et al. have pointed out that it is equivalent to an Ising model with a Hamiltonian that effectively counts the number of edge rule violations. Therefore, finding the assignment of factions that minimizes energy of the equivalent Ising system yields an estimate of the imbalance in the network. Based on the Ising model equivalence of the signed-social network balance problem, we have used the D-Wave 2X quantum annealing computer to explore some aspects of signed social networks. Because connectivity in the D-Wave computer is limited to its particular native topology, arbitrary networks cannot be represented directly. Rather, they must be “embedded” using a technique in which multiple qubits are chained together with special weights to simulate a collection of nodes with the required connectivity. This limits the size of a fully connected network in the D-Wave to about 50 simulated nodes, using all of the approximately 1150 qubits in the machine. In order to keep within this limitation, while exploring a problem of potential social relevance, we constructed time series of historical network snapshots from Stanford’s Mapping Militants Project, where nodes represent militant organizations, and edges represent either alliances or rivalries between organizations. We constructed two series from different theaters – Iraq and Syria – spanning timelines from about 2000 to 2016, each with networks whose maximum size was in the 20-30 node range. Computationally, our experience suggests D-Wave technology is promising, providing fast, nearly constant scaling of computational effort in the main part of the calculation that relies on the quantum annealing cycle. However, the cost of embedding an arbitrary network of interest in the D-Wave native topology scales poorly. If the embedding cost can be amortized relative to the annealing cycle, it may be possible to gain a substantial advantage over classical computing methods, provided a large enough network can be accommodated by partitioning into subnetworks or some similar strategy. In terms of our application to networks of militant organizations, we found a rise in network imbalance in the Syrian theater that appears to correspond roughly with the entrance of the Islamic State into a milieu already populated with other groups, a phenomenon we plan to explore in more detail. In these very preliminary results, we also noticed that during at least one period where both the size and imbalance of the network increased substantially, the imbalance per edge seemed to remain fairly steady. This may suggest some adaptive behavior among the participating factions, which may also warrant further exploration.

97 MATHEMATICS AND COMPUTING↗

Uranium particle age dating, aggregation, and model age best estimators

We present important aspects of uranium particle age dating by Large-Geometry Secondary Ion Mass Spectrometry (LG-SIMS) that can introduce bias and increase model age uncertainties, especially for small, young, and/or low-enriched particles. This metrology is important for applications related to International Nuclear Safeguards. We explore influential factors related to model age estimation, including the effects of evolving surface chemistry on inter-element measurements of particles (e.g., Th and U), detector background, and aggregation methods using simulated and actual particle samples. We introduce a new model age estimator, called “mid68”, that supplements 95% confidence intervals, providing a “best estimate” and uncertainty about the most likely age. The mid68 estimator can be calculated using the Feldman and Cousins method or Bayesian methods and provides a value with a symmetric uncertainty that can be used for calculations and approximate aggregation of processed model age values when the raw data and correction factors are not available. For particles yielding low 230 Th counts amidst nonzero detector background, their underlying model age probability distributions are asymmetric, so the mid68 estimator provides additional robust information regarding the underlying model age likelihood. This study provides a comprehensive and timely examination of critical aspects of uranium particle age dating as more laboratories establish particle chronometry capabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Wedge sampling for computing clustering coefficients and triangle counts on large graphs

Graphs are used to model interactions in a variety of contexts, and there is a growing need to quickly assess the structure of such graphs. Some of the most useful graph metrics are based on triangles, such as those measuring social cohesion. Despite the importance of these triadic measures, algorithms to compute them can be extremely expensive. We discuss the method of wedge sampling. This versatile technique allows for the fast and accurate approximation of various types of clustering coefficients and triangle counts. Furthermore, these techniques are extensible to counting directed triangles in digraphs. Our methods come with provable and practical time-approximation tradeoffs for all computations. We provide extensive results that show our methods are orders of magnitude faster than the state of the art, while providing nearly the accuracy of full enumeration.

97 MATHEMATICS AND COMPUTING↗

Bounds on galaxy stochasticity from halo occupation distribution modeling

The joint probability distribution of matter overdensity and galaxy counts in cells is a powerful probe of cosmology, and the extent to which variance in galaxy counts at fixed matter density deviates from Poisson shot noise is not fully understood. The lack of informed bounds on this stochasticity is currently the limiting factor in constraining cosmology with the galaxy–matter probability distribution function (PDF). We investigate stochasticity in the conditional distribution of galaxy counts along lines of sight with fixed matter density, and we present a halo occupation distribution (HOD)-based approach for obtaining plausible ranges for stochasticity parameters. To probe the high-dimensional space of possible galaxy–matter connections, we derive a set of HODs that conserve the galaxies’ linear bias and number density to produce RED M A G I C-like galaxy catalogs within the A BACUS S UMMIT suite of N -body simulations. We study the impact of individual HOD parameters and cosmology on stochasticity and perform a Monte Carlo search in HOD parameter space subject to the constraints on bias and density. In mock catalogs generated by the selected HODs, shot noise in galaxy counts spans both sub-Poisson and super-Poisson values, ranging from 80% to 133% of Poisson variance for cells with mean matter density. Nearly all of the derived HODs show a positive relationship between local matter density and stochasticity. For galaxy catalogs with higher stochasticity, modeling galaxy bias to second order is required for an accurate description of the conditional PDF of galaxy counts at fixed matter density. The presence of galaxy assembly bias also substantially extends the range of stochasticity in the super-Poisson direction. This HOD-based approach leverages degrees of freedom in the galaxy–halo connection to obtain informed bounds on nuisance model parameters and can be adapted to study other parametrizations of shot noise in galaxy counts, in particular to motivate prior ranges on stochasticity for cosmological analyses.

Britt, Dylan (ORCID:000000019905601X)↗

Haffner 16 Redux: Revisiting a Young Cluster in the Outer Galaxy

Images and spectra recorded with the Gemini Multi-Object Spectrograph on Gemini South are used to investigate the stellar content of the open cluster Haffner 16. The (i{sup ′},g{sup ′}−i{sup ′}) color–magnitude diagram (CMD) constructed from these data extends over 10 mag in i{sup ′}, sampling the cluster main sequence (MS) and 5 mag of the pre-MS (PMS). The fraction of unresolved equal mass binaries among PMS stars is estimated to be 0.6 ± 0.1. The isochrones do not track the PMS on the CMD, in the sense that the PMS has a shallower slope on the CMD than predicted by the models. Still, a dip in star counts, which is associated with the relaxation of PMS stars onto the MS, is identified near i{sup ′}=17. The depth and brightness of this feature—as well as the morphology of the cluster MS on the CMD—are matched by models with a slightly sub-solar metallicity that have an age of ∼20 Myr and a distance modulus of 12.3 ± 0.2. A light profile of Haffner 16 is constructed in the W1 filter (λ{sub cen}=3.4 μm), which suggests that the cluster is surrounded by a diffuse stellar halo. Spectra of candidate cluster MS and PMS stars selected according to location on the CMD are presented. The spectra show characteristics that are suggestive of a sub-solar metallicity. Hα emission is common among objects on the PMS locus on the CMD near i{sup ′}=18. It is suggested that the location of the Haffner 16 PMS on the CMD is affected by large-scale cool spot activity, likely induced by rapid stellar rotation.

79 ASTRONOMY AND ASTROPHYSICS↗

Implementation of a Portable HPGe for Field Contamination Assay

In this work, using MCNP to construct a detector model based initially on x-ray images of a portable high purity germanium (HPGe) detector followed by normalizing covering material values to also agree with check source responses, a validation of the model was attained. By calibrating the detector parameters using large count spectra, rigorous reproducibility is attained for high activity measurements but does not prevent deviations from normality in error distributions at the very low count events where spectral peaks are not always identifiable. The resulting model was created to allow operational assay of contamination over large areal distributions that could not otherwise be measured, such as the exhaust shaft at the Waste Isolation Pilot Plant (WIPP). Results indicate that contamination levels of activity in the exhaust shaft can be assayed to within a factor of 2. Lastly, detection limits are evaluated to be well below the contamination levels, which would constitute a legal environmental release if unfiltered ventilation of the underground facility were used.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

RGG: Reactor geometry (and mesh) generator

The reactor geometry (and mesh) generator RGG takes advantage of information about repeated structures in both assembly and core lattices to simplify the creation of geometry and mesh. It is released as open source software as a part of the MeshKit mesh generation library. The methodology operates in three stages. First, assembly geometry models of various types are generated by a tool called AssyGen. Next, the assembly model or models are meshed by using MeshKit tools or the CUBIT mesh generation tool-kit, optionally based on a journal file output by AssyGen. After one or more assembly model meshes have been constructed, a tool called CoreGen uses a copy/move/merge process to arrange the model meshes into a core model. In this paper, we present the current state of tools and new features in RGG. We also discuss the parallel-enabled CoreGen, which in several cases achieves super-linear speedups since the problems fit in available RAM at higher processor counts. Several RGG applications - 1/6 VHTR model, 1/4 PWR reactor core, and a full-core model for Monju - are reported. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗