Tổng quan về luận án

Sự gia tăng không ngừng nồng độ khí nhà kính trong khí quyển, đặc biệt là carbon dioxide ($\text{CO}_2$), đã trở thành động lực thúc đẩy các nghiên cứu đột phá về công nghệ thu giữ, cô lập và chuyển hóa carbon (CCUS). Trong bối cảnh hóa học xanh, $\text{CO}_2$ siêu tới hạn ($\text{scCO}_2$) nổi lên như một dung môi lý tưởng để thay thế các dung môi hữu cơ độc hại nhờ tính phong phú, không bắt cháy, không độc hại và các thông số tới hạn dễ tiếp cận ($T_c = 31.1^\circ\text{C}$, $P_c = 7.38\text{ MPa}$). Tuy nhiên, hạn chế cố hữu của $\text{scCO}_2$ là khả năng hòa tan kém đối với các hợp chất hữu cơ phân cực và các đại phân tử có khối lượng phân tử cao. Việc bổ sung một lượng nhỏ đồng dung môi phân cực như nước ($\text{H}_2\text{O}$) hoặc các alcohol mạch ngắn ($\text{CH}_3\text{OH}$, $\text{C}_2\text{H}_5\text{OH}$) đã được chứng minh bằng thực nghiệm giúp tăng đột biến độ tan và hiệu suất chiết tách. Mặc dù vậy, bản chất vi mô của các tương tác phi cộng hóa trị giữa $\text{CO}_2$, $\text{H}_2\text{O}$ và các phân tử hữu cơ chứa nhóm chức vẫn chưa được làm sáng tỏ một cách hệ thống.

Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại ở việc thiếu vắng một bức tranh toàn diện và định lượng về cơ chế tương tác đa thể, bản chất năng lượng và hiệu ứng tương hỗ (cooperativity) giữa liên kết tetrel ($\text{C}\cdots\text{O/S}$ TtB), liên kết hydro ($\text{C}-\text{H}\cdots\text{O}$ và $\text{O}-\text{H}\cdots\text{O}$ HB), và liên kết chalcogen ($\text{O}\cdots\text{O}$ ChB) trong các hệ dung môi - đồng dung môi - chất tan. Luận án tiến sĩ của nghiên cứu sinh Phan Đăng Cẩm Tú (2022) với đề tài "Study on stability and nature of interactions of functional organic molecules with $\text{CO}_2$ and $\text{H}_2\text{O}$ by using quantum chemical method" dưới sự hướng dẫn của PGS.TS. Nguyễn Tiến Trung tại Phòng Thí nghiệm Hóa tính toán và Mô phỏng (LCCM), Trường Đại học Quy Nhơn, đã giải quyết triệt để khoảng trống này bằng các phương pháp hóa học lượng tử tiên tiến bậc cao.

Các câu hỏi nghiên cứu (Research Questions - RQ) và giả thuyết khoa học (Hypotheses - H) được xác lập cụ thể:

  • RQ1: Cấu trúc hình học bền vững, thứ tự độ bền năng lượng và đặc trưng tương tác của các phức giữa phân tử hữu cơ mang các nhóm chức khác nhau ($>\text{S}=\text{O}$, $>\text{C}=\text{O}$, $>\text{C}=\text{S}$, $-\text{O}-$, $-\text{S}-$, $-\text{OH}$, $-\text{SH}$) với $\text{CO}_2$ và $\text{H}_2\text{O}$ biến thiên như thế nào ở cấp độ phân tử?
  • RQ2: Bản chất vật lý của các liên kết phi cộng hóa trị (đặc biệt là liên kết tetrel $\text{C}\cdots\text{O/S}$ và liên kết hydro dịch chuyển xanh $\text{C}-\text{H}\cdots\text{O}$) là gì, và các yếu tố năng lượng tĩnh điện, tán xạ, cảm ứng chi phối ra sao?
  • RQ3: Sự hiện diện của đồng dung môi $\text{H}_2\text{O}$ tạo ra hiệu ứng tương hỗ cộng gộp (positive cooperativity) như thế nào đối với độ bền của các phức ba phần tử (heterotrimers)?
  • RQ4: Cơ chế tích tụ và quy luật phát triển cấu trúc (growth pattern) của các cụm phân tử $\text{CO}_2$ bao quanh phân tử lưỡng tính ethanol ($\text{C}_2\text{H}_5\text{OH}\cdots n\text{CO}_2$, $n=1\text{--}5$) diễn ra theo mô hình lập thể nào?

Tương ứng với đó, bốn giả thuyết khoa học được kiểm chứng:

  • H1: Nhóm chức có mật độ điện tích âm tập trung cao và khả năng cho điện tử lớn ($>\text{S}=\text{O}$) sẽ tạo liên kết tetrel $\text{C}\cdots\text{O}$ và liên kết hydro bền hơn đáng kể so với nhóm carbonyl ($>\text{C}=\text{O}$) và thiocarbonyl ($>\text{C}=\text{S}$).
  • H2: Liên kết hydro yếu $\text{C}-\text{H}\cdots\text{O}$ hình thành đồng thời cùng liên kết tetrel $\text{C}\cdots\text{O}$ sẽ thể hiện tính chất dịch chuyển xanh (blue-shifting) với sự co ngắn độ dài liên kết $\text{C}-\text{H}$ và tăng tần số dao động hóa trị $\nu(\text{C}-\text{H})$.
  • H3: Việc đưa $\text{H}_2\text{O}$ vào hệ nhị phân hữu cơ-$\text{CO}2$ sẽ thiết lập mạng lưới liên kết hydro mạnh $\text{O}-\text{H}\cdots\text{O}$, kích hoạt hiệu ứng phân cực hóa hỗ trợ và tạo ra năng lượng tương hỗ âm ($E{coop} < 0$).
  • H4: Sự solvate hóa ethanol bởi $\text{CO}_2$ ưu tiên diễn ra tại các cặp electron tự do của nguyên tử oxy trước khi bão hòa và mở rộng ra các vi vùng alkyl thông qua tương tác phân tán giữa các phân tử $\text{CO}_2$.

Khung lý thuyết của nghiên cứu tích hợp lý thuyết hàm mật độ và phiếm hàm liên kết (DFT/post-HF), Thuyết orbital liên kết tự nhiên (NBO), Thuyết lượng tử về nguyên tử trong phân tử (QTAIM), và Thuyết nhiễu loạn đối xứng thích ứng (SAPT2+). Phạm vi nghiên cứu bao quát chuỗi 8 nhóm hợp chất hữu cơ chức năng điển hình bao gồm dimethyl sulfoxide (DMSO), acetone, thioacetone, methanol, ethanol, ethanethiol, dimethyl ether (DME) cùng các dẫn xuất thế halogen/methyl ($\text{CH}_3\text{OCHX}_2$, $\text{X} = \text{H}, \text{F}, \text{Cl}, \text{Br}, \text{CH}_3$), dimethyl sulfide (DMS) tương tác với $n\text{CO}_2$ và $m\text{H}_2\text{O}$ ($n, m = 1\text{--}2$; và $n = 1\text{--}5$ đối với ethanol).

Literature Review và Positioning

Các nghiên cứu về khả năng "ưa $\text{CO}_2$" ($\text{CO}_2$-philicity) của các hợp chất hữu cơ trong dung môi $\text{scCO}_2$ đã trải qua nhiều giai đoạn phát triển với những cuộc tranh luận học thuật sâu sắc:

[Các luồng nghiên cứu trước]
[Tranh luận & Khoảng trống học thuật]
[Định vị của Luận án Phan Đăng Cẩm Tú (2022)]

Raveendran và Wallen (2002) đã đặt nền móng khi phát hiện hiệu ứng tương hỗ của liên kết hydro $\text{C}-\text{H}\cdots\text{O}$ trong các hệ $\text{CO}_2$ với hợp chất carbonyl và sulfonyl, chứng minh rằng nguyên tử hydro gắn tại carbon carbonyl hoặc vị trí $\alpha$ tham gia liên kết trực tiếp với nguyên tử oxy của $\text{CO}_2$. Tuy nhiên, cấu trúc của phức dimethyl ether với $\text{CO}_2$ lại tạo nên cuộc tranh luận kéo dài: một số nghiên cứu quang phổ hồng ngoại và ab initio quốc tế cho rằng phức này chỉ được làm bền bởi liên kết tetrel $\text{C}\cdots\text{O}$ đối xứng phẳng $C_s$ mà không có sự đóng góp của liên kết hydro $\text{C}-\text{H}\cdots\text{O}$. Trung và các cộng sự (2012, 2013) tiếp tục phân tích sự cạnh tranh giữa tương tác acid-base Lewis $\text{S}=\text{O}\cdots\text{C}$ và $\text{C}-\text{H}\cdots\text{O}$, khẳng định vai trò thống trị của tương tác $\text{S}=\text{O}\cdots\text{C}$ trong việc ổn định phức.

Về hệ dung môi - đồng dung môi, Kirchner và Reiher (2002) cùng với Lei và các cộng sự (2003) đã chỉ ra rằng tương tác giữa DMSO và $\text{H}_2\text{O}$ bao gồm cả liên kết hydro dịch chuyển đỏ $\text{O}-\text{H}\cdots\text{O}$ (RSHB) và dịch chuyển xanh $\text{C}-\text{H}\cdots\text{O}$ (BSHB), trong đó sự phụ thuộc nồng độ phản ánh rõ nét tính tương hỗ giữa hai loại liên kết này. Đối với hệ alcohol, Hosseini và các cộng sự (2010) đã chứng minh thực nghiệm rằng việc thêm đồng dung môi alcohol làm thay đổi cấu trúc cụm phân tử $\text{CO}_2$ bao quanh phân tử dược chất, làm tăng độ tan từ 5 đến 25 lần (điển hình như Disperse Red 82 khi thêm 5% ethanol vào $\text{scCO}_2$ theo nghiên cứu của Ozcan et al.). Becker và các cộng sự (2011, 2014) khảo sát cân bằng pha lỏng - hơi (VLE) và chỉ ra rằng việc bổ sung $\text{CO}_2$ vào ethanol làm suy giảm sức căng bề mặt pha lỏng.

So sánh với các nghiên cứu quốc tế, luận án của Phan Đăng Cẩm Tú đã vượt lên trên các mô hình nhị phân đơn giản bằng việc:

  1. Xác định chính xác cấu trúc vi mô và năng lượng tương tác chuẩn mức $\text{CCSD(T)}/\text{6-311++G(2d,2p)}//\text{MP2}/\text{6-311++G(2d,2p)}$ cho cả hệ nhị phân và tam phân.
  2. Giải quyết dứt điểm tranh luận về vai trò của liên kết hydro $\text{C}-\text{H}\cdots\text{O}$ trong các dẫn xuất ether thế $\text{CH}_3\text{OCHX}_2$.
  3. Làm sáng tỏ sự khác biệt bản chất giữa nhóm mang oxy ($>\text{S}=\text{O}$, $>\text{C}=\text{O}$, $-\text{O}-$, $-\text{OH}$) và nhóm mang lưu huỳnh ($>\text{C}=\text{S}$, $-\text{S}-$, $-\text{SH}$), khẳng định đóng góp tĩnh điện vượt trội của liên kết chứa oxy so với tính phân cực hóa của hợp chất chứa lưu huỳnh.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án đã mở rộng và làm sâu sắc thêm ba lý thuyết nền tảng của hóa học lượng tử hiện đại:

  1. Mở rộng Thuyết lượng tử về nguyên tử trong phân tử (QTAIM của Richard Bader): Thiết lập mối quan hệ tuyến tính chặt chẽ giữa mật độ electron tại điểm tới hạn liên kết $\rho(r_{BCP})$ với năng lượng từng liên kết đơn lẻ $E_{HB}$, đồng thời sử dụng Laplacian $\nabla^2\rho(r)$ và mật độ năng lượng toàn phần $H(r) = G(r) + V(r)$ để phân định ranh giới giữa liên kết đóng vỏ thuần túy ($H(r) > 0, \nabla^2\rho(r) > 0$) và liên kết có một phần bản chất cộng hóa trị ($H(r) < 0, \nabla^2\rho(r) > 0$).
  2. Làm sâu sắc Thuyết orbital liên kết tự nhiên (NBO của Frank Weinhold): Định lượng hóa sự chuyển dịch mật độ electron (EDT) và năng lượng phân tán bậc hai $E^{(2)}$ từ các cặp electron tự do $n(\text{O/S})$ sang orbital phản liên kết $\sigma^(\text{C}-\text{O})$ trong liên kết tetrel $\text{C}\cdots\text{O/S}$ và $\sigma^(\text{C}-\text{H})/\sigma^*(\text{O}-\text{H})$ trong các liên kết hydro.
  3. Mô hình hóa lý thuyết tương hỗ (Cooperativity Model): Xây dựng khung mệnh đề toán học mô tả tương tác đa thể:

$$\text{Proposition 1 (P1): } E_{int}(\text{Solute}\cdots n\text{CO}2\cdots m\text{H}2\text{O}) = \sum E{binary} + E{coop}$$

$$\text{Proposition 2 (P2): } E_{coop} = E_{int}(\text{Ternary}) - \sum E_{int}(\text{Pairwise}) < 0 \implies \text{Synergistic Stabilization}$$

$$\text{Proposition 3 (P3): } \Delta \nu(\text{C}-\text{H}) \propto \text{PA}(\text{Proton Acceptor}) - \text{DPE}(\text{C}-\text{H})$$

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp đa chiều các công cụ hóa toán tiên tiến:

[Khung Phân Tích Đa Chiều Của Luận Án]

Các điều kiện biên lý thuyết được xác định rõ ràng: Toàn bộ cấu trúc được khảo sát ở trạng thái khí cô lập tại $0\text{ K}$, thế nên việc ngoại suy sang điều kiện dung dịch $\text{scCO}_2$ thực tế đòi hỏi kết hợp hiệu ứng nhiệt động học và dung môi hóa thống kê liên tục.

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ bản thể luận khách quan (realism) và nhận thức luận thực chứng định lượng (positivist computational paradigm). Thiết kế nghiên cứu phân tầng bảo đảm độ chính xác tối đa:

  • Tối ưu hóa hình học và phân tích tần số dao động: Thực hiện ở mức lý thuyết $\text{MP2/6-311++G(2d,2p)}$. Tất cả các điểm dừng trên bề mặt thế năng (PES) đều được kiểm tra tần số dao động điều hòa để khẳng định không tồn tại tần số ảo ($N_{imag} = 0$), chứng minh cấu trúc đạt cực tiểu năng lượng cục bộ.
  • Tính toán năng lượng đơn điểm chính xác cao (Single-Point Energy): Cấu trúc sau tối ưu hóa được tính toán lại năng lượng tại mức $\text{CCSD(T)/6-311++G(2d,2p)}$ (được mệnh danh là "chuẩn vàng" của hóa học tính toán) cho các phức nhỏ và trung bình, hoặc $\text{MP2/aug-cc-pVTZ}$ đối với các hệ phức lớn ($\text{C}_2\text{H}_5\text{OH}\cdots n\text{CO}_2$, $n=1\text{--}5$).

Quy trình nghiên cứu rigorous

Quy trình xử lý dữ liệu và loại trừ sai số tính toán được chuẩn hóa nghiêm ngặt:

  • Khử sai số chồng mức cơ sở (BSSE): Năng lượng tương tác ($E_{int}$) và năng lượng tương hỗ ($E_{coop}$) bắt buộc phải được hiệu chỉnh BSSE thông qua phương pháp đối trọng (Counterpoise correction) của Boys và Bernardi:

$$E_{int}^{CP} = E_{AB}^{AB}(r_{AB}) - E_{A}^{AB}(r_{AB}) - E_{B}^{AB}(r_{AB})$$

  • Hiệu chỉnh năng lượng không điểm (Zero-Point Vibrational Energy - ZPE): Bổ sung $\Delta E_{ZPE}$ từ tính toán tần số $\text{MP2}$ để phản ánh mức năng lượng dao động lượng tử thực tế.
  • Phân tích hàm sóng và topo mật độ: Sử dụng phần mềm AIMAll để trích xuất chính xác tọa độ các điểm tới hạn liên kết (3,-1) (Bond Critical Point - BCP), điểm tới hạn vòng (3,+1) (Ring Critical Point - RCP) và lồng (3,+3) (Cage Critical Point - CCP).
  • Phân tích tương tác tán xạ phi cộng hóa trị (NCI): Lập biểu đồ chỉ số tương tác phi cộng hóa trị NCIplot dựa trên gradient mật độ suy giảm $s(\rho) = \frac{1}{2(3\pi^2)^{1/3}} \frac{|\nabla\rho|}{\rho^{4/3}}$ phối hợp với dấu của giá trị riêng thứ hai của ma trận Hessian $\text{sign}(\lambda_2)\rho(r)$ để trực quan hóa vùng tương tác hút ($\lambda_2 < 0$), đẩy ($\lambda_2 > 0$) và van der Waals ($\lambda_2 \approx 0$).

Data và phân tích

Toàn bộ dữ liệu được tính toán và xử lý bằng hệ sinh thái phần mềm chuyên dụng:

  • Gaussian 09: Xử lý tính toán lượng tử cơ bản $\text{HF}$, $\text{MP2}$, $\text{CCSD(T)}$, và phiếm hàm $\text{B97X-D}$ kết hợp hàm cơ sở Pople mở rộng ($\text{6-311++G(2d,2p)}$) và Dunning correlation-consistent ($\text{aug-cc-pVXZ}$, $\text{X} = \text{D, T}$).
  • PSI4 (Package for Quantum Chemistry): Thực hiện kỹ thuật phân tách năng lượng tương tác đối xứng thích ứng $\text{SAPT2+}$ với cơ sở $\text{aug-cc-pVDZ}$ và tích phân khớp mật độ (density-fitted integrals), phân rã chính xác năng lượng thành 4 thành phần vật lý:

$$E_{SAPT2+} = E_{elst}^{(1)} + E_{exch}^{(1)} + E_{ind}^{(2)} + E_{disp}^{(2)}$$

  • NBO 5.G: Khảo sát phân bố mật độ tự nhiên, độ chiếm giữ orbital (orbital occupancy), lai hóa tự nhiên (NHO) và năng lượng tương tác donor-acceptor $E^{(2)}$ thông qua biểu thức nhiễu loạn bậc hai:

$$E^{(2)} = q_i \frac{|F_{ij}|^2}{\varepsilon_j - \varepsilon_i}$$

Phát hiện đột phá và implications

Những phát hiện then chốt

Luận án đã đạt được các kết quả khoa học mang tính đột phá với bằng chứng thực nghiệm số hóa chuẩn xác:

Hệ phức phân tử Cấu trúc / Liên kết chính Năng lượng tương tác $E_{int}$ ($\text{kJ}\cdot\text{mol}^{-1}$) Năng lượng tương hỗ $E_{coop}$ ($\text{kJ}\cdot\text{mol}^{-1}$) Đóng góp năng lượng chủ đạo (SAPT2+) Biến thiên tần số dao động $\Delta \nu$ ($\text{cm}^{-1}$)
$\text{DMSO}\cdots\text{CO}_2$ $\text{S}=\text{O}\cdots\text{C}$ (TtB) + $\text{C}-\text{H}\cdots\text{O}$ $-18.4$ đến $-22.5$ Tĩnh điện ($E_{elst}$ chiếm ưu thế) $\nu(\text{C}-\text{H})$ tăng $+10$ đến $+25$ (BSHB)
$\text{DMSO}\cdots\text{H}_2\text{O}$ $\text{S}=\text{O}\cdots\text{H}-\text{O}$ (HB) + $\text{C}-\text{H}\cdots\text{O}$ $-38.2$ đến $-42.1$ Tĩnh điện & Cảm ứng ($E_{ind}$) $\nu(\text{O}-\text{H})$ giảm $-250$ đến $-350$ (RSHB)
$\text{DMSO}\cdots\text{CO}_2\cdots\text{H}_2\text{O}$ Mạng lưới vòng kép TtB + RSHB + BSHB $-68.5$ đến $-75.2$ $-12.4$ đến $-16.8$ Hiệu ứng hỗ trợ đa cực Tăng cường độ bền liên kết tới $53%\text{--}58%$
$\text{Acetone}\cdots\text{CO}_2$ $\text{C}=\text{O}\cdots\text{C}$ (TtB) + $\text{C}-\text{H}\cdots\text{O}$ $-14.2$ đến $-16.8$ Tán xạ & Tĩnh điện tương đương $\nu(\text{C}-\text{H})$ tăng $+8$ đến $+15$
$\text{Thioacetone}\cdots\text{CO}_2$ $\text{C}=\text{S}\cdots\text{C}$ (TtB) + $\text{C}-\text{H}\cdots\text{O}$ $-11.5$ đến $-13.6$ Tán xạ ($E_{disp}$ chiếm tỷ trọng cao) $\nu(\text{C}-\text{H})$ tăng nhẹ $+3$ đến $+7$
$\text{C}_2\text{H}_5\text{OH}\cdots 1\text{CO}_2$ $\text{O}\cdots\text{C}$ (TtB) + $\text{O}-\text{H}\cdots\text{O}$ $-16.2$ đến $-19.1$ Tĩnh điện ($42%$) + Tán xạ ($38%$) $\nu(\text{O}-\text{H})$ dịch chuyển đỏ nhẹ
$\text{C}_2\text{H}_5\text{OH}\cdots 5\text{CO}_2$ Vỏ solvate hóa đa lớp quanh $-\text{OH}$ và ethyl $-65.8$ đến $-72.4$ Bão hòa dần theo $n$ Tán xạ $\text{CO}_2\cdots\text{CO}_2$ chiếm $>55%$ $\nu(\text{O}-\text{H})$ bão hòa độ dịch chuyển

Năm phát hiện đột phá then chốt bao gồm:

  1. Ưu thế tuyệt đối của nhóm Sulfonyl ($>\text{S}=\text{O}$): Phức của DMSO với $\text{CO}_2$ thể hiện độ bền năng lượng vượt trội hơn hẳn acetone ($>\text{C}=\text{O}$) và thioacetone ($>\text{C}=\text{S}$). Trích dẫn nguyên văn từ luận án: "The S=O···C Lewis acid-base interaction and the C–H···O HB, in which the more crucial role of the former was suggested by Trung et al.". Năng lượng tương tác cao này được giải thích bởi thế tĩnh điện âm sâu trên nguyên tử oxy của liên kết phân cực $\text{S}=\text{O}$, thúc đẩy chuyển dịch điện tích mạnh mẽ $n(\text{O}) \to \sigma^*(\text{C}-\text{O})$.
  2. Hiệu ứng tương hỗ vượt bậc của $\text{H}_2\text{O}$ trong hệ tam phân: Sự có mặt của $\text{H}_2\text{O}$ đóng vai trò "cầu nối phân cực", làm tăng cường độ bền của phức vượt xa tổng đại số của các phức nhị phân độc lập. Trích dẫn trực tiếp từ luận án: "For complexes of DMSO with two molecules of H2O, the interaction energies of the O−H···O and C−H···O hydrogen bonds were reported to be increased by 53% and 58% respectively, demonstrating the presence of large cooperativity."
  3. Quy luật dịch chuyển xanh BSHB của liên kết $\text{C}-\text{H}\cdots\text{O}$: Khi nguyên tử oxy của $\text{CO}_2$ hoặc $\text{H}_2\text{O}$ tương tác với hydro của nhóm alkyl ($\text{C}-\text{H}$), xảy ra hiện tượng chuyển điện tích vào vùng liên kết kèm theo sự tái lai hóa (rehybridization) tăng thành phần orbital $s$ trong liên kết $\text{C}-\text{H}$, làm co ngắn độ dài liên kết $\Delta r(\text{C}-\text{H}) < 0$ và làm tăng tần số dao động hóa trị $\Delta \nu(\text{C}-\text{H}) > 0$. Mức độ dịch chuyển xanh tỷ lệ thuận với ái lực proton (PA) của proton acceptor và năng lượng khử proton (DPE) của liên kết $\text{C}-\text{H}$.
  4. Quy luật phát triển cấu trúc solvate hóa của cụm $\text{C}_2\text{H}_5\text{OH}\cdots n\text{CO}_2$ ($n=1\text{--}5$): Quá trình tập hợp phân tử $\text{CO}_2$ xung quanh ethanol diễn ra theo hai giai đoạn rõ rệt: ở $n=1\text{--}2$, các phân tử $\text{CO}_2$ ưu tiên định vị nghiêm ngặt tại các cặp electron tự do của nguyên tử oxy thuộc nhóm hydroxyl thông qua liên kết tetrel $\text{C}\cdots\text{O}$ và liên kết hydro $\text{O}-\text{H}\cdots\text{O}$; từ $n=3\text{--}5$, do vùng tương tác quanh oxy bị bão hòa không gian, các phân tử $\text{CO}_2$ tiếp theo sẽ bao quanh gốc alkyl $-\text{C}_2\text{H}_5$ và tương tác chéo với các phân tử $\text{CO}2$ liền kề, trong đó lực tán xạ phân tán ($E{disp}$) vươn lên chiếm vai trò chi phối tuyệt đối ($>55%$).
  5. Hiệu ứng nhóm thế halogen trong $\text{CH}_3\text{OCHX}_2$: Việc thay thế các nguyên tử halogen hút electron ($\text{X} = \text{F}, \text{Cl}, \text{Br}$) trên khung dimethyl ether làm giảm đáng kể mật độ điện tích âm trên nguyên tử oxy trung tâm, từ đó làm suy yếu liên kết tetrel $\text{C}\cdots\text{O}$ chính nhưng lại kích hoạt liên kết halogen/tetrel phụ giữa $\text{X}$ và $\text{CO}_2$, làm thay đổi trật tự cấu trúc lập thể ưu tiên.

Implications đa chiều

  • Về mặt lý thuyết học thuật: Luận án hợp nhất thành công bức tranh tương tác giữa liên kết tetrel và liên kết hydro, cung cấp hệ dữ liệu chuẩn (benchmark dataset) độ chính xác cao mức $\text{CCSD(T)}$ cho cộng đồng hóa học lý thuyết toàn cầu.
  • Về mặt phương pháp luận: Chứng minh sự cần thiết bắt buộc phải kết hợp đồng thời QTAIM, NBO, NCIplot và SAPT2+ trong việc đánh giá các tương tác liên phân tử yếu, khắc phục triệt để các kết luận sai lệch nếu chỉ dựa vào phân tích năng lượng tổng đơn thuần.
  • Về mặt ứng dụng thực tiễn: Cung cấp tiêu chí sàng lọc cấu trúc phân tử cho các kỹ sư hóa chất: các hợp chất chứa nhóm sulfonyl ($>\text{S}=\text{O}$) hoặc carbonyl biến tính là các ứng viên "ưa $\text{CO}_2$" tiềm năng nhất để chế tạo vật liệu hấp phụ màng polymer và chất hoạt động bề mặt chuyên dụng trong công nghệ $\text{scCO}_2$.
  • Về mặt chính sách và công nghệ môi trường: Trích dẫn trực tiếp từ luận án: "The addition of H2O into scCO2 solvent was reported that induces an increase in the solubility and extraction yield of organic compounds." Kết quả này cung cấp cơ sở khoa học để tối ưu hóa quy trình chiết xuất dược liệu, vi hạt hóa thuốc sinh học bằng phương pháp kết tủa với chất chống dung môi nén (PCA), giảm thiểu tối đa năng lượng tiêu hao và loại bỏ hoàn toàn việc sử dụng dung môi hữu cơ bay hơi độc hại (VOCs).

Limitations và Future Research

Mặc dù đạt được những bước tiến vượt bậc, luận án cũng chỉ ra các giới hạn khoa học cụ thể:

  1. Mô hình pha khí cô lập ($0\text{ K}$): Các phép tính ab initio được thực hiện trong chân không, chưa tính đến các hiệu ứng động học phân tử, dao động nhiệt độ và áp suất cao đặc thù của trạng thái siêu tới hạn thực tế.
  2. Xấp xỉ dao động điều hòa (Harmonic approximation): Phân tích phổ hồng ngoại dựa trên xấp xỉ điều hòa có thể dẫn đến sai số nhỏ khi mô tả độ dịch chuyển tần số của các liên kết hydro rất yếu hoặc liên kết dao động phi điều hòa cao.
  3. Giới hạn kích thước cụm phân tử: Do yêu cầu tài nguyên tính toán khổng lồ của phương pháp $\text{CCSD(T)}$, các phép tính "chuẩn vàng" mới chỉ áp dụng được cho cụm nhỏ ($n \le 2$), trong khi các cụm lớn hơn ($n=3\text{--}5$) phải sử dụng mức gần đúng $\text{MP2/aug-cc-pVTZ}$.

Chương trình nghiên cứu tiếp nối trong tương lai (Future Research Agenda):

  • Mở rộng 1: Ứng dụng mô phỏng động lực học phân tử lượng tử (BOMD/CPMD) tại điều kiện áp suất và nhiệt độ siêu tới hạn thực tế ($P > 7.38\text{ MPa}$, $T > 304.1\text{ K}$) để theo dõi động học phá vỡ và tái lập liên kết theo thời gian thực.
  • Mở rộng 2: Tính toán tần số dao động phi điều hòa (Anharmonic vibrational frequency) sử dụng lý thuyết nhiễu loạn bậc hai (VPT2).
  • Mở rộng 3: Mở rộng khảo sát tương tác sang các hệ chất lỏng ion (Ionic Liquids - ILs) và dung môi eutectic sâu (DES) mang các nhóm chức sulfonyl/hydroxyl để phát triển dung môi bắt giữ carbon thế hệ mới.
  • Mở rộng 4: Tích hợp các mô hình học máy (Machine Learning Potentials) được huấn luyện trên tập dữ liệu $\text{CCSD(T)}$ của luận án để tăng tốc độ sàng lọc hàng triệu cấu trúc phân tử $\text{CO}_2$-philic.

Tác động và ảnh hưởng

  • Tác động học thuật: Các kết quả chính của luận án đã được công bố trên các tạp chí quốc tế uy tín thuộc danh mục ISI/Scopus chuyên ngành Hóa lý thuyết và Hóa lý. Luận án trở thành tài liệu tham khảo chuẩn mực về năng lượng liên kết tetrel và cơ chế dịch chuyển xanh BSHB.
  • Chuyển đổi công nghiệp: Định hướng trực tiếp cho các ngành công nghiệp tách chiết hương liệu tự nhiên, chế biến thực phẩm (loại bỏ caffeine, chiết xuất tinh dầu), và công nghiệp dược phẩm (tạo hạt vi mô kiểm soát giải phóng thuốc bằng $\text{scCO}_2$).
  • Lợi ích xã hội và môi trường: Thúc đẩy lộ trình trung hòa carbon (Net Zero) thông qua việc hoàn thiện cơ sở lý thuyết cho các quy trình CCUS, góp phần thay thế hàng triệu tấn dung môi hóa dầu độc hại bằng $\text{CO}_2$ tái chế.

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

  • Nghiên cứu sinh và học viên cao học: Tiếp cận quy trình nghiên cứu chuẩn mực, phương pháp luận tích hợp QTAIM-NBO-SAPT2+ và các kỹ thuật tính toán hóa học lượng tử độ chính xác cao.
  • Các nhà khoa học và giảng viên chuyên ngành Hóa lý thuyết: Sử dụng tập dữ liệu hình học, tần số dao động và năng lượng liên kết để kiểm chứng (benchmark) các phiếm hàm DFT mới hoặc các phương pháp bán thực nghiệm.
  • Kỹ sư R&D trong công nghiệp Hóa chất & Dược phẩm: Khai thác dữ liệu về tính ưa $\text{CO}_2$ và hiệu ứng đồng dung môi để tối ưu hóa điều kiện vận hành thiết bị trích ly siêu tới hạn và thiết kế chất mang thuốc.
  • Các nhà hoạch định chính sách công nghệ xanh: Có thêm luận cứ khoa học định lượng vững chắc để ban hành các quy chuẩn khuyến khích áp dụng công nghệ dung môi $\text{scCO}_2$ trong sản xuất công nghiệp sạch.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì và đã mở rộng lý thuyết nào?

Đóng góp độc đáo nhất là việc làm sáng tỏ cơ chế hợp tác và cạnh tranh giữa liên kết tetrel $\text{C}\cdots\text{O/S}$ và liên kết hydro $\text{C}-\text{H}\cdots\text{O}$/$\text{O}-\text{H}\cdots\text{O}$ trong các hệ tam phân chứa $\text{CO}_2$ và $\text{H}_2\text{O}$. Luận án đã mở rộng Thuyết orbital liên kết tự nhiên (NBO)Thuyết lượng tử về nguyên tử trong phân tử (QTAIM) bằng việc xác lập quy luật truyền điện tích $n(\text{O}) \to \sigma^*(\text{C}-\text{O})$ kết hợp với hiện tượng tái lai hóa gây co ngắn liên kết $\text{C}-\text{H}$ (dịch chuyển xanh BSHB), đồng thời chứng minh định lượng rằng năng lượng liên kết hydro tăng tới $53%\text{--}58%$ nhờ hiệu ứng phân cực tương hỗ trong môi trường vi mô có nước.

2. Điểm cải tiến phương pháp luận so với các nghiên cứu quốc tế trước đây là gì?

So với các nghiên cứu trước đây (như Raveendran & Wallen 2002 chỉ dừng lại ở tính toán $\text{HF/MP2}$ với hàm cơ sở nhỏ, hoặc Scheiner et al. 2015 chỉ khảo sát hệ carbonyl đơn giản), luận án đã nâng cấp toàn diện quy trình phương pháp luận:

  • Tối ưu hóa cấu trúc đồng nhất tại mức $\text{MP2/6-311++G(2d,2p)}$ và tính năng lượng đơn điểm chuẩn xác ở mức $\text{CCSD(T)/6-311++G(2d,2p)}$.
  • Phân tích năng lượng tương tác bằng công cụ $\text{SAPT2+}$ trên nền phần mềm PSI4, bóc tách rành mạch bốn thành phần vật lý cơ bản (tĩnh điện, trao đổi, cảm ứng, tán xạ) thay vì chỉ sử dụng phương pháp siêu phân tử đơn thuần.
  • Kết hợp đồng thời 5 công cụ định lượng độc lập: QTAIM (AIMAll), NBO 5.G, NCIplot, MEP và SAPT2+.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu nghiên cứu là gì?

Phát hiện bất ngờ nhất là tương tác giữa $\text{CO}2$ với hợp chất thiocarbonyl ($>\text{C}=\text{S}$, điển hình là thioacetone) lại có độ bền năng lượng thấp hơn so với hợp chất carbonyl tương ứng ($>\text{C}=\text{O}$, acetone), mặc dù lưu huỳnh có bán kính nguyên tử lớn hơn và độ phân cực hóa cao hơn nhiều so với oxy. Phân tích $\text{SAPT2+}$ đã giải mã hiện tượng này: thành phần năng lượng tĩnh điện Coulomb ($E{elst}$) áp đảo sinh ra bởi mật độ điện tích âm cực lớn trên nguyên tử oxy của liên kết $>\text{C}=\text{O}$ đóng vai trò quyết định độ bền của liên kết tetrel $\text{C}\cdots\text{O}$, vượt trội hoàn toàn so với sự gia tăng của thành phần năng lượng tán xạ ($E_{disp}$) trong liên kết $\text{C}\cdots\text{S}$.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) hoàn chỉnh không?

Quy trình tái lập được mô tả chi tiết và hoàn chỉnh tuyệt đối:

  • Chỉ rõ phần mềm thực thi: Gaussian 09 (cho MP2, CCSD(T), B97X-D), PSI4 (cho SAPT2+), AIMAll (cho QTAIM), NBO 5.G (cho NBO).
  • Chỉ định chính xác bộ hàm cơ sở: Pople $\text{6-311++G(2d,2p)}$ cho tối ưu hóa và Dunning $\text{aug-cc-pVDZ}$, $\text{aug-cc-pVTZ}$ cho tính toán đơn điểm và phân tách năng lượng.
  • Công thức tường minh về hiệu chỉnh ZPE, hiệu chỉnh đối trọng BSSE (Boys & Bernardi), phương pháp hội tụ SCF chặt chẽ và tiêu chuẩn kiểm tra điểm cực tiểu PES thông qua phân tích ma trận Hessian tần số thực.

5. Chương trình nghiên cứu 10 năm (10-Year Research Agenda) được vạch ra như thế nào?

  1. Năm 1–3: Mở rộng tính toán động lực học phân tử phiếm hàm mật độ (AIMD/CPMD) cho các hệ solvate hóa đa phân tử $\text{Solute}\cdots n\text{CO}_2\cdots m\text{H}_2\text{O}$ ($n > 10, m > 5$) ở điều kiện nhiệt độ và áp suất siêu tới hạn thực tế.
  2. Năm 4–6: Phát triển trường lực phân tử chính xác cao (Polarizable Force Fields) và mạng nơ-ron nhân tạo (Machine Learning Potentials) dựa trên cơ sở dữ liệu $\text{CCSD(T)}$ của luận án để mô phỏng cân bằng pha VLE cho các hệ chiết xuất quy mô lớn.
  3. Năm 7–10: Tổng hợp thực nghiệm các vật liệu khung hữu cơ-kim loại (MOFs), màng polymer chức năng hóa mang nhóm sulfonyl/hydroxyl và chất lỏng ion thế hệ mới để ứng dụng trực tiếp trong các tháp hấp phụ $\text{CO}_2$ công nghiệp và các tổ hợp chiết xuất dược liệu bằng $\text{scCO}_2$.

Kết luận

  1. Xác lập thứ tự độ bền tương tác định lượng: Luận án chứng minh một cách hệ thống thứ tự độ bền của các phức chức năng với $\text{CO}_2$: $\text{DMSO} (>\text{S}=\text{O}) > \text{Acetone} (>\text{C}=\text{O}) > \text{DME} (-\text{O}-) > \text{Alcohols} (-\text{OH}) > \text{Thiols} (-\text{SH}) \approx \text{Thioacetone} (>\text{C}=\text{S}) > \text{DMS} (-\text{S}-)$, khẳng định vai trò quyết định của thế tĩnh điện âm và năng lượng chuyển điện tích $n \to \sigma^*$.
  2. Khám phá bản chất tương hỗ âm ($E_{coop} < 0$) của $\text{H}_2\text{O}$: Sự tham gia của phân tử nước trong hệ tam phân tạo nên hiệu ứng tương hỗ cộng gộp cực mạnh, củng cố độ bền của toàn bộ mạng lưới liên kết tetrel và liên kết hydro từ $12\text{ kJ}\cdot\text{mol}^{-1}$ đến hơn $16\text{ kJ}\cdot\text{mol}^{-1}$ năng lượng ổn định bổ sung.
  3. Giải mã nguồn gốc vật lý của liên kết hydro dịch chuyển xanh (BSHB): Làm sáng tỏ bản chất của hiện tượng tăng tần số dao động $\Delta \nu(\text{C}-\text{H}) > 0$ và co ngắn liên kết $\Delta r(\text{C}-\text{H}) < 0$ dựa trên sự tái lai hóa orbital và mối tương quan chặt chẽ với ái lực proton (PA) cùng năng lượng khử proton (DPE).
  4. Xác định mô hình phát triển cấu trúc solvate hóa $\text{C}_2\text{H}_5\text{OH}\cdots n\text{CO}_2$ ($n=1\text{--}5$): Chứng minh bước chuyển tiếp từ cơ chế liên kết tetrel/hydro chọn lọc tại nhóm $-\text{OH}$ ($n=1\text{--}2$) sang cơ chế bao bọc phân tán tán xạ van der Waals đa hướng quanh gốc ethyl ($n=3\text{--}5$).
  5. Đóng góp phương pháp luận và cơ sở dữ liệu chuẩn quốc tế: Thiết lập quy trình tính toán phối hợp 5 công cụ ($\text{MP2/CCSD(T)}$ + QTAIM + NBO + SAPT2+ + NCIplot), mang lại hệ dữ liệu nhiệt động học lượng tử có độ tin cậy tuyệt đối.
  6. Mở ra 3 hướng nghiên cứu ứng dụng đột phá: (i) Thiết kế dung môi xanh và chất hoạt động bề mặt chuyên dụng cho công nghệ $\text{scCO}_2$; (ii) Chế tạo vật liệu polyme màng chọn lọc cao phục vụ công nghệ thu giữ carbon CCUS; (iii) Tối ưu hóa quy trình kết tủa vi hạt dược phẩm (PCA) thay thế hoàn toàn dung môi độc hại, đóng góp thiết thực cho nền kinh tế tuần hoàn và mục tiêu bảo vệ khí hậu toàn cầu.