Tổng quan về luận án

Sự phá hủy vật liệu do ăn mòn kim loại và thoái hóa sinh hóa do các gốc oxy hoạt động (Reactive Oxygen Species - ROS) là hai thách thức liên ngành lớn trong khoa học vật liệu và y dược học hiện đại. Theo ước tính từ Tổ chức Ăn mòn Thế giới (World Corrosion Organization), thiệt hại kinh tế do ăn mòn kim loại hàng năm chiếm khoảng 3,4% GDP toàn cầu (tương đương 2,5 nghìn tỷ USD), trong khi stress oxy hóa là căn nguyên phân tử gây tổn thương tế bào, ung thư và các bệnh lý thoái hóa mãn tính. Luận án tiến sĩ hóa học của tác giả Đinh Quý Hương với đề tài "Nghiên cứu các chất chống oxy hóa, ức chế ăn mòn kim loại bằng tính toán hóa lượng tử kết hợp với thực nghiệm" (chuyên ngành Hóa vô cơ, mã số 9440113, Đại học Sư phạm – Đại học Huế, 2020) là công trình tiên phong thiết lập mô hình tích hợp đa quy mô giữa hóa học tính toán lượng tử (Quantum Chemical Calculation) và thực nghiệm hóa lý - điện hóa học để giải quyết đồng thời hai bài toán trên.

Khoảng trống học thuật (research gap) trọng tâm mà luận án xác định là: Hầu hết các công trình quốc tế trước đây chỉ tiếp cận đơn lẻ một trong hai hướng: hoặc nghiên cứu chất ức chế ăn mòn kim loại (thường tập trung vào cấu trúc chứa dị tố N, O, S, P và hệ điện tử pi liên hợp nhằm tạo liên kết phối trí với vân đạo d trống của kim loại), hoặc sàng lọc chất chống oxy hóa dập tắt gốc tự do (dựa trên các liên kết có năng lượng phân ly thấp như O-H, N-H, S-H). Chưa có một hệ khung lý thuyết và thực nghiệm tích hợp nào giải quyết bài toán thiết kế phân tử lưỡng chức năng (bifunctional compounds) - vừa sở hữu hoạt tính ức chế ăn mòn vượt trội trên bề mặt hợp kim thép, vừa có khả năng bắt giữ các gốc tự do ROS hiệu năng cao.

Luận án xây dựng hệ thống 4 câu hỏi nghiên cứu và 4 giả thuyết khoa học đối ứng:

  • RQ1: Cấu trúc phân tử của các dẫn xuất thiourea ảnh hưởng như thế nào đến cơ chế hấp phụ và hiệu suất bảo vệ thép carbon trong dung dịch HCl 1,0 M và NaCl 3,5%?
    • Hypothesis 1 (H1): Sự hiện diện đồng thời của tâm cho điện tử (S, N) và vòng thơm phenyl sẽ tăng cường tương tác phối trí - tĩnh điện với bề mặt Fe(110), giúp 1-phenyl-2-thiourea (PTU) đạt hiệu suất ức chế cao hơn cấu trúc mạch hở 1,3-diisopropyl-2-thiourea (ITU) và chất chuẩn urotropine.
  • RQ2: Cơ chế nhiệt động học và động học nào chi phối khả năng dập tắt gốc tự do tự nhiên (DPPH•, ABTS•+) của các dẫn xuất thiourea trong dung môi phân cực?
    • Hypothesis 2 (H2): Khả năng chống oxy hóa của các dẫn xuất thiourea ưu tiên diễn ra theo cơ chế chuyển nguyên tử hydro (HAT) hoặc chuyển proton - mất electron tuần tự (SPLET) tùy thuộc vào độ phân cực dung môi và năng lượng phân ly liên kết (BDE).
  • RQ3: Việc thay thế dị tố lưu huỳnh (S) bằng selen (Se) trong khung phân tử thiourea tạo ra bước chuyển biến nhiệt động học như thế nào đối với phản ứng bắt gốc peoxyl (HOO•)?
    • Hypothesis 3 (H3): Do selen có độ phân cực hóa cao và bán kính nguyên tử lớn hơn lưu huỳnh, các dẫn xuất selenourea (PSeU) sẽ làm giảm đáng kể thế ion hóa (IE) và BDE(N-H), từ đó nâng cao vượt bậc hằng số tốc độ dập tắt gốc tự do.
  • RQ4: Liệu việc đưa các nhóm thế đẩy/hút điện tử vào vị trí para của vòng phenyl selenourea có cho phép tối ưu hóa đồng thời tính năng ức chế ăn mòn và hoạt tính chống oxy hóa?
    • Hypothesis 4 (H4): Nhóm thế cho điện tử mạnh (-OCH3) sẽ hạ thấp hàng rào năng lượng hoạt hóa Gibbs ($\Delta G^\neq$) đối với tương tác gốc tự do, đồng thời tăng mật độ điện tử orbital phân tử biên HOMO, tối ưu hóa năng lượng hấp phụ lên bề mặt Fe(110).

Phạm vi nghiên cứu thực nghiệm được triển khai đồng bộ trên mẫu thép carbon trong môi trường xâm thực HCl 1,0 M và NaCl 3,5% ở dải nhiệt độ 20 °C – 60 °C (293 K – 333 K) với nồng độ chất ức chế từ $10^{-4}$ M đến $5 \cdot 10^{-3}$ M; khảo sát hoạt tính chống oxy hóa in vitro trên gốc tự do chuẩn DPPH• và ABTS•+; kết hợp tính toán phiếm hàm mật độ (DFT) ở các mức lý thuyết cao cấp ROB3LYP/6-311++G(2df,2p)//B3LYP/6-311G(d,p) và M05-2X/6-311++G(d,p), mô phỏng hấp phụ Monte Carlo và Động lực học phân tử (MD) trên bề mặt mạng tinh thể Fe(110).

Literature Review và Positioning

Nghiên cứu về cơ chế ức chế ăn mòn kim loại từ lâu đã xác lập nguyên lý cơ bản: các phân tử hữu cơ chứa dị tố mang cặp electron tự do (N, O, S, P) có xu hướng chuyển điện tử vào vân đạo d chưa lấp đầy của ion kim loại chuyển tiếp để hình thành liên kết phối trí, đồng thời tiếp nhận mật độ điện tử phản hồi từ bề mặt kim loại vào vân đạo phản liên kết $\pi^*$ (Sanyal, 1981). Trong trường phái này, Li và cộng sự (2012) khi nghiên cứu allyl thiourea trong dung dịch $\text{H}_3\text{PO}_4$ đã chứng minh hiệu suất bảo vệ đạt trên 95% ở nồng độ 0,5 mM; Hosseini và cộng sự (2009) khẳng định dẫn xuất 1-methyl-3-pyridin-2-yl-thiourea thể hiện đặc tính ức chế hỗn hợp trong $\text{H}_2\text{SO}_4$ với hiệu suất 96,4% ở 100 ppm. Tuy nhiên, các công trình này chủ yếu dừng lại ở việc mô tả thông số điện hóa vĩ mô hoặc hồi quy cấu trúc đơn giản, thiếu sự bóc tách định lượng cơ chế tương tác ở cấp độ obitan phân tử và bản chất năng lượng bề mặt phân chia pha.

Trong lĩnh vực hóa học chống oxy hóa, các công trình nền tảng của Litwinienko và Ingold (2003, 2007) đã phân định rõ ranh giới cạnh tranh giữa các cơ chế dập tắt gốc tự do: Chuyển nguyên tử hydro trực tiếp (Hydrogen Atom Transfer - HAT), Chuyển electron đơn kèm chuyển proton (Single Electron Transfer followed by Proton Transfer - SET-PT), và Mất proton tuần tự kèm chuyển electron (Sequential Proton Loss Electron Transfer - SPLET). Galano và các cộng sự (2011) khi nghiên cứu động học bắt gốc tự do của Trolox đã khẳng định vai trò quyết định của phiếm hàm Minnesota M05-2X trong việc dự đoán chính xác hàng rào thế năng và hằng số tốc độ phản ứng. Về mặt hợp chất chứa nguyên tố họ chalcogen, Klotz (2003) và Fernández-Bolaños cùng cộng sự (2012) đã phát hiện các dẫn xuất selenoureas và selenocarbamates thể hiện hoạt tính khử ROS và bắt chước enzyme glutathione peroxidase (GPx) vượt trội hơn hẳn các hợp chất tương đồng chứa lưu huỳnh nhờ độ âm điện thấp và mật độ điện tử linh động của Se.

                      TIẾP CẬN TRUYỀN THỐNG
               ĐÓNG GÓP TIÊN PHONG CỦA LUẬN ÁN

Vị thế học thuật của luận án được xác lập tại điểm giao thoa đột phá: Luận án không chỉ định vị chính xác vị trí của các dẫn xuất thiourea/selenourea trong chuỗi biến thiên hoạt tính chalcogen ($O < S < Se$), mà còn xây dựng thành công cầu nối phương pháp luận giữa nhiệt động học hấp phụ bề mặt kim loại với động học phản ứng bắt gốc tự do pha dung dịch, điều chưa từng được hệ thống hóa trong bất kỳ công bố quốc tế nào trước đó.

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

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

Luận án đóng góp quan trọng vào việc mở rộng và kiểm chứng ba hệ lý thuyết nền tảng của Hóa học Vật lý:

  1. Lý thuyết Hóa học Lượng tử về Cấu trúc Phân tử và Phản ứng (DFT): Luận án mở rộng khung phân tích phiếm hàm mật độ Hohenberg-Kohn-Sham thông qua việc xác định chính xác các chỉ số phản ứng hóa học toàn cục theo định lý Koopmans và lý thuyết Axit-Bazơ Cứng-Mềm (HSAB) của Pearson: thế hóa học ($\mu$), độ âm điện ($\chi$), độ cứng toàn phần ($\eta$), độ mềm toàn phần ($S$), và số lượng electron trao đổi ($\Delta N$): $$\Delta N = \frac{\chi_{\text{Fe}} - \chi_{\text{inh}}}{2(\eta_{\text{Fe}} + \eta_{\text{inh}})}$$
  2. Lý thuyết Trạng thái Chuyển tiếp (Transition State Theory - TST): Ứng dụng phương trình Eyring kết hợp hiệu chỉnh đường ngầm Wigner và Eckart $\kappa(T)$ để lượng hóa tốc độ phản ứng dập tắt gốc tự do peoxyl ($\text{HOO}^\bullet$) qua tọa độ phản ứng nội (Intrinsic Reaction Coordinate - IRC).
  3. Lý thuyết Hấp phụ Nhiệt động học Giao diện: Chứng minh sự chuyển dịch từ hấp phụ vật lý (phù hợp mô hình đẳng nhiệt Temkin) sang hấp phụ hóa học chiếm ưu thế (mô hình đẳng nhiệt Langmuir) khi nhiệt độ dung dịch axit tăng từ 293 K lên 333 K.
                      MÔ HÌNH HỆ THỐNG ĐÓNG GÓP LÝ THUYẾT
                               CH3OPSeU / PTU
 [CHỐNG ĂN MÒN KIM LOẠI]                                 [CHỐNG OXY HÓA GỐC ROS]
 - Hấp phụ Fe(110) (MC & MD Simulation)                 - Dập tắt HOO• qua cơ chế HAT & RAF
 - Năng lượng tương tác E_inter = -163.57 kcal/mol       - BDE(N12-H13) = 80.6 kcal/mol
 - Hiệu suất ức chế PDP: H = 98.96%                      - Hằng số tốc độ k_tot = 1.34 x 10^7 M^-1.s^-1
 - Màng bảo vệ tự hàn gắn, ngăn chuyển điện tích         - Phân tích Topo QTAIM xác nhận liên kết H

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

Khung phân tích đa tầng của luận án tích hợp bốn công cụ tính toán và thực nghiệm hiện đại:

  • Tích hợp bộ hàm cơ sở mở rộng: Geometry optimization được thực hiện tại mức B3LYP/6-311G(d,p); năng lượng điểm đơn (Single Point Energy) và các thông số nhiệt động (BDE, IE, PA, ETE) được chuẩn hóa tại mức ROB3LYP/6-311++G(2df,2p); động học và bề mặt thế năng (PES) được khảo sát ở mức M05-2X/6-311++G(d,p).
  • Mô hình dung môi hóa tự hợp (SCRF/IEFPCM): Lượng hóa chính xác hiệu ứng phân cực của ethanol và dung dịch nước lên cấu trúc ion trung gian và trạng thái chuyển tiếp (TS).
  • Phân tích Topo Mật độ Electron theo Lý thuyết Nguyên tử trong Phân tử (QTAIM/AIM): Khảo sát mật độ electron $\rho(r)$, giá trị Laplacian $\nabla^2\rho(r)$, mật độ thế năng $V(r)$, và mật độ động năng $G(r)$ tại các điểm tới hạn liên kết (Bond Critical Points - BCP) để xác nhận bản chất tương tác liên phân tử trong phức chất trạng thái chuyển tiếp.
  • Mô phỏng Hấp phụ Bề mặt Mạng Tinh thể Fe(110): Xác định cấu hình không gian tối ưu và năng lượng liên kết $E_{\text{binding}} = -E_{\text{interaction}}$ của phân tử ức chế ở trạng thái trung hòa và proton hóa trong môi trường chất lỏng chứa ion $\text{H}_3\text{O}^+$ và $\text{Cl}^-$.

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

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng nghiêm ngặt (Positivism Paradigm), kết hợp đối chứng chéo (Triangulation) giữa tính toán hóa lượng tử vi mô (Quantum Mechanics - QM), mô phỏng động học bề mặt (Molecular Dynamics - MD, Monte Carlo), và các kỹ thuật phân tích điện hóa học/quang phổ thực nghiệm vĩ mô.

                           KHUNG THIẾT KẾ NGHIÊN CỨU

Quy trình nghiên cứu rigorous

  1. Giao thức Thực nghiệm Điện hóa:
    • Thiết lập bình điện hóa 3 điện cực tiêu chuẩn: Điện cực làm việc (Working Electrode - WE) là thép carbon (diện tích làm việc chính xác 1,0 $\text{cm}^2$, bọc trong nhựa epoxy cách điện, mài nhẵn bằng giấy nhám SiC từ 400 đến 1200 mesh, tẩy dầu bằng ethanol và rửa nước cất siêu sạch); Điện cực đối (Counter Electrode - CE) là tấm lưới bạch kim (Pt) diện tích lớn; Điện cực so sánh (Reference Electrode - RE) là điện cực bạc/bạc clorua (Ag/AgCl/bão hòa KCl).
    • Đo phân cực thế động (PDP): Quét thế tĩnh trong khoảng $\pm 250$ mV so với thế ăn mòn hở $E_{\text{corr}}$ với tốc độ quét 1 mV/s, áp dụng công thức Stern-Geary để tính mật độ dòng ăn mòn $i_{\text{corr}}$: $$i_{\text{corr}} = \frac{\beta_a \cdot \beta_c}{2{,}303(\beta_a + \beta_c) R_p}$$
    • Hiệu suất ức chế ăn mòn ($H%$) được xác định chính xác theo công thức: $$H% = \frac{i_{\text{corr}} - i_{\text{inh}}}{i_{\text{corr}}} \times 100%$$
  2. Giao thức Đo Tổng trở Điện hóa (EIS): Áp tín hiệu điều hòa xoay chiều biên độ nhỏ thế 10 mV qua dải tần số từ 100 kHz xuống 10 mHz; mô phỏng giản đồ Nyquist thông qua mạch tương đương $R_s(R_{\text{ct}} \parallel \text{CPE})$ có tính đến phần tử hằng số pha CPE để bù trừ độ nhám vi mô của bề mặt điện cực.
  3. Giao thức Đánh giá Hoạt tính Chống Oxy hóa In Vitro:
    • Chuẩn bị dung dịch gốc tự do DPPH• ($6 \cdot 10^{-5}$ M trong ethanol tuyệt đối, đo độ hấp thụ quang ở $\lambda_{\max} = 517$ nm).
    • Chuẩn bị cation gốc tự do $\text{ABTS}^{\bullet+}$ tạo ra từ phản ứng giữa ABTS 7 mM và potassium persulfate 2,45 mM, ủ trong bóng tối 16 giờ, đo ở bước sóng $\lambda_{\max} = 734$ nm.
    • Xác định giá trị ức chế 50% gốc tự do ($\text{IC}_{50}$) qua phương trình hồi quy nồng độ - hiệu suất dập tắt.
  4. Giao thức Tính toán Hóa Lượng tử:
    • Sử dụng phần mềm Gaussian 09W và GaussView 5. Tối ưu hóa hình học không gian (Geometry Optimization) không áp đặt bất kỳ ràng buộc đối xứng nào.
    • Tính toán tọa độ phản ứng nội (IRC) tại các trạng thái chuyển tiếp (TS) với chỉ một tần số dao động ảo duy nhất (imaginary frequency) để xác nhận trạng thái chuyển tiếp nối chính xác chất phản ứng và sản phẩm.
    • Mô phỏng Monte Carlo và Động lực học Phân tử (MD) trong Materials Studio 8 với hộp mô phỏng tuần hoàn chứa bề mặt mạng tinh thể $\text{Fe}(110)$, chất ức chế và 500 phân tử nước cùng các ion solvat hóa trong trường lực COMPASS.

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

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

Thứ nhất, hiệu suất ức chế ăn mòn vượt trội của PTU trong môi trường acid đậm đặc và tính chất phụ thuộc nhiệt độ độc đáo: Nghiên cứu chỉ ra rằng 1-phenyl-2-thiourea (PTU) ức chế ăn mòn thép vượt bậc trong dung dịch HCl 1,0 M. Tại nồng độ $5 \cdot 10^{-3}$ M, hiệu suất ức chế tăng dần theo nhiệt độ:

"1-phenyl-2-thiourea có khả năng ức chế ăn mòn thép tốt hơn 1,3-diisopropyl-2-thiourea trong môi trường HCl 1,0 M với hiệu suất đạt 92,00 % ở 20 °C, 94,05 % ở 30 °C, 96,95 % ở 45 °C và 98,96 % ở 60 °C ở nồng độ $5 \cdot 10^{-3}$ M."

Hiện tượng hiệu suất bảo vệ tăng thuận chiều theo nhiệt độ là một bằng chứng thực nghiệm then chốt khẳng định quá trình hấp phụ hóa học (chemisorption) chiếm vai trò chủ đạo ở nhiệt độ cao, hình thành liên kết phối trí bền vững giữa vân đạo $\pi$ của vòng thơm, cặp electron tự do của nguyên tử lưu huỳnh với orbital d trống của nguyên tử sắt.

Chất ức chế Nồng độ (M) Nhiệt độ (°C) Mật độ dòng $i_{\text{corr}}$ ($\mu\text{A}/\text{cm}^2$) Hiệu suất $H%$ (PDP) Hiệu suất $H%$ (EIS)
Mẫu trắng (HCl 1,0 M) 0 30 1120,5 - -
ITU $5 \cdot 10^{-3}$ 30 185,2 83,47% 82,10%
PTU $1 \cdot 10^{-4}$ 30 310,4 72,30% 71,50%
PTU $1 \cdot 10^{-3}$ 30 125,6 88,79% 87,90%
PTU $5 \cdot 10^{-3}$ 20 89,6 92,00% 91,40%
PTU $5 \cdot 10^{-3}$ 30 66,7 94,05% 93,80%
PTU $5 \cdot 10^{-3}$ 45 34,2 96,95% 96,20%
PTU $5 \cdot 10^{-3}$ 60 11,6 98,96% 98,50%
Urotropine (Chuẩn) $5 \cdot 10^{-3}$ 30 240,8 78,51% 77,20%

Thứ hai, ưu thế tuyệt đối của PTU so với chất ức chế công nghiệp truyền thống Urotropine: Luận án xác nhận:

"Khả năng ức chế ăn mòn thép của 1-phenyl-2-thiourea tốt hơn urotropine trong cả môi trường HCl 1,0 M và NaCl 3,5 % khi xét ở cùng điều kiện nồng độ và nhiệt độ."

Trong môi trường muối NaCl 3,5%, hiệu suất ức chế của PTU đạt 78,4% so với mức 54,2% của Urotropine ở cùng nồng độ $5 \cdot 10^{-3}$ M. Ảnh chụp kính hiển vi điện tử quét (SEM) sau 24 giờ ngâm mẫu thép khẳng định bề mặt kim loại khi có mặt PTU giữ được độ phẳng mịn nguyên vẹn, các vết rỗ và mảng ăn mòn sâu bị triệt tiêu hoàn toàn so với mẫu đối chứng ngâm trong HCl 1,0 M thuần túy.

Thứ ba, sự vượt trội về hoạt tính chống oxy hóa của dẫn xuất Selenourea so với Thiourea: Luận án chứng minh:

"1-phenyl-2-selenourea thể hiện khả năng chống oxy hóa tốt hơn 1-phenyl-2-thiourea. Trong các dẫn xuất của selenourea, các dẫn xuất chứa nhóm đẩy electron cho khả năng chống oxy hóa tốt hơn các dẫn xuất chứa nhóm hút electron."

Tính toán năng lượng phân ly liên kết BDE(N12-H13) và thế ion hóa (IE) ở mức lý thuyết ROB3LYP/6-311++G(2df,2p) chỉ ra rằng PSeU có giá trị BDE(N-H) thấp hơn PTU xấp xỉ 4,2 kcal/mol trong pha khí và 3,8 kcal/mol trong dung môi ethanol. Khi gắn nhóm thế cho điện tử mạnh $-\text{OCH}_3$ tại vị trí para (hợp chất $\text{CH}_3\text{OPSeU}$), giá trị thế ion hóa giảm mạnh từ 7,42 eV xuống còn 7,01 eV, thúc đẩy phản ứng diễn ra theo cả hai kênh HAT và chuyển electron đơn (SET).

Thứ tư, sàng lọc và thiết kế phân tử lưỡng chức năng $\text{CH}_3\text{OPSeU}$: Luận án đã xác lập phát hiện mang tính đột phá: 1-(4-methoxyphenyl)-2-selenourea ($\text{CH}3\text{OPSeU}$) là ứng viên hoàn hảo nhất với năng lượng liên kết hấp phụ trên bề mặt $\text{Fe}(110)$ trong môi trường acid đạt $E{\text{binding}} = -163{,}57$ kcal/mol (mạnh hơn đáng kể so với PTU là $-142{,}18$ kcal/mol). Đồng thời, tính toán động học dập tắt gốc $\text{HOO}^\bullet$ qua Eyringpy tại mức M05-2X/6-311++G(d,p) cho thấy hằng số tốc độ tổng cộng $k_{\text{tot}}$ đạt $1{,}34 \cdot 10^7\ \text{M}^{-1}\text{s}^{-1}$, chứng minh khả năng dập tắt gốc tự do cực kỳ nhanh trong dung dịch.

          BỀ MẶT THẾ NĂNG (PES) PHẢN ỨNG DẬP TẮT GỐC HOO•
   Năng lượng tự do Gibbs (kcal/mol)
                               Tọa độ phản ứng nội (IRC)

Implications đa chiều

  • Về mặt Lý thuyết: Luận án cung cấp một hệ thống dữ liệu lượng tử chuẩn xác về nhiệt động học liên kết (BDE, PA, IE, ETE) cho họ hợp chất selenourea thế para, đóng góp vào ngân hàng dữ liệu nhiệt hóa học quốc tế và củng cố phương pháp tiếp cận đa phiếm hàm (B3LYP kết hợp M05-2X) trong hóa học tính toán.
  • Về mặt Phương pháp luận: Xây dựng quy trình tích hợp khép kín giữa mô phỏng Monte Carlo, lý thuyết topo QTAIM và thực nghiệm phân cực điện hóa/quang phổ UV-Vis. Quy trình này có thể chuyển giao trực tiếp để nghiên cứu các hệ chất ức chế xanh chiết xuất từ tự nhiên hoặc vật liệu bán dẫn hữu cơ.
  • Về mặt Ứng dụng Thực tiễn: Hợp chất $\text{CH}_3\text{OPSeU}$ và dẫn xuất thiourea PTU là các ứng viên thương mại hóa tiềm năng cao trong công nghệ tẩy rửa gỉ sét đường ống công nghiệp bằng acid (acid pickling inhibitors) và phụ gia chống lão hóa, chống oxy hóa trong công nghiệp chế biến dầu nhờn và polyme sinh học.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án vẫn tồn tại một số giới hạn học thuật cần được mở rộng:

  1. Giới hạn về Bề mặt Tinh thể Kim loại: Các mô phỏng Monte Carlo và Động lực học phân tử mới chỉ tập trung vào mặt phẳng tinh thể $\text{Fe}(110)$ đơn tinh thể lý tưởng. Trong thực tế, thép carbon công nghiệp tồn tại nhiều pha đa tinh thể phức tạp ($\text{Fe}(100)$, $\text{Fe}(111)$, ranh giới hạt, và pha tạp chất carbide/cementite $\text{Fe}_3\text{C}$).
  2. Khảo sát In Vitro đối với Gốc Tự do: Các thử nghiệm thực nghiệm chống oxy hóa mới chỉ dừng lại ở các gốc tự do mô hình bền là DPPH• và $\text{ABTS}^{\bullet+}$. Hoạt tính bắt các gốc tự do oxy sinh học nguy hiểm trực tiếp như gốc hydroxyl ($\text{HO}^\bullet$), superoxide ($\text{O}_2^{\bullet-}$), và singlet oxygen ($^1\text{O}_2$) mới chỉ được tiếp cận qua mô hình tính toán lượng tử lý thuyết đối với $\text{HOO}^\bullet$.
  3. Độc tính Sinh thái của Hợp chất Chứa Selen: Luận án chưa thực hiện các thử nghiệm độc tính tế bào (cytotoxicity assays in vitro) và khả năng phân hủy sinh học của dẫn xuất selenourea trong môi trường tự nhiên.

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

  • Triển khai mô phỏng Động lực học lượng tử Ab Initio (AIMD) để khảo sát trực tiếp quá trình phân ly liên kết và truyền điện tử động lực học giữa phân tử ức chế và bề mặt hợp kim trong dung dịch điện phân theo thời gian thực (real-time scale fs/ps).
  • Tổng hợp hóa học xanh các dẫn xuất nano-selen gắn khung polymer sinh học (chitosan, alginate) nhằm tăng độ tan, giảm độc tính và tối ưu hóa tính năng bảo vệ kép.
  • Mở rộng thử nghiệm ức chế ăn mòn trong điều kiện công nghiệp khắc nghiệt (nhiệt độ cao $> 100\ ^\circ\text{C}$, áp suất cao, dòng chảy xoáy đa pha chứa khí xâm thực $\text{H}_2\text{S}$ và $\text{CO}_2$).

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

Công trình luận án mang lại những tác động sâu rộng trên cả phương diện học thuật lẫn kinh tế - xã hội:

                            SƠ ĐỒ TÁC ĐỘNG ĐA TẦNG
  • Tác động Học thuật: Cung cấp hệ phương pháp luận mẫu mực kết hợp QM/MD với thực nghiệm điện hóa, tạo tiền đề cho các nhóm nghiên cứu trong nước và quốc tế trích dẫn và phát triển trong các lĩnh vực liên ngành Hóa Lý thuyết - Khoa học Bề mặt - Hóa Dược.
  • Tác động Kinh tế Công nghiệp: Việc ứng dụng các dẫn xuất thiourea (PTU) với hiệu suất bảo vệ đạt đến 98,96% ở 60 °C giúp các nhà máy luyện kim và xử lý bề mặt kim loại giảm thiểu lượng tiêu hao acid tẩy gỉ, kéo dài tuổi thọ thiết bị, tiết kiệm hàng triệu USD chi phí bảo dưỡng định kỳ.
  • Tác động Xã hội và Môi trường: Góp phần giảm thiểu phát tán ion kim loại nặng ($\text{Fe}^{2+}/\text{Fe}^{3+}$) vào môi trường nước thải công nghiệp nhờ khả năng kìm hãm hiệu quả phản ứng hòa tan anốt của thép.

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

  • Nghiên cứu sinh & Nhà khoa học trẻ (PhD/Postdoc Researchers): Tiếp cận một khung phương pháp luận hoàn chỉnh từ thiết lập phần mềm tính toán (Gaussian, Materials Studio, AIM2000, Eyringpy) đến các kỹ thuật điện hóa chuyên sâu (PDP, EIS).
  • Các Giáo sư & Chuyên gia Hóa học Lý thuyết / Điện hóa: Sử dụng bộ dữ liệu nhiệt động, hằng số tốc độ phản ứng, và các thông số topo electron để hiệu chuẩn và phát triển các mô hình tính toán mới.
  • Kỹ sư R&D trong Ngành Luyện kim & Hóa chất: Ứng dụng trực tiếp công thức chất ức chế gốc thiourea/selenourea trong sản xuất dung dịch ức chế ăn mòn, sơn bảo vệ kim loại, và chất chống oxy hóa cho dầu động cơ.
  • Các Nhà hoạch định Chính sách & Cơ quan Quản lý Môi trường: Có thêm cơ sở khoa học để xây dựng tiêu chuẩn kỹ thuật về việc tái sử dụng và kiểm soát hóa chất ức chế ăn mòn thân thiện với môi trường trong công nghiệp.

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 mở rộng Lý thuyết Axit-Bazơ Cứng-Mềm (HSAB) của Pearson và Lý thuyết Mật độ Electron Phân tử Biên (Frontier Molecular Orbital Theory) sang trạng thái hấp phụ của các hợp chất chứa liên kết chalcogen ($C=S$, $C=Se$) trên bề mặt tinh thể kim loại chuyển tiếp $\text{Fe}(110)$. Luận án chứng minh rằng nguyên tử Se trong khung selenourea có độ mềm hóa học cao hơn S ($S_{\text{PSeU}} > S_{\text{PTU}}$), giúp tương tác cho - nhận điện tử giữa HOMO của chất ức chế và LUMO của mạng sắt diễn ra mạnh mẽ hơn, giảm thiểu năng lượng liên kết hấp phụ $E_{\text{binding}}$ xuống mức $-163{,}57$ kcal/mol.

2. Sự đổi mới về phương pháp luận nghiên cứu khi so sánh với ít nhất hai công trình quốc tế đi trước?

So với nghiên cứu của Li và cộng sự (2012) (chỉ sử dụng phương pháp đường cong phân cực để khảo sát allyl thiourea) và công trình của Hosseini và cộng sự (2009) (chỉ dừng lại ở mức tính toán bán thực nghiệm AM1/PM3), luận án của Đinh Quý Hương đã tạo nên bước đột phá vượt bậc:

  • Tích hợp đồng thời phiếm hàm Minnesota cấp cao M05-2X/6-311++G(d,p) với phần mềm Eyringpy để giải quyết động học dập tắt gốc tự do ở cấp độ chuyển động hạt nhân kèm hiệu chỉnh đường ngầm Wigner/Eckart.
  • Ứng dụng lý thuyết lượng tử nguyên tử trong phân tử (QTAIM) qua phần mềm AIM2000 để giải mã bản chất liên kết Hydro và tương tác nội phân tử tại các điểm tới hạn liên kết (BCP) của trạng thái chuyển tiếp (TS).

3. Phát hiện nào gây bất ngờ nhất (counter-intuitive) và giải thích cơ chế lý thuyết?

Phát hiện bất ngờ nhất là: Hiệu suất ức chế ăn mòn thép của 1-phenyl-2-thiourea trong môi trường acid HCl 1,0 M tăng mạnh khi nhiệt độ tăng (từ 92,00% ở 20 °C lên tới 98,96% ở 60 °C). Theo quy luật nhiệt động học thông thường đối với các chất ức chế hấp phụ vật lý, nhiệt độ tăng sẽ làm tăng động năng phân tử, thúc đẩy quá trình giải hấp phụ (desorption), dẫn đến suy giảm hiệu suất bảo vệ. Tuy nhiên, với PTU, nhiệt độ cao cung cấp năng lượng vượt qua hàng rào hoạt hóa để phân tử chuyển hóa hoàn toàn sang cơ chế hấp phụ hóa học (chemisorption), trong đó các liên kết phối trí bền vững được hình thành trực tiếp giữa S, N và các orbital d trống của Fe, tạo thành lớp màng bảo vệ thụ động sít chặt và bền vững nhiệt.

4. Luận án có cung cấp giao thức tái lập (replication protocol) chi tiết không?

Có. Luận án cung cấp chi tiết toàn bộ thông số kỹ thuật phục vụ việc tái lập độc lập:

  • Tọa độ không gian nguyên tử ($XYZ$ coordinates) của tất cả các cấu trúc tối ưu hóa hình học, trạng thái chuyển tiếp (TS), chất trung gian (Inter1, Inter2).
  • Thiết lập từ khóa tính toán chính xác trong Gaussian 09W (# opt freq=noraman rob3lyp/6-311++g(2df,2p) scrf=(iefpcm,solvent=ethanol)).
  • Giao thức xử lý điện cực thép (chuẩn bị bề mặt, đánh bóng, kích thước $1{,}0\ \text{cm}^2$, dải thế phân cực, tốc độ quét thế $1\ \text{mV/s}$, nồng độ hóa chất tinh khiết phân tích).

5. Khung chương trình nghiên cứu 10 năm (10-year research agenda) được định hình như thế nào?

Khung nghiên cứu 10 năm tập trung vào 3 trụ cột:

  1. Năm 1 - 3: Mở rộng thư viện phân tử selenourea gắn khung dị vòng đa nhân (benzothiazole, indole) và kiểm tra hoạt tính dập tắt các gốc tự do ROS trên mô hình tế bào sống in vitro.
  2. Năm 4 - 6: Ứng dụng kỹ thuật máy học (Machine Learning / QSAR) trên tập dữ liệu lớn các thông số hóa lượng tử để dự đoán tự động hiệu suất ức chế ăn mòn của hàng nghìn cấu trúc mới trước khi tổng hợp thực nghiệm.
  3. Năm 7 - 10: Phát triển công nghệ vi bọc (microencapsulation) các hợp chất $\text{CH}_3\text{OPSeU}$ vào lớp phủ thông minh tự hàn gắn (self-healing smart coatings) ứng dụng cho kết cấu thép công trình biển và giàn khoan dầu khí ngoài khơi.

Kết luận

Công trình luận án tiến sĩ của Đinh Quý Hương đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp học thuật cốt lõi:

  1. Khẳng định tính vượt trội của PTU: Chứng minh 1-phenyl-2-thiourea đạt hiệu suất ức chế ăn mòn thép 98,96% ở 60 °C trong môi trường HCl 1,0 M, vượt xa chất chuẩn công nghiệp Urotropine và dẫn xuất mạch hở ITU.
  2. Làm sáng tỏ cơ chế nhiệt động học: Thiết lập chính xác mô hình hấp phụ hóa học của các dẫn xuất thiourea trên bề mặt $\text{Fe}(110)$ tuân theo đẳng nhiệt hấp phụ Langmuir và Temkin.
  3. Giải mã động học chống oxy hóa chalcogen: Xác lập quy luật biến thiên hoạt tính dập tắt gốc tự do của họ selenourea theo cơ chế HAT và SPLET, khẳng định vai trò của nhóm thế cho điện tử para trong việc hạ thấp BDE và thế ion hóa IE.
  4. Phát hiện cấu trúc lưỡng chức năng đột phá: Thiết kế và chứng minh thành công phân tử 1-(4-methoxyphenyl)-2-selenourea ($\text{CH}_3\text{OPSeU}$) sở hữu đồng thời năng lượng liên kết hấp phụ kim loại cực đại ($-163{,}57$ kcal/mol) và hằng số tốc độ dập tắt gốc $\text{HOO}^\bullet$ đạt $1{,}34 \cdot 10^7\ \text{M}^{-1}\text{s}^{-1}$.
  5. Chuẩn hóa phương pháp luận liên ngành: Xây dựng thành công quy trình kết hợp khép kín giữa tính toán lượng tử vi mô (DFT, QTAIM, TST) và thực nghiệm hóa lý vĩ mô (PDP, EIS, SEM, UV-Vis).

Công trình mở ra các hướng nghiên cứu chuyển giao công nghệ vật liệu mới, đóng góp thiết thực vào sự phát triển của ngành Hóa học Vô cơ - Hóa Lý thuyết tại Việt Nam và hội nhập sâu rộng với dòng chảy học thuật quốc tế.