All calculations were performed using Gaussian 16 software. Geometry, frequency, and IRC analysis calculations were conducted at the M06-2X/aug-cc-pVTZ level of theory. The heats of reaction have been evaluated at the highest level, CCSD(T)/CBS(T,Q,5)//M06-2x/aug-cc-pvTz.