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 775 records · Page 43

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↗

Draft Function Allocation Framework and Preliminary Technical Basis for Advanced SMR Concepts of Operations

This report presents preliminary research results from the investigation in to the development of new models and guidance for concepts of operations (ConOps) in advanced small modular reactor (aSMR) designs. In support of this objective, three important research areas were included: operating principles of multi-modular plants, functional allocation models and strategies that would affect the development of new, non-traditional concept of operations, and the requirements for human performance, based upon work domain analysis and current regulatory requirements. As part of the approach for this report, we outline potential functions, including the theoretical and operational foundations for the development of a new functional allocation model and the identification of specific regulatory requirements that will influence the development of future concept of operations. The report also highlights changes in research strategy prompted by confirmation of the importance of applying the work domain analysis methodology to a reference aSMR design. It is described how this methodology will enrich the findings from this phase of the project in the subsequent phases and help in identification of metrics and focused studies for the determination of human performance criteria that can be used to support the design process.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Complexity Reduction in Large Quantum Systems: Fragment Identification and Population Analysis via a Local Optimized Minimal Basis

We present, within Kohn-Sham Density Functional Theory calculations, a quantitative method to identify and assess the partitioning of a large quantum mechanical system into fragments. We then introduce a simple and efficient formalism (which can be written as generalization of other well-known population analyses) to extract, from first principles, electrostatic multipoles for these fragments. The corresponding fragment multipoles can in this way be seen as reliable (pseudo-) observables. By applying our formalism within the code BigDFT, we show that the usage of a minimal set of in-situ optimized basis functions is of utmost importance for having at the same time a proper fragment definition and an accurate description of the electronic structure. With this approach it becomes possible to simplify the modeling of environmental fragments by a set of multipoles, without notable loss of precision in the description of the active quantum mechanical region. Furthermore, this leads to a considerable reduction of the degrees of freedom by an effective coarsegraining approach, eventually also paving the way towards efficient QM/QM and QM/MM methods coupling together different levels of accuracy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modulation of hemoglobin dynamics by an allosteric effector

Hemoglobin (Hb) is an extensively studied paradigm of proteins that alter their function in response to allosteric effectors. Models of its action have been used as prototypes for structure-function relationships in many proteins, and models for the molecular basis of its function have been deeply studied and extensively argued. Recent reports suggest that dynamics may play an important role in its function. Relatively little is known about the slow, correlated motions of hemoglobin subunits in various structural states because experimental and computational strategies for their characterization are challenging. Allosteric effectors such as inositol hexaphosphate (IHP) bind to both deoxy-Hb and HbCO, albeit at different sites, leading to a lowered oxygen affinity. The manner in which these effectors impact oxygen binding is unclear and may involve changes in structure, dynamics or both. Here we use neutron spin echo (NSE) measurements accompanied by wideangle x-ray scattering (WAXS) to show that binding of IHP to HbCO results in an increase in the rate of coordinated motions of Hb subunits relative to one another with little if any change in large scale structure. This increase of large-scale dynamics seems to be coupled with a decrease in the average magnitude of higher frequency modes of individual residues. Furthermore, these observations indicate that enhanced dynamic motions contribute to the functional changes induced by IHP and suggest that they may be responsible for the lowered oxygen affinity triggered by these effectors.

59 BASIC BIOLOGICAL SCIENCES↗

Determination of a Density Functional Tight Binding Model with an Extended Basis Set and Three-Body Repulsion for Carbon Under Extreme Pressures and Temperatures

In this paper, we report here on development of a density functional tight binding (DFTB) simulation approach for carbon under extreme pressures and temperatures that includes an expanded basis set and an environmentally dependent repulsive energy. We find that including d-orbital interactions in the DFTB Hamiltonian improves determination of the electronic states at high pressure–temperature conditions, compared to standard DFTB implementations that utilize s- and p-orbitals only for carbon. We then determine a three-body repulsive energy through fitting to diamond, BC8, and simple cubic cold compression curve data, as well pressures from metallic liquid configurations from density functional theory (DFT) simulations. Our new model (DFTB-p3b) yields approximately 2 orders of magnitude increase in computational efficiency over standard DFT while retaining its accuracy for condensed phases of carbon under a wide range of conditions, including the metallic liquid phase at conditions up to 2000 GPa and 30 000 K. Finally, our results provide a straightforward method by which DFTB can be extended to studies of covalently bonded materials under extremely high pressures and temperatures such as the interiors of planets and other large celestial bodies.

36 MATERIALS SCIENCE↗

Exchange correlation potentials from full configuration interaction in a Slater orbital basis

Ryabinkin–Kohut–Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, full configuration interaction computations in the Slater orbital basis are employed to give quality input to RKS, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. SlaterRKS is shown to be an efficient algorithm to arrive at exchange-correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. As a result, examples of weakly and strongly correlated molecular systems will demonstrate the main features of SlaterRKS.

74 ATOMIC AND MOLECULAR PHYSICS↗

Influence of the complete basis set approximation, tight weighted-core, and diffuse functions on the DLPNO-CCSD(T1) atomization energies of neutral H,C,O-compounds

The impact of complete basis set extrapolation schemes (CBS), diffuse functions, and tight weighted core functions on enthalpies of formation predicted via the DLPNO-CCSD(T1) reduced Feller Peterson-Dixon approach has been examined for neutral H,C,O-compounds. All tested three-point (TZ/QZ/5Z) extrapolation schemes result in mean unsigned deviation (MUD) below 2 kJ mol -1 relative to the experiment. Here, the two-point QZ/5Z and TZ/QZ CBS 1/$l^3_{max}$ extrapolation schemes are inferior to their inverse power counterpart (1/($l_{max}$ + 1/2) 4 ) by 1.3 and 4.3 kJ mol -1 . The CBS extrapolated frozen core atomization energies are insensitive (within 1 kJ mol -1 ) to augmentation of the basis set with tight weighted core functions. The core-valence correlation effects converge already at triple-ζ, although double- ζ/triple- ζ CBS extrapolation performs better and is recommended. The effect of diffuse function augmentation converges slowly, and cannot be reproduced with double- ζ or triple- ζ calculations as these are plagued with basis set superposition and incompleteness errors.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multicopy Single-Stranded DNA Directs Intestinal Colonization of Enteric Pathogens

Multicopy single-stranded DNAs (msDNAs) are hybrid RNA-DNA molecules encoded on retroelements called retrons and produced by the action of retron reverse transcriptases. Retrons are widespread in bacteria but the natural function of msDNA has remained elusive despite 30 years of study. The major roadblock to elucidation of the function of these unique molecules has been the lack of any identifiable phenotypes for mutants unable to make msDNA. We report that msDNA of the zoonotic pathogen Salmonella Typhimurium is necessary for colonization of the intestine. Similarly, we observed a defect in intestinal persistence in an enteropathogenic E. coli mutant lacking its retron reverse transcriptase. Under anaerobic conditions in the absence of msDNA, proteins of central anaerobic metabolism needed for Salmonella colonization of the intestine are dysregulated. We show that the msDNA-deficient mutant can utilize nitrate but not other alternate electron acceptors in anaerobic conditions. Consistent with the availability of nitrate in the inflamed gut, a neutrophilic inflammatory response partially rescued the ability of a mutant lacking msDNA to colonize the intestine. These findings together indicate that the mechanistic basis of msDNA function during Salmonella colonization of the intestine is proper production of proteins needed for anaerobic metabolism. We further conclude that a natural function of msDNA is to regulate protein abundance, the first attributable function for any msDNA. Our data provide novel insight into the function of this mysterious molecule that likely represents a new class of regulatory molecules.

Environmental Molecular Sciences Laboratory↗

State-of-the-art ab initio potential energy curve for the krypton atom pair and thermophysical properties of dilute krypton gas

A new reference krypton-krypton interatomic potential energy curve was developed by means of quantum-chemical ab initio calculations for 36 interatomic separations. Highly accurate values for the interaction energies at the complete basis set limit were obtained using the coupled-cluster method with single, double, and perturbative triple excitations as well as t-aug-cc-pV5Z and t-aug-cc-pV6Z basis sets including mid-bond functions, with the 6Z basis set being newly constructed for this study. Higher orders of coupled-cluster terms were considered in a successive scheme up to full quadruple excitations. Core-core and core-valence correlation effects were included. Furthermore, relativistic effects were studied not only at a scalar relativistic level using second-order direct perturbation theory, but also utilizing full four-component and Gaunt-effect computations. An analytical pair potential function was fitted to the interaction energies, which is characterized by a depth of 200.88 K with an estimated standard uncertainty of 0.51 K. Thermophysical properties of low-density krypton were calculated for temperatures up to 5000 K. Second and third virial coefficients were obtained from statistical thermodynamics. Viscosity and thermal conductivity as well as the self-diffusion coefficient were computed using the kinetic theory of gases. The theoretical results are compared with experimental data and with results for other pair potential functions from the literature, especially with those calculated from the recently developed ab initio potential of Waldrop et al. [J. Chem. Phys. 142, 204307 (2015)]. Highly accurate experimental viscosity data indicate that both the present ab initio pair potential and the one of Waldrop et al. can be regarded as reference potentials, even though the quantum-chemical methods and basis sets differ. However, the uncertainties of the present potential and of the derived properties are estimated to be considerably lower.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Quaternionic Harmonic Analysis of Texture

QHAT uses various functions and data structures native to MATLAB to analyze crystallographic texture information using harmonic functions on the space of rotations represented as normalized quaternions. These harmonic functions generalize the spherical harmonics in three dimensions, and form the basis for the irreducible representations of the four-dimensional rotation group. This allows the basis of harmonic functions to be reduced to linearly independent combinations that satisfy the crystal and sample symmetry point groups.

Mason, J.↗

The Electrolyte Genome project: A big data approach in battery materials discovery

We present a high-throughput infrastructure for the automated calculation of molecular properties with a focus on battery electrolytes. The infrastructure is largely open-source and handles both practical aspects (input file generation, output file parsing, and information management) as well as more complex problems (structure matching, salt complex generation, and failure recovery). Using this infrastructure, we have computed the ionization potential (IP) and electron affinities (EA) of 4830 molecules relevant to battery electrolytes (encompassing almost 55,000 quantum mechanics calculations) at the B3LYP/6-31+G(*) level. We describe automated workflows for computing redox potential, dissociation constant, and salt-molecule binding complex structure generation. We present routines for automatic recovery from calculation errors, which brings the failure rate from 9.2% to 0.8% for the QChem DFT code. Automated algorithms to check duplication between two arbitrary molecules and structures are described. We present benchmark data on basis sets and functionals on the G2-97 test set; one finding is that a IP/EA calculation method that combines PBE geometry optimization and B3LYP energy evaluation requires less computational cost and yields nearly identical results as compared to a full B3LYP calculation, and could be suitable for the calculation of large molecules. Our data indicates that among the 8 functionals tested, XYGJ-OS and B3LYP are the two best functionals to predict IP/EA with an RMSE of 0.12 and 0.27 eV, respectively. Application of our automated workflow to a large set of quinoxaline derivative molecules shows that functional group effect and substitution position effect can be separated for IP/EA of quinoxaline derivatives, and the most sensitive position is different for IP and EA. Published by Elsevier B.V

Battery↗

Structural and Functional Analysis of JMJD2D Reveals Molecular Basis for Site-Specific Demethylation among JMJD2 Demethylases

We found that JMJD2 lysine demethylases (KDMs) participate in diverse genomic processes. Most JMJD2 homologs display dual selectivity toward H3K9me3 and H3K36me3, with the exception of JMJD2D, which is specific for H3K9me3. Here, we report the crystal structures of the JMJD2D•2-oxoglutarate•H3K9me3 ternary complex and JMJD2D apoenzyme. Utilizing structural alignments with JMJD2A, molecular docking, and kinetic analysis with an array of histone peptide substrates, we elucidate the specific signatures that permit efficient recognition of H3K9me3 by JMJD2A and JMJD2D, and the residues in JMJD2D that occlude H3K36me3 demethylation. Surprisingly, these results reveal that JMJD2A and JMJD2D exhibit subtle yet important differences in H3K9me3 recognition, despite the overall similarity in the substrate-binding conformation. Further, we show that H3T11 phosphorylation abrogates demethylation by JMJD2 KDMs. These studies reveal the molecular basis for JMJD2 site specificity and provide a framework for structure-based design of selective inhibitors of JMJD2 KDMs implicated in disease.

Krishnan, Swathi↗