Dataset Information

Critical Assessment of the Performance of Density Functional Methods for Several Atomic and Molecular Properties.

ABSTRACT: The reliable prediction of molecular properties is a vital task of computational chemistry. In recent years, density functional theory (DFT) has become a popular method for calculating molecular properties for a vast array of systems varying in size from small organic molecules to large biological compounds such as proteins. In this work we assess the ability of many DFT methods to accurately determine atomic and molecular properties for small molecules containing elements commonly found in proteins, DNA, and RNA. These properties include bond lengths, bond angles, ground state vibrational frequencies, electron affinities, ionization potentials, heats of formation, hydrogen bond interaction energies, conformational energies, and reaction barrier heights. Calculations are carried out with the 3-21G*, 6-31G*, 3-21+G*, 6-31+G*, 6-31++G*, cc-pVxZ, and aug-cc-pVxZ (x=D,T) basis sets, while bond distance and bond angle calculations are also done using the cc-pVQZ and aug-cc-pVQZ basis sets. Members of the popular functional classes, namely, LSDA, GGA, meta-GGA, hybrid-GGA, and hybrid-meta-GGA, are considered in this work. For the purpose of comparison, Hartree-Fock (HF) and second order many-body perturbation (MP2) methods are also assessed in terms of their ability to determine these physical properties. Ultimately, it is observed that the split valence bases of the 6-31G variety provide accuracies similar to those of the more computationally expensive Dunning type basis sets. Another conclusion from this survey is that the hybrid-meta-GGA functionals are typically among the most accurate functionals for all of the properties examined in this work.

SUBMITTER: Riley KE

PROVIDER: S-EPMC2581803 | biostudies-literature | 2007

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Critical Assessment of the Performance of Density Functional Methods for Several Atomic and Molecular Properties.

Riley Kevin E KE Op't Holt Bryan T BT Merz Kenneth M KM

Journal of chemical theory and computation 20070101 2

The reliable prediction of molecular properties is a vital task of computational chemistry. In recent years, density functional theory (DFT) has become a popular method for calculating molecular properties for a vast array of systems varying in size from small organic molecules to large biological compounds such as proteins. In this work we assess the ability of many DFT methods to accurately determine atomic and molecular properties for small molecules containing elements commonly found in prot ...[more]

PMID: 19002267

Similar Datasets

Project description:How to describe nondynamic electron correlation is still a major challenge to density functional theory (DFT). Recent models designed particularly for this problem, such as Becke'05 (B05) and Perdew-Staroverov-Tao-Scuseria (PSTS) functionals employ the exact-exchange density, the efficient calculation of which is technically quite challenging. We have recently implemented self-consistently the B05 functional based on an efficient resolution-identity (RI) technique. In this study, we report a self-consistent RI implementation of the PSTS functional. In contrast to its original implementation, our version brings no limitation on the choice of the basis set. We have also implemented the Mori-Sanchez-Cohen-Yang-2 (MCY2) functional, another recent DFT method that includes full exact exchange. The performance of PSTS, B05, and MCY2 is validated on thermochemistry, reaction barriers, and dissociation energy curves, with an emphasis on nondynamic correlation effects in the discussion. All three methods perform rather well in general, B05 and MCY2 being on average somewhat better than PSTS. We include also results with other functionals that represent various aspects of the development in this field in recent years, including B3LYP, M06-HF, M06-2X, ωB97X, and TPSSh. The performance of the heavy-parameterized functionals M06-2X and ωB97X is on average better than that of B05, MCY2, and PSTS for standard thermodynamic properties and reactions, while the latter functionals do better in hydrogen abstraction reactions and dissociation processes. In particular, B05 is found to be the only functional that yields qualitatively correct dissociation curves for two-center symmetric radicals like He(2)(+). Finally, we compare the performance of all these functionals on a strongly correlated exemplary case system, the NO dimer. Only PSTS, B05, and MCY2 describe the system qualitatively correctly. Overall, this new type of functionals show good promise of overcoming some of the difficulties DFT encounters for systems with strong nondynamic correlation.

Project description:Quantum chemical calculations are important for elucidating light-capturing mechanisms in photobiological systems. The time-dependent density functional theory (TDDFT) has become a popular methodology because of its balance between accuracy and computational scaling, despite its problems in describing, for example, charge transfer states. As a step toward systematically understanding the performance of TDDFT calculations on biomolecular systems, we study here 17 commonly used density functionals, including seven long-range separated functionals, and compare the obtained results with excitation energies calculated at the approximate second order coupled-cluster theory level (CC2). The benchmarking set includes the first five singlet excited states of 11 chemical analogues of biochromophores from the green fluorescent protein, rhodopsin/bacteriorhodopsin (Rh/bR), and the photoactive yellow protein. We find that commonly used pure density functionals such as BP86, PBE, M11-L, and hybrid functionals with 20-25% of Hartree-Fock (HF) exchange (B3LYP, PBE0) have a tendency to consistently underestimate vertical excitation energies (VEEs) relative to the CC2 values, whereas hybrid density functionals with around 50% HF exchange such as BHLYP, PBE50, and M06-2X and long-range corrected functionals such as CAM-B3LYP, ωPBE, ωPBEh, ωB97X, ωB97XD, BNL, and M11 overestimate the VEEs. We observe that calculations using the CAM-B3LYP and ωPBEh functionals with 65% and 100% long-range HF exchange, respectively, lead to an overestimation of the VEEs by 0.2-0.3 eV for the benchmarking set. To reduce the systematic error, we introduce here two new empirical functionals, CAMh-B3LYP and ωhPBE0, for which we adjusted the long-range HF exchange to 50%. The introduced parameterization reduces the mean signed average (MSA) deviation to 0.07 eV and the root mean square (rms) deviation to 0.17 eV as compared to the CC2 values. In the present study, TDDFT calculations using the aug-def2-TZVP basis sets, the best performing functionals relative to CC2 are ωhPBE0 (rms = 0.17, MSA = 0.06 eV); CAMh-B3LYP (rms = 0.16, MSA = 0.07 eV); and PBE0 (rms = 0.23, MSA = -0.14 eV). For the popular range-separated CAM-B3LYP functional, we obtain an rms value of 0.31 eV and an MSA value of 0.25 eV, which can be compared with the rms and MSA values of 0.37 and -0.31 eV, respectively, as obtained at the B3LYP level.

Dataset Information

Critical Assessment of the Performance of Density Functional Methods for Several Atomic and Molecular Properties.

Publications

Critical Assessment of the Performance of Density Functional Methods for Several Atomic and Molecular Properties.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure