Structures of stable intermediates and transition states were obtained by standard structure optimization techniques. Frequency calculations were carried out to ensure transition states have only one imaginary frequency. IRC calculations were performed for all transition states to ensure that they connect relevant reaction intermediates. Energy decomposition analysis based on absolutely localized molecular orbitals was performed using the Q-Chem software package (version 5.3.1) .