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 793 records · Page 44

Development of a Technical Basis and Guidance for Advanced SMR Function Allocation

This report presents the results from three key activities for FY13 that influence the definition of new concepts of operations for advanced Small Modular Reactors (AdvSMR: a) the development of a framework for the analysis of the functional environmental, and structural attributes, b) the effect that new technologies and operational concepts would have on the way functions are allocated to humans or machines or combinations of the two, and c) the relationship between new concepts of operations, new function allocations, and human performance requirements.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Draft Function Allocation Framework and Preliminary Technical Basis for Advanced SMR Concepts of Operations

This report presents preliminary research results from the investigation into the development of new models and guidance for Concepts of Operations in advanced small modular reactor (AdvSMR) designs. AdvSMRs are nuclear power plants (NPPs), but unlike conventional large NPPs that are constructed on site, AdvSMRs systems and components will be fabricated in a factory and then assembled on site. AdvSMRs will also use advanced digital instrumentation and control systems, and make greater use of automation. Some AdvSMR designs also propose to be operated in a multi-unit configuration with a single central control room as a way to be more cost-competitive with existing NPPs. These differences from conventional NPPs not only pose technical and operational challenges, but they will undoubtedly also have regulatory compliance implications, especially with respect to staffing requirements and safety standards.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Variational treatment of electron–polyatomic-molecule scattering calculations using adaptive overset grids

In this paper, the complex Kohn variational method for electron–polyatomic-molecule scattering is formulated using an overset-grid representation of the scattering wave function. The overset grid consists of a central grid and multiple dense atom-centered subgrids that allow the simultaneous spherical expansions of the wave function about multiple centers. Scattering boundary conditions are enforced by using a basis formed by the repeated application of the free-particle Green's function and potential G ^ 0 + V ^ on the overset grid in a Born-Arnoldi solution of the working equations. The theory is shown to be equivalent to a specific Padé approximant to the T matrix and has rapid convergence properties, in both the number of numerical basis functions employed and the number of partial waves employed in the spherical expansions. The method is demonstrated in calculations on methane and CF 4 in the static-exchange approximation and compared in detail with calculations performed with the numerical Schwinger variational approach based on single-center expansions. Finally, an efficient procedure for operating with the free-particle Green's function and exchange operators (to which no approximation is made) is also described.

74 ATOMIC AND MOLECULAR PHYSICS↗

Heterogeneous Mixtures of Dictionary Functions to Approximate Subspace Invariance in Koopman Operators: Why Deep Koopman Operators Work

Abstract Koopman operators model nonlinear dynamics as a linear dynamic system acting on a nonlinear function as the state. This nonstandard state is often called a Koopman observable and is usually approximated numerically by a superposition of functions drawn from a dictionary . In a widely used algorithm, extended dynamic mode decomposition (EDMD), the dictionary functions are drawn from a fixed class of functions. Deep learning combined with EDMD has been used to learn novel dictionary functions in an algorithm called deep dynamic mode decomposition (deepDMD). The learned representation both (1) accurately models and (2) scales well with the dimension of the original nonlinear system. In this paper, we analyze the learned dictionaries from deepDMD and explore the theoretical basis for their strong performance. We explore State-Inclusive Logistic Lifting (SILL) dictionary functions to approximate Koopman observables. Error analysis of these dictionary functions show they satisfy a property of subspace approximation, which we define as uniform finite approximate closure. Typically, a Koopman dictionary’s nonlinear functions are homogeneous. In this paper, we discover that structured mixing of heterogeneous dictionary functions drawn from different classes of nonlinear functions achieve the same accuracy and dimensional scaling as the deep-learning-based deepDMD algorithm Yeung et al. ( In: 2019 American Control Conference (ACC), 2019). We specifically show this by building a heterogeneous dictionary comprised of SILL functions and conjunctive radial basis functions (RBFs). This mixed dictionary achieves similar accuracy and dimensional scaling to deepDMD with an order of magnitude reduction in parameters, while maintaining geometric interpretability. These results strengthen the viability of dictionary-based Koopman models to solving high-dimensional nonlinear learning problems.

Johnson, Charles A.↗

Wave function continuity and the diagonal Born-Oppenheimer correction at conical intersections

We demonstrate that though exact in principle, the expansion of the total molecular wave function as a sum over adiabatic Born-Oppenheimer (BO) vibronic states makes inclusion of the second-derivative nonadiabatic energy term near conical intersections practically problematic. In order to construct a well-behaved molecular wave function that has density at a conical intersection, the individual BO vibronic states in the summation must be discontinuous. When the second-derivative nonadiabatic terms are added to the Hamiltonian, singularities in the diagonal BO corrections (DBOCs) of the individual BO states arise from these discontinuities. In contrast to the well-known singularities in the first-derivative couplings at conical intersections, these singularities are non-integrable, resulting in undefined DBOC matrix elements. Though these singularities suggest that the exact molecular wave function may not have density at the conical intersection point, there is no physical basis for this constraint. Instead, the singularities are artifacts of the chosen basis of discontinuous functions. We also demonstrate that continuity of the total molecular wave function does not require continuity of the individual adiabatic nuclear wave functions. We classify nonadiabatic molecular dynamics methods according to the constraints placed on wave function continuity and analyze their formal properties. Based on our analysis, it is recommended that the DBOC be neglected when employing mixed quantum-classical methods and certain approximate quantum dynamical methods in the adiabatic representation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Structural basis for enzymatic terminal C–H bond functionalization of alkanes

Alkane monooxygenase (AlkB) is a widely occurring integral membrane metalloenzyme that catalyzes the initial step in the functionalization of recalcitrant alkanes with high terminal selectivity. AlkB enables diverse microorganisms to use alkanes as their sole carbon and energy source. Here we present the 48.6-kDa cryo-electron microscopy structure of a natural fusion from Fontimonas thermophila between AlkB and its electron donor AlkG at 2.76 Å resolution. The AlkB portion contains six transmembrane helices with an alkane entry tunnel within its transmembrane domain. A dodecane substrate is oriented by hydrophobic tunnel-lining residues to present a terminal C–H bond toward a diiron active site. AlkG, an [Fe–4S] rubredoxin, docks via electrostatic interactions and sequentially transfers electrons to the diiron center. The archetypal structural complex presented reveals the basis for terminal C–H selectivity and functionalization within this broadly distributed evolutionary class of enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

A three-body model of the {sup 11}B nucleus

The binding energy and the rms charge and mass radii have been calculated in terms of the single-channel three-body {sup 4}He{sup 4}He{sup 3}H model of the {sup 11}B nucleus with an expansion of the three-body wave function in a nonorthogonal Gaussian basis. Parameters of the wave function are presented and convergence of the three-body energy depending on the number of expansion terms is demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A projection-free method for representing plane-wave DFT results in an atom-centered basis

Plane wave density functional theory (DFT) is a powerful tool for gaining accurate, atomic level insight into bulk and surface structures. Yet, the delocalized nature of the plane wave basis set hinders the application of many powerful post-computation analysis approaches, many of which rely on localized atom-centered basis sets. Traditionally, this gap has been bridged via projection-based techniques from a plane wave to atom-centered basis. We instead propose an alternative projection-free approach utilizing direct calculation of matrix elements of the converged plane wave DFT Hamiltonian in an atom-centered basis. This projection-free approach yields a number of compelling advantages, including strict orthonormality of the resulting bands without artificial band mixing and access to the Hamiltonian matrix elements, while faithfully preserving the underlying DFT band structure. The resulting atomic orbital representation of the Kohn-Sham wavefunction and Hamiltonian provides a gateway to a wide variety of analysis approaches. We demonstrate the utility of the approach for a diverse set of chemical systems and example analysis approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Edge-Functionalized Graphene Nanoribbon Encapsulation To Enhance Stability and Control Kinetics of Hydrogen Storage Materials

Hydrogen is a long-term clean energy carrier that enables completely carbon-free energy production. However, practical implementation of hydrogen fuel technologies is restricted because of lack of safe and high-performing storage materials. Here, we report Mg nanocrystals encapsulated by narrow, bottom-up synthesized graphene nanoribbons (GNRs) as environmentally stable and high-capacity hydrogen storage materials. As an encapsulation medium, GNRs offer similar functionalities as reduced graphene oxide to protect the encapsulated Mg nanocrystals from extensive oxidation, while allowing penetrations of hydrogen. In addition, the GNRs can be edge functionalized to tune the (de-)hydrogenation kinetics, in particular for the processes occurred at the GNR–Mg interfaces. In this work, four different types of edge-functional groups were introduced into GNRs with the goal of comparing their cycling performances because of edge functionalization. On the basis of detailed kinetic analysis coupled with first-principles calculations, we propose that edge-functional groups can contribute to the reduction of kinetic barriers for surface hydrogen reactions at the interface with the GNR by stabilizing surface defects. Furthermore, the GNR–Mg composite exhibited higher hydrogen storage capacity (7.1 wt % of hydrogen based on the total composite) compared with the current literature while demonstrating long-term air stability. This work suggests that the rational design of edge-functional groups in graphene derivatives can provide a general and novel paradigm for simultaneous encapsulation and hydrogen storage catalysis in simple metal or complex metal nanocrystals.

08 HYDROGEN↗

Coupled-cluster, Möller Plesset (MP2), Density Fitted Local MP2, and Density Functional Theory Examination of the Energetic and Structural Features of Hydrophobic Solvation: Water and Pentane

The interaction potentials between immiscible polar and non-polar solvents are a major driving force behind the formation of liquid:liquid interfaces. In this work, the interaction energy of water–pentane dimer has been determined using coupled-cluster theory with single double (triple) excitations [CCSD(T)], 2nd order Möller Plesset perturbation theory (MP2), density fitted local MP2 (DF-LMP2), as well as density functional theory using a wide variety of density functionals and several different basis sets. The M05-2X exchange correlation functionals exhibit excellent agreement with CCSD(T) and DF-LMP2 after taking into account basis set superposition error. The gas phase water–pentane interaction energy is found to be quite sensitive to the specific pentane isomer (2,2- dimethylpropane vs. n-pentane) and relative orientation of the monomeric constituents. Subsequent solution phase cluster calculations of 2,2-dimethylpropane and n-pentane solvated by water indicate a positive free energy of solvation that is in good agreement with available experimental data. Structural parameters are quite sensitive to the density functional employed and reflect differences in the two-body interaction energy calculated by each method. In contrast, cluster calculations of pentane solvation of H2O solute are found to be inadequate for describing the organic solvent, likely due to limitations associated with the functionals employed (B3LYP, BHandH, and M05-2X).

Environmental Molecular Sciences Laboratory↗

Effective empirical corrections for basis set superposition error in the def2-SVPD basis: gCP and DFT-C

With the aim of mitigating the basis set error in density functional theory (DFT) calculations employing local basis sets, we herein develop two empirical corrections for basis set superposition error (BSSE) in the def2-SVPD basis, a basis which - when stripped of BSSE - is capable of providing near-complete-basis DFT results for non-covalent interactions. Specifically, we adapt the existing pairwise geometrical counterpoise (gCP) approach to the def2-SVPD basis, and we develop a beyond-pairwise approach, DFT-C, which we parameterize across a small set of intermolecular interactions. Both gCP and DFT-C are evaluated against the traditional Boys-Bernardi counterpoise correction across a set of 3402 non-covalent binding energies and isomerization energies. We find that the DFT-C method represents a significant improvement over gCP, particularly for non-covalently-interacting molecular clusters. Furthermore, DFT-C is transferable among density functionals and can be combined with existing functionals - such as B97M-V - to recover large-basis results at a fraction of the cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Linear-scaling implementation of the direct random-phase approximation

We report the linear-scaling implementation of the direct random-phase approximation (dRPA) for closed-shell molecular systems. As a bonus, linear-scaling algorithms are also presented for the second-order screened exchange extension of dRPA as well as for the second-order Møller–Plesset (MP2) method and its spin-scaled variants. Our approach is based on an incremental scheme which is an extension of our previous local correlation method [Rolik et al., J. Chem. Phys. 139, 094105 (2013)]. The approach extensively uses local natural orbitals to reduce the size of the molecular orbital basis of local correlation domains. In addition, we also demonstrate that using natural auxiliary functions [M. Kállay, J. Chem. Phys. 141, 244113 (2014)], the size of the auxiliary basis of the domains and thus that of the three-center Coulomb integral lists can be reduced by an order of magnitude, which results in significant savings in computation time. The new approach is validated by extensive test calculations for energies and energy differences. Our benchmark calculations also demonstrate that the new method enables dRPA calculations for molecules with more than 1000 atoms and 10 000 basis functions on a single processor.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessment of enzyme active site positioning and tests of catalytic mechanisms through X-ray–derived conformational ensembles

How enzymes achieve their enormous rate enhancements remains a central question in biology, and our understanding to date has impacted drug development, influenced enzyme design, and deepened our appreciation of evolutionary processes. While enzymes position catalytic and reactant groups in active sites, physics requires that atoms undergo constant motion. Numerous proposals have invoked positioning or motions as central for enzyme function, but a scarcity of experimental data has limited our understanding of positioning and motion, their relative importance, and their changes through the enzyme’s reaction cycle. To examine positioning and motions and test catalytic proposals, we collected “room temperature” X-ray crystallography data for Pseudomonas putida ketosteroid isomerase (KSI), and we obtained conformational ensembles for this and a homologous KSI from multiple PDB crystal structures. Here, ensemble analyses indicated limited change through KSI’s reaction cycle. Active site positioning was on the 1- to 1.5-Å scale, and was not exceptional compared to noncatalytic groups. The KSI ensembles provided evidence against catalytic proposals invoking oxyanion hole geometric discrimination between the ground state and transition state or highly precise general base positioning. Instead, increasing or decreasing positioning of KSI’s general base reduced catalysis, suggesting optimized Ångstrom-scale conformational heterogeneity that allows KSI to efficiently catalyze multiple reaction steps. Ensemble analyses of surrounding groups for WT and mutant KSIs provided insights into the forces and interactions that allow and limit active-site motions. Most generally, this ensemble perspective extends traditional structure–function relationships, providing the basis for a new era of “ensemble–function” interrogation of enzymes.

59 BASIC BIOLOGICAL SCIENCES↗

3D spherical functional expansion tallies in Serpent 2 Monte Carlo code

This work extends the application of functional expansion tallies to 3D spherical geometries. The 3D Zernike polynomials are set as an orthonormal polynomials basis for the functional reconstruction. The study describes the construction of the complete set of polynomials, a natural expansion of the spherical harmonics polynomials where 3D Zernike moments can be evaluated as a linear combination of the geometrical moments. The 3D Zernike polynomials formulation and the computational approach implemented in Serpent 2 are presented and tested through the Godiva model from the ICSBEP criticality benchmark test cases. The implementation results are in agreement with a reference solution described in a fine-resolution mesh, enhancing also the performance and memory demand. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Helical Organic and Inorganic Polymers

Despite being a staple of synthetic plastics and biomolecules, helical polymers are scarcely studied with Gaussian-basis-set ab initio electron-correlated methods on an equal footing with molecules. This article introduces an ab initio second-order many-body Green’s-function [MBGF(2)] method with nondiagonal, frequency-dependent Dyson self-energy for infinite helical polymers using screw-axis-symmetry-adapted Gaussian-spherical-harmonics basis functions. Together with the Gaussian-basis-set density-functional theory for energies, analytical atomic forces, translational-period force, and helical-angle force, it can compute correlated energy, quasiparticle energy bands, structures, and vibrational frequencies of an infinite helical polymer, which smoothly converge at the corresponding oligomer results. These methods can handle incommensurable structures, which have an infinite translational period and are hard to characterize by any other method, just as efficiently as commensurable structures. We apply them to polyethylene (2/1 helix), polyacetylene (Peierls’ system), and polytetrafluoroethylene (13/6 helix) to establish the quantitative accuracy of MBGF(2)/cc-pVDZ in simulating their (angle-resolved) ultraviolet photoelectron spectra, and of B3LYP/cc-pVDZ or 6-31G** in reproducing their structures, infrared and Raman band positions, phonon dispersions, and (coherent and incoherent) inelastic neutron scattering spectra. We then predict the same properties for infinitely catenated chains of nitrogen or oxygen and discuss their possible metastable existence under ambient conditions. They include planar zigzag polyazene (N2)x (Peierls’ system), 11/3-helical isotactic polyazane (NH)x, 9/4-helical isotactic polyfluoroazane (NF)x, and 7/2-helical polyoxane (O)x as potential high-energy-density materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Simulations of Attosecond Metallization in Quartz and Diamond Probed with Inner-Shell Transient Absorption Spectroscopy

When dielectrics are hit with intense infrared (IR) laser pulses, transient metalization can occur. The initial attosecond dynamics behind this metallization are not entirely understood. Therefore, simulations are needed to understand this process and to help interpret experimental observations of it, such as with attosecond transient absorption (ATA). In this paper, we present first-principles simulations of ATA based on bulk-mimicking clusters and real-time time-dependent density functional theory (RT-TDDFT), with Koopmans-tuned range-separated hybrid functionals and Gaussian basis sets. Our method gives good agreement with the experiment for the breakdown threshold in silica and diamond. This breakdown voltage corresponds to a Keldysh parameter of approximately one and thus involves a transition to a regime where the dynamics are driven by tunneling. Pumping at an amplitude just below this value causes a mixture of multiphoton and tunneling excitations across the band gap to occur. The computed extreme ultraviolet and X-ray attosecond transient spectra also agree well with the experiment and show a decrease in optical density due to the transient population of the conduction band from the IR field. First-principles approaches such as this are valuable for interpreting the complicated modulations in a spectrum and for guiding future attosecond experiments on solids.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Momentum-space Argonne V18 interaction

This paper gives a momentum-space representation of the Argonne V18 potential as an expansion in products of spin-isospin operators with scalar coefficient functions of the momentum transfer. Two representations of the scalar coefficient functions for the strong part of the interaction are given. One is as an expansion in an orthonormal basis of rational functions and the other as an expansion in Chebyshev polynomials on different intervals. Both provide practical and efficient representations for computing the momentum-space potential that do not require integration or interpolation. Programs based on both expansions are available as supplementary material. Analytic expressions are given for the scalar coefficient functions of the Fourier transform of the electromagnetic part of the Argonne V18. A simple method for computing the partial-wave projections of these interactions from the operator expressions is also given.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Basis convergence of range-separated density-functional theory

Range-separated density-functional theory (DFT) is an alternative approach to Kohn-Sham density-functional theory. The strategy of range-separated density-functional theory consists in separating the Coulomb electron-electron interaction into long-range and short-range components and treating the long-range part by an explicit many-body wave-function method and the short-range part by a density-functional approximation. Among the advantages of using many-body methods for the long-range part of the electron-electron interaction is that they are much less sensitive to the one-electron atomic basis compared to the case of the standard Coulomb interaction. Here, we provide a detailed study of the basis convergence of range-separated density-functional theory. We study the convergence of the partial-wave expansion of the long-range wave function near the electron-electron coalescence. We show that the rate of convergence is exponential with respect to the maximal angular momentum L for the long-range wave function, whereas it is polynomial for the case of the Coulomb interaction. We also study the convergence of the long-range second-order Møller-Plesset correlation energy of four systems (He, Ne, N{sub 2}, and H{sub 2}O) with cardinal number X of the Dunning basis sets cc − p(C)V XZ and find that the error in the correlation energy is best fitted by an exponential in X. This leads us to propose a three-point complete-basis-set extrapolation scheme for range-separated density-functional theory based on an exponential formula.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗