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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Data-Centric Approach to Capture Non-Polynomial Nonlinear Dynamics

We propose an analytical construction of observable functions in the extended dynamic mode decomposition (EDMD) algorithm. EDMD is a numerical method for approximating the spectral properties of the Koopman operator. The choice of observable functions is fundamental for applying EDMD to nonlinear problems arising in systems and control. Existing methods either start from a set of dictionary functions and look for the subset that best fits the underlying nonlinear dynamics or rely on machine learning algorithms to “learn” observable functions. Conversely, in this paper, we start from the dynamical system model and lift it through the Lie derivatives, rendering it into a polynomial form. This proposed transformation into a polynomial form is exact and provides an adequate set of observable functions. The strength of the proposed approach is its applicability to a broader class of nonlinear dynamical systems, particularly those with nonpolynomial functions and compositions thereof. Moreover, it retains the physical interpretability of the underlying dynamical system and can be readily integrated into existing numerical libraries. We demonstrate the proposed approach with an application to electric power systems. The modeled system consists of a single generator connected to an infinite bus, where nonlinear terms include sine and cosine functions. The results demonstrate the effectiveness of the proposed procedure in off-attractor nonlinear dynamics for estimation and prediction; the observable functions obtained from the proposed construction outperform methods that use dictionary functions comprising monomials or radial basis functions.

extended dynamic mode decomposition↗
DOE OSTI2016

A Nonorthogonal State-Interaction Approach for Matrix Product State Wave Functions

We present a state-interaction approach for matrix product state (MPS) wave functions in a nonorthogonal molecular orbital basis. Our approach allows us to calculate, for example, transition and spin–orbit coupling matrix elements between arbitrary electronic states, provided that they share the same one-electron basis functions and size of the active orbital space, respectively. The key element is the transformation of the MPS wave functions of different states from a nonorthogonal to a biorthonormal molecular orbital basis representation, by exploiting a sequence of nonunitary transformations, following a proposal by Malmqvist [Int. J. Quantum Chem.1986, 30, 479]. This is well-known for traditional wave function parametrizations but has not yet been exploited for MPS wave functions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2011

Generalized variational principle for excited states using nodes of trial functions

The familiar variational principle provides an upper bound to the ground-state energy of a given Hamiltonian. This allows one to optimize a trial wave function by minimizing the expectation value of the energy. This approach is also trivially generalized to excited states, so that given a trial wave function of a certain symmetry, one can compute an upper bound to the lowest-energy level of that symmetry. In order to generalize further and build an upper bound of an arbitrary excited state of the desired symmetry, a linear combination of basis functions is generally used to generate an orthogonal set of trial functions, all bounding their respective states. However, sometimes a compact wave-function form is sought, and a basis-set expansion is not desirable or possible. Here we present an alternative generalization of the variational principle to excited states that does not require explicit orthogonalization to lower-energy states. It is valid for one-dimensional systems and, with additional information, to at least some n-dimensional systems. This generalized variational principle exploits information about the nodal structure of the trial wave function, giving an upper bound to the exact energy without the need to build a linear combination of basis functions. To illustrate the theorem we apply it to a nontrivial example: the 1s2s {sup 1} S excited state of the helium atom.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Multiresolution quantum chemistry in multiwavelet bases: excited states from time-dependent Hartree–Fock and density functional theory via linear response

Using the fully numerical method for time-dependent Hartree–Fock and density functional theory (TD-HF/DFT) with the Tamm–Dancoff (TD) approximation we use a multiresolution analysis (MRA) approach to present our findings. From a reformulation with effective use of the density matrix operator, we obtain a general form of the HF/DFT linear response equation in the first quantization formalism. It can be readily rewritten as an integral equation with the bound-state Helmholtz (BSH) kernel for the Green's function. The MRA implementation of the resultant equation permits excited state calculations without virtual orbitals. Moreover, the integral equation is efficiently and adaptively solved using a numerical multiresolution solver with multiwavelet bases. Our implementation of the TD-HF/DFT methods is applied for calculating the excitation energies of H 2 , Be, N 2 , H 2 O, and C 2 H 4 molecules. The numerical errors of the calculated excitation energies converge in proportion to the residuals of the equation in the molecular orbitals and response functions. The energies of the excited states at a variety of length scales ranging from short-range valence excitations to long-range Rydberg-type ones are consistently accurate. It is shown that the multiresolution calculations yield the correct exponential asymptotic tails for the response functions, whereas those computed with Gaussian basis functions are too diffuse or decay too rapidly. Finally, we introduce a simple asymptotic correction to the local spin-density approximation (LSDA) so that in the TDDFT calculations, the excited states are correctly bound.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Flexible nuclear screening approximation to the two-electron spin–orbit coupling based on ab initio parameterization

The derivation, implementation, and validation of a new approximation to the two-electron spin–orbit coupling (SOC) terms is reported. The approximation, referred to as flexible nuclear screening spin–orbit, is based on the effective one-electron spin–orbit operator and accounts for two-electron SOC effects by screening nuclear charges. A highly flexible scheme for the nuclear screening is developed, mainly using parameterization based on ab initio atomic SOC calculations. Tabulated screening parameters are provided for contracted and primitive Gaussian-type basis functions of the ANO-RCC basis set for elements from H to Cm. The strategy for their adaptation to any other Gaussian basis set is presented and validated. A model to correct for the effect of splitting of transition metal d orbitals on their SOC matrix elements is introduced. The method is applied to a representative set of molecules, and compared to exact treatment and other approximative approaches at the same level of relativistic theory. The calculated SOC matrix elements are in very good agreement with their “exact” values; deviation below 1% is observed on average. The presented approximation is considered to be generally applicable, simple to implement, highly efficient, and accurate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2011

Reducing sequencing complexity in dynamical quantum error suppression by Walsh modulation

We study dynamical error suppression from the perspective of reducing sequencing complexity, with an eye toward facilitating the development of efficient semiautonomous quantum-coherent systems. To this end, we focus on digital sequences where all interpulse time periods are integer multiples of a minimum clock period and compatibility with digital classical control circuitry is intrinsic. We use so-called Walsh functions as a unifying mathematical framework; the Walsh functions are an orthonormal set of basis functions which may be associated directly with the control propagator for a digital modulation scheme. Using this insight, we characterize the suite of resulting Walsh dynamical decoupling sequences--including both familiar and novel control sequences--and identify the number of periodic square-wave (Rademacher) functions required to generate the associated Walsh function as the key determinant of the error-suppressing features. We also show how Walsh modulation may be employed for the protection of certain nontrivial logic gates. Based on these insights, we identify Walsh modulation as a digital-efficient approach for physical-layer error suppression.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2012

An Adaptive Wavelet Stochastic Collocation Method for Irregular Solutions of Stochastic Partial Differential Equations

Accurate predictive simulations of complex real world applications require numerical approximations to first, oppose the curse of dimensionality and second, converge quickly in the presence of steep gradients, sharp transitions, bifurcations or finite discontinuities in high-dimensional parameter spaces. In this paper we present a novel multi-dimensional multi-resolution adaptive (MdMrA) sparse grid stochastic collocation method, that utilizes hierarchical multiscale piecewise Riesz basis functions constructed from interpolating wavelets. The basis for our non-intrusive method forms a stable multiscale splitting and thus, optimal adaptation is achieved. Error estimates and numerical examples will used to compare the efficiency of the method with several other techniques.

adaptive sparse grids↗
NASA NTRS1995

The Schwinger Variational Method

Variational methods have proven invaluable in theoretical physics and chemistry, both for bound state problems and for the study of collision phenomena. The application of the Schwinger variational (SV) method to e-molecule collisions and molecular photoionization has been reviewed previously. The present chapter discusses the implementation of the SV method as applied to e-molecule collisions. Since this is not a review of cross section data, cross sections are presented only to server as illustrative examples. In the SV method, the correct boundary condition is automatically incorporated through the use of Green's function. Thus SV calculations can employ basis functions with arbitrary boundary conditions. The iterative Schwinger method has been used extensively to study molecular photoionization. For e-molecule collisions, it is used at the static exchange level to study elastic scattering and coupled with the distorted wave approximation to study electronically inelastic scattering.

Huo, Winifred M.↗
NASA NTRS2003

Electron Collisional Excitation Rates for OI USING the B-Spline R-Matrix Approach

The B-spline R-matrix approach has been used to calculate electron collisional excitation strengths and rates for transitions between the 3P, 1D, and 1S states of ground configuration and from these states to the states of the excited 2s(sup 2)2p(sup 3)ns (n = 3-5), 2s(sup 2)2p(sup 3)np (n = 3-4), 2s(sup 2)2p(sup 3)nd (n = 3-4), 2s(sup 2)2p(sup 3)4f, and 2s2p(sup 5) configurations. The nonorthogonal orbitals are used for an accurate description of both the target wave functions and the R-matrix basis functions. The thermally averaged collision strengths are obtained from the collision strengths by integrating over a Maxwellian velocity distribution of electron energies, and these are tabulated over a temperature range from 1000 to 60,000 K. The parametric functions of scaled energy have also been obtained to represent collision strengths over a wide energy range or thermally averaged collision strengths at any desired temperature.

Zatsarinny, O.↗
NASA NTRS2006

New Accurate Oscillator Strengths and Electron Excitation Collision Strengths for N1

The nonorthogonal orbitals technique in a multiconfiguration Hartree-Fock approach is used to calculate oscillator strengths and transition probabilities of N(I) lines. The relativistic effects are allowed by means of Breit-Pauli operators. The length and velocity forms of oscillator strengths show good agreement for most transitions. The B-spline R-matrix with pseudostates approach has been used to calculate electron excitation collision strengths and rates. The nonorthogonal orbitals are used for an accurate description of both target wave functions and the R-matrix basis functions. The 24 spectroscopic bound and autoionizing states together with 15 pseudostates are included in the close-coupling expansion. The collision strengths for transitions between fine-structure levels are calculated by transforming the LS-coupled K-matrices to K-matrices in an intermediate coupling scheme. Thermally averaged collision strengths have been determined by integrating collision strength over a Maxwellian distribution of electron energies over a temperature range suitable for the modeling of astrophysical plasmas. The oscillator strengths and thermally averaged collision strengths are presented for transitions between the fine-structure levels of the 2s(sup 2)p(sup 3) (sup 4)S(sup 0), (sup 2)D(sup 0), (sup 2)P(sup 0), 2s2p(sup 4) (sup 4)P, 2s(sup 2)2p(sup 2)3s (sup 4)P, and (sup 2)P terms and from these levels to the levels of the 2s(sup 2)2p(sup 2)3p (sup 2)S(sup 0), (sup 4)D(sup 0), (sup 4)P(sup 0), (sup 4)S(sup 0), (sup 2)D(sup 0), (sup 2)P(sup 0),2s(sup 2)2p(sup 2)3s(sup 2)D, 2s(sup 2)2p(sup 2)4s(sup 4)P, (sup 2)P, 2s(sup 2)2p(sup 2)3d(sup 2)P, (sup 4)F,(sup 2)F,(sup 4)P, (sup 4)D, and (sup 2)D terms. Thermally averaged collision strengths are tabulated over a temperature range from 500 to 50,000 K.

Tayal, S. S.↗
NASA NTRS1994

Scattering from perfectly conducting and resistive strips on a grounded dielectric slab

The scattering properties of perfectly conducting and resistive strips are predicted for strips which are located on a dielectric slab backed by a perfectly conducting ground plane. The spectral domain Green's function is used to relate the currents and fields on the strip, and the resulting integral equation is solved using the method of moments. Both TE and TM strips are examined using piecewise linear and pulse subdomain basis functions, respectively, to model the current on the strip. Calculated results are compared with results measured at the NASA Langley Research Center.

PULSE SUBDOMAIN BASIS FUNCTION↗
DOE OSTI2014

Accurate adiabatic correction in the hydrogen molecule

A new formalism for the accurate treatment of adiabatic effects in the hydrogen molecule is presented, in which the electronic wave function is expanded in the James-Coolidge basis functions. Systematic increase in the size of the basis set permits estimation of the accuracy. Numerical results for the adiabatic correction to the Born-Oppenheimer interaction energy reveal a relative precision of 10{sup −12} at an arbitrary internuclear distance. Such calculations have been performed for 88 internuclear distances in the range of 0 < R ⩽ 12 bohrs to construct the adiabatic correction potential and to solve the nuclear Schrödinger equation. Finally, the adiabatic correction to the dissociation energies of all rovibrational levels in H{sub 2}, HD, HT, D{sub 2}, DT, and T{sub 2} has been determined. For the ground state of H{sub 2} the estimated precision is 3 × 10{sup −7} cm{sup −1}, which is almost three orders of magnitude higher than that of the best previous result. The achieved accuracy removes the adiabatic contribution from the overall error budget of the present day theoretical predictions for the rovibrational levels.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

The benzene⋯naphthalene complex: A more challenging system than the benzene dimer for newly developed computational methods

High-level coupled cluster singles, doubles, and perturbative triples [CCSD(T)] computations with up to the aug-cc-pVQZ basis set (1924 basis functions) and various extrapolations toward the complete basis set (CBS) limit are presented for the sandwich, T-shaped, and parallel-displaced benzene⋯naphthalene complex. Using the CCSD(T)/CBS interaction energies as a benchmark, the performance of some newly developed wave function and density functional theory methods has been evaluated. The best performing methods were found to be the dispersion-corrected PBE0 functional (PBE0-D3) and spin-component scaled zeroth-order symmetry-adapted perturbation theory (SCS-SAPT0). The success of SCS-SAPT0 is very encouraging because it provides one method for energy component analysis of π-stacked complexes with 200 atoms or more. Most newly developed methods do, however, overestimate the interaction energies. The results of energy component analysis show that interaction energies are overestimated mainly due to the overestimation of dispersion energy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Strong Correlation DMRG and DFT

This project developed new ways to improve computer simulations of materials where electrons interact strongly with each other, a challenge for today’s most widely used method, density functional theory (DFT). We used an exact numerical method, the density matrix renormalization group (DMRG), to create highly accurate reference results for simple model systems, and used these to test DFT, prove when it will converge, and even train machine-learned functionals. We also invented new kinds of localized basis functions (“gausslets” and “multi-sliced gausslets”) and a “sliced-basis” approach that make high-accuracy simulations faster and more practical. These methods were applied to extended hydrogen systems, enabling the direct derivation of accurate low-energy models from first-principles calculations. We also introduced a new formalism, Conditional-Probability DFT, which could bypass traditional approximations. The tools and results from this work, including open-source software releases, will help scientists design and understand complex quantum materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

A meshless Galerkin method for non-local diffusion using localized kernel bases

Here, we introduce a meshless method for solving both continuous and discrete variational formulations of a volume constrained, nonlocal diffusion problem. We use the discrete solution to approximate the continuous solution. Our method is nonconforming and uses a localized Lagrange basis that is constructed out of radial basis functions. By verifying that certain inf-sup conditions hold, we demonstrate that both the continuous and discrete problems are well-posed, and also present numerical and theoretical results for the convergence behavior of the method. The stiffness matrix is assembled by a special quadrature routine unique to the localized basis. Combining the quadrature method with the localized basis produces a well-conditioned, symmetric matrix. This then is used to find the discretized solution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

A Massively Parallel Implementation of the CCSD(T) Method Using the Resolution-of-the-Identity Approximation and a Hybrid Distributed/Shared Memory Parallelization Model

In this work, a parallel algorithm is described for the coupled-cluster singles and doubles method augmented with a perturbative correction for triple excitations [CCSD(T)] using the resolution-of-the-identity (RI) approximation for two-electron repulsion integrals (ERIs). The algorithm bypasses the storage of four-center ERIs by adopting an integral-direct strategy. The CCSD amplitude equations are given in a compact quasi-linear form by factorizing them in terms of amplitude-dressed three-center intermediates. A hybrid MPI/OpenMP parallelization scheme is employed, which uses the OpenMP-based shared memory model for intranode parallelization and the MPI-based distributed memory model for internode parallelization. Parallel efficiency has been optimized for all terms in the CCSD amplitude equations. Two different algorithms have been implemented for the rate-limiting terms in the CCSD amplitude equations that entail and -scaling computational costs, where N O and N V denote the number of correlated occupied and virtual orbitals, respectively. One of the algorithms assembles the four-center ERIs requiring N V 4 and N O 2 N V 2 -scaling memory costs in a distributed manner on a number of MPI ranks, while the other algorithm completely bypasses the assembling of quartic memory-scaling ERIs and thus largely reduces the memory demand. It is demonstrated that the former memory-expensive algorithm is faster on a few hundred cores, while the latter memory-economic algorithm shows a better strong scaling in the limit of a few thousand cores. The program is shown to exhibit a near-linear scaling, in particular for the compute-intensive triples correction step, on up to 8000 cores. The performance of the program is demonstrated via calculations involving molecules with 24–51 atoms and up to 1624 atomic basis functions. As the first application, the complete basis set (CBS) limit for the interaction energy of the π-stacked uracil dimer from the S66 data set has been investigated. This work reports the first calculation of the interaction energy at the CCSD(T)/aug-cc-pVQZ level without local orbital approximation. The CBS limit for the CCSD correlation contribution to the interaction energy was found to be -8.01 kcal/mol, which agrees very well with the value -7.99 kcal/mol reported by Schmitz, Hättig, and Tew [ Phys. Chem. Chem. Phys. 2014 , 16 , 22167-22178]. The CBS limit for the total interaction energy was estimated to be -9.64 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Non-perturbative calculation of molecular magnetic properties within current-density functional theory

We present a novel implementation of Kohn–Sham density-functional theory utilizing London atomic orbitals as basis functions. External magnetic fields are treated non-perturbatively, which enable the study of both magnetic response properties and the effects of strong fields, using either standard density functionals or current-density functionals—the implementation is the first fully self-consistent implementation of the latter for molecules. Pilot applications are presented for the finite-field calculation of molecular magnetizabilities, hypermagnetizabilities, and nuclear magnetic resonance shielding constants, focusing on the impact of current-density functionals on the accuracy of the results. Existing current-density functionals based on the gauge-invariant vorticity are tested and found to be sensitive to numerical details of their implementation. Furthermore, when appropriately regularized, the resulting magnetic properties show no improvement over standard density-functional results. An advantage of the present implementation is the ability to apply density-functional theory to molecules in very strong magnetic fields, where the perturbative approach breaks down. Comparison with high accuracy full-configuration-interaction results show that the inadequacies of current-density approximations are exacerbated with increasing magnetic field strength. Standard density-functionals remain well behaved but fail to deliver high accuracy. The need for improved current-dependent density-functionals, and how they may be tested using the presented implementation, is discussed in light of our findings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗