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 55 records · Page 3

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 reversible 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 HNO, 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 double bond H, F, and Cl) species. Hence the motivation for this study was to determine a very accurate (Delta)H(sub f, sup o) 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.↗

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↗

The Structure and Stability of Bn(+) Clusters

The geometries of B+n clusters for n less than 14 have been optimized using density functional theory with the B3LYP functional. The most stable structure for each cluster is planar or quasi-planar. The B3LYP fragmentation energies are calibrated using coupled cluster theory. Overall, our corrected fragmentation energies are in reasonable agreement with experiment. Our results are compared with previous theoretical results.

Ricca, Alessandra↗

Theoretical investigations of the structures and binding energies of Be(sub n) and Mg(sub n) (n = 3-5) clusters

Researchers determined the equilibrium geometries and binding energies of Be and Mg trimers, tetramers and pentamers using single and double excitation coupled cluster (CCSD) and complete active space self-consistent-field (CASSCF) multireference configuration interaction (MRCI) wave functions in conjunction with extended atomic basis sets. Best estimates of the cluster binding energies are 24, 83 and 110 kcal/mole for Be3, Be4 and Be5; and 9, 31 and 41 kcal/mole for Mg3, Mg4 and Mg5, respectively. A comparison of the MRCI and CCSD results shows that even the best single-reference approach (limited to single and double excitations) is not capable of quantitative accuracy in determining the binding energies of Be and Mg clusters.

Lee, Timothy J.↗

Theoretical investigations of the structures and binding energies of Be(n) and Mg(n) (n = 3-5) clusters

The equilibrium geometries and binding energies of Be and Mg trimers, tetramers and pentamers have been determined using single and double excitation coupled cluster (CCSD) and complete active space self-consistent-field (CASSCF) multireference configuration interaction (MRCI) wave functions in conjunction with extended atomic basis sets. The best estimates of the cluster binding energies are 24, 83, and 110 kcal/mol for Be3, Be4, and Be5; and 9, 31, and 41 kcal/mol for Mg3, Mg4, and Mg5, respectively. A comparison of the MRCI and CCSD results shows that even the best single-reference approach (limited to single and double excitations) is not capable of quantitative accuracy in determining the binding energies of Be and Mg clusters.

Lee, Timothy J.↗

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

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

Pendulum Motion in Main Parachute Clusters

The coupled dynamics of a cluster of parachutes to a payload are notoriously difficult to predict. Often the payload is designed to be insensitive to the range of attitude and rates that might occur, but spacecraft generally do not have the mass and volume budgeted for this robust of a design. The National Aeronautics and Space Administration (NASA) Orion Capsule Parachute Assembly System (CPAS) implements a cluster of three mains for landing. During testing of the Engineering Development Unit (EDU) design, it was discovered that with a cluster of two mains (a fault tolerance required for human rating) the capsule coupled to the parachute cluster could get into a limit cycle pendulum motion which would exceed the spacecraft landing capability. This pendulum phenomenon could not be predicted with the existing models and simulations. A three phased effort has been undertaken to understand the consequence of the pendulum motion observed, and explore potential design changes that would mitigate this phenomenon. This paper will review the early analysis that was performed of the pendulum motion observed during EDU testing, summarize the analysis ongoing to understand the root cause of the pendulum phenomenon, and discuss the modeling and testing that is being pursued to identify design changes that would mitigate the risk.

Ray, Eric S.↗

Kramers' Restricted Closed Shell CCSD Theory

A Kramers' restricted version of the closed shell coupled cluster singles doubles theory is presented. The theory may be used in conjunction with 2 or 4-component relativistic reference wavefunctions. The intrinsic treatment of the spin-orbit coupling doubles the number of independent quantities (amplitudes and integrals) relative to a spin-independent formalism. The number of operations required to evaluate the equations is four times larger than in the optimal spin-independent closed shell formalism.

Visscher, Lucas↗

Vibrations in small Mg clusters

The equilibrium geometries, binding energies, and harmonic frequencies of Mg3 and Mg4 have been determined using large atomic natural orbital basis sets in conjunction with high levels of electron correlation. The correlation treatments comprise multireference configuration interaction (MRCI), singles and doubles coupled-cluster (CCSD) theory and the CCSD(T) extension that includes a perturbational estimate of connected triple excitations. As with previous studies of small Be clusters, the CCSD(T) method is found to reproduce the MRCI results with a remarkable degree of accuracy and at a fraction of the computational cost. Using the CCSD(T) method, full quartic force fields for Mg3 and Mg4 have been determined and anharmonic analyses have been performed using second-order perturbation theory. Vibration-rotation interaction constants and centrifugal distortion constants have also been determined.

Lee, Timothy J.↗

Quantum mechanical calculations to chemical accuracy

The accuracy of current molecular-structure calculations is illustrated with examples of quantum mechanical solutions for chemical problems. Two approaches are considered: (1) the coupled-cluster singles and doubles (CCSD) with a perturbational estimate of the contribution of connected triple excitations, or CCDS(T); and (2) the multireference configuration-interaction (MRCI) approach to the correlation problem. The MRCI approach gains greater applicability by means of size-extensive modifications such as the averaged-coupled pair functional approach. The examples of solutions to chemical problems include those for C-H bond energies, the vibrational frequencies of O3, identifying the ground state of Al2 and Si2, and the Lewis-Rayleigh afterglow and the Hermann IR system of N2. Accurate molecular-wave functions can be derived from a combination of basis-set saturation studies and full configuration-interaction calculations.

Bauschlicher, Charles W., Jr.↗

Cr2 Revisited

Abstract Cr2 is studied at the coupled cluster singles and doubles level of theory, including a perturbational estimate of the connected triples [CCSD(T)]. Both spin restricted, RCCSD(T), and spin unrestricted, UCCSD(T), approaches are used. The UCCSD(T) bond length is much longer than the experimental value, while the RCCSD(T) bond length is shorter. This difference is due to the inability of the spin restricted Hartree-Fock reference to describe the antiferromagnetic coupling at larger r values. In contrast with the results obtained using density functional theory in conjunction with older functionals, the results obtained with recently developed functionals are in agreement with the UCCSD(T) results.

Bauschlicher, Charles W., Jr.↗

Theoretical Study of the First Transition Row Oxides and Sulfides

The first transition row oxides and sulfides are studied using several different levels of theory. The calculations show the bonding mechanism in the sulfides and oxides to be very similar. For the oxides, accurate experimental data allow the theoretical methods to be calibrated. The same level of theory is used to study the sulfides where there is far less experimental information. For ScO through MnO and CuO the coupled cluster singles and doubles technique including a perturbational estimate of the unliked triple excitations [CCSD(T)] yields spectroscopic constants ((tau)e, (omega)e, and D0) in good agreement with experiment. The triple excitations are found to be very important in achieving this accuracy. For FeO to NiO, the single determinant self-consistent-field (SCF) approach yields pi orbitals that are localized on the metal or oxygen. This appears to cause problems for the single reference techniques; this is discussed in detail for NiO. The complete-active-space SCF/internally contracted averaged coupled pair functional approach (CASSCF/ICACPF) works well for FeO to NiO. The calculation of accurate dipole moments is found to be very difficult.

Bauschlicher, Charles W., Jt.↗

The Anharmonic Force Field of Ethylene, C2H4, by Means of Accurate Ab Initio Calculations

The quartic force field of ethylene, C2H4, has been calculated ab initio using augmented coupled cluster, CCSD(T), methods and correlation consistent basis sets of spdf quality. For the C-12 isotopomers C2H4, C2H3D, H2CCD2, cis-C2H2D2, trans-C2H2D2, C2HD3, and C2D4, all fundamentals could be reproduced to better than 10 per centimeter, except for three cases of severe Fermi type 1 resonance. The problem with these three bands is identified as a systematic overestimate of the Kiij Fermi resonance constants by a factor of two or more; if this is corrected for, the predicted fundamentals come into excellent agreement with experiment. No such systematic overestimate is seen for Fermi type 2 resonances. Our computed harmonic frequencies suggest a thorough revision of the accepted experimentally derived values. Our computed and empirically corrected re geometry differs substantially from experimentally derived values: both the predicted rz geometry and the ground-state rotational constants are, however, in excellent agreement with experiment, suggesting revision of the older values. Anharmonicity constants agree well with experiment for stretches, but differ substantially for stretch-bend interaction constants, due to equality constraints in the experimental analysis that do not hold. Improved criteria for detecting Fermi and Coriolis resonances are proposed and found to work well, contrary to the established method based on harmonic frequency differences that fails to detect several important resonances for C2H4 and its isotopomers. Surprisingly good results are obtained with a small spd basis at the CCSD(T) level. The well-documented strong basis set effect on the v8 out-of-plane motion is present to a much lesser extent when correlation-optimized polarization functions are used. Complete sets of anharmonic, rovibrational coupling, and centrifugal distortion constants for the isotopomers are available as supplementary material to the paper.

Martin, Jan M. L.↗

Ordering of the O-O stretching vibrational frequencies in ozone

The ordering of nu1 and nu3 for O3 is incorrectly predicted by most theoretical methods, including some very high level methods. The first systematic electron correlation method based on one-reference configuration to solve this problem is the coupled cluster single and double excitation method. However, a relatively large basis set, triple zeta plus double polarization is required. Comparison with other theoretical methods is made.

Scuseria, Gustavo E.↗

A high level Ab initio study of the anionic hydrogen-bonded complexes FH-CN-, FH-NC-, H2O-CN- and H2O-NC-

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗

A high-level ab initio study of the anionic hydrogen-bonded complexes FH-CN(-), FH-NC(-), H2O-CN(-), and H2O-NC(-)

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗

The polarizabilities of neon

The static dipole polarizability, second hyperpolarizability, quadrupole polarizability and dipole-dipole-quadrupole hyperpolarizability for the Ne atom are computed, using large Gaussian basis sets and coupled cluster wavefunctions. Correlation effects are substantial, ranging from 15 percent for the dipole polarizability to 40 percent for the second hyperpolarizability. It is confirmed by explicit calculation earlier hypotheses that g-type and higher angular momentum functions and core correlation effects make almost no contribution to these properties in Ne. Triple excitations, on the other hand, are very important, accounting for as much as 25 percent of the correlation contribution to the second hyperpolarizability. The best estimate of the second hyperpolarizability is 119 + or - 4 au, in good agreement with the recent calculations of Maroulis and Thakkar (1989) and the latest experimental result of Shelton (1989).

Taylor, Peter R.↗