Search NASA⌕ Search

SEARCH · Search NASA

Results for “Basis 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,009 records · Page 56

Half-metallicity and stability of the rock salt BaC and SrC (111) surfaces: A density functional study

The electronic structure and magnetic properties of relaxed (111) surfaces of the alkaline-earth monocarbides BaC and SrC in the stable rock salt structure, are calculated on the basis of first principle density functional theory within the framework of self-consistent field plane wave pseudo-potential method, using the generalized gradient approximation for the exchange-correlation functional. The results of this study reveal that the C-terminated (111) surfaces retain the bulk half-metallic property in both BaC and SrC. The half-metallicity of the C-terminated BaC surface is found to be more robust compared to the bulk BaC due to the larger half-metallic energy gap. In contrast, the half-metallic energy gap of the C-terminated SrC surface is found to be smaller than that of the bulk. The Ba-terminated surface of BaC and the Sr-terminated surface of SrC, however, lose their bulk half-metallicity due to the formation of surface states in the majority spin band gap. The calculations also show that the atomic magnetic moments at the half-metallic C-terminated surfaces in both BaC and SrC increase considerably with respect to the corresponding bulk values, which is explained in terms of an increase in the number of unpaired 2p electrons of the carbon atom at the surface. We also discuss the stability of the surfaces via the calculated bulk formation energies. The bulk formation energies for both BaC and SrC in the rock salt structure are found to be positive, which indicate that the surfaces are not stable at normal pressure and temperature conditions, and non-equilibrium growth techniques may be required for the realization of BaC and SrC thin films.

74 ATOMIC AND MOLECULAR PHYSICS↗

Structure-Function Analysis of Friedreich's Ataxia Mutants Reveals Determinants of Frataxin Binding and Activation of the Fe-S Assembly Complex

Friedreich's ataxia (FRDA) is a progressive neurodegenerative disease associated with the loss of function of the protein frataxin (FXN) that results from low FXN levels due to a GAA triplet repeat expansion or, occasionally, from missense mutations in the FXN gene. Here biochemical and structural properties of FXN variants, including three FRDA missense mutations (N146K, Q148R, and R165C) and three related mutants (N146A, Q148G, and Q153A), were determined in an effort to understand the structural basis for the loss of function. In vitro assays revealed that although the three FRDA missense mutations exhibited similar losses of cysteine desulfurase and Fe-S cluster assembly activities, the causes for these activation defects were distinct. The R165C variant exhibited a k cat /K M higher than that of native FXN but weak binding to the NFS1, ISD11, and ISCU2 (SDU) complex, whereas the Q148R variant exhibited the lowest k cat /K M of the six tested FXN variants and only a modest binding deficiency. The order of the FXN binding affinities for the SDU Fe-S assembly complex was as follows: FXN > Q148R > N146A > Q148G > N146K > Q153A > R165C. Four different classes of FXN variants were identified on the basis of their biochemical properties. Together, these structure-function studies reveal determinants for the binding and allosteric activation of the Fe-S assembly complex and provide insight into how FRDA missense mutations are functionally compromised.

59 BASIC BIOLOGICAL SCIENCES↗

Approaching the complete-basis limit with a truncated many-body expansion

High-accuracy electronic structure calculations with correlated wave functions demand the use of large basis sets and complete-basis extrapolation, but the accuracy of fragment-based quantum chemistry methods has most often been evaluated using double-ζ basis sets, with errors evaluated relative to a supersystem calculation using the same basis set. Here, we examine the convergence towards the basis-set limit of two- and three-body expansions of the energy, for water clusters and ion–water clusters, focusing on calculations at the level of second-order Møller-Plesset perturbation theory (MP2). Several different corrections for basis-set superposition error (BSSE), each consistent with a truncated many-body expansion, are examined as well. We present a careful analysis of how the interplay of errors (from all sources) influences the accuracy of the results. We conclude that fragment-based methods often benefit from error cancellation wherein BSSE offsets both incompleteness of the basis set as well as higher-order many-body effects that are neglected in a truncated many-body expansion. An n-body counterpoise correction facilitates smooth extrapolation to the MP2 basis-set limit, and at n = 3 affords accurate results while requiring calculations in subsystems no larger than trimers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

EXPERIMENTAL AND THEORETICAL STUDY ON THE INFRARED SPECTROSCOPY OF ASTROPHYSICALLY RELEVANT POLYCYCLIC AROMATIC HYDROCARBON DERIVATIVES 2- AND 9-VINYLANTHRACENE

We propose to evaluate the contribution of polycyclic aromatic hydrocarbon (PAH) molecules that carry side groups to the mid-infrared emission spectra. Within this framework, the IR absorption spectra of 2-vinylanthracene (2-VA) and 9-vinylanthracene (9-VA) were measured in Ar matrices at 12 K and in CsI and polyethylene pellets at room temperature. The laboratory spectra were analyzed with the support of simulations based on the density functional theory. For each PAH molecule, eight IR spectra were computed by combining the B3LYP functional with as many different basis sets, namely, 4-31G, 4-31G(d), 6-31G, 6-311G, 6-31G(d), 6-31G(d,p), 6-31+G(d,p), and 6-31++G(d,p). The comparison of the theoretical spectra with the laboratory data allowed us to determine the most suitable combinations for modeling the IR spectra of neutral PAH molecules that carry a vinyl side group. It was concluded from the examples of 2- and 9-VA that the optimum basis set is 6-31G unless a steric interaction has to be taken into account, in which case the optimum basis set is 6-31G(d). Thus, in the presence of such an interaction, the use of d-type polarization functions is recommended. We discuss the possibility for neutral vinyl-substituted PAHs to contribute to the mid-infrared emission spectra and find that their specific features do not match with the mid-infrared aromatic emission bands.

79 ASTRONOMY AND ASTROPHYSICS↗

Computing Legacy Software Behavior to Understand Functionality and Security Properties: An IBM/370 Demonstration

Organizations maintaining mainframe legacy software can benefit from code modernization and incorporation of security capabilities to address the current threat environment. Oak Ridge National Laboratory is developing the Hyperion system to compute the behavior of software as a means to gain understanding of software functionality and security properties. Computation of functionality is critical to revealing security attributes, which are in fact specialized functional behaviors of software. Oak Ridge is collaborating with MITRE Corporation to conduct a demonstration project to compute behavior of legacy IBM Assembly Language code for a federal agency. The ultimate goal is to understand functionality and security vulnerabilities as a basis for code modernization. This paper reports on the first phase, to define functional semantics for IBM Assembly instructions and conduct behavior computation experiments.

reverse engineering↗

ROVIBRATIONAL QUENCHING RATE COEFFICIENTS OF HD IN COLLISIONS WITH He

Along with H{sub 2}, HD has been found to play an important role in the cooling of the primordial gas for the formation of the first stars and galaxies. It has also been observed in a variety of cool molecular astrophysical environments. The rate of cooling by HD molecules requires knowledge of collisional rate coefficients with the primary impactors, H, He, and H{sub 2}. To improve knowledge of the collisional properties of HD, we present rate coefficients for the He-HD collision system over a range of collision energies from 10{sup -5} to 5 Multiplication-Sign 10{sup 3} cm{sup -1}. Fully quantum mechanical scattering calculations were performed for initial HD rovibrational states of j = 0 and 1 for v = 0-17 which utilized accurate diatom rovibrational wave functions. Rate coefficients of all {Delta}v = 0, -1, and -2 transitions are reported. Significant discrepancies with previous calculations, which adopted a small basis and harmonic HD wave functions for excited vibrational levels, were found for the highest previously considered vibrational state of v = 3. Applications of the He-HD rate coefficients in various astrophysical environments are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

How Accurate Are the Minnesota Density Functionals for Noncovalent Interactions, Isomerization Energies, Thermochemistry, and Barrier Heights Involving Molecules Composed of Main-Group Elements?

The 14 Minnesota density functionals published between the years 2005 and early 2016 are benchmarked on a comprehensive database of 4986 data points (84 data sets) involving molecules composed of main-group elements. The database includes noncovalent interactions, isomerization energies, thermochemistry, and barrier heights, as well as equilibrium bond lengths and equilibrium binding energies of noncovalent dimers. Additionally, the sensitivity of the Minnesota density functionals to the choice of basis set and integration grid is explored for both noncovalent interactions and thermochemistry. By and large, the main strength of the hybrid Minnesota density functionals is that the best ones provide very good performance for thermochemistry (e.g., M06-2X), barrier heights (e.g., M08-HX, M08-SO, MN15), and systems heavily characterized by self-interaction error (e.g., M06-2X, M08-HX, M08-SO, MN15), while the main weakness is that none of them are state-of-the-art for the full spectrum of noncovalent interactions and isomerization energies (although M06-2X is recommended from the 10 hybrid Minnesota functionals). Similarly, the main strength of the local Minnesota density functionals is that the best ones provide very good performance for thermochemistry (e.g., MN15-L), barrier heights (e.g., MN12-L), and systems heavily characterized by self-interaction error (e.g., MN12-L and MN15-L), while the main weakness is that none of them are state-of-the-art for the full spectrum of noncovalent interactions and isomerization energies (although M06-L is clearly the best from the four local Minnesota functionals). Finally, as an overall guide, M06-2X and MN15 are perhaps the most broadly useful hybrid Minnesota functionals, while M06-L and MN15-L are perhaps the most broadly useful local Minnesota functionals, although each has different strengths and weaknesses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Massless and massive photons within basis light-front quantization

We compute the structure function, transverse momentum dependent parton distributions (TMDs), and generalized parton distributions (GPDs) for the photon from the light-front quantum electrodynamics (QED) Hamiltonian, determined for its constituent bare photon and electron-positron Fock components within the basis light-front quantization framework. While the structure function and GPDs align well with expectations, some oscillatory behavior in the TMDs emerges as a finite-basis artifact. We also compute the GPD and the corresponding structure function for the longitudinal polarization mode of the massive photon. We further study the photon mass dependence of the calculated observables. After performing nonperturbative renormalization, we obtain a good quality description of the perturbative QED properties of the photon, as well as good agreement with the experimental data for the massless photon’s structure function.

Ab initio calculations↗

Draft Methodology for Cybersecurity Analysis for Adoption of Wireless Technology in Nuclear Power Plants

The report is delivered as DOE-NE Cybersecurity Program Milestone M2CT-22IN1104014-Industry cybersecurity guidance adoption status. The milestone is specifically to document the process and the status of the efforts towards developing guidance in support of industry’s adoption of wireless communication technologies. To transform the operation of nuclear power plants, including domestic light-water reactors, advanced reactors, microreactors, and fission battery, for grid and non-grid applications, secure and resilient wireless capabilities are required. In addition, to realize new operational concepts such as autonomous operation and remote monitoring, advanced sensor, and instrumentation, wireless technology is essential. However, industry implementation is low, and there is no technical basis for how to understand and address potential risks for wireless communications for critical plant functions. A methodology with a technical basis for implementing secure wireless communication is crucial. This wireless adoption methodology is intended to assist a licensee in identifying an appropriate technological approach for securing wireless communications in nuclear power plant. However, this methodology does not provide guidance for addressing wireless design criteria or for addressing all the cybersecurity commitments in licensees’ cybersecurity plans (CSPs) including addressing security impact analysis that a licensee must perform before an update to a Critical Digital Asset (CDA). This methodology guides a licensee through a process of evaluating a proposed wireless implementation and drafting a plant change notification for an example use case. The example is a generic/nonproprietary version of the first test case of the methodology, which, if successfully implemented, is part of the consideration for protecting other functions’ use of wireless in a cybersecure manner.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Introducing GPU Acceleration into the Python-Based Simulations of Chemistry Framework

We introduce the first version of GPU4P Y SCF, a module that provides GPU acceleration of methods in P Y SCF. As a core functionality, this provides a GPU implementation of two-electron repulsion integrals (ERIs) for contracted basis sets comprising up to g functions using the Rys quadrature. As an illustration of how this can accelerate a quantum chemistry workflow, we describe how to use the ERIs efficiently in the integral-direct Hartree–Fock build and nuclear gradient construction. Benchmark calculations show a significant speedup of 2 orders of magnitude with respect to the multithreaded CPU Hartree–Fock code of P Y SCF and the performance comparable to other open-source GPU-accelerated quantum chemical packages, including GAMESS and QUICK, on a single NVIDIA A100 GPU.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Fundamental Vibrational Analyses of the HCN Monomer, Dimer and Associated Isotopologues

In this work we provide high level ab initio treatments of the structures, vibrational frequencies, and electronic energies of the HCN monomer and dimer systems along with several isotopologues. Additionally, the plethora of information related to this system within the literature is summarized and serves as a basis for comparison with the results of this paper. The geometry of the dimer and monomer are reported at the all electroncoupled-cluster singles, doubles, and perturbative triples level of theory [AE-CCSD(T)] with the correlation consistent quadruple-zeta quality basis sets with extra core functions (cc-pCVQZ) from Dunning. The theoretical geometries and electronic structures are further analyzed through the use of the Natural Bond Orbital (NBO) method and Natural Resonance Theory (NRT). At the AE-CCSD(T)/cc-pCVQZ level of theory, the full cubic with semi-diagonal quartic force field for nine dimer and four monomer isotopologues (the parent isotopologue along with 15N, 13C, and D derivatives) were obtained to treat the anharmonicity of the vibrations via second order vibrational perturbation theory (VPT2). Lastly, the enthalpy change associated with the formation of the dimer from two monomer units was determined using the focal point analysis. Computations including coupled-cluster through perturbative quadruples as well as basis sets up to six-zeta quality, including core functions (cc-pCVXZ, X=D,T,Q,5,6) were used to extrapolate to the AE-CCSDT(Q)/CBS energy associated with this hydrogen-bond forming process. After appending anharmonic zero-point vibrational, relativistic, and diagonal Born-Oppenheimer corrections, we report a value of -3.93 kcal mol -1 for the enthalpy of formation. To our knowledge, each set of results (geometries, vibrational frequencies, and energetics) reported in this study represents the highest-level and most reliable theoretical predictions reported for this system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab Initio Thermodynamics of Iridium Surface Oxidation and Oxygen Evolution Reaction

Iridium-based materials are considered as state-of-the-art electrocatalysts for oxygen evolution reaction (OER), however, their stability and catalytic activity greatly depend on surface-state changes induced by electrochemical cycling. To better understand the behavior of the low-index Ir surfaces in an electrochemical environment, we perform a systematic thermodynamic analysis by means of the density functional theory (DFT) calculations. On the basis of computed surface energies of the Ir (111), (110) and (100) facets as a function of applied electrode potential and coverage of adsorbed water species we determine stability maps and predict equilibrium shapes of Ir nanoparticles. Our calculations also show that metastable oxide precursors formed at the initial stages of Ir surface oxidation are responsible for enhanced catalytic activity toward OER as compared to metal surfaces covered by oxygen adsorbates and thick-oxide films. Such enhancement occurs not only due to the modified thermodynamic stability of OER intermediates, but also because thin-oxide layers may display the more energetically favorable I2M (interaction of two M–O units) rather than WNA (water nucleophilic attack) OER mechanism.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Properties of a potential energy matrix in oscillator basis

Matrix elements of potential energy are examined in detail. We consider a model problem — a particle in a central potential. The most popular forms of central potential are taken up, namely, square-well potential, Gaussian, Yukawa and exponential potentials. We study eigenvalues and eigenfunctions of the potential energy matrix constructed with oscillator functions. It is demonstrated that eigenvalues coincide with the potential energy in coordinate space at some specific discrete points. We establish approximate values for these points. It is also shown that the eigenfunctions of the potential energy matrix are the expansion coefficients of the spherical Bessel functions in a harmonic oscillator basis. We also demonstrate a close relation between the separable approximation and L{sup 2} basis (J-matrix) method for the quantum theory of scattering.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A NEW TOOL FOR IMAGE ANALYSIS BASED ON CHEBYSHEV RATIONAL FUNCTIONS: CHEF FUNCTIONS

We introduce a new approach to the modeling of the light distribution of galaxies, an orthonormal polar basis formed by a combination of Chebyshev rational functions and Fourier polynomials that we call CHEF functions, or CHEFs. We have developed an orthonormalization process to apply this basis to pixelized images, and implemented the method as a Python pipeline. The new basis displays remarkable flexibility, being able to accurately fit all kinds of galaxy shapes, including irregulars, spirals, ellipticals, highly compact, and highly elongated galaxies. It does this while using fewer components than similar methods, as shapelets, and without producing artifacts, due to the efficiency of the rational Chebyshev polynomials to fit quickly decaying functions like galaxy profiles. The method is linear and very stable, and therefore is capable of processing large numbers of galaxies in a fast and automated way. Due to the high quality of the fits in the central parts of the galaxies, and the efficiency of the CHEF basis modeling galaxy profiles up to very large distances, the method provides highly accurate estimates of total galaxy fluxes and ellipticities. Future papers will explore in more detail the application of the method to perform multiband photometry, morphological classification, and weak shear measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

First-Principles Modeling of Non-Covalent Interactions in Supramolecular Systems: The Role of Many-Body Effects

Supramolecular host–guest systems play an important role for a wide range of applications in chemistry and biology. The prediction of the stability of host–guest complexes represents a great challenge to first-principles calculations due to an interplay of a wide variety of covalent and noncovalent interactions in these systems. In particular, van der Waals (vdW) dispersion interactions frequently play a prominent role in determining the structure, stability, and function of supramolecular systems. On the basis of the widely used benchmark case of the buckyball catcher complex (C60@C60H28), we assess the feasibility of computing the binding energy of supramolecular host–guest complexes from first principles. Large-scale diffusion Monte Carlo (DMC) calculations are carried out to accurately determine the binding energy for the C60@C60H28 complex (26 ± 2 kcal/mol). On the basis of the DMC reference, we assess the accuracy of widely used and efficient density-functional theory (DFT) methods with dispersion interactions. The inclusion of vdW dispersion interactions in DFT leads to a large stabilization of the C60@C60H28 complex. However, DFT methods including pairwise vdW interactions overestimate the stability of this complex by 9–17 kcal/mol compared to the DMC reference and the extrapolated experimental data. A significant part of this overestimation (9 kcal/mol) stems from the lack of dynamical dielectric screening effects in the description of the molecular polarizability in pairwise dispersion energy approaches. The remaining overstabilization arises from the isotropic treatment of atomic polarizability tensors and the lack of many-body dispersion interactions. A further assessment of a different supramolecular system – glycine anhydride interacting with an amide macrocycle – demonstrates that both the dynamical screening and the many-body dispersion energy are complex contributions that are very sensitive to the underlying molecular geometry and type of bonding. We discuss the required improvements in theoretical methods for achieving “chemical accuracy” in the first-principles modeling of supramolecular systems.

Chemistry↗

Tunable-fidelity wave functions for the ab initio description of scattering and reactions

Here, the no-core shell model (NCSM) is an ab initio method that solves the nuclear many-body problem by expanding the many-particle wave function into a (typically) harmonic oscillator basis and minimizing the energy to obtain the expansion coefficients. Extensions of the NCSM, such as its coupling with microscopic-cluster basis states, further allow for an ab initio treatment of light-ion nuclear reactions of interest for both astrophysics and nuclear technology applications. A downside of the method is the exponential scaling of the basis size with increasing number of nucleons and excitation quanta, which limits its applicability to mass A ≲ 16 nuclei, except for variants where the basis is further down-selected via some truncation scheme. We consider a basis selection method for the NCSM that was first introduced in the context of the large-scale shell model and captures the essential degrees of freedom of the nuclear wave function leading to a favorable complexity scaling for calculations and enabling ab initio reaction calculations in sd-shell nuclei. The particle configurations within the NCSM basis are ordered based on their contribution to the first moment of the Hamiltonian matrix that results from the projection onto the many-body basis. The truncation scheme then consists in retaining only the lowest-first-moment configurations, which typically contain only few many-body basis states (Slater determinants). As the energy threshold above which configurations are disregarded is increased, the size of the basis becomes an almost-continuous variable, allowing for tunable fidelity in the obtained wave functions. The resulting wave functions can then be used directly in ab initio reaction calculations. We present calculations for 7 Li and n + 12 C scattering using nucleon-nucleon interactions derived from chiral effective field theory and softened using the similarity renormalization group method. The obtained energy levels invariably demonstrate exponential convergence with the size of the basis, and we find improved convergence in scattering calculations. To demonstrate the possibilities enabled by the approach, we also present a first calculation for the scattering of neutrons from 24 Mg. The method presented in this work appears promising for future studies of nuclei with mass A > 16, opening multiple future research directions impacting both nuclear astrophysics and nuclear technology applications.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Min-Max Spaces and Complexity Reduction in Min-Max Expansions

Idempotent methods have been found to be extremely helpful in the numerical solution of certain classes of nonlinear control problems. In those methods, one uses the fact that the value function lies in the space of semiconvex functions (in the case of maximizing controllers), and approximates this value using a truncated max-plus basis expansion. In some classes, the value function is actually convex, and then one specifically approximates with suprema (i.e., max-plus sums) of affine functions. Note that the space of convex functions is a max-plus linear space, or moduloid. In extending those concepts to game problems, one finds a different function space, and different algebra, to be appropriate. Here we consider functions which may be represented using infima (i.e., min-max sums) of max-plus affine functions. It is natural to refer to the class of functions so represented as the min-max linear space (or moduloid) of max-plus hypo-convex functions. We examine this space, the associated notion of duality and min-max basis expansions. In using these methods for solution of control problems, and now games, a critical step is complexity-reduction. In particular, one needs to find reduced-complexity expansions which approximate the function as well as possible. We obtain a solution to this complexity-reduction problem in the case of min-max expansions.

97 MATHEMATICS AND COMPUTING↗

Hybrid preconditioning for iterative diagonalization of ill-conditioned generalized eigenvalue problems in electronic structure calculations

The iterative diagonalization of a sequence of large ill-conditioned generalized eigenvalue problems is a computational bottleneck in quantum mechanical methods employing a nonorthogonal basis for ab initio electronic structure calculations. We propose a hybrid preconditioning scheme to effectively combine global and locally accelerated preconditioners for rapid iterative diagonalization of such eigenvalue problems. In partition-of-unity finite-element (PUFE) pseudopotential density-functional calculations, employing a nonorthogonal basis, we show that the hybrid preconditioned block steepest descent method is a cost-effective eigensolver, outperforming current state-of-the-art global preconditioning schemes, and comparably efficient for the ill-conditioned generalized eigenvalue problems produced by PUFE as the locally optimal block preconditioned conjugate-gradient method for the well-conditioned standard eigenvalue problems produced by planewave methods.

74 ATOMIC AND MOLECULAR PHYSICS↗