Engineering topics

Taylor, Peter R.

Publications and source records attributed to Taylor, Peter R..

87 records · Page 5

NASA NTRS1987

Accurate ab initio calculations for the ground states of N2, O2 and F2

Spectroscopic constants and dissociation energies for the ground states of N2, O2, and F2 determined at the CAS SCF MRCI correlation level are in excellent agreement with experiment when very large primitive valence and polarization one-particle Gaussian basis sets are employed. The dissociation energy (De) for N2 is larger than experiment unless the 2s electrons are correlated. The basis set requirement for an accurate determination of De is found to increase with the degree of multiple-bond character in the molecule.

Langhoff, Stephen R.↗

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

A full CI treatment of the 1A1, 1B1, and 3B1 states of SiH2

Full CI calculations are presented for the 1A1, 3B1, and 1B1 states of SiH2 at their respective equilibrium geometries and at geometries with the SiH bonds stretched. These results are compared with those obtained from single-reference and multireference CI calculations. The computed Te values agree well with the full CI results, provided that the effects of higher-than-double excitations are accounted for either by the Davidson correction or by a multireference approach. When the SiH bonds are stretched, the single-reference methods are not sufficiently flexible, and only CASSCF/MRCI achieves chemical accuracy (i.e., agrees with the full CI to 1 kcal/mol). Overall, the accuracy of the various approximate methods is very similar to that found for H2O, NH2, and CH2.

Bauschlicher, Charles W., Jr.↗

NASA NTRS1987

Full CI studies of the collinear transition state for the reaction F + H2 yields HF + H

Full CI calculations on the collinear transition state for the reaction F + H2 yields HF + H are reported. The full CI results are compared with those obtained from single-reference and multireference CI calculations and from single-reference CPF calculations. In general, only those methods which attempt to account for the effects of higher excitations, such as CPF or CI plus the Davidson correction, yield a transition-state location and barrier height in good agreement with the full CI. In an extended basis, the effect of higher excitations is estimated to lower the barrier by as much as 1.5 kcal/mol; such an effect would essentially eliminate the present discrepancy between theory and experiment.

Bauschlicher, Charles W., Jr.↗

NASA NTRS1987

Strategies for obtaining the maximum performance from current supercomputers

The exploitation of special features of the current generation of supercomputers is discussed with particular reference to computational quantum chemistry. Modification of algorithms in the light of the very large memory available on the CRAY 2 is described, and multi-processing of programs is considered for both the CRAY XMP and CRAY 2. The overhead in multi-processing can be reduced to a small percentage even for steps that perform extensive IO.

Taylor, Peter R.↗

NASA NTRS1987

A simple interpretation of the Fe2(-) photoelectron spectrum

The photoelectron spectrum of Fe2(-) can be simply interpreted in terms of electron detachment from the mildly bonding 4s sigma sup * sub u orbital of a (4s sigma sub s)(sup 2) (4s sigma sup * sub u)(sup 2) (3d)(sup 13) anion. This interpretation implies a (4s sigma sub g)(sup 2) (4s sigma sup * sub u)(sup1) (3d)(sup13) configuration for the ground state of Fe2, correlating with one ground state (4s2 3d6) and one excited state (4s1 3d7) Fe atom. A comparison of the bond length and vibrational frequency of Fe2 to values for transition metal dimers containing single 4s-4s bonds is suggestive of 3d-3d bonding in this molecule. The results of preliminary full-valence configuration interaction calculations provide stron support for the proposed Fe2 and Fe2(-) configurations.

Leopold, Doreen G.↗

NASA NTRS1987

Theoretical studies of the potential surface for the F - H2 greater than HF + H reaction

The F + H2 yields HF + H potential energy hypersurface was studied in the saddle point and entrance channel regions. Using a large (5s 5p 3d 2f 1g/4s 3p 2d) atomic natural orbital basis set, a classical barrier height of 1.86 kcal/mole was obtained at the CASSCF/multireference CI level (MRCI) after correcting for basis set superposition error and including a Davidson correction (+Q) for higher excitations. Based upon an analysis of the computed results, the true classical barrier is estimated to be about 1.4 kcal/mole. The location of the bottleneck on the lowest vibrationally adiabatic potential curve was also computed and the translational energy threshold determined from a one-dimensional tunneling calculation. Using the difference between the calculated and experimental threshold to adjust the classical barrier height on the computed surface yields a classical barrier in the range of 1.0 to 1.5 kcal/mole. Combining the results of the direct estimates of the classical barrier height with the empirical values obtained from the approximation calculations of the dynamical threshold, it is predicted that the true classical barrier height is 1.4 + or - 0.4 kcal/mole. Arguments are presented in favor of including the relatively large +Q correction obtained when nine electrons are correlated at the CASSCF/MRCI level.

Bauschlicher, Charles W., Jr.↗

NASA NTRS1987

The 2D Rydberg series in Al I

High quality ab initio electonic structure calculations were performed on the 2D Rydberg series in Al I. The configuration 3s3p2(2D) is shown to contribute substantially to the lowest four 2D Rydberg states. The same configuration also contributes substantially to a 2D state embedded in the ionization continuum. Computed oscillator strengths for the first six members of the 2D Rydberg transitions are given: these should be of substantially high accuracy than currently available values.

Taylor, Peter R.↗

NASA NTRS1987

On the bond distance in methane

The equilibrium bond distance in methane was optimized using coupled-pair functional and contracted CI wave functions, and a Gaussian basis that includes g-type functions on carbon and d-type functions on hydrogen. The resulting bond distance, when corrected for core-valence correlation effects, agrees with the experimental value of 2.052 a(0) to within the experimental uncertainty of 0.002 a(0). The main source of error in the best previous studies, which showed discrepancies with experiment of 0.007 a(0) is shown to be basis set incompleteness. In particular, it is important that the basis set be close to saturation, at least for the lower angular quantum numbers.

Bowen-Jenkins, Philippa↗

NASA NTRS1987

Theoretical approaches to metal chemistry

Theoretical calculations on metals ranging from very accurate ab initio studies of diatomic and triatomic systems to model studies of larger clusters are presented. Recent improvements in the representation of one-particle and n-particle spaces have made possible the prediction that Al2 has a 3Pi(u) ground state, even though the 3Sigma(-)g state lies within 200/cm. Results suggest that cluster geometry varies dramatically with cluster size, and that rather large clusters are required before the bulk structure becomes optimal. Al cluster studies show that three-body terms are needed for quantitative agreement with ab initio calculations.

Bauschlicher, Charles W.↗

NASA NTRS1987

Theoretical determination of the ground state of N2(2+)

The low-lying 1Sigma(g)(+), 3Pi(u), and 3Sigma(g)(-) states of N2(2+) in the immediate region of the minima in their respective potential curves are investigated. Spectroscopic constants for the 1Sigma(u)(+) state are computed for comparison with experiment. Approaches to computing the entire potential curves of many states N2(2+) are calibrated.

Taylor, Peter R.↗

NASA NTRS1987

Theoretical study of the dissociation energy and the red and violet band systems of CN

The dissociation energy (D sub O) of CN is determined to be 7.65 + or - 0.06 eV. This corresponds to delta H sub f (CN) = 105.3 + or - 1.5 kcal/mole, in excellent agreement with Engleman and Rouse (1975), but considerably larger than the recent value deduced from shock-tube studies by Colket (1984). The result is obtained not only from extensive ab initio MRCI calculations using a very large Gaussian basis set, but also from extrapolation of the directly computed value by comparison of computed and experimental results fo NO, C2, and N2. As an additional calibration of the methods, the D sub O value for CN was computed from the corresponding value for CN(-) using the experimental electron affinity data. The lifetime of the nu prime = 0 level of the violet (B 2 sigma + yields X 2 sigma +) system was computed to be 62.4 ns, in good agreement with both experiment and previous calculations. Lifetimes for the red (A 2 pi yields X 2 sigma +) system decrease with increasing nu prime, which is consistent both with the recent experiment and calculations. While the computed lifetimes are significantly longer that those obtained from the experiment, they are shorter than those deduced from an analysis of the solar spectrum. However the D sub O and f (sub OO) are consistent with Lambert's model for the solar spectrum.

Bauschlicher, Charles W., Jr.↗

NASA NTRS1986

A full CI treatment of the 1A1-3B1 separation in methylene

The accuracy of recent theoretical computations of the total energies and the adiabatic separation of the 1A1 and 3B1 states of CH2 is investigated on the basis of complete CI calculations using the double-zeta basis sets of Dunning (1970). The results are presented in a table and characterized in detail. The errors in the separation values are found to range from less than 0.01 kcal/mol for a CASSCF/MRSDCI calculation, to 0.38 kcal/mol for a Davidson-corrected SCF/SDCI calculation, to 14.17 kcal/mol for an uncorrected SCF calculation.

Bauschlicher, Charles W., Jr.↗

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

Accurate ab initio calculations which demonstrate a 3 Pi u ground state for Al2

The spectroscopic parameters and separations between the three low-lying X 3 Pi u, A 3 Sigma g -, and a 1 Sigma g + states of Al2 are studied as a function of both the one-particle and n-particle basis set. Approximate correlation treatments are calibrated against full Cl calculations correlating the six valence electrons in a double-zeta plus two d-function basis set. Since the CASSCF/MRCI 3 Pi u to 3 Sigma g - separation is in excellent agreement wtih the FCI value, the MRCI calculations were carried out in an extended (20s13p6d4f)/(6s5p3d2f) gaussian basis. Including a small correction for relativistic effects, the best estimate is that 3 Sigma g - state lies 174/cm above the 3 Pi u ground state. The 1 Sigma g + state lies at least 2000/cm higher in energy. At the CPF level, inclusion of 2s and 2p correlation has little effect on D sub e, reduces T sub e by only 26/cm, and shortens the bond lengths by about 0.02 a sub o. Further strong support for a 3 Pi u ground state comes from the experimental absorption spectra, since both observed transitions can be convincingly assigned as 3 Pi u yields 3 Pi g. The (2) 3 Pi g state is observed to be sensitive to the level of correlation treatment, and to have its minimum shifted to shorter rho values, such that the strongest experimental absorption peak probably corresponds to the 0 yields 2 transition.

Bauschlicher, Charles W., Jr.↗