Search NASA⌕ Search

SEARCH · Search NASA

Results for “density 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 145 records · Page 8

Software For Multivariate Bayesian Classification

PHD general-purpose classifier computer program. Uses Bayesian methods to classify vectors of real numbers, based on combination of statistical techniques that include multivariate density estimation, Parzen density kernels, and EM (Expectation Maximization) algorithm. By means of simple graphical interface, user trains classifier to recognize two or more classes of data and then use it to identify new data. Written in ANSI C for Unix systems and optimized for online classification applications. Embedded in another program, or runs by itself using simple graphical-user-interface. Online help files makes program easy to use.

Saul, Ronald↗

ESTIMATES OF DENSITIES AND FILLING FACTORS FROM A COOLING TIME ANALYSIS OF SOLAR MICROFLARES OBSERVED WITH RHESSI

We use more than 4500 microflares from the RHESSI microflare data set to estimate electron densities and volumetric filling factors of microflare loops using a cooling time analysis. We show that if the filling factor is assumed to be unity, the calculated conductive cooling times are much shorter than the observed flare decay times, which in turn are much shorter than the calculated radiative cooling times. This is likely unphysical, but the contradiction can be resolved by assuming that the radiative and conductive cooling times are comparable, which is valid when the flare loop temperature is a maximum and when external heating can be ignored. We find that resultant radiative and conductive cooling times are comparable to observed decay times, which has been used as an assumption in some previous studies. The inferred electron densities have a mean value of 10{sup 11.6} cm{sup -3} and filling factors have a mean of 10{sup -3.7}. The filling factors are lower and densities are higher than previous estimates for large flares, but are similar to those found for two microflares by Moore et al.

79 ASTRONOMY AND ASTROPHYSICS↗

High-dimensional and permutation invariant anomaly detection

Methods for anomaly detection of new physics processes are often limited to low-dimensional spaces due to the difficulty of learning high-dimensional probability densities. Particularly at the constituent level, incorporating desirable properties such as permutation invariance and variable-length inputs becomes difficult within popular density estimation methods. In this work, we introduce a permutation-invariant density estimator for particle physics data based on diffusion models, specifically designed to handle variable-length inputs. We demonstrate the efficacy of our methodology by utilizing the learned density as a permutation-invariant anomaly detection score, effectively identifying jets with low likelihood under the background-only hypothesis. To validate our density estimation method, we investigate the ratio of learned densities and compare to those obtained by a supervised classification algorithm.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Estimates of Ta strength at ultrahigh pressures and strain rates using thin-film graded-density impactors

In this paper, we present results from an experimental technique used to estimate the strength of Ta at extreme pressures (150 GPa) and strain rates (10 7 s -1 ). A graded-density impactor (GDI) was fabricated using sputter deposition to produce an approximately 40-μm-thick film containing alternating layers of Al and Cu. The thicknesses of the respective layers are adjusted to give an effective density gradient through the film. The GDIs were launched with a 2-stage light gas gun, and shock-ramp-release velocity profiles were measured over timescales of ~10 ns. Results are presented for the direct impact of the film onto LiF windows, which allows for a dynamic characterization of the GDI, as well as from impact onto thin (~40μm) sputtered Ta samples backed by a LiF window. The measurements were coupled with mesoscale numerical simulations to infer the strength of Ta, and the results agree well with other high-pressure platforms, particularly when strain-rate, microstructural, and thermodynamic-path differences are considered.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Robust Estimation of Electron Density From Anatomic Magnetic Resonance Imaging of the Brain Using a Unifying Multi-Atlas Approach

Purpose: To develop a reliable method to estimate electron density based on anatomic magnetic resonance imaging (MRI) of the brain. Methods and Materials: We proposed a unifying multi-atlas approach for electron density estimation based on standard T1- and T2-weighted MRI. First, a composite atlas was constructed through a voxelwise matching process using multiple atlases, with the goal of mitigating effects of inherent anatomic variations between patients. Next we computed for each voxel 2 kinds of conditional probabilities: (1) electron density given its image intensity on T1- and T2-weighted MR images; and (2) electron density given its spatial location in a reference anatomy, obtained by deformable image registration. These were combined into a unifying posterior probability density function using the Bayesian formalism, which provided the optimal estimates for electron density. We evaluated the method on 10 patients using leave-one-patient-out cross-validation. Receiver operating characteristic analyses for detecting different tissue types were performed. Results: The proposed method significantly reduced the errors in electron density estimation, with a mean absolute Hounsfield unit error of 119, compared with 140 and 144 (P<.0001) using conventional T1-weighted intensity and geometry-based approaches, respectively. For detection of bony anatomy, the proposed method achieved an 89% area under the curve, 86% sensitivity, 88% specificity, and 90% accuracy, which improved upon intensity and geometry-based approaches (area under the curve: 79% and 80%, respectively). Conclusion: The proposed multi-atlas approach provides robust electron density estimation and bone detection based on anatomic MRI. If validated on a larger population, our work could enable the use of MRI as a primary modality for radiation treatment planning.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Dwarf galaxies with ionizing radiation feedback. II. Spatially resolved star formation relation

We investigate the spatially resolved star formation relation using a galactic disk formed in a comprehensive high-resolution (3.8 pc) simulation. Our new implementation of stellar feedback includes ionizing radiation as well as supernova explosions, and we handle ionizing radiation by solving the radiative transfer equation rather than by a subgrid model. Photoheating by stellar radiation stabilizes gas against Jeans fragmentation, reducing the star formation rate (SFR). Because we have self-consistently calculated the location of ionized gas, we are able to make simulated, spatially resolved observations of star formation tracers, such as Hα emission. We can also observe how stellar feedback manifests itself in the correlation between ionized and molecular gas. Applying our techniques to the disk in a galactic halo of 2.3 × 10{sup 11} M {sub ☉}, we find that the correlation between SFR density (estimated from mock Hα emission) and H{sub 2} density shows large scatter, especially at high resolutions of ≲75 pc that are comparable to the size of giant molecular clouds (GMCs). This is because an aperture of GMC size captures only particular stages of GMC evolution and because Hα traces hot gas around star-forming regions and is displaced from the H{sub 2} peaks themselves. By examining the evolving environment around star clusters, we speculate that the breakdown of the traditional star formation laws of the Kennicutt-Schmidt type at small scales is further aided by a combination of stars drifting from their birthplaces and molecular clouds being dispersed via stellar feedback.

79 ASTRONOMY AND ASTROPHYSICS↗

Gradient-based stochastic estimation of the density matrix

Fast estimation of the single-particle density matrix is key to many applications in quantum chemistry and condensed matter physics. The best numerical methods leverage the fact that the density matrix elements f(H) ij decay rapidly with distance rij between orbitals. This decay is usually exponential. However, for the special case of metals at zero temperature, algebraic decay of the density matrix appears and poses a significant numerical challenge. Here, we introduce a gradient-based probing method to estimate all local density matrix elements at a computational cost that scales linearly with system size. For zero-temperature metals, the stochastic error scales like S -(d+2)/2d , where d is the dimension and S is a prefactor to the computational cost. The convergence becomes exponential if the system is at finite temperature or is insulating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Estimation of Enceladus Plume Density Using Cassini Flight Data

The Cassini spacecraft was launched on October 15, 1997 by a Titan 4B launch vehicle. After an interplanetary cruise of almost seven years, it arrived at Saturn on June 30, 2004. In 2005, Cassini completed three flybys of Enceladus, a small, icy satellite of Saturn. Observations made during these flybys confirmed the existence of water vapor plumes in the south polar region of Enceladus. Five additional low-altitude flybys of Enceladus were successfully executed in 2008-9 to better characterize these watery plumes. During some of these Enceladus flybys, the spacecraft attitude was controlled by a set of three reaction wheels. When the disturbance torque imparted on the spacecraft was predicted to exceed the control authority of the reaction wheels, thrusters were used to control the spacecraft attitude. Using telemetry data of reaction wheel rates or thruster on-times collected from four low-altitude Enceladus flybys (in 2008-10), one can reconstruct the time histories of the Enceladus plume jet density. The 1 sigma uncertainty of the estimated density is 5.9-6.7% (depending on the density estimation methodology employed). These plume density estimates could be used to confirm measurements made by other onboard science instruments and to support the modeling of Enceladus plume jets.

Wang, Eric K.↗

Application of the Redlich-Kister expansion for estimating the density of molten fluoride psuedo-ternary salt systems of nuclear industry interest

The development of next-generation molten salt reactors relies on accurate knowledge of the thermophysical properties of the candidate coolant and fueled molten salts. These thermophysical properties include density, viscosity, thermal conductivity, and heat capacity. Because of difficulties in measuring thermophysical properties of molten salts, there are many gaps in the current state of thermophysical property knowledge of these salts, particularly those that contain actinides or beryllium. Therefore, leveraging modeling techniques to estimate unknown molten salt thermophysical properties and guide future experimental measurements has high value for the nuclear industry. Here, the densities of molten fluoride pseudo-ternary salt systems, which are of interest to the nuclear industry, were estimated using Redlich-Kister expansion and Muggianu interpolation techniques. The pseudo-ternary systems considered for estimation in this study were NaF-LiF-ZrF 4 , LiF-BeF 2 -ZrF 4 , LiF-BeF 2 -ThF 4 , NaF-LiF-BeF 2 , NaF-KF-BeF 2 , NaF-ZrF 4 -UF 4 , and NaF-BeF 2 -UF 4 . This Redlich-Kister estimation approach accounts for nonideal mixing behavior based on pseudo-binary subsystem interaction parameters determined from experimentally measured pseudo-binary system density data sets. The Redlich-Kister estimation was compared with the method of additive molar volumes, which assumes ideal mixing. Additionally, the Redlich-Kister approach was used to determine previously unknown binary and ternary interaction parameters based on experimentally measured density data sets for select pseudo-ternary salt systems. The results of this study show improvement in density estimation using the Redlich-Kister approach for all systems considered compared with estimation by additive molar volumes. Furthermore, this analysis allowed for the estimation of nonideal density behavior in experimentally unstudied ZrF 4 -UF 4 and BeF 2 -UF 4 , as well as the quantification of ternary interaction in NaF-LiF-ZrF 4 , NaF-BeF 2 -UF 4 , and NaF-ZrF 4 -UF 4 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Functional Data Analysis in NTCP Modeling: A New Method to Explore the Radiation Dose-Volume Effects

Purpose/Objective(s): To describe a novel method to explore radiation dose-volume effects. Functional data analysis is used to investigate the information contained in differential dose-volume histograms. The method is applied to the normal tissue complication probability modeling of rectal bleeding (RB) for patients irradiated in the prostatic bed by 3-dimensional conformal radiation therapy. Methods and Materials: Kernel density estimation was used to estimate the individual probability density functions from each of the 141 rectum differential dose-volume histograms. Functional principal component analysis was performed on the estimated probability density functions to explore the variation modes in the dose distribution. The functional principal components were then tested for association with RB using logistic regression adapted to functional covariates (FLR). For comparison, 3 other normal tissue complication probability models were considered: the Lyman-Kutcher-Burman model, logistic model based on standard dosimetric parameters (LM), and logistic model based on multivariate principal component analysis (PCA). Results: The incidence rate of grade ≥2 RB was 14%. V{sub 65Gy} was the most predictive factor for the LM (P=.058). The best fit for the Lyman-Kutcher-Burman model was obtained with n=0.12, m = 0.17, and TD50 = 72.6 Gy. In PCA and FLR, the components that describe the interdependence between the relative volumes exposed at intermediate and high doses were the most correlated to the complication. The FLR parameter function leads to a better understanding of the volume effect by including the treatment specificity in the delivered mechanistic information. For RB grade ≥2, patients with advanced age are significantly at risk (odds ratio, 1.123; 95% confidence interval, 1.03-1.22), and the fits of the LM, PCA, and functional principal component analysis models are significantly improved by including this clinical factor. Conclusion: Functional data analysis provides an attractive method for flexibly estimating the dose-volume effect for normal tissues in external radiation therapy.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

PICS: SIMULATIONS OF STRONG GRAVITATIONAL LENSING IN GALAXY CLUSTERS

Gravitational lensing has become one of the most powerful tools available for investigating the “dark side” of the universe. Cosmological strong gravitational lensing, in particular, probes the properties of the dense cores of dark matter halos over decades in mass and offers the opportunity to study the distant universe at flux levels and spatial resolutions otherwise unavailable. Studies of strongly lensed variable sources offer even further scientific opportunities. One of the challenges in realizing the potential of strong lensing is to understand the statistical context of both the individual systems that receive extensive follow-up study, as well as that of the larger samples of strong lenses that are now emerging from survey efforts. Motivated by these challenges, we have developed an image simulation pipeline, Pipeline for Images of Cosmological Strong lensing (PICS), to generate realistic strong gravitational lensing signals from group- and cluster-scale lenses. PICS uses a low-noise and unbiased density estimator based on (resampled) Delaunay Tessellations to calculate the density field; lensed images are produced by ray-tracing images of actual galaxies from deep Hubble Space Telescope observations. Other galaxies, similarly sampled, are added to fill in the light cone. The pipeline further adds cluster member galaxies and foreground stars into the lensed images. The entire image ensemble is then observed using a realistic point-spread function that includes appropriate detector artifacts for bright stars. Noise is further added, including such non-Gaussian elements as noise window-paning from mosaiced observations, residual bad pixels, and cosmic rays. The aim is to produce simulated images that appear identical—to the eye (expert or otherwise)—to real observations in various imaging surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Mixture density network estimation of continuous variable maximum likelihood using discrete training samples

Abstract Mixture density networks (MDNs) can be used to generate posterior density functions of model parameters $$\varvec{\theta }$$ θ given a set of observables $${\mathbf {x}}$$ x . In some applications, training data are available only for discrete values of a continuous parameter $$\varvec{\theta }$$ θ . In such situations, a number of performance-limiting issues arise which can result in biased estimates. We demonstrate the usage of MDNs for parameter estimation, discuss the origins of the biases, and propose a corrective method for each issue.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Notos - a galaxy tool to analyze CpN observed expected ratios for inferring DNA methylation types

Background: DNA methylation patterns store epigenetic information in the vast majority of eukaryotic species. The relatively high costs and technical challenges associated with the detection of DNA methylation however have created a bias in the number of methylation studies towards model organisms. Consequently, it remains challenging to infer kingdom-wide general rules about the functions and evolutionary conservation of DNA methylation. Methylated cytosine is often found in specific CpN dinucleotides, and the frequency distributions of, for instance, CpG observed/expected (CpG o/e) ratios have been used to infer DNA methylation types based on higher mutability of methylated CpG. Results: Predominantly model-based approaches essentially founded on mixtures of Gaussian distributions are currently used to investigate questions related to the number and position of modes of CpG o/e ratios. These approaches require the selection of an appropriate criterion for determining the best model and will fail if empirical distributions are complex or even merely moderately skewed. We use a kernel density estimation (KDE) based technique for robust and precise characterization of complex CpN o/e distributions without a priori assumptions about the underlying distributions. Conclusions: We show that KDE delivers robust descriptions of CpN o/e distributions. For straightforward processing, we have developed a Galaxy tool, called Notos and available at the ToolShed, that calculates these ratios of input FASTA files and fits a density to their empirical distribution. Based on the estimated density the number and shape of modes of the distribution is determined, providing a rational for the prediction of the number and the types of different methylation classes. Notos is written in R and Perl.

59 BASIC BIOLOGICAL SCIENCES↗

The circumstellar nature of the metallic features in a hot DA white dwarf

A new co-added IUE echelle spectrum of the bright DA white dwarf CD -38 deg 10980, together with a newly determined radial velocity for this star, indicate that the sharp lined Si and C absorption features seen in the UV are clearly circumstellar in origin. Absorption in both excited and ground state transitions occurs at a velocity displaced by -12.1 +/- 2.0 km.s with respect to the photospheric velocity. Weak features due to the Si IV doublet are seen at a velocity intermediate between that of the circumstellar features and the photosphere. First time estimates of column densities for excited and ground states of C II, Si II, and Si III are derived. These quantities are used with electron density estimates derived from these species to determine the location and physical conditions of the circumstellar gas in the vicinity of CD -38 deg 10980. If collisional excitation alone is responsible for the excited levels of Si III observed in CD -38 deg 10980, then electron densities in the circumstellar gas must exceed 10(exp 9)/cu cm. Substantially lower electron densities are possible if the circumstellar gas is located near enough to the star so that photoexcitation is the dominant process responsible for the excited lines seen in the UV. Strong limits are placed on the photospheric abundance of Si and C in the star itself. These limits are in sharp contrast to the theoretical predictions of radiative levitation in which Si, but not C, is expected in the photosphere of a white dwarf such as CD -38 deg 10980. The interstellar line of sight to CD -38 deg 10980 is also investigated.

Holberg, J. B.↗

Two-chord interferometry using 3.39 {mu}m He-Ne laser on a flux-coil-generated FRC

A two-chord {lambda}{sub IR}{approx}3.39 {mu}m He-Ne laser interferometer system was developed for a flux-coil-generated field-reversed configuration to estimate the electron density and the total temperature of the field-reversed configuration (FRC) plasma. This two-chord heterodyne interferometer system consists of a single {approx}2 mW infrared He-Ne laser, a visible ({lambda}{sub vis}{approx}632.8 nm) He-Ne laser for the alignment, a 40 MHz acousto-optic modulator, photodetectors, and quadrature phase detectors. Initial measurement was performed and the measured average electron densities were 2-10x10{sup 19} m{sup -3} at two different radial positions in the midplane. A time shift in density was observed as the FRC expands radially. The time evolution of the line-averaged density agrees with the density estimated from the in situ internal magnetic probes, based on a rigid-rotor profile model.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Estimation of the number of wild pigs found in the United States

Based on a compilation of three estimation approaches, the total nationwide population of wild pigs in the United States numbers approximately 6.3 million animals, with that total estimate ranging from 4.4 up to 11.3 million animals. The majority of these numbers (99 percent), which were encompassed by ten states (i.e., Alabama, Arkansas, California, Florida, Georgia, Louisiana, Mississippi, Oklahoma, South Carolina and Texas), were based on defined estimation methodologies (e.g., density estimates correlated to the total potential suitable wild pig habitat statewide, statewide harvest percentages, statewide agency surveys regarding wild pig distribution and numbers). In contrast to the pre-1990 estimates, none of these more recent efforts, collectively encompassing 99 percent of the total, were based solely on anecdotal information or speculation. To that end, one can defensibly state that the wild pigs found in the United States number in the millions of animals, with the nationwide population estimated to arguably vary from about four million up to about eleven million individuals.

59 BASIC BIOLOGICAL SCIENCES↗