Search NASA⌕ Search

SEARCH · Search NASA

Results for “Basis function”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,063 records · Page 59

The dissociation energy of N2

The requirements for very accurate ab initio quantum chemical prediction of dissociation energies are examined using a detailed investigation of the nitrogen molecule. Although agreement with experiment to within 1 kcal/mol is not achieved even with the most elaborate multireference CI (configuration interaction) wave functions and largest basis sets currently feasible, it is possible to obtain agreement to within about 2 kcal/mol, or 1 percent of the dissociation energy. At this level it is necessary to account for core-valence correlation effects and to include up to h-type functions in the basis. The effect of i-type functions, the use of different reference configuration spaces, and basis set superposition error were also investigated. After discussing these results, the remaining sources of error in our best calculations are examined.

Almlof, Jan↗

Curcumin inhibits development and cell adhesion in Dictyostelium discoideum: Implications for YakA signaling and GST enzyme function

The molecular basis for nutraceutical properties of the polyphenol curcumin (Curcuma longa, Turmeric) is complex, affecting multiple factors that regulate cell signaling and homeostasis. Here, we report the effect of curcumin on cellular and developmental mechanisms in the eukaryotic model, Dictyostelium discoideum. Dictyostelium proliferation was inhibited in the presence of curcumin, which also suppressed the prestarvation marker, discoidin I, members of the yakA-mediated developmental signaling pathway, and expression of the extracellular matrix/cell adhesion proteins (DdCAD and csA). This resulted in delayed chemotaxis, adhesion, and development of the organism. In contrast to the inhibitory effects on developmental genes, curcumin induced gstA gene expression, overall GST activity, and generated production of reactive oxygen species. These studies expand our knowledge of developmental and biochemical signaling influenced by curcumin, and lends greater consideration of GST enzyme function in eukaryotic cell signaling, development, and differentiation.

60 APPLIED LIFE SCIENCES↗

An Estimate of the Size and Shape of Sunspot Cycle 24 Based on its Early Cycle Behavior using the Hathaway-Wilson-Reichmann Shape-Fitting Function

On the basis of 12-month moving averages (12-mma) of monthly mean sunspot number (R), sunspot cycle 24 had its minimum amplitude (Rm = 1.7) in December 2008. At 12 mo past minimum, R measured 8.3, and at 18 mo past minimum, it measured 16.4. Thus far, the maximum month-to-month rate of rise in 12-mma values of monthly mean sunspot number (AR(t) max) has been 1.7, having occurred at elapsed times past minimum amplitude (t) of 14 and 15 mo. Compared to other sunspot cycles of the modern era, cycle 24?s Rm and AR(t) max (as observed so far) are the smallest on record, suggesting that it likely will be a slow-rising, long-period sunspot cycle of below average maximum amplitude (RM). Supporting this view is the now observed relative strength of cycle 24?s geomagnetic minimum amplitude as measured using the 12-mma value of the aa-geomagnetic index (aam = 8.4), which also is the smallest on record, having occurred at t equals 8 and 9 mo. From the method of Ohl (the inferred preferential association between RM and aam), one predicts RM = 55 +/- 17 (the ?1 se prediction interval) for cycle 24. Furthermore, from the Waldmeier effect (the inferred preferential association between the ascent duration (ASC) and RM) one predicts an ASC longer than 48 mo for cycle 24; hence, maximum amplitude occurrence should be after December 2012. Application of the Hathaway-Wilson-Reichmann shape-fitting function, using an RM = 70 and ASC = 56 mo, is found to adequately fit the early sunspot number growth of cycle 24.

Wilson, Robert M.↗

Communication: Practical and rigorous reduction of the many-electron quantum mechanical Coulomb problem to O(N{sup 2/3}) storage

It is tacitly accepted that, for practical basis sets consisting of N functions, solution of the two-electron Coulomb problem in quantum mechanics requires storage of O(N{sup 4}) integrals in the small N limit. For localized functions, in the large N limit, or for planewaves, due to closure, the storage can be reduced to O(N{sup 2}) integrals. Here, it is shown that the storage can be further reduced to O(N{sup 2/3}) for separable basis functions. A practical algorithm, that uses standard one-dimensional Gaussian-quadrature sums, is demonstrated. The resulting algorithm allows for the simultaneous storage, or fast reconstruction, of any two-electron Coulomb integral required for a many-electron calculation on processors with limited memory and disk space. For example, for calculations involving a basis of 9171 planewaves, the memory required to effectively store all Coulomb integrals decreases from 2.8 Gbytes to less than 2.4 Mbytes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Improving scalability of electronic structure code for molecular simulations in the presence of environment

A scalable density functional electronic code with Gaussian basis set, called UTEP-NRLMOL, is developed to perform simulations of molecular systems in the presence of the environment with particular attention to the memory requirements. In the electronic structure calculations, the memory and computation time are proportional to the number of atoms. Memory requirements for density functional calculations scale as N*N, where N is the number of atoms. While the recent advances in HPC offer platforms with large numbers of cores, the limited amount of memory available on a given node and poor scalability of the electronic structure codes hinder their efficient usage of these platforms. We have introduced new scaling and parallelization paradigms using MPI-3 shared-memory functionality combined with usage of sparse algebra and storage of matrices in sparse format. This extends the range of applicability of the UTEP-NRLMOL code to large systems over 10,000 atoms, or using up to 67,000 basis functions, and making use of HPC architectures using over 6,000 processors utilizing all available cores. We have also interfaced code with effective fragment potential and polarizable continuum model libraries. The code was used in simulations of several applications which are published in reputed scientific journals.

74 ATOMIC AND MOLECULAR PHYSICS↗

The study of basis sets for the calculation of the structure and dynamics of the benzene-Kr complex

An ab initio intermolecular potential energy surface (PES) has been constructed for the benzene-krypton (BKr) van der Waals (vdW) complex. The interaction energy has been calculated at the coupled cluster level of theory with single, double, and perturbatively included triple excitations using different basis sets. As a result, a few analytical PESs of the complex have been determined. They allowed a prediction of the complex structure and its vibrational vdW states. The vibrational energy level pattern exhibits a distinct polyad structure. Comparison of the equilibrium structure, the dipole moment, and vibrational levels of BKr with their experimental counterparts has allowed us to design an optimal basis set composed of a small Dunning’s basis set for the benzene monomer, a larger effective core potential adapted basis set for Kr and additional midbond functions. Such a basis set yields vibrational energy levels that agree very well with the experimental ones as well as with those calculated from the available empirical PES derived from the microwave spectra of the BKr complex. The basis proposed can be applied to larger complexes including Kr because of a reasonable computational cost and accurate results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Density functional theory study of the interaction of vinyl radical, ethyne, and ethene with benzene, aimed to define an affordable computational level to investigate stability trends in large van der Waals complexes

Our purpose is to identify a computational level sufficiently dependable and affordable to assess trends in the interaction of a variety of radical or closed shell unsaturated hydro-carbons A adsorbed on soot platelet models B. These systems, of environmental interest, would unavoidably have rather large sizes, thus prompting to explore in this paper the performances of relatively low-level computational methods and compare them with higher-level reference results. To this end, the interaction of three complexes between non-polar species, vinyl radical, ethyne, or ethene (A) with benzene (B) is studied, since these species, involved themselves in growth processes of polycyclic aromatic hydrocarbons (PAHs) and soot particles, are small enough to allow high-level reference calculations of the interaction energy ΔE{sub AB}. Counterpoise-corrected interaction energies ΔE{sub AB} are used at all stages. (1) Density Functional Theory (DFT) unconstrained optimizations of the A−B complexes are carried out, using the B3LYP-D, ωB97X-D, and M06-2X functionals, with six basis sets: 6-31G(d), 6-311 (2d,p), and 6-311++G(3df,3pd); aug-cc-pVDZ and aug-cc-pVTZ; N07T. (2) Then, unconstrained optimizations by Møller-Plesset second order Perturbation Theory (MP2), with each basis set, allow subsequent single point Coupled Cluster Singles Doubles and perturbative estimate of the Triples energy computations with the same basis sets [CCSD(T)//MP2]. (3) Based on an additivity assumption of (i) the estimated MP2 energy at the complete basis set limit [E{sub MP2/CBS}] and (ii) the higher-order correlation energy effects in passing from MP2 to CCSD(T) at the aug-cc-pVTZ basis set, ΔE{sub CC-MP}, a CCSD(T)/CBS estimate is obtained and taken as a computational energy reference. At DFT, variations in ΔE{sub AB} with basis set are not large for the title molecules, and the three functionals perform rather satisfactorily even with rather small basis sets [6-31G(d) and N07T], exhibiting deviation from the computational reference of less than 1 kcal mol{sup −1}. The zero-point vibrational energy corrected estimates Δ(E{sub AB}+ZPE), obtained with the three functionals and the 6-31G(d) and N07T basis sets, are compared with experimental D{sub 0} measures, when available. In particular, this comparison is finally extended to the naphthalene and coronene dimers and to three π−π associations of different PAHs (R, made by 10, 16, or 24 C atoms) and P (80 C atoms)

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The static response function in Kohn-Sham theory: An appropriate basis for its matrix representation in case of finite AO basis sets

The role of the static Kohn-Sham (KS) response function describing the response of the electron density to a change of the local KS potential is discussed in both the theory of the optimized effective potential (OEP) and the so-called inverse Kohn-Sham problem involving the task to find the local KS potential for a given electron density. In a general discussion of the integral equation to be solved in both cases, it is argued that a unique solution of this equation can be found even in case of finite atomic orbital basis sets. It is shown how a matrix representation of the response function can be obtained if the exchange-correlation potential is expanded in terms of a Schmidt-orthogonalized basis comprising orbitals products of occupied and virtual orbitals. The viability of this approach in both OEP theory and the inverse KS problem is illustrated by numerical examples.

97 MATHEMATICS AND COMPUTING↗

Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets

We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (~1.2 × 10 -3 E h / a 0 ) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A comprehensive analysis of molecule-intrinsic quasi-atomic, bonding, and correlating orbitals. I. Hartree-Fock wave functions

Through a basis-set-independent web of localizing orbital-transformations, the electronic wave function of a molecule is expressed in terms of a set of orbitals that reveal the atomic structure and the bonding pattern of a molecule. The analysis is based on resolving the valence orbital space in terms of an internal space, which has minimal basis set dimensions, and an external space. In the internal space, oriented quasi-atomic orbitals and split-localized molecular orbitals are determined by new, fast localization methods. The density matrix between the oriented quasi-atomic orbitals as well as the locations of the split-localized orbitals exhibit atomic populations and inter-atomic bonding patterns. A correlation-adapted quasi-atomic basis is determined in the external orbital space. The general formulations are specified in detail for Hartree-Fock wave functions. Applications to specific molecules exemplify the general scheme.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Theoretical Benchmark of the Geometric and Optical Properties for 3d Transition Metal Nanoclusters via Density Functional Theory

Understanding structure–property relationships in atomically precise metal nanoclusters is vital in finding selective and tunable catalysts. Here, in this study, density functional theory (DFT) was used to benchmark seven exchange correlation functionals at different basis sets for 17 atomically precise nanoclusters against experimentally determined geometries, band gaps, and optical gaps. The set contains both monometallic and bimetallic clusters that possess at least two types of 3d transition metals (specifically, Cu, Ni, Fe, or Co). The benchmark highlights that PBE0 is a good functional to use regardless of the basis set, and Minnesota functionals do well with respect to specific metals. Further, while long-range corrected functionals overestimate band and optical gaps, they model absorption features better than the other considered functionals. The study additionally looks at the photoinduced hydrogen evolution reaction (HER) and the CO 2 reduction mechanism on nanoclusters reported from the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Ab initio structural and vibrational properties of GaAs diamondoids and nanocrystals

Gallium arsenide diamondoids structural and vibrational properties are investigated using density functional theory at the PBE/6-31(d) level and basis including polarization functions. Variation of energy gap as these diamondoids increase in size is seen to follow confinement theory for diamondoids having nearly equiaxed dimensions. Density of energy states transforms from nearly single levels to band structure as we reach larger diamondoids. Bonds of surface hydrogen with As atoms are relatively localized and shorter than that bonded to Ga atoms. Ga-As bonds have a distribution range of values due to surface reconstruction and effect of bonding to hydrogen atoms. Experimental bulk Ga-As bond length (2.45 Å) is within this distribution range. Tetrahedral and dihedral angles approach values of bulk as we go to higher diamondoids. Optical-phonon energy of larger diamondoids stabilizes at 0.037 eV (297 cm{sup -1}) compared to experimental 0.035 eV (285.2 cm{sup -1}). Ga-As force constant reaches 1.7 mDyne/Å which is comparable to Ga-Ge force constant (1.74 mDyne/Å). Hydrogen related vibrations are nearly constant and serve as a fingerprint of GaAs diamondoids while Ga-As vibrations vary with size of diamondoids.

36 MATERIALS SCIENCE↗

Exploring metastable states in UO 2 using hybrid functionals and dynamical mean field theory

A detailed exploration of the f-atomic orbital occupancy space for UO2 is performed using a first principles approach based on density functional theory (DFT), employing a full hybrid functional within a systematic basis set. Specifically, the PBE0 functional is combined with an occupancy biasing scheme implemented in a wavelet-based algorithm which is adapted to large supercells. The results are compared with previous DFT + U calculations reported in the literature, while dynamical mean field theory is also performed to provide a further base for comparison. This work shows that the computational complexity of the energy landscape of a correlated f-electron oxide is much richer than has previously been demonstrated. The resulting calculations provide evidence of the existence of multiple previously unexplored metastable electronic states of UO 2 , including those with energies which are lower than previously reported ground states.

36 MATERIALS SCIENCE↗

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree–Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Upadhyay, Shiv (ORCID:0000000285010501)↗

Asymptotic method for determining the amplitude for three-particle breakup: Neutron-deuteron scattering

The process of neutron-deuteron scattering at energies above the deuteron-breakup threshold is described within the three-body formalism of Faddeev equations. Use is made of the method of solving Faddeev equations in configuration space on the basis of expanding wave-function components in the asymptotic region in bases of eigenfunctions of specially chosen operators. Asymptotically, wave-function components are represented in the form of an expansion in an orthonormalized basis of functions depending on the hyperangle. This basis makes it possible to orthogonalize the contributions of elastic-scattering and breakup channels. The proposed method permits determining scattering and breakup parameters from the asymptotic representation of the wave function without reconstructing it over the entire configuration space. The scattering and breakup amplitudes for states of total spin S = 1/2 and 3/2 were obtained for the s-wave Faddeev equation.

74 ATOMIC AND MOLECULAR PHYSICS↗

A Density Functional Tight Binding Model with an Extended Basis Set and Three-Body Repulsion for Hydrogen under Extreme Thermodynamic Conditions

Here, we present a new DFTB-p3b density functional tight binding model for hydrogen at extremely high pressures and temperatures, which includes a polarizable basis set ( p ) and a three-body environmentally dependent repulsive potential ( 3b ). We find that use of an extended basis set is necessary under dissociated liquid conditions to account for the substantial p-orbital character of the electronic states around the Fermi energy. The repulsive energy is determined through comparison to cold curve pressures computed from density functional theory (DFT) for the hexagonal close-packed solid, as well as pressures from thermally equilibrated DFT-MD simulations of the liquid phase. In particular, we observe improved agreement in our DFTB-p3b model with previous theoretical and experimental results for the shock Hugoniot of hydrogen up to 100 GPa and 25000 K, compared to a standard DFTB model using pairwise interactions and an s-orbital basis set, only. In conclusion, the DFTB-p3b approach discussed here provides a general method to extend the DFTB method for a wide variety of materials over a significantly larger range of thermodynamic conditions than previously possible.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterizing visible and invisible cell wall mutant phenotypes

About 10% of a plant's genome is devoted to generating the protein machinery to synthesize, remodel, and deconstruct the cell wall. High-throughput genome sequencing technologies have enabled a reasonably complete inventory of wall-related genes that can be assembled into families of common evolutionary origin. Assigning function to each gene family member has been aided immensely by identification of mutants with visible phenotypes or by chemical and spectroscopic analysis of mutants with ‘invisible’ phenotypes of modified cell wall composition and architecture that do not otherwise affect plant growth or development. This review connects the inference of gene function on the basis of deviation from the wild type in genetic functional analyses to insights provided by modern analytical techniques that have brought us ever closer to elucidating the sequence structures of the major polysaccharide components of the plant cell wall.

59 BASIC BIOLOGICAL SCIENCES↗

Hyperspherical Coulomb spheroidal basis in the Coulomb three-body problem

A hyperspherical Coulomb spheroidal (HSCS) representation is proposed for the Coulomb three-body problem. This is a new expansion in the set of well-known Coulomb spheroidal functions. The orthogonality of Coulomb spheroidal functions on a constant-hyperradius surface {rho} = const rather than on a constant-internuclear-distance surface R = const, as in the traditional Born-Oppenheimer approach, is a distinguishing feature of the proposed approach. Owing to this, the HSCS representation proves to be consistent with the asymptotic conditions for the scattering problem at energies below the threshold for three-body breakup: only a finite number of radial functions do not vanish in the limit of {rho}{yields}{infinity}, with the result that the formulation of the scattering problem becomes substantially simpler. In the proposed approach, the HSCS basis functions are considerably simpler than those in the well-known adiabatic hyperspherical representation, which is also consistent with the asymptotic conditions. Specifically, the HSCS basis functions are completely factorized. Therefore, there arise no problems associated with avoided crossings of adiabatic hyperspherical terms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗