URVA and Local Mode Analysis of an Iridium Pincer Complex Efﬁciently Catalyzing the Hydrogenation of Carbon Dioxide

: The catalytic effects of iridium pincer complexes for the hydrogenation of carbon dioxide were investigated with the Uniﬁed Reaction Valley Approach (URVA), exploring the reaction mechanism along the reaction path traced out by the reacting species on the potential energy surface. Further details were obtained with the Local Mode Analysis performed at all stationary points, complemented by the Natural Bond Orbital and Bader’s Quantum Atoms in Molecules analyses. Each of the ﬁve reaction paths forming the catalytic cycle were calculated at the DFT level complemented with DLPNO-CCSD(T) single point calculations at the stationary points. For comparison, the non-catalytic reaction was also investigated. URVA curvature proﬁles identiﬁed all important chemical events taking place in the non-catalyzed reaction and in the ﬁve reactions forming the catalytic cycle, and their contribution to the activation energy was disclosed. The non-catalytic reaction has a large unfavorable activation energy of 76.3 kcal/mol, predominately caused by HH bond cleave in the H 2 reactant. As shown by our study, the main function of the iridium pincer catalyst is to split up the one–step non-catalytic reaction into an energy efﬁcient multistep cycle, where HH bond cleavage is replaced by the cleavage of a weaker IrH bond with a small contribution to the activation energy. The dissociation of the ﬁnal product from the catalyst requires the cleavage of an IrO bond, which is also weak, and contributes only to a minor extent to the activation energy. This, in summary, leads to the substantial lowering of the overall activation barrier by about 50 kcal/mol for the catalyzed reaction. We hope that this study inspires the community to add URVA to their repertoire for the investigation of catalysis reactions.

The Nozaki group identified six-coordinate Ir(III)-PNP trihydride (PNP = 2,6-bis(diisopropylphosphinomethyl)pyridine) as a highly active catalyst with high turn over frequencies (TOF) and turnover number (TON), reflecting the stability of the active site for the CO 2 hydrogenation [41,43,44].The original mechanism of the CO 2 hydrogenation to formic acid was reported to include three major steps, (i) CO 2 addition to the catalyst, (ii) formic acid dissociation, and (iii) H 2 addition in order to regenerate the catalyst [43].In the following years, the mechanism was extended including, e.g., basic conditions leading to a more complex catalytic cycle and a larger number of intermediates confirmed by DFT calculations [44][45][46].Because the insertion of CO 2 into the metal hydride is the crucial step of the catalysis [47,48], possible CO 2 insertion into a five-coordinate Iridium(III) dihydride complex as alternative was experimentally investigated [49], showing the formation of k 2 -bound formate monohydride products by selective electro-catalytic reduction of CO 2 .The mechanism of the CO 2 insertion into five-coordinate Iridium(III) dihydrides and four-coordinate Iridium(I) hydrides was also analyzed theoretically [50], exploring the different reactivity of these complexes with different iridium oxidation states.The catalytic reduction of CO 2 to methane with silanes using Brookhart's cationic Ir(III) pincer complex was investigated theoretically [51], focusing on the transfer of silane hydrogens to CO 2 and the formation of silylformate, bis(silyl)acetal, methoxysilane and methane.The importance of a hydrogen bond donor in the secondary coordination sphere of the Ir(III)-PNP catalyst was analyzed theoretically [52], finding that the isolation of the formate product from the reaction is straight forward and therefore makes this reaction favorable.CO 2 hydrogenation was also investigated using a series of modifications of the Iridium complexes bearing sophisticated NˆN-bidentate ligands [53] showing importance of hydroxy groups as proton-responsive substituents in CO 2 hydrogenation reactions.Other modifications of the iridium complexes included amide-based ligands, showing that the electron-donating effects of an anionic nitrogen atom and the presence of an OH group near the metal center improves the catalytic activity for CO 2 hydrogenation [54].Although these studies provide valuable insights, a clear understanding of the mechanistic details of each reaction step of the CO 2 hydrogenation catalytic cycle and how the catalyst actually lowers the activation barrier has been missing so far.
Therefore, we investigated in this study all major chemical events taking place during the CO 2 hydrogenation to formic acid, for the original catalytic cycle proposed by Nozaki [43], using as computational tools the Unified Reaction Valley Approach (URVA) and the Local Modes Analysis (LMA) developed in our group, which are both based on vibrational spectroscopy.URVA requires the accurate and precise determination of the reaction path starting at the transition state (TS) and moving far down into both the entrance and exit channels, a time consuming enterprise.Therefore, we simplified the original Ir(III)-PNP complex in our study by replacing the isopropyl groups with hydrogen atoms.The catalytic cycle involves the five reaction steps R1-R5 shown in Figure 1.For comparison we also included the non-catalyzed reaction R0.Reaction movies illustrating the movements of all atoms of the complexes along the reaction path together with the Cartesian coordinates of the stationary points for each of the five reaction paths R1-R5 of the catalytic cycle and the non-catalyzed reaction R0, both investigated in the gas phase, are presented in the Supplementary Materials.

Materials and Methods
The theoretical background of URVA is provided in a comprehensive review article [55], therefore only some highlights are summarized.URVA explores the reaction complex (RC, i.e., the union of reacting molecules) along the path it traces out on the potential energy surface (PES) starting from the TS down into the entrance and exit channels, described by a large amplitude motion [56][57][58][59] as well as in the close vicinity of the reaction path, the so-called reaction valley being spanned by the vibrations perpendicular to the reaction path.[60,61].The major focus of URVA is on the curvature of the reaction path.Since the reaction path is a curved line in N vib dimensional space with N vib = 3N L internal coordinates (L = 6 for a non-linear N-atomic RC and 5 for an N-atomic linear RC), its direction and curvature can be derived with the Frenet-Serret formalism [62].The reaction path direction at a path point s is given by the unit vector h(s): where the derivative of the mass-weighted reaction coordinate x(s) with regard to s is the normalized mass-weighted gradient vector g(s) ⌘ g( x(s)) = M 1/2 g(s) and M is a diagonal matrix of atomic masses.c(s) is the normalization constant being equal to the length of the gradient vector k g(s)k.The curvature vector k(s) is given by [63,64] where fx (s) is the mass-weighted Hessian matrix in Cartesian coordinates x.Any electronic structure change, e.g., bond breaking/forming processes, but also rehybridization, charge polarization and transfer, etc. leads to distinct curvature peaks, which are directly reflected in the scalar reaction path curvature k(s) = ||k(s)|| in contrast to regions of minimal electronic change, which are reflected by curvature minima as sketched in Figure 2 showing the scalar curvature as a function of the reaction path parameter s.We call the region from one curvature minimum to the next embedding a curvature peak a reaction phase.Each chemical reaction has a unique pattern of curvature maxima and minima and as such a different number of reaction phases leading to a specific reaction profile, the so-called fingerprint of the reaction [55].Deeper insights into the nature of a chemical event reflected by a curvature peak are obtained from the decomposition of the scalar reaction curvature into individual components [65] such bond lengths, bond angles, and dihedral angles or other components such as puckering coordinates [66] denoted by the colored lines in Figure 2. Interesting to note is that even for complex catalytic reactions a curvature peak is most likely composed of a few components, which facilitates the analysis of the mechanistic details considerably.The sign of a component indicates if the parameter in question supports the chemical event (positive sign) or if it resists the chemical change (negative sign) [65].For a detailed mathematical derivation, its connection to the work of Miller, Handy and Adams [67] and Kato and Morokuma [68] on the reaction path Hamiltonian, and recent advances of URVA, interested readers are referred to Refs [55,69].Applications of URVA can be found, e.g., in Refs [69][70][71][72][73][74][75][76][77][78][79][80][81][82][83][84].Similar as URVA, LMA is based on vibrational spectroscopy and it refines the use of normal vibrational force constants and frequencies derived from normal vibrational modes to characterize chemical bonds and/or weak chemical interactions.Normal vibrational modes of a polyatomic molecule are generally delocalized, as already stated by Wilson in 1941 via his proof that the associated normal mode coordinates are a linear combination of internal coordinates q or Cartesian coordinates x [85].Therefore, normal mode stretching frequencies and associated stretching force constants are of limited use as individual bond strength descriptors.Konkoli and Cremer solved this problem via the transformation of normal vibrational modes into their local mode counterparts [86][87][88][89][90] (for details the reader may refer to two comprehensive review articles [91,92]).A local vibrational mode a n is defined as implying the important result that all what is needed for the local mode analysis are the diagonal normal mode force constant matrix K in normal mode coordinates and the normal mode vectors d n in internal coordinates, which can be obtained from a standard vibrational frequency calculation via the Wilson GF formalism [93][94][95], a routine part of most modern quantum chemistry packages.
Once the local mode vector a n is known, one can define molecular properties corresponding to this local motion, such as local mode force constant, local mass and local frequency [86][87][88][89][90].
The corresponding local force constant k a n of local mode n (superscript a denotes an adiabatically relaxed, i.e., local mode) [96] can be expressed as As has been shown by Zou et al., local mode stretching force constants reflect the intrinsic strength of a chemical bond and/or weak chemical interaction [97].
The associated local mass is defined as where G nn is a diagonal element of Wilson's inverse kinetic energy matrix G matrix and the Wilson B matrix provides the important relationship between internal and Cartesian coordinates via the first derivatives of the internal coordinates q n (n = 1, 2, 3 . . .N vib ) with regard to the Cartesian coordinates x i (i = 1, 2, 3 . . .3N) [93], From the local mode force constant k a and local mode mass m a the local mode frequency w n can be calculated with c being the speed of light.In addition, the local mode infrared intensity has been defined which can be related to bond dipole moments [98].
It is convenient to transform local mode stretching force constants k a into relative bond strength order (BSO) via a power relationship BSO = A ⇤ (k a ) B being being based on the generalized Badger rule [83,139].The dimensionless parameters A and B are obtained from two reference molecules with known BSO and k a values for the bond in question and the requirement that for a zero k a the BSO is zero too.We used different reference molecules for different bond types investigated in our study, for the Ir-ligand bonds we used Mayer's bond orders [140][141][142].Details of the reference bonds and reference molecules are summarized in Table 1.The resulting A and B parameters for each bond are documented in the corresponding plots of the results and discussion part. 2Based on the CH 3 OH and CH 2 O references. 3 Based on the on the HF and F 2 H references.
Bader's quantum theory of atoms-in-molecules (QTAIM) characterizes of chemical bonds based on the topological features of the total electron density r(r) [143][144][145][146].According to QTAIM a chemical bond between two atoms A and B is characterized by a bond path and a bond critical point r.The covalent character of the bond AB can be determined via the Cremer-Kraka criterion of covalent bonding [147][148][149].According to this criterion, a chemical bond has covalent character when the energy density H r , i.e., the energy density H(r) taken at the at the bond critical point r is negative, while a positive value of H r indicates an electrostatic character.H(r) = G(r) + V(r); the positive kinetic energy density G(r) describes electron density depletion and the negative potential energy density V(r) describes electron density accumulation.Therefore, a chemical bond has covalent character when the potential energy density dominates.
We used in this study as a reaction path the intrinsic coordinate (IRC) path of Fukui [150,151] combined with the improved reaction path following procedure of Hratchian and Kraka, enabling the tracking of a chemical reaction far out into the entrance and exit channel [151].We employed an IRC step size of 0.03 amu 1/2 Bohr and an ultrafine grid for numerical DFT integrations.IRC calculations and the investigation of all stationary points were performed with the B3LYP density functional [152][153][154] utilizing Pople's 6-31G(d,p) basis set [155][156][157] and the Stuttgart-Dresden (SDD) effective core potential for iridium [158].The energetics of the investigated chemical reactions was additionally recalculated via single energy point calculations with the coupled cluster DLPNO-CCSD(T) method [159] utilizing the def2-TZVP basis set [160] and the def2-ECP pseudo potential for iridium [158].The DLPNO-CCSD(T) calculations were based on the DFT geometries and applying thermochemical corrections from the DFT frequency calculations.DFT calculations were performed with Gaussian program package [161], and DLPNO-CCSD(T) calculations with ORCA [162].URVA was carried out using the program pURVA [163], and LMA analysis was performed using the LModeA program [164].The Natural Bond Orbital (NBO) analysis [165][166][167][168] was applied for calculating NBO atomic charges along the reaction path.The QTAIM analysis was performed with the AIMALL program package [169].

Energetics
Table 2 presents the activation and reaction energies E a , E R and activation and reaction enthalpies H a , H R of reactions R0-R5 calculated at the B3LYP/6-31G(d,p)/SDD(Ir) and DLPNO-CCSD(T)/def2-TZVP/def2-ECP(Ir) levels of theory.The energy and enthalpy values in Table 2 are taken for each reaction relative to the energy of first point on the IRC (i.e., first point in the entrance channel corresponding to a van der Waals complex).The following discussion is based on the DLPNO-CCSD(T) results.According to Table 2, the activation energy E a of the non-catalyzed reaction with a value of 76.3 kcal/mol is considerably larger then that for all reactions of the catalytic cycle, ranging between 1.7 kcal/mol for reaction R4, and 23.9 kcal/mol for reaction R2. Figure 3 displays the overall energetics of the catalytic cycle being composed of reactions R1-R5.All energies are relative to the energy of the van der Waals complex formed between the catalyst and CO 2 1 in the entrance channel.All energies in Figure 3 are based on the DLPNO-CCSD(T) calculations performed at the end points of the DFT IRC calculations and for simplicity they have been calculated by adding the activation and reaction energy of the next reaction step of the cycle, to the product energy of the previous reaction step.
Table 2.The activation and reaction energy E a , E R and the activation and reaction enthalpy H a , H R of the reactions R0-R5 investigated in this study 1 .The thermochemistry corrections of the DLPNO-CCSD(T) method is based on the DFT optimal geometry and frequency calculations.The B3LYP/6-31G(d,p)/SDD(Ir) and DLPNO-CCSD(T)/def2-TZVP/def2-ECP(Ir) levels of theory.According to Table 2 and Figure 3, reactions R1-CO 2 addition (E a = 17.6 kcal/mol), R2-hydrogen transfer (E a = 23.9kcal/mol), and R5-catalyst regeneration (E a = 23.5 kcal/mol), have the largest activation energies, i.e., they are the decisive steps of the overall catalytic cycle.The values of the activation enthalpies H a for reactions R1 (18.9 kcal/mol), R2 (22.7 kcal/mol), and R5 (21.3 kcal/mol) show a similar trend as seen from the activation energies E a of these reactions, which indicates small thermal energy contributions to the overall reaction energies.Although information about the energetics is useful, it does not disclose the actual mechanism taking place in each reaction step forming the catalytic cycle.This is explored in the next section.

Reaction Mechanism
Reaction R0 Figure 4a presents the energy profile of the non-catalyzed reaction R0 along the reaction path.The decomposition of the reaction curvature into selected components is presented in Figure 4b, selected geometrical parameters in Figure 4c and atomic NBO charges along the reaction path are presented in Figure 4d.According to Figure 4b, the reaction starts with initial formation of the O a H a bond in reaction phase 2 (green line), with a resisting contribution of the O a C a O b bond angle (purple line).Bending of CO 2 molecule is also reflected in Figure 4c.The O a C a O b bond angle starts to decrease in phase 2. As shown in Figure 4a, both events the initial O a H a bond formation and the O a C a O b bending, contribute 6.8 kcal/mol to the activation energy in this reaction phase.As revealed in Figure 4b, breaking of the H a H b bond takes place in phase 4, which is confirmed by a strong supporting contribution of the H a H b component (blue line).HH bond cleavage is also reflected by the change of H a H b bond length (see Figure 4c), which starts to elongate in phase 4, making H a positively charged (see Figure 4d).H a H b bond breaking is supported by the start of C a H b bond formation (red line in Figure 4b), which is also confirmed by changes of C a and H b atomic charges (see Figure 4d).As shown in Figure 4a, the chemical events happening in phase 4, in particular the energy demanding breakage of the hydrogen bond, contributing 40.8 kcal/mol to the activation energy is the major cause of the high barrier.Finalization of O a H a and C a H b bonds takes place in phases 5 and 6 after TS (see Figure 4b),which is also reflected by O a H a and C a H b bond length changes (Figure 4c) and the NBO charge changes (Figure 4d).Obviously, the catalyzed reactions circumvents direct HH bond breakage, which is explored next.
Reaction R1 Properties of reaction R1 of the catalytic cycle along the reaction path are presented in Figure 5a-d   Reaction R2 Reaction R2 of the catalytic cycle is characterized by the H atom transfer from the C b carbon atom of the catalyst to the O a atom of CO 2 .This is a necessary step for formation of the formic acid product, which goes along with elongation of the intermediate IrO b bond.Figure 6a-d show the corresponding URVA analysis.According to Figure 6b, C b H a bond cleavage (orange line) takes place in phase 5, strongly supporting, and in phase 6, strongly resisting.The C b H a distance starts to increase in phase 5, however it takes until the end of phase 10 to reach its final distance of 3.1 Å (see Figure 6c).
The formation of the new O a H a bond starts in phase 4 (red line), as depicted in Figure 6b and is finalized in phase 6, also reflected by the corresponding change of the O a H a distance shown in Figure 6c.What is intersting to note is that the changes of O a H a and C b H a distances complement each other.A secondary contribution to the mechanism is the change of C a O a from double to single bond character (purple line in Figure 6b) in phases 5 and 6.Overall, changes in the atomic charges are only marginal (see Figure 6d).According to Figure 6a, the largest contributions to the activation energy are i) O a H a bond formation and elongation of the IrO b bond in phase 4 (3.8 kcal/mol) and ii) hydrogen atom transfer accompanied by the change of the CO bond from double to single bond in phase 5 (9.8 kcal/mol).Reaction R3 In reaction R3 of the catalytic cycle, formic acid is released, accompanied by a cleavage of the intermediate IrO b bond (already weakened in R2), reducing the metal coordination number from 6 to 5, which leads to a reorganization the Ir coordination sphere, in particular a change in the H c and H d ligand positions.The URVA results of R3 are summarized in Figure 7a-d.The curvature profile 7b shows that IrO b bond cleavage (blue line) starts in phase 2, with a supporting contribution.This is also reflected by a steep increase in the IrO b distance starting in that phase, (see Figure 7c) and an increase of the Ir charge (see Figure 7d).The next stages of IrO b bond cleavage occur in phases 3 and 4 (see IrO b ) first resisting and then supporting (phase 5, after the TS) with the final IrO b distance of 4.1 Å reached at the end of the reaction (see Figure 7c).
According to Figure 7b, there is also a contribution from the NH a hydrogen bond cleavage (pink line) in phase 3 (the NH a hydrogen bond is not shown in Figure 1 for simplicity), which starts in phase 2 and is finalized in phase 6.The cleavage of the NH a bond is also reflected in Figure 7c showing a steady increase of the NH bond distance from  Reaction R4 Figure 8a-d illustrate the change of the reaction properties along the reaction path for reaction R4 of the catalytic cycle, which involves the addition of diatomic hydrogen to the coordination sphere of the metal via h 2 -H 2 coordination; the first part of restoring the catalyst, which implies an increase of the Ir coordination number from 5 to 6.According to Figure 8b the IrH a (green line) and IrH b (blue line) bond components contribute to the reaction curvature in phases 2, 3, and 5.As depicted in Figure 8c IrH a and IrH b distances start to decrease in phase 2. Although the starting H a distance is 4.2 Å, considerably larger than the H b distance of 3.45 Å, both bonds are finalized simultaneously at the end of phase 5.There is also a minor contribution of the H a H b component (red line) to the reaction curvature in phase 5 (see Figure 8b), reflecting changes of the H a H b bond in this phase, where the h 2 -H 2 complexation to the metal is finalized.The change of the metal coordination sphere in this reaction phase leads to an increase of negative charge on Ir atomic atom, as shown in Figure 8d.The activation energy of reaction R4 is relatively small, (DFT result: E a = 1.2 kcal/mol, Figure 8a), mainly resulting from electronic structure reorganization of the IrH a and IrH b bonds taking place in phase 3. Reaction R5 In reaction R5 the initial catalyst is recovered, which implies H a H b bond cleavage of the h 2 -H 2 coordinated H 2 ligand, IrH a bond cleavage, migration of H a to C b and C b H a bond formation.The corresponding URVA analysis is illustrated in Figure 9a-d.Figure 9b shows that IrH a bond cleavage (orange line) starts in phase 4 and proceeds through the following phases until phase 10.As reflected in Figure 9b, the initial IrH a distance of 1.8 Å starts to increase in phase 4 until its final value of 3.5 Å is reached at the end of the reaction.H a H b bond cleavage (red line) starts at the border between phases 3 and 4, dominates phase 6 shortly before the TS with a large supporting contribution and proceeds parallel to IrH a bond cleavage until phase 10.It is interesting to note that elongation of the H a H b bond starts in phase 6 (see Figure 9c) which parallels charge separation, H a migrating to the C b carbon atom to form a CH bond becomes more positively charged and H b forming a stronger IrH bond, more negatively charged (see Figure 9d).According to Figure 9b, C b H a bond formation (green line) starts in phase 3 and is finalized in phase 7, after the TS and is revealed by a large curvature contribution.In phase 7 the C b H a bond length reached its final value, as shown in Figure 9c.According to Figure 9a, the largest energy contributions to the activation energy occur in reaction phases 4 (5.5 kcal/mol) and 5 (7.1 kcal/mol), corresponding to IrH a bond cleavage and strengthening of the IrH b bond (see Figure 9b).In addition, it is important to note that H a H f b bond cleavage contributes only 4.2 kcal/mol to the activation energy in the reaction phase 6, assisted by the charge polarization of both H atoms.

Chemical Bond Analysis
The URVA analysis was complemented with LMA and QTAIM bond analyses for all stationary points of reactions R0-R5 in order provide further insides into the mechanism and to confirm the URVA analysis.Table 3 presents bond length R, local mode force constant k a and related bond strength order BSO and energy density H r for those bonds of reactants RE, transition states TS and products PR which change their character during the reaction.In order to limit the number of entries reported in this table, we assumed that the product of the previous catalyzed reaction is the same as the reactant of the next reaction in the catalytic cycle, although the IRC calculations did finish at slightly different geometries.and 4.805 mDyn/Å for RE, TS and PR, respectively), its length is decreasing (3.128, 1.447, and 1.108 Å for RE, TS and PR, respectively), and its covalent character is increasing ( 0.0666 and 0.2896 Hr/Bohr 3 , for TS and PR, respectively).Similarly, the strength of the O a H a bond is increasing (0.017, 1.686, and 8.005 mDyn/Å for RE, TS and PR, respectively), its length is decreasing (2.959, 1.308, and 0.969 Å for RE, TS and PR, respectively), and its covalent character is increasing (0.0010, 0.1021, and 0.5918 Hr/Bohr 3 for RE, TS and PR, respectively).Finally, the H a H b bond is cleaved in this reaction, which is confirmed by its decreasing strength (5.898, 1.210, and 0.469 mDyn/Å for RE, TS and PR, respectively), by its increasing length (0.743 1.030, and 2.120 Å for RE, TS and PR, respectively), and by its decreasing covalent character ( 0.3364 and 0.0861 Hr/Bohr 3 for RE and TS, respectively).In summary, the LMA and QTAIM analyses taken at the stationary points are completely in line with the URVA analysis.3).The C a O a bond of CO 2 keeps more or less its double character during reaction R1, although its strength is somewhat decreasing as reflected by the k a values (16.098, 9.679, and 10.562 mDyn/Å for RE, TS and PR, respectively), increasing bond length (1.170, 1.242, and 1.240 Å for RE, TS and PR, respectively) and decreasing covalent character ( 0.7449, 0.6444, 0.6432 Hr/Bohr 3 for RE, TS and PR, respectively).In contrast, the C a O b bond of CO 2 is changing its character from a double to a single bond upon CO 2 addition to the metal.The strength of the C a O b bond is decreasing (16.039, 8.344, and 8.253 mDyn/Å for RE, TS and PR, respectively), the bond length is increasing (1.170 1.268, and 1.283 Å for RE, TS and PR, respectively) and the covalent character is decreasing ( 0.7436, 0.6083, and 0.5787 Hr/Bohr 3 for RE, TS and PR, respectively).The changes of the C a O b bond are a consequences of the IrO b bond formation in this reaction.The strength of the IrO b bond is consequently increasing (k a = 0.049, 0.415, and 1.184 mDyn/Å for RE, TS and PR, respectively) and its length is decreasing (4.262, 2.893, and 2.222 Å for RE, TS and PR, respectively).We observed only for the PR IrO b bond a bond critical with a small covalent character ( 0.0038 Hr/Bohr 3 ).The reaction R1 also involves a H atom transfer from the catalyst to the incoming CO 2 ligand, where the IrH b bond of the catalyst is cleaved, and the new C a H b bond is formed.Consequently, the strength of the IrH b bond is decreasing (k a = 1.983 and 0.415 mDyn/Å for RE and TS, respectively), its length is increasing (1.679, 2.359, and 4.079 Å for RE, TS and PR, respectively), and its covalent character is decreasing ( 0.0574 and 0.0003 Hr/Bohr 3 for RE and TS, respectively).The formation of the C a H b bond is reflected by increasing bond strength (0.041, 3.090, and 4.692 mDyn/Å for RE, TS and PR, respectively), decreasing bond length (2.769 1.153, and 1.111 Å for RE, TS and PR, respectively) and increasing covalent character (0.0014, 0.2259, and 0.2828 Hr/Bohr 3 for RE, TS and PR, respectively).In conclusion, we observe that the results of the LMA and QTAIM analyses of the reaction R1 are in support of the results of the URVA analysis, similarly as for the reaction R0.
Reaction R2 Reaction R2 involves H atom transfer from the catalyst to the CO 2 ligand and subsequent OH bond formation implying the cleavage of a catalyst CH bond.There are six bonds involved in this process, C a O a , C a O b , C a H b , C b H a , O a H a , and IrO b (see Table 3).As the H a atom is transferred from C b to the CO 2 ligand, the strength of the C b H a bond is decreasing (k a = 5.237, 0.138, and 0.117 mDyn/Å for RE, TS and PR, respectively), and its length is consequently increasing (1.094, 2.623, and 3.084 Å for RE, TS and PR, respectively).We observed only a bond critical point for the RE C b H a bond with a relatively large covalent character ( 0.2896 Hr/Bohr 3 ).The change of O a H a bond properties are going into the opposite direction.The O a H a bond strength is increasing (k a = 0.029, 6.893, and 4.053 mDyn/Å for RE, TS and PR, respectively), the bond length is decreasing (2.403, 0.985, and 1.011 Å for RE, TS and PR, respectively), and the covalent character is increasing (0.0020, 0.5652, and 0.4934 Hr/Bohr 3 for RE, TS and PR, respectively).It is interesting to note that at the TS the O a H a bond is stronger (6.893 mDyn/Å), shorter (0.985 Å) and more covalent ( 0.5652 Hr/Bohr 3 ) than in the PR.As the reaction proceeds, the C a O a bond is changing its character from double bond to a single bond, which is reflected in the RE, TS, and PR bond properties.The strength of the C a O a bond is decreasing (k a = 10.562,7.086, and 7.134 mDyn/Å for RE, TS and PR, respectively), its length is increasing (1.240, 1.318, and 1.311 Å for RE, TS and PR, respectively), and its covalent character is decreasing ( 0.6432, 0.5196, and 0.5320 Hr/Bohr 3    3).The IrO b bond strength is decreasing (k a = 0.635, 0.076, and 0.050 mDyn/Å for RE, TS and PR, respectively), the bond becomes longer (2.333, 3.476, and 4.085 Å for RE, TS and PR, respectively) and less covalent ( 0.0002 and 0.0004 Hr/Bohr 3 for RE and TS , respectively; no bond critical point was found for PR).Cleavage of the IrO b bond affects the properties of the C a O b bond, which becomes stronger (k a = 11.316,12.883, and 13.069 mDyn/Å for RE, TS and PR, respectively), shorter (1.229, 1.210, and 1.209 Å for RE, TS and PR, respectively), and more covalent ( 0.6671, 0.6978, and 0.6967 Hr/Bohr 3 for RE, TS and PR, respectively).In contrast, the C a O a bond becomes weaker (k a = 7.134, 6.260, and 6.280 mDyn/Å for RE, TS and PR, respectively), longer (1.311, 1.337, and 1.338 Å for RE, TS and PR, respectively) and less covalent ( 0.5320, 0.4924, and 0.4907 Hr/Bohr 3 for RE, TS and PR, respectively), although the overall changes of C a O a bond properties are smaller than those for the C a O b bond.Change of the metal coordination number from 6 to 5 during this reaction implies a change of the H c and H d ligand positions in the metal coordination sphere, which affects the IrH c and IrH d bond properties.The IrH c bond becomes weaker (3.139, 3.129, and 2.758 mDyn/Å for RE, TS and PR, respectively), longer (1.564, 1.548, and 1.576 Å for RE, TS and PR, respectively), and less covalent ( 0.0981, 0.1078, and 0.0949 Hr/Bohr 3 for RE, TS and PR, respectively).However, although the IrH d bond becomes somewhat weaker (k a = 2.833, 2.737, and 2.728 mDyn/Å for RE, TS and PR, respectively), it becomes shorter (1.590, 1.586, and 1.578 Å for RE, TS and PR, respectively), and slightly more covalent ( 0.0874, 0.0900, and 0.0939 Hr/Bohr 3 for RE, TS and PR, respectively).The values of the IrH c and IrH d PR bond properties are similar (k a = 2.758 and 2.728 mDyn/Å, R = 1.576 and 1.578 Å, and H r = 0.0949 and 0.0939 Hr/Bohr 3 , for IrH c and IrH d , respectively), which is consistent with their similar positions in the metal coordination sphere.There is an additional NH bond involved in reaction R3 with hydrogen bond character N• • • H a O a , which is cleaved during the reaction.The strength of the NH a bond continuously decreases (k a = 0.219, 0.129, and 0.029 mDyn/Å for RE, TS and PR, respectively), the length increases (1.747, 2.031, and 3.225 Å for RE, TS and PR, respectively), and the covalent character decreases ( 0.0040 and 0.0007 Hr/Bohr 3 for RE and TS R, respectively, no bond critical point was found for PR.In conclusion, the results of the LMA and QTAIM analyses are consistent with the results of the URVA analysis. Reaction R4 Reaction R4 involves the h 2 -H 2 coordination of H 2 to the metal center, the first step of catalyst recovery.The coordination number of the metal is changed back from 5 to 6 and two new bonds are formed between the metal and H 2 , i.e., IrH a and IrH b , changing the properties of the H a H b bond from the incoming H 2 ligand.H 2 addition also changes properties of the IrH c and IrH d bonds caused by the change of the H c and H d ligand positions in the coordination sphere.Five bonds are involved in this process, H a H b , IrH a , IrH b , IrH c and IrH d (see Table 3).The strength of the IrH a bond is increasing (k a = 0.004, 0.071, and 0.730 mDyn/Å) for RE, TS and PR, respectively) and its length is decreasing (4.201, 3.210, and 1.829 Å for RE, TS and PR, respectively).No bond critical point could be found for any of the three stationary points.Similarly, the strength of the IrH b bond is increasing as well (k a = 0.004, 0.050, and 0.807 mDyn/Å for RE, TS and PR, respectively), the length is decreasing (3.461 2.597, and 1.787 mDyn/Å for RE, TS and PR, respectively), and its covalent character is increasing (0.0004, 0.0001, and 0.0275 Hr/Bohr 3 for RE, TS and PR, respectively).The H 2 addition to the metal changes properties of the H a H b bond, which becomes weaker (k a = 5.744, 5.371, and 1.826 mDyn/Å for RE, TS and PR, respectively), longer (0.745 0.750, and 0.842 Å for RE, TS and PR, respectively) and less covalent ( 0.3330, 0.3225, and 0.2020 Hr/Bohr 3 for RE, TS and PR, respectively).This is important for the final step of the catalytic cycle, in which the H a H b bond has to be cleaved.The H c and H d atoms are moving into new positions, affecting the IrH c bond which becomes slightly stronger (k a = 2.758, 3.318, and 2.783 mDyn/Å for RE, TS and PR, respectively), slightly longer (1.576, 1.545, and 1.589 Å for RE, TS and PR, respectively), and less covalent ( 0.0949, 0.1091, and 0.0889 Hr/Bohr 3 for RE, TS and PR, respectively).Similar changes were found for the IrH d bond.The bond strength is increasing (k a = 2.728, 2.760, and 2.762 mDyn/Å for RE, TS and PR, respectively), the bond length is increasing (1.578, 1.595, and 1.604 Å for RE, TS and PR, respectively), and the covalent character is decreasing ( 0.0939, 0.0859, and 0.0829 Hr/Bohr 3 for RE, TS and PR, respectively).These results fully support the URVA analysis.

Reaction R5
The last reaction of the catalytic cycle involves the dissociation of the H a H b bond from the H 2 ligand, the transfer of the H a atom from the metal center to the C b carbon atom of the catalyst accompanied by the cleavage of the IrH a bond, and the strengthening of the two IrH b and C b H a bonds, regenerating the original catalyst.Six bonds are involved in this process, H a H b , IrH a , IrH b , IrH c , IrH d and C b H a (see Table 3).The strength of the already weakened H a H b bond decreases further (k a = 1.826, 0.707, and 0.037 mDyn/Å for RE, TS and PR, respectively), its length increases (0.842, 0.986, and 3.182 Å for RE, TS and PR, respectively), and its covalent character decreases ( 0.2020 and 0.0907 Hr/Bohr 3 for RE and TS, respectively, no bond critical point was found for PR).The strength of the IrH a bond is decreasing (k a = 0.730, 0.355, and 0.107 mDyn/Å or for RE, TS and PR, respectively) and its length is increasing (1.829, 2.335, and 3.588 Å for RE, TS and PR, respectively) without bond critical points in this reaction.C b H a bond formation is confirmed by increasing bond strength (k a = 0.105, 0.932, and 5.118 mDyn/Å for RE, TS and PR, respectively)), decreasing bond length (3.481, 1.580, and 1.098 Å for RE, TS and PR, respectively)) and increasing covalent character ( 0.0371 and 0.2792 Hr/Bohr 3 for TS and PR, respectively, no bond critical point was found for RE).Strengthening of the IrH b bond is reflected by increasing k a values ( 0.807, 0.567, and 1.987 mDyn/Å or RE, TS and PR, respectively), decreasing bond length (1.787, 1.909, and 1.679 Å or RE, TS and PR, respectively), and increasing covalent character ( 0.0275, 0.0202, and 0.0581 Hr/Bohr 3 or RE, TS and PR, respectively).The strength of the IrH c bond is decreasing (k a = 2.783, 2.737, and 1.988 mDyn/Å for RE, TS and PR, respectively), the length is increasing (1.589, 1.591, and 1.679 Å for RE, TS and PR, respectively), and the covalent character is decreasing ( 0.0889, 0.0863, and 0.0581 Hr/Bohr In conclusion, the LMA and QTAIM analyses performed at all stationary points of reactions R0-R5 deepen the insights obtained with URVA by adding bond specific details.In the last section, overall trends of the LMA and QTAIM bond properties based on Table 3 data are presented to round up the results and discussion part.

Comparison of Relative Bond Strength Orders BSO
Local mode force constants k a can be transformed into relative bond strength orders BSO, which asses bond strength in a more commonly used language referring to single and multiple bond character of a bond, and in this way offering a more convenient way to compare bond strengths in aseries of molecules.BSO values were obtained described in the computational section, based on the reference molecules summarized in Table 1. Figure 10a-d 1).The smallest BSO values are in the range of 1.2 indicating single bond character, with the TS values in between.The weakest CO bond is found again for the PR of the non-catalyzed reaction R0.
Figure 10b shows the BSO values for C a H b and C b H a bonds of the catalyst, which are changing in reactions R0-R2, and R5.For example, the dissociating C b H a bond of reaction R2 starts with an RE BSO value of 0.821, which is close to a single bond bond strength and drops down in the PR to a value of 0.228, indicating a weak CH interaction, already shown in the TS with a BSO value of 0.241.The BSO values for the H a H b bonds are presented in Figure 10c for reactions R0, R4, and R5.The HH bond of the RE for both non-catalyzed reaction R0 and catalyzed reaction R4 are the strongest with BSO values in the range of 0.85, which are close to the BSO value of 0.855 of the HH bond of the free H 2 molecule (see Table 1.The weakest HH interaction is found for the RE of reaction R5 with a BSO value of 0.575, facilitating the cleavage of this bond.There is one OH formed in reactions R0 and R2.The small RE BSO value of 0.119 reflects that there is only a weak interaction in the RE van der Waals complex.The BSO value increases to 0.560 at the TS and the PR BSO of 0.947 indicates OH single bond character upon OH bond formation.Figure 11a-b  Corresponding BSO values are presented in Figure 11a.For example, in reaction R1 the IrO b bond is formed, starting with a small RE BSO value of 0.107, which increases to 0.352 for the TS and to 0.633 for the PR.It is interesting to note that the overall BSO values of all IrO bonds in this study are smaller than 1, indicating relatively weak bond strength character.Figure 11b shows the BSO values of the IrH bonds, such as IrH a , IrH b , IrH c , and IrH d , which are formed, cleaved, and changed in reactions R1, and R3-R5.For example, in reaction R3, upon dissociation of the HCOOH ligand, the IrH c bond is only changing slightly its character, with a BSO values of 0.772 for the RE, 0.770 for the TS, and 0.710 for the PR.

Correlations between Bond Properties and Local Mode Force Constants
Figures 12 and 13 present the correlation between bond strength expressed by the local mode force constants and bond length, which according to the Badger rule [139,170] leads to a power relationship.Figure 12a shows the correlation for CO bonds of reactions R0-R3.According to Figure 12a, stronger CO bonds are connected to shorter CO bonds.The R 2 value of 0.9858 confirms a significant correlation between these two properties.larger divergencies are observed for the C a O a bond at TS of reaction R0, and C a O a and C a O b bonds at the TSs of reaction R1. Figure 12b shows a similar correlation between bond length and bond strength for the CH bonds in reactions R0-R2, and R5.The R 2 factor of this correlation has a somewhat lower value of 0.9310, corresponding to a good correlation.The divergency occurs for the C a H b bond at the TS of reaction R0.The correlation between bond strength and bond length for the HH bonds are presented in Figure 12c, showing the results for the H a H b bonds in reactions R0, R4, and R5.This correlation has an R 2 value of 0.8801, with outliers for PRs of reactions R0 and R4 and the RE of reaction R5. Figure 12d shows a significant correlation between bond strength and bond length for the O a H a bonds of reactions R0 and R2.The value of the R 2 factor of this correlation is 0.9898 with no substantial outlier.The bond strength bond versus length correlation for the IrO b bonds in reactions R1-R3 is presented in Figure 13a.The R 2 factor of this correlation has a moderate value 0.9545, with the largest outlier for the TSs of reactions R1 and R2. Figure 13b shows the correlation between bond strength and bond length for the IrH bonds of reactions R1, and R3-R5.The R 2 value of 0.9099 indicates a moderate correlation, smaller than that found for the IrO bonds.Largest outliers are found for the IrH b TS bond of reaction R1, the IrH a and IrH b TS bonds of reactions R5, and for the PR IrH a and IrH c bonds of reaction R4.In summary, we observed a satisfactory correlation between bond strength and bond length, in line with the popular rule that the stronger bond is also the shorter bond.However, this is not always the case as demonstrated by the increasing number of exceptions from this simple rule, reported in the literature [87,88,90,171,172].The outliers in the correlations presented in this work are mostly related to the TS stationary points.At the TS, in particular for bond breaking/forming processes, significant electronic structure rearrangements take place, e.g., leading to a stage between two different hybridization forms, which can cause an unusual bonding situation that cannot be captured by such an empirical correlation.
Energy density versus Local mode force constant As outlined above, the energy density H r is convenient measure of the covalent character of a chemical bond or weak chemical interaction based on the electron density.Overall, we observed in our study that the stronger bonds have more covalent character, expressed by a more negative value of H r .Figures 14 and 15 present the correlations between H r and k a based on a linear relationship.Figure 14a shows the correlation for the C a O a and C a O b bonds of reactions R0-R3.The R 2 factor has a value of 0.9411, with a number of small outliers mostly for TSs.The bond strength-energy density correlation for the C a H b and C b H a bonds of reactions R0-R2, is presented in Figure 14b.The R 2 value of 0.9758 reflects a good correlation, with the largest outlier for the TS C a H b bond of reaction R1. Figure 14c shows the same correlation for the H a H b bonds of reactions R0, R4, and R5, with a R 2 value of 0.9438.The largest outlier was found for the PR H a H b bond of reaction R4.The correlation between energy density and bond strength for the O a H a bonds of reactions R0 and R2 is shown in Figure 14d.The R 2 factor for this reaction, has a value of 0.9256 reflecting a moderate correlation with the largest outlier for the PR O a H a bond of reaction R2. Figure 15a presents the correlation between H r and k a for the IrO b bonds in reactions R1-R3.As obvious from the R 2 factor of 0.7816, there is only some trend between these two properties.One also has to consider the scarcity of the data points.It is interesting to note that the PR IrO b bond of reaction R2 has an H r close to zero, indicating a more electrostatic character of this bond.A similar situation occurs for the TS IrO b bond of reaction R3.A better correlation between energy density and bond strength is observed for the IrH bonds of reactions R1 and R3-R5 shown in Figure 15b.The R 2 factor of that correlation has a value of 0.9829, with some outlier such as the TS IrH b bond of the reaction R1, or the TS IrH c bond of reaction R3.
In summary, overall, we observed a good correlation between energy density H r and bond strength expressed by the local mode force constant k a , which indicates that the strength of the bonds investigated in this work is predominantly based on their covalent character.However, a caveat is appropriate.H r reflects bonding at just one point of the electron density, namely the bond critical point, whereas the k a as a second order property picks up the electronic environment [91].As such it is more sensitive in the case of complex bonding situation as, e.g., found for the IrO bonds showing only a trend but not a significant correlation between these two properties.

Conclusions and Outlook
We investigated in this study the catalytic cycle (reactions R1-R5) of carbon dioxide hydrogenation catalyzed with an iridium pincer model complex, which reduces the unfavorable activation energy of 76 kcal/mol for the non-catalyzed gas phase reaction R0 by more than 100%.As a computational tool we used URVA to explore the mechanistic features of each reaction step in all detail, shedding new light on which are the most energy consuming events in the non-catalyzed reaction, and how are they avoided in the catalyzed reaction.URVA was complemented by LMA and QTAIM analyses at all stationary points of each reaction.We used the DFT level of theory for URVA, LMA, and QTAIM, whereas for the energetics a dual approach (i.e., single point DLPNO-CCSD(T) energy calculations for all stationary points based on the DFT geometries) was applied.The most important findings of our study showing for the first time how the catalyst effectively manages to keep energy demand at a balance are summarized in the following list:

•
According to the URVA analysis, the most energy consuming process of the noncatalyzed reaction R0 is the cleavage of the HH bond of the H 2 reactant, which takes place before TS, and as such contributes to the activation energy with 40.8 kcal/mol.

•
In order to avoid direct HH bond cleavage, the catalyst divides the one-step noncatalytic reaction into a multistep catalytic cycle; CO 2 addition to the catalyst, H atom transfer between catalyst and CO 2 ligand, release of HCOOH product, addition of H 2 , reorganization of the coordination sphere of the catalyst to achieve its original form.Each of these decisive steps could be clearly identified in the characteristic curvature profiles for the first time.

•
The expensive cleavage of the HH bond in the non-catalyzed reaction, is replaced in the catalyzed reaction by H transfer (reaction R2) requiring the cleavage of an IrH bond with a significantly smaller contribution of 9.8 kcal/mol to the activation energy, which is revealed by both energy and curvature profiles.

•
The dissociation of the final product from the catalyst (reaction R3) is characterized by the cleavage of an IrO bond and an intermediate NH hydrogen bond, and according to our URVA analysis, both contribute to the activation energy with a moderate amount of 12.3 kcal/mol.

•
As unravelled by the URVA curvature profiles most of the events related to the reorganization of the catalyst to restore its original form (reaction R5) occur after the TS, i.e., they do not contribute to the activation energy.
In summary, the URVA analysis complimented with the LMA and QTAIM methods provided new comprehensive mechanistic details of all reaction steps forming the catalytic cycle of the hydrogenation of carbon dioxide, which hopefully will stimulate and inspire the community working on this important topic.Experiments suggest that the catalytic activity of iridium pincer complexes for the hydrogenation of carbon dioxide can be considerably increased under basic reaction conditions [43,44].A follow-up URVA study is planned to explore this.

Supplementary Materials:
The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/inorganics10120234/s1,containing reaction movies R0-R5 and the Cartesian coordinates of the stationary point, i.e, the reactant complex, TS and product complex for all six reactions R0-R5.
Funding: This research was funded by the National Science Foundation NSF, grant CHE 2102461.

Figure 2 .
Figure 2. Schematic representation of the reaction profile for a model reaction defined as the scalar curvature k(s) (solid black line) given as a function of the reaction parameter s.Curvature minima M and curvature maxima K are shown.The location of the TS is denoted by a dotted blue line.Start and end of each reaction phase are denoted by dashed black lines.In addition the decomposition of the scalar curvature k(s) into four components is shown, blue, red, green and purple colored lines.

Figure 3 .
Figure 3. Energy diagram for the catalyzed reactions R1-R5, relative to the initial reactant van der Waals complex of 1.For simplicity, the energies have been calculated by adding the activation and reaction energy of the next reaction step of the cycle, to the product energy of the previous reaction step.DLPNO-CCSD(T)/def2-TZVP/def2-ECP(Ir) level of theory performed at the end points of the DFT IRC calculations with B3LYP/6-31G(d,p)/SDD(Ir).
. Reaction R1 involves addition of CO 2 to the coordination sphere of iridium and hydrogen atom transfer from Ir to one of the CO 2 carbon atoms.According to Figure 5b the reaction starts in phase 3 with the initial stage of the IrO b bond formation (blue line) first resisting, followed by the next stage in phase 4 where it becomes supporting, whereas the final step of IrO b bond formation occurs after TS in phase 9.The different stages of IrO bond formation are in line with changes in IrO b bond lengths shown in Figure 5c.In phase 3 the IrO b bond distance starts to decrease, reaching its final value in phase 9. Figure 5d presents the corresponding changes of Ir and O b NBO charges taking place between reaction phases 4 and 9.The accompanying transfer of the H atom from Ir to the C a carbon of CO 2 implies IrH b bond cleavage and C a H b bond formation.According to Figure 5b, the transfer process (red line) starts at the end of phase 4, where the formation of the new C a H b bond (green line) is still resisting, changing into a large supporting contribution in phase 5.

Figure 4 .
Figure 4. Properties of the gas phase reaction R0 along the reaction path.(a) Energy profile, energy contributions of each reaction phase to the activation energy are indicated by italic numbers; (b) decomposition of the reaction curvature into selected components; (c) change of selected bond lengths along s; (d) change of selected NBO atomic charges along s.Positions of the curvature minima are shown as dashed vertical lines and are labeled as M1, M2, and so on.The position of the TS is indicated by a dotted line.Reaction phases are indicated by blue numbers.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.The next stage of H b transfer from Ir to C a occurs in phase 5, being finalized far out in the exit channel in phases 8-10.Figure 5c,d confirm these processes.Largest C a H b bond length changes are found in phases 8 and 9 after the TS, as well as the changes in the C a and H b NBO charges.Figure 5b also reveals the involvement of the O a C a O b bond angle (purple line) in the reaction mechanism, with largest contribution between phases 4 and 9, in line with the OCO bond angle changes from 180 to 120 occurring between phases 4 and 9 (see Figure 5c) which is facilitated by the coordination of C b to the metal atom.According to Figure 5a, the two largest energy contributions to TS energy are 4.6 kcal/mol in phase 4, caused by IrO b and C a H b bond formation, together with changes of the O a C a O b bond angle, and 6.2 kcal/mol in phase 8, connected with IrH b bond cleavage and IrO b bond formation, which in summary is far less than the direct cleavage of the HH bond in reaction R0.

Figure 5 .
Figure 5. Properties of reaction R1 of the catalytic cycle along the reaction path.(a) Energy profile, energy contributions of each reaction phase to the activation energy are indicated by italic numbers; (b) decomposition of the reaction curvature into selected components; (c) change of selected bond lengths along s; (d) change of selected NBO atomic charges along s.Positions of the curvature minima are shown as dashed vertical lines and are labeled as M1, M2, and so on.The position of the TS is indicated by a dotted line.Reaction phases are indicated by blue numbers.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Figure 6 .
Figure 6.Properties of reaction R2 of the catalytic cycle along the reaction path.(a) Energy profile, energy contributions of each reaction phase to the activation energy are indicated by italic numbers; (b) decomposition of the reaction curvature into selected components; (c) change of selected bond lengths along s; (d) change of selected NBO atomic charges along s.Positions of the curvature minima are shown as dashed vertical lines and are labeled as M1, M2, and so on.The position of the TS is indicated by a dotted line.Reaction phases are indicated by blue numbers.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.
1.75 Å to 3.2 Å.The change of the metal coordination sphere being related to the change of the H c and H d ligand positions is reflected by IrH c (orange line) and IrH d (red line) curvature components in phases 3 and 4. The contribution of the IrH d component is larger, indicating a stronger electronic structure change of the metal caused by displacement of this equatorial ligand.As shown in Figure 7a, the largest energy contribution of 7.0 kcal/mol to the activation energy results from the first stage of IrO b bond cleavage occurring in phase 2. The second stage of IrO b bond cleavage in phase 3, together with IrH c and IrH d ligand reorganization and NH a hydrogen bond cleavage contribute 5.3 kcal/mol to the activation energy.

Figure 7 .
Figure 7. Properties of reaction R3 of the catalytic cycle along the reaction path.(a) Energy profile, energy contributions of each reaction phase to the activation energy are indicated by italic numbers; (b) decomposition of the reaction curvature into selected components; (c) change of selected bond lengths along s; (d) change of selected NBO atomic charges along s.Positions of the curvature minima are shown as dashed vertical lines and are labeled as M1, M2, and so on.The position of the TS is indicated by a dotted line.Reaction phases are indicated by blue numbers.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Figure 8 .
Figure 8. Properties of reaction R4 of the catalytic cycle along the reaction path.(a) Energy profile, energy contributions of each reaction phase to the activation energy are indicated by italic numbers; (b) decomposition of the reaction curvature into selected components; (c) change of selected bond lengths along s; (d) change of selected NBO atomic charges along s.Positions of the curvature minima are shown as dashed vertical lines and are labeled as M1, M2, and so on.The position of the TS is indicated by a dotted line.Reaction phases are indicated by blue numbers.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Figure 9 .
Figure 9. Properties of reaction R5 of the catalytic cycle along the reaction path.(a) Energy profile, energy contributions of each reaction phase to the activation energy are indicated by italic numbers; (b) decomposition of the reaction curvature into selected components; (c) change of selected bond lengths along s; (d) change of selected NBO atomic charges along s.Positions of the curvature minima are shown as dashed vertical lines and are labeled as M1, M2, and so on.The position of the TS is indicated by a dotted line.Reaction phases are indicated by blue numbers.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.
Reaction R1 Reaction R1 involves addition of CO 2 to the metal center forming the IrO b bond, H atom transfer from the metal to the CO 2 ligand via cleavage of the IrH b bond, and the formation of the C a H b bond.Five bonds are involved in this process, namely C a O a , C a O b , C a H b , IrO b , and IrH b (see Table Reaction R3 Reaction R3 involves the dissociation of the HCOOH product from the metal coordination sphere via IrO b bond cleavage, affecting also the properties of both C a O a and C a O b bonds.In addition, this process changes the metal coordination numbers from 6 to 5, which is accompanied by the change of H c and H d atomic positions and the change of IrH c and IrH d bond properties.As such, six bonds are involved, C a O a , C a O b , IrO b , IrH c , IrH d , and NH a (see Table show BSO values as a function of the local mode force constants k a , for the CO, CH, HH, and OH bonds and Figure 11a-b for IrO and IrH bonds.There are two types of CO bonds, C a O a and C a O b , which are changing their bond strengths in reactions R0-R3.According to Figure 10a, largest CO BSO values are in the range of 2.2 indicating strong double bond character.The strongest strongest CO bonds were found for the RE of the non-catalyzed reaction R0 with BSO values close to the BSO of 2.271 of the CO bond in the CO 2 molecule (see Table present plots of the BSO values for the IrO and IrH bonds.The IrO b bonds change in reactions R1-R3.

Figure 10 .
Figure 10.BSO values as a function of the local mode force constants k a , calculated from the power relationship BSO = A ⇤ (k a ) B ; for details, see text and Table 1.(a) CO bonds; (b) CH bonds; (c) HH bonds; (d) OH bonds.Labels R0C a O a R, R0C a O a T, and R0C a O a P indicate reactant RE (red color), TS (blue color), and product PR (green color) of the bond involving the C a and O a atoms of reaction R0.Similarly notation for the reactions R1-R5.For the bond label description see the text.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Figure 11 .
Figure 11.BSO values as a function of the local mode force constants k a , calculated from the power relationship BSO = A ⇤ (k a ) B ; for details, see text and Table 1.(a) IrO bonds; (b) IrH bonds.Labels R1IrO b R, R1IrO b T, and R1IrO b O indicate reactant RE (red color), TS (blue color), and product PR (green color) of the bond involving the Ir and O b atoms of thereaction R1.Similarly notation for reactions R2-R5.For the bond label description see the text.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Figure 12 .
Figure 12.Correlation between bond lengths and local mode force constants.(a) CO bonds; (b) CH bonds; (c) HH bonds; (d) OH bonds.Labels such as R0C a O a R, R0C a O a T, and R0C a O a P indicate reactant RE (red color), TS (blue color), and product PR (green color) C a O a bond in reaction R0.Similarly notation for reactions R1-R5.For the bond label description see the text.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Figure 13 .
Figure 13.Correlation between bond length and the local mode force constants.(a) IrO bon; (b) IrH bonds.Labels such as R1IrO b R, R1IrO b T, and R1IrO b O indicate reactant RE (red color), TS (blue color), and product PR (green color) IrO b bond in reaction R1.Similarly notation for reactions R2-R5.For the bond label description see the text.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Figure 14 .
Figure 14.Correlation between energy density H r and local mode force constant k a .(a) CO bonds; (b) CH bonds; (c) HH bonds; (d) OH bonds.Labels such as R0C a O a R, R0C a O a T, and R0C a O a P indicate reactant RE (red color), TS (blue color), and product PR (green color) C a O a in reaction R0.Similar notation for reactions R1-R5.For the bond label description see the text.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Figure 15 .
Figure 15.Correlation between energy density H r and local mode force constant k a .(a) IrO bonds; (b) IrH bonds.label such as R1IrO b R, R1IrO b T, and R1IrO b O indicate reactant RE (red color), TS (blue color), and product PR (green color) IrO b bond in R1.Similar notation for reactions R2-R5.For the bond label description see the text.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Table 1 .
Bond length R, local mode force constant k a , bond strength order BSO, and energy density H r evaluated at the bond critical point r(r) of the reference bonds and reference molecules used in our study.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.

Table 3 .
Bond length R, local mode force constant k a , bond strength order BSO, and energy density H r of the selected bonds for reactant RE, transition state TS, and product PR of reactions R0-R5.For the bond label description, see the text.B3LYP/6-31G(d,p)/SDD(Ir) level of theory.Reaction R0 Reaction R0 includes HH bond cleavage, formation of a new CH and OH bond, and the change of CO bond character in CO 2 .There are five bonds involved in this reaction, namely C a O a , C a O b , C a H b , O a H a , and H a H b (see Table 3).The C a O O b bond of CO 2 keeps its double character during this reaction, however, its strength is slightly decreasing (16.402, 14.830, and 13.994 mDyn/Å, for RE, TS and PR, respectively), the bond length is slightly increasing (1.169, 1.176, and 1.198 Å for RE, TS and PR, respectively), and the covalent character is slightly decreasing ( 0.7597, 0.7484, and 0.7315 Hr/Bohr 3 for RE, TS and PR, respectively).Two bonds are formed in this reaction, namely C a H b and O a H a , while the H a H b bond is cleaved.Consequently, the strength of the C a H b bond is increasing (0.019, 0.865, RE, TS and PR, respectively) and the decreasing covalent character ( 0.7594, 0.6537, and 0.4670, Hr/Bohr 3 for RE, TS and PR, respectively).In contrast, the C a for RE, TS and PR, respectively).C a O b bond properties are going in the opposite direction, as the character of this bond is changing from a single to a double bond.The strength of the C a O b bond is increasing (k a = 8.253, 11.827, and 11.316 mDyn/Å for RE, TS and PR, respectively), its length is decreasing (1.283 1.224, and 1.229 Å RE, TS and PR, respectively), and its covalent character is increasing ( 0.5787, 0.6731, and 0.6671 Hr/Bohr 3 for RE, TS and PR, respectively).The change of the C a O b bond properties, affects in turn the properties of the IrO b bond.The IrO b bond becomes weaker (k a = 1.184, 0.573, and 0.635 mDyn/Å for RE, TS and PR, respectively), longer (2.222, 2.351, and 2.333 Å for RE, TS and PR, respectively), and less covalent ( 0.0038, 0.0002, and 0.0002 Hr/Bohr 3 for RE, TS and PR, respectively).The strength of the C a H b bond increases only slightly, as reflected by its k a values (4.692, 5.302, and 5.270 mDyn/Å for RE, TS and PR, respectively), its bond length (1.111, 1.096, and 1.097 Å for RE, TS and PR, respectively), and its covalent character ( 0.2828, 0.3133, and 0.3115 Hr/Bohr 3 for RE, TS and PR, respectively).Again, these results are in line with the URVA analysis.
3for RE, TS and PR, respectively) whereas the properties of the IrH d bond are changing in the opposite direction.The strength of this bond increases (k a = 2.762, 2.842, and 2.914 mDyn/Å for RE, TS and PR, respectively), the length decreases (1.604, 1.592, and 1.585 Å for RE, TS and PR, respectively), and the covalent character increases ( 0.0829, 0.0869, and 0.0894 Hr/Bohr 3 RE, TS and PR, respectively).Again, these results fully support the URVA analysis.