Search NASA⌕ Search

SEARCH · Search NASA

Results for “Basis sets”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Ab initio characterization of ClOOH - Implications for atmospheric chemistry

The equilibrium structure, dipole moment, harmonic vibrational frequencies, and infrared intensities of ClOOH are determined using the CCSD(T) (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations) electronic structure method in conjunction with a TZ2P (triple xi plus double polarization) basis set. The heat of formation of CIOOH is determined (using two different isodesmic reactions) to be +1.5 +/- 1 kcal/mol at 0 K or +0.2 +/- 1 kcal/mol at 298.15 K. Using the computed heat of formation, we examined the stability of ClOOH with respect to the ClO + OH, ClOO + H, and HOO + Cl dissociation limits. Since ClOOH is found to be quite stable, it is argued that the chemistry of ClOOH should be included in any accurate modeling of the stratosphere.

Lee, Timothy J.↗
NASA NTRS1992

A global ab initio potential for HCN/HNC, exact vibrational energies, and comparison to experiment

An ab initio (i.e., from first principles) calculation of vibrational energies of HCN and HNC is reported. The vibrational calculations were done with a new potential derived from a fit to 1124 ab initio electronic energies which were calculated using the highly accurate CCSD(T) coupled-cluster method in conjunction with a large atomic natural orbital basis set. The properties of this potential are presented, and the vibrational calculations are compared to experiment for 54 vibrational transitions, 39 of which are for zero total angular momentum, J = 0, and 15 of which are for J = 1. The level of agreement with experiment is unprecedented for a triatomic with two nonhydrogen atoms, and demonstrates the capability of the latest computational methods to give reliable predictions on a strongly bound triatomic molecule at very high levels of vibrational excitation.

Bentley, Joseph A.↗
NASA NTRS1993

Calculation of molecular excitation rates

State-to-state collisional excitation rates for interstellar molecules observed by radio astronomers continue to be required to interpret observed line intensities in terms of local temperatures and densities. A problem of particular interest is collisional excitation of water which is important for modeling the observed interstellar masers. In earlier work supported by a different NASA Grant, excitation of water in collisions with He atoms was studied; after many years of successively more refined calculations that problem now seems to be well understood, and discrepancies with earlier experimental data for related (pressure broadening) phenomena are believed to reflect experimental errors. Because of interstellar abundances, excitation by H2, the dominant interstellar species, is much more important than excitation by He, although it has been argued that rates for excitation by these are similar. Under the current grant theoretical study of this problem has begun which is greatly complicated by the additional degrees of freedom which must be included both in determining the interaction potential and also in the molecular scattering calculation. We have now computed the interaction forces for nearly a thousand molecular geometries and are close to having an acceptable global fit to these points which is necessary for the molecular dynamics calculations. Also, extensive modifications have been made to the molecular scattering code, MOLSCAT. These included coding the rotational basis sets and coupling matrix elements required for collisions of an asymmetric top with a linear rotor. A new method for numerical solution of the coupled equations has been incorporated. Because of the long-ranged nature of the water-hydrogen interaction it is necessary to integrate the equations to rather large intermolecular separations, and the integration methods previously available in MOLSCAT are not ideal for such cases. However, the method used by Alexander in his HIBRIDON code is particularly suited for such cases. We have obtained this code and incorporated that part which solves the coupled differential equations as an option in the MOLSCAT program.

Flynn, George↗
NASA NTRS1993

Theoretical research program to study chemical reactions in AOTV bow shock tubes

The main focus was the development, implementation, and calibration of methods for performing molecular electronic structure calculations to high accuracy. These various methods were then applied to a number of chemical reactions and species of interest to NASA, notably in the area of combustion chemistry. Among the development work undertaken was a collaborative effort to develop a program to efficiently predict molecular structures and vibrational frequencies using energy derivatives. Another major development effort involved the design of new atomic basis sets for use in chemical studies: these sets were considerably more accurate than those previously in use. Much effort was also devoted to calibrating methods for computing accurate molecular wave functions, including the first reliable calibrations for realistic molecules using full CI results. A wide variety of application calculations were undertaken. One area of interest was the spectroscopy and thermochemistry of small molecules, including establishing small molecule binding energies to an accuracy rivaling, or even on occasion surpassing, the experiment. Such binding energies are essential input to modeling chemical reaction processes, such as combustion. Studies of large molecules and processes important in both hydrogen and hydrocarbon combustion chemistry were also carried out. Finally, some effort was devoted to the structure and spectroscopy of small metal clusters, with applications to materials science problems.

Taylor, Peter R.↗
NASA NTRS1994

(H2)2 mole-fraction altitude profile in the atmosphere of Jupiter: A computational study

The mole fraction x(sub 2) of (H2)2 in equilibrium mixture with H2 under the atmospheric conditions of Jupiter is evaluated from the dimerization equilibrium constant calculated by quantum-chemical treatments and also from the Lennard-Jones potential. The treatments are of an ab initio type with the second and fourth order Moller-Plesset perturbation techniques and a basis set superposition error evaluation. The computed dimerization equilibrium constant is combined with observed height profiles of temperature and pressure. In six treatments considered it is found that the mole fraction decreases with increasing height. Various approximations suggest the dimeric mole fraction at the Jupiter 1 atm pressure level between 0.04 and 1.06%.

Slanina, Zdenek↗
NASA NTRS1993

Applicability of PM3 to transphosphorylation reaction path: Toward designing a minimal ribozyme

A growing body of evidence shows that RNA can catalyze many of the reactions necessary both for replication of genetic material and the possible transition into the modern protein-based world. However, contemporary ribozymes are too large to have self-assembled from a prebiotic oligonucleotide pool. Still, it is likely that the major features of the earliest ribozymes have been preserved as molecular fossils in the catalytic RNA of today. Therefore, the search for a minimal ribozyme has been aimed at finding the necessary structural features of a modern ribozyme (Beaudry and Joyce, 1990). Both a three-dimensional model and quantum chemical calculations are required to quantitatively determine the effects of structural features of the ribozyme on the reaction it catalyzes. Using this model, quantum chemical calculations must be performed to determine quantitatively the effects of structural features on catalysis. Previous studies of the reaction path have been conducted at the ab initio level, but these methods are limited to small models due to enormous computational requirements. Semiempirical methods have been applied to large systems in the past; however, the accuracy of these methods depends largely on a simple model of the ribozyme-catalyzed reaction, or hydrolysis of phosphoric acid. We find that the results are qualitatively similar to ab initio results using large basis sets. Therefore, PM3 is suitable for studying the reaction path of the ribozyme-catalyzed reaction.

Manchester, John I.↗
NASA NTRS1994

An ab initio study of HCuCO

HCuCO is studied using a large Gaussian basis set at the coupled cluster singles and doubles level of theory, including a perturbational estimate of the connected triples (CCSD(T)). In contrast with CuCO, HCuCO is linear. The Cu-CO bond in HCuCO is significantly stronger than in CuCO. These differences between HCuCO and CuCO are discussed in terms of theCu-H bond polarizing the Cu 4s electron away from the CO.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1994

The Fraternal Twins of Quartet O4(+)

Eleven stationary geometries of quartet O4(+) have been studied by ab initio methods. The geometries were optimized at the Complete Active Space Self-Consistent Field (CASSCF) level of theory and the energies were calculated by the multiconfigurational second order pertubation method (CASPT2), using Double-Zeta Plus polarization (DZP), Triple-Zeta Plus Double Polarization (TZ2P), average Atomic Natural Orbital (ANO) (5s4p2d) and average ANO (6s5p3d2f) basis sets. The rectangular and trans-planar structures are found to be the most stable, with an energy barrier to conversion between the two at the threshold of dissociation. Both have a delocalized hole and are stable relative to separated 02 and 02(+)by 11.0 and 11.5 kcal/mol for the rectangular and the 2 trans-planar structure, respectively, compared with the experimentally deduced energy in the range of 9.2 to 10.8 kcal/mol. The adiabatic ionization potentials of 04 and 02 are computed to be 11.67 and 12.21 eV, while experimental values are 11.66 and 12.07 eV, respectively. The vibrational frequencies have been computed for all degrees of freedom at the CASSCF level of theory. Symmetry breaking is found to be a particular problem in the computation of the antisymmetric stretch frequency for the delocalized structures at the CASSCF level of theory. Attempts to rectify these problems using the restricted active space self-consistent field (RASSCF) method leads to additional difficulties, but further analysis yields insight into the symmetry breaking and problems with earlier calculations. Finally, a nonorthogonal configuration interaction (CI) calculation based on the interaction of localized CASSCF wave functions using the Complete Active Space State Interation (CASSI) method leads to a balanced treatment of the antisymmetric stretch which is free from symmetry breaking. The study explains the four most prominent absorption frequencies observed in the partially unassigned IR spectrum of O4(+) isolated in solid neon as the antisymmetric OO stretch, and the combination band of the symmetric and antisymmetric 00 stretch of both the rectangular and trans-planar structures.

Lindh, Roland↗
NASA NTRS1996

Use of Density Functional Method to Study Molecular Vibrations

Use of a scale factor reduces the errors in computed vibrational frequencies (errors due to anharmonicities, basis set deficiencies, and approximate treatment of electron correlation). Additional errors come from modelling solids with small molecules. Nevertheless, we hope that the present calculations add to the understanding of surface IR measurements of the heat tiles.

Chong, Delano P.↗
NASA NTRS1994

Coupled-Cluster Theory Employing Approximate Integrals: An Approach to Avoid the Input/Output and Storage Bottlenecks

By representing orbital products in an expansion basis, certain classes of two-electron integrals are approximated for use in CCSD(T) calculations (singles and doubles coupled-cluster plus a perturbational estimate of the effects of connected triple excitations). This leads to a very large reduction in disk storage and input/output requirements, with usually only a modest increase in computational effort. The new procedure will allow very large CCSD(T) calculations to be undertaken, limited only by available processor time. Using the molecular basis as the expansion basis, explicit numerical comparisons of equilibrium geometries, harmonic frequencies, and energy differences indicate that the error due to the use of approximate integrals is less than the error associated with truncation of the molecular basis set.

Rendell, Alistair P.↗
NASA NTRS1996

AB Initio Characterization of MgCCH, MgCCH(+), and MgC2, and Pathways to their Formation in the Interstellar Medium

A study of Mg-bearing compounds has been performed in order to determine molecular properties which are critical for planning new astronomical searches and laboratory studies. The primary focus of the work is on MgCCH, MgCCH(+), and the isomers of MgC2. Only MgCCH has been identified in laboratory studies. Additional calculations have been carried out on MgH, MgNC, MgCN, and their cations in an effort to evaluate pathways to the formation of MgCCH and MgCCH(+) in the InterStellar Medium (ISM) or in circumstellar envelopes. Correlated ab initio methods and correlation-consistent basis sets have been employed. Properties including structures, rotational constants, dipole moments, and harmonic frequencies are reported. A transition state between linear MgCC and cyclic MgC2 has been characterized and was found to yield a minimal barrier (approx. 0.5 kcal/mole), indicating easy interconversion to the cyclic form. Direct reactions in the ISM between Mg or Mg(+) and HCCH are precluded by energetic considerations, but a number of ion- molecule or neutral-neutral exchange reactions between CCH and various Mg-containing species offer plausible pathways to MgCCH or MgCCH(+). Weakly bound MgH may react with CCH to form MgCCH, but MgH has not been detected. Both MgNC and MgCN have been observed, but reactions with CCH are slightly endothermic by 1-3 kcal/mole. Although MgH(+), MgNC(+), and MgCN(+) have not been detected, their reactions with CCH to form MgCCH(+) are all exothermic. With only a small barrier separating linear MgCC and cyclic MgC2, the dissociative recombination of MgCCH(+) with an electron is expected to yield cyclic MgC2, and regenerate Mg and CCH. New astronomical searches for MgCCH, MgCCH(+), cyclic MgC2, MgNC(+), and MgCN(+) will provide further insight into organo-magnesium astrochemistry.

Woon, David E.↗
NASA NTRS1996

Ab initio Potential-Energy Surfaces and Electron-Spin-Exchange Cross Sections for H-O2 Interactions

Accurate quartet- and doublet-state potential-energy surfaces for the interaction of a hydrogen atom and an oxygen molecule in their ground states have been determined from an ab initio calculation using large-basis sets and the internally contracted multireference configuration interaction method. These potential surfaces have been used to calculate the H-O2 electron-spin-exchange cross section; the square root of the cross section (in a(sub 0)), not taking into account inelastic effects, can be obtained approximately from the expressions 2.390E(sup -1/6) and 5.266-0.708 log10(E) at low and high collision energies E (in E(sub h)), respectively. These functional forms, as well as the oscillatory structure of the cross section found at high energies, are expected from the nature of the interaction energy. The mean cross section (the cross section averaged over a Maxwellian velocity distribution) agrees reasonably well with the results of measurements.

Stallcop, James R.↗
NASA NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO rev. reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X=H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate(DELTA)H(sup o)(sub f) value for HNO using state of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

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 NTRS1995

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO reversible reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HNO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) value for HNO using state-of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +/- 0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗
NASA NTRS1995

Characterization of the Minimum Energy Paths for the Ring Closure Reactions of C4H3 with Acetylene

The ring closure reaction of C4H3 with acetylene to give phenyl radical is one proposed mechanism for the formation of the first aromatic ring in hydrocarbon combustion. There are two low-lying isomers of C4H3; 1-dehydro-buta-l-ene-3-yne (n-C4H3) and 2-dehydro-buta-l-ene-3-yne (iso-C4H3). It has been proposed that only n-C4H3 reacts with acetylene to give phenyl radical, and since iso-C4H3 is more stable than n-C4H3, formation of phenyl radical by this mechanism is unlikely. We report restricted Hartree-Fock (RHF) plus singles and doubles configuration interaction calculations with a Davidson's correction (RHF+1+2+Q) using the Dunning correlation consistent polarized valence double zeta basis set (cc-pVDZ) for stationary point structures along the reaction pathway for the reactions of n-C4H3 and iso-C4H3 with acetylene. n-C4H3 plus acetylene (9.4) has a small entrance channel barrier (17.7) (all energetics in parentheses are in kcal/mol with respect to iso-C4H3 plus acetylene) and the subsequent closure steps leading to phenyl radical (-91.9) are downhill with respect to the entrance channel barrier. Iso-C4H3 Plus acetylene also has an entrance channel barrier (14.9) and there is a downhill pathway to 1-dehydro-fulvene (-55.0). 1-dehydro-fulvene can rearrange to 6-dehydro-fulvene (-60.3) by a 1,3-hydrogen shift over a barrier (4.0), which is still below the entrance channel barrier, from which rearrangement to phenyl radical can occur by a downhill pathway. Thus, both n-C4H3 and iso-C4H3 can react with acetylene to give phenyl radical with small barriers.

Walch, Stephen P.↗
NASA NTRS1986

Benchmark Full Configuration-Interaction Calculations on HF and NH2

Full configuration-interaction (FCI) calculations are performed at selected geometries for the (sup 1)Sigma(+) state of HCF and the (sup 2)B(sub 1) and (sup 3)A(sub 1) states NH2 using both DZ and DZP Gaussian basis sets. Higher excitations become more important when the bonds are stretched and the SCF reference becomes a poorer zeroth-order description of the wave function. The CASSCF-MRCI procedure gives excellent agreement with the FCI potentials, especially when corrected with a multireference analog of the Davidson correction.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1986

On the Bonding in Fe2(CO)9

The bonding in Fe2(CO)9 is analyzed using an self consistend field (SCF) wave function for a large basis set. There is no direct Fe-Fe metal-metal bond. The bridging CO's hold the two Fe(CO)3 fragments together by a sigma donation into the empty Fe-Fe d pi orbital and metal donation from the d pi* orbital into the CO 2pi* orbital. The bonding of the terminal CO is similar to that in Ni(CO)4 and the equatorial groups in Fe(CO)5.

Bauschlicher, Charles W., Jr.↗