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 73 records · Page 4

Effective core potential parameters for first- and second-row atoms

An improved effective core potential (ECP) technique is described and used to give ECP parameters for the atoms of the first two rows of the periodic table. A given basis set is parametrized which allows for a direct comparison with all-electron calculations. Extensive test calculations on first- and second-row molecules using the ECP have been performed, giving excellent agreement with the all-electron results at the SCF, CASSCF, and CI levels of treatment. Correlating and diffuse functions may be added without modifying the ECP parameters. The present ECP descriptions result in CPU time reductions of the order of 50 percent in addition to reduced disk storage.

Pettersson, Lars G. M.↗

An initio calculations on C2, Si2, and SiC

Full-configuration-interaction (FCI) calculations have been performed for the Z 1Sigma(g)(+), a 3Pi(u), and b 3Sigma(g)(-) states of C2 and the X 3Sigma(g)(-) and A 3Pi(u) states of Si2 in valence double-zeta plus polarization one-particle basis sets. By carrying out the CASSCF multireference singles plus doubles CI (MRCI) procedure in a very large Gaussian basis set and correcting for errors with the FCI, the triplet-triplet separation in C2 is reproduced to within 47/cm. Analogous calculations on Si2 definitively determine the ground state to be 3Sigma(g)(-); the 3Pi(u) state is only 440 + or - 100/cm higher in energy. For the SiC molecule, which has never been observed spectroscopically but is believed to be an important component in carbon stars and to be present in interstellar space, the A 3Sigma(-) state is nearly 4000/cm above the X 3Pi ground state.

Bauschlicher, Charles W., Jr.↗

Theoretical study of the scandium and yttrium halides

The X1Sigma(+) ground states and a3Delta excited states of the diatomic halides of Sc and Y are characterized theoretically, using the SDCI coupled-pair functional method and the state-averaged CASSCF method to determine the spectroscopic constants and related properties. The techniques employed are discussed, and the results are presented in extensive tables. The dissociation energies are given as D0 = 6.00 eV for ScF, 4.55 eV for ScCl, 3.90 eV for ScBr, 6.72 eV for YF, 5.36 eV for YCl, and 4.74 eV for YBr.

Langhoff, Stephen R.↗

The computed spectrum of AlC

The low-lying (up to about 40,000/cm) doublet and quartet states of AlC have been studied at the CASSCF and MRCI levels of theory. The bonding in the X4Sigma(-) state is shown to involve three one-electron bonds. The first excited state, a 2Pi, lies about 8700/cm higher and like the X4Sigma(-) state is derived from ground state Al and C, but has a two-electron pi bond leading to a much shorter bond length.

Bauschlicher, Charles W., Jr.↗

Theoretical characterization of the minimum energy path for hydrogen atom addition to N2 - Implications for the unimolecular lifetime of HN2

Results are reported from CASSCF externally contracted CI ab initio computations of the minimum-energy path for the addition of H to N2. The theoretical basis and numerical implementation of the computations are outlined, and the results are presented in extensive tables and graphs and characterized in detail. The zero-point-corrected barrier for HN2 dissociation is estimated as 8.5 kcal/mol, and the lifetime of the lowest-lying quasi-bound vibrational state of HN2 is found to be between 88 psec and 5.8 nsec (making experimental observation of this species very difficult).

Walch, Stephen P.↗

Theoretical study of the far-infrared A 3Sigma(-)g - X 3Pi(u) transition in Al2

The transition between the X 3Pi(u) ground state of Al2 and the first excited state A 3Sigma(-)g is investigated theoretically by means of CASSCF computations. The Einstein coefficients are determined on the basis of the potentials published by Bauschlicher et al. (1989) and listed in a table. The problems involved in devising an experiment to verify the theoretical results are briefly considered.

Bauschlicher, Charles W., Jr.↗

Theoretical characterization of the potential energy surface for H + O2 yields HO2(asterisk) yields HO + O. II - The potential for H atom exchange in HO2

The results of CASSCF multireference contracted CI calculations with large ANO basis sets are presented for the exchange region of the HO2 potential-energy surface. The saddle point for H atom exchange is about 13 kcal/mol below the energy of H + O2; therefore, this region of the surface should be accessible during H + O2 recombination and methathesis reactions.

Walch, Stephen P.↗

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