Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional approximation”

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,675 records · Page 93

Response function of modulated grid Faraday cup plasma instruments

Modulated grid Faraday cup plasma analyzers are a very useful tool for making in situ measurements of space plasmas. One of their great attributes is that their simplicity permits their angular response function to be calculated theoretically. An expression is derived for this response function by computing the trajectories of the charged particles inside the cup. The Voyager plasma science experiment is used as a specific example. Two approximations to the rigorous response function useful for data analysis are discussed. Multisensor analysis of solar wind data indicates that the formulas represent the true cup response function for all angles of incidence with a maximum error of only a few percent.

Barnett, A.↗

Communication: Near-locality of exchange and correlation density functionals for 1- and 2-electron systems

The uniform electron gas and the hydrogen atom play fundamental roles in condensed matter physics and quantum chemistry. The former has an infinite number of electrons uniformly distributed over the neutralizing positively charged background, and the latter only one electron bound to the proton. The uniform electron gas was used to derive the local spin density approximation to the exchange-correlation functional that undergirds the development of the Kohn-Sham density functional theory. We show here that the ground-state exchange-correlation energies of the hydrogen atom and many other 1- and 2-electron systems are modeled surprisingly well by a different local spin density approximation (LSDA0). LSDA0 is constructed to satisfy exact constraints but agrees surprisingly well with the exact results for a uniform two-electron density in a finite, curved three-dimensional space. We also apply LSDA0 to excited or noded 1-electron densities, where it works less well. Furthermore, we show that the localization of the exact exchange hole for a 1- or 2-electron ground state can be measured by the ratio of the exact exchange energy to its optimal lower bound.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Approximate formulation of redistribution in the Ly-alpha, Ly-beta, H-alpha system

Simple approximate formulas are given for the coupled redistribution of Ly-alpha, Ly-beta, and H-alpha, by using well-defined approximations to an essentially exact formulation. These formulas incorporate all the essential physics including Raman scattering, lower state radiative decay, and correlated terms representing emission during a collision which must be retained in order that the emission coefficients are properly behaved in the line wings. Approximate expressions for the appropriate line broadening parameters are collected. Finally, practical expressions for the source functions are given. These are formulated through newly introduced nonimpact redistribution functions, which are shown to be reasonably approximated by existing (ordinary and generalized) redistribution functions.

Cooper, J.↗

A simple parameterization for quality factor as a function of linear energy transfer

The paper presents a simple analytic approximation of the radiation quality factor (Q) as a function of linear energy transfer for use in radiation protection calculations. The paper also presents estimated quality factors in water for protons over a broad range of incident energies. It is shown that the quality factors are less than unity for all proton energies greater than 13 MeV.

Townsend, Lawrence W.↗

Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. I. An efficient and simple linear scaling local MP2 method that uses an intermediate basis of pair natural orbitals

In this work, a systematic infrastructure is described that formalizes concepts implicit in previous work and greatly simplifies computer implementation of reduced-scaling electronic structure methods. The key concept is sparse representation of tensors using chains of sparse maps between two index sets. Sparse map representation can be viewed as a generalization of compressed sparse row, a common representation of a sparse matrix, to tensor data. By combining few elementary operations on sparse maps (inversion, chaining, intersection, etc.), complex algorithms can be developed, illustrated here by a linear-scaling transformation of three-center Coulomb integrals based on our compact code library that implements sparse maps and operations on them. The sparsity of the three-center integrals arises from spatial locality of the basis functions and domain density fitting approximation. A novel feature of our approach is the use of differential overlap integrals computed in linear-scaling fashion for screening products of basis functions. Finally, a robust linear scaling domain based local pair natural orbital second-order Möller-Plesset (DLPNO-MP2) method is described based on the sparse map infrastructure that only depends on a minimal number of cutoff parameters that can be systematically tightened to approach 100% of the canonical MP2 correlation energy. With default truncation thresholds, DLPNO-MP2 recovers more than 99.9% of the canonical resolution of the identity MP2 (RI-MP2) energy while still showing a very early crossover with respect to the computational effort. Based on extensive benchmark calculations, relative energies are reproduced with an error of typically <0.2 kcal/mol. The efficiency of the local MP2 (LMP2) method can be drastically improved by carrying out the LMP2 iterations in a basis of pair natural orbitals. While the present work focuses on local electron correlation, it is of much broader applicability to computation with sparse tensors in quantum chemistry and beyond.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extending Parsimonious Bayesian Inference

Parsimonious Bayesian inference is a theoretical framework for efficient data assimilation that seeks to balance increased consistency between predictions and training data against corresponding increases in model complexity. Within this framework, over-training is understood as optimization that encodes excessive information within model parameters while only achieving small improvements between predictions and training data. This project aims to develop practical methods of limiting excess model information during optimization. One key observation is that practical heuristics for parsimonious learning in high-dimensions must balance expressivity, i.e. the ability of the model to capture diverse predictions with only a few non-zero parameters, against discoverability, i.e. the ability to train the model with gradient-based optimization and drive parameters to low information states. As such, we developed logical activation functions that are able to adaptively approximate arbitrary truth tables that define Boolean logic operations within a probabilistic framework. These functions have demonstrated the ability to learn exclusive disjunction (XOR) and conditioned disjunction (if [condition] then [result_if_true] else [result_if_false]) within a single layer of a neural network. To efficiently exploit these activation functions to drive parsimonious learning required several other advances within the domain of variational inference. The most efficient form of complexity suppression is structured sparsification, driving most model parameters to zero while achieving the structural coherence among nonzeros needed for bandwidth reduction. Such models are not only far more efficient at suppressing information-theoretic complexity, they also reduce the other forms of complexity (computations, communication, storage, and the number of dependencies needed to evaluate predictions). Aiming to support enhanced sparsification, this project examined new approaches to high-dimensional variational inference that allow us to calibrate and control parameter uncertainty during optimization. By identifying which parameters can sustain sparsifying perturbations with little impact on prediction quality, we can develop better pruning strategies by framing them as approximate Bayesian inference. These advances also open paths to mitigate concerns with deploying advanced learning methods in resource-constrained environments, such as running models on power-limited or communication-limited devices.

97 MATHEMATICS AND COMPUTING↗

The response function of modulated grid Faraday cup plasma instruments

Modulated grid Faraday cup plasma analyzers are a very useful tool for making in situ measurements of space plasmas. One of their great attributes is that their simplicity permits their angular response function to be calculated theoretically. An expression is derived for this response function by computing the trajectories of the charged particles inside the cup. The Voyager Plasma Science (PLS) experiment is used as a specific example. Two approximations to the rigorous response function useful for data analysis are discussed. The theoretical formulas were tested by multi-sensor analysis of solar wind data. The tests indicate that the formulas represent the true cup response function for all angles of incidence with a maximum error of only a few percent.

Barnett, A.↗

Seasonal Variability of Middle Latitude Ozone in the Lowermost Stratosphere Derived from Probability Distribution Functions

We present a study of the distribution of ozone in the lowermost stratosphere with the goal of characterizing the observed variability. The air in the lowermost stratosphere is divided into two population groups based on Ertel's potential vorticity at 300 hPa. High (low) potential vorticity at 300 hPa indicates that the tropopause is low (high), and the identification of these two groups is made to account for the dynamic variability. Conditional probability distribution functions are used to define the statistics of the ozone distribution from both observations and a three-dimensional model simulation using winds from the Goddard Earth Observing System Data Assimilation System for transport. Ozone data sets include ozonesonde observations from northern midlatitude stations (1991-96) and midlatitude observations made by the Halogen Occultation Experiment (HALOE) on the Upper Atmosphere Research Satellite (UARS) (1994- 1998). The conditional probability distribution functions are calculated at a series of potential temperature surfaces spanning the domain from the midlatitude tropopause to surfaces higher than the mean tropical tropopause (approximately 380K). The probability distribution functions are similar for the two data sources, despite differences in horizontal and vertical resolution and spatial and temporal sampling. Comparisons with the model demonstrate that the model maintains a mix of air in the lowermost stratosphere similar to the observations. The model also simulates a realistic annual cycle. Results show that during summer, much of the observed variability is explained by the height of the tropopause. During the winter and spring, when the tropopause fluctuations are larger, less of the variability is explained by tropopause height. This suggests that more mixing occurs during these seasons. During all seasons, there is a transition zone near the tropopause that contains air characteristic of both the troposphere and the stratosphere. The relevance of the results to the assessment of the environmental impact of aircraft effluence is also discussed.

Rood, Richard B.↗

The Unsteady Lift of a Finite Wing

Unsteady lift function for wings of finite aspect ratio have been calculated by approximate methods involving corrections of the aerodynamic inertia and of the angle of the infinite wing. The starting lift of the finite wing is found to be only slightly less than that of the infinite wing; whereas the final lift may be considerably less. The calculations indicate that the distribution of lift near the start is similar to the final distribution. Both the indicia and the oscillating lift functions are given. Approximate operational equivalents of the functions have been devised to facilitate the calculation of lift under various conditions of motion.

Jones, Robert T↗

Kinetic theory of turbulence for parallel propagation revisited: Low-to-intermediate frequency regime

A previous paper [P. H. Yoon, “Kinetic theory of turbulence for parallel propagation revisited: Formal results,” Phys. Plasmas 22, 082309 (2015)] revisited the second-order nonlinear kinetic theory for turbulence propagating in directions parallel/anti-parallel to the ambient magnetic field, in which the original work according to Yoon and Fang [Phys. Plasmas 15, 122312 (2008)] was refined, following the paper by Gaelzer et al. [Phys. Plasmas 22, 032310 (2015)]. The main finding involved the dimensional correction pertaining to discrete-particle effects in Yoon and Fang's theory. However, the final result was presented in terms of formal linear and nonlinear susceptibility response functions. In the present paper, the formal equations are explicitly written down for the case of low-to-intermediate frequency regime by making use of approximate forms for the response functions. The resulting equations are sufficiently concrete so that they can readily be solved by numerical means or analyzed by theoretical means. The derived set of equations describe nonlinear interactions of quasi-parallel modes whose frequency range covers the Alfvén wave range to ion-cyclotron mode, but is sufficiently lower than the electron cyclotron mode. The application of the present formalism may range from the nonlinear evolution of whistler anisotropy instability in the high-beta regime, and the nonlinear interaction of electrons with whistler-range turbulence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Small-angle approximation to the transfer of narrow laser beams in anisotropic scattering media

The broadening and the signal power detected of a laser beam traversing an anisotropic scattering medium were examined using the small-angle approximation to the radiative transfer equation in which photons suffering large-angle deflections are neglected. To obtain tractable answers, simple Gaussian and non-Gaussian functions for the scattering phase functions are assumed. Two other approximate approaches employed in the field to further simplify the small-angle approximation solutions are described, and the results obtained by one of them are compared with those obtained using small-angle approximation. An exact method for obtaining the contribution of each higher order scattering to the radiance field is examined but no results are presented.

Box, M. A.↗

Direct non-Born-Oppenheimer variational calculations of all bound vibrational states corresponding to the first rotational excitation of D{sub 2} performed with explicitly correlated all-particle Gaussian functions

Direct variational calculations where the Born-Oppenheimer approximation is not assumed are done for all rovibrational states of the D{sub 2} molecule corresponding to first excited rotational level (the N = 1 states). All-particle explicitly correlated Gaussian basis functions are used in the calculations. The exponential parameters of the Gaussians are optimized with the aid of analytically calculated energy gradient determined with respect to these parameters. The results allow to determine the ortho-para spin isomerization energies as a function of the vibrational quantum number.

74 ATOMIC AND MOLECULAR PHYSICS↗

Quantum Effects on the Capacitance of Graphene-Based Electrodes

We recently measured quantum capacitance for electric double layers (EDL) at electrolyte/graphene interfaces. However, the importance of quantum capacitance in realistic carbon electrodes is not clear. Toward understanding that from a theoretical perspective, here we studied the quantum capacitance and total capacitance of graphene electrodes as a function of the number of graphene layers. The quantum capacitance was obtained from electronic density functional theory based on fixed band approximation with an implicit solvation model, while the EDL capacitances were from classical density functional theory. We found that quantum capacitance plays a dominant role in total capacitance of the single-layer graphene both in aqueous and ionic-liquid electrolytes but the contribution decreases as the number of graphene layers increases. Moreover, the total integral capacitance roughly levels off and is dominated by the EDL capacitance beyond about four graphene layers. Finally, because many porous carbons have nanopores with stacked graphene layers at the surface, this research provides a good estimate of the effect of quantum capacitance on their electrochemical performance.

36 MATERIALS SCIENCE↗

Parity violation in nuclear magnetic resonance frequencies of chiral tetrahedral tungsten complexes NWXYZ (X, Y, Z = H, F, Cl, Br or I)

Density functional theory within the two-component quasi-relativistic zeroth-order regular approximation (ZORA) is used to predict parity violation shifts in {sup 183}W nuclear magnetic resonance shielding tensors of chiral, tetrahedrally bonded tungsten complexes of the form NWXYZ (X, Y, Z = H, F, Cl, Br or I), as well as for the heavier systems NWHAtF and NWH(117)F for comparison. The calculations reveal that sub-mHz accuracy is required to detect such tiny effects in this class of compounds, and that parity violation effects are very sensitive to the choice of ligands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electronic structure of antifluorite Cu{sub 2}X (X = S, Se, Te) within the modified Becke-Johnson potential plus an on-site Coulomb U

The traditional photon absorbers Cu{sub 2−x}X (X = S, Se, and Te) have regained significant research attention in the search of earth-abundant photovoltaic materials. These moderate- and narrow-gap materials have also been shown to exhibit excellent thermoelectric properties recently. However, semimetallic band structures with inverted band orderings are predicted for antifluorite structure Cu{sub 2}X using density functional theory with the local density approximation or the generalized gradient approximation. We find that semiconducting band structures and normal band orderings can be obtained using the modified Becke-Johnson potential plus an on-site Coulomb U (the mBJ+U approach), which is consistent with our earlier finding for diamond-like Cu-based multinary semiconductors [Y. Zhang, J. Zhang, W. Gao, T. A. Abtew, Y. Wang, P. Zhang, and W. Zhang, J. Chem. Phys. 139, 184706 (2013)]. The trend of the chemical bonding of Cu{sub 2}X is analyzed, which shows that the positions of the valence band maximum and conduction band minimum are strongly affected by the inter-site pd and intra-site sp hybridizations, respectively. The calculated gaps of Cu{sub 2}S and Cu{sub 2}Se still seem to be underestimated compared with experimental results. We also discuss the effects of different structural phases and Cu disordering and deficiency on the bandgaps of these materials.

36 MATERIALS SCIENCE↗

Platinum atomic wire encapsulated in gold nanotubes: A first principle study

The nanotubes of gold incorporated with platinum atomic wire have been investigated by means of firstprinciples density functional theory with plane wave pseudopotential approximation. The structure with zig-zag chain of Pt atoms in side gold is found to be 0.73 eV lower in energy in comparison to straight chain of platinum atoms. The Fermi level of the composite tube was consisting of d-orbitals of Pt atoms. Further interaction of oxygen with these tubes reveals that while tube with zig-zag Pt prefers dissociative adsorption of oxygen molecule, the gold tube with linear Pt wire favors molecular adsorption.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

First-order nonadiabatic coupling matrix elements between excited states: A Lagrangian formulation at the CIS, RPA, TD-HF, and TD-DFT levels

Analytic expressions for the first-order nonadiabatic coupling matrix elements between electronically excited states are first formulated exactly via both time-independent equation of motion and time-dependent response theory, and are then approximated at the configuration interaction singles, particle-hole/particle-particle random phase approximation, and time-dependent density functional theory/Hartree-Fock levels of theory. Note that, to get the Pulay terms arising from the derivatives of basis functions, the standard response theory designed for electronic perturbations has to be extended to nuclear derivatives. The results are further recast into a Lagrangian form that is similar to that for excited-state energy gradients and allows to use atomic orbital based direct algorithms for large molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗