In the calculation of reaction energies, all possible "A↔B type" isodesmic reactions were identified among the QM9 molecules, and their corresponding indices were saved in a separate CSV file named "reactions.csv". Subsequently, the energy changes associated with each reaction were computed using the same level of theory as for the energy calculations of individual molecules and stored in additional CSV files.