Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio calculations”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Energetics of the protonation of CO - Implications for the observation of HOC(+) in dense interstellar clouds

A number of molecular species on the H3CO(+) energy hypersurface is examined. Ab initio molecular orbital theory is used to determine the structures and relative energetics of the two isomers of HCO(+) and HOC(+) together with the affinity of CO for protonation at either end. The proton affinities of H2 and H2CO are also examined. The calculations are performed using large basis sets and include the effects of electron correlation. The calculated vibrational frequencies are used to correct for zero point energy differences. The results show that the proton affinities of H2 and CO to form HOC(+) are within 1 kcal of each other. The calculations demonstrate that there is no thermodynamic driving force to form HOC(+) in collisions of H3(+) with CO, and that the formation of HCO(+) in such collisions is very exoergic. A plausible mechanism is suggested to explain the differences observed between the laboratory and the interstellar medium.

Dixon, D.↗
NASA NTRS2015

Empirical InfraRed Line Lists for Five SO2 Isotopologues: 32/33/34/36S16O2 and 32S18O2

Using the latest published, empirically refined potential energy surface (PES) Ames-1 and purely ab initio CCSD(T)/aug-cc-pVQZ dipole moment surface (DMS), we have computed Infrared line lists for five symmetric isotopologues of sulfur dioxide: 32S16O2 (626), 33S16O2 (636), 34S16O2 (646), 36S16O2 (666), and 32S18O2 (828). The line lists are based on J=0-80 rovibrational variational calculations with E’ ≤ 8000 cm-1. The 34S16O2 and 33S16O2 line lists are compared to the experiment-based models in the HIgh-resolution TRANsmission molecular absorption database (HITRAN2012, http://www.cfa.harvard.edu/hitran/) and the Cologne Database for Molecular Spectroscopy, CDMS (http://www.astro.uni-koeln.de/cdms/). The accuracy for computed 646 band origins is similar to what has been reported for the main isotopologue, i.e. 0.01 – 0.03 cm-1 for bands up to 5500 cm-1. For rovibrational transitions, the 646 line position and intensity deviation patterns are much simpler and more self-consistent than those of the main isotopologue 626. The discrepancies are mainly found for higher Ka/J transitions. 626 and 646 exhibit comparable line position and intensity agreement for lower Ka/J transitions. The line position deviations for the 636 purely rotational band are parallel to those of 626 and 646, while its line intensity deviations do not show branching patterns as we found in the 626 and 646 cases. Predictions for the other minor isotopologues are expected to exhibit similar accuracy. These line lists are accurate enough to provide alternatives for missing bands of 626 and the minor isotopologues. It may significantly facilitate the laboratory spectroscopic measurement and analysis, as well as to identify these isotopologues in various astrophysical environments.

Sulfur Dioxide↗
NASA NTRS1987

Experimental and theoretical investigation of the far-infrared spectrum of H2-He mixtures

New measurements of the rototranslational collision-induced absorption spectra in the far-infrared (FIR) frequency region of gaseous mixtures of hydrogen and helium are reported. The frequency range extends up to 1500/cm at the temperature of 195 K, and 1700/cm at 296 K. The uncertainties of the measurement are assessed. Comparison with previous measurements shows substantial agreement at frequencies where the measurements overlap. For comparison with the measurement, spectra are computed from the fundamental theory using a state of the art ab initio dipole function and an isotropic potential surface as input. The observed agreement suggests that, for important applications, such as the analysis of the FIR spectra of the outer planets, the spectra can be obtained reliably as a function of frequency and temperature from quantum calculations of the kind presented here.

Birnbaum, George↗
NASA NTRS1992

Theoretical study of hydrogen and nitrogen interactions - N-H transport cross sections and collision integrals

Potential energy curves have been calculated for the X 3Sigma(-) and 5Sigma(-) states of NH using a complete-active-space self-consistent-field/multireference configuration-interaction procedure and extensive Gaussian basis sets. An accurate dipole moment function is computed for the X 3Sigma(-) state. Transport cross sections have been calculated for the collisions of hydrogen and nitrogen atoms in their ground states. The mean transport cross sections are tabulated for collision energies in the range 0.0001-1 E(h); the mean collision integrals are listed for temperatures between 100 and 100,000 K. The ab initio energies for the 5 Sigma(-) state of NH are found to be consistent with the results of correlation studies and agree well with the prediction from combining relations using the energies of the highest spin states of H2 and N2.

Stallcop, James R.↗
NASA NTRS1978

Effective core potentials for the cadmium and mercury atoms

Ab initio effective core potentials have been obtained for the cadmium and mercury atoms by the methods of Kahn et al. (1976). Both two and twelve valence electron representations of Cd and Hg were tested for various atom state-configurations by comparison with all-electron calculations. The generated potentials were used to obtain the equilibrium bond distances and molecular binding energies for the dichloride and dimethyl compounds of both atoms from single and optimum-double configuration self-consistent field calculations.

Basch, H.↗
NASA NTRS1993

Statistical modeling of capture, association, and exit-channel dynamics in the CH3(+)/CH3CN system

The ion-molecule CH3(+) + CH3CN reaction is presently modeled by means of a master equation treatment incorporating weak collisions. The parameter required for the Rice-Ramsberger-Kassel-Markus (RRKM) treatment is derived from an ab initio investigation of the energy surface in question, and a means of including the capture rate coefficients in the RRKM approach is developed in which only the hindered dipole rotation is coupled in the reaction coordinate at large separations. Unimolecular rate coefficients for the (CH3CNCH3+)-activated complex are calculated for all product channels in the 300-600 K temperature range.

Smith, S. C.↗
NASA NTRS1996

Optimized Structures and Proton Affinities of Fluorinated Dimethyl Ethers: An Ab Initio Study

Ab initio methods have been used to investigate the proton affinity and the geometry changes upon protonation for the molecules (CH3)2O, (CH2F)2O, (CHF2)2O, and (CF3)2O. Geometry optimizations were performed at the MP2/3-2 I G level, and the resulting geometries were used for single-point energy MP2/6-31G calculations. The proton affinity calculated for (CH3)2O was 7 Kjoule/mole from the experimental value, within the desired variance of +/- 8Kjoule/mole for G2 theory, suggesting that the methodology used in this study is adequate for energy difference considerations. For (CF3)20, the calculated proton affinity of 602 Kjoule/mole suggests that perfluorinated ether molecules do not act as Lewis bases under normal circumstances; e.g. degradation of commercial lubricants in tribological applications.

Orgel, Victoria B.↗
NASA NTRSDate not supplied

Towards Accurate Predictions of Martensitic Transition Temperatures for Shape Memory Alloys from Ab Initio Simulations

Experimentally, NiTi undergoes a single martensitic phase transition around 341 K from the lowtemperature (T) monoclinic B19’ phase (P21/m) to the high-temperature cubic B2 phase (Pm3m). Theoretically, an orthorhombic B33 (Cmcm) has also been proposed as the T=0 ground state structure, although this phase has never been observed. Accurate predictions of martensitic transition temperatures (MTT) have remained elusive in part due to several well-known theoretical complexities of these systems including low temperature instabilities of the B2 phase. Recently, we proposed a rigorous thermodynamic integration approach based on ab initio simulations to resolve many of these difficulties [1,2]. However, an unsatisfying overprediction of the MTT relative to experiment (by ~100 K) means a fully quantitative theory is still lacking. In this work, we report several new developments to our method that bring first principles theory and experiment much closer into agreement. Our calculations indicate that phonon free energies at low temperature stabilizes B19’ over B33, rationalizing B19’ as the ground state down to T=0. We also find that accurate computations of the electronic free energy, i.e. the change in energy and the appearance of electronic configurational entropy due to finite temperature, is crucial to obtain accurate MTT. Incorporating these corrections results in an MTT prediction of 365 K for binary NiTi, which is in very close agreement with experiment. Our theoretical approach is expected to be a broadly applicable and predictive theory for MTT of SMAs.

Zhigang Wu↗
NASA NTRS1984

Quantum chemical studies of a model for peptide bond formation. 3. Role of magnesium cation in formation of amide and water from ammonia and glycine

The SN2 reaction between glycine and ammonia molecules with magnesium cation Mg2+ as a catalyst has been studied as a model reaction for Mg(2+)-catalyzed peptide bond formation using the ab initio Hartree-Fock molecular orbital method. As in previous studies of the uncatalyzed and amine-catalyzed reactions between glycine and ammonia, two reaction mechanisms have been examined, i.e., a two-step and a concerted reaction. The stationary points of each reaction including intermediate and transition states have been identified and free energies calculated for all geometry-optimized reaction species to determine the thermodynamics and kinetics of each reaction. Substantial decreases in free energies of activation were found for both reaction mechanisms in the Mg(2+)-catalyzed amide bond formation compared with those in the uncatalyzed and amine-catalyzed amide bond formation. The catalytic effect of the Mg2+ cation is to stabilize both the transition states and intermediate, and it is attributed to the neutralization of the developing negative charge on the electrophile and formation of a conformationally flexible nonplanar five-membered chelate ring structure.

Magnesium/chemistry↗
NASA NTRS1987

A high-precision ab initio determination of the equilibrium geometry and force field of HOC(+)

The results of an ab initio molecular orbital investigation of the isoformyl cation, HOC(+), shape are reported. The effects of expanding the basis set to near the Hartree-Fock limit and of electron correlation were examined, and the results indicate that near the Hartree-Fock limit the HOC(+) is linear. An analytic potential function is presented, from which the calculated rotational energies are only 0.03 percent different from the experimental values. This represents a nearly two orders of magnitude reduction in error from earlier work.

Defrees, D. J.↗
NASA NTRS1992

Ab initio predictions on the rotational spectra of carbon-chain carbene molecules

We predict rotational constants for the carbon-chain molecules H2C=(C=)nC, n=3-8, using ab initio computations, observed values for the earlier members in the series, H2CCC and H2CCCC with n=1 and 2, and empirical geometry corrections derived from comparison of computation and experiment on related molecules. H2CCC and H2CCCC have already been observed by radioastronomy; higher members in the series, because of their large dipole moments, which we have calculated, are candidates for astronomical searches. Our predictions can guide searches and assist in both astronomical and laboratory detection.

NASA Discipline Number 52-10↗
NASA NTRSDate not supplied

Metal Hydroxide and Oxyhydroxide Gas Database

In a high temperature combustion situation such as a gas turbine or rocket engine, combustion products include water vapor. This water vapor reacts with oxides in the hot-gas-path to form volatile hydroxides and oxyhydroxides. The oxides may be structural oxide components or protective oxides on alloys or non-oxide ceramics. In either case, this enhanced volatilization leads to degradation of components. Reliable thermochemical data for many metal hydroxides and oxyhydroxides has not been available. We review a database put together at NASA over many years based on experimental and ab initio computations for thermochemical data of Al, Si, Ti, Cr, Mn, Y, Zr, La, Ta, Gd, and Yb hydroxides and oxyhydroxides. Experimental methods for measuring thermodynamic quantities of metal hydroxides and oxyhydroxides are limited since the required high temperatures and high pressures preclude many conventional methods. Nonetheless, transpiration has proven a valuable tool in our laboratories for the Si-O-H, Cr-O-H, and Ti-O-H systems. Data are analyzed with second and third law methods. These data can be leveraged with quantum chemistry to provide reliable free energy functions for the third law calculation and comparison to the derived enthalpies of formation. Both wave function and DFT methods can be utilized for deriving the thermodynamic functions for the various metal hydroxides and oxyhydroxides. We have found that, in general, reasonable molecular geometries and spectroscopic data can be derived with a relatively low level of theory, i.e. DFT with the B3LYP functional. Accurate calculations of the enthalpy of formation require higher levels of theory, such as the Coupled Cluster, Singles, Doubles, and perturbative Triples CCSD(T) approach with complete basis sets (CBS). Also, corrections for hindered rotations and anharmonicity are necessary. For some symmetries, the anharmonic correction has proven to be a challenge and we have found that this is only treated properly by certain codes. The database lists the fHo(298), So(298), and Cp in the form Cp = A + BT + CT2 + D/T + E/T2, where A, B, C, D, and E are the fitting constants. Thus, it can be easily used in any of the common computational thermochemistry codes. To date, the database contains 51 gaseous species of the formula MOx(OH)y.

Metal Hydroxides↗
NASA NTRS1994

Spectral Moments of Collision-Induced Absorption of CO2 Pairs: The Role of the Intermolecular Potential

In this paper we examine the role of the anisotropy of the intermolecular potential in the rototranslational collision-induced absorption of the CO2 pairs. Using newly developed formulas that include the effects of anisotropy of the potential to all orders, we calculate the two lowest spectral moments gamma(prime), and alpha(prime), for four different classes of C02 pair potentials and compare the results with the experimental values. We assumed only multipolar induction in the process of forming the induced dipole, with the second-order contributions included. Using a site-site LJ and a site-site semi-ab initio intermolecular potentials we were able to reproduce the experimental values of gamma(prime), and alpha(prime) moments over entire temperature range from 230 to 330 K. Also, the role of an electrostatic interaction between two C02 molecules and its impact on the spectral moments is thoroughly investigated. An isotropic core with a point quadrupole centered at each molecule is shown to be an inadequate representation of the C02-CO2 potential. Additionally, we show the results obtained with the first- and second-order perturbation theory to be more than twice too small.

Gruszka, Marcin↗
NASA NTRS1991

Use of scaled external correlation, a double many-body expansion, and variational transition state theory to calibrate a potential energy surface for FH2

New ab initio results and a double many-body expansion formalism have been used to parameterize a new FH2 potential energy surface with improved properties near the saddle point and in the region of long-range attraction. The functional form of the new surface includes dispersion forces by a double many-body expansion. Stationary point properties for the new surface are calculated along with the product-valley barrier maxima of vibrationally adiabatic potential curves for F + H2 - HF(nu-prime = 3) + H, F + HD - HF(nu-prime = 3) + D, and F + D2 - DF(nu-prime = 4) + D. The new surface should prove useful for studying the effect on dynamics of a low, early barrier with a wide, flat bend potential.

Lynch, Gillian C.↗
NASA NTRS1995

Ab Initio Studies of Fluorine and Chlorine Oxide and Nitrogen Oxide Species of Interest in Stratospheric Chemistry

The quality of fundamental vibrational frequencies determined using the CCSD(T) method (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations) is shown to be very good, usually predicting band centers to within plus or minus 8 per centimeter. This approach is applied to several molecules of interest in atmospheric chemistry, such as HNO, cis-FONO, cis-ClONO, and ClOOH. The HNO molecule displays a large and unusual anharmonicity in the H-N stretch. For the calculation of ultraviolet (UV) spectra, the linear response CCSD (LRCCSD) approach (which is equivalent to EOM-CCSD) has been shown to yield vertical excitation energies that are accurate to approximately equal to 0.1 eV for singly excited electronic states. This method together with more approximate methods is used to examine the UV spectra of several molecules important in stratospheric chemistry, including HOCl, Cl2O, ClONO2, HONO2, ClOOCl, ClOOH, and HOOH.

Lee, Timothy J.↗
NASA NTRS1995

Accurate ab Initio Quartic Force Fields, Vibrational Frequencies, and Heats of Formation for FCN, FNC, ClCN, and ClNC

The XCN and XNC (X = F, Cl) isomers have been investigated using the CCSD(T) method in conjunction with correlation consistent basis sets. Equilibrium geometries, harmonic frequencies, anharmonic constants, fundamental frequencies, and heats of formation have been evaluated. Agreement with experiment for the fundamental frequencies is very good, even for nu(sub 2), for CICN, which is subject to a strong Fermi resonance with 2nu(sub 3). It is also shown that a second-order perturbation theory approach to solving the nuclear Schroedinger equation gives results in excellent agreement with essentially exact variational calculations. This is true even for nu(sub 2) of ClCN, provided that near-singular terms are eliminated from the perturbation theory formulas and the appropriate Fermi interaction energy matrix is then diagonalized. A band at 615/cm, tentatively assigned as the Cl-N stretch in ClNC in matrix isolation experiments, is shown not to be due to ClNC. Accurate atomization energies are determined and are used to evaluate accurate heats of formation (3.1 +/- 1.5, 33.2 +/- 1.5, 72.6 +/- 1.5, and 75.9 +/- 1.5 kcal/mol for FCN, ClCN, FNC, and ClNC, respectively). It is expected that the theoretical heats of formation for FCN, FNC, and ClNC are the most accurate available.

Lee, Timothy J.↗
NASA NTRS1989

The effect of ionization on the infrared absorption spectra of PAHs: A preliminary report

The emission lines observed in many interstellar IR sources at 3.28, 6.2, 7.7, 8.7, and 11.3 microns are theorized to originate from polycyclic aromatic hydrocarbons (PAHs). These assignments are based on analyses of lab IR spectra of neutral PAHs. However, it is likely that in the interstellar medium that PAHs are ionized, i.e., are positively charged. Besides, as pointed out by Allamandola et al., although the IR emission band spectrum resembles what one might expect from a mixture of PAHs, it does not match in details such as frequency, band profile, or relative intensities predicted from the absorption spectra of any known PAH molecule. One source of more information to test the PAH theory is ab initio molecular orbital theory. It can be used to compute, from first principles, the geometries, vibrational frequencies, and vibrational intensities for model PAH compounds which are difficult to study in the lab. The Gaussian 86 computer program was used to determine the effect of ionization on the infrared absorption spectra of several small PAHs: naphthalene and anthracene. A preliminary report is presented of the results of these calculations.

Defrees, Doug J.↗
NASA NTRS1989

The dissociation energy of He2(+)

Ab initio potential energy curves for He2(+) are computed using full CI wavefunctions in conjuction with large Gaussian and Slater basis sets. The computed curves yield rotational-vibrational transitions with errors of less than 1/cm. They also yield quasi-bound levels with energies and lifetimes in reasonable agreement with translational spectroscopy measurements. The best D(e) is estimated to be 2.470 eV. Green's function Monte Carlo calculations yield a very similar value, 2.466 + or - 0.005.

Bauschlicher, Charles W., Jr.↗