Theoretical study on reactions of ground state boron atom with ethylene(c 2 H 4 ), allene (H 2 CCCH 2 ), and methylacetylen (CH 3 CCH) C. H. Huang 1, C. H. Kao 1, H. L. Sun 1, A. H. H. Chang 1, and R. I. Kaiser 2 (1) Department of Chemistry, National Dong Hwa University, Shoufeng, Hualien, Taiwan, Fax: 011-886-38633570, air0319@yahoo.com.tw, (2) Department of Chemistry, University of Hawaii at Manoa, Honolulu, HI 96822
Abstract The reactions of ground-state boron atom with prototype unsaturated hydrocarbons : ethylene, allene, and methylacetylene, respectively, are investigated theoretically to explore probable routes and dynamics. Both abstraction and insertion channels on C 2 H 4 B and C 3 H 4 B adiabatic doublet ground-state potential energy surfaces for reactions with ethylene, and isomerpair,allene and methylacetylene, respectively, are Characterized by utilizing the unrestricted B3LYP/6-311G(d,p) level of theory and the CCSD(T)/cc-pVTZ calculations. With facilitation of RRKM and variational RRKM rate constants at collision energies of 0-10 kcal/mol, the most probable paths, thus reaction mechanisms, are determined. The corresponding rate equations are then solved such that the evolutions of concentrations of collision complexes, intermediates, and dissociation products versus time are obtained. As a result, the final products and yields are identified-assuming the energy randomization is complete.
Strategy Ab initio calculations on doublet C 2 H 4 B and CH 3 CCCHB and H 2 CCCH 2 B ground state surfaces Reaction paths for each collision complex Capturing cross-sections (σ cap 's) of forming all collision complexes Unimolecular rate constants Most probable paths (reaction mechanism) Solve rate equations Product yields
Theoretical methods Ab initio electronic structure calculation for reaction paths B3LYP/6-311G(d,p) optimized geometry, harmonic frequencies CCSD(T)/cc-pVTZ energy RRKM and variational RRKM rate constant -- For reaction A * k A, where A*: energized reactant P A P A : transition state P : product A RRKM rate constant: σ W ( E E ) k( E) = h ρ( E) σ where : symmetry factor : number of state of W A : density of state of A* ρ -- For barrierless reactions, ie. simple bond breaking reaction : variational RRKM, the geometry where W is the transition state R = 0 A + A P
methods Capturing cross-section σ cap -- For long-range intermolecular potential of a bimolecular reaction, A+B P: V ( R σ cap C ) = 6, where R : distance between centers of mass of two reactants: A-B R R C 1 1 3 3 ( E) = 3π ( ), or σ cap C ------ Langevin model 4E -- now there are 2 or 3 collision complexes: C(1) V1( R1 ) = 6 R C(2) V2( R2 ) =, R 6 R M σ c1 σ c2 M 1 2 1 ( C(1) ) 1 ( C(2) ), R 3 3 2 1 : reaction coordinate of forming collision complex C : reaction coordinate of forming collision complex C 2 1 Solve rate equations concentation evolutions product yields
B( 2 P) + C 2 H 4 2 collision complexes
B( 2 P) + H 2 CCCH 2 3 collision complexes
B( 2 P) + CH 3 CCH 2 collision complexes
B( 2 P) + C 2 H 4 C1 paths and most probable paths(highlighted)
B( 2 P) + C 2 H 4 C2 paths and most probable paths(highlighted)
B( 2 P) + C 2 H 4 reaction mechanism ( most probable paths )
C1 rate equations based on reaction mechanism: d[c1] = -(k + k + k )[c1] + k [i1] + k [c2] 1 3 2-1 -3 dt d[c2] = -k -3[c2] + k3[c1] dt d[i1] = -(k + k )[i1] + k [c1] + k [i6] 5-1 1-5 dt d[i6] = k [i1] - k [i6] 5-5 dt d[i2] = -(k + k + k )[i2] + k [c1] + k [i4] + k 6 7 19 2-6 dt d[i4] = -(k + k )[i4] + k [i5] + k [i2] 10-6 -10 6 dt d[i5] = -(k + k + k )[i5] + k [i9] + k [i4] 20 14-10 -18 10 dt -19 [i9] d[i9] = -(k -18 + k dt d[p4] = k 7[i2] dt d[p9] = k [i5] 20 dt d[p14] = k14[i5] dt -19 )[i9] + k 19 [i2]
C1 evolution
C1 evolution
C2 evolution
C2 evolution
product yields: B( 2 P) + C 2 H 4 B + (p4) + H 1 : (p9) + H 2.71 : (p14) + H 1.56
B( 2 P) + H 2 CCCH 2 i5 paths and most probable paths(highlighted)
B( 2 P) + H 2 CCCH 2 i8 paths and most probable paths(highlighted)
B( 2 P) + H 2 CCCH 2 C1a paths and most probable paths(highlighted)
B( 2 P) + H 2 CCCH 2 reaction mechanism ( most probable paths )
i5 evolution
i8 evolution
C1a evolution
product yields: B( 2 P) + H 2 CCCH 2 B + (p13) + H 1 (p24) + H 0.29 :
B( 2 P) + CH 3 CCH C1 paths and most probable paths(highlighted)
B( 2 P) + CH 3 CCH C2 paths and most probable paths(highlighted)
B( 2 P) + CH 3 CCH reaction mechanism ( most probable paths )
C1 evolution
C2 evolution
product yields: B( 2 P) + CH 3 CCH B + (p1) + H
summary Barrierless B + C 2 H 4, H 2 CCCH 2,CH 3 CCH reactions have been investigated theoretically by combining ab initio calculation, RRKM and variational RRKM theory, and Langevin model. Reaction paths, most probable paths (reaction mechanisms), product yields are predicted.