C-H versus C-O Addition: A DFT Study of the Catalytic Cleavage of the β -O-4 Ether Linkage in Lignin by Iridium and Cobalt Pincer Complexes

: The potential energy surfaces of the reactions involved in the catalytic cleavage of 2 - phenoxy-1 -phenylethanol, a model of the β -O-4 linkage in lignin, by ( ipr PCP)-Ir, ( ipr PCOP)-Ir, ( ipr PCP)- Co and ( ipr PCOP)-Co complexes have been studied using the M06/6-311G**/LANL2TZ level of theory. Both iridium and cobalt are found to be active towards the cleavage of the β -O-4 linkage, with rate constants of 44.7 s − 1 and 5.1 × 10 6 s − 1 , respectively. The iridium catalysts prefer the ‘initial C-H addition’ pathway, showing a kinetic preference of 16.8 kcal mol − 1 over the ‘direct C-O insertion’ pathway, while the cobalt catalysts prefer the ‘direct C-O insertion’ route which is kinetically favored by 15.7 kcal mol − 1 over the ‘initial C-H addition’ pathway. A two-state reactivity occurs along the preferred pathway for the cobalt-catalyzed reaction.


Introduction
The mono-anionic ECE-type pincer systems are very common due to their comparative ease of synthesis.Known to be planar with a single aromatic backbone, they possess a central aryl anionic carbon 'C' and ortho substituents bearing side arm donor groups 'E', where E is N or P.These pincer ligands make up an attractive niche of ligand frameworks that provide a suitable balance between stability and reactivity and have become widely used in homogenous catalysis [1].For example, the complex (PCP)Ir (PCP = κ3-C 6 H 3 -2,6-[CH 2 P(t-Bu) 2 ] 2 ) and its derivatives are known as the most effective alkane dehydrogenation catalysts developed to date [2][3][4].Pincer complexes have also been shown to be catalytically active in the activation of C-X bonds (X = C, O, N, F etc.).In 2008, Chirik et al. reported the cleavage of C-O bonds in alkyl-substituted esters via binuclear oxidative addition using ( ipr PDI)Fe(N 2 ) 2 (Scheme 1) [5].Other noteworthy examples of sp 3 C-O bond cleavage include the rearrangement reactions of methyl aryl ethers by a tris(pyrazolyl)borate iridium complex reported by Carmona and Paneque [6][7][8], as well as reports by Ozerov and Grubbs on the cleavage of the tert-butyl-oxygen and benzyl-oxygen bonds of methyl ethers by a PNP-pincer iridium complex [9].Recently, Choi et al. [10] reported the oxidative addition of C(sp 3 -O) bonds of methyl esters, methyl tosylate and methyl aryl ethers using (PCP)Ir (PCP = κ3-C 6 H 3 -2,6-[CH 2 P(t-Bu) 2 ] 2 ).
These reports highlight some of the progress that has been made in addressing the challenges associated with the catalytic activation of sp 3 C-O bonds, one of the primary linkages in organic chemistry.The activation of these types of C-O bonds is also a focus area of interest for the catalytic valorization of lignin to deoxygenated fuels and commercial chemicals.Lignin, which makes up 40% of the weight of biomass, is mostly comprised of alkyl aryl ether bonds that are difficult to break due to their strength and stability.Current research focuses on understanding how the β-O-4 linkage, which is the most abundant linkage in lignin, can be cleaved [11].In 2014, Haibach et al. reported the dehydroaryloxylation of alkyl aryl ethers using PCP-type iridium complexes [12], making it the second report of this type of cleavage reaction being undertaken in a fully atomeconomic fashion, following in investigations carried out by Bergman et al. [13] A previous study conducted by Haibach et al. on the formation of ether C-O bonds through the same PCP-type iridium-catalyzed olefin hydroaryloxylation [14] led to their 2014 investigation after they discovered that the catalyst could also carry out the reverse reaction, i.e., 'ether dehydroaryloxylation'.From the optimization studies carried out by Haibach et al. [12], the use of less sterically bulky precursors ( ipr PCP)Ir A and ( ipr PCOP)Ir B (Scheme 2) was found to be beneficial as they led to much higher conversions.These reports highlight some of the progress that has been made in addressing the challenges associated with the catalytic activation of sp 3 C-O bonds, one of the primary linkages in organic chemistry.The activation of these types of C-O bonds is also a focus area of interest for the catalytic valorization of lignin to deoxygenated fuels and commercial chemicals.Lignin, which makes up 40% of the weight of biomass, is mostly comprised of alkyl aryl ether bonds that are difficult to break due to their strength and stability.Current research focuses on understanding how the β-O-4 linkage, which is the most abundant linkage in lignin, can be cleaved [11].In 2014, Haibach et al. reported the dehydroaryloxylation of alkyl aryl ethers using PCP-type iridium complexes [12], making it the second report of this type of cleavage reaction being undertaken in a fully atom-economic fashion, following in investigations carried out by Bergman et al. [13] A previous study conducted by Haibach et al. on the formation of ether C-O bonds through the same PCPtype iridium-catalyzed olefin hydroaryloxylation [14] led to their 2014 investigation after they discovered that the catalyst could also carry out the reverse reaction, i.e., 'ether dehydroaryloxylation'.From the optimization studies carried out by Haibach et al. [12], the use of less sterically bulky precursors ( ipr PCP)Ir A and ( ipr PCOP)Ir B (Scheme 2) was found to be beneficial as they led to much higher conversions.These reports highlight some of the progress that has been made in addressing the challenges associated with the catalytic activation of sp 3 C-O bonds, one of the primary linkages in organic chemistry.The activation of these types of C-O bonds is also a focus area of interest for the catalytic valorization of lignin to deoxygenated fuels and commercial chemicals.Lignin, which makes up 40% of the weight of biomass, is mostly comprised of alkyl aryl ether bonds that are difficult to break due to their strength and stability.Current research focuses on understanding how the β-O-4 linkage, which is the most abundant linkage in lignin, can be cleaved [11].In 2014, Haibach et al. reported the dehydroaryloxylation of alkyl aryl ethers using PCP-type iridium complexes [12], making it the second report of this type of cleavage reaction being undertaken in a fully atom-economic fashion, following in investigations carried out by Bergman et al. [13] A previous study conducted by Haibach et al. on the formation of ether C-O bonds through the same PCPtype iridium-catalyzed olefin hydroaryloxylation [14] led to their 2014 investigation after they discovered that the catalyst could also carry out the reverse reaction, i.e., 'ether dehydroaryloxylation'.From the optimization studies carried out by Haibach et al. [12], the use of less sterically bulky precursors ( ipr PCP)Ir A and ( ipr PCOP)Ir B (Scheme 2) was found to be beneficial as they led to much higher conversions.[10].Aryl C-H insertion was also observed to be kinetically favored in the reaction of RuH 2 CO(PPH 3 ) with aryl ethers [15].However, DFT studies backed by experimental findings reported by Chmely et al. showed the kinetic preference of the direct C-O activation over C-H activation in the cleavage of an alkyl aryl ether using a Ru-xantphos catalyst [16].Liu and Wilson also reported a DFT study on the reaction between PCP-type group 8 complexes with a model of the β-O-4 linkage in lignin, where the reaction proceeds via direct C-O bond cleavage [17].Consequently, in the effort to develop an efficient pincer-type catalyst for the cleavage of the alkyl-aryl ether bonds in lignin, the question of whether the reaction proceeds via initial C-H bond activation or via direct C-O bond insertion needs to be resolved.
Inspired by the work by Haibach and coworkers [11], here we report an exploratory mechanistic study on the catalytic activity of ( ipr PCP)Ir and ( ipr PCOP)Ir towards the cleavage of the β-O-4 linkage, 2-phenoxy-1-phenyl-ethanol (Scheme 3).We have investigated the energetics (kinetics and thermodynamics) of the reaction along three proposed pathways (Schemes 1-3) to determine whether the C-O bond cleavage proceeds via initial C-H activation or via direct C-O bond insertion.For comparison, we have also investigated the catalytic activity of ( ipr PCP)Co and ( ipr PCOP)Co for the cleavage reaction and compared the energetics with that of the iridium-catalyzed reaction.Cobalt-based pincer-type complexes are a relatively cheap and abundant alternative for catalytic applications, which have been applied successfully in a variety of dehydrogenation and hydrogenation reactions [18].
netically favored in the reaction of RuH2CO(PPH3) with aryl ethers [15].However, DFT studies backed by experimental findings reported by Chmely et al. showed the kinetic preference of the direct C-O activation over C-H activation in the cleavage of an alkyl aryl ether using a Ru-xantphos catalyst [16].Liu and Wilson also reported a DFT study on the reaction between PCP-type group 8 complexes with a model of the β-O-4 linkage in lignin, where the reaction proceeds via direct C-O bond cleavage [17].Consequently, in the effort to develop an efficient pincer-type catalyst for the cleavage of the alkyl-aryl ether bonds in lignin, the question of whether the reaction proceeds via initial C-H bond activation or via direct C-O bond insertion needs to be resolved.
Inspired by the work by Haibach and coworkers [11], here we report an exploratory mechanistic study on the catalytic activity of ( ipr PCP)Ir and ( ipr PCOP)Ir towards the cleavage of the β-O-4 linkage, 2-phenoxy-1-phenyl-ethanol (Scheme 3).We have investigated the energetics (kinetics and thermodynamics) of the reaction along three proposed pathways (Schemes 1-3) to determine whether the C-O bond cleavage proceeds via initial C-H activation or via direct C-O bond insertion.For comparison, we have also investigated the catalytic activity of ( ipr PCP)Co and ( ipr PCOP)Co for the cleavage reaction and compared the energetics with that of the iridium-catalyzed reaction.Cobalt-based pincer-type complexes are a relatively cheap and abundant alternative for catalytic applications, which have been applied successfully in a variety of dehydrogenation and hydrogenation reactions [18].The second mechanism as proposed by Liu and Wilson [17] begins with the coordination of the ether oxygen to the metal center to form an adduct.Direct C-O bond addition follows to form the five-coordinate intermediate B3, followed by the α-C-H bond addition to form B4. Finally, the phenol and enolate, which rearrange to acetophenone, are released to regenerate the Ir-pincer active species, as shown in Scheme 5.
The third proposed reaction pathway, shown in Scheme 6, is similar to the mechanism proposed by Wu et al. [19] in their experimental study of the cleavage of a β-O-4 model compound using a Ru-xantphos catalyst.It begins with dehydrogenation of the substrate 2-phenoxy-1-phenyl-ethanol 1 to its corresponding ketone motif C2.C2 undergoes rearrangement via TSC1 to a slightly more stable adduct with a three-membered ring formed between the ketone carbon, the ketone oxygen and the Ir center.This is followed by the direct C-O bond insertion to form C4, which undergoes hydrogenation using the H 2 generated earlier.The H 2 molecule, which binds to the metal to form a dihydrogen complex C5, is used in the subsequent reductive elimination of phenol and acetophenone and the Ir-pincer catalyst is regenerated to repeat the catalytic cycle.The second mechanism as proposed by Liu and Wilson [17] begins with the coordination of the ether oxygen to the metal center to form an adduct.Direct C-O bond addition follows to form the five-coordinate intermediate B3, followed by the α-C-H bond addition to form B4. Finally, the phenol and enolate, which rearrange to acetophenone, are released to regenerate the Ir-pincer active species, as shown in Scheme 5.The third proposed reaction pathway, shown in Scheme 6, is similar to the mechanism proposed by Wu et al. [19] in their experimental study of the cleavage of a β-O-4 model compound using a Ru-xantphos catalyst.It begins with dehydrogenation of the substrate 2-phenoxy-1-phenyl-ethanol 1 to its corresponding ketone motif C2.C2 undergoes rearrangement via TSC1 to a slightly more stable adduct with a three-membered ring formed between the ketone carbon, the ketone oxygen and the Ir center.This is fol- goes rearrangement via TSC1 to a slightly more stable adduct with a three-membered ring formed between the ketone carbon, the ketone oxygen and the Ir center.This is followed by the direct C-O bond insertion to form C4, which undergoes hydrogenation using the H2 generated earlier.The H2 molecule, which binds to the metal to form a dihydrogen complex C5, is used in the subsequent reductive elimination of phenol and acetophenone and the Ir-pincer catalyst is regenerated to repeat the catalytic cycle.Scheme 6. Proposed mechanism for the cleavage 2-phenoxy-1-phenyl-ethanol by pincer complexes along pathway C.

The Reaction of ( ipr PCP)Ir with 2-Phenoxy-1-phenylethanol 1
The optimized geometries of the stationary points, as well as the Gibbs free energy profile for the reaction of ( ipr PCP)Ir with 2-phenoxy-1-phenylethanol 1 along path A are shown in Figures 1 and 2. All stationary points in the catalytic cycle were sought on both the singlet and triplet state potential energy surface (PES), but triplet-state energies were found to be highly unstable, with the reactants being 104 kcal mol −1 less stable on the triplet surface than on the singlet surface (Table 1).

The Reaction of ( ipr PCP)Ir with 2-Phenoxy-1-phenylethanol 1
The optimized geometries of the stationary points, as well as the Gibbs free energy profile for the reaction of ( ipr PCP)Ir with 2-phenoxy-1-phenylethanol 1 along path A are shown in Figures 1 and 2. All stationary points in the catalytic cycle were sought on both the singlet and triplet state potential energy surface (PES), but triplet-state energies were found to be highly unstable, with the reactants being 104 kcal mol −1 less stable on the triplet surface than on the singlet surface (Table 1).The first step in pathway A, which is the C-H addition, proceeds with an activation energy barrier of 11.4 kcal mol −1 , leading to a five-coordinate complex A2.Both the hydrogen being extracted by the metal and the β-O-4 group sit axially to the pincer backbone in TSA1, and the transition state is a three-membered structure involving the α-carbon, the Table 1.Comparison of singlet-state and triplet-state energies of the stationary points involved in the reaction of ( ipr PCP)Ir with 1 along path A. The first step in pathway A, which is the C-H addition, proceeds with an activation energy barrier of 11.4 kcal mol −1 , leading to a five-coordinate complex A2.Both the hydrogen being extracted by the metal and the β-O-4 group sit axially to the pincer backbone in TSA1, and the transition state is a three-membered structure involving the α-carbon, the corresponding hydrogen and the metal.The Ir-C bond in the three-membered structure transitions from a bond length of 2.35 Å in TSA1 to 2.12 Å in A2 (d Ir-C = 2.12 Å) and both the abstracted hydrogen and the β-O-4 compound retain their axial position in A2.On the triplet PES, the activation energy barrier for the C-H addition is 34.1 kcal mol −1 and the resulting triplet intermediate is 21.1 kcal mol −1 less stable than the singlet state intermediate.The next step, which is the insertion of the ether C-O bond for subsequent cleavage to the metal center, proceeds via a four-membered transition state TSA2 involving the metal, the α-carbon, the ether oxygen and the ether carbon, with an activation energy barrier of 11.2 kcal mol −1 .The bond between the iridium and the ether carbon transitions from having no significant interaction (d Ir-C = 3.15 Å) in A2 to 2.34 Å in A3, with the O-Ir-C angle moving from 63.9 • in A2 to 92.2 • in TSA2 to 53.1 • in A3.Triplet state structures for TSA2 and A3 were found not to exist on the reaction surface.In the six-coordinate intermediate ( ipr PCP)Ir(H)(OPh)(enolate) A3, the enolate which was proposed to form a π-complex with the metal center [14] preferably binds to the metal center via the methylene carbon in an equatorial fashion while the phenolate sits axially to the pincer backbone.The two possible orders of product release were considered, and it was revealed that the thermodynamic product is the enolate, which rearranges to the more stable ketone while the kinetic product is the phenol.On the singlet surface, the formation of A4i mixed with the free enolate 6, which can rearrange to the more stable ketone 7, is calculated to be the major resting state in the catalytic cycle with a free energy of −25.6 kcal mol −1 relative to the reactants, whereas A4i undergoes reductive elimination to give phenol with an activation barrier of 26.8 kcal mol −1 ; this step is endergonic by 20.6 kcal mol −1 .In contrast, the release of phenol first proceeds with an activation barrier of 16.4 kcal mol −1 and the subsequent release of the catalyst is endergonic by 16.4 kcal mol −1 .An exhaustive search for the transition state from A3 and A4ii that leads to the bond cleavage of the enolate and the release of the catalyst was not fruitful.The energies for both orders of product release show that subsequent removal of either the phenol or the enolate, which rearranges to the more stable ketone, to give back the catalyst is kinetically unfavorable.This means that the likely predominant products are A4i+6+7 (the thermodynamic product) and A4ii+5 (the kinetic product).These findings are consistent with reports from experimental work carried out by Kundu et al.where the intermediate prior to the release of phenol was isolated instead of the phenol and the catalyst [14].On the triplet PES, the formation of the A4i is less stable than the singlet product by 83.5 kcal mol −1 while formation of A4ii is less stable by 118.1 kcal mol −1 .Although some of the minima (intermediates) and transition states in the catalytic cycle exist both in the singlet state and in the triplet state, the singlet state structures would deplete the triplet state reaction surface and no multiple-state reactivity is observed in the catalytic cycle.

Singlet State/Kcal mol
The optimized geometries of the stationary points, as well as the Gibbs free energy profile for the reaction of ( ipr PCP)Ir with 2-phenoxy-1-phenylethanol 1 along pathway B are shown in Figures 3 and 4. Triplet-state energies obtained along this pathway were also found to be less stable than the singlet-state energies (Table 2).The reaction begins with the ether oxygen coordinating axially to the metal center, which is exergonic by −6.9 kcal mol −1 , and a bond is formed between the Ir and O to give B2, which was not found to exist on the triplet surface.Subsequent addition of the β-carbon to Ir and cleavage of the ether C-O bond takes place with an activation energy barrier of 34.1 kcal mol −1 .
Here, the C-O bond insertion proceeds via a three-membered transition state TSB1, leading to the five-coordinate intermediate B3, which is exergonic by 10.5 kcal mol −1 .The Ir-O bond length transitions from a bond length of 2.32 Å in B2 to 2.23 Å in B3, while the phenolate remains axial to the pincer backbone and the enolate sits in an equatorial fashion.The Ir-C(H 2 ) bond also moves from having no significant interaction (d Ir-C = 3.32 Å) in B2 to a covalent bond (d Ir-C = 2.08 Å) in B3.The intermediate B3 was also located on the triplet PES, but was found to be less stable by 19.6 kcal mol −1 , while the transition state leading to B3 was found not to exist on the triplet surface.Subsequent H abstraction by Ir from the α-carbon proceeds via a four-membered early transition state TSB2 with an activation energybarrier of 22.5 kcal mol −1 , leading to a six-Ir-coordinated complex B4i.In the intermediate B4i, the enolate sits axially to the aromatic backbone of the pincer ligand, forming a π-complex with the iridium center and moving the phenolate and the abstracted hydrogen to the equatorial position.Neither TSB2 nor B4 could be located on the triplet potential energy surface.In order to determine the preferred order of product release, a different isomer of B4i which would allow for the release of the phenol first was located on the PES, but this isomer B4ii was found to be less stable by 12.4 kcal mol −1 than B4i.The corresponding activation barrier for the initial release of the phenol from the predominant isomer of B4i is 34.3 kcal mol −1 .Initial release of the enolate proceeds without barrier via the bond cleavage of the η 2 -Ir-C bond in a slightly endergonic fashion (∆G • = 1.0 kcal mol −1 ), followed by the release of phenol which proceeds with an activation barrier of 26.8 kcal mol −1 and is endergonic by 20.5 kcal mol −1 .This shows that the preferred order of product release is the initial release of the enolate followed by the release of the phenol.However, the energetics also show that the predominant product is B8i+6+7, as subsequent release of the phenol and the catalyst is shown to be kinetically less favored.Similar to pathway A, the formation of B8 and the transition state structure for the release of phenol are highly unstable on the triplet PES, and the major product formed is the enolate 6 which rearranges to the acetophenone 7 (essentially 7), as shown in Figure 4.No multiple state reactivity is observed along path B.
The optimized geometries of the stationary points, as well as the Gibbs free energy profile for the reaction of ( ipr PCP)Ir with 2-phenoxy-1-phenylethanol 1 along pathway C are shown in Figures 5 and 6.The reaction begins with the substrate binding to the catalyst via the ketone oxygen after undergoing dehydrogenative equilibrium in an exergonic fashion to give the adduct C2.C2 is located both on the singlet and triplet surface, with the triplet state structure being less stable by 12.9 kcal mol −1 (Table 3).The complex C2, which is proposed to undergo direct C-O bond cleavage, rather undergoes initial rearrangement to form a π-complex across the ketone C-O bond to the Ir center, allowing the ether oxygen of the substrate to sit in a much closer proximity to the metal center to allow for subsequent C-O addition.This rearrangement is thermoneutral and proceeds via TSC1, which could only be located on the triplet PES with an activation barrier of 31.5 kcal mol −1 .Although the activation barrier for the formation of C3 would be lower on the triplet surface, the more stable singlet-state C2 would deplete the surface of the triplet-state C2, making the formation of C3 proceed from singlet-state C2 via triplet-state TSC1 to singlet-state C3.The intermediate C3 is also located on the triplet surface, but is found to be less stable by 30.4 kcal mol −1 than the singlet-state C3.Subsequent C-O bond cleavage follows via TSC2, which involves the ether oxygen, the ether carbon, α-carbon and the Ir, with an activation barrier of 26.7 kcal mol −1 to give a five-membered Ir-complex C5, containing the phenoxide and what should have been a π-complex between the C-C bond and the Ir, but which is optimized as a C-bound ketone.The formation of C4 is exergonic by 21.4 kcal mol −1 on the singlet PES while the triplet state C4 is less stable by 20.4 kcal mol −1 .Coordination of initially released H 2 is exergonic by 3.2 kcal mol −1 leading to an η 2 -H 2 complex C5, shown to be the major resting state in the kinetically preferred catalytic cycle.In C5, which is only located to the singlet PES, the H 2 ligand binds axially to the aromatic backbone of the pincer ligand while both the phenoxide and ketone are arranged equatorial to it, allowing for reductive elimination of either the phenol or the acetophenone.The release of the acetophenone first proceeds via TSC3i with a barrier of 24.1 kcal mol −1 , followed by the release of the phenol via TSC4i with a barrier of 26.7 kcal mol −1 , while the barriers for release of the phenol first, followed by subsequent release of the acetophenone, proceeding via TSC3ii and TSC4ii, are 19.2 kcal mol −1 and 4.5 kcal mol −1 , respectively.The formation of C6i+7 from the initial release of the acetophenone is seen as the major resting state (∆G = −41.8kcal mol −1 ) in the reaction cycle, also showing that the acetophenone 7 would be the predominant product obtained, while the phenol 5 is the kinetic product.The release of the catalyst is also seen to be kinetically feasible if the pathway that follows the initial release of the phenol is considered.The η 2 -H 2 complex C5, or the transition states TS3i and TS3ii, could not be located on the triplet PES, while the intermediates C6i and C6ii, as well as the transition states TS4i and TS4ii that were found on the triplet PES, were highly unstable.A two-state reactivity is observed along this pathway (path C) with the transition state for the rearrangement of C2 to C3 located only in the triplet state.However, the other stationary points on the catalytic cycle, which exist both on the singlet surface as well as on the triplet surface, show the singlet state structures to be the ground state structures of the system.Tables 4-6 show comparisons in activation energies, energies of formation and apparent activation energies (δE) for the entire cycles calculated for the three pathways and for the two different ligand combinations.For pathway A, where the ether bond cleavage pro-ceeds via the initial C-H addition, the δE was calculated to be 15.2 kcal mol −1 , while the δE for pathway B and pathway C, where ether bond cleavage proceeds via direct C-O bond addition, were calculated to be 32.0 kcal mol −1 and 27.0 kcal mol −1 , respectively.Taking into account preferred product outcomes, the δE for pathway A reduces to 11.1 kcal mol −1 when the phenol is the desired product, while the δE for pathway C increases to 46.7 kcal mol −1 when the acetophenone is the desired product.From these apparent activation energies, the selectivity for the initial C-H addition pathway for this type of cleavage reaction, as observed in studies carried out by Choi et al. [14] and Haibach et al. [13], has been demonstrated.The activation barrier for the ether bond cleavage is also lower in pathway A (Ea = 16.4 kcal mol −1 ) than in pathways B (Ea = 34.1 kcal mol −1 ) and C (Ea = 26.7 kcal mol −1 ).For all three pathways, the phenol is observed as the kinetic product, while the enolate/acetophenone is observed as the thermodynamic product, and release of the catalyst becomes more kinetically feasible when phenol is the desired product.It should also be noted that for pathways B and C, in which the reaction proceeds via direct C-O bond addition, C-O bond cleavage in pathway C, which involves the presence of the ketone functionality of the substrate, is kinetically favored by 7.4 kcal mol −1 over pathway B. These results support both experimental [12] and theoretical [20] findings, which have reported that the presence of the ketone functionality facilitates the C-O bond cleavage and that the bond dissociation energy of the ketone is far less than its alcohol counterpart.

Comparison between the Two Ligands ipr PCP and ipr PCOP
The differences in δE between the ( ipr PCP)Ir-catalyzed reaction and the ( ipr PCOP)Ircatalyzed process are well within the margin of error for the three pathways (∆δE pathA = 1.5 kcal mol −1 , ∆δE pathB = 1.8 kcal mol −1 and ∆δE pathC = 3.8 kcal mol −1 ).Similar to the ( ipr PCP)Ir-catalyzed reaction, the observed trend in δE for the reaction along all three pathways is path B > path C > path A. Thus, there is also a kinetic preference for the pathway which begins with the C-H addition for the ( ipr PCOP)Ir-catalyzed β-O-4 compound in lignin.

The Reaction of ( ipr PCP)Co with 2-Phenoxy-1-phenylehtanol
Figures 7-12 show the optimized geometries as well as the Gibbs free energy profiles for the catalytic dehydroaryloxylation of 2-phenoxy-1-phenylehtanol 1 using ( ipr PCP)Co along pathways A, B and C, respectively.All stationary points have been sought on both the singlet and triplet potential energy surface as shown in Tables 7-9, with the triplet state reactants being 77.2 kcal mol −1 less stable than the singlet state reactants.Along pathway A, the ( ipr PCP)Co-catalyzed cleavage of 1 proceeds initially via C-H addition through TSA1 with an activation barrier of 13.9 kcal mol −1 to form A2, whereas C-O bond insertion follows via TSA2 with an activation barrier of 17.9 kcal mol −1 .The abstracted hydrogen sits axially to the aromatic backbone of the pincer ligand in TSA1 but shifts to an equatorial position during the C-O bond cleavage, while the cleaved phenolate replaces it in the axial position in TSA2.It is found that the six-membered intermediate is not formed after the C-O bond cleavage occurs, as the enolate is released in this step with a reaction energy of 27.9 kcal mol              In the ( ipr PCP)Co-catalyzed reaction of 1 along pathway B, the ether oxygen forms a covalent bond with the cobalt complex to give an adduct B2 in an exergonic reaction (∆G • = −7.3kcal mol −1 ).Direct ether cleavage via C-O bond addition follows with an activation energy barrier of 27.2 kcal mol −1 , leading to the O-bound and C-bound complex B3.Here, a two-state reactivity is observed, as the triplet state B3 is found to be more stable than the singlet state B3 by 9.8 kcal mol −1 .The ether-carbon transitions from a position axial to the aromatic backbone in TSB1 to an equatorial position in B3, while the etheroxygen sits axially to the aromatic backbone in both TSB1 and B3.Subsequent hydrogen abstraction via C-H addition to the metal center proceeds from the more stable triplet-state B3/t via singlet-state TSB2 with an activation energy barrier of 23.3 kcal mol −1 to form a singlet-state B4 which is endergonic by 5.1 kcal mol −1 .Similar to pathway A, cleavage of the enolate occurs in this step leaving only the phenolate and the hydride bound to the metal center, as seen in intermediate B4 which then undergoes reductive elimination of the phenol and regeneration of the catalyst with an activation energy barrier of 26.3 kcal mol −1 .The enolate 6, and the more stable acetophenone 7, to which it rearranges, are also shown to be the major products formed in this reaction.Apart from the intermediate B3, all other stationary points along the reaction cycle are found to be unstable on the triplet PES.
In the ( ipr PCP)Co-catalyzed reaction of 1 along pathway C, all but one of the stationary points in the catalytic cycle exhibit multiple spin ground states, where more than one of those stationary points with lower energies belong to the higher spin state and a two-state reactivity is observed as shown in Figure 12.The coordination of the ketone substrate to the metal center to form the adduct C2 is exergonic by −14.1 kcal mol −1 on the triplet surface and −9.9 kcal mol −1 on the singlet surface.This causes subsequent rearrangement of the adduct to form C3 containing the η 2 -bound ketone to proceed via TSC1 with an activation energy barrier of 15.7 kcal mol −1 instead of a lower barrier of 11.5 kcal mol −1 .Similar to what was observed in the ( ipr PCP)Ir-catalyzed reaction, the ketone-carbon transitions from having no significant interaction with the cobalt center in C2 (d Ir-C = 2.95 Å) to forming a covalent bond with the metal center in C3 (d Ir-C = 1.99 Å), and moves from occupying an axial position relative to the aromatic backbone of the pincer ligand in C2 to an equatorial position in C3.The bond distance between the ether-oxygen and the metal center also reduces from 4.37 Å to 2.18 Å after the rearrangement of C2 to C3, bringing the ether oxygen in close proximity for subsequent C-O bond cleavage.Ether bond cleavage via TSC2/t is kinetically more favored on the triplet surface with an activation energy barrier of 15.1 kcal mol −1 (Ea [singlet state] = 33.2kcal mol −1 ), leading to the formation of the triplet state O-bound and C-bound ketone C4/t, which is exergonic by 28.5 kcal mol −1 and more stable than the singlet-state C4 by 4.2 kcal mol −1 .The phenolate occupies the axial position relative to the aromatic backbone of the pincer ligand, while the ketone occupies the equatorial position.However, both the phenolate and the ketone move equatorially to the pincer ligand after hydrogenation occurs, putting the η 2 -H 2 in the axial position in C5.The hydrogenation step is endergonic by 18.2 kcal mol −1 on the triplet surface and endergonic by 3.7 kcal mol −1 on the singlet surface.Similar to the observed trend in the ( ipr PCP)Ir-catalyzed reaction, the phenol is the kinetic product while the acetophenone is the thermodynamic product.Reductive elimination of the phenol first proceeds with an activation barrier of 14.4 kcal mol −1 with the released product being formed at −12.4 kcal mol −1 , while initial release of the acetophenone proceeds with activation barrier of 15.7 kcal mol −1 and the released product is formed at −39.2 kcal mol −1 .The intermediate C5, and the transition states TSC3i and TSC3ii are all found to be more stable on the singlet surface.Subsequent release of the acetophenone from C6ii, which proceeds with an activation energy barrier of 3.1 kcal mol −1 on the singlet surface, and is lower than the triplet state activation energy barrier by 26.4 kcal mol −1 , leads to the release of the catalyst, which is exergonic by 9.4 kcal mol −1 , while subsequent release of the phenol from C6i proceeds with an activation energy barrier of 26.3 kcal mol −1 which is less kinetically favored for the release of the catalyst.Hence, the likely observed product from this sequence of product release would be C6i+7.show the comparison in activation energies, energies of formation and apparent activation energies (δE) for the entire cycles, calculated for the three pathways and two different ligand combinations.The trend observed here in the activation energies and apparent activation energies (δE) for the entire cycle is different from the ( ipr PCP)Ir-catalyzed process.The activation barrier for the C-O bond addition is seen to be lower in pathway C (Ea = 15.1 kcal mol −1 ) than pathways A (Ea = 17.9 kcal mol −1 ) and B (Ea = 27.2 kcal mol −1 ), whereas, ultimately, the calculated apparent activation energies (δE) for pathway C are shown to be kinetically favored by 15.7 kcal mol −1 over pathway A and by 17.3 kcal mol −1 over pathway B. These energetics show that the kinetically favored pathway for the ( ipr PCP)Co-catalyzed cleavage of a β-O-4 compound is pathway C, which involves the ketone motif of the substrate and begins with the C-O bond addition, even when the desired product is acetophenone (δE = 16.7 kcal mol −1 ).The differences in the apparent activation energy (δE) between the ( ipr PCP)Cocatalyzed process and the ( ipr PCOP)Co-catalyzed process are also well within the margin of error for all three pathways (∆δE pathA = 2.1 kcal mol −1 , ∆δE pathB = 3.0 kcal mol −1 , ∆δE pathC = 5.5 kcal mol −1 ).When we compare the δE values for the ( ipr PCOP)Co-catalyzed reaction among all three pathways, the observed trend is path B > path A > path C.These values also show that pathway C is the kinetically favored pathway for the ( ipr PCOP)Cocatalyzed cleavage of the β-O-4 linkage in lignin, and that the PCOP ligand displays comparatively similar catalytic activity to the PCP ligand.

Iridium vs. Cobalt
Along pathway A, the ( ipr PCP)Ir-catalyzed process is kinetically favored by 8.8 kcal mol −1 over the ( ipr PCP)Co-catalyzed reaction, while along pathways B and C, the Co-catalyzed reaction is kinetically favored over their Ir-catalyzed counterpart by 6.4 kcal mol −1 and 18.7 kcal mol −1 , respectively.From the obtained data, the use of the ( ipr PCP)Co catalyst for the cleavage of alkyl aryl ethers is a comparable alternative to the ( ipr PCP)Ir catalyst, with the Co being more active, although the ( ipr PCP)Co-catalyzed process reaction would preferably proceed via direct C-O insertion, while the ( ipr PCP)Ir-catalyzed process would rather proceed via direct C-H addition.

Calculation of Rate Constants for the Cleavage Reactions
The apparent activation energies are used to calculate the rate constants of the reactions using the equation derived from Erying's transition state theory [21], where ∆G • in the equation is replaced by δE.
The rate constant k(T) is also defined as the turnover frequency (TOF) by Kozuch and Shaik [22] based on the energetic span approximation.Taking the temperature as 298.15K, the computed rate constants/TOF are shown in Table 13 below.It is apparent from these rate constant values that the fastest reaction is the cleavage of the β-O-4 linkage with the ( ipr PCP)Co catalyst via the pathway which involves the ketone motif of the substrate and proceeds via direct C-O bond insertion (Path C), followed by the use of the ( ipr PCP)Ir catalyst via the direct C-H bond addition mechanism (Path A).The fact that more intermediates and transition states occur in the triplet state along pathway C for the cleavage reaction using ( ipr PCP)Co greatly affects the kinetics of the reaction, causing it to proceed much faster.

Computational Details
All calculations were carried out using the Spartan 14 Molecular Modelling program [23] and the Gaussian 09 package [24].The structures and energies of all the stationary points along the reaction pathways were computed using the M06 functional.The reference energy and the energy differences for each stationary point in the "free energy landscape" figures correspond to the reactants; these are the Ir/Co catalyst and 2-phenoxy-1-phenylethanol (I).The atoms hydrogen and argon were described with the 6-311G** basis set [25], while the metals Ir and Co were modeled with the quasi-relativistic pseudopotential of triple-ζ quality, LANL2TZ [26,27].The M06 functional is a hybrid metagradient-corrected functional (meta-GGA) with 27% of Hartree-Fock exchange, which has been benchmarked and is recommended for the study of organometallic thermochemistry as well as non-covalent interactions [28,29].All stationary points were verified by full frequency calculations.Minima (reactants, intermediates and products) were shown to have a Hessian matrix whose eigenvalues are all positive, leading to vibrational frequencies that are real, while transition states were shown to have a Hessian matrix with all positive eigenvalues, except for a single negative eigenvalue characterized by a vibration along the reaction coordinate.
The starting geometries of the molecular systems were constructed using Spartan's graphical model builder and minimized interactively using the sybyl force field.Input structures for transition state (TS) optimizations were obtained after a potential energy surface (PES) scan was performed for each reaction step.This calculation was also useful in showing that the minima obtained do in fact correspond to the reactant and product for that step.The PES scan gives an approximate transition state structure, which is then submitted for transition state calculations using the synchronous Transit-Guided Quasi-Newton (STQN) Method developed by Schlegel and coworkers [30].General effects of a surrounding solvent were evaluated using a polarizable continuum model derived with an integral equation formalism (IEF-PCM) [31] of p-xylene solvation, which is the solvent used in experiments.However, the calculations that included solvent effects showed differences of just <2 kcal/mol in the calculated barrier heights with gas phase calculations, and hence only gas phase calculations are reported below.
The apparent activatis.nenergy, also known as the rate-limiting activation energy, of the catalytic cycle is obtained from the energetic span model as defined by Kozuch and Shaik [22], which is an extension of the Curtin-Hammett principle.Two terms appear in this model that determine the turnover frequency (TOF) of a catalytic cycle: the TDTS (TOF-determining transition state) and the TDI (TOF-determining intermediate).These species define the energetic span within the cyclic constraints and in turn measure the kinetics of the catalytic cycle.The energetic span, δE, is thus defined as δE = T TDTS − T TDI if TDTS appears after TDI T TDTS − T TDI + ∆G r if TDTS appears before TDI (2) To ascertain the suitability of the M06 functional to give reliable energetics, some of the steps of the cycle were computed at the M06-2X-D3 level for comparison.The M06-2X-D3 optimized geometries are very similar to those of the M06 geometries, and the differences in activation barriers range from 0 to 5 kcal mol −1 (Supporting Information, Tables S1-S3).

Conclusions
We have reported the kinetics and thermodynamics of the cleavage of a model of the β-O-4 linkage in lignin (2-phenoxy-1-phenylethanol) using iridium and cobalt pincer complexes.The observed selectivity for the initial C-H addition mechanism (Path A) for the cleavage reaction in experiment [11,13] is confirmed, as the initial C-H addition pathway is kinetically favored by 16.8 kcal mol −1 over the direct C-O insertion pathway in the ( ipr PCP)Ir-catalyzed reaction.
The cobalt catalysts show a kinetic preference for the mechanism that involves an initial dehydrogenative equilibrium of the substrate with its ketone motif, and then proceeds via direct C-O insertion (Path C).This pathway is found to be kinetically favored by 15.7 kcal mol −1 over the initial C-H addition pathway.
The cleavage reaction is kinetically controlled, having a rate constant of 44.7 s −1 when the ( ipr PCP)Ir catalyst is used and a rate constant of 5.1 × 10 6 s −1 when the ( ipr PCP)Co is used.The cobalt-pincer catalyst has been shown to be more active for the cleavage reaction, having a two-state reactivity along the catalytic cycle, which greatly reduces its activation barrier.
It is observed from the calculated energies that the kinetic product is the phenol, while the acetophenone is the thermodynamic product for both the iridium-and cobalt-catalyzed reactions.The catalyst can be recovered easily when cobalt is used and when the phenol is the desired product.Some of the intermediates and the transition states located in this work display multiple-spin ground states and a two-state reactivity is observed in the ( ipr PCP)Cocatalyzed reaction.
In summary, this study has shown that PCP-type iridium and cobalt catalysts are active towards the cleavage of alkyl aryl ethers, such as the β-O-4 linkage in lignin, and it has provided molecular level insights into the preferred mechanism for the cleavage reactions.This new knowledge will aid in the design of an atom-economic catalyst for the cleavage of alkyl aryl ethers, which are widely used in chemistry, as well as the valorization of lignin towards the production of value-added aromatic compounds.

Scheme 2 . 2 .
Scheme 2. Dehydroaryloxylation of aryl ethers using PCP-type iridium complexes.Study carried out by Haibach et al.However, it is still not understood from a mechanistic point of view how these bonds are cleaved.In studies carried out by Choi et al., both DFT calculations and experimental work have shown that the C-O bond cleavage proceeded via initial C-H bond activation Scheme 2. Dehydroaryloxylation of aryl ethers using PCP-type iridium complexes.Study carried out by Haibach et al.However, it is still not understood from a mechanistic point of view how these bonds are cleaved.In studies carried out by Choi et al., both DFT calculations and experimental work have shown that the C-O bond cleavage proceeded via initial C-H bond activation rather than the direct C-O insertion [10].Aryl C-H insertion was also observed to be kinetically favored in the reaction of RuH 2 CO(PPH 3 ) with aryl ethers [15].However, DFT studies backed by experimental findings reported by Chmely et al. showed the kinetic preference of the direct C-O activation over C-H activation in the cleavage of an alkyl aryl ether using a Ru-xantphos catalyst [16].Liu and Wilson also reported a DFT study on the reaction between PCP-type group 8 complexes with a model of the β-O-4 linkage in lignin, where the reaction proceeds via direct C-O bond cleavage [17].Consequently, in the effort

Scheme 3 .
Scheme 3. Chemical equation for the catalytic cleavage of the β-O-4 linkage, 2-phenoxy-1-phenylethanol using PCP-type iridium complexes and the chemical structures of the catalysts used.Scheme 4 shows the mechanistic pathway as proposed by Haibach et al. (route A).The initial step is the addition of the α-C-H bond to the metal center to yield A2, followed by the C-O insertion step to form the six-coordinate Ir intermediate A3.Loss of the phenol 5 and enolate 6, which rearranges to the more stable acetophenone 7, regenerates the Irpincer active species and completes the catalytic cycle.

Scheme 3 .
Scheme 3. Chemical equation for the catalytic cleavage of the β-O-4 linkage, 2-phenoxy-1-phenylethanol using PCP-type iridium complexes and the chemical structures of the catalysts used.Scheme 4 shows the mechanistic pathway as proposed by Haibach et al. (route A).The initial step is the addition of the α-C-H bond to the metal center to yield A2, followed by the C-O insertion step to form the six-coordinate Ir intermediate A3.Loss of the phenol 5 and enolate 6, which rearranges to the more stable acetophenone 7, regenerates the Ir-pincer active species and completes the catalytic cycle.The second mechanism as proposed by Liu and Wilson[17] begins with the coordination of the ether oxygen to the metal center to form an adduct.Direct C-O bond addition follows to form the five-coordinate intermediate B3, followed by the α-C-H bond addition to form B4. Finally, the phenol and enolate, which rearrange to acetophenone, are released to regenerate the Ir-pincer active species, as shown in Scheme 5.The third proposed reaction pathway, shown in Scheme 6, is similar to the mechanism proposed by Wu et al.[19] in their experimental study of the cleavage of a β-O-4 model compound using a Ru-xantphos catalyst.It begins with dehydrogenation of the substrate 2-phenoxy-1-phenyl-ethanol 1 to its corresponding ketone motif C2.C2 undergoes rearrangement via TSC1 to a slightly more stable adduct with a three-membered ring formed between the ketone carbon, the ketone oxygen and the Ir center.This is followed by the direct C-O bond insertion to form C4, which undergoes hydrogenation using the H 2 generated earlier.The H 2 molecule, which binds to the metal to form a dihydrogen complex

Figure 1 .
Figure 1.Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Ir with 2-phenoxy-1-phenylethanol 1 along pathway A. Distances in Å.

Figure 1 .
Figure 1.Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Ir with 2-phenoxy-1-phenylethanol 1 along pathway A. Distances in Å.

Figure 2 .
Figure 2. Free energy landscape of the overall catalytic cycle for the reaction of ( ipr PCP)Ir with 1 along pathway A.

Table 1 .
Comparison of singlet-state and triplet-state energies of the stationary points involved in the reaction of ( ipr PCP)Ir with 1 along path A.

Figure 2 .
Figure 2. Free energy landscape of the overall catalytic cycle for the reaction of ( ipr PCP)Ir with 1 along pathway A.

Catalysts 2023 , 24 Figure 3 .
Figure 3. Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Ir with 2-phenoxy-1phenylethanol 1 along pathway B. Distances in Å.

Figure 4 .
Figure 4. Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Ir with 1 along pathway B.

Figure 4 .
Figure 4. Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Ir with 1 along pathway B.

Figure 4 .
Figure 4. Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Ir with 1 along pathway B.

Figure 5 .
Figure 5. Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Ir with 2-phenoxy-1-phenylethanol 1 along pathway C. Distances in Å.

Figure 5 .
Figure 5. Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Ir with 2-phenoxy-1-phenylethanol 1 along pathway C. Distances in Å.

Figure 6 .Table 3 .
Figure 6.Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Ir with 1 along pathway C. Triplet state energies shown in blue.

Figure 6 .
Figure 6.Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Ir with 1 along pathway C. Triplet state energies shown in blue.

Figure 7 .
Figure 7. Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Co with 2-ph noxy-1-phenylethanol 1 along pathway A. Distances in Å.

Figure 7 .
Figure 7. Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Co with 2-phenoxy-1-phenylethanol 1 along pathway A. Distances in Å.

Figure 7 .
Figure 7. Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Co with 2-phenoxy-1-phenylethanol 1 along pathway A. Distances in Å.

Figure 8 .
Figure 8. Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Co with 1 along pathway A. Triplet state energies shown in blue.Figure 8. Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Co with 1 along pathway A.

Figure 8 . 24 Figure 9 .
Figure 8. Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Co with 1 along pathway A. Triplet state energies shown in blue.Figure 8. Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Co with 1 along pathway A. Catalysts 2023, 13, x FOR PEER REVIEW 16 of 24

Figure 9 .
Figure 9. Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Co with 2-phenoxy-1-phenylethanol 1 along pathway B. Distances in Å.

Figure 9 .
Figure 9. Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Co with 2-phenoxy-1-phenylethanol 1 along pathway B. Distances in Å.

Figure 10 .
Figure 10.Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Co with 1 along pathway B. Triplet state energies shown in blue.Figure 10.Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Co with 1 along pathway B. Triplet state energies shown in blue.

Figure 10 . 24 Figure 11 .
Figure 10.Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Co with 1 along pathway B. Triplet state energies shown in blue.Figure 10.Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Co with 1 along pathway B. Triplet state energies shown in blue.Catalysts 2023, 13, x FOR PEER REVIEW 17 of 24

Figure 11 .
Figure 11.Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Co with 2-phenoxy-1-phenylethanol 1 along pathway C. Distances in Å.Figure 11.Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Co with 2-phenoxy-1-phenylethanol 1 along pathway C. Distances in Å.

Figure 11 .
Figure 11.Optimized geometries of the main stationary points involved in the reaction of ( ipr PCP)Co with 2-phenoxy-1-phenylethanol 1 along pathway C. Distances in Å.

Figure 12 .Table 7 . 1 Figure 12 .
Figure 12.Free energy landscape for the overall catalytic cycle for the reaction of ( ipr PCP)Co with 1 along pathway C. Triplet state energies shown in blue.Table 7.Comparison of singlet-state and triplet-state energies of the stationary points involved in the reaction of ( ipr PCP)Co with 1 along pathway A. Singlet State/Kcal mol −1 Triplet State/Kcal mol −1

Table 2 .
Comparison of singlet-state and triplet-state energies of the stationary points involved in the reaction of ( ipr PCP)Ir with 1 along path B.

Table 3 .
Comparison of singlet-state and triplet-state energies of the stationary points involved in the reaction of ( ipr PCP)Ir with 1 along path C.

Table 4 .
Activation energies and reaction energies for the reaction of 1 with [IrPCP] and [IrPOCP] along path A. Energies in kcal mol −1 .

Table 5 .
Activation energies and reaction energies for the reaction of 1 with [IrPCP] and [IrPOCP] along path B. Energies in kcal mol −1 .

Table 6 .
Activation energies and reaction energies for the reaction of 1 with [IrPCP] and [IrPOCP] along path C. Energies in kcal mol −1 .

Table 7 .
Comparison of singlet-state and triplet-state energies of the stationary points involved in the reaction of ipr PCP)Co with 1 along pathway A.Singlet State/Kcal mol−1Triplet State/Kcal mol−1

Table 8 .
Comparison of singlet-state and triplet-state energies of the stationary points involved in the reaction of ( ipr PCP)Co with 1 along pathway B.

Table 9 .
Comparison of singlet-state and triplet-state energies of the stationary points involved in the reaction of ( ipr PCP)Co with 1 along pathway C.

Table 10 .
Activation energies and reaction energies for the reaction of 1

Table 11 .
Activation energies and reaction energies for the reaction of 1 with [CoPCP] and [CoPOCP] along pathway B. Energies in kcal mol −1 .

Table 12 .
Activation energies and reaction energies for the reaction of 1 with [CoPCP] and [CoPOCP] along pathway C. Energies in kcal mol −1 .

Table 13 .
Rate constants for the ( ipr