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

Exploring mechanisms of a tropospheric archetype: CH{sub 3}O{sub 2} + NO

Methylperoxy radical (CH{sub 3}O{sub 2}) and nitric oxide (NO) contribute to the propagation of photochemical smog in the troposphere via the production of methoxy radical (CH{sub 3}O) and nitrogen dioxide (NO{sub 2}). This reaction system also furnishes trace quantities of methyl nitrate (CH{sub 3}ONO{sub 2}), a sink for reactive NO{sub x} species. Here, the CH{sub 3}O{sub 2} + NO reaction is examined with highly reliable coupled-cluster methods. Specifically, equilibrium geometries for the reactants, products, intermediates, and transition states of the ground-state potential energy surface are characterized. Relative reaction enthalpies at 0 K (ΔH{sub 0K}) are reported; these values are comprised of electronic energies extrapolated to the complete basis set limit of CCSDT(Q) and zero-point vibrational energies computed at CCSD(T)/cc-pVTZ. A two-part mechanism involving CH{sub 3}O and NO{sub 2} production followed by radical recombination to CH{sub 3}ONO{sub 2} is determined to be the primary channel for formation of CH{sub 3}ONO{sub 2} under tropospheric conditions. Constrained optimizations of the reaction paths at CCSD(T)/cc-pVTZ suggest that the homolytic bond dissociations involved in this reaction path are barrierless.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Ultraviolet photodissociation of gas-phase iron pentacarbonyl probed with ultrafast infrared spectroscopy

It is well known that ultraviolet photoexcitation of iron pentacarbonyl results in rapid loss of carbonyl ligands leading to the formation of coordinatively unsaturated iron carbonyl compounds. Here we employ ultrafast mid-infrared transient absorption spectroscopy to probe the photodissociation dynamics of gas-phase iron pentacarbonyl following ultraviolet excitation at 265 and 199 nm. After photoexcitation at 265 nm, our results show evidence for sequential dissociation of iron pentacarbonyl to form iron tricarbonyl via a short-lived iron tetracarbonyl intermediate. Photodissociation at 199 nm results in the prompt production of Fe(CO) 3 within 0.25 ps via several energetically accessible pathways. An additional 15 ps time constant extracted from the data is tentatively assigned to intersystem crossing to the triplet manifold of iron tricarbonyl or iron dicarbonyl. Mechanisms for formation of iron tetracarbonyl, iron tricarbonyl, and iron dicarbonyl are proposed and theoretically validated with one-dimensional cuts through the potential energy surface as well as bond dissociation energies. Ground state calculations are computed at the CCSD(T) level of theory and excited states are computed with EOM-EE-CCSD(dT).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Enhancing the accuracy of XPS calculations: Exploring hybrid basis set schemes for CVS-EOMIP-CCSD calculations

Reliable computational methodologies and basis sets for modeling x-ray spectra are essential for extracting and interpreting electronic and structural information from experimental x-ray spectra. In particular, the trade-off between numerical accuracy and computational cost due to the size of the basis set is a major challenge, since molecular orbitals undergo extreme relaxation in the core-hole state. To gain clarity on the changes in electronic structure induced by the formation of a core-hole, the use of sufficiently flexible basis for expanding the orbitals, particularly for the core region, has been shown to be essential. This work focuses on the refinement of core-hole ionized state calculations using the equation-of-motion coupled cluster family of methods through an extensive analysis on the effectiveness of “hybrid” and mixed basis sets. In this investigation, we utilize the CVS-EOMIP-CCSD method in combination and construct hybrid basis sets piecewise from readily available Dunning’s correlation consistent basis sets in order to calculate x-ray ionization energies (IEs) for a set of small gas phase molecules. Our results provide insights into the impact of basis sets on the CVS-EOMIP-CCSD calculations of K-edge IEs of first-row p-block elements. Furthermore, these insights enable us to understand more about the basis set dependence of the core IEs computed and allow us to establish a protocol for deriving reliable and cost-effective theoretical estimates for computing IEs of small molecules containing such elements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Communication: Helium nanodroplet isolation and rovibrational spectroscopy of hydroxymethylene

Hydroxymethylene (HCOH) and its d{sub 1}-isotopologue (HCOD) are isolated in low temperature helium nanodroplets following pyrolysis of glyoxylic acid. Transitions identified in the infrared spectrum are assigned exclusively to the trans-conformation based on previously reported anharmonic frequency computations [P. R. Schreiner, H. P. Reisenauer, F. C. Pickard, A. C. Simmonett, W. D. Allen, E. Mátyus, and A. G. Császár, Nature 453, 906 (2008); L. Koziol, Y. M. Wang, B. J. Braams, J. M. Bowman, and A. I. Krylov, J. Chem. Phys. 128, 204310 (2008)]. For the OH(D) and CH stretches, a- and b-type transitions are observed, and when taken in conjunction with CCSD(T)/cc-pVTZ computations, lower limits to the vibrational band origins are determined. The relative intensities of the a- and b-type transitions provide the orientation of the transition dipole moment in the inertial frame. The He nanodroplet data are in excellent agreement with anharmonic frequency computations reported here and elsewhere, confirming an appreciable Ar-matrix shift of the OH and OD stretches and strong anharmonic resonance interactions in the high-frequency stretch regions of the mid-infrared.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗
DOE OSTI2024

Exploring the Tl 2H 2 potential energy surface: A comparative analysis with group 13 systems and experiment

Abstract Thallium chemistry is experiencing unprecedented importance. Therefore, it is valuable to characterize some of the simplest thallium compounds. Stationary points along the singlet and triplet TlH potential energy surface have been characterized. Stationary point geometries were optimized with the CCSD(T)/aug‐cc‐pwCVQZ‐PP method. Harmonic vibrational frequencies were computed at the same level of theory while anharmonic vibrational frequencies were computed at the CCSD(T)/aug‐cc‐pwCVTZ‐PP level of theory. Final energetics were obtained with the CCSDT(Q) method. Basis sets up to augmented quintuple‐zeta cardinality (aug‐cc‐pwCV5Z‐PP) were employed to obtain energetics in order to extrapolate to the complete basis set limits using the focal point approach. Zero‐point vibrational energy corrections were appended to the extrapolated energies in order to determine relative energies at 0 K. It was found that the planar dibridged isomer lies lowest in energy while the linear structure lies highest in energy. The results were compared to other group 13 MH (M = B, Al, Ga, In, and Tl) theoretical studies and some interesting variations are found. With respect to experiment, incompatibilities exist.

Chemistry↗
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.↗
DOE OSTI2014

Accurate molecular structures and infrared spectra of trans-2,3-dideuterooxirane, methyloxirane, and trans-2,3-dimethyloxirane

Oxirane derivatives are the most used benchmarks for chiroptical spectroscopies in view of their small size and relative rigidity. The molecular structure, vibrational harmonic and anharmonic frequencies, and infrared intensities of the ground electronic states are analyzed in this paper. Equilibrium structure and harmonic force fields have been evaluated by means of high-level quantum-chemical calculations at the coupled-cluster level including single and double excitations together with a perturbative treatment of triples (CCSD(T)). Extrapolation to the complete basis-set limit as well as core-correlation effects have also been taken into account. Anharmonic contributions have been computed at the CCSD(T)/cc-pVTZ level for trans-2,3-dideuterooxirane. These data can serve as references to evaluate the accuracy of less expensive computational approaches rooted in the density functional theory (DFT). The latter have been used within hybrid CC/DFT approaches, which have been applied to simulate fully anharmonic infrared (IR) spectra. Finally, the best theoretical estimates of the equilibrium structures and vibrational wavenumbers are compared to the most accurate experimental data and show in all cases very good agreement, i.e., within 0.001 Å, 0.1 deg, 10 cm{sup −1}, and 0.5 km mol{sup −1}, for bond lengths, angles, wavenumbers, and IR intensities, respectively.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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.↗
DOE OSTI2020

The atmospheric importance of methylamine additions to Criegee intermediates

We report that Criegee intermediates are important targets for study in atmospheric chemistry because of their capacity to oxidize airborne species. Among these species, ammonia has received critical attention for its presence in polluted agricultural or industrial areas and its role in forming particulate matter and condensation nuclei. Although methylamine has been given less attention than ammonia, both theoretical and experimental studies have demonstrated that the additional methyl substitution on the ammonia derivatives increases the rate constants for some systems. This suggests that the methylamine addition to Criegee intermediates could be more significant to atmospheric processes. In this work, geometries are optimized at the DF-CCSD(T)/ANO1 level for the methylamine addition reactions to the simplest Criegee intermediate and the anti- and syn-methylated Criegee intermediates. Energies for each stationary point were computed at the CCSD(T)/CBS level with corrections from the CCSDT(Q) method. Rate constants are obtained for each reaction using canonical transition state theory. Although methylamine addition proved to be a more favorable reaction relative to ammonia addition, the significantly lower concentration of atmospheric methylamine limits the prevalence of these reactions, even in the most optimal conditions. It is unlikely that the methylamine addition to Criegee intermediates will contribute significantly to the consumption of Criegee intermediates in the atmosphere.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Equation-of-motion coupled-cluster method with double electron-attaching operators: Theory, implementation, and benchmarks

We report a production-level implementation of the equation-of-motion (EOM) coupled-cluster (CC) method with double electron-attaching (DEA) EOM operators of 2 p and 3 p 1 h types, EOM-DEA-CCSD. This ansatz, suitable for treating electronic structure patterns that can be described as two-electrons-in-many orbitals, represents a useful addition to the EOM-CC family of methods. We analyze the performance of EOM-DEA-CCSD for energy differences and molecular properties. By considering reduced quantities, such as state and transition one-particle density matrices, we compare EOM-DEA-CCSD wave functions with wave functions computed by other EOM-CCSD methods. The benchmarks illustrate that EOM-DEA-CCSD is capable of treating diradicals, bond-breaking, and some types of conical intersections.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Optimizing Tensor Contractions in CCSD(T) for Efficient Execution on GPUs

Tensor contractions are higher dimensional analogs of matrix multiplications, used in many computational contexts such as high order models in quantum chemistry, deep learning, nite element methods etc. In contrast to the wide availability of high-performance libraries for matrix multiplication on GPUs, the same is not true for tensor contractions. In this paper, we address the optimization of a set of symmetrized tensor contractions that form the computational bottleneck in the CCSD(T) coupled-cluster method in computational chemistry suites like NWChem. Some of the challenges in optimizing tensor contractions that arise in practice from the variety of dimensionalities and shapes for tensors include effective mapping of the high-dimensional iteration space to threads, choice of data buffering in shared-memory and registers, and tile sizes for multi-level tiling. Furthermore, in the case of symmetrized tensor contractions in CCSD(T), it is also a challenge to fuse contractions to reduce data movement cost by exploiting reuse of intermediate tensors. In this paper, we developed an efficient GPU implementation of the tensor contractions in CCSD(T) using shared-memory buffering, register tiling, loop fusion and register transpose. Experimental results demonstrate significant improvement over the current state-of-the-art.

Kim, Jinsung↗
DOE OSTI2022

How to renormalize coupled cluster theory

Coupled cluster theory is an attractive tool to solve the quantum many-body problem because its singles and doubles (CCSD) approximation is computationally affordable and yields about 90% of the correlation energy. Capturing the remaining 10%, e.g., via including triples, is numerically expensive. In this work, we assume that short-range three-body correlations dominate and—following Lepage (arXiv:nucl-th/9706029)—that their effects can be included within CCSD by renormalizing the three-body contact interaction. We renormalize this contact in 16 O and obtain systematically improved CCSD results for 24 O, 20–34 Ne, 40,48 Ca, 78 Ni, 90 Zr, and 100 Sn.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Three-dimensional potential energy surfaces of ArNO ($\tilde{X}$ 2 Π)

Until now, the potential energy surfaces (PESs) of the ArNO complex found in the literature were two-dimensional, with the NO interatomic distance being fixed. In this work, we present the first accurate three-dimensional ground state X 2 Π PESs (both A ' and A ") of ArNO computed at the CCSD(T)/CBS level of theory. The equilibrium geometries and the well depths ( D e ) are compared to several other electronic structure methods. We found that using the multireference method, MRCI-F12 makes the surfaces much shallower (by 25%) and the depth of the surfaces does not agree with experimental data. The explicitly correlated coupled-cluster method underestimates the well depth as well. Analytic representations for both A ' and A " surfaces were fit to 4380 ab initio points to within 2.71 cm -1 .A three-dimensional Numerov propagator method in Delves coordinates is used to compute the bound state spectrum up to J tot = 6.5. The recommended dissociation energies are D 0 = 97.2 cm -1 for the adiabatic ground state and D e = 133.7 (128.1) cm -1 for A ' ( A ")

36 MATERIALS SCIENCE↗
DOE OSTI2023

Anion-Activated Bases and Nucleophiles Characterized by Photoelectron Spectroscopy

In this study, negative ion photoelectron spectra at 20 K along with ab initio (CCSD(T)) and M06-2X density function theory calculations are reported for a series of six basic and nucleophilic pyridine derivatives with an anionic substituent (i.e., 3- and 4-PyrBX 3 –, where X = F, 4-t-BuC 6 H 4 , 4- MeOC 6 H 4 , and 3,5-(MeO) 2 C 6 H 3 ). Vertical detachment energies spanning from 4.50 – 5.85 eV are observed and well reproduced by M06-2X/aug-cc-pVTZ and CCSD(T)/maug-cc-pVTZ computations. Surprisingly, the VDEs are found to correlate with the S N 2 reactivity of the PPh 4 + salts of the anions with 1-iodooctane. Spin densities of the neutral structures formed from the anionic species are also reported and indicate where the excess electron density in the anions resides and likely cation-interaction sites in their corresponding salts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae

Here, we present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5–50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called “thermalized” rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

The ethyl radical in superfluid helium nanodroplets: Rovibrational spectroscopy and ab initio computations

The ethyl radical has been isolated and spectroscopically characterized in 4He nanodroplets. The band origins of the five CH stretch fundamentals are shifted by < 2 cm –1 from those reported for the gas phase species. The symmetric CH 2 stretching band (v 1 ) is rotationally resolved, revealing nuclear spin statistical weights predicted by G 12 permutation-inversion group theory. A permanent electric dipole moment of 0.28 (2) D is obtained via the Stark spectrum of the v 1 band. The four other CH stretch fundamental bands are significantly broadened in He droplets and lack rotational fine structure. This broadening is attributed to symmetry dependent vibration-to-vibration relaxation facilitated by the He droplet environment. In addition to the five fundamentals, three a 1 ' overtone/combination bands are observed, and each of these have resolved rotational substructure. As a result, these are assigned to the 2v 12 , v 4 + v 6 , and 2v 6 bands through comparisons to anharmonic frequency computations at the CCSD(T)/cc-pVTZ level of theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

Theoretical spectroscopic characterization at low temperatures of S-methyl thioformate and O-methyl thioformate

Highly correlated ab initio methods are employed to determine spectroscopic properties at low temperatures of two S-analogs of methyl formate: S-methyl thioformate CH{sub 3}-S-CHO (MSCHO) and O-methyl thioformate CH{sub 3}-O-CHS (MOCHS). Both species are detectable and they are expected to play an important role in Astrochemistry. Molecular properties are compared with those of the O-analog, methyl formate. Both isomers present two conformers cis and trans. cis-CH{sub 3}-S-CHO represents the most stable structure lying 4372.2 cm{sup −1} below cis-CH{sub 3}-O-CHS. The energy difference between the cis and trans forms is drastically lower for MSCHO (1134 cm{sup −1}) than for MOCHS (1963.6 cm{sup −1}). Harmonic and anharmonic fundamentals and the corresponding intensities, as well as the rotational constants for the ground vibrational and first excited torsional states and the centrifugal distortions constants, are provided. Low torsional energy levels have been obtained by solving variationally a two dimensional Hamiltonian expressed in terms of the two torsional degrees of freedom. The corresponding 2D potential energy surfaces have been computed at the CCSD(T)/aug-cc-pVTZ level of theory. The methyl torsional barriers V{sub 3}(cis) are determined to be 139.7 cm{sup −1} (CH{sub 3}-S-CHO) and 670.4 cm{sup −1} (CH{sub 3}-O-CHS). The A/E splitting of ground torsional state has been estimated to be 0.438 cm{sup −1} for CH{sub 3}-S-CHO and negligible for CH{sub 3}-O-CHS.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗