Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Effective Potential Energies and Transport Cross Sections for Interactions of Hydrogen and Nitrogen

The interaction energies for N2-He and N2-H2 are calculated by accurate ab initio methods. The virial coefficient and differential scattering cross section for N2-H2 are calculated; the theoretical results are compared with experimental data. The transport collision integrals for N2-H2 and N2-N2 interactions are calculated and tabulated; the results yield transport coefficients that compare well with measured data. Transport coefficients are found to be determined accurately from the interaction energies for a specific configuration of the molecule formed from the interaction partners. Comparisons with results of measurement and accurate calculations demonstrate that the transport properties of complex molecular interactions can be determined rapidly and fairly accurately from the interaction energies of simpler system using combination rules for the short-range parameters of effective interaction energies and the coefficients for the long-range forces. The coefficients for a two-parameter temperature expansion of diffusion and viscosity are tabulated for a realistic universal potential energy that is based primarily on the results of very accurate calculations of the He-He interaction energy.

Stallcop, James R.↗
NASA NTRS2001

Effective Potential Energies and Transport Cross Sections for Atom-Molecule Interactions of Nitrogen and Nitrogen

The potential energy surfaces for H2-N and N2-N interactions are calculated by accurate ab initio methods and applied to determine transport data. The results confirm that an effective potential energy for accurately determining transport properties can be calculated using a single orientation. A simple method is developed to determine the dispersion coefficients of effective potential energies Effective potential energies required for O2-O collisions are determ=ined. The H2-N, N2-N, O2-H, and O2-O collision integrals are calculated and tabulated for a large range of temperatures. The theoretical values of the N2-N and O2-O diffusion coefficients compare well with measured data available at room temperature.

Stallcop, James R.↗
NASA NTRS1995

An Extended Ab Initio and Theoretical Thermodynamics Studies of the Bergman Reaction and the Energy Splitting of the Singlet Ortho-, Meta-, and Para-Benzynes

The autoaromatization of (Z)-hex-3-ene-1,5-diyne to the singlet biradical para-benzyne has been reinvestigated by state of the art ab initio methods. Previous CCSD(T)/6-31G(d,p) and CASPT2[0]/ANO[C(5s4p2d1f)/H(3s2p)] calculations estimated the the reaction heat at 298 K to be 8-10 and 4.9 plus or minus 3.2 kcal/mol, respectively. Recent NO- and oxygen-dependent trapping experiments and collision-induced dissociation threshold energy experiments estimate the heat of reaction to be 8.5 plus or minus 1.0 at 470 K (recomputed to 9.5 plus or minus 1.0 at 298 K) and 8.4 plus or minus 3.0 kcal/mol at 298 K, respectively. New theoretical estimates at 298 K predict the values at the basis set limit for the CCSD(T) and CASPT2(g1) methods to be 12.7 plus or minus 2.0 and 5.4 plus or minus 2.0 kcal/mol, respectively. The experimentally predicted electronic contribution to the heat of activation is 28.6 kcal/mol. This can be compared with 25.5 and 29.8 kcal/mol from the CASPT2[g1] and the CCSD(T) methods, respectively. The new study has in particular improved on the one-particle basis set for the CCSD(T) method as compared to earlier studies. For the CASPT2 investigation the better suited CASPT2[g1] approximation is utilized. The original CASPT2 method, CASPT2[0], systematically favors open shell systems relative to closed shell systems. This was previously corrected empirically. The study shows that the energy difference between CCSD(T) and CASPT2[g1] at the basis set limit is estimated to be 7 plus or minus 2 kcal/mol. The study also demonstrates that the estimated heat of reaction is very sensitive to the quality of the basis set.

Lindh, Roland↗
NASA NTRS1996

Quantum Mechanical Determination of Potential Energy Surfaces for TiO and H2O

We discuss current ab initio methods for determining potential energy surfaces, in relation to the TiO and H2O molecules, both of which make important contributions to the opacity of oxygen-rich stars. For the TiO molecule we discuss the determination of the radiative lifetimes of the excited states and band oscillator strengths for both the triplet and singlet band systems. While the theoretical radiative lifetimes for TiO agree well with recent measurements, the band oscillator strengths differ significantly from those currently employed in opacity calculations. For the H2O molecule we discuss the current results for the potential energy and dipole moment ground state surfaces generated at NASA Ames. We show that it is necessary to account for such effects as core-valence Correlation energy to generate a PES of near spectroscopic accuracy. We also describe how we solve the ro-vibrational problem to obtain the line positions and intensities that are needed for opacity sampling.

Langhoff, Stephen R.↗
NASA NTRS1996

Towards the Synthesis and Identification of Tetrahedral N4

Tetrahedral (Td)-N4 has attracted attention as a high energy density material candidate. Although computational studies have shown Td-N4 to be a stable structure on the ground-state potential energy surface, it has yet to be experimentally observed. This study investigates possible synthetic pathways via excited state surfaces. Preliminary results include a comparison of vertical excitation energies from Td-N4 using CIS, CIS(D), CIS(3), LRCCSD, CASSCF, CASPT2, and MRCI ab initio methods. Previous studies have predicted weak IR vibrational intensities. Therefore, alternative detection schemes of Td-N4 are also investigated, such as UV and vibrational Raman spectroscopy.

Dateo, Christopher E.↗
NASA NTRS2002

Gas-Phase Reaction Pathways and Rate Coefficients for the Dichlorosilane-Hydrogen and Trichlorosilane-Hydrogen Systems

As part of NASA Ames Research Center's Integrated Process Team on Device/Process Modeling and Nanotechnology our goal is to create/contribute to a gas-phase chemical database for use in modeling microelectronics devices. In particular, we use ab initio methods to determine chemical reaction pathways and to evaluate reaction rate coefficients. Our initial studies concern reactions involved in the dichlorosilane-hydrogen (SiCl2H2--H2) and trichlorosilane-hydrogen (SiCl2H-H2) systems. Reactant, saddle point (transition state), and product geometries and their vibrational harmonic frequencies are determined using the complete-active-space self-consistent-field (CASSCF) electronic structure method with the correlation consistent polarized valence double-zeta basis set (cc-pVDZ). Reaction pathways are constructed by following the imaginary frequency mode of the saddle point to both the reactant and product. Accurate energetics are determined using the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations (CCSD(T)) extrapolated to the complete basis set limit. Using the data from the electronic structure calculations, reaction rate coefficients are obtained using conventional and variational transition state and RRKM theories.

Dateo, Christopher E.↗
NASA NTRS2017

The Effect of the Spin-Forbidden Co((sup 1) Sigma plus) plus O((sup 3) P) Yields CO2 (1 Sigma (sub G) plus) Recombination Reaction on Afterbody Heating of Mars Entry Vehicles

Vibrationally excited CO2, formed by two-body recombination from CO((sup 1) sigma plus) and O((sup 3) P) in the wake behind spacecraft entering the Martian atmosphere reaction, is potentially responsible for the higher than anticipated radiative heating of the backshell, compared to pre-flight predictions. This process involves a spin-forbidden transition of the transient triplet CO2 molecule to the longer-lived singlet. To accurately predict the singlet-triplet transition probability and estimate the thermal rate coefficient of the recombination reaction, ab initio methods were used to compute the first singlet and three lowest triplet CO2 potential energy surfaces and the spin-orbit coupling matrix elements between these states. Analytical fits to these four potential energy surfaces were generated for surface hopping trajectory calculations, using Tully's fewest switches surface hopping algorithm. Preliminary results for the trajectory calculations are presented. The calculated probability of a CO((sup 1) sigma plus) and O((sup 3) P) collision leading to singlet CO2 formation is on the order of 10 (sup -4). The predicted flowfield conditions for various Mars entry scenarios predict temperatures in the range of 1000 degrees Kelvin - 4000 degrees Kelvin and pressures in the range of 300-2500 pascals at the shoulder and in the wake, which is consistent with a heavy-particle collision frequency of 10 (sup 6) to 10 (sup 7) per second. Owing to this low collision frequency, it is likely that CO((sup 1) sigma plus) molecules formed by this mechanism will mostly be frozen in a highly nonequilibrium rovibrational energy state until they relax by photoemission.

CO2 recombination↗
NASA NTRSDate not supplied

Hierarchical screening for Li-based solid electrolytes using fast, interpretable machine-learned potentials

Li-based solid-state electrolyte materials enable safer, all-solid-state batteries but the computational search for candidates with favorable stability and Li-ion conductivity is challenging due to the size of the search space and the cost of evaluating transport properties with ab initio methods. The prohibitive cost of high-throughput screening with DFT has lead to the development of surrogate models using geometric analysis, empirical potentials, and descriptors for ionic transport. Here, I will discuss a hierarchical screening approach for identifying promising materials using a combination of density functional theory, bond-valence methods, and machine learning potentials generated with the Ultra-Fast Force Fields (UF3) framework. We show how the inexpensive bond-valence method can be used to guide the generation of training samples for machine learning, in addition to filtering candidates. Finally, we apply the hierarchical workflow to screen for ionic conductivity across a database of Li-containing compounds.

Materials discovery↗
NASA NTRSDate not supplied

Hierarchical Screening for Li-Based Solid Electrolytes Using Fast, Interpretable Machine-Learned Potentials

Li-based solid-state electrolyte materials enable safer, all-solid-state batteries but the computational search for candidates with favorable stability and Li-ion conductivity is challenging due to the size of the search space and the cost of evaluating transport properties with ab initio methods. The prohibitive cost of high-throughput screening with DFT has lead to the development of surrogate models using geometric analysis, empirical potentials, and descriptors for ionic transport. Here, I will discuss a hierarchical screening approach for identifying promising materials using a combination of density functional theory, bond-valence methods, and machine learning potentials generated with the Ultra-Fast Force Fields (UF3) framework. We show how the inexpensive bond-valence method can be used to guide the generation of training samples for machine learning, in addition to filtering candidates.

Materials discovery↗
NASA NTRSDate not supplied

First-Principles Investigation of Phase Transitions in RuNb, RuTa, and Ru-Nb-Ta Ultrahigh Temperature Shape Memory Alloys

Ultrahigh temperature shape memory alloys (UHT-SMAs) have transition temperatures above 600 C, and they have found applications for sensing and actuating devices in aerospace industry. Among very few such UHT-SMAs currently known are Ru-based alloys such as RuNb and RuTa, whose martensite structures and phase transitions are totally different from those of NiTi-based SMAs and were poorly understood. In this work, we carried out a systematical study of RuNb using first-principles total energy calculations and molecular dynamics (MD) simulations. The transition paths and mechanisms in cubic → tetragonal →monoclinic transitions are revealed. and the transition sequence and martensitic transition temperatures are determined (MTTs) by evaluating the Gibbs free energies using thermodynamic integration. The calculated MTTs are in very good agreement with the experimental data. We found that the monoclinic phase at the second transition has the 𝑃21/m symmetry instead of experimentally identified 𝑃2/m. Our calculations demonstrate that RuTa has very similar phase transitions to those of RuTa. Furthermore, we studied the Ru0.5Nb0.25Ta0.25 ternary. Our results can explain the measured significant drop of MTT (~ 200 C) in the second transition for Ru0.5Nb0.25Ta0.25 compared with those of RuNb and RuTa, while in the first transition its MTT is between MTTs of RuNb and RuTa. The insights gained by this study and the verified ab initio methods for accurate MTT calculations can be applied to fast screen and quantitatively design novel UHT-SMAs having similar properties with desirable MTTs and much reduced cost.

Zhigang Wu↗
DOE OSTI2016

DEVELOPMENT OF INTERATOMIC POTENTIALS IN TUNGSTEN-RHENIUM SYSTEMS

Reference data are generated using the ab initio method to fit interatomic potentials for the W-Re system. The reference data include single phases of W and Re, strained structures, slabs, systems containing several concentrations of vacancies, systems containing various types of interstitial defects, melt structures, structures in the σ and χ phases, and structures containing several concentrations of solid solutions of Re in bcc W and W in hcp Re. Future work will start the fitting iterations.

Density Functional Theory↗
DOE OSTI2026

Computing neutrino cross sections from Euclidian responses

Energy integrated neutrino cross sections are integrals of nuclear responses weighted with kinematic prefactors. We decompose the prefactors into a limited set of functions of energy transfer and show the relevant integrals are the moments of the responses, and integrals weighted with $1/(a+ω)^n$ with $n\leq 2$. These can be directly obtained from the Euclidean response, avoiding the need for inversion of the Laplace transform. As a proof of concept we study the procedure with toy-model responses for the quasielastic peak. We show that the different contributions can be straightforwardly organized in terms of relative importance, and how flux-averaged cross sections can be obtained. Using a realistic model for the response and numerical uncertainty we show that it is feasible to obtain the required integrals from the Euclidean response, with large uncertainties only for the third moment. Due to kinematic restrictions, the integrals contain contributions from the unphysical region for neutrino scattering, coming from high-momentum nucleons. We show that (in the absence of two-body currents) robust corrections for this contamination are obtained from the single-nucleon momentum distribution. These results present an opportunity to compute certain neutrino cross sections with ab-initio methods with controlled uncertainties.

Nikolakopoulos, A. [Washington U., Seattle]↗
DOE OSTI2026

Computing neutrino cross sections from Euclidian responses

Energy integrated neutrino cross sections are integrals of nuclear responses weighted with kinematic prefactors. We decompose the prefactors into a limited set of functions of energy transfer and show the relevant integrals are the moments of the responses, and integrals weighted with $1/(a+ω)^n$ with $n\leq 2$. These can be directly obtained from the Euclidean response, avoiding the need for inversion of the Laplace transform. As a proof of concept we study the procedure with toy-model responses for the quasielastic peak. We show that the different contributions can be straightforwardly organized in terms of relative importance, and how flux-averaged cross sections can be obtained. Using a realistic model for the response and numerical uncertainty we show that it is feasible to obtain the required integrals from the Euclidean response, with large uncertainties only for the third moment. Due to kinematic restrictions, the integrals contain contributions from the unphysical region for neutrino scattering, coming from high-momentum nucleons. We show that (in the absence of two-body currents) robust corrections for this contamination are obtained from the single-nucleon momentum distribution. These results present an opportunity to compute certain neutrino cross sections with ab-initio methods with controlled uncertainties.

Nikolakopoulos, A. [Washington U., Seattle]↗
DOE OSTI2019

The comparison of semiempirical and ab initio molecular modeling methods in activity and property evaluation of selected antimicrobial sulfonamides

A group of sulfonamides exhibiting antimicrobial activity was analyzed with the application of quantitative structure–activity/property relationships method. The purpose of this study was to show the common and differentiating characteristics of the analyzed chemical structures alike physicochemically as well as pharmacologically based on the quantum-chemical calculations both in vacuo and in the aquatic environment together with their microbiological activity and chromatographic retention data. The semiempirical and ab initio level of in silico molecular modeling was performed for calculations of statistically significant molecular descriptors to compare obtained results. The relationship between the structure and biological activity and physicochemical parameters data was able to class and describe analyzed molecules and the applied chemometric approaches (principal component analysis, factor analysis, and multiple regression analysis) revealed the influential features of the tested structures responsible for the antimicrobial activity.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Toward Generalized Tensor Algebra for ab initio Quantum Chemistry Methods

The widespread use of tensor operations in describing electronic structure calculations has motivated the design of software frameworks for productive development of scalable optimized tensor-based electronic structure methods. Whereas prior work focused on Cartesian abstractions for dense tensors, we present an algebra to specify and perform tensor operations on a larger class of block-sparse tensors. We illustrate the use of this framework in expressing real-world computational chemistry calculations beyond the reach of existing frameworks.

Mutlu, Erdal↗
DOE OSTI2019

Derivation of Lanthanide Series Crystal Field Parameters From First Principles.

Two series of lanthanide complexes were chosen to analyze the trends of the magnetic properties along the series: the highly symmetric LnZn 16 (picHA) 16 series with Ln = Tb, Dy, Ho, Er, Yb and picHA=picoline hydroxamic, and the [Ln(DPA) 3 ](C 3 H 5 N 2 ) 3 · 3H 2 O series with DPA = 2,6-dipicolinic acid and Ln = Ce–Yb and the approximate threefold symmetry. The first series present a compressed coordination sphere of 8 oxygen atoms environment while in the 2nd series, the coordination sphere is formed by an elongated coordination sphere formed by 6 oxygen atoms. Crystal field parameters (CFPs) are deduced from ab initio calculations using two methods: the AILFT (Ab Initio Ligand Field Theory) which determines the parameters at the orbital level, and the ITO (Irreducible Technique Operator) decomposition which treats the problem at the many-electron level. It is shown that the CFPs are transferable from one derivative to the other, within a given series. The sign of the 2nd order parameter $B_0^2$ differs in the two series reflecting the difference of environment. It is shown that the use of strength parameter $S$ permits a easy comparison between complexes. Furthermore, in both cases, the parameters are found to decrease in magnitude along the series and this decrease is imputed to covalent effects.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗