In order to gain a better understanding of the reaction mechanism and the origin of stereocontrol of the current pincer Ir-catalyzed asymmetric C-H insertion, density functional theory (DFT) calculations were performed by employment of Gaussian 09 program [
30]. The geometry optimization was carried out by using the M06–2X [
31,
32] functional with the basis set of 6–31G(d, p) for C, N, Cl, O, H and SDD for Ir [
33,
34] and accounting for the dichloromethane solvent effect by employing the IEF-PCM [
35,
36] solvation model (M06–2X/6–31G(d, p)+SDD//IEF-PCM
DCM). As shown in
Fig. 4, the reaction initiates with dissociation of the water from the catalyst
1g, which requires only 1.4 kcal/mol energy, giving rise to the formation of the 16-electron intermediate
trans-
5 with the two chloride ligands located
trans to each other. Then the resulting
trans-
5 intermediate would interact with the 3-diazooxindole to generate the iridium carbenoid complex along with extrusion of nitrogen. The energy barrier associated with this process (through the transition state
8-TS) is found to be 25.9 kcal/mol, which is difficult to attain under the optimized reaction conditions. Meanwhile,
trans-
5 can rearrange to its
cis counterpart
cis-
5 through
6-TS with a much lower energy barrier of 17.5 kcal/mol. Furthermore, the interaction of
cis-
5 intermediate with the 3-diazooxindole proceeds very readily
via the
7-TS because the
7-TS has a lower energy than
cis-
5, furnishing the iridium carbenoid intermediate
9 with the carbene ligand coordinated at the axial position. Subsequently, the carbenoid intermediate
9 undergoes the C-H insertion step with 1, 4-cyclohexadiene to give the product
3a. It was found that the reaction did not proceed
via a three-centered transition state which would yield directly the insertion product
3a. Instead, a stepwise process involving hydride transfer followed by C-C bond formation was proposed computationally (
Fig. 4). This is quite different from the related Ir and Rh catalyzed asymmetric C-H insertion with cyclohexadiene, where a concerted or concerted asynchronous mechanism involving a three-centered transition state was well supported by computational analysis [
37,
38]. In the hydride transfer step,
Si-face approach of the cyclodiene to the intermediate
9 through the
10-TS has an energy barrier of 15.3 kcal/mol, whereas a higher barrier of 20.2 kcal/mol (
11-TS) is present for the corresponding
Re-face approach. The higher energy of the
11-TS is related to the steric hindrance of upward 4
S-phenyl substituent on the imidazoline ligand, which shields the
Re-face of the carbenoid intermediate, thereby making the
Si-face approach of the cyclodiene preferentially and leading to the formation of (
S)-isomer of the product with high enantioselectivity. In addition, noncovalent interaction (NCI) [
39,
40] was used to analyze the transition states
10-TS and
11-TS (
Fig. 5). It is found that there are five stronger interactions (C–H…Cl, C–H…
π, C–H…Cl and two
π…
π) in
10-TS. While only four interactions (C–H…
π, C–H…Cl, N–H…
π and
π…
π) exist in
11-TS. The observation also indicates that
10-TS is more stable than
11-TS. The lower energy
10-TS results in generation of the zwitterionic intermediate
12. From this intermediate, the reaction continues and the C-C bond is formed to afford the product (
S)
−3a via the
14-TS with an energy barrier of 10.7 kcal/mol.