skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


Title: Analytic Gradients for Equation-of-Motion Coupled Cluster with Single, Double, and Perturbative Triple Excitations
Understanding the process of molecular photoexcitation is crucial in various fields, including drug development, materials science, photovoltaics, and more. The electronic vertical excitation energy is a critical property, for example in determining the singlet–triplet gap of chromophores. However, a full understanding of excited-state processes requires additional explorations of the excited-state potential energy surface and electronic properties, which is greatly aided by the availability of analytic energy gradients. Owing to its robust high accuracy over a wide range of chemical problems, equation-of-motion coupled cluster with single and double excitations (EOM-CCSD) is a powerful method for predicting excited-state properties, and the implementation of analytic gradients of many EOM-CCSD variants (excitation energies, ionization potentials, electron attachment energies, etc.) along with numerous successful applications highlights the flexibility of the method. In specific cases where a higher level of accuracy is needed or in more complex electronic structures, the inclusion of triple excitations becomes essential, for example, in the EOM-CCSD* approach of Saeh and Stanton. In this work, we derive and implement for the first time the analytic gradients of EOMEE-CCSD*, which also provides a template for analytic gradients of related excited-state methods with perturbative triple excitations. The capabilities of analytic EOMEE-CCSD* gradients are illustrated by several representative examples.  more » « less
Award ID(s):
2143725
PAR ID:
10581540
Author(s) / Creator(s):
;
Publisher / Repository:
ACS
Date Published:
Journal Name:
Journal of Chemical Theory and Computation
ISSN:
1549-9618
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, The Journal of Chemical Physics)
    The equation-of-motion coupled-cluster singles and doubles method with double electron attachment (EOM-DEA-CCSD) is capable of computing reliable energies, wave functions, and first-order properties of excited states in diradicals and polyenes that have a significant doubly excited character with respect to the ground state, without the need for including the computationally expensive triple excitations. Here, we extend the capabilities of the EOM-DEA-CCSD method to the calculations of a multiphoton property, two-photon absorption (2PA) cross sections. Closed-form expressions for the 2PA cross sections are derived within the expectation-value approach using response wave functions. We analyze the performance of this new implementation by comparing the EOM-DEA-CCSD energies and 2PA cross sections with those computed using the CC3 quadratic response theory approach. As benchmark systems, we consider transitions to the states with doubly excited character in twisted ethene and in polyenes, for which EOM-EE-CCSD (EOM-CCSD for excitation energies) performs poorly. The EOM-DEA-CCSD 2PA cross sections are comparable with the CC3 results for twisted ethene; however, the discrepancies between the two methods are large for hexatriene. The observed trends are explained by configurational analysis of the 2PA channels. 
    more » « less
  2. ; ; (, International Journal of Mass Spectrometry)
    Modeling charge migration resulting from the coherent superposition of cation ground and excited states requires information about the potential energy surfaces of the relevant cation states. Since these states are often of the same electronic symmetry as the ground state of the cation, conventional single reference methods such as coupled cluster cannot be used for the excited states. The EOMCCSD-IP (equation of motion coupled cluster with single and double excitations and ionization) is a convenient and reliable “black-box” method that can be used for the ground and excited states of cations, yielding results of CCSD (coupled cluster with singles and double excitation) quality. Charge migration in haloacetylene cations arises from the superposition of the X and A states of HCCX+ (X = F, Cl, Br and I). The geometries, ionization potentials and vibrational frequencies have been calculated by CCSD/cc-pVTZ for neutral HCCX and the X state of HCCX+ and by EOM CCSD-IP/cc-pVTZ for the X and A states of HCCX+. The results agree very well with each other and with experiment. The very good agreement between CCSD and EOMCCSD-IP for the X states demonstrates that EOMCCSD-IP is a suitable method for calculating the structure and properties of ground and excited states for the HCCX cations. 
    more » « less
  3. ; ; ; ; (, Physical Chemistry Chemical Physics)
    We present a novel approach for computing resonant inelastic X-ray scattering (RIXS) cross sections within the equation-of-motion coupled-cluster (EOM-CC) framework. The approach is based on recasting the sum-over-states expressions for RIXS moments into closed-form expressions by using damped response theory. Damped response formalism allows one to circumvent problems of divergent behavior of response equations in the resonant regime. However, the convergence of response equations in the X-ray frequency range is often erratic due to the electronically metastable ( i.e. , resonant) nature of the virtual core-excited states embedded in the valence ionization continuum. We circumvent this problematic behavior by extending the core–valence separation (CVS) scheme, which decouples the valence-excited and core-excited configurations of the excitation manifold, into the response domain. The accuracy of the CVS-enabled damped response theory, implemented within the EOM-EE-CCSD (EOM-CC for excitation energies with single and double excitations) framework, is assessed by comparison against standard damped EOM-EE-CCSD response calculations. The capabilities of the new approach are illustrated by calculations of RIXS cross sections for benzene and benzene radical cation. 
    more » « less
  4. ; ; ; (, The Journal of Chemical Physics)
    It is demonstrated that a double hybrid density functional approximation, ωB88PTPSS, that incorporates equipartition of density functional theory and the non-local correlation, however with a meta-generalized gradient approximation correlation functional, as well as with the range-separated exchange of ωB2PLYP, provides accurate excitation energies for conventional systems, as well as correct prescription of negative singlet–triplet gaps for non-conventional systems with inverted gaps, without any necessity for parametric scaling of the same-spin and opposite-spin non-local correlation energies. Examined over “safe” excitations of the QUESTDB set, ωB88PTPSS performs quite well for open-shell systems, correctly and fairly accurately [relative to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) reference] predicts negative gaps for 50 systems with inverted singlet–triplet gaps, and is one of the leading performers for intramolecular charge-transfer excitations and achieves near-second-order approximate coupled cluster (CC2) and second-order algebraic diagrammatic construction quality for the Q1 and Q2 subsets. Subsequently, we tested ωB88PTPSS on two sets of real-life examples from recent computational chemistry literature–the low energy bands of chlorophyll a (Chl a) and a set of thermally activated delayed fluorescence (TADF) systems. For Chl a, ωB88PTPSS qualitatively and quantitatively achieves DLPNO-STEOM-CCSD-level performance and provides excellent agreement with experiment. For TADF systems, ωB88PTPSS agrees quite well with spin-component-scaled CC2 (SCS-CC2) excitation energies, as well as experimental values, for the gaps between the S1 and T1 excited states. 
    more » « less
  5. ; ; (, The Journal of Chemical Physics)
    We calculate bandgaps of 12 inorganic semiconductors and insulators composed of atoms from the first three rows of the Periodic Table using periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). Our calculations are performed with atom-centered triple-zeta basis sets and up to 64 k-points in the Brillouin zone. We analyze the convergence behavior with respect to the number of orbitals and number of k-points sampled using composite corrections and extrapolations to produce our final values. When accounting for electron–phonon corrections to experimental bandgaps, we find that EOM-CCSD has a mean signed error of −0.12 eV and a mean absolute error of 0.42 eV; the largest outliers are C (error of −0.93 eV), BP (−1.00 eV), and LiH (+0.78 eV). Surprisingly, we find that the more affordable partitioned EOM-MP2 theory performs as well as EOM-CCSD. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email