In the present study, YARP used Gaussian 16 as the reference quantum chemistry engine for the DFT calculations associated with the Berny optimizations and IRC calculations. Calculations were performed at the B3LYP/6-31G level of theory during network exploration, while reaction energies and reaction barriers were refined at the B3LYP-D3/6-311G(d,p) level of theory for validation and comparison with periodic calculations. The GSM calculations were performed by the pyGSM package using eleven images, fixed reactant and product geometries, and other default hyperparameters. All GFN2-xTB calculations were performed with the xTB program (version 6.2.3). Periodic DFT calculations were performed using Vienna Ab-initio Simulation Package (VASP, 5.4.1), where plane-wave basis sets expanded the Kohn-Sham orbitals, and the Kohn-Sham equations were solved self-consistently. The BEEF-VdW exchange-correlation functional with pro-jector augmented wave (PAW) pseudopotentials was employed. A Monkhorst-Pack k -sampling was used, and a k point grid of 2 × 2 × 1 was applied. A cutoff energy of 400 eV and a forceconvergence criterion of 20 meV Å-1 for energy local minima were used. The climbing image nudged-elastic-band (CI-NEB) method was used as a first step to locate transition states. Seven images were used in each NEB calculation as generated by the Image Dependent Pair Potential (IDPP) tool. Following each NEB calculation, Lanczos diagonalization was used to identify the transition state with a greater accuracy. The force-convergence criterion of a transition state optimization was 20 meV Å-1 . All energies are reported with respect to the ground state energy of Ga-ethyl plus a gaseous ethylene molecule.