The Tricarbonylhydridocobalt-Based Hydroformylation Reaction

process results in the formation of the catalytically active species tricarbon- ..... as phosphines and retarded by IT acceptors such as carbon monoxi...
0 downloads 0 Views 2MB Size
5

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

The Tricarbonylhydridocobalt-Based Hydroformylation Reaction A Theoretical Study

Tom Ziegler and Louis Versluis Department of Chemistry, University of Calgary, Calgary, Alberta T2N 1N4, Canada

Density functional calculations have been carried out on the mechanism proposed by Heck and Breslow for the hydroformylation process based on tricarbonylhydridocobalt [HCo(CO) ]. Geometries and relative energies were determined for intermediates involved in each elementary step. Reaction profiles were further traced by an approximate linear transit procedure. 3

THE

O X O O R H Y D R O F O R M Y L A T I O N R E A C T I O N , discovered in 1938

by Roe-

len, is used on a large industrial scale (1-3) to convert olefins and synthesis gas into aldehydes. The process employs homogeneous catalysts based on cobalt (1-3) or rhodium (4). The most commonly used precatalyst is H C o ( C O ) , which is generated in situ from the hydrogénation of C o ( C O ) by H . A mechanism for the cobalt-based hydroformylation process was first proposed by Heck and Breslow (5) in 1961 (Scheme I). The catalytic cycle in Scheme I consists of a number of elementary reaction steps (a-e) that we will discuss. The emphasis of the investigation lies on the identification of the equilibrium geometries and the relative energies of stable intermediates in Scheme I. The energy profile of the reaction paths connecting the intermediates will be modeled by an approximate linear transit procedure (6, 7). 4

2

2

0065-2393/92/0230-0()75$06.00/0 © 1992 American Chemical Society

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

8

76

H O M O G E N E O U S TRANSITION M E T A L C A T A L Y Z E D REACTIONS

HCo(CO)

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

4

Scheme I.

Computational Details The calculations were all based on the Hartree-Fock-Slater (HFS) model as implemented by Baerends and co-workers (8, 9). With it we used the latest version of the fully vectorized H F S - L C A O - S T O program developed by Ravenek (10). The bonding energies were calculated by the generalized transition-state method (II, 12) in conjunction with Becke's (13) nonlocal exchange corrections, as well as corrections to allow for correlations between electrons of different spins (14). The numerical integration scheme employed in this work was formulated by Becke (15). A n uncontracted triple-ζ Slater-type orbitals (STO) basis set (16) was used for the 3s, 3p, 3d, 4s, and 4p shells of cobalt. The 2s and 2p shells of C and Ο and the Is shell of Η were described by a double-ζ S T O basis set (16, 17), which was extended with one polarization function (2p on Η and 3d on C and O). The electrons in the lower energy shells on C o , C , and Ο were considered as core electrons and treated by the frozen-core approxi­ mation method according to Baerends and co-workers (8, 9). All molecular structures were optimized within the C symmetry group. The geometry optimizations were carried out according to the algorithm developed by Versluis and Ziegler (18). s

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

5.

ZiEGLER & VERSLUIS

The Hydroformylation

77

Reaction

Application of Density Functional Methods to Organometallic Substances Calculations on metal carbonyls (19), binuclear metal complexes (20), alkyl and hydride complexes (21-23), and complexes containing metal-ligand bonds for a number of different ligands (24) have shown that the approximate density functional method employed here affords metal-ligand and metal-metal bond energies of nearly chemical accuracy ( ± 5 kcal mol" ). Approximate density functional methods have also been tested in connection with vibrational frequencies (25), conformational energies (6, 7), trip­ let-singlet separations (26), and transition-state structures (27). More than 50 molecular structures optimized by approximate density functional theory have been compared with experiment (18). The agreement between exper­ iment and approximate density functional theory is excellent in most cases. The present method has also been applied to a study of C - H activation by late transition metals (28), as well as organosilane polymerization (29) and halogen abstraction by metal carbonyls (30).

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

1

Dissociation of CO from HCo(CO) HCo(CO)

4

To Form the Catalyst

3

The initial key step (a of Scheme I) in the hydroformylation process is rep­ resented by the dissociation of a C O ligand from H C o ( C O ) . The dissociation process results in the formation of the catalytically active species tricarbonylhydridocobalt [HCo(CO) ] (eq 1). 4

3

HCo(CO)

4

HCo(CO)

3

+ CO

(1)

The dissociation is assumed to take place prior to, or in concert with, the complexation of an olefin leading to the olefinic ττ complex H C o ( C O ) ( η olefin) (step b of Scheme I). The coordinatively unsaturated 16-electron species, H C o ( C O ) , has been identified by matrix isolation techniques (31), but its structure is unknown. We explored (6, 7) possible structures for H C o ( C O ) . O u r investigation was confined to singlet states of H C o ( C O ) because the carbonyl dissociation process of eq 1 probably takes place on the singlet surface. The precatalyst H C o ( C O ) has (6, 7) a trigonal bipyramidal (TBP) groundstate conformation with the hydride in an axial position. Dissociation of either an equatorial or axial C O ligand from H C o ( C O ) will thus give rise to l a and l b , respectively. 3

2

3

3

3

4

4

In our calculations both l a and l b constitute local minima on the singlet surface of H C o ( C O ) . We calculate H C o ( C O ) with the butterfly shape (la) 3

3

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

78

H O M O G E N E O U S TRANSITION M E T A L C A T A L Y Z E D REACTIONS

1b

to be 38 kj mol lower in energy than the trigonal-shaped H C o ( C O ) species, l b . The equatorial ( H C o ( C O ) - * la) and axial [ H C o ( C O ) lb] C O dissociation energies were found to be 186 and 224 kj mol" , respectively. The C O dissociation energies for the d complex H C o ( C O ) are not known experimentally. However, our calculated values compare well with an experimental (32) (singlet) dissociation energy of 183 kj mol" in the d complex F e ( C O ) . The energy required to convert the precatalyst H C o ( C O ) into the active coordinatively unsaturated 16-electron species H C o ( C O ) of structure l a (186 kj m o l ) is substantial. We shall show that step a is the most energetically demanding of the steps in Scheme I. 1

3

4

4

1

8

4

1

8

5

4

3

1

Some experimental evidence for the existence of H C o ( C O ) in different configurations was given by Sweany and Russell (33, 34), who inferred on the basis of results from the photolysis of H C o ( C O ) in an argon matrix that H C o ( C O ) forms two isomers consistent with structures l a and l b . 3

4

3

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

5.

Z I E G L E R & VERSLUIS

The Hydroformylation

79

Reaction

Olefin Insertion into the Co-Η Bond The active catalyst H C o ( C O ) of conformation l a combines in step b of Scheme I with olefin to generate a ir-olefin complex where C H is coor­ dinated to the vacant equatorial site and the C - C bond is placed either in the equatorial plane (2a) or parallel to the C o - Η bond (2b). Conformation 2a was, as one might expect (35), calculated (7, 36) to be more stable than 2b. However, the difference is only 20 kj mol , and the C o - C H bond energy in 2a was estimated to be 70 kj mol" . 3

2

l

4

2

4

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

1

According to step c of Scheme I, the olefin will undergo a migratory insertion into the C o - H bond after its complexation, and thus form an ethyl complex. HCo(CO) en -C H ) - » Co(CO) C H 3

2

2

4

3

2

5

(2)

Only 2b has the proper relative orientation of ethylene and the hydride for the insertion. Conformation 2a must, as a consequence, rearrange to 2b before the process in eq 2 can take place. The reaction profile for the insertion is shown in Figure 1. The insertion process 2b —» 3a is exothermic by 8 kj mol and has a small activation barrier of 5 kj m o l . The calculated exothermicity and modest activation for the process is in agreement with the experimental observation. Thus, the migration of a hydride to a coordinated olefin group is observed experimentally to be very facile (37). In fact, the hydride-olefin insertion reaction has, with a few exceptions (38), rarely been directly observed. As a consequence, metal complexes containing both hy­ dride and olefin are scarce. 1

1

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

80

H O M O G E N E O U S TRANSITION M E T A L C A T A L Y Z E D REACTIONS

Reaction Coordinate Figure

1. Energy

profile

of the hydride

HCo(CO)^f\ -C H ). 2

2

migration

to the ethylene

The energy zero refers to structure

4

group

in

3a.

ο

3a A hydrogen atom bound to carbon can interact weakly with a metal center. We shall in the following refer to such an interaction as agostic. The optimized structure, 3a, for the resulting ethyl complex exhibits a clear agostic interaction between a β-hydrogen and the vacant metal center. However, under catalytic conditions, with P = 200-300 atm, the coordinatively unsaturated complex 3a will coordinate a C O ligand to form the saturated complex, C H C o ( C O ) , with the ethyl group in the equatorial position (3b). The coordinatively saturated ethyl complex of conformation 3b can subsequently rearrange to the more stable (6) conformation 3c, in which the ethyl group is in the axial position, by a Berry pseudorotation for which the activation energy is predicted (39) to be low. ( : o

2

5

4

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

5.

ZiECLER & VERSLUIS

The Hydroformylation Reaction

3b

81

3c

We might conclude that the olefin insertion into the C o - Η bond (steps b and c of Scheme I) is very facile. Thus the olefin insertion should not constitute a bottleneck in the hydroformylation process.

Migratory Insertion of Alkyl into the Co-CO Bond We modeled (6, 7) the migratory insertion process (step d of Scheme I) with C H C o ( C O ) (4) rather than C H C o ( C O ) . 3

4

2

5

4

RCo(CO) -> RC(0)Co(CO) 4

3

101

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

(3)

82

H O M O G E N E O U S TRANSITION M E T A L C A T A L Y Z E D

REACTIONS

The process in eq 3 could in principle proceed by an insertion of a cis-CO into the C o - C H bond. This insertion would produce the coordinatively unsaturated complex, 5a, with the acyl group in an axial position. Alternatively, the methyl group might migrate to a cis-earbonyl and thus form complex 5b with the acyl group in an equatorial position. Perhaps not surprisingly, we find that the energy profile for the C O insertion, 4 —> 5a, into

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

3

5a

5b

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

5.

Z I E G L E R & VERSLUIS

The Hydroformyfotion

83

Reaction

the C 0 - C H 3 bond (Figure 2a) has a prohibitively high activation barrier of 200 kj mol" . C O insertion, as a consequence, cannot be a viable mechanism for the process in eq 3. 1

The migration of C H to the cis~CO ligand, 4 —> 5b, was calculated (Figure 2b) to have an endothermicity, Δ Η , of 71 kj mol" and a very modest activation barrier, Δ £ * , of only 9 kj mol" . Thus the C H migration, 4 —» 5b, seems to be favored as the mechanism for the process in eq 3. The calculated reaction enthalpy and activation barrier for 4 —» 5b compare well with an earlier study (40) on the C H - » C O migration in C H M n ( C O ) . In that study we found ΔΗ to be 75 kj mol" and Δ £ * to be 11 kj mol" . O u r findings are also in agreement with a recent kinetic study by Roe (41), who found the rate constant of the methyl back migration of C H C ( 0 ) C o ( C O ) to be considerably larger than the rate constant for the corresponding forward reaction. The structures in Figure 2a illustrate how the methyl group can slide almost parallel along the cis C - C o bond onto the cis carbonyl carbon while the remaining C o ( C O ) framework stays almost unchanged. The 9 kj mol calculated for Δ £ * in the present study is an upper bound (6, 7) to the actual value. We can thus conclude that the methyl migration, 4 —» 5b, should proceed with a rather modest activation barrier. 3

1

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

1

3

3

3

1

5

1

3

3

3

1

The 1,2 shift reaction of an alkyl group in which a metal-alky 1 system is converted into a metal-acyl complex, 4 —> 5b, is well documented for a variety of alkyl complexes. The corresponding 1,2 shift reaction, 6a —» 6b, involving Η rather than alkyl, has proven to be rather elusive. The 1,2 hydride shift reaction was inferred in earlier work (42-44) as an elementary reaction step. In spite of considerable efforts it has been detected with certainty only in a few cases (45-47). It is now widely accepted that the hydride migration, in contrast to the alkyl migration, is thermodynamically unfavorable, at least for middle to late transition metals. We studied (6, 7) the 1,2 shift reaction 6a —> 6b, which represents a 1,2 shift of a hydride in H C o ( C O ) with C symmetry. The formyl structures 6b do not represent a local energy minima on the H F S energy surface. Thus, any attempt to optimize 6b resulted in a back migration of the formyl hy­ drogen to the parent hydrido metal complexes 6a. Our findings indicate that the formyl complex 6b is kinetically unstable with respect to the parent hydrido complex 6a. That is, the decarbonylation reaction (6b —» 6a) should have at most a minimal activation barrier. Our findings can be reconciled with the experimental observation that most neutral metal formyl complexes decarbonylate readily to the corresponding hydrido complexes (48-51). The decomposition is believed to occur by a back migration of the formyl hy­ drogen to the metal center under loss of a ligand of the coordinatively saturated formyl complex. We calculate (6, 7) the process in eq 4 to have an exothermicity of AH = -69 kj mol" . 4

3 v

1

H C ( 0 ) C o ( C O ) -> H - C o ( C O ) 4

4

+ CO

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

(4)

84

H O M O G E N E O U S TRANSITION M E T A L C A T A L Y Z E D REACTION

J

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

260

Reaction Coordinate Figure 2. a: Energy profile for the insertion of CO into the C0-CH3 bond, 4 —• 5a. b: Energy profile for the migration of CH to C O , 4 -> 5b. The zero point refers to 4 in both plots. 3

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

5.

Z I E G L E R & VERSLUIS

The Hydroformylation Reaction

H

H

CO

OC-Ci

CO

C ο

ο

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

6a

85

CO CO

C ο

6b

Acyl complexes are, in contrast to their formyl counterparts, well-known molecules. We calculate the corresponding decarbonylation process of C H C ( 0 ) C o ( C O ) to be endothermic with Δ Η - 20 kj mol" . Thus the acyl complex is thermodynamically stable, although the formyl analog is unstable and decomposes to H C o ( C O ) according to eq 4. The higher exothermicity of the decarbonylation process for the formyl compared to the acyl can largely be ascribed to the higher bond strength (21, 22, 52) of C o - Η compared to C o - C H . For middle to late transition metals D ( M - H ) is - 2 4 0 kj m o l whereas the corresponding D ( M - C H ) bond strength is only ~160 kj m o l . The isoelectronic d formyl complex HC(0)Fe(CO) ~ has been found (53) to decompose slowly to H F e ( C O ) " . The kinetic stability can probably be as­ cribed to the stabilization of the M - C O bond through increased back-bond­ ing interactions of the carbonyl ligands in the charged species. However, the overall reaction is thermodynamically favorable with a reported exo­ thermicity of Δ Η = -43 ± 30 kj mol" (54). 3

1

4

4

1

3

-1

3

8

4

4

1

H 'lnduced Aldehyde Elimination 2

The last step in the catalytic cycle of the hydroformylation process, e of Scheme I, is the reaction of the acyl intermediate with H . This reaction results in the formation of the desired aldehyde molecule and the regen­ eration of the catalyst H C o ( C O ) . 2

3

The aldehyde product can be formed from the acyl intermediate by several possible routes. Heck and Breslow (5) proposed a mechanism in which the coordinatively unsaturated acyl complex undergoes first an oxi­ dative addition of H to afford a dihydro acyl compound (eq 5a) followed by an irreversible reductive elimination of an aldehyde molecule (eq 5b). 2

R(0)CCo(CO)

3

+ H -> R(0)CCo(CO) (H) 2

R(0)CCo(CO) (H) 3

2

3

(5a)

2

R C H O + HCo(CO)

3

(5b)

This type of process has been inferred for numerous catalytic and stoichio­ metric systems (37). As a consequence, these reactions have been studied

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

86

H O M O G E N E O U S TRANSITION M E T A L C A T A L Y Z E D REACTIONS

extensively by both experimental (55-65) and theoretical (66-75) techniques. However, in the cobalt-based hydroformylation process it has not been es­ tablished that the product formation proceeds in fact via the oxidative ad­ dition-reductive elimination mechanism. Some experimental observations indicate (76-79) that the acyl complex might react with H C o ( C O ) , and thereby form an aldehyde molecule and a binuclear cobalt compound (eq 6a).

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

4

CO

R(0)CCo(CO)

3

+ HCo(CO)

Co (CO) 2

8

> R C H O + Co (CO)

4

+ H

2

2

> 2HCo(CO)

8

4

(6a) (6b)

In a subsequent reaction C o ( C O ) is then (according to eq 6b) transformed back to the mononuclear hydrido-cobalt complex. Experimental studies (5) of the stoichiometric reactions (eqs 6a and 6b) revealed that the process is very facile. However, under catalytic conditions the overall concentration of cobalt species is low in comparison to the reactants. That is, the probability for reaction between two cobalt complexes is small (77). For electron-poor systems such as early transition metal complexes (80-83), the hydrogenolysis of the M - C bond proceeds via a mechanism in which an incoming H molecule initially forms a η adduct with the metal complex. This reaction is followed by the concerted cleavage of the hydrogen bond and the formation of H - M and H - C bonds by way of a four-center intermediary structure such as that illustrated by structure 7. 2

8

2

2

Η

Η

/

\/°

M

C

7

This reaction mode thus omits the oxidative addition-reductive elimi­ nation mechanism. Such a process could also be envisioned for the cobalt system studied here and was speculated upon by earlier workers (2, 48). In this section we will concentrate our efforts on the study of those cobalt complexes that can result from an interaction of the acyl intermediates with an incoming H molecule. Furthermore, attention will be given to the ox­ idative addition process of H to the metal fragment as well as to the hy­ drogenolysis reaction through a four-center intermediary structure. 2

2

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

5.

Z I E G L E R & VERSLUIS

The Hydroformylation

87

Reaction

The two most stable products from the interaction between an acyl complex and H are the dihydrogen complexes 8a and 8b with H in the 2

2

equatorial position. Both represent energy minima on the H F S energy sur­ face; 8b is about 19 kj mol" above 8a in energy.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

1

8a

8b

The mechanism by Heck and Breslow (5) suggests that H adds oxidatively to the unsaturated acyl complex to form a dihydride, 8c. The dihydride complex of configuration 8c was found to be 25 kj m o l higher in energy than the η adduct 8a. Thus the H complex is more stable than the product of the oxidative addition. Complex 8c is also 6 kj mol" higher in energy than the η compound 8b. This result is somewhat surprising because a number of d complexes containing phosphine ligands are known to add H readily; this addition results in the formation of dihydrides (66-74). Ther­ modynamically stable η - Η complexes have been prepared (85-91). T h e first of these complexes was the d compound W[P(i-Pr) ] (CO) (H ) (85), for which the H - H distance was found to be 0.75 ± 0.16 Â. 2

1

2

2

1

2

8

2

2

2

6

3

2

3

2

A theoretical study (92) on the related model systems W ( P H ) ( C O ) ( H ) and W ( P H ) ( H ) revealed that the η - Η complex is stabilized by the iracceptor C O ligands, which lower the energy levels of the metal d orbitals 3

3

5

2

2

2

2

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

3

2

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

88

H O M O G E N E O U S TRANSITION M E T A L C A T A L Y Z E D REACTIONS

8c w i t h respect to the c o r r e s p o n d i n g orbitals i n W ( P H ) ( H ) . A s a conse­ q u e n c e , the c a p a b i l i t y of the m e t a l fragment to donate electrons into the a n t i b o n d i n g σ * o r b i t a l of the H m o l e c u l e has d i m i n i s h e d . I n t u r n , this r e s t r i c t i o n prevents the system from u n d e r g o i n g oxidative a d d i t i o n . It has also b e e n f o u n d i n e x p e r i m e n t a l studies i n v o l v i n g d a n d d complexes that the oxidative a d d i t i o n of Η is facilitated b y electron-releasing ligands such as p h o s p h i n e s a n d r e t a r d e d b y IT acceptors s u c h as c a r b o n m o n o x i d e (37). 3

5

2

2

8

1 0

2

F i g u r e 3 displays the e n e r g y profile for the oxidative a d d i t i o n reaction 8a —» 8c. W e calculated an activation energy Δ Ε * of 77 k j m o l " . T h i s v a l u e is m a r k e d l y larger than the activation energies that have b e e n d e t e r m i n e d theoretically for the oxidative a d d i t i o n of H to transition m e t a l complexes c o n t a i n i n g p h o s p h i n e ligands (49). T h i s difference can b e a t t r i b u t e d to the stabilization of the η - Η adduct b y the ττ-acceptor C O ligands. T h e e n e r g y c u r v e m o d e l e d i n F i g u r e 3 ascends steeply d u r i n g the early stages of the reaction. T h e slope indicates that the activation e n e r g y arises largely f r o m the i n i t i a l elongation o f the H - H b o n d distance. 1

2

2

2

T h e fact that the oxidative a d d i t i o n process, 8a —» 8c, exhibits a sizable activation b a r r i e r l e d us to consider an alternative m e c h a n i s m for the last p r o d u c t - f o r m i n g step (e of S c h e m e I) i n the h y d r o f o r m y l a t i o n c y c l e . T h e

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

5.

Z I E G L E R & VERSLUIS

The Hydroformyhtion

89

Reaction

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

100 H

04

Reaction Coordinate Figure

3. Energy

profile for

the oxidative

addition

of H to

CH C(0)Co(CO) .

2

3

3

alternative mechanism assumes that one hydrogen atom of the η - Η com­ plex shifts directly toward the acyl group and thereby forms an aldehyde molecule and regenerates the catalyst H C o ( C O ) . This type of reaction mech­ anism, in which a direct oxidative addition of Η is avoided, has, to our knowledge, never been investigated theoretically for the cobalt-based hy­ droformylation process. 2

2

3

2

The approximate profile for the reaction 9a —> 9b —• 9c is given in Figure 4. The first part of the profile connecting 9a and 9b was obtained by changing the internal coordinates of 9a into those of 9b in a linear and stepwise fashion. A total of six steps was used in the transit. We find the reaction proceeds with a minimal activation barrier (Figure 4). The activation energy for the first part, 9a —» 9b, is therefore essentially equivalent to the reaction enthalpy Δ£ for the reaction 9a —* 9b, which was calculated to be 83 kj mol" . We have also traced the energy surface for the reaction 9b —> 9c by using a similar five-step linear transit procedure. The step 9b —• 9c, which is exo­ thermic with a reaction enthalpy of Δ£ = -26 kj mol" , has a negligible activation energy. Finally, the system, 9c, is lowered in energy by 69 kj m o l as the adduct 9c breaks up into acetaldehyde and H C o ( C O ) in its ground-state conformation la. 1

1

1

3

C H

r £

Ο

1

C

ο 9a

c^ I r c ο 9b

3

_

C

Co—Η

9c

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

90

H O M O G E N E O U S TRANSITION M E T A L C A T A L Y Z E D REACTIONS

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

100+

Reaction Coordinate Figure

Concluding

4. Energy

profile

for

the reaction

9a - » 9b —» 9c.

Rental

We carried out calculations on the elementary steps of the hydroformylation process. Our calculations indicate that the last product-forming step (e of Scheme I, in which aldehyde elimination is induced) has the highest acti­ vation energy with Δ Η * ~ 80-90 kj mol" . Step d was calculated to have a similar activation energy (ΔΗ* ~ 70-80 kj mol" ). The insertion of olefin into the C o - Η bond (steps b-c) was, on the other hand, estimated to have a negligible activation energy of ΔΗ* ~ 10-20 kj mol" . We finally found that the formation of the active catalyst H C o ( C O ) from H C o ( C O ) by C O dissociation (step a) requires 185 kj mol" . More detailed studies of steps a (6), b (93), c (6), d (94), and e (94) of Scheme I can be found elsewhere. 1

1

1

3

4

1

Acknowledgment This investigation was supported by the Natural Sciences and Engineering Research Council of Canada (NSERC). We also acknowledge access to the Cyber-205 installations at the University of Calgary. We are thankful to E . J . Baerends and W. Ravenek from the Free University of Amsterdam for a copy of their latest vectorized version of the H F S - L C A O program system. We would also like to thank Pieter Vernooijs for help with the installation of the program system.

References 1. Heck, R. F. Adv. Organomet. Chem. 1966, 4, 243. 2. Orchin, M.; Rupilius, W. Catal. Rev. 1972, 6, 85. 3. Orchin, M. Acc. Chem. Res. 1981, 14, 25.

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

5.

ZIEGLER & VERSLUIS The Hydroformyhtion Reaction

91

4. Pino, P.; Piacenti, F.; Bianchi, M. In Organic Synthesis via Metal Carbonyls, Vol. II; Wender, I.; Pino, P., Eds.; Wiley: New York, 1977; pp 43-135. 5. Heck, R. F.; Breslow, D. S. J. Am. Chem. Soc. 1961, 83, 4023. 6. Versluis, L.; Ziegler, T.; Baerends, E. J.; Ravenek, W. J. Am. Chem. Soc. 1989, 111, 2018. 7. Versluis, L. Ph.D. Thesis, University of Calgary, 1989. 8. Baerends, E. J.; Ellis, D. E.; Ros, P. Chem. Phys. 1973, 2, 41. 9. Baerends, E. J., Ph.D. Thesis, Vrije Universiteit, Amsterdam, 1975. 10. Ravenek, W. In Algorithms and Applications on Vector and Parallel Comput te Riele, H. J. J.; Dekker, Th. J.; van de Vorst, Η. Α., Eds.; Elsevier: Amsterdam, 1987. 11. Ziegler, T.; Rauk, A. Theor. Chim. Acta 1977, 46, 1. 12. Ziegler, T., Ph.D. Thesis, University of Calgary, 1978. 13. Becke, A. D. J. Chem. Phys. 1986, 84, 4524. 14. Stoll, H.; Golka, E.; Preuss, H. Theor. Chim. Acta 1980, 55, 29. 15. Becke, A. D. J. Chem. Phys. 1988, 88, 2547. 16. Snijders, G. J.; Baerends, E. J.; Vernooijs, P. At. Nucl. Data Tables 1982, 26, 483. 17. Vernooijs, P.; Snijders, G. J.; Baerends, E. J. "Slater Type Basis Functions for the Whole Periodic System"; internal report, Free University of Amsterdam: Netherlands, 1981. 18. Versluis, L.; Ziegler, T. J. Chem. Phys. 1988, 88, 322. 19. Ziegler, T.; Tschinke, V.; Ursenbach, C. J. Am. Chem. Soc. 1987, 109, 4825. 20. Ziegler, T.; Tschinke, V.; Becke, A. D. Polyhedron 1987, 6, 685. 21. Ziegler, T.; Tschinke, V.; Becke, A. D. J. Am. Chem. Soc. 1987, 109, 1351. 22. Ziegler, T.; Cheng, W.; Baerends, E. J.; Ravenek, W. Inorg. Chem. 1988, 27, 3458. 23. Ziegler, T.; Tschinke, V.; Baerends, E. J.; Snijders, J. G.; Ravenek, W.J.Phys. Chem. 1989, 93, 3050. 24. Ziegler, T.; Tschinke, V.; Versluis, L.; Baerends, E. J.; Ravenek, W. Polyhedron 1988, 7, 1625. 25. Fan, L.; Versluis, L.; Ziegler, T.; Baerends, E. J.; Ravenek, W. Int. J. Quantum Chem. 1988, S22, 173. 26. Ziegler, T.; Rauk, Α.; Baerends, E. J. Theor. Chim. Acta 1977, 43, 261. 27. Fan, L.; Ziegler, T. J. Chem. Phys. 1990, 92, 897. 28. Ziegler, T.; Tschinke, V.; Fan, L.; Becke, A. D. J. Am. Chem. Soc. 1989, 111, 9177. 29. Harrod, J. F.; Ziegler, T.; Tschinke, V. Organometallics 1990, 9, 897. 30. Master, A. P.; Soerensen, T. S.; Ziegler, T. Organometallics 1989, 8, 1088. 31. Wermer, P.; Ault, B. S.; Orchin, M. J. Organomet. Chem. 1978, 162, 189. 32. Lewis, Κ. E.; Golden, D. M.; Smith, G. P. J. Am. Chem. Soc. 1984, 106, 3906. 33. Sweany, R. L.; Russell, F. N. Organometallics 1988, 7, 719. 34. Sweany, R. L. Organometallics 1989, 8, 175. 35. Rossi, A. R.; Hoffmann, R. Inorg. Chem. 1975, 14, 365. 36. Versluis, L.; Ziegler, T. J. Am. Chem. Soc. 1990, 112, 6763. 37. Collman, J. P.; Hegedus, L. S. In Principles and Applications of Organotransition Metal Chemistry; Kelly, Α., Ed.; University Science Books: Mill Valley, CA, 1980. 38. Werner, H.; Feser, R. Angew. Chem., Int. Ed. Engl. 1979, 18, 157. 39. Koga, N.; Jin, S. Q.; Morokuma, K. J. Am. Chem. Soc. 1988, 110, 3417. 40. Ziegler, T.; Versluis, L.; Tschinke, V. J. Am. Chem. Soc. 1986, 108, 612. 41. Roe, D. C. Organometallics 1987, 6, 942. 42. Dombek, B. D. Ann. Ν. Y. Acad. Sci. 1983, 415, 176.

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

92

H O M O G E N E O U S TRANSITION M E T A L C A T A L Y Z E D REACTIONS

43. Dombek, B. D. In Adv.Catal.1983, 32, 325. 44. Byers, Β. H.; Brown, T. L. J. Am. Chem. Soc. 1977, 99, 2527. 45. Wayland, Β. B.; Duttahmed, Α.; Woods, Β. Α.; Pierce, R.J.Chem. Soc., Chem. Commun. 1983, 143-144. 46. Fagan, P. J.; Manriquers, J. M.; Marks, T. J.; Vollmer, S. H.; Day, C. S.; Day, V. W. J. Am. Chem. Soc. 1981, 103, 2206. 47. Moloy, K. G.; Marks, T. J. J. Am. Chem. Soc. 1984, 106, 7051. 48. Tam, W.; Wong, W. K.; Gladysz, J. A. J. Am. Chem. Soc. 1979, 101, 1589. 49. Collins, T. J.; Roper, W. R. J. Chem. Soc., Chem. Commun. 1976, 1044. 50. Brown, K. L.; Clark, G. R.; Headford, C. E. L.; Marsden, K.; Roper, W. R. J. Am. Chem. Soc. 1979, 101, 503. 51. Chaudret, Β. N.; Cole-Hamilton, D. J.; Nohr, R. S.; Wilinnson, G.J.Chem. Soc., Dalton Trans. 1977, 1546. 52. Martinho Simoes, J. Α.; Beauchamp, J. L. Chem. Rev. 1990, 90, 629. 53. Collman, J. P.; Winter, S. R. J. Am. Chem. Soc. 1973, 95, 4089. 54. Lane, K. R.; Sallans, L.; Squires, R. R. Organometallics 1985, 4, 408. 55. Collman, J. P.; Roper, W. R. Adv. Organomet. Chem. 1968, 7, 53. 56. James, B. R. In Homogeneous Hydrogenation; Wiley: New York, 1973. 57. Chock, P. B.; Halpern, J. J. Am. Chem. Soc. 1966, 88, 3511. 58. Halpern, J. Acc. Chem. Res. 1970, 3, 386. 59. Lappert, M. F.; Lednor, P. W. Adv. Organomet. Chem. 1976, 14, 345. 60. Crabtree, R. H.; Morehouse, S. M. Adv. Organomet. Chem. 1982, 21, 4210. 61. Harrod, J. F.; Yorke, W. J. J. Inorg. Chem. (USSR) 1981, 20, 1156. 62. Norton, J. Acc. Chem. Res. 1979, 12, 139. 63. Kemmitt, R. D. W.; Smith, M. A. R. Inorg. React. Mech. 1976, 4, 3119. 64. Abis, L.; Sen, Α.; Halpern, J. J. Am. Chem. Soc. 1978, 100, 2915. 65. Halpern, J. Acc. Chem. Res. 1982, 15, 332. 66. Low, J. J.; Goddard, W. Α.,IIIOrganometallics 1986, 5, 609. 67. Low, J. J.; Goddard, W. Α.,IIIJ.Am. Chem. Soc. 1984, 106, 6228. 68. Noell, J. O.; Hay, P. J. J. Am. Chem. Soc. 1982, 104, 4578. 69. Obara, S.; Kitaura, K.; Morokuma, K. J. Am. Chem. Soc. 1984, 106, 7482. 70. Dedieu, Α.; Strich, A. Inorg. Chem. 1979, 18, 2940. 71. Sevin, A. Nouv. J. Chim. 1981, 5, 233. 72. Balazs, A. C.; Johnson, K. H.; Whitesides, G. M. Inorg. Chem. 1982, 21, 2162. 73. Lauher, J. W.; Hoffmann, R. J. Am. Chem. Soc. 1976, 98, 1729. 74. Tatsumi, K.; Hoffmann, R.; Yamamoto, Α.; Stille, J. K. Bull. Chem. Soc. Jpn. 1981, 54, 1857. 75. Rabaâ, H.; Saillard, J.-Y.; Hoffmann, R. J. Am. Chem. Soc. 1986, 108, 4327. 76. Alemdarogly, Ν. H.; Penninger, J. L. M.; Oltay, E. Monatsh. Chem. 1976, 107, 1153. 77. Kovás, I.; Ungváry, F.; Markó, L. Organometallics 1986, 5, 209. 78. Ungváry, F.; Markó, L. Organometallics 1983, 2, 1608. 79. Hoff, C. D.; Ungváry, F.; King, R. B.; Markó, L. J. Am. Chem. Soc. 1985, 107, 666. 80. Lin, Z.; Marks, J. J. Am. Chem. Soc. 1987, 109, 7979. 81. Steigerwald, M. L.; Goddard, W. Α.,IIIJ.Am. Chem.Soc.1984, 106, 308. 82. Jeske, G.; Lauke, H.; Mauermann, H.; Schumann, H.; Marks, T. J. J. Am. Chem.Soc.1985, 107, 8111. 83. Watson, P. J.; Parshall, G. W. Acc. Chem. Res. 1985, 18, 51. 84. Brintzinger, Η. H. J. Organomet. Chem. 1979, 171, 337. 85. Kubas, G. J.; Ryan, R. R.; Swanson, Β. I.; Vergamini, P. J.; Wasserman, H. J. Am. Chem. Soc. 1984, 106, 451.

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.

Downloaded by UNIV OF CALIFORNIA SAN DIEGO on June 9, 2015 | http://pubs.acs.org Publication Date: March 1, 1992 | doi: 10.1021/ba-1992-0230.ch005

5.

Z I E G L E R & VERSLUIS

The Hydroformyhtion

Reaction

93

86. Kubas, G. J.; Unkefer, G. J.; Swanson, B.; Fukishima, E. J. Am. Chem. Soc. 1986, 108, 7000. 87. Wasserman, H. J.; Kubas, G. J.; Ryan, R. R. J. Am. Chem. Soc. 1986, 108, 2294. 88. Upmacis, R. K.; Gadd, G. E.; Poliakoff, M.; Simpson, M. B.; Turner, J. J.; Whyman, R.; Simpson, A. F. J. Chem. Soc., Chem. Commun. 1985, 27. 89. Sweany, R. L. J. Am. Chem. Soc. 1985, 107, 2374. 90. Church, S. P.; Grevels, F. W.; Hermann, H.; Schaffner, K. J. Chem. Soc., Chem. Commun. 1985, 30. 91. Crabtree, R. H.; Lavin, M. J. Chem. Soc., Chem. Commun. 1985, 794. 92. Hay, P. J. J. Am. Chem. Soc. 1987, 109, 705. 93. Versluis, L.; Ziegler, T.; Fan, L. Inorg. Chem. 1990, 29, 5340. 94. Versluis, L.; Ziegler, T. Organometallics 1990, 9, 2985. RECEIVED for review October 19, 1990. ACCEPTED revised manuscript June 12, 1991.

In Homogeneous Transition Metal Catalyzed Reactions; Moser, W., et al.; Advances in Chemistry; American Chemical Society: Washington, DC, 1992.