Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled cluster method”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Projected coupled cluster theory: Optimization of cluster amplitudes in the presence of symmetry projection

Methods which aim at universal applicability must be able to describe both weak and strong electronic correlation with equal facility. Such methods are in short supply. The combination of symmetry projection for strong correlation and coupled cluster theory for weak correlation offers tantalizing promise to account for both on an equal footing. In order to do so, however, the coupled cluster portion of the wave function must be optimized in the presence of the symmetry projection. Finally, this paper discusses how this may be accomplished, and shows the importance of doing so for both the Hubbard model Hamiltonian and the molecular Hamiltonian, all with a computational scaling comparable to that of traditional coupled cluster theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

A new ab initio potential energy surface and infrared spectra for the Ar–CS{sub 2} complex

We report a new three-dimensional potential energy surface for Ar–CS{sub 2} involving the Q{sub 3} normal mode for the υ{sub 3} antisymmetric stretching vibration of the CS{sub 2} molecule. The potential energies were calculated using the supermolecular method at the coupled-cluster singles and doubles level with noniterative inclusion of connected triples, using augmented correlation-consistent quadruple-zeta basis set plus midpoint bond functions. Two vibrationally averaged potentials with CS{sub 2} at both the ground (υ = 0) and the first excited (υ = 1)υ{sub 3} vibrational states were generated from the integration of the three-dimensional potential over the Q{sub 3} coordinate. Each potential was found to have a T-shaped global minimum and two equivalent linear local minima. The radial discrete variable representation /angular finite basis representation method and the Lanczos algorithm were applied to calculate the rovibrational energy levels. The calculated band origin shift of the complex (0.0622 cm{sup −1}) is very close to the observed one (0.0671 cm{sup −1}). The predicted infrared spectra and spectroscopic parameters based on the two averaged potentials are in excellent agreement with the available experimental data.

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

Probing ultrafast ππ*/nπ* internal conversion in organic chromophores via K-edge resonant absorption

Many photoinduced processes including photosynthesis and human vision happen in organic molecules and involve coupled femtosecond dynamics of nuclei and electrons. Organic molecules with heteroatoms often possess an important excited-state relaxation channel from an optically allowed ππ* to a dark nπ* state. The ππ*/nπ* internal conversion is difficult to investigate, as most spectroscopic methods are not exclusively sensitive to changes in the excited-state electronic structure. Here, we report achieving the required sensitivity by exploiting the element and site specificity of near-edge soft X-ray absorption spectroscopy. As a hole forms in the n orbital during ππ*/nπ* internal conversion, the absorption spectrum at the heteroatom K-edge exhibits an additional resonance. We demonstrate the concept using the nucleobase thymine at the oxygen K-edge, and unambiguously show that ππ*/nπ* internal conversion takes place within (60 ± 30) fs. Furthermore, high-level-coupled cluster calculations confirm the method’s impressive electronic structure sensitivity for excited-state investigations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Fermi–Löwdin orbital self-interaction correction of adsorption energies on transition metal ions

Density functional theory (DFT)-based descriptions of the adsorption of small molecules on transition metal ions are prone to self-interaction errors. Here, we show that such errors lead to a large over-estimation of adsorption energies of small molecules on Cu + , Zn + , Zn 2+ , and Mn + in local spin density approximation (LSDA) and Perdew, Burke, Ernzerhof (PBE) generalized gradient approximation calculations compared to reference values computed using the coupled-cluster with single, doubles, and perturbative triple excitations method. These errors are significantly reduced by removing self-interaction using the Perdew–Zunger self-interaction correction (PZ-SIC) in the Fermi–Löwdin Orbital (FLO) SIC framework. In the case of FLO-PBE, typical errors are reduced to less than 0.1 eV. Furthermore, analysis of the results using DFT energies evaluated on self-interaction-corrected densities [DFT(@FLO)] indicates that the density-driven contributions to the FLO-DFT adsorption energy corrections are roughly the same size in DFT = LSDA and PBE, but the total corrections due to removing self-interaction are larger in LSDA.

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

Anion photoelectron spectroscopy and chemical bonding of ThS 2 – and ThSO –

Anion photoelectron spectra of ThSO – and ThS 2 – were recorded using the third (355 nm) harmonic of an Nd-YAG laser; these provided the measured vertical detachment energies of each anion. The experiments are supported by extensive coupled cluster calculations on ThSO, ThSO – , ThS 2 , and ThS 2 – , as well as the oxygen congeners ThO 2 and ThO 2 – . The ab initio calculations, which included complete basis set extrapolations, spin–orbit effects using four-component coupled cluster, and higher-order correlation contributions through CCSDT(Q), yielded an adiabatic electron affinity for ThO 2 – that was within 0.02 eV of the previously determined experimental value. The singly occupied molecular orbital (SOMO) in all three anions corresponds primarily to the 7s orbital on Th. Successive substitution of S for each O in ThO 2 leads to larger electron affinities and smaller bond angles in the neutral molecules, but larger angles in the anions. As demonstrated by Franck–Condon simulations of the spectra using the CCSD(T) spectroscopic constants, substitution of O by S significantly complicates the resulting detachment spectra due to the lower vibrational frequencies in the sulfur species. Altogether the calculated vertical detachment energies are in very good agreement with the experiment. In addition to the adiabatic electron affinities of each species, atomization energies and heats of formation have also been determined via the FPD approach with expected uncertainties of 1–2 kcal/mol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1994

Electron Correlation in 4-Component Relativistic Calculations

The full 4-component Dirac-Coulomb equation can nowadays be used in molecular calculations, The first step in solving this relativistic many-electron equation usually consists of solving the closed or open-shell Diarc-Fock equations. Like in non-relativistic calculations the outcome does not account for the effects of electron correlation. This can in principle be remedied by developing relativistic variants of electron correlation methods like Configuration Interaction or Coupled Cluster. In this talk the differences and similarities of such relativistic approaches as compared to non-relativistic methods will be reviewed. Results of Configuration Interaction calculations on the PtH molecule and on the MeF(sub 6, sup 2-) (Me= Co, Rh, Ir) complexes will be presented to give an impression of the kind of results that currently can be obtained.

Visscher, Luuk↗
DOE OSTI2025

Slimmer Geminals For Accurate F12 Electronic Structure Models

The Slater-type F12 geminal length scales originally tuned for the second-order Mo̷ller-Plesset F12 method are too large for higher-order F12 methods formulated using the SP (diagonal fixed-coefficient spin-adapted) F12 ansatz. The new geminal parameters reported herein reduce the basis set incompleteness errors (BSIEs) of absolute coupled-cluster singles and doubles F12 correlation energies by a significant─and increase with the cardinal number of the basis─margin. The effect of geminal reoptimization is especially pronounced for the cc-pVXZ-F12 basis sets (specifically designed for use with F12 methods) relative to their conventional aug-cc-pVXZ counterparts. The BSIEs of relative energies are less affected, but substantial reductions can be obtained, especially for atomization energies and ionization potentials with the cc-pVXZ-F12 basis sets. The new geminal parameters are therefore recommended for all applications of high-order F12 methods, such as coupled-cluster F12 methods and transcorrelated F12 methods.

Powell, Samuel R. [Virginia Polytechnic Inst. and ↗
DOE OSTI2014

Second-order many-body perturbation and coupled-cluster singles and doubles study of ice VIII

The structure, equation of state, IR, Raman, and inelastic neutron scattering (INS) spectra of high-pressure, proton-ordered phase VIII of ice are calculated by the second-order many-body perturbation and coupled-cluster singles and doubles methods. Nearly all the observed features of the pressure-dependence of the structures and spectra are reproduced computationally up to 60 GPa insofar as the anharmonic effects can be neglected. The calculations display no sign of the hypothetical isostructural transition in 2–3 GPa to phase VIII{sup ′}, the existence of which has been a matter of controversy for over a decade, while they do not contradict the interpretation of the spectral anomaly at 10–14 GPa as a precursor of the VIII-X phase transition. The calculated INS spectra correct a systematic error in the peak positions of the observed spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

High-accuracy coupled cluster calculations of atomic properties

The four-component Fock-space coupled cluster and intermediate Hamiltonian methods are implemented to evaluate atomic properties. The latter include the spectra of nobelium and lawrencium (elements 102 and 103) in the range 20000-30000 cm{sup −1}, the polarizabilities of elements 112-114 and 118, required for estimating their adsorption enthalpies on surfaces used to separate them in accelerators, and the nuclear quadrupole moments of some heavy atoms. The calculations on superheavy elements are supported by the very good agreement with experiment obtained for the lighter homologues.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2015

Two-photon absorption cross sections within equation-of-motion coupled-cluster formalism using resolution-of-the-identity and Cholesky decomposition representations: Theory, implementation, and benchmarks

The equation-of-motion coupled-cluster (EOM-CC) methods provide a robust description of electronically excited states and their properties. Here, we present a formalism for two-photon absorption (2PA) cross sections for the equation-of-motion for excitation energies CC with single and double substitutions (EOM-CC for electronically excited states with single and double substitutions) wave functions. Rather than the response theory formulation, we employ the expectation-value approach which is commonly used within EOM-CC, configuration interaction, and algebraic diagrammatic construction frameworks. In addition to canonical implementation, we also exploit resolution-of-the-identity (RI) and Cholesky decomposition (CD) for the electron-repulsion integrals to reduce memory requirements and to increase parallel efficiency. The new methods are benchmarked against the CCSD and CC3 response theories for several small molecules. We found that the expectation-value 2PA cross sections are within 5% from the quadratic response CCSD values. The RI and CD approximations lead to small errors relative to the canonical implementation (less than 4%) while affording computational savings. RI/CD successfully address the well-known issue of large basis set requirements for 2PA cross sections calculations. The capabilities of the new code are illustrated by calculations of the 2PA cross sections for model chromophores of the photoactive yellow and green fluorescent proteins.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Characterizing metastable states beyond energies and lifetimes: Dyson orbitals and transition dipole moments

The theoretical description of electronic resonances is extended beyond calculations of energies and lifetimes. We present the formalism for calculating Dyson orbitals and transition dipole moments within the equation-of-motion coupled-cluster singles and doubles method for electron-attached states augmented by a complex absorbing potential (CAP-EOM-EA-CCSD). The capabilities of the new methodology are illustrated by calculations of Dyson orbitals of various transient anions. We also present calculations of transition dipole moments between transient and stable anionic states as well as between different transient states. Dyson orbitals characterize the differences between the initial neutral and final electron-attached states without invoking the mean-field approximation. By extending the molecular-orbital description to correlated many-electron wave functions, they deliver qualitative insights into the character of resonance states. Dyson orbitals and transition moments are also needed for calculating experimental observables such as spectra and cross sections. Physically meaningful results for those quantities are obtained only in the framework of non-Hermitian quantum mechanics, e.g., in the presence of a complex absorbing potential (CAP), when studying resonances. We investigate the dependence of Dyson orbitals and transition moments on the CAP strength and illustrate how Dyson orbitals help understand the properties of metastable species and how they are affected by replacing the usual scalar product by the so-called c-product.

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

Coupled cluster valence bond theory for open-shell systems with application to very long range strong correlation in a polycarbene dimer

Here, the Coupled Cluster Valence Bond (CCVB) method, previously presented for closed-shell (CS) systems, is extended to open-shell (OS) systems. The theoretical development is based on embedding the basic OS CCVB wavefunction in a fictitious singlet super-system. This approach reveals that the OS CCVB amplitude equations are quite similar to those of CS CCVB, and thus that OS CCVB requires the same level of computational effort as CS CCVB, which is an inexpensive method. We present qualitatively correct CCVB potential energy curves for all low-lying spin states of P 2 and Mn 2 + . CCVB is successfully applied to the low-lying spin states of some model linear polycarbenes, systems that appear to be a hindrance to standard density functionals. We examine an octa-carbene dimer in a side-by-side orientation, which, in the monomer dissociation limit, exhibits maximal strong correlation over the length of the polycarbene.

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

Correlation consistent basis sets for actinides. II. The atoms Ac and Np–Lr

New correlation consistent basis sets optimized using the all-electron third-order Douglas-Kroll-Hess (DKH3) scalar relativistic Hamiltonian are reported for the actinide elements Ac and Np through Lr. These complete the series of sets reported previously for Th–U [K. A. Peterson, J. Chem. Phys. 142, 074105 (2015); M. Vasiliu et al., J. Phys. Chem. A 119, 11422 (2015)]. The new sets range in size from double- to quadruple-zeta and encompass both those optimized for valence (6s6p5f7s6d) and outer-core electron correlations (valence + 5s5p5d). The final sets have been contracted for both the DKH3 and eXact 2-component (X2C) Hamiltonians, yielding cc-pVnZ-DK3/cc-pVnZ-X2C sets for valence correlation and cc-pwCVnZ-DK3/cc-pwCVnZ-X2C sets for outer-core correlation (n = D, T, Q in each case). In order to test the effectiveness of the new basis sets, both atomic and molecular benchmark calculations have been carried out. In the first case, the first three atomic ionization potentials (IPs) of all the actinide elements Ac–Lr have been calculated using the Feller-Peterson-Dixon (FPD) composite approach, primarily with the multireference configuration interaction (MRCI) method. Excellent convergence towards the respective complete basis set (CBS) limits is achieved with the new sets, leading to good agreement with experiment, where these exist, after accurately accounting for spin-orbit effects using the 4-component Dirac-Hartree-Fock method. For a molecular test, the IP and atomization energy (AE) of PuO 2 have been calculated also using the FPD method but using a coupled cluster approach with spin-orbit coupling accounted for using the 4-component MRCI. The present calculations yield an IP 0 for PuO 2 of 159.8 kcal/mol, which is in excellent agreement with the experimental electron transfer bracketing value of 162 ± 3 kcal/mol. Likewise, the calculated 0 K AE of 305.6 kcal/mol is in very good agreement with the currently accepted experimental value of 303.1 ± 5 kcal/mol. The ground state of PuO 2 is predicted to be the 5 Σ$+\atop{0g}$ state.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2017

Excitation energies from particle-particle random phase approximation with accurate optimized effective potentials

The optimized effective potential (OEP) that gives accurate Kohn-Sham (KS) orbitals and orbital energies can be obtained from a given reference electron density. These OEP-KS orbitals and orbital energies are used here for calculating electronic excited states with the particle-particle random phase approximation (pp-RPA). Our calculations allow the examination of pp-RPA excitation energies with the exact KS density functional theory (DFT). Various input densities are investigated. Specifically, the excitation energies using the OEP with the electron densities from the coupled-cluster singles and doubles method display the lowest mean absolute error from the reference data for the low-lying excited states. This study probes into the theoretical limit of the pp-RPA excitation energies with the exact KS-DFT orbitals and orbital energies. We believe that higher-order correlation contributions beyond the pp-RPA bare Coulomb kernel are needed in order to achieve even higher accuracy in excitation energy calculations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree–Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Upadhyay, Shiv (ORCID:0000000285010501)↗
DOE OSTI2017

Spin-projected generalized Hartree-Fock method as a polynomial of particle-hole excitations

The past several years have seen renewed interest in the use of symmetry-projected Hartree-Fock for the description of strong correlations. Unfortunately, these symmetry-projected mean-field methods do not adequately account for dynamic correlation. Presumably, this shortcoming could be addressed if one could combine symmetry-projected Hartree-Fock with a many-body method such as coupled-cluster theory, but this is by no means straightforward because the two techniques are formulated in very different ways. However, we have recently shown that the singlet S 2 -projected unrestricted Hartree-Fock wave function can in fact be written in a coupled-cluster-like wave function. That is, the spin-projected unrestricted Hartree-Fock wave function can be written as a polynomial of a double-excitation operator acting on some closed-shell reference determinant. Here, we extend this result and show that the spin-projected generalized Hartree-Fock wave function (which has both S 2 and S z projection) is likewise a polynomial of low-order excitation operators acting on a closed-shell determinant and provide a closed-form expression for the resulting polynomial coefficients. The spin projection of the generalized Hartree-Fock wave function introduces connected triple and quadruple excitations which are absent when spin-projecting an unrestricted Hartree-Fock determinant. We include a few preliminary applications of the combination of this spin-projected Hartree-Fock and coupled-cluster theory to the Hubbard Hamiltonian and comment on generalizations of the methodology. Results here are not for production level, but a similarity-transformed theory that combines the two offers the promise of being accurate for both weak and strong correlation, and may offer significant improvements in the intermediate correlation regime where neither projected Hartree-Fock nor coupled cluster is particularly accurate.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Open sd-shell nuclei from first principles

We extend the ab initio coupled-cluster effective interaction (CCEI) method to open-shell nuclei with protons and neutrons in the valence space, and compute binding energies and excited states of isotopes of neon and magnesium. We employ a nucleon-nucleon and three-nucleon interaction from chiral effective field theory evolved to a lower cutoff via a similarity renormalization group transformation. We find good agreement with experiment for binding energies and spectra, while charge radii of neon isotopes are underestimated. For the deformed nuclei 20 Ne and 24 Mg we reproduce rotational bands and electric quadrupole transitions within uncertainties estimated from an effective field theory for deformed nuclei, thereby demonstrating that collective phenomena in sd-shell nuclei emerge from complex ab initio calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Nuclear Computational Low Energy Initiative (NUCLEI)

This is the final report for University of Washington for the NUCLEI SciDAC-3. The NUCLEI -project, as defined by the scope of work, will develop, implement and run codes for large-scale computations of many topics in low-energy nuclear physics. Physics to be studied include the properties of nuclei and nuclear decays, nuclear structure and reactions, and the properties of nuclear matter. The computational techniques to be used include Quantum Monte Carlo, Configuration Interaction, Coupled Cluster, and Density Functional methods. The research program will emphasize areas of high interest to current and possible future DOE nuclear physics facilities, including ATLAS and FRIB (nuclear structure and reactions, and nuclear astrophysics), TJNAF (neutron distributions in nuclei, few body systems, and electroweak processes), NIF (thermonuclear reactions), MAJORANA and FNPB (neutrino-less double-beta decay and physics beyond the Standard Model), and LANSCE (fission studies).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗