Engineering Papers⌕ Search

Engineering topics

Lee, Timothy J.

Publications and source records attributed to Lee, Timothy J..

At least 109 records · Page 6

Ab initio multireference study of the BN molecule

The lowest 1Sigma(+) and 3Pi states of the BN molecule are studied using multireference configuration interaction (MRCI) and averaged coupled-pair functional (ACPF) methods and large atomic natural orbital (ANO) basis sets, as well as several coupled cluster methods. Our calculations strongly support a 3Pi ground state, but the a1Sigma(+) state lies only 381 +/- 100/cm higher. The a1Sigma(+) state wave function exhibits strong multireference character and, consequently, the predictions of the perturbationally-based single-reference CCSD(T) coupled cluster method are not as reliable in this case as the multireference results. The theoretical predictions for the spectroscopic constants of BN are in good agreement with experiment for the Chi3Pi state, but strongly suggest a misassignment of the fundamental vibrational frequency for the a1Sigma(+) state.

Martin, J. M. L.↗

An extensive ab initio study of the structures, vibrational spectra, quadratic force fields, and relative energetics of three isomers of Cl2O2

Quantum mechanical computational methods are employed for an ab initio investigation of: (1) the molecular properties of the lowest isomers of the ClO dimer; and (2) predicted molecular and thermochemical properties. Techniques employed include electron correlation and particularly singles and doubles coupled-cluster (CCSD) theory with or without perturbational estimates of the effects of connected triple excitations. The isomers ClOClO and ClClO2 are found to have higher energies than the ClOOCl isomer, and the theoretical vibrational frequencies of the isomers are well correlated with experimental data. Experimental values of the heat of formation for the isomers are also compared with calculations based on an isodesmic reaction with Cl2O, H2O, and HOOH.

Lee, Timothy J.↗

A global ab initio potential for HCN/HNC, exact vibrational energies, and comparison to experiment

An ab initio (i.e., from first principles) calculation of vibrational energies of HCN and HNC is reported. The vibrational calculations were done with a new potential derived from a fit to 1124 ab initio electronic energies which were calculated using the highly accurate CCSD(T) coupled-cluster method in conjunction with a large atomic natural orbital basis set. The properties of this potential are presented, and the vibrational calculations are compared to experiment for 54 vibrational transitions, 39 of which are for zero total angular momentum, J = 0, and 15 of which are for J = 1. The level of agreement with experiment is unprecedented for a triatomic with two nonhydrogen atoms, and demonstrates the capability of the latest computational methods to give reliable predictions on a strongly bound triatomic molecule at very high levels of vibrational excitation.

Bentley, Joseph A.↗

FONO - A difficult case for theory and experiment

High levels of ab initio theory are used to investigate the equilibrium structures, vibrational spectra, and relative energetics of FNO2, cis-FONO, and trans-FONO isomers. FNO2 is determined to be 36.9 +/- 2.5 kcal/mol (0 K) more stable than cis-FONO, which is more stable than trans-FONO by 2.5 +/- 1.0 kcal/mol (0 K). The molecular structure of cis-FONO is shown to have typical F-O and central O-N single bond distances, indicating that cis-FONO has a structure analogous to HONO. The computed vibrational spectrum of cis-FONO is shown to be consistent with IR matrix isolation experiments, including isotopic shifts. The experimentally deduced structure of FONO is shown to be incorrect because one of the vibrational bands included in the normal coordinate analysis is either a combination band or an overtone.

Lee, Timothy J.↗

The vibrational frequencies of difluoroethyne

Ab initio coupled-cluster calculations with single and double excitations and with a perturbational treatment of connected triple excitations are reported for difluoroethyne using large basis sets. The results for the transbending mode nu-4 are extremely sensitive to electron correlation and basis set effects. The best theoretical and experimental estimates for the fundamental vibrational frequencies are in excellent agreement.

Breidung, Juergen↗

Comparison of the Brueckner and coupled-cluster approaches to electron correlation

The coupled-cluster and Brueckner electron correlation methods are used to determine the equilibrium structures, dipole moments, harmonic vibrational frequencies, and IR intensities of NH3, FON, Be3, BeC2, and BeO2. The singles and doubles coupled-cluster and Brueckner doubles methods are employed, and the corresponding methods that include a perturbational estimate of connected triple excitations are also investigated. The T sub 1 diagnostic is found to provide a good indication of the magnitude of the difference between the results obtained with the coupled-cluster and Brueckner methods. For NH3, the T sub 1 diagnostic is small and so the differences between results obtained from coupled-cluster and Brueckner theories are quite small. It is found for all of the molecules under consideration that inclusion of the contribution from connected triple excitations is more important than the differences between the Brueckner and coupled-cluster correlation methods.

Lee, Timothy J.↗

The structures, binding energies and vibrational frequencies of Ca3 and Ca4: An application of the CCSD(T) method

The Ca3 and Ca4 metallic clusters have been investigated using state-of-the-art ab initio quantum mechanical methods. Large atomic natural orbital basis sets have been used in conjunction with the singles and doubles coupled-cluster (CCSD) method, a coupled-cluster method that includes a perturbational estimate of connected triple excitations, denoted CCSD(T), and the multireference configuration interaction (MRCI) method. The equilibrium geometries, binding energies and harmonic vibrational frequencies have been determined with each of the methods so that the accuracy of the coupled-cluster methods may be assessed. Since the CCSD(T) method reproduces the MRCI results very well, cubic and quartic force fields of Ca3 and Ca4 have been determined using this approach and used to evaluate the fundamental vibrational frequencies. The infrared intensities of both the e' mode of Ca3 and the t2 mode of Ca4 are found to be small. The results obtained in this study are compared and contrasted with those from our earlier studies on small Be and Mg clusters.

Lee, Timothy J.↗

Connected triple excitations in coupled-cluster calculations of hyperpolarizabilities: Neon

We have calculated the second hyperpolarizability gamma of neon using the CCSD(T) method. The accuracy of the CCSD(T) approach has been established by explicit comparison with the single, double and triple excitation coupled-cluster (CCSDT) method using extended basis sets that are known to be adequate for the description of gamma. Our best estimate for gamma(sub 0) of 110 +/- 3 a.u. is in good agreement with other recent theoretical values and with Shelton's recent experimental estimate of 108 +/- 2 a.u. Comparison of the MP2 and CCSD(T) hyperpolarizability values indicates that MP2 gives a very good description of the electron correlation contribution to gamma(sub 0). We have combined MP2 frequency-dependent corrections with the CCSD(T) gamma(sub 0) to yield values of gamma(-2 omega;omega,omega,0) and gamma(exp K)(-omega;omega,0,0).

Rice, Julia E.↗

The calculation of the vibrational frequencies of CuCO+, NiCO and CuCH3

The vibrational frequencies of CuCO+, NiCO, and CuCH3 are computed at both the self-consistent-field (SCF) and coupled-cluster singles and doubles, with a perturbational estimate of the triples (CCSD(T)), levels of theory. Overall the SCF frequencies are in qualitative agreement with the CCSD(T) results. The calculations show that the frequencies for the electrostatically bound CuCO+ are quite different from the datively bound NiCO. The CCSD(T) calculations support a previous suggestion that the C-H stretches in CuCH3 have been misassigned.

Sodupe, M.↗

Proton affinity of methyl nitrate - Less than proton affinity of nitric acid

Several state-of-the-art ab initio quantum mechanical methods were used to investigate the equilibrium structure, dipole moments, harmonic vibrational frequencies, and IR intensities of methyl nitrate, methanol, and several structures of protonated methyl nitrate, using the same theoretical methods as in an earlier study (Lee and Rice, 1992) of nitric acid. The ab initio results for methyl nitrate and methanol were found to be in good agreement with available experimental data. The proton affinity (PA) of methyl nitrate was calculated to be 176.9 +/-5 kcal/mol, in excellent agreement with the experimental value 176 kcal/mol obtained by Attina et al. (1987) and less than the PA value of nitric acid. An explanation of the discrepancy of the present results with those of an earlier study on protonated nitric acid is proposed.

Lee, Timothy J.↗

Ab initio study of the molecular structure and vibrational spectrum of nitric acid and its protonated forms

The equilibrium structures, harmonic vibrational frequencies, IR intensities, and relative energetics of HNO3 and its protonated form H2NO3+ were investigated using double-zeta plus polarization and triple-zeta plus polarization basis sets in conjunction with high-level ab initio methods. The latter include second-order Moller-Plesset perturbation theory, the single and double excitation coupled cluster (CCSD) methods, a perturbational estimate of the effects of connected triple excitations (CCSD(T)), and the self-consistent field. To determine accurate energy differences CCSD(T) energies were computed using large atomic natural orbital basis sets. Four different isomers of H2NO3+ were considered. The lowest energy form of protonated nitric acid was found to correspond to a complex between H2O and NO2+, which is consistent with earlier theoretical and experimental studies.

Lee, Timothy J.↗

The energy separation between the classical and nonclassical isomers of protonated acetylene - An extensive study in one- and n-particle space saturation

The energy separation between the classical and nonclassical forms of protonated acetylene has been reinvestigated in light of the recent experimentally deduced lower bound to this value of 6.0 kcal/mol. The objective of the present study is to use state-of-the-art ab initio quantum mechanical methods to establish this energy difference to within chemical accuracy (i.e., about 1 kcal/mol). The one-particle basis sets include up to g-type functions and the electron correlation methods include single and double excitation coupled-cluster (CCSD), the CCSD(T) extension, multireference configuration interaction, and the averaged coupled-pair functional methods. A correction for zero-point vibrational energies has also been included, yielding a best estimate for the energy difference between the classical and nonclassical forms of 3.7 + or - 1.3 kcal/mol.

Lindh, Roland↗

An efficient formulation and implementation of the analytic energy gradient method to the single and double excitation coupled-cluster wave function - Application to Cl2O2

The analytic energy gradient for the single and double excitation coupled-cluster (CCSD) wave function has been reformulated and implemented in a new set of programs. The reformulated set of gradient equations have a smaller computational cost than any previously published. The iterative solution of the linear equations and the construction of the effective density matrices are fully vectorized, being based on matrix multiplications. The new method has been used to investigate the Cl2O2 molecule, which has recently been postulated as an important intermediate in the destruction of ozone in the stratosphere. In addition to reporting computational timings, the CCSD equilibrium geometries, harmonic vibrational frequencies, infrared intensities, and relative energetics of three isomers of Cl2O2 are presented.

Rendell, Alistair P.↗

Analytic gradients for coupled-cluster energies that include noniterative connected triple excitations - Application to cis- and trans-HONO

An efficient formulation of the analytic energy gradient for the single and double excitation coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations is presented. The formulation has a small computational cost, and the algebraic manipulations may be applied generally to the analytic gradient of Moller-Plesset perturbation theory energies. The new formulation has been implemented in an efficient set of programs that utilize highly vectorized algorithms and has been used to investigate the equilibrium structures, harmonic vibrational frequencies, IR intensities, and energy separation of cis- and trans-HONO.

Lee, Timothy J.↗

A parallel vectorized implementation of triple excitations in CCSD(T) - Application to the binding energies of the AlH3, AlH2F, AlHF2 and AlF3 dimers

An efficient method for various noniterative estimates of connected triple excitations in coupled-cluster theory is outlined and related to a similar expression occurring in Moller-Plesset perturbation theory. The method is highly vectorized and capable of utilizing multiple processors on a shared-memory machine, leading to computational rates in excess of one billion floating-point operations per second on four processors of a CRAY Y-MP. Using the new procedure, the binding energies of the D(2h) diborane-type dimers of AlH3, AlH2F, AlHF2, and AlF3 have been determined to be 32, 40, 20, and 47 kcal/mol, respectively. For Al2F6, the correlation procedure includes 232 molecular orbitals and over 1.5 x 10 to the 6th single and double coupled-cluster amplitudes, effectively accounting for over 2 x 10 to the 9th connected triple excitations.

Rendell, Alistair P.↗

The determination of accurate dipole polarizabilities alpha and gamma for the noble gases

Accurate static dipole polarizabilities alpha and gamma of the noble gases He through Xe were determined using wave functions of similar quality for each system. Good agreement with experimental data for the static polarizability gamma was obtained for Ne and Xe, but not for Ar and Kr. Calculations suggest that the experimental values for these latter ions are too low.

Rice, Julia E.↗

The structure and energetics of the HCN-HNC transition state

The optimum geometries and quadratic force constants of HCN, HNC and the transition state connecting them have been determined at the single- and double-excitation coupled-cluster (CCSD) and CCSD(T) levels of theory. Energy differences were evaluated using the CCSD and CCSD(T) methods in conjunction with large atomic natural orbital basis sets containing g-type basis functions on the heavy atoms and f-type functions on hydrogen. The most reliable structure obtained for the transition state has bond distances of 1.194, 1.188, and 1.389 A for r(CN), r(CH), and r(NH), respectively. Including a correction for zero-point vibrational energies, the transition state is predicted to be 44.6 + or - 1.0 kcal/mol above the HCN isomer, while HNC is predicted to be 14.4 + or - 1.0 kcal/mol above HCN. The latter value is in excellent agreement with the most recent experimental determination (14.8 + or - 2.0 kcal/mol).

Lee, Timothy J.↗

Theoretical characterization of tetrahedral N4

Results are reported from ab initio theoretical calculations of the structure and energetics of tetrahedral N4. The theoretical basis and implementation of the computations are reviewed, and the results are presented in tables. It is found that a tetrahedral N4 molecule has energy 186 kcal/mole greater than that of two N2 molecules, with a 61-kcal/mole barrier between the two species. Brief consideration is given to the possible use of N4 (should a process for synthesizing it be developed) as a high-energy-density fuel for supersonic transport aircraft.

Lee, Timothy J.↗