Easy to use DFT approach for computational pKa determination of carboxylic acids

In pKa computational determination, the challenge in exploring and fostering new methodologies and approaches goes in parallel with the amelioration of computational performances. In this paper a "ready to use methodology" has been compared to other strategies, such as the re-shaping in solvation cavity (Bondi radius re-shaping), wanting to assess its reliability in predicting the pKa of a broad list of carboxylic acids. Thus, the functionals B3LYP and CAM-B3LYP have been selected, using SMD as continuum solvation model. Exploiting our previous results, two water molecules were made explicit on the reaction centre. Data show that our model (CAM-B3LYP/2H2 O) is capable to accurately predict pKa, leading to mean absolute error (MAE) values lower than 0.5. Noteworthy, good results were achieved in computing the pKa of substituents bearing nitro and cyano groups. Focusing on B3LYP, eventually remarkable outputs were obtained only when Bondi correction was applied to the complex with two water molecules. Hence, massive outcomes were obtained in foreseeing the trichloro and trifluoro acetic acid pKa. These findings demonstrated that no complex level of theory nor external factor is required to accurately predict carboxylic acids pKa, with MAE well below 0.5 units.