Search NASA⌕ Search

SEARCH · Search NASA

Results for “Coupled cluster theory”

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 235 records · Page 13

Alkali Metal Control of Triplet-Mediated C–H Activation in Iron-Mediated Coupling of Dinitrogen and Benzene

Direct coupling of N 2 with abundant feedstocks like benzene to form N-containing organic compounds is a promising strategy for N 2 fixation pathways. The challenge of coupling N 2 activation and C–H bond oxidative addition was recently solved by introducing a reversible benzene C–H bond activation process mediated by a low-valent Fe(0) complex, which gave an organometallic product that could couple with partially reduced N 2 . Interestingly, the energetics of the C–H oxidative addition/reductive elimination step depends on the choice of alkali metal. However, the reason why the alkali metal influences the C–H bond activation remained elusive. In this work, we present a comprehensive study on this Fe(0)-mediated reversible C–H activation. Through density functional theory combined with high-level coupled cluster calculations, we discovered that the intermediate-spin triplet (S = 1) controls the energy of the transition state for C–H cleavage, while the high-spin quintet (S = 2) controls the position of the equilibrium. Na + drives the equilibrium toward oxidative addition due to an electrostatic effect, while K + and Rb + are dominated by a steric effect that favors the iron(0) species. The key role played by nonbonding interactions in the Fe-mediated C–H activation provides a conceptual model for alkali control over organometallic transformations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electric dipole polarizability of 58 Ni

The electric dipole strength distribution in 58 Ni between 6 and 20 MeV has been determined from proton inelastic scattering experiments at very forward angles at RCNP, Osaka. The experimental data are rather well reproduced by quasiparticle random-phase approximation calculations including vibration coupling, despite a mild dependence on the adopted Skyrme interaction. They allow an estimate of the experimentally inaccessible high-energy contribution above 20 MeV, leading to an electric dipole polarizability α D ⁡ ( 58 Ni) = 3.48 (31)⁢ fm 3 . This serves as a test case for recent extensions of coupled-cluster calculations with chiral effective field theory interactions to nuclei with two nucleons on top of a closed-shell system.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Two-Neutron Halo is Unveiled in F 29

Here, we report the measurement of reaction cross sections ( σ R ex ) of F 27 , 29 with a carbon target at RIKEN. The unexpectedly large σ R ex and derived matter radius identify F 29 as the heaviest two-neutron Borromean halo to date. The halo is attributed to neutrons occupying the 2 p 3 / 2 orbital, thereby vanishing the shell closure associated with the neutron number N = 20 . The findings are explained by state-of-the-art shell model calculations. Coupled-cluster computations based on effective field theories of the strong nuclear force describe the matter radius of F 27 but are challenged for F 29 .

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

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

An ab initio study of HCuCO

HCuCO is studied using a large Gaussian basis set at the coupled cluster singles and doubles level of theory, including a perturbational estimate of the connected triples (CCSD(T)). In contrast with CuCO, HCuCO is linear. The Cu-CO bond in HCuCO is significantly stronger than in CuCO. These differences between HCuCO and CuCO are discussed in terms of theCu-H bond polarizing the Cu 4s electron away from the CO.

Bauschlicher, Charles W., Jr.↗

Mechanism of superconductivity in the Hubbard model at intermediate interaction strength

We study the fluctuations responsible for pairing in the d -wave superconducting state of the two-dimensional Hubbard model at intermediate coupling within a cluster dynamical mean-field theory with a numerically exact quantum impurity solver. By analyzing how momentum- and frequency-dependent fluctuations generate the d -wave superconducting state in different representations, we identify antiferromagnetic fluctuations as the pairing glue of superconductivity in both the underdoped and the overdoped regime. Nevertheless, in the intermediate coupling regime, the predominant magnetic fluctuations may differ significantly from those described by conventional spin fluctuation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Real-time coupled-cluster approach for the cumulant Green's function

Green’s function methods within many-body perturbation theory provide a general framework for treating electronic correlations in excited states. Here we investigate the cumulant form of the one-electron Green’s function based on the coupled-cluster equation of motion approach in an extension of our previous study. The approach yields a non-perturbative expression for the cumulant in terms of the solution to a set of coupled first order, non-linear differential equations. The method thereby adds non-linear corrections to traditional cumulant methods linear in the self energy. The approach is applied to the core-hole Green’s function and illustrated for a number of small molecular systems. For these systems we find that the non-linear contributions lead to significant improvements both for quasiparticle properties such as core-level binding energies, as well as the satellites corresponding to inelastic losses observed in photoemission spectra.

Vila, Fernando D.↗

Quadratic Unitary Coupled-Cluster Singles and Doubles Scheme: Efficient Implementation, Benchmark Study, and Formulation of an Extended Version

An efficient implementation of the quadratic unitary coupled-cluster singles and doubles (qUCCSD) scheme for calculations of electronic ground and excited states using an unrestricted molecular spin–orbital formulation and an efficient tensor contraction library is reported. The accuracy of the qUCCSD scheme and the efficiency of the present implementation are demonstrated using extensive benchmark calculations of excitation energies and an application to S o → S 1 vertical excitation energies for cis- and trans-4a,4b-dihydrotriphenylene. The qUCCSD scheme has been shown to provide improved excitation energies compared with the UCC3 scheme formulated based on perturbation theory. Here, a UCC truncation scheme that can provide excitation energies correct through the fourth order is also presented to further improve the accuracy of the qUCCSD scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The behavior of methane–water mixtures under elevated pressures from simulations using many-body potentials

Non-polarizable empirical potentials have been proven to be incapable of capturing the mixing of methane–water mixtures at elevated pressures. Although density functional theory-based ab initio simulations may circumvent this discrepancy, they are limited in terms of the relevant time and length scales associated with mixing phenomena. Here, we show that the many-body MB-nrg potential, designed to reproduce methane–water interactions with coupled cluster accuracy, successfully captures this phenomenon up to 3 GPa and 500 K with varying methane concentrations. Two-phase simulations and long time scales that are required to fully capture the mixing, affordable due to the speed and accuracy of the MBX software, are assessed. Constructing the methane–water equation of state across the phase diagram shows that the stable mixtures are denser than the sum of their parts at a given pressure and temperature. As a result, we find that many-body polarization plays a central role, enhancing the induced dipole moments of methane by 0.20 D during mixing under pressure. Overall, the mixed system adopts a denser state, which involves a significant enthalpic driving force as elucidated by a systematic many-body energy decomposition analysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗

Dynamics of Density Fluctuations in Atomic Nuclei

We study the spatiotemporal patterns of density fluctuations in 16,24O and 48Ca using nuclear interactions from chiral effective field theory and the time-dependent coupled-cluster method. We find that two-particle–two-hole excitations generate small-amplitude fluctuations that are fast, short-ranged, and of stochastic character.

Bonaiti, Francesca [ORNL] (ORCID:0000000239261609)↗

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

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

The structure and energetics of Cr(CO)6 and Cr(CO)5

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. 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. Calculations using larger basis sets reduce 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. In the largest basis set, the total CO binding energy of Cr(CO)6 is estimated to be 140 kcal/mol at the CCSD(T) level of theory, or about 86 percent of the experimental value. The remaining discrepancy between the experimental and theoretical value 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 set (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.↗

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

Structure and energetics of Cr(CO)6 and Cr(CO)5

The geometric structures and energetics of Cr(CO)6 and Cr(CO)5 are determined at the modified coupled-pair functional, single and double excitation coupled-cluster (CCSD), and CCSD(T) levels of theory. For Cr(CO)6, the structure and force constants for the totally symmetric representation are in good agreement with experimental data once basis set constants are taken into account. In the largest basis set at the CCSD(T) level of theory, the total binding energy of CR(CO)6 is estimated at around 140 kcal/mol, or about 86 percent of the experimental value. In contrast, the first bond energy of Cr(CO)6 is very well described at the CCSD(T) level of theory, with the best estimated value of 38 kcal/mol being within the experimental uncertainty.

Barnes, Leslie A.↗

Theoretical Characterization of the NNO...O3 and NNO...O Complexes

The bonding mechanisms and bond strengths of complexes between two molecular and/or atomic species will be discussed. The implications of complex formation for stratospheric chemistry will be presented. Results of theoretical calculations on NNO ... O3 and NNO ...O will be used as specific examples to illustrate the principles of complex formation. The equilibrium structures and binding energies of these complexes will be presented as well as shifts in other molecular properties due to complexation, such as vibrational frequencies. Thermal dissociation rates will be presented. The theoretical methodology used in this study consists of high-level coupled-cluster calculations and more approximate density functional theory approaches, together with large one-particle basis sets.

Lee, Timothy J.↗