Tổng quan nghiên cứu

Liên kết hydrogen đóng vai trò cốt lõi trong việc duy trì cấu trúc bậc cao của hơn 90% các phân tử sinh học quan trọng như DNA, RNA và protein. Trong khi liên kết hydrogen truyền thống kiểu chuyển dời đỏ đã được khám phá từ năm 1931 với cơ chế tĩnh điện kinh điển, dạng liên kết hydrogen chuyển dời xanh với sự tham gia của liên kết Csp2-H vẫn còn nhiều bí ẩn chưa được giải mã trọn vẹn. Luận văn thạc sĩ chuyên ngành Hóa lý thuyết và hóa lý được thực hiện tại Trường Đại học Quy Nhơn vào năm 2022 tập trung giải quyết vấn đề tương tác phi cộng hóa trị đặc biệt này. Mục tiêu nghiên cứu cụ thể là khảo sát cấu trúc không gian, đánh giá độ bền năng lượng và cơ chế chuyển dịch phổ của 24 hệ dimer RCHZ...RCHZ được tạo bởi 6 nhóm thế hóa học và 4 nguyên tố chalcogen gồm oxy, lưu huỳnh, selen và telu. Về mặt phạm vi, đề tài áp dụng các thuật toán hóa học lượng tử tiên tiến để mô phỏng tương tác trên không gian ba chiều trong điều kiện pha khí chuẩn. Kết quả từ công trình mang ý nghĩa học thuật vượt trội khi cung cấp hệ thống dữ liệu toàn diện với độ chính xác cao về 24 điểm cực tiểu thế năng, 48 thông số enthalpy tách proton và ái lực proton. Việc làm sáng tỏ quy luật tương tác giữa các tâm chalcogen nặng không chỉ lấp đầy khoảng trống dữ liệu thực nghiệm mà còn đóng góp các giá trị định lượng chuẩn xác nhằm định hướng cho quá trình mô hình hóa cấu trúc phân tử thuốc và kỹ thuật tinh thể trong giai đoạn tiếp theo.

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

Khung lý thuyết áp dụng

Luận văn vận dụng phương trình sóng Schrödinger kết hợp phép gần đúng Born-Oppenheimer nhằm phân lập chuyển động của electron trong trường thế hạt nhân cố định. Thuyết nhiễu loạn bậc hai Møller-Plesset và thuyết chùm tương tác kích thích đơn, đôi kết hợp kích thích ba chuẩn hóa được sử dụng để xử lý năng lượng tương quan electron, vốn chiếm tới khoảng 80% đến 90% độ chính xác của tương tác yếu. Bốn khái niệm cốt lõi được tích hợp xuyên suốt mô hình gồm: năng lượng tương tác hiệu chỉnh sai số chồng chất bộ cơ sở, điểm tới hạn liên kết trong thuyết nguyên tử trong phân tử, ái lực proton và enthalpy tách proton. Đặc biệt, thuyết nguyên tử trong phân tử của Bader cùng phân tích orbital liên kết tự nhiên đóng vai trò khung phân tích hình học topo và bản chất truyền điện tích nội phân tử. Các mô hình lý thuyết này cho phép mô tả chính xác sự phân bố mật độ electron và định vị chính xác vị trí không gian của các liên kết phi cộng hóa trị.

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

Nghiên cứu thu thập dữ liệu tính toán từ tập hợp mẫu gồm 24 cấu hình dimer và 24 phân tử monomer được chọn lọc theo nguyên tắc biến đổi tuần hoàn hệ thống từ 6 nhóm thế gồm H, F, Cl, Br, CH3, NH2 kết hợp 4 tâm chalcogen là O, S, Se, Te. Quy trình chọn mẫu đảm bảo tính đại diện cho cả nhóm hút electron và đẩy electron nhằm phản ánh đầy đủ phổ biến thiên mật độ điện tích. Phương pháp nhiễu loạn bậc hai với bộ hàm 6-311++G(3df,2pd) và aug-cc-pVDZ-PP được chọn để tối ưu hóa hình học vì tối ưu được thời gian tính toán trong khi vẫn đảm bảo độ hội tụ cao cho các nguyên tử nặng. Năng lượng điểm đơn được chuẩn hóa ở mức chùm tương tác CCSD(T) nhằm triệt tiêu sai số năng lượng dao động điểm không và hiệu chỉnh chồng chất bộ cơ sở theo kỹ thuật đối trọng Boys-Bernardi. Toàn bộ chuỗi tính toán mô phỏng được thực hiện trên hệ thống phần mềm Gaussian 16 và AIMALL trong khung thời gian nghiên cứu liên tục kéo dài 12 tháng tại phòng thí nghiệm chuyên ngành.

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

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

Thứ nhất, toàn bộ 24 dimer đều tự hình thành cấu trúc vòng 6 cạnh phẳng khép kín rất bền vững nhờ cặp liên kết hydrogen đôi Csp2-H...Z đối xứng. Thứ hai, khoảng cách liên kết đo được giữa proton H và tâm nhận Z luôn ngắn hơn tổng bán kính van der Waals từ 8% đến 13%, cụ thể khoảng cách O...H dao động từ 2,357 Å đến 2,497 Å so với chuẩn 2,72 Å; S...H đạt 2,667 Å đến 2,858 Å so với chuẩn 3,00 Å; Se...H từ 2,764 Å đến 2,977 Å; và Te...H từ 2,911 Å đến 3,149 Å. Thứ ba, giá trị ái lực proton tại tâm chalcogen tăng vọt theo thứ tự nhóm thế, đạt đỉnh cao nhất ở phân tử NH2CHTe với mức 909,5 kJ/mol, cao hơn tới 41,6% so với FCHO chỉ đạt 642,2 kJ/mol. Thứ tư, độ phân cực của liên kết Csp2-H tăng dần khi chuyển dịch từ O sang Te, khiến enthalpy tách proton giảm mạnh từ mức 1684,9 kJ/mol ở HCHO xuống còn 1540,8 kJ/mol ở BrCHTe.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự bền hóa cấu trúc vòng bắt nguồn từ sự cộng hưởng tương tác tĩnh điện và hiệu ứng chuyển dịch mật độ electron từ đôi electron tự do của tâm Z sang orbital phản liên kết sigma của Csp2-H. Dữ liệu tính toán có thể được trực quan hóa qua biểu đồ phân tán 2 chiều biểu diễn tương quan tuyến tính giữa mật độ electron tại điểm tới hạn liên kết và khoảng cách tương tác, cùng bảng đối chiếu ma trận năng lượng tương tác theo từng nhóm thế. So sánh với các công bố trước đây trên hệ Csp3-H vốn có độ co rút liên kết nhỏ, hệ Csp2-H trong nghiên cứu này ghi nhận mức chuyển dời xanh tần số dao động hóa trị rất rõ nét với độ tăng dao động vượt mốc 80 cm-1 đến 96 cm-1 ở các hệ aldehyde mang nhóm hút electron mạnh. Kết quả này khẳng định vai trò quyết định của nhóm thế trong việc điều biến bản chất liên kết từ chuyển dời đỏ truyền thống sang chuyển dời xanh hiện đại.

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

Để phát triển và ứng dụng tối đa kết quả hóa học lượng tử vào thực tiễn, luận văn đưa ra 4 khuyến nghị then chốt:

  1. Mở rộng khảo sát mô phỏng sang các hệ dị vòng phức tạp chứa liên kết Csp2-H trong vòng 6 đến 12 tháng tới do các nhóm nghiên cứu hóa dược thực hiện, nhằm nâng cao độ chính xác dự đoán liên kết thụ thể lên trên 95%.
  2. Tích hợp bộ thông số nhiệt động của 24 dimer vào các phần mềm sàng lọc ảo trong thời hạn 18 tháng bởi các chuyên gia tin sinh học, giúp giảm thiểu ít nhất 30% chi phí thử nghiệm thực nghiệm in vitro đối với các phối tử chứa chalcogen.
  3. Triển khai phương pháp nhiễu loạn bậc cao và phiếm hàm mật độ điều chỉnh tán xạ cho các hợp chất cơ kim chứa selen và telu trong giai đoạn 2023 đến 2025 do các nhà hóa học tính toán chủ trì, hướng tới mục tiêu hiệu chỉnh sai số năng lượng dưới 1,5 kJ/mol.
  4. Ứng dụng quy luật điều biến nhóm thế hút và đẩy electron để thiết kế các vật liệu mạng tinh thể hữu cơ xốp mới trong chu kỳ 2 năm do các viện nghiên cứu vật liệu tiên tiến tiến hành, phấn đấu tăng độ bền nhiệt của khung cấu trúc lên thêm 15% đến 20%.

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

Luận văn là tài liệu tham khảo giá trị cho 4 nhóm đối tượng trọng tâm:

  1. Các nhà nghiên cứu hóa học tính toán và lượng tử: Tiếp cận quy trình chuẩn hóa kết hợp phương pháp MP2 và CCSD(T) với hơn 50 cấu trúc tối ưu để ứng dụng trực tiếp vào việc tính toán tương tác phi cộng hóa trị chính xác.
  2. Các kỹ sư thiết kế thuốc và hóa dược phẩm: Sử dụng bộ 48 chỉ số về ái lực proton và enthalpy tách proton để mô phỏng và tối ưu hóa lực liên kết giữa phân tử thuốc tiềm năng với các enzyme đích.
  3. Giảng viên và học viên cao học chuyên ngành hóa học: Khai thác tài liệu làm bài giảng chuyên đề về thuyết nguyên tử trong phân tử và cơ chế chuyển dời xanh với 7 bảng số liệu chi tiết.
  4. Các nhà khoa học vật liệu tinh thể siêu phân tử: Ứng dụng dữ liệu tương tác của 4 nguyên tố chalcogen để định hướng tổng hợp mạng lưới polymer phối trí có cấu trúc ổn định và kiểm soát hướng sắp xếp không gian.

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

Liên kết hydrogen chuyển dời xanh khác biệt gì so với liên kết hydrogen cổ điển?

Khác với liên kết cổ điển làm kéo dài liên kết và giảm tần số dao động, liên kết hydrogen chuyển dời xanh làm co rút độ dài liên kết Csp2-H và làm gia tăng tần số dao động hóa trị. Trong thực tế tính toán của luận văn, tần số dao động có thể tăng từ 80 cm-1 đến trên 90 cm-1 khi hình thành cấu trúc dimer bền.

Tại sao nghiên cứu lại chọn khảo sát 4 nguyên tố nhóm chalcogen từ O đến Te?

Các nguyên tố O, S, Se, Te có bán kính nguyên tử và độ âm điện biến thiên tuần hoàn rất rõ rệt. Việc khảo sát đồng thời 4 nguyên tố giúp làm sáng tỏ sự biến đổi của ái lực proton từ 642,2 kJ/mol đến 909,5 kJ/mol, tạo cơ sở so sánh toàn diện giữa nguyên tố nhẹ và nguyên tố nặng.

Thuyết nguyên tử trong phân tử đóng vai trò gì trong việc xác định liên kết?

Thuyết nguyên tử trong phân tử cung cấp các tiêu chuẩn định lượng thông qua điểm tới hạn liên kết. Mật độ electron tại điểm tới hạn trong khoảng 0,002 đến 0,035 au cùng dấu dương của toán tử Laplacian là bằng chứng vững chắc khẳng định sự tồn tại của liên kết hydrogen dạng vòng trong cả 24 dimer.

Nhóm thế hút electron ảnh hưởng thế nào đến độ bền của liên kết Csp2-H...Z?

Các nhóm thế hút electron mạnh như F, Cl, Br làm tăng độ phân cực của liên kết Csp2-H, làm giảm enthalpy tách proton xuống mức khoảng 1505 kJ/mol đến 1583 kJ/mol. Hiện tượng này giúp tăng cường khả năng cho proton và làm cho liên kết hydrogen trong dimer trở nên bền vững hơn.

Phương pháp tính toán nào được sử dụng để loại bỏ sai số chồng chất bộ cơ sở?

Luận văn áp dụng kỹ thuật hiệu chỉnh đối trọng Counterpoise của Boys và Bernardi ở mức lý thuyết chùm tương tác CCSD(T). Kỹ thuật này giúp loại bỏ hoàn toàn hiện tượng năng lượng tương tác bị ước tính quá cao do sự dùng chung hàm cơ sở giữa hai monomer trong không gian tương tác.

Kết luận

Tổng kết các kết quả trọng tâm của luận văn:

  • Xác lập cấu trúc hình học tối ưu và khẳng định tính bền vững của 24 hệ dimer dạng vòng 6 cạnh khép kín.
  • Đo đạc chính xác khoảng cách liên kết H...Z ngắn hơn tổng bán kính van der Waals từ 8% đến 13%.
  • Định lượng toàn diện giá trị ái lực proton trong khoảng 642,2 đến 909,5 kJ/mol và enthalpy tách proton từ 1505,0 đến 1684,9 kJ/mol.
  • Giải mã cơ chế chuyển dời xanh tần số dao động dựa trên thuyết nguyên tử trong phân tử và phân tích orbital liên kết tự nhiên.
  • Hoàn thiện hệ thống dữ liệu nhiệt động học ở mức lý thuyết chùm tương tác CCSD(T) có độ chuẩn xác cao.

Đóng góp then chốt của công trình là đặt nền móng lý thuyết vững chắc cho các nghiên cứu tương tác chalcogen phi cộng hóa trị. Trong giai đoạn 12 tháng tiếp theo, các mô hình này cần được mở rộng sang hệ dung môi và phức chất sinh học. Độc giả và các nhóm nghiên cứu quan tâm hãy khai thác toàn văn công trình để ứng dụng trọn vẹn bộ cơ sở dữ liệu khoa học quý giá này.