Search

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. [ABSTRACT FROM AUTHOR]

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-1≤Na=M, and so can avoid these errors. This desideratum is a natural generalization to all M of the more familiar one of one-electron self-interaction freedom. The intent of this paper is not to advocate for any functional, but to understand what is required for a functional to be M-electron self-interaction-free and thus correct even for highly stretched bonds. The original Perdew-Zunger self-interaction correction (SIC) and our scaled-down variant of it are exactly one- and nearly two-electron self-interaction-free, but only the former is nearly so for atoms with M>2. 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. [ABSTRACT FROM AUTHOR]

Semilocal density functional approximations for the exchange-correlation energy can improperly dissociate a neutral molecule XY (Y≠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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

While most molecules and solids are spin-unpolarized, most chemically-active atoms are partly spin-polarized. As a result, the errors of the spin-dependence of a density functional are much more troublesome for atomization energies than they are for typical reaction or formation energies. This observation explains why the atomization energy errors of approximate functionals do not correlate with their other errors, and why the errors of atomization energies for a given functional can be radically reduced by fitting the energies of the atoms. We present an illustrative example from the recent nonempirical construction of the SCAN meta-generalized gradient approximation. [ABSTRACT FROM AUTHOR]

In the assessment of density functional approximations for the exchange‐correlation energy, great weight is usually given to the accuracy of molecular atomization energies, or the enthalpies of formation constructed from calculated atomization energies. Here we show that a recent nonempirical functional, the meta‐generalized gradient approximation of Tao, Perdew, Staroverov, and Scuseria, achieves remarkably accurate atomization energies even for the larger organic molecules of the G3‐3 test set. But we also present strong evidence that most of the error of previous nonempirical functionals resides in the energy of the free atom and so cancels out of typical reaction energies. Finally, we suggest that enthalpies of formation calculated without any reference to the free atoms would provide a fairer assessment of the performance of approximate density functionals. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2005 [ABSTRACT FROM AUTHOR]

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-1≤Na=M, and so can avoid these errors. This desideratum is a natural generalization to all M of the more familiar one of one-electron self-interaction freedom. The intent of this paper is not to advocate for any functional, but to understand what is required for a functional to be M-electron self-interaction-free and thus correct even for highly stretched bonds. The original Perdew-Zunger self-interaction correction (SIC) and our scaled-down variant of it are exactly one- and nearly two-electron self-interaction-free, but only the former is nearly so for atoms with M>2. 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Abstract: 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. [Copyright &y& Elsevier]

Abstract: Eighteen empirical force fields and the semi-empirical quantum method PM3CARB-1 were compared for studying β-cellobiose, α-maltose, and α-galabiose [α-d-Galp-(1→4)-α-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. [Copyright &y& Elsevier]