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.

SourceResearch connectionFrom yearThrough year

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

At least 1,243 records · Page 69

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

A Parallelized Oxidation-Driven Surface Recession Framework in DSMC Code, SPARTA

Spacecrafts rely on ablative thermal protection systems (TPS) made of composites consisting of a carbon-based reinforcement and a polymeric matrix. These materials are designed to withstand high-temperature oxidation and surface recession during re-entry into the Earth's atmosphere. However, ablation occurs due to a complex interplay of thermal, mechanical, and chemical factors, making it challenging to determine the individual impact of each on the TPS's overall degradation. In this study, we have developed an ablation model that can leverage a finite rate carbon oxidation model to predict material recession and surface states more accurately. Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al can perform both gas-surface and pure-surface reactions and is now adapted to tally surface data on a per grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also briefly discuss some unique challenges associated with parallelizing this dynamic surface state and geometry. Finally, we analyze the performance of this parallelized implicit chemistry model with simple 2D and 3D benchmark cases by producing surface state statistics, area changes over time, and visualization across a range of surface temperatures and processors with and without load-balancing.

DSMC↗

NASA NTRSDate not supplied

NASA CARA Tools

The NASA Conjunction Assessment Risk Analysis (CARA) team has created and posted to GitHub a set of MATLAB software tools for collision assessment (CA) analysis. The software is provided to the public under a NASA Open-Source Software Agreement and is maintained/updated by the CARA Analysis team. Capabilities provided within the toolset include algorithms in assessing probability of collision (Pc), collision consequence, covariance realism, Orbit Determination (OD) quality assessment, and single covariance max Pc. Additionally, a set of utilities are provided for astrodynamics topics closely related to CA, such as: Conjunction Data Message (CDM) parsing, covariance transformations, coordinate system transformations, etc. The repository is regularly updated on a monthly basis and new functionality will be added as algorithms are approved for public release. This presentation will provide a brief introduction to the toolset, highlight some of the key components which can be used in operations, and will present information on algorithms that are in the release pipeline and can be expected soon. This software will enable operators to independently perform CA calculations using the same core algorithms that CARA uses in operations.

L. Baars↗

NASA NTRSDate not supplied

NASA CARA Tools

Luis Baars↗

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↗

NASA NTRS1977

Symbolic Laplace transforms of special functions

A MACSYMA implementation of the Laplace transform for special functions is described. The generalized hypergeometric functions are used as a basis for the representation of approximately fifty special functions. Only a relatively small number of formulas that generally involve generalized hypergeometric functions are utilized for the integration stage. A sample of actual examples and their timing is provided.

Avgoustis, Y.↗

NASA NTRS1998

The Heats of Formation of GaCl3 and its Fragments

The heats of formation of GaC13 and its fragments are computed. The geometries and frequencies are obtained at the B3LYP level. The CCSD(T) approach is used to solve the correlation problem. The effect of Ga 3d correlation is studied, and found to affect the bond energies by up to 1 kcal/mol. Both basis set extrapolation and bond functions are considered as ways to approach the basis set limit. Spin-orbit and scalar relativistic effects are also considered.

Bauschlicher, Charles W., Jr.↗

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↗

NASA NTRS1998

The Dissociation Energies of He2, HeH, and ArH; A Bond Function Study

The bond energies and bond lengths are determined for He2, HeH, and ArH at the CCSD(T) level using both atom-centered basis sets and those that include bond functions. The addition of bond functions dramatically improves the rate of convergence of the results with respect to the size of the atom-centered basis set; with bond functions, triple zeta atom-centered basis set, outperform quintuple zeta basis sets without bond functions. The addition of bond functions also reduces the number of diffuse functions that must be added to the atom-centered sets. Employing bond functions appear to offer a very cost effective method of computing the interaction between weakly bound systems, especially for He.

Bauschlicher, Charles W., Jr.↗

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↗

NASA NTRS1986

Transfer function concept for ultrasonic characterization of material microstructures

The approach given depends on treating material microstructures as elastomechanical filters that have analytically definable transfer functions. These transfer functions can be defined in terms of the frequency dependence of the ultrasonic attenuation coefficient. The transfer function concept provides a basis for synthesizing expressions that characterize polycrystalline materials relative to microstructural factors such as mean grain size, grain-size distribution functions, and grain boundary energy transmission. Although the approach is nonrigorous, it leads to a rational basis for combining the previously mentioned diverse and fragmented equations for ultrasonic attenuation coefficients.

Vary, A.↗

NASA NTRS1991

Improved techniques for outgoing wave variational principle calculations of converged state-to-state transition probabilities for chemical reactions

Improved techniques and well-optimized basis sets are presented for application of the outgoing wave variational principle to calculate converged quantum mechanical reaction probabilities. They are illustrated with calculations for the reactions D + H2 yields HD + H with total angular momentum J = 3 and F + H2 yields HF + H with J = 0 and 3. The optimization involves the choice of distortion potential, the grid for calculating half-integrated Green's functions, the placement, width, and number of primitive distributed Gaussians, and the computationally most efficient partition between dynamically adapted and primitive basis functions. Benchmark calculations with 224-1064 channels are presented.

Mielke, Steven L.↗