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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

Chaos control of the brushless direct current motor using adaptive dynamic surface control based on neural network with the minimum weights

This paper investigates chaos control for the brushless DC motor (BLDCM) system by adaptive dynamic surface approach based on neural network with the minimum weights. The BLDCM system contains parameter perturbation, chaotic behavior, and uncertainty. With the help of radial basis function (RBF) neural network to approximate the unknown nonlinear functions, the adaptive law is established to overcome uncertainty of the control gain. By introducing the RBF neural network and adaptive technology into the dynamic surface control design, a robust chaos control scheme is developed. It is proved that the proposed control approach can guarantee that all signals in the closed-loop system are globally uniformly bounded, and the tracking error converges to a small neighborhood of the origin. Simulation results are provided to show that the proposed approach works well in suppressing chaos and parameter perturbation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Analytical gradients of the state-average complete active space self-consistent field method with density fitting

An efficient implementation of the state-averaged complete active space self-consistent field (SA-CASSCF) gradients employing density fitting (DF) is presented. The DF allows a reduction both in scaling and prefactors of the different steps involved. The performance of the algorithm is demonstrated on a set of molecules ranging up to an iron-Heme b complex which with its 79 atoms and 811 basis functions is to our knowledge the largest SA-CASSCF gradient computed. For smaller systems where the conventional code could still be used as a reference, both the linear response calculation and the gradient formation showed a clear timing reduction and the overall cost of a geometry optimization is typically reduced by more than one order of magnitude while the accuracy loss is negligible.

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

A multi-agent quantum Monte Carlo model for charge transport: Application to organic field-effect transistors

We have developed a multi-agent quantum Monte Carlo model to describe the spatial dynamics of multiple majority charge carriers during conduction of electric current in the channel of organic field-effect transistors. The charge carriers are treated by a neglect of diatomic differential overlap Hamiltonian using a lattice of hydrogen-like basis functions. The local ionization energy and local electron affinity defined previously map the bulk structure of the transistor channel to external potentials for the simulations of electron- and hole-conduction, respectively. The model is designed without a specific charge-transport mechanism like hopping- or band-transport in mind and does not arbitrarily localize charge. An electrode model allows dynamic injection and depletion of charge carriers according to source-drain voltage. The field-effect is modeled by using the source-gate voltage in a Metropolis-like acceptance criterion. Although the current cannot be calculated because the simulations have no time axis, using the number of Monte Carlo moves as pseudo-time gives results that resemble experimental I/V curves.

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

Restricted and unrestricted non-Hermitian Hartree-Fock: Theory, practical considerations, and applications to metastable molecular anions

This work describes the implementation and applications of non-Hermitian self-consistent field (NH-SCF) theory with complex basis functions for the ab initio computation of positions and widths of shape resonances in molecules. We utilize both the restricted open-shell and the previously unexplored spin-unrestricted variants to compute Siegert energies of several anionic shape resonances in small diatomic and polyatomic molecules including carbon tetrafluoride which has been the subject of several recent experimental studies. The computation of general molecular properties from a non-Hermitian wavefunction is discussed, and a density-based analysis is applied to the {sup 2}B{sub 1} shape resonance in formaldehyde. Spin-unrestricted NH-SCF is used to compute a complex potential energy surface for the carbon monoxide anion which correctly describes dissociation.

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

Uniform electron gases. III. Low-density gases on three-dimensional spheres

By combining variational Monte Carlo (VMC) and complete-basis-set limit Hartree-Fock (HF) calculations, we have obtained near-exact correlation energies for low-density same-spin electrons on a three-dimensional sphere (3-sphere), i.e., the surface of a four-dimensional ball. In the VMC calculations, we compare the efficacies of two types of one-electron basis functions for these strongly correlated systems and analyze the energy convergence with respect to the quality of the Jastrow factor. The HF calculations employ spherical Gaussian functions (SGFs) which are the curved-space analogs of Cartesian Gaussian functions. At low densities, the electrons become relatively localized into Wigner crystals, and the natural SGF centers are found by solving the Thomson problem (i.e., the minimum-energy arrangement of n point charges) on the 3-sphere for various values of n. We have found 11 special values of n whose Thomson sites are equivalent. Three of these are the vertices of four-dimensional Platonic solids — the hyper-tetrahedron (n = 5), the hyper-octahedron (n = 8), and the 24-cell (n = 24) — and a fourth is a highly symmetric structure (n = 13) which has not previously been reported. By calculating the harmonic frequencies of the electrons around their equilibrium positions, we also find the first-order vibrational corrections to the Thomson energy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Chaos control for the output-constrained system by using adaptive dynamic surface technology and application to the brushless DC motor

In this paper, chaos control is proposed for the output- constrained system with uncertain control gain and time delay and is applied to the brushless DC motor. Using the dynamic surface technology, the controller overcomes the repetitive differentiation of backstepping and boundedness hypothesis of pre-determined control gain by incorporating radial basis function neural network and adaptive technology. The tangent barrier Lyapunov function is employed for time-delay chaotic system to prevent constraint violation. It is proved that the proposed control approach can guarantee asymptotically stable in the sense of uniformly ultimate boundedness without constraint violation. Finally, the effectiveness of the proposed approach is demonstrated on the brushless DC motor example.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Communication: An efficient approach to compute state-specific nuclear gradients for a generic state-averaged multi-configuration self consistent field wavefunction

We present a new, very efficient semi-numerical approach for the computation of state-specific nuclear gradients of a generic state-averaged multi-configuration self consistent field wavefunction. Our approach eliminates the costly coupled-perturbed multi-configuration Hartree-Fock step as well as the associated integral transformation stage. The details of the implementation within the Firefly quantum chemistry package are discussed and several sample applications are given. The new approach is routinely applicable to geometry optimization of molecular systems with 1000+ basis functions using a standalone multi-core workstation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Development of massively parallel quantum chemistry program SMASH

A massively parallel program for quantum chemistry calculations SMASH was released under the Apache License 2.0 in September 2014. The SMASH program is written in the Fortran90/95 language with MPI and OpenMP standards for parallelization. Frequently used routines, such as one- and two-electron integral calculations, are modularized to make program developments simple. The speed-up of the B3LYP energy calculation for (C{sub 150}H{sub 30}){sub 2} with the cc-pVDZ basis set (4500 basis functions) was 50,499 on 98,304 cores of the K computer.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Sparse maps—A systematic infrastructure for reduced-scaling electronic structure methods. II. Linear scaling domain based pair natural orbital coupled cluster theory

Domain based local pair natural orbital coupled cluster theory with single-, double-, and perturbative triple excitations (DLPNO-CCSD(T)) is a highly efficient local correlation method. It is known to be accurate and robust and can be used in a black box fashion in order to obtain coupled cluster quality total energies for large molecules with several hundred atoms. While previous implementations showed near linear scaling up to a few hundred atoms, several nonlinear scaling steps limited the applicability of the method for very large systems. In this work, these limitations are overcome and a linear scaling DLPNO-CCSD(T) method for closed shell systems is reported. The new implementation is based on the concept of sparse maps that was introduced in Part I of this series [P. Pinski, C. Riplinger, E. F. Valeev, and F. Neese, J. Chem. Phys. 143, 034108 (2015)]. Using the sparse map infrastructure, all essential computational steps (integral transformation and storage, initial guess, pair natural orbital construction, amplitude iterations, triples correction) are achieved in a linear scaling fashion. In addition, a number of additional algorithmic improvements are reported that lead to significant speedups of the method. The new, linear-scaling DLPNO-CCSD(T) implementation typically is 7 times faster than the previous implementation and consumes 4 times less disk space for large three-dimensional systems. For linear systems, the performance gains and memory savings are substantially larger. Calculations with more than 20 000 basis functions and 1000 atoms are reported in this work. In all cases, the time required for the coupled cluster step is comparable to or lower than for the preceding Hartree-Fock calculation, even if this is carried out with the efficient resolution-of-the-identity and chain-of-spheres approximations. The new implementation even reduces the error in absolute correlation energies by about a factor of two, compared to the already accurate previous implementation.

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

Communication: An effective linear-scaling atomic-orbital reformulation of the random-phase approximation using a contracted double-Laplace transformation

An atomic-orbital (AO) reformulation of the random-phase approximation (RPA) correlation energy is presented allowing to reduce the steep computational scaling to linear, so that large systems can be studied on simple desktop computers with fully numerically controlled accuracy. Our AO-RPA formulation introduces a contracted double-Laplace transform and employs the overlap-metric resolution-of-the-identity. First timings of our pilot code illustrate the reduced scaling with systems comprising up to 1262 atoms and 10 090 basis functions. .

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A reactant-coordinate-based wave packet method for full-dimensional state-to-state quantum dynamics of tetra-atomic reactions: Application to both the abstraction and exchange channels in the H + H 2 O reaction

Here, an efficient and accurate wave packet method is proposed for the calculation of the state-to-state S -matrix elements in bimolecular reactions involving four atoms. This approach propagates an initial state specific wave packet in reactant Jacobi coordinates. The projection in product channels is carried out on projection planes, which have one less degree of freedom, by transforming both the time-dependent wave packet and final product states into a set of intermediate coordinates. This reactant-coordinate-based method is more efficient than product-coordinate-based methods because it typically requires a smaller number of basis functions or grid points and allows the determination of S -matrix elements for multiple product channels from a single propagation. This method is demonstrated in calculating the ( J tot = 0) state-to-state S -matrix elements for both the abstraction and exchange channels of the H + H 2 O reaction.

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

SparseMaps—A systematic infrastructure for reduced-scaling electronic structure methods. III. Linear-scaling multireference domain-based pair natural orbital N-electron valence perturbation theory

Multi-reference (MR) electronic structure methods, such as MR configuration interaction or MR perturbation theory, can provide reliable energies and properties for many molecular phenomena like bond breaking, excited states, transition states or magnetic properties of transition metal complexes and clusters. However, owing to their inherent complexity, most MR methods are still too computationally expensive for large systems. Therefore the development of more computationally attractive MR approaches is necessary to enable routine application for large-scale chemical systems. Among the state-of-the-art MR methods, second-order N-electron valence state perturbation theory (NEVPT2) is an efficient, size-consistent, and intruder-state-free method. However, there are still two important bottlenecks in practical applications of NEVPT2 to large systems: (a) the high computational cost of NEVPT2 for large molecules, even with moderate active spaces and (b) the prohibitive cost for treating large active spaces. In this work, we address problem (a) by developing a linear scaling “partially contracted” NEVPT2 method. This development uses the idea of domain-based local pair natural orbitals (DLPNOs) to form a highly efficient algorithm. As shown previously in the framework of single-reference methods, the DLPNO concept leads to an enormous reduction in computational effort while at the same time providing high accuracy (approaching 99.9% of the correlation energy), robustness, and black-box character. In the DLPNO approach, the virtual space is spanned by pair natural orbitals that are expanded in terms of projected atomic orbitals in large orbital domains, while the inactive space is spanned by localized orbitals. The active orbitals are left untouched. Our implementation features a highly efficient “electron pair prescreening” that skips the negligible inactive pairs. The surviving pairs are treated using the partially contracted NEVPT2 formalism. A detailed comparison between the partial and strong contraction schemes is made, with conclusions that discourage the strong contraction scheme as a basis for local correlation methods due to its non-invariance with respect to rotations in the inactive and external subspaces. A minimal set of conservatively chosen truncation thresholds controls the accuracy of the method. With the default thresholds, about 99.9% of the canonical partially contracted NEVPT2 correlation energy is recovered while the crossover of the computational cost with the already very efficient canonical method occurs reasonably early; in linear chain type compounds at a chain length of around 80 atoms. Calculations are reported for systems with more than 300 atoms and 5400 basis functions.

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

Numerical evaluation of two-center integrals over Slater type orbitals

Slater Type Orbitals (STOs) which one of the types of exponential type orbitals (ETOs) are used usually as basis functions in the multicenter molecular integrals to better understand physical and chemical properties of matter. In this work, we develop algorithms for two-center overlap and two-center two-electron hybrid and Coulomb integrals which are calculated with help of translation method for STOs and some auxiliary functions by V. Magnasco’s group. We use Mathematica programming language to produce algorithms for these calculations. Numerical results for some quantum numbers are presented in the tables. Consequently, we compare our obtained numerical results with the other known literature results and other details of evaluation method are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Communication: A difference density picture for the self-consistent field ansatz

We formulate self-consistent field (SCF) theory in terms of an interaction picture where the working variable is the difference density matrix between the true system and a corresponding superposition of atomic densities. As the difference density matrix directly represents the electronic deformations inherent in chemical bonding, this “difference self-consistent field (dSCF)” picture provides a number of significant conceptual and computational advantages. We show that this allows for a stable and efficient dSCF iterative procedure with wholly single-precision Coulomb and exchange matrix builds. We also show that the dSCF iterative procedure can be performed with aggressive screening of the pair space. These approximations are tested and found to be accurate for systems with up to 1860 atoms and >10 000 basis functions, providing for immediate overall speedups of up to 70% in the heavily optimized TERACHEM SCF implementation.

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

Atomic orbital-based SOS-MP2 with tensor hypercontraction. I. GPU-based tensor construction and exploiting sparsity

We present a tensor hypercontracted (THC) scaled opposite spin second order Møller-Plesset perturbation theory (SOS-MP2) method. By using THC, we reduce the formal scaling of SOS-MP2 with respect to molecular size from quartic to cubic. We achieve further efficiency by exploiting sparsity in the atomic orbitals and using graphical processing units (GPUs) to accelerate integral construction and matrix multiplication. The practical scaling of GPU-accelerated atomic orbital-based THC-SOS-MP2 calculations is found to be N{sup 2.6} for reference data sets of water clusters and alanine polypeptides containing up to 1600 basis functions. The errors in correlation energy with respect to density-fitting-SOS-MP2 are less than 0.5 kcal/mol for all systems tested (up to 162 atoms).

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

The hydrogen tunneling splitting in malonaldehyde: A full-dimensional time-independent quantum mechanical method

The accurate time-independent quantum dynamics calculations on the ground-state tunneling splitting of malonaldehyde in full dimensionality are reported for the first time. This is achieved with an efficient method developed by us. In our method, the basis functions are customized for the hydrogen transfer process which has the effect of greatly reducing the size of the final Hamiltonian matrix, and the Lanczos method and parallel strategy are used to further overcome the memory and central processing unit time bottlenecks. The obtained ground-state tunneling splitting of 24.5 cm{sup −1} is in excellent agreement with the benchmark value of 23.8 cm{sup −1} computed with the full-dimensional, multi-configurational time-dependent Hartree approach on the same potential energy surface, and we estimate that our reported value has an uncertainty of less than 0.5 cm{sup −1}. Moreover, the role of various vibrational modes strongly coupled to the hydrogen transfer process is revealed.

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

Computation of resistive instabilities by matched asymptotic expansions

Here, we present a method for determining the linear resistive magnetohydrodynamic (MHD) stability of an axisymmetric toroidal plasma, based on the method of matched asymptotic expansions. The plasma is partitioned into a set of ideal MHD outer regions, connected through resistive MHD inner regions about singular layers where q = m/n, with m and n toroidal mode numbers, respectively, and q the safety factor. The outer regions satisfy the ideal MHD equations with zero-frequency, which are identical to the Euler-Lagrange equations for minimizing the potential energy delta W. The solutions to these equations go to infinity at the singular surfaces. The inner regions satisfy the equations of motion of resistive MHD with a finite eigenvalue, resolving the singularity. Both outer and inner regions are solved numerically by newly developed singular Galerkin methods, using specialized basis functions. These solutions are matched asymptotically, providing a complex dispersion relation which is solved for global eigenvalues and eigenfunctions in full toroidal geometry. The dispersion relation may have multiple complex unstable roots, which are found by advanced root-finding methods. These methods are much faster and more robust than the previous numerical methods. The new methods are applicable to more challenging high-pressure and strongly shaped plasma equilibria and generalizable to more realistic inner region dynamics. In the thermonuclear regime, where the outer and inner regions overlap, they are also much faster and more accurate than the straight-through methods, which treat the resistive MHD equations in the whole plasma volume.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2017

CC2 oscillator strengths within the local framework for calculating excitation energies (LoFEx)

In a recent work [P. Baudin and K. Kristensen, J. Chem. Phys. 144, 224106 (2016)], we introduced a local framework for calculating excitation energies (LoFEx), based on second-order approximated coupled cluster (CC2) linear-response theory. LoFEx is a black-box method in which a reduced excitation orbital space (XOS) is optimized to provide coupled cluster (CC) excitation energies at a reduced computational cost. In this article, we present an extension of the LoFEx algorithm to the calculation of CC2 oscillator strengths. Two different strategies are suggested, in which the size of the XOS is determined based on the excitation energy or the oscillator strength of the targeted transitions. The two strategies are applied to a set of medium-sized organic molecules in order to assess both the accuracy and the computational cost of the methods. The results show that CC2 excitation energies and oscillator strengths can be calculated at a reduced computational cost, provided that the targeted transitions are local compared to the size of the molecule. Finally, to illustrate the potential of LoFEx for large molecules, both strategies have been successfully applied to the lowest transition of the bivalirudin molecule (4255 basis functions) and compared with time-dependent density functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗