Search NASA⌕ Search

SEARCH · Search NASA

Results for “Value function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,189 records · Page 66

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗

Including frameworks of public health ethics in computational modelling of infectious disease interventions

Decisions on public health interventions to control infectious diseases are often informed by computational models. Interpreting the predicted outcomes of a public health decision requires not only high-quality modelling but also an ethical framework for assessing the benefits and harms associated with different options. The design and specification of ethical frameworks matured independently of computational modelling, so many values recognized as important for ethical decision-making are missing from computational models. We demonstrate a proof-of-concept approach to incorporate multiple public health values into the evaluation of a simple computational model for vaccination against a pathogen such as SARS-CoV-2. By examining a bounded space of alternative prioritizations of three values relevant to public health ethics (aggregate clinical burden, equity in clinical burden, equity in adverse effects from vaccination), we identify value trade-offs, where the outcomes of optimal strategies differ depending on the ethical framework. This work demonstrates an approach to incorporating diverse values into decision criteria used to evaluate outcomes of models of infectious disease interventions.

"Mathematical Biology"↗

Starting Vortex Identified as Key to Unsteady Ejector Performance

Unsteady ejectors are currently under investigation for use in some pulse-detonation-engine-based propulsion systems. Experimental measurements made in the past, and recently at the NASA Glenn Research Center, have demonstrated that thrust augmentation can be enhanced considerably when the driver is unsteady. In ejector systems, thrust augmentation is defined as = T(sup Total)/T(sup j), where T(sup Total) is the total thrust of the combined ejector and driving jet and T(sup j) is the thrust due to the driving jet alone. There are three images in this figure, one for each of the named thrust sources. The images are color contours of measured instantaneous vorticity. Each image is an ensemble average of at least 150 phase-locked measurements. The flow is from right to left, and the shape and location of each driver is shown on the far right of each image. The emitted vortex is a clearly defined "doughnut" of highly vortical (spinning) flow. In these planar images, the vortex appears as two distorted circles, one above, and one below the axis of symmetry. Because they are spinning in the opposite direction, the two circles have vorticity of opposite sign and thus are different colors. There is also a rectangle shown in each image. Its width represents the ejector diameter that was found experimentally to yield the highest thrust augmentation. It is apparent that the optimal ejector diameter is that which just "captures" the vortex: that is, the diameter bounding the outermost edge of the vortex structure. The exact mechanism behind the enhanced performance is unclear; however, it is believed to be related to the powerful vortex emitted with each pulse of the unsteady driver. As such, particle imaging velocimetry (PIV) measurements were obtained for three unsteady drivers: a pulsejet, a resonance tube, and a speaker-driven jet. All the drivers were tested with ejectors, and all exhibited performance enhancement over similarly sized steady drivers. The characteristic starting vortices of each driver are shown in these images. The images are color contours of measured instantaneous vorticity. Each image is an ensemble average of at least 150 phase-locked measurements. The flow is from right to left. The shape and location of each driver is shown on the far right of each image. The rectangle shown in each image represents the ejector diameter that was found experimentally to yield the highest thrust augmentation. It is apparent that the optimal ejector diameter is that which just "captures" the vortex: that is, the diameter bounding the outermost edge of the vortex structure. Although not shown, it was observed that the emitted vortex spread as it traveled downstream. The spreading rate for the pulsejet is shown as the dashed lines in the top image. A tapered ejector was fabricated that matched this shape. When tested, the ejector demonstrated superior performance to all those previously tested at Glenn (which were essentially of straight, cylindrical form), achieving a remarkable thrust augmentation of 2. The measured thrust augmentation is shown as a function of ejector length. Also shown are the thrust augmentation values achieved with the straight, cylindrical ejectors of varying diameters. Here, thrust augmentation is plotted as a function of ejector length for several families of ejector diameters. It can be seen that large thrust augmentation values are indeed obtained and that they are sensitive to both ejector length and diameter, particularly the latter. Five curves are shown. Four correspond to straight ejector diameters of 2.2, 3.0, 4.0, and 6.0 in. The fifth curve corresponds to the tapered ejector contoured to bound the emitted vortex. For each curve, there are several data points corresponding to different lengths. The largest value of thrust augmentation is 2.0 for the tapered ejector and 1.81 for the straight ejectors. Regardless of their diameters, all the ejectors trend toward peak performance at a particular leng. That the cross-sectional dimensions of optimal ejectors scaled precisely with the vortex dimensions on three separate pulsed thrust sources demonstrates that the action of the vortex is responsible for the enhanced ejector performance. The result also suggests that, in the absence of a complete understanding of the entrainment and augmentation mechanisms, methods of characterizing starting vortices may be useful for correlating and predicting unsteady ejector performance.

Paxson, Daniel E.↗

A light weight regularization for wave function parameter gradients in quantum Monte Carlo

The parameter derivative of the expectation value of the energy, ∂E/∂p, is a key ingredient in variational Monte Carlo (VMC) wave function optimization methods. In some cases, a naïve estimate of this derivative suffers from an infinite variance, which inhibits the efficiency of optimization methods that rely on a stable estimate of the derivative. In this work, we derive a simple regularization of the naïve estimator, which is trivial to implement in existing VMC codes, has finite variance, and a negligible bias, which can be extrapolated to zero bias with no extra cost. We use this estimator to construct an unbiased, finite variance estimation of ∂E/∂p for a multi-Slater–Jastrow trial wave function on the LiH molecule and in the optimization of a multi-Slater–Jastrow trial wave function on the CuO molecule. This regularized estimator is a simple and efficient estimator of ∂E/∂p for VMC optimization techniques.

36 MATERIALS SCIENCE↗

Assessing the accuracy of compound formation energies with quantum Monte Carlo

Accurately predicting the formation energy of a compound, which describes its thermodynamic stability, is a key challenge in materials physics. Here, we employ a many-body quantum Monte Carlo (QMC) with single-reference trial functions to compute the formation energy of two electronically disparate compounds, the intermetallic VPt2 and the semiconductor CuI, for which standard density functional theory (DFT) predictions using both the Perdew-Burke-Ernzerhof (PBE) and the strongly constrained and appropriately normed (SCAN) density functional approximations deviate markedly from available experimental values. For VPt2, we find an agreement between QMC, SCAN, and PBE0 estimates, which therefore remain in disagreement with the much less exothermic experimental value. For CuI, the QMC result agrees with neither SCAN nor PBE pointing towards DFT exchange-correlation biases, likely related to the localized Cu 3d electrons. Compared to the behavior of some density functional approximations within DFT, spin-averaged QMC exhibits a smaller but still appreciable deviation when compared to experiment. The QMC result is slightly improved by incorporating spin-orbit corrections for CuI and solid I 2 , so that experiment and theory are brought into imperfect but reasonable agreement within about 120 meV/atom.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Oscillation characteristics of neutrino in the model with three sterile neutrinos for analysis of the anomalies on small distances

In the framework of the model with three sterile neutrinos, the transition probabilities for different flavours of neutrino are calculated and the graphical dependences are obtained, in particular, for the appearance probability of electron neutrino and antineutrino in the muon neutrino and antineutrino jets as a function of distance and other model parameters at their acceptable values and at the neutrino energy less than 50 MeV, as well as a function of a ratio of distance to the neutrino energy. The theoretical results obtained can be used for analysis of the neutrino data related to the anomalies on small distances.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

External calibration of SIR-B Imagery with area-extended and point targets

Data collected by the Shuttle Imaging Radar-B (SIR-B) on two ascending orbits over a site in west-central Illinois were used to establish end-to-end transfer functions for conversion of the digital numbers on the eight-bit image to values of the radar backscattering coefficient, sigma(0). The transfer function for each data take was defined by the SIR response to an array of six calibrated point targets of known radar cross section and to a large number of area-extended targets, also with a known cross section, as measured by externally calibrated, truck-mounted scatterometers. In spite of variance in the SIR transmit power, end-to-end transfer functions generated from either point targets or area-extended targets were found to yield close agreement for a given pass, indicating that targets of known radar cross section can be used to provide transfer functions for orbital SAR. Although the use of area-extended targets yielded the most satisfactory results, the use of point targets is far simpler and may be the preferred alternative for many situations.

Dobson, M. C.↗

Experimental determination of the valence band offsets of ZnGeN 2 and (ZnGe) 0.94 Ga 0.12 N 2 with GaN

A predicted type-II staggered band alignment with an approximately 1.4 eV valence band offset at the ZnGeN 2 /GaN heterointerface has inspired novel band-engineered III-N/ZnGeN 2 heterostructure-based device designs for applications in high performance optoelectronics. We report on the determination of the valence band offset between metalorganic chemical vapor deposition grown (ZnGe) 1- x Ga 2 x N 2 , for x = 0 and 0.06, and GaN using x-ray photoemission spectroscopy. The valence band of ZnGeN 2 was found to lie 1.45–1.65 eV above that of GaN. This result agrees well with the value predicted by first-principles density functional theory calculations using the local density approximation for the potential profile and quasiparticle self-consistent GW calculations of the band edge states relative to the potential. For (ZnGe) 0.94 Ga 0.12 N 2 the value was determined to be 1.29 eV, ~10%–20% lower than that of ZnGeN 2 . The experimental determination of the large band offset between ZnGeN 2 and GaN provides promising alternative solutions to address challenges faced with pure III-nitride-based structures and devices.

Physics↗

Turbulent boundary layer separation over a rearward facing ramp and its control through mechanical excitation

A vane which oscillates about a fixed point at the inlet to a two-dimensional, 20-deg rearward-facing ramp has been noted to effectively delay the separation of a turbulent boundary layer. A parametric study has been conducted to determine the variation of the maximum values of the static pressure recovery, measured close to the ramp exit as a function of the oscillation frequency of the vane. The resulting curves show a peak value in the maximum static pressure recovery as a function of the vane's oscillation frequency. For the excited cases, corona anemometer measurements in the inner wall regions of the boundary layer along the surface of the ramp show a large range of unsteadiness in the local velocities. Finally, a qualitative explanation of the phenomena of delay of separation is discussed which appears to be supported by the surface static pressure and flow field measurements obtained.

Mckinzie, Daniel J., Jr.↗

Numerical simulations of turbulent trapping in the weak beam-plasma instability

Numerical simulations of the weak beam-plasma instability were done in the turbulent regime where small-scale trapping is a dominant feature of the instability, a regime with behavior not predicted by quasi-linear theory. The results of the simulations were compared with those of a specific model of the turbulence, the so-called 'turbulent trapping' model, which gives precise formulas for the particle correlation functions, and predicts a growth rate well enhanced over the quasi-linear value. It was found that the model gives accurate predictions for the correlation functions. On the other hand, while growth rates were enhanced over the quasi-linear values, the enhancements observed are smaller than expected from the quantitative predictions of the model.

Theilhaber, K.↗

A simple method for obtaining a three-dimensional proton distribution function from Voyager plasma data

The main sensor of the Vogager plasma experiment consists of a cluster of three, modulated-grid Faraday cups whose normals are arranged symmetrically about the symmetry axis of the cluster at an angle of 20 degrees to that axis. In interplanetary space, each cup explores the positive ion distribution by accepting particles from contiguous slices in velocity space. The slices are narrow in the direction of the normal to the modulating grid but are broad in planes parallel to that grid. The resulting three sets of measurements can be combined to yield the three-dimensional distribution function in the following way: the distribution function is assumed to be gyrotropic. For each value of speed in a frame of reference moving with the bulk velocity of the solar wind, the variation of the distribution function with angle from the field direction is represented by a series of Legendre polynomials. Effects such as double-streaming and heat flow can be well represented by using only the first three terms of the series which are fully specified by the measurements. Examples of the use of this method in the analysis of Voyager data are shown.

Olbert, S.↗

Frequency functions of Venus nightside ion densities

The variability of the total and major ion densities in the nightside Venus ionosphere are studied by deriving an experimental frequency function for these quantities. The frequency function N(i) derived for the total ion density is found to be insensitive to division of the parent set into two approximately equal subsets representing dawn and dusk hemispheres, an upper and lower altitude interval, or two solar zenith angle intervals. The N(i) function is reasonably well peaked near the median value, and 78 percent of the normalized values fall within a factor of 2.8 of the median value. The frequency of occurrence decreases more rapidly on the high side of the median than on the low side. The frequency functions for the constituent ions are similar to that for N(i).

Knudsen, William C.↗

A Search for Correlations between Turbulence and Star Formation in LITTLE THINGS Dwarf Irregular Galaxies

Turbulence has the potential for creating gas density enhancements that initiate cloud and star formation (SF), and it can be generated locally by SF. To study the connection between turbulence and SF, we looked for relationships between SF traced by FUV images, and gas turbulence traced by kinetic energy density (KED) and velocity dispersion (v {sub disp}) in the LITTLE THINGS sample of nearby dIrr galaxies. We performed 2D cross-correlations between FUV and KED images, measured cross-correlations in annuli to produce correlation coefficients as a function of radius, and determined the cumulative distribution function of the cross-correlation value. We also plotted on a pixel-by-pixel basis the locally excess KED, v {sub disp}, and H i mass surface density, Σ{sub HI}, as determined from the respective values with the radial profiles subtracted, versus the excess SF rate density Σ{sub SFR}, for all regions with positive excess Σ{sub SFR}. We found that Σ{sub SFR} and KED are poorly correlated. The excess KED associated with SF implies a ∼0.5% efficiency for supernova energy to pump local H i turbulence on the scale of the resolution here, which is a factor of ∼2 too small for all of the turbulence on a galactic scale. The excess v {sub disp} in SF regions is also small, only ∼0.37 km s{sup −1}. The local excess in Σ{sub HI} corresponding to an excess in Σ{sub SFR} is consistent with a H i consumption time of ∼1.6 Gyr in the inner parts of the galaxies. The similarity between this timescale and the consumption time for CO implies that CO-dark molecular gas has comparable mass to H i in the inner disks.

74 ATOMIC AND MOLECULAR PHYSICS↗

Systematic characterization of gene function in the photosynthetic alga Chlamydomonas reinhardtii

Most genes in photosynthetic organisms remain functionally uncharacterized. Here, using a barcoded mutant library of the model eukaryotic alga Chlamydomonas reinhardtii, we determined the phenotypes of more than 58,000 mutants under more than 121 different environmental growth conditions and chemical treatments. A total of 59% of genes are represented by at least one mutant that showed a phenotype, providing clues to the functions of thousands of genes. Mutant phenotypic profiles place uncharacterized genes into functional pathways such as DNA repair, photosynthesis, the CO 2 -concentrating mechanism and ciliogenesis. We illustrate the value of this resource by validating phenotypes and gene functions, including three new components of an actin cytoskeleton defense pathway. The data also inform phenotype discovery in land plants; mutants in Arabidopsis thaliana genes exhibit phenotypes similar to those we observed in their Chlamydomonas homologs. We anticipate that this resource will guide the functional characterization of genes across the tree of life.

59 BASIC BIOLOGICAL SCIENCES↗

Evaluation of tests for vestibular function

Vestibular testing methodologies are examined in terms of their clinical value and applicability to the assessment of astronaut candidates and human responses to microgravity. Classical and novel testing approaches are examined in terms of nine functions such as otolith organs, vestibular neural processing, and ocular motor systems. The tests reviewed include the basic vestibular-function test battery, postural stability testing, visual-vestibular interaction, reactions to off-axis rotation, and vestibular evoked responses. The basic test battery is shown to be essential for astronaut and inertial-environment baseline-data evaluations, whereas the newer and research tests are only recommended in special cases. Stimulus rearrangement and motion-sickness tests are indicated for testing adaptation to inertial environments. Postural stability and visual-vestibular interaction are recommended for evaluating one's suitability for being an astronaut.

Source record↗

(U) A Linear Response Model Predicts Reactivity From a Density Profile

We tested the ability to predict the system reactivity, described by alpha, given a density profile using a simple linear system response. We generated a suite of 1-dimensional density profiles that consisted of nominal density, a discontinuity, and a decay. These profiles were prescribed a functional form and the mass was conserved in all cases. From these density profiles, we calculated the alpha value of the 3-dimensional system.We calculated a linear response function given a training set of the 1-dimensional density profiles, and the system reactivity described by alpha. We tested the robustness of the response function using the remaining test data. Our results showed very good agreement between the predicted and calculated test values, where the distribution of alpha differences was centered about zero and had a standard deviation of 0.005 gens/shake. The predicted and calculated alpha values did not significantly differ (t=-0.0009 p<0.99). We used Singular Value Decomposition (SVD) to reduce the matrix rank by retaining95% of the cumulative singular value contributions. This reduced the matrix rank by 91.7%. We generated the linear response matrix and calculated the difference between the predicted and calculated alpha values. Using the reduced order matrix, we showed good agreement between the predicted and calculated alpha values where the distribution of differences was centered near zero, the standard deviation was 0.006 gens/shake, and the statistical t-test showed good agreement (t=0.02, p<0.98). These results show a linear relationship between a series of 1-dimensional density profiles,where the mass was conserved, and the system reactivity. The next steps of this work will be to investigate the linear response using 2-dimensional density profiles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Position error propagation in the simplex strapdown navigation system

The results of an analysis of the effects of deterministic error sources on position error in the simplex strapdown navigation system were documented. Improving the long term accuracy of the system was addressed in two phases: understanding and controlling the error within the system, and defining methods of damping the net system error through the use of an external reference velocity or position. Review of the flight and ground data revealed error containing the Schuler frequency as well as non-repeatable trends. The only unbounded terms are those involving gyro bias and azimuth error coupled with velocity. All forms of Schuler-periodic position error were found to be sufficiently large to require update or damping capability unless the source coefficients can be limited to values less than those used in this analysis for misalignment and gyro and accelerometer bias. The first-order effects of the deterministic error sources were determined with a simple error propagator which provided plots of error time functions in response to various source error values.

Source record↗