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

Use of Nanostructure-Initiator Mass Spectrometry to Deduce Selectivity of Reaction in Glycoside Hydrolases

Chemically synthesized nanostructure-initiator mass spectrometry (NIMS) probes derivatized with tetrasaccharides were used to study the reactivity of representative Clostridium thermocellum β-glucosidase, endoglucanases, and cellobiohydrolase. Diagnostic patterns for reactions of these different classes of enzymes were observed. Results show sequential removal of glucose by the β-glucosidase and a progressive increase in specificity of reaction from endoglucanases to cellobiohydrolase. Time-dependent reactions of these polysaccharide-selective enzymes were modeled by numerical integration, which provides a quantitative basis to make functional distinctions among a continuum of naturally evolved catalytic properties. Consequently, our method, which combines automated protein translation with high-sensitivity and time-dependent detection of multiple products, provides a new approach to annotate glycoside hydrolase phylogenetic trees with functional measurements.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A Study of the Radiation Tolerance of CVD Diamond to 70 MeV Protons, Fast Neutrons and 200 MeV Pions

We measured the radiation tolerance of commercially available diamonds grown by the Chemical Vapor Deposition process by measuring the charge created by a 120 GeV hadron beam in a 50 μm pitch strip detector fabricated on each diamond sample before and after irradiation. We irradiated one group of samples with 70 MeV protons, a second group of samples with fast reactor neutrons (defined as energy greater than 0.1 MeV), and a third group of samples with 200 MeV pions, in steps, to (8.8±0.9) × 1015 protons/cm2, (1.43±0.14) × 1016 neutrons/cm2, and (6.5±1.4) × 1014 pions/cm2, respectively. By observing the charge induced due to the separation of electron–hole pairs created by the passage of the hadron beam through each sample, on an event-by-event basis, as a function of irradiation fluence, we conclude all datasets can be described by a first-order damage equation and independently calculate the damage constant for 70 MeV protons, fast reactor neutrons, and 200 MeV pions. We find the damage constant for diamond irradiated with 70 MeV protons to be 1.62±0.07(stat)±0.16(syst)× 10−18 cm2/(p μm), the damage constant for diamond irradiated with fast reactor neutrons to be 2.65±0.13(stat)±0.18(syst)× 10−18 cm2/(n μm), and the damage constant for diamond irradiated with 200 MeV pions to be 2.0±0.2(stat)±0.5(syst)× 10−18 cm2/(π μm). The damage constants from this measurement were analyzed together with our previously published 24 GeV proton irradiation and 800 MeV proton irradiation damage constant data to derive the first comprehensive set of relative damage constants for Chemical Vapor Deposition diamond. We find 70 MeV protons are 2.60 ± 0.29 times more damaging than 24 GeV protons, fast reactor neutrons are 4.3 ± 0.4 times more damaging than 24 GeV protons, and 200 MeV pions are 3.2 ± 0.8 more damaging than 24 GeV protons. We also observe the measured data can be described by a universal damage curve for all proton, neutron, and pion irradiations we performed of Chemical Vapor Deposition diamond. Finally, we confirm the spatial uniformity of the collected charge increases with fluence for polycrystalline Chemical Vapor Deposition diamond, and this effect can also be described by a universal curve.

36 MATERIALS SCIENCE↗
DOE OSTI2018

An Improved Formulation for Hybridizable Discontinuous Galerkin Fluid-Structure Interaction Modeling with Reduced Computational Expense

Here, this work presents two computational efficiency improvements for the hybridizable discontinuous Galerkin (HDG) fluid-structure interaction (FSI) model presented by Sheldon et al. A new formulation for the solid is presented that eliminates the global displacement, resulting in the velocity being the only global solid variable. This necessitates a change to the solid-mesh displacement coupling, which is accounted for by coupling the local solid displacement to the global mesh displacement. Additionally, the mesh basis and test functions are restricted to linear polynomials, rather than being equal-order with the fluid and solid. This change increases the computational efficiency dynamically, with greater benefit the higher order the computation, when compared to an equal-order formulation. These two improvements result in a 50% reduction in the number of global degrees of freedom for high-order simulations for both the fluid and solid domains, as well as an approximately 50% reduction in the number of local fluid domain degrees of freedom for high-order simulations. Finally, the new, more efficient formulation is compared against that from Sheldon et al. and negligible change of accuracy is found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Crystal Structure of the Human Symplekin-Ssu72-CTD Phosphopeptide Complex

Symplekin (Pta1 in yeast) is a scaffold in the large protein complex that is required for 3'-end cleavage and polyadenylation of eukaryotic messenger RNA precursors (pre-mRNAs); it also participates in transcription initiation and termination by RNA polymerase II (Pol II). Symplekin mediates interactions between many different proteins in this machinery, although the molecular basis for its function is not known. Here we report the crystal structure at 2.4 {angstrom} resolution of the amino-terminal domain (residues 30-340) of human symplekin in a ternary complex with the Pol II carboxy-terminal domain (CTD) Ser5 phosphatase Ssu72 and a CTD Ser5 phosphopeptide. The N-terminal domain of symplekin has the ARM or HEAT fold, with seven pairs of antiparallel {alpha}-helices arranged in the shape of an arc. The structure of Ssu72 has some similarity to that of low-molecular-mass phosphotyrosine protein phosphatase, although Ssu72 has a unique active-site landscape as well as extra structural features at the C terminus that are important for interaction with symplekin. Ssu72 is bound to the concave face of symplekin, and engineered mutations in this interface can abolish interactions between the two proteins. The CTD peptide is bound in the active site of Ssu72, with the pSer5-Pro6 peptide bond in the cis configuration, which contrasts with all other known CTD peptide conformations. Although the active site of Ssu72 is about 25 {angstrom} from the interface with symplekin, we found that the symplekin N-terminal domain stimulates Ssu72 CTD phosphatase activity in vitro. Furthermore, the N-terminal domain of symplekin inhibits polyadenylation in vitro, but only when coupled to transcription. Because catalytically active Ssu72 overcomes this inhibition, our results show a role for mammalian Ssu72 in transcription-coupled pre-mRNA 3'-end processing.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Monte Carlo Explicitly Correlated Second-Order Many-Body Green’s Function Calculations of Semiconductor Band Gaps

A systematically converging series of ab initio, post-density-functional, size-consistent, electron-correlated approximations is desired for predictive computing of felectronic band structures of insulating, semiconducting, and metallic solids. A series that meets all of these desiderata (except the applicability to metals) is ab initio many-body Green's function theory based on Gaussian-type-orbital (GTO) basis sets. Here, its leading-order approximation, the second-order Green's function (GF2) method in the diagonal and frequency-independent approximations with the aug-cc-pVDZ basis set, is applied to the fundamental band gaps of three semiconductors (diamond, silicon, and silicon carbide in the zincblende structure) using cluster models. Corrections are made to the basis-set-incompleteness errors by the explicit-correlation (F12) ansatz (GF2-F12) for the valence band edges. The crystals are modeled as surface-passivated clusters of increasing sizes, whose wave functions are expanded by up to 2709 GTO basis functions. Immense computational costs of these calculations are overcome by the highly scalable stochastic algorithm of the Monte Carlo GF2-F12 method, whose operation cost per state increases only as a cubic power of system size, which has a tiny memory footprint and easily achieves near-perfect parallel efficiency on thousands of CPUs or on hundreds of GPUs. The correlated, F12-corrected highest-occupied and lowest-unoccupied molecular-orbital energy (HOMO-LUMO) gap is 5.78 ± 0.07 eV for C 87 H 76 as compared with the experimental value of the fundamental (indirect) band gap of bulk diamond at 5.48 eV. The correlated, F12-corrected HOMO-LUMO gaps for Si 75 H 76 and Si 32 C 43 H 76 are 2.56 ± 0.15 eV and 3.50 ± 0.12 eV, respectively, which are expected to decrease further with increasing cluster sizes. As a result, the experimental fundamental (indirect) band gaps of bulk silicon and silicon carbide are 1.17 eV and 2.42 eV, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Refining perovskite structures to pair distribution function data using collective Glazer modes as a basis

Structural modelling of octahedral tilts in perovskites is typically carried out using the symmetry constraints of the resulting space group. In most cases, this introduces more degrees of freedom than those strictly necessary to describe only the octahedral tilts. It can therefore be a challenge to disentangle the octahedral tilts from other structural distortions such as cation displacements and octahedral distortions. This paper reports the development of constraints for modelling pure octahedral tilts and implementation of the constraints in diffpy-CMI, a powerful package to analyse pair distribution function (PDF) data. The model in the program allows features in the PDF that come from rigid tilts to be separated from non-rigid relaxations, providing an intuitive picture of the tilting. The model has many fewer refinable variables than the unconstrained space group fits and provides robust and stable refinements of the tilt components. It further demonstrates the use of the model on the canonical tilted perovskite CaTiO 3 which has the known Glazer tilt system α + β – β – . The Glazer model fits comparably to the corresponding space-group model Pnma below r = 14 Å and becomes progressively worse than the space-group model at higher r due to non-rigid distortions in the real material.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Accurate Electron Affinities and Orbital Energies of Anions from a Nonempirically Tuned Range-Separated Density Functional Theory Approach

The treatment of atomic anions with Kohn–Sham density functional theory (DFT) has long been controversial because the highest occupied molecular orbital (HOMO) energy, E HOMO , is often calculated to be positive with most approximate density functionals. We assess the accuracy of orbital energies and electron affinities for all three rows of elements in the periodic table (H–Ar) using a variety of theoretical approaches and customized basis sets. Among all of the theoretical methods studied here, we find that a nonempirically tuned range-separated approach (constructed to satisfy DFT-Koopmans’ theorem for the anionic electron system) provides the best accuracy for a variety of basis sets, even for small basis sets where most functionals typically fail. Previous approaches to solve this conundrum of positive E HOMO values have utilized non-self-consistent methods; however, electronic properties, such as electronic couplings/gradients (which require a self-consistent potential and energy), become ill-defined with these approaches. In contrast, the nonempirically tuned range-separated procedure used here yields well-defined electronic couplings/gradients and correct E HOMO values because both the potential and resulting electronic energy are computed self-consistently. Orbital energies and electron affinities are further analyzed in the context of the electronic energy as a function of electronic number (including fractional numbers of electrons) to provide a stringent assessment of self-interaction errors for these complex anion systems.

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

Generalization of the Kohn-Sham system that can represent arbitrary one-electron density matrices

Density functional theory is currently the most widely applied method in electronic structure theory. The Kohn-Sham method, based on a fictitious system of noninteracting particles, is the workhorse of the theory. The particular form of the Kohn-Sham wave function admits only idempotent one-electron density matrices whereas wave functions of correlated electrons in post-Hartree-Fock methods invariably have fractional occupation numbers. Here we show that by generalizing the orbital concept and introducing a suitable dot product as well as a probability density, a noninteracting system can be chosen that can represent the one-electron density matrix of any system, even one with fractional occupation numbers. This fictitious system ensures that the exact electron density is accessible within density functional theory. It can also serve as the basis for reduced density matrix functional theory. Moreover, to aid the analysis of the results the orbitals may be assigned energies from a mean-field Hamiltonian. This produces energy levels that are akin to Hartree-Fock orbital energies such that conventional analyses based on Koopmans' theorem are available. Lastly, this system is convenient in formalisms that depend on creation and annihilation operators as they are trivially applied to single-determinant wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Amorphous transparent conducting oxides in context: Work function survey, trends, and facile modification

The work functions of various amorphous and crystalline transparent conducting oxides (TCO5) were measured using Kelvin probe. The films, made by pulsed laser deposition, exhibited varying work functions dependent on the composition and deposition parameters. Tin oxide showed the largest work functions of the oxides measured, while zinc oxide showed the lowest. Binary and ternary combinations of the basis TCOs showed intermediate work functions dependent on the endpoint components. Amorphous TCO5, important in OPV and other technological applications, exhibited similar work functions to their crystalline counterparts. UV/ozone treatment of TCOs temporarily increased the work function, consistent with proposed defect mechanisms associated with near-surface changes in carrier content and Fermi level. Finally, a method for facile adjustment of the work function of commercial TCOs by atomic layer deposition (ALD) capping layers was presented, illustrated by the growth of zinc oxide layers on commercial crystalline ITO films.

Yeh, T. C.↗
DOE OSTI2021

Elimination of the linearization error in APW/LAPW basis set: Dirac-Kohn-Sham equations

Here, a detailed account of the implementation of equations of the relativistic density functional theory (RDFT) using basis sets of APW/LAPW type with flexible extensions provided by local orbitals is given. Earlier discoveries of the importance of the high derivative local orbital (HDLO) extension of the APW/LAPW basis set for enhancing the accuracy of DFT calculations are confirmed using a fully relativistic approach and α – U as an example. High energy local orbitals (HELO's), however indispensable for GW calculations, are considerably less efficient in enhancing the accuracy of DFT applications. It is shown that a simplified approach to the relativistic effects, namely considering them only inside the muffin-tin (MT) spheres, produces basically identical results (as compared to fully relativistic approach) for the electronic free energy of the five materials considered in this work. By comparing the effect of the simplified approach on the electronic free energy with its effect on the electronic kinetic energy we conclude that the insensitivity of the free energy to the way we describe the relativistic effects in the interstitial region is related to the variational property of this quantity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Density functional study of the electronic structure of dye-functionalized fullerenes and their model donor-acceptor complexes containing P3HT

We study the electronic structure of C{sub 60} fullerenes functionalized with a thiophene-diketo-pyrrolopyrrole-thiophene based chromophore using density functional theory combined with large polarized basis sets. As the attached chromophore has electron donor character, the functionalization of the fullerene leads to a donor-acceptor (DA) system. We examine in detail the effect of the linker and the addition site on the electronic structure of the functionalized fullerenes. We further study the electronic structure of these DA complexes with a focus on the charge transfer excitations. Finally, we examine the interface of the functionalized fullerenes with the widely used poly(3-hexylthiophene-2,5-diyl) (P3HT) donor. Our results show that all functionalized fullerenes with an exception of the C{sub 60}-pyrrolidine [6,6], where the pyrrolidine is attached at a [6,6] site, have larger electron affinities relative to the pristine C{sub 60} fullerene. We also estimate the quasi-particle gap, lowest charge transfer excitation energy, and the exciton binding energies of the functionalized fullerene-P3MT model systems. Results show that the exciton binding energies in these model complexes are slightly smaller compared to a similarly prepared phenyl-C{sub 61}-butyric acid methyl ester (PCBM)-P3MT complex.

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

Accurate potential energy, dipole moment curves, and lifetimes of vibrational states of heteronuclear alkali dimers

We calculate the potential energy curves, the permanent dipole moment curves, and the lifetimes of the ground and excited vibrational states of the heteronuclear alkali dimers XY (X, Y = Li, Na, K, Rb, Cs) in the X{sup 1}Σ{sup +} electronic state using the coupled cluster with singles doubles and triples method. All-electron quadruple-ζ basis sets with additional core functions are used for Li and Na, and small-core relativistic effective core potentials with quadruple-ζ quality basis sets are used for K, Rb, and Cs. The inclusion of the coupled cluster non-perturbative triple excitations is shown to be crucial for obtaining the accurate potential energy curves. A large one-electron basis set with additional core functions is needed for the accurate prediction of permanent dipole moments. The dissociation energies are overestimated by only 14 cm{sup −1} for LiNa and by no more than 114 cm{sup −1} for the other molecules. The discrepancies between the experimental and calculated harmonic vibrational frequencies are less than 1.7 cm{sup −1}, and the discrepancies for the anharmonic correction are less than 0.1 cm{sup −1}. We show that correlation between atomic electronegativity differences and permanent dipole moment of heteronuclear alkali dimers is not perfect. To obtain the vibrational energies and wave functions the vibrational Schrödinger equation is solved with the B-spline basis set method. The transition dipole moments between all vibrational states, the Einstein coefficients, and the lifetimes of the vibrational states are calculated. We analyze the decay rates of the vibrational states in terms of spontaneous emission, and stimulated emission and absorption induced by black body radiation. In all studied heteronuclear alkali dimers the ground vibrational states have much longer lifetimes than any excited states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Systematic Evaluation of Counterpoise Correction in Density Functional Theory

A widespread belief persists that the Boys–Bernardi function counterpoise (CP) procedure “overcorrects” supramolecular interaction energies for the effects of basis-set superposition error. To the extent that this is true for correlated wave function methods, it is usually an artifact of low-quality basis sets. The question has not been considered systematically in the context of density functional theory, however, where basis-set convergence is generally less problematic. We present a systematic assessment of the CP procedure for a representative set of functionals and basis sets, considering both benchmark data sets of small dimers and larger supramolecular complexes. The latter include layered composite polymers with ~150 atoms and ligand–protein models with ~300 atoms. Provided that CP correction is used, we find that intermolecular interaction energies of nearly complete-basis quality can be obtained using only double-ζ basis sets. Furthermore, this is less expensive as compared to triple-ζ basis sets without CP correction. CP-corrected interaction energies are less sensitive to the presence of diffuse basis functions as compared to uncorrected energies, which is important because diffuse functions are expensive and often numerically problematic for large systems. Our results upend the conventional wisdom that CP “overcorrects” for basis-set incompleteness. In small basis sets, CP correction is mandatory in order to demonstrate that the results do not rest on error cancellation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Structural basis for full-spectrum inhibition of translational functions on a tRNA synthetase

The polyketide natural product borrelidin displays antibacterial, antifungal, antimalarial, anticancer, insecticidal and herbicidal activities through the selective inhibition of threonyl-tRNA synthetase (ThrRS). How borrelidin simultaneously attenuates bacterial growth and suppresses a variety of infections in plants and animals is not known. Here we show, using X-ray crystal structures and functional analyses, that a single molecule of borrelidin simultaneously occupies four distinct subsites within the catalytic domain of bacterial and human ThrRSs. These include the three substrate-binding sites for amino acid, ATP and tRNA associated with aminoacylation, and a fourth ‘orthogonal’ subsite created as a consequence of binding. Thus, borrelidin competes with all three aminoacylation substrates, providing a potent and redundant mechanism to inhibit ThrRS during protein synthesis. These results highlight a surprising natural design to achieve the quadrivalent inhibition of translation through a highly conserved family of enzymes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Native state volume fluctuations in proteins as a mechanism for dynamic allostery

Allostery enables tight regulation of protein function in the cellular environment. While existing models of allostery are firmly rooted in the current structure-function paradigm, the mechanistic basis for allostery in the absence of structural change remains unclear. In this study, we show that a typical globular protein is able to undergo significant changes in volume under native conditions while exhibiting no additional changes in protein structure. These native state volume fluctuations were found to correlate with changes in internal motions that were previously recognized as a source of allosteric entropy. This finding offers a novel mechanistic basis for allostery in the absence of canonical structural change. As a result, the unexpected observation that function can be derived from expanded, low density protein states has broad implications for our understanding of allostery and suggests that the general concept of the native state be expanded to allow for more variable physical dimensions with looser packing.

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

Resolution-of-identity stochastic time-dependent configuration interaction for dissipative electron dynamics in strong fields

In this contribution, we introduce a method for simulating dissipative, ultrafast many-electron dynamics in intense laser fields. The method is based on the norm-conserving stochastic unraveling of the dissipative Liouville-von Neumann equation in its Lindblad form. The N-electron wave functions sampling the density matrix are represented in the basis of singly excited configuration state functions. The interaction with an external laser field is treated variationally and the response of the electronic density is included to all orders in this basis. The coupling to an external environment is included via relaxation operators inducing transition between the configuration state functions. Single electron ionization is represented by irreversible transition operators from the ionizing states to an auxiliary continuum state. The method finds its efficiency in the representation of the operators in the interaction picture, where the resolution-of-identity is used to reduce the size of the Hamiltonian eigenstate basis. The zeroth-order eigenstates can be obtained either at the configuration interaction singles level or from a time-dependent density functional theory reference calculation. The latter offers an alternative to explicitly time-dependent density functional theory which has the advantage of remaining strictly valid for strong field excitations while improving the description of the correlation as compared to configuration interaction singles. The method is tested on a well-characterized toy system, the excitation of the low-lying charge transfer state in LiCN.

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

Galaxy bispectrum in the spherical Fourier-Bessel basis

The bispectrum, the three-point correlation in Fourier space, is a crucial statistic for studying many effects targeted by the next-generation galaxy surveys, such as primordial non-Gaussianity (PNG) and general relativistic (GR) effects on large scales. In this work we develop a formalism for the bispectrum in the spherical Fourier-Bessel (SFB) basis—a natural basis for computing correlation functions on the curved sky, as it diagonalizes the Laplacian operator in spherical coordinates. Working in the SFB basis allows for line-of-sight effects such as redshift space distortions and GR to be accounted for exactly, i.e., without having to resort to perturbative expansions to go beyond the plane-parallel approximation. Only analytic results for the SFB bispectrum exist in the literature given the intensive computations needed. We numerically calculate the SFB bispectrum for the first time, enabled by a few techniques: We implement a template decomposition of the redshift-space kernel Z 2 into Legendre polynomials, and separately treat the PNG and velocity-divergence terms. We derive an identity to integrate a product of three spherical harmonics connected by a Dirac delta function as a simple sum and use it to investigate the limit of a homogeneous and isotropic Universe. Furthermore, we present a formalism for convolving the signal with separable window functions and use a toy spherically symmetric window to demonstrate the computation and give insights into the properties of the observed bispectrum signal. While our implementation remains computationally challenging, it is a step toward a feasible full extraction of information on large scales via a SFB bispectrum analysis.

79 ASTRONOMY AND ASTROPHYSICS↗