Charge-transfer separability and size-extensivity in the equation-of-motion coupled cluster method: EOM-CCx

Perspective on Coupled-cluster Theory. The evolution toward simplicity in quantum chemistry

Coupled-cluster theory has revolutionized quantum chemistry. It has provided the framework to effectively solve the problem of electron correlation, the main focus of the field for over 60 years. This has enabled ab initio quantum chemistry to provide predictive quality results for most quantities of interest that are obtainable from first-principle calculations. The best that one can do in a basis is the 'full CI,' the exact solution of the non-relativistic Schrödinger equation or, if need be, the relativistic Dirac equation. With due regard to converging the basis set and adequate consideration of higher clusters and relativity in a calculation, virtually predictive results can be obtained. But in addition to its numerical performance, coupled-cluster theory also offers a conceptually new, many-body foundation for the theory that should be appreciated by all practitioners. The latter is emphasized in this perspective, leading to the 'evolution toward simplicity' in the title. The ultimate theory will benefit from the several features that are uniquely exact in coupled-cluster theory and its equation-of-motion (EOM-CC) extensions.

Benchmark Study of the Electronic States of the LiRb Molecule: Ab Initio Calculations with the Fock Space Coupled Cluster Approach

Accurate potential energy curves (PECs) are determined for the twenty-two electronic states of LiRb. In contrast to previous studies, the applied approach relies on the first principle calculations involving correlation among all electrons. The current methodology is founded on the multireference coupled cluster (CC) scheme constructed within the Fock space (FS) formalism, specifically for the (2,0) sector. The FS methodology is established within the framework of the intermediate Hamiltonian formalism and offers an intruder-free, efficient computational scheme. This method has a distinctive feature that, when applied to the doubly ionized system, provides the characteristics of the neutral case. This proves especially beneficial when investigating PECs in situations where a closed-shell molecule dissociates into open-shell fragments, yet its double positive ion forms closed-shell species. In every instance, we successfully computed continuous PECs spanning the entire range of interatomic distances, from the equilibrium to the dissociation limit. Moreover, the spectroscopic characteristic of various electronic states is presented, including relativistic effects. Relativistic corrections included at the third-order Douglas-Kroll level have a non-negligible effect on the accuracy of the determined spectroscopic constants.

Excited-State-Specific Pseudoprojected Coupled-Cluster Theory

We present an excited-state-specific coupled-cluster approach in which both the molecular orbitals and cluster amplitudes are optimized for an individual excited state. The theory is formulated via a pseudoprojection of the traditional coupled-cluster wavefunction that allows correlation effects to be introduced atop an excited-state mean field starting point. The approach shares much in common with ground-state CCSD, including size extensivity and an N6 cost scaling. Preliminary numerical tests show that, when augmented with N5 cost perturbative corrections for key terms, the method can improve over excited-state-specific second-order perturbation theory in valence, charge transfer, and Rydberg states.

Accurate Prediction of Vertical Ionization Potentials and Electron Affinities from Spin-Component Scaled CC2 and ADC(2) Models

The CC2 and ADC(2) wave function models and their spin-component scaled modifications are adopted for predicting vertical ionization potentials (VIPs) and electron affinities (VEAs). The ionic solutions are obtained as electronic excitations in the continuum orbital formalism, making possible the use of existing, widespread quantum chemistry codes with minimal modifications, in full consistency with the treatment of charge transfer excitations. The performance of different variants is evaluated via benchmark calculations on various sets from previous works, containing small- and medium-sized systems, including the nucleobases. It is shown that with the spin-scaled approximate methods, in particular the scaled opposite-spin variant of the ADC(2) method, the accuracy of EOM-CCSD is achievable at a fraction of the computational cost, also outperforming many common electron propagator approaches.

Intermediate Hamiltonian Fock-space coupled-cluster theory for excitation energies, double ionization potentials, and double electron attachments with spin–orbit coupling

The intermediate Hamiltonian Fock-space coupled-cluster methods at the singles and doubles level (IHFSCCSD) for excitation energies in the (1p, 1h) sector, double ionization potentials in the (0p, 2h) sector, and double electron attachments in the (2p, 0h) sector of the Fock space are implemented based on the CCSD method with spin–orbit coupling (SOC) included in the post-Hartree–Fock treatment using a closed-shell reference in this work. The active space is chosen to contain those orbitals that have the largest contribution to principal ionized or electron-attached states obtained from the equation-of-motion coupled-cluster calculations. Both time-reversal symmetry and spatial symmetry are exploited in the implementation. Our results show that the accuracy of IHFSCCSD results is closely related to the active space, and the sufficiency of the active space can be assessed from the percentage of transitions within the active space. In addition, unreasonable results may be encountered when the ionized or electron-attached states with a somewhat larger contribution from double excitations are included to determine the active space and cluster operators in the (0p, 1h) or (1p, 0h) sector of the Fock space. A larger active space may be required to describe SO splitting reliably than that in the scalar-relativistic calculations in some cases. The IHFSCCSD method with SOC developed in this work can provide reliable results for heavy-element systems when a sufficient active space built upon the principal ionization potential/electron affinity states is adopted.

Decay Processes in Cationic Alkali Metals in Microsolvated Clusters: A Complex Absorbing Potential Based Equation-of-Motion Coupled Cluster Investigation

We have employed the highly accurate complex absorbing potential based ionization potential equation-of-motion coupled cluster singles and doubles (CAP-IP-EOM-CCSD) method to study the various intermolecular decay processes in ionized metals (Li⁺, Na⁺, K⁺) microsolvated by water molecules. For the Li atom, the electron is ionized from the 1s subshell. However, for Na and K atoms, the electron is ionized from 2s and both 2s and 2p subshells, respectively. We have investigated decay processes for the Li⁺-(H₂O)_n (n = 1-3) systems, as well as Na⁺-(H₂O)_n (n = 1, 2), and K⁺-H₂O. The lithium cation in water can decay only via electron transfer mediated decay (ETMD) as there are no valence electrons in lithium. We have investigated how the various decay processes change in the presence of different alkali metal atoms and how the increasing number of water molecules play a significant role in the decay of microsolvated systems. To see the effect of the environment, we have studied Li⁺-NH₃ in comparison to Li⁺-H₂O. In the case of Na⁺-H₂O, we have studied the impact of bond distance on the decay width. The effect of polarization on decay width was checked for the X⁺-H₂O (X = Li, Na) systems. We used the PCM model to study the polarization effect. We have compared our results with existing theoretical and experimental results wherever available in the literature.

Equation of motion coupled-cluster study of core excitation spectra II: Beyond the dipole approximation

We present the time-independent (TI) and time-dependent (TD) equation of motion coupled-cluster (EOM-CC) oscillator strengths not limited to those obtained by the dipole approximation. For the conventional TI-EOM-CC, we implement all the terms in the multipole expansion through second order that contributes to the oscillator strength. These include contributions such as magnetic dipole, electric quadrupole, electric octupole, and magnetic quadrupole. In TD-EOM-CC, we only include the quadrupole moment contributions. This augments our previous work [Y. C. Park, A. Perera, and R. J. Bartlett, J. Chem. Phys. 151, 164117 (2019)]. The inclusion of the quadrupole contributions (and all the other contributions through second order in the case of TI-EOM-CCSD) enables us to obtain the intensities for the pre-edge transitions in the metal K-edge spectra, which are dipole inactive. The TI-EOM-CCSD and TD-EOM-CCSD spectra of Ti⁴⁺ atoms are used to showcase the implementation of the second-order oscillator strengths. The origin of 1s → e and 1s → t₂ in core spectra from iron tetrachloride and titanium tetrachloride is discussed and compared with the experiment.

An efficient Fock space multi-reference coupled cluster method based on natural orbitals: Theory, implementation, and benchmark

We present a natural orbital-based implementation of the intermediate Hamiltonian Fock space coupled-cluster method for the (1, 1) sector of Fock space. The use of natural orbitals significantly reduces the computational cost and can automatically choose an appropriate set of active orbitals. The new method retains the charge transfer separability of the original intermediate Hamiltonian Fock space coupled-cluster method and gives excellent performance for valence, Rydberg, and charge-transfer excited states. It offers significant computational advantages over the popular equation of motion coupled cluster method for excited states dominated by single excitations.

The Devil's Triangle of Kohn-Sham density functional theory and excited states

Exchange-correlation (XC) functionals from Density Functional Theory (DFT) developed under the rigorous arguments of Correlated Orbital Theory (COT) address the Devil's Triangle of prominent errors in Kohn-Sham (KS) DFT. At the foundation of this triangle lie the incorrect one-particle spectrum, the lack of integer discontinuity, and the self-interaction error. At the top level, these failures manifest themselves in incorrect charge transfer and Rydberg excitation energies, along with poor activation barriers. Accordingly, the Quantum Theory Project (QTP) XC functionals have been created to address several of the long-term issues encountered in KS theory and its Time Dependent DFT (TDDFT) variant for electronic excitations. Recognizing that COT starts with a correct one-particle spectrum, a condition imposed on the QTP functionals by means of minimum parameterization, the question that arises is how does this affect the electronically excited states? Among up to 28 XC functionals considered, the QTP family provides one of the smallest mean absolute deviations for charge-transfer excitations while also showing excellent results for Rydberg states. However, there is some room for improvement in the case of excitation energies to valence states, which are systematically underestimated by all functionals investigated. An alternative path for better treatment of excitation energies, mainly for valence states, is offered by using orbital energies from QTP functionals, especially by CAM-QTP-02 and LC-QTP. In this case, the deviations from the reference data can be reduced approximately by half.

Electronic states of NaLi molecule: Benchmark results with Fock space coupled cluster approach

Accurate potential energy curves (PECs) are obtained for 20 lowest lying electronic states of the NaLi molecule. The computational scheme used here is based on the multireference coupled cluster theory formulated in the (2,0) sector of the Fock space. The latter sector provides the description of states obtained by attachment of two electrons to the reference system. This makes it possible to adopt the doubly ionized NaLi⁺² molecule as a Fermi vacuum. The latter has a very concrete advantage in calculations of the PECs since it dissociates into closed shell fragments (NaLi⁺² → Na⁺ + Li⁺); hence, the restricted Hartree-Fock method can be used within the whole range of interatomic distances. Computed PECs and spectroscopic constants stay very close to the experimental values (if the latter are available) with the accuracy exceeding the other theoretical approaches including those based on the effective core polarization potentials. Relativistic corrections included at the infinite-order two-component level have a non-negligible effect on the accuracy of computed excitation and dissociation energies with contributions up to 50 cm^-1.

Equation-of-motion coupled-cluster theory for double electron attachment with spin-orbit coupling

We report implementation of the equation-of-motion coupled-cluster (EOM-CC) method for double electron-attachment (DEA) with spin-orbit coupling (SOC) at the CC singles and doubles (CCSD) level using a closed-shell reference in this work. The DEA operator employed in this work contains two-particle and three-particle one-hole excitations, and SOC is included in post-Hartree-Fock treatment. Time-reversal symmetry and spatial symmetry are exploited to reduce computational cost. The EOM-DEA-CCSD method with SOC allows us to investigate SOC effects of systems with two-unpaired electrons. According to our results on atoms, double ionization potentials (DIPs), excitation energies (EEs), and SO splittings of low-lying states are calculated reliably using the EOM-DEA-CCSD method with SOC. Its accuracy is usually higher than that of EOM-CCSD for EEs or DIPs if the same target can be reached from single excitations by choosing a proper closed-shell reference. However, performance of the EOM-DEA-CCSD method with SOC on molecules is not as good as that for atoms. Bond lengths for the ground and the several lowest excited states of GaH, InH, and TlH are underestimated pronouncedly, although reasonable EEs are obtained, and splittings of the ³Σ^- state from the π² configuration are calculated to be too small with EOM-DEA-CCSD.

A confinement induced spectroscopic study of noble gas atoms using equation of motion architecture: Encapsulation within fullerene's voids

A relativistic study of spectroscopic properties of the endohedral fullerenes Ng@C₆₀^q (where Ng = He, Ne and q=0,±1,±2 are the charges) associated with the C₆₀ molecule has been done using the equation of motion coupled cluster (EOM-CC) methodology. Specific properties estimated are the transition energies, dipole oscillator strengths, and transition probabilities for the low-lying excitations 1s²(¹S₀) → 1snp (¹P₁) (n = 2, 3, 4) for He@C₆₀^q and 1s²2s²2p⁶ (¹S₀) → 1s²2s²2p⁵ns∕nd (¹P₁) (n = 3, 4) for Ne@C₆₀^q, which have been compared with those for the isolated atom to depict the confinement effect of the host molecule on the encapsulated atom. This is accomplished by introducing an effective potential to the atomic Hamiltonian induced by the fullerene moiety and its charge. The EOM-CC results have been compared with those estimated with the random phase approximation (and configuration interaction singles) to understand the effect of electron correlation under such confinement. The systematic and interesting behavior of the properties is highlighted indicating the effect of fullerene cage potential on the redistribution of electron density of the guest atom.