Influence of Electronic Correlations on Electron-Phonon Interactions of Molecular Systems with the GW and Coupled Cluster Methods.

Electron-phonon interactions are of great importance to a variety of physical phenomena, and their accurate description is an important goal for first-principles calculations. Isolated examples of materials and molecular systems have emerged where electron-phonon coupling is enhanced over density functional theory (DFT) when using the Green's-function-based ab initio GW method, which provides a more accurate description of electronic correlations. It is, however, unclear how general this enhancement is and how employing high-end quantum chemistry methods, which further improve the description of electronic correlations, might further alter electron-phonon interactions over GW or DFT. Here, we address these questions by computing the renormalization of the highest occupied molecular orbital energies of Thiel's set of organic molecules by harmonic vibrations using DFT, GW, and equation-of-motion coupled-cluster calculations. We find that, depending on the amount of exact exchange included in the DFT starting point, GW can increase the magnitude of the electron-phonon coupling across Thiel's set of molecules by an average factor of 1.1-1.8 compared to the underlying DFT, while equation-of-motion coupled-cluster leads to an increase of 1.4-2. The electron-phonon coupling predicted with the ab initio GW method is generally in much closer agreement to coupled cluster values compared to DFT, establishing GW as a promising route for accurately computing electron-phonon phenomena in molecules and beyond at a much lower computational cost than higher-end quantum chemistry techniques.

Fully Dynamic G3W2 Self-Energy for Finite Systems: Formulas and Benchmark.

Over the years, Hedin's GW self-energy has been proven to be a rather accurate and simple approximation to evaluate electronic quasiparticle energies in solids and in molecules. Attempts to improve over the simple GW approximation, the so-called vertex corrections, have been constantly proposed in the literature. Here, we derive, analyze, and benchmark the complete second-order term in the screened Coulomb interaction W for finite systems. This self-energy named G3W2 contains all the possible time orderings that combine 3 Green's functions G and 2 dynamic W. We present the analytic formula and its imaginary frequency counterpart, with the latter allowing us to treat larger molecules. The accuracy of the G3W2 self-energy is evaluated on well-established benchmarks (GW100, Acceptor 24, and Core 65) for valence and core quasiparticle energies. Its link with the simpler static approximation, named SOSEX for static screened second-order exchange, is analyzed, which leads us to propose a more consistent approximation named 2SOSEX. In the end, we find that neither the G3W2 self-energy nor any of the investigated approximations to it improve over one-shot G0W0 with a good starting point. Only quasi-particle self-consistent GW HOMO energies are slightly improved by addition of the G3W2 self-energy correction. We show that this is due to the self-consistent update of the screened Coulomb interaction, leading to an overall sign change of the vertex correction to the frontier quasiparticle energies.

RPA, an Accurate and Fast Method for the Computation of Static Nonlinear Optical Properties.

The accurate computation of static nonlinear optical properties (SNLOPs) in large polymers requires accounting for electronic correlation effects with a reasonable computational cost. The Random Phase Approximation (RPA) used in the adiabatic connection fluctuation theorem is known to be a reliable and cost-effective method to render electronic correlation effects when combined with density-fitting techniques and integration over imaginary frequencies. We explore the ability of the RPA energy expression to predict SNLOPs by evaluating RPA electronic energies in the presence of finite electric fields to obtain (using the finite difference method) static polarizabilities and hyperpolarizabilities. We show that the RPA based on hybrid functional self-consistent field calculations yields accurate SNLOPs as the best-tuned double-hybrid functionals developed today, with the additional advantage that the RPA avoids any system-specific adjustment.

An optimally tuned range-separated hybrid starting point for ab initio GW plus Bethe-Salpeter equation calculations of molecules.

The ab initio GW plus Bethe-Salpeter equation (GW-BSE, where G is the one particle Green's function and W is the screened Coulomb interaction) approach has emerged as a leading method for predicting excitations in both solids and molecules with a predictive power contingent upon several factors. Among these factors are the (1) generalized Kohn-Sham eigensystem used to construct the GW self-energy and to solve the BSE and (2) the efficacy and suitability of the Tamm-Dancoff approximation. Here, we present a detailed benchmark study of low-lying singlet excitations from a generalized Kohn-Sham (gKS) starting point based on an optimally tuned range-separated hybrid (OTRSH) functional. We show that the use of this gKS starting point with one-shot G0W0 and G0W0-BSE leads to the lowest mean absolute errors (MAEs) and mean signed errors (MSEs), with respect to high-accuracy reference values, demonstrated in the literature thus far for the ionization potentials of the GW100 benchmark set and for low-lying neutral excitations of Thiel's set molecules in the gas phase, without the need for self-consistency. The MSEs and MAEs of one-shot G0W0-BSE@OTRSH excitation energies are comparable to or lower than those obtained with other functional starting points after self-consistency. Additionally, we compare these results with linear-response time-dependent density functional theory (TDDFT) calculations and find GW-BSE to be superior to TDDFT when calculations are based on the same exchange-correlation functional. This work demonstrates tuned range-separated hybrids used in combination with GW and GW-BSE can greatly suppress starting point dependence for molecules, leading to accuracy similar to that for higher-order wavefunction-based theories for molecules without the need for costlier iterations to self-consistency.

Ab Initio Prediction of a Negative Barkas Coefficient for Slow Protons and Antiprotons in LiF.

We report the ab initio prediction of a negative Barkas coefficient in lithium fluoride (LiF) insulator at low velocity (v<0.25 a.u., E_{kin}∼2 keV). The electronic stopping power of protons in LiF has been extensively studied both experimentally and theoretically because of a controversial threshold effect. While our time-dependent density-functional theory simulations confirm the presence of a velocity threshold below which the proton stopping power vanishes, our calculations demonstrate that the antiprotons do not experience such a threshold. The combination of those two contrasting behaviors gives rise to an unprecedented negative Barkas effect: the stopping power of antiprotons is larger than that of protons. We identify that the slow antiproton at close encounter destabilizes a p orbital of the F^{-} anion pointing toward the antiproton. This particular orbital becomes highly polarizable and hence contributes much to the stopping power.

Coupling Natural Orbital Functional Theory and Many-Body Perturbation Theory by Using Nondynamically Correlated Canonical Orbitals.

We develop a new family of electronic structure methods for capturing at the same time the dynamic and nondynamic correlation effects. We combine the natural orbital functional theory (NOFT) and many-body perturbation theory (MBPT) through a canonicalization procedure applied to the natural orbitals to gain access to any MBPT approximation. We study three different scenarios: corrections based on second-order Møller-Plesset (MP2), random-phase approximation (RPA), and coupled-cluster singles doubles (CCSD). Several chemical problems involving different types of electron correlation in singlet and multiplet spin states have been considered. Our numerical tests reveal that RPA-based and CCSD-based corrections provide similar relative errors in molecular dissociation energies (De) to the results obtained using a MP2 correction. With respect to the MP2 case, the CCSD-based correction improves the prediction, while the RPA-based correction reduces the computational cost.

Improved One-Shot Total Energies from the Linearized GW Density Matrix.

The linearized GW density matrix (γGW) is an efficient method to improve the static portion of the self-energy compared to that of ordinary perturbative GW while keeping the single-shot simplicity of the calculation. Previous work has shown that γGW gives an improved Fock operator and total energy components that approach the self-consistent GW quality. Here, we test γGW for dimer dissociation for the first time by studying N2, LiH, and Be2. We also calculate a set of self-consistent GW results in identical basis sets for a direct and consistent comparison. γGW approaches self-consistent GW total energies for a starting point based on a high amount of exact exchange. We also compare the accuracy of different total energy functionals, which differ when evaluated with a non-self-consistent density or density matrix. While the errors in total energies among different functionals and starting points are small, the individual energy components show noticeable errors when compared to reference data. The energy component errors of γGW are smaller than functionals of the density and we suggest that the linearized GW density matrix is a route to improving total energy evaluations in the adiabatic connection framework.

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.

Extrapolating Unconverged GW Energies up to the Complete Basis Set Limit with Linear Regression.

The GW approximation to the electronic self-energy is now a well-recognized approach to obtain the electron quasiparticle energies of molecules and, in particular, their ionization potential and electron affinity. Though much faster than the corresponding wavefunction methods, the GW energies are still affected by slow convergence with respect to the basis completeness. This limitation hinders a wider application of the GW approach. Here, we show that we can reach the complete basis set limit for the cumbersome GW calculations solely based on fast preliminary calculations with an unconverged basis set. We introduce a linear model that correlates the molecular orbital characteristics and the basis convergence error for a large database of approximately 600 states in 104 organic molecules that contain H, C, O, N, F, P, S, and Cl. The model employs molecular-orbital-based non-linear descriptors that encode efficiently the chemical space offering outstanding transferability. Using a low number of descriptors (17) the performance of this extrapolation procedure is superior to that of the earlier more physically motivated approaches. The predictive power of the method is finally demonstrated for a selection of large acceptor molecules.

ABINIT: Overview and focus on selected capabilities.

abinit is probably the first electronic-structure package to have been released under an open-source license about 20 years ago. It implements density functional theory, density-functional perturbation theory (DFPT), many-body perturbation theory (GW approximation and Bethe-Salpeter equation), and more specific or advanced formalisms, such as dynamical mean-field theory (DMFT) and the "temperature-dependent effective potential" approach for anharmonic effects. Relying on planewaves for the representation of wavefunctions, density, and other space-dependent quantities, with pseudopotentials or projector-augmented waves (PAWs), it is well suited for the study of periodic materials, although nanostructures and molecules can be treated with the supercell technique. The present article starts with a brief description of the project, a summary of the theories upon which abinit relies, and a list of the associated capabilities. It then focuses on selected capabilities that might not be present in the majority of electronic structure packages either among planewave codes or, in general, treatment of strongly correlated materials using DMFT; materials under finite electric fields; properties at nuclei (electric field gradient, Mössbauer shifts, and orbital magnetization); positron annihilation; Raman intensities and electro-optic effect; and DFPT calculations of response to strain perturbation (elastic constants and piezoelectricity), spatial dispersion (flexoelectricity), electronic mobility, temperature dependence of the gap, and spin-magnetic-field perturbation. The abinit DFPT implementation is very general, including systems with van der Waals interaction or with noncollinear magnetism. Community projects are also described: generation of pseudopotential and PAW datasets, high-throughput calculations (databases of phonon band structure, second-harmonic generation, and GW computations of bandgaps), and the library libpaw. abinit has strong links with many other software projects that are briefly mentioned.

Assessment of the Linearized GW Density Matrix for Molecules.

The GW approximation is well-known for the calculation of high-quality ionization potentials and electron affinities in solids and molecules. Recently, it has been identified that the density matrix that is obtained from the contraction of the GW Green's function allows one to include Feynman diagrams that are significant for the ionization potentials. However, the Green's function contains much more information than the mere quasi-particle energies. Here, we test and assess the quality of this so-called linearized GW density matrix for several molecular properties. First of all, we extend the original formulation to perform the linearization starting from any self-consistent mean-field approximation, and not only from Hartree-Fock. We then demonstrate the reliability and the stability of the linearized GW density matrix to evaluate the total energy out of a non-self-consistent GW calculation. Based on a comprehensive benchmark of 34 molecules, we compare the quality of the electronic density, Hartree energy, exchange energy, and the Fock operator expectation values against other well-established techniques. In particular, we show that the obtained linearized GW densities markedly differ from those calculated within the widespread quasi-particle self-consistent GW approximation.

Photoelectron spectra of copper oxide cluster anions from first principles methods.

We present results and analyses for the photoelectron spectra of small copper oxide cluster anions (CuO-, Cu O2- , Cu O3- , and Cu2O-). The spectra are computed using various techniques, including density functional theory (DFT) with semi-local, global hybrid, and optimally tuned range-separated hybrid functionals, as well as many-body perturbation theory within the GW approximation based on various DFT starting points. The results are compared with each other and with the available experimental data. We conclude that as in many metal-organic systems, self-interaction errors are a major issue that is mitigated by hybrid functionals. However, these need to be balanced against a strong role of non-dynamical correlation-especially in smaller, more symmetric systems-where errors are alleviated by semi-local functionals. The relative importance of the two phenomena, including practical ways of balancing the two constraints, is discussed in detail.

Optical Properties of Saturated and Unsaturated Carbonyl Defects in Polyethylene.

Polyethylene (PE), one of the simplest and most used aliphatic polymers, is generally provided with a number of additives, in particular antioxidants, because of its tendency to get oxidized. Carbonyl defects, a product of the oxidation of PE, are occurring in various forms, in particular saturated ones, known as ketones, where a C═O double bond substitutes a CH2 group, and various unsaturated ones, i.e., with further missing hydrogens. Many experimental investigations of the optical properties in the visible/UV range mainly attribute the photoluminescence of PE to one specific kind of unsaturated carbonyls, following analogies to the emission spectra of similar small molecules. However, the reason why saturated carbonyls should not be optically detected is not clear. We investigated the optical properties of PE with and without carbonyl defects using perturbative GW and the Bethe-Salpeter equation in order to take into account excitonic effects. We discuss the calculated excitonic states in comparison with experimental absorption/emission energies and the stability of both saturated and unsaturated carbonyl defects. We conclude that the unsaturated defects are indeed the best candidate for the luminescence of oxidized PE, and the reason is mainly due to oscillator strengths.

Partial Molar Volumes of Aqua Ions from First Principles.

Partial molar volumes of ions in water solution are calculated through pressures obtained from ab initio molecular dynamics simulations. The correct definition of pressure in charged systems subject to periodic boundary conditions requires access to the variation of the electrostatic potential upon a change of volume. We develop a scheme for calculating such a variation in liquid systems by setting up an interface between regions of different density. This also allows us to determine the absolute deformation potentials for the band edges of liquid water. With the properly defined pressures, we obtain partial molar volumes of a series of aqua ions in very good agreement with experimental values.

An assessment of low-lying excitation energies and triplet instabilities of organic molecules with an ab initio Bethe-Salpeter equation approach and the Tamm-Dancoff approximation.

The accurate prediction of singlet and triplet excitation energies is an area of intense research of significant fundamental interest and critical for many applications. Most calculations of singlet and triplet energies use time-dependent density functional theory (TDDFT) in conjunction with an approximate exchange-correlation functional. In this work, we examine and critically assess an alternative method for predicting low-lying neutral excitations with similar computational cost, the ab initio Bethe-Salpeter equation (BSE) approach, and compare results against high-accuracy wavefunction-based methods. We consider singlet and triplet excitations of 27 prototypical organic molecules, including members of Thiel's set, the acene series, and several aromatic hydrocarbons exhibiting charge-transfer-like excitations. Analogous to its impact in TDDFT, we find that the Tamm-Dancoff approximation (TDA) overcomes triplet instabilities in the BSE approach, improving both triplet and singlet energetics relative to higher level theories. Finally, we find that BSE-TDA calculations built on effective DFT starting points, such as those utilizing optimally tuned range-separated hybrid functionals, can yield accurate singlet and triplet excitation energies for gas-phase organic molecules.

Benchmarking the GW Approximation and Bethe-Salpeter Equation for Groups IB and IIB Atoms and Monoxides.

Energies from the GW approximation and the Bethe-Salpeter equation (BSE) are benchmarked against the excitation energies of transition-metal (Cu, Zn, Ag, and Cd) single atoms and monoxide anions. We demonstrate that best estimates of GW quasiparticle energies at the complete basis set limit should be obtained via extrapolation or closure relations, while numerically converged GW-BSE eigenvalues can be obtained on a finite basis set. Calculations using real-space wave functions and pseudopotentials are shown to give best-estimate GW energies that agree (up to the extrapolation error) with calculations using all-electron Gaussian basis sets. We benchmark the effects of a vertex approximation (ΓLDA) and the mean-field starting point in GW and the BSE, performing computations using a real-space, transition-space basis and scalar-relativistic pseudopotentials. While no variant of GW improves on perturbative G0W0 at predicting ionization energies, G0W0ΓLDA-BSE computations give excellent agreement with experimental absorption spectra as long as off-diagonal self-energy terms are included. We also present G0W0 quasiparticle energies for the CuO-, ZnO-, AgO-, and CdO- anions, in comparison to available anion photoelectron spectra.

Optimized virtual orbital subspace for faster GW calculations in localized basis.

The popularity of the GW approximation to the self-energy to access the quasiparticle energies of molecules is constantly increasing. As the other methods addressing the electronic correlation, the GW self-energy unfortunately shows a very slow convergence with respect to the basis complexity, which precludes the calculation of accurate quasiparticle energies for large molecules. Here we propose a method to mitigate this issue that relies on two steps: (i) the definition of a reduced virtual orbital subspace, thanks to a much smaller basis set; (ii) the account of the remainder through the simpler one-ring approximation to the self-energy. We assess the quality of the corrected quasiparticle energies for simple molecules, and finally we show an application to large graphene chunks to demonstrate the numerical efficiency of the scheme.

Evaluating the GW Approximation with CCSD(T) for Charged Excitations Across the Oligoacenes.

Charged excitations of the oligoacene family of molecules, relevant for astrophysics and technological applications, are widely studied and therefore provide an excellent system for benchmarking theoretical methods. In this work, we evaluate the performance of many-body perturbation theory within the GW approximation relative to new high-quality CCSD(T) reference data for charged excitations of the acenes. We compare GW calculations with a number of hybrid density functional theory starting points and with eigenvalue self-consistency. Special focus is given to elucidating the trend of GW-predicted excitations with molecule length increasing from benzene to hexacene. We find that GW calculations with starting points based on an optimally tuned range-separated hybrid (OTRSH) density functional and eigenvalue self-consistency can yield quantitative ionization potentials for the acenes. However, for larger acenes, the predicted electron affinities can deviate considerably from reference values. Our work paves the way for predictive and cost-effective GW calculations of charged excitations of molecules and identifies certain limitations of current GW methods used in practice for larger molecules.