Tổng quan luận án

Tính cấp thiết và khoảng trống nghiên cứu

Cùng với sự phát triển của nền công nghiệp, nhu cầu năng lượng ngày càng gia tăng trong khi các nguồn nhiên liệu hóa thạch truyền thống như than đá và dầu mỏ dần cạn kiệt. Việc đốt cháy nhiên liệu hóa thạch giải phóng lượng lớn khí thải carbon monoxide (CO) và carbon dioxide ($\text{CO}_2$), gây ô nhiễm môi trường và ảnh hưởng xấu đến sức khỏe con người. Việc chuyển hóa các khí thải chứa carbon thành nhiên liệu sạch thông qua hỗn hợp syngas ($\text{CO} + \text{H}_2$) đóng vai trò quan trọng nhằm quay vòng chu trình carbon và tận dụng nguồn nguyên liệu sẵn có từ sinh khối, khí tự nhiên, than đá hoặc rác thải hữu cơ.

Trong các hướng chuyển hóa syngas, quá trình tổng hợp alcohol mạch cao (Higher Alcohol Synthesis - HAS, chứa từ hai nguyên tử carbon trở lên) có ý nghĩa thực tiễn lớn do các sản phẩm này sở hữu phân tử khối cao, áp suất hơi thấp, nhiệt độ sôi cao, khả năng hòa tan tốt với hydrocarbon và khả năng chịu nước cao, phù hợp làm nhiên liệu trực tiếp hoặc phụ gia tăng chỉ số octane. Tuy nhiên, liên kết $\text{C}\equiv\text{O}$ có độ bền nhiệt động rất cao, khiến quá trình đòi hỏi năng lượng hoạt hóa lớn. Các hệ xúc tác kim loại chuyển tiếp đơn kim loại thường gặp hạn chế về độ chọn lọc: xúc tác trên cơ sở Ni và Co có xu hướng thúc đẩy phản ứng methane hóa tạo $\text{CH}_4$, trong khi xúc tác Cu có hoạt tính hydro hóa tạo alcohol nhưng cần tối ưu hóa để tăng hiệu suất tạo mạch $\text{C}_2^+$.

Khoảng trống nghiên cứu lớn được tác giả chỉ ra là: phần lớn các nghiên cứu lý thuyết trước đây chỉ tập trung vào các hệ xúc tác kim loại đơn lẻ (Ni, Co, Cu). Các nghiên cứu tính toán lượng tử về cơ chế phản ứng syngas trên các hệ xúc tác đa thành phần phức tạp — đặc biệt là các cụm (cluster) lưỡng kim loại phân tán trên các chất mang oxide và carbon — còn rất hạn chế. Do đó, việc nghiên cứu có hệ thống ở cấp độ phân tử nhằm làm rõ vai trò tương tác giữa các tâm kim loại và ảnh hưởng của chất mang là một đòi hỏi cấp thiết.

Mục tiêu nghiên cứu

Mục đích của luận án là sử dụng các phương pháp hóa học tính toán để khảo sát cơ chế phản ứng hydrogen hóa CO trên các hệ xúc tác cluster kim loại chuyển tiếp Ni, Cu, Co, các hệ xúc tác lưỡng kim loại Ni-Cu, Co-Cu phân tán trên các chất mang than hoạt tính (AC), magnesium oxide (MgO) và nhôm oxide ($\text{Al}_2\text{O}_3$); so sánh và làm rõ vai trò của các tâm xúc tác đơn/lưỡng kim loại; đồng thời làm rõ vai trò của các chất mang trong việc định hướng phản ứng và độ chọn lọc sản phẩm.

Nhiệm vụ nghiên cứu

Luận án xác lập hai nhiệm vụ nghiên cứu chính:

a) Nghiên cứu tài liệu và xây dựng tổng quan:

  1. Tổng hợp cơ sở lý thuyết hóa học lượng tử, nhiệt động lực học và động hóa học liên quan; hệ thống hóa phương pháp phiếm hàm mật độ (DFT), phương pháp dải đàn hồi leo dốc (CI-NEB), phương pháp động lực phân tử (MD) và mô phỏng Monte Carlo.
  2. Đánh giá tình hình nghiên cứu thực nghiệm và lý thuyết về phản ứng chuyển hóa syngas trên các hệ xúc tác trong nước và trên thế giới, xác định các vấn đề còn tồn đọng.

b) Thực hiện tính toán cơ chế phản ứng hydrogen hóa CO trên các hệ xúc tác:

  1. Xây dựng mô hình và tối ưu hóa cấu trúc hình học của các phân tử $\text{CO}$, $\text{H}_2$, các cluster $\text{Ni}_4$, $\text{Cu}_4$, $\text{Co}_4$, $\text{Ni}_2\text{Cu}_2$, $\text{Co}_2\text{Cu}_2$, các chất mang $\text{AC}$, $\text{MgO}$, $\text{Al}_2\text{O}_3$ và các hệ xúc tác mang kim loại tương ứng.
  2. Khảo sát quá trình hấp phụ và hoạt hóa của $\text{CO}$ và $\text{H}2$ trên các hệ xúc tác: tính toán năng lượng hấp phụ ($E{ads}$), biến thiên năng lượng ($\Delta E$), phân bố điện tích ($Q$), độ dài liên kết ($d$), bậc liên kết Mayer ($B$) và mật độ trạng thái (DOS) để xác định bản chất hấp phụ vật lý hay hóa học.
  3. Khảo sát các đường phản ứng chuyển hóa $\text{CO}$ tạo thành methanol ($\text{CH}_3\text{OH}$), ethanol ($\text{CH}_3\text{CH}_2\text{OH}$), methane ($\text{CH}_4$), formaldehyde ($\text{HCHO}$), acetaldehyde ($\text{CH}_3\text{CHO}$)...; xác định trạng thái chuyển tiếp (TS), sản phẩm trung gian (IS) và xây dựng bề mặt thế năng (PES).
  4. So sánh khả năng hoạt động, năng lượng hoạt hóa ($E_a$) và độ chọn lọc sản phẩm giữa các hệ xúc tác để đề xuất hệ xúc tác tối ưu.

Đối tượng và phạm vi nghiên cứu

  • Chất phản ứng ban đầu: Các phân tử $\text{CO}$ và $\text{H}_2$.
  • Sản phẩm khảo sát: Methane ($\text{CH}_4$), methanol ($\text{CH}_3\text{OH}$), ethanol ($\text{C}_2\text{H}_5\text{OH}$), formaldehyde ($\text{HCHO}$), formic acid ($\text{HCOOH}$), acetaldehyde ($\text{CH}_3\text{CHO}$)...
  • Mô hình cluster kim loại: Cụm đơn kim loại $\text{Ni}_4$, $\text{Cu}_4$, $\text{Co}_4$ và cụm lưỡng kim loại $\text{NiCu}$, $\text{Ni}_2\text{Cu}_2$, $\text{Co}_2\text{Cu}_2$.
  • Chất mang: Than hoạt tính (AC), magnesium oxide $\text{MgO}(200)$ và nhôm oxide $\text{Al}_2\text{O}_3(104)$.
  • Phạm vi hệ khảo sát: 7 hệ xúc tác cụ thể gồm $\text{NiCu/AC}$, $\text{Ni}_2\text{Cu}_2/\text{AC}$, $\text{Ni}_2\text{Cu}_2/\text{MgO}$, $\text{Co}_2\text{Cu}_2/\text{MgO}$, $\text{Co}_4/\text{Al}_2\text{O}_3$, $\text{Cu}_4/\text{Al}_2\text{O}_3$ và $\text{Co}_2\text{Cu}_2/\text{Al}_2\text{O}_3$.

Tổng quan tài liệu và vị trí của luận án

Các hướng nghiên cứu trước đây

Tác giả đã điểm lại lịch sử và tình hình nghiên cứu thực nghiệm cũng như lý thuyết về chuyển hóa syngas:

  1. Lịch sử phát hiện và xúc tác nhiệt độ cao/áp suất cao:

    • Frans Fischer và Hans Tropsch (1923): Lần đầu tiên phát hiện phản ứng tổng hợp alcohol mạch cao trực tiếp từ syngas qua quá trình Fischer-Tropsch (FT).
    • Natta và cộng sự: Nghiên cứu hệ xúc tác $\text{Zn/Cr}$ tổng hợp methanol và alcohol mạch cao ở $400 - 450^\circ\text{C}$ và $10 - 25{,}5\text{ MPa}$. Khi bổ sung kim loại kiềm ($\text{K}$, $\text{Rb}$, $\text{Cs}$), tốc độ và độ chọn lọc tạo isobutanol tăng rõ rệt; điều kiện công nghệ tối ưu được xác định tại $440^\circ\text{C}$, $10{,}3\text{ MPa}$, tỷ số mol $\text{H}_2/\text{CO} = 1{,}0$ và tốc độ không gian thể tích $\text{GHSV} = 12,000\text{ h}^{-1}$ trên xúc tác $\text{Zn/Cr}$ biến tính $1%\text{ K}$.

  2. Xúc tác nhiệt độ thấp và xúc tác Fischer-Tropsch biến tính:

    • Lewis và Frolich (1928): Phát hiện hệ xúc tác tổng hợp methanol ở áp suất và nhiệt độ thấp $\text{Cu/ZnO/Al}_2\text{O}_3$ (hoạt động ở $250 - 300^\circ\text{C}$, $2 - 10\text{ MPa}$). Khi có thêm $\text{K}$, xúc tác cho sản phẩm chọn lọc alcohol cao hơn $\text{C}_1$, tuy nhiên sản phẩm chính vẫn là methanol.
    • Zhigang và cộng sự: Chỉ ra rằng biến tính hệ xúc tác FT bằng $\text{Ni}$, $\text{Fe}$ hoặc $\text{Co}$ làm tăng độ chọn lọc alcohol mạch thẳng bậc cao ở điều kiện $4 - 10\text{ MPa}$, $200 - 350^\circ\text{C}$, tỷ số $\text{H}_2/\text{CO}$ từ $1$ đến $2$.
    • Các hệ xúc tác gốc Rhodium: Hoạt động ở điều kiện ôn hòa ($150 - 250^\circ\text{C}$, $0{,}1 - 2{,}5\text{ MPa}$). Khi thêm các oxide $\text{MnO}_x$, $\text{TiO}_2$, $\text{ZrO}_2$, độ chuyển hóa CO tăng từ 10 đến 50 lần và tăng độ chọn lọc đối với các sản phẩm $\text{C}_2$ chứa oxygen. Biến tính bởi $\text{Fe}$ hoặc $\text{Mo}$ thúc đẩy hình thành sản phẩm mạch nhánh.
    • Các hệ xúc tác gốc Molybdenum ($\text{MoS}_2$, $\beta\text{-Mo}_2\text{C}$, $\text{MoP}$): Phát triển bởi Dow Chemical và Union Carbide, có ưu điểm bền với ngộ độc lưu huỳnh và có hoạt tính cao với phản ứng chuyển dịch khí than nước (water-gas shift). Biến tính bằng kim loại nhóm VIII ($\text{Ni}$, $\text{Co}$) hoặc oxide kiềm giúp tăng độ chọn lọc alcohol mạch thẳng bậc cao.

  3. Chất mang xúc tác:

    • Các oxide truyền thống ($\text{Al}_2\text{O}_3$, $\text{SiO}_2$, $\text{CeO}_2$, $\text{ZrO}_2$): Bề mặt acid của $\text{Al}_2\text{O}_3$ có thể cản trở sự tạo thành alcohol nếu không được biến tính trung hòa.
    • Than hoạt tính (AC) và ống nano carbon (CNTs): Bề mặt riêng lớn, tính chất trung tính, bền hóa học, tương tác vừa phải với pha hoạt tính và không có tính chọn lọc đối với hydrocarbon.

  4. Cơ chế phản ứng chuyển hóa syngas:

    • Cơ chế chèn $\text{CO}$ ($\text{CO}$ insertion): Phân ly $\text{CO}$, hydrogen hóa tạo các mảnh $\text{C}_1$ ($\text{CH}_x^$), chèn $\text{CO}^$ vào nhóm alkyl để tạo cấu trúc acyl ($\text{CH}_3\text{CO}^*$), sau đó hydrogen hóa tiếp tạo acetaldehyde và ethanol.
    • Cơ chế chèn $\text{CH}_2$ ($\text{CH}_2$ insertion): Ghép nối liên tiếp các nhóm methylene để tăng mạch hydrocarbon hoặc alcohol.
    • Cơ chế hydroxycarbene: Chèn $\text{CH}_2$ vào phối tử $\text{HCOH}^$ hấp phụ tạo thành alkylhydroxycarbene ($\text{CH}_3\text{-CH-OH}^$).

Khoảng trống luận án lựa chọn giải quyết

Mặc dù thực nghiệm đã ghi nhận vai trò của các chất trợ xúc tác và chất mang, cơ chế vi mô chi tiết về sự phối hợp giữa hai tâm kim loại khác nhau (như Ni-Cu, Co-Cu) và tương tác của chúng với các bề mặt chất mang (AC, MgO, $\text{Al}_2\text{O}_3$) vẫn chưa được làm sáng tỏ đầy đủ. Luận án giải quyết khoảng trống này bằng cách mô phỏng lượng tử toàn diện các bước hấp phụ, phân ly, chèn mạch và hydrogen hóa, xác định các trạng thái chuyển tiếp và rào cản năng lượng của từng giai đoạn phản ứng.

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

Khung lý thuyết và mô hình hóa học lượng tử

Luận án áp dụng nền tảng Cơ học lượng tử và Hóa học tính toán với các mô hình và xấp xỉ cơ bản:

  • Phương trình Schrödinger dừng: $\hat{H}\Psi = E\Psi$, kết hợp xấp xỉ Born-Oppenheimer nhằm tách rời chuyển động của electron khỏi hạt nhân đứng yên.
  • Hàm sóng nhiều electron: Sử dụng định thức Slater nhằm thỏa mãn nguyên lý loại trừ Pauli và tính phản đối xứng của hàm sóng electron.
  • Thuyết phiếm hàm mật độ (DFT): Dựa trên hai định lý Hohenberg-Kohn:
    • Định lý thứ nhất: Thế năng ngoài $v(\mathbf{r})$ và năng lượng trạng thái cơ bản được xác định duy nhất bởi mật độ electron $\rho(\mathbf{r})$.
    • Định lý thứ hai: Năng lượng toàn phần $E[\rho]$ thỏa mãn nguyên lý biến phân, đạt cực tiểu tại mật độ electron trạng thái cơ bản $E[\rho] \ge E_0$.
  • Hệ phương trình Kohn-Sham: Chuyển bài toán tương tác nhiều hạt về hệ các phương trình tự hợp một hạt không tương tác, biểu diễn động năng $T_s[\rho]$, năng lượng đẩy tĩnh điện cổ điển Coulomb $J[\rho]$, và nhóm năng lượng trao đổi - tương quan $E_{xc}[\rho]$.
  • Phiếm hàm trao đổi - tương quan: Sử dụng xấp xỉ gradient tổng quát (GGA) với phiếm hàm Perdew-Burke-Ernzerhof (PBE), kết hợp đối chiếu với các phiếm hàm PW91, BLYP, B3LYP và xấp xỉ mật độ địa phương (LDA/HFS).
  • Bộ hàm cơ sở: Áp dụng tổ hợp tuyến tính các orbital nguyên tử (LCAO-MO) với các bộ hàm Gauss (GTO, PGTO, CGTO) và bộ hàm Slater (STO).

Phương pháp xác định trạng thái chuyển tiếp và đường phản ứng

  • Bề mặt thế năng (PES): Khảo sát các điểm cực tiểu địa phương (ứng với chất phản ứng RA, chất trung gian IS, sản phẩm PR) và các điểm yên ngựa bậc một (ứng với trạng thái chuyển tiếp TS, nơi ma trận đạo hàm bậc hai Hessian có duy nhất một trị riêng âm).
  • Phương pháp Climbing Image Nudged Elastic Band (CI-NEB): Dùng chuỗi hình ảnh liên kết đàn hồi (mặc định 7 ảnh trung gian) giữa cấu trúc đầu và cuối để tìm đường phản ứng có năng lượng tối thiểu (MEP) và định vị chính xác năng lượng của điểm yên ngựa mà không cần tính ma trận đạo hàm bậc hai toàn phần cho mọi cấu trúc lớn.
  • Tọa độ phản ứng nội tại (IRC): Kiểm tra và xác nhận cấu trúc TS kết nối chính xác chất phản ứng với sản phẩm.

Các đại lượng và công cụ phân tích dữ liệu

Ký hiệu đại lượng Tên gọi và ý nghĩa Đơn vị / Phương pháp
$E_{ads}$ Năng lượng hấp phụ phân tử trên bề mặt $\text{kJ/mol}$
$E_{rel}$ Năng lượng tương đối dọc theo đường phản ứng $\text{kJ/mol}$
$E_a$ Năng lượng hoạt hóa của từng giai đoạn phản ứng $\text{kJ/mol}$
$\Delta E$ Biến thiên năng lượng phản ứng (nhiệt phản ứng) $\text{kJ/mol}$
$d$ Độ dài liên kết giữa các nguyên tử Å (Angstrom)
$B$ Bậc liên kết hóa học Chỉ số bậc liên kết Mayer
$Q$ Điện tích nguyên tử từng tâm Phân tích điện tích Hirshfeld
$N_{ue}$ Số electron độc thân trong hệ cluster Electron
$\text{DOS}$ Mật độ trạng thái điện tử (spin $\alpha$, spin $\beta$) Trạng thái / eV
$\text{HOMO / LUMO}$ Orbital phân tử bị chiếm cao nhất / còn trống thấp nhất Biểu diễn đẳng giá trị mật độ electron

Nội dung chính theo từng chương

Chương 1. Cơ sở lý thuyết

Chương 1 (độ dài 22 trang) trình bày toàn diện các nguyên lý toán học và vật lý của hóa học lượng tử ứng dụng:

  • Thiết lập phương trình Schrödinger, cấu trúc toán tử Hamilton toàn phần gồm động năng hạt nhân, động năng electron, thế tương tác hút electron-hạt nhân, tương tác đẩy electron-electron và tương tác đẩy giữa các hạt nhân.
  • Trình bày chi tiết xấp xỉ Born-Oppenheimer, cấu trúc tích Hartree và định thức Slater cho hệ vỏ kín và vỏ mở.
  • Tổng quan các phương pháp tính toán: từ bán kinh nghiệm (Hückel mở rộng, CNDO, INDO, MINDO, NDDO), ab-initio (HF, MPn, CI, CC, CCSD(T)) đến lý thuyết phiếm hàm mật độ (DFT).
  • Phân tích mô hình Thomas-Fermi và hạn chế khi bỏ qua tương quan trao đổi; dẫn dắt hai định lý nền tảng Hohenberg-Kohn và phương pháp Kohn-Sham; chi tiết hóa các phiếm hàm LDA (phương pháp Thomas-Fermi-Dirac, Slater HFS) và GGA (PBE, PW91).
  • Trình bày lý thuyết bộ hàm cơ sở hóa trị, tối thiểu, mở rộng; lý thuyết bề mặt thế năng PES, ma trận Hessian, điểm cực tiểu và điểm yên ngựa bậc một.
  • Phân tích nguyên lý thuật toán tìm kiếm trạng thái chuyển tiếp: phương pháp dải đàn hồi NEB, kỹ thuật leo dải CI-NEB, phương pháp CPR, Drag, DHS, Dimer và phương pháp phát triển chuỗi GSM.

Chương 2. Tổng quan hệ chất nghiên cứu

Chương 2 (độ dài 10 trang) hệ thống hóa các tư liệu thực nghiệm và lý thuyết về phản ứng chuyển hóa syngas:

  • Phân loại 4 nhóm xúc tác chính cho quá trình tổng hợp methanol và alcohol mạch cao: xúc tác biến tính methanol (Zn/Cr, $\text{Cu/ZnO/Al}_2\text{O}_3$), xúc tác Fischer-Tropsch biến tính (Cu-Co, Cu-Fe, Cu-Ni), xúc tác gốc Rhodium và xúc tác gốc Molybdenum ($\text{MoS}_2$, $\beta\text{-Mo}_2\text{C}$).
  • Đánh giá đặc tính của các chất mang xúc tác: than hoạt tính (AC) có diện tích bề mặt lớn, tính trung tính, hạn chế tạo hydrocarbon nhưng có vi mao quản dễ tạo coke; ống nano carbon (CNTs) có cấu trúc mao quản trung bình đồng nhất; các oxide kim loại $\text{MgO}$ (tính base) và $\text{Al}_2\text{O}_3$ (tính acid).
  • Tổng hợp mạng lưới 10 phương trình phản ứng nhiệt động học của quá trình chuyển hóa syngas tạo methanol, phản ứng chuyển dịch khí than nước, phản ứng tạo ethanol, alcohol cao, ketone, aldehyde, ester và phản ứng phụ tạo methane ($\text{CO} + 3\text{H}_2 \leftrightarrow \text{CH}_4 + \text{H}_2\text{O}$).
  • Phân tích chi tiết 3 cơ chế phản ứng chính: cơ chế chèn $\text{CO}$, cơ chế chèn $\text{CH}_2$, cơ chế hydroxycarbene và vai trò của các tiểu phân trung gian $\text{C}_1$ ($\text{CHO}$, $\text{HCOH}$, $\text{CH}_2\text{OH}$, $\text{CH}_3$, $\text{CH}_2$, $\text{CH}$, $\text{COH}$).

Chương 3. Kết quả nghiên cứu và thảo luận

Chương 3 (độ dài 98 trang, gồm 29 bảng số liệu và 31 hình vẽ minh họa) là chương nội dung thực nghiệm tính toán chính của luận án, được chia thành các phần nghiên cứu cụ thể:

  1. Phản ứng hydrogen hóa CO trên hệ xúc tác Ni-Cu và $\text{Ni}_2\text{Cu}_2$ trên chất mang than hoạt tính (AC):

    • Khảo sát các cấu trúc hấp phụ của cluster $\text{NiCu}$ và $\text{Ni}2\text{Cu}2$ trên bề mặt graphene/than hoạt tính (AC). Tính toán năng lượng tối ưu hóa ($E{opt}$), năng lượng liên kết trung bình ($E_b$), số electron độc thân ($N{ue}$), và mật độ trạng thái DOS phân lớp $3d, 4s$ của Ni và Cu.
    • Khảo sát quá trình hấp phụ phân tử và hấp phụ phân ly của $\text{H}_2$ và $\text{CO}$ trên $\text{NiCu/AC}$ và $\text{Ni}_2\text{Cu}2/\text{AC}$. Phân tích mật độ trạng thái DOS spin $\alpha$ cho các orbital $2p$ của C, O và $1s$ của H; đánh giá sự suy giảm năng lượng liên kết $E{lk}(\text{C-O})$ khi có mặt xúc tác so với phân tử tự do.
    • Xây dựng bề mặt thế năng cho các đường phản ứng hydrogen hóa $\text{CO}$ tạo $\text{CH}_3\text{OH}$, $\text{CH}_4$, và đường chèn $\text{CO}^$ vào $\text{CH}_3^$ / $\text{CH}_2^$ tạo $\text{CH}_3\text{CO}^$, hydrogen hóa tiếp thành $\text{CH}_3\text{CHO}$ và $\text{CH}_3\text{CH}_2\text{OH}$.

  2. Phản ứng hydrogen hóa CO trên hệ xúc tác $\text{Ni}_2\text{Cu}_2$ trên chất mang magnesium oxide (MgO):

    • Xây dựng mô hình bề mặt tinh thể $\text{MgO}(200)$ dựa trên giản đồ XRD tính toán. Tối ưu hóa cấu trúc $\text{Ni}_2\text{Cu}_2/\text{MgO}$, xác định độ dài liên kết và điện tích tương tác tại bề mặt phân cách.
    • Khảo sát sự hấp phụ của $\text{CO}$ và $\text{H}_2$, phân tích hình ảnh orbital HOMO-LUMO của hệ $(\text{CO}^* + \text{H}_2^)/\text{Ni}_2\text{Cu}_2/\text{MgO}$ trước khi chuyển hóa $\text{CO}^$ thành $\text{CHO}^*$.
    • Tính toán rào cản năng lượng hoạt hóa ($E_a$) và biến thiên năng lượng ($\Delta E$) của các nhánh phản ứng: phân ly $\text{CH}_x\text{O}^$, tạo $\text{CH}_3\text{OH}$, $\text{CH}_4$ và chèn $\text{CO}$ tạo $\text{CH}_3\text{CO}^$ dẫn đến hình thành ethanol.

  3. Phản ứng hydrogen hóa CO trên hệ xúc tác $\text{Co}_2\text{Cu}_2$ trên chất mang MgO:

    • Khảo sát các cấu trúc hình học tối ưu của cluster $\text{Co}_2\text{Cu}_2$ trên $\text{MgO}(200)$, tính toán năng lượng hấp phụ $\text{H}_2$, $\text{CO}$ và phân tích DOS spin $\alpha$ của các phân lớp $2s, 2p$ của nguyên tử C và O trong cấu trúc hấp phụ $\text{CO}^*\text{-Co}_2\text{Cu}_2/\text{MgO}$.
    • Khảo sát các đường phản ứng chuyển hóa $\text{CO}$ qua các tâm Co-Cu tạo $\text{HCHO}$, $\text{CH}_3\text{OH}$, $\text{CH}_3^$, $\text{CH}_4$, $\text{CH}_3\text{CHO}$ và $\text{CH}_3\text{CH}_2\text{OH}$; phân tích DOS của nhóm $\text{CO}$ trong phối tử $\text{CH}_3\text{CO}^$.

  4. Phản ứng hydrogen hóa CO trên hệ xúc tác đơn kim loại $\text{Co}_4$ và $\text{Cu}_4$ trên chất mang $\text{Al}_2\text{O}_3$:

    • Xây dựng mô hình mặt tinh thể $\text{Al}2\text{O}3(104)$ từ dữ liệu XRD tính toán. Xác định thông số tối ưu $E{opt}$, $E_b$, $N{ue}$, tổng điện tích $Q$ Hirshfeld và tổng bậc liên kết Mayer giữa cụm kim loại $M_4$ với các nguyên tử O của bề mặt chất mang $\text{Al}_2\text{O}_3$.
    • Khảo sát độc lập quá trình hấp phụ $\text{H}_2$, $\text{CO}$ và tính toán các đường phản ứng hydrogen hóa toàn phần $\text{CO}$ trên $\text{Co}_4/\text{Al}_2\text{O}_3$ và $\text{Cu}_4/\text{Al}_2\text{O}_3$, chỉ rõ đặc tính ưu tiên tạo methane trên tâm Co thuần túy và ưu tiên tạo alcohol trên tâm Cu.

  5. Phản ứng hydrogen hóa CO trên hệ xúc tác lưỡng kim loại $\text{Co}_2\text{Cu}_2$ trên chất mang $\text{Al}_2\text{O}_3$ và so sánh tổng hợp:

    • Xác định cấu trúc tối ưu của $\text{Co}_2\text{Cu}_2/\text{Al}_2\text{O}_3$, các trạng thái hấp phụ $\text{H}_2$, $\text{CO}$, sơ đồ chuyển hóa $\text{CO}^* \rightarrow \text{C}^, \text{CHO}^, \text{COH}^$, phân tích DOS spin $\alpha$ của hai nguyên tử $\text{H}^$ và phân lớp $2s, 2p$ của C, O.
    • Tính toán toàn diện năng lượng tương đối $E_{rel}$, năng lượng hoạt hóa $E_a$ và nhiệt phản ứng $\Delta E$ cho các giai đoạn chuyển hóa từ $\text{CO}^$ thành $\text{CH}_2^$, $\text{CH}_3^*$, $\text{CH}_3\text{OH}$, và quá trình ghép nối chèn $\text{CO}$ tạo $\text{CH}_3\text{CHO}$, $\text{CH}_3\text{CH}_2\text{OH}$.
    • So sánh tổng hợp 7 hệ xúc tác đã khảo sát, khẳng định vị trí lưỡng kim loại Co-Cu kết hợp chất mang $\text{Al}_2\text{O}_3$ mang lại hiệu quả vượt trội về nhiệt động và động học.
Hệ xúc tác khảo sát Chất mang Tâm hoạt tính chính Xu hướng phản ứng ưu tiên Sản phẩm chính ghi nhận
$\text{NiCu/AC}$ Than hoạt tính Ni, Cu, Ni-Cu Hydrogen hóa $\text{CO}$, chèn $\text{CO}^*$ vào alkyl $\text{CH}_4$, $\text{CH}_3\text{OH}$, $\text{C}_2\text{H}_5\text{OH}$
$\text{Ni}_2\text{Cu}_2/\text{AC}$ Than hoạt tính $\text{Ni}_2\text{Cu}_2$ Giảm $E_a$ chèn $\text{CO}$, ức chế methane hóa $\text{CH}_3\text{OH}$, $\text{CH}_3\text{CHO}$, $\text{C}_2\text{H}_5\text{OH}$
$\text{Ni}_2\text{Cu}_2/\text{MgO}$ $\text{MgO}(200)$ Tâm phân cách $\text{Ni-Cu/MgO}$ Hoạt hóa $\text{CO}$ mạnh, chuyển hóa $\text{CHO}^*$ $\text{CH}_3\text{OH}$, $\text{CH}_3\text{CO}^*$, $\text{C}_2\text{H}_5\text{OH}$
$\text{Co}_2\text{Cu}_2/\text{MgO}$ $\text{MgO}(200)$ Tâm phân cách $\text{Co-Cu/MgO}$ Hydrogen hóa tạo $\text{CH}_3^*$, chèn $\text{CO}$ tạo acyl $\text{CH}_3\text{CHO}$, $\text{C}_2\text{H}_5\text{OH}$
$\text{Co}_4/\text{Al}_2\text{O}_3$ $\text{Al}_2\text{O}_3(104)$ Tâm Co thuần Phân ly sâu liên kết $\text{C-O}$, methane hóa mạnh $\text{CH}_4$, hydrocarbon
$\text{Cu}_4/\text{Al}_2\text{O}_3$ $\text{Al}_2\text{O}_3(104)$ Tâm Cu thuần Hydrogen hóa không phân ly $\text{C-O}$ $\text{CH}_3\text{OH}$
$\text{Co}_2\text{Cu}_2/\text{Al}_2\text{O}_3$ $\text{Al}_2\text{O}_3(104)$ Tâm phối hợp $\text{Co-Cu/Al}_2\text{O}_3$ Đa số phản ứng có $E_a$ nhỏ, $\Delta E < 0$, ưu tiên ethanol $\text{CH}_3\text{CH}_2\text{OH}$ (chọn lọc cao nhất)

Kết quả và những đóng góp mới

Những đóng góp mới về mặt lý luận và khoa học

  1. Xây dựng bức tranh bề mặt thế năng hoàn chỉnh: Luận án đã mô phỏng chi tiết các giai đoạn hấp phụ, hoạt hóa $\text{CO}$ và $\text{H}_2$, cũng như cơ chế phản ứng hydrogen hóa $\text{CO}$ tạo thành nhiều sản phẩm khác nhau ($\text{CH}_3\text{OH}$, $\text{CH}_4$, $\text{HCHO}$, $\text{CH}_3\text{CHO}$, $\text{CH}_3\text{CH}_2\text{OH}$) trên 7 hệ xúc tác đơn và lưỡng kim loại: $\text{NiCu/AC}$, $\text{Ni}_2\text{Cu}_2/\text{AC}$, $\text{Ni}_2\text{Cu}_2/\text{MgO}$, $\text{Co}_2\text{Cu}_2/\text{MgO}$, $\text{Co}_4/\text{Al}_2\text{O}_3$, $\text{Cu}_4/\text{Al}_2\text{O}_3$ và $\text{Co}_2\text{Cu}_2/\text{Al}_2\text{O}_3$.
  2. Làm rõ hiệu ứng hiệp đồng lưỡng kim loại: Các kết quả tính toán chứng minh rằng sự kết hợp giữa hai tâm kim loại (Ni-Cu hoặc Co-Cu) có tác dụng trực tiếp làm giảm năng lượng hoạt hóa ($E_a$) của bước chèn $\text{CO}^$ và các phản ứng hydrogen hóa tạo gốc methyl ($\text{CH}_3^$), dẫn đến sự hình thành ưu tiên các sản phẩm $\text{C}_2$ chứa oxygen (đặc biệt là ethanol) như sản phẩm chính trên các vị trí này.
  3. Xác định các chất trung gian then chốt: Luận án đã làm rõ các tiểu phân trung gian mang tính quyết định đến độ chọn lọc tạo ethanol là $\text{CH}_3\text{O}^$, $\text{CH}_2\text{OH}^$, $\text{CH}_3^$ và $\text{CH}_3\text{CO}^$. Sự cạnh tranh giữa khả năng hydrogen hóa và phân ly của $\text{CH}_3\text{O}^$ và $\text{CH}_2\text{OH}^$ quyết định tỷ lệ giữa methanol và ethanol. Việc gia tăng diện tích tiếp xúc của các vị trí lưỡng kim loại giúp làm suy yếu lực hấp phụ $\text{CO}$ quá mức và ngăn chặn phản ứng methane hóa tạo $\text{CH}_4$.
  4. Xác định hệ xúc tác tối ưu: Luận án chỉ ra rằng hệ xúc tác $\text{Co}_2\text{Cu}_2/\text{Al}_2\text{O}_3$ là hệ tiềm năng và thuận lợi nhất cho quá trình tổng hợp ethanol từ syngas. Đa số các giai đoạn phản ứng trên hệ này đều có rào cản năng lượng hoạt hóa $E_a$ nhỏ và biến thiên năng lượng $\Delta E$ âm (tỏa nhiệt thuận lợi về mặt nhiệt động học), đồng thời chỉ rõ vai trò tương tác điện tử tích cực của chất mang $\text{Al}_2\text{O}_3$.

Ý nghĩa thực tiễn và kiến nghị giải pháp

  • Kết quả nghiên cứu cung cấp cơ sở dữ liệu lý thuyết đáng tin cậy để định hướng thiết kế và tổng hợp thực nghiệm các hệ xúc tác lưỡng kim loại mới có hiệu suất và độ chọn lọc cao cho quá trình sản xuất ethanol và alcohol mạch cao từ khí thải carbon và khí syngas.
  • Góp phần thúc đẩy phát triển công nghệ tái chế carbon, sản xuất nhiên liệu lỏng sạch thay thế nhiên liệu hóa thạch, hỗ trợ giải quyết đồng thời bài toán kinh tế năng lượng và giảm thiểu ô nhiễm môi trường.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của đề tài

  • Nghiên cứu được thực hiện dựa trên các mô hình cụm kim loại hữu hạn (cluster từ 2 đến 4 nguyên tử: $\text{NiCu}$, $\text{Ni}_2\text{Cu}_2$, $\text{Co}_2\text{Cu}_2$, $\text{Co}_4$, $\text{Cu}_4$) trên các phiến cắt bề mặt tinh thể lý tưởng ($\text{MgO}(200)$, $\text{Al}_2\text{O}_3(104)$ và mô hình than hoạt tính đơn lớp), chưa phản ánh toàn bộ tính đa dạng cấu trúc bề mặt thực tế (như các bậc thang, khuyết tật mạng phức tạp hay kích thước hạt nano thực tế).
  • Các tính toán hóa học lượng tử DFT khảo sát phản ứng chủ yếu ở điều kiện không tuyệt đối ($0\text{ K}$) hoặc thông qua hiệu chỉnh nhiệt động cơ bản, chưa xét hết các biến số động học phức tạp trong môi trường dòng chảy công nghiệp ở áp suất và nhiệt độ cao.

Kiến nghị các nghiên cứu tiếp theo

  • Mở rộng tính toán mô phỏng trên các cluster kim loại có kích thước lớn hơn hoặc mô hình định kỳ bề mặt (periodic slab models) với tỷ lệ hợp kim Ni-Cu và Co-Cu khác nhau.
  • Khảo sát sâu hơn ảnh hưởng của các chất trợ xúc tác kim loại kiềm (K, Cs, Rb) khi bổ sung vào hệ lưỡng kim loại mang trên oxide.
  • Ứng dụng phương pháp động lực học phân tử lượng tử (ab-initio Molecular Dynamics - AIMD) kết hợp mô phỏng Monte Carlo để khảo sát động thái phản ứng tại điều kiện nhiệt độ và áp suất vận hành thực tế.
  • Tiến hành thực nghiệm tổng hợp và đánh giá hoạt tính xúc tác trên hệ $\text{Co-Cu/Al}_2\text{O}_3$ và $\text{Ni-Cu/AC}$ dựa trên các thông số cấu trúc tối ưu mà luận án đã đề xuất.

Giá trị tham khảo

Luận án là tài liệu chuyên khảo có giá trị trong các lĩnh vực:

  1. Nghiên cứu sinh, học viên cao học và giảng viên ngành Hóa lý thuyết và Hóa lý: Cung cấp quy trình phương pháp luận chuẩn xác trong việc ứng dụng DFT và thuật toán CI-NEB để xây dựng bề mặt thế năng, tối ưu cấu trúc bề mặt xúc tác dị thể, xác định trạng thái chuyển tiếp và tính toán các đại lượng hóa học lượng tử ($E_{ads}$, $E_a$, DOS, Hirshfeld charge, Mayer bond order).
  2. Các nhà nghiên cứu trong lĩnh vực Xúc tác - Hấp phụ và Công nghệ Hóa học: Cung cấp bộ thông số cấu trúc, năng lượng liên kết và cơ chế phản ứng chi tiết của 7 hệ xúc tác kim loại chuyển tiếp trên các chất mang $\text{AC}$, $\text{MgO}$, $\text{Al}_2\text{O}_3$, hỗ trợ trực tiếp cho việc lựa chọn thành phần vật liệu xúc tác.
  3. Các nhà phát triển công nghệ năng lượng và vật liệu mới: Cung cấp luận cứ khoa học về khả năng thay thế xúc tác kim loại quý (Rhodium) bằng xúc tác lưỡng kim loại chuyển tiếp chi phí thấp ($\text{Co-Cu}$, $\text{Ni-Cu}$) cho quá trình tổng hợp nhiên liệu lỏng sạch từ khí syngas.

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

1. Tại sao hệ xúc tác lưỡng kim loại $\text{Co}_2\text{Cu}_2/\text{Al}_2\text{O}_3$ lại được đánh giá là tiềm năng nhất cho phản ứng tổng hợp ethanol?

Tính toán hóa học lượng tử DFT cho thấy trên hệ $\text{Co}_2\text{Cu}_2/\text{Al}_2\text{O}_3$, đa số các bước phản ứng chuyển hóa $\text{CO}$ và $\text{H}_2$ tạo các mảnh trung gian $\text{CH}_2^$, $\text{CH}_3^$ và bước chèn $\text{CO}^$ tạo nhóm acyl $\text{CH}_3\text{CO}^$ đều có năng lượng hoạt hóa ($E_a$) thấp và biến thiên năng lượng ($\Delta E$) mang giá trị âm (tỏa nhiệt). Sự phối hợp giữa tâm Co (có khả năng hoạt hóa, phân ly $\text{CO}$ và hydrogen hóa) và tâm Cu (thúc đẩy chèn $\text{CO}$ và tạo alcohol) cùng tương tác phân bố điện tích với chất mang $\text{Al}_2\text{O}_3$ giúp tối ưu hóa rào cản động học và nhiệt động học theo hướng tạo ethanol thay vì methane.

2. Các chất trung gian phản ứng nào đóng vai trò quyết định độ chọn lọc ethanol trong luận án?

Luận án xác định được 4 chất trung gian then chốt quyết định độ chọn lọc ethanol là $\text{CH}_3\text{O}^$, $\text{CH}_2\text{OH}^$, $\text{CH}_3^$ và $\text{CH}_3\text{CO}^$. Trong đó, sự cạnh tranh giữa khả năng hydrogen hóa tiếp tục hay phân ly của $\text{CH}_3\text{O}^$ và $\text{CH}_2\text{OH}^$ quyết định độ chọn lọc giữa methanol và các sản phẩm mạch cao, còn sự hình thành bền vững của nhóm acyl $\text{CH}_3\text{CO}^$ qua bước chèn $\text{CO}^$ là tiền đề trực tiếp để hydrogen hóa tạo thành acetaldehyde và ethanol.

3. Sự khác biệt cơ bản về hướng tạo sản phẩm giữa xúc tác đơn kim loại Co, Cu và xúc tác lưỡng kim loại Co-Cu là gì?

Trên xúc tác đơn kim loại $\text{Co}_4/\text{Al}_2\text{O}_3$, tâm Co có hoạt tính phân ly liên kết $\text{C-O}$ rất mạnh, dẫn đến phản ứng hydrogen hóa sâu tạo thành methane ($\text{CH}_4$) là sản phẩm chủ yếu. Ngược lại, trên xúc tác đơn kim loại $\text{Cu}_4/\text{Al}_2\text{O}_3$, tâm Cu không phân ly mạnh liên kết $\text{C-O}$ mà chủ yếu hydrogen hóa không phân ly tạo thành methanol ($\text{CH}_3\text{OH}$). Khi kết hợp thành hệ lưỡng kim loại $\text{Co}_2\text{Cu}_2$, bề mặt phân cách giữa hai kim loại làm giảm năng lượng hoạt hóa của phản ứng chèn $\text{CO}^$ vào nhóm methyl $\text{CH}_3^$, ức chế quá trình methane hóa và chuyển hướng ưu tiên chọn lọc tạo ethanol.

4. Phương pháp CI-NEB có ưu thế gì trong việc khảo sát cơ chế phản ứng trên bề mặt xúc tác trong luận án?

Phương pháp Climbing Image Nudged Elastic Band (CI-NEB) cho phép xác định chính xác cấu trúc và năng lượng của trạng thái chuyển tiếp (TS) dọc theo đường phản ứng có năng lượng tối thiểu (MEP) chỉ dựa trên cấu trúc ban đầu đã biết (chất phản ứng) và cấu trúc cuối (sản phẩm). Phương pháp này giúp điểm ảnh có năng lượng cao nhất leo lên đỉnh bề mặt thế năng, tránh việc phải tính ma trận đạo hàm bậc hai Hessian toàn phần vô cùng tốn kém tài nguyên tính toán đối với các hệ mô hình bề mặt lớn chứa hàng trăm nguyên tử như $\text{M}_4/\text{MgO}$ hay $\text{M}_4/\text{Al}_2\text{O}_3$.

5. Chất mang than hoạt tính (AC) thể hiện đặc tính gì khi phân tán các cụm kim loại Ni-Cu?

Than hoạt tính (AC) có diện tích bề mặt riêng lớn, tính chất hóa học trung tính và không có tính chọn lọc đối với hydrocarbon. Kết quả tính toán cho thấy các cluster $\text{NiCu}$ và $\text{Ni}_2\text{Cu}_2$ hấp phụ ổn định trên bề mặt AC thông qua tương tác orbital, làm giảm năng lượng liên kết $\text{C-O}$ của phân tử $\text{CO}$ bị hấp phụ và tạo điều kiện thuận lợi cho phản ứng chèn $\text{CO}^*$ tạo các hợp chất $\text{C}_2$ chứa oxygen.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Bình Long đã ứng dụng phương pháp phiếm hàm mật độ (DFT) kết hợp kỹ thuật CI-NEB để làm sáng tỏ cơ chế phản ứng hydrogen hóa $\text{CO}$ trên 7 hệ xúc tác đơn và lưỡng kim loại Ni, Cu, Co phân tán trên các chất mang $\text{AC}$, $\text{MgO}$ và $\text{Al}_2\text{O}_3$. Công trình chứng minh rõ hiệu ứng hiệp đồng lưỡng kim loại trong việc giảm năng lượng hoạt hóa của bước chèn $\text{CO}^*$, ức chế quá trình methane hóa và nâng cao độ chọn lọc tổng hợp ethanol. Trong các hệ khảo sát, $\text{Co}_2\text{Cu}_2/\text{Al}_2\text{O}_3$ được xác định là hệ xúc tác tiềm năng và thuận lợi nhất về cả phương diện động học lẫn nhiệt động học, cung cấp cơ sở khoa học định lượng hỗ trợ thiết kế các vật liệu xúc tác chuyển hóa syngas trong thực tiễn.