Issue 19, 2020

Rapid and accurate molecular deprotonation energies from quantum alchemy


We assess the applicability of alchemical perturbation density functional theory (APDFT) for quickly and accurately estimating deprotonation energies. We have considered all possible single and double deprotonations in one hundred small organic molecules drawn at random from QM9 [Ramakrishnan et al., JCTC, 2015]. Numerical evidence is presented for 5160 deprotonated species at both HF/def2-TZVP and CCSD/6-31G* levels of theory. We show that the perturbation expansion formalism of APDFT quickly converges to reliable results: using CCSD electron densities and derivatives, regular Hartree–Fock calculations are outperformed at the second or third order for ranking all possible doubly or singly deprotonated molecules, respectively. CCSD single deprotonation energies are reproduced within 1.4 kcal mol−1 on average within third order APDFT. We introduce a hybrid approach where the computational cost of APDFT is reduced even further by mixing first order terms at a higher level of theory (CCSD) with higher order terms at a lower level of theory only (HF). We find that this approach reaches 2 kcal mol−1 accuracy in absolute deprotonation energies compared to CCSD at 2% of the computational cost of third order APDFT.

Graphical abstract: Rapid and accurate molecular deprotonation energies from quantum alchemy

Supplementary files

Article information

Article type
29 Nov 2019
19 Dec 2019
First published
19 Dec 2019
This article is Open Access
Creative Commons BY license

Phys. Chem. Chem. Phys., 2020,22, 10519-10525

Rapid and accurate molecular deprotonation energies from quantum alchemy

G. F. von Rudorff and O. A. von Lilienfeld, Phys. Chem. Chem. Phys., 2020, 22, 10519 DOI: 10.1039/C9CP06471K

This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. You can use material from this article in other publications without requesting further permissions from the RSC, provided that the correct acknowledgement is given.

Read more about how to correctly acknowledge RSC content.

Social activity