Second-order Screened Exchange Research Articles

Assessment of the Second-Order Statically Screened Exchange Correction to the Random Phase Approximation for Correlation Energies.

With increasing interelectronic distance, the screening of the electron–electron interaction by the presence of other electrons becomes the dominant source of electron correlation. This effect is described by the random phase approximation (RPA) which is therefore a promising method for the calculation of weak interactions. The success of the RPA relies on the cancellation of errors, which can be traced back to the violation of the crossing symmetry of the 4-point vertex, leading to strongly overestimated total correlation energies. By the addition of second-order screened exchange (SOSEX) to the correlation energy, this issue is substantially reduced. In the adiabatic connection (AC) SOSEX formalism, one of the two electron–electron interaction lines in the second-order exchange term is dynamically screened (SOSEX(W, vc)). A related SOSEX expression in which both electron–electron interaction lines are statically screened (SOSEX(W(0), W(0))) is obtained from the G3W2 contribution to the electronic self-energy. In contrast to SOSEX(W, vc), the evaluation of this correlation energy expression does not require an expensive numerical frequency integration and is therefore advantageous from a computational perspective. We compare the accuracy of the statically screened variant to RPA and RPA+SOSEX(W, vc) for a wide range of chemical reactions. While both methods fail for barrier heights, SOSEX(W(0), W(0)) agrees very well with SOSEX(W, vc) for charged excitations and noncovalent interactions where they lead to major improvements over RPA.

Staggered Mesh Method for Correlation Energy Calculations of Solids: Random Phase Approximation in Direct Ring Coupled Cluster Doubles and Adiabatic Connection Formalisms

We propose a staggered mesh method for correlation energy calculations of periodic systems under the random phase approximation (RPA), which generalizes the recently developed staggered mesh method for periodic second order Møller-Plesset perturbation theory (MP2) calculations [Xing; Li; Lin J. Chem. Theory Comput. 2021]. Compared to standard RPA calculations, the staggered mesh method introduces negligible additional computational cost. It avoids a significant portion of the finite-size error and can be asymptotically advantageous for quasi-1D systems and certain quasi-2D and 3D systems with high symmetries. We demonstrate the applicability of the method using two different formalisms: the direct ring coupled cluster doubles (drCCD) theory, and the adiabatic-connection (AC) fluctuation-dissipation theory. In the drCCD formalism, the second order screened exchange (SOSEX) correction can also be readily obtained using the staggered mesh method. In the AC formalism, the staggered mesh method naturally avoids the need of performing "head/wing" corrections to the dielectric operator. The effectiveness of the staggered mesh method for insulating systems is theoretically justified by investigating the finite-size error of each individual perturbative term in the RPA correlation energy, expanded as an infinite series of terms associated with ring diagrams. As a side contribution, our analysis provides proof that the finite-size error of each perturbative term of standard RPA and SOSEX calculations scales as , where Nk is the number of grid points in a Monkhorst-Pack mesh.

The GW Miracle in Many-Body Perturbation Theory for the Ionization Potential of Molecules.

We use the GW100 benchmark set to systematically judge the quality of several perturbation theories against high-level quantum chemistry methods. First of all, we revisit the reference CCSD(T) ionization potentials for this popular benchmark set and establish a revised set of CCSD(T) results. Then, for all of these 100 molecules, we calculate the HOMO energy within second and third-order perturbation theory (PT2 and PT3), and, GW as post-Hartree-Fock methods. We found GW to be the most accurate of these three approximations for the ionization potential, by far. Going beyond GW by adding more diagrams is a tedious and dangerous activity: We tried to complement GW with second-order exchange (SOX), with second-order screened exchange (SOSEX), with interacting electron-hole pairs (W TDHF), and with a GW density-matrix (γ GW ). Only the γ GW result has a positive impact. Finally using an improved hybrid functional for the non-interacting Green’s function, considering it as a cheap way to approximate self-consistency, the accuracy of the simplest GW approximation improves even more. We conclude that GW is a miracle: Its subtle balance makes GW both accurate and fast.

Reducing the Many-Electron Self-Interaction Error in the Second-Order Screened Exchange Method.

Preserving the beneficial properties of the second-order screened exchange (SOSEX) method, such as its freedom from one-body self-correlation error, and its seamless treatment of long-range dispersion interaction, we construct new nonlocal functionals by down-scaling the higher-order terms in the SOSEX series to reduce the many-body self-correlation error in molecular systems. Our down-scaled SOSEX (dsSOSEX) and scaled equation SOSEX (seSOSEX) approaches deliver considerably more accurate noncovalent interaction energies, reaction energies, and barrier heights than the original SOSEX method. The dsSOSEX approach improves on the description of medium- and long-range correlation, for example, in isogyric processes, while the seSOSEX approach improves on the description of processes with short- and medium-range rearrangement in the electronic structure, such as atomization. The He2+ potential energy curve shows that the many-body self-correlation error is slightly smaller in the dsSOSEX method than in the seSOSEX approach. Furthermore, the H2 potential energy curve reveals that the static correlation error can also be reduced at medium range by the dsSOSEX and seSOSEX schemes. As our dsSOSEX method provides a better description of medium- and long-range correlation, it generally leads to more accurate energetics than the seSOSEX method; hence it should be preferred over the latter. We compare our results to other empirical and nonempirical direct random phase approximation (dRPA) and beyond-dRPA methods and also discuss the possibilities for further development.

Description of noncovalent interactions involving π-system with high precision: An assessment of RPA, MP2, and DFT-D methods.

Efficient approaches with high precision are essential for understanding the formation and stability of noncovalent interaction complexes. Here, 21 noncovalent interaction complexes involving π-system are selected and grouped in three subsets according to ETS-NOCV method: dispersion-dominated, electrostatic-dominated, and mixed. We mainly focus on examining the performance of random-phase approximation (RPA) on these π systems. The tested RPA-based method includes standard RPA and its variants including the related single excitations (SEs), renormalized single excitations (rSEs), second-order screened exchange (SOSEX), and the renormalized second-order perturbation theory (rPT2). The routine second-order Møller-Plesset perturbation theory (MP2) and three popular DFT-D functionals (M06-2X-D3, ωB97XD, and PBE-D3(BJ)) are also assessed for comparison. In this work, besides the calculation of interaction energies at Dunning-type aug-cc-pVDZ and aug-cc-pVTZ basis set, we also present a larger database of interaction energies calculated using MP2 and RPA methods with Dunning-type aug-cc-pVQZ basis set. An accurate CCSD(T)/CBS scheme is used to provide benchmark database. In addition to the high-level results, we also provide potential energy surfaces (PES) of different interaction type. Among all the tested methods, MP2 has a satisfactory performance on electrostatic-dominated and mixed-type systems, except for dispersion-dominated systems. DFT-D functionals, especially ωB97XD functional, has a balanced performance across all the tested systems. Importantly, for RPA-based methods, the calculation accuracy can be dramatically improved by taking into account SE or exchange effects, especially in the mixed complexes. We conclude that rPT2 among all the test RPA-based methods gives an overall satisfactory performance across different interaction types. © 2019 Wiley Periodicals, Inc.

Screened Exchange Corrections to the Random Phase Approximation from Many-Body Perturbation Theory.

The random phase approximation (RPA) systematically overestimates the magnitude of the correlation energy and generally underestimates cohesive energies. This originates in part from the complete lack of exchange terms that would otherwise cancel Pauli exclusion principle violating (EPV) contributions. The uncanceled EPV contributions also manifest themselves in form of an unphysical negative pair density of spin-parallel electrons close to electron-electron coalescence. We follow considerations of many-body perturbation theory to propose an exchange correction that corrects the largest set of EPV contributions, while having the lowest possible computational complexity. The proposed method exchanges adjacent particle/hole pairs in the RPA diagrams, considerably improving the pair density of spin-parallel electrons close to coalescence in the uniform electron gas (UEG). The accuracy of the correlation energy is comparable to other variants of second-order screened exchange (SOSEX) corrections although it is slightly more accurate for the spin-polarized UEG. Its computational complexity scales as [Formula: see text] or [Formula: see text] in orbital space or real space, respectively. Its memory requirement scales as [Formula: see text].

Performance and Scope of Perturbative Corrections to Random-Phase Approximation Energies

It has been suspected since the early days of the random-phase approximation (RPA) that corrections to RPA correlation energies result mostly from short-range correlation effects and are thus amenable to perturbation theory. Here we test this hypothesis by analyzing formal and numerical results for the most common beyond-RPA perturbative corrections, including the bare second-order exchange (SOX), second-order screened exchange (SOSEX), and approximate exchange kernel (AXK) methods. Our analysis is facilitated by efficient and robust algorithms based on the resolution-of-the-identity (RI) approximation and numerical frequency integration, which enable benchmark beyond-RPA calculations on medium- and large-size molecules with size-independent accuracy. The AXK method systematically improves upon RPA, SOX, and SOSEX for reaction barrier heights, reaction energies, and noncovalent interaction energies of main-group compounds. The improved accuracy of AXK compared with SOX and SOSEX is attributed to stronger screening of bare SOX in AXK. For reactions involving transition-metal compounds, particularly 3d transition-metal dimers, the AXK correction is too small and can even have the wrong sign. These observations are rationalized by a measure α̅ of the effective coupling strength for beyond-RPA correlation. When the effective coupling strength increases beyond a critical α̅ value of approximately 0.5, the RPA errors increase rapidly and perturbative corrections become unreliable. Thus, perturbation theory can systematically correct RPA but only for systems and properties qualitatively well captured by RPA, as indicated by small α̅ values.

Beyond the random phase approximation with a local exchange vertex

With the aim of constructing an electronic structure approach that systematically goes beyond the GW and random phase approximation (RPA) we introduce a vertex correction based on the exact-exchange (EXX) potential of time-dependent density functional theory. The EXX vertex function is constrained to be local but is expected to capture similar physics as the Hartree-Fock vertex. With the EXX vertex we then consistently unify different beyond-RPA approaches such as the various re-summations of RPA with exchange (RPAx) and the second order screened exchange (SOSEX) approximation. The theoretical analysis is supported by numerical studies on the hydrogen dimer and the electron gas, and we discuss the importance of including the vertex correction in both the screened interaction and the self energy. Finally, we give details on our implementation within the plane-wave pseudo potential framework and demonstrate the excellent performance of the different RPAx methods in describing the energetics of hydrogen and van der Waals bonds.

Short-range second order screened exchange correction to RPA correlation energies.

Direct random phase approximation (RPA) correlation energies have become increasingly popular as a post-Kohn-Sham correction, due to significant improvements over DFT calculations for properties such as long-range dispersion effects, which are problematic in conventional density functional theory. On the other hand, RPA still has various weaknesses, such as unsatisfactory results for non-isogyric processes. This can in parts be attributed to the self-correlation present in RPA correlation energies, leading to significant self-interaction errors. Therefore a variety of schemes have been devised to include exchange in the calculation of RPA correlation energies in order to correct this shortcoming. One of the most popular RPA plus exchange schemes is the second order screened exchange (SOSEX) correction. RPA + SOSEX delivers more accurate absolute correlation energies and also improves upon RPA for non-isogyric processes. On the other hand, RPA + SOSEX barrier heights are worse than those obtained from plain RPA calculations. To combine the benefits of RPA correlation energies and the SOSEX correction, we introduce a short-range RPA + SOSEX correction. Proof of concept calculations and benchmarks showing the advantages of our method are presented.

Random phase approximation with second-order screened exchange for current-carrying atomic states.

The direct random phase approximation (RPA) and RPA with second-order screened exchange (SOSEX) have been implemented with complex orbitals as a basis for treating open-shell atoms. Both RPA and RPA+SOSEX are natural implicit current density functionals because the paramagnetic current density implicitly is included through the use of complex orbitals. We confirm that inclusion of the SOSEX correction improves the total energy accuracy substantially compared to RPA, especially for smaller-Z atoms. Computational complexity makes post self-consistent-field (post-SCF) evaluation of RPA-type expressions commonplace, so orbital basis origins and properties become important. Sizable differences are found in correlation energies, total atomic energies, and ionization energies for RPA-type functionals evaluated in the post-SCF fashion with orbital sets obtained from different schemes. Reference orbitals from Kohn-Sham calculations with semi-local functionals are more suitable for RPA+SOSEX to generate accurate total energies, but reference orbitals from exact exchange (non-local) yield essentially energetically degenerate open-shell atom ground states. RPA+SOSEX correlation combined with exact exchange calculated from a hybrid reference orbital set (half the exchange calculated from exact-exchange orbitals, the other half of the exchange from orbitals optimized for the Perdew-Burke-Ernzerhof (PBE) exchange functional) gives the best overall performance. Numerical results show that the RPA-like functional with SOSEX correction can be used as a practical implicit current density functional when current effects should be included.

Dielectric Matrix Formulation of Correlation Energies in the Random Phase Approximation: Inclusion of Exchange Effects.

Starting from the general expression for the ground state correlation energy in the adiabatic-connection fluctuation-dissipation theorem (ACFDT) framework, it is shown that the dielectric matrix formulation, which is usually applied to calculate the direct random phase approximation (dRPA) correlation energy, can be used for alternative RPA expressions including exchange effects. Within this famework, the ACFDT analog of the second order screened exchange (SOSEX) approximation leads to a logarithmic formula for the correlation energy similar to the direct RPA expression. Alternatively, the contribution of the exchange can be included in the kernel used to evaluate the response functions. In this case, the use of an approximate kernel is crucial to simplify the formalism and to obtain a correlation energy in logarithmic form. Technical details of the implementation of these methods are discussed, and it is shown that one can take advantage of density fitting or Cholesky decomposition techniques to improve the computational efficiency; a discussion on the numerical quadrature made on the frequency variable is also provided. A series of test calculations on atomic correlation energies and molecular reaction energies shows that exchange effects are instrumental for improvement over direct RPA results.

Accurate Ionization Potentials and Electron Affinities of Acceptor Molecules III: A Benchmark of GW Methods

The performance of different GW methods is assessed for a set of 24 organic acceptors. Errors are evaluated with respect to coupled cluster singles, doubles, and perturbative triples [CCSD(T)] reference data for the vertical ionization potentials (IPs) and electron affinities (EAs), extrapolated to the complete basis set limit. Additional comparisons are made to experimental data, where available. We consider fully self-consistent GW (scGW), partial self-consistency in the Green's function (scGW0), non-self-consistent G0W0 based on several mean-field starting points, and a "beyond GW" second-order screened exchange (SOSEX) correction to G0W0. We also describe the implementation of the self-consistent Coulomb hole with screened exchange method (COHSEX), which serves as one of the mean-field starting points. The best performers overall are G0W0+SOSEX and G0W0 based on an IP-tuned long-range corrected hybrid functional with the former being more accurate for EAs and the latter for IPs. Both provide a balanced treatment of localized vs delocalized states and valence spectra in good agreement with photoemission spectroscopy (PES) experiments.

Accurate, precise, and efficient theoretical methods to calculate anion-π interaction energies in model structures.

A correct description of the anion-π interaction is essential for the design of selective anion receptors and channels and important for advances in the field of supramolecular chemistry. However, it is challenging to do accurate, precise, and efficient calculations of this interaction, which are lacking in the literature. In this article, by testing sets of 20 binary anion-π complexes of fluoride, chloride, bromide, nitrate, or carbonate ions with hexafluorobenzene, 1,3,5-trifluorobenzene, 2,4,6-trifluoro-1,3,5-triazine, or 1,3,5-triazine and 30 ternary π-anion-π' sandwich complexes composed from the same monomers, we suggest domain-based local-pair natural orbital coupled cluster energies extrapolated to the complete basis-set limit as reference values. We give a detailed explanation of the origin of anion-π interactions, using the permanent quadrupole moments, static dipole polarizabilities, and electrostatic potential maps. We use symmetry-adapted perturbation theory (SAPT) to calculate the components of the anion-π interaction energies. We examine the performance of the direct random phase approximation (dRPA), the second-order screened exchange (SOSEX), local-pair natural-orbital (LPNO) coupled electron pair approximation (CEPA), and several dispersion-corrected density functionals (including generalized gradient approximation (GGA), meta-GGA, and double hybrid density functional). The LPNO-CEPA/1 results show the best agreement with the reference results. The dRPA method is only slightly less accurate and precise than the LPNO-CEPA/1, but it is considerably more efficient (6-17 times faster) for the binary complexes studied in this paper. For 30 ternary π-anion-π' sandwich complexes, we give dRPA interaction energies as reference values. The double hybrid functionals are much more efficient but less accurate and precise than dRPA. The dispersion-corrected double hybrid PWPB95-D3(BJ) and B2PLYP-D3(BJ) functionals perform better than the GGA and meta-GGA functionals for the present test set.

Renormalized second-order perturbation theory for the electron correlation energy: Concept, implementation, and benchmarks

We present a renormalized second-order perturbation theory (rPT2), based on a Kohn-Sham (KS) reference state, for the electron correlation energy that includes the random-phase approximation (RPA), second-order screened exchange (SOSEX), and renormalized single excitations (rSE). These three terms all involve a summation of certain types of diagrams to infinite order, and can be viewed as "renormalization" of the 2nd-order direct, exchange, and single excitation (SE) terms of Rayleigh-Schr\"odinger perturbation theory based on an KS reference. In this work we establish the concept of rPT2 and present the numerical details of our SOSEX and rSE implementations. A preliminary version of rPT2, in which the renormalized SE (rSE) contribution was treated approximately, has already been benchmarked for molecular atomization energies and chemical reaction barrier heights and shows a well balanced performance [Paier et al, New J. Phys. 14, 043002 (2012)]. In this work, we present a refined version of rPT2, in which we evaluate the rSE series of diagrams rigorously. We then extend the benchmark studies to non-covalent interactions, including the rare-gas dimers, and the S22 and S66 test sets. Despite some remaining shortcomings, we conclude that rPT2 gives an overall satisfactory performance across different chemical environments, and is a promising step towards a generally applicable electronic structure approach.

Assessment of correlation energies based on the random-phase approximation

The random-phase approximation to the ground state correlation energy (RPA) in combination with exact exchange (EX) has brought the Kohn–Sham (KS) density functional theory one step closer towards a universal, ‘general purpose first-principles method’. In an effort to systematically assess the influence of several correlation energy contributions beyond RPA, this paper presents dissociation energies of small molecules and solids, activation energies for hydrogen transfer and non-hydrogen transfer reactions, as well as reaction energies for a number of common test sets. We benchmark EX + RPA and several flavors of energy functionals going beyond it: second-order screened exchange (SOSEX), single-excitation (SE) corrections, renormalized single-excitation (rSE) corrections and their combinations. Both the SE correction and the SOSEX contribution to the correlation energy significantly improve on the notorious tendency of EX + RPA to underbind. Surprisingly, activation energies obtained using EX + RPA based on a KS reference alone are remarkably accurate. RPA + SOSEX + rSE provides an equal level of accuracy for reaction as well as activation energies and overall gives the most balanced performance, because of which it can be applied to a wide range of systems and chemical reactions.

Correlation Energy Expressions from the Adiabatic-Connection Fluctuation-Dissipation Theorem Approach.

We explore several random phase approximation (RPA) correlation energy variants within the adiabatic-connection fluctuation-dissipation theorem approach. These variants differ in the way the exchange interactions are treated. One of these variants, named dRPA-II, is original to this work and closely resembles the second-order screened exchange (SOSEX) method. We discuss and clarify the connections among different RPA formulations. We derive the spin-adapted forms of all the variants for closed-shell systems and test them on a few atomic and molecular systems with and without range separation of the electron-electron interaction.