Harmonic Vibrational Frequencies: An Evaluation of Hartree−Fock, Møller−Plesset, Quadratic Configuration Interaction, Density Functional Theory, and Semiempirical Scale Factors
10 Oct 1996-The Journal of Physical Chemistry (American Chemical Society)-Vol. 100, Iss: 41, pp 16502-16513
TL;DR: In this paper, scaling factors for fundamental vibrational frequencies, low-frequency vibrations, zero-point vibrational energies (ZPVE), and thermal contributions to enthalpy and entropy from harmonic frequencies determined at 19 levels of theory have been derived through a least-squares approach.
Abstract: Scaling factors for obtaining fundamental vibrational frequencies, low-frequency vibrations, zero-point vibrational energies (ZPVE), and thermal contributions to enthalpy and entropy from harmonic frequencies determined at 19 levels of theory have been derived through a least-squares approach. Semiempirical methods (AM1 and PM3), conventional uncorrelated and correlated ab initio molecular orbital procedures [Hartree−Fock (HF), Moller−Plesset (MP2), and quadratic configuration interaction including single and double substitutions (QCISD)], and several variants of density functional theory (DFT: B-LYP, B-P86, B3-LYP, B3-P86, and B3-PW91) have been examined in conjunction with the 3-21G, 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-311G(d,p), and 6-311G(df,p) basis sets. The scaling factors for the theoretical harmonic vibrational frequencies were determined by a comparison with the corresponding experimental fundamentals utilizing a total of 1066 individual vibrations. Scaling factors suitable for low-frequency vib...
Citations
More filters
TL;DR: The M06-2X meta-exchange correlation function is proposed in this paper, which is parametrized including both transition metals and nonmetals, and is a high-non-locality functional with double the amount of nonlocal exchange.
Abstract: We present two new hybrid meta exchange- correlation functionals, called M06 and M06-2X. The M06 functional is parametrized including both transition metals and nonmetals, whereas the M06-2X functional is a high-nonlocality functional with double the amount of nonlocal exchange (2X), and it is parametrized only for nonmetals.The functionals, along with the previously published M06-L local functional and the M06-HF full-Hartree–Fock functionals, constitute the M06 suite of complementary functionals. We assess these four functionals by comparing their performance to that of 12 other functionals and Hartree–Fock theory for 403 energetic data in 29 diverse databases, including ten databases for thermochemistry, four databases for kinetics, eight databases for noncovalent interactions, three databases for transition metal bonding, one database for metal atom excitation energies, and three databases for molecular excitation energies. We also illustrate the performance of these 17 methods for three databases containing 40 bond lengths and for databases containing 38 vibrational frequencies and 15 vibrational zero point energies. We recommend the M06-2X functional for applications involving main-group thermochemistry, kinetics, noncovalent interactions, and electronic excitation energies to valence and Rydberg states. We recommend the M06 functional for application in organometallic and inorganometallic chemistry and for noncovalent interactions.
22,326 citations
TL;DR: An extension of the CHARMM force field to drug‐like molecules is presented, making it possible to perform “all‐CHARMM” simulations on drug‐target interactions thereby extending the utility ofCHARMM force fields to medicinally relevant systems.
Abstract: The widely used CHARMM additive all-atom force field includes parameters for proteins, nucleic acids, lipids, and carbohydrates. In the present article, an extension of the CHARMM force field to drug-like molecules is presented. The resulting CHARMM General Force Field (CGenFF) covers a wide range of chemical groups present in biomolecules and drug-like molecules, including a large number of heterocyclic scaffolds. The parametrization philosophy behind the force field focuses on quality at the expense of transferability, with the implementation concentrating on an extensible force field. Statistics related to the quality of the parametrization with a focus on experimental validation are presented. Additionally, the parametrization procedure, described fully in the present article in the context of the model systems, pyrrolidine, and 3-phenoxymethylpyrrolidine will allow users to readily extend the force field to chemical groups that are not explicitly covered in the force field as well as add functional groups to and link together molecules already available in the force field. CGenFF thus makes it possible to perform "all-CHARMM" simulations on drug-target interactions thereby extending the utility of CHARMM force fields to medicinally relevant systems.
4,553 citations
TL;DR: The new local density functional, called M06-L, is designed to capture the main dependence of the exchange-correlation energy on local spin density, spin density gradient, and spin kinetic energy density, and it is parametrized to satisfy the uniform-electron-gas limit.
Abstract: We present a new local density functional, called M06-L, for main-group and transition element thermochemistry, thermochemical kinetics, and noncovalent interactions. The functional is designed to capture the main dependence of the exchange-correlation energy on local spin density, spin density gradient, and spin kinetic energy density, and it is parametrized to satisfy the uniform-electron-gas limit and to have good performance for both main-group chemistry and transition metal chemistry. The M06-L functional and 14 other functionals have been comparatively assessed against 22 energetic databases. Among the tested functionals, which include the popular B3LYP, BLYP, and BP86 functionals as well as our previous M05 functional, the M06-L functional gives the best overall performance for a combination of main-group thermochemistry, thermochemical kinetics, and organometallic, inorganometallic, biological, and noncovalent interactions. It also does very well for predicting geometries and vibrational frequencies. Because of the computational advantages of local functionals, the present functional should be very useful for many applications in chemistry, especially for simulations on moderate-sized and large systems and when long time scales must be addressed. © 2006 American Institute of Physics. DOI: 10.1063/1.2370993
4,154 citations
TL;DR: This chapter discusses the development of DFT as a tool for Calculating Atomic andMolecular Properties and its applications, as well as some of the fundamental and Computational aspects.
Abstract: I. Introduction: Conceptual vs Fundamental andComputational Aspects of DFT1793II. Fundamental and Computational Aspects of DFT 1795A. The Basics of DFT: The Hohenberg−KohnTheorems1795B. DFT as a Tool for Calculating Atomic andMolecular Properties: The Kohn−ShamEquations1796C. Electronic Chemical Potential andElectronegativity: Bridging Computational andConceptual DFT1797III. DFT-Based Concepts and Principles 1798A. General Scheme: Nalewajski’s ChargeSensitivity Analysis1798B. Concepts and Their Calculation 18001. Electronegativity and the ElectronicChemical Potential18002. Global Hardness and Softness 18023. The Electronic Fukui Function, LocalSoftness, and Softness Kernel18074. Local Hardness and Hardness Kernel 18135. The Molecular Shape FunctionsSimilarity 18146. The Nuclear Fukui Function and ItsDerivatives18167. Spin-Polarized Generalizations 18198. Solvent Effects 18209. Time Evolution of Reactivity Indices 1821C. Principles 18221. Sanderson’s Electronegativity EqualizationPrinciple18222. Pearson’s Hard and Soft Acids andBases Principle18253. The Maximum Hardness Principle 1829IV. Applications 1833A. Atoms and Functional Groups 1833B. Molecular Properties 18381. Dipole Moment, Hardness, Softness, andRelated Properties18382. Conformation 18403. Aromaticity 1840C. Reactivity 18421. Introduction 18422. Comparison of Intramolecular ReactivitySequences18443. Comparison of Intermolecular ReactivitySequences18494. Excited States 1857D. Clusters and Catalysis 1858V. Conclusions 1860VI. Glossary of Most Important Symbols andAcronyms1860VII. Acknowledgments 1861VIII. Note Added in Proof 1862IX. References 1865
3,890 citations
Book•
01 Sep 2001
TL;DR: A Chemist's Guide to Density Functional Theory should be an invaluable source of insight and knowledge for many chemists using DFT approaches to solve chemical problems.
Abstract: "Chemists familiar with conventional quantum mechanics will applaud and benefit greatly from this particularly instructive, thorough and clearly written exposition of density functional theory: its basis, concepts, terms, implementation, and performance in diverse applications. Users of DFT for structure, energy, and molecular property computations, as well as reaction mechanism studies, are guided to the optimum choices of the most effective methods. Well done!" Paul von RaguE Schleyer "A conspicuous hole in the computational chemist's library is nicely filled by this book, which provides a wide-ranging and pragmatic view of the subject.[...It] should justifiably become the favorite text on the subject for practioneers who aim to use DFT to solve chemical problems." J. F. Stanton, J. Am. Chem. Soc. "The authors' aim is to guide the chemist through basic theoretical and related technical aspects of DFT at an easy-to-understand theoretical level. They succeed admirably." P. C. H. Mitchell, Appl. Organomet. Chem. "The authors have done an excellent service to the chemical community. [...] A Chemist's Guide to Density Functional Theory is exactly what the title suggests. It should be an invaluable source of insight and knowledge for many chemists using DFT approaches to solve chemical problems." M. Kaupp, Angew. Chem.
3,550 citations
References
More filters
TL;DR: Numerical calculations on a number of atoms, positive ions, and molecules, of both open- and closed-shell type, show that density-functional formulas for the correlation energy and correlation potential give correlation energies within a few percent.
Abstract: A correlation-energy formula due to Colle and Salvetti [Theor. Chim. Acta 37, 329 (1975)], in which the correlation energy density is expressed in terms of the electron density and a Laplacian of the second-order Hartree-Fock density matrix, is restated as a formula involving the density and local kinetic-energy density. On insertion of gradient expansions for the local kinetic-energy density, density-functional formulas for the correlation energy and correlation potential are then obtained. Through numerical calculations on a number of atoms, positive ions, and molecules, of both open- and closed-shell type, it is demonstrated that these formulas, like the original Colle-Salvetti formulas, give correlation energies within a few percent.
84,646 citations
TL;DR: This work reports a gradient-corrected exchange-energy functional, containing only one parameter, that fits the exact Hartree-Fock exchange energies of a wide variety of atomic systems with remarkable accuracy, surpassing the performance of previous functionals containing two parameters or more.
Abstract: Current gradient-corrected density-functional approximations for the exchange energies of atomic and molecular systems fail to reproduce the correct 1/r asymptotic behavior of the exchange-energy density. Here we report a gradient-corrected exchange-energy functional with the proper asymptotic limit. Our functional, containing only one parameter, fits the exact Hartree-Fock exchange energies of a wide variety of atomic systems with remarkable accuracy, surpassing the performance of previous functionals containing two parameters or more.
45,683 citations
TL;DR: A simple analytic representation of the correlation energy for a uniform electron gas, as a function of density parameter and relative spin polarization \ensuremath{\zeta}, which confirms the practical accuracy of the VWN and PZ representations and eliminates some minor problems.
Abstract: We propose a simple analytic representation of the correlation energy ${\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}$ for a uniform electron gas, as a function of density parameter ${\mathit{r}}_{\mathit{s}}$ and relative spin polarization \ensuremath{\zeta}. Within the random-phase approximation (RPA), this representation allows for the ${\mathit{r}}_{\mathit{s}}^{\mathrm{\ensuremath{-}}3/4}$ behavior as ${\mathit{r}}_{\mathit{s}}$\ensuremath{\rightarrow}\ensuremath{\infty}. Close agreement with numerical RPA values for ${\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}$(${\mathit{r}}_{\mathit{s}}$,0), ${\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}$(${\mathit{r}}_{\mathit{s}}$,1), and the spin stiffness ${\mathrm{\ensuremath{\alpha}}}_{\mathit{c}}$(${\mathit{r}}_{\mathit{s}}$)=${\mathrm{\ensuremath{\partial}}}^{2}$${\mathrm{\ensuremath{\varepsilon}}}_{\mathit{c}}$(${\mathit{r}}_{\mathit{s}}$, \ensuremath{\zeta}=0)/\ensuremath{\delta}${\mathrm{\ensuremath{\zeta}}}^{2}$, and recovery of the correct ${\mathit{r}}_{\mathit{s}}$ln${\mathit{r}}_{\mathit{s}}$ term for ${\mathit{r}}_{\mathit{s}}$\ensuremath{\rightarrow}0, indicate the appropriateness of the chosen analytic form. Beyond RPA, different parameters for the same analytic form are found by fitting to the Green's-function Monte Carlo data of Ceperley and Alder [Phys. Rev. Lett. 45, 566 (1980)], taking into account data uncertainties that have been ignored in earlier fits by Vosko, Wilk, and Nusair (VWN) [Can. J. Phys. 58, 1200 (1980)] or by Perdew and Zunger (PZ) [Phys. Rev. B 23, 5048 (1981)]. While we confirm the practical accuracy of the VWN and PZ representations, we eliminate some minor problems with these forms. We study the \ensuremath{\zeta}-dependent coefficients in the high- and low-density expansions, and the ${\mathit{r}}_{\mathit{s}}$-dependent spin susceptibility. We also present a conjecture for the exact low-density limit. The correlation potential ${\mathrm{\ensuremath{\mu}}}_{\mathit{c}}^{\mathrm{\ensuremath{\sigma}}}$(${\mathit{r}}_{\mathit{s}}$,\ensuremath{\zeta}) is evaluated for use in self-consistent density-functional calculations.
21,353 citations
TL;DR: In this article, the unpolarized absorption and circular dichroism spectra of the fundamental vibrational transitions of the chiral molecule, 4-methyl-2-oxetanone, are calculated ab initio using DFT, MP2, and SCF methodologies and a 5S4P2D/3S2P (TZ2P) basis set.
Abstract: : The unpolarized absorption and circular dichroism spectra of the fundamental vibrational transitions of the chiral molecule, 4-methyl-2-oxetanone, are calculated ab initio. Harmonic force fields are obtained using Density Functional Theory (DFT), MP2, and SCF methodologies and a 5S4P2D/3S2P (TZ2P) basis set. DFT calculations use the Local Spin Density Approximation (LSDA), BLYP, and Becke3LYP (B3LYP) density functionals. Mid-IR spectra predicted using LSDA, BLYP, and B3LYP force fields are of significantly different quality, the B3LYP force field yielding spectra in clearly superior, and overall excellent, agreement with experiment. The MP2 force field yields spectra in slightly worse agreement with experiment than the B3LYP force field. The SCF force field yields spectra in poor agreement with experiment.The basis set dependence of B3LYP force fields is also explored: the 6-31G* and TZ2P basis sets give very similar results while the 3-21G basis set yields spectra in substantially worse agreements with experiment. jg
17,871 citations
TL;DR: A way is found to visualize and understand the nonlocality of exchange and correlation, its origins, and its physical effects as well as significant interconfigurational and interterm errors remain.
Abstract: Generalized gradient approximations (GGA's) seek to improve upon the accuracy of the local-spin-density (LSD) approximation in electronic-structure calculations. Perdew and Wang have developed a GGA based on real-space cutoff of the spurious long-range components of the second-order gradient expansion for the exchange-correlation hole. We have found that this density functional performs well in numerical tests for a variety of systems: (1) Total energies of 30 atoms are highly accurate. (2) Ionization energies and electron affinities are improved in a statistical sense, although significant interconfigurational and interterm errors remain. (3) Accurate atomization energies are found for seven hydrocarbon molecules, with a rms error per bond of 0.1 eV, compared with 0.7 eV for the LSD approximation and 2.4 eV for the Hartree-Fock approximation. (4) For atoms and molecules, there is a cancellation of error between density functionals for exchange and correlation, which is most striking whenever the Hartree-Fock result is furthest from experiment. (5) The surprising LSD underestimation of the lattice constants of Li and Na by 3--4 % is corrected, and the magnetic ground state of solid Fe is restored. (6) The work function, surface energy (neglecting the long-range contribution), and curvature energy of a metallic surface are all slightly reduced in comparison with LSD. Taking account of the positive long-range contribution, we find surface and curvature energies in good agreement with experimental or exact values. Finally, a way is found to visualize and understand the nonlocality of exchange and correlation, its origins, and its physical effects.
17,848 citations