Quantum Chemical Topology Group, University of Manchester, UK
E-mail: paul.popelier@manchester.ac.uk
The Relative Energy Gradient (REG) Method or
How to Compute Chemical Insight
Quantum Chemical Topology Group, University of Manchester, UK
E-mail: paul.popelier@manchester.ac.uk
Much chemical insight ultimately comes down to finding out which fragment of a total system behaves like the total system. What is needed is a geometrical change exposing a chemical phenomenon (e.g. rotation barrier) at hand thereby generating an energy profile governed by a control coordinate.
A simple example is that of the water dimer. It is said that this system is held together by a hydrogen bond, and the distance between the two waters can serve as a control coordinate. In this case the fragment consists of two atoms: H···O making up the hydrogen bond. However, from a quantum mechanical point of view, each atom in the dimer interacts with any other atom. Thus, the view that the hydrogen bond by itself governs the energetic stability of the water dimer thus needs rigorous scrutiny.
We propose a solution for this general and important problem by presenting the Relative Energy Gradient (REG) method1. This automatable and unbiased method operates on energy contributions calculated by the Interaction Quantum Atoms (IQA) method, which belongs to the real-space partitioning framework of Quantum Chemical Topology2, 3.
The REG method has been implemented4 in the in-house program REG.py5, and also in the Amsterdam Density Functional (ADF) package. REG is able to explain many chemical phenomena. To name some: the gauche effect6, the torsional barrier in biphenyl7 (which was long controversial, see Figure), the “arrow-pushing” scheme behind an enzymatic reaction (peptide hydrolysis in the HIV-1 Protease active site)8, halogen-alkane nucleophilic substitution (SN2) reactions9, and the factors controlling the nature of complementary hydrogen-bonded networks as found in nucleic acid base pairs10.
References
[1] J. C. R. Thacker and P. L. A. Popelier, Theor.Chem.Acc., 2017, 136, 86.
[2] P. L. A. Popelier, in Challenges and Advances in Computational Chemistry and Physics dedicated to "Applications of Topological Methods in Molecular Chemistry", eds. R. Chauvin, C. Lepetit, E. Alikhani and B. Silvi, Springer, Switzerland, 2016, ch. 2, pp. 23-52.
[3] P. L. A. Popelier, in The Chemical Bond - 100 years old and getting stronger, ed. M. Mingos, Springer, Switzerland, 2016, ch. 2, pp. 71-117.
[4] F. Falcioni and P. L. A. Popelier, J.Chem.Inf.Mod., 2023, 63, 4312-4327.
[5] F. Falcioni, L. J. Duarte and P. L. A. Popelier, REG.py (Version 0.1), https://github.com/popelier-group/REG).
[6] J. C. R. Thacker and P. L. A. Popelier, J.Phys.Chem.A, 2018, 122, 1439−1450.
[7] P. L. A. Popelier, P. I. Maxwell, J. C. R. Thacker and I. Alkorta, Theor.Chem.Accs., 2019, 138:12, 1-16.
[8] J. C. R. Thacker, M. A. Vincent and P. L. A. Popelier, Chem.Eur.J., 2018, 14, 11200-11210.
[9] I. Alkorta, J. C. R. Thacker and P. L. A. Popelier, J.Comp.Chem. , 2018, 39, 546–556.
[10] O. J. Backhouse, J. C. R. Thacker and P. L. A. Popelier, ChemPhysChem, 2019, 20, 555-564.