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]