The bond-based approach proposed in this article lets us rely on the linear regression method. We relied on the QM9 data to build the dataset of chemical reactions. The QM9 dataset is a subset of the GDB-17 chemical universe and contains atomization energies and other properties of 133,885 small organic molecules with up to 9 heavy atoms. The geometry of the molecules of the QM9 dataset is optimized at the B3LYP/6-31G(2df,p) level of theory.