Role of Electronic and Steric Effects on Ruthenium Catalysts with Bulky NHC Ligands and Relationship with the Z-Selectivity in Oleﬁn Metathesis

: Recently, sterically demanding N-heterocyclic cyclometalated ruthenium were reported as efﬁcient Z-selective catalysts for cross-metathesis, showing a different reactivity in the function of the auxiliary ligand and the bulky ligand. To understand the origin of this behavior, we carried out density functional (M06-L) calculations to explore the reaction mechanism and insight from the energetic contributions into the determinant step. We emphasize the differences that occur when the 2,6-diisopropylphenyl (Dipp) and 2,6-diisopentylphenyl (Dipep) are employed. The results show that the barrier energies, ∆ G ‡ , increase when the bulky ligand is greater, using nitrate as an auxiliary ligand, while the opposite behavior is obtained when pivalate is the auxiliary ligand. This tendency has its origin in the low reorganization energy and the less steric hindrance (%V bur ) obtained in catalysts that involve nitrate ligand and Dipep group. Moreover, by scrutinizing the energy decomposition analysis (EDA), it is found that the electronic contributions are also dominant and are not uniquely the steric effects that control the Z-selectivity.

Other Keitz's results showed higher activity and stability compared with substituted carboxylate when the anionic group is substituted by nitrate (Ru-2) [17].Several attempts to improve the catalytic activity, selectivity, and stability have been realized through the design of new catalysts such as those described in Figure 2.For instance, the substitution of the O-bidentate by an S-bidentate anionic group, catecholate type, (Ru-3) [18] and the chlorine addition to the catecholate group (Ru-4) [19], both without the bulky ligand Adm, gave origin to catalysts with high selectivity (>98%) in the CM of allylic alcohols.Also, it is important to highlight the challenging proposal to include chelating laterals with a smaller volume than the Adm group (Ru-5) [20].Recently, the ab initio design proposed and synthesized the first S-monodentate catalyst (Ru-6) [21], reaching up to 80% in prototypical homocoupling metathesis reactions.Thus, the nature of the anionic group, different electron acceptor residues, and the size of the chelate ligand would influence the origin of the Z-selectivity, suggesting that the steric effect is not uniquely responsible.Other Keitz's results showed higher activity and stability compared with substituted carboxylate when the anionic group is substituted by nitrate (Ru-2) [17].Several attempts to improve the catalytic activity, selectivity, and stability have been realized through the design of new catalysts such as those described in Figure 2.For instance, the substitution of the O-bidentate by an S-bidentate anionic group, catecholate type, (Ru-3) [18] and the chlorine addition to the catecholate group (Ru-4) [19], both without the bulky ligand Adm, gave origin to catalysts with high selectivity (>98%) in the CM of allylic alcohols.Also, it is important to highlight the challenging proposal to include chelating laterals with a smaller volume than the Adm group (Ru-5) [20].Recently, the ab initio design proposed and synthesized the first S-monodentate catalyst (Ru-6) [21], reaching up to 80% in prototypical homocoupling metathesis reactions.Thus, the nature of the anionic group, different electron acceptor residues, and the size of the chelate ligand would influence the origin of the Z-selectivity, suggesting that the steric effect is not uniquely responsible.Recently, Grubbs reported the preparation of sterically demanding N-heterocyclic cyclometalated ruthenium catalysts with applications in CM between acrylamides and terminal olefins [22].Historically, the acrylamide-olefin metathesis reaction (Figure 3a) using Ru-1 and Ru-2 was difficult due to the high thermodynamic stability towards Eselectivity (>25:1) and the long reaction times exceeding 48 h [23].To overcome this diffi-  Other Keitz's results showed higher activity and stability compared with substituted carboxylate when the anionic group is substituted by nitrate (Ru-2) [17].Several attempts to improve the catalytic activity, selectivity, and stability have been realized through the design of new catalysts such as those described in Figure 2.For instance, the substitution of the O-bidentate by an S-bidentate anionic group, catecholate type, (Ru-3) [18] and the chlorine addition to the catecholate group (Ru-4) [19], both without the bulky ligand Adm, gave origin to catalysts with high selectivity (>98%) in the CM of allylic alcohols.Also, it is important to highlight the challenging proposal to include chelating laterals with a smaller volume than the Adm group (Ru-5) [20].Recently, the ab initio design proposed and synthesized the first S-monodentate catalyst (Ru-6) [21], reaching up to 80% in prototypical homocoupling metathesis reactions.Thus, the nature of the anionic group, different electron acceptor residues, and the size of the chelate ligand would influence the origin of the Z-selectivity, suggesting that the steric effect is not uniquely responsible.Recently, Grubbs reported the preparation of sterically demanding N-heterocyclic cyclometalated ruthenium catalysts with applications in CM between acrylamides and terminal olefins [22].Historically, the acrylamide-olefin metathesis reaction (Figure 3a) using Ru-1 and Ru-2 was difficult due to the high thermodynamic stability towards Eselectivity (>25:1) and the long reaction times exceeding 48 h [23].To overcome this difficulty, proposed catalysts were derived from the incorporation of an even bulkier N-substituent than the Mes, such as the 2,6-diisopropylphenyl (Dipp) [24,25] and 2,6-diisopentylphenyl (Dipep) [26] groups, as is described in Figure 3b for the complexes Ru-7 to Ru-8 and Ru-9 to Ru-10, respectively.In this work, when Ru-7 was employed, a Z-amide and Z-olefin were obtained with a yield of 27% and 63%, respectively; if the anionic group was Recently, Grubbs reported the preparation of sterically demanding N-heterocyclic cyclometalated ruthenium catalysts with applications in CM between acrylamides and terminal olefins [22].Historically, the acrylamide-olefin metathesis reaction (Figure 3a) using Ru-1 and Ru-2 was difficult due to the high thermodynamic stability towards Eselectivity (>25:1) and the long reaction times exceeding 48 h [23].To overcome this difficulty, proposed catalysts were derived from the incorporation of an even bulkier N-substituent than the Mes, such as the 2,6-diisopropylphenyl (Dipp) [24,25] and 2,6diisopentylphenyl (Dipep) [26] groups, as is described in Figure 3b for the complexes Ru-7 to Ru-8 and Ru-9 to Ru-10, respectively.In this work, when Ru-7 was employed, a Z-amide and Z-olefin were obtained with a yield of 27% and 63%, respectively; if the anionic group was exchanged to carboxylate (Ru-8), the yield increased remarkably up to 49% to the Z-amide and decreased to 11% to the olefin.In turn, this limit was exceeded by the Ru-10 system, which produced the maximum yield (80%) for Z-amide and a new increase of 27% for Z-olefin.Thus, the reactivity of the CM can be significantly influenced by the combination of the auxiliary ligand nature and the bulkier N-substituent NHC ligand.
Hence, gaining knowledge of these contributions can facilitate a broader comprehension of olefin metathesis.
exchanged to carboxylate (Ru-8), the yield increased remarkably up to 49% to the Z-amide and decreased to 11% to the olefin.In turn, this limit was exceeded by the Ru-10 system, which produced the maximum yield (80%) for Z-amide and a new increase of 27% for Zolefin.Thus, the reactivity of the CM can be significantly influenced by the combination of the auxiliary ligand nature and the bulkier N-substituent NHC ligand.Hence, gaining knowledge of these contributions can facilitate a broader comprehension of olefin metathesis.The experimental and theoretical description of the reaction mechanism in olefin metathesis has been a challenge [27][28][29][30][31][32][33].The generally accepted mechanism, depicted in Figure 4, begins with the activation of the 16-electron pre-catalyst through the dissociation reaction of the Ru-O bond originating from the isopropyl, resulting in the formation of a 14-electron Ru-alkylidene active intermediate .
The formation of the active intermediate gives rise to the next stage, which entails the coordination of the olefinic substrate in two coordination modes: a) trans to NHC (called bottom pathway or b) trans to bidentate anionic group (called side pathway).In both pathways, metallocyclobutane (MC) is produced as a result of [2+2] cycloaddition (TSCA) being surmounted.In the following step of the isomerization of the MC intermediate, a [2+2] cycloreversion (TSCR) gives rise to the decoordination of the new olefin.These pathways are crucial for the E-/Z-selectivity.If methyl substituents are located behind (B) or in front of (F) to the metallocycle plane on the bottom pathway, an E-2-butene olefin is obtained, whereas if they are located up (U) and down (D) on the side pathway, a Z-2-butene alkene is produced.Figure 4 depicts the main transition states of each substage, along with their respective abbreviation.
The origin of Z-selectivity in olefin metathesis has been the focus of a variety of studies in recent years.In this field, Liu [32] reported that the side preference for Ru-1 can be explained by steric and electronic factors of the different ligands that influence the metal center's stability.According to the results, the cis products are driven by the side mechanism due to a kinetic preference.The Gibbs energy difference between the most favorable mechanism (side) and the less favorable pathway (bottom) was reported at 10.4 kcal/mol (∆∆ ).Within this mechanism, the substituents originating from the olefin acquire a spatial orientation opposite to that of the bulky Mes.In terms of the electronic effect, the transition states of the "bottom" mechanism are destabilized due to the backdonation from π* orbitals of both NHC and alkylidene to two different types of d orbitals on Ruthenium, whereas in the "side" mechanism, the retro donation involves the same d orbital, resulting in the stability.Wang and Dang [34] found comparable results for Ru-2.The experimental and theoretical description of the reaction mechanism in olefin metathesis has been a challenge [27][28][29][30][31][32][33].The generally accepted mechanism, depicted in Figure 4, begins with the activation of the 16-electron pre-catalyst through the dissociation reaction of the Ru-O bond originating from the isopropyl, resulting in the formation of a 14-electron Ru-alkylidene active intermediate .
The formation of the active intermediate gives rise to the next stage, which entails the coordination of the olefinic substrate in two coordination modes: a) trans to NHC (called bottom pathway or b) trans to bidentate anionic group (called side pathway).In both pathways, metallocyclobutane (MC) is produced as a result of [2+2] cycloaddition (TSCA) being surmounted.In the following step of the isomerization of the MC intermediate, a [2+2] cycloreversion (TSCR) gives rise to the decoordination of the new olefin.These pathways are crucial for the E-/Zselectivity.If methyl substituents are located behind (B) or in front of (F) to the metallocycle plane on the bottom pathway, an E-2-butene olefin is obtained, whereas if they are located up (U) and down (D) on the side pathway, a Z-2-butene alkene is produced.Figure 4 depicts the main transition states of each substage, along with their respective abbreviation.
The origin of Z-selectivity in olefin metathesis has been the focus of a variety of studies in recent years.In this field, Liu [32] reported that the side preference for Ru-1 can be explained by steric and electronic factors of the different ligands that influence the metal center's stability.According to the results, the cis products are driven by the side mechanism due to a kinetic preference.The Gibbs energy difference between the most favorable mechanism (side) and the less favorable pathway (bottom) was reported at 10.4 kcal/mol (∆∆G Bottom−Side ).Within this mechanism, the substituents originating from the olefin acquire a spatial orientation opposite to that of the bulky Mes.In terms of the electronic effect, the transition states of the "bottom" mechanism are destabilized due to the back-donation from π* orbitals of both NHC and alkylidene to two different types of d orbitals on Ruthenium, whereas in the "side" mechanism, the retro donation involves the same d orbital, resulting in the stability.Wang and Dang [34] found comparable results for Ru-2.
In addition, Hoveyda and co-workers have conducted an analysis of other catalysts, such as the Ru-3 [18].The researchers studied four possible reaction pathways, including styrene as substrate, and determined the steric effect by analyzing the spatial orientation of the carbene substituent in relation to the NHC.In the most favorable pathway, the geometry of the transition state for the [2+2] cycloaddition reveals that the Ru carbene substituent is oriented away from the NHC, exhibiting a barrier of 3.2 kcal/mol.During the exploration of the subsequent, more stable route, a steric hindrance between the species was identified, with an energetic value of 6.4 kcal/mol.Therefore, this work highlighted that the presence of a significant steric repulsion in the metallocycle plays a crucial role in promoting the preferential formation of Z-selectivity.
In addition, Hoveyda and co-workers have conducted an analysis of other catalysts, such as the Ru-3 [18].The researchers studied four possible reaction pathways, including styrene as substrate, and determined the steric effect by analyzing the spatial orientation of the carbene substituent in relation to the NHC.In the most favorable pathway, the geometry of the transition state for the [2+2] cycloaddition reveals that the Ru carbene substituent is oriented away from the NHC, exhibiting a barrier of 3.2 kcal/mol.During the exploration of the subsequent, more stable route, a steric hindrance between the species was identified, with an energetic value of 6.4 kcal/mol.Therefore, this work highlighted that the presence of a significant steric repulsion in the metallocycle plays a crucial role in promoting the preferential formation of Z-selectivity.
Also, a modified Ru-3 catalyst [37] with Dipp ligand was analyzed, suggesting that this substituent would lead to a greater steric effect, thereby favoring a suitable spatial position in the MCy for Z-selectivity.Martinez [38] studied the initial stage of the Ru-3 and ethylene-based olefin metathesis to delve deeper into this subject.In the initial stage, by altering the dihedral angle around styrene, eight structures with distinct positions and orientations were discovered.The energetic values of the structures evaluated in the initiation stage range from 0 to 22 kcal/mol.
The most important result is the dissociation Ru-S bond, which suggests that the Smonodentate promotes the subsequent stages of olefin metathesis.In the mechanism, they evaluated bottom and side pathways and quantified the steric and electronic properties using the activation strain model.They discovered that the Z-isomer can be formed by kinetic control and suggest that the metallocycle can be stabilized by the electronic effect of the active species and a reduction in the olefin´s strain.These works demonstrate the relevance of understanding the factors that govern Z-selectivity in olefin metathesis.
As a result of these findings, it is essential to identify, comprehend, and investigate the steric and electronic effects that control the Z-selectivity in ruthenium catalysts.Experiments have shown that changing the auxiliary ligand and increasing the steric volume around the central may have a considerable impact on selectivity [22].In this direction, we explore the ethane cross-metathesis reaction mechanism and understand the energetic contributions (steric and electronic) to the rate-determining step, emphasizing the differences between 2,6-diisopropylphenyl (Dipp) and 2,6-diisopentylphenyl (Dipep).Therefore, the dominance of steric-electronic factors in the activation energies of the Z-selective Also, a modified Ru-3 catalyst [37] with Dipp ligand was analyzed, suggesting that this substituent would lead to a greater steric effect, thereby favoring a suitable spatial position in the MCy for Z-selectivity.Martinez [38] studied the initial stage of the Ru-3 and ethylene-based olefin metathesis to delve deeper into this subject.In the initial stage, by altering the dihedral angle around styrene, eight structures with distinct positions and orientations were discovered.The energetic values of the structures evaluated in the initiation stage range from 0 to 22 kcal/mol.
The most important result is the dissociation Ru-S bond, which suggests that the S-monodentate promotes the subsequent stages of olefin metathesis.In the mechanism, they evaluated bottom and side pathways and quantified the steric and electronic properties using the activation strain model.They discovered that the Z-isomer can be formed by kinetic control and suggest that the metallocycle can be stabilized by the electronic effect of the active species and a reduction in the olefin s strain.These works demonstrate the relevance of understanding the factors that govern Z-selectivity in olefin metathesis.
As a result of these findings, it is essential to identify, comprehend, and investigate the steric and electronic effects that control the Z-selectivity in ruthenium catalysts.Experiments have shown that changing the auxiliary ligand and increasing the steric volume around the central may have a considerable impact on selectivity [22].In this direction, we explore the ethane cross-metathesis reaction mechanism and understand the energetic contributions (steric and electronic) to the rate-determining step, emphasizing the differences between 2,6-diisopropylphenyl (Dipp) and 2,6-diisopentylphenyl (Dipep).Therefore, the dominance of steric-electronic factors in the activation energies of the Z-selective represents a new limit to be overcome, allowing advancement in the design of more specialized catalysts.

Bottom and Side Mechanism Employing the Keitz-Grubbs-Substituted, Ru-1 and Ru-2 Catalysts
Figures S1 and S2 depict the energy profiles for the CM reaction with Ru-1 and Ru-2 catalysts, respectively.In both instances, the Side-DD pathway is observed to be more stable than the Bottom-BF.The limiting stage or Ru-1 is the TSCR with barrier energies of 20.3 kcal/mol and 32.4 kcal/mol, respectively.Consequently, the barrier Gibbs energy ∆∆G Bottom-Side difference for Ru-1 and Ru-2 are 12.1 kcal/mol and 18.4 kcal/mol, respectively.Thus, our methodology can reproduce the main Z-product formation in these catalysts as observed experimentally.In addition, these results are consistent with the energy gap ∆∆G Bottom-Side closest to 10 kcal/mol obtained by Liu et al., [32] which demon-strates that a kinetic control is essential for Z-selectivity.Figures S3 and S4 depict the optimized intermediate and transition structures for the side pathway of each catalyst.
In addition, an energy decomposition analysis (EDA) was conducted.Table 1 describes the reorganization and interaction contribution energies of the limiting stages for Ru-1 and Ru-2.In general, the reorganization contribution ( ∆E reorg t in Ru-1 is significantly higher in Bottom-BF than Side-DD, with respective values of 50.9 kcal/mol and 41.5 kcal/mol.Similar results were discovered for Ru-2.Thus, the side pathway requires 7-8 kcal/mol less energy than the bottom pathway to reorganize.Regarding the interaction contribution ∆E int , the obtained energy values indicate that the Side-DD mechanism is more stable than in Bottom-BF.Comparing these contributions energies in both catalysts (∆E reorg t versus ∆E int ), we find that the interaction energy is greater in the side mechanism than the reorganization energy.Consequently, it is reported that in the limiting stage of the Side-DD mechanism, the electronic effects predominate over the steric ones, whereas in Bottom-BF, the converse is true.Thus, the electronic contributions are accountable for the Z-selectivity's stability.Thus, our methodology can reproduce the main Z-product formation in these catalysts as observed experimentally.In addition, these results are consistent with the energy gap ∆∆GBottom-Side closest to 10 kcal/mol obtained by Liu et al., [32] which demonstrates that a kinetic control is essential for Z-selectivity.Figures S3 and S4 depict the optimized intermediate and transition structures for the side pathway of each catalyst.
In addition, an energy decomposition analysis (EDA) was conducted.Table 1 describes the reorganization and interaction contribution energies of the limiting stages for Ru-1 and Ru-2.In general, the reorganization contribution (∆ in Ru-1 is significantly higher in Bottom-BF than Side-DD, with respective values of 50.9 kcal/mol and 41.5 kcal/mol.Similar results were discovered for Ru-2.Thus, the side pathway requires 7-8 kcal/mol less energy than the bottom pathway to reorganize.Regarding the interaction contribution ∆ , the obtained energy values indicate that the Side-DD mechanism is more stable than in Bottom-BF.Comparing these contributions energies in both catalysts (∆  ∆ ), we find that the interaction energy is greater in the side mechanism than the reorganization energy.Consequently, it is reported that in the limiting stage of the Side-DD mechanism, the electronic effects predominate over the steric ones, whereas in Bottom-BF, the converse is true.Thus, the electronic contributions are accountable for the Z-selectivity's stability.In Tables 2 and 3, the energies of intermediates and transition states associated with the Bottom and Side pathways of the CM utilizing Ru-7 to Ru-10 catalysts are listed.The Side-DD demonstrates low energy barriers, indicating that thermal equilibrium could be readily achieved in all catalysts.In particular, the TSCR with Gibbs energies (∆G ‡ ) between 16. stable than the Bottom-BF.The limiting stage or Ru-1 is the TSCR with barrier energies of 20.3 kcal/mol and 32.4 kcal/mol, respectively.Consequently, the barrier Gibbs energy ∆∆GBottom-Side difference for Ru-1 and Ru-2 are 12.1 kcal/mol and 18.4 kcal/mol, respectively.Thus, our methodology can reproduce the main Z-product formation in these catalysts as observed experimentally.In addition, these results are consistent with the energy gap ∆∆GBottom-Side closest to 10 kcal/mol obtained by Liu et al., [32] which demonstrates that a kinetic control is essential for Z-selectivity.Figures S3 and S4 depict the optimized intermediate and transition structures for the side pathway of each catalyst.
In addition, an energy decomposition analysis (EDA) was conducted.Table 1 describes the reorganization and interaction contribution energies of the limiting stages for Ru-1 and Ru-2.In general, the reorganization contribution (∆ in Ru-1 is significantly higher in Bottom-BF than Side-DD, with respective values of 50.9 kcal/mol and 41.5 kcal/mol.Similar results were discovered for Ru-2.Thus, the side pathway requires 7-8 kcal/mol less energy than the bottom pathway to reorganize.Regarding the interaction contribution ∆ , the obtained energy values indicate that the Side-DD mechanism is more stable than in Bottom-BF.Comparing these contributions energies in both catalysts (∆  ∆ ), we find that the interaction energy is greater in the side mechanism than the reorganization energy.Consequently, it is reported that in the limiting stage of the Side-DD mechanism, the electronic effects predominate over the steric ones, whereas in Bottom-BF, the converse is true.Thus, the electronic contributions are accountable for the Z-selectivity's stability.In Tables 2 and 3, the energies of intermediates and transition states associated with the Bottom and Side pathways of the CM utilizing Ru-7 to Ru-10 catalysts are listed.The Side-DD demonstrates low energy barriers, indicating that thermal equilibrium could be readily achieved in all catalysts.In particular, the TSCR with Gibbs energies (∆G ‡ ) between 16.In Tables 2 and 3, the energies of intermediates and transition states associated with the Bottom and Side pathways of the CM utilizing Ru-7 to Ru-10 catalysts are listed.The Side-DD demonstrates low energy barriers, indicating that thermal equilibrium could be readily achieved in all catalysts.In particular, the TSCR with Gibbs energies (∆G ‡ ) between 16.4-17.4kcal/mol is the limiting stage.Optimized structures are depicted in Figure S5-S8.For the Bottom-BF, the limiting stage is the TSCA in the Ru-8 and Ru-10 with ∆G ‡ values of 28.6 kcal/mol and 33.4 kcal/mol, respectively.In Ru-7 and Ru-9, the TSCR has a ∆G ‡ of 33.3 kcal/mol.Consequently, these results are comparable to those of Ru-1 and Ru-2, indicating that Ru-7 to Ru-10 are active and stable catalysts for the production of Z-alkenes.
In addition, Table 4 displays the difference in the barrier Gibbs energy between the two mechanisms ∆∆G Bottom-Side for Ru-7 to Ru-10.The minimum value reported is 11.2 kcal/mol for the Ru-8 catalyst, while the utmost is 19.3 kcal/mol for the Ru-10 catalyst.These results demonstrate the Side preference of these catalysts, which is consistent with what has been reported for comparable complexes [32]  Figure 5 depicts the profile energy of cross metathesis utilizing Ru-1, Ru-7, and Ru-9 to compare the role of bulky ligands (Mes, Dipp, and Dipep) when a carboxylate group is used as an auxiliary ligand.Ru-7 catalysts are the most stable thermodynamically in the initial stages of the mechanism, followed by Ru-1 and Ru-9.However, when the MC-II is formed and [2+2] cycloreversion occurs, the stability order changes to Ru-9 > Ru-7 > Ru-1 (Dipep > Dipp > Mes), showing barrier energies of 18.4 kcal/mol for Ru-1, 17.4 kcal/mol for Ru-7 and 16.8 kcal/mol for Ru-9.
The energetic values corresponding to the steric and electronic contributions to the determinant step are displayed in Table 5.The binding energy (∆E = ) for Ru-1, Ru-7, and Ru-9 is −5.7 kcal/mol, −6.9 kcal/mol, and −6.6 kcal/mol, respectively.The results indicated that ∆E = is favored using bulky ligands.Additionally, reorganization (∆E reorg t ) and interaction ( ∆E int ) energies exhibit a trend like the obtained in the binding energy, demonstrating that the volume of the ligand has a significant impact on the bond energy contributions.Specifically, it can be observed that Ru-9 has a greater 14e-fragment reorganization (∆E reorg 14e ) than Ru-7 and Ru-1.Regarding ∆E int, Ru-1 and Ru-7 catalysts exhibit comparable values, −47.0 kcal/mol and -47.6 kcal/mol, respectively, whereas Ru-9 is −50.9 kcal/mol.This result indicates that Ru-9 (Dipep) has the lowest interaction energy.Furthermore, if we compare the total reorganization energy contribution to that of the interaction, we find that the former is less than the latter.Therefore, the electronic contributions are more prominent than the steric contributions at the mechanism's limiting stage.The energetic values corresponding to the steric and electronic contributions to the determinant step are displayed in Table 5.The binding energy (∆E ≠ ) for Ru-1, Ru-7, and Ru-9 is −5.7 kcal/mol, −6.9 kcal/mol, and −6.6 kcal/mol, respectively.The results indicated that ∆E ≠ is favored using bulky ligands.Additionally, reorganization (∆ ) and interaction (∆ energies exhibit a trend like the obtained in the binding energy, demonstrating that the volume of the ligand has a significant impact on the bond energy contributions.Specifically, it can be observed that Ru-9 has a greater 14e-fragment reorganization (∆ ) than Ru-7 and Ru-1.Regarding ∆ , Ru-1 and Ru-7 catalysts exhibit comparable values, −47.0 kcal/mol and -47.6 kcal/mol, respectively, whereas Ru-9 is −50.9 kcal/mol.This result indicates that Ru-9 (Dipep) has the lowest interaction energy.Furthermore, if we compare the total reorganization energy contribution to that of the interaction, we find that the former is less than the latter.Therefore, the electronic contributions are more prominent than the steric contributions at the mechanism's limiting stage.Figure 6 depicts the profile energy of cross-metathesis using Ru-2, Ru-8, and Ru-10.This study compares the role of bulky ligands (Mes, Dipp, and Dipep) when the nitrate group is employed as an auxiliary ligand.In general, the thermodynamic stability of the pathways involving the catalysts Ru-8 and Ru-10 is greater than that of Ru-2.
Nevertheless, when the determinant step was analyzed in detail, ∆ ‡ values of 16.0 kcal/mol, 17.4 kcal/mol, and 16.4 kcal/mol were determined for Ru-2, Ru-8, and Ru-10, respectively.This trend is the inverse of what was observed for Ru-1 (Mes), Ru-7 (Dipp),  Figure 6 depicts the profile energy of cross-metathesis using Ru-2, Ru-8, and Ru-10.This study compares the role of bulky ligands (Mes, Dipp, and Dipep) when the nitrate group is employed as an auxiliary ligand.In general, the thermodynamic stability of the pathways involving the catalysts Ru-8 and Ru-10 is greater than that of Ru-2.
Nevertheless, when the determinant step was analyzed in detail, ∆G ‡ values of 16.0 kcal/mol, 17.4 kcal/mol, and 16.4 kcal/mol were determined for Ru-2, Ru-8, and Ru-10, respectively.This trend is the inverse of what was observed for Ru-1 (Mes), Ru-7 (Dipp), and Ru-9 (Dipep), in which the highest barrier energy was reached for the catalyst with less bulky ligand Ru-1.
Table 6 displays the results of the EDA with Ru-2, Ru-8 and Ru-10 catalysts.The bond energy (∆E = ) is −6.2 kcal/mol, −7.5 kcal/mol, and −13.3 kcal/mol, respectively, which increased in correlation with the steric demand of the lateral NHC.Regarding the contribution of reorganization (∆E reorg t ) and interaction energy (∆E int ), the latter is greater than the former, similar to Ru-1, Ru-7, and Ru-9, emphasizing the importance of the electronic effect over the steric effect in the determinant stage.Specifically, the interaction energy for Ru-10 is 8.4 and 7.0 kcal/mol more stable than Ru-8 and Ru-2, respectively.On the other hand, the values of ∆E reorg 14e show that this contribution decreases with the increase of the bulky ligand, which is opposite to Ru-1, Ru-7, and Ru-9.Specifically, the Ru-10 catalyst exhibits a decrease of 2.4 and 1.4 kcal/mol in comparison to the Ru-2 and Ru-8 catalysts, respectively.This lower reorganization energy for the 14e-fragment is compensated by an increase of 2.8 kcal/mol in the reorganization energy of the olefin.Therefore, despite the fact that the barrier energies for Ru-9 and Ru-10 are similar and are controlled by the interaction energy, the steric effects are different due to the nature of the auxiliary ligand.
and Ru-9 (Dipep), in which the highest barrier energy was reached for the catalyst with less bulky ligand Ru-1. ) and interaction energy (∆ , the latter is greater than the former, similar to Ru-1, Ru-7, and Ru-9, emphasizing the importance of the electronic effect over the steric effect in the determinant stage.Specifically, the interaction energy for Ru-10 is 8.4 and 7.0 kcal/mol more stable than Ru-8 and Ru-2, respectively.On the other hand, the values of ∆ show that this contribution decreases with the increase of the bulky ligand, which is opposite to Ru-1, Ru-7, and Ru-9.Specifically, the Ru-10 catalyst exhibits a decrease of 2.4 and 1.4 kcal/mol in comparison to the Ru-2 and Ru-8 catalysts, respectively.This lower reorganization energy for the 14e-fragment is compensated by an increase of 2.8 kcal/mol in the reorganization energy of the olefin.Therefore, despite the fact that the barrier energies for Ru-9 and Ru-10 are similar and are controlled by the interaction energy, the steric effects are different due to the nature of the auxiliary ligand.To gain a better understanding of why Ru-10 has a lower steric effect than Ru-8, the optimized structures of the transition state were analyzed.According to Figure 7, it can be noted that in the Ru-8 catalyst, the Dipp side arm is located at 3.646 Å from the olefin, forming an angle dihedral of Ru-C1-N-C2 of 20.7° and an angle C1-Ru-C(olefin) of −6.9°.However, this distance is greater for the Ru-10 complex, being 3.720 Å and the dihedral  To gain a better understanding of why Ru-10 has a lower steric effect than Ru-8, the optimized structures of the transition state were analyzed.According to Figure 7, it can be noted that in the Ru-8 catalyst, the Dipp side arm is located at 3.646 Å from the olefin, forming an angle dihedral of Ru-C1-N-C2 of 20.7 • and an angle C1-Ru-C(olefin) of −6.9 • .However, this distance is greater for the Ru-10 complex, being 3.720 Å and the dihedral angles and angle are 35.7 • and −11 • , respectively.These results show that in Ru-10, both the side arm and the olefin adopt favorable orientations for sterically decongesting the cycloreversion reaction.Therefore, the repulsive forces between the side arm and the olefin are stronger with the Dipep compared to Dipp groups (Dipp <<<< Dipep).

Discussion
According to the preceding, a distinct behavior is observed when pivalate and nitrate are used as auxiliary ligands.For instance, the barrier energies ∆G ‡ increase when the bulky ligand is greater if the auxiliary ligand is a nitrate and decrease when the auxiliary ligand is pivalate.This result is consistent with the report by Grubbs and co-workers [22].Ru-7 (nitrate ligand) produces more Z-olefin than Ru-8 (pivalate ligand).Moreover, the interaction energies between Ru-1 and Ru-2, as well as Ru-9 and Ru-10, exhibit greater favorability in Ru-2 and Ru-10, respectively, owing to the inclusion of the nitrate ligand.
These differences can be explained by the fact that Ru-10 exhibited the lowest reorganization energy ∆E reorg 14e when compared to Ru-2 and Ru-8, suggesting that in the presence of the nitrate ligand, a rearrangement occurs more easily despite the presence of the bulky ligands.Nevertheless, for the catalysts Ru-1, Ru-7, and Ru-9 with carboxylate as an auxiliary ligand, a higher reorganization energy ∆E reorg 14e is found when there are bulky ligands, as should be expected.Thus, the structural reorganization as energy that controls the reactivity in olefin metathesis has also been observed in Mo catalysts [39], in 1st-and 2nd-generation Grubbs catalysts [30] and currently, in RCM utilizing 2nd-Grubbs catalysts, where the much larger distortion disfavors the olefin coordination [40].Table 7 presents a summary of these results.This nitrate behavior difference can be explained by a strong electronic interaction generated between an oxygen coming from the auxiliary group and one of the substituted hydrogens of the auxiliary olefin (O---H).Figure 8 illustrates the optimized structure of the limiting stage of Ru-9 and Ru-10 catalysts to demonstrate the existence of this effect.The O---H interaction lengths in the Ru-9 and Ru-10 catalyst are 2.619 Å and 2.584 Å, respectively.Therefore, the system containing nitrate ligand presents a distance smaller O---H than the carboxylate ligand, which demands low reorganization energy and demonstrates the possible existence of a "nitrate effect".Importantly, this analysis was also carried out and observed with the other catalysts Ru-1 versus Ru-2 and Ru-7 versus Ru-8, as shown in Figures S9 and S10 of the Supplementary Material (SM), respectively.
Catalysts 2023, 13, x FOR PEER REVIEW 10 This nitrate behavior difference can be explained by a strong electronic intera generated between an oxygen coming from the auxiliary group and one of the substi hydrogens of the auxiliary olefin (O---H).Figure 8 illustrates the optimized structu the limiting stage of Ru-9 and Ru-10 catalysts to demonstrate the existence of this e The O---H interaction lengths in the Ru-9 and Ru-10 catalyst are 2.619 Å and 2.584 Å spectively.Therefore, the system containing nitrate ligand presents a distance small --H than the carboxylate ligand, which demands low reorganization energy and de strates the possible existence of a "nitrate effect".Importantly, this analysis was also ried out and observed with the other catalysts Ru-1 versus Ru-2 and Ru-7 versus Ru shown in Figures S9 and S10 of the Supplementary Material (SM), respectively.Moreover, another outstanding result is the fact that the interaction energy is sli higher than the reorganization energy, suggesting that both contributions are cruci understanding the reactivity of these systems.We quantify the steric congestion in selective ruthenium catalysts using the steric maps of Cavallo and co-workers [41] (F 9).The steric maps disclose that these catalysts have high steric hindrance (%Vbur), w it reached a minimum limit of 68.3% with the Ru-2 catalyst and a maximum of 75.9%Moreover, another outstanding result is the fact that the interaction energy is slightly higher than the reorganization energy, suggesting that both contributions are crucial for understanding the reactivity of these systems.We quantify the steric congestion into Zselective ruthenium catalysts using the steric maps of Cavallo and co-workers [41] (Figure 9).The steric maps disclose that these catalysts have high steric hindrance (%V bur ), where it reached a minimum limit of 68.3% with the Ru-2 catalyst and a maximum of 75.9% with the Ru-9 catalyst.These percentages are substantially greater than those obtained with the Hoveyda-Grubbs catalyst (Figure S11 and Tables S1 and S2 of the supplementary material).The results show that if the carboxylate ligand is used as an anionic group, an increase in the steric hindrance (%V bur ) of 72, 74.8, and 75.9% for Ru-1, Ru-7, and Ru-9, respectively.In contrast, when the nitrate ligand is present, the %V bur decreases to 68.3, 70.6, and 72.0% in Ru-2, Ru-8, and Ru-10, respectively.Additionally, it is essential to analyze the free volume (%V free ).It can be associated with the free space required for the formation of the new cis-olefin (Tables S3-S6 of the supplementary material).The results show a %V free of 28, 31.7,25.2, 29.4,24.1 and 28.0 for Ru-1, Ru-2, Ru-7, Ru-8, Ru-9 and Ru-10, respectively.These results show that the larger free volume occurs in Ru-2 and Ru-8.Thus, the use of nitrate and Dipp as auxiliary and bulky ligands in a Ru catalyst can result in a higher free volume.Consequently, nitrate ligands allow a larger free volume than carboxylate, which is associated with the catalytic superiority of the nitrate complexes over the substituted tertiary butyls.To complement the discussion, an analysis of the charge distribution in the determining step, the [2+2] cycloreversion, was performed for each of the studied catalysts Ru-1, Ru-2; Ru-7, Ru-8, and Ru-9, Ru-10, with Mes, Dipp, and Dippep as the NHC ligand arms, respectively.Table 8 displays the results for each system considering the charge associated with the central atom ( ), to anionic auxiliary ligand ( ), to the NHC ligand ( ), to the adamant group ( ), the carbene bond ( ), and the olefin ( ).The results reveal that the auxiliary ligand and the carbene bond are electron acceptors, whereas the rest of the molecule has an electron donor character.Specifically, the charge distribution in  is observed to be greater with nitrate than with carboxylate ligands.The charge associated with the carbene bond  exhibits the opposite trend, with greater electron acceptor character in catalysts containing carboxylate.In addition, the ruthenium atom is a strong charge donor, whereas NHC, adamantly and the olefin are weak charge donors.Importantly, the charge distribution is independent of the steric influence that contributes to the arm NHC ligand (L = Mes, Dipp, Dipep).However, as was observed in Ru-2 (0.108), Ru-8 (0.081), and Ru-10 (0.027), the donor character in the olefin ( ) decreases with the presence of a bulky ligand.Therefore, the exceptional charac- To complement the discussion, an analysis of the charge distribution in the determining step, the [2+2] cycloreversion, was performed for each of the studied catalysts Ru-1, Ru-2; Ru-7, Ru-8, and Ru-9, Ru-10, with Mes, Dipp, and Dippep as the NHC ligand arms, respectively.Table 8 displays the results for each system considering the charge associated with the central atom (Q Ru ), to anionic auxiliary ligand (Q Anionic ), to the NHC ligand (Q NHC ), to the adamant group (Q Adm ), the carbene bond (Q carbene ), and the olefin (Q ole f in ).The results reveal that the auxiliary ligand and the carbene bond are electron acceptors, whereas the rest of the molecule has an electron donor character.Specifically, the charge distribution in Q Anionic is observed to be greater with nitrate than with carboxylate ligands.The charge associated with the carbene bond Q carbene exhibits the opposite trend, with greater electron acceptor character in catalysts containing carboxylate.In addition, the ruthenium atom is a strong charge donor, whereas NHC, adamantly and the olefin are weak charge donors.Importantly, the charge distribution is independent of the steric influence that contributes to the arm NHC ligand (L = Mes, Dipp, Dipep).However, as was observed in Ru-2 (0.108), Ru-8 (0.081), and Ru-10 (0.027), the donor character in the olefin (Q ole f in ) decreases with the presence of a bulky ligand.Therefore, the exceptional characteristics of catalysts with nitrate ligand, such as the highest electron acceptor behavior and the lowest reorganization energy, are a result of a combination of electronic and steric effects.

Models
The studied catalysts were divided into two groups, as shown in Figure 10.The first group consists of the active species (14-Int), so-called Keitz-Grubbs-substituted complexes.These complexes are characterized by the presence of an NHC with a Mes substituent and bidentate anionic ligands, such as tert-butyl (Ru-1) and nitrate (Ru-2).The second group includes the active species of complexes reported by Grubbs et al.

Models
The studied catalysts were divided into two groups, as shown in Figure 10.The first group consists of the active species (14-Int), so-called Keitz-Grubbs-substituted complexes.These complexes are characterized by the presence of an NHC with a Mes substituent and bidentate anionic ligands, such as tert-butyl (Ru-1) and nitrate (Ru-2).The second group includes the active species of complexes reported by Grubbs et al.The active catalytic species (14-Int) of the catalysts were chosen based on the findings reported by Núñez-Zarur and co-workers [42].According to his research, the formation of this intermediate in the second-generation Hoveyda-Grubbs precatalysts occurs by a dissociative mechanism, primarily driven by the kinetic preference for the elimination of the alkoxy group.Thus, with Hoveyda-Grubbs catalysts, the initiation step involves (a) the loss of the chelating alkoxy group and (b) a cross-metathesis process in which the initial substituent of the active species is exchanged for the substituent of the reactive alkene.The preceding was depicted in Figure S12 of the Supplementary material.

Computational Details
Density Functional Theory (DFT) was used as it is implemented in Gaussian-16 software.All molecular structures were optimized in the gas phase using the meta-GGA-type functional M06-L.This function has been previously studied for the Grubbs catalysts because it includes attractive non-covalent interactions that are determinants for an accurate The active catalytic species (14-Int) of the catalysts were chosen based on the findings reported by Núñez-Zarur and co-workers [42].According to his research, the formation of this intermediate in the second-generation Hoveyda-Grubbs precatalysts occurs by a dissociative mechanism, primarily driven by the kinetic preference for the elimination of the alkoxy group.Thus, with Hoveyda-Grubbs catalysts, the initiation step involves (a) the loss of the chelating alkoxy group and (b) a cross-metathesis process in which the initial substituent of the active species is exchanged for the substituent of the reactive alkene.The preceding was depicted in Figure S12 of the Supplementary material.

Computational Details
Density Functional Theory (DFT) was used as it is implemented in Gaussian-16 software.All molecular structures were optimized in the gas phase using the meta-GGAtype functional M06-L.This function has been previously studied for the Grubbs catalysts because it includes attractive non-covalent interactions that are determinants for an accurate energy description [43,44].The electronic configuration of the ruthenium atom combines the MWB28 [45,46] pseudopotential, which replaces the inner electrons with a non-local effective potential, and the (8s7p6d)/[6s5p3d] basis set to describe the valence electrons.The 6-31 + G(d, p) basis set was utilized to represent carbon, hydrogen, oxygen, and nitrogen atoms.To ensure the nature of each structure as an intermediate or transition state on the potential energy surfaces (PES), a vibrational analysis of all stationary points was performed by calculating the Hessian matrix in internal coordinates (second derivative of the energy with respect to position).These frequencies were used to determine the thermal contributions to enthalpy (H) and Gibbs free energy (G) utilizing the harmonic oscillator and rigid rotor model implemented in Gaussian-16 and assuming a temperature of 298 K at standard pressure (1 atm).The effects of the tetrahydrofuran (THF) solvent were included with the density-based solvation model (SMD) [45,46] as energetic corrections from single-point calculations.In summary, the reported Gibbs free energy for each species in the PES was obtained by G = G opt + G solv .

Theoretical Background 4.3.1. Energy Decomposition Analysis (EDA)
The energy decomposition analysis allows us to understand the changes of a chemical process through the equation proposed by Houk [47] to quantitatively assess binding energy (E = ).This model has been successfully used for olefin metathesis reactions in Molybdenum and Ruthenium complexes [39,48].The binding energy of the determinant step into the mechanism may be comprehended by considering two distinct contributions.These can be attributed to the rearrangement ∆E reorg total and the interaction energy ∆E int between the metal fragment and the olefin.The former is related to the steric effect, while the latter is associated with the electronic influence.The energy partitioning from the energy barrier (E = ) is defined as ∆E reorg total analyzes the changes in the geometry of the active species ∆E reorg 14 and of the olefin ∆E reorg ole f in allowing the degree of structural distortion between the fundamental state and the distorted state to be quantified so that: ∆E reorg total = ∆E reorg 14 + ∆E reorg ole f in = E 14 e−atcomplex − E 14 e−reactive + E ole f in at complex − E ole f in The interaction energy is then obtained using the results of ∆E = and ∆E reorg total :

Topographic Maps
The application of topographic steric maps has been extremely useful in the design of transition metal complexes [25,38,49,50].The computational tool developed by Carvallo et al. [41] provides a picture of the interaction surface between the catalyst and the substrates that are formed by the ligands in the complex.The steric hindrance of the above catalysts in the limiting stage has been represented using topographic maps through a comparative analysis of the buried volume (%V bur ); the preceding was interpreted [51].The orientation of the map was adjusted to avoid steric interference from NHC side arms.The Ru atom was defined to be located at the origin, while the Ru-alkylidene resides on the z+ axis, and the bulk-type ligand lies in the xy plane.(Figure S13) The radius of the sphere surrounding the center of the metal was set to 5.0 Å [52,53].(Figure S11) For the other atoms, the Bondi radii were scaled by 1.17, and a mesh of 0.1 Å was selected to scan the sphere for buried voxels.H atoms were included in the calculations.

Conclusions
This study investigates the impact of bulky ligands (L = Mes, Dipp, and Dippep) on novel sterically demanding N-heterocyclic cyclometalated ruthenium complexes.Nitrate and pivalate are employed as auxiliary ligands to understand the underlying factors contributing to the Z-selectivity in cross-metathesis reactions.Our results reveal the Side preference in the reaction mechanism indicated by the difference in Gibbs energy barriers ∆∆G Bottom-Side ranging from 11.2-19.3kcal/mol.The [2+2] cycloreversion was identified as the rate-determinant step, showing comparable energy values in all the catalysts.However, it has been shown that the activation energies ∆G ‡ exhibit an increase as the size of the bulky ligand increases in catalysts, including nitrate as an auxiliary ligand.Conversely, the energies display a decrease when pivalate is employed as the auxiliary ligand.This difference can be attributed to the relatively low reorganization energy required for a rearrangement in the coordination environment.This is further evidenced by the short O---H distance between the oxygen atom of the auxiliary ligand and the hydrogen atom of the olefin.Moreover, the results also explain that the rate-determinant step is not just influenced by steric effects but also by electronic factors coming from the electron acceptor character of the auxiliary ligand and the donor characteristics of the ruthenium atom.

Supplementary Materials:
The following supplementary material information can be downloaded at: https://www.mdpi.com/article/10.3390/catal13091305/s1,Figures S1 and S2: Gibbs profile energies of catalysts Ru-1 and Ru-2.Figures S3-S8: Geometry of the optimized structures associated to Ru-1, Ru-2, Ru-7, Ru-8, Ru-9 and Ru-10.Figure S9: Optimized geometry of the TSCR of the Ru-1 and Ru-2 catalysts.Figure S10: Optimized geometry of the TSCR of the Ru-7 and Ru-8 catalysts.Figure S11: Radius of the sphere built around Ru-1. Figure S12: Dissociation mechanism in Hoveyda-Grubbs catalysts.Figure S13: Visualization of the orientation used in steric maps with the Ru-1 catalyst.Table S1: Values of %V free and %Buried for Ru-1; Table S2: Values of %V free and %Buried for Ru-2; Table S3: Values of %V free and %Buried for Ru-7; Table S4: Values of %V free and %Buried for Ru-8; Table S5: Values of %V free and %Buried for Ru-9; Table S6: Values of %V free and %Buried for Ru-10; Cartesian coordinates of the structure optimized by DFT.
4-17.4 kcal/mol is the limiting stage.Optimized structures are depicted in Figure S5-S8.For the Bottom-BF, the limiting stage is the TSCA in the Ru-8 and Ru-10 with ∆G ‡ values of 28.6 kcal/mol and 33.4 kcal/mol, respectively.In Ru-7 and Ru-9, the TSCR has a ∆G ‡ of 33.3 kcal/mol.Consequently, these results are comparable to those of Ru-1 and Ru-

2. 2 .
4-17.4 kcal/mol is the limiting stage.Optimized structures are depicted in Figure S5-S8.For the Bottom-BF, the limiting stage is the TSCA in the Ru-8 and Ru-10 with ∆G ‡ values of 28.6 kcal/mol and 33.4 kcal/mol, respectively.In Ru-7 and Ru-9, the TSCR has a ∆G ‡ of 33.3 kcal/mol.Consequently, these results are comparable to those of Ru-1 and Ru-Bottom and Side DFT Mechanism Using Novel Ruthenium Catalysts, Ru-7 to Ru-10

Figure 5 .
Figure 5. Profile of Gibbs energy for the CM reaction of propene with Ru-1, Ru-7, and Ru-9 via Side-DD.Values are in kcal/mol.

Figure 5 .
Figure 5. Profile of Gibbs energy for the CM reaction of propene with Ru-1, Ru-7, and Ru-9 via Side-DD.Values are in kcal/mol.

Figure 6 .
Figure 6.Profile of Gibbs energy for the CM reaction of propene with Ru-2, Ru-8, and Ru-10 via Side-DD.Values are in kcal/mol.Table6displays the results of the EDA with Ru-2, Ru-8 and Ru-10 catalysts.The bond energy (∆E ≠ ) is −6.2 kcal/mol, −7.5 kcal/mol, and −13.3 kcal/mol, respectively, which increased in correlation with the steric demand of the lateral NHC.Regarding the contribution of reorganization (∆ ) and interaction energy (∆ , the latter is greater than

Figure 6 .
Figure 6.Profile of Gibbs energy for the CM reaction of propene with Ru-2, Ru-8, and Ru-10 via Side-DD.Values are in kcal/mol.

Catalysts 2023 ,
13, x FOR PEER REVIEW 9 o angles and angle are 35.7°and −11°, respectively.These results show that in Ru-10, b the side arm and the olefin adopt favorable orientations for sterically decongesting cycloreversion reaction.Therefore, the repulsive forces between the side arm and the ole are stronger with the Dipep compared to Dipp groups (Dipp <<<< Dipep).

Figure 7 .
Figure 7. Optimized geometry of the TSCR for Ru-8 and Ru-10 in the mechanism Side-DD.Distan are given in Å.

Figure 7 .
Figure 7. Optimized geometry of the TSCR for Ru-8 and Ru-10 in the mechanism Side-DD.Distances are given in Å.

Figure 8 .
Figure 8. Optimized geometry of the TSCR for Ru-9 and Ru-10 in the mechanism Side-DD.Dist are expressed in Å.

Figure 8 .
Figure 8. Optimized geometry of the TSCR for Ru-9 and Ru-10 in the mechanism Side-DD.Distances are expressed in Å.

Figure 9 .
Figure 9. Topographic steric maps of ruthenium catalysts in the determinant stage of the reaction mechanism.

Figure 9 .
Figure 9. Topographic steric maps of ruthenium catalysts in the determinant stage of the reaction mechanism.
in 2020.These systems incorporate the bulky lateral arms of Dipp and Dipep.The catalysts are denoted as Ru-7 (L= Dipp, t-Bu-C); Ru-8 (L = Dipp, O-N+); Ru-9 (L = Dipep, t-Bu-C); and Ru-10 (L = Dipep, O-N+).It is important to highlight that Ru-10 is a proposed catalyst because computational and experimental efforts have shown that the inclusion of the nitrate group in conventional catalysts has significantly improved Z-selectivity [17].Catalysts 2023, 13, x FOR PEER REVIEW 12 of 16 in 2020.These systems incorporate the bulky lateral arms of Dipp and Dipep.The catalysts are denoted as Ru-7 (L= Dipp, t-Bu-C); Ru-8 (L = Dipp, O-N+); Ru-9 (L = Dipep, t-Bu-C); and Ru-10 (L = Dipep, O-N+).It is important to highlight that Ru-10 is a proposed catalyst because computational and experimental efforts have shown that the inclusion of the nitrate group in conventional catalysts has significantly improved Z-selectivity [17].

Figure 10 .
Figure 10.Active species used in this work.

Figure 10 .
Figure 10.Active species used in this work.

Table 1 .
∆E = , ∆E reorg t and ∆E int for the determinant step of the CM reaction of propene using the Keitz-Grubbs-substituted complexes, Ru-1 and Ru-2.Values are in kcal/mol.stablethan the Bottom-BF.The limiting stage or Ru-1 is the TSCR with barrier energies of 20.3 kcal/mol and 32.4 kcal/mol, respectively.Consequently, the barrier Gibbs energy ∆∆GBottom-Side difference for Ru-1 and Ru-2 are 12.1 kcal/mol and 18.4 kcal/mol, respectively.

Table 1 .
∆E ≠ , ∆ and ∆ for the determinant step of the CM reaction of propene using the Keitz-Grubbs-substituted complexes, Ru-1 and Ru-2.Values are in kcal/mol.

Table 1 .
∆E ≠ , ∆ and ∆ for the determinant step of the CM reaction of propene using the Keitz-Grubbs-substituted complexes, Ru-1 and Ru-2.Values are in kcal/mol.

Table 2 .
Gibbs energy values (kcal/mol) of the intermediates and transition states connected to the Bottom-BF mechanism of Ru-7 to Ru-10 catalysts.

Table 3 .
Gibbs energy values (kcal/mol) of the intermediates and transition states connected to the Side-DD mechanism of Ru-7 to Ru-10 catalysts.

Table 4 .
Difference in the barrier Gibbs energy ∆∆G Bottom-Side between Ru-7 and Ru-10 catalysts.Values are in kcal/mol.
2.3.Side Mechanism and Analysis of the Electronic and Steric Role in Recent Z-Selective Catalysts (Ru-1, Ru-7, and Ru-9) Containing a Carboxylate as Anionic Group

Table 5 .
∆E ≠ , ∆ , and ∆ for the determinant step of the CM reaction of propene using Ru-1, Ru-7, and Ru-9.Values are in kcal/mol.

Table 5 .
∆E = , ∆E reorg t , and ∆E int for the determinant step of the CM reaction of propene using Ru-1, Ru-7, and Ru-9.Values are in kcal/mol.

Table 6 .
∆E ≠ , ∆ , and ∆ for the determinant step of the CM reaction of propene using Ru-2, Ru-8, and Ru-10.Values are in kcal/mol.

Table 6 .
∆E = , ∆E reorg t , and ∆E int for the determinant step of the CM reaction of propene using Ru-2, Ru-8, and Ru-10.Values are in kcal/mol.

Table 7 .
Summary of the most significant results.Values are in kcal/mol.

Table 8 .
Atomic charges (Q) associated with the stabilizing ligands and the metal fragment.