Spurious Oscillations Caused by Density Functional Approximations: Who is to Blame? Exchange or Correlation?

We analyze the varying susceptibilities of different density functional approximations (DFAs) to present spurious oscillations on the profiles of several vibrational properties. Among other problems, these spurious oscillations cause significant errors in harmonic and anharmonic IR and Raman frequencies and intensities. This work hinges on a judicious strategy to dissect the exchange and correlation components of DFAs and pinpoint the origins of these oscillations. We identify spurious oscillations in derivatives of all energy components with respect to nuclear displacements, including those energy terms that do not involve numerical integrations. These indirect spurious oscillations are attributed to suboptimal electron densities resulting from a self-consistent field procedure using a DFA that exhibits direct spurious oscillations. Direct oscillations stem from inaccurate numerical integration of the exchange and correlation energy density functionals. A thorough analysis of direct spurious oscillations reveals that only a handful of exchange and correlation components are insensitive to spurious oscillations, giving rise to three families of functionals, BH&H, LSDA, and BLYP. Among the functionals in these families, we encounter four widespread DFAs: BLYP, B3LYP, LC-BLYP, and CAM-B3LYP. Certain DFAs like PBE appear less sensitive to spurious oscillations due to compensatory cancellations between their energy components. Additionally, we found non-negligible but small oscillations in PBE and TPSS, which could be safely employed provided a sufficiently large integration grid is used in the calculations. These findings hint at the key components of current approximations to be improved and emphasize the necessity to develop accurate DFAs suitable for studying molecular spectroscopies.

All-Purpose Measure of Electron Correlation for Multireference Diagnostics

We present an analytical relationship between two natural orbital occupancy-based indices, I N D ¯ and INDmax, and two established electron correlation metrics: the leading term of a configuration interaction expansion, c0, and the D2 diagnostic. Numerical validation revealed that I N D ¯ and INDmax can effectively substitute for c0 and D2, respectively. These indices offer three distinct advantages: (i) they are universally applicable across all electronic structure methods, (ii) their interpretation is more intuitive, and (iii) they can be readily incorporated into the development of hybrid electronic structure methods. Additionally, we draw a distinction between correlation measures and correlation diagnostics, establishing MP2 and CCSD numerical thresholds for INDmax, which are to be used as a multireference diagnostic. Our findings further demonstrate that establishing thresholds for other electronic structure methods can be easily accomplished using small data sets.

How Reliable Are Modern Density Functional Approximations to Simulate Vibrational Spectroscopies?

We show that properties of molecules with low-frequency modes calculated with density functional approximations (DFAs) suffer from spurious oscillations along the nuclear displacement coordinate due to numerical integration errors. Occasionally, the problem can be alleviated using extensive integration grids that compromise the favorable cost-accuracy ratio of DFAs. Since spurious oscillations are difficult to predict or identify, DFAs are exposed to severe performance errors in IR and Raman intensities and frequencies or vibrational contributions to any molecular property. Using Fourier spectral analysis and digital signal processing techniques, we identify and quantify the error due to these oscillations for 45 widely used DFAs. LC-BLYP and BH&H are revealed as the only functionals showing robustness against the spurious oscillations of various energy, dipole moment, and polarizability derivatives with respect to a nuclear displacement coordinate. Given the ubiquitous nature of molecules with low-frequency modes, we warrant caution in using modern DFAs to simulate vibrational spectroscopies.

Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package

This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.

Role of Dispersion Interactions in Endohedral TM@(ZnS)12 Structures

II-VI semiconducting materials are gaining attention due to their optoelectronic properties. Moreover, the addition of transition metals, TMs, might give them magnetic properties. The location and distance of the TM are crucial in determining such magnetic properties. In this work, we focus on small hollow (ZnS)₁₂ nanoclusters doped with TMs. Because (ZnS)₁₂ is a cage-like spheroid, the cavity inside the structure allows for the design of endohedral compounds resembling those of C₆₀. Previous studies theoretically predicted that the first-row TM(ZnS)₁₂ endohedral compounds were thermodynamically unstable compared to the surface compounds, where the TM atom is located at the surface of the cluster. The transition states connecting both structure families were calculated, and the estimated lifetimes of these compounds were predicted to be markedly small. However, in such works dispersion effects were not taken into account. Here, in order to check for the influence of dispersion on the possible stabilization of the desired TM(ZnS)₁₂ endohedrally doped clusters, several functionals are tested and compare to MP2. It is found that the dispersion effects play a very important role in determining the location of the metals, especially in those TMs with the 4s3d shell half-filled or completely filled. In addition, a complete family of TM doped (ZnS)₁₂ nanoclusters is explored using ab initio molecular dynamics simulations and local minima optimizations that could guide the experimental synthesis of such compounds. From the magnetic point of view, the Cr(⁷S)@(ZnS)₁₂ compound is the most interesting case, since the endohedral isomer is predicted to be the global minimum. Moreover, molecular dynamics simulations show that when the Cr atom is located at the surface of the cluster, it spontaneously migrates toward the center of the cavity at room temperature.

How Many Electrons Does a Molecular Electride Hold?

Electrides are very peculiar ionic compounds where electrons occupy the anionic positions. In a crystal lattice, these isolated electrons often form channels or surfaces, furnishing electrides with many traits with promising technological applications. Despite their huge potential, thus far, only a few stable electrides have been produced because of the intricate synthesis they entail. Due to the difficulty in assessing the presence of isolated electrons, the characterization of electrides also poses some serious challenges. In fact, their properties are expected to depend on the arrangement of these electrons in the molecule. Among the criteria that we can use to characterize electrides, the presence of a non-nuclear attractor (NNA) of the electron density is both the rarest and the most salient feature. Therefore, a correct description of the NNA is crucial to determine the properties of electrides. In this paper, we analyze the NNA and the surrounding region of nine molecular electrides to determine the number of isolated electrons held in the electride. We have seen that the correct description of a molecular electride hinges on the electronic structure method employed for the analyses. In particular, one should employ a basis set with sufficient flexibility to describe the region close to the NNA and a density functional approximation that does not suffer from large delocalization errors. Finally, we have classified these nine molecular electrides according to the most likely number of electrons that we can find in the NNA. We believe this classification highlights the strength of the electride character and will prove useful in designing new electrides.

Singling Out Dynamic and Nondynamic Correlation

The correlation part of the pair density is separated into two components, one of them being predominant at short electronic ranges and the other at long ranges. The analysis of the intracular part of these components permits to classify molecular systems according to the prevailing correlation: dynamic or nondynamic. The study of the long-range asymptotics reveals the key component of the pair density that is responsible for the description of London dispersion forces and a universal decay with the interelectronic distance. The natural range-separation, the identification of the dispersion forces, and the kind of predominant correlation type that arise from this analysis are expected to be important assets in the development of new electronic structure methods in wave function, density, and reduced density-matrix functional theories.

Nine questions on energy decomposition analysis

The paper collects the answers of the authors to the following questions: Is the lack of precision in the definition of many chemical concepts one of the reasons for the coexistence of many partition schemes? Does the adoption of a given partition scheme imply a set of more precise definitions of the underlying chemical concepts? How can one use the results of a partition scheme to improve the clarity of definitions of concepts? Are partition schemes subject to scientific Darwinism? If so, what is the influence of a community's sociological pressure in the "natural selection" process? To what extent does/can/should investigated systems influence the choice of a particular partition scheme? Do we need more focused chemical validation of Energy Decomposition Analysis (EDA) methodology and descriptors/terms in general? Is there any interest in developing common benchmarks and test sets for cross-validation of methods? Is it possible to contemplate a unified partition scheme (let us call it the "standard model" of partitioning), that is proper for all applications in chemistry, in the foreseeable future or even in principle? In the end, science is about experiments and the real world. Can one, therefore, use any experiment or experimental data be used to favor one partition scheme over another? © 2019 Wiley Periodicals, Inc.

The Coulomb Hole of the Ne Atom

We analyze the Coulomb hole of Ne from highly-accurate CISD wave functions obtained from optimized even-tempered basis sets. Using a two-fold extrapolation procedure we obtain highly accurate results that recover 97 % of the correlation energy. We confirm the existence of a shoulder in the short-range region of the Coulomb hole of the Ne atom, which is due to an internal reorganization of the K-shell caused by electron correlation of the core electrons. The feature is very sensitive to the quality of the basis set in the core region and it is not exclusive to Ne, being also present in most of second-row atoms, thus confirming that it is due to K-shell correlation effects.

Comprehensive benchmarking of density matrix functional approximations

The energy usually serves as a yardstick in assessing the performance of approximate methods in computational chemistry. After all, these methods are mostly used for the calculation of the electronic energy of chemical systems. However, computational methods should be also aimed at reproducing other properties, such strategy leading to more robust approximations with a wider range of applicability. In this study, we suggest a battery of ten tests with the aim to analyze density matrix functional approximations (DMFAs), including several properties that the exact functional should satisfy. The tests are performed on a model system with varying electron correlation, carrying a very small computational effort. Our results not only put forward a complete and exhaustive benchmark test for DMFAs, currently lacking, but also reveal serious deficiencies of existing approximations that lead to important clues in the construction of more robust DMFAs.

Electron correlation effects in third-order densities

The electronic energy of a system of fermions can be obtained from the second-order reduced density matrix through the contracted Schrödinger equation or its anti-Hermitian counterpart. Both energy expressions depend on the third-order reduced density matrix (3-RDM) which is usually approximated from lower-order densities. The accuracy of these methods depends critically on the set of N-representability conditions enforced in the calculation and the quality of the approximate 3-RDM. There are no benchmark studies including most 3-RDM approximations and, thus far, no assessment of the deterioration of the approximations with correlation effects has been performed. In this paper we introduce a series of tests to assess the performance of 3-RDM approximations in a model system with varying electron correlation effects, the three-electron harmonium atom. The results of this work put forward several limitations of the currently most used 3-RDM approximations for systems with important electron correlation effects.

Separation of dynamic and nondynamic correlation

The account of electron correlation and its efficient separation into dynamic and nondynamic parts plays a key role in the development of computational methods. In this paper we suggest a physically-sound matrix formulation to split electron correlation into dynamic and nondynamic parts using the two-particle cumulant matrix and a measure of the deviation from idempotency of the first-order density matrix. These matrices are applied to a two-electron model, giving rise to a simplified electron correlation index that (i) depends only on natural orbitals and their occupancies, (ii) can be straightforwardly decomposed into orbital contributions and (iii) splits into dynamic and nondynamic correlation parts that (iv) admit a local version. These expressions are shown to account for dynamic and nondynamic correlation in a variety of systems containing different electron correlation regimes, thus providing the first separation of dynamic and nondynamic correlation using solely natural orbital occupancies.

H4: A challenging system for natural orbital functional approximations

The correct description of nondynamic correlation by electronic structure methods not belonging to the multireference family is a challenging issue. The transition of D(2h) to D(4h) symmetry in H4 molecule is among the most simple archetypal examples to illustrate the consequences of missing nondynamic correlation effects. The resurgence of interest in density matrix functional methods has brought several new methods including the family of Piris Natural Orbital Functionals (PNOF). In this work, we compare PNOF5 and PNOF6, which include nondynamic electron correlation effects to some extent, with other standard ab initio methods in the H4 D(4h)/D(2h) potential energy surface (PES). Thus far, the wrongful behavior of single-reference methods at the D(2h)-D(4h) transition of H4 has been attributed to wrong account of nondynamic correlation effects, whereas in geminal-based approaches, it has been assigned to a wrong coupling of spins and the localized nature of the orbitals. We will show that actually interpair nondynamic correlation is the key to a cusp-free qualitatively correct description of H4 PES. By introducing interpair nondynamic correlation, PNOF6 is shown to avoid cusps and provide the correct smooth PES features at distances close to the equilibrium, total and local spin properties along with the correct electron delocalization, as reflected by natural orbitals and multicenter delocalization indices.

A computational study of radical initiated protein backbone homolytic dissociation on all natural amino acids

Hydroxyl radical (˙OH) is known to be one of the most reactive species. The attack of this radical onto the backbone of all natural amino acids is investigated.