Tổng quan về luận án

Sự phá hủy vật liệu kim loại do ăn mòn điện hóa và sự thoái hóa tế bào do stress oxy hóa gây ra bởi các gốc oxy hoạt động (Reactive Oxygen Species - ROS) là hai thách thức khoa học và kinh tế mang tính toàn cầu. Theo thống kê của Tổ chức Ăn mòn Thế giới, tổn thất kinh tế do ăn mòn vật liệu ước tính chiếm khoảng 3,4% GDP toàn cầu (tương đương 2,5 nghìn tỷ USD). Song song đó, sự tích tụ các gốc tự do như peroxyl ($\text{ROO}^\bullet$), hydroxyl ($^\bullet\text{OH}$), và superoxide ($\text{O}_2^{\bullet-}$) là nguyên nhân gốc rễ gây đột biến cấu trúc sinh học và thoái hóa polyme. Luận án tiến sĩ hóa học 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" của nghiên cứu sinh Đinh Quý Hương (chuyên ngành Hóa vô cơ, mã số 9440113, thực hiện tại Trường Đại học Sư phạm – Đại học Huế dưới sự hướng dẫn khoa học của PGS.TS. Phạm Cẩm Nam và PGS.TS. Trần Dương) đã tiên phong giải quyết bài toán giao thoa này.

Điểm nghẽn nghiên cứu (research gap) cốt lõi được luận án định vị là: Các công trình trên thế giới trước đây thường tách biệt hoàn toàn giữa việc khảo sát chất ức chế ăn mòn kim loại (thường tập trung vào các dị vòng chứa N, S, O) và sàng lọc chất dập tắt gốc tự do chống oxy hóa (thường tập trung vào các hợp chất phenolic, hydroxylamine). Cho đến thời điểm luận án được triển khai, chưa có một nghiên cứu tích hợp có hệ thống nào thiết lập mối liên hệ định lượng giữa cấu trúc điện tử, cơ chế hóa học phân tử và tính năng kép: vừa ức chế ăn mòn thép trong môi trường xâm thực ($\text{HCl}\text{ 1,0 M}$ và $\text{NaCl}\text{ 3,5}%$) vừa chống oxy hóa dập tắt gốc tự do ở mức độ vi mô.

Luận án thiết lập hệ thống câu hỏi nghiên cứu và giả thuyết khoa học chặt chẽ:

  1. Câu hỏi nghiên cứu 1 (RQ1): Mối tương quan định lượng giữa cấu trúc phân tử, các tham số lượng tử (HOMO, LUMO, năng lượng ion hóa $\text{IE}$, năng lượng phân ly liên kết $\text{BDE}$) và hiệu suất ức chế ăn mòn thép cacbon trong dung dịch $\text{HCl}\text{ 1,0 M}$ và $\text{NaCl}\text{ 3,5}%$ của các dẫn xuất thiourea diễn ra theo quy luật nào?
    • Giả thuyết 1 (H1): Sự hiện diện của vòng phenyl liên hợp ($\pi$-electrons) và dị tố mang cặp electron tự do (N, S) trong 1-phenyl-2-thiourea (PTU) tạo mật độ điện tích chuyển dịch hiệu quả vào obitan d trống của sắt, đem lại hiệu suất bảo vệ bề mặt cao hơn cấu trúc mạch nhánh không liên hợp như 1,3-diisopropyl-2-thiourea (ITU).
  2. Câu hỏi nghiên cứu 2 (RQ2): Cơ chế chuyển giao electron/nguyên tử (HAT, SET-PT, SPLET, RAF) nào chi phối hoạt tính dập tắt gốc tự do của dẫn xuất thiourea và selenourea trong các pha môi trường khác nhau?
    • Giả thuyết 2 (H2): Việc thay thế nguyên tử lưu huỳnh bằng nguyên tố selen (Se) có độ âm điện thấp hơn, bán kính nguyên tử lớn hơn và giàu mật độ electron hơn sẽ làm giảm đáng kể rào cản năng lượng $\text{BDE}(\text{N}-\text{H})$ và $\text{IE}$, thúc đẩy phản ứng bắt gốc tự do qua cơ chế chuyển nguyên tử hydro (HAT) và chuyển electron đơn lẻ (SET).
  3. Câu hỏi nghiên cứu 3 (RQ3): Liệu có thể thiết kế hợp lý một phân tử đa chức năng mới thỏa mãn đồng thời tính năng ức chế ăn mòn kim loại xuất sắc và hoạt tính chống oxy hóa vượt trội bằng phương pháp mô hình hóa tính toán lượng tử kết hợp thực nghiệm?
    • Giả thuyết 3 (H3): Hợp chất 1-(4-methoxyphenyl)-2-selenourea ($\text{CH}_3\text{OPSeU}$) tích hợp nhóm đẩy electron $-\text{OCH}_3$ ở vị trí para sẽ tối ưu hóa cấu hình hấp phụ phẳng trên bề mặt $\text{Fe}(110)$ và hạ thấp rào cản Gibbs hoạt hóa ($\Delta G^\neq$) trong phản ứng dập tắt gốc hydroperoxyl ($\text{HOO}^\bullet$).

Khung lý thuyết của công trình tích hợp Lý thuyết Phiếm hàm Mật độ (Density Functional Theory - DFT) của Hohenberg-Kohn và Kohn-Sham, Lý thuyết Trạng thái Chuyển tiếp (Transition State Theory - TST), Thuyết Tương tác Biên obitan Phân tử (FMO Theory của Fukui), cùng Mô hình Nhiệt động học Hấp phụ Langmuir/Temkin và Động học Điện hóa Stern-Geary. Phạm vi nghiên cứu bao quát từ tính toán hóa lượng tử pha khí/pha dung môi phân cực liên tục ($\text{IEFPCM}$) trên các cụm nguyên tử mô phỏng bề mặt mạng tinh thể $\text{Fe}(110)$, kết hợp thực nghiệm phân cực thế động (PDP), phổ tổng trở điện hóa (EIS), kính hiển vi điện tử quét (SEM), đến các phép thử quang phổ dập tắt gốc tự do $\text{DPPH}^\bullet$ và $\text{ABTS}^{\bullet+}$.

Literature Review và Positioning

Nghiên cứu về chất ức chế ăn mòn kim loại và chất chống oxy hóa đã trải qua nhiều giai đoạn phát triển với những luồng học thuật chuyên biệt. Trong lĩnh vực ăn mòn kim loại, dòng nghiên cứu kinh điển của Sanyal (1981) và Hackerman (1987) đã xác lập vai trò của các phân tử hữu cơ chứa dị tố (N, O, S, P) trong việc chia sẻ cặp electron tự do hoặc electron $\pi$ vào obitan d trống của kim loại để hình thành liên kết phối trí bề mặt. Tiếp nối nền tảng đó, Li et al. (2012) khi nghiên cứu dẫn xuất allyl thiourea trong dung dịch $\text{H}_3\text{PO}_4$ đã chứng minh hiệu suất ức chế đạt trên $95%$ ở nồng độ $0,5\text{ mM}$. Tương tự, Hosseini et al. (2010) công bố hợp chất 1-methyl-3-pyridin-2-yl-thiourea thể hiện đặc tính ức chế hỗn hợp trong dung dịch $\text{H}_2\text{SO}_4$ với hiệu suất bảo vệ $96,4%$ ở nồng độ $100\text{ ppm}$. Tại Việt Nam, Trương Thị Thảo (2012) và Lục Văn Thụ (2012) đã tập trung khảo sát các cao chiết tự nhiên (chè, thuốc lá, polyphenol, caffein) để ức chế ăn mòn thép carbon trong môi trường axit. Tuy nhiên, phần lớn các công trình này chỉ dừng lại ở mô tả thông số vĩ mô của quá trình điện hóa hoặc sử dụng các phép tính bán thực nghiệm (như AM1 của Vũ Minh Tân et al., 2008) với độ chính xác hóa học chưa cao đối với cấu trúc ranh giới phân tử.

Trong lĩnh vực chống oxy hóa, các nghiên cứu cơ chế nhiệt động học và động học phát triển mạnh mẽ nhờ sự đóng góp của Blois (1958) với phép thử $\text{DPPH}^\bullet$, Re et al. (1999) với phép thử $\text{ABTS}^{\bullet+}$, và đặc biệt là Ingold & Litwinienko (2003, 2007) khi chuẩn hóa các con đường phản ứng: Chuyển nguyên tử hydro (Hydrogen Atom Transfer - HAT), Chuyển electron đơn lẻ kèm chuyển proton (Single Electron Transfer followed by Proton Transfer - SET-PT), Mất proton tuần tự kèm chuyển electron (Sequential Proton Loss Electron Transfer - SPLET), và Tạo sản phẩm cộng gốc tự do (Radical Adduct Formation - RAF). Về mặt hóa học tính toán, Alvarez-Idaboy et al. (2001) và Galano et al. (2011) đã chứng minh vai trò quyết định của phiếm hàm khuếch tán $\text{M05-2X}$ trong việc tính toán rào cản năng lượng hoạt hóa và hằng số tốc độ phản ứng dập tắt gốc tự do của các phân tử phenolic như Trolox. Đối với hợp chất chứa selen, Klotz (2003) và Fernández-Bolaños et al. (2012) phát hiện các dẫn xuất selenourea, selenocarbamate có hoạt tính bắt gốc tự do vượt trội với giá trị $\text{IC}_{50}$ đạt $19 - 46\text{ }\mu\text{M}$ trong phép thử $\text{DPPH}^\bullet$, nhờ tính ái nhân và khả năng cho electron ưu việt của nguyên tử Se so với S và O.

Tranh luận học thuật quốc tế nổi bật nằm ở câu hỏi: Liệu cơ chế bảo vệ kim loại của dẫn xuất thiourea bị chi phối bởi hấp phụ vật lý tĩnh điện hay hấp phụ hóa học tạo liên kết phối trí; và trong phản ứng chống oxy hóa, dung môi phân cực điều hướng phản ứng theo cơ chế HAT hay SPLET? Luận án của Đinh Quý Hương đã định vị chính xác ở giao điểm này, thiết lập một quy trình chuẩn hóa từ tính toán lượng tử mức cao ($\text{ROB3LYP/6-311++G(2df,2p)//B3LYP/6-311G(d,p)}$ và $\text{M05-2X/6-311++G(d,p)}$), mô phỏng động lực học phân tử/Monte Carlo bề mặt $\text{Fe}(110)$, đối chiếu trực tiếp với các thông số phân cực điện hóa thực nghiệm và phổ bắt gốc tự do quang học.

Đó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 thuyết tương tác bề mặt kim loại - phối tử và lý thuyết động học dập tắt gốc tự do qua các đóng góp cụ thể:

  1. Làm sâu sắc thuyết Hấp phụ Phối trí Kim loại: Luận án mở rộng mô hình tương tác biên của Fukui và Sanyal, chứng minh rằng mật độ điện tích Mulliken trên dị tố lưu huỳnh ($q_{\text{S}} = -0,266$ e ở PTU) kết hợp với hệ electron $\pi$ định cư trên vòng phenyl đóng vai trò trung tâm cho-nhận electron kép với mạng tinh thể sắt: Cặp electron tự do cho vào orbital $3d$ trống của $\text{Fe}$, đồng thời các orbital phản liên kết rỗng ($\pi^*$) tiếp nhận electron phản phối trí từ bề mặt kim loại, tạo liên kết hấp phụ hóa học bền vững.
  2. Chuẩn hóa ranh giới cơ chế chống oxy hóa HAT và SET-PT: Luận án xác định năng lượng phân ly liên kết $\text{BDE}(\text{N}-\text{H})$ là thông số nhiệt động học quyết định cơ chế dập tắt gốc trong pha khí và dung môi kém phân cực, trong khi năng lượng ion hóa $\text{IE}$ và ái lực proton $\text{PA}$ điều khiển phản ứng trong dung môi phân cực có khả năng sonvat hóa proton/electron.
  3. Mô hình hóa hiệu ứng thế nhân dị tố (S $\rightarrow$ Se) và nhóm thế para: Cung cấp bằng chứng lý thuyết định lượng rằng việc thay thế S bằng Se làm giảm mạnh giá trị $\text{BDE}(\text{N}-\text{H})$ từ $82,4\text{ kcal}\cdot\text{mol}^{-1}$ (ở PTU) xuống còn $75,8\text{ kcal}\cdot\text{mol}^{-1}$ (ở PSeU), đồng thời hạ thấp thế ion hóa $\text{IE}$, giải thích triệt để sự gia tăng vượt bậc của hoạt tính bắt gốc oxy hóa tự do.

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

Khung phân tích của luận án là sự hợp nhất của bốn trụ cột lý thuyết: (1) Thuyết Phiếm hàm Mật độ lượng tử (Conceptual DFT), (2) Lý thuyết Nguyên tử trong Phân tử (Quantum Theory of Atoms in Molecules - QTAIM của Bader), (3) Thuyết Động học Trạng thái Chuyển tiếp (TST kết hợp hiệu chỉnh đường ngầm Wigner qua phần mềm Eyringpy), và (4) Nhiệt động lực học Bề mặt Điện hóa (Adsorption Isotherms).

Phương pháp tiếp cận phân tích diễn ra qua 3 tầng nấc nghiêm ngặt:

  • Tầng 1 - Cấu trúc điện tử đơn phân tử: Xác định các tham số lượng tử toàn cục gồm năng lượng $E_{\text{HOMO}}$, $E_{\text{LUMO}}$, khoảng cách năng lượng $\Delta E = E_{\text{LUMO}} - E_{\text{HOMO}}$, độ cứng hóa học $\eta = (\text{IE}-\text{EA})/2$, độ mềm $S = 1/(2\eta)$, độ âm điện $\chi = (\text{IE}+\text{EA})/2$, thế hóa học $\mu = -\chi$, moment lưỡng cực $\mu_{\text{D}}$, và tỷ lệ electron chuyển dịch $\Delta N = (\Phi - \chi_{\text{inh}})/[2(\eta_{\text{Fe}} + \eta_{\text{inh}})]$ với hàm công thoát kim loại $\Phi_{\text{Fe}(110)} = 4,82\text{ eV}$.
  • Tầng 2 - Cấu trúc trạng thái chuyển tiếp và bề mặt: Khảo sát Tọa độ Phản ứng Nội (Intrinsic Reaction Coordinate - IRC) để tìm kiếm trạng thái chuyển tiếp (TS), tính toán rào cản năng lượng Gibbs hoạt hóa $\Delta G^\neq$ và hằng số tốc độ phản ứng dập tắt gốc tự do $k$, kết hợp tính toán topo mật độ electron $\rho(r)$, ma trận Hessian $\nabla^2\rho(r)$, mật độ thế năng $V(r)$ và động năng $G(r)$ tại các điểm tới hạn liên kết BCP (Bond Critical Points).
  • Tầng 3 - Điều kiện biên: Xác lập rõ ràng trong môi trường axit mạnh ($\text{HCl}\text{ 1,0 M}$), môi trường muối trung tính ($\text{NaCl}\text{ 3,5}%$), và dung môi ethanol đối với phản ứng gốc tự do.

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

Thiết kế nghiên cứu

Luận án tuân thủ triết lý nghiên cứu thực chứng nghiêm ngặt (positivism paradigm), kết hợp phương pháp luận tính toán mô phỏng lý thuyết từ nguyên lý ban đầu (ab initio) và thực nghiệm điện hóa - hóa sinh đa chiều (triangulation methodology). Thiết kế nghiên cứu tích hợp đa cấp độ: từ mức độ quỹ đạo electron hạ nano-mét ($10^{-10}\text{ m}$), cấu trúc tập hợp phân tử dung môi bề mặt ($10^{-9}\text{ m}$), đến bề mặt vi mô ăn mòn điện cực ($10^{-6}\text{ m}$) và mẫu đo thực tế ($10^{-2}\text{ m}$).

Quy trình nghiên cứu rigorous

  1. Giao thức Thực nghiệm Điện hóa:
    • Mẫu thép cacbon thấp chuẩn bị với diện tích làm việc chính xác $1,0\text{ cm}^2$, đúc trong nhựa epoxy, đánh bóng bằng giấy nhám SiC từ cỡ hạt 400 đến 1200 mesh, rửa sạch bằng nước cất hai lần, tẩy nhờn bằng axeton và sấy khô trước khi nhúng vào cell điện hóa 3 điện cực: Điện cực làm việc (WE - Thép), Điện cực đối trợ (CE - Lưới Bạch kim $\text{Pt}$), và Điện cực so sánh (RE - Bạc/Bạc Clorua $\text{Ag/AgCl}$ hoặc Calomel bão hòa).
    • Đo phân cực thế động (PDP): Quét thế trong dải $\pm 250\text{ mV}$ so với thế ăn mòn hở $E_{\text{corr}}$, tốc độ quét $1\text{ mV}\cdot\text{s}^{-1}$. Ngoại suy các đường dốc Tafel anốt ($\beta_{\text{a}}$) và catốt ($\beta_{\text{c}}$) để xác định mật độ dòng ăn mòn $i_{\text{corr}}$.
    • Phổ tổng trở điện hóa (EIS): Kích thích bằng tín hiệu điện thế xoay chiều hình sin biên độ $10\text{ mV}$ trong dải tần số từ $100\text{ kHz}$ đến $10\text{ mHz}$. Mô hình hóa giản đồ Nyquist và Bode bằng mạch tương đương chứa phần tử pha không đổi (Constant Phase Element - CPE).
    • Phân tích hình thái bề mặt: Sử dụng Kính hiển vi điện tử quét (SEM) với điện thế gia tốc chùm electron để đánh giá độ gồ ghề và mức độ tổn hại bề mặt thép sau $24\text{ giờ}$ ngâm trong dung dịch có và không có chất ức chế.
  2. Giao thức Thực nghiệm Chống Oxy hóa:
    • Thử nghiệm $\text{DPPH}^\bullet$: Trộn dung dịch mẫu thử ở các nồng độ khác nhau với dung dịch $\text{DPPH}^\bullet\text{ }0,1\text{ mM}$ trong ethanol, ủ tối ở $30^\circ\text{C}$ trong $30\text{ phút}$, đo mật độ quang tại bước sóng $\lambda = 517\text{ nm}$ bằng máy quang phổ UV-Vis. Xác định phần trăm bắt gốc $\text{SA}{\text{DPPH}}%$ và ngoại suy giá trị $\text{IC}{50}$.
    • Thử nghiệm $\text{ABTS}^{\bullet+}$: Tạo cation gốc tự do $\text{ABTS}^{\bullet+}$ bằng phản ứng giữa dung dịch $\text{ABTS}\text{ 7,0 mM}$ với $\text{K}_2\text{S}_2\text{O}8\text{ 2,45 mM}$, ủ tối $16\text{ giờ}$, pha loãng đến độ hấp thụ $A{734} = 0,70 \pm 0,02$, đo sự suy giảm quang phổ tại bước sóng $\lambda = 734\text{ nm}$ sau khi bổ sung chất chống oxy hóa.
  3. Quy trình Tính toán Hóa lượng tử và Mô phỏng Phân tử:
    • Tối ưu hóa cấu trúc hình học và tính tần số dao động điểm không ($\text{ZPE}$) ở mức lý thuyết $\text{B3LYP/6-311G(d,p)}$. Năng lượng điểm đơn hiệu chỉnh chính xác cao được tính ở mức $\text{ROB3LYP/6-311++G(2df,2p)}$. Hiệu ứng dung môi được tính toán thông qua mô hình phân cực liên tục $\text{IEFPCM}$.
    • Khảo sát động học phản ứng dập tắt gốc $\text{HOO}^\bullet$ sử dụng phiếm hàm lai hóa $\text{M05-2X}$ với bộ hàm $\text{6-311++G(d,p)}$. Trạng thái chuyển tiếp được kiểm tra nghiêm ngặt chỉ có duy nhất một tần số dao động ảo (imaginary frequency) và được xác nhận qua phân tích đường tọa độ phản ứng nội (IRC). Hằng số tốc độ phản ứng $k$ được hiệu chỉnh đường ngầm lượng tử bằng hệ số $\kappa(T)$ thông qua phần mềm Eyringpy.
    • Mô phỏng Monte Carlo và Động lực học phân tử (MD) được thực hiện trên phần mềm Materials Studio 8 với mặt phẳng tinh thể $\text{Fe}(110)$, hộp mô phỏng chứa các phân tử chất ức chế, nước và ion $\text{H}3\text{O}^+$, $\text{Cl}^-$ để tính toán năng lượng hấp phụ $E{\text{ads}}$ và năng lượng liên kết $E_{\text{binding}}$.

Data và phân tích

Hợp chất khảo sát Nồng độ ($M$) Môi trường ăn mòn Nhiệt độ ($^\circ\text{C}$) Mật độ dòng $i_{\text{corr}}$ ($\mu\text{A}\cdot\text{cm}^{-2}$) Hiệu suất ức chế $H%$ Cơ chế hấp phụ $\Delta G^\circ_{\text{ads}}$ ($\text{kJ}\cdot\text{mol}^{-1}$)
Mẫu trắng (Blank) 0 $\text{HCl}\text{ 1,0 M}$ 30 1120,5 - - -
ITU $5\cdot 10^{-3}$ $\text{HCl}\text{ 1,0 M}$ 30 185,3 83,46% Langmuir -28,45
PTU $5\cdot 10^{-3}$ $\text{HCl}\text{ 1,0 M}$ 20 89,6 92,00% Langmuir/Temkin -34,12
PTU $5\cdot 10^{-3}$ $\text{HCl}\text{ 1,0 M}$ 30 66,7 94,05% Langmuir/Temkin -35,80
PTU $5\cdot 10^{-3}$ $\text{HCl}\text{ 1,0 M}$ 45 34,2 96,95% Langmuir/Temkin -37,92
PTU $5\cdot 10^{-3}$ $\text{HCl}\text{ 1,0 M}$ 60 11,7 98,96% Langmuir/Temkin -41,05
Urotropine $5\cdot 10^{-3}$ $\text{HCl}\text{ 1,0 M}$ 30 240,8 78,51% Langmuir -25,60
PTU $5\cdot 10^{-3}$ $\text{NaCl}\text{ 3,5}%$ 30 42,1 76,20% Hỗn hợp -24,18
Urotropine $5\cdot 10^{-3}$ $\text{NaCl}\text{ 3,5}%$ 30 78,5 55,65% Hỗn hợp -19,85

Số liệu thực nghiệm khẳng định độ tin cậy thống kê cao với hệ số tương quan tuyến tính $R^2 > 0,99$ cho các phương trình đẳng nhiệt hấp phụ Langmuir.

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

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

  1. Hiệu năng ức chế ăn mòn tăng thuận theo nhiệt độ của PTU: Trái ngược với quy luật suy giảm hấp phụ vật lý thông thường ở các chất hữu cơ khi tăng nhiệt độ, hiệu suất ức chế của PTU trong $\text{HCl}\text{ 1,0 M}$ tăng lũy tiến từ $92,00%$ ($20^\circ\text{C}$) lên $94,05%$ ($30^\circ\text{C}$), $96,95%$ ($45^\circ\text{C}$) và đạt đỉnh $98,96%$ ở $60^\circ\text{C}$ tại nồng độ $5\cdot 10^{-3}\text{ M}$. Biến thiên năng lượng tự do Gibbs hấp phụ tiêu chuẩn $\Delta G^\circ_{\text{ads}}$ dao động từ $-34,12$ đến $-41,05\text{ kJ}\cdot\text{mol}^{-1}$ chứng minh quá trình chuyển hóa dứt khoát từ hấp phụ tĩnh điện sang hấp phụ hóa học tạo liên kết phối trí với bề mặt kim loại.
  2. Sự vượt trội của 1-phenyl-2-thiourea so với 1,3-diisopropyl-2-thiourea và Urotropine: PTU thể hiện khả năng bảo vệ thép vượt trội hơn ITU và chất ức chế truyền thống urotropine trong cả môi trường axit ($\text{HCl}\text{ 1,0 M}$) lẫn môi trường muối ($\text{NaCl}\text{ 3,5}%$). Nguyên nhân là do cấu trúc phẳng của vòng phenyl trong PTU giúp phân tử nằm áp sát song song với bề mặt tinh thể $\text{Fe}(110)$, tăng diện tích che phủ bề mặt $\theta$, tối đa hóa tương tác phân tán $\pi - d$.
  3. Đột phá hoạt tính chống oxy hóa khi chuyển vị từ Lưu huỳnh sang Selen: Hợp chất 1-phenyl-2-selenourea (PSeU) thể hiện khả năng quét gốc tự do $\text{DPPH}^\bullet$ và $\text{ABTS}^{\bullet+}$ vượt bậc so với chất tương đồng chứa lưu huỳnh PTU. Dữ liệu lượng tử chỉ ra rằng năng lượng phân ly liên kết $\text{BDE}(\text{N}{12}-\text{H}{13})$ của PSeU giảm xuống chỉ còn $75,8\text{ kcal}\cdot\text{mol}^{-1}$ (thấp hơn $6,6\text{ kcal}\cdot\text{mol}^{-1}$ so với PTU), giải phóng nguyên tử H dễ dàng hơn cho các gốc ROS.
  4. Quy luật thế para và phát hiện phân tử kép tối ưu $\text{CH}_3\text{OPSeU}$: Đối với các dẫn xuất selenourea thế ở vị trí para, các nhóm đẩy electron ($-\text{OCH}_3$, $-\text{CH}_3$) làm giảm đáng kể năng lượng ion hóa $\text{IE}$ và $\text{BDE}$, trong khi các nhóm hút electron ($-\text{NO}_2$, $-\text{CF}3$) gây hiệu ứng ngược lại. Hợp chất 1-(4-methoxyphenyl)-2-selenourea ($\text{CH}3\text{OPSeU}$) được chứng minh là phân tử đa năng xuất sắc nhất: Rào cản hoạt hóa phản ứng dập tắt gốc $\text{HOO}^\bullet$ qua cơ chế HAT giảm xuống chỉ còn $\Delta G^\neq = 7,32\text{ kcal}\cdot\text{mol}^{-1}$, hằng số tốc độ $k{\text{HAT}} = 4,28\cdot 10^7\text{ M}^{-1}\cdot\text{s}^{-1}$, đồng thời duy trì năng lượng hấp phụ cực mạnh trên bề mặt $\text{Fe}(110)$ ($E{\text{ads}} = -142,6\text{ kcal}\cdot\text{mol}^{-1}$).

Implications đa chiều

  • Về mặt Lý thuyết: Mở rộng thành công biên giới học thuật giữa hóa học lượng tử tính toán và điện hóa học ăn mòn, xác lập phương pháp luận dự đoán hoạt tính ức chế/chống oxy hóa thông qua các chỉ số descriptors mà không cần qua quy trình tổng hợp mù quáng.
  • Về mặt Phương pháp luận: Chuẩn hóa quy trình tính toán phối hợp $\text{ROB3LYP}$ (tính nhiệt động) và $\text{M05-2X}$ (tính động học) kết hợp hiệu chỉnh đường ngầm Wigner, có khả năng chuyển giao áp dụng cho hàng loạt hệ phân tử dị vòng khác.
  • Về mặt Thực tiễn và Ứng dụng Công nghiệp: Cung cấp giải pháp phụ gia ức chế ăn mòn thế hệ mới cho các hệ thống đường ống dẫn dầu mỏ, dung dịch tẩy rửa cáu cặn công nghiệp bằng axit $\text{HCl}$, hệ thống trao đổi nhiệt và kết cấu thép ven biển, giảm thiểu phụ thuộc vào các chất ức chế độc hại truyền thống như cromat hay hợp chất chứa kim loại nặng.
  • Về mặt Chính sách và Tiêu chuẩn Kỹ thuật: Đóng góp luận cứ khoa học để hoàn thiện tiêu chuẩn kiểm soát chất ức chế thân thiện môi trường trong ngành hóa dầu và luyện kim Việt Nam.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn nghiên cứu:

  1. Mô hình mô phỏng tinh thể bề mặt: Việc mô phỏng lượng tử bề mặt $\text{Fe}(110)$ sử dụng cụm tinh thể hữu hạn và màng dung dịch mỏng, chưa phản ánh đầy đủ khuyết tật mạng tinh thể thực tế (dislocations, grain boundaries) và dòng chảy thủy động lực học có áp suất.
  2. Môi trường pH kiềm và vi sinh vật: Nghiên cứu tập trung chủ yếu vào môi trường axit mạnh ($\text{HCl}\text{ 1,0 M}$) và muối ($\text{NaCl}\text{ 3,5}%$), chưa khảo sát chi tiết môi trường kiềm hóa hoặc môi trường ăn mòn sinh học (Microbiologically Influenced Corrosion - MIC) do vi khuẩn khử sunfat (SRB).
  3. Độ bền lâu dài và độc tính in vivo của Selen hữu cơ: Dù các hợp chất selenourea thể hiện hoạt tính hóa học xuất sắc, độ bền thủy phân dài hạn và hồ sơ độc tính sinh học toàn diện trên mô sống động vật cần được đánh giá lâm sàng sâu hơn.

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

  • Phát triển mô phỏng Động lực học Phân tử Phản ứng (ReaxFF MD) để mô tả chi tiết sự hình thành và đứt gãy liên kết theo thời gian thực trên quy mô hàng trăm nghìn nguyên tử.
  • Tổng hợp xanh và thử nghiệm bán công nghiệp dẫn xuất $\text{CH}_3\text{OPSeU}$ làm phụ gia chống ăn mòn cho dầu bôi trơn động cơ và sơn phủ bảo vệ công trình biển.
  • Đánh giá khả năng ức chế stress oxy hóa nội bào và bảo vệ tế bào thần kinh chống thoái hóa sinh học của các hợp chất selen hữu cơ.

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

  • Tác động Học thuật: Công trình đóng góp nền tảng dữ liệu hóa lượng tử và điện hóa chất lượng cao, tạo nguồn trích dẫn học thuật quan trọng cho các nghiên cứu thuộc nhóm ngành Hóa lý thuyết, Vật liệu Chống ăn mòn và Hóa Dược hữu cơ.
  • Tác động Ngành Công nghiệp: Giúp các doanh nghiệp sản xuất hóa chất phụ gia, lọc hóa dầu, đóng tàu biển và kết cấu công trình ngoài khơi định hình hướng phát triển sản phẩm ức chế ăn mòn hiệu năng cao, bền nhiệt độ và giảm thiểu tiêu hao kim loại.
  • Tác động Xã hội và Môi trường: Giảm thiểu thiệt hại kinh tế do sự cố ăn mòn gây rò rỉ hóa chất, bảo vệ tài nguyên sắt thép quốc gia, hướng tới mục tiêu phát triển bền vững và sản xuất hóa học xanh.

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

  • Nghiên cứu sinh và Giới Học thuật: Tiếp cận phương pháp luận tích hợp DFT, QTAIM, Eyringpy và kỹ thuật điện hóa để mở rộng cho các hệ phối tử - kim loại phức tạp.
  • Kỹ sư R&D Công nghiệp Vật liệu & Luyện kim: Nhận được các công thức phân tử tối ưu ($\text{PTU}$, $\text{CH}_3\text{OPSeU}$) với thông số liều lượng chính xác ($5\cdot 10^{-3}\text{ M}$) để phối chế chất ức chế ăn mòn.
  • Các nhà Hóa Dược và Công nghệ Sinh học: Khai thác cấu trúc dẫn xuất selenourea trong việc phát triển các hoạt chất quét gốc tự do ROS bảo vệ sức khỏe con người.

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ì? Trả lời: Đóng góp độc đáo nhất là việc mở rộng Thuyết Hấp phụ Phối trí Kim loại và Lý thuyết Tương tác Biên obitan (FMO) sang hệ thống phân tử đa chức năng: Luận án đã làm sáng tỏ cơ chế cho-nhận electron kép giữa hệ orbital $\pi$ liên hợp/cặp electron tự do của phối tử thiourea/selenourea với orbital $3d$ của bề mặt $\text{Fe}(110)$, giải thích bản chất bước nhảy từ hấp phụ vật lý sang hấp phụ hóa học tự phát khi tăng nhiệt độ ($\Delta G^\circ_{\text{ads}}$ âm hơn ở nhiệt độ cao).
  2. Tính đổi mới trong phương pháp luận so với các nghiên cứu trước đây? Trả lời: Khác với nghiên cứu bán thực nghiệm của Vũ Minh Tân et al. (2008) hoặc các khảo sát điện hóa đơn thuần của Li et al. (2012) và Hosseini et al. (2010), luận án thiết lập quy trình tính toán đa thang bậc nghiêm ngặt: Kết hợp $\text{ROB3LYP/6-311++G(2df,2p)}$ cho nhiệt động học, $\text{M05-2X/6-311++G(d,p)}$ kèm phần mềm Eyringpy để tính hằng số tốc độ động học $k$ có hiệu chỉnh đường ngầm Wigner, mô phỏng Monte Carlo/MD bề mặt tinh thể thực trên Materials Studio 8, và đối chứng đa biến với thực nghiệm PDP, EIS (Nyquist), SEM và UV-Vis ($\text{DPPH}/\text{ABTS}$).
  3. Phát hiện nào gây bất ngờ nhất và bằng chứng dữ liệu hỗ trợ? Trả lời: Phát hiện bất ngờ nhất là hiệu suất ức chế ăn mòn của 1-phenyl-2-thiourea (PTU) trong dung dịch $\text{HCl}\text{ 1,0 M}$ tỷ lệ thuận với nhiệt độ môi trường, tăng từ $92,00%$ tại $20^\circ\text{C}$ lên $98,96%$ tại $60^\circ\text{C}$ ở nồng độ $5\cdot 10^{-3}\text{ M}$. Dữ liệu nhiệt động học $\Delta H^\circ_{\text{ads}} > 0$ và $\Delta G^\circ_{\text{ads}}$ đạt $-41,05\text{ kJ}\cdot\text{mol}^{-1}$ ở $60^\circ\text{C}$ cung cấp bằng chứng thực nghiệm và lý thuyết không thể bác bỏ về hiện tượng hấp phụ hóa học thu nhiệt bền vững.
  4. Giao thức tái lặp nghiên cứu (Replication Protocol) có được cung cấp đầy đủ? Trả lời: Luận án cung cấp chi tiết toàn bộ tọa độ hình học phân tử tối ưu hóa, cơ sở dữ liệu hàm sóng, bộ thông số phiếm hàm, điều kiện biên mô phỏng dung môi $\text{IEFPCM}$, chuẩn độ nồng độ dung dịch thử nghiệm điện hóa ($\text{HCl}\text{ 1,0 M}$, $\text{NaCl}\text{ 3,5}%$), thông số quét phân cực ($1\text{ mV}\cdot\text{s}^{-1}$, $\pm 250\text{ mV}$), giải thuật khớp phổ tổng trở qua mạch tương đương, và quy trình pha chế gốc tự do $\text{DPPH}^\bullet/\text{ABTS}^{\bullet+}$, bảo đảm khả năng tái lập độc lập $100%$ trong các phòng thí nghiệm quốc tế.
  5. Chương trình nghiên cứu 10 năm được định hình ra sao? Trả lời: Mở rộng mô hình tính toán lượng tử lên hệ thống màng bảo vệ đa phân tử tự lắp ráp (Self-Assembled Monolayers - SAMs), thử nghiệm các phối tử phức cơ kim chứa selen trong môi trường dòng chảy áp suất cao của giếng khoan dầu khí, và tích hợp trí tuệ nhân tạo (Machine Learning Descriptors) để sàng lọc hàng triệu dẫn xuất dị vòng đa năng trong thời gian thực.

Kết luận

  1. Khẳng định 1-phenyl-2-thiourea (PTU) là chất ức chế ăn mòn thép xuất sắc trong môi trường $\text{HCl}\text{ 1,0 M}$, đạt hiệu suất bảo vệ $98,96%$ tại $60^\circ\text{C}$ ở nồng độ $5\cdot 10^{-3}\text{ M}$, vượt trội so với chất ức chế mạch nhánh ITU và chất ức chế chuẩn công nghiệp Urotropine.
  2. Xác định quá trình hấp phụ của PTU trên thép tuân theo quy luật đẳng nhiệt Langmuir và Temkin, có bản chất là hấp phụ hóa học thu nhiệt với biến thiên năng lượng tự do Gibbs $\Delta G^\circ_{\text{ads}}$ đạt $-41,05\text{ kJ}\cdot\text{mol}^{-1}$.
  3. Chứng minh hiệu ứng chuyển dịch dị tố từ Lưu huỳnh sang Selen (PSeU) làm giảm mạnh năng lượng liên kết $\text{BDE}(\text{N}-\text{H})$ xuống $75,8\text{ kcal}\cdot\text{mol}^{-1}$ và hạ thấp năng lượng ion hóa $\text{IE}$, nâng cao vượt bậc hoạt tính quét gốc tự do $\text{DPPH}^\bullet$ và $\text{ABTS}^{\bullet+}$ theo cơ chế chuyển nguyên tử hydro (HAT).
  4. Thiết lập quy luật ảnh hưởng của nhóm thế para trên nhân thơm của dẫn xuất selenourea: Nhóm đẩy electron ($-\text{OCH}_3$) thúc đẩy mạnh mẽ khả năng cho electron và nguyên tử hydro, trong khi nhóm hút electron làm suy giảm hoạt tính.
  5. Thiết kế và sàng lọc thành công phân tử đột phá 1-(4-methoxyphenyl)-2-selenourea ($\text{CH}_3\text{OPSeU}$) sở hữu tính năng kép hoàn hảo: Hằng số tốc độ dập tắt gốc $\text{HOO}^\bullet$ đạt $k = 4,28\cdot 10^7\text{ M}^{-1}\cdot\text{s}^{-1}$ và năng lượng hấp phụ bảo vệ bề mặt $\text{Fe}(110)$ đạt $-142,6\text{ kcal}\cdot\text{mol}^{-1}$.
  6. Xây dựng thành công khung phương pháp luận kết hợp chuẩn xác giữa tính toán hóa lượng tử ($\text{DFT/MD}$) và các phương pháp thực nghiệm hiện đại ($\text{PDP/EIS/SEM/UV-Vis}$), đóng góp một di sản học thuật quan trọng cho ngành Hóa học Vô cơ và Khoa học Vật liệu Chống ăn mòn.