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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

On The Choice of a Zeroth-Order Hamiltonian for Second-Order Perturbation Theory with A CASSCF Reference Function

A new approach to perturbation theory based on a CASSCF reference function has been developed. The key to the approach is the definition of the zeroth order Hamiltonian, H(sub 0), which includes the full CI Hamiltonian for the active space. In the the inactive and secondary spaces, operators may be chosen which reduce to the usual Moller-Plesset or Epstein-Nesbet forms in the limit of a null active space. These operators are diagonal in the orbital indices and permit the block-diagonalization of H(sub 0). The reference is an eigenfunction of Ho without N-particle projection. H(sub 0) automatically incorporates denominator shifts in the style of those appearing in recent open-shell perturbation theories. Comparative results are presented for a few test cases.

Dyall, Kenneth G.↗
NASA NTRS2001

Theoretical Studies of Possible Synthetic Routes for the High Energy Density Material Td N4: Excited Electronic States

Vertical electronic excitation energies for single states have been computed for the high energy density material (HEDM) Td N4 in order to assess possible synthetic routes that originate from excited electronic states of N2 molecules. Several ab initio theoretical approaches have been used, including complete active space self-consistent field (CASSCF), state averaged CASSCF (SA-CASSCF), singles configuration interaction (CIS), CIS with second-order and third-order correlation corrections [CIS(D)) and CIS(3)], and linear response singles and doubles coupled-cluster (LRCCSD), which is the highest level of theory employed. Standard double zeta polarized (DZP) and triple zeta double polarized (TZ2P) one-particle basis sets were used. The CASSCF calculations are found to overestimate the excitation energies, while the SA-CASSCF approach rectifies this error to some extent, but not completely. The accuracy of the CIS calculations varied depending on the particular state, while the CIS(D), CIS(3), and LRCCSD results are in generally good agreement. Based on the LRCCSD calculations, the lowest six excited singlet states are 9.35(l(sup)T1), 10.01(l(sup)T2), 10.04(1(sup)A2), 10.07(1(sup)E), 10.12(2(sup)T1), and 10.42(2(sup)T2) eV above the ground state, respectively. Comparison of these excited state energies with the energies of possible excited states of N2+N2 fragments, leads us to propose that the most likely synthetic route for Td N4 involving this mechanism arises from combination of two bound quintet states of N2.

Lee, Timothy J.↗
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 NTRS2007

Theoretical Studies of Routes to Synthesis of Tetrahedral N4

A paper [Chem. Phys. Lett. 345, 295 (2001)] describes theoretical studies of excited electronic states of nitrogen molecules, with a view toward utilizing those states in synthesizing tetrahedral N4, or Td N4 a metastable substance under consideration as a high-energy-density rocket fuel. Several ab initio theoretical approaches were followed in these studies, including complete active space self-consistent field (CASSCF), state-averaged CASSCF (SA-CASSCF), singles configuration interaction (CIS), CIS with second-order and third-order correlation corrections [CIS(D) and CIS(3)], and linear response singles and doubles coupled-cluster (LRCCSD). Standard double zeta polarized and triple zeta double polarized one-particle basis sets were used. The CASSCF calculations overestimated the excitation energies, while SACASSCF calculations partly corrected these overestimates. The accuracy of the CIS calculations varied, depending on the particular state, while the CIS(D), CIS(3), and LRCCSD results were in generally good agreement. The energies of the lowest six excited singlet states of Td N4 as calculated by the LRCCSD were compared with the energies of possible excited states of N2 + N2 fragments, leading to the conclusion that the most likely route for synthesis of Td N4 would involve a combination of two bound quintet states of N2.

Lee, Timothy J.↗
NASA NTRS1988

Full configuration-interaction study of the ionic-neutral curve crossing in LiF

Full configuration-interaction (FCI) calculations are used to assess the relative ability of methods for truncating the n-particle expansion in describing the ionic-neutral curve crossing between the two lowest Sigma(+) states of LiF. While the FCI calculations yield a smooth dipole moment function, MRCI calculations based on CASSCF orbitals optimized for the lowest state at all r values yield a discontinuous dipole moment function. However, when the orbitals are optimized using a state-averaged CASSCF procedure, with equal weights for the ionic and neutral solutions, both the CASSCF and MRCI dipole moment functions are smooth and in reasonable agreement with the FCI. No single-reference-based method is found to work satisfactorily. Potential curves for the lowest two Sigma(+) states are determined in both the adiabatic and diabatic representations using a large atomic natural orbit Gaussian basis set and a state-averaged CASSCF/MRCI treatment of electron correlation.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1995

The Ozonide Anion: A Theoretical Study

Complete active space self-consistent field (CASSCF) and CASSCF second order perturbation theory (CASPT2) methods have been used for the geometry optimization and calculation of harmonic and fundamental frequencies of the ozonide ion O3(-) and the ozonide lithium complex Li(+)O3(-). For O3(-) harmonic frequencies have also been obtained using the coupled-cluster method, CCSD(T). Infrared intensities are computed from dipole moment derivatives at the CASSCF level. The predicted equilibrium geometry for O3(-) is: Roo = 1.361 Angstroms and angle ooo = 115.4 degrees, and the fundamental frequencies are: nu(sub 1) = 989 per centimeter, nu(sub 2) 556 per centimeter, nu(sub 3) = 870 per centimeter (experimental values are: Roo =1.36 plus or minus 0.02 Angstroms, angle (ooo) = 111.8 plus or minus 2.0 degrees, nu(sub 1) = 975(50) per centimeter, nu(sub 2) = 550(50) per centimeter, nu(sub 3) = 880(50) per centimeter). Corresponding data for the lithium ozonide complex have also been obtained. The presented data contradict the previous interpretation of the IR and Raman spectrum obtained after deposition of ozone in N2, argon, or neon matrices with atomic beams of alkali metals. The presence of the lithium cation raises the asymmetric stretch frequency to about 940 per centimeter, which is contradictory to assumptions made in the assignments of the matrix spectra. Calculations made in a dielectric medium strongly suggest that the effect of the matrix on the IR spectrum is small for O3(-) itself. The dissociation and atomization energies of O3(-) are found to be in agreement with experiment.

Borowski, Poitr↗
NASA NTRS1986

On the electron affinity of the oxygen atom

The electron affinity (EA) of oxygen is computed to be 1.287 eV, using 2p electron full configuration-interaction (CI) wave functions expanded in a 6s5p3d2f Slater-type orbital basis. The best complete active space self-consistent field - multireference CI (CASSCF-MRCI) result including only 2p correlation is 1.263 eV. However, inclusion of 2s intrashell and 2s2p intershell correlation increases the computed EA to 1.290 at the CASSCF-MRCI level. At the full CI basis set limit, the 2s contribution to the electron affinity is estimated to be as large as 0.1 eV. This study clearly establishes the synergistic effect between the higher excitations and basis set completeness on the electron affinity when the 2s electrons are correlated.

Bauschlicher, C. W., Jr.↗
NASA NTRS1986

On the electron affinity of the oxygen atom

The electron affinity of oxygen is computed to be 1.287 eV, at the full CI level using a 6s5p3d 2f Slater-type orbital basis and correlating only the 2p electrons. The best CASSCF-MRCI result including only 2p correlation is 1.263 eV. However, inclusion of 2s intrashell and 2s2p intershell coorelation increases the computed EA to 1.290 eV at the CASSCF-MRCI level. At the full CI basis set limit, the 2s contribution to the electron affinity is estimated to be as large as 0.1 eV. The higher excitation contribution to the electron affinity is found to increase substantially with basis set completeness, especially when the 2s electrons are correlated. Relativistic effects are shown to make a small (less than 0.01 eV) change in the EA.

Bauschlicher, C. W., Jr.↗
NASA NTRS1987

Full CI benchmark calculations for several states of the same symmetry

Full CI (FCI) wave functions are used to compute energies for several electronic states of the same symmetry for SiH2, CH2, and CH2(+). It is found that CASSCF/multireference CI wave functions yield results very similar to FCI, irrespective of whether the CASSCF MOs are optimized independently for each state or using an average of the CASSCF energies for all desired states. The ionization potentials and excitation energies obtained from the FCI calculations should help calibrate methods (such as Green's function approaches, equations of motion and propagator methods, and cluster expansions) in which energy differences are computed directly.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1987

Theoretical study of the dipole moment function of OH(X 2Pi)

A complete theoretical electric dipole moment function mu(r) is computed for OH(X 2Pi) at the full configuration-interaction (FCI) level for DZP + diffuse Gaussian basis set. The energy and dipole moment are also determined with 2d functions on oxygen at seven r values. The following methods are compared to the FCI for both the oxygen 1d and 2d basis sets: single-reference singles-plus-doubles configuration interaction (SDCI), SDCI with Davidson's correction, the coupled pair functional approach (CPF) and a modified form of CPF, CASSCF, and CASSCF-MRSDCI, and the external contracted CI(CCI) method. The dipole moments are evaluated both as an expectation value and by finite-field methods. The results support the superiority of evaluating the dipole moment as an energy derivative rather than as an expectation value.

Langhoff, Stephen R.↗
NASA NTRS1985

Theoretical studies of transition metal dimers

The CASSCF approach was used to perform the MCSCF calculations for a number of transition metal dimers, including the Sc2, Ti2, Cr2, Cu2, TiV, Y2, Nb2, and Mo2 molecules; in addition, CASSCF/CI calculations were carried out for Sc2, Ti2, Cu2, and Y2. The CASSCF procedure is shown to provide a consistent set of calculations for these molecules, from which trends and a simple qualitative picture of the electronic structure may be derived. In particular, the calculations confirmed the ground states of the Sc2 and the TiV, and led to predictions for other molecules in this series. In addition to specific predictions, the study provides a simple qualitative picture of the bonding in these dimers.

Walch, Stephen P.↗
NASA NTRS1991

The computation of C-C and N-N bond dissociation energies for singly, doubly, and triply bonded systems

The bond dissociation energies (De) of C2H2, C2H4, C2H6, N2, N2H2, and N2H4 are studied at various levels of correlation treatment. The convergence of De with respect to the one-particle basis is studied at the single-reference modified coupled-pair-functional (MCPF)level. At all levels of correlation treatment, the errors in the bond dissociation energies increase with the degree of multiple bond character. The multireference configuration-interaction (MRCI) De values, corrected for an estimate of higher excitations, are in excellent agreement with those determined using the size-extensive averaged-coupled-pair-functional (ACPF) method. The full-valence complete-active-space self-consistent-field (CASSCF)/MRCI calculations are reproduced very well by MRCI calculations based on a CASSCF calculation that includes in the active space only those electrons involved in the C-C or N-N bonds.

Langhoff, Stephen R.↗
NASA NTRS1989

The computation of C-C and N-N bond dissociation energies for singly, doubly, and triply bonded systems

The bond dissociation energies (D sub e) of C2H2, C2H4, C2H6, N2, N2H2, and N2H4 are studied at various levels of correlation treatment. The convergence of D sub e with respect to the one particle basis is studied at the single reference modified coupled-pair functional (MCPF) level. At all levels of correlation treatment, the errors in the bond dissociation energies increase with the degree of multiple bond character. The multireference configuration interaction (MRCI) D sub e values, corrected for an estimate of higher excitations, are in excellent agreement with those determined using the size extensive averaged coupled pair functional (ACPF) method. It was found that the full valence complete active space self consistent field (CASSCF)/MRCI calculations are reproduced very well by MRCI calculations based on a CASSCF calculation that includes in the active space only those electrons involved in the C-C or N-N bonds. To achieve chemical accuracy (1 kcal/mole) for the D sub e values of the doubly bonded species C2H4 and N2H2 requires one particle basis sets including up through h angular momentum functions (l = 5) and a multireference treatment of electron correlation: still higher levels of calculation are required to achieve chemical accuracy for the triply bonded species C2H2 and N2.

Langhoff, Stephen R.↗
NASA NTRS1994

Theoretical Study of the Nitric Oxide epsilon and 1100 A Bands

Potential-energy curves, transition-moment functions, Einstein Alpha coefficients and radiative lifetimes associated with the Rydberg state of nitric oxide are calculated in large Gaussian basis sets employing internally contracted MRCI from state-averaged CASSCF wave functions. Computationally, a strong sensitivity of the results to the choice of the CASSCF active space is seen, and careful examination of this behavior gives insight into matters related to the treatment of weakly occupied orbitals, which is a consideration in any study Rydberg excited states. The epsilon system transition-moment function and radiative lifetimes are in disagreement with another calculation in the literature, but reasons for the discrepancies are examined, and the present results are determined to be more accurate.

Sheehy, Jeffrey A.↗
NASA NTRS1995

The MCH2+ Systems: Do ScCH2+ and TiCH2+ Have C(sub s) or C(sub 2v) Symmetry and A Comparison of the B3LYP Method to Other Approaches

ScCH2(+) and TiCH2(+) are found to have C(sub s) symmetry, while the remaining systems have the expected C(sub 2v) symmetry. The C(sub s) symmetry structure is favored for ScCH2(+) and TiCH2(+) because it allows donation from one of the CH bonds into an empty 3d orbital. This distortion stabilizes the systems by less than 2 kcal/mol at the B3LYP, CASSCF, and CCSD(T) levels of theory. The SCF and B3LYP approaches find the C2(sub v) structure to be a saddle point, while the CASSCF approach finds the C(sub 2v) structure to be a stationary point. For V-Cu there are no empty 3d orbitals of the correct symmetry, and therefore the distortion to C(sub s), is unfavorable. The B3LYP binding energies are in good agreement with experiment and our previous best estimates.

Ricca, Alessandra↗
NASA NTRS1983

CASSCF/CI calculations for first row transition metal hydrides - The TiH(4-phi), VH(5-delta), CrH(6-sigma-plus), MnH(7-sigma-plus), FeH(4,6-delta) and NiH(2-delta) states

Calculations are performed for the predicted ground states of TiH(4-phi), VH(5-delta), CrH(6-sigma-plus), MnH(7-sigma-plus), Fett(4,6-delta) and NiH(2-delta). For FeH both the 6-delta and 4-delta states are studied, since both are likely candidates for the ground state. The ground state symmetries are predicted based on a combination of atomic coupling arguments and coupling of 4s(2)3d(n) and 4s(1)3d(n+1) terms in the molecular system. Electron correlation is included by a CASSCF/CI (SD) treatment. The CASSCF includes near-degeneracy effects, while correlation of the 3d electrons in included at the CI level.

Walch, S. P.↗
NASA NTRS1983

On the nature of the bonding in Cu2 - An ab initio viewpoint

The bonding in Cu2 has been calculatively investigated. A CASSCF calculation is performed where the two 4s electrons are distributed in all possible ways among the 4s and 4p orbitals. The CAS results have a 3d population of 9.99. The corresponding orbitals between the CAS core orbitals at r(e) and the atomic SCF core orbitals show an overlap of 13.9992 out of a possible 14, with all of the change in sigma space. Thus, the CAS results show virtually no change in the d orbitals relative to the atom and can be considered as showing only a single bond. The reduction in 3d population, the decrease in r(e), and the increase in D(e) upon correlating the 3d shell indicates that the 3d electrons are slightly involved in the bonding, but there is no indication of the formation of a triple bond.

Bauschlicher, C. W., Jr.↗