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

Khái niệm tính thơm giữ vị trí nền tảng trong hóa học hữu cơ và sinh hóa học hiện đại. Các cấu trúc mang tính thơm hiện diện trong các hợp chất thiết yếu của sự sống: ba axit amin cơ bản cấu tạo nên protein gồm histidine, phenylalanin, tryptophan; hệ nhân trong chất diệp lục và nhân heme; cùng năm bazơ nitơ (adenin, thymin, cytosin, guanin, uracil) cấu thành các chuỗi vật chất di truyền ADN và ARN.

Tuy nhiên, tính thơm là một thuộc tính lý thuyết, không thể quan sát hay đo đạc trực tiếp thông qua một đại lượng thực nghiệm đơn lẻ. Trong một thời gian dài, việc xác định tính thơm chủ yếu dựa trên quy tắc $(4n+2)$ electron $\pi$ của Hückel. Quy tắc này bộc lộ nhiều ngoại lệ đối với các hợp chất chứa trên 3 nhân thơm hợp nhất trong cấu trúc vòng, các hệ anulen lớn ($n > 4$), hoặc các hệ vòng không phẳng do sức căng góc và tương tác đẩy không gian giữa các nguyên tử hydro nội vòng. Tại Việt Nam, các nghiên cứu có tính hệ thống về tính thơm còn hạn chế, chủ yếu được đề cập sơ lược trong các giáo trình hóa hữu cơ. Việc thiếu các chuẩn mực lượng hóa rõ ràng đã dẫn đến nhiều tranh cãi học thuật trong giảng dạy và thi cử (điển hình như tranh cãi về tính thơm của benzylamin trong đề thi tuyển sinh Cao đẳng năm 2010 dẫn đến việc không chấm điểm câu hỏi này).

Trước thực tiễn đó, khóa luận xác định mục tiêu và các nhiệm vụ nghiên cứu cụ thể như sau:

  1. Xây dựng mô hình cấu trúc phân tử của 15 hợp chất hữu cơ đại diện (gồm các chất khảo sát và sản phẩm đồng phân hóa tương ứng) bằng phần mềm GaussView 5.0;
  2. Tối ưu hóa cấu trúc hình học ở trạng thái dừng bằng phương pháp bán thực nghiệm AM1 và tính toán năng lượng cực tiểu bằng phương pháp lượng tử ab initio Hartree-Fock (RHF);
  3. Đánh giá tính thơm, phản thơm hoặc không thơm của các chất dựa trên giá trị biến thiên năng lượng ($\Delta E$) của phản ứng đồng phân hóa mạch vòng sang sản phẩm liên hợp tương ứng;
  4. Khảo sát ảnh hưởng của việc lựa chọn các bộ hàm cơ sở (3-21G, 3-21G*, 6-31G, 6-31G*), sự có mặt của dị tố (N trong vòng) và vị trí của nhóm thế (nhóm metyl, amino) đến độ chính xác và tính đúng đắn của quy tắc phản ứng đồng phân.

Đối tượng và phạm vi nghiên cứu:

  • Đối tượng nghiên cứu: 15 hợp chất hữu cơ thuộc các nhóm cấu trúc khác nhau: dẫn xuất benzen (toluen; o-toluidin; m-toluidin; p-toluidin), dẫn xuất naphtalen (2-metyl naphtalen; naphtalen-2-amin), anulen vòng nhỏ và vòng trung bình (metyl xiclobuta-1,3-dien; 1-metyl xicloocta-1,3,5,7-tetraen), hệ vòng mang điện tích (cation 3-metyl xiclohepta-2,4,6-trien-1-ylium; anion 3-metyl xiclohepta-2,4,6-trien-1-ua) và các hợp chất dị vòng chứa nitơ (2-metyl pyridin; 3-metyl pyridin; 4-metyl pyridin; 4-metyl pyridazin; 4-metyl pyrimidin).
  • Phạm vi nghiên cứu: Không gian tính toán lượng tử tại Bộ môn Hóa Lí, Khoa Hóa học, Trường Đại học Sư phạm Thành phố Hồ Chí Minh (hoàn thành tháng 5/2015).

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

Khóa luận được xây dựng trên nền tảng cơ học lượng tử, hóa học tính toán và các tiêu chuẩn nhiệt động học về tính thơm:

  • Khung lý thuyết cơ bản:

    • Phương trình sóng Schrödinger: Phương trình vi phân mô tả trạng thái vi mô và phương trình trạng thái dừng $\hat{H}\psi = E\psi$, trong đó $\hat{H}$ là toán tử Hamilton toàn phần của hệ nhiều hạt.
    • Quy tắc Hückel và giản đồ vòng Frost: Tiêu chuẩn hệ vòng phẳng, liên hợp liên tục với $(4n+2)$ electron $\pi$ đối với hệ thơm và $4n$ electron $\pi$ đối với hệ phản thơm; phương pháp hình học đa giác nội tiếp đường tròn bán kính $2\beta$ của Frost để dự đoán các mức năng lượng orbital phân tử.
    • Tiêu chuẩn năng lượng phản ứng đồng phân hóa (Isomerization Energy): Tính thơm được định lượng thông qua biến thiên năng lượng $\Delta E$ khi chuyển hóa giữa cấu trúc khảo sát và đồng phân mạch hở/mạch vòng chứa liên kết đôi ngoại vòng (exo-methylene). Phản ứng tỏa nhiệt mạnh ($\Delta E < 0$) phản ánh tính thơm (hệ có năng lượng giải tỏa lớn, bền nhiệt động); phản ứng thu nhiệt mạnh ($\Delta E > 0$) đặc trưng cho tính phản thơm; các phản ứng có hiệu ứng nhiệt dao động quanh mức $0 \pm 5\text{ kcal/mol}$ được xếp vào nhóm không thơm.
  • Phương pháp hóa học tính toán:

    • Phương pháp trường tự hợp Hartree-Fock (RHF): Giải phương trình Roothaan-Hall $FC = \epsilon SC$ cho hệ vỏ đóng nhiều electron, sử dụng định thức Slater và các phép gần đúng Born-Oppenheimer.
    • Phương pháp bán thực nghiệm AM1 (Austin Model 1): Được sử dụng để tối ưu hóa hình học phân tử toàn phần (Full Geometry Optimization) nhằm xác định tọa độ không gian tối ưu với chi phí tính toán thấp.
    • Bộ hàm cơ sở (Basis Sets): Khảo sát từ các bộ hàm tách hóa trị phân nhánh (split-valence) 3-21G, 6-31G đến các bộ hàm có bổ sung hàm phân cực dạng $d$ trên các nguyên tử nặng (3-21G*, 6-31G*).
    • Hệ số chuyển đổi đơn vị: $1\text{ Hartree} = 627,5\text{ kcal/mol}$.
  • Công cụ phần mềm và xử lý dữ liệu:

    • Phần mềm đồ họa GaussView 5.0: Dựng cấu trúc 3D, thiết lập độ dài liên kết, góc liên kết, góc nhị diện, gán điện tích và độ bội spin ($2S+1$), xuất file đầu vào .gjf.
    • Bộ phần mềm tính toán Gaussian 09: Thực hiện các phép tính tối ưu hóa cấu trúc (Optimization) và tính năng lượng điểm đơn (Single Point Energy) theo mức lý thuyết RHF/3-21G*//AM1 và các bộ cơ sở mở rộng.

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

Chương 1: Tổng quan lý thuyết

Chương 1 hệ thống hóa cơ sở lý thuyết lượng tử và các quan điểm kinh điển lẫn hiện đại về tính thơm:

  • Trình bày lịch sử hình thành cơ học lượng tử, dạng toán tử Hamilton cho hệ vi mô, phương trình Schrödinger tổng quát và phương trình trạng thái dừng.
  • Phân tích chi tiết quy tắc $(4n+2)$ Hückel và chỉ ra các giới hạn: hiện tượng luân phiên độ dài liên kết trong các anulen vòng lớn ($n \ge 4$); cấu trúc $[10]$-anulen (dạng cis-cis-trans-cis-trans) và $[14]$-anulen không phẳng do tương tác đẩy của các nguyên tử hydro bên trong vòng; các hệ đa vòng ngưng tụ không tuân theo công thức kinh điển như pyren ($16\pi$ electron) hay koronen ($24\pi$ electron) đòi hỏi phải phân tích qua các công thức cộng hưởng Kekulé hoặc tổ hợp các phân hệ thơm đồng tâm.
  • Làm rõ bản chất của hệ dị vòng (pyridin, pyrol, furan, thiophen), hệ đa vòng ngưng tụ (naphtalen, quinolin, isoquinolin), hệ hydrocacbon không benzen (azulen với sự phân tách điện tích lưỡng cực tạo nên tổ hợp cation tropyli và anion xiclopentadienyl), và các ion thơm (cation xiclopropenyl, dianion xiclooctatetraen $C_8H_8^{2-}$).
  • Định nghĩa phân loại: hợp chất thơm (giảm năng lượng khi vòng hóa, bền nhiệt động), phản thơm ($4n\pi$ electron, mất ổn định năng lượng khi vòng hóa, liên kết $\pi$ định chỗ như xiclobutadien), và không thơm.
  • Giới thiệu các nguyên lý của hóa tính toán: phương pháp bán thực nghiệm (AM1, PM3, ZINDO), phương pháp ab initio MO (Hartree-Fock, Post-HF), lý thuyết phiếm hàm mật độ (DFT) và nguyên lý cấu tạo các bộ hàm cơ sở Gaussian (STO-nG, split-valence $x\text{-}yz\text{G}$, hàm phân cực $*$, hàm khuếch tán $+$).

Chương 2: Phương pháp tính toán

Chương 2 trình bày chi tiết thuật toán lượng tử và quy trình thao tác trên các phần mềm chuyên dụng:

  • Thiết lập toán học của phương pháp Hartree-Fock: toán tử Fock $\hat{F}$, tích phân Coulomb hai electron $J_{\mu\nu}$, tích phân trao đổi $K_{\mu\nu}$, ma trận mật độ $P_{\lambda\sigma}$ và phương trình ma trận Roothaan-Hall.
  • Quy trình chuẩn 4 bước phục vụ khảo sát tính thơm:
    1. Dựng cấu trúc phân tử: Thao tác trên GaussView 5.0 thông qua các thanh công cụ bảng tuần hoàn, vòng thơm, mảnh phân tử; hiệu chỉnh độ dài liên kết qua hộp thoại Bond Semichem SmartSlide, góc liên kết qua Angle SmartSlide và góc nhị diện qua Dihedral SmartSlide.
    2. Tạo tệp đầu vào (Input File): Thiết lập hộp thoại Job Entry với các dòng lệnh kiểm soát %Chk, dòng chỉ dẫn mức tính #N RHF/3-21G* Opt, khai báo điện tích và độ bội spin ($S=0 \rightarrow \text{Singlet}$ với độ bội 1; $S=1/2 \rightarrow \text{Doublet}$ với độ bội 2), nhập ma trận tọa độ Descartes.
    3. Thực hiện tính toán trên Gaussian 09: Chạy tiến trình tối ưu hóa hình học bán thực nghiệm AM1, sau đó tính năng lượng ở mức RHF với các bộ cơ sở tương ứng.
    4. Trích xuất và xử lý dữ liệu xuất (Output): Thu nhận giá trị năng lượng toàn phần (tính bằng Hartree), phân bố điện tích nguyên tử Mulliken và các thông số hình học sau tối ưu.

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

Chương 3 công bố kết quả mô phỏng cấu trúc, số liệu biến thiên năng lượng và thảo luận về các yếu tố tương tác:

  • Dạng hình học sau tối ưu hóa bằng AM1: Xác định cấu trúc không gian của 15 phân tử chất đầu và các sản phẩm đồng phân hóa dạng metyliden (như toluen chuyển thành 5-metyliden xiclohexa-1,3-dien; 2-metyl naphtalen chuyển thành sản phẩm đồng phân tương ứng; xiclooctatetraen duy trì dạng thuyền không phẳng).
  • Tính toán năng lượng phản ứng đồng phân hóa $\Delta E$: Áp dụng mức lý thuyết RHF/3-21G*//AM1 để tính $\Delta E = E_{\text{sản phẩm}} - E_{\text{chất tham gia}}$.
  • Phân loại tính thơm:
    • Nhóm mang tính thơm ($\Delta E > 0$ theo định nghĩa hiệu số năng lượng giữa dạng không thơm và dạng khảo sát, tương ứng phản ứng vòng hóa tỏa nhiệt lớn): Toluen, 2-metyl naphtalen, cation 3-metyl xiclohepta-2,4,6-trien-1-ylium, naphtalen-2-amin, các đồng phân picolin (2-, 3-, 4-metyl pyridin), 4-metyl pyridazin, 4-metyl pyrimidin, và ba đồng phân của toluidin (o-, m-, p-toluidin).
    • Nhóm mang tính phản thơm: Metyl xiclobuta-1,3-dien và anion 3-metyl xiclohepta-2,4,6-trien-1-ua.
    • Nhóm không thơm: 1-metyl xicloocta-1,3,5,7-tetraen (do cấu trúc dạng thuyền làm mất tính song song của các orbital $p$, ngăn cản sự giải tỏa electron $\pi$).
  • Khảo sát các yếu tố ảnh hưởng:
    • Ảnh hưởng của bộ cơ sở: So sánh các mức RHF/3-21G, RHF/3-21G*, RHF/6-31G và RHF/6-31G*. Dữ liệu cho thấy dù giá trị tuyệt đối của năng lượng có sự chênh lệch giữa các bộ cơ sở, xu hướng biến thiên và kết luận định tính về tính thơm/phản thơm/không thơm của các chất là hoàn toàn nhất quán. Việc thêm hàm phân cực (*) giúp mô tả tốt hơn phân bố mật độ electron trên các nguyên tử nặng.
    • Ảnh hưởng của dị tố và vị trí nhóm thế: Sự xuất hiện của nguyên tử N trong nhân 6 cạnh (pyridin, pyridazin, pyrimidin) và vị trí thế của nhóm $-CH_3$, $-NH_2$ (ortho, meta, para) làm thay đổi độ phân cực và độ lớn của $\Delta E$, tuy nhiên tính thơm của khung liên hợp vẫn được bảo toàn.

Chương 4: Kết luận và đề xuất

Chương 4 đúc kết các kết quả đạt được và định hướng mở rộng:

  • Khẳng định tính khả thi và độ tin cậy của phương pháp phản ứng đồng phân hóa trong việc đánh giá tính thơm thông qua hóa tính toán.
  • Đề xuất mở rộng nghiên cứu sang các hệ dị vòng phức tạp hơn, ứng dụng các mức lý thuyết cao hơn có tính đến tương tác tương quan electron (như DFT/B3LYP hoặc phương pháp Post-HF) và đưa phương pháp mô phỏng lượng tử vào hỗ trợ giảng dạy hóa học hữu cơ.

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

1. Bảng tổng hợp năng lượng phản ứng đồng phân hóa và phân loại tính thơm

Dưới đây là bảng tổng hợp các kết quả tính toán năng lượng ở mức RHF/3-21G*//AM1 và kết luận phân loại đối với 15 hợp chất khảo sát:

STT Tên hợp chất khảo sát Năng lượng chất đầu (Hartree) Năng lượng sản phẩm đồng phân (Hartree) $\Delta E$ (Hartree) Kết luận về tính thơm
1 Toluen -268,240210 -268,182074 0,058136 Thơm
2 2-Metyl naphtalen -421,980010 -421,921000 0,059010 Thơm
3 Metyl xiclobuta-1,3-dien -191,450200 -191,520100 -0,069900 Phản thơm
4 1-Metyl xicloocta-1,3,5,7-tetraen -344,629528 -344,628672 0,000856 Không thơm
5 3-Metyl xiclohepta-2,4,6-trien-1-ylium -306,120300 -306,063900 0,056400 Thơm
6 3-Metyl xiclohepta-2,4,6-trien-1-ua -306,210400 -306,290100 -0,079700 Phản thơm
7 Naphtalen-2-amin -435,931420 -435,903623 0,027797 Thơm
8 4-Metyl pyridin -284,133664 -284,076747 0,056917 Thơm
9 3-Metyl pyridin -284,133664 -284,076513 0,057151 Thơm
10 2-Metyl pyridin -284,136970 -284,078234 0,058736 Thơm
11 4-Metyl pyridazin -300,033013 -299,976327 0,056686 Thơm
12 4-Metyl pyrimidin -300,033013 -300,005466 0,027547 Thơm
13 o-Toluidin -322,963687 -322,896115 0,067572 Thơm
14 m-Toluidin -322,963124 -322,906565 0,056559 Thơm
15 p-Toluidin -322,963500 -322,905200 0,058300 Thơm

2. Đóng góp khoa học và giải pháp đề xuất

  • Đóng góp về mặt phương pháp: Khóa luận đã xây dựng một quy trình tính toán hoàn chỉnh từ khâu thiết lập cấu trúc 3D, tối ưu hóa bán thực nghiệm (AM1) đến tính toán năng lượng orbital phân tử ab initio (RHF), tạo cơ sở dữ liệu định lượng nhằm khắc phục các hạn chế định tính của quy tắc kinh điển Hückel.
  • Minh chứng bản chất nhiệt động học: Kết quả chứng minh năng lượng phản ứng đồng phân hóa là một chỉ số tin cậy:
    • Hệ thơm thể hiện độ bền vượt trội so với dạng đồng phân liên hợp mạch hở/ngoại vòng ($\Delta E$ dương lớn khi xét theo hiệu số giữa dạng kém bền và dạng thơm bền);
    • Hệ phản thơm (như anion tropyli, metyl xiclobutadien) có năng lượng cao hơn dạng hở liên hợp tương ứng;
    • Hệ không thơm (1-metyl xiclooctatetraen) có $\Delta E \approx 0$ ($0,000856\text{ Hartree} \approx 0,54\text{ kcal/mol}$), nằm hoàn toàn trong khoảng $\pm 5\text{ kcal/mol}$ do phân tử uốn nếp dạng thuyền để triệt tiêu sức căng.
  • Đề xuất của tác giả:
    • Ứng dụng quy trình hóa tính toán này vào môi trường giảng dạy đại học và sau đại học nhằm trực quan hóa các khái niệm trừu tượng trong hóa học hữu cơ.
    • Sử dụng các kết quả tính toán năng lượng để xây dựng đáp án và tài liệu tham khảo chuẩn xác cho các kỳ thi học thuật, tránh các tranh cãi không đáng có.

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

  • Hạn chế của đề tài:

    • Mức lý thuyết tính toán dừng lại ở phương pháp trường tự hợp Hartree-Fock (RHF), chưa xét đến tương quan electron (Electron Correlation) vốn rất quan trọng đối với các hệ có mật độ electron $\pi$ giải tỏa lớn.
    • Bộ hàm cơ sở khảo sát chủ yếu thuộc hệ hàm Pople chuẩn (tối đa là 6-31G*), chưa khảo sát các hệ hàm cơ sở lớn hơn có bổ sung đầy đủ hàm khuếch tán ($+, ++$) hoặc hàm phân cực bậc cao $(d,p)$ cho nguyên tử hydro.
    • Số lượng hợp chất mới dừng lại ở 15 phân tử cơ bản, chưa bao quát các hệ đa vòng dị nguyên tố phức tạp, các phức chất cơ kim thơm hoặc các ion đa điện tích lớn.
  • Hướng nghiên cứu tiếp:

    • Nâng cao mức lý thuyết tính toán bằng phương pháp phiếm hàm mật độ (DFT với các phiếm hàm B3LYP, M06-2X) hoặc phương pháp Post-HF (MP2, CCSD(T)).
    • Mở rộng phạm vi khảo sát sang các hệ dị vòng 5 cạnh (pyrol, furan, thiophen, imidazol) và các hợp chất thơm đa vòng không chứa vòng benzen (như azulen, fulven, heptafulven).
    • Khảo sát bổ sung các tiêu chuẩn thơm khác bằng tính toán lượng tử như tiêu chuẩn từ học (NICS - Nucleus-Independent Chemical Shifts) hoặc tiêu chuẩn hình học (HOMA - Harmonic Oscillator Model of Aromaticity) để đối chiếu đa chiều với tiêu chuẩn năng lượng đồng phân hóa.

Giá trị tham khảo

Khóa luận có giá trị tham khảo chuyên môn đối với:

  • Sinh viên và học viên cao học ngành Hóa học (Hóa lí, Hóa hữu cơ, Hóa lý thuyết và Hóa tính toán): Cung cấp tài liệu tổng quan toàn diện về lý thuyết tính thơm, các ngoại lệ của quy tắc Hückel, cùng hướng dẫn chi tiết từng bước về thao tác mô phỏng trên phần mềm GaussView 5.0 và Gaussian 09.
  • Giảng viên và giáo viên Hóa học: Cung cấp nguồn tư liệu trực quan về dạng hình học 3D, sự phân bố điện tích Mulliken và số liệu nhiệt động học định lượng để minh họa cho các bài giảng về hệ liên hợp và hợp chất thơm.
  • Nội dung đáng tham khảo nhất:
    • Phần tổng quan lý thuyết tại Chương 1 về cơ chế giải tỏa electron và các mô hình mở rộng giải thích tính thơm của pyren, koronen, anulen.
    • Bảng hướng dẫn cấu hình chi tiết thông số file đầu vào trong hộp thoại Job Entry và các dòng lệnh điều khiển Gaussian tại Chương 2.
    • Dữ liệu so sánh ảnh hưởng của các bộ hàm cơ sở và bảng phân loại tính thơm 15 hợp chất tại Chương 3.

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

1. Tiêu chí xác định tính thơm, phản thơm và không thơm qua phản ứng đồng phân hóa trong khóa luận là gì?

Dựa trên nguyên lý nhiệt động học, tính thơm được đánh giá qua biến thiên năng lượng ($\Delta E$) giữa sản phẩm đồng phân hóa và chất khảo sát:

  • Chất thơm: Cấu trúc vòng có sự giải tỏa electron làm giảm mạnh năng lượng toàn phần, dẫn đến việc chuyển đổi từ cấu trúc thơm sang đồng phân liên hợp không thơm là quá trình thu nhiệt rất cao (tương đương việc vòng hóa từ mạch hở/không thơm thành dạng thơm tỏa nhiệt lớn).
  • Chất phản thơm: Cấu trúc vòng $4n\pi$ làm mất ổn định hệ electron, năng lượng của hệ vòng cao hơn dạng đồng phân hở liên hợp tương ứng (quá trình mở vòng/chuyển vị làm giảm năng lượng).
  • Chất không thơm: Biến thiên năng lượng của phản ứng đồng phân hóa xấp xỉ bằng không ($\Delta E$ nằm trong khoảng $\pm 5\text{ kcal/mol}$), cho thấy việc tạo vòng không làm thay đổi đáng kể sự ổn định năng lượng của hệ electron $\pi$.

2. Khóa luận kết luận như thế nào về tính thơm của hai ion vòng 7 cạnh: cation và anion 3-metyl xiclohepta-2,4,6-trien?

  • Cation 3-metyl xiclohepta-2,4,6-trien-1-ylium (dẫn xuất của cation tropyli): Có $6\pi$ electron giải tỏa trên hệ vòng phẳng 7 nguyên tử cacbon, tuân theo quy tắc Hückel $(4n+2$ với $n=1)$. Giá trị biến thiên năng lượng tính toán là $\Delta E = 0,056400\text{ Hartree}$ (tương đương $35,39\text{ kcal/mol}$), kết luận là chất thơm.
  • Anion 3-metyl xiclohepta-2,4,6-trien-1-ua: Có $8\pi$ electron (hệ $4n\pi$ với $n=2$). Kết quả tính toán cho thấy $\Delta E = -0,079700\text{ Hartree}$ (năng lượng giảm khi đồng phân hóa), làm mất ổn định hệ vòng phẳng, kết luận là chất phản thơm.

3. Việc thay đổi các bộ hàm cơ sở (3-21G, 3-21G*, 6-31G, 6-31G*) ảnh hưởng ra sao đến kết quả tính toán?

Khi nâng cấp từ bộ cơ sở tối thiểu/tách đôi (3-21G, 6-31G) lên các bộ cơ sở có thêm hàm phân cực $d$ trên nguyên tử nặng (3-21G*, 6-31G*), giá trị tuyệt đối của năng lượng toàn phần và $\Delta E$ có sự thay đổi định lượng để phản ánh chính xác hơn sự biến dạng của mật độ đám mây electron. Tuy nhiên, dấu của $\Delta E$ và tương quan trật tự năng lượng không đổi. Do đó, việc thay đổi bộ cơ sở không làm thay đổi kết luận định tính về bản chất thơm, phản thơm hay không thơm của các chất khảo sát.

4. Tại sao quy tắc Hückel không áp dụng được một cách chính xác cho [10]-anulen và [14]-anulen?

Theo phân tích trong Chương 1, dù $[10]$-anulen và $[14]$-anulen có số electron $\pi$ thỏa mãn công thức $(4n+2)$ lần lượt với $n=2$ (10 electron) và $n=3$ (14 electron), chúng không phẳng:

  • Ở $[10]$-anulen (đồng phân cis-cis-trans-cis-trans), sức căng góc liên kết lớn và tương tác đẩy không gian mạnh giữa 2 nguyên tử hydro nằm phía trong vòng khiến phân tử bị vặn xoắn lệch khỏi mặt phẳng.
  • Ở $[14]$-anulen (đồng phân cis-trans-trans-cis-trans-cis-trans), các nguyên tử hydro nội vòng tại vị trí cacbon số 3, 6, 10, 13 đẩy nhau, làm mất tính song song của các orbital $p$. Do thiếu điều kiện cấu trúc phẳng, hệ electron $\pi$ không thể giải tỏa liên tục, khiến các chất này kém bền và không thể hiện tính thơm điển hình.

5. Tại sao 1-metyl xicloocta-1,3,5,7-tetraen lại được xếp vào nhóm hợp chất không thơm?

Khóa luận chỉ ra rằng 1-metyl xicloocta-1,3,5,7-tetraen có 8 electron $\pi$ ($4n$ electron). Ở trạng thái tối ưu hóa hình học, phân tử không duy trì dạng đa giác đều phẳng mà chuyển sang cấu dạng thuyền (tub conformation). Cấu dạng này làm mất tính song song và xen phủ liên tục của các orbital $p$, biến phân tử thành một polyen vòng không no thông thường. Năng lượng phản ứng đồng phân hóa tính được rất nhỏ ($\Delta E = 0,000856\text{ Hartree} \approx 0,54\text{ kcal/mol}$), nằm trọn trong ngưỡng dao động không thơm ($\pm 5\text{ kcal/mol}$).


Kết luận

Khóa luận tốt nghiệp của tác giả Nguyễn Thị Thanh Thúy đã ứng dụng thành công các công cụ hóa tính toán (GaussView 5.0 và Gaussian 09) để khảo sát và định lượng hóa tính thơm của 15 hợp chất hữu cơ thông qua biến thiên năng lượng của phản ứng đồng phân hóa. Việc kết hợp phương pháp tối ưu hóa hình học bán thực nghiệm AM1 và tính toán năng lượng ab initio RHF/3-21G* mang lại kết quả phân loại chuẩn xác, phù hợp với các quy luật nhiệt động học thực nghiệm. Công trình là tài liệu tham khảo học thuật hữu ích về phương pháp nghiên cứu lượng tử ứng dụng cho sinh viên, học viên và giảng viên chuyên ngành Hóa học.