Tổng quan nghiên cứu

Nghiên cứu về tương tác không cộng hóa trị, đặc biệt là liên kết hydro, đóng vai trò nền tảng trong hóa học hiện đại với hơn 90% các quá trình sinh hóa và cấu trúc phân tử sinh học chịu sự chi phối trực tiếp từ dạng liên kết này. Luận văn thạc sĩ chuyên ngành Hóa lý thuyết và hóa lý của tác giả Nguyễn Thị Thanh Cúc, thực hiện tại Trường Đại học Quy Nhơn vào năm 2021 dưới sự hướng dẫn của PGS. Nguyễn Tiến Trung, tập trung làm sáng tỏ bản chất của liên kết hydro chuyển dời xanh và tương tác chalcogen mới. Đối tượng nghiên cứu bao gồm 12 hệ phức cấu trúc được hình thành giữa monomer CH3CHZ (với Z lần lượt là các nguyên tố chalcogen: oxy, lưu huỳnh, selen, telur) và các phân tử nước nH2O (với n = 1, 2, 3).

Vấn đề then chốt được đặt ra là giải quyết các tranh luận lý thuyết kéo dài nhiều thập kỷ về cơ chế của liên kết hydro chuyển dời xanh tại vị trí Csp2–H, đồng thời đánh giá vai trò của các nguyên tử chalcogen nặng như Se và Te vốn rất hiếm khi được khảo sát có hệ thống. Mục tiêu cụ thể của công trình là xác định cấu trúc hình học cân bằng, định lượng độ bền tương tác thông qua mức lý thuyết lượng tử cao cấp CCSD(T)/6-311++G(3df,2pd)//MP2/6-311++G(3df,2pd), phân tích bản chất liên kết bằng thuyết AIM, NBO, SAPT2+ và làm rõ hiệu ứng cộng kết khi tăng số lượng phân tử nước từ 1 lên 3. Nghiên cứu thực hiện trong phạm vi hệ pha khí tại phòng thí nghiệm Hóa học tính toán và Mô phỏng trong năm 2021. Kết quả đạt được mang ý nghĩa khoa học sâu sắc, bổ sung dữ liệu nhiệt động học chuẩn xác với năng lượng tương tác hiệu chỉnh dao động từ khoảng -15 kJ/mol đến hơn -60 kJ/mol, tạo tiền đề định lượng vững chắc cho việc thiết kế cấu trúc phân tử và mô phỏng phản ứng sinh hóa.

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

Khung lý thuyết áp dụng

Công trình vận dụng toàn diện hệ thống lý thuyết hóa học lượng tử hiện đại nhằm mô tả bản chất tương tác giữa các phân tử:

  • Thuyết nhiễu luận Møller-Plesset và phương pháp chùm tương tác (Coupled Cluster): Sử dụng phương pháp nhiễu loạn bậc hai MP2 để thu hồi khoảng 80% đến 90% năng lượng tương quan electron khi tối ưu hóa hình học, kết hợp phương pháp CCSD(T) được xem là "tiêu chuẩn vàng" để tính toán năng lượng điểm đơn có độ chính xác cao.
  • Thuyết phiếm hàm mật độ (DFT) kết hợp hiệu chỉnh phân tán: Áp dụng phiếm hàm ωB97X-D nhằm xử lý tương tác siêu liên hợp và hiệu ứng phân tán tầm xa.
  • Thuyết nguyên tử trong phân tử (AIM) của Richard Bader: Khảo sát các điểm tới hạn liên kết (BCP), mật độ electron ρ(r), toán tử Laplacian ∇²ρ(r), mật độ động năng G(r) và mật độ thế năng V(r) để xác định bản chất liên kết.
  • Thuyết nhiễu loạn phù hợp đối xứng (SAPT2+): Tách tổng năng lượng tương tác thành 4 thành phần vật lý tường minh: năng lượng tĩnh điện (Eelst), năng lượng trao đổi (Eexch), năng lượng cảm ứng (Eind) và năng lượng phân tán (Edisp).
  • Thuyết orbital liên kết tự nhiên (NBO) và phân tích vùng tương tác không cộng hóa trị (NCIPLOT): Định lượng sự chuyển dịch mật độ electron (EDT) giữa các orbital liên kết và phản liên kết σ*, giải thích hiện tượng tái lai hóa trên nguyên tử cacbon.

Các khái niệm cốt lõi bao gồm liên kết hydro chuyển dời đỏ truyền thống, liên kết hydro chuyển dời xanh (Blue-Shifting Hydrogen Bond - BSHB), enthalpy tách proton (DPE), ái lực proton (PA) và sai số do chồng chất bộ cơ sở (BSSE) được hiệu chỉnh theo sơ đồ Boys - Bernardi.

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

Tập dữ liệu nghiên cứu bao gồm 4 monomer aldehyde và dẫn xuất chalcogenaldehyde tương ứng (CH3CHO, CH3CHS, CH3CHSe, CH3CHTe) cùng các cụm phức chất với 1, 2 và 3 phân tử nước. Toàn bộ cấu trúc được tối ưu hóa hoàn toàn ở mức lý thuyết MP2 sử dụng bộ hàm cơ sở mở rộng 6-311++G(3df,2pd) tích hợp đầy đủ hàm khuếch tán và hàm phân cực bậc cao.

Lý do lựa chọn phương pháp phân tích này xuất phát từ yêu cầu khắt khe về độ chính xác khi mô tả các liên kết hydro yếu và liên kết chuyển dời xanh, nơi mà các phương pháp DFT thông thường hoặc Hartree-Fock thường đánh giá sai lệch năng lượng tương quan. Năng lượng dao động điểm không (ZPE) và sai số BSSE được tính toán cẩn trọng để đảm bảo giá trị nhiệt động thu được có độ tin cậy tương đương dữ liệu thực nghiệm phổ hồng ngoại (IR) và phổ cộng hưởng từ hạt nhân (NMR). Quá trình tính toán phân tích SAPT2+ được tiến hành trên phần mềm Psi4 với bộ cơ sở def2-TZVPP, trong khi phân tích mật độ topo electron được thực hiện qua phần mềm AIMAll và NBO 5.0 trong suốt tiến trình nghiên cứu năm 2021.

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

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

Nghiên cứu đã ghi nhận 4 phát hiện quan trọng có giá trị khoa học vượt bậc:

  1. Quy luật hình thành liên kết hydro dị thể O–H···Z (Z=O, S, Se, Te): Tất cả các phức đều hình thành liên kết hydro dị thể bền vững giữa nhóm O–H của nước và nguyên tử chalcogen Z. Độ bền liên kết giảm dần theo thứ tự từ O > S > Se > Te, tương ứng với sự suy giảm ái lực proton (PA) của các nguyên tử trung tâm từ oxy xuống telur trong bảng tuần hoàn.
  2. Xác nhận bản chất chuyển dời xanh của liên kết Csp2–H···O: Khi tạo phức, liên kết Csp2–H thể hiện sự co ngắn độ dài liên kết từ 0,002 Å đến 0,012 Å và tần số dao động hóa trị tăng vọt (chuyển dịch xanh) với biên độ từ khoảng 15 cm⁻¹ lên tới hơn 85 cm⁻¹ so với monomer tự do.
  3. Hiệu ứng cộng kết vượt trội khi gia tăng số phân tử nước: Khi số lượng phân tử nước tăng từ n = 1 lên n = 3, năng lượng tương tác toàn phần của hệ phức tăng từ khoảng -18 kJ/mol lên vượt mức -65 kJ/mol. Hiệu ứng cộng kết giữa các liên kết hydro vòng làm tăng độ bền chung của hệ thêm khoảng 15% đến 28% so với tổng năng lượng của từng liên kết đơn lẻ cấu thành.
  4. Phân tích topo mật độ electron tại điểm tới hạn (BCP): Kết quả AIM chỉ ra mật độ electron ρ(r) tại các điểm BCP của liên kết Csp2–H···O nằm trong khoảng 0,010 đến 0,025 au với Laplacian mang giá trị dương (∇²ρ(r) > 0), chứng minh đây là các tương tác đóng vỏ mang bản chất tĩnh điện ưu thế.

Thảo luận kết quả

Cơ chế chuyển dời xanh của liên kết Csp2–H được làm sáng tỏ thông qua sự phối hợp giữa phân tích NBO và SAPT2+. Khi phức hình thành, sự chuyển dịch mật độ electron ngoại phân tử từ đôi điện tử tự do của oxy trong nước không đi thẳng vào orbital phản liên kết σ*(C–H) mà được định hướng vào các vùng lân cận, kích hoạt quá trình tái lai hóa làm gia tăng tỷ lệ phần trăm đặc tính orbital s (%s) trong liên kết C–H. Sự gia tăng đặc tính s làm liên kết C–H trở nên ngắn hơn và bền hơn, trực tiếp dẫn đến hiện tượng dịch chuyển tần số dao động về vùng sóng ngắn (blue-shift).

Dữ liệu nghiên cứu có thể được trực quan hóa rất sinh động thông qua biểu đồ phân tán 2D thể hiện tương quan tuyến tính chặt chẽ giữa mật độ electron tại BCP và năng lượng liên kết riêng phần EHB, cũng như biểu đồ cột chồng mô tả tỷ lệ phần trăm đóng góp của các thành phần năng lượng SAPT2+. Trong đó, năng lượng tĩnh điện đóng góp chi phối từ 50% đến 65% lực hút làm bền, tiếp theo là năng lượng cảm ứng chiếm 20% đến 30% và năng lượng phân tán chiếm khoảng 15% đến 20%. So sánh với các công bố trước đây trên hệ formaldehyde và acid formic ghi nhận độ chuyển dịch xanh từ 81 đến 96 cm⁻¹, kết quả của luận văn thể hiện sự tương thích hoàn hảo, đồng thời mở rộng biên giới hiểu biết sang các dị tố Se và Te vốn chưa từng được định lượng chi tiết.

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

Nhằm phát huy tối đa các phát hiện lý thuyết của luận văn vào thực tiễn nghiên cứu hóa học và công nghệ vật liệu, 4 nhóm giải pháp chiến lược được đề xuất:

  1. Mở rộng quy mô mô phỏng động lực học phân tử lượng tử (QM/MM): Các nhóm nghiên cứu hóa lý tính toán cần triển khai mô phỏng solvat hóa mở rộng với số lượng phân tử nước từ 4 đến 20 phân tử trong giai đoạn 2022-2024, nhằm tiệm cận môi trường dung dịch thực tế và kiểm chứng tính bền của hiệu ứng cộng kết với sai số nhiệt động mục tiêu dưới 1,0 kJ/mol.
  2. Ứng dụng quy luật chuyển dời xanh trong thiết kế xúc tác hữu cơ: Các nhà hóa học tổng hợp hữu cơ nên ứng dụng cơ chế hoạt hóa liên kết Csp2–H thông qua tương tác hydro chuyển dời xanh để phát triển các tâm xúc tác chọn lọc lập thể mới, hướng đến mục tiêu gia tăng hiệu suất phản ứng tổng hợp từ 15% đến 25% trong các quy trình hóa học xanh từ năm 2023.
  3. Xây dựng cơ sở dữ liệu phổ chuẩn cho hợp chất chứa chalcogen nặng: Các phòng thí nghiệm phân tích quang phổ cần phối hợp thiết lập thư viện dữ liệu đối soát phổ hồng ngoại IR và phổ NMR tính toán cho hơn 50 hợp chất hữu cơ chứa selen và telur vào năm 2024, hỗ trợ nhận diện nhanh các cấu trúc trung gian trong hóa dược.
  4. Chuẩn hóa quy trình tính toán năng lượng tương tác không cộng hóa trị: Các viện nghiên cứu và trường đại học cần tích hợp quy trình phối hợp MP2, CCSD(T), AIM và SAPT2+ vào các chương trình đào tạo sau đại học từ học kỳ năm 2023, nâng cao tỷ lệ công bố quốc tế chuẩn ISI/Scopus của học viên lên trên 40%.

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

Công trình luận văn là tài liệu tham khảo chuyên sâu đặc biệt giá trị cho 4 nhóm đối tượng sau:

  • Học viên cao học và nghiên cứu sinh chuyên ngành Hóa lý thuyết, Hóa lý, Hóa học tính toán: Tiếp cận phương pháp luận chuẩn mực về tối ưu hóa hình học phân tử, xử lý sai số BSSE và áp dụng thành thạo các gói phần mềm hiện đại như Gaussian, Psi4, AIMAll trong việc xử lý dữ liệu phức chất.
  • Giảng viên và chuyên gia nghiên cứu cấu trúc phân tử: Sử dụng tài liệu như một nguồn tham khảo bài giảng chi tiết về thuyết tương tác liên phân tử, cơ chế liên kết hydro chuyển dời xanh và hóa học các hợp chất nhóm chalcogen.
  • Nhà khoa học phát triển dược phẩm và hóa sinh học cấu trúc: Khai thác các dữ liệu về độ bền liên kết O–H···Se và O–H···Te để mô phỏng tương tác giữa các phối tử chứa selen/telur với các thụ thể protein đích trong cơ thể.
  • Kỹ sư công nghệ vật liệu tinh thể và hóa học siêu phân tử: Ứng dụng các thông số năng lượng tương tác định lượng để định hướng thiết kế mạng lưới liên kết hydro đa chiều trong vật liệu tự lắp ráp và khung hữu cơ kim loại.

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

1. Liên kết hydro chuyển dời xanh (Blue-shifting hydrogen bond) có điểm gì khác biệt so với liên kết hydro truyền thống? Liên kết hydro truyền thống (chuyển dời đỏ) làm giãn dài liên kết cho proton A–H và giảm tần số dao động hóa trị. Ngược lại, liên kết hydro chuyển dời xanh làm co ngắn độ dài liên kết C–H từ 0,002 đến 0,012 Å và làm tăng tần số dao động từ khoảng 15 cm⁻¹ lên tới hơn 85 cm⁻¹ do hiện tượng tái lai hóa orbital.

2. Tại sao nghiên cứu cần khảo sát các nguyên tố chalcogen nặng như Selen (Se) và Telur (Te)? Các nguyên tố Se và Te có bán kính nguyên tử lớn, độ âm điện thấp và mật độ điện tích khuếch tán, giữ vai trò quan trọng trong các enzyme sinh học nhưng hiếm khi được khảo sát. Việc khảo sát giúp kiểm chứng tính tổng quát của thuyết liên kết hydro và định lượng ái lực proton thực tế.

3. Hiệu ứng cộng kết (cooperativity) thể hiện vai trò gì trong hệ phức CH3CHZ···nH2O? Khi tăng số phân tử nước từ 1 lên 3, các liên kết hydro tạo thành hệ thống vòng khép kín hỗ trợ lẫn nhau. Hiệu ứng này làm năng lượng tương tác toàn phần tăng phi tuyến tính, đóng góp thêm từ 15% đến 28% vào tổng độ bền nhiệt động học của toàn hệ phức.

4. Vì sao cần kết hợp nhiều phương pháp tính toán như MP2, CCSD(T), AIM và SAPT2+? Mỗi phương pháp cung cấp một khía cạnh bổ trợ: MP2 tối ưu hình học tin cậy, CCSD(T) cung cấp năng lượng chuẩn xác, AIM định lượng mật độ electron tại điểm tới hạn BCP, và SAPT2+ phân rã tường minh các thành phần lực tĩnh điện, phân tán và cảm ứng.

5. Kết quả của luận văn có thể ứng dụng trực tiếp vào nghiên cứu thực nghiệm như thế nào? Dữ liệu dịch chuyển tần số dao động Δν và sự thay đổi độ dài liên kết Δr cung cấp thang đo chuẩn xác để các nhà thực nghiệm đối chiếu, giải đoán phổ hồng ngoại (FT-IR) và phổ Raman của các hợp chất trung gian trong phản ứng hóa học.

Kết luận

  • Luận văn đã mô phỏng và tối ưu hóa thành công cấu trúc không gian của 12 hệ phức chất CH3CHZ···nH2O (Z=O, S, Se, Te; n=1, 2, 3) tại mức lý thuyết chất lượng cao MP2/6-311++G(3df,2pd).
  • Chứng minh vững chắc sự tồn tại của liên kết hydro chuyển dời xanh tại vị trí Csp2–H···O với mức tăng tần số dao động lên đến hơn 85 cm⁻¹ và co ngắn độ dài liên kết C–H tương ứng.
  • Định lượng độ suy giảm độ bền của liên kết hydro dị thể O–H···Z theo trật tự tuần hoàn O > S > Se > Te, phù hợp hoàn toàn với quy luật biến thiên ái lực proton và mật độ điện tích topo AIM.
  • Xác lập vai trò chi phối của lực hút tĩnh điện (50-65%) thông qua phân tích SAPT2+ và làm rõ cơ chế tăng cường độ bền từ 15-28% nhờ hiệu ứng cộng kết đa phân tử nước.
  • Thiết lập lộ trình mở rộng nghiên cứu mô phỏng động lực học lượng tử trong dung môi thực giai đoạn 2022-2025, mở ra hướng ứng dụng đột phá trong thiết kế xúc tác và công nghệ hóa dược.

Quý độc giả, nghiên cứu sinh và các nhà khoa học quan tâm đến lĩnh vực hóa học tính toán có thể khai thác toàn bộ dữ liệu tọa độ hình học và ma trận năng lượng chi tiết của công trình để phục vụ cho các nghiên cứu chuyên sâu tiếp theo.