Tổng quan về luận án

Hiện tượng ăn mòn kim loại và sự thoái hóa vật liệu do các tác nhân oxy hóa là những thách thức kỹ thuật và kinh tế nghiêm trọng hàng đầu trên quy mô toàn cầu. Theo báo cáo thống kê từ Tổ chức Ăn mòn Thế giới (World Corrosion Organization), tổng thiệt hại kinh tế trực tiếp do ăn mòn gây ra ước tính chiếm khoảng 3,4% tổng sản phẩm quốc nội (GDP) toàn cầu, tương đương 2,5 nghìn tỷ USD mỗi năm; riêng tại Hoa Kỳ, con số này đã vượt ngưỡng 1,1 nghìn tỷ USD (chiếm xấp xỉ 6,2% GDP). Song song với đó, quá trình hủy hoại sinh học và biến tính vật liệu gây nên bởi các gốc oxy hoạt động (Reactive Oxygen Species - ROS như $\text{HO}^\bullet$, $\text{HOO}^\bullet$, $\text{O}_2^{\bullet-}$) đòi hỏi các giải pháp can thiệp ở cấp độ phân tử. Luận án tiến sĩ hóa học của nghiên cứu sinh Đ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 Huế, 2020; dưới sự hướng dẫn khoa học của PGS. Trần Dương và PGS. Phạm Cẩm Nam) đã giải quyết trọn vẹn điểm nghẽn học thuật này bằng việc tích hợp mô hình hóa lượng tử đa thang đo với kiểm chứng thực nghiệm hiện đại nhằm phát triển các cấu trúc phân tử lưỡng chức năng (vừa ức chế ăn mòn kim loại, vừa dập tắt gốc tự do).

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực trạng các công trình quốc tế trước đây (điển hình như Li et al., 2011; Hosseini et al., 2009; Madkour et al., 2016; Klotz, 2003) thường khảo sát độc lập hai mảng: hoặc tối ưu hóa chất ức chế ăn mòn, hoặc nghiên cứu cơ chế chống oxy hóa sinh hóa. Việc thiếu vắng một khung lý thuyết thống nhất kết hợp mô phỏng hóa học tính toán (Computational Quantum Chemistry) với thực nghiệm điện hóa - quang phổ để thiết kế các dẫn xuất thiourea và selenourea đa năng đã tạo ra sự đứt gãy trong việc giải mã bản chất tương tác liên diện (interfacial interactions) và truyền điện tích.

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ẽ:

  • Câu hỏi nghiên cứu 1 (RQ1): Cấu trúc hình học và mật độ electron của các dị tố (N, S, Se) cùng vòng thơm trong dẫn xuất thiourea ảnh hưởng như thế nào đến khả năng hấp phụ và ức chế ăn mòn thép trong môi trường axit ($\text{HCl}\text{ 1,0 M}$) và môi trường muối ($\text{NaCl}\text{ 3,5%}$)?
  • Câu hỏi nghiên cứu 2 (RQ2): Động học và nhiệt động học của quá trình dập tắt gốc tự do ($\text{DPPH}^\bullet, \text{ABTS}^{\bullet+}, \text{HOO}^\bullet$) bởi dẫn xuất thiourea và selenourea tuân theo cơ chế ưu thế nào ($\text{HAT}, \text{SET-PT}, \text{SPLET}$ hay $\text{RAF}$)?
  • Câu hỏi nghiên cứu 3 (RQ3): Mối tương quan định lượng giữa cấu trúc - hoạt tính (QSAR) và các thông số Hammett hiệu chỉnh ($\sigma_p^+$) tác động ra sao đến năng lượng phân ly liên kết ($\text{BDE}$) và năng lượng ion hóa ($\text{IE}$)?
  • Câu hỏi nghiên cứu 4 (RQ4): Có thể thiết kế thành công một phân tử tối ưu tích hợp đồng thời hiệu suất bảo vệ bề mặt kim loại vượt trội và hoạt tính chống oxy hóa vượt bậc dựa trên khung phân tử selenourea hay không?

Hệ thống giả thuyết tương ứng:

  • Giả thuyết 1 (H1): Sự phối trí của cặp electron không liên kết từ dị tố S/Se và hệ electron $\pi$ liên hợp vào obitan d trống của nguyên tử Fe trên bề mặt $\text{Fe}(110)$ hình thành liên kết chemisorption bền vững, làm tăng hiệu suất bảo vệ khi nhiệt độ tăng.
  • Giả thuyết 2 (H2): Dẫn xuất selenourea có năng lượng ion hóa ($\text{IE}$) và $\text{BDE}(\text{N-H})$ thấp hơn dẫn xuất thiourea tương ứng, tạo ra hằng số tốc độ dập tắt gốc tự do cao hơn qua kênh chuyển nguyên tử hydro ($\text{HAT}$).
  • Giả thuyết 3 (H3): Các nhóm thế đẩy electron tại vị trí para của vòng phenyl sẽ làm giảm đáng kể rào cản năng lượng hoạt hóa ($\Delta G^\ddagger$) trong phản ứng với gốc hydroperoxyl ($\text{HOO}^\bullet$).
  • Giả thuyết 4 (H4): Phân tử 1-(4-methoxyphenyl)-2-selenourea ($\text{CH}_3\text{OPSeU}$) thể hiện năng lượng liên kết hấp phụ âm vượt trội trên bề mặt sắt đồng thời sở hữu hoạt tính dập tắt ROS cao nhất trong chuỗi nghiên cứu.

Khung lý thuyết nền tảng của công trình bao gồm Lý thuyết Phiếm hàm Mật độ (Density Functional Theory - DFT với các phiếm hàm Becke 3-parameter Lee-Yang-Parr - B3LYP, Restricted Open-shell B3LYP - ROB3LYP, và Minnesota M05-2X), Lý thuyết Trạng thái Chuyển tiếp (Transition State Theory - TST/Eyring), Lý thuyết Lớp điện kép và Động học Phân cực Điện hóa Stern-Geary, cùng Lý thuyết Nguyên tử trong Phân tử (Quantum Theory of Atoms in Molecules - QTAIM của Richard Bader).

Phạm vi nghiên cứu bao quát hệ mẫu thép cacbon thấp trong các dung dịch điện ly ăn mòn ($\text{HCl}\text{ 1,0 M}$ và $\text{NaCl}\text{ 3,5%}$) tại các dải nhiệt độ $20^\circ\text{C}, 30^\circ\text{C}, 45^\circ\text{C}, 60^\circ\text{C}$ với nồng độ chất ức chế từ $1 \times 10^{-4}\text{ M}$ đến $5 \times 10^{-3}\text{ M}$; khảo sát thực nghiệm quang phổ $\text{UV-Vis}$ với gốc $\text{DPPH}^\bullet$ và cation gốc $\text{ABTS}^{\bullet+}$; mô phỏng lượng tử pha khí và pha dung môi (ethanol, nước) bằng mô hình phân cực liên tục IEFPCM.

Literature Review và Positioning

Tổng quan y văn chỉ ra hai dòng nghiên cứu học thuật chủ đạo:

Dòng nghiên cứu ức chế ăn mòn kim loại: Khởi xướng từ các nghiên cứu nền tảng về chất ức chế hữu cơ chứa dị tố (Sanyal, 1981), các tác giả quốc tế đã chứng minh các hợp chất chứa đồng thời nitơ và lưu huỳnh mang lại hiệu quả vượt trội so với các hợp chất đơn lẻ (Madkour et al., 2016; Li et al., 2011; Hosseini et al., 2009). Điển hình, Hosseini et al. (2009) đã khảo sát 1-methyl-3-pyridin-2-yl-thiourea trong $\text{H}_2\text{SO}_4$, ghi nhận hiệu suất đạt 96,4% ở nồng độ $100\text{ ppm}$. Li et al. (2011) công bố allyl thiourea đạt hiệu suất 95% trong $\text{H}_3\text{PO}_4$ ở nồng độ $0,5\text{ mM}$. Tại Việt Nam, Trương Thị Thảo (2012) và Lục Văn Thụ (2012) nghiên cứu dịch chiết tự nhiên chứa polyphenol và cafein, song các hệ này thiếu tính ổn định hóa học và khó chuẩn hóa cơ chế lượng tử.

Dòng nghiên cứu chất chống oxy hóa: Phát triển từ phương pháp đo quang phổ gốc tự do $\text{DPPH}^\bullet$ của Blois (1958) và $\text{ABTS}^{\bullet+}$ của Re et al. (1999). Về mặt cơ chế, Ingold & Litwinienko (2003, 2007) đã hệ thống hóa các con đường: Chuyển nguyên tử Hydro ($\text{Hydrogen Atom Transfer - HAT}$), Chuyển electron đơn kèm chuyển proton ($\text{Single Electron Transfer followed by Proton Transfer - SETPT}$), và Mất proton tuần tự kèm chuyển electron ($\text{Sequential Proton Loss Electron Transfer - SPLET}$). Về hợp chất cơ-selen (organoselenium), Klotz (2003) và Fernández-Bolaños et al. (2012) chứng minh các dẫn xuất selenourea có khả năng bẫy gốc tự do mạnh mẽ với giá trị $\text{IC}_{50}$ đạt $19 - 46\text{ }\mu\text{M}$, mô phỏng hoạt tính enzyme glutathione peroxidase.

Tranh luận học thuật cốt lõi tồn tại ở hai khía cạnh:

  1. Bản chất hấp phụ nhiệt động: Quan điểm truyền thống cho rằng chất ức chế hữu cơ chủ yếu hấp phụ vật lý (physisorption) tĩnh điện, do đó hiệu suất bảo vệ sẽ giảm khi nhiệt độ tăng do sự giải hấp. Ngược lại, trường phái hóa học bề mặt hiện đại chứng minh các dẫn xuất có obitan liên hợp trải qua quá trình hấp phụ hóa học (chemisorption) thu nhiệt, tạo liên kết phối trí với kim loại, dẫn đến hiệu suất tăng theo nhiệt độ.
  2. Cơ chế phản ứng gốc tự do trong môi trường phân cực: Tồn tại sự bất đồng giữa vai trò của cơ chế $\text{HAT}$ thuần túy và cơ chế ion hóa phân ly $\text{SPLET}$ khi dung môi thay đổi từ không phân cực sang dung môi phân cực như ethanol hay nước.

Vị thế nghiên cứu của luận án được xác lập độc bản thông qua việc:

  • Hợp nhất phương pháp điện hóa thực nghiệm (PDP, EIS), hiển vi bề mặt (SEM) với hóa lượng tử nâng cao ($\text{ROB3LYP/6-311++G(2df,2p)//B3LYP/6-311G(d,p)}$ và $\text{M05-2X/6-311++G(d,p)}$).
  • So sánh định lượng trực tiếp với hai nghiên cứu quốc tế tiêu chuẩn: nghiên cứu của Madkour et al. (2016) trên 5 dẫn xuất phenylthiourea trong $\text{HCl}\text{ 2,0 M}$ (hiệu suất cao nhất 85,3% ở nồng độ $11 \times 10^{-3}\text{ M}$) và nghiên cứu của Alberto et al. (2013) về động học Trolox, từ đó chứng minh các dẫn xuất selenourea thế para trong luận án đạt tốc độ dập tắt gốc tự do vượt bậ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 Obitan Phân tử Biên (Frontier Molecular Orbital Theory) của Fukui và nguyên lý Độ cứng/Độ mềm Cặp Axit-Bazơ (HSAB) của Pearson trong động học ăn mòn kim loại:

  • Thiết lập mô hình chuyển dịch điện tích liên diện: Lượng tử hóa phân số electron trao đổi ($\Delta N$) từ phân tử ức chế sang vùng dẫn kim loại sắt: $$\Delta N = \frac{\chi_{\text{Fe}} - \chi_{\text{inh}}}{2(\eta_{\text{Fe}} + \eta_{\text{inh}})}$$
  • Mở rộng lý thuyết SPLET của Litwinienko & Ingold sang hệ dị tố chu kỳ cao (Selenium): Khẳng định sự thay thế đồng đẳng từ lưu huỳnh sang selen làm biến đổi căn bản năng lượng ion hóa ($\text{IE}$) và ái lực proton ($\text{PA}$), biến selenourea thành chất cho electron và hydro siêu việt.

Mô hình lý thuyết bao gồm 4 mệnh đề chuẩn hóa:

  • Mệnh đề 1 (P1): Giá trị $E_{\text{HOMO}}$ càng cao và khoảng cách năng lượng $\Delta E = E_{\text{LUMO}} - E_{\text{HOMO}}$ càng hẹp, khả năng nhường electron vào obitan d trống của sắt càng mạnh, tỷ lệ thuận với điện trở chuyển điện tích ($R_{\text{ct}}$).
  • Mệnh đề 2 (P2): Sự hiện diện của vòng benzen định hướng không gian phẳng làm tăng mật độ electron $\pi$, tạo tương tác đa tâm với bề mặt tinh thể $\text{Fe}(110)$.
  • Mệnh đề 3 (P3): Năng lượng phân ly liên kết $\text{BDE}(\text{N-H})$ là thông số nhiệt động học quyết định tốc độ phản ứng $\text{HAT}$; liên kết $\text{N-H}$ càng yếu, hằng số tốc độ $k_{\text{HAT}}$ càng lớn.
  • Mệnh đề 4 (P4): Nhóm thế đẩy electron ($-\text{OCH}_3$) tại vị trí para làm bền hóa gốc tự do trung gian thông qua hiệu ứng liên hợp, làm giảm năng lượng tự do hoạt hóa $\Delta G^\ddagger$.

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

Khung phân tích của luận án tích hợp bốn trụ cột lý thuyết:

  1. Hóa lượng tử cấu trúc tĩnh: Khảo sát các chỉ số phản ứng toàn cục ($E_{\text{HOMO}}, E_{\text{LUMO}}, \Delta E, \mu, \eta, S, \omega, \mu_D$) và điện tích cục bộ Mulliken.
  2. Nhiệt động học và Động học trạng thái chuyển tiếp: Tính toán chính xác $\text{BDE}, \text{IE}, \text{PA}, \text{ETE}, \Delta G^\ddagger$ kết hợp hiệu chỉnh đường ngầm Wigner/Eckart qua phần mềm Eyringpy.
  3. Mô phỏng Động lực học Phân tử (MD) & Monte Carlo: Định lượng năng lượng hấp phụ ($E_{\text{ads}}$) và cấu hình không gian của chất ức chế trên mạng tinh thể $\text{Fe}(110)$ tiếp xúc dung dịch axit bằng Materials Studio.
  4. Phân tích Tô-pô Mật độ Electron (QTAIM): Giải mã bản chất liên kết hidro và tương tác phối trí tại các điểm tới hạn liên kết (Bond Critical Points - BCP) thông qua mật độ electron $\rho(r)$, toán tử Laplacian $\nabla^2\rho(r)$ và mật độ năng lượng toàn phần $H(r) = G(r) + V(r)$.

Điều kiện biên xác lập: Áp dụng cho bề mặt kim loại nền sắt/thép cacbon thấp; dung dịch điện ly có tính xâm thực cao ($\text{Cl}^-$); nhiệt độ giới hạn từ $293,15\text{ K}$ đến $333,15\text{ K}$; dung môi mô phỏng là ethanol và nước.

                    ┌────────────────────────────────────────────────────────┐
                    │     KHUNG PHÂN TÍCH LƯỢNG TỬ - THỰC NGHIỆM ĐA TẦNG     │
                    └──────────────────────────┬─────────────────────────────┘
                                               │
               ┌───────────────────────────────┴───────────────────────────────┐
               ▼                                                               ▼
┌──────────────────────────────┐                               ┌──────────────────────────────┐
│     HÓA LƯỢNG TỬ & MÔ PHỎNG  │                               │    THỰC NGHIỆM ĐIỆN HÓA &    │
│            (DFT / MD)        │                               │       QUANG PHỔ HỌC          │
├──────────────────────────────┤                               ├──────────────────────────────┤
│ • Gaussian 09W / GaussView 5 │                               │ • Đường cong phân cực (PDP)  │
│   - B3LYP/6-311G(d,p)        │                               │   - Xác định icorr, Ecorr    │
│   - ROB3LYP/6-311++G(2df,2p) │                               │   - Độ dốc Tafel βa, βc      │
│   - M05-2X/6-311++G(d,p)     │                               │ • Tổng trở điện hóa (EIS)    │
│ • Materials Studio (MD/MC)   │   Tương quan đối sánh thực    │   - Điện trở Rct, tụ điện CPE│
│   - Hấp phụ trên Fe(110)     │ ◄───────────────────────────► │ • Hiển vi quét SEM           │
│ • AIM2000 (Bader QTAIM)      │      nghiệm - lý thuyết       │   - Hình thái học bề mặt     │
│   - Điểm tới hạn BCP, ρ(r)   │                               │ • Quang phổ UV-Vis           │
│ • Eyringpy: TST Kinetics     │                               │   - Đo IC50 với DPPH•, ABTS•+│
└──────────────────────────────┘                               └──────────────────────────────┘
               │                                                               │
               └───────────────────────────────┬───────────────────────────────┘
                                               ▼
                    ┌────────────────────────────────────────────────────────┐
                    │      THIẾT KẾ PHÂN TỬ LƯỠNG CHỨC NĂNG ĐỘT PHÁ          │
                    │   1-(4-methoxyphenyl)-2-selenourea (CH3OPSeU):         │
                    │    - Hiệu suất ức chế ăn mòn Fe(110) vượt bậc          │
                    │    - Dập tắt gốc tự do HOO• siêu tốc (kHAT cực đại)    │
                    └────────────────────────────────────────────────────────┘

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

Thiết kế nghiên cứu

Nghiên cứu vận hành theo thế giới quan Hiện thực Phê phán (Critical Realism) và Thực chứng (Positivism), đặt trọng tâm vào việc giải thích các quy luật nhân quả vi mô (cấu trúc orbital phân tử) chi phối các biểu hiện vĩ mô (tốc độ ăn mòn điện hóa và tốc độ dập tắt gốc tự do). Thiết kế hỗn hợp đa phương pháp (multi-method design) kết hợp đồng thời ba tầng mức:

  • Tầng vi mô: Tính toán lượng tử cấu trúc điện tử bằng các phiếm hàm DFT chuẩn xác cao.
  • Tầng trung mô (meso-scale): Mô phỏng Monte Carlo và Động lực học phân tử xác định năng lượng tương tác $E_{\text{inter}}$ trên giao diện $\text{Fe}(110)/\text{dung dịch}$.
  • Tầng vĩ mô: Đo đạc các thông số động học dòng ăn mòn điện hóa và độ hấp thụ quang học.

Quy trình nghiên cứu rigorous

  1. Chuẩn bị mẫu và đo điện hóa:
    • Điện cực làm việc (WE): Thép cacbon chế tạo dạng trụ tròn, đúc trong nhựa epoxy với diện tích bề mặt lộ sáng chính xác $1,0\text{ cm}^2$. Bề mặt kim loại được mài cơ học bằng giấy nhám SiC từ độ hạt 400 đến 1500 mesh, đánh bóng, tẩy dầu mỡ bằng axeton và rửa bằng nước cất hai lần.
    • Hệ đo điện hóa 3 điện cực: Điện cực so sánh Calomel bão hòa (SCE) hoặc $\text{Ag/AgCl}$, điện cực phụ trợ bằng lưới Bạch kim ($\text{Pt}$), bình điện phân hai lớp điều nhiệt.
    • Đo đường cong 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ở mạch ($E_{\text{corr}}$) với tốc độ quét $1\text{ mV/s}$.
    • Đo phổ tổng trở điện hóa (EIS): Biên độ tín hiệu kích thích xoay chiều $10\text{ mV}$ (peak-to-peak) trong dải tần số từ $100\text{ kHz}$ đến $10\text{ mHz}$.
  2. Khảo sát quang phổ hoạt tính chống oxy hóa:
    • Thử nghiệm $\text{DPPH}^\bullet$: Trộn dung dịch chất thử ở các nồng độ khác nhau với dung dịch $\text{DPPH}^\bullet$ trong ethanol, ủ tối $30\text{ phút}$, đo độ giảm mật độ quang tại bước sóng cực đại $\lambda_{\text{max}} = 517\text{ nm}$ trên máy quang phổ UV-Vis.
    • Thử nghiệm $\text{ABTS}^{\bullet+}$: Tạo cation gốc tự do bằng phản ứng giữa $\text{ABTS}$ và $\text{K}_2\text{S}_2\text{O}8$, pha loãng đạt độ hấp thụ chuẩn tại $\lambda{\text{max}} = 734\text{ nm}$, ghi nhận độ khử màu sau khi thêm chất nghiên cứu.
  3. Quy trình tính toán hóa lượng tử:
    • Tối ưu hóa cấu trúc hình học ở trạng thái dừng (stationary points) và phân tích tần số dao động (harmonic vibrational frequency) tại mức lý thuyết $\text{B3LYP/6-311G(d,p)}$ trên phần mềm Gaussian 09W. Xác nhận cấu trúc đáy không có tần số ảo, cấu trúc trạng thái chuyển tiếp ($\text{TS}$) có duy nhất 1 tần số ảo.
    • Phân tích tọa độ phản ứng nội (Intrinsic Reaction Coordinate - IRC) để kiểm chứng $\text{TS}$ kết nối chính xác chất đầu và sản phẩm.
    • Tính toán năng lượng đơn điểm chính xác cao tại mức $\text{ROB3LYP/6-311++G(2df,2p)}$ và $\text{M05-2X/6-311++G(d,p)}$ có hiệu chỉnh năng lượng điểm không ($\text{ZPE}$) và hiệu chỉnh nhiệt động học enthalpy/Gibbs free energy tại $298,15\text{ K}, 1\text{ atm}$.

Data và phân tích

Phân tích dữ liệu thực nghiệm và lý thuyết đạt độ tin cậy thống kê cao:

  • Dữ liệu đường cong phân cực được ngoại suy Tafel tuyến tính tự động để xác định mật độ dòng ăn mòn ($i_{\text{corr}}$), thế ăn mòn ($E_{\text{corr}}$), và độ dốc anốt/catốt ($\beta_a, \beta_c$). Hiệu suất ức chế $H%$ tính theo công thức: $$H% = \frac{i_{\text{corr}} - i_{\text{inh}}}{i_{\text{corr}}} \times 100%$$
  • Phổ EIS được khớp (fitting) vào mạch điện tương đương $R_s(R_{\text{ct}}Q_{\text{dl}})$ sử dụng phần tử pha không đổi ($\text{CPE} - Q$) để bù trừ cho độ nhám và tính không đồng nhất của bề mặt điện cực, với hệ số tương quan $R^2 > 0,999$.
  • Đường đẳng nhiệt hấp phụ Langmuir và Temkin được xử lý hồi quy tuyến tính với hệ số xác định $R^2 > 0,99$, cho phép rút ra các thông số nhiệt động then chốt: hằng số cân bằng hấp phụ ($K_{\text{ads}}$), biến thiên năng lượng tự do hấp phụ tiêu chuẩn ($\Delta G^\circ_{\text{ads}}$), enthalpy hấp phụ ($\Delta H^\circ_{\text{ads}}$) và entropy hấp phụ ($\Delta S^\circ_{\text{ads}}$).

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

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

Luận án đã đạt được 5 phát hiện mang tính đột phá khoa học:

  1. Hiệu suất ức chế vượt trội của 1-phenyl-2-thiourea (PTU): Trong môi trường $\text{HCl}\text{ 1,0 M}$, PTU ức chế ăn mòn thép cacbon vượt trội hoàn toàn so với dẫn xuất béo 1,3-diisopropyl-2-thiourea (ITU). Tại nồng độ $5 \times 10^{-3}\text{ M}$, hiệu suất ức chế của PTU tăng đơn điệu 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 $\text{HCl}\text{ 1,0 M}$ với hiệu suất đạt 92,00 % ở $20^\circ\text{C}$, 94,05 % ở $30^\circ\text{C}$, 96,95 % ở $45^\circ\text{C}$ và 98,96 % ở $60^\circ\text{C}$ ở nồng độ $5 \times 10^{-3}\text{ M}$." (Trích dẫn trực tiếp từ luận án). Đặc tính tăng hiệu suất theo nhiệt độ kèm giá trị $\Delta G^\circ_{\text{ads}} < -20\text{ kJ/mol}$ khẳng định quá trình hấp phụ mang bản chất hấp phụ hóa học tự phát (chemisorption).

  2. Cơ chế ức chế hỗn hợp và ưu thế môi trường: Đo điện hóa PDP chứng minh PTU làm dịch chuyển cả hai nhánh phân cực anốt và catốt, hoạt động như một chất ức chế hỗn hợp (mixed-type inhibitor). Giá trị điện trở chuyển điện tích $R_{\text{ct}}$ tăng mạnh từ $12,4\text{ }\Omega\cdot\text{cm}^2$ (mẫu trắng) lên $210,5\text{ }\Omega\cdot\text{cm}^2$ khi có mặt PTU ở $30^\circ\text{C}$. Đáng chú ý, khả năng bảo vệ của PTU trong môi trường axit ($\text{HCl}\text{ 1,0 M}$) cao hơn đáng kể so với môi trường muối ($\text{NaCl}\text{ 3,5%}$) do trong dung dịch axit, phân tử PTU bị proton hóa ($\text{pPTU}$) kết hợp hiệp đồng với ion $\text{Cl}^-$ hấp phụ trước trên bề mặt sắt tích điện dương, tạo lớp màng bảo vệ kép bền vững. Ảnh chụp SEM hiển thị bề mặt thép ngâm 24 giờ có PTU phẳng mịn, gần như không bị rỗ ăn mòn so với bề mặt bị phá hủy nặng nề của mẫu trắng.

  3. Vượt trội so với chất ức chế công nghiệp chuẩn: So sánh trực tiếp ở cùng điều kiện nồng độ ($5 \times 10^{-3}\text{ M}$) và nhiệt độ ($30^\circ\text{C}$), PTU cho hiệu suất bảo vệ đạt 94,05% trong $\text{HCl}\text{ 1,0 M}$ và 78,2% trong $\text{NaCl}\text{ 3,5%}$, vượt trội hoàn toàn so với Urotropine (chỉ đạt lần lượt 72,3% và 51,4%).

  4. Động học bắt gốc tự do và hiệu ứng dị tố Selen: Thực nghiệm $\text{DPPH}^\bullet$ và $\text{ABTS}^{\bullet+}$ cho thấy PTU có hoạt tính chống oxy hóa cao hơn ITU ($\text{IC}_{50}$ của PTU thấp hơn ITU 2,8 lần). Tính toán hóa lượng tử $\text{ROB3LYP/6-311++G(2df,2p)}$ giải thích tường tận: năng lượng phân ly liên kết $\text{BDE}(\text{N-H})$ của PTU thấp hơn ITU ($82,4\text{ kcal/mol}$ so với $91,2\text{ kcal/mol}$) nhờ gốc tự do tạo thành sau khi mất H được giải tỏa cộng hưởng trên hệ vòng phenyl. Đặc biệt, việc thay thế nguyên tử S bằng Se trong 1-phenyl-2-selenourea (PSeU) làm giảm mạnh năng lượng ion hóa $\text{IE}$ từ $7,82\text{ eV}$ xuống $7,15\text{ eV}$, đưa PSeU trở thành tác nhân chống oxy hóa mạnh hơn PTU.

  5. Thiết kế phân tử lưỡng chức năng $\text{CH}_3\text{OPSeU}$: Khảo sát chuỗi dẫn xuất selenourea thế para ($\text{R-C}_6\text{H}_4\text{-NH-C(=Se)-NH}_2$ với $\text{R} = -\text{OCH}_3, -\text{CH}_3, -\text{H}, -\text{Cl}, -\text{NO}_2$) chỉ ra mối tương quan tuyến tính chặt chẽ giữa hằng số Hammett $\sigma_p^+$ với $\text{IE}$ và $\text{BDE}$. Hợp chất 1-(4-methoxyphenyl)-2-selenourea mang nhóm $-\text{OCH}3$ đẩy electron mạnh nhất sở hữu $\text{BDE}(\text{N-H})$ thấp nhất ($76,8\text{ kcal/mol}$), rào cản năng lượng hoạt hóa $\Delta G^\ddagger$ với gốc $\text{HOO}^\bullet$ qua cơ chế $\text{HAT}$ chỉ $6,42\text{ kcal/mol}$ (tương ứng hằng số tốc độ $k{\text{HAT}} = 4,18 \times 10^7\text{ M}^{-1}\text{s}^{-1}$ ở $298,15\text{ K}$). Đồng thời, mô phỏng Monte Carlo trên mặt tinh thể $\text{Fe}(110)$ xác nhận phân tử $\text{CH}3\text{OPSeU}$ hấp phụ song song hoàn hảo với năng lượng hấp phụ cực tiểu $E{\text{ads}} = -148,6\text{ kcal/mol}$.

  HIỆU SUẤT ỨC CHẾ ĂN MÒN CỦA PTU TRONG HCl 1.0 M (Nồng độ 5.10^-3 M)
  ┌─────────────────────────────────────────────────────────────────┐
  │ 100% ┤                                                98.96%   │
  │      │                                      96.95%    ┌────┐   │
  │  95% ┤                            94.05%    ┌────┐    │    │   │
  │      │                  92.00%    ┌────┐    │    │    │    │   │
  │  90% ┤                  ┌────┐    │    │    │    │    │    │   │
  │      │                  │    │    │    │    │    │    │    │   │
  │  85% ┤                  │    │    │    │    │    │    │    │   │
  │      │  Mẫu trắng (0%)  │    │    │    │    │    │    │    │   │
  │   0% ┴──────┬───────────┴────┴────┴────┴────┴────┴────┴────┴───┤
  │            0 °C          20 °C     30 °C     45 °C     60 °C    │
  │                       Nhiệt độ môi trường phản ứng             │
  └─────────────────────────────────────────────────────────────────┘

Implications đa chiều

  • Về mặt lý thuyết: Cung cấp bằng chứng thực nghiệm và số liệu lượng tử chuẩn xác chứng minh cơ chế hấp phụ hóa học thu nhiệt của thiourea thơm; hoàn thiện mô hình tương quan Hammett-DFT trong dự đoán hoạt tính chống oxy hóa của các hợp chất cơ-selen.
  • Về mặt phương pháp luận: Chuẩn hóa quy trình tích hợp giữa tính toán hóa lượng tử đa cấu hình (DFT/Eyringpy/QTAIM) với mô phỏng giao diện kim loại (Materials Studio) và thử nghiệm điện hóa thực, thiết lập tiêu chuẩn nghiên cứu mới cho lĩnh vực thiết kế vật liệu bảo vệ bề mặt.
  • Về mặt ứng dụng thực tiễn: Định hướng quy trình tổng hợp chất ức chế ăn mòn thế hệ mới thay thế các chất độc hại truyền thống (như hợp chất cromat, amin bậc cao độc tính); mở ra ứng dụng phụ gia đa chức năng vừa chống oxy hóa dầu bôi trơn vừa chống ăn mòn đường ống dẫn trong công nghiệp lọc hóa dầu và hàng hải.

Limitations và Future Research

Luận án thẳng thắn thừa nhận 4 hạn chế học thuật:

  1. Mô hình hóa bề mặt lý tưởng hóa: Mô phỏng Monte Carlo/Động lực học phân tử thực hiện trên mặt đơn tinh thể phẳng $\text{Fe}(110)$, chưa xét đến độ gồ ghề thực tế, các biên hạt (grain boundaries), và sự hiện diện của các pha tạp chất cacbit ($\text{Fe}_3\text{C}$) trong thép cacbon đa tinh thể.
  2. Xử lý dung môi liên tục: Mô hình IEFPCM chưa mô tả được chi tiết các liên kết hidro vi mô cục bộ giữa phân tử dung môi với các tâm phản ứng dị tố của chất ức chế.
  3. Giới hạn thử nghiệm động lực học dòng chảy: Thử nghiệm điện hóa thực hiện trong điều kiện dung dịch tĩnh, chưa khảo sát điều kiện dòng chảy rối (turbulent flow) và ứng suất cắt trong đường ống công nghiệp.
  4. Độc tính sinh thái và chi phí tổng hợp của dẫn xuất Selen: Mặc dù selenourea thể hiện hoạt tính vượt bậc, tính bền oxy hóa trong không khí lâu dài và chi phí nguyên liệu tổng hợp hợp chất cơ-selen quy mô công nghiệp cần được đánh giá bổ sung.

Chương trình nghiên cứu tương lai (Future Agenda):

  • Phát triển mô hình cơ lượng tử đa tỷ lệ (QM/MM) kết hợp trường lực phản ứng ReaxFF để mô phỏng động học phá hủy màng ức chế dưới dòng chảy ăn mòn xâm thực.
  • Ứng dụng thuật toán Máy học (Machine Learning) và mạng nơ-ron tích chập trên tập dữ liệu cấu trúc DFT để sàng lọc tự động hàng nghìn dẫn xuất selenourea tiềm năng.
  • Triển khai thử nghiệm thực tế tại hiện trường trên các hệ thống tháp giải nhiệt công nghiệp và bồn chứa xăng dầu.
  • Khảo sát độc tính tế bào (in vitro cytotoxicity) của hợp chất $\text{CH}_3\text{OPSeU}$ nhằm mở rộng ứng dụng trong bảo quản dược - thực phẩm.

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

  • Tác động học thuật: Dự báo tạo xung lực trích dẫn cao trong các tạp chí quốc tế chuyên ngành (Q1/Q2) thuộc hệ thống ISI/Scopus như Corrosion Science, Electrochimica Acta, Journal of Molecular Liquids, và Food Chemistry.
  • Chuyển đổi công nghiệp: Cung cấp công thức phụ gia ức chế ăn mòn cho ngành luyện kim, chế tạo máy và công nghiệp khai thác dầu khí biển. Việc áp dụng dẫn xuất thiourea/selenourea với liều lượng cực nhỏ ($5\text{ mM}$) cho phép nâng cao tuổi thọ thiết bị trao đổi nhiệt thêm $30 - 50%$, giảm chi phí bảo dưỡng định kỳ.
  • Lợi ích xã hội và môi trường: Góp phần giảm thiểu tiêu hao kim loại, ngăn ngừa rò rỉ hóa chất và ô nhiễm môi trường do bục vỡ đường ống công nghiệp; giảm thiểu việc sử dụng các hóa chất chống gỉ gốc kim loại nặng độc hại.

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

  • Nghiên cứu sinh & Nhà nghiên cứu Hóa học/Vật liệu: Tiếp cận phương pháp luận mẫu mực kết hợp DFT với thực nghiệm phân cực và tổng trở điện hóa.
  • Giảng viên & Nhà khoa học cao cấp: Tài liệu tham khảo giá trị về nhiệt động học hấp phụ, động học trạng thái chuyển tiếp và hóa học các hợp chất cơ-chalcogen (S, Se).
  • Bộ phận R&D doanh nghiệp Luyện kim, Sơn mạ & Hóa chất: Nắm bắt công nghệ cốt lõi để pha chế các chế phẩm tẩy rửa kim loại trong môi trường axit ($\text{pickling inhibitors}$) và dung dịch làm mát chống ăn mòn.
  • Nhà hoạch định chính sách & Quản lý môi trường: Căn cứ khoa học để xây dựng tiêu chuẩn kỹ thuật về chất chống ăn mòn an toàn, thân thiện môi trường trong công nghiệp nặng.

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?

Luận án đã mở rộng Thuyết Obitan Phân tử Biên (FMO) và Mô hình SPLET (Sequential Proton Loss Electron Transfer) của Litwinienko & Ingold sang hệ dị tố cơ-selen. Điểm độc đáo nhất là chứng minh sự thay thế nguyên tử lưu huỳnh bằng selen trong khung phân tử thiourea không chỉ làm tăng khả năng cho electron từ obitan HOMO lên vùng d của kim loại sắt, mà còn làm suy yếu liên kết $\text{N-H}$ liền kề thông qua hiệu ứng trường lượng tử, biến phân tử thành chất dập tắt gốc tự do $\text{HOO}^\bullet$ siêu tốc qua đồng thời hai cơ chế $\text{HAT}$ và $\text{RAF}$.

2. Đổi mới phương pháp luận thể hiện như thế nào so với hai nghiên cứu quốc tế điển hình?

So với nghiên cứu của Madkour et al. (2016) (chỉ dùng thực nghiệm điện hóa cổ điển khảo sát phenylthiourea) và nghiên cứu của Alberto et al. (2013) (chỉ tính toán DFT cho phân tử Trolox), luận án đã:

  • Thiết lập quy trình khép kín: Tính toán vi mô ($\text{DFT B3LYP/M05-2X}$) $\rightarrow$ Mô phỏng trung mô (Monte Carlo/MD trên tinh thể $\text{Fe}(110)$) $\rightarrow$ Đo đạc thực nghiệm vĩ mô (PDP, EIS, SEM, quang phổ UV-Vis).
  • Áp dụng mức lý thuyết hiệu chỉnh tương quan cao $\text{ROB3LYP/6-311++G(2df,2p)}$ kết hợp phân tích QTAIM để giải mã chính xác bản chất liên kết hydro và rào cản năng lượng hoạt hóa.

3. Phát hiện bất ngờ nhất được hỗ trợ bởi dữ liệu là gì?

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) tăng tỉ lệ thuận với nhiệt độ trong $\text{HCl}\text{ 1,0 M}$ (từ 92,00% ở $20^\circ\text{C}$ lên 98,96% ở $60^\circ\text{C}$). Kết quả này phản trực giác so với các chất ức chế hữu cơ thông thường (thường bị giải hấp khi tăng nhiệt độ) và là bằng chứng đanh thép chứng minh PTU đã trải qua quá trình hấp phụ hóa học thu nhiệt mạnh mẽ trên bề mặt thép, kích hoạt bởi sự phối trí electron $\pi$ của vòng thơm.

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

Hoàn toàn có. Luận án cung cấp chi tiết:

  • Thông số chế tạo điện cực thép ($1\text{ cm}^2$), quy trình xử lý cơ học bề mặt và nồng độ pha chế chính xác dung dịch chất ức chế ($1 \times 10^{-4}\text{ M}$ đến $5 \times 10^{-3}\text{ M}$).
  • Điều kiện đo đạc điện hóa (tốc độ quét $1\text{ mV/s}$, biên độ EIS $10\text{ mV}$, dải tần $100\text{ kHz} - 10\text{ mHz}$).
  • Bộ từ khóa và mức lý thuyết tính toán trên Gaussian 09W (# opt freq b3lyp/6-311g(d,p), # rob3lyp/6-311++g(2df,2p) scrf=(iefpcm,solvent=ethanol)), tọa độ không gian nguyên tử và thông số mô phỏng Monte Carlo trên Materials Studio.

5. Chương trình nghiên cứu 10 năm được phác thảo ra sao?

Chương trình 10 năm tập trung vào 4 giai đoạn:

  • Năm 1-3: Phát triển thư viện chất ức chế cơ-selen đa vòng và tối ưu hóa quy trình tổng hợp xanh (green synthesis).
  • Năm 4-5: Ứng dụng trí tuệ nhân tạo (AI/Machine Learning) dự đoán hoạt tính ức chế và độ bền màng tự lắp ghép (SAMs).
  • Năm 6-8: Nghiên cứu công nghệ đóng gói nano (nano-encapsulation) tạo lớp phủ thông minh tự phục hồi (self-healing coatings).
  • Năm 9-10: Thử nghiệm pilot quy mô công nghiệp trong ngành vận tải biển và chế biến dầu khí.

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 6 đóng góp cốt lõi:

  1. Xác định 1-phenyl-2-thiourea (PTU) là chất ức chế ăn mòn thép cacbon xuất sắc trong $\text{HCl}\text{ 1,0 M}$, đạt hiệu suất bảo vệ $98,96%$ ở $60^\circ\text{C}$ ($5 \times 10^{-3}\text{ M}$), vượt trội hoàn toàn so với ITU và chất công nghiệp Urotropine.
  2. Làm sáng tỏ cơ chế ức chế hỗn hợp của dẫn xuất thiourea, chứng minh bản chất hấp phụ hóa học tuân theo phương trình đẳng nhiệt Langmuir và Temkin.
  3. Giải mã vai trò quyết định của môi trường axit trong việc thúc đẩy cơ chế hấp phụ cầu nối ion giữa phân tử proton hóa và bề mặt kim loại qua trung gian anion $\text{Cl}^-$.
  4. Khẳng định ưu thế vượt trội của dẫn xuất selenourea so với thiourea trong khả năng bẫy gốc tự do $\text{DPPH}^\bullet, \text{ABTS}^{\bullet+}, \text{HOO}^\bullet$, thiết lập định luật định lượng tương quan giữa hằng số Hammett $\sigma_p^+$ với $\text{IE}$ và $\text{BDE}$.
  5. Thiết kế thành công phân tử đột phá 1-(4-methoxyphenyl)-2-selenourea ($\text{CH}3\text{OPSeU}$) tích hợp hoàn hảo hai chức năng: bám dính ức chế ăn mòn bề mặt $\text{Fe}(110)$ cực mạnh ($E{\text{ads}} = -148,6\text{ kcal/mol}$) và dập tắt gốc $\text{HOO}^\bullet$ siêu tốc qua cơ chế $\text{HAT}$ ($k_{\text{HAT}} = 4,18 \times 10^7\text{ M}^{-1}\text{s}^{-1}$).
  6. Thiết lập một chuẩn mực nghiên cứu mới kết hợp chặt chẽ giữa hóa lượng tử độ chính xác cao và thực nghiệm điện hóa - quang phổ học, mở ra hướng đi đầy tiềm năng cho việc phát triển các vật liệu thông minh bảo vệ kim loại và chống lão hóa sinh hóa trong tương lai.