Increasing the applicability of density functional theory. II. Correlation potentials from the random phase approximation and beyond

Density functional theory (DFT) results are mistrusted at times due to the presence of an unknown exchange correlation functional, with no practical way to guarantee convergence to the right answer. The use of a known exchange correlation functional based on wave-function theory helps to alleviate such mistrust. The exchange correlation functionals can be written exactly in terms of the density-density response function using the adiabatic-connection and fluctuation-dissipation framework. The random phase approximation (RPA) is the simplest approximation for the density-density response function. Since the correlation functional obtained from RPA is equivalent to the direct ring coupled cluster doubles (ring-CCD) correlation functional, meaning only Coulomb interactions are included, one can bracket RPA between many body perturbation theory (MBPT)-2 and CCD with the latter having all ring, ladder, and exchange contributions. Using an optimized effective potential strategy, we obtain correlation potentials corresponding to MBPT-2, RPA (ring-CCD), linear-CCD, and CCD. Using the suitable choice of the unperturbed Hamiltonian, Kohn-Sham self-consistent calculations are performed. The spatial behavior of the resulting potentials, total energies, and the HOMO eigenvalues are compared with the exact values for spherical atoms. Further, we demonstrate that the self-consistent eigenvalues obtained from these consistent potentials used in ab initio dft approximate all principal ionization potentials as demanded by ionization potential theorem.     

Quantum‐based models of charge‐dependent potential energy surfaces: Three‐state models

Quantum‐based models of how potential energies depend on charge are developed from a three‐state model, at the level of neglecting state‐to‐state overlap. The energy as a function of charge is defined as proposed previously (Valone and Atlas, J Chem Phys 2004, 120, 7262). With this definition, addition of a third state smooths the derivatives of the energy model with respect to charge at integer values of charge that are in the interior of the allowed charge range. These derivatives are related to the chemical potential. At the dissociation limit, this model converges to established limits. Another dependence is proposed that uses two different charges simultaneously. The concepts are illustrated, with calculations on an OH molecule. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Third-order corrections to random-phase approximation correlation energies

Several random-phase approximation (RPA) correlation methods were compared in third order of perturbation theory. While all of the considered approaches are exact in second order of perturbation theory, it is found that their corresponding third-order correlation energy contributions strongly differ from the exact third-order correlation energy contribution due to missing interactions of the particle-particle-hole-hole type. Thus a simple correction method is derived which makes the different RPA methods also exact to third-order of perturbation theory. By studying the reaction energies of 16 chemical reactions for 21 small organic molecules and intermolecular interaction energies of 23 intermolecular complexes comprising weakly bound and hydrogen-bridged systems, it is found that the third-order correlation energy correction considerably improves the accuracy of RPA methods if compared to coupled-cluster singles doubles with perturbative triples as a reference.     

Homolytic molecular dissociation in natural orbital functional theory

The dissociation of diatomic molecules of the 14-electron isoelectronic series N(2), O(2)(2+), CO, CN(-) and NO(+) is examined using the Piris natural orbital functional. It is found that the method describes correctly the dissociation limit yielding an integer number of electrons on the dissociated atoms, in contrast to the fractional charges obtained when using the variational two-particle reduced density matrix method under the D, Q and G positivity necessary N-representability conditions. The chemistry of the considered systems is discussed in terms of their dipole moments, natural orbital occupations and bond orders as well as atomic Mulliken populations at the dissociation limit. The values obtained agree well with accurate multiconfigurational wave function based CASSCF results and the available experimental data.     

Density functionals that are one- and two- are not always many-electron self-interaction-free, as shown for H2+, He2+, LiH+, and Ne2+

The common density functionals for the exchange-correlation energy make serious self-interaction errors in the molecular dissociation limit when real or spurious noninteger electron numbers N are found on the dissociation products. An "M-electron self-interaction-free" functional for positive integer M is one that produces a realistic linear variation of total energy with N in the range of M-12. Thus all these SIC's produce an exact binding energy curve for H2+, and an accurate one for He2+, but only the unscaled Perdew-Zunger SIC produces an accurate one for Ne2+, where there are more than two electrons on each fragment Ne+0.5. We also discuss LiH+, which is relatively free from self-interaction errors. We suggest that the ability of the original and unscaled Perdew-Zunger SIC to be nearly M-electron self-interaction-free for atoms of all M stems in part from its formal resemblance to the Hartree-Fock theory, with which it shares a sum rule on the exchange-correlation hole of an open system.     

Wavelength-decomposition-based embedded cluster density approximation for systems with nonlocal electron correlation

Local correlation methods rely on the assumption that electron correlation is nearsighted. In this work, we develop a method to alleviate this assumption. This new method is demonstrated by calculating the random phase approximation (RPA) correlation energies in several one-dimensional model systems. In this new method, the first step is to approximately decompose the RPA correlation energy to the nearsighted and farsighted components based on the wavelength decomposition of electron correlation developed by Langreth and Perdew. The short-wavelength (SW) component of the RPA correlation energy is then considered to be nearsighted, and the long-wavelength (LW) component of the RPA correlation energy is considered to be farsighted. The SW RPA correlation energy is calculated using a recently developed local correlation method: the embedded cluster density approximation (ECDA). The LW RPA correlation energy is calculated globally based on the system's Kohn-Sham orbitals. This new method is termed λ-ECDA, where λ indicates the wavelength decomposition. The performance of λ-ECDA is examined on a one-dimensional model system: a H₂₄ chain, in which the RPA correlation energy is highly nonlocal. In this model system, a softened Coulomb interaction is used to describe the electron-electron and electron-ion interactions, and slightly stronger nuclear charges (1.2e ) are assigned to the pseudo-H atoms. Bond stretching energies, RPA correlation potentials, and Kohn-Sham eigenvalues predicted by λ-ECDA are in good agreement with the benchmarks when the clusters are made reasonably large. We find that the LW RPA correlation energy is critical for obtaining accurate prediction of the RPA correlation potential, even though the LW RPA correlation energy contributes to only a few percent of the total RPA correlation energy.     

Fractional spin in reduced density-matrix functional theory

We study the behavior of different functionals of the one-body reduced density matrix (1RDM) for systems with fractional z-component of the total spin. We define these systems as ensembles of integer spin states. It is shown that, similarly to density functional theory, the error in the dissociation of diatomic molecules is directly related to the deviation from constancy of the atomic total energies as functions of the fractional spin. However, several functionals of the 1RDM show a size inconsistency which leads to additional errors. We also investigate the difference between a direct evaluation of the energy of an ensemble of integer-spin systems and a direct minimization of the energy of a fractional-spin system.     

Excitation energies from range-separated time-dependent density and density matrix functional theory

Time-dependent density functional theory (TD-DFT) in the adiabatic formulation exhibits known failures when applied to predicting excitation energies. One of them is the lack of the doubly excited configurations. On the other hand, the time-dependent theory based on a one-electron reduced density matrix functional (time-dependent density matrix functional theory, TD-DMFT) has proven accurate in determining single and double excitations of H(2) molecule if the exact functional is employed in the adiabatic approximation. We propose a new approach for computing excited state energies that relies on functionals of electron density and one-electron reduced density matrix, where the latter is applied in the long-range region of electron-electron interactions. A similar approach has been recently successfully employed in predicting ground state potential energy curves of diatomic molecules even in the dissociation limit, where static correlation effects are dominating. In the paper, a time-dependent functional theory based on the range-separation of electronic interaction operator is rigorously formulated. To turn the approach into a practical scheme the adiabatic approximation is proposed for the short- and long-range components of the coupling matrix present in the linear response equations. In the end, the problem of finding excitation energies is turned into an eigenproblem for a symmetric matrix. Assignment of obtained excitations is discussed and it is shown how to identify double excitations from the analysis of approximate transition density matrix elements. The proposed method used with the short-range local density approximation (srLDA) and the long-range Buijse-Baerends density matrix functional (lrBB) is applied to H(2) molecule (at equilibrium geometry and in the dissociation limit) and to Be atom. The method accounts for double excitations in the investigated systems but, unfortunately, the accuracy of some of them is poor. The quality of the other excitations is in general much better than that offered by TD-DFT-LDA or TD-DMFT-BB approximations if the range-separation parameter is properly chosen. The latter remains an open problem.     

Electron correlation methods based on the random phase approximation

In the past decade, the random phase approximation (RPA) has emerged as a promising post-Kohn–Sham method to treat electron correlation in molecules, surfaces, and solids. In this review, we explain how RPA arises naturally as a zero-order approximation from the adiabatic connection and the fluctuation-dissipation theorem in a density functional context. This is contrasted to RPA with exchange (RPAX) in a post-Hartree–Fock context. In both methods, RPA and RPAX, the correlation energy may be expressed as a sum over zero-point energies of harmonic oscillators representing collective electronic excitations, consistent with the physical picture originally proposed by Bohm and Pines. The extra factor 1/2 in the RPAX case is rigorously derived. Approaches beyond RPA are briefly summarized. We also review computational strategies implementing RPA. The combination of auxiliary expansions and imaginary frequency integration methods has lead to recent progress in this field, making RPA calculations affordable for systems with over 100 atoms. Finally, we summarize benchmark applications of RPA to various molecular and solid-state properties, including relative energies of conformers, reaction energies involving weak and covalent interactions, diatomic potential energy curves, ionization potentials and electron affinities, surface adsorption energies, bulk cohesive energies and lattice constants. RPA barrier heights for an extended benchmark set are presented. RPA is an order of magnitude more accurate than semi-local functionals such as B3LYP for non-covalent interactions rivaling the best empirically parametrized methods. Larger but systematic errors are observed for processes that do not conserve the number of electron pairs, such as atomization and ionization.     

Optimized effective potential from a correlated wave function: optimized effective potential-generalized valence bond (OEP-GVB)

The optimized effective potential (OEP) method allows orbital-dependent functionals to be used in density functional theory. Traditionally the orbital-dependent functional of interest has been the Hartree-Fock energy, leading to exact exchange density functional theory. Here we present results that use a generalized valence-bond (GVB) wave function, a multiconfigurational wave function that includes static correlation and dissociates to the proper limits. We demonstrate the effectiveness of the OEP-GVB method by showing the dissociation of H2 and the excitation spectrum of He.     

A simple but fully nonlocal correction to the random phase approximation

The random phase approximation (RPA) stands on the top rung of the ladder of ground-state density functional approximations. The simple or direct RPA has been found to predict accurately many isoelectronic energy differences. A nonempirical local or semilocal correction to this direct RPA leaves isoelectronic energy differences almost unchanged, while improving total energies, ionization energies, etc., but fails to correct the RPA underestimation of molecular atomization energies. Direct RPA and its semilocal correction may miss part of the middle-range multicenter nonlocality of the correlation energy in a molecule. Here we propose a fully nonlocal, hybrid-functional-like addition to the semilocal correction. The added full nonlocality is important in molecules, but not in atoms. Under uniform-density scaling, this fully nonlocal correction scales like the second-order-exchange contribution to the correlation energy, an important part of the correction to direct RPA, and like the semilocal correction itself. For the atomization energies of ten molecules, and with the help of one fit parameter, it performs much better than the elaborate second-order screened exchange correction.     

Exact expressions for ensemble functionals from particle number dependence

Some properties of exact ensemble density functionals can be determined by examining the particle number dependence of ground state ensemble density matrices for systems where the integer ground state energies satisfy a convexity condition. The results include the observation that the integral of the product of the functional derivative and Fukui function of functionals that can be expressed as the trace of an operator is particle number independent for particle numbers between successive integers and the integral itself is equal to the difference between functionals evaluated at successive integer particle numbers. Expressions that must be satisfied by 2nd and higher order functional derivatives are formulated and equations that must be satisfied point by point in space are derived. Using the analytic Hooke's atom model, it is shown that commonly used correlation functional approximations do not bear any resemblance to a spatially dependent expression derived from the exact second order functional derivative of the correlation functional. It is also shown that two expressions for the mutual Coulomb energy are not equal when approximate exchange and correlation functionals are used.     

Spurious fractional charge on dissociated atoms: pervasive and resilient self-interaction error of common density functionals

Semilocal density functional approximations for the exchange-correlation energy can improperly dissociate a neutral molecule XY (Y not =X) to fractionally charged fragments X(+q)...Y(-q) with an energy significantly lower than X0...Y0. For example, NaCl can dissociate to Na(+0.4)...Cl(-0.4). Generally, q is positive when the lowest-unoccupied orbital energy of atom Y0 lies below the highest-occupied orbital energy of atom X0. The first 24 open sp-shell atoms of the Periodic Table can form 276 distinct unlike pairs XY, and in the local spin density approximation 174 of these display fractional-charge dissociation. Finding these lowest-energy solutions with standard quantum chemistry codes, however, requires special care. Self-interaction-corrected (SIC) semilocal approximations are exact for one-electron systems and also reduce the spurious fractional charge q. The original SIC of Perdew and Zunger typically reduces q to 0. A scaled-down SIC with better equilibrium properties sometimes fails to reduce q all the way to 0. The desideratum of "many-electron self-interaction freedom" is introduced as a generalization of the one-electron concept.     

Microphase separation in polyelectrolytic diblock copolymer melt: weak segregation limit

The authors present a generalized theory of microphase separation for charged-neutral diblock copolymer melt. The stability limit of the disordered phase for salt-free melt has been calculated using random phase approximation (RPA) and self-consistent-field theory (SCFT). Explicit analytical free energy expressions for different classical ordered microstructures (lamellar, cylinder, and sphere) are presented. The authors demonstrate that the chemical mismatch required for the onset of microphase separation (chi*N) in charged-neutral diblock melt is higher and the period of ordered microstructures is lower than those for the corresponding neutral-neutral diblock system. Theoretical predictions on the period of ordered structures in terms of Coulomb electrostatic interaction strength, chain length, block length, and chemical mismatch between blocks are presented. SCFT has been used to go beyond the stability limit, where electrostatic potential and charge distribution are calculated self-consistently. Stability limits calculated using RPA are in perfect agreement with the corresponding SCFT calculations. Limiting laws for the stability limit and the period of ordered structures are presented and comparisons are made with an earlier theory. Also, transition boundaries between different morphologies have been investigated.     

Dissociative double photoionization of N2 using synchrotron radiation: appearance energy of the N2+ dication

Photoionization cross sections for the production of the doubly charged ion N2+ from N2 have been measured by means of synchrotron radiation in the photon energy range from 50 to 110 eV. The appearance energy for N2+ has been determined as 55.2+/-0.2 eV, i.e., about 1.3 eV higher than the spectroscopic dissociation limit leading to the charge asymmetric dissociation channel N2+(2P)+N(4S) at 53.9 eV. The onset of a second threshold at 59.9+/-0.2 eV is detected and the energy dependence of photoion intensities near the threshold regions is interpreted in terms of the Wannier theory. The production of the N2+ dication is discussed in terms of direct and indirect mechanisms for dissociative charge asymmetric photoionization and by comparison with the potential energy curves of the intermediate N(2)2+ dication. Experimental evidences for the opening of the Coulomb explosion channel N2++N+ at high photon energies are provided by measuring the kinetic energy release spectra of N2+ fragments at selected photon energies.     

Basis set convergence of molecular correlation energy differences within the random phase approximation

The basis set convergence of energy differences obtained from the random phase approximation (RPA) to the correlation energy is investigated for a wide range of molecular interactions. For dispersion bound systems the basis set incompleteness error is most pronounced, as shown for the S22 benchmark [P. Jurecka et al., Phys. Chem. Chem. Phys. 8, 1985 (2006)]. The use of very large basis sets (> quintuple-zeta) or extrapolation to the complete basis set (CBS) limit is necessary to obtain a reliable estimate of the binding energy for these systems. Counterpoise corrected results converge to the same CBS limit, but counterpoise correction without extrapolation is insufficient. Core-valence correlations do not play a significant role. For medium- and short-range correlation, quadruple-zeta results are essentially converged, as demonstrated for relative alkane conformer energies, reaction energies dominated by intramolecular dispersion, isomerization energies, and reaction energies of small organic molecules. Except for weakly bound systems, diffuse augmentation almost universally slows down basis set convergence. For most RPA applications, quadruple-zeta valence basis sets offer a good balance between accuracy and efficiency.     

Increasing physical constraints and improving performances in a parameter-free GGA functional

In this Letter we present an evolution of the TCA functional [V. Tognetti, P. Cortona, C. Adamo, J. Chem. Phys. 128 (2008) 034101] that gives a vanishing correlation energy for hydrogenoid atoms. This feature of the exact functional is incorporated in the TCA approximation at the generalized-gradient level, without kinetic energy density or higher order density derivatives dependence. A significant improvement of atomic and atomization energies, activation barriers and ionization potentials is found when the new correlation is coupled with modified PBE exchange functionals, while structural parameters are close in accuracy to those provided by the original PBE functional.     

Open-shell reduced density matrix functional theory

Open-shell reduced density matrix functional theory is established by investigating the domain of the exact functional. For spin states that are the ground state, a particularly simple set is found to be the domain. It cannot be generalized to other spin states. A number of conditions satisfied by the exact density matrix functional is formulated and tested for approximate functionals. The exact functional does not suffer from fractional spin error, which is the source of the static correlation error in dissociated molecules. We prove that a simple approximation (called the Buijse-Baerends functional, Mu?ller or square root functional) has a non-positive fractional spin error. In the case of the H atom the error is zero. Numerical results for a few atoms are given for approximate density and density matrix functionals as well as a recently developed range-separated combination of both.