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.     

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.     

On the equivalence of ring-coupled cluster and adiabatic connection fluctuation-dissipation theorem random phase approximation correlation energy expressions

The correlation energy in the direct random phase approximation (dRPA) can be written, among other possibilities, either in terms of the interaction strength averaged correlation density matrix, or in terms of the coupled cluster doubles amplitudes obtained in the direct ring approximation (drCCD). Although the corresponding dRPA correlation density matrix on the one hand, and the drCCD amplitude matrix on the other hand, differ significantly, they yield identical energies. Similarly, the analogous RPA and rCCD correlation energies calculated from antisymmetrized two-electron integrals are identical to each other despite very different underlying working equations. In the present communication, a direct correspondence between amplitudes and densities is established and investigated with perturbation theory arguments. Our analysis also sheds some light on the properties of recently proposed RPA/rCCD variants which use antisymmetrized integrals in part of the equations and nonantisymmetrized integrals in others.     

Long-range-corrected hybrids using a range-separated Perdew-Burke-Ernzerhof functional and random phase approximation correlation

We build on methods combining a short-range density functional approximation with a long-range random phase approximation [B. G. Janesko, T. M. Henderson, and G. E. Scuseria, J. Chem. Phys. 130, 081105 (2009)] or second-order screened exchange [J. Paier et al., J. Chem. Phys. 132, 094103 (2010)] by replacing the range-separated local density approximation functional with a range-separated generalized gradient approximation functional in the short range. We present benchmark results that show a marked improvement in the thermodynamic tests over the previous local density approximation-based methods while retaining those methods' excellent performance in van der Waals interactions.     

Theoretical investigation of carrier envelope phase dependence of all-optical poling with the third- and higher-order nonlinear processes

Dependence of efficiency in all-optical poling with nonlinear processes, up to eighth order, is considered. The explicit form of the nonlinear susceptibility that is responsible for the poling is derived, which shows both CEP and phase mismatch dependence. On the basis of an analysis of pulse propagation in a nonlinear material, it is shown that one can identify the order of nonlinearity that is relevant to the poling process, relying on current technology of CEP stabilization and thin-film growth.     

Dynamic polarizability of helium: A random phase approximation calculation

The random phase approximation (RPA ) or time-dependent Hartree–Fock approximation (TDHF ) is reconsidered for the calculation of the dynamic polarizability for atoms. An integral equation which admits a simple numerical treatment is established. The asymptotic approximation for the electron propagator is tested for its applicability by means of comparisons with earlier results.     

Ionization energy of atoms obtained from GW self-energy or from random phase approximation total energies

A systematic evaluation of the ionization energy within the GW approximation is carried out for the first row atoms, from H to Ar. We describe a Gaussian basis implementation of the GW approximation, which does not resort to any further technical approximation, besides the choice of the basis set for the electronic wavefunctions. Different approaches to the GW approximation have been implemented and tested, for example, the standard perturbative approach based on a prior mean-field calculation (Hartree-Fock GW@HF or density-functional theory GW@DFT) or the recently developed quasiparticle self-consistent method (QSGW). The highest occupied molecular orbital energies of atoms obtained from both GW@HF and QSGW are in excellent agreement with the experimental ionization energy. The lowest unoccupied molecular orbital energies of the singly charged cation yield a noticeably worse estimate of the ionization energy. The best agreement with respect to experiment is obtained from the total energy differences within the random phase approximation functional, which is the total energy corresponding to the GW self-energy. We conclude with a discussion about the slight concave behavior upon number electron change of the GW approximation and its consequences upon the quality of the orbital energies.     

Short-range exchange and correlation energy density functionals: beyond the local-density approximation

We propose approximations which go beyond the local-density approximation for the short-range exchange and correlation density functionals appearing in a multideterminantal extension of the Kohn-Sham scheme. A first approximation consists of defining locally the range of the interaction in the correlation functional. Another approximation, more conventional, is based on a gradient expansion of the short-range exchange-correlation functional. Finally, we also test a short-range generalized-gradient approximation by extending the Perdew-Burke-Ernzerhof exchange-correlation functional to short-range interactions.     

Characteristics of the consistent ground state of the random phase approximation

Detailed considerations of the ground-state vector derived for the random phase approximation obtained earlier in a generator coordinate representation employing the unitary group parameter space reveal a particular correlated pair structure. The results of explicit calculations of ground-state averages are discussed.     

Justifiability of the ZDO approximation in terms of a power series expansion

The Zero Differential Overlap Approximation cannot be justified for all-valence-electron calculations in terms of a power series expansion of the overlap matrix, because the expansion diverges.     

The principle-quantum-number (and the radial-quantum-number) expansion of the correlation energy of two-electron atoms

The second-order correlation energy of two-electron ions is studied in terms of an expansion in minimal approximations to the first-order natural orbitals (NOs). The non-linear parameters of these NOs are determined by minimization of the second-order energy. An approximation to the total second-order correlation energy is obtained as a sum of increments e(lp), depending on the angular quantum number l and the radial quantum number p. (Either l or p can be eliminated in favor of the principal quantum number n = l + p.) Closed expressions for these energy increments are derived. For fixed p the increments go as (l + 1)(-5). This is consistent with the behavior of the exact partial wave increments (that depend on the parameter l only) as (l + 1/2)(-4). While the partial wave increments correspond to a summation of e(lp) over p, other partial summations of the two-parameter increments lead to either the principal-quantum-number expansion (PQNE) with energy increments approximately n(-4), or the radial-quantum-number expansion, with a less transparent convergence pattern. Unfortunately these partial summations can neither be done in closed form nor from the asymptotic expansion, but some insight is obtained from a numerical summation. The hope to find a rigorous derivation of the PQNE has not been fulfilled.     

Calculation of molecular polarizabilities for large systems within the random phase approximation

We examine the static-field molecular polarizability from a sum over uncoupled Hartree—Fock states (SOS), the Tamm—Dancoff approximation (TDA), and the random phase approximation (RPA). An efficient algorithm for the inversion of the TDA or RPA matrix is outlined, which avoids matrix diagonalization and explicit construction of matrix elements over states, allowing for rapid calculation of molecular polarizabilities. The extension of the method is straightforward; third-order hyperpolarizability is developed as an example. Test cases are reported for molecules represented by an intermediate neglect of differential overlap (INDO) wavefunction.     

The first and second derivative matrices in the random phase approximation scheme by using the Lagrangian technique

We have presented the explicit formulas for first and second derivatives of A and B matrices, appearing in the random phase approximation (RPA), with the aid of Lagrangian technique. Owing to the 2n + 1 rule, the Lagrangian approach is more efficient than the conventional approach to evaluate the higher‐order matrix elements. We have confirmed the validity of our formulation by demonstrating the geometry optimization of the first‐excited singlet states of formaldehyde, ethylene, and 1‐amino‐3‐propenal molecules. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

On the role of higher-order correlation effects on the induction interactions between closed-shell molecules

High-order correlation contributions to the second-order induction energy were studied for various representative van der Waals complexes. It was found that the induction energy obtained by the truncation of the relaxed M?ller-Plesset expansion in the second or third order is in most cases quite close to the induction energy computed with the coupled-cluster method (restricted to single and double excitations). Also, the effect of triples excitations on this perturbation term is usually small. However, given an oscillatory behaviour of the M?ller-Plesset induction corrections, the coupled-cluster method seems to be better suited to a reliable calculation of the induction energy. The sources of the remaining differences between the interaction energies computed by symmetry-adapted perturbation theory and those computed by the supermolecule coupled-cluster method (restricted to single, double, and noniterative triple excitations) are examined. It has been found that they can be attributed to the higher-order correlation terms in the second-order dispersion and exchange-induction corrections.     

Connected moments expansion calculations of the correlation energy in small molecules

The second-order connected moments expansion (CMX(2)) approach to calculation of the correlation energy is tested numerically on several closed-shell di- and tri-atomic molecules. Benchmark computations performed within 6–31G⁷ basis set reveal that CMX(2) usually recovers more than 50% of the MP3 correlation energy and improves the SCF molecular geometries at a cost comparable to the MP3 calculations.     

Dispersion and repulsion contributions to the solvation energy: Refinements to a simple computational model in the continuum approximation

A computational method for the evaluation of dispersion and repulsion contributions to the solvation energy is here presented in a formulation which makes use of a continuous distribution of the solvent, without introducing additional assumptions (e.g., local isotropy in the solvent distribution). The analysis is addressed to compare the relative importance of the various components of the dispersion energy (n = 6, 8, 10) and of the repulsion term, to compare several molecular indicators (molecular surface and volume, number of electrons) which may be put in relation to the dispersion-repulsion energy, and to define simplified computational strategies. The numerical examples refer to saturated hydrocarbons in water, treated with the homogeneous approximation of the distribution function which for this type of solution appears to be acceptable.     

Resolution of identity approach for the Kohn-Sham correlation energy within the exact-exchange random-phase approximation

Two related methods to calculate the Kohn-Sham correlation energy within the framework of the adiabatic-connection fluctuation-dissipation theorem are presented. The required coupling-strength-dependent density-density response functions are calculated within exact-exchange time-dependent density-functional theory, i.e., within time-dependent density-functional response theory using the full frequency-dependent exchange kernel in addition to the Coulomb kernel. The resulting resolution-of-identity exact-exchange random-phase approximation (RI-EXXRPA) methods in contrast to previous EXXRPA methods employ an auxiliary basis set (RI basis set) to improve the computational efficiency, in particular, to reduce the formal scaling of the computational effort with respect to the system size N from N(6) to N(5). Moreover, the presented RI-EXXRPA methods, in contrast to previous ones, do not treat products of occupied times unoccupied orbitals as if they were linearly independent. Finally, terms neglected in previous EXXRPA methods can be included, which leads to a method designated RI-EXXRPA+, while the method without these extra terms is simply referred to as RI-EXXRPA. Both EXXRPA methods are shown to yield total energies, reaction energies of small molecules, and binding energies of noncovalently bonded dimers of a quality that is similar and in some cases even better than that obtained with quantum chemistry methods such as Mo?ller-Plesset perturbation theory of second order (MP2) or with the coupled cluster singles doubles method. In contrast to MP2 and to conventional density-functional methods, the presented RI-EXXRPA methods are able to treat static correlation.