Engineering Papers⌕ Search

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 91 records · Page 5

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.↗

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 vinylacetylene 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 vinylacetylene.

Walch, Stephen P.↗

Characterization of the Minimum Energy Path for the Reaction of Singlet Methylene with N2: The Role of Singlet Methylene in Prompt NO

We report calculations of the minimum energy pathways connecting CH2 + N2 to diazomethane and diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to HCN2+H, and of diazomethane to CH2N+N. The calculations use Complete Active Space Self-Consistent Field (CASSCF) derivative methods to characterize the stationary points and Internally Contracted Configuration Interaction (ICCI) to determine the energetics. The calculations suggest a potential new source of prompt NO from the reaction CH2 with N2 to give diazirine, and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN2, which leads to HCN+N(S-4) on the previously studied CH+N2 surface. The calculations also predict accurate 0 K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, respectively.

Walch, Stephen P.↗

Characterization of the Minimum Energy Path for the Reaction of Singlet Methylene with N2: The Role of Singlet Methylene in Prompt NO

We report calculations of the minimum energy pathways connecting (1)CH2+N2 to diazomethane and diazirine, for the rearrangement of diazirine to diazomethane, for the dissociation of diazirine to HCN2+H, and of diazomethane to CH2N+N. The calculations use complete active space self-consistent field (CASSCF) derivative methods to characterize the stationary points and internally contracted configuration interaction (ICCI) to determine the energetics. The calculations suggest a potential new source of prompt NO from the reaction of (1)CH2 with N2 to give diazirine, and subsequent reaction of diazirine with hydrogen abstracters to form doublet HCN2, which leads to HCN+N(S-4) on the previously studied CH+N2 Surface. The calculations also predict accurate 0 K heats of formation of 77.7 kcal/mol and 68.0 kcal/mol for diazirine and diazomethane, respectively.

Walch, Stephen P.↗

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.↗

Theoretical Study of the Electric Dipole Moment Function of the CIO Molecule

The potential energy function and electric dipole moment function (EDMF) are computed for CIO Chi(sup 2)Pi using several different techniques to include electron correlation. The EDMF is used to compute Einstein coefficients, vibrational lifetimes, and dipole moments in higher vibrational levels. Remaining questions concerning the position of the maximum of the EDMF may be resolved through experimental measurement of dipole moments of higher vibrational levels. The band strength of the 1-0 fundamental transition is computed to be 12 +/- 2 /sq cm atm in good agreement with three experimental values, but larger than a recent value of 5 /sq cm atm determined from infrared heterodyne spectroscopy. The theoretical methods used include SCF, CASSCF, multireference singles plus doubles configuration interaction (MRCI) and contracted CI, coupled pair functional (CPF), and a modified version of the CPF method. The results obtained using the different methods are critically compared.

Pettersson, Lars G. M.↗

Computed Minimum Energy Pathway for Isomerization in Ketene

A minimum energy pathway for interchange of the CH groups in ketene via a C2v structure has been obtained using complete active space self consistent field (CASSCF) derivative methods with a polarized valence double zeta basis set to define the reaction pathway followed by multi-reference internally contracted configuration interaction (ICCI) calculations with a [3s3p2d/3s2p] basis set to determine the energetics. Qualitatively, the C2v structure is found to be a shallow minimum on the potential energy surface separated from ketene by a small barrier (0.2 kcal/mol), a second minimum, and a larger barrier (3.0 kcal/mol). The minimum energy pathway leading from the C2v minimum to ketene starts by simultaneous rotation of the farther CH group out of the plane and away from the oxygen followed by increase of the CCO angle and subsequent 1,2-migration of the H of the nearer CH group toward the carbon of the farther CH group.

Walch, Stephen P.↗

The Theoretical Transition Probabilities Between the B(sup 3)Pi(sub g) and the A(sup 3)Sigma(Sup +, sub u), W(sup 3)Delta(sub u), B'(sup 3)Sigma(sup -, sub u) States of N2

The electronic transition moment functions between the B(sup 3)Pi(sub g) and the A(sup 3)Sigma(sup +, sub u), W(sup 3)Delta(sub u), B'(sup 3)Sigma(sup -, sub u) states of N2 are studied using the internally contracted multireference configuration interaction (ICMRCI) method based upon complete active space SCF (CASSCF) reference wave-functions. The dependence of the moments on both the one and n-particle basis sets has been investigated in detail. The calculated radiative lifetimes for the vibrational levels of B(sup 3)Pi(sub g) are in excellent agreement with the most recent measurement of Euler and Pipkin (1983)

Thuemmel, Helmar T.↗

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.↗

On the Reaction of FNO2 with CH3, t-butyl, and C13H21

Theoretical studies are reported for the reaction of FN02 with the radicals CH3, t-butyl, and C13H21, which are templates for the radical site of a hydrogenated diamond surface. All structures axe fully optimized using density functional theory (DFT) based on the B3LYP functional. Calibration calculations axe performed for CH3 + FNO2 using the coupled cluster approach, the the internally contracted multireference configuration. interaction method, and second order perturbation theory based upon complete active space SCF (CASSCF) reference wave function. These calibration calculations support the B3LYP approach for the calculation of bond energies, but show the B3LYP barrier is too low. Combining the calibration calculations with the larger clusters yields our best estimate of a barrier of about 10 kcal/mol. for the reaction of FNO2 with a radical site on hydrogenated diamond.

Thuemmel, H. T.↗

Comparison of DFT with Traditional Methods for the Calculation of Vibrational Frequencies and Bond Energies

The vibrational frequencies of MO2 are computed at many levels of theory, including HF, B3LYP, BP86, CASSCF, MP2, and CCSD(T). The computed results are compared with the available experimental results. Most of the methods fail for at least one state of the systems considered. The accuracy of the results and the origin of the observed failures are discussed. The B3LYP bond energies are compared with traditional methods for a variety of systems, ranging from FeCOn+ to SiCln and its positive ions. The cases where B3LYP differs from the traditional methods are discussed.

Bauschlicher, Charles W.↗

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.↗