Giới thiệu dự án

Nghiên cứu cơ chế phản ứng hữu cơ ở cấp độ phân tử đóng vai trò then chốt trong công nghiệp tổng hợp hóa chất, dược phẩm và giảng dạy chuyên sâu. Trong hóa học thực nghiệm truyền thống, việc theo dõi các trạng thái chuyển tiếp (Transition State - TS) và tiểu phân trung gian không bền (như carbocation) có thời gian sống cực ngắn (thường $< 10^{-12}\text{ s}$) gặp nhiều rào cản kỹ thuật và tiêu tốn chi phí phân tích quang phổ phức tạp. Theo các thống kê từ ngành hóa học tính toán hiện đại, hơn $80%$ các công trình nghiên cứu cơ chế phản ứng tiên tiến trên thế giới hiện nay đều kết hợp mô phỏng hóa học lượng tử (Quantum Chemistry) để tối ưu hóa lộ trình phản ứng trước khi thử nghiệm thực địa.

Đồ án khóa luận "Ả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" tập trung giải quyết bài toán định lượng hóa ảnh hưởng cấu tạo phân tử (mạch carbon, nhóm thế, hiệu ứng không gian và điện tử) đến khả năng phản ứng cộng ái điện tử ($A_E$) của alkene với các tác nhân hydro halide ($HX$) và nước ($H_2O$).

                      Giai đoạn 1 (Chậm - Quyết định tốc độ)             Giai đoạn 2 (Nhanh)
Alkene + H-X  ───────────────────────────────────────────> [Carbocation] ────────────────────> Sản phẩm no
               TS1 (Ea = 80.89 kcal/mol, v = -1106.71i)                   + X⁻ (Nucleophile)

Mục tiêu nghiên cứu cụ thể

  1. Mô hình hóa và tối ưu hóa cấu trúc: Khảo sát hình học không gian trạng thái cơ bản của các alkene ($C_3 - C_4$: propene, but-1-ene, 2-methylpropene) và tác nhân electrophile ($HBr, H_2O$) ở mức lý thuyết phiếm hàm mật độ (DFT).
  2. Khảo sát orbital phân tử biên (FMO): Xác định chính xác các mức năng lượng $E_{\text{HOMO}}$, $E_{\text{LUMO}}$ và khe năng lượng $\Delta E$ nhằm dự đoán hướng và vị trí tấn công ái điện tử.
  3. Định vị trạng thái chuyển tiếp (TS) và sản phẩm trung gian: Tìm kiếm cấu trúc $TS(1)$ với một tần số dao động ảo duy nhất ($\nu < 0$) và xác định cấu trúc carbocation trung gian.
  4. Định lượng nhiệt động học và động học: Tính toán biến thiên Enthalpy ($\Delta H^\circ_{298}$), Entropy ($\Delta S^\circ_{298}$), năng lượng tự do Gibbs ($\Delta G^\circ_{298}$) và năng lượng hoạt hóa ($E_a$).
  5. Chứng minh quy tắc Markovnikov: Giải thích định lượng độ bền của carbocation bậc I, II, III thông qua sự giải tỏa điện tích Mulliken và hiệu ứng cảm ứng $+I$.

Phương pháp tiếp cận và kết quả kỳ vọng

  • Phương pháp tiếp cận: Ứng dụng gói phần mềm hóa học lượng tử Gaussian 09W kết hợp giao diện đồ họa GaussView 6.0, sử dụng phương pháp phiếm hàm mật độ lai hóa DFT/B3LYP cùng hệ hàm cơ sở phân cực $6-31\text{G}(d)$.
  • Kết quả đo lường được: Xác định năng lượng hoạt hóa giai đoạn quyết định tốc độ $E_a(1) \approx 80.89\text{ kcal/mol}$; định vị $TS(1)$ có tần số ảo $\nu = -1106.71i\text{ cm}^{-1}$; chứng minh phản ứng giữa propene và $HBr$ tỏa nhiệt với $\Delta H^\circ_{298} = -23.37\text{ kcal/mol}$ và tự diễn ra với $\Delta G^\circ_{298} = -12.79\text{ kcal/mol}$.
  • Phạm vi và giới hạn: Nghiên cứu khảo sát các alkene mạch ngắn ($C_3 - C_4$) trong pha khí chuẩn ($298.15\text{ K}$, $1\text{ atm}$); không mở rộng sang các hệ dị vòng hoặc tác nhân $HI$ do giới hạn hệ hàm không có thế lõi hiệu dụng (ECP).

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Phương pháp nghiên cứu Ưu điểm Nhược điểm Chi phí & Rủi ro
Thực nghiệm hóa học truyền thống (Wet Lab) Cung cấp sản phẩm thực tế, đo đạc trực tiếp hiệu suất tổng hợp. Không cô lập được trạng thái chuyển tiếp ($TS$); khó đo đạc mật độ điện tích nguyên tử tức thời. Chi phí hóa chất cao; rủi ro độc hại môi trường phòng thí nghiệm.
Bán thực nghiệm (AM1, PM3, MNDO) Tốc độ tính toán nhanh; tốn ít tài nguyên phần cứng. Độ chính xác cấu trúc và năng lượng thấp; phụ thuộc nhiều vào bộ tham số thực nghiệm hóa. Chi phí thấp nhưng sai số năng lượng hoạt hóa lớn ($> 15\text{ kcal/mol}$).
DFT-B3LYP/6-31G(d) (Giải pháp đồ án) Độ chính xác cao; mô tả tốt tương tác orbital và phân bố điện tích; tối ưu hóa $TS$ chuẩn xác. Thời gian tính toán lớn hơn bán thực nghiệm; yêu cầu cấu hình máy trạm chuyên dụng. Tiết kiệm $100%$ chi phí hóa chất thử nghiệm; độ tin cậy tương đương phổ thực nghiệm.

Yêu cầu hệ thống theo mô hình MoSCoW

  • Must have: Tối ưu hóa cấu trúc hình học ($Opt$); tính toán tần số dao động ($Freq$); xác định $HOMO/LUMO$; kiểm tra $TS$ với $1$ tần số ảo.
  • Should have: Tính toán các đại lượng nhiệt động học ($\Delta H, \Delta S, \Delta G$); hiệu chỉnh năng lượng dao động điểm không ($ZPE$) với hệ số quy chuẩn $0.9804$.
  • Could have: Khảo sát mô hình dung môi phân cực liên tục (PCM, IPCM, SCI-PCM) cho các phản ứng trong môi trường phân cực.
  • Won't have (lần này): Khảo sát các alkene đa vòng phức tạp hoặc phản ứng quang hóa kích thích nhiều photon.

Thiết kế hệ thống

flowchart TD
    A["Cấu trúc 3D sơ bộ (GaussView 6.0)"] --> B["Tạo file cấu hình đầu vào (*.gjf)"]
    B --> C["Tối ưu hóa hình học: #p B3LYP/6-31G(d) Opt"]
    C --> D{"Kiểm tra hội tụ (Forces & Displacements)"}
    D -- "Chưa hội tụ" --> C
    D -- "Hội tụ" --> E["Tính toán tần số: #p B3LYP/6-31G(d) Freq"]
    E --> F{"Kiểm tra số tần số ảo (NIMAG)"}
    F -- "NIMAG = 0 (Cấu trúc bền)" --> G["Trích xuất E_ZPE, HOMO, LUMO, Điện tích"]
    F -- "NIMAG = 1 (Trạng thái chuyển tiếp)" --> H["Xác nhận TS1: v = -1106.71i cm⁻¹"]
    F -- "NIMAG > 1" --> I["Hiệu chỉnh ma trận hằng số lực (Opt=CalcFC)"]
    I --> C
    G --> J["Phân tích Nhiệt động & Động học"]
    H --> J

Công nghệ và công cụ sử dụng

  • Phần mềm tính toán lượng tử: Gaussian 09W (Revision D.01).
  • Giao diện mô hình hóa: GaussView 6.0.16.
  • Mức lý thuyết (Level of Theory): Density Functional Theory (DFT) với phiếm hàm 3 thông số Becke kết hợp tương quan Lee-Yang-Parr (B3LYP).
  • Hệ hàm cơ sở (Basis Set): Pople Split-Valence Polarized $6-31\text{G}(d)$.
  • Tiêu chuẩn hội tụ trường tự hợp (SCF): $10^{-8}\text{ a.u.}$; điều kiện dừng tối ưu hình học: Maximum Force $< 0.000450\text{ a.u.}$, RMS Force $< 0.000300\text{ a.u.}$, Maximum Displacement $< 0.001800\text{ a.u.}$, RMS Displacement $< 0.001200\text{ a.u.}$

Phương pháp luận (Methodology)

Quy trình thực thi nghiên cứu được chuẩn hóa qua 4 giai đoạn logic:

  1. Giai đoạn 1 - Chuẩn bị mô hình phân tử: Vẽ khung liên kết nguyên tử trong GaussView, thiết lập các góc liên kết và độ dài ban đầu dựa trên dữ liệu hóa lập thể chuẩn.
  2. Giai đoạn 2 - Khảo sát đơn điểm và tối ưu hóa: Thực thi tính toán lượng tử để đạt cấu trúc cực tiểu năng lượng (Energy Minimum) trên mặt thế năng (PES).
  3. Giai đoạn 3 - Định vị trạng thái chuyển tiếp STQN: Sử dụng thuật toán Opt=QST2 với đầu vào gồm hai cấu trúc: chất đầu ($Reactant$) và chất trung gian ($Carbocation$) để dò tìm tọa độ $TS$.
  4. Giai đoạn 4 - Hiệu chỉnh và xử lý số liệu nhiệt động: Hiệu chỉnh $ZPE$ theo tỷ lệ $\lambda = 0.9804$, chuyển đổi đơn vị từ Hartree/particle sang $\text{kcal/mol}$ ($1\text{ Hartree} = 627.509391\text{ kcal/mol}$).

Triển khai và kết quả thực nghiệm

Quy trình tính toán lượng tử

Ví dụ file đầu vào Gaussian (.gjf) cho tối ưu hóa và tính tần số của phân tử propene ($CH_2=CH-CH_3$):

%chk=propene_opt.chk
%mem=4GB
%nprocshared=4
#p opt freq b3lyp/6-31g(d) pop=full gfinput gfprint

Propene Optimization and Frequency Calculation at B3LYP/6-31G(d)

0 1
 C                  0.00000000    0.00000000    0.00000000
 C                  1.33334000    0.00000000    0.00000000
 C                  2.15200000    1.25600000    0.00000000
 H                 -0.54000000    0.93500000    0.00000000
 H                 -0.54000000   -0.93500000    0.00000000
 H                  1.87334000   -0.93500000    0.00000000
 H                  3.21200000    1.01500000    0.00000000
 H                  1.92100000    1.86000000    0.88000000
 H                  1.92100000    1.86000000   -0.88000000

Dữ liệu output kiểm tra hội tụ và năng lượng dao động điểm không ($ZPE$) thu được từ Gaussian log:

 SCF Done:  E(RB3LYP) =  -117.90711310     A.U. after    11 cycles
 Convg  =    0.3421D-08             -V/T =  2.0048
 Zero-point correction=                           0.080081 (Hartree/Particle)
 Zero-point vibrational energy=                   50.25142 (Kcal/Mol)
 Sum of electronic and zero-point Energies=       -117.827032
 Sum of electronic and thermal Enthalpies=        -117.822455
 Sum of electronic and thermal Free Energies=     -117.854972

Kiểm định và đánh giá cấu trúc

1. Thông số hình học các alkene tối ưu hóa

Độ dài liên kết đôi $C=C$ phản ánh mật độ electron và khả năng tham gia phản ứng cộng:

  • Propene ($CH_2=CH-CH_3$): $d(C=C) = 1.33334\text{ \AA}$
  • But-1-ene ($CH_2=CH-CH_2-CH_3$): $d(C=C) = 1.33346\text{ \AA}$
  • 2-Methylpropene ($CH_2=C(CH_3)_2$): $d(C=C) = 1.33662\text{ \AA}$

Liên kết $C=C$ của 2-methylpropene giãn dài nhất ($1.33662\text{ \AA}$) do chịu tác động đẩy electron ($+I$) mạnh mẽ từ hai nhóm methyl ($-CH_3$), làm tăng mật độ electron tại vùng liên kết $\pi$ và tạo điều kiện thuận lợi cho tác nhân electrophile tấn công.

2. Xác thực trạng thái chuyển tiếp giai đoạn 1 ($TS1$)

Trong phản ứng giữa propene và $HBr$, cấu trúc $TS1$ xác định được có các thông số:

  • Độ dài liên kết $C=C$ tại $TS1$: Giãn từ $1.33334\text{ \AA}$ lên $1.41106\text{ \AA}$, minh chứng cho sự chuyển đổi từ liên kết đôi $sp^2$ sang liên kết đơn $sp^3$.
  • Góc liên kết: $\angle HCH = 115.41^\circ$, $\angle HCC = 117.24^\circ$ (tiến gần cấu trúc tam giác phẳng của carbocation).
  • Tần số dao động ảo: $\nu = -1106.71i\text{ cm}^{-1}$ (dao động kéo dãn liên kết $H-C$ và phân cắt $H-Br$).
       H                H
        \              /
         C ═══ C      +   H ─── Br
        /       \
       H         CH₃
   [Chất đầu: d(C=C)=1.333 Å]
                │
                ▼ (Tấn công chậm, Ea = 80.89 kcal/mol)
       H ... H ... Br
        \   :
         C ─-─ C ⁺ ── H        [TS1: d(C=C)=1.411 Å, v = -1106.71i cm⁻¹]
        /       \
       H         CH₃
                │
                ▼ (Hình thành tiểu phân trung gian)
       H         H
        \       / ⁺
     H ── C ─── C ── CH₃   +   Br ⁻
        /       
       H
   [Carbocation Propyl bậc II: Phẳng sp², Góc CCC=121.32°]
                │
                ▼ (Tấn công nhanh giai đoạn 2)
       H    Br   H
        \   │   /
     H ── C ─ C ─ CH₃
        /   │   \
       H    H    H
   [2-Bromopropane: Sản phẩm chính Markovnikov]

Kết quả đạt được và thông số nhiệt động

Bảng tổng hợp năng lượng orbital biên FMO (giai đoạn 1)

Hợp chất Hệ số HOMO (a.u.) Hệ số LUMO (a.u.) $\Delta E(\text{HOMO}{\text{Alkene}} - \text{LUMO}{HX})$ $\Delta E(\text{HOMO}{HX} - \text{LUMO}{\text{Alkene}})$ Hướng tương tác ưu tiên
$HBr$ $0.30115$ $0.00645$ — — —
$C_3H_6$ (Propene) $0.24980$ $0.02829$ $0.24335$ $0.27286$ $\text{HOMO}{C_3H_6} \rightarrow \text{LUMO}{HBr}$

Chênh lệch năng lượng giữa $\text{HOMO}(C_3H_6)$ và $\text{LUMO}(HBr)$ nhỏ hơn ($0.24335\text{ a.u.}$ so với $0.27286\text{ a.u.}$), khẳng định electron $\pi$ từ orbital $HOMO$ của propene sẽ chủ động xen phủ và chuyển dịch sang orbital phản liên kết $LUMO$ của $HBr$.

Bảng thông số nhiệt động phản ứng tổng: $C_3H_6 + HBr \rightarrow C_3H_7Br$

Đại lượng Giá trị tính toán DFT Đơn vị Ý nghĩa hóa học
Biến thiên Enthalpy ($\Delta H^\circ_{298}$) $-23.3716$ $\text{kcal/mol}$ Phản ứng tỏa nhiệt mạnh ($\Delta H < 0$).
Biến thiên Entropy ($\Delta S^\circ_{298}$) $-35.4790$ $\text{cal/mol}\cdot\text{K}$ Giảm độ hỗn loạn do $2$ phân tử kết hợp thành $1$ phân tử no.
Năng lượng tự do Gibbs ($\Delta G^\circ_{298}$) $-12.7935$ $\text{kcal/mol}$ Phản ứng tự phát về mặt nhiệt động học ở $298\text{ K}$.
Năng lượng hoạt hóa ($E_a(1)$) $80.8903$ $\text{kcal/mol}$ Rào cản thế năng giai đoạn tạo carbocation (quyết định tốc độ).

Mật độ điện tích Mulliken trên C trung tâm của carbocation

Cấu trúc Carbocation Bậc Carbocation Mật độ điện tích C trung tâm Độ bền tương đối
Propyl ($CH_3-\overset{+}{C}H-CH_3$) Bậc II $-0.090$ Khá bền
Sec-butyl ($CH_3-\overset{+}{C}H-CH_2CH_3$) Bậc II $-0.089$ Bền
Tert-butyl ($(CH_3)_3\overset{+}{C}$) Bậc III $+0.141$ Cực kỳ bền (giải tỏa điện tích cao)

Theo Tiên đề Hammond (Hammond's Postulate), do giai đoạn 1 thu nhiệt mạnh ($\Delta H^\ddagger = 48.01\text{ kcal/mol}$), cấu trúc trạng thái chuyển tiếp $TS1$ mang nhiều đặc tính hình học tương đồng với carbocation trung gian. Carbocation bậc III ($(CH_3)_3C^+$) có độ giải tỏa điện tích tốt nhất nhờ $3$ nhóm methyl đẩy electron, giải thích nguyên nhân tốc độ phản ứng cộng của 2-methylpropene vượt trội so với propene và but-1-ene.

Đổi mới và đóng góp khoa học

Cải tiến kỹ thuật nổi bật

  • Định lượng hóa thuyết cấu trúc: Chuyển đổi các giải thích định tính kinh điển (hiệu ứng cảm ứng $+I$, siêu liên hợp) thành các số liệu định lượng cụ thể: mật độ điện tích Mulliken, độ dài liên kết chính xác đến $10^{-5}\text{ \AA}$, và khe năng lượng orbital $HOMO-LUMO$.
  • Mô phỏng tường minh trạng thái chuyển tiếp: Dựng thành công tọa độ phản ứng không dừng $TS1$ với một tần số dao động ảo duy nhất ($\nu = -1106.71i\text{ cm}^{-1}$), điều mà thực nghiệm quang phổ thông thường không thể ghi nhận trực tiếp.
  • Tiết kiệm chi phí nghiên cứu: Giảm thiểu $100%$ lượng dung môi độc hại ($CH_2Cl_2, CCl_4$) và hóa chất halogenua tiêu tốn trong giai đoạn thử nghiệm sơ bộ.

So sánh với các giải pháp hiện hành

Tiêu chí Khảo sát thực nghiệm phòng Lab Bán thực nghiệm PM3 Mô hình DFT-B3LYP/6-31G(d) (Đồ án)
Chi phí triển khai Rất cao ($10 - 50$ triệu VNĐ/loạt phản ứng) Rất thấp (chạy trên PC thường) Thấp (đầu tư máy trạm tính toán một lần)
Thời gian khảo sát $2 - 4$ tuần tổng hợp và đo phổ $1 - 2$ phút $2 - 6$ giờ/cấu trúc
Độ chính xác $\Delta H, \Delta G$ Phụ thuộc độ tinh khiết mẫu Sai số $> 20%$ Sai số $< 5%$ so với thực nghiệm NIST
Quan sát trạng thái chuyển tiếp Không thể thực hiện trực tiếp Độ tin cậy hình học thấp Rõ nét hình học và mode dao động ảo

Ứng dụng thực tế và triển khai

Kịch bản ứng dụng trong thực tế

  1. Giảng dạy và bồi dưỡng chuyên sâu: Ứng dụng mô hình hóa trực quan $3D$ để giảng dạy bản chất phản ứng cộng $A_E$, cơ chế định hướng Markovnikov cho học sinh chuyên hóa và sinh viên ngành Hóa học sư phạm/hóa dược.
  2. Thiết kế quy trình tổng hợp công nghiệp: Dự đoán hiệu suất nhiệt động ($\Delta G$) và năng lượng rào cản ($E_a$) của các dẫn xuất alkyl halide trước khi mở rộng quy mô phản ứng (scale-up) trong công nghiệp hóa dầu.
  3. Sàng lọc xúc tác: Định hướng bổ sung xúc tác acid Lewis để hạ thấp năng lượng hoạt hóa $E_a(1) = 80.89\text{ kcal/mol}$ xuống ngưỡng phản ứng tối ưu ở nhiệt độ phòng.

Cấu hình triển khai hệ thống tính toán

  • Phần cứng tối thiểu: CPU 4 Cores (Intel Core i5/AMD Ryzen 5 trở lên), RAM 8GB DDR4, Ổ cứng SSD 100GB trống.
  • Phần cứng khuyến nghị cho mở rộng: Workstation dual Intel Xeon / AMD EPYC (16-32 Cores), RAM 64GB ECC, lưu trữ NVMe RAID 0.
  • Môi trường phần mềm: HĐH Linux Ubuntu 20.04 LTS / Windows 10 Pro 64-bit; Gaussian 09W/16; GaussView 6.0.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Hệ hàm cơ sở: Hệ hàm $6-31\text{G}(d)$ chưa bao gồm hàm khuếch tán ($+$) cho anion halogenua lớn như $I^-$ và chưa tích hợp thế lõi hiệu dụng (ECP) cho các nguyên tố chu kỳ cao.
  • Môi trường phản ứng: Mô hình khảo sát chính đặt trong pha khí, chưa đánh giá đầy đủ năng lượng solvat hóa phức tạp của các dung môi phân cực mạnh qua mô hình SMD/PCM.

Hướng nghiên cứu mở rộng

  • Nâng cấp mức lý thuyết lên hệ hàm cơ sở mở rộng $\text{B3LYP}/6-311++\text{G}(d,p)$ hoặc phiếm hàm tầm xa $\omega\text{B97X-D}$ để mô tả chính xác tương tác phân tán Van der Waals.
  • Mở rộng khảo sát phản ứng cộng electrophile với các hệ liên hợp phức tạp (dien, alkyne) và xúc tác acid Bronsted trong dung môi phân cực.

Đối tượng hưởng lợi

  • Sinh viên & Học viên cao học: Tiếp cận tài liệu chuẩn hóa về quy trình tính toán lượng tử (Gaussian input/output, xử lý $ZPE$, phân tích FMO).
  • Giáo viên & Giảng viên Hóa học: Sở hữu bộ tư liệu hình ảnh và dữ liệu năng lượng trực quan phục vụ bài giảng cơ chế phản ứng hữu cơ chuyên sâu.
  • Kỹ sư R&D Hóa dược: Ứng dụng quy trình mô phỏng rút ngắn thời gian thử nghiệm phản ứng halogen hóa và hydrat hóa trong phòng Lab.
  • Nhà nghiên cứu hóa lý: Kế thừa bộ dữ liệu nhiệt động ($\Delta H, \Delta S, \Delta G, E_a$) chính xác làm mốc tham chiếu cho các nghiên cứu alkene mở rộng.

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

1. Yêu cầu phần cứng tối thiểu để chạy mô phỏng Gaussian cho các alkene này là gì?

Hệ thống cần tối thiểu CPU 4 nhân thực, 8GB RAM và khoảng 50GB không gian đĩa trống cho các file checkpoint (.chk) và scratch files (.rwf). Khuyến nghị sử dụng hệ thống Linux để tối ưu hóa hiệu năng tính toán đa luồng OpenMP.

2. Làm thế nào để xác nhận chắc chắn một cấu trúc tìm được là trạng thái chuyển tiếp (TS)?

Cấu trúc trạng thái chuyển tiếp hợp lệ bắt buộc phải thỏa mãn $2$ tiêu chí: (1) Kiểm tra tần số dao động ($Freq$) chỉ xuất hiện duy nhất $1$ giá trị tần số ảo ($\text{NIMAG} = 1$, ví dụ $\nu = -1106.71i\text{ cm}^{-1}$); (2) Vector dao động của tần số ảo đó phải mô tả đúng chuyển động phá vỡ liên kết cũ và hình thành liên kết mới dọc theo tọa độ phản ứng (Intrinsic Reaction Coordinate - IRC).

3. Tại sao cần nhân hệ số hiệu chỉnh 0.9804 cho năng lượng dao động điểm không (ZPE)?

Do các phương pháp lượng tử (bao gồm B3LYP) sử dụng phép gần đúng dao động điều hòa (Harmonic Approximation) nên thường đánh giá cao hơn tần số dao động thực tế. Hệ số kinh nghiệm $\lambda = 0.9804$ giúp bù trừ sai số phi điều hòa và tương quan electron ở mức DFT/6-31G(d).

4. Tại sao đồ án không thực hiện tính toán với tác nhân hydro iodide (HI)?

Nguyên tử Iodine ($I$) có số hiệu nguyên tử lớn ($Z = 53$) với nhiều electron lõi. Hệ hàm Pople chuẩn $6-31\text{G}(d)$ không hỗ trợ đầy đủ cho Iodine. Để tính toán $HI$, cần sử dụng hệ hàm cơ sở thế lõi hiệu dụng như LanL2DZ hoặc Def2-TZVP.

5. Chi phí đầu tư phần mềm và thời gian hoàn vốn (ROI) khi ứng dụng hóa học lượng tử?

So với chi phí vận hành phòng lab hữu cơ tiêu chuẩn (từ $200 - 500$ triệu VNĐ/năm cho hóa chất tinh khiết và bảo dưỡng thiết bị đo phổ), việc đầu tư máy trạm tính toán lượng tử (khoảng $30 - 50$ triệu VNĐ) giúp hoàn vốn sau dưới $6$ tháng thông qua việc loại bỏ các thử nghiệm sai hỏng.

Kết luận

Khóa luận tốt nghiệp "Ả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" đã ứng dụng thành công phương pháp phiếm hàm mật độ DFT-B3LYP/6-31G(d) trên bộ phần mềm Gaussian 09W và GaussView 6.0. Đề tài đã giải quyết trọn vẹn bài toán mô hình hóa hình học, giải thích định lượng cơ chế phản ứng hai giai đoạn, định vị chính xác trạng thái chuyển tiếp $TS1$ ($\nu = -1106.71i\text{ cm}^{-1}$, $E_a = 80.89\text{ kcal/mol}$), và chứng minh tính ưu việt của quy tắc Markovnikov thông qua mật độ điện tích và thuyết orbital phân tử biên FMO. Kết quả nghiên cứu là nền tảng vững chắc kết nối lý thuyết lượng tử với thực hành sư phạm và nghiên cứu hóa học ứng dụng hiện đại.