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
Theory, implementation, and disappointing results for two-photon absorption cross sections within the doubly electron-attached equation-of-motion coupled-cluster framework
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
- Award ID(s):
- 2154482
- PAR ID:
- 10420493
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 158
- Issue:
- 5
- ISSN:
- 0021-9606
- Page Range / eLocation ID:
- 054102
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Heteroaromatic species are commonly found in complex gaseous mixtures, from tobacco smoke to petroleum and asphaltene combustion products. At high temperatures, C–H bond rupture produces various dehydro radical isomers. We have used the spin–flip formulation of equation-of-motion coupled cluster theory with single and double substitutions (EOM-SF-CCSD) to characterize the energies and wave functions of the lowest lying singlet and triplet states of the diradical (2,3), (2,4), (2,5), and (3,4) di-dehydro isomers of pyrrole, furan, and thiophene. In all cases, these diradicals are minima on the broken-symmetry ωB97X-D/cc-pVDZ potential energy surface. In most cases, the diradical geometries distort to enhance through-space or through-bond coupling in the singlet states and to avoid Coulombic or exchange repulsion in the triplet states. EOM-SF-CCSD results indicate that all diradical isomers are two-configurational, closed shell singlet states. The only exceptions to this are for (2,3) and (2,4) thiophene and (2,3) pyrrole, which each contain more than two configurations. In all cases, the leading term in the multiconfigurational diradical wave function doubly occupies the symmetric radical σ orbital, indicative of either through-space or 1,3 through-bond coupling. We utilized the nucleus-independent chemical shift (NICS) approach to qualitatively assess aromaticity and find that this property varies and may be related to the energetic splittings in these diradical isomers.more » « less
-
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
-
We report on the implementation of Dyson orbitals within the recently introduced frozen-core (fc) core–valence separated (CVS) equation-of-motion (EOM) coupled-cluster singles and doubles (CCSD) method, which enables efficient and reliable characterization of core-level states. The ionization potential (IP) variant of fc-CVS-EOM-CCSD, in which the EOM target states have one electron less than the reference, gives access to core-ionized states thus enabling modeling of X-ray photoelectron spectra (XPS) and its time-resolved variant (TR-XPS). Dyson orbitals are reduced quantities that can be interpreted as correlated states of the ejected/attached electron; they enter the expressions of various experimentally relevant quantities. In the context of photoelectron spectroscopy, Dyson orbitals can be used to estimate the strengths of photoionization transitions. We illustrate the utility of Dyson orbitals and fc-CVS-EOM-IP-CCSD by calculating XPS of the ground state of adenine and TR-XPS of the excited states of uracil.more » « less
-
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
An official website of the United States government

