Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled-cluster methods”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

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

A Doubles Correction to Electronic Excited States from Configuration Interaction in the Space of Single Substitutions

A perturbative correction to the method of configuration interaction with single substitutions (CIS) is presented. This CIS(D) correction approximately introduces the effect of double substitutions which are absent in CIS excited states. CIS(D) is a second-order perturbation expansion of the coupled-cluster excited state method, restricted to single and double substitutions, in a series in which CIS is zeroth order, and the first-order correction vanishes. CIS (D) excitation energies are size consistent and the calculational complexity scales with the fifth power of molecular size, akin to second-order Moller-Plesset theory for the ground state. Calculations on singlet excited states of ethylene, formaldehyde, acetaldehyde, butadiene and benzene show that CIS (D) is a uniform improvement over CIS. CIS(D) appears to be a promising method for examining excited states of large molecules, where more accurate methods are not feasible.

Head-Gordon, Martin↗
NASA NTRS1990

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

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

Accurate Quartic Force Fields and Vibrational Frequencies for HCN and HNC

The quartic force fields of HCN and HNC are determined using atomic natural orbital one-particle basis sets of spdf/spd and spdfg/spdf quality in conjunction with the CCSD(T) electron correlation method (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations). The HCN force field is in good agreement with a recent experimentally derived force field and also with the force field recently computed by Wong and Bacskay. On the basis of the good agreement obtained for HCN, it is argued that the ab initio quartic force field for HNC is superior to a prior force field derived from experiment. The harmonic frequencies of HNC are predicted to be 3822 +/- 10,472 +/- 5, and 2051 +/-10/cm for omega1, omega2, and omega3, respectively; the experimentally derived values are above these values and fall outside the estimated uncertainties. Using the quartic force field, spectroscopic constants are predicted for HNC based on a vibrational second-order perturbation theory analysis. It is also asserted that the gas-phase fundamental v(sub 3) for HNC is slightly lower than the matrix isolation value. The range of validity of the quartic force fields is investigated by comparison of variational vibrational energies computed with the quartic force fields to those obtained from our recently reported global HCN/HNC potential energy surface and also to experimental data.

Lee, Timothy J.↗
NASA NTRS1993

Accurate Quartic Force Fields and Vibrational Frequencies for HCN and HNC

The quartic force fields of HCN and HNC are determined using atomic natural orbital one-particle basis sets of spdf/spd and spdfg/spdf quality in conjunction with the CCSD(T) electron correlation method (singles and doubles coupled-cluster theory plus a perturbation estimate of the effects of connected triple excitations). The HCN force field is in good agreement with a recent experimentally derived force field and also with the force field recently computed by Wong and Bacskay. On the basis of the good agreement obtained for HCN, it is argued that the ab initio quartic force field for HNC is superior to a prior force field derived from experiment. The harmonic frequencies of HNC are predicted to be 3822 +/- 10, 472 +/- 5, and 2051 +/- 10 cm(exp -1) for omega(sub 1), omega(sub 2), and omega(sub 3), respectively; the experimentally derived values are above these values and fall outside the estimated uncertainties. Using the quartic force field, spectroscopic constants are predicted for HNC based on a vibrational second-order perturbation theory analysis. It is also asserted that the gas-phase fundamental nu(sub 3) for HNC is slightly lower than the matrix isolation value. The range of validity of the quartic force fields is investigated by comparison of variational vibrational energies computed with the quartic force fields to those obtained from our recently reported global HCN/HNC potential energy surface and also to experimental data.

Lee, Timothy J.↗
NASA NTRS1989

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

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↗
NASA NTRS1991

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

The Molecular Structure of cis-FONO

The molecular structure of cis-FONO has been determined with the CCSD(T) correlation method using an spdf quality basis set. In agreement with previous coupled-cluster calculations but in disagreement with density functional theory, cis-FONO is found to exhibit normal bond distances. The quadratic and cubic force fields of cis-FONO have also been determined in order to evaluate the effect of vibrational averaging on the molecular geometry. Vibrational averaging is found to increase bond distances, as expected, but it does not affect the qualitative nature of the bonding. The CCSD(T)/spdf harmonic frequencies of cis-FONO support our previous assertion that a band observed at 1200 /cm is a combination band (upsilon(sub 3) + upsilon(sub 4)), and not a fundamental.

Lee, Timothy J.↗
NASA NTRS1993

Ab initio characterization of ClOOH - Implications for atmospheric chemistry

The equilibrium structure, dipole moment, harmonic vibrational frequencies, and infrared intensities of ClOOH are determined using the CCSD(T) (singles and doubles coupled-cluster theory plus a perturbational estimate of the effects of connected triple excitations) electronic structure method in conjunction with a TZ2P (triple xi plus double polarization) basis set. The heat of formation of CIOOH is determined (using two different isodesmic reactions) to be +1.5 +/- 1 kcal/mol at 0 K or +0.2 +/- 1 kcal/mol at 298.15 K. Using the computed heat of formation, we examined the stability of ClOOH with respect to the ClO + OH, ClOO + H, and HOO + Cl dissociation limits. Since ClOOH is found to be quite stable, it is argued that the chemistry of ClOOH should be included in any accurate modeling of the stratosphere.

Lee, Timothy J.↗
NASA NTRS1991

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

Triple and quadruple excitation contributions to the binding in Be clusters: Calibration calculations on Be3

The contribution of connected triple and quadruple excitations to the binding in Be3 is investigated by comparing various coupled-cluster (CC) and truncated configuration interaction (CI) treatments with multireference CI (MRCI) and full CI(FCI) calculations. The CC method with single and double excitations (CCSD) produces results that differ substantially from more elaborate treatments, but most extensions to CCSD that account approximately for connected triple excitations perform very well. In constrast, good agreement with FCI for Be2 can be achieved only with the highest level CC and MRCI methods.

Watts, John D.↗
NASA NTRS1990

Triple and quadruple excitation contributions to the binding in Be clusters - Calibration calculations on Be3

The contribution of connected triple and quadruple excitations to the binding in Be3 is investigated by comparing various coupled-cluster (CC) and truncated configuration-interaction (CI) treatments with multireference CI (MRCI) and full CI (FCI) calculations. The CC method with single and double excitations (CCSD) produces results that differ substantially from more elaborate treatments, but most extensions to CCSD that account approximately for connected triple excitations perform very well. In contrast, good agreement with CFI for Be2 can be achieved only with the highest level CC and MRCI methods.

Watts, John D.↗
NASA NTRS1994

A Coupled-Cluster Study of the Molecular Structure, Vibrational Spectrum, and Relative Energies of the XCN and XNC (X=F, Cl) Isomers

The XCN and XNC (X=F, Cl) isomers have been investigated using the CCSD and CCSD(T) methods in conjunction with a TZ2P basis set. Equilibrium geometries, dipole moments, harmonic frequencies, IR intensities and relative energies have been evaluated. The CCSD(T) geometries and vibrational frequencies for the XCN isomers are in good agreement with the available experimental data. The CCSD(T) results for FCN and FNC are in good agreement with the CEPA calculations of Botshwina et al., with the exception of the energy difference, which the CEPA method underestimates by about 1.2 kcal/mol. FCN and CICN are shown to be lower in energy than the FNC and ClNC isomers by 69.511.0 and 42.711.0 kcal/mol (0 K), respectively.

Lee, Timothy J.↗
NASA NTRS1993

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

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

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