2. COMPUTATIONAL DETAILS
All structures were optimized at the B3LYP-D3[34-36]/BSI level (BSI designates the basis set combination of LanL2DZ[37-39] for Rh atom, and 6-31G(d,p) for main group elements) in the gas phase. Harmonic vibrational frequencies were also calculated at the same level of theory to identify all stationary points as minima (zero imaginary frequencies) or transition states (one imaginary frequency). Intrinsic coordinate reaction (IRC)[40-42] calculations were carried out to examine the connectivity of a transition state with its backward and forward minima when necessary. The energetic results were then further refined by single-point calculations at the M06[43,44]/BSII level with solvation effects accounted for by the SMD[45-49]solvent model using ethanol or 1,4-dioxane as solvent according to the experimental conditions, where BSII denotes the basis set combination of SDD[50-52] for Rh atom, and 6-311++G(d,p) for main group elements. Natural bond orbital (NBO) analyses were performed at the B3LYP/BSII level on selected systems with the NBO code included in Gaussian 09.[53] In all of the figures that contain energy diagrams, calculated relative Gibbs free energies are presented. For reference, relative enthalpic energies are also given in parentheses. All the energies are given in kcal/mol. Unless otherwise stated, Gibbs free energies are used for the discussion of reaction mechanisms. All the calculations were performed with the Gaussian 09 software package.[53]