Search

Are the vertical ionization energies from a bound electronic system, initially in its ground state, equal to minus the corresponding exact Kohn-Sham orbital energies of density functional theory (DFT)? This is known to be true for the first or lowest vertical ionization energy. We show that the correction from time-dependent DFT arises from the continuum and need not vanish. Recent work compared the experimental photoemission thresholds of the molecules Cu 2 O - , CuO - , CuO 2 - , and CuO 3 - with minus the corresponding orbital energies from a generalized gradient approximation (GGA) and its global and range-separated hybrids with exact exchange, finding striking differences which were attributed to self-interaction error, strong correlation, or both. Here, we extend that work to include the local spin density approximation (LSDA), its Perdew-Zunger self-interaction correction with Fermi-Löwdin localized orbitals (LSDA-SIC), a quasi-self-consistent locally scaled-down version of LSDA-SIC (QLSIC), and the Quantum Theory Project QTP02 range-separated hybrid functional, all but LSDA implemented in a generalized Kohn-Sham approach. QTP02 impressively yields a near equality for many sp-bonded molecules. However, for the copper oxide anions studied here, none of the tested methods reproduces the experimental photoemission thresholds.

We analyzed various possibilities to improve upon the SCAN meta-generalized gradient approximation density functional obeying all known properties of the exact functional that can be satisfied at this level of approximation. We examined the necessity of locally satisfying a strongly tightened lower bound for the exchange energy density in single-orbital regions, the nature of the error cancellation between the exchange and correlation parts in two-electron regions, and the effect of the fourth-order term in the gradient expansion of the correlation energy density. We have concluded that the functional can be modified to separately reproduce the exchange and correlation energies of the helium atom by locally releasing the strongly tightened lower bound for the exchange energy density in single-orbital regions, but this leads to an unbalanced improvement in the single-orbital electron densities. Therefore, we decided to keep the F X ≤ 1.174 exact condition for any single-orbital density, where F X is the exchange enhancement factor. However, we observed a general improvement in the single-orbital electron densities by revising the correlation functional form to follow the second-order gradient expansion in a wider range. Our new revSCAN functional provides more-accurate atomization energies for the systems with multireference character, compared to the SCAN functional. The nonlocal VV10 dispersion-corrected revSCAN functional yields more-accurate noncovalent interaction energies than the VV10-corrected SCAN functional. Furthermore, its global hybrid version with 25% of exact exchange, called revSCAN0, generally performs better than the similar SCAN0 for reaction barrier heights. Here, we also analyzed the possibility of the construction of a local hybrid from the SCAN exchange and a specific locally bounded nonconventional exact exchange energy density. We predict compatibility problems since this nonconventional exact exchange energy density does not really obey the strongly tightened lower bound for the exchange energy density in single-orbital regions.

Since its formal introduction, density functional theory has achieved many successes in the fields of molecular and solid-state chemistry. According to its central theorems, the ground state of a many-electron system is fully described by its electron density, and the exact functional minimizes the energy at the exact electron density. For many years of density functional development, it was assumed that the improvements in the energy are accompanied by the improvements in the density, and the approximations approach the exact functional. In a recent analysis ( Medvedev et al. Science 2017 , 355 , 49 - 52 .), it has been pointed out for 14 first row (Be-Ne) atoms and cations with 2, 4, or 10 electrons that the nowadays popular flexible but physically less rigorous approximate density functionals may provide large errors in the calculated electron densities despite the accurate energies. Although far-reaching conclusions have been drawn in this work, the methodology used by the authors may need improvements. Most importantly, their benchmark set was biased toward small atomic cations with compressed, high electron densities. In our paper, we construct a molecular test set with chemically relevant densities and analyze the performance of several density functional approximations including the less-investigated double hybrids. We apply an intensive error measure for the density, its gradient, and its Laplacian and examine how the errors in the density propagate into the semilocal exchange-correlation energy. While we have confirmed the broad conclusions of Medvedev et al., our different way of analyzing the data has led to conclusions that differ in detail. Finally, seeking for a rationale behind the global hybrid or double hybrid methods from the density's point of view, we also analyze the role of the exact exchange and second-order perturbative correlation mixing in PBE-based global hybrid and double hybrid functional forms.

Recently, we have constructed a dual-hybrid direct random phase approximation method, called dRPA75, and demonstrated its good performance on reaction energies, barrier heights, and noncovalent interactions of main-group elements. However, this method has also shown significant but quite systematic errors in the computed atomization energies. In this paper, we suggest a constrained spin-component scaling formalism for the dRPA75 method (SCS-dRPA75) in order to overcome the large error in the computed atomization energies, preserving the good performance of this method on spin-unpolarized systems at the same time. The SCS-dRPA75 method with the aug-cc-pVTZ basis set results in an average error lower than 1.5 kcal mol -1 for the entire n-homodesmotic hierarchy of hydrocarbon reactions (RC0-RC5 test sets). The overall performance of this method is better than the related direct random phase approximation-based double-hybrid PWRB95 method on open-shell systems of main-group elements (from the GMTKN30 database) and comparable to the best O(N 4 )-scaling opposite-spin second-order perturbation theory-based double-hybrid methods like PWPB95-D3 and to the O(N 5 )-scaling RPAX2@PBEx method, which also includes exchange interactions. Furthermore, it gives well-balanced performance on many types of barrier heights similarly to the best O(N 5 )-scaling second-order perturbation theory-based or spin-component scaled second-order perturbation theory-based double-hybrid methods such as XYG3 or DSD-PBEhB95. Finally, we show that the SCS-dRPA75 method has reduced self-interaction and delocalization errors compared to the parent dRPA75 method and a slightly smaller static correlation error than the related PWRB95 method.

Without extensive fitting, accurate prediction of transition metal chemistry is a challenge for semilocal and hybrid density funcitonals. The Random Phase Approximation (RPA) has been shown to yield superior results to semilocal functionals for main group thermochemistry, but much less is known about its performance for transition metals. We have therefore analyzed the behavior of reaction energies, barrier heights, and ligand dissociation energies obtained with RPA and compare our results to several semilocal and hybrid functionals. Particular attention is paid to the reference determinant dependence of RPA. We find that typically the results do not vary much between semilocal or hybrid functionals as a reference, as long as the fraction of exact exchange (EXX) mixing in the hybrid functional is small. For large fractions of EXX mixing, however, the Hartree-Fock-like nature of the determinant can severely degrade the performance. Overall, RPA systematically reduces the errors of semilocal functionals and delivers excellent performance from a single reference determinant for inherently multireference reactions. The behavior of dual hybrids that combine RPA correlation with a hybrid exchange energy was also explored, but ultimately did not lead to a systematic improvement compared to traditional RPA for these systems. We rationalize this conclusion by decomposing the contributions to the reaction energies, and briefly discuss the possible implications for double-hybrid functionals based on RPA. The correlation between EXX mixing and spin-symmetry breaking is also discussed.

In water clusters, there is a delicate balance of van der Waals interactions and hydrogen bonds. Although semilocal and nonlocal density functional approximations have been recently routinely applied to water in various phases, the accurate description of hydrogen bonds remains a challenge. The most popular density functional approaches fail to predict the correct ordering of the energies of water clusters. To illustrate the required accuracy, the CCSD(T) complete basis set extrapolated dissociation energy difference between the two lowest energy hexamer structures is 0.06 kcal mol(-1) per monomer. In this work, we assessed interaction energies in neutral and ionic water clusters with various density functionals with or without van der Waals correction. Generally, van der Waals approximations play a significant role in clusters with increasing size, while hybrid functionals improve the description of hydrogen bonds. Despite these general trends, none of the tested density functional approximations with or without van der Waals correction and exact exchange mixing can lead to a uniform performance for neutral and ionic water clusters. The recently constructed dual-hybrid dRPA75 approximation is a successful combination of exact and semilocal exchange, and nonlocal correlation in its energy, while utilizing a high fraction of exact exchange. We have shown that the dRPA75 method has a systematic error, which can be efficiently compensated for by the aug-cc-pVTZ basis set for small- and medium-sized water clusters.

The direct random phase approximation (dRPA) combined with Kohn-Sham reference orbitals is among the most promising tools in computational chemistry and applicable in many areas of chemistry and physics. The reason for this is that it scales as N(4) with the system size, which is a considerable advantage over the accurate ab initio wave function methods like standard coupled-cluster. dRPA also yields a considerably more accurate description of thermodynamic and electronic properties than standard density-functional theory methods. It is also able to describe strong static electron correlation effects even in large systems with a small or vanishing band gap missed by common single-reference methods. However, dRPA has several flaws due to its self-correlation error. In order to obtain accurate and precise reaction energies, barriers and noncovalent intra- and intermolecular interactions, we construct a new dual-hybrid dRPA (hybridization of exact and semilocal exchange in both the energy and the orbitals) and test the performance of this new functional on isogyric, isodesmic, hypohomodesmotic, homodesmotic, and hyperhomodesmotic reaction classes. We also use a test set of 14 Diels-Alder reactions, six atomization energies (AE6), 38 hydrocarbon atomization energies, and 100 reaction barrier heights (DBH24, HT-BH38, and NHT-BH38). For noncovalent complexes, we use the NCCE31 and S22 test sets. To test the intramolecular interactions, we use a set of alkane, cysteine, phenylalanine-glycine-glycine tripeptide, and monosaccharide conformers. We also discuss the delocalization and static correlation errors. We show that a universally accurate description of chemical properties can be provided by a large, 75% exact exchange mixing both in the calculation of the reference orbitals and the final energy.

The direct random phase approximation (dRPA) is a promising way to obtain improvements upon the standard semilocal density functional results in many aspects of computational chemistry. In this paper, we address the slow convergence of the calculated dRPA correlation energy with the increase of the quality and size of the popular Gaussian-type Dunning's correlation consistent aug-cc-pVXZ split valence atomic basis set family. The cardinal number X controls the size of the basis set, and we use X = 3-6 in this study. It is known that even the very expensive X = 6 basis sets lead to large errors for the dRPA correlation energy, and thus complete basis set extrapolation is necessary. We study the basis set convergence of the dRPA correlation energies on a set of 65 hydrocarbon isomers from CH4 to C6H6. We calculate the iterative density fitted dRPA correlation energies using an efficient algorithm based on the CC-like form of the equations using the self-consistent HF orbitals. We test the popular inverse cubic, the optimized exponential, and inverse power formulas for complete basis set extrapolation. We have found that the optimized inverse power based extrapolation delivers the best energies. Further analysis showed that the optimal exponent depends on the molecular structure, and the most efficient two-point energy extrapolations that use X = 3 and 4 can be improved considerably by considering the atomic composition and hybridization states of the atoms in the molecules. Our results also show that the optimized exponents that yield accurate X = 3 and 4 extrapolated dRPA energies for atoms or small molecules might be inaccurate for larger molecules.

We assess the performance of the semilocal PBE functional; its global hybrid variants; the highly parametrized empirical M06-2X and M08-SO; the range separated rCAM-B3LYP and MCY3; the atom-pairwise or nonlocal dispersion corrected semilocal PBE and TPSS; the dispersion corrected range-separated ωB97X-D; the dispersion corrected double hybrids such as PWPB95-D3; the direct random phase approximation, dRPA, with Hartree-Fock, Perdew-Burke-Ernzerhof, and Perdew-Burke-Ernzerhof hybrid reference orbitals and the RPAX2 method based on a Perdew-Burke-Ernzerhof exchange reference orbitals for the Diels-Alder, DARC; and self-interaction error sensitive, SIE11, reaction energy test sets with large, augmented correlation consistent valence basis sets. The dRPA energies for the DARC test set are extrapolated to the complete basis set limit. CCSD(T)/CBS energies were used as a reference. The standard global hybrid functionals show general improvements over the typical endothermic energy error of semilocal functionals, but despite the increased accuracy the precision of the methods increases only slightly, and thus all reaction energies are simply shifted into the exothermic direction. Dispersion corrections give mixed results for the DARC test set. Vydrov-Van Voorhis 10 correction to the reaction energies gives superior quality results compared to the too-small D3 correction. Functionals parametrized for energies of noncovalent interactions like M08-SO give reasonable results without any dispersion correction. The dRPA method that seamlessly and theoretically correctly includes noncovalent interaction energies gives excellent results with properly chosen reference orbitals. As the results for the SIE11 test set and H2(+) dissociation show that the dRPA methods suffer from delocalization error, good reaction energies for the DARC test set from a given method do not prove that the method is free from delocalization error. The RPAX2 method shows good performance for the DARC, the SIE11 test sets, and for the H2(+) and H2 potential energy curves showing no one-electron self-interaction error and reduced static correlation errors at the same time. We also suggest simplified DARC6 and SIE9 test sets for future benchmarking.

A correct description of the anion-π interaction is essential for the design of selective anion receptors and channels and important for advances in the field of supramolecular chemistry. However, it is challenging to do accurate, precise, and efficient calculations of this interaction, which are lacking in the literature. In this article, by testing sets of 20 binary anion-π complexes of fluoride, chloride, bromide, nitrate, or carbonate ions with hexafluorobenzene, 1,3,5-trifluorobenzene, 2,4,6-trifluoro-1,3,5-triazine, or 1,3,5-triazine and 30 ternary π-anion-π' sandwich complexes composed from the same monomers, we suggest domain-based local-pair natural orbital coupled cluster energies extrapolated to the complete basis-set limit as reference values. We give a detailed explanation of the origin of anion-π interactions, using the permanent quadrupole moments, static dipole polarizabilities, and electrostatic potential maps. We use symmetry-adapted perturbation theory (SAPT) to calculate the components of the anion-π interaction energies. We examine the performance of the direct random phase approximation (dRPA), the second-order screened exchange (SOSEX), local-pair natural-orbital (LPNO) coupled electron pair approximation (CEPA), and several dispersion-corrected density functionals (including generalized gradient approximation (GGA), meta-GGA, and double hybrid density functional). The LPNO-CEPA/1 results show the best agreement with the reference results. The dRPA method is only slightly less accurate and precise than the LPNO-CEPA/1, but it is considerably more efficient (6-17 times faster) for the binary complexes studied in this paper. For 30 ternary π-anion-π' sandwich complexes, we give dRPA interaction energies as reference values. The double hybrid functionals are much more efficient but less accurate and precise than dRPA. The dispersion-corrected double hybrid PWPB95-D3(BJ) and B2PLYP-D3(BJ) functionals perform better than the GGA and meta-GGA functionals for the present test set.

Computationally efficient semilocal approximations of density functional theory at the level of the local spin density approximation (LSDA) or generalized gradient approximation (GGA) poorly describe weak interactions. We show improved descriptions for weak bonds (without loss of accuracy for strong ones) from a newly developed semilocal meta-GGA (MGGA), by applying it to molecules, surfaces, and solids. We argue that this improvement comes from using the right MGGA dimensionless ingredient to recognize all types of orbital overlap.

Among the computationally efficient semilocal density functionals for the exchange-correlation energy, meta-generalized-gradient approximations (meta-GGAs) are potentially the most accurate. Here, we assess the performance of three new meta-GGAs (revised Tao-Perdew-Staroverov-Scuseria or revTPSS, regularized revTPSS or regTPSS, and meta-GGA made simple or MGGA_MS), within and beyond their "comfort zones," on Grimme's big test set of main-group molecular energetics (thermochemistry, kinetics, and noncovalent interactions). We compare them against the standard Perdew-Burke-Ernzerhof (PBE) GGA, TPSS, and Minnesota M06L meta-GGAs, and Becke-3-Lee-Yang-Parr (B3LYP) hybrid of GGA with exact exchange. The overall performance of these three new meta-GGA functionals is similar. However, dramatic differences occur for different test sets. For example, M06L and MGGA_MS perform best for the test sets that contain noncovalent interactions. For the 14 Diels-Alder reaction energies in the "difficult" DARC subset, the mean absolute error ranges from 3 kcal mol(-1) (MGGA_MS) to 15 kcal mol(-1) (B3LYP), while for some other reaction subsets the order of accuracy is reversed; more generally, the tested new semilocal functionals outperform the standard B3LYP for ring reactions. Some overall improvement is found from long-range dispersion corrections for revTPSS and regTPSS but not for MGGA_MS. Formal and universality criteria for the functionals are also discussed.

The van der Waals coefficients between quasispherical nanostructures can be modeled accurately and analytically by those of classical solid spheres (for nanoclusters) or spherical shells (for fullerenes) of uniform valence electron density, with the true static dipole polarizability. Here, we derive analytically and confirm numerically from this model the size dependencies of the van der Waals coefficients of all orders, showing, for example, that the asymptotic dependence for C(6) is the expected n(2) for pairs of nanoclusters A(n)-A(n), each containing n atoms, but n(2.75) for pairs of single-walled fullerenes C(n)-C(n). Large fullerenes are argued to have much larger polarizabilities and dispersion coefficients than those predicted by either the standard atom pair-potential model or widely used nonlocal van der Waals correlation energy functionals.

Fullerene molecules such as C(60) are large nearly spherical shells of carbon atoms. Pairs of such molecules have a strong long-range van der Waals attraction that can produce scattering or binding into molecular crystals. A simplified classical-electrodynamics model for a fullerene is a spherical metal shell, with uniform electron density confined between outer and inner radii (just as a simplified model for a nearly spherical metallic nanocluster is a solid metal sphere or filled shell). For the spherical-shell model, the exact dynamic multipole polarizabilities are all known analytically. From them, we can derive exact analytic expressions for the van der Waals coefficients of all orders between two spherical metal shells. The shells can be identical or different, and hollow or filled. To connect the model to a real fullerene, we input the static dipole polarizability, valence electron number and estimated shell thickness t of the real molecule. Our prediction for the leading van der Waals coefficient C(6) between two C(60) molecules ((1.30 ± 0.22) × 10(5) hartree bohr(6)) agrees well with a prediction for the real molecule from time-dependent density functional theory. Our prediction is remarkably insensitive to t. Future work might include the prediction of higher-order (e.g. C(8) and C(10)) coefficients for C(60), applications to other fullerenes or nearly spherical metal clusters, etc. We also make general observations about the van der Waals coefficients.

We have improved the revised Tao-Perdew-Staroverov-Scuseria (revTPSS) meta-generalized gradient approximation (GGA) in order to remove the order of limits anomaly in its exchange energy. The revTPSS meta-GGA recovers the second-order gradient expansion for a wide range of densities and therefore provides excellent atomization energies and lattice constants. For other properties of materials, however, even the revTPSS does not give the desired accuracy. The revTPSS does not perform as well as expected for the energy differences between different geometries for the same molecular formula and for the related nonbarrier height chemical reaction energies. The same order of limits problem might lead to inaccurate energy differences between different crystal structures and to inaccurate cohesive energies of insulating solids. Here we show a possible way to remove the order of limits anomaly with a weighted difference of the revTPSS exchange between the slowly varying and iso-orbitals (one- or two-electron) limits. We show that the new regularized (regTPSS) gives atomization energies comparable to revTPSS and preserves the accurate lattice constants as well. For other properties, the regTPSS gives at least the same performance as the revTPSS or TPSS meta-GGAs.

Adiabatic Φ/ψ maps for cellobiose were prepared with B3LYP density functional theory. A mixed basis set was used for minimization, followed with 6-31+G(d) single-point calculations, with and without SMD continuum solvation. Different arrangements of the exocyclic groups (38 starting geometries) were considered for each Φ/ψ point. The vacuum calculations agreed with earlier computational and experimental results on the preferred gas phase conformation (anti-Φ(H), syn-ψ(H)), and the results from the solvated calculations were consistent with the (syn Φ(H)/ψ(H) conformations from condensed phases (crystals or solutions). Results from related studies were compared, and there is substantial dependence on the solvation model as well as arrangements of exocyclic groups. New stabilizing interactions were revealed by Atoms-In-Molecules theory., (Published by Elsevier Ltd.)

Correlated ab initio wave function calculations have been performed, using nonrelativistic frozen core MP2 complete basis set extrapolation model chemistry. The calculations have been made for three test sets of gas-phase saccharide conformations to provide reference values for their relative energies. The remaining correlation effects are estimated from frozen core coupled-cluster singles and doubles [CCSD(T)] calculations. The test sets consist of 15 conformers of α- and β-d-allopyranose, 15 of 3,6-anhydro-4-O-methyl-d-galactitol, and four of β-d-glucopyranose. For each set, conformational energies varied by about 7 kcal/mol. These benchmark quality relative conformational energies are used to re-evaluate the performance of the best density functional methods for conformational analyses of saccharides. Our results show that the B3PW91 and PBE0 relative energies are systematically better than the B3LYP and M05-2X results. Overall, the functionals based on the exact constraints perform better for the relative energies of monosaccharide conformers than the empirically fitted functionals.

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.

There is current interest in the random phase approximation (RPA), a "fifth-rung" density functional for the exchange-correlation energy. RPA has full exact exchange and constructs the correlation with the help of the unoccupied Kohn-Sham orbitals. In many cases (uniform electron gas, jellium surface, and free atom), the correction to RPA is a short-ranged effect that is captured by a local spin density approximation (LSDA) or a generalized gradient approximation (GGA). Nonempirical density functionals for the correction to RPA were constructed earlier at the LSDA and GGA levels (RPA+), but they are constructed here at the fully nonlocal level (RPA++), using the van der Waals density functional (vdW-DF) of Langreth, Lundqvist, and collaborators. While they make important and helpful corrections to RPA total and ionization energies of free atoms, they correct the RPA atomization energies of molecules by only about 1 kcal/mol. Thus, it is puzzling that RPA atomization energies are, on average, about 10 kcal/mol lower than those of accurate values from experiment. We find here that a hybrid of 50% Perdew-Burke-Ernzerhof GGA with 50% RPA+ yields atomization energies much more accurate than either one does alone. This suggests a solution to the puzzle: While the proper correction to RPA is short-ranged in some systems, its contribution to the correlation hole can spread out in a molecule with multiple atomic centers, canceling part of the spread of the exact exchange hole (more so than in RPA or RPA+), making the true exchange-correlation hole more localized than in RPA or RPA+. This effect is not captured even by the vdW-DF nonlocality, but it requires the different kind of full nonlocality present in a hybrid functional.

Eighteen empirical force fields and the semi-empirical quantum method PM3CARB-1 were compared for studying beta-cellobiose, alpha-maltose, and alpha-galabiose [alpha-D-Galp-(1-->4)-alpha-D-Galp]. For each disaccharide, the energies of 54 conformers with differing hydroxymethyl, hydroxyl, and glycosidic linkage orientations were minimized by the different methods, some at two dielectric constants. By comparing these results and the available crystal structure data and/or higher level density functional theory results, it was concluded that the newer parameterizations for force fields (GROMOS, GLYCAM06, OPLS-2005 and CSFF) give results that are reasonably similar to each other, whereas the older parameterizations for Amber, CHARMM or OPLS were more divergent. However, MM3, an older force field, gave energy and geometry values comparable to those of the newer parameterizations, but with less sensitivity to dielectric constant values. These systems worked better than MM2 variants, which were still acceptable. PM3CARB-1 also gave adequate results in terms of linkage and exocyclic torsion angles. GROMOS, GLYCAM06, and MM3 appear to be the best choices, closely followed by MM4, CSFF, and OPLS-2005. With GLYCAM06 and to a lesser extent, CSFF, and OPLS-2005, a number of the conformers that were stable with MM3 changed to other forms.

Semilocal density functionals for the exchange-correlation energy are needed for large electronic systems. The Tao-Perdew-Staroverov-Scuseria (TPSS) meta-generalized gradient approximation (meta-GGA) is semilocal and usefully accurate, but predicts too-long lattice constants. Recent "GGA's for solids" yield good lattice constants but poor atomization energies of molecules. We show that the construction principle for one of them (restoring the density gradient expansion for exchange over a wide range of densities) can be used to construct a "revised TPSS" meta-GGA with accurate lattice constants, surface energies, and atomization energies for ordinary matter.

Some fundamental issues in ground-state density functional theory are discussed without equations: (1) The standard Hohenberg-Kohn and Kohn-Sham theorems were proven for a Hamiltonian that is not quite exact for real atoms, molecules, and solids. (2) The density functional for the exchange-correlation energy, which must be approximated, arises from the tendency of electrons to avoid one another as they move through the electron density. (3) In the absence of a magnetic field, either spin densities or total electron density can be used, although the former choice is better for approximations. (4) "Spin contamination" of the determinant of Kohn-Sham orbitals for an open-shell system is not wrong but right. (5) Only to the extent that symmetries of the interacting wave function are reflected in the spin densities should those symmetries be respected by the Kohn-Sham noninteracting or determinantal wave function. Functionals below the highest level of approximations should however sometimes break even those symmetries, for good physical reasons. (6) Simple and commonly used semilocal (lower-level) approximations for the exchange-correlation energy as a functional of the density can be accurate for closed systems near equilibrium and yet fail for open systems of fluctuating electron number. (7) The exact Kohn-Sham noninteracting state need not be a single determinant, but common approximations can fail when it is not. (8) Over an open system of fluctuating electron number, connected to another such system by stretched bonds, semilocal approximations make the exchange-correlation energy and hole-density sum rule too negative. (9) The gap in the exact Kohn-Sham band structure of a crystal underestimates the real fundamental gap but may approximate the first exciton energy in the large-gap limit. (10) Density functional theory is not really a mean-field theory, although it looks like one. The exact functional includes strong correlation, and semilocal approximations often overestimate the strength of static correlation through their semilocal exchange contributions. (11) Only under rare conditions can excited states arise directly from a ground-state theory.

Correlated ab initio wave function calculations using MP2/aug-cc-pVTZ model chemistry have been performed for three test sets of gas phase saccharide conformations to provide reference values for their relative energies. The test sets consist of 15 conformers of α- and β-d-allopyranose, 15 of 3,6-anhydro-4-O-methyl-d-galactitol, and four of β-d-glucopyranose. For each set, conformational energies varied by about 7 kcal/mol. Results obtained with the Hartree-Fock method, with pure density functional approximations (DFAs) like LSDA, PBEsol, PBE, and TPSS and with hybrid DFAs like B3PW91, B3LYP, PBEh, and M05-2X, were then compared to the reference and local MP2 relative energies. Basis sets included 6-31G*, 6-31G**, 6-31+G*, 6-31+G**, 6-311+G**, 6-311++G**, cc-pVTZ(-f), cc-pVTZ, and aug-cc-pVTZ(-f). The smallest basis set that gives good DFA relative energies is 6-31+G**, and more converged results can be obtained with 6-311+G**. The optimized geometries obtained from a smaller basis set, 6-31+G*, were useful for subsequent single point energy calculations with larger basis sets. The best agreement with MP2 was shown by M05-2X, but only when using a dense DFT grid. The popular B3LYP functional is not the best for saccharide conformational studies. The B3PW91 functional gives systematically better results, but other hybrid functionals like PBEh or TPSSh are even better. Overall, the nonempirical PBE GGA and TPSS meta-GGA functionals also performed better than B3LYP.

A new, regularized gradient expansion (RGE) approximation density functional (i.e., a generalized gradient approximation or GGA that recovers the second-order gradient expansion for exchange in the slowly varying limit) was designed in an attempt to obtain good solid-state and molecular properties at the same time from a single GGA. We assess the performance of this functional for molecular atomization energies, solid lattice constants, and jellium surface energies. We compare the performance of this functional to the modified Perdew-Burke-Ernzerhof generalized gradient approximation (PBEsol GGA), the original PBE GGA, and the Tao-Perdew-Staroverov-Scuseria (TPSS) meta-GGA.

Serious and systematic errors with popular density functionals occur for isodesmic stabilization energies of n-alkanes, isomerization, and dimerization energies of hydrocarbons and geometries of sterically overcrowded aromatic systems. These functionals are too biased toward the correct description of free atoms. Changing two parameters within the Perdew-Burke-Ernzerhof approximation leads to a new nonempirical functional, PBEsol, that improves the description of large organic systems.

Popular modern generalized gradient approximations are biased toward the description of free-atom energies. Restoration of the first-principles gradient expansion for exchange over a wide range of density gradients eliminates this bias. We introduce a revised Perdew-Burke-Ernzerhof generalized gradient approximation that improves equilibrium properties of densely packed solids and their surfaces.

The Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation for the exchange-correlation energy functional has two nonempirical constructions, based on satisfaction of universal exact constraints on the hole density or on the energy. We show here that, by identifying one possible free parameter in exchange and a second in correlation, we can continue to satisfy these constraints while diminishing the gradient dependence almost to zero (i.e., almost recovering the local spin density approximation or LSDA). This points out the important role played by the Perdew-Wang 1991 nonempirical hole construction in shaping PBE and later constructions. Only the undiminished PBE is good for atoms and molecules, for reasons we present, but a somewhat diminished PBE could be useful for solids; in particular, the surface energies of solids could be improved. Even for atoms and molecules, a strongly diminished PBE works well when combined with a scaled-down self-interaction correction (although perhaps not significantly better than LSDA). This shows that the undiminished gradient dependence of PBE and related functionals works somewhat like a scaled-down self-interaction correction to LSDA.

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.

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.

Semilocal density functional approximations (DFAs) for the exchange-correlation energy suffer from self-interaction error, which is believed to be the cause of many of the failures of common DFAs, such as poor description of charge transfer and transition states of chemical reactions. The standard self-interaction correction (SIC) of Perdew and Zunger mends some of these failures but spoils such essential properties as thermochemistry and equilibrium bond lengths. The Perdew-Zunger SIC seems to overcorrect many-electron systems. In this paper, we propose a modified SIC, which is scaled down in many-electron regions. The new SIC has an improved performance for many molecular properties, including total energies, atomization energies, barrier heights of chemical reactions, ionization potentials, electron affinities, and bond lengths. The local spin-density approximation (LSDA) benefits from SIC more than higher-level functionals do. The scaled-down SIC has only one adjustable parameter. Rationalization of the optimal value of this parameter enables us to construct an almost-nonempirical version of the scaled-down SIC-LSDA, which is significantly better than uncorrected LSDA and even better than the uncorrected generalized gradient approximation. We present an analysis of the formal properties of the scaled-down SIC and define possible directions for further improvements. In particular, we find that exactness for all one-electron densities does not guarantee correct asymptotics for the exchange-correlation potential of a many-electron system.

Binding or potential energy curves have been calculated for the ground-state diatomics H(2)(+), He(2)(+), LiH(+), H(2), N(2), and C(2), for the transition state H(3), and for the triplet first excited state of H(2) using the nonempirical density functionals from the first three rungs of a ladder of approximations: the local spin density (LSD) approximation, the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and the Tao-Perdew-Staroverov-Scuseria (TPSS) meta GGA. Good binding energy curves in agreement with coupled cluster or configuration interaction calculations are found from the PBE GGA and especially from the TPSS meta GGA. Expected exceptions are the symmetric radicals H(2)(+) and He(2)(+), where the functionals suffer from self-interaction error, and the exotically bonded C(2). Although the energy barrier for the reaction H(2) + H --> H + H(2) is better in PBE than in TPSS, the transition state H(3) is a more properly positioned and curved saddle point of the energy surface in TPSS. The triplet first excited state of H(2) obeys the Aufbau principle and thus is one of the exceptional excited states that are computable in principle from the ground-state functional. The PBE GGA and TPSS meta GGA are useful not only for chemical applications but also for the construction of higher-rung nonempirical functionals that can further improve the binding energy curves.

Binding energy curves have been calculated for the ground-state rare-gas diatomics Ne(2) and Ar(2) and for the alkaline-earth diatomic Be(2) using the nonempirical density functionals from the first three rungs of a ladder of approximations: the local spin density (LSD) approximation, the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and the Tao-Perdew-Staroverov-Scuseria (TPSS) meta-GGA. Binding energy curves in reasonable agreement with those constructed from experiment are found from PBE and TPSS, which incorporate inhomogeneity corrections that satisfy the Lieb-Oxford bound and so describe the short-range part of the van der Waals interaction. At large internuclear separation, these functionals produce an exponentially decaying attraction in place of the correct long-range -C(6)/R(6). Basis-set and exchange-only effects are also discussed.

The accuracy of the Perdew-Burke-Ernzerhof and Tao-Perdew-Staroverov-Scuseria density functionals for describing noncovalent interaction energies in small water clusters is studied by testing 11 basis sets on a reduced test set proposed by Dahlke and Truhlar (J. Phys. Chem. B 2005, 109, 15677). We have also tested variants of the Perdew-Burke-Ernzerhof functional and the Becke98 hybrid functional. While moderate basis sets give converged density functional theory results for covalent dissociation energies, this is not true for noncovalent interaction energies. Our results show that density functionals give converged interaction energies with aug-cc-pVTZ and aug-cc-pVQZ basis sets. Gradual simplification of the basis set introduces an increasing overbinding effect. The best agreement with the high-level result was obtained by the Perdew-Burke-Ernzerhof functional at the basis set limit. The converged Tao-Perdew-Staroverov-Scuseria interaction energies show a systematic underbinding effect that can be compensated by a somewhat systematic overbinding basis set effect of smaller basis sets such as 6-31+G(d,2p). The inclusion of the diffuse functions in the oxygen basis set is very important, while the inclusion of the f functions practically does not influence the results. Improvement can be obtained by adding more hydrogen p functions to the 6-31+G basis set.

We present the case for the nonempirical construction of density functional approximations for the exchange-correlation energy by the traditional method of "constraint satisfaction" without fitting to data sets, and present evidence that this approach has been successful on the first three rungs of "Jacob's ladder" of density functional approximations [local spin-density approximation (LSD), generalized gradient approximation (GGA), and meta-GGA]. We expect that this approach will also prove successful on the fourth and fifth rungs (hyper-GGA or hybrid and generalized random-phase approximation). In particular, we argue for the theoretical and practical importance of recovering the correct uniform density limit, which many semiempirical functionals fail to do. Among the beyond-LSD functionals now available to users, we recommend the nonempirical Perdew-Burke-Ernzerhof (PBE) GGA and the nonempirical Tao-Perdew-Staroverov-Scuseria (TPSS) meta-GGA, and their one-parameter hybrids with exact exchange. TPSS improvement over PBE is dramatic for atomization energies of molecules and surface energies of solids, and small or moderate for other properties. TPSS is now or soon will be available in standard codes such as GAUSSIAN, TURBOMOLE, NWCHEM, ADF, WIEN, VASP, etc. We also discuss old and new ideas to eliminate the self-interaction error that plagues the functionals on the first three rungs of the ladder, bring up other related issues, and close with a list of "do's and don't's" for software developers and users.

For accurate thermochemical tests of electronic structure theory, accurate true anharmonic zero-point vibrational energies ZPVE(true) are needed. We discuss several possibilities to extract this information for molecules from density functional or wave function calculations and/or available experimental data: (1) Empirical universal scaling of density-functional-calculated harmonic ZPVE(harm)s, where we find that polyatomics require smaller scaling factors than diatomics. (2) Direct density-functional calculation by anharmonic second-order perturbation theory PT2. (3) Weighted averages of harmonic ZPVE(harm) and fundamental ZPVE(fund) (from fundamental vibrational transition frequencies), with weights (3/4, 1/4) for diatomics and (5/8,3/8) for polyatomics. (4) Experimental correction of the PT2 harmonic contribution, i.e., the estimate ZPVE(true)PT2 + (ZPVE(fund)expt - ZPVE(fund)PT2) for ZPVE(true). The (5/8,3/8) average of method 3 and the additive correction of method 4 have been proposed here. For our database of experimental ZPVE(true), consisting of 27 diatomics and 8 polyatomics, we find that methods 1 and 2, applied to the popular B3LYP and the nonempirical PBE and TPSS functionals and their one-parameter hybrids, yield polyatomic errors on the order of 0.1 kcal/mol. Larger errors are expected for molecules larger than those in our database. Method 3 yields errors on the order of 0.02 kcal/mol, but requires very accurate (e.g., experimental, coupled cluster, or best-performing density functional) input harmonic ZPVE(harm). Method 4 is the best-founded one that meets the requirements of high accuracy and practicality, requiring as experimental input only the highly accurate and widely available ZPVE(fund)expt and producing errors on the order of 0.05 kcal/mol that are relatively independent of functional and basis set. As a part of our study, we also test the ability of the density functionals to predict accurate equilibrium bond lengths and angles for a data set of 21 mostly polyatomic molecules (since all calculated ZPVEs are evaluated at the correspondingly calculated molecular geometries).

Relative stabilities of rotamers of the N-acetyl-O-(2-acetamido-2-deoxy-alpha-D-galactopyranosyl)-L-seryl-N'-methyl amide (1) and eleven analogous molecules containing beta-galactose, alpha- and beta-mannose, alpha- and beta-glucose, and L-threonine were calculated to learn whether they could explain the natural preference for 1 in linkages between the carbohydrate and protein in glycoproteins. The lowest energy rotamers of four O-glycoside models of serine diamide were identified with a Monte Carlo search coupled with molecular mechanics (MM2*). These rotamers were further optimized with an ab initio level of theory (HF/6-31G(d)). Subsequently, B3LYP/6-31 + G(d) single point energies were calculated for the most stable HF structures. The most favorable interactions are present in 1 and its glucose analogue. The monosaccharide for the carbohydrate antenna is anchored to the serine residue with an AcNH...O=C-NHMe hydrogen bond in the most stable rotamers. The mannose analogue and the beta-anomers are considerably less stable according to the MM2* and especially to the ab inito energy values. The three analogues have HF/6-31 G(d) energies which are 4-6 kcal mol-1 higher; the single point B3LYP/6-31 + G(d)//HF/6-31 G(d) calculations yield preferences of 3-5 kcal mol-1 for 1. The most stable L-threonine analogues show a behaviour very similarly to the corresponding serine analogues. The ZPE and thermal correction components of the calculated delta H298 and delta G298 values are relatively small (< 0.4 kcal mol-1). However, the T delta S298 term can be as large as 2.6 kcal mol-1. The entropy terms stabilize the alpha-anomers relative to beta-anomers, and ManNAc relative to GalNAc. The largest stabilization effect is observed for one of the rotamers of the alpha-anomer of ManNAc.

A brief survey is given on the synthesis and structure elucidation for dihydro-1,2,4-triazines. The relative stability of nine possible dihydro-1,2,4-triazines and three dihydrotriazinium cations is studied at HF, MP2, generalized gradient approximation DFT, and CBS-4 levels of theory. The structural consequences of the inclusion of the electron correlation are also given. We attempt to rationalize the experimental findings using high-quality theoretical results. The quantum chemical calculations support that the most stable isomer is the 2,5-dihydro-1,2,4-triazine and all the other relatively stable isomers have been experimentally identified correctly. Several experimental papers report structures that have been proved to be nonexistent. These structures have energy that is too high according to the best-quality calculations.