Extension of the B3LYP–dispersion-correcting potential approach to the accurate treatment of both inter- and intra-molecular interactions
Tóm tắt
We recently demonstrated that dispersion-correcting potentials (DCPs), which are atom-centered Gaussian-type functions that were developed for use with B3LYP (Torres and DiLabio in J Phys Chem Lett 3:1738–1744, 2012), greatly improved the ability of the underlying functional to predict non-covalent interactions. However, the recent application of the B3LYP–DCP approach to study the β-scission of the cumyloxyl radical led to a calculated barrier height that was over-estimated by ca. 8 kcal/mol. We demonstrate in the present work that the source of this error arises from the previously developed carbon atom DCPs, which erroneously alters the electron density in the C–C covalent-bonding region. In this work, we developed a new C-DCP with a form that was expected to less strongly influence the electron density in the covalent bonding region. Tests of the new C-DCP, in conjunction with previously published H-, N-, and O-DCPs, with B3LYP–DCP/6-31+G(2d,2p) on the S66, S22B, HSG-A, and HC12 databases of non-covalently interacting dimers showed that it is one of the most accurate methods available for treating intermolecular interactions, giving mean absolute errors (MAEs) of 0.19, 0.27, 0.16, and 0.18 kcal/mol, respectively. Additional testing on the S12L database of very large complexation systems gave an MAE of 2.6 kcal/mol, demonstrating that the B3LYP–DCP/6-31+G(2d,2p) approach to be one of the best-performing and most feasible methods for treating large systems containing significant non-covalent interactions. Finally, we showed that the modeling of C–C-making/C–C-breaking chemistry is well predicted using the newly developed DCPs. In addition to predicting a barrier height for the β-scission of the cumyloxyl radical, that is, within 1.7 kcal/mol of the high-level value, application of B3LYP–DCP/6-31+G(2d,2p) to 10 databases that include reaction barrier heights and energies, isomerization energies, and relative conformation energies gives performance that is among the best of all available dispersion-corrected density-functional theory approaches.
Tài liệu tham khảo
Pauling L, Corey RB, Branson HR (1951) The structure of proteins: two hydrogen-bonded helical configurations of the polypeptide chain. Proc Natl Acad Sci USA 37:205–211
Astbury WT (1933) Some problems in the X-ray analysis of the structure of animal hairs and other protein fibres. Trans Faraday Soc 29:193–205
Wu X, Vargas MC, Nayak S, Lotrich V, Scoles G (2001) Towards extending the applicability of density functional theory to weakly bound systems. J Chem Phys 115:8748–8757
Wu Q, Yang W (2002) Empirical correction to density functional theory for van der Waals interactions. J Chem Phys 116:515–524
Grimme S (2004) Accurate description of van der Waals complexes by density functional theory including empirical corrections. J Comput Chem 25:1463–1473
Grimme S (2006) Semiempirical GGA-type density functional constructed with a long-range dispersion correction. J Comput Chem 27:1787–1799
Grimme S, Antony J, Ehrlich S, Krieg H (2010) A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H–Pu. J Chem Phys 132:154104
Johnson ER, Becke AD (2005) A post-Hartree–Fock model of intermolecular interactions. J Chem Phys 123:024101
Becke AD, Arabi AA, Kannemann FO (2010) Nonempirical density-functional theory for van der Waals interactions. Can J Chem 88:1057–1062
Otero-de-la-Roza A, Johnson ER (2012) A benchmark for non-covalent interactions in solids. J Chem Phys 137:054103
Tkatchenko A, Scheffler M (2009) Accurate molecular van der Waals interactions from ground-state electron density and free-atom reference data. Phys Rev Lett 102:073005
Marom N, Tkatchenko A, Rossi M, Gobre VV, Hod O, Scheffler M, Kronik L (2011) Dispersion interactions with density-functional theory: benchmarking semiempirical and interatomic pairwise corrected density functionals. J Chem Theory Comput 7:3944–3951
Zhao Y, Truhlar DG (2008) The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor Chem Acc 120:215–241
Zhao Y, Truhlar DG (2008) Exploring the Limit of accuracy of the global hybrid meta density functional for main-group thermochemistry, kinetics, and noncovalent interactions. J Chem Theory Comput 4:1849–1868
DiLabio GA (2008) Accurate treatment of van der Waals interactions using standard density functional theory with effective core-type potentials: application to carbon-containing dimers. Chem Phys Lett 455:348–353
Mackie ID, DiLabio GA (2008) Interactions in large, polyaromatic hydrocarbons dimers: application of density functional theory with dispersion corrections. J Phys Chem A 112:10968–10976
Johnson ER, DiLabio GA (2009) Theoretical study of dispersion binding of hydrocarbon molecules to hydrogen-terminated Silicon(100)-2 × 1. J Phys Chem C 113:5681–5689
Mackie ID, DiLabio GA (2010) Accurate dispersion interactions from standard density-functional theory methods with small basis sets. Phys Chem Chem Phys 12:6092–6098
Torres E, DiLabio GA (2012) A (nearly) universally applicable method for modeling noncovalent interactions using B3LYP. J Phys Chem Lett 3:1738–1744
von Lilienfeld OA, Tavernelli I, Rothlisberger U (2004) Optimization of effective atom centered potentials for London dispersion forces in density functional theory. Phys Rev Lett 93:153004
von Lilienfeld OA, Tavernelli I, Rothlisberger U (2005) Performance of optimized atom-centered potentials for weakly bonded systems using density functional theory. Phys Rev B 71:195119
Lin I-C, Coutinho-Neto MD, Felsenheimer C, von Lilienfeld OA, Tavernelli I, Rothlisberger U (2007) Library of dispersion-corrected atom-centered potentials for generalized gradient approximation functionals: elements H, C, N, O, He, Ne, Ar, and Kr. Phys Rev B 75:205131
Christiansen PA, Lee YS, Pitzer KS (1979) Improved ab initio effective core potentials for molecular calculations. J Chem Phys 71:4445–4450
Becke AD (1988) Density-functional exchange-energy approximation with correct asymptotic behavior. Phys Rev A 38:3098–3100
Lee C, Yang W, Parr RG (1988) Development of the Colle–Salvetti correlation-energy formula into a functional of the electron density. Phys Rev B 37:785–789
D’Alfonso C, Bietti M, DiLabio GA, Lanzalunga O, Salamone M (2013) Reactions of the phthalimide N-oxyl radical (PINO) with activated phenols: the contribution of pi-stacking interactions to hydrogen atom transfer rates. J Org Chem 78:1026–1037
Salamone M, DiLabio GA, Bietti M (2012) Reactions of the cumyloxyl and benzyloxyl radicals with strong hydrogen bond acceptors. Large enhancements in hydrogen abstraction reactivity determined by substrate/radical hydrogen bonding. J Org Chem 77:10479–10487
Salamone M, DiLabio GA, Bietti M (2011) Hydrogen atom abstraction selectivity in the reactions of alkylamines with the benzyloxyl and cumyloxyl radicals. The importance of structure and of substrate radical hydrogen bonding. J Am Chem Soc 133:16625–16634
Torres E, DiLabio GA (2013) Density-functional theory with dispersion-correcting potentials for methane: bridging the efficiency and accuracy gap between high-level wave function and classical molecular mechanics methods. J Chem Theory Comput. doi:10.1021/ct4003114
Boys SF, Bernardi F (1970) The calculation of small molecular interactions by the differences of separate total energies. Some procedures with reduced errors. Mol Phys 19:553–566
Sherrill CD, Takatani T, Hohenstein EG (2009) An assessment of theoretical methods for nonbonded interactions: comparison to complete basis set limit coupled-cluster potential energy curves for the benzene dimer, the methane dimer, benzene–methane, and benzene–H2S. J Phys Chem A 113:10146–10159
Mackie ID, DiLabio GA (2011) Approximations to complete basis set-extrapolated, highly correlated non-covalent interaction energies. J Chem Phys 135:134318
Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Mennucci B, Petersson GA, Nakatsuji H, Caricato M, Li X, Hratchian HP, Izmaylov AF, Bloino J, Zheng G, Sonnenberg JL, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Vreven T, Montgomery JA Jr, Peralta JE, Ogliaro F, Bearpark M, Heyd JJ, Brothers E, Kudin KN, Staroverov VN, Kobayashi R, Normand J, Raghavachari K, Rendell A, Burant JC, Iyengar SS, Tomasi J, Cossi M, Rega N, Millam JM, Klene M, Knox JE, Cross JB, Bakken V, Adamo C, Jaramillo J, Gomperts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Martin RL, Morokuma K, Zakrzewski VG, Voth GA, Salvador P, Dannenberg JJ, Dapprich S, Daniels AD, Farkas Ö, Foresman JB, Ortiz JV, Cioslowski J, Fox DJ (2009) Gaussian 09, Revision D.01. Gaussian, Inc., Wallingford, CT
Řezáč J, Riley KE, Hobza P (2011) Extensions of the S66 data set: more accurate interaction energies and angular-displaced nonequilibrium geometries. J Chem Theory Comput 7:3466–3470
Jurečka P, Šponer J, Černý J, Hobza P (2006) Benchmark database of accurate (MP2 and CCSD(T) complete basis set limit) interaction energies of small model complexes, DNA base pairs, and amino acid pairs. Phys Chem Chem Phys 8:1985–1993 (The “B” designation in the set name refers to revised BEs as computed by Marshall et al. in reference 39)
Marshall MS, Burns LA, Sherrill CD (2011) Basis set convergence of the coupled-cluster correction, \(\delta_{\text{MP2}}^{\text{CCSD(T)}}\): best practices for benchmarking non-covalent interactions and the attendant revision of the S22, NBC10, HBC6, and HSG databases. J Chem Phys 135:194102
Granatier J, Pitoňák M, Hobza P (2012) Accuracy of several wave function and density functional theory methods for description of noncovalent interaction of saturated and unsaturated hydrocarbon dimers. J Chem Theor Comput 8:2282–2292
Grimme S (2012) Supramolecular binding thermodynamics by dispersion-corrected density functional theory. Chem Euro J 18:9955–9964
Goerigk L, Grimme S (2010) A general database for main group thermochemistry, kinetics, and noncovalent interactions—assessment of common and reparameterized (meta-)GGA density functionals. J Chem Theory Comput 6:107–126
Goerigk L, Grimme S (2011) Efficient and accurate double-hybrid-meta-GGA density functionals—evaluation with the extended GMTKN30 database for general main group thermochemistry, kinetics, and noncovalent interactions. J Chem Theory Comput 7:291–309
Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Montgomery JA Jr, Vreven T, Kudin KN, Burant JC, Millam JM, Iyengar SS, Tomasi J, Barone V, Mennucci B, Cossi M, Scalmani G, Rega N, Petersson GA, Nakatsuji H, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Klene M, Li X, Knox JE, Hratchian HP, Cross JB, Adamo C, Jaramillo J, Gomperts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Ayala PY, Morokuma K, Voth GA, Salvador P, Dannenberg JJ, Zakrzewski VG, Dapprich S, Daniels AD, Strain MC, Farkas O, Malick DK, Rabuck AD, Raghavachari K, Foresman JB, Ortiz JV, Cui Q, Baboul AG, Clifford S, Cioslowski J, Stefanov BB, Liu G, Liashenko A, Piskorz P, Komaromi I, Martin RL, Fox DJ, Keith T, Al-Laham MA, Peng CY, Nanayakkara A, Challacombe M, Gill PMW, Johnson B, Chen W, Wong MW, Gonzalez C, Pople JA (2004) Gaussian 03, Revision D.01. Gaussian Inc., Pittsburgh, PA
Avila DV, Brown CE, Ingold KU, Lusztyk J (1993) Solvent effects on the competitive β-scission and hydrogen atom abstraction reactions of the cumyloxyl radical. Resolution of a long-standing problem. J Am Chem Soc 115:466–470
Montgomery JA Jr, Frisch MJ, Ochterski JW, Petersson GA (2000) A complete basis set model chemistry. VII. Use of the minimum population localization method. J Chem Phys 112:6532–6542
Johnson ER, Clarkin OJ, DiLabio GA (2003) Density functional theory based model calculations for accurate bond dissociation enthalpies. 3. a single approach for X–H, X–X, and X–Y (X, Y) C, N, O, S, halogen) bonds. J Phys Chem A 107:9953–9963
DiLabio GA, Johnson ER, Otero-de-la-Roza A (2013) Performance of conventional and dispersion-corrected density-functional theory methods for hydrogen bonding interaction energies. Phys Chem Chem Phys 15:12821–12828
Goerigk L, Kruse H, Grimme S (2011) Benchmarking density functional methods against the S66 and S66 × 8 datasets for noncovalent interactions. ChemPhysChem 12:3421–3433
Faver JC, Benson ML, He XA, Roberts BP, Wang B, Marshall MS, Kennedy MR, Sherrill CD, Merz KM (2011) Formal estimation of errors in computed absolute interaction energies of protein–ligand complexes. J Chem Theory Comput 7:790–797
Burns LA, Vázquez-Mayagoitia Á, Sumpter BG, Sherrill CD (2011) Density-functional approaches to noncovalent interactions: a comparison of dispersion corrections (DFT-D), exchange-hole dipole moment (XDM) theory, and specialized functionals. J Chem Phys 134:084107
Perdew JP, Burke K, Ernzerhof M (1996) Generalized gradient approximation made simple. Phys Rev Lett 77:3865–3868
Risthaus T, Grimme S (2013) Benchmarking of London dispersion-accounting density functional theory methods on very large molecular complexes. J Chem Theory Comput 9:1580–1591
Karton A, Tarnopolsky A, Lamére JF, Schatz GC, Martin JML (2008) Highly accurate first-principles benchmark data sets for the parametrization and thermochemical kinetics. J Phys Chem A 112:12868–12886
Guner V, Khuong KS, Leach AG, Lee PS, Bartberger MD, Houk KN (2003) A standard set of pericyclic reactions of hydrocarbons for the benchmarking of computational methods: the performance of ab initio, density functional, CASSCF, CASPT2, and CBS-QB3 methods for the prediction of activation barriers, reaction energetics, and transition state geometries. J Phys Chem A 107:11445–11459
Ess DH, Houk KN (2005) Activation energies of pericyclic reactions: performance of DFT, MP2, and CBS-QB3 methods for the prediction of activation barriers and reaction energetics of 1, 3-dipolar cycloadditions, and revised activation enthalpies for a standard set of hydrocarbon pericyclic reactions. J Phys Chem A 109:9542–9553
Dinadayalane TC, Vijaya R, Smitha A, Narahari Sastry G (2002) Diels–alder reactivity of butadiene and cyclic five-membered dienes ((CH)4X, X = CH2, SiH2, O, NH, PH, and S) with ethylene: a benchmark study. J Phys Chem A 106:1627–1633
Johnson ER, Mori-Sánchez P, Cohen AJ, Yang W (2008) Delocalization errors in density functionals and implications for main-group thermochemistry. J Chem Phys 129:204112
Krieg H, Grimme S (2010) Thermochemical benchmarking of hydrocarbon bond separation reaction energies: Jacob’s ladder is not reversed! Mol Phys 108:2655–2666
Grimme S, Steinmetz M, Korth M (2007) How to compute isomerization energies of organic molecules with quantum chemical methods. J Org Chem 72:2118–2126
Huenerbein R, Schirmer B, Moellmann J, Grimme S (2010) Effects of London dispersion on the isomerization reactions of large organic molecules: a density functional benchmark study. Phys Chem Chem Phys 12:6940–6948
Řeha D, Valdes H, Vondrášek J, Hobza P, Abu-Riziq A, Crews B, de Vries MS (2005) Structure and IR spectrum of phenylalanyl–glycyl–glycine tripeptide in the gas-phase: IR/UV experiments, Ab initio quantum chemical calculations, and molecular dynamic simulations. Chem Euro J 11:6803–6817
Gruzman D, Karton A, Martin JML (2009) Performance of Ab initio and density functional methods for conformational equilibria of CnH2n+2 alkane isomers (n = 4–8). J Phys Chem A 113:11974–11983
Csonka GI, French AD, Johnson GP, Stortz CA (2009) Evaluation of density functionals and basis sets for carbohydrates. J Chem Theory Comput 5:679–692