Engineering PapersSearch

SEARCH · Engineering Papers

Results for “CASSCF”

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 55 records · Page 3

Theoretical investigations of the structures and binding energies of Be(sub n) and Mg(sub n) (n = 3-5) clusters

Researchers determined the equilibrium geometries and binding energies of Be and Mg trimers, tetramers and pentamers using single and double excitation coupled cluster (CCSD) and complete active space self-consistent-field (CASSCF) multireference configuration interaction (MRCI) wave functions in conjunction with extended atomic basis sets. Best estimates of the cluster binding energies are 24, 83 and 110 kcal/mole for Be3, Be4 and Be5; and 9, 31 and 41 kcal/mole for Mg3, Mg4 and Mg5, respectively. A comparison of the MRCI and CCSD results shows that even the best single-reference approach (limited to single and double excitations) is not capable of quantitative accuracy in determining the binding energies of Be and Mg clusters.

Lee, Timothy J.

Theoretical investigations of the structures and binding energies of Be(n) and Mg(n) (n = 3-5) clusters

The equilibrium geometries and binding energies of Be and Mg trimers, tetramers and pentamers have been determined using single and double excitation coupled cluster (CCSD) and complete active space self-consistent-field (CASSCF) multireference configuration interaction (MRCI) wave functions in conjunction with extended atomic basis sets. The best estimates of the cluster binding energies are 24, 83, and 110 kcal/mol for Be3, Be4, and Be5; and 9, 31, and 41 kcal/mol for Mg3, Mg4, and Mg5, respectively. A comparison of the MRCI and CCSD results shows that even the best single-reference approach (limited to single and double excitations) is not capable of quantitative accuracy in determining the binding energies of Be and Mg clusters.

Lee, Timothy J.

An ab initio study of the C3(+) cation using multireference methods

The energy difference between the linear 2 sigma(sup +, sub u) and cyclic 2B(sub 2) structures of C3(+) has been investigated using large (5s3p2d1f) basis sets and multireference electron correlation treatments, including complete active space self consistent fields (CASSCF), multireference configuration interaction (MRCI), and averaged coupled-pair functional (ACPF) methods, as well as the single-reference quadratic configuration interaction (QCISD(T)) method. Our best estimate, including a correction for basis set incompleteness, is that the linear form lies above the cyclic from by 5.2(+1.5 to -1.0) kcal/mol. The 2 sigma(sup +, sub u) state is probably not a transition state, but a local minimum. Reliable computation of the cyclic/linear energy difference in C3(+) is extremely demanding of the electron correlation treatment used: of the single-reference methods previously considered, CCSD(T) and QCISD(T) perform best. The MRCI + Q(0.01)/(4s2p1d) energy separation of 1.68 kcal/mol should provide a comparison standard for other electron correlation methods applied to this system.

Taylor, Peter R.

Theoretical study of the low-lying bound states of O2

It is demonstrated that a complete-active-space self-consistent-field (CASSCF) (2p)/MRCI + Q (multireference configuration interaction with a Davidson correction) description in a (13s8p6d 4f2g)/((5s4p3d 2f1g) atomic natural orbits (ANO) basis set supplemented with diffuse functions provides a quantitative description of the six lowest states of O2. The calculated potentials are within 0.05 eV (1.2 kilocal/mol) of accurate experimental results. The importance of substantially expanding the primitive basis set has been investigated, and it is demonstrated that such expansions yield insignificant improvement in the spectroscopic constants. Potential energy curves have also been reported for the weakly bound states of O2. The 5Pi(g) state is estimated to have a D(e) of 0.16 +/- 0.03 eV. The upper bound of D(e) is found to be sufficiently large that the importance of this state as a precursor for the formation of O2 (b 1Sigma(t)(+)) and O(1S) should be reconsidered.

Partridge, Harry

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

Effort continued through this period to refine and expand the SIRIUS/ABACUS program package for CASSCF and RASSCF second derivatives. A new approach to computing the Gaussian integral derivatives that require much of the time in gradient and Hessian calculations was devised. Several different studies were undertaken in the area of application calculations. These include a study of proton transfer in the HF trimer, which provides an analog of rearrangement reactions, and the extension of our previous work on Be and Mg clusters to Ca clusters. In addition, a very accurate investigation of the lowest-lying potential curves of the O2 molecule was completed. These curves are essential for evaluating different models of the terrestrial atmosphere nightglow. The most effort this year was devoted to a large scale investigation of stationary points on the C4H4 surface, and the thermochemistry of acetylene/acetylene reaction.

Taylor, Peter

The vibrational spectrum of H2O3: An ab initio investigation

Theoretically determined frequencies and absorption intensities are reported for the vibrational spectrum of the covalent HOOOH and hydrogen bonded HO---HOO intermediates that may form in the reaction of the hydroxyl and hydroperoxyl radicals. Basis sets of DZP quality, augmented by diffuse and second sets of polarization functions have been used with CASSCF wave functions. The calculated harmonic vibrational frequencies of HOOOH have been corrected with empirical factors and presented in the form of a 'stick' spectrum. The oxygen backbone vibrations, predicted to occur at 519, 760, and 870 cm(exp -1), are well separated from most interferences, and may be the most useful for the species' identification. In the case of the hydrogen bonded isomer, emphasis has been placed upon prediction of the shifts in the intramolecular vibrational frequencies that take place upon formation of the complex. In particular, the HO stretch and HOO bend of HO2 are predicted to have shifts of -59 and 53 cm(exp -1), respectively, which should facilitate their identification. It is also noted that the antisymmetric stretching frequency of the oxygen backbone in HOOOH exhibits a strong sensitivity to the degree of electron correlation, such as has been previously observed for the same mode in ozone.

Jackels, Charles F.

Theoretical characterization of the potential energy surface for NH + NO

The potential energy surface (PES) for NH + NO was characterized using complete active space self-consistent field (CASSCF) gradient calculations to determine the stationary point geometries and frequencies followed by CASSCF/internally contracted configuration interaction (CCI) calculations to refine the energetics. The present results are in qualitative accord with the BAC-MP4 calculations, but there are differences as large as 8 kcal/mol in the detailed energetics. Addition of NH to NO on a (2)A' surface, which correlated with N2 + OH or H + N2O products, involves barriers of 3.2 kcal/mol (trans) and 6.3 kcal/mol (cis). Experimental evidence for these barriers is found in earlier works. The (2)A' surface has no barrier to addition, but does not correlate with products. Surface crossings between the barrierless (2)A' surface and the (2)A' surface may be important. Production of N2 + OH products is predicted to occur via a planar saddle point of (2)A' symmetry. This is in accord with the preferential formation of II(A') lambda doublet levels of OH in earlier experiments. Addition of NH (1)delta to NO is found to occur on an excited state surface and is predicted to lead to N2O product as observed in earlier works.

Walch, Stephen P.

Theoretical characterization of the reaction CH3 +OH yields CH3OH yeilds products: The (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO channels

The potential energy surface (PES) for the CH3OH system has been characterized for the (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO product channels using complete-active-space self-consistent-field (CASSCF) gradient calculations to determine the stationary point geometries and frequencies followed by CASSCF/internally contracted configuration-interaction (CCI) calculations to refine the energetics. The (1)CH2 + H2O channel is found to have no barrier. The long range interaction is dominated by the dipole-dipole term, which orients the respective dipole moments parallel to each other but pointing in opposite directions. At shorter separations there is a dative bond structure in which a water lone pair donates into the empty a" orbital of CH2. Subsequent insertion of CH2 into an OH bond of water have barriers located at -5.2 kcal/mol and 1.7 kcal/mol, respectively, with respect to CH3 + OH. From comparison of the computed energetics of the reactants and products to known thermochemical data it is estimated that the computed PES is accurate to plus or minus 2 kcal/mol.

Walch, Stephen P.

Computed barrier heights for H + CH2O yields CH3O yields CH2OH

The barrier heights (including zero-point effects) for H + CH2O yields CH3O and CH3O yields CH2OH have been computed using complete active space self consistent field (CASSCF)/gradient calculations to define the stationary point geometries and harmonic frequencies and internally contracted configuration-interaction (CCI) to refine the energetics. The computed barrier heights are 5.6 kcal/mol and 30.1 kcal/mol, respectively. The former barrier height compares favorably to an experimental activation energy of 5.2 kcal/mol.

Walch, Stephen P.

Theoretical characterization of the potential energy surface for NH + NO

The potential energy surface for NH + NO was characterized using complete active space self-consistent field (CASSCF) gradient calculation to determine the stationary point geometries and frequencies followed by CASSCF/internally contracted configuration interaction calculations to refine the energetics. The present results are in qualitative accord with the BAC-MP4 calculations, but there are differences as large as 8 kcal/mol in the detailed energetics.

Walch, Stephen P.

Computed barrier heights for H + CH2O yields CH3O yields CH2OH

The barrier heights (including zero-point effects) for H + CH2O yields CH3O and CH3O yields CH2OH have been computed using complete active space self consistent field (CASSCF)/gradient calculations to define the stationary point geometries and harmonic frequencies and internally contracted configuration-interaction (CCI) to refine the energetics. The computed barrier heights are 5.6 kcal/mol and 30.1 kcal/mol, respectively. The former barrier height compares favorably to an experimental activation energy of 5.2 kcal/mol.

Walch, Stephen P.

Theoretical characterization of the reaction CH3 + OH yields CH3OH yields products - The (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO channels

The potential energy surface (PES) for the CH3OH system has been characterized for the (1)CH2 + H2O, H2 + HCOH, and H2 + H2CO product channels using complete-active-space self-consistent-field (CASSCF) gradient calculations to determine the stationary point geometries and frequencies followed by CASSCF/internally contracted configuration-interaction (CCI) calculations to refine the energetics. The (1)CH2 + H2O channel is found to have no barrier. The long range interaction is dominated by the dipole-dipole term, which orients the respective dipole moments parallel to each other but pointing in opposite directions. At shorter separations there is a dative bond structure in which a water lone pair donates into the empty 'a' orbital of CH2. Subsequent insertion of CH2 into an OH bond of water have barriers located at -5.2 kcal/mol and 1.7 kcal/mol, respectively, with respect to CH3 + OH. From comparison of the computed energetics of the reactants and products to known thermochemical data it is estimated that the computed PES is accurate to plus or minus 2 kcal/mol.

Walch, Stephen P.

The protonation of N2O reexamined - A case study on the reliability of various electron correlation methods for minima and transition states

The protonation of N2O and the intramolecular proton transfer in N2OH(+) are studied using various basis sets and a variety of methods, including second-order many-body perturbation theory (MP2), singles and doubles coupled cluster (CCSD), the augmented coupled cluster (CCSD/T/), and complete active space self-consistent field (CASSCF) methods. For geometries, MP2 leads to serious errors even for HNNO(+); for the transition state, only CCSD/T/ produces a reliable geometry due to serious nondynamical correlation effects. The proton affinity at 298.15 K is estimated at 137.6 kcal/mol, in close agreement with recent experimental determinations of 137.3 +/- 1 kcal/mol.

Martin, J. M. L.

Theoretical characterization of the reaction NH2 + O yields products

The potential energy surface for NH2+O has been characterized using complete active space self-consistent field (CASSCF)/derivative calculations to determine stationary point geometries and frequencies followed by internally contracted configuration interaction (ICCI) calculations to determine the energetics. The calculations predict a NO bond strength of 85.8 kcal/mol for NH2O. The barrier for isomerization of NH2O to trans-HNOH is predicted to be 48.0 kcal/mol and the barriers for H+HNO forming NH2O and NHOH are predicted to be 2.1 and 8.3 kcal/mol, respectively (all corrected for zero-point energy). The computed heats of formation for NH2O and cis- and trans-HNOH are in good agreement with the present results. The barrier for H + HNO yields H2 + NO is computed to be about 0.3 kcal/mol.

Walch, Stephen P.

Theoretical characterization of the reaction NH2+NO - products

The potential energy surface (PES) for NH2+NO is characterized using complete active space self-consistent-field (CASSCF)/derivative methods to locate the stationary points, followed by internally contracted configuration interaction (ICCI) calculations to determine the energetics. Production of N2+H2O is found to involve a complex mechanism, which, however, has no barrier with respect to NH2+NO. This pathway is exothermic by 124.5 kcal/mol. Production of NH2+OH can occur with no barrier other than the exothermicity from any of three isomers of HNNOH.

Walch, Stephen P.

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(C0)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule 04(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in 04(+) has also been studied using the Brueckner coupled cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.

Computational studies of metal-metal and metal-ligand interactions

Accurate calculations on the bond length, dipole moment, and harmonic frequency of CO are presented, using large basis sets and high levels of electron correlation. The geometric structure, force constants and binding energies of Cr(CO)6 and Cr(CO)5 are computed using large basis sets and high levels of electron correlation. The molecule O4(+) is studied, using large basis sets and high levels of electron correlation, including the CASSCF, CASSI,and CASPT2 methods. Binding energies, geometries and frequencies are computed. Symmetry breaking is a particular problem for the antisymmetric stretch, which is addressed using the CASSI method. The symmetry breaking problem in O4(+) has also been studied using the Brueckner coupled-cluster method. This gives results in good agreement with CASSI. A multi-region numerical integration scheme is investigated for use in Density Functional Calculations. This scheme is found to give comparable results to a widely used scheme based on the Euler-Maclaurin technique.

Barnes, Leslie A.

Characterization of the Minimum Energy Paths and Energetics for the Reaction of Vinylidene with Acetylene

The reaction of vinylidene (CH2C) with acetylene may be an initiating reaction in soot formation. We report minimum energy paths and accurate energetics for a pathway leading to vinyl-acetylene and for a number of isomers of C4H4. The calculations use complete active space self-consistent field (CASSCF) derivative methods to characterize the stationary points and internally contacted configuration interaction (ICCI) and/or coupled cluster singles and doubles with a perturbational estimate of triple excitations (CCSD(T)) to determine the energetics. We find an entrance channel barrier of about 5 kcal/mol for the addition of vinylidene to acetylene, but no barriers above reactants for the reaction pathway leading to vinyl-acetylene.

Walch, Stephen P.