Tổng quan luận án
Công trình luận án tiến sĩ hóa học của tác giả Phạm Văn Tiến với đề tài "Nghiên cứu lý thuyết cơ chế phản ứng và động học của gốc propargyl với một số phân tử trong pha khí" được thực hiện tại Bộ môn Hóa lý thuyết và Hóa lý, Trung tâm Khoa học tính toán thuộc Khoa Hóa học, Trường Đại học Sư phạm Hà Nội, dưới sự hướng dẫn của PGS. Nguyễn Thị Minh Huệ và PGS. Lê Kim Long.
Tính cấp thiết và khoảng trống nghiên cứu
Quá trình đốt cháy nhiên liệu hóa thạch trong các hoạt động sản xuất công nghiệp và giao thông vận tải là nguồn phát thải chính sinh ra các khí như COx, NOx, hơi nước và nhiều gốc hydrocacbon trung gian. Trong số các sản phẩm trung gian này, gốc propargyl ($C_3H_3$) giữ vai trò quan trọng trong việc hình thành các hợp chất thơm đầu tiên (như benzen $C_6H_6$, gốc phenyl $C_6H_5$) – vốn được coi là chất nền cho sự hình thành hydrocacbon thơm đa vòng (PAHs) và muội than gây ô nhiễm môi trường. Đồng thời, gốc $C_3H_3$ tham gia vào quá trình biến đổi hóa học các oxit nitơ ($NO_x$), phản ứng với oxy nguyên tử và oxy phân tử phát xạ năng lượng làm tăng đặc tính xanh trong phát quang hóa học của sự cháy. Mặc dù các phương pháp thực nghiệm như phổ khối lượng hay phổ laser có thể theo dõi sự tồn tại của các gốc hydrocacbon, nhưng việc đo đạc thực nghiệm gặp khó khăn trong việc mô tả chi tiết bản chất vi mô, cấu trúc hình học của trạng thái chuyển tiếp và cơ chế phản ứng ở mức độ phân tử. Việc ứng dụng hóa học lượng tử và động hóa học tính toán là giải pháp nhằm làm sáng tỏ bề mặt thế năng, thông số nhiệt động và hằng số tốc độ của các phản ứng liên quan đến gốc propargyl.
Mục tiêu nghiên cứu
- Sử dụng phương pháp phiếm hàm mật độ (DFT) và phương pháp tương tác chùm $CCSD(T)$ với bộ hàm cơ sở 6-311++G(3df,2p) để nghiên cứu cơ chế phản ứng và các thông số nhiệt động của các phản ứng trong pha khí.
- Xây dựng bề mặt thế năng (PES) của các phản ứng giữa gốc propargyl với 9 chất thường gặp trong phản ứng cháy: $H$, $H_2$, $O_2$, $OH$, $H_2O$, $NO$, $NH_3$, $CH_3$ và $C_3H_3$. Cấu trúc hình học được tối ưu hóa ở mức B3LYP/6-311++G(3df,2p) và năng lượng điểm đơn được tính ở mức CCSD(T)/6-311++G(3df,2p).
- Tính các thông số nhiệt động học gồm biến thiên entanpi ($\Delta H^\circ_{298}$), biến thiên entropi ($\Delta S^\circ_{298}$) và biến thiên năng lượng tự do Gibbs ($\Delta G^\circ_{298}$) cho các đường phản ứng theo mức lý thuyết CCSD(T)/6-311++G(3df,2p).
- Tính toán động học cho các phản ứng, xác định hằng số tốc độ cho từng sản phẩm và hằng số tốc độ tổng cộng ($k_{tot}$) ở áp suất tiêu chuẩn 760 Torr trong khoảng nhiệt độ từ 300 K đến 2000 K (hoặc điều chỉnh phù hợp với điều kiện thực nghiệm đã công bố).
- Dự đoán và giải thích các kết quả thực nghiệm dựa trên cơ sở tính toán lý thuyết hóa học lượng tử.
Đối tượng và phạm vi nghiên cứu
- Đối tượng nghiên cứu: Bề mặt thế năng, thông số nhiệt động học, cơ chế phản ứng và hằng số tốc độ phản ứng giữa gốc propargyl ($C_3H_3$) với 9 tác nhân: nguyên tử $H$, phân tử $H_2$, phân tử $O_2$, gốc $OH$, phân tử $H_2O$, phân tử $NO$, phân tử $NH_3$, gốc $CH_3$ và gốc $C_3H_3$.
- Phạm vi nghiên cứu: Các phản ứng khảo sát hoàn toàn trong pha khí; áp suất khảo sát chuẩn là 1 atm (760 Torr); khoảng nhiệt độ khảo sát động học từ 300 K đến 2000 K.
Tổng quan tài liệu và vị trí của luận án
Các hướng nghiên cứu trong nước
Luận án điểm lại các công trình nghiên cứu lý thuyết trong nước về cơ chế và động học của các gốc hydrocacbon trong pha khí:
- Trần Quốc Trị (2011): Khảo sát thông số nhiệt động học và đường phản ứng của gốc tự do etinyl ($C_2H$) trong 9 hệ phản ứng với $H_2S$, $PH_3$, $HCN$, $CH_3CN$, $C_2H_5CN$, $C_3H_7CN$, $CH_3OH$, $HCHO$, $HCOOH$ bằng phương pháp hóa học lượng tử.
- Nguyễn Hữu Thọ (2012): Nghiên cứu thông số nhiệt động học và đường phản ứng của gốc tự do metylidin ($CH$) trong 9 hệ phản ứng với $H_2O$, cis-$CH_3OH$, trans-$CH_3OH$, $H_2CO$, cis-$HCOH$, trans-$HCOH$, $HCN$, $CH_3CN$, $C_6H_4$.
- Nguyễn Trọng Nghĩa (2014): Nghiên cứu lý thuyết cơ chế và động học phản ứng của axit fulminic ($HCNO$) với 9 chất ($H$, $F$, $OH$, $CH_3$, $SH$, $NH_2$, $SiH_3$, $NH_2$, $PH_2$), tính hằng số tốc độ phản ứng cho hai hệ $HCNO + H$ và $HCNO + OH$.
Các hướng nghiên cứu quốc tế
- Eugene V. Slagle và Vadim D. Knyazev (2003, Đại học Catholic, Hoa Kỳ): Thực nghiệm phản ứng giữa hai gốc propargyl trong khoảng nhiệt độ 500 - 1000 K, xác định các sản phẩm $C_6H_5 + H$, $C_6H_6$ (mạch thẳng và vòng benzen), $C_8H_6 + CH_3$, $C_9H_8 + H$, $C_6H_6 + C_3H_3$, trong đó benzen là sản phẩm chính; xác định biểu thức hằng số tốc độ $k(T) = 4,49 \cdot 10^{-13} \exp(-128/T) \text{ cm}^3 \text{molecule}^{-1} \text{s}^{-1}$, tại 500 K đạt $(3,3 \pm 0,35) \cdot 10^{-11} \text{ cm}^3 \text{molecule}^{-1} \text{s}^{-1}$ và tại 1000 K đạt $(1,2 \pm 0,14) \cdot 10^{-10} \text{ cm}^3 \text{molecule}^{-1} \text{s}^{-1}$.
- Khoa Hóa, Trường Đại học Seoul, Hàn Quốc (2006): Thực nghiệm phản ứng $C_3H_3 + O(^3P)$ bằng phổ phát xạ tử ngoại chân không kết hợp ab initio, thu được sản phẩm chính là $H(^2S) + C_3H_2O$ (propynal) với năng lượng dịch chuyển trung bình $5,09 \pm 0,36 \text{ kcal/mol}$, cùng các kênh phụ tạo $C_2H_3 + CO$, $C_3H_2 + OH$, $C_2H_2 + CHO$, $C_2H_2O + CH$, $C_2HO + CH_2$, $CH_2O + C_2H$.
- Slagle và Gutman (2012): Nghiên cứu thực nghiệm hệ $C_3H_3 + O_2$ ở dải nhiệt độ 500 - 900 K, đo được các thông số nhiệt động của cân bằng tạo $C_3H_3O_2$ ($\Delta H^\circ_{298} = -79 \pm 6 \text{ kJ/mol}$, $\Delta S^\circ_{298} = -131 \pm 12 \text{ J/mol}\cdot\text{K}$); ở nhiệt độ cao phân hủy tạo $H_2C_2O + HCO$ với hằng số tốc độ $k(T) = 5,0 \cdot 10^{-13} \exp(-E_a/RT) \text{ cm}^3 \text{molecule}^{-1} \text{s}^{-1}$.
- Antonio Fernández-Ramos, James A. Klippenstein và Donald G. Truhlar: Nghiên cứu 9 đường thực nghiệm của hệ $C_3H_3 + H$ dẫn tới 4 sản phẩm.
Khoảng trống nghiên cứu luận án giải quyết
Phần lớn các nghiên cứu trước đây chỉ khảo sát thực nghiệm hoặc tính toán rời rạc trên một vài phản ứng đơn lẻ của gốc propargyl với $O_2$, $O(^3P)$ hoặc $C_3H_3$. Luận án lựa chọn giải quyết khoảng trống về một hệ thống dữ liệu toàn diện, có độ chính xác cao về bề mặt thế năng, cơ chế chi tiết, thông số nhiệt động chuẩn ($\Delta H^\circ_{298}$, $\Delta S^\circ_{298}$, $\Delta G^\circ_{298}$) và động học phụ thuộc nhiệt độ cho 9 hệ phản ứng của gốc propargyl với các tác nhân đại diện trong phản ứng cháy và hóa học khí quyển.
Cơ sở lý thuyết và phương pháp nghiên cứu
Khung lý thuyết sử dụng
- Cơ học lượng tử và hóa học lượng tử: Ứng dụng phương trình Schrödinger cho trạng thái dừng: $\hat{H}\Psi = E\Psi$; xấp xỉ Born - Oppenheimer (tách biệt chuyển động hạt nhân và electron); nguyên lý biến phân và lý thuyết nhiễu loạn.
- Phương pháp phiếm hàm mật độ (DFT): Sử dụng hàm mật độ electron $\rho(x,y,z)$ để xác định năng lượng hệ thay cho hàm sóng $3N$ biến số; sử dụng phiếm hàm lai hóa trao đổi - tương quan B3LYP (kết hợp các phiếm hàm Becke 88, Hartree-Fock, Lee-Yang-Parr và Vosko-Wilk-Nusair).
- Phương pháp tương tác chùm Coupled-Cluster: Phương pháp $CCSD(T)$ (Coupled Cluster với các kích thích đơn, đôi và bổ chính kích thích ba qua nhiễu loạn) dựa trên toán tử hàm mũ $\hat{T} = \hat{T}_1 + \hat{T}_2 + \hat{T}_3$.
- Lý thuyết bề mặt thế năng (PES): Xác định các điểm dừng trên PES thỏa mãn $\partial E/\partial q_i = 0$. Điểm cực tiểu có ma trận Hessian dương ($\partial^2 E/\partial q_i \partial q_j > 0$ với mọi tọa độ $q$); trạng thái chuyển tiếp (TS) là điểm yên ngựa bậc một có đúng một giá trị riêng âm, tương ứng với một tần số dao động ảo dọc theo tọa độ phản ứng nội tại (IRC).
- Lý thuyết động hóa học:
- Định luật Arrhenius: $k(T) = A \exp(-E_a/RT)$.
- Thuyết va chạm lưỡng phân tử (W. Lewis, 1918): $k = P \cdot Z_0 \cdot \exp(-E/RT)$.
- Thuyết phức chất hoạt động (TST - Eyring, Polanyi, 1935): $k_r = (k_B T / h) \cdot C_0^{\Delta n^\ddagger} \exp(\Delta S^\ddagger/R) \exp(-\Delta H^\ddagger/RT)$.
- Thuyết trạng thái chuyển tiếp biến phân (VTST): Tiếp cận vi cổ điển xác định vị trí trạng thái chuyển tiếp tại điểm có tổng trạng thái $G^\ddagger(E, J)$ đạt cực tiểu ($dk(E)/dq^\ddagger = 0$).
- Thuyết Rice-Ramsperger-Kassel-Marcus (RRKM) kết hợp phương trình Master (Master Equation): Tính toán hằng số tốc độ cho các phản ứng đơn phân tử và lưỡng phân tử đa bước qua các mức năng lượng $k(E,J)$ và xác suất va chạm truyền năng lượng $P_{ij}$.
Bảng dữ liệu kiểm chứng lý thuyết VTST trong luận án
Để minh họa hiệu quả của lý thuyết VTST đối với phản ứng không có rào cản năng lượng truyền thống, luận án dẫn số liệu tính toán phản ứng chuẩn $H + CH_3 \rightarrow CH_4$:
| Nhiệt độ $T$ (K) |
Chiều dài liên kết $H-CH_3$ ($R$, Å) |
Tần số dao động $\nu$ ($\text{cm}^{-1}$) |
Hiệu năng lượng $E_0$ (kcal/mol) |
Hằng số tốc độ $k$ ($10^8 \text{ m}^3\text{mol}^{-1}\text{s}^{-1}$) |
| 200 |
3,54 |
95 |
-0,51 |
1,24 |
| 400 |
3,33 |
138 |
-1,05 |
1,56 |
| 600 |
3,19 |
173 |
-1,62 |
1,70 |
| 800 |
3,09 |
204 |
-2,22 |
1,78 |
| 1000 |
3,00 |
233 |
-2,84 |
1,83 |
Quy trình tính toán thực tế và công cụ phần mềm
- Tối ưu hóa cấu trúc hình học và tần số dao động: Mức lý thuyết B3LYP/6-311++G(3df,2p). Áp dụng phương pháp cấu hình không hạn chế (UB3LYP) cho các hệ vỏ mở có electron độc thân; giá trị độ nhiễm spin $\langle S^2 \rangle$ được kiểm soát trong khoảng 0,76 - 0,79 (nằm trong giới hạn chuẩn 0,75 - 1,0 đối với hệ doublet).
- Tính toán năng lượng điểm đơn: Mức lý thuyết CCSD(T)/6-311++G(3df,2p). Năng lượng tổng cộng được hiệu chỉnh với năng lượng dao động điểm không (ZPE) từ mức B3LYP nhân hệ số hiệu chỉnh kinh nghiệm 0,97:
$$E_{total} = E_{CCSD(T)} + 0,97 \cdot ZPE_{B3LYP}$$
- Tính toán nhiệt động học: Tính biến thiên entanpi ($\Delta H^\circ_{298}$), entropi ($\Delta S^\circ_{298}$) và thế đẳng áp Gibbs ($\Delta G^\circ_{298}$) cho từng hướng phản ứng dựa trên hiệu chỉnh nhiệt $H_{corr}$ và $G_{corr}$.
- Tính toán động học: Sử dụng phần mềm Chemrate để tính theo TST và phần mềm Variflex để tính toán theo VTST và RRKM kết hợp giải phương trình Master; biểu diễn đường năng lượng tối thiểu (MEP) không có TS bằng hàm thế Morse: $V(R) = D_e [1 - \exp(-\beta(R - R_e))]^2$.
- Bộ phần mềm tính toán: Gaussian 09, GaussView, ChemCraft, ChemOffice, Chemrate, Variflex.
Nội dung chính theo từng chương
Bố cục cấu trúc của luận án
Luận án có độ dài 146 trang, bao gồm 29 bảng số liệu, 31 hình vẽ, 87 tài liệu tham khảo và 13 trang phụ lục:
- Phần mở đầu: 7 trang.
- Chương 1: Cơ sở lý thuyết (16 trang).
- Chương 2: Tổng quan về gốc propargyl và phương pháp tính (10 trang).
- Chương 3: Kết quả nghiên cứu và thảo luận (103 trang).
- Kết luận chung và kiến nghị: 3 trang.
- Danh mục công trình công bố của tác giả: 1 trang (12 công trình).
- Tài liệu tham khảo: 6 trang.
Mở đầu (7 tr.) -> Chương 1 (16 tr.) -> Chương 2 (10 tr.) -> Chương 3 (103 tr.) -> Kết luận & Kiến nghị (3 tr.)
Chương 1: Cơ sở lý thuyết phương pháp tính hóa học lượng tử, bề mặt thế năng và động hóa học
Chương này trình bày có hệ thống các nguyên lý cơ bản của hóa học lượng tử và động hóa học:
- Phương pháp hàm sóng và DFT: Trình bày chi tiết toán tử Hamilton, phép gần đúng Hartree-Fock (HF), phương pháp tương tác cấu hình (CI), phương pháp tương tác chùm (CC, CCSD, CCSD(T)), các bộ hàm cơ sở Slater (STO) và Gauss (GTO), các phiếm hàm mật độ (BLYP, B3LYP, BHandHLYP).
- Lý thuyết bề mặt thế năng: Định nghĩa toán học về PES $E = f(q_1, q_2, \dots, q_n)$, phân biệt điểm cực tiểu, điểm yên ngựa bậc một (TS) và đường tọa độ phản ứng thực IRC; ý nghĩa của tần số dao động thực và tần số dao động ảo.
- Lý thuyết động học: Thiết lập công thức định luật Arrhenius, thuyết va chạm Lewis, thuyết trạng thái chuyển tiếp Eyring-Polanyi, thuyết trạng thái chuyển tiếp cải biến VTST và lý thuyết phản ứng đơn phân tử RRKM kết hợp giải phương trình vi phân Master Equation.
Chương 2: Tổng quan về gốc propargyl và phương pháp tính
Chương 2 phân tích đặc điểm hóa học của gốc propargyl và quy chuẩn tính toán:
- Đặc tính của gốc propargyl ($C_3H_3$): Là gốc hydrocacbon nhỏ với 3 nguyên tử C và 3 nguyên tử H. Tồn tại 2 dạng đồng phân hỗ biến, trong đó dạng cấu trúc thứ nhất ($H_2C=C=\dot{C}H \leftrightarrow H_2\dot{C}-C\equiv CH$) hoạt động hơn và chiếm tỷ lệ 65% trong pha khí. Gốc có hai trung tâm phản ứng chính: nguyên tử carbon chứa electron độc thân và hệ liên kết $\pi$.
- Cơ chế phản ứng chính của gốc cacbo tự do: Phản ứng cộng hợp, phản ứng tách, phản ứng chèn, phản ứng chuyển vị H nội phân tử hoặc gian phân tử, phản ứng cắt chuyển H và phản ứng cắt mạch C-C.
- Nguồn gốc hình thành: Phân tích các kênh hình thành gốc propargyl từ phản ứng nhiệt phân, quang phân etilen, axetilen hoặc phản ứng giữa nguyên tử Clo với propin ($Cl + C_3H_4 \rightarrow C_3H_3 + HCl$) sau khi quang phân $(CClO)_2$.
- Chiến lược tính toán: Thiết lập cấu hình tính toán CCSD(T)//B3LYP/6-311++G(3df,2p), áp dụng phương pháp phiếm hàm mật độ không hạn chế UB3LYP nhằm loại trừ sai số do nhiễm spin, tính toán nhiệt động ở 298 K và động học trong dải 300 - 2000 K ở 1 atm.
Chương 3: Kết quả nghiên cứu và thảo luận
Chương 3 là phần thực nghiệm tính toán trọng tâm của luận án (103 trang), khảo sát toàn diện 9 hệ phản ứng giữa gốc propargyl với $H$, $H_2$, $O_2$, $OH$, $H_2O$, $NO$, $NH_3$, $CH_3$ và $C_3H_3$.
1. Hệ phản ứng giữa gốc propargyl và nguyên tử Hydro ($C_3H_3 + H$)
Tác giả khảo sát chi tiết PES ở trạng thái spin singlet, phát hiện 6 trạng thái trung gian ($I_1$ đến $I_6$), 16 trạng thái chuyển tiếp (TS) và 11 kênh sản phẩm ($P_1$ đến $P_{11}$):
| Ký hiệu |
Cấu trúc sản phẩm |
Năng lượng tương quan (kcal/mol) |
Đặc điểm cơ chế và nhiệt động |
| $P_1$ |
$C_3H_4$ (propyne/allene) |
- |
Tạo thành từ $I_1, I_2, I_3$ qua các rào cản TS tương ứng |
| $P_2$ |
$C_2H_2 + CH_2$ |
+12,41 |
Đi từ $I_3$ qua trạng thái chuyển tiếp $T_{3/P2}$ (20,95 kcal/mol) |
| $P_3$ |
$CH_3 + CCH$ |
+35,35 |
Hình thành trực tiếp từ $I_2$ không qua TS, phản ứng thu nhiệt |
| $P_4$ |
$H_2CCC + H_2$ |
-3,22 |
Đi từ $I_1$ qua $T_{1/P4}$ (5,93 kcal/mol), thuận lợi nhất về mặt năng lượng |
| $P_5$ |
cyc-$C_3H_2 + H_2$ |
-16,72 |
Đi từ $I_5$ qua $T_{5/P5}$ (19,86 kcal/mol), bền vững nhất về mặt nhiệt động |
| $P_6$ |
$C_2H_3 + CH$ |
- |
Đi từ $I_4$ bằng cách cắt đứt mạch hoặc trực tiếp từ $I_3$ |
| $P_7$ |
$H_2CC + CH_2$ |
- |
Đi từ $I_4$ bằng phản ứng cắt mạch |
| $P_8$ |
$H_3CC + CH$ |
+124,01 |
Tạo trực tiếp từ $I_2$, mức năng lượng cao nhất, khó tạo thành nhất |
| $P_9$ |
cyc-$C_3H_3 + H$ |
+40,31 |
Tạo trực tiếp từ $I_6$ bằng sự tách H; H đóng vai trò chất xúc tác |
| $P_{10}$ |
cyc-$C_3H_3 + H$ |
- |
Tạo trực tiếp từ $I_5$ bằng sự tách H; H đóng vai trò chất xúc tác |
| $P_{11}$ |
$HCCCH + H_2$ |
-5,36 |
Phản ứng tách H trực tiếp qua $T_{0/P11}$ (12,34 kcal/mol), $\Delta H^\circ_{298} = -10,95 \text{ kcal/mol}$ |
- Phản ứng tách H trực tiếp: Nguyên tử H tấn công vào đầu $C_1$ của propargyl qua $T_{0/P11}$ (hàng rào 12,34 kcal/mol). Tại $T_{0/P11}$, khoảng cách liên kết $C_1-H$ bị kéo dài ra 1,30 Å, khoảng cách $H\cdots H$ hình thành là 0,96 Å (rút ngắn về 0,74 Å trong $P_{11}$). Góc liên kết $HC_1C_2$ mở rộng từ 120,87° lên 128,41°, góc $HC_3C_2$ giảm từ 180° xuống 128,29°, góc $C_1C_2C_3$ uốn cong về 164,61°.
- Phản ứng cộng hợp không có rào cản:
- Cộng vào $C_3$ tạo $I_1$ ($H_2C=C=CH_2$, propadiene): năng lượng tương quan -86,48 kcal/mol.
- Cộng vào $C_1$ tạo $I_2$ ($H_3C-C\equiv CH$, propyne): năng lượng tương quan -86,46 kcal/mol (-87,46 kcal/mol).
- Cộng vào $C_2$ tạo $I_3$ ($H_2C=CH\dot{C}H$): năng lượng tương quan -26,37 kcal/mol (kém bền hơn do mất tính đối xứng và còn 2 electron chưa ghép đôi).
- So sánh với thực nghiệm: Nghiên cứu thực nghiệm trước đây của Antonio Fernández-Ramos và cộng sự chỉ phát hiện 4 sản phẩm ($P_1, P_4, P_5, P_{11}$). Luận án đã xác định thêm 7 hướng phản ứng tạo sản phẩm ($P_2, P_3, P_6, P_7, P_8, P_9, P_{10}$) và giải thích tính chất kém bền nhiệt động của các sản phẩm này.
2. Khảo sát các hệ phản ứng còn lại
Luận án thực hiện tính toán PES, thông số nhiệt động và hằng số tốc độ cho 8 hệ phản ứng tiếp theo:
- Hệ $C_3H_3 + H_2$: Khảo sát các đường phản ứng cộng và tách H, lập bảng thông số nhiệt động và đồ thị hằng số tốc độ theo nhiệt độ.
- Hệ $C_3H_3 + O_2$: So sánh $\Delta H^\circ_{298}$ của 11 đường phản ứng, xác định các giá trị $\Delta S^\circ_{298}$ và $\Delta G^\circ_{298}$, phân tích kênh tạo $H_2C_2O + HCO$ và hằng số tốc độ tổng cộng trong dải 300 - 2000 K.
- Hệ $C_3H_3 + OH$: So sánh $\Delta H^\circ_{298}$ của 8 đường phản ứng, tính toán hằng số tốc độ tạo thành các chất trung gian $IS_1, IS_2$ và các sản phẩm tách nước, cộng nhóm hydroxyl.
- Hệ $C_3H_3 + H_2O$: So sánh $\Delta H^\circ_{298}$ của 11 đường phản ứng, phân tích cơ chế cộng liên kết O-H của nước vào liên kết $\pi$ của propargyl tạo $H_2CCHCHOH$ hoặc $H_3CCCHOH$.
- Hệ $C_3H_3 + NO$: So sánh $\Delta H^\circ_{298}$ của 14 đường phản ứng, xác định hằng số tốc độ phản ứng $C_3H_3 + NO \rightarrow IS_1$ và các sản phẩm phân hủy chứa nitơ.
- Hệ $C_3H_3 + NH_3$: Xác định các thông số nhiệt động ở 298 K, xây dựng bề mặt thế năng rút gọn và tính hằng số tốc độ $k(T)$.
- Hệ $C_3H_3 + CH_3$: So sánh $\Delta H^\circ_{298}$ của 16 đường phản ứng, tính biến thiên $\Delta S^\circ_{298}$, $\Delta G^\circ_{298}$ và hằng số tốc độ tạo $IS_1, IS_2$.
- Hệ $C_3H_3 + C_3H_3$: Khảo sát phản ứng tái kết hợp của hai gốc propargyl tạo các trạng thái trung gian mạch thẳng và đóng vòng benzen ($C_6H_6$), gốc phenyl ($C_6H_5 + H$), tính hằng số tốc độ $k(T)$ và so sánh với biểu thức thực nghiệm của Slagle và Knyazev.
Kết quả và những đóng góp mới
- Xây dựng bề mặt thế năng toàn diện: Đã thiết lập chi tiết bề mặt thế năng (PES) cho 9 hệ phản ứng của gốc propargyl với 9 tác nhân trong pha khí ($H$, $H_2$, $O_2$, $OH$, $H_2O$, $NO$, $NH_3$, $CH_3$, $C_3H_3$) ở mức lý thuyết $CCSD(T)//B3LYP/6-311++G(3df,2p)$.
- Xác định hệ thống thông số nhiệt động học chuẩn: Đã tính toán đầy đủ các giá trị $\Delta H^\circ_{298}$, $\Delta S^\circ_{298}$ và $\Delta G^\circ_{298}$ cho tất cả các đường phản ứng trong cả 9 hệ; so sánh, đối chiếu hiệu ứng nhiệt tính toán với các số liệu thực nghiệm đã công bố của nhiều nhóm tác giả quốc tế.
- Tính toán động hóa học và thiết lập phương trình Arrhenius: Đã tính hằng số tốc độ của các đường phản ứng ưu tiên và hằng số tốc độ tổng cộng ($k_{tot}$) của từng hệ phản ứng trong khoảng nhiệt độ từ 300 K đến 2000 K (ở áp suất chuẩn 760 Torr) dựa trên lý thuyết TST, VTST và RRKM kết hợp giải phương trình Master; xây dựng các biểu thức giải tích mô tả sự phụ thuộc của hằng số tốc độ vào nhiệt độ dạng $k(T) = A \cdot T^n \exp(-E_a/RT)$.
- Làm sáng tỏ cơ chế phản ứng ở mức độ phân tử: Giải thích cơ chế chuyển dịch electron, sự chuyển vị hydro, sự cạnh tranh giữa các kênh phản ứng cộng hợp và phản ứng tách, đồng thời chỉ ra các kênh sản phẩm mới mà thực nghiệm chưa phát hiện do rào cản năng lượng cao.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế được ghi nhận từ phạm vi nghiên cứu
- Các tính toán động học và nhiệt động học trong luận án được thực hiện hoàn toàn đối với các phản ứng đơn phân tử và lưỡng phân tử trong pha khí, chưa khảo sát ảnh hưởng của hiệu ứng dung môi hoặc pha lỏng.
- Dải áp suất khảo sát chủ yếu tập trung tại áp suất thường (760 Torr / 1 atm); chưa mở rộng khảo sát ở các dải áp suất cực cao hoặc áp suất siêu thấp trong điều kiện đặc thù của thiên văn học liên sao.
Hướng nghiên cứu tiếp theo
- Mở rộng phạm vi nghiên cứu lý thuyết đối với phản ứng của gốc propargyl với các phân tử chứa các dị tố khác như lưu huỳnh ($S$), phospho ($P$), hoặc các gốc hydrocacbon chưa bão hòa phức tạp hơn.
- Khảo sát sâu hơn cơ chế động học hình thành các vòng thơm đa vòng (PAHs) bậc cao từ các sản phẩm trung gian $C_6H_5$ và $C_6H_6$ trong môi trường nhiệt độ và áp suất biến đổi liên tục của buồng đốt động cơ.
Giá trị tham khảo
- Đối tượng tham khảo:
- Nghiên cứu sinh, học viên cao học và giảng viên: Thuộc các chuyên ngành Hóa lý thuyết và Hóa lý, Hóa học tính toán, Hóa học môi trường, Kỹ thuật Hóa học.
- Kỹ sư và nhà nghiên cứu công nghệ năng lượng: Các nhóm nghiên cứu về cơ chế cháy, mô phỏng phản ứng đốt nhiên liệu hóa thạch, cơ chế giảm thiểu phát thải khí nhà kính và kiểm soát bụi mịn/muội than.
- Các nội dung có giá trị tham khảo chuyên sâu:
- Quy trình phối hợp đa phương pháp tính toán lượng tử: tối ưu hóa cấu trúc ở mức B3LYP/6-311++G(3df,2p) và hiệu chỉnh năng lượng chính xác cao ở mức CCSD(T)/6-311++G(3df,2p).
- Phương pháp xử lý hệ vỏ mở và kiểm soát độ nhiễm spin trong các gốc tự do.
- Phương pháp ứng dụng phần mềm Chemrate và Variflex để tính toán hằng số tốc độ vi cổ điển theo lý thuyết VTST và RRKM giải phương trình Master.
Câu hỏi thường gặp
1. Luận án sử dụng tổ hợp phương pháp và bộ hàm cơ sở nào để nghiên cứu cấu trúc và năng lượng của các chất?
Tác giả sử dụng phương pháp phiếm hàm mật độ không hạn chế UB3LYP kết hợp bộ hàm cơ sở 6-311++G(3df,2p) để tối ưu hóa cấu trúc hình học và tính tần số dao động. Năng lượng điểm đơn được tính ở mức tương tác chùm cấp cao CCSD(T) với cùng bộ hàm cơ sở 6-311++G(3df,2p) và được hiệu chỉnh năng lượng dao động điểm không (ZPE) với hệ số kinh nghiệm 0,97.
2. Trong hệ phản ứng $C_3H_3 + H$, sản phẩm nào bền vững nhất về mặt nhiệt động và sản phẩm nào thuận lợi nhất về mặt năng lượng?
- Bền vững nhất về mặt nhiệt động: Sản phẩm $P_5$ (cyc-$C_3H_2 + H_2$) có năng lượng tương quan thấp nhất toàn hệ (-16,72 kcal/mol).
- Thuận lợi nhất về mặt năng lượng (động học): Sản phẩm $P_4$ ($H_2CCC + H_2$) hình thành từ chất trung gian $I_1$ qua trạng thái chuyển tiếp $T_{1/P4}$ có hàng rào năng lượng thấp nhất chỉ 5,93 kcal/mol.
3. Phương pháp thực nghiệm nào được tác giả tổng quan để điều chế gốc propargyl trong thực tế?
Gốc propargyl trong thực nghiệm được tạo ra bằng phương pháp quang phân (photolysis): dùng tia laser có bước sóng 248 nm để phân tách oxalyl chloride $(CClO)_2$ thành nguyên tử Clo và khí $CO$, sau đó nguyên tử Clo phản ứng với phân tử propin ($C_3H_4$) để tạo ra gốc propargyl ($C_3H_3$) và phân tử $HCl$.
4. Việc tính toán hằng số tốc độ phản ứng theo các lý thuyết động học trong luận án sử dụng các công cụ phần mềm nào?
Tác giả sử dụng phần mềm Gaussian 09 để tính toán cấu trúc lượng tử và bề mặt thế năng; phần mềm Chemrate để tính hằng số tốc độ theo thuyết trạng thái chuyển tiếp (TST); và phần mềm Variflex để giải phương trình Master, tính toán hằng số tốc độ theo thuyết trạng thái chuyển tiếp cải biến (VTST) và thuyết RRKM.
5. So với nghiên cứu thực nghiệm của Antonio Fernández-Ramos và cộng sự trên hệ $C_3H_3 + H$, luận án đã bổ sung những kết quả mới nào?
Nghiên cứu thực nghiệm trước đây chỉ phát hiện 4 kênh sản phẩm ($P_1, P_4, P_5, P_{11}$). Luận án đã xây dựng bề mặt thế năng chi tiết xác định toàn bộ 11 kênh sản phẩm, bổ sung thêm 7 sản phẩm ($P_2, P_3, P_6, P_7, P_8, P_9, P_{10}$) và giải thích nguyên nhân các sản phẩm này khó quan sát được trong thực nghiệm do năng lượng rào cản cao và kém bền nhiệt động.
Kết luận
Luận án của tác giả Phạm Văn Tiến đã hoàn thành việc nghiên cứu lý thuyết có hệ thống về cơ chế phản ứng, nhiệt động học và động học của gốc propargyl ($C_3H_3$) với 9 tác nhân trong pha khí bằng tổ hợp phương pháp CCSD(T)//B3LYP/6-311++G(3df,2p). Công trình đã thiết lập các bề mặt thế năng chi tiết, xác định cấu trúc hình học của các chất trung gian và trạng thái chuyển tiếp, tính toán các đại lượng $\Delta H^\circ_{298}$, $\Delta S^\circ_{298}$, $\Delta G^\circ_{298}$ và hằng số tốc độ phản ứng $k(T)$ trong dải nhiệt độ 300 - 2000 K theo các lý thuyết TST, VTST và RRKM. Kết quả tính toán cung cấp cơ sở dữ liệu định lượng phục vụ mô phỏng quá trình cháy của nhiên liệu hóa thạch và hóa học khí quyển, được công bố qua 12 bài báo trên các tạp chí chuyên ngành.
