Tổng quan nghiên cứu

Trong hóa học lập thể hiện đại, sự phân bố cấu trạng phân tử đóng vai trò then chốt quyết định độ bền nhiệt động, tính chất hóa lý và hoạt tính sinh học của vật chất. Theo nguyên lý lập thể kinh điển, phân tử etan có rào thế cản quay khoảng 2,8 kcal/mol, trong khi n-butan thể hiện độ chênh lệch thế năng 0,8 kcal/mol nghiêng về dạng đối lệch (Anti). Cấu trạng Anti với góc nhị diện 180 độ thường bền nhất do triệt tiêu tối đa lực đẩy không gian Van der Waals và tương tác tĩnh điện giữa các nhóm thế. Dẫu vậy, nhiều khảo sát thực nghiệm ghi nhận nghịch lý khi cấu trạng bán lệch (Gauche) với góc nhị diện 60 độ lại chiếm tỉ lệ vượt trội, tạo nên hiệu ứng bán lệch (Gauche effect). Luận văn đi sâu phân tích bản chất cơ chế này bằng phương pháp hóa học lượng tử trên 7 dẫn xuất 1,2-thế của etan, bao gồm các hợp chất đối xứng và bất đối xứng chứa nhóm thế halogen hoặc hydroxyl. Phạm vi nghiên cứu được khảo sát tại 3 mốc nhiệt độ tiêu chuẩn gồm 0 độ C, 25 độ C và 100 độ C, đồng thời đánh giá sự chuyển dịch cấu trạng trong 2 dung môi đặc trưng là nước (hằng số điện môi khoảng 78,4) và cacbon tetraclorua (hằng số điện môi khoảng 2,2). Mục tiêu trọng tâm là tối ưu hóa hình học phân tử, tính toán chính xác năng lượng tự do Gibbs và xác định hằng số cân bằng nhiệt động để giải mã bản chất lực tương tác vi mô, cung cấp cơ sở dữ liệu quan trọng phục vụ công tác thiết kế phân tử trong hóa dược và tổng hợp hữu cơ.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng của thuyết obitan phân tử (MO) kết hợp lý thuyết trường tự hợp Hartree-Fock (HF) theo phép gần đúng Born-Oppenheimer. Nhằm xử lý hiệu ứng tương quan electron phức tạp mà phương pháp Hartree-Fock thuần túy bỏ qua, lý thuyết hàm mật độ (DFT) với phiếm hàm lai hóa ba tham số B3LYP của Becke được lựa chọn làm công cụ tính toán cốt lõi. Mô hình B3LYP phối hợp hài hòa năng lượng trao đổi chính xác từ thuyết Hartree-Fock với hàm tương quan gradient do Lee, Yang và Parr đề xuất. Ba khái niệm trung tâm được vận dụng xuyên suốt gồm: góc nhị diện trong công thức chiếu Newmann, tương tác siêu liên hợp donor-acceptor giữa obitan liên kết sigma và obitan phản liên kết sigma sao, cùng hệ thức nhiệt động học xác định hằng số cân bằng hóa học qua biến thiên năng lượng tự do Gibbs.

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

Cỡ mẫu nghiên cứu gồm 7 hợp chất dẫn xuất 1,2-thế tiêu biểu của etan, được phân chia thành hai nhóm chính: nhóm dihalogen đối xứng (1,2-difloroetan, 1,2-dicloroetan, 1,2-dibromoetan) và nhóm chứa chức hydroxyl (etylen glycol, 2-floroetanol, 2-cloroetanol, 2-bromoetanol). Phương pháp chọn mẫu có chủ đích được áp dụng nhằm bao quát quy luật biến thiên liên tục theo độ âm điện (từ flo 3,98 đến brom 2,96 theo thang Pauling), bán kính Van der Waals và khả năng hình thành liên kết hydro nội phân tử. Toàn bộ quá trình tính toán nhiệt động và tối ưu hóa hình học không gian được thực hiện trên phần mềm Gaussian 98W với các tập cơ sở phân tách hóa trị phân cực 6-31G(d) và khuếch tán 6-31+G(d). Lý do lựa chọn tổ hợp phương pháp DFT/B3LYP cùng các tập cơ sở này là khả năng mô tả chính xác tương tác liên kết hydro yếu và phân tán lưỡng cực, đồng thời tối ưu hóa chi phí và thời gian tính toán. Lộ trình tính toán và phân tích dữ liệu được triển khai đồng bộ tại phòng tính toán Hóa lý thuộc Trường Đại học Sư phạm Thành phố Hồ Chí Minh.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Thứ nhất, đối với dãy dihalogenoetan trong pha khí tại 25 độ C, 1,2-dicloroetan và 1,2-dibromoetan ưu tiên cấu trạng Anti với tỉ lệ chiếm từ 75% đến 82%, trong khi 1,2-difloroetan thể hiện hiệu ứng bán lệch rõ rệt với dạng Gauche chiếm ưu thế khoảng 58% đến 65%. Thứ hai, ở các dẫn xuất chứa nhóm hydroxyl như etylen glycol và 2-floroetanol, cấu trạng Gauche chiếm ưu thế áp đảo từ 88% đến 94% nhờ sự xuất hiện của liên kết hydro nội phân tử tạo vòng giả 5 cạnh bền vững có năng lượng liên kết ước tính 2,5 đến 3,8 kcal/mol. Thứ ba, môi trường dung môi phân cực nước làm gia tăng tỉ lệ cấu trạng Gauche thêm khoảng 15% đến 22% so với pha khí và dung môi kém phân cực cacbon tetraclorua, do mômen lưỡng cực dạng Gauche (dao động từ 2,1 đến 3,2 Debye) tương tác solvat hóa mạnh hơn dạng Anti có mômen lưỡng cực xấp xỉ 0 Debye. Khi nhiệt độ tăng từ 0 độ C lên 100 độ C, tỉ lệ dạng Anti tăng nhẹ từ 3% đến 6% do hiệu ứng entropy thúc đẩy chuyển động quay nhiệt hỗn loạn.

Thảo luận kết quả

Bản chất của hiệu ứng bán lệch trong 1,2-difloroetan bắt nguồn từ sự xen phủ obitan siêu liên hợp giữa liên kết sigma C-H và obitan phản liên kết sigma sao C-F có năng lượng thấp, bù trừ hoàn toàn lực đẩy tĩnh điện giữa hai nguyên tử flo. Dữ liệu tính toán được biểu diễn trực quan qua biểu đồ đường cong thế năng quét góc quay 360 độ quanh trục liên kết cacbon-cacbon, xác định rõ hai cực tiểu thế năng tại 60 độ và 180 độ cùng các đỉnh rào cản năng lượng từ 3,5 đến 9,2 kcal/mol tại 0 độ và 120 độ. Đồng thời, bảng đối chiếu năng lượng tự do Gibbs và mômen lưỡng cực giúp minh chứng tính ổn định của các cấu trúc. Kết quả này hoàn toàn tương thích với các công bố thực nghiệm của David O'Hagan năm 2000 trên dẫn xuất fluoroamid và nghiên cứu của Allen năm 2004 trên 2-floroetanol, khẳng định tính chuẩn xác và độ tin cậy cao của mô hình cơ lượng tử đã thiết lập.

Đề xuất và khuyến nghị

Thứ nhất, mở rộng mô hình tính toán lượng tử sang các hệ phân tử dược chất đa chức phức tạp và các dẫn xuất peptit tự nhiên, hướng tới mục tiêu phân tích cấu trạng ưu đãi cho hơn 50 hợp chất có hoạt tính sinh học trong giai đoạn 2026-2027 do các viện nghiên cứu hóa dược chủ trì.

Thứ hai, nâng cấp phương pháp cơ lượng tử bằng cách áp dụng các tập cơ sở mở rộng cao cấp hơn như aug-cc-pVTZ kết hợp phương pháp tương quan điện tử CCSD(T) hoặc CBS-QB3, nhằm giảm sai số năng lượng tự do Gibbs xuống dưới ngưỡng 0,1 kcal/mol trong vòng 12 tháng tới.

Thứ ba, tích hợp mô hình solvat hóa kết hợp giữa trường điện môi liên tục PCM và các phân tử dung môi tường minh xung quanh vị trí tạo liên kết hydro, nâng cao độ tương thích với thực nghiệm quang phổ cộng hưởng từ hạt nhân NMR và phổ hồng ngoại IR lên trên 95% trước năm 2028.

Thứ tư, chuẩn hóa giáo trình thực hành hóa tính toán tại các trường đại học sư phạm và bách khoa, đảm bảo 100% sinh viên và học viên cao học chuyên ngành Hóa lý được đào tạo chuyên sâu về kỹ năng mô phỏng trên các phần mềm Gaussian và ORCA trước năm học tới.

Đối tượng nên tham khảo luận văn

Nhóm thứ nhất là giảng viên và nghiên cứu viên Hóa lý, Hóa hữu cơ và Hóa học tính toán, sử dụng luận văn làm tài liệu chuyên khảo để kiểm chứng tương tác xen phủ obitan và hoàn thiện các bài giảng chuyên đề hóa lượng tử.

Nhóm thứ hai là học viên cao học và nghiên cứu sinh ngành Hóa học, tham khảo quy trình chuẩn mực về tối ưu hóa hình học phân tử, xử lý dữ liệu nhiệt động và phân tích rào cản năng lượng quay từ phần mềm Gaussian.

Nhóm thứ ba là chuyên gia nghiên cứu phát triển R&D trong ngành dược phẩm và vật liệu mới, khai thác quy luật cấu trạng ưu đãi để thiết kế cấu trúc phân tử thuốc có độ gắn kết tối ưu với thụ thể sinh học.

Nhóm thứ tư là sinh viên chuyên ngành Sư phạm Hóa học và Hóa học ứng dụng, củng cố kiến thức nền tảng về hóa học lập thể, đồng phân quay và các tương tác phi cổ điển chi phối độ bền phân tử.

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

Hiệu ứng bán lệch là gì và nguyên nhân vi mô nào chi phối hiện tượng này?

Hiệu ứng bán lệch là hiện tượng cấu trạng Gauche với góc nhị diện khoảng 60 độ trở nên bền hơn và chiếm tỉ lệ cao hơn cấu trạng Anti trong hỗn hợp cân bằng. Hiện tượng này xuất phát từ tương tác siêu liên hợp bền hóa giữa obitan liên kết sigma C-H và obitan phản liên kết sigma sao C-X, kết hợp với liên kết hydro nội phân tử trong các dẫn xuất chứa nhóm hydroxyl.

Tại sao phương pháp DFT với phiếm hàm B3LYP lại được lựa chọn để khảo sát các dẫn xuất etan?

Phiếm hàm lai hóa B3LYP kết hợp năng lượng trao đổi Hartree-Fock với hàm tương quan gradient, mang lại độ chính xác cao khi xác định cấu trúc hình học và năng lượng dao động của phân tử. Phương pháp này mô tả chính xác tương tác Van der Waals và liên kết hydro yếu với sai số dưới 1,0 kcal/mol, đồng thời tối ưu hóa thời gian tính toán.

Môi trường dung môi phân cực tác động thế nào đến sự chuyển dịch cân bằng cấu trạng?

Dung môi có hằng số điện môi cao như nước hạ thấp thế năng của cấu trạng Gauche thông qua tương tác solvat hóa lưỡng cực mạnh, làm tăng tỉ lệ dạng Gauche thêm khoảng 15% đến 22% so với pha khí. Trong khi đó, dung môi không phân cực như cacbon tetraclorua hầu như không làm biến đổi tương quan năng lượng gốc giữa hai cấu trạng.

Liên kết hydro nội phân tử đóng vai trò gì trong việc ổn định cấu trạng Gauche của 2-floroetanol?

Trong phân tử 2-floroetanol, cấu trạng Gauche đưa nguyên tử hydro linh động của nhóm OH lại gần nguyên tử flo âm điện, tạo nên liên kết hydro nội phân tử dạng vòng 5 cạnh. Tương tác này đóng góp năng lượng bền hóa khoảng 2,8 kcal/mol, giúp cấu trạng Gauche chiếm tỉ lệ áp đảo trên 90% ở nhiệt độ phòng.

Sự thay đổi nhiệt độ từ 0 độ C đến 100 độ C ảnh hưởng ra sao đến tỉ lệ phân bố cấu trạng?

Khi nhiệt độ tăng từ 0 độ C lên 100 độ C, số hạng entropy trong phương trình Gibbs tăng lên, kích thích phân tử vượt qua rào cản năng lượng quay từ 3,0 đến 4,5 kcal/mol. Do đó, tỉ lệ cấu trạng Anti vốn kém bền hơn về mặt thế năng trong hệ có hiệu ứng bán lệch sẽ tăng nhẹ từ 3% đến 6% trong hỗn hợp.

Kết luận

Luận văn đã hoàn thành xuất sắc việc khảo sát hiệu ứng bán lệch trên 7 dẫn xuất của etan bằng phần mềm Gaussian 98W với 5 kết luận trọng tâm:

  • Xác lập quy luật chuyển dịch ưu thế cấu trạng từ Anti sang Gauche khi thay thế nhóm halogen từ brom, clo sang flo.
  • Làm sáng tỏ vai trò quyết định của liên kết hydro nội phân tử trong etylen glycol và 2-floroetanol với tỉ lệ dạng Gauche đạt trên 88%.
  • Chứng minh cơ chế solvat hóa của dung môi nước làm tăng độ bền dạng Gauche thêm khoảng 15% đến 22% nhờ mômen lưỡng cực cao.
  • Xây dựng hệ thống dữ liệu nhiệt động học định lượng mối quan hệ giữa năng lượng tự do Gibbs và hằng số cân bằng tại 3 mốc nhiệt độ.
  • Đóng góp nguồn tài liệu tham khảo cơ lượng tử chuẩn mực cho nghiên cứu hóa học lập thể và hóa học tính toán tại Việt Nam.

Trong giai đoạn 2026-2028, các nghiên cứu tiếp theo nên mở rộng sang mô phỏng động lực học phân tử trên các hệ đa chức phức tạp. Quý độc giả và các nhà khoa học hãy ứng dụng quy trình tính toán Gaussian để tối ưu hóa việc nghiên cứu cấu trúc phân tử và phát triển các hợp chất mới.