Search NASA⌕ Search

SEARCH · Search NASA

Results for “CCSD Computations”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

The Successive OH Binding Energies of Sc(OH)n+ for n=1-3

The geometries of Sc(OH)n+, for n = 1-3, have been optimized using density functional theory, in conjunction with the B3LYP hybrid functional. The zero-point energies are computed at the same level of theory. The successive OH bond energies have been computed at the CCSD(T) level for ScOH+ and Sc(OH)2+. The computed result for ScOD+ is in excellent agreement with the recent experiment of Armentrout and co-workers. There is a dramatic drop for the third OH, because Sc+ has only two valence electrons and therefore the bonding changes when the third OH is added. The difference between the B3LYP and CCSD(T) OH binding energies for the first two OH groups is discussed.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1997

Comparison of DFT with Traditional Methods for the Calculation of Vibrational Frequencies and Bond Energies

The vibrational frequencies of MO2 are computed at many levels of theory, including HF, B3LYP, BP86, CASSCF, MP2, and CCSD(T). The computed results are compared with the available experimental results. Most of the methods fail for at least one state of the systems considered. The accuracy of the results and the origin of the observed failures are discussed. The B3LYP bond energies are compared with traditional methods for a variety of systems, ranging from FeCOn+ to SiCln and its positive ions. The cases where B3LYP differs from the traditional methods are discussed.

Bauschlicher, Charles W.↗
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 NTRS1990

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

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

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

The Structure and Stability of B(n)H(+) Clusters

The geometries of the B(n)H(+) clusters for n=1-13 have been optimized at the B3LYP level of theory. Excluding B2H(+) all clusters are low-spin coupled. The structures appear to be related to the bare B(+)(n) clusters in two ways: 1) an H atom is added to the B(+)(n) cluster or 2) a BH is added to the B(+)(n-1) cluster. The B(+)(n)-H binding energies are computed at the B3LYP level of theory, and calibrated n using the CCSD(T) level of theory. The computed results tend to fall between the experimental best estimates and the experimental lower bounds. The B(+)(n)-H binding energies show an inverse correspondence with the B+(n-1)-B values) while the B(+)(n-1)-B and B(+)(n-1)-BH results tend to parallel each other.

Ricca, Alessandra↗
NASA NTRS1997

Positive Anharmonicities: The Oxonide Anion as an Example

An accurate ab initio quartic force field for the ozonide anion has been determined at the singles and doubles coupled-cluster level of theory that includes a perturbational estimate of the effects of connected triple excitations, denoted CCSD(T), using the augmented valence triple-zeta correlation consistent one-particle basis set. Convergence of the harmonic frequencies with respect to the one-particle basis set has been demonstrated by computing quadratic force fields at the CCSD(T) level using augmented valence double-zeta and quadruple-zeta basis sets. Fundamental vibrational frequencies have been determined via second-order ro-vibrational perturbation theory and also using exact variational methods. Agreement is very good, and they both show that the antisymmetric O-O stretch, v 3, possesses a positive anharmonic correction (that is, the fundamental frequency is predicted to be higher in energy than the harmonic frequency). Comparison of the O_3 and O3 quartic force fields shows that the positive anharmonic correction is the result of a particularly large F3333 symmetry internal coordinate force constant. The reasonableness of this force constant has been tested by computing both F33 and F3333 at the CCSD, CCSD(T), and CASPT2 levels of theory. A discussion of known positive anharmonicities for stretching modes is presented.

Lee, Timothy J.↗
NASA NTRS1995

The Ozonide Anion: A Theoretical Study

Complete active space self-consistent field (CASSCF) and CASSCF second order perturbation theory (CASPT2) methods have been used for the geometry optimization and calculation of harmonic and fundamental frequencies of the ozonide ion O3(-) and the ozonide lithium complex Li(+)O3(-). For O3(-) harmonic frequencies have also been obtained using the coupled-cluster method, CCSD(T). Infrared intensities are computed from dipole moment derivatives at the CASSCF level. The predicted equilibrium geometry for O3(-) is: Roo = 1.361 Angstroms and angle ooo = 115.4 degrees, and the fundamental frequencies are: nu(sub 1) = 989 per centimeter, nu(sub 2) 556 per centimeter, nu(sub 3) = 870 per centimeter (experimental values are: Roo =1.36 plus or minus 0.02 Angstroms, angle (ooo) = 111.8 plus or minus 2.0 degrees, nu(sub 1) = 975(50) per centimeter, nu(sub 2) = 550(50) per centimeter, nu(sub 3) = 880(50) per centimeter). Corresponding data for the lithium ozonide complex have also been obtained. The presented data contradict the previous interpretation of the IR and Raman spectrum obtained after deposition of ozone in N2, argon, or neon matrices with atomic beams of alkali metals. The presence of the lithium cation raises the asymmetric stretch frequency to about 940 per centimeter, which is contradictory to assumptions made in the assignments of the matrix spectra. Calculations made in a dielectric medium strongly suggest that the effect of the matrix on the IR spectrum is small for O3(-) itself. The dissociation and atomization energies of O3(-) are found to be in agreement with experiment.

Borowski, Poitr↗
NASA NTRS2003

Anharmonic Vibrational Spectroscopy of the F-(H20)n, complexes, n=1,2

We report anharmonic vibrational spectra (fundamentals, first overtones) for the F-(H(sub 2)O) and F-(H(sub 2)O)2 clusters computed at the MP2 and CCSD(T) levels of theory with basis sets of triple zeta quality. Anharmonic corrections were estimated via the correlation-corrected vibrational self-consistent field (CC-VSCF) method. The CC-VSCF anharmonic spectra obtained on the potential energy surfaces evaluated at the CCSD(T) level of theory are the first ones reported at a correlated level beyond MP2. We have found that the average basis set effect (TZP vs. aug-cc-pVTZ) is on the order of 30-40 cm(exp -1), whereas the effects of different levels of electron correlation [MP2 vs. CCSD(T)] are smaller, 20-30 cm(exp -1). However, the basis set effect is much larger in the case of the H-bonded O-H stretch of the F-(H(sub 2)O) cluster amounting to 100 cm(exp -1) for the fundamentals and 200 cm (exp -1) for the first overtones. Our calculations are in agreement with the limited available set of experimental data for the F-(H(sub 2)O) and F-(H(sub 2)O)2 systems and provide additional information that can guide further experimental studies.

Chaban, Galina M.↗
NASA NTRS1992

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

Ames-2021 CO 2 Dipole Moment Surface and IR Line Lists: Toward 0.1% Uncertainty for CO 2 IR Intensities

A highly accurate CO 2 ab initio dipole moment surface (DMS), Ames-2021, is reported along with 12 C 16 O 2 infrared (IR) intensity comparisons approaching a 1–4‰ level of agreement and uncertainty. The Ames-2021 DMS was accurately fitted from CCSD(T) finite-field dipoles computed with the aug-cc-pVXZ (X = T, Q, 5) basis for C atom and the d-aug-cc-pVXZ (X = T, Q, 5) basis for O atoms, and extrapolated to the one particle basis set limit. Fitting σ(rms) is 3.8 × 10^(–7) au for 4443 geometries below 15 000 per cm . The corresponding IR intensity, S Ames-2021 , are computed using the Ames-2 potential energy surface (PES), which is the best PES available for CO 2 . Compared to high accuracy IR studies for 2001i–00001 and 3001i–00001 bands, S Ames-2021 matches NIST experiment-based intensities [S NIST-HIT16 or S HIT20 ] to −1.0 ± 1.3‰, or matches DLR experiment-based intensities [S DLR-HIT16/UCL/Ames ] to 1.9 ± 3.7‰. This indicates the systematic deviations and uncertainties have been significantly reduced in S Ames-2021 . The S UCL2015 (or S HITRAN2016 ) have larger deviations (vs S DLR ) and uncertainties (vs S DLR , S NIST ) which are attributed to the less accurate Ames-1 PES adopted in UCL-296 line list calculation. The S Ames-2021 intensity of 12 C 16 O 2 and 13 C 16 O 2 is utilized to derive new absolute 13 C/ 12 C ratios for Vienna PeeDee Belemnite (VPDB) with uncertainty reduced by 1/3 or 2/3. Further evaluation of S Ames-2021 are carried out on those CO 2 bands discussed in the HITRAN2020 update paper. Consistent improvements and better accuracies are found in band-by-band analysis, except for those bands strongly affected by Coriolis couplings, or very weak bands measured with relatively larger experimental uncertainties. The Ames-2021 296 K IR line lists are generated for 13 CO 2 isotopologues, with 18,000 per and S 296K > 1E–31) cm/molecule cutoff and then combined with CDSD line positions (except 14 C 16 O 2 ). The Ames-2021 DMS and 296 K IR line lists represent a major improvement over previous CO 2 theoretical IR intensity studies, including Ames-2016, UCL-296, and recent UCL DMS 2021 update. A real 1 permille level of agreement and uncertainty will definitely require both more accurate PES and more accurate DMS.

Chemical calculations↗
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 NTRS2002

An Accurate Quartic Force Field, Fundamental Frequencies, and Binding Energy for the High Energy Density Material T(d)N4

The CCSD(T) method has been used to compute a highly accurate quartic force field and fundamental frequencies for all N-14 and N-15 isotopomers of the high energy density material T(sub d)N(sub 4). The computed fundamental frequencies show beyond doubt that the bands observed in a matrix isolation experiment by Radziszewski and coworkers are not due to different isotopomers of T(sub d)N(sub 4). The most sophisticated thermochemical calculations to date yield a N(sub 4) -> 2N(sub 2) heat of reaction of 182.22 +/- 0.5 kcal/mol at 0 K (180.64 +/- 0.5 at 298 K). It is hoped that the data reported herein will aid in the ultimate detection of T(sub d)N(sub 4).

Lee, Timothy J.↗
NASA NTRS1994

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

The Effect of Approximating Some Molecular Integrals in Coupled-Cluster Calculations: Fundamental Frequencies and Rovibrational Spectroscopic Constants of Cyclopropenylidene

The singles and doubles coupled-cluster method that includes a perturbational estimate of connected triple excitations, denoted CCSD(T), has been used, in conjunction with approximate integral techniques, to compute highly accurate rovibrational spectroscopic constants of cyclopropenylidene, C3H2. The approximate integral technique was proposed in 1994 by Rendell and Lee in order to avoid disk storage and input/output bottlenecks, and today it will also significantly aid in the development of algorithms for distributed memory, massively parallel computer architectures. It is shown in this study that use of approximate integrals does not impact the accuracy of CCSD(T) calculations. In addition, the most accurate spectroscopic data yet for C3H2 is presented based on a CCSD(T)/cc-pVQZ quartic force field that is modified to include the effects of core-valence electron correlation. Cyclopropenylidene is of great astronomical and astrobiological interest because it is the smallest aromatic ringed compound to be positively identified in the interstellar medium, and is thus involved in the prebiotic processing of carbon and hydrogen. The singles and doubles coupled-cluster method that includes a perturbational estimate of

Lee, Timothy J.↗
NASA NTRS2022

Energy Landscape, and Structural and Spectroscopic Characterization of Diazirine and Its Cyclic Isomers

Identifying new nitrogenated hydrocarbon molecules in the interstellar medium (ISM) is challenging due to the lack of comprehensive spectroscopic data from experiments. In this computational work, we focus on investigating the structures, relative energies, spectroscopic constants, and energy landscape of the cyclic isomers of diazirine (c-CH2N2) using ab initio quantum chemical methods. Density functional theory (DFT) methods and coupled cluster theory with singles and doubles including perturbative triples [CCSD(T)] and CCSD(T) with the explicitly correlated F12 b correction [CCSD(T)-F12b] were employed for this purpose along with large correlation consistent cc-pVTZ, cc-pVQZ, and cc-pV5Z basis sets. Harmonic vibrational frequencies, infrared vibrational intensities, rotational constants, and dipole moments are reported. An harmonic vibrational fundamentals along with centrifugal distortion constants, and vibration-rotation interaction constants are also reported for all the cyclic isomers. The energies computed with the CCSD(T) and CCSD(T)-F12b methods were extrapolated to the one-particle complete basis set (CBS) limit following a three-point formula. At the CCSD(T)-F12b/CBS level of theory the 3,3H-diazirine (c-CH2N2) is the lowest energy cyclic isomer followed by 1,3H-diazirine, (E)-1,2H-Diazirine, and (Z)-1,2H-diazirine which are 20.1, 47.8, and 51.3kcal mol-1above the 3,3H-Diazirine, respectively. Accurate structures and spectroscopic constants that are reported here could be useful for future identification of these cyclic nitrogenated organic molecules in the interstellar medium or circumstellar disks.

Energy Landscape↗
NASA NTRS2004

Vibrational Spectroscopy of Ions and Radicals Present in the Interstellar Medium and in Planetary Atmospheres: A Theoretical Study

Anharmonic vibrational frequencies and intensities are calculated for OH(H2O)n and H(H2O)n radicals (that form on icy particles of the interstellar medium), HCO radical (the main intermediate in the synthesis of organic molecules in space), NH2(-) and C2H(-) anions, H5(+) cation, and other systems relevant to interstellar chemistry. In addition to pure ions and radicals, their complexes with water are studied to assess the effects of water environment on infrared spectra. The calculations are performed using the correlation-corrected vibrational self-consistent field (CC-VSCF) method with ab initio potential surfaces at the MP2 and CCSD(T) levels. Fundamental, overtone, and combination excitations are computed. The results are in good agreement with available experimental data and provide reliable predictions for vibrational excitations not yet measured in laboratory experiments. The data should be useful for interpretation of astronomically observed spectra and identification of ions and radicals present in the interstellar medium and in planetary atmospheres.

Chaban, Galina M.↗