Search

Computational study of electronic resonances is still a very challenging topic, with the phenomenon of dissociative electron attachment (DEA) being one of the multiple features worth investigating. Recently, we extended the charge stabilization method from energies to properties of conceptual density functional theory and applied this to metastable anionic states of ethene and chlorinated ethene derivatives to study the DEA mechanism present in these compounds. We now present an extension to spatial functions, namely, the electronic Fukui function and the electron localization function. The results of our analysis show that extrapolated spatial functions are relevant and useful for more precise localization of the unbound electron. Furthermore, we report for the first time the combination of the electron localization function with Berlin's binding function for these challenging electronic states. This promising methodology allows for accurate predictions of when and where DEA will happen in the molecules studied and provides more insight into the process. [ABSTRACT FROM AUTHOR]

When a vacancy is created in an inner-valence orbital of a dimer of atoms or molecules, the resulting species can undergo interatomic/intermolecular Coulombic decay (ICD): the hole is filled through a relaxation process that leads to a doubly ionized cluster with two positively charged atoms or molecules. Since they are subject to electronic decay, inner-valence ionized states are not bound states but electronic resonances whose transient nature can only be described with special quantum-chemical methods. In this work, we explore the capacity of equation-of-motion coupled-cluster theory with two techniques from non-Hermitian quantum mechanics, complex basis functions and Feshbach–Fano projection with a plane wave description of the outgoing electron, to describe ICD. To this end, we compute the decay rates of several dimers: Ne2, NeAr, NeMg, and (HF)2, among which the energy of the outgoing electron varies between 0.3 and 16 eV. We observe that both methods deliver better results when the outgoing electron is fast, but the characteristic R−6 distance dependence of the ICD width is captured much better with complex basis functions. [ABSTRACT FROM AUTHOR]

We present an ab initio computational study of the Auger electron spectrum of benzene. Auger electron spectroscopy exploits the Auger–Meitner effect, and although it is established as an analytic technique, the theoretical modeling of molecular Auger spectra from first principles remains challenging. Here, we use coupled-cluster theory and equation-of-motion coupled-cluster theory combined with two approaches to describe the decaying nature of core-ionized states: (i) Feshbach–Fano resonance theory and (ii) the method of complex basis functions. The spectra computed with these two approaches are in excellent agreement with each other and also agree well with experimental Auger spectra of benzene. The Auger spectrum of benzene features two well-resolved peaks at Auger electron energies above 260 eV, which correspond to final states with two electrons removed from the 1e1g and 3e2g highest occupied molecular orbitals. At lower Auger electron energies, the spectrum is less well resolved, and the peaks comprise multiple final states of the benzene dication. In line with theoretical considerations, singlet decay channels contribute more to the total Auger intensity than the corresponding triplet decay channels. [ABSTRACT FROM AUTHOR]

The charge stabilization method has often been used before for obtaining energies of temporary anions. Herein, we combine this method for the first time with conceptual density functional theory (DFT) and quantum theory of atoms in molecules by extending it to the study of nuclear Fukui functions, atom-condensed electronic Fukui functions, and bond critical points. This is applied to temporary anions of ethene and chlorinated ethene compounds, which are known to undergo dissociative electron attachment (DEA). It appears that the method is able to detect multiple valence resonance states in the same molecule, namely, a Π and a Σ state. The obtained nuclear and atom-condensed electronic Fukui functions are interpreted as nuclear forces and electron distributions, respectively, and show clear differences between the Π and Σ states. This enables a more profound characterization and understanding of how the DEA process proceeds. The conclusions are in line with findings from earlier publications, proving that the combination of conceptual DFT with the charge stabilization method yields reasonable results at rather low computational cost. [ABSTRACT FROM AUTHOR]

The emission of an Auger electron is the predominant relaxation mechanism of core-vacant states in molecules composed of light nuclei. In this non-radiative decay process, one valence electron fills the core vacancy, while a second valence electron is emitted into the ionization continuum. Because of this coupling to the continuum, core-vacant states represent electronic resonances that can be tackled with standard quantum-chemical methods only if they are approximated as bound states, meaning that Auger decay is neglected. Here, we present an approach to compute Auger decay rates of core-vacant states from coupled-cluster and equation-of-motion coupled-cluster wave functions combined with complex scaling of the Hamiltonian or, alternatively, complex-scaled basis functions. Through energy decomposition analysis, we illustrate how complex-scaled methods are capable of describing the coupling to the ionization continuum without the need to model the wave function of the Auger electron explicitly. In addition, we introduce in this work several approaches for the determination of partial decay widths and Auger branching ratios from complex-scaled coupled-cluster wave functions. We demonstrate the capabilities of our new approach by computations on core-ionized states of neon, water, dinitrogen, and benzene. Coupled-cluster and equation-of-motion coupled-cluster theory in the singles and doubles approximation both deliver excellent results for total decay widths, whereas we find partial widths more straightforward to evaluate with the former method. [ABSTRACT FROM AUTHOR]

An up-to-date overview of the CFOUR program system is given. After providing a brief outline of the evolution of the program since its inception in 1989, a comprehensive presentation is given of its well-known capabilities for high-level coupled-cluster theory and its application to molecular properties. Subsequent to this generally well-known background information, much of the remaining content focuses on lesser-known capabilities of CFOUR, most of which have become available to the public only recently or will become available in the near future. Each of these new features is illustrated by a representative example, with additional discussion targeted to educating users as to classes of applications that are now enabled by these capabilities. Finally, some speculation about future directions is given, and the mode of distribution and support for CFOUR are outlined. [ABSTRACT FROM AUTHOR]

We study the performance of the resolution-of-the-identity (RI) approximation for complex basis functions that we recently introduced [M. Hernández Vera and T.-C. Jagau, J. Chem. Phys. 151, 111101 (2019)] for second-order Møller–Plesset (MP2) perturbation theory as well as for the Coulomb and exchange contributions in Hartree–Fock theory. The sensitivity of this new RI-MP2 method toward the basis set and the auxiliary basis set is investigated, and computation times are analyzed. We show that the auxiliary basis set can be chosen purely real, that is, no complex-scaled functions need to be included. This approximation enables a further speedup of the method without compromising accuracy. We illustrate the application range of our implementation by computing static-field ionization rates of several polyacenes up to pentacene (C22H18) at the RI-MP2 level of theory. Pronounced anisotropies are observed for the ionization rates of these molecules. [ABSTRACT FROM AUTHOR]

A generalization of the Schwarz bound employed to reduce the scaling of quantum-chemical calculations is introduced in the context of non-Hermitian methods employing complex-scaled basis functions. Non-Hermitian methods offer a treatment of molecular metastable states in terms of L2-integrable wave functions with complex energies, but until now, an efficient upper bound for the resulting electron-repulsion integrals has been unavailable due to the complications from non-Hermiticity. Our newly formulated bound allows us to inexpensively and rigorously estimate the sparsity in the complex-scaled two-electron integral tensor, providing the basis for efficient integral screening procedures. We have incorporated a screening algorithm based on the new Schwarz bound into the state-of-the-art complex basis function integral code by White, Head-Gordon, and McCurdy [J. Chem. Phys. 142, 054103 (2015)]. The effectiveness of the screening is demonstrated through non-Hermitian Hartree-Fock calculations of the static field ionization of the 2-pyridoxine 2-aminopyridine molecular complex. [ABSTRACT FROM AUTHOR]

A resolution-of-the-identity (RI) approximation for two-electron integrals over Gaussian basis functions with a complex-scaled exponent is presented. Such functions are used in non-Hermitian quantum mechanics to represent electronic resonances by L2 integrable wave functions with complex energies. We have implemented this new RI approximation for second-order Møller-Plesset perturbation (MP2) theory as well as for the Coulomb and exchange contributions in Hartree-Fock (HF) theory. We discuss the differences to the standard RI approximation of Hermitian quantum mechanics and demonstrate the utility of the non-Hermitian RI-MP2 and RI-HF methods by computations of the orientation-dependent ionization rates of CO, C6H6, and C10H8 in static electric fields. Our results illustrate that RI-MP2 correctly describes correlation effects in molecular electronic resonances while the computational cost is low enough to allow for investigations of medium-sized molecules. [ABSTRACT FROM AUTHOR]

Photoelectron spectra of 63CuF2− are reported at wavelengths 310 nm, 346.6 nm, and 350.1 nm, obtained via velocity map imaging. The photoelectron angular distributions allow for the unambiguous assignment of a 2Σg+ neutral CuF2 ground state. Vibrational analysis of the direct detachment transitions in the spectra enables accurate determination of the anion and neutral bond length difference (0.073 Å), adiabatic electron affinity of CuF2 (3.494 eV) and symmetric stretching (500 cm−1, anion, and 630 cm−1, neutral) and antisymmetric stretching (610 cm−1, anion, and 782 cm−1 neutral) frequencies of the ground electronic states. Strongly photon energy dependent intensities are also observed for select transitions. Equation-of-motion coupled-cluster singles and doubles calculations augmented by a complex absorbing potential reveal a metastable 1Πg anion state which is optically accessible due to Renner-Teller coupling. Mediation of the detachment process by this state allows measurement of the bending frequencies (177 cm−1, anion, and 200 cm−1, neutral) completing the inventory of experimentally measured vibrational properties of the ground electronic states. [ABSTRACT FROM AUTHOR]

Ionization rates and Stark shifts of H2, CO, O2, H2O, and CH4 in static electric fields have been computed with coupled-cluster methods in a basis set of atom-centered Gaussian functions with a complex-scaled exponent. Consideration of electron correlation is found to be of great importance even for a qualitatively correct description of the dependence of ionization rates and Stark shifts on the strength and orientation of the external field. The analysis of the second moments of the molecular charge distribution suggests a simple criterion for distinguishing tunnel and barrier suppression ionization in polyatomic molecules. [ABSTRACT FROM AUTHOR]

The impact of residual electron correlation beyond the equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) approximation on positions and widths of electronic resonances is investigated. To establish a method that accomplishes this task in an economical manner, several approaches proposed for the approximate treatment of triple excitations are reviewed with respect to their performance in the electron attachment (EA) variant of EOM-CC theory. The recently introduced EOM-CCSD(T)(a)* method [D. A. Matthews and J. F. Stanton, J. Chem. Phys. 145, 124102 (2016)], which includes non-iterative corrections to the reference and the target states, reliably reproduces vertical attachment energies from EOM-EA-CC calculations with single, double, and full triple excitations in contrast to schemes in which non-iterative corrections are applied only to the target states. Applications of EOM-EA-CCSD(T)(a)* augmented by a complex absorbing potential (CAP) to several temporary anions illustrate that shape resonances are well described by EOM-EA-CCSD, but that residual electron correlation often makes a non-negligible impact on their positions and widths. The positions of Feshbach resonances, on the other hand, are significantly improved when going from CAP-EOM-EA-CCSD to CAP-EOM-EA-CCSD(T)(a)*, but the correct energetic order of the relevant electronic states is still not achieved. [ABSTRACT FROM AUTHOR]

The general theory of analytic energy gradients is presented for the complex absorbing potential equation-of-motion coupled-cluster (CAP-EOM-CC) method together with an implementation within the singles and doubles approximation. Expressions for theCAP-EOM-CCenergy gradient are derived based on a Lagrangian formalism with a special focus on the extra terms arising from the presence of the CAP. Our implementation allows for locating minima on high-dimensional complex-valued potential energy surfaces and thus enables geometry optimizations of resonance states of polyatomic molecules. The applicability of our CAP-EOM-CC gradients is illustrated by computations of the equilibrium structures and adiabatic electron affinities of the temporary anions of formaldehyde, formic acid, and ethylene. The results are compared to those obtained from standard EOM-CC calculations and the advantages of CAP methods are emphasized. [ABSTRACT FROM AUTHOR]

The theory and implementation of the complex-scaled coupled-cluster method with singles and doubles excitations (cs-CCSD) for studying resonances induced by static electric fields are presented. Within this framework, Stark shifts and ionization rates are obtained directly from the real and imaginary parts of the complex energy. The method is applied to the ground states of hydrogen, helium, lithium, beryllium, neon, argon, and carbon at varying field strengths. Complex-scaled Hartree-Fock, second-order many-body perturbation theory, and CCSD results are reported and analyzed with a focus on the impact of electron correlation on the ionization process. cs-CCSD calculations with suitably augmented standard Gaussian basis sets are found to deliver accurate strong-field ionization rates over a range of six orders of magnitude. The field-induced resonances are characterized beyond energy and ionization rate through their dipole moments, second moments, as well as Dyson orbitals and comparisons are drawn to autoionizing and autodetaching resonances. Marked differences are found between the tunneling and above-barrier regimes allowing for a clear distinction of the two mechanisms. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

A production-level implementation of equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) for electron attachment and excitation energies augmented by a complex absorbing potential (CAP) is presented. The new method enables the treatment of metastable states within the EOM-CC formalism in a similar manner as bound states. The numeric performance of the method and the sensitivity of resonance positions and lifetimes to the CAP parameters and the choice of oneelectron basis set are investigated. A protocol for studying molecular shape resonances based on the use of standard basis sets and a universal criterion for choosing the CAP parameters are presented. Our results for a variety of π* shape resonances of small to medium-size molecules demonstrate that CAP-augmented EOM-CCSD is competitive relative to other theoretical approaches for the treatment of resonances and is often able to reproduce experimental results. [ABSTRACT FROM AUTHOR]

The general theory required for the calculation of analytic third energy derivatives at the coupled-cluster level of theory is presented and connected to preceding special formulations for hyperpolarizabilities and polarizability gradients. Based on our theory, we have implemented a scheme for calculating the dipole Hessian matrix in a fully analytical manner within the coupled-cluster singles and doubles approximation. The dipole Hessian matrix is the second geometrical derivative of the dipole moment and thus a third derivative of the energy. It plays a crucial role in IR spectroscopy when taking into account anharmonic effects and is also essential for computing vibrational corrections to dipole moments. The superior accuracy of the analytic evaluation of third energy derivatives as compared to numerical differentiation schemes is demonstrated in some pilot calculations. [ABSTRACT FROM AUTHOR]

Analytic gradients for the state-specific multireference coupled-cluster method suggested by Mahapatra et al. [Mol. Phys. 94, 157 (1998)] (Mk-MRCC) are reported within the singles and doubles approximation using two-configurational self-consistent field (TCSCF) orbitals. The present implementation extends our previous work on Mk-MRCC gradients [E. Prochnow et al., J. Chem. Phys. 131, 064109 (2009)] which is based on restricted Hartree–Fock orbitals and consequently the main focus of the present paper is on the treatment of orbital relaxation at the TCSCF level using coupled-perturbed TCSCF theory. Geometry optimizations on m-arynes and nitrenes are presented to illustrate the influence of the orbitals on the computed equilibrium structures. The results are compared to those obtained at the single-reference coupled-cluster singles and doubles and at the Mk-MRCC singles and doubles level of theory when using restricted Hartree–Fock orbitals. [ABSTRACT FROM AUTHOR]

The formalism of response theory is applied to derive expressions for static and dynamic polarizabilities within the state-specific multireference coupled-cluster theory suggested by Mukherjee and co-workers (Mk-MRCC) [J. Chem. Phys. 110, 6171 (1998)]. We show that the redundancy problem inherent to Mk-MRCC theory gives rise to spurious poles in the Mk-MRCC response functions, which hampers the reliable calculation of dynamic polarizabilities. Furthermore, we demonstrate that in the case of a symmetry-breaking perturbation a working response theory is obtained only if certain internal excitations are included in the responses of the cluster amplitudes. Exemplary calculations within the singles and doubles approximation (Mk-MRCCSD) are carried out on aryne compounds to illustrate the impact of a multireference ansatz on the polarizability. [ABSTRACT FROM AUTHOR]

The recently presented linear-response function for Mukherjee's multireference coupled-cluster method (Mk-MRCC) [T.-C. Jagau and J. Gauss, J. Chem. Phys. 137, 044115 (2012)] is employed to determine vertical excitation energies within the singles and doubles approximation (Mk-MRCCSD-LR) for ozone as well as for o-benzyne, m-benzyne, and p-benzyne, which display increasing multireference character in their ground states. In order to assess the impact of a multireference ground-state wavefunction on excitation energies, we compare all our results to those obtained at the single-reference coupled-cluster level of theory within the singles and doubles as well as within the singles, doubles, and triples approximation. Special attention is paid to the artificial splitting of certain excited states which arises from the redundancy intrinsic to Mk-MRCC theory and hinders the straightforward application of the Mk-MRCC-LR method. [ABSTRACT FROM AUTHOR]

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.