Bối cảnh và vấn đề nghiên cứu

Trong nghiên cứu hóa học hữu cơ, phản ứng cộng electrophile ($A_E$) vào liên kết đôi C=C của alkene là một trong những chuyển hóa cơ bản đã được khảo sát thực nghiệm qua nhiều công trình. Tuy nhiên, các phương pháp thực nghiệm truyền thống gặp nhiều giới hạn khi cần quan sát trực tiếp các cấu trúc vi mô có thời gian tồn tại rất ngắn như trạng thái chuyển tiếp (transition state - TS) hay các tiểu phân trung gian (carbocation). Sự phát triển của hóa học lượng tử và các phần mềm tính toán hiện đại như Gaussian và GaussView cho phép mô phỏng bề mặt thế năng, tối ưu hóa cấu trúc hình học và tính toán các thông số nhiệt động - động học nhằm kiểm chứng, làm sáng tỏ bản chất của các quy luật thực nghiệm.

Khóa luận tốt nghiệp ngành Sư phạm Hóa học của tác giả Võ Thị Mỹ Huyền (Trường Đại học An Giang, thực hiện dưới sự hướng dẫn của TS. Phạm Phát Tân, bảo vệ tháng 05/2022) tập trung khảo sát đề tài: "Ảnh hưởng của cấu trúc đến hoạt tính hoá học của alkene trong phản ứng cộng electrophile".

Mục tiêu và nhiệm vụ nghiên cứu của khóa luận được xác định cụ thể như sau:

  1. Nghiên cứu sự ảnh hưởng cấu trúc của chất nền (các alkene) và tác nhân electrophile đến hoạt tính hóa học trong phản ứng cộng $A_E$.
  2. Xây dựng cấu trúc ban đầu của các phân tử alkene, tác nhân electrophile, trạng thái chuyển tiếp và sản phẩm trung gian trên phần mềm GaussView 6.0.
  3. Thực hiện tính toán lượng tử và xử lý số liệu bằng phần mềm Gaussian 09W.
  4. Lập bảng so sánh, thống kê độ dài liên kết, năng lượng liên kết, mật độ điện tích nguyên tử, mức năng lượng orbital phân tử và các đại lượng nhiệt động ($\Delta H$, $\Delta S$, $\Delta G$, $E_a$, $k$) để mô tả, giải thích và đối chiếu với thực nghiệm.

Đối tượng nghiên cứu là phản ứng cộng của alkene với tác nhân electrophile. Phạm vi nghiên cứu giới hạn ở mô hình tính toán hóa học lượng tử cho các alkene mạch ngắn thể khí gồm propene ($CH_2=CH-CH_3$), but-1-ene ($CH_2=CH-CH_2-CH_3$) và 2-methylpropene ($CH_2=C(CH_3)_2$); các tác nhân phản ứng được khảo sát là $HBr$ và $H_2O$ trong điều kiện pha khí (không sử dụng mô hình dung môi). Nghiên cứu không khảo sát tác nhân $HI$ do hạn chế của hệ hàm cơ sở 6-31G(d).

Cơ sở lý thuyết và phương pháp

Khóa luận vận dụng kết hợp các cơ sở lý thuyết hóa hữu cơ kinh điển và lý thuyết hóa học lượng tử hiện đại:

  • Cơ chế phản ứng cộng electrophile ($A_E$): Cơ chế hai giai đoạn nối tiếp. Giai đoạn 1 là giai đoạn chậm, quyết định tốc độ phản ứng, trong đó liên kết $\pi$ của alkene tương tác với tác nhân electrophile ($E^+$) để tạo thành carbocation trung gian. Giai đoạn 2 diễn ra nhanh khi tác nhân nucleophile ($Y^-$) tấn công vào carbocation phẳng để tạo thành sản phẩm no. Động học phản ứng tuân theo phương trình bậc hai: $v = k[alkene][EY]$.
  • Tiên đề Hammond: Sử dụng để dự đoán cấu trúc của trạng thái chuyển tiếp. Đối với giai đoạn thu nhiệt và có năng lượng hoạt hóa cao, cấu trúc của trạng thái chuyển tiếp có nhiều nét tương đồng với sản phẩm trung gian (carbocation).
  • Thuyết Orbital phân tử biên (Frontier Molecular Orbital - FMO): Khung lý thuyết của Kenichi Fukui (Nobel Hóa học 1982) giải thích tương tác giữa orbital phân tử bị chiếm có mức năng lượng cao nhất (HOMO) của tác chất cho electron và orbital phân tử trống có mức năng lượng thấp nhất (LUMO) của tác chất nhận electron. Cặp tương tác nào có độ chênh lệch năng lượng ($\Delta E = |E_{HOMO} - E_{LUMO}|$) nhỏ hơn sẽ là tương tác ưu đãi quyết định chiều hướng phản ứng.

Phương pháp nghiên cứu tính toán cụ thể gồm:

  • Mức lý thuyết tính toán: Phương pháp phiếm hàm mật độ (Density Functional Theory - DFT) sử dụng phiếm hàm lai hóa ba thông số Becke kết hợp hiệu chỉnh tương quan Lee-Yang-Parr (B3LYP).
  • Hệ hàm cơ sở: Hệ hàm phân cực 6-31G(d), bổ sung orbital hàm $d$ cho các nguyên tử nặng nhằm mô tả độ biến dạng không gian của đám mây electron.
  • Quy trình tối ưu hóa cấu trúc: Dùng từ khóa Opt với 4 tiêu chuẩn hội tụ: Maximum Force (ngưỡng 0.000450), RMS Force (ngưỡng 0.000300), Maximum Displacement (ngưỡng 0.001800), và RMS Displacement (ngưỡng 0.001200).
  • Tính toán tần số dao động: Dùng từ khóa Freq tại cùng mức lý thuyết B3LYP/6-31G(d) để xác định năng lượng dao động điểm không (Zero-Point Vibrational Energy - $E_{ZPV}$) và kiểm tra các điểm dừng trên bề mặt thế năng. Trạng thái chuyển tiếp (TS) được xác nhận khi có duy nhất một tần số dao động ảo ($i$).
  • Tìm kiếm trạng thái chuyển tiếp: Ứng dụng phương pháp STQN với từ khóa Opt=QST2 hoặc Opt=(calcfc, ts) kết hợp đọc ma trận lực Opt=ReadFC.
  • Phân tích mật độ điện tích và orbital: Sử dụng lệnh phân tích quần cư Pop=Full để trích xuất hệ số orbital và điện tích Mulliken trên từng nguyên tử.
  • Xử lý số liệu nhiệt động: Năng lượng tổng cộng được xác định theo biểu thức $E_{tổng} = E_{SCF} + 0.9804 \times E_{ZPV}$ (với 0.9804 là hệ số hiệu chỉnh kinh nghiệm cho mức B3LYP/6-31G(d)). Hệ số quy đổi năng lượng sử dụng: $1\text{ Hartree} = 627.509391\text{ kcal/mol} = 2622.989254\text{ kJ/mol}$.

Nội dung chính theo từng chương

Chương 1: Giới thiệu

Chương mở đầu trình bày tính cấp thiết của việc áp dụng hóa học lượng tử tính toán trong việc giải thích cấu trúc và cơ chế phản ứng hữu cơ, nêu rõ vai trò của bộ công cụ Gaussian 09W và GaussView 6.0. Tác giả xác lập mục đích nghiên cứu, đối tượng, phạm vi khảo sát phản ứng cộng $A_E$ trên hệ alkene $C_3 - C_4$, đồng thời nêu rõ các đóng góp học thuật và sư phạm của đề tài.

Chương 2: Tổng quan về đề tài nghiên cứu

Chương 2 hệ thống hóa cơ sở lý thuyết tổng quan về phản ứng cộng, phân loại các phản ứng cộng ($A_R, A_E, A_{Nu}$, cộng đồng bộ) và các loại tác nhân electrophile. Tác giả phân tích chi tiết cơ chế phản ứng cộng $A_E$ vào alkene qua hai giai đoạn, giản đồ năng lượng, động học bậc hai, và hóa học lập thể (tính chọn lọc lập thể dạng cộng anti ưu tiên hơn syn qua các bảng số liệu thực nghiệm của phản ứng halogen hóa). Chương này cũng tổng hợp số liệu thực nghiệm về tốc độ cộng $Br_2$ tương đối của một số alkene và trình bày nguyên lý của thuyết Orbital phân tử biên (FMO), nhấn mạnh tương tác 2 điện tử giữa HOMO và LUMO.

Chương 3: Phương pháp nghiên cứu

Nội dung chương 3 mô tả chi tiết chức năng của Gaussian 09W và GaussView 6.0. Tác giả trình bày bản chất toán học của phương pháp Hartree-Fock (HF), lý thuyết phiếm hàm mật độ (DFT-B3LYP), các hệ hàm cơ sở từ tối thiểu (STO-nG) đến phân cực (6-31G(d)) và khuếch tán. Ngoài ra, chương này giới thiệu tổng quan các mô hình dung dịch liên tục (Onsager, PCM, IPCM, SCI-PCM), các bước tối ưu hóa hình học, thuật toán tìm kiếm trạng thái chuyển tiếp (QST2, CalcFC, ReadFC), quy tắc xác nhận tần số ảo, bảng hệ số hiệu chỉnh năng lượng $E_{ZPV}$ và các công thức tính toán biến thiên entanpy ($\Delta H^\circ_{298}$), entropy ($\Delta S^\circ_{298}$), năng lượng tự do Gibbs ($\Delta G^\circ_{298}$), năng lượng hoạt hóa ($E_a$) và hằng số tốc độ ($k$).

Chương 4: Kết quả và thảo luận

Chương 4 là phần nội dung thực nghiệm tính toán trọng tâm của khóa luận, gồm 8 mục nhỏ:

  1. Khảo sát trạng thái bền của các alkene và HX: Tối ưu hóa cấu trúc propene, but-1-ene, $HBr$ và $H_2O$ ở pha khí.
  2. Nhận định khả năng phản ứng qua mật độ điện tích của phân tử: Phân tích điện tích Mulliken trên propene ($C_1 = -0.333; C_2 = -0.041$) và but-1-ene ($C_2 = -0.037$) cho thấy $C_2$ dương điện hơn, làm sáng tỏ tính bất đối xứng điện tử. Phân tích hình thái HOMO-LUMO cho thấy mật độ điện tử tập trung lớn ở liên kết đôi $C=C$.
  3. Khảo sát phản ứng cộng $A_E$ giữa propene và HBr: Bảng so sánh hệ số orbital chỉ ra tương tác ưu đãi diễn ra giữa HOMO($C_3H_6$) và LUMO($HBr$) do có độ chênh lệch năng lượng thấp hơn (0.24335 so với 0.27286). Bảng tính toán nhiệt động học cho thấy phản ứng tổng quát tỏa nhiệt và tự diễn biến:
Hợp chất / Đại lượng Entanpy (Hartree/Particle) Entropy ($\text{cal}\cdot\text{mol}^{-1}\cdot\text{K}^{-1}$) Năng lượng tự do ($\text{cal}\cdot\text{mol}^{-1}$)
$C_3H_6$ -117.822455 63.189 -73953.5368
$HBr$ -2572.292071 47.457 -1614151.58
$C_3H_7Br$ -2690.151771 75.167 -1688117.91
Biến thiên phản ứng -23.37158727 kcal/mol -35.479 cal/(mol.K) -12.7935234 kcal/mol
  1. Trạng thái chuyển tiếp giai đoạn 1: Tối ưu hóa thành công cấu trúc $TS(1)$ với liên kết $C=C$ dãn dài từ $1.33334\text{ Å}$ lên $1.41106\text{ Å}$, tần số dao động ảo duy nhất đạt $\nu = -1106.71i\text{ cm}^{-1}$, góc liên kết $\widehat{HCH} = 115.41039^\circ$, $\widehat{HCC} = 117.23722^\circ$. Cấu trúc trung gian carbocation propyl hình thành có dạng tam giác phẳng lai hóa $sp^2$, độ dài liên kết $C-C$ là $1.49268\text{ Å}$, các góc $\widehat{HCC} = 118.43649^\circ$ và $\widehat{CCC} = 121.32246^\circ$, điện tích tại $C$ trung tâm là $-0.090$.
  2. Các giá trị nhiệt động trong giai đoạn 1: Giai đoạn 1 có $\Delta H_1 = +48.011371\text{ kcal/mol}$ (thu nhiệt), $\Delta S_1 = -0.3600\text{ cal/(mol.K)}$, $\Delta G_1 = +48.118470\text{ kcal/mol}$ và năng lượng hoạt hóa $E_{a(1)} = 80.89029873\text{ kcal/mol}$. Phù hợp với tiên đề Hammond, cấu trúc của $TS(1)$ nghiêng hẳn về phía carbocation trung gian.
  3. Khảo sát giai đoạn 2: Giai đoạn phản ứng giữa carbocation propyl và $Br^-$. Tương tác orbital chi phối là giữa HOMO($C_3H_7^+$) và LUMO($Br^-$) với $\Delta E = 0.04486$. Cấu trúc phẳng $sp^2$ cho phép $Br^-$ tấn công đối xứng từ hai phía.
  4. Ảnh hưởng của gốc hydrocarbon đến phản ứng cộng $A_E$: So sánh cấu trúc 3 alkene cho thấy chiều dài liên kết $C=C$ tăng dần:
Phân tử Độ dài liên kết $C=C$ ($\text{Å}$)
$CH_2=CH-CH_3$ 1.33334
$CH_2=CH-CH_2-CH_3$ 1.33346
$CH_2=C(CH_3)_2$ 1.33662

Liên kết $C=C$ càng dài thì độ phân cực càng lớn, làm tăng khả năng phản ứng. Mật độ điện tích dương trên nguyên tử carbon trung tâm của các carbocation tương ứng: $CH_3-C^+H-CH_3$ ($-0.090$), $CH_3-C^+H-CH_2-CH_3$ ($-0.089$), và $(CH_3)_3C^+$ ($+0.141$). Carbocation tert-butyl có điện tích dương lớn nhất nhưng do hiệu ứng cản trở không gian của 3 nhóm methyl nên mức chênh lệch orbital trong giai đoạn 2 lớn hơn ($\Delta E = 0.05705$ so với $0.04486$ của propyl). 8. Ảnh hưởng của tác nhân electrophile: Phân tích mức năng lượng HOMO-LUMO chỉ ra trong giai đoạn 1 của phản ứng với $HBr$ hay $H_2O$, tương tác chủ đạo luôn là HOMO của alkene tấn công vào LUMO của tác nhân mang proton $H^+$.

Chương 5: Kết luận và kiến nghị

Tóm tắt các thành tựu tính toán hình học, nhiệt động học, động học và phân tích FMO; đồng thời nêu các kiến nghị mở rộng về hệ hàm, đối tượng alkene lớn hơn, và khảo sát hiệu ứng dung môi.

Kết quả và đóng góp

Khóa luận đã trích xuất được hệ thống dữ liệu lượng tử chi tiết, làm sáng tỏ cơ chế phản ứng cộng electrophile:

  • Thông số hình học và tần số trạng thái chuyển tiếp: Xác định chính xác cấu trúc trạng thái chuyển tiếp giai đoạn 1 ($TS1$) trong phản ứng giữa propene và $HBr$ với tần số ảo duy nhất $\nu = -1106.71i\text{ cm}^{-1}$. Liên kết $C=C$ dãn dài từ $1.33334\text{ Å}$ ở chất phản ứng lên $1.41106\text{ Å}$ ở $TS1$, và đạt $1.49268\text{ Å}$ ở carbocation trung gian. Góc liên kết tại carbocation chuyển dịch về cấu trúc tam giác phẳng $sp^2$ ($\widehat{HCC} = 118.43649^\circ$, $\widehat{CCC} = 121.32246^\circ$).
  • Số liệu nhiệt động và động học:
    • Phản ứng tổng thể giữa propene và $HBr$: $\Delta H^\circ_{298} = -23.37158727\text{ kcal/mol}$ (tỏa nhiệt); $\Delta S^\circ_{298} = -35.479\text{ cal/(mol.K)}$; $\Delta G^\circ_{298} = -12.7935234\text{ kcal/mol}$ (phản ứng tự diễn biến một chiều).
    • Giai đoạn 1 (tạo carbocation): $\Delta H_1 = +48.011371\text{ kcal/mol}$ (thu nhiệt); $\Delta G_1 = +48.118470\text{ kcal/mol}$; năng lượng hoạt hóa $E_{a(1)} = 80.89029873\text{ kcal/mol}$. Năng lượng hoạt hóa cao chứng minh giai đoạn 1 quyết định tốc độ phản ứng.
  • Quy luật cấu trúc - hoạt tính: Độ dài liên kết $C=C$ tăng dần theo mức độ thế alkyl: propene ($1.33334\text{ Å}$) < but-1-ene ($1.33346\text{ Å}$) < 2-methylpropene ($1.33662\text{ Å}$), chứng minh nhóm đẩy electron làm tăng độ phân cực liên kết đôi và làm giảm năng lượng hoạt hóa.
  • Đóng góp của đề tài: Ứng dụng thành công phần mềm Gaussian 09W và GaussView 6.0 vào mô phỏng cơ chế phản ứng hữu cơ ở bậc đại học; cung cấp dữ liệu vi mô định lượng để giải thích quy tắc cộng Markovnikov và tiên đề Hammond; xây dựng nguồn tài liệu tham khảo có giá trị lý thuyết phục vụ giảng dạy chuyên đề hóa hữu cơ và bồi dưỡng học sinh giỏi hóa học ở trường phổ thông.

Hạn chế và hướng nghiên cứu tiếp

Khóa luận tự chỉ ra một số hạn chế và định hướng phát triển:

  • Hạn chế:
    • Chỉ khảo sát ở mức lý thuyết B3LYP với hệ hàm cơ sở trung bình 6-31G(d), chưa sử dụng các hệ hàm khuếch tán hoặc hệ hàm cao hơn để nâng cao độ chuẩn xác tuyệt đối.
    • Chưa tính toán được cho tác nhân $HI$ do hệ hàm 6-31G(d) không hỗ trợ tối ưu hóa nguyên tử iod mà cần dùng hệ hàm giả thế lõi như LanL2DZ.
    • Toàn bộ mô hình phản ứng mới thực hiện ở pha khí, chưa xét đến ảnh hưởng của hiệu ứng solvat hóa trong dung môi phân cực.
  • Hướng nghiên cứu tiếp:
    • Nâng cấp hệ hàm cơ sở cao hơn và áp dụng các mô hình dung môi liên tục (như PCM, SMD).
    • Mở rộng khảo sát cho các alkene có mạch carbon dài hơn ($>C_4$), các dẫn xuất phức tạp hơn và các tác nhân electrophile khác ($Cl_2, Br_2, BrI$, hệ xúc tác acid).
    • Khảo sát thêm các cơ chế cạnh tranh khác như phản ứng cộng nucleophile ($A_{Nu}$) hoặc phản ứng tách ($E1, E2$).

Giá trị tham khảo

Bản khóa luận là tài liệu tham khảo hữu ích cho:

  • Sinh viên và giảng viên ngành Sư phạm Hóa học, Hóa hữu cơ, Hóa lý thuyết: Cung cấp quy trình chi tiết thiết lập các file đầu vào (Job Entry, Route Section), cách sử dụng các từ khóa (Opt=QST2, calcfc, ts, Pop=Full, Freq) và phương pháp trích xuất dữ liệu lượng tử từ file ngõ ra (.out/.log) của Gaussian.
  • Giáo viên trường THPT chuyên: Cung cấp cơ sở lý thuyết lượng tử (mật độ điện tích Mulliken, tương tác FMO, tiên đề Hammond) để giải thích sâu sắc bản chất cơ chế phản ứng cộng vào hợp chất không no trong công tác bồi dưỡng học sinh giỏi quốc gia, quốc tế.
  • Nội dung đáng tham khảo nhất: Quy trình tìm kiếm và xác thực trạng thái chuyển tiếp qua tần số ảo duy nhất tại Chương 3 và số liệu phân tích tương tác HOMO-LUMO giai đoạn 1, giai đoạn 2 tại Chương 4.

Câu hỏi thường gặp

1. Tiêu chí lượng tử nào trong khóa luận khẳng định cấu trúc TS(1) của phản ứng giữa propene và HBr đã được định vị chính xác?
Trạng thái chuyển tiếp $TS(1)$ được xác nhận khi tối ưu hóa thành công với từ khóa Opt=(calcfc, ts) và qua phép tính tần số dao động (Freq) xuất hiện duy nhất một tần số ảo có giá trị $\nu = -1106.71i\text{ cm}^{-1}$, đồng thời góc liên kết $\widehat{HCH}$ ($115.41^\circ$) và $\widehat{HCC}$ ($117.24^\circ$) chuyển dịch tương ứng về trạng thái lai hóa $sp^2$.

2. Cặp tương tác orbital biên (FMO) nào chi phối chiều hướng phản ứng ở giai đoạn 1 giữa propene và HBr?
Theo kết quả tính toán tại Bảng 7, độ chênh lệch hệ số giữa HOMO($C_3H_6$) và LUMO($HBr$) là $0.24335$, nhỏ hơn chênh lệch giữa HOMO($HBr$) và LUMO($C_3H_6$) là $0.27286$. Do đó, tương tác ưu đãi quyết định phản ứng là orbital HOMO của propene tấn công vào orbital LUMO của $HBr$.

3. Số liệu nhiệt động nào chứng minh giai đoạn 1 là giai đoạn quyết định tốc độ phản ứng tổng?
Giai đoạn 1 đòi hỏi vượt qua rào cản năng lượng hoạt hóa rất cao $E_{a(1)} = 80.89029873\text{ kcal/mol}$ với biến thiên entanpy mang giá trị dương ($\Delta H_1 = +48.011371\text{ kcal/mol}$, phản ứng thu nhiệt mạnh), trong khi phản ứng tổng thể lại tỏa nhiệt ($\Delta H^\circ_{298} = -23.37158727\text{ kcal/mol}$).

4. Vì sao carbocation tert-butyl có điện tích dương lớn nhất nhưng mức chênh lệch hệ số FMO ở giai đoạn 2 lại lớn hơn carbocation propyl?
Điện tích Mulliken tại carbon trung tâm của $(CH_3)_3C^+$ đạt $+0.141$ (dương nhất so với $-0.090$ của propyl). Tuy nhiên, mức chênh lệch orbital $\Delta E(HOMO - LUMO)$ của carbocation tert-butyl với $Br^-$ đạt $0.05705$, cao hơn mức $0.04486$ của carbocation propyl do sự cản trở không gian của 3 nhóm methyl xung quanh trung tâm phản ứng.

5. Vì sao nghiên cứu không khảo sát tác nhân hydro iodide (HI)?
Tác giả nêu rõ nguyên nhân là do hệ hàm cơ sở phân cực 6-31G(d) sử dụng thống nhất trong toàn bộ đề tài không hỗ trợ tối ưu hóa cho nguyên tử iod ($I$) vì iod nằm ở chu kỳ lớn, đòi hỏi phải sử dụng các hệ hàm giả thế lõi hiệu dụng như LanL2DZ.

Kết luận

Khóa luận tốt nghiệp của tác giả Võ Thị Mỹ Huyền đã ứng dụng thành công phương pháp tính toán hóa học lượng tử DFT/B3LYP/6-31G(d) trên bộ phần mềm Gaussian 09W và GaussView 6.0 để giải bài toán cơ chế phản ứng cộng electrophile của các alkene mạch ngắn. Các kết quả tính toán về cấu trúc hình học, tần số dao động ảo của trạng thái chuyển tiếp, tương tác orbital biên HOMO-LUMO và các thông số nhiệt động - động học đã kiểm chứng và giải thích thuyết phục các quy luật thực nghiệm kinh điển như quy tắc Markovnikov và tiên đề Hammond. Đây là công trình khóa luận có phương pháp luận chặt chẽ, số liệu minh bạch, mang lại giá trị học thuật và tư liệu sư phạm thiết thực cho công tác nghiên cứu, giảng dạy hóa học hữu cơ.