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 37 records · Page 2

Coupled-Cluster Theory Employing Approximate Integrals: An Approach to Avoid the Input/Output and Storage Bottlenecks

By representing orbital products in an expansion basis, certain classes of two-electron integrals are approximated for use in CCSD(T) calculations (singles and doubles coupled-cluster plus a perturbational estimate of the effects of connected triple excitations). This leads to a very large reduction in disk storage and input/output requirements, with usually only a modest increase in computational effort. The new procedure will allow very large CCSD(T) calculations to be undertaken, limited only by available processor time. Using the molecular basis as the expansion basis, explicit numerical comparisons of equilibrium geometries, harmonic frequencies, and energy differences indicate that the error due to the use of approximate integrals is less than the error associated with truncation of the molecular basis set.

Rendell, Alistair P.↗

A Coupled-Cluster Study of the Molecular Structure, Vibrational Spectrum, and Heats of Formation of XONO2 (X=H, F, Cl)

The equilibrium structures, harmonic vibrational frequencies, dipole moments, and IR intensities of nitric acid, fluorine nitrate, and chlorine nitrate have been investigated by using the singles and doubles coupled cluster method that also includes a perturbational estimate of the effects of connected triple excitations, CCSD(T). A standard triple-zeta double-polarized basis set was utilized. The equilibrium geometries and vibrational spectra of HONO2 and ClONO2 are shown to be in excellent agreement with the available experimental data. The ab initio vibrational spectrum of FONO2 is also shown to be in excellent agreement with experiment. Unlike the FOOF and FON molecules, but similar to the cis- and trans-FONO molecules, FONO2 is shown to possess normal bond distances. The bonding in FONO2 is shown to be more similar to that in ClONO2 than that in HONO2, although there are still significant differences, especially in the partial atomic charges as deduced from Mulliken populations. This causes FONO2 to possess almost no dipole moment, which is very different to both HONO2 and ClONO2. By using large atomic natural orbital basis sets, CCSD(T) energies are computed for four isodesmic reactions in order to determine an accurate heat of formation for FONO2. Our best estimate for DELTA H(raised circle) (sub f,298) is 3.1 plus or minus 2.0 kcal/mol, indicating that the F-ONO2 bond energy is 31.3 kcal/mol.

Lee, Timothy J.↗

Recent Advances in Open-Shell Perturbation Theory and Coupled-Cluster Theory

Comparisons of various recently developed open-shell RHF perturbation theories will be presented. Among the aspects considered are spin-contamination, computational cost, and quality of numerical results. In addition, a new approach to avoid the disk storage and I/O bottlenecks in large scale coupled-cluster calculations will be discussed.

Lee, Timothy J.↗

Symmetry breaking in O4(+): An application of the Brueckner coupled-cluster method

A recent calculation of the antisymmetric stretch frequency for the rectangular structure of quartet O4(+) using the singles and doubles quadratic configuration interaction method with a perturbational estimate of connected triple excitations (QCISD(T)) method gave a value of 3710 cm(exp -1). This anomalous frequency is shown to be a consequence of symmetry breaking effects, which occur even though the QCISD(T) solution derived from a delocalized SCF reference function lies energetically well below the two localized (symmetry-broken) solutions at the equilibrium geometry. The symmetry breaking is almost eliminated at the CCSD level of theory, but the small remaining symmetry breaking effects are magnified at the CCSD(T) level of theory so that the antisymmetric stretch frequency is still significantly in error. The use of the Brueckner coupled cluster method, however, leads to a symmetrical solution which is free of symmetry breaking effects, with an antisymmetric stretch frequency of 1322 cm(exp -1), in good agreement with our earlier calculations using the complete active space self consistent field/complete active space state interaction (CASSCF/CASSI) method.

Barnes, Leslie A.↗

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

The Convergence of the Coupled Cluster Approach for MgO

The X(exp 1)Sigma(+) state of MgO shows significant differences between the CCSD(T) and CCSDT results and even some significant differences between the CCSDT and CCSDTQ results. The CCSDT(Q) results for this state do not agree well with the CCSDTQ results, as has been found for other systems. Unlike the X state, the state, which is well described by a single configuration, shows much less dependence on level of correlation treatment. Despite the slow convergence of the valence treatment of the X state, the effect of core correlation converges rapidly with level of correlation treatment.

Bauschlicher, Charles W.↗

The protonation of N2O reexamined - A case study on the reliability of various electron correlation methods for minima and transition states

The protonation of N2O and the intramolecular proton transfer in N2OH(+) are studied using various basis sets and a variety of methods, including second-order many-body perturbation theory (MP2), singles and doubles coupled cluster (CCSD), the augmented coupled cluster (CCSD/T/), and complete active space self-consistent field (CASSCF) methods. For geometries, MP2 leads to serious errors even for HNNO(+); for the transition state, only CCSD/T/ produces a reliable geometry due to serious nondynamical correlation effects. The proton affinity at 298.15 K is estimated at 137.6 kcal/mol, in close agreement with recent experimental determinations of 137.3 +/- 1 kcal/mol.

Martin, J. M. L.↗

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

Computational Thermochemistry of Jet Fuels and Rocket Propellants

The design of new high-energy density molecules as candidates for jet and rocket fuels is an important goal of modern chemical thermodynamics. The NASA Glenn Research Center is home to a database of thermodynamic data for over 2000 compounds related to this goal, in the form of least-squares fits of heat capacities, enthalpies, and entropies as functions of temperature over the range of 300 - 6000 K. The chemical equilibrium with applications (CEA) program written and maintained by researchers at NASA Glenn over the last fifty years, makes use of this database for modeling the performance of potential rocket propellants. During its long history, the NASA Glenn database has been developed based on experimental results and data published in the scientific literature such as the standard JANAF tables. The recent development of efficient computational techniques based on quantum chemical methods provides an alternative source of information for expansion of such databases. For example, it is now possible to model dissociation or combustion reactions of small molecules to high accuracy using techniques such as coupled cluster theory or density functional theory. Unfortunately, the current applicability of reliable computational models is limited to relatively small molecules containing only around a dozen (non-hydrogen) atoms. We propose to extend the applicability of coupled cluster theory- often referred to as the 'gold standard' of quantum chemical methods- to molecules containing 30-50 non-hydrogen atoms. The centerpiece of this work is the concept of local correlation, in which the description of the electron interactions- known as electron correlation effects- are reduced to only their most important localized components. Such an advance has the potential to greatly expand the current reach of computational thermochemistry and thus to have a significant impact on the theoretical study of jet and rocket propellants.

Crawford, T. Daniel↗

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

Singlet Excited States of Cl and Br Molecules: New Theories Applied to the -XO and -XO2 (X=C1 and Br) Chromophores

Electronic excitation energies are determined using single-reference based theories derived from response equations involving perturbation theory and coupled-cluster theory. These methods are applied to the singlet manifold of excited electronic states of the HClO, HBrO, HOClO, HOBrO, HClO2, and HBrO2 molecules. The reliability of the various perturbation theory approaches is assessed by comparison to the linear-response singles and doubles coupled-cluster (LRCCSD) method. The excitation energies for the Y-XO compounds are compared and contrasted for Y=H and HO, and X=Cl and Br. A similar comparison is performed for the H-XO2 compounds.

Lee, Timothy J.↗

Concerted hydrogen atom exchange between three HF molecules

We have investigated the termolecular reaction involving concerted hydrogen exchange between three HF molecules, with particular emphasis on the effects of correlation at the various stationary points along the reaction. Using an extended basis, we have located the geometries of the stable hydrogen-bonded trimer, which is of C(sub 3h) symmetry, and the transition state for hydrogen exchange, which is of D(sub 3h) symmetry. The energies of the exchange reation were then evaluated at the correlated level, using a large atomic natural orbital basis and correlating all valence electrons. Several correlation treatments were used, namely, configration interaction with single and double excitations, coupled-pair functional, and coupled-cluster methods. We are thus able to measure the effect of accounting for size-extensivity. Zero-point corrections to the correlated level energetics were determined using analytic second derivative techniques at the SCF level. Our best calculations, which include the effects of connected triple excitations in the coupled-cluster procedure, indicate that the trimer is bound by 9 +/- 1 kcal/mol relative to three separate monomers, in excellent agreement with previous estimates. The barrier to concerted hydrogen exchange is 15 kcal/mol above the trimer, or only 4.7 kcal/mol above three separated monomers. Thus the barrier to hydrogen exchange between HF molecules via this termolecular process is very low.

Komornicki, Andrew↗

A diagnostic for determining the quality of single-reference electron correlation methods

It was recently proposed that the Euclidian norm of the t sub 1 vector of the coupled cluster wave function (normalized by the number of electrons included in the correlation procedure) could be used to determine whether a single-reference-based electron correlation procedure is appropriate. This diagnostic, T sub 1, is defined for use with self consistent field molecular orbitals and is invariant to the same orbital rotations as the coupled cluster energy. T sub 1 is investigated for several different chemical systems which exhibit a range of multireference behavior, and is shown to be an excellent measure of the importance of nondynamical electron correlation and is far superior to C sub 0 from a singles and doubles configuration interaction wave function. It is further suggested that when the aim is to recover a large fraction of the dynamical electron correlation energy, a large T sub 1 (i.e., greater than 0.02) probably indicates the need for a multireference electron correlation procedure.

Lee, Timothy J.↗

A diagnostic for determining the quality of single-reference electron correlation methods

It was recently proposed that the Euclidian norm of the t(sub 1) vector of the coupled cluster wave function (normalized by the number of electrons included in the correlation procedure) could be used to determine whether a single-reference-based electron correlation procedure is appopriate. This diagnostic, T(sub 1) is defined for use with self-consistent-field molecular orbitals and is invariant to the same orbital rotations as the coupled cluster energy. T(sub 1) is investigated for several different chemical systems which exhibit a range of multireference behavior, and is shown to be an excellent measure of the importance of non-dynamical electron correlation and is far superior to C(sub 0) from a singles and doubles configuration interaction wave function. It is further suggested that when the aim is to recover a large fraction of the dynamical electron correlation energy, a large T(sub 1) (i.e., greater than 0.02) probably indicates the need for a multireference electron correlation procedure.

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

Initiating Molecular Growth in the Interstellar Medium via Dimeric Complexes of Observed Ions and Molecules

A feasible initiation step for particle growth in the interstellar medium (ISM) is simulated by means of ab quantum chemistry methods. The systems studied are dimer ions formed by pairing nitrogen containing small molecules known to exist in the ISM with ions of unsaturated hydrocarbons or vice versa. Complexation energies, structures of ensuing complexes and electronic excitation spectra of the encounter complexes are estimated using various quantum chemistry methods. Moller-Plesset perturbation theory (MP2, Z-averaged perturbation theory (ZAP2), coupled cluster singles and doubles with perturbative triples corrections (CCSD(T)), and density functional theory (DFT) methods (B3LYP) were employed along with the correlation consistent cc-pVTZ and aug-cc-pVTZ basis sets. Two types of complexes are predicted. One type of complex has electrostatic binding with moderate (7-20 kcal per mol) binding energies, that are nonetheless significantly stronger than typical van der Waals interactions between molecules of this size. The other type of complex develops strong covalent bonds between the fragments. Cyclic isomers of the nitrogen containing complexes are produced very easily by ion-molecule reactions. Some of these complexes show intense ultraviolet visible spectra for electronic transitions with large oscillator strengths at the B3LYP, omegaB97, and equations of motion coupled cluster (EOM-CCSD) levels. The open shell nitrogen containing carbonaceous complexes especially exhibit a large oscillator strength electronic transition in the visible region of the electromagnetic spectrum.

Bera, Partha P.↗

The Heat of Formation of HNO

The HNO molecule is of interest in both combustion and atmospheric chemistry. For example, Guadagnini et al. have recently presented ab initio potential energy surfaces for the three lowest lying electronic states of HNO and then used these in examining several chemical reactions that take place in the combustion of nitrogen containing fuels and in the oxidation of atmospheric nitrogen. We have previously studied the ground state potential energy surface (i.e., stationary points along the HNO rev. reaction HON path), vibrational spectrum (using an accurate quartic force field), zero-point energy, and bonding of HNO using coupled-cluster ab initio methods. HNO is also very interesting because of the unique nature of its bonding characteristics. That is, the potential energy surface is very flat along the H-N bonding coordinate thereby giving unusual harmonic and fundamental vibrational frequencies, and the H-N bond energy is rather weak in comparison to other H-N bond energies. In fact, using experimental heats of formation for HO, H, and NO, the H- bond energy is computed to be only 49.9 kcal/ mol (298 K). However, ab initio calculations of isodesmic reaction energies involving HNO, FNO, ClNO, and several other molecules have shown that there is an inconsistency in the experimental heats of formation of the XNO (X=H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate(DELTA)H(sup o)(sub f) value for HNO using state of-the-art ab initio methods. Based on many recent studies it is evident that the singles and doubles coupled-cluster method that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), in conjunction with large one-particle basis sets should be reliable to better than +0.8 kcal/mol for this quantity. The computational methodology is described in the next section followed by our results and discussion. Conclusions are presented in the final section.

Lee, Timothy J.↗