On the Accuracy of the Direct Method to Calculate pK(a) from Electronic Structure Calculations

Autor(es)

Dutra, Felipe Ribeiro ^[1] ; Silva, Cleuton de Souza ^[2] ; Custodio, Rogerio ^[1]

Afiliação do(s) autor(es)

^[1] Univ Estadual Campinas, Inst Quim, BR-13083970 Campinas, SP - Brazil

^[2] Univ Fed Amazonas, Inst Ciencias Exatas & Tecnol, BR-69100021 Itacoatiara, Amazonas - Brazil

• Tipo de documento: Artigo Científico
• Fonte: Journal of Physical Chemistry A; v. 125, n. 1, p. 65-73, JAN 14 2021.
• Citações Web of Science: 0

Resumo

The direct method (HA((soln)) reversible arrow A((soln))(-) + H-(so(ln))+) for calculating pK(a) of monoprotic acids is as efficient as thermodynamic cycles. A selective adjustment of proton free energy in solution was used with experimental pK(a) data. The procedure was analyzed at different levels of theory. The solvent was described by the solvation model density (SMD) model, including or not explicit water molecules, and three training sets were tested. The best performance under any condition was obtained by the G4CEP method with a mean absolute error close to 0.5 units of pK(a) and an uncertainty around +/- 1 unit of pK(a) for any training set including or excluding explicit solvent molecules. PM6 and AM1 performed very well with average absolute errors below 0.75 units of pK(a) but with uncertainties up to +/- 2 units of pK(a,) using only the SMD solvent model. Density functional theory (DFT) results were highly dependent on the basis functions and explicit water molecules. The best performance was observed for the local spin density approximation (LSDA) functional in almost all calculations and under certain conditions, as high as those obtained by G4CEP. Basis set complexity and explicit solvent molecules were important factors to control DFT calculations. The training set molecules should consider the diversity of compounds. (AU)

• Processo FAPESP: 13/08293-7 - CECC - Centro de Engenharia e Ciências Computacionais
• Beneficiário: Munir Salomao Skaf
• Modalidade de apoio: Auxílio à Pesquisa - Centros de Pesquisa, Inovação e Difusão - CEPIDs
• Processo FAPESP: 17/11485-6 - Desenvolvimentos computacionais e teóricos baseados em métodos ab initio e na Teoria do Funcional de Densidade
• Beneficiário: Rogério Custodio
• Modalidade de apoio: Auxílio à Pesquisa - Regular