Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled cluster”

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

Vibrational frequencies for Be3 and Be4

The harmonic frequencies of Be3 and Be4, evaluated by multireference configuration-interaction (MRCI) and singles and doubles coupled-cluster (CCSD) methods, are reported, as are IR intensities evaluated by using the double harmonic approximation. A comparison of the results shows that the CCSD method augmented with a correction for connected triple excitations, CCSD(T), can be a viable alternative to MRCI for determining the ground-state properties of small Be clusters (other than the dimer). Results of normal mode analyses indicate that the anharmonic vibrational frequencies of Be3 and Be4 are reasonably consistent at the CCSD, CCSD(T), or MRCI levels of treatment.

Rendell, Alistair P.↗

The vibrational frequencies of ozone

The singles and doubles coupled-cluster (CCSD) method is used, together with the recently proposed CCSD(T) and CCSDT-1 methods, to determine the equilibrium structure and harmonic frequencies of ozone at various levels, on the basis of large one-particle basis sets. It is demonstrated that f-type functions are required in order to obtain quantitative accuracy for the equilibrium structure and harmonic frequencies. The CCSD(T) level of theory, in conjunction with a (4s 3p 2d 1f) atomic natural orbital set, yields an equilibrium geometry and harmonic frequencies that are in excellent agreement with experimental results; an effort has accordingly been made to determine the cubic and quartic force constants, and evaluate the anharmonic corrections to the fundamental vibrational frequencies.

Lee, Timothy J.↗

Theoretical study of polarizabilities and hyperpolarizabilities of Ne, HF, F(-), and OH(-)

Theoretical calculations are presented for the polarizability parameters through the fourth power of an applied static electric field for Ne and HF and through the second power in the field for the negative ions F(-) and OH(-). Large Gaussian basis sets are employed and electron correlation is incorporated using the modified coupled pair functional (MCPF) method and for Ne and HF using coupled-cluster single- and double-excitation (CCSD) methods with two approaches of including the effect of triple excitations. The MCPF method is found to perform well in relation to the more rigorous CCSD method.

Chong, Delano P.↗

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

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

The static dipole polarizabilities alpha and gamma for the noble gases helium through xenon were determined using large flexible one-particle basis sets in conjunction with high-level treatments of electron correlation. The electron correlation methods include single and double excitation coupled-cluster theory (CCSD), an extension of CCSD that includes a perturbational estimate of connected triple excitations, CCSD(T), and second order perturbation theory (MP2). The computed alpha and gamma values are estimated to be accurate to within a few percent. Agreement with experimental data for the static hyperpolarizability gamma is good for neon and xenon, but for argon and krypton the differences are larger than the combined theoretical and experimental uncertainties. Based on our calculations, we suggest that the experimental value of gamma for argon is too low; adjusting this value would bring the experimental value of gamma for krypton into better agreement with our computed result. The MP2 values for the polarizabilities of neon, argon, krypton and zenon are in reasonabe agreement with the CCSD and CCSD(T) values, suggesting that this less expensive method may be useful in studies of polarizabilities for larger systems.

Rice, Julia E.↗

Computational studies of metal-metal and metal-ligand interactions

The geometric structure of Cr(CO)6 is optimized at the modified coupled-pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. A detailed comparison of the properties of free CO is therefore given, at both the MCPF and CCSD/CCSD(T) levels of treatment, using a variety of basis sets. With very large one-particle basis sets, the SSCD(T) method gives excellent results for the bond distance, dipole moment and harmonic frequency of free CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used here and in a previous study. Calculations using larger basis sets reduced the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. The remaining discrepancy between the experimental and theoretical total binding energy of Cr(CO)6 is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive se (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.

Barnes, Leslie A.↗

Interconversion of diborane(4) isomers

Highly correlated electronic structure computations using many-body perturbation theory and coupled-cluster gradient techniques are used to study the reaction pathway that links the two forms (C2u and D2d) of diborane(4). The results obtained indicate that a low-energy pathway exists for interconversion of the two low-lying isomers of diborane(4). The proposed mechanism consists of a single concerted but nonsynchronous rotation of the BH2 groups. The pathway first follows an idealized reaction coordinate which preserves C2 symmetry, but then bifurcates at a branch point, leading to two equivalent transition states which lack nontrivial elements of symmetry.

Stanton, John F.↗

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

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↗

The structure and energetics of Cr(CO)6 and Cr(CO)5

The geometric structure of Cr(CO)6 is optimized at the modified coupled pair functional (MCPF), single and double excitation coupled-cluster (CCSD) and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), and the force constants for the totally symmetric representation are determined. The geometry of Cr(CO)5 is partially optimized at the MCPF, CCSD, and CCSD(T) levels of theory. Comparison with experimental data shows that the CCSD(T) method gives the best results for the structures and force constants, and that remaining errors are probably due to deficiencies in the one-particle basis sets used for CO. The total binding energies of Cr(CO)6 and Cr(CO)5 are also determined at the MCPF, CCSD, and CCSD(T) levels of theory. The CCSD(T) method gives a much larger total binding energy than either the MCPF or CCSD methods. An analysis of the basis set superposition error (BSSE) at the MCPF level of treatment points out limitations in the one-particle basis used. Calculations using larger basis sets reduce the BSSE, but the total binding energy of Cr(CO)6 is still significantly smaller than the experimental value, although the first CO bond dissociation energy of Cr(CO)6 is well described. An investigation of 3s3p correlation reveals only a small effect. In the largest basis set, the total CO binding energy of Cr(CO)6 is estimated to be 140 kcal/mol at the CCSD(T) level of theory, or about 86 percent of the experimental value. The remaining discrepancy between the experimental and theoretical value is probably due to limitations in the one-particle basis, rather than limitations in the correlation treatment. In particular an additional d function and an f function on each C and O are needed to obtain quantitative results. This is underscored by the fact that even using a very large primitive set (1042 primitive functions contracted to 300 basis functions), the superposition error for the total binding energy of Cr(CO)6 is 22 kcal/mol at the MCPF level of treatment.

Barnes, Leslie A.↗

On the performance of large Gaussian basis sets for the computation of total atomization energies

The total atomization energies of a number of molecules have been computed using an augmented coupled-cluster method and (5s4p3d2f1g) and 4s3p2d1f) atomic natural orbital (ANO) basis sets, as well as the correlation consistent valence triple zeta plus polarization (cc-pVTZ) correlation consistent valence quadrupole zeta plus polarization (cc-pVQZ) basis sets. The performance of ANO and correlation consistent basis sets is comparable throughout, although the latter can result in significant CPU time savings. Whereas the inclusion of g functions has significant effects on the computed Sigma D(e) values, chemical accuracy is still not reached for molecules involving multiple bonds. A Gaussian-1 (G) type correction lowers the error, but not much beyond the accuracy of the G1 model itself. Using separate corrections for sigma bonds, pi bonds, and valence pairs brings down the mean absolute error to less than 1 kcal/mol for the spdf basis sets, and about 0.5 kcal/mol for the spdfg basis sets. Some conclusions on the success of the Gaussian-1 and Gaussian-2 models are drawn.

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

An accurate ab initio quartic force field for ammonia

The quartic force field of ammonia is computed using basis sets of spdf/spd and spdfg/spdf quality and an augmented coupled cluster method. After correcting for Fermi resonance, the computed fundamentals and nu 4 overtones agree on average to better than 3/cm with the experimental ones except for nu 2. The discrepancy for nu 2 is principally due to higher-order anharmonicity effects. The computed omega 1, omega 3, and omega 4 confirm the recent experimental determination by Lehmann and Coy (1988) but are associated with smaller error bars. The discrepancy between the computed and experimental omega 2 is far outside the expected error range, which is also attributed to higher-order anharmonicity effects not accounted for in the experimental determination. Spectroscopic constants are predicted for a number of symmetric and asymmetric top isotopomers of NH3.

Martin, J. M. L.↗

Bond length, dipole moment, and harmonic frequency of CO

A detailed comparison of some properties of CO is given, at the modified coupled-pair functional, single and double excitation coupled-cluster (CCSD), and CCSD(T) levels of theory (including a perturbational estimate for connected triple excitations), using a variety of basis sets. With very large one-particle basis sets, the CCSD(T) method gives excellent results for the bond distance, dipole moment, and harmonic frequency of CO. In a (6s 5p 4d 3f 2g 1h) + (1s 1p 1d) basis set, the bond distance is about 0.005a0 too large, the dipole moment about 0.005 a.u. too small, and the frequency about 6/cm too small, when compared with experimental results.

Barnes, Leslie A.↗

Structure and energetics of Cr(CO)6 and Cr(CO)5

The geometric structures and energetics of Cr(CO)6 and Cr(CO)5 are determined at the modified coupled-pair functional, single and double excitation coupled-cluster (CCSD), and CCSD(T) levels of theory. For Cr(CO)6, the structure and force constants for the totally symmetric representation are in good agreement with experimental data once basis set constants are taken into account. In the largest basis set at the CCSD(T) level of theory, the total binding energy of CR(CO)6 is estimated at around 140 kcal/mol, or about 86 percent of the experimental value. In contrast, the first bond energy of Cr(CO)6 is very well described at the CCSD(T) level of theory, with the best estimated value of 38 kcal/mol being within the experimental uncertainty.

Barnes, Leslie A.↗

Accurate ab initio quartic force fields for borane and BeH2

The quartic force fields of BH3 and BeH2 have been computed ab initio using an augmented coupled cluster (CCSD(T)) method and basis sets of spdf and spdfg quality. For BH3, the computed spectroscopic constants are in very good agreement with recent experimental data, and definitively confirm misassignments in some older work, in agreement with recent ab initio studies. Using the computed spectroscopic constants, the rovibrational partition function for both molecules has been constructed using a modified direct numerical summation algorithm, and JANAF-style thermochemical tables are presented.

Martin, J. M. L.↗

An open-shell restricted Hartree-Fock perturbation theory based on symmetric spin orbitals

A new open-shell perturbation theory is formulated in terms of symmetric spin orbitals. Only one set of spatial orbitals is required, thereby reducing the number of independent coefficients in the perturbed wavefunctions. For second order, the computational cost is shown to be similar to a closed-shell calculation. This formalism is therefore more efficient than the recently developed RMP, ROMP or RMP-MBPT theories. The perturbation theory described herein was designed to have a close correspondence with our recently proposed coupled-cluster theory based on symmetric spin orbitals. The first-order wavefunction contains contributions from only doubly excited determinants. Equilibrium structures and vibrational frequencies determined from second-order perturbation theory are presented for OH, NH, CH, 02, NH2 and CH2.

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