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.     

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.     

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.     

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.     

Accuracy of basis-set extrapolation schemes for DFT-RPA correlation energies in molecular calculations

We construct a reference benchmark set for atomic and molecular random phase approximation (RPA) correlation energies in a density functional theory framework at the complete basis-set limit. This set is used to evaluate the accuracy of some popular extrapolation schemes for RPA all-electron molecular calculations. The results indicate that for absolute energies, accurate results, clearly outperforming raw data, are achievable with two-point extrapolation schemes based on quintuple- and sextuple-zeta basis sets. Moreover, we show that results in good agreement with the benchmark can also be obtained by using a semiempirical extrapolation procedure based on quadruple- and quintuple-zeta basis sets. Finally, we analyze the performance of different extrapolation schemes for atomization energies.     

Noniterative Doubles Corrections to the Random Phase and Higher Random Phase Approximations: Singlet and Triplet Excitation Energies

The second-order noniterative doubles-corrected random phase approximation (RPA) method has been extended to triplet excitation energies and the doubles-corrected higher RPA method as well as a shifted version for calculating singlet and triplet excitation energies are presented here for the first time. A benchmark set consisting of 20 molecules with a total of 117 singlet and 71 triplet excited states has been used to test the performance of the new methods by comparison with previous results obtained with the second-order polarization propagator approximation (SOPPA) and the third order approximate coupled cluster singles, doubles and triples model CC3. In general, the second-order doubles corrections to RPA and HRPA significantly reduce both the mean deviation as well as the standard deviation of the errors compared to the CC3 results. The accuracy of the new methods approaches the accuracy of the SOPPA method while using only 10–60% of the calculation time. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.     

Density-functional theory with effective potential expressed as a direct mapping of the external potential: applications to atomization energies and ionization potentials

In this paper the authors further develop and apply the direct-mapping density functional theory to calculations of the atomization energies and ionization potentials. Single-particle orbitals are determined by solving the Kohn-Sham [Phys. Rev. A. 140, 1133 (1965)] equations with a local effective potential expressed in terms of the external potential. A two-parametric form of the effective potential for molecules is proposed and equations for optimization of the parameters are derived using the exchange-only approximation. Orbital-dependent correlation functional is derived from the second-order perturbation theory in its Moller-Plesset-type zeroth-order approximation based on the Kohn-Sham orbitals and orbital energies. The total atomization energies and ionization potentials computed with the second-order perturbation theory were found to be in agreement with experimental values and benchmark results obtained with ab initio wave mechanics methods.     

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.     

Random-phase-approximation-based correlation energy functionals: benchmark results for atoms

The random phase approximation for the correlation energy functional of the density functional theory has recently attracted renewed interest. Formulated in terms of the Kohn-Sham orbitals and eigenvalues, it promises to resolve some of the fundamental limitations of the local density and generalized gradient approximations, as, for instance, their inability to account for dispersion forces. First results for atoms, however, indicate that the random phase approximation overestimates correlation effects as much as the orbital-dependent functional obtained by a second order perturbation expansion on the basis of the Kohn-Sham Hamiltonian. In this contribution, three simple extensions of the random phase approximation are examined; (a) its augmentation by a local density approximation for short-range correlation, (b) its combination with the second order exchange term, and (c) its combination with a partial resummation of the perturbation series including the second order exchange. It is found that the ground state and correlation energies as well as the ionization potentials resulting from the extensions (a) and (c) for closed subshell atoms are clearly superior to those obtained with the unmodified random phase approximation. Quite some effort is made to ensure highly converged data, so that the results may serve as benchmark data. The numerical techniques developed in this context, in particular, for the inherent frequency integration, should also be useful for applications of random phase approximation-type functionals to more complex systems.     

Transport properties of boron and aluminum

Transport cross sections and collision integrals are tabulated for a wide range of energies and temperatures for the interactions B–B and Al–Al. For aluminum, a semiclassical approximation was used to determine the scattering phase shifts from which the transport cross sections were calculated. For boron, the smaller reduced mass and the deep potential wells required the phase shifts at lower energies to be determined from a numerical solution of the time-independent Schroedinger equation; the semiclassical approximation was used at higher energies where the two methods agree. The variations of the collision integrals for viscosity and diffusion are presented graphically as a function of temperature. The results are applied to estimate the transport properties of gallium. Received: 9 July 1999 / Accepted: 18 August 1999 / Published online: 2 November 1999     

The influence of electron correlation on reaction energies. The dimerization energies of BH3 and LiH

Results of rigorous computations employing extended Gaussian-type basis sets are reported for BH₃, B₂H₆, LiH, and Li₂H₂ in their respective equilibrium geometries. The dimerization energy of BH₃ is calculated as −20.7 kcal/mol within the Hartree-Fock approximation and as −36.6 kcal/mol if electron correlation is included. The corresponding results for the dimerization of LiH are −47.3 kcal/mol and −48.3 kcal/mol. Partitioning of the correlation energy contributions allows to attribute the effect of electron correlation to the increase of next neighbour bond interactions on the dimerization of BH₃ and LiH. The difficulties of accurate computations of reaction energies are discussed in detail.     

Fluctuation-dissipation theorem density-functional theory

Using the fluctuation-dissipation theorem (FDT) in the context of density-functional theory (DFT), one can derive an exact expression for the ground-state correlation energy in terms of the frequency-dependent density response function. When combined with time-dependent density-functional theory, a new class of density functionals results that use approximations to the exchange-correlation kernel fxc as input. This FDT-DFT scheme holds promise to solve two of the most distressing problems of conventional Kohn-Sham DFT: (i) It leads to correlation energy functionals compatible with exact exchange, and (ii) it naturally includes dispersion. The price is a moderately expensive O(N6) scaling of computational cost and a slower basis set convergence. These general features of FDT-DFT have all been recognized previously. In this paper, we present the first benchmark results for a set of molecules using FDT-DFT beyond the random-phase approximation (RPA)-that is, the first such results with fxc not equal to 0. We show that kernels derived from the adiabatic local-density approximation and other semilocal functionals suffer from an "ultraviolet catastrophe," producing a pair density that diverges at small interparticle distance. Nevertheless, dispersion interactions can be treated accurately if hybrid functionals are employed, as is demonstrated for He2 and HeNe. We outline constraints that future approximations to fxc should satisfy and discuss the prospects of FDT-DFT.     

Cyanovinyl radical: an illustration of the poor performance of unrestricted perturbation theory and density functional theory procedures in calculating radical stabilization energies

Stabilization energies for the 1-cyanovinyl radical (CH₂=CCN) have been calculated using a variety of conventional ab initio (M?ller–Plesset, quadratic configuration interaction and coupled-cluster) and density functional theory (B-LYP, B3-LYP) procedures, as well as with a range of compound methods. Compared with a high-level benchmark value (that predicts a stabilization energy of 17.1 kJ mol⁻¹), UMP2 and UMP4 give the wrong sign and magnitude of the stabilization energy (both methods predicting desta- bilization instead of stabilization), while B-LYP and B3-LYP overestimate the degree of stabilization. The RMP2, RMP4, QCISD(T) and CCSD(T) techniques, and several, but not all, variants of G2 and CBS theories give radical stabilization energies in good agreement with the benchmark value. Received: 15 June 1998 / Accepted: 19 August 1998 /  Published online: 15 February 1999     

Assessment of the Van Voorhis–Scuseria exchange–correlation functional for predicting excitation energies using time-dependent density functional theory

We assess the performance of the Van Voorhis–Scuseria exchange–correlation functional (VSXC), a kinetic-energy-density-dependent exchange–correlation functional recently developed in our group, for calculating vertical excitation energies using time-dependent density functional theory in a benchmark set of molecules. Overall, VSXC performs very well, with accuracy similar to that of hybrid functionals such as the hybrid Perdew–Burke–Ernzerhof functional and Becke's three parameter hybrid method with the Lee, Yang, and Parr correlation functional, which contain a portion of Hartree–Fock exchange. Received: 29 December 1999 / Accepted: 5 June 2000 / Published online: 11 September 2000     

Relative Cross Sections for L- and M-Shell Ionization by Electron Impact

 Results from measurements and calculations of relative L- and M-shell ionization cross sections by electron impact are presented. Measurements were performed for elements Te, Au and Bi on an electron microprobe with specimens consisting of extremely thin films of the studied element deposited on thin, self-supporting, carbon layers. The relative variation of the ionization cross section was obtained by counting the number of characteristic X-rays from the considered element and shell, for varying incident electron energies, from the ionization energy up to 40 keV. Measured data were corrected to account for the energy-dependent spread of the electron beam within the active film and for the ionization due to the electrons backscattered from the carbon layer, using Monte Carlo simulation. Cross sections were evaluated in the Born approximation using an optical-data model with numerically evaluated dipole photoelectric cross sections. Calculated ionization cross section were converted to vacancy production cross sections, which can be directly compared with our experimental data.     

Electron propagator calculations show that alkyl substituents alter porphyrin ionization energies

The effect of alkyl substituents on the four lowest vertical ionization energies of porphyrins is determined with ab initio electron propagator calculations on porphine and octamethylporphyrin. With the use of the partial third-order approximation, predicted ionization energies are in close agreement with recent photoelectron spectra. These data and the associated Dyson orbitals, which describe changes in electronic structure that accompany photoionization, enable assignment of photoelectron spectra and determination of alkyl-induced shifts. Hyperconjugation is most evident in the Dyson orbitals associated with the third and fourth ionization energies of octamethylporphyrin and is least prominent in the Dyson orbital of the second ionization energy. There is a positive correlation between the shift in an ionization energy produced by alkyl substitution and the degree of hyperconjugation in the associated Dyson orbital. Alkyl substitutions, therefore, may be employed to adjust the ionization energies of porphyrins and, consequently, their reactivity patterns that depend on charge-transfer capabilities and disposition to electrophilic attack.     

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.     

On the Stability of Cyclophane Derivates Using a Molecular Fragmentation Method

A molecular fragmentation method is used to study the stability of cyclophane derivates by decomposing the molecular energy into the molecular strain and intramolecular interaction energies. The molecular strain energies obtained by utilising the fragmentation method are in good agreement with existing experimental data. The intramolecular interaction energies calculated as the difference between the supermolecular energy and the bonded fragment energies are repulsive in the cyclophanes studied. The nature of this interaction is studied for groups of systematically extended doubled layered paracyclophane systems using the random‐phase approximation (RPA), two recently developed extensions to the RPA and standard density functional theory (DFT) methods including dispersion corrections. Upon a systematic increase in conjugation the strongly repulsive intramolecular interaction energy reduces and thus leads to an increase in the stability. Finally, existing experimental and theoretical estimates of the molecular strain are compared with the results of this work.     

Correlations of the physicochemical parameters of azaindoles with their reactions

Using the AM1 semiempirical quantum method the enthalpies of formation, ionization energies, electron affinities, energy differences between highest occupied and lowest unoccupied orbitals, atomic charges, bond orders, and dipole moments have been calculated for 4-, 5-, 6-, and 7-azaindoles. A correlation has been built up between the calculated physicochemical parameters and the Hammett para-substituent and inductive constants. The 1H to 7H proton transfer in 7-azaindole has been quantitatively described. __________ Translated from Khimiya Geterotsiklicheskikh Soedinenii, No. 7, pp. 1062–1072, July, 2006.     

Calculated vertical ionization energies of the common α-amino acids in the gas phase and in solution

The vertical ionization energies of the low-lying conformers of the α-amino acids found in proteins have been calculated. Geometry optimizations were first performed at the B3LYP/6-311G(d,p) level of theory, and then reoptimized at the MP2/6-311G(d,p) level of theory. Vertical ionization energies were then computed by three methods, electron propagator in the partial third-order (P3) approximation, Outer-Valence-Green's Functions, and by evaluating the difference in the total energy between the cation radical and the neutral amino acid in the geometry of the neutral species. When available, the results are compared to the experimental vertical ionization energies. The vertical ionization energies calculated using the MP2/P3 method gave the best overall agreement with the experimental results. Next, the ionization energies in solution are calculated for the zwitterionic forms of the α-amino acids by using IEFPCM methods. To obtain the vertical ionization energy in solution, it is necessary to use the nonequilibrium polarizable continuum model (NEPCM), the results of which are reported here for the α-amino acids.