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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

What Is the Fractional Laplacian? A Comparative Review with New Results

The fractional Laplacian in $\mathbb{R}^{d}$, which we write as (–Δ) α/2 with, α ϵ (0, 2) has multiple equivalent characterizations. Furthermore, in bounded domains, boundary conditions must be incorporated in these characterizations in mathematically distinct ways, and there is currently no consensus in the literature as to which definition of the fractional Laplacian in bounded domains is most appropriate for a given application. The Riesz (or integral) definition, for example, admits a nonlocal boundary condition, where the value of a function must be prescribed on the entire exterior of the domain in order to compute its fractional Laplacian. Yet, the spectral definition requires only the standard local boundary condition. These differences, among others, lead us to ask the question: “What is the fractional Laplacian?” Beginning from first principles, we compare several commonly used definitions of the fractional Laplacian theoretically, through their stochastic interpretations as well as their analytical properties. Next, we present quantitative comparisons using a sample of state-of-the-art methods. We discuss recent advances on nonzero boundary conditions and present new methods to discretize such boundary value problems: radial basis function collocation (for the Riesz fractional Laplacian) and nonharmonic lifting (for the spectral fractional Laplacian). In our numerical studies, we aim to compare different definitions on bounded domains using a collection of benchmark problems. We consider the fractional Poisson equation with both zero and nonzero boundary conditions, where the fractional Laplacian is defined according to the Riesz definition, the spectral definition, the directional definition, and the horizon-based nonlocal definition. We verify the accuracy of the numerical methods used in the approximations for each operator, and we focus on identifying differences in the boundary behaviors of solutions to equations posed with these different definitions. Through our efforts, we aim to further engage the research community in open problems and assist practitioners in identifying the most appropriate definition and computational approach to use for their mathematical models in addressing anomalous transport in diverse applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

Analysis and comparison of CVS-ADC approaches up to third order for the calculation of core-excited states

The extended second order algebraic-diagrammatic construction (ADC(2)-x) scheme for the polarization operator in combination with core-valence separation (CVS) approximation is well known to be a powerful quantum chemical method for the calculation of core-excited states and the description of X-ray absorption spectra. For the first time, the implementation and results of the third order approach CVS-ADC(3) are reported. Therefore, the CVS approximation has been applied to the ADC(3) working equations and the resulting terms have been implemented efficiently in the adcman program. By treating the α and β spins separately from each other, the unrestricted variant CVS-UADC(3) for the treatment of open-shell systems has been implemented as well. The performance and accuracy of the CVS-ADC(3) method are demonstrated with respect to a set of small and middle-sized organic molecules. Therefore, the results obtained at the CVS-ADC(3) level are compared with CVS-ADC(2)-x values as well as experimental data by calculating complete basis set limits. The influence of basis sets is further investigated by employing a large set of different basis sets. Besides the accuracy of core-excitation energies and oscillator strengths, the importance of cartesian basis functions and the treatment of orbital relaxation effects are analyzed in this work as well as computational timings. It turns out that at the CVS-ADC(3) level, the results are not further improved compared to CVS-ADC(2)-x and experimental data, because the fortuitous error compensation inherent in the CVS-ADC(2)-x approach is broken. While CVS-ADC(3) overestimates the core excitation energies on average by 0.61% ± 0.31%, CVS-ADC(2)-x provides an averaged underestimation of −0.22% ± 0.12%. Eventually, the best agreement with experiments can be achieved using the CVS-ADC(2)-x method in combination with a diffuse cartesian basis set at least at the triple-ζ level.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Porting Fragmentation Methods to Graphical Processing Units Using an OpenMP Application Programming Interface: Offloading the Fock Build for Low Angular Momentum Functions

Here, a framework to offload four-index two-electron repulsion integrals to graphical processing units (GPUs) using OpenMP is discussed. The method has been applied to the Fock build for low angular momentum s and p functions in both the restricted Hartree–Fock (RHF) and in the effective fragment molecular orbital (EFMO) framework. Benchmark calculations for the GPU code for the pure RHF method show an increasing speedup relative to the existing OpenMP CPU code in GAMESS from 1.04 to 52× for clusters of 70–569 water molecules. The parallel efficiency on 24 NVIDIA V100 GPU boards also increases when increasing the system size: from 75 to 94% for water clusters that contain 303–1120 molecules. In the EFMO framework, the GPU Fock build shows a high linear scalability up to 4608 V100s with a parallel efficiency of 96% for calculations on a solvated mesoporous silica nanoparticle system with ~67,000 basis functions.

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

Analysis of nonlinear optical properties in donor–acceptor materials

Time-dependent density functional theory has been used to calculate nonlinear optical (NLO) properties, including the first and second hyperpolarizabilities as well as the two-photon absorption cross-section, for the donor-acceptor molecules p-nitroaniline and dimethylamino nitrostilbene, and for respective materials attached to a gold dimer. The CAMB3LYP, B3LYP, PBE0, and PBE exchange-correlation functionals all had fair but variable performance when compared to higher-level theory and to experiment. The CAMB3LYP functional had the best performance on these compounds of the functionals tested. However, our comprehensive analysis has shown that quantitative prediction of hyperpolarizabilities is still a challenge, hampered by inadequate functionals, basis sets, and solvation models, requiring further experimental characterization. Attachment of the Au{sub 2}S group to molecules already known for their relatively large NLO properties was found to further enhance the response. While our calculations show a modest enhancement for the first hyperpolarizability, the enhancement of the second hyperpolarizability is predicted to be more than an order of magnitude.

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

Electronic switching by metastable polarization states in BiFe O 3 thin films

Here, we present an approach to control resistive switching in metal-ferroelectric contacts using a radially symmetric electric field. In ferroelectrics with significant polarization along the corresponding field lines, the field above a critical threshold will induce polarization discontinuity, a corresponding nanoscale volume of space charge, and a conducting junction. We demonstrate this principle using nanoscale polarization switching of a conventional (001)-oriented thin film of BiFeO 3 . Without any optimization, the conducting state created in this regime of resistive switching exhibits local currents of ~1–10nA, approaching the ~100nA threshold required for device implementation [J. Jiang et al., Nat. Mater. 17, 49 (2018)]. The corresponding electronic function is that of a volatile resistive switch, which is directly compatible with neuristor functionality that encodes the functioning basis of an axon [M. D. Pickett et al., Nat. Mater. 12, 114 (2013)]. Phase-field modeling further reveals that in the strongly charged local configuration, BiFeO 3 locally undergoes a rhombohedral-tetragonal (R-T) phase transition, in part due to substantial piezoelectric expansion of the lattice. The estimated local charge density can be as high as ~10 21 cm -3 , which would be extremely difficult to achieve by conventional doping approaches without altering other material properties. Therefore, this method for creating stable and reproducible strongly charged ferroelectric junctions enables more systematic studies of their physical properties, such as the possibility of structural and electronic phase transitions, and it can lead to new ferroelectric devices for advanced information functions.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Calculation of {sup 15}N NMR chemical shifts of amines in the framework of the density functional theory

The {sup 15}N NMR magnetic shielding constants (followed by recalculation into chemical shifts) in a representative series of amines were calculated in the framework of the density functional theory. The results were compared with experiment. The accuracy factors of this calculation (functional, basis set, solvent effects, and multistandard) were studied. Taking into account all the above factors leads to a noticeable decrease in the average absolute error in the calculation of the {sup 15}N NMR chemical shifts (from 13 to 4 ppm) in a range of their changing in the studies series of compounds of ∼60 ppm (which is 6—7% in relative units).

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

Parsimonious Potential Energy Surface Expansions Using Dictionary Learning with Multipass Greedy Selection

Potential energy surfaces fit with basis set expansions have been shown to provide accurate representations of electronic energies and have enabled a variety of high-accuracy dynamics, kinetics, and spectroscopy applications. The number of terms in these expansions scales poorly with system size, a drawback that challenges their use for systems with more than similar to 10 atoms. A solution is presented here using dictionary learning. Subsets of the full set of conventional basis functions are optimized using a newly developed multipass greedy regression method inspired by forward and backward selection methods from the statistics, signal processing, and machine learning literatures. Here, the optimized representations have accuracies comparable to the full set but are 1 or more orders of magnitude smaller, and notably, the number of terms in the optimized multipass greedy expansions scales approximately linearly with the number of atoms.

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

Prediction of {sup 1}P Rydberg energy levels of beryllium based on calculations with explicitly correlated Gaussians

Benchmark variational calculations are performed for the seven lowest 1s{sup 2}2s np ({sup 1}P), n = 2…8, states of the beryllium atom. The calculations explicitly include the effect of finite mass of {sup 9}Be nucleus and account perturbatively for the mass-velocity, Darwin, and spin-spin relativistic corrections. The wave functions of the states are expanded in terms of all-electron explicitly correlated Gaussian functions. Basis sets of up to 12 500 optimized Gaussians are used. The maximum discrepancy between the calculated nonrelativistic and experimental energies of 1s{sup 2}2s np ({sup 1}P) →1s{sup 2}2s{sup 2} ({sup 1}S) transition is about 12 cm{sup −1}. The inclusion of the relativistic corrections reduces the discrepancy to bellow 0.8 cm{sup −1}.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A coupled discontinuous Galerkin-Finite Volume framework for solving gas dynamics over embedded geometries

Herein, we present a computational framework for solving the equations of inviscid gas dynamics using structured grids with embedded geometries. The novelty of the proposed approach is the use of high-order discontinuous Galerkin (dG) schemes and a shock-capturing Finite Volume (FV) scheme coupled via an hp adaptive mesh refinement (hp-AMR) strategy that offers high-order accurate resolution of the embedded geometries. The hp-AMR strategy is based on a multi-level block-structured domain partition in which each level is represented by block-structured Cartesian grids and the embedded geometry is represented implicitly by a level set function. The intersection of the embedded geometry with the grids produces the implicitly-defined mesh that consists of a collection of regular rectangular cells plus a relatively small number of irregular curved elements in the vicinity of the embedded boundaries. High-order quadrature rules for implicitly-defined domains enable high-order accuracy resolution of the curved elements with a cell-merging strategy to address the small-cell problem. The hp-AMR algorithm treats the system with a second-order finite volume scheme at the finest level to dynamically track the evolution of solution discontinuities while using dG schemes at coarser levels to provide high-order accuracy in smooth regions of the flow. On the dG levels, the methodology supports different orders of basis functions on different levels. The space-discretized governing equations are then advanced explicitly in time using high-order Runge-Kutta algorithms. Numerical tests are presented for two-dimensional and three-dimensional problems involving an ideal gas. The results are compared with both analytical solutions and experimental observations and demonstrate that the framework provides high-order accuracy for smooth flows and accurately captures solution discontinuities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2011

Autoionization resonance states of two-electron atomic systems with finite spherical confinement

We investigate the lowest-lying S-wave resonant states of two-electron atoms confined by a spherical quantum cavity under the framework of the stabilization method. Hylleraas-type wave functions (basis length N = 444) taking the correlation effects between all the charged particles into account are used in the present paper. The finite oscillator potential is used to represent the confinement potential. We present the resonant parameters (energies and widths) of the quantum-confined two-electron atoms with different depths and various ranges of the potentials.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Magnetizability and rotational g tensors for density fitted local second-order Møller-Plesset perturbation theory using gauge-including atomic orbitals

In this paper, we present theory and implementation of an efficient program for calculating magnetizabilities and rotational g tensors of closed-shell molecules at the level of local second-order Møller-Plesset perturbation theory (MP2) using London orbitals. Density fitting is employed to factorize the electron repulsion integrals with ordinary Gaussians as fitting functions. The presented program for the calculation of magnetizabilities and rotational g tensors is based on a previous implementation of NMR shielding tensors reported by S. Loibl and M. Schütz [J. Chem. Phys. 137, 084107 (2012)]. Extensive test calculations show (i) that the errors introduced by density fitting are negligible, and (ii) that the errors of the local approximation are still rather small, although larger than for nuclear magnetic resonance (NMR) shielding tensors. Electron correlation effects for magnetizabilities are tiny for most of the molecules considered here. MP2 appears to overestimate the correlation contribution of magnetizabilities such that it does not constitute an improvement over Hartree-Fock (when comparing to higher-order methods like CCSD(T)). For rotational g tensors the situation is different and MP2 provides a significant improvement in accuracy over Hartree-Fock. The computational performance of the new program was tested for two extended systems, the larger comprising about 2200 basis functions. It turns out that a magnetizability (or rotational g tensor) calculation takes about 1.5 times longer than a corresponding NMR shielding tensor calculation.

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

Selection of polynomial chaos bases via Bayesian model uncertainty methods with applications to sparse approximation of PDEs with stochastic inputs

Generalized polynomial chaos (gPC) expansions allow us to represent the solution of a stochastic system using a series of polynomial chaos basis functions. The number of gPC terms increases dramatically as the dimension of the random input variables increases. When the number of the gPC terms is larger than that of the available samples, a scenario that often occurs when the corresponding deterministic solver is computationally expensive, evaluation of the gPC expansion can be inaccurate due to over-fitting. We propose a fully Bayesian approach that allows for global recovery of the stochastic solutions, in both spatial and random domains, by coupling Bayesian model uncertainty and regularization regression methods. It allows the evaluation of the PC coefficients on a grid of spatial points, via (1) the Bayesian model average (BMA) or (2) the median probability model, and their construction as spatial functions on the spatial domain via spline interpolation. The former accounts for the model uncertainty and provides Bayes-optimal predictions; while the latter provides a sparse representation of the stochastic solutions by evaluating the expansion on a subset of dominating gPC bases. Moreover, the proposed methods quantify the importance of the gPC bases in the probabilistic sense through inclusion probabilities. We design a Markov chain Monte Carlo (MCMC) sampler that evaluates all the unknown quantities without the need of ad-hoc techniques. The proposed methods are suitable for, but not restricted to, problems whose stochastic solutions are sparse in the stochastic space with respect to the gPC bases while the deterministic solver involved is expensive. We demonstrate the accuracy and performance of the proposed methods and make comparisons with other approaches on solving elliptic SPDEs with 1-, 14- and 40-random dimensions.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Electronic and thermoelectric properties of InN studied using ab initio density functional theory and Boltzmann transport calculations

The thermoelectric properties of indium nitride in the most stable wurtzite phase (w-InN) as a function of electron and hole concentrations and temperature were studied by solving the semiclassical Boltzmann transport equations in conjunction with ab initio electronic structure calculations, within Density Functional Theory. Based on maximally localized Wannier function basis set and the ab initio band energies, results for the Seebeck coefficient are presented and compared with available experimental data for n-type as well as p-type systems. Also, theoretical results for electric conductivity and power factor are presented. Most cases showed good agreement between the calculated properties and experimental data for w-InN unintentionally and p-type doped with magnesium. Our predictions for temperature and concentration dependences of electrical conductivity and power factor revealed a promising use of InN for intermediate and high temperature thermoelectric applications. The rigid band approach and constant scattering time approximation were utilized in the calculations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Efficient calculation of nuclear spin-rotation constants from auxiliary density functional theory

The computation of the spin-rotation tensor within the framework of auxiliary density functional theory (ADFT) in combination with the gauge including atomic orbital (GIAO) scheme, to treat the gauge origin problem, is presented. For the spin-rotation tensor, the calculation of the magnetic shielding tensor represents the most demanding computational task. Employing the ADFT-GIAO methodology, the central processing unit time for the magnetic shielding tensor calculation can be dramatically reduced. In this work, the quality of spin-rotation constants obtained with the ADFT-GIAO methodology is compared with available experimental data as well as with other theoretical results at the Hartree-Fock and coupled-cluster level of theory. It is found that the agreement between the ADFT-GIAO results and the experiment is good and very similar to the ones obtained by the coupled-cluster single-doubles-perturbative triples-GIAO methodology. With the improved computational performance achieved, the computation of the spin-rotation tensors of large systems or along Born-Oppenheimer molecular dynamics trajectories becomes feasible in reasonable times. Three models of carbon fullerenes containing hundreds of atoms and thousands of basis functions are used for benchmarking the performance. Furthermore, a theoretical study of temperature effects on the structure and spin-rotation tensor of the H{sup 12}C–{sup 12}CH–DF complex is presented. Here, the temperature dependency of the spin-rotation tensor of the fluorine nucleus can be used to identify experimentally the so far unknown bent isomer of this complex. To the best of our knowledge this is the first time that temperature effects on the spin-rotation tensor are investigated.

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

Excited-state electronic structure of molecules using many-body Green’s functions: Quasiparticles and electron–hole excitations with VOTCA-XTP

Here, we present the open-source VOTCA-XTP software for the calculation of the excited-state electronic structure of molecules using many-body Green’s function theory in the GW approximation with the Bethe–Salpeter equation (BSE). This work provides a summary of the underlying theory and discusses the details of its implementation based on Gaussian orbitals, including resolution-of-identity techniques and different approaches to the frequency integration of the self-energy or acceleration by offloading compute-intensive matrix operations using graphics processing units in a hybrid OpenMP/Cuda scheme. A distinctive feature of VOTCA–XTP is the capability to couple the calculation of electronic excitations to a classical polarizable environment on an atomistic level in a coupled quantum- and molecular-mechanics (QM/MM) scheme, where a complex morphology can be imported from Molecular Dynamics simulations. The capabilities and limitations of the GW–BSE implementation are illustrated with two examples. First, we study the dependence of optically active electron–hole excitations in a series of diketopyrrolopyrrole-based oligomers on molecular-architecture modifications and the number of repeat units. Second, we use the GW–BSE/MM setup to investigate the effect of polarization on localized and intermolecular charge-transfer excited states in morphologies of low-donor content rubrene–fullerene mixtures. These showcases demonstrate that our implementation currently allows us to treat systems with up to 2500 basis functions on regular shared-memory workstations, providing accurate descriptions of quasiparticle and coupled electron–hole excited states of various characters on an equal footing.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Unbiased particle conformation extraction from scattering spectra using orthonormal basis expansions

A strategy is outlined for quantitatively evaluating the particle density profiles from small-angle scattering spectra of dilute solutions. The approach employs an orthonormal basis function expansion method, enabling the determination of characteristic mass distributions in self-assembled structures without the need for a specific structural model. Through computational benchmarking, the efficacy of this approach is validated by effectively reconstructing the density profile of soft-ball systems with varying fuzziness from their scattering signatures. Further, the feasibility of the method is demonstrated by fitting small-angle neutron scattering data obtained from Pluronic L64 micelles at different temperatures. This proposed approach is both simple and analytical, eliminating the requirement for a presumptive structural model in scattering analysis. The new method could therefore facilitate quantitative descriptions of complex nanoscopic structures inherent to numerous soft-matter systems using small-angle scattering techniques.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Multinode Multi-GPU Two-Electron Integrals: Code Generation Using the Regent Language

The computation of two-electron repulsion integrals (ERIs) is often the most expensive step of integral-direct self-consistent field methods. Formally it scales as O(N 4 ), where N is the number of Gaussian basis functions used to represent the molecular wave function. In practice, this scaling can be reduced to O(N 2 ) or less by neglecting small integrals with screening methods. The contributions of the ERIs to the Fock matrix are of Coulomb (J) and exchange (K) type and require separate algorithms to compute matrix elements efficiently. We previously implemented highly efficient GPU-accelerated J-matrix and K-matrix algorithms in the electronic structure code TeraChem. Although these implementations supported the use of multiple GPUs on a node, they did not support the use of multiple nodes. This presents a key bottleneck to cutting-edge ab initio simulations of large systems, e.g., excited state dynamics of photoactive proteins. We present our implementation of multinode multi-GPU J- and K-matrix algorithms in TeraChem using the Regent programming language. Regent directly supports distributed computation in a task-based model and can generate code for a variety of architectures, including NVIDIA GPUs. We demonstrate multinode scaling up to 45 GPUs (3 nodes) and benchmark against hand-coded TeraChem integral code. Finally, we also outline our metaprogrammed Regent implementation, which enables flexible code generation for integrals of different angular momenta.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗