Can Density Functional Theory Be Trusted for High-Order Electric Properties? The Case of Hydrogen-Bonded Complexes

This work reports on an extensive assessment of the performance of a wide palette of density functional approximations in predicting the (high-order) electric properties of hydrogen-bonded complexes. To this end, we compute the electronic and vibrational contributions to the electric polarizability...

Descripción completa

Detalles Bibliográficos
Autores: Zaleśny, Robert, Medved, Miroslav, Sitkiewicz, Sebastian P., Matito i Gras, Eduard, Luis Luis, Josep Maria
Tipo de recurso: artículo
Estado:Versión aceptada para publicación
Fecha de publicación:2019
País:España
Institución:Varias* (Consorci de Biblioteques Universitáries de Catalunya, Centre de Serveis Científics i Acadèmics de Catalunya)
Repositorio:Recercat. Dipósit de la Recerca de Catalunya
OAI Identifier:oai:recercat.cat:10256/28569
Acceso en línea:http://hdl.handle.net/10256/28569
https://hdl.handle.net/10256/28569
Access Level:acceso abierto
Palabra clave:Funcional de densitat, Teoria del
Density functionals
Descripción
Sumario:This work reports on an extensive assessment of the performance of a wide palette of density functional approximations in predicting the (high-order) electric properties of hydrogen-bonded complexes. To this end, we compute the electronic and vibrational contributions to the electric polarizability and the first and second hyperpolarizabilities, using the CCSD(T)/aug-cc-pVTZ level of theory as reference. For all the studied properties, the average absolute errors below 20% can only be obtained using the CAM-B3LYP functional, while LC-BLYP and MN15 are shown to be only slightly less accurate (average absolute errors not exceeding 30%). Among Minnesota density functionals, i.e., M06, M06-2X, and MN15, we only recommend the latter one, which quite accurately predicts the electronic and vibrational (hyper)polarizabilities. We also analyze the optimal tuning of the range-separation parameter μ for the LC-BLYP functional, finding that this approach does not bring any systematic improvement in the predictions of electronic and vibrational (hyper)polarizabilities and the accuracy of computed properties is largely system-dependent. Finally, we report huge errors in predicting the vibrational second hyperpolarizability by ωB97X, M06, and M06-2X functionals. Based on the explicit evaluation of anharmonic terms contributing to the second hyperpolarizability, this failure is traced down to a poor determination of third- and fourth-order energy derivatives with respect to normal modes. These results reveal serious flaws of some density functional approximations and suggest caution in selecting the appropriate functional to calculate not only electronic and vibrational (hyper)polarizabilities but also other molecular properties that contain vibrational anharmonic contributions