“…SPE calculations were carried out using the 6-311+G(d,p) basis set at B3LYP, M06-2X, and ωB97X-D functionals. 35 The polarizable continuum model 36–38 (PCM) implemented in the Gaussian 16 program package was used to incorporate the solvation (THF) effect at B3LYP, M06-2X, and ωB97X-D theory levels, the last two being hybrid long-range separated functionals, which take care of empirically-corrected dispersion interactions for studying structures with non-covalent interactions of the catalyst C 60 -substrate type. In separate sets of calculations, relative free energies were obtained using the B3LYP-D3 level of theory to incorporate dispersion correction for all the intermediates in Scheme 1A and B, the details of which are provided in Tables S1 and S2 (see ESI†).…”