Gábor I. Csonka

Restoring the density-gradient expansion for exchange in solids and surfaces.

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.Successful modern generalized gradient approximations (GGAs) are biased toward atomic energies. Restoration of the first-principles gradient expansion for the exchange energy over a wide range of density gradients eliminates this bias. With many collaborators, I introduce PBEsol, a revised Perdew-Burke-Ernzerhof GGA that improves equilibrium properties of densely-packed solids and their surfaces.

Assessing the performance of recent density functionals for bulk solids

We assess the performance of recent density functionals for the exchange-correlation energy of a nonmolecular solid, by applying accurate calculations with the GAUSSIAN, BAND, and VASP codes to a test set of 24 solid metals and nonmetals. The functionals tested are the modified Perdew-Burke-Ernzerhof generalized gradient approximation PBEsol GGA, the second-order GGA SOGGA, and the Armiento-Mattsson 2005 AM05 GGA. For completeness, we also test more standard functionals: the local density approximation, the original PBE GGA, and the Tao-Perdew-Staroverov-Scuseria meta-GGA. We find that the recent density functionals for solids reach a high accuracy for bulk properties lattice constant and bulk modulus. For the cohesive energy, PBE is better than PBEsol overall, as expected, but PBEsol is actually better for the alkali metals and alkali halides. For fair comparison of calculated and experimental results, we consider the zeropoint phonon and finite-temperature effects ignored by many workers. We show how GAUSSIAN basis sets and inaccurate experimental reference data may affect the rating of the quality of the functionals. The results show that PBEsol and AM05 perform somewhat differently from each other for alkali metal, alkaline-earth metal, and alkali halide crystals where the maximum value of the reduced density gradient is about 2, but perform very similarly for most of the other solids where it is often about 1. Our explanation for this is consistent with the importance of exchange-correlation nonlocality in regions of core-valence overlap.

Evaluation of Density Functionals and Basis Sets for Carbohydrates

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.

Exchange and correlation in open systems of fluctuating electron number

While the exact total energy of a separated open system varies linearly as a function of average electron number between adjacent integers, the energy predicted by semilocal density-functional approximations is concave up and the exact-exchange-only or Hartree-Fock energy is concave down. As a result, semilocal density functionals fail for separated open systems of fluctuating electron number, as in stretched molecular ions A{sub 2}{sup +} and in solid transition-metal oxides. We develop an exact-exchange theory and an exchange-hole sum rule that explain these failures and we propose a way to correct them via a local hybrid functional.

Proper basis set for quantum mechanical studies of potential energy surfaces of carbohydrates

Abstract In the current paper we address the optimal selection of basis set for carbohydrates. The following basis sets are used with B3LYP method: 6-31G(d), 6-31G(d,p) 5d, 6-31+G(d,p), 6-311+G(d,p), 6-311++G(d,p), 6-311+G(2d,2p), and 6-311++G(2d,2p). The following molecules were used for the comparison: water, conformers of α- l -fucose, and β- d -glucose. The 6-31G(d) or 6-31G(d,p) basis set in combination with B3LYP functional provide unacceptably poor results for carbohydrates. (Although the HF/6-31G(d) results are quite good.) The introduction of the diffuse functions on heavy atoms is necessary for good results if B3LYP functional is used (suggested methods are B3LYP/6-31+G(d,p) or 6-311+G(d,p), the latter is closer to basis set limit). The introduction of the diffuse functions (++) on hydrogen atoms is not necessary if B3LYP or other density functional method is used for carbohydrates

THE ORIGIN OF THE PROBLEMS WITH THE PM3 CORE REPULSION FUNCTION

Abstract We investigate the cause of failures of PM3 for H⋯H, O⋯H and N⋯H interactions. We show that the actual choice of parameters for the Gaussian correction functions leads to spurious oscillations in the interatomic core repulsion functions and in the potential energy surfaces. The positions of the minima of the Gaussian correction functions considerably influence the positions of the minima on the potential energy hypersurface of weak interactions. The PM3 potential energy curve of water dimer as a function of the O⋯H distance clearly illustrates the problem. The minimum at 1.8 A and the wide shoulder at 2.1 A correspond to a minimum and a maximum of the O⋯H Gaussian correction function, respectively.

Comparison of different force fields for the study of disaccharides

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.

Scaling down the Perdew-Zunger self-interaction correction in many-electron regions

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.

Density functionals that recognize covalent, metallic, and weak bonds.

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.

Relative stability of 1C4 and 4C1 chair forms of β-d-glucose: a density functional study

Abstract The method and basis set dependence of the relative energies of the 1 C 4 and 4 C 1 chair forms of β- d -glucose were calculated for two selected, low-energy hydroxyl rotamers at various levels of generalized gradient approximation density functional theory (GGA-DFT). The GGA-DFT and MP2 methods provide similar energetic differences for β- d -glucose conformers. Addition of the diffuse functions to a double-zeta quality basis set and inclusion of the HF exchange into the DFT functionals improve the agreement between the DFT and the best composite estimates of the energetic differences. The GGA- or hybrid-DFT methods reproduce the geometrical consequences of correlation effects correctly for glucose.