Đang tải... (xem toàn văn)
Động học phản ứng giữa gốc propargyl với phân tử nước và gốc hyđroxyl được nghiên cứu bằng lý thuyết trạng thái chuyển tiếp (TST), lý thuyết trạng thái chuyển tiếp biến phân (VTST) và RRKM trên cơ sở cơ chế của cả 2 phản ứng đã được nghiên cứu bằng phương pháp phiếm hàm mật độ (DFT).
Tạp chí Khoa học Cơng nghệ 140 (2020) 050-055 Nghiên cứu động học phản ứng gốc propargyl (C3H3) với phân tử nước (H2O) gốc hyđroxyl (OH) pha khí Study on Kinetics of the Reactions between Propargyl Radical (C3H3) and Water Molecule (H2O) and Hydroxyl Radical (OH) in the Gas Phase Phạm Văn Tiến* , Nguyễn Ngọc Tuệ Trường Đại học Bách khoa Hà Nội - Số 1, Đại Cồ Việt, Hai Bà Trưng, Hà Nội Đến Tòa soạn: 26-11-2018; chấp nhận đăng: 20-01-2020 Tóm tắt Động học phản ứng gốc propargyl với phân tử nước gốc hyđroxyl nghiên cứu lý thuyết trạng thái chuyển tiếp (TST), lý thuyết trạng thái chuyển tiếp biến phân (VTST) RRKM sở chế phản ứng nghiên cứu phương pháp phiếm hàm mật độ (DFT) Kết tính tốn -34 -1 -1 hệ phản ứng C3H3 + H2O khó xảy điều kiện thường (k = 4,10x10 cm molecule s 300 K, 1atm), hệ phản ứng C3H3 + OH xảy nhanh với số tốc độ -11 0,35 -1 -1 k(T) = 1,39x10 T exp(-26,42/T) cm molecule s phụ thuộc vào nhiệt độ khoảng 300-2100 K Ngồi ra, kết tính số tốc độ hệ phản ứng thứ hai phù hợp với kết thực nghiệm Hansen Miller (500-2100 K) kết Eiteneer Frenklach (1100-2100 K) Từ khóa: Động học phản ứng, gốc propargyl, TST, nước, gốc hyđroxyl Abstract Kinetics of the reactions between propargyl radicals with water molecule and hydroxyl radical were investigated by TST, VST and RRKM theories on the basis of the mechanisms of both reactions researched by means of density functional theory (DFT) The calculated results show that the C3H3 + H2O reaction is -34 -1 -1 unlikely to occur at room temperature (k = 4,10x10 cm molecule s at 300 K, 1atm), while the C3H3 + -11 0,35 -1 -1 OH system reacts very fast with a rate constant k(T) = 1,39x10 T exp(-26,42/T) cm molecule s depending on the temperature in the range 300-2100 K In addition, the results of the second reaction are in good agreement with the experimental results of Hansen and Miller (500-2100 K) and Eiteneer, B and Frenklach (1100-2100 K) Keywords: Kinetics, propargyl radical, TST, water, hydroxyl radical oxi hóa khí cháy [5] Vì vậy, phản ứng gốc OH quan tâm nghiên cứu lý thuyết thực nghiệm nhiều năm qua [6] Phản ứng OH với C3H3, mắt xích trung gian quan trọng phản ứng cháy nhiên liệu hóa thạch có mặt khí nghiên cứu [7] Eiteneer Frenklach đo số tốc độ phản ứng C3H3 + OH HCO + C2H3 khoảng nhiệt độ 1150 - 2130 K, áp suất 0,9 - 1,9 atm, khí trơ Ar, kết cho thấy phản ứng xảy nhanh với số tốc độ k = 2,21x10-11 cm3 molecule-1 s-1 [8] Hansen Miller [9] nghiên cứu lý thuyết cho phản ứng C3H3 + OH CO + C2H4 C3H3 + OH HCO + C2H3 khoảng nhiệt độ 500 - 2100K, tính số tốc độ phản ứng tương ứng 4,98x10-11 cm3 molecule-1 s-1 8,3x1011 cm3 molecule-1 s-1 Tuy nhiên chưa có nghiên cứu đầy đủ cho toàn đường phản ứng xảy hệ C3H3 + OH Vì vậy, cần có khảo sát chi tiết, đầy đủ động học hệ phản ứng Mở đầu Nước hợp chất phổ biến khí hình thành hầu hết trình đốt cháy hợp chất hữu Đã có nhiều nghiên cứu nước với gốc tự như: CCl2, CH, C2H, CCl, C2F [13], chưa có cơng trình cơng bố kết nghiên cứu lý thuyết thực nghiệm phản ứng gốc C3H3 với H2O Do vậy, việc nghiên cứu động học phản ứng gốc C3H3 với H2O mở rộng hiểu biết khả phản ứng C3H3, mà giải thích hình thành sản phẩm q trình chuyển hóa gốc propargyl Gốc tự hyđroxyl (OH) gốc đóng vai trò quan trọng phản ứng cháy hóa học khí Trong khí quyển, OH ảnh hưởng đến hình thành phá hủy tầng ozon ảnh hưởng tới nhiệt độ toàn cầu [4] OH tác nhân Địa liên hệ: Tel.: (+84) 917 604 068 Email: tien.phamvan@hust.edu.vn 50 Tạp chí Khoa học Công nghệ 140 (2020) 050-055 Phương pháp tính Kết thảo luận Việc tính động học phản ứng C3H3 với chất khác (H2O OH) pha khí thực áp suất thường (1atm) khoảng nhiệt độ từ 300 đến 2100 K theo lý thuyết TST, VTST RRKM [10] Động học phản ứng bao gồm: xác định số tốc độ hướng đầu vào, số tốc độ tạo thành sản phẩm số tốc độ chung hệ (ktot) Việc tính số tốc độ (k) theo TST sử dụng phần mềm Chemrate [11] Việc tính k theo VTST RRKM thực phần mềm Variflex, có giải phương trình master liên quan đến trình có nhiều bước chuyển tiếp [12] Đối với q trình khơng có TS, đường lượng cực tiểu (MEP) biểu diễn sau làm khớp với hàm Morse: E(R)=De [1-e-β(R-Re ) ]2 Trong 3.1 Động học hệ phản ứng C3H3 + H2O Để tính số tốc độ cho hệ phản ứng C3H3 + H2O, sử dụng đường phản ứng có hàng rào lượng thấp so với đường phản ứng lại, đường thứ qua T0/P7 tạo sản phẩm P7, đường lại qua T0/P6 sinh P6 Từ bề mặt (hình 1) nhận thấy, hai đường phản ứng qua trạng thái chuyển tiếp T0/P6 T0/P7 Do đó, để tính số tốc độ cho đường phản ứng này, mơ hình tính tốn thuyết phức chất hoạt động thuyết RRKM sử dụng thông qua phần mềm Chemrate [11] Các thông số LennardJones cho chất phản ứng cần sử dụng q trình tính tốn số tốc độ hai đường phản ứng nói lấy từ tài liệu [13], C3H3 = 4,76 Å; /K = 252 K H2O = 2,605; /K = 572,4 K Các thơng số cấu trúc hình học (góc liên kết, độ dài liên kết, góc nhị diện) tần số dao động chất phản ứng, trạng thái chuyển tiếp, trạng thái trung gian sản phẩm liên quan đến đường phản ứng nói lấy từ file liệu đầu chất tối ưu hóa mức B3LYP/6311++G(3df,2p) Trong đó, tần số ảo T0/P6 T0/P7 lược bỏ để bảo toàn tổng trạng thái phức chất hoạt động đó, Re giá trị cân độ dài liên kết (R) cấu trúc trạng thái trung gian (IS); β số; De lượng đứt gãy liên kết khơng tính đến giá trị lượng dao động điểm không (ZPE) Hằng số tốc độ phản ứng theo RRKM tính sau: k T hQr Qv E0 J 1 G E, J e E k T k E, J / J 0 dE b (1) Trong đó: h số Planck; Qr, Qv hàm phân bố lượng quay, dao động; E0 lượng hoạt hóa phản ứng; J số lượng tử; tần số va chạm; G‡ tổng trạng thái giai đoạn chuyển tiếp; k(E,J) số tốc độ phản ứng vi cổ điển Giá trị số tốc độ phản ứng đường phản ứng nói (k1, k2), số tốc độ tổng (ktot = k1 + k2, cm3 molecule-1 s-1) nhiệt độ khác trình bày bảng Trong k1 số tốc độ đường phản ứng thứ qua trạng thái chuyển tiếp T0/P7 (27,83 kcal/mol), k2 số tốc độ ứng với đường phản ứng thứ hai qua T0/P6 (30,70 kcal/mol) Phương trình master tập hợp phương trình vi phân phức tạp, mơ tả q trình tiến triển theo thời gian hàm tổng mật độ va chạm hệ: m d i (t ) i Pij i (t ) i (t ) (ki1 ki ) i (t ) (2) dt j 1 Bảng Các giá trị số tốc độ, k [cm3 molecule-1 s-1], phụ thuộc vào nhiệt độ phản ứng C3H3 + H2O Trong đó, i biểu diễn tốc độ đo hai chất phản ứng; m số hạt chọn cho đóng góp hạt thứ m không đáng kể số tốc độ phản ứng lưỡng phân tử; tần số va chạm; ki1(E) ki2(E) số tốc độ vi cổ điển riêng cho phân hủy phân ly trở lại; P ij xác suất chuyển lượng từ hạt j đến hạt i sau va chạm Pij = Aj exp[−α(Ej − Ei)], ji (3) Trong đó, thơng số liên quan đến hệ số chuyển lượng ; Aj số chuẩn hóa đạt từ điều kiện chuẩn hóa T (K) k1 k2 ktot 300 500 700 900 1100 1300 1500 1700 1900 2100 4,09x10-34 1,37x10-25 7,65x10-22 1,08x10-19 2,85x10-18 2,98x10-17 1,78x10-16 7,32x10-16 2,33x10-15 6,15x10-15 1,00x10-36 2,26x10-27 2,90x10-23 6,59x10-21 2,36x10-19 3,09x10-18 2,18x10-17 1,03x10-16 3,63x10-16 1,04x10-15 4,10x10-34 1,39x10-25 7,94x10-22 1,15x10-19 3,09x10-18 3,29x10-17 2,00x10-16 8,35x10-16 2,69x10-15 7,19x10-15 Các số tốc độ sử dụng để làm khớp với biểu thức tính số tốc độ k phụ thuộc vào nhiệt độ T khoảng nhiệt độ từ 300 đến 2100 K theo công thức định luật Arrhenius [14] 51 Tạp chí Khoa học Cơng nghệ 140 (2020) 050-055 Hình 1: Bề mặt chi tiết hệ phản ứng C3H3 với H2O Năng lượng (kcal/mol) tính mức CCSD(T)/6-311++G(3df,2p)//B3LYP/6-311++G(3df,2p) + ZPE -20 Chính thế, đóng góp đường phản ứng thứ hai vào số tốc độ chung hệ không đáng kể vùng nhiệt độ thấp (dưới 500 K) Do đó, khảo sát hệ vùng nhiệt độ khơng cần tính đến đường phản ứng nói Còn vùng nhiệt độ cao (trên 500 K) đóng góp đường phản ứng đáng kể việc tạo sản phẩm P6 P7 -1 -1 Log(k, cm molecule s ) phản ứng k1 qua trạng thái chuyển tiếp T0/P7 (27,83 kcal/mol) thấp hẳn so với đường thứ hai qua trạng thái chuyển tiếp T0/P6 (30,7 kcal/mol) k1: RA > 0/P7 > P7 k2: RA > 0/P6 > P6 ktot = k1 + k2 -15 -25 -30 -35 3.2 Động học hệ phản ứng C3H3 + OH 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 3.2.1 Hằng số tốc độ trình C3H3 + OH IS1 C3H3 + OH IS2 1000/T (K) Hình Sự phụ thuộc số tốc độ k(T) vào nhiệt độ T hệ C3H3 + H2O Bề mặt phản ứng cho thấy trình đầu vào C3H3 + OH CH2=C=CH-OH (IS1) C3H3 + OH HO-CH2-C≡CH (IS2) trạng thái chuyển tiếp Để giải khó khăn này, phương pháp cải biến áp dụng sở lượng tự cực đại (G‡max) Trước tiên quét bề mặt công nhóm OH vào vị trí C hai đầu mạch C3H3 để tìm đường lượng cực tiểu phương pháp khác Khoảng cách C1-O trạng thái trung gian IS1 thay đổi từ vị trí cân 1,38Å đến khoảng cách phân ly 7,38Å với bước nhảy 0,1Å Tương tự vậy, khoảng cách C3-O trạng thái trung gian IS2 Sự phụ thuộc k vào T đường phản ứng biểu diễn hình Kết bảng cho thấy, hai đường phản ứng có số tốc độ nhỏ nhiệt độ thường, chứng tỏ hệ phản ứng khó xảy điều kiện này, điều hồn tồn phù hợp với phân tích lượng Khi nhiệt độ tăng, đóng góp đường tăng lên, nhiệt độ cao giúp chất phản ứng dễ dàng việc vượt qua hàng rào lượng T0/P6 T0/P7 Ngoài ra, đường phản ứng thứ có số tốc độ k1 lớn so với số tốc độ k2 đường thứ hai, điều hồn tồn phù hợp đường 52 Tạp chí Khoa học Cơng nghệ 140 (2020) 050-055 thay đổi từ vị trí cân 1,43Å đến khoảng cách phân ly 7,43Å với bước nhảy 0,1Å Khi tính tốn số tốc độ cho giai đoạn này, cấu trúc ứng với G‡max sử dụng làm trạng thái chuyển tiếp bước C3H3 + OH IS1 C3H3 + OH IS2 Do giá trị G‡max đặt cấu trúc tính được, nên tần số dao động momen quán tính ngoại suy cách tuyến tính điểm Các thông số Lennard-Jone C3H3: σ = 4,76Å; ε/K = 252K giá trị L-J OH: σ = 2,75Å; ε/K = 80 K [13] Các thông số LJ tất trạng thái trung gian tính dựa kết quét bề mặt va chạm khí Ar với chất trung gian Các thơng số hình học khác tối ưu khoảng cách C1-O C3-O Đối với cấu trúc tối ưu theo phương pháp B3LYP/6311++G(3df,2p), tính 3N - tần số dao động gradient tương ứng Để đạt độ tin cậy cao hơn, lượng điểm đơn tính theo phương pháp CCSD(T)/6-311++G(3df,2p) cho cấu trúc hình học tối ưu với điểm xung quanh vị trí có lượng tự cực đại Ngồi ra, Morse, E(R) = De[1-exp(β(R-Re))]2, sử dụng để làm khớp với tính mức B3LYP/6-311++G(3df,2p) Đường lượng cực tiểu biểu diễn trình đứt gãy khơng có trạng thái chuyển tiếp CH2CHCHO C2H3 + HCO đạt việc tính tốn đường cong dọc theo tọa độ phản ứng liên kết C-C trạng thái chuyển tiếp IS5 Trong đó, độ dài liên kết C-C kéo dài từ giá trị cân 1,47Å đến 8,47Å với bước nhảy 0,1Å thơng số khác tối ưu hồn toàn Đường cong đứt gãy làm khớp với Morse sử dụng để tính gần đường lượng cực tiểu trạng thái chuyển tiếp cải biến tính tốn tốc độ phản ứng Các thông số thu sau: De = 32051,5 cm-1; β = 1,159Å-1; Re = 1,47Å Trong đó, R khoảng cách C1-O C3-O C C3H3 với O OH; Re = 1,38Å IS1 1,43Å IS2 (giá trị cân R) Đối với IS1, giá trị β xác định 1,346Å-1 De = 30966,36 cm-1 lượng đứt gãy liên kết C-O kcal/mol, khơng tính đến lượng điểm không (ZPE); tương tự IS2, giá trị β De tương ứng 1,419 Å-1 29958,61 cm-1 Từ Morse, momen quán tính tần số dao động tính được, cho phép xác định G Chúng tơi tìm kiếm giá trị cực đại G‡max nhiệt độ khác từ 300 đến 2100 K Sử dụng giá trị cực đại G‡max, thu số tốc độ VTST, k01 k02 tương ứng cho kết hợp C3H3 + OH CH2=C=CH-OH (IS1) C3H3 + OH HO-CH2-C≡CH (IS2) Kết bảng Các đường phản ứng liên quan đến dịch chuyển tách nguyên tử H liên kết C-C tính đến hiệu ứng đường hầm Eckart việc tính tốn số tốc độ phản ứng Hằng số tốc độ đường phản ứng (k1, k2, k3, k4, k5, k6) số tốc độ tổng (ktot = k1 + k2 + k3 + k4 + k5 + k6) nhiệt độ khác trình bày bảng 3.2.4b Trong đó, k1 số tốc độ đường RA → IS1 → IS5 → PR4; k2: RA → IS1 → IS5 → PR8; k3: RA → IS1 → IS5 → PR5; k4: RA → IS1 → PR6; k5: RA → IS1 → PR1 k6: RA → IS2 → PR2 3.2.2 Hằng số tốc độ cho chế phản ứng C3H3 + OH → sản phẩm Như phân tích trên, đường phản ứng tạo sản phẩm PR1, PR2, PR4, PR5, PR6 PR8 đường phản ứng thuận lợi mặt lượng so với đường lại Chính vậy, đường phản ứng sử dụng để tính số tốc độ phản ứng (hình 3) Kết bảng 3.2.4b cho thấy, nhiệt độ thường đường phản ứng tạo sản phẩm PR1, PR2 PR6 có số tốc độ phản ứng nhỏ so với đường phản ứng tạo sản phẩm PR4, PR5 PR8 Điều phù hợp với bề mặt chúng phải trải qua trạng thái chuyển tiếp cao -2,78, -5,99 -6,62 kcal/mol so với đường k1 (-19,68 kcal/mol) Khi nhiệt độ phản ứng tăng, đóng góp đường phản ứng tăng lên, nhiệt độ cao giúp chất phản ứng dễ dàng việc vượt qua hàng rào lượng Đường phản ứng k1 đường phản ứng k2 có số tốc độ cao đường lại, phù hợp với bề mặt năng, k1 đường có hàng rào lượng thấp nhất, k2 đường tạo sản phẩm PR8 không thông qua trạng thái chuyển tiếp giai đoạn cuối Hình Các đường phản ứng hệ C3H3 + OH Năng lượng (kcal/mol) tính mức CCSD(T)/6-311++G(3df,2p)//B3LYP/6311++G(3df,2p) + ZPE Đồng thời k1 k2 có giá trị gần nhau, phù hợp với kết phân tích khả cạnh tranh 53 Tạp chí Khoa học Cơng nghệ 140 (2020) 050-055 sản phẩm PR4 PR8 Mặt khác, số tốc độ trung bình k1 k2 thu khoảng nhiệt độ 500 - 2100 K có giá trị tương ứng 6,66x10-11 6,62x10-11 cm3 molecule-1 s-1 phù hợp tốt so với giá trị tính Hansen Miller [9], tương ứng 4,8x10-11 8,3x10-11 cm3 molecule-1 s-1 Hơn số tốc độ k2 khoảng nhiệt độ 1100 2100 K có giá trị 6,92x10-11 cm3 molecule-1 s-1 phù hợp tốt với giá trị đo thực nghiệm Eiteneer, B Frenklach, M 2,21x10-11 cm3 molecule-1 s-1 [8] Điều cho thấy kết tính tốn lý thuyết đáng tin cậy dạng: k(T) = molecule-1 s-1 1,39x10-11T0,35exp(-26,42/T) cm3 -11 -1 -1 Log(k, cm molecule s ) -10 -12 k1: RA > PR4 k2: RA > PR8 k3: RA > PR5 k4: RA > PR6 k5: RA > PR1 k6: RA > PR2 ktot Để khảo sát phụ thuộc số tốc độ phản ứng k vào nhiệt độ T, phương trình cải biến Arrhenius sử dụng dạng k(T) = A.Tn.exp(-B/T), A, B n thừa số cần xác định, với B = Ea/R Bằng việc sử dụng phần mềm Origin với liệu đầu vào giá trị k nhiệt độ T tương ứng, sau làm khớp với đường cong động học k(T), giá trị đầu thu sau: A = 1,39x10-11; B = 26,42 n = 0,35; biểu thức số tổng cộng ktot phụ thuộc vào T có -13 -14 -15 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 1000(K)/T Hình Sự phụ thuộc số tốc độ k(T) vào nhiệt độ hệ phản ứng C3H3 + OH Bảng Giá trị cực đại G‡max ứng với độ dài liên kết C-O số tốc độ phản ứng cho C3H3 + OH → IS1 C3H3 + OH → IS2 tính theo lý thuyết VTST T (K) 300 500 700 900 1100 1300 1500 1700 1900 2100 d (C1-O) (Å) 2,17 2,15 2,12 2,08 2,03 1,98 1,94 1,91 1,87 1,82 G‡max (kcal/mol) -2,33 -3,43 -4,45 -5,42 -6,36 -7,20 -8,03 -8,81 -9,59 -10,4 k01 (cm molecule-1s-1) 5,20x10-10 5,48x10-10 5,94x10-10 6,45x10-10 6,97x10-10 7,31x10-10 7,67x10-10 7,98x10-10 8,34x10-10 8,76x10-10 G‡max (kcal/mol) -2,30 -3,38 -4,41 -5,42 -6,32 -7,17 -7,97 -8,77 -9,44 -10,2 d (C3-O) (Å) 2,21 2,19 2,16 2,12 2,07 2,02 1,97 1,94 1,90 1,85 k02 (cm3molecule-1s-1) 4,91x10-10 5,20x10-10 5,75x10-10 6,43x10-10 6,85x10-10 7,23x10-10 7,51x10-10 7,89x10-10 8,01x10-10 8,42x10-10 Bảng Hằng số tốc độ, k [cm3 molecule-1 s-1], số đường phản ứng hệ C3H3 + OH T (K) 300 500 700 900 1100 1300 1500 1700 1900 2100 k1 k2 k3 k4 k5 k6 ktot 5,27x10-11 5,63x10-11 6,10x10-11 6,48x10-11 6,72x10-11 6,86x10-11 6,97x10-11 7,03x10-11 7,07x10-11 7,12x10-11 5,21x10-11 5,56x10-11 6,06x10-11 6,41x10-11 6,66x10-11 6,83x10-11 6,94x10-11 7,00x10-11 7,03x10-11 7,05x10-11 5,33x10-12 7,24x10-12 1,12x10-11 1,61x10-11 2,03x10-11 2,42x10-11 2,83x10-11 3,25x10-11 3,59x10-11 4,18x10-11 3,28x10-14 5,35x10-13 2,23x10-12 3,47x10-12 3,89x10-12 4,02x10-12 4,41x10-12 4,62x10-12 4,86x10-12 5,12x10-12 3,09x10-14 5,22x10-13 2,08x10-12 3,39x10-12 3,81x10-12 3,99x10-12 4,37x10-12 4,55x10-12 4,80x10-12 5,06x10-12 2,77x10-15 8,41x10-14 4,36x10-13 8,58x10-13 1,15x10-12 1,36x10-12 1,57x10-12 1,78x10-12 2,02x10-12 2,16x10-12 1,10x10-10 1,20x10-10 1,38x10-10 1,53x10-10 1,63x10-10 1,70x10-10 1,78x10-10 1,84x10-10 1,89x10-10 1,96x10-10 54 Tạp chí Khoa học Công nghệ 140 (2020) 050-055 Kết luận and propyne flames J Combust Flame, 156, (2009) 2153-2164 Bằng việc sử dụng lý thuyết TST, VTST RRKM, số tốc độ cho đường phản ứng hệ C3H3 + H2O C3H3 + OH tính tốn điều kiện p=1 atm T=300-2100 K Kết tính động học cho thấy, hệ phản ứng gốc propargyl với phân tử nước xảy chậm, nhiệt độ thường (T = 300 K) số tốc độ k = 4,10x10-34 cm3 molecule-1 s-1 Kết hoàn toàn phù hợp với phân tích bề mặt đường phản ứng hệ C3H3 + H2O phải trải qua hàng rào lượng cao, thấp 27,83 kcal/mol (hình 1) Trong đó, phản ứng gốc propargyl với gốc hyđroxyl lại xảy nhanh (k = 1,10x10-10 cm3 molecule-1 s-1, T = 300 K) Ngồi kết tính động học cho hệ phản ứng thứ hai phù hợp với kết thực nghiệm nhóm tác giả Hansen Miller (500 - 2100 K) [9] Eiteneer, B Frenklach, M (1100 - 2100 K) [8] [10] Wardlaw, D M.; Marcus, R A RRKM Reaction Theory for Transition States of any Looseness Chem Phys Lett., 110, (2013) 230-234 [11] Vladimir, M V B.; Tsang, W.; Zachariah, M R.; Knyazev, V D.; Sean, M W ChemRate, MD 20899 NIST, Gaithersburg, USA, RRKM/Master Equation Modeling (2011) [12] Klippenstein, S J.; Wagner, A F.; Dunbar, R C.; Wardlaw, D M.; Robertson, S H VARIFLEX Argonne National Laboratory, Argonne, IL (1999) [13] Jasper, A W.; Miller, J A Lennard-Jones parameters for combustion and chemical kinetics modeling from full-dimensional intermolecular potentials Combust and Flame, (2013) 1-30 [14] Steinfeld, J I.; Francisco, J S.; Hase, W L Chemical Kinetics and Dynamics Prentice-Hall, Inc., New Jersey, USA (1989) Tài liệu tham khảo [1] Carl, S A.; Nguyen, H M T.; Elsamra, R MI.; Nguyen, M T.; Peeters, J Pulsed laser photolysis and quantum chemical-statistical rate study of the reaction of the ethynyl radical with water vapor J Chem Phys., 122, (2005) 114307 [2] Liu, Y Z.; Zhang, Z Q.; Pei, L S.; Chen, Y.; Chen, C X Reaction kinetic studies of CCl2 with C2H2 and H2O molecules Chem Phys., 303, (2004) 255 - 263 [3] Zabarnick, S.; Fleming, J W.; Lin, M C Temperature dependence of CH radical reactions with H2O and CH2O Symp Int Combust Proc., 21, (1998) 713-719 [4] Miller J A Concluding Remarks Faraday Discuss., 119, (2001) 461 [5] Cheung R.; Li, K F.; Wang, S.; Pongetti, T J.; Cageao, R P.; Sander, S P.; and Yung, Y L Atmospheric hydroxyl radical (OH) abundances from ground-based ultraviolet solar spectra: an improved retrieval method Applied Optics, 47, (2008) 62776284 [6] Frost, M J.; Sharkey, P.; Smith, I W M Reaction between hydroxyl (deuteroxyl) radicals and carbon monoxide at temperatures down to 80 K: experiment and theory J Phys Chem., 97, (1993) 12254 [7] Dong, F.; Wang, S.; Kong, F Reaction of propargyl with oxygen J Phys Chem A, 107, (2003) 9374 9379 [8] Eiteneer, B.; Frenklach, M Experimental and Modeling Study of Shock-Tube Oxidatioon of Acetylene Int J Chem Kinet., 35, (2003) 391-414 [9] Hansen, N.; Miller, J A.; Westmoreland, P R.; Kasper, T.; Kohse-Hoinghaus, K.; Wang, J.; Cool, T A Isomer-specific combustion chemistry in allene 55 Tạp chí Khoa học Cơng nghệ 140 (2020) 050-055 56 ... K Kết tính động học cho thấy, hệ phản ứng gốc propargyl với phân tử nước xảy chậm, nhiệt độ thường (T = 300 K) số tốc độ k = 4,10x10-34 cm3 molecule-1 s-1 Kết hoàn toàn phù hợp với phân tích bề... thuộc k vào T đường phản ứng biểu diễn hình Kết bảng cho thấy, hai đường phản ứng có số tốc độ nhỏ nhiệt độ thường, chứng tỏ hệ phản ứng khó xảy điều kiện này, điều hồn tồn phù hợp với phân tích... cho chế phản ứng C3H3 + OH → sản phẩm Như phân tích trên, đường phản ứng tạo sản phẩm PR1, PR2, PR4, PR5, PR6 PR8 đường phản ứng thuận lợi mặt lượng so với đường lại Chính vậy, đường phản ứng sử