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,117 records · Page 62
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Helical Organic and Inorganic Polymers

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

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

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

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

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

Momentum-space Argonne V18 interaction

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

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2015

Basis convergence of range-separated density-functional theory

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

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Structural insights into E1 recognition and the ubiquitin-conjugating activity of the E2 enzyme Cdc34

Ubiquitin (Ub) signaling requires the sequential interactions and activities of three enzymes, E1, E2, and E3. Cdc34 is an E2 that plays a key role in regulating cell cycle progression and requires unique structural elements to function. The molecular basis by which Cdc34 engages its E1 and the structural mechanisms by which its unique C-terminal extension functions in Cdc34 activity are unknown. Here, we present crystal structures of Cdc34 alone and in complex with E1, and a Cdc34~Ub thioester mimetic that represents the product of Uba1-Cdc34 Ub transthiolation. These structures reveal conformational changes in Uba1 and Cdc34 and a unique binding mode that are required for transthiolation. The Cdc34~Ub structure reveals contacts between the Cdc34 C-terminal extension and Ub that stabilize Cdc34~Ub in a closed conformation and are critical for Ub discharge. Altogether, our structural, biochemical, and cell-based studies provide insights into the molecular mechanisms by which Cdc34 function in cells.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2015

Vibrationally resolved NEXAFS at C and N K-edges of pyridine, 2-fluoropyridine and 2,6-difluoropyridine: A combined experimental and theoretical assessment

In the present work, the near edge X-ray absorption spectroscopy (NEXAFS) spectra at both C and N K-edges of pyridine, 2-fluoropyridine, and 2,6-difluoropyridine have been studied both experimentally and theoretically. From an electronic point of view, both transition potential density functional theory and time-dependent density functional theory approaches lead to reliable results provided that suitable basis sets and density functionals are employed. In this connection, the global hybrid B3LYP functional in conjunction with the EPR-III basis set appears particularly suitable after constant scaling of the band positions. For the N K-edge, vertical energies obtained at these levels and broadened by symmetric Gaussian distributions provide spectra in reasonable agreement with the experiment. Vibronic contributions further modulate the band-shapes leading to a better agreement with the experimental results, but are not strictly necessary for semi-quantitative investigations. On the other hand, vibronic contributions are responsible for strong intensity redistribution in the NEXAFS C K-edge spectra, and their inclusion is thus mandatory for a proper description of experiments. In this connection, the simple vertical gradient model is particularly appealing in view of its sufficient reliability and low computational cost. For more quantitative results, the more refined vertical Hessian approach can be employed, and its effectiveness has been improved thanks to a new least-squares fitting approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Accurate ab initio calculations which demonstrate a 3 Pi u ground state for Al2

The spectroscopic parameters and separations between the three low-lying X 3 Pi u, A 3 Sigma g -, and a 1 Sigma g + states of Al2 are studied as a function of both the one-particle and n-particle basis set. Approximate correlation treatments are calibrated against full Cl calculations correlating the six valence electrons in a double-zeta plus two d-function basis set. Since the CASSCF/MRCI 3 Pi u to 3 Sigma g - separation is in excellent agreement wtih the FCI value, the MRCI calculations were carried out in an extended (20s13p6d4f)/(6s5p3d2f) gaussian basis. Including a small correction for relativistic effects, the best estimate is that 3 Sigma g - state lies 174/cm above the 3 Pi u ground state. The 1 Sigma g + state lies at least 2000/cm higher in energy. At the CPF level, inclusion of 2s and 2p correlation has little effect on D sub e, reduces T sub e by only 26/cm, and shortens the bond lengths by about 0.02 a sub o. Further strong support for a 3 Pi u ground state comes from the experimental absorption spectra, since both observed transitions can be convincingly assigned as 3 Pi u yields 3 Pi g. The (2) 3 Pi g state is observed to be sensitive to the level of correlation treatment, and to have its minimum shifted to shorter rho values, such that the strongest experimental absorption peak probably corresponds to the 0 yields 2 transition.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2014

An efficient extrapolation to the (T)/CBS limit

We extrapolate to the perturbative triples (T)/complete basis set (CBS) limit using double ζ basis sets without polarization functions (Wesleyan-1-Triples-2ζ or “Wes1T-2Z”) and triple ζ basis sets with a single level of polarization functions (Wesleyan-1-Triples-3ζ or “Wes1T-3Z”). These basis sets were optimized for 102 species representing the first two rows of the Periodic Table. The species include the entire set of neutral atoms, positive and negative atomic ions, as well as several homonuclear diatomic molecules, hydrides, rare gas dimers, polar molecules, such as oxides and fluorides, and a few transition states. The extrapolated Wes1T-(2,3)Z triples energies agree with (T)/CBS benchmarks to within ±0.65 mE{sub h}, while the rms deviations of comparable model chemistries W1, CBS-APNO, and CBS-QB3 for the same test set are ±0.23 mE{sub h}, ±2.37 mE{sub h}, and ±5.80 mE{sub h}, respectively. The Wes1T-(2,3)Z triples calculation time for the largest hydrocarbon in the G2/97 test set, C{sub 6}H{sub 5}Me{sup +}, is reduced by a factor of 25 when compared to W1. The cost-effectiveness of the Wes1T-(2,3)Z extrapolation validates the usefulness of the Wes1T-2Z and Wes1T-3Z basis sets which are now available for a more efficient extrapolation of the (T) component of any composite model chemistry.

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

An efficient basis set representation for calculating electrons in molecules

The method of McCurdy, Baertschy, and Rescigno, is generalised to obtain a straightforward, surprisingly accurate, and scalable numerical representation for calculating the electronic wave functions of molecules. It uses a basis set of product sinc functions arrayed on a Cartesian grid, and yields 1 kcal/mol precision for valence transition energies with a grid resolution of approximately 0.1 bohr. The Coulomb matrix elements are replaced with matrix elements obtained from the kinetic energy operator. A resolution-of-the-identity approximation renders the primitive one- and two-electron matrix elements diagonal; in other words, the Coulomb operator is local with respect to the grid indices. The calculation of contracted two-electron matrix elements among orbitals requires only O( N log (N)) multiplication operations, not O( N 4 ), where N is the number of basis functions; N = n 3 on cubic grids. The representation not only is numerically expedient, but also produces energies and properties superior to those calculated variationally. Absolute energies, absorption cross sections, transition energies, and ionisation potentials are reported for 1- (He + , H + 2 ), 2- (H 2 , He), 10- (CH 4 ), and 56-electron (C 8 H 8 ) systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

The importance of diffuse f functions for transition metals

The importance of a diffuse f function for determining the dissociation energy (De) of Cu2 and the 3F-3D energy separation in Ni atom is investigated. It is found that the diffuse f contributes at most 0.05 eV to the De of Cu2 when added to a basis containing tight f functions and a flexibly contracted d basis. The diffuse f function is found to decrease the 3F-3D separation in Ni, but by substantially less than the tight f functions.

Langhoff, S. R.↗
DOE OSTI2019

First-principles study of the infrared spectrum in liquid water from a systematically improved description of H-bond network

An accurate ab initio theory of the H-bond structure of liquid water requires a high-level exchange correlation approximation from density functional theory. Based on the liquid structures modeled by ab initio molecular dynamics by using maximally localized Wannier functions as a basis, we study the infrared spectrum of water within the canonical ensemble. In particular, we employ both the Perdew-Burke-Ernzerhof (PBE) functional within the generalized gradient approximation (GGA) and the state-of-the-art meta-GGA level approximation provided by the strongly constrained and appropriately normed (SCAN) functional. In this work, we demonstrate that the SCAN functional improves not only the water structure but also the theoretical infrared spectrum of water. Our analyses show that the improvement in the stretching and bending bands can be mainly attributed to better descriptions of directional H bonding and the covalency at the inter- and intramolecular levels, respectively. On the other hand, better agreements in libration and hindered translation bands are due to the improved dynamics of the H-bond network enabled by a less structured liquid in the experimental direction. The spectrum predicted by SCAN shows much better agreement with experimental data than the conventionally widely adopted PBE functional at the GGA level.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Nuclear-electronic orbital reduced explicitly correlated Hartree-Fock approach: Restricted basis sets and open-shell systems

The nuclear electronic orbital (NEO) reduced explicitly correlated Hartree-Fock (RXCHF) approach couples select electronic orbitals to the nuclear orbital via Gaussian-type geminal functions. This approach is extended to enable the use of a restricted basis set for the explicitly correlated electronic orbitals and an open-shell treatment for the other electronic orbitals. The working equations are derived and the implementation is discussed for both extensions. The RXCHF method with a restricted basis set is applied to HCN and FHF{sup −} and is shown to agree quantitatively with results from RXCHF calculations with a full basis set. The number of many-particle integrals that must be calculated for these two molecules is reduced by over an order of magnitude with essentially no loss in accuracy, and the reduction factor will increase substantially for larger systems. Typically, the computational cost of RXCHF calculations with restricted basis sets will scale in terms of the number of basis functions centered on the quantum nucleus and the covalently bonded neighbor(s). In addition, the RXCHF method with an odd number of electrons that are not explicitly correlated to the nuclear orbital is implemented using a restricted open-shell formalism for these electrons. This method is applied to HCN{sup +}, and the nuclear densities are in qualitative agreement with grid-based calculations. Future work will focus on the significance of nonadiabatic effects in molecular systems and the further enhancement of the NEO-RXCHF approach to accurately describe such effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

STEP and STEPSPL: Computer programs for aerodynamic model structure determination and parameter estimation

The successful parametric modeling of the aerodynamics for an airplane operating at high angles of attack or sideslip is performed in two phases. First the aerodynamic model structure must be determined and second the associated aerodynamic parameters (stability and control derivatives) must be estimated for that model. The purpose of this paper is to document two versions of a stepwise regression computer program which were developed for the determination of airplane aerodynamic model structure and to provide two examples of their use on computer generated data. References are provided for the application of the programs to real flight data. The two computer programs that are the subject of this report, STEP and STEPSPL, are written in FORTRAN IV (ANSI l966) compatible with a CDC FTN4 compiler. Both programs are adaptations of a standard forward stepwise regression algorithm. The purpose of the adaptation is to facilitate the selection of a adequate mathematical model of the aerodynamic force and moment coefficients of an airplane from flight test data. The major difference between STEP and STEPSPL is in the basis for the model. The basis for the model in STEP is the standard polynomial Taylor's series expansion of the aerodynamic function about some steady-state trim condition. Program STEPSPL utilizes a set of spline basis functions.

Batterson, J. G.↗
DOE OSTI2018

Modeling the diffusion of heat energy within composites of homogeneous materials using the uncertainty principle

The goal of this paper is to develop a highly accurate and efficient numerical method for the solution of a time-dependent partial differential equation with a piecewise constant coefficient, on a finite interval with periodic boundary conditions. The resulting algorithm can be used, for example, to model the diffusion of heat energy in one space dimension, in the case where the spatial domain represents a medium consisting of two homogeneous materials. The resulting model has, to our knowledge, not yet been solved in closed form through analytical methods, and is difficult to solve using existing numerical methods, thus suggesting an alternative approach. The approach presented in this paper is to represent the solution as a linear combination of wave functions that change frequencies at the interfaces between different materials. It is demonstrated through numerical experiments that using the Uncertainty Principle to construct a basis of such functions, in conjunction with a spectral method, a mathematical model for heat diffusion through different materials can be solved much more efficiently than with conventional time-stepping methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗
DOE OSTI2024

Specialization Restricts the Evolutionary Paths Available to Yeast Sugar Transporters

Functional innovation at the protein level is a key source of evolutionary novelties. The constraints on functional innovations are likely to be highly specific in different proteins, which are shaped by their unique histories and the extent of global epistasis that arises from their structures and biochemistries. These contextual nuances in the sequence–function relationship have implications both for a basic understanding of the evolutionary process and for engineering proteins with desirable properties. Here, we have investigated the molecular basis of novel function in a model member of an ancient, conserved, and biotechnologically relevant protein family. These Major Facilitator Superfamily sugar porters are a functionally diverse group of proteins that are thought to be highly plastic and evolvable. By dissecting a recent evolutionary innovation in an α-glucoside transporter from the yeast Saccharomyces eubayanus, we show that the ability to transport a novel substrate requires high-order interactions between many protein regions and numerous specific residues proximal to the transport channel. To reconcile the functional diversity of this family with the constrained evolution of this model protein, we generated new, state-of-the-art genome annotations for 332 Saccharomycotina yeast species spanning ~400 My of evolution. By integrating phylogenetic and phenotypic analyses across these species, we show that the model yeast α-glucoside transporters likely evolved from a multifunctional ancestor and became subfunctionalized. The accumulation of additive and epistatic substitutions likely entrenched this subfunction, which made the simultaneous acquisition of multiple interacting substitutions the only reasonably accessible path to novelty.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Reproducibility of calculations on Li species with correlation-consistent basis sets

Here we have found that two different basis set series, both labeled cc-pVxZ, are widely available and widely used for calculation on Li compounds such that calculations that should be the same are different, causing a reproducibility problem. Here we study the magnitude of the disagreement of these nominally-the-same basis sets for wave function and density functional calculations of Li clusters. We find deviations among the basis sets of as large as 22 millihartrees. In most cases one obtains lower energies with the noncanonical choices. Therefore, it is important to specify which version of the basis set is used for calculations.

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

HVEM structures and mutants reveal distinct functions of binding to LIGHT and BTLA/CD160

HVEM is a TNF (tumor necrosis factor) receptor contributing to a broad range of immune functions involving diverse cell types. It interacts with a TNF ligand, LIGHT, and immunoglobulin (Ig) superfamily members BTLA and CD160. Assessing the functional impact of HVEM binding to specific ligands in different settings has been complicated by the multiple interactions of HVEM and HVEM binding partners. To dissect the molecular basis for multiple functions, we determined crystal structures that reveal the distinct HVEM surfaces that engage LIGHT or BTLA/CD160, including the human HVEM–LIGHT–CD160 ternary complex, with HVEM interacting simultaneously with both binding partners. Based on these structures, we generated mouse HVEM mutants that selectively recognized either the TNF or Ig ligands in vitro. Knockin mice expressing these muteins maintain expression of all the proteins in the HVEM network, yet they demonstrate selective functions for LIGHT in the clearance of bacteria in the intestine and for the Ig ligands in the amelioration of liver inflammation.

60 APPLIED LIFE SCIENCES↗