Luận án: Hydrogen hóa CO2 xúc tác Ni5/MgO, AC bằng DFT - Văn Thị Minh Huệ

Tổng quan về luận án

Sự gia tăng nồng độ khí nhà kính trong khí quyển, đặc biệt là carbon dioxide ($CO_2$) phát sinh từ các hoạt động công nghiệp và sử dụng nhiên liệu hóa thạch, đang đặt ra thách thức sinh thái và khí hậu toàn cầu. Trước bối cảnh cạn kiệt nguồn tài nguyên hóa thạch, chiến lược thu hồi, cố định và tái sử dụng $CO_2$ nhằm chuyển hóa thành nguồn năng lượng mới và nguyên liệu hóa học có giá trị cao đóng vai trò then chốt. Trong số các lộ trình chuyển hóa, phản ứng hydrogen hóa $CO_2$ tạo methane ($CH_4$)—được biết đến là phản ứng Sabatier ($CO_2 + 4H_2 \rightarrow CH_4 + 2H_2O$)—giữ vị trí đặc biệt quan trọng nhờ khả năng tích hợp vào công nghệ "Power-to-Gas", biến đổi điện năng tái tạo dư thừa thành năng lượng hóa học dự trữ, điển hình như dự án thương mại công suất 1,0 MW vận hành tại Copenhagen từ năm 2016.

Mặc dù các kim loại quý như Ru, Rh, Pd thể hiện hoạt tính và độ chọn lọc cao, rào cản chi phí lớn khiến xúc tác kim loại chuyển tiếp trên cơ sở nickel (Ni) trở thành lựa chọn tối ưu về mặt kinh tế kỹ thuật. Tuy nhiên, phân tử $CO_2$ sở hữu độ bền nhiệt động rất cao với năng lượng liên kết $C=O$ lên tới $1680,8\text{ kJ}\cdot\text{mol}^{-1}$ và độ dài liên kết $1,17\text{ \AA}$, gây khó khăn lớn cho quá trình hoạt hóa phân tử. Hơn nữa, xúc tác Ni thường đối mặt với nguy cơ bị oxy hóa và tạo phức carbonyl độc hại $Ni(CO)_4/Ni(CO)_5$ gây mất hoạt tính. Khoảng trống nghiên cứu (research gap) cốt lõi tồn tại trong y văn học thuật quốc tế hiện nay là: cơ chế phản ứng ở cấp độ phân tử trên các cụm nguyên tử kích thước nanomet (nanocluster) vẫn chưa sáng tỏ; vai trò tương tác điện tử giữa cụm kim loại và chất mang (oxide bazơ như $MgO$ hay vật liệu carbon xốp $AC$) trong việc điều biến đường phản ứng ưu tiên (qua trung gian $CO$ hay qua trung gian formate $HCOO$) vẫn còn nhiều mâu thuẫn giữa các công trình thực nghiệm và tính toán.

Để giải quyết triệt để khoảng trống này, luận án của nghiên cứu sinh Văn Thị Minh Huệ (chuyên ngành Hóa lí thuyết và Hóa lí, Trường Đại học Sư phạm Hà Nội, dưới sự hướng dẫn của PGS. Nguyễn Ngọc Hà) đặt ra 4 câu hỏi nghiên cứu và các giả thuyết khoa học tương ứng:

• Câu hỏi nghiên cứu 1 (RQ1): Cấu trúc hình học, độ bền nhiệt động và trạng thái spin cơ bản của cluster $Ni_5$ tự do biến đổi như thế nào dưới tác động của hiệu ứng lượng tử?
  Giả thuyết 1 (H1): Cluster $Ni_5$ tồn tại ở dạng cấu trúc 3 chiều lưỡng tháp tam giác bị biến dạng phi đối xứng theo hiệu ứng Jahn-Teller nhằm giải tỏa sự suy biến năng lượng orbital.
• Câu hỏi nghiên cứu 2 (RQ2): Bản chất của quá trình hấp phụ và hoạt hóa $H_2$, $CO_2$ trên các tâm kim loại $Ni_5$ tự do và $Ni_5$ gắn trên chất mang $MgO(200)$, $AC$ là hấp phụ vật lý hay hấp phụ hóa học?
  Giả thuyết 2 (H2): Phân tử $H_2$ bị hấp phụ phân ly tự phát, trong khi $CO_2$ nhận điện tích từ orbital d của Ni vào orbital phản liên kết $2\pi^*$ làm suy yếu và hoạt hóa mạnh liên kết $C=O$.
• Câu hỏi nghiên cứu 3 (RQ3): Giữa các sản phẩm chứa 1 carbon ($CH_4, CH_3OH, HCHO, HCOOH$), đường phản ứng nào là ưu tiên nhất về mặt nhiệt động học và động học?
  Giả thuyết 3 (H3): Phản ứng ưu tiên tạo thành $CH_4$ thông qua con đường phân cắt liên kết tạo chất trung gian $CO$ (CO route) thay vì con đường qua trung gian formate.
• Câu hỏi nghiên cứu 4 (RQ4): Chất mang $MgO$ và $AC$ đóng vai trò như thế nào trong việc thay đổi cấu trúc electron của cluster kim loại và điều chỉnh hàng rào năng lượng hoạt hóa ($E_a$)?
  Giả thuyết 4 (H4): Chất mang không chỉ phân tán tâm hoạt động mà còn gây dịch chuyển mật độ điện tích bề mặt, làm biến đổi thế năng của các trạng thái chuyển tiếp.

Khung lý thuyết của nghiên cứu được xây dựng vững chắc trên nền tảng Lý thuyết Phiếm hàm Mật độ (Density Functional Theory - DFT), sử dụng phương pháp dải đàn hồi cải tiến CI-NEB (Climbing Image Nudged Elastic Band) để định vị trạng thái chuyển tiếp (Transition State - TS). Phạm vi nghiên cứu bao quát 03 hệ xúc tác gồm $Ni_5$, $Ni_5/MgO(200)$ và $Ni_5/AC$ với quy mô hệ mô phỏng lên tới 95 nguyên tử trong hộp kích thước $25\times 25\times 25\text{ \AA}^3$. Toàn bộ luận án gồm 150 trang, 28 bảng số liệu, 42 hình vẽ đồ thị, 130 tài liệu tham khảo, khảo sát toàn diện 33 đường phản ứng với 93 bước cơ bản. Các kết quả đột phá của luận án đã được công bố trong 06 công trình khoa học chuyên ngành, bao gồm 01 bài báo trên tạp chí quốc tế thuộc danh mục ISI, 01 báo cáo tại hội thảo quốc tế và 04 công trình tại các tạp chí/hội nghị chuyên ngành quốc gia.

Literature Review và Positioning

Quá trình chuyển hóa $CO_2$ thành nhiên liệu và hóa chất hữu cơ đã hình thành hai luồng quan điểm lý thuyết và thực nghiệm lớn trong cộng đồng khoa học quốc tế:

Luồng nghiên cứu thứ nhất tập trung vào quá trình tổng hợp methanol ($CH_3OH$) và cơ chế qua chất trung gian formate ($HCOO$). Tiêu biểu là công trình của Grabow và Mavrikakis (2011) nghiên cứu phản ứng trên bề mặt xúc tác $Cu(111)$ bằng tính toán DFT-PW91 kết hợp mô hình vi động học (microkinetic), đề xuất cơ chế thuận lợi qua chuỗi: $CO_2 \rightarrow HCOO \rightarrow HCOOH \rightarrow CH_3O_2 \rightarrow CH_2O \rightarrow CH_3O \rightarrow CH_3OH$. Đồng thuận với quan điểm này, Zhao và cộng sự (2011) trên hệ $Cu(111)$ xác định hàng rào năng lượng hình thành cấu trúc trung gian $HCOO$ ($E_a = 0,6\text{ eV}$) thấp hơn đáng kể so với cấu trúc $COOH$ ($E_a = 1,17\text{ eV}$). Về mặt thực nghiệm, Gonzalez và các cộng sự (2017) sử dụng kỹ thuật phổ DRIFTS chứng minh $CO_2$ hấp phụ tạo bicarbonate và formate trên xúc tác có nhóm hydroxyl bề mặt trước khi hydrogen hóa tạo methane.

Luồng nghiên cứu thứ hai khẳng định cơ chế phản ứng đi qua chất trung gian carbon monoxide ($CO$) và sự phân cắt trực tiếp liên kết $C-O$. Nghiên cứu của Ren và cộng sự (2015) trên bề mặt tinh thể $Ni(111)$ bằng DFT chỉ ra đường phản ứng ưu tiên nhất là $CO_2 \rightarrow CO + O \rightarrow C + O + 4H \rightarrow CH_2 + 2H \rightarrow CH_3 + H \rightarrow CH_4$. Tương tự, nghiên cứu lý thuyết bán kinh nghiệm ASED-MO của Cnudde và cộng sự cũng như nghiên cứu thực nghiệm AP-XPS của Weatherup và cộng sự (2014) xác nhận sự phân hủy $CO^$ ($E_a = 2,97\text{ eV}$) hoặc sự kết hợp $2CO^ \rightarrow C^* + CO_{2(g)}$ ($E_a = 1,93\text{ eV}$) là giai đoạn mấu chốt. Bằng kỹ thuật phổ hồng ngoại operando (operando IR), Henriques và cộng sự trên xúc tác zeolite NiUSY đã chứng minh $CO$ là chất trung gian thực sự và phản ứng hydrogen hóa/phân ly $CO$ chính là giai đoạn quyết định tốc độ phản ứng (Rate-Determining Step - RDS).

Tranh luận học thuật càng trở nên phức tạp khi xem xét ảnh hưởng của cấu trúc bề mặt và môi trường chất mang. Vesselli và cộng sự (2013) khi kết hợp thực nghiệm TPD, XPS, HREELS với tính toán DFT trên $Ni(110)$ đã chứng minh sự hình thành $HCOO$ không thuận lợi bằng sự tạo thành trung gian $CO$. Trong khi đó, Pan và cộng sự (2014) khi khảo sát cluster $Ni_4/\gamma-Al_2O_3$ chỉ ra rằng trên chất mang khô, sự tạo thành $CO$ có $E_a = 2,13\text{ eV}$ (kém thuận lợi động học hơn $HCOO$ với $E_a = 1,25\text{ eV}$), nhưng trên chất mang ướt có nhóm -OH bề mặt, sự hình thành $CO$ lại chiếm ưu thế tuyệt đối ($E_a = 0,69\text{ eV}$ so với $2,32\text{ eV}$ của $HCOO$).

Vị trí của luận án được định vị tại ranh giới tiên phong của chuyên ngành: thay vì khảo sát các bề mặt tinh thể kim loại khối lý tưởng (như $Ni(111), Ni(110)$) vốn đã bão hòa trong y văn, luận án tập trung khảo sát hệ nanocluster kim loại kích thước siêu nhỏ ($Ni_5$) ở trạng thái tự do và trạng thái tương tác với hai loại chất mang có bản chất hóa học trái ngược—chất mang oxide kim loại kiềm thổ có tính bazơ ($MgO(200)$) và chất mang carbon vô định hình dẫn điện ($AC$). Luận án cung cấp bức tranh cơ chế định lượng đầy đủ đầu tiên với 33 đường phản ứng song song, giải quyết triệt để tranh cãi về con đường trung gian $CO$ so với $HCOO$ trên xúc tác cụm nguyên tử phân tán.

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

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

Luận án mang lại những đóng góp nền tảng cho lý thuyết hóa học lượng tử và xúc tác dị thể:

• Mở rộng lý thuyết liên kết kim loại chuyển tiếp và định lý Jahn-Teller (H. Jahn & E. Teller, 1937): Luận án chứng minh rằng trong cluster $Ni_5$, cấu hình không gian bền nhất là dạng lưỡng tháp tam giác không đối xứng ($Ni_5(8)$) với năng lượng liên kết $E_b = 3,60\text{ eV}$ và tổng bậc liên kết Mayer đạt 5,1, vượt trội hơn dạng tháp vuông ($Ni_5(7)$, $E_b = 3,38\text{ eV}$, $E_{rel} = 1,10\text{ eV}$). Cấu trúc này bị biến dạng hình học nhằm giải tỏa sự suy biến năng lượng orbital d, làm giảm số electron độc thân từ 10 (của 5 nguyên tử Ni tự do cấu hình $3d^8 4s^2$) xuống còn đúng 8 electron độc thân, phản ánh sự ghép đôi tạo liên kết cộng hóa trị $Ni-Ni$.
• Bổ sung lý thuyết chuyển dịch điện tích trong hoạt hóa phân tử nhỏ: Luận án làm sáng tỏ cơ chế cho-nhận electron theo mô hình Dewar-Chatt-Duncanson mở rộng, trong đó orbital d của cluster Ni đóng vai trò chất cho electron vào orbital phản liên kết LUMO $2\pi^*$ của phân tử $CO_2$. Sự chuyển dịch điện tích này làm giảm bậc liên kết $C=O$, kéo dãn độ dài liên kết và bẻ cong góc phân tử, tạo tiền đề động học thuận lợi cho bước cắt đứt liên kết tiếp theo.
• Xác lập mô hình tương tác điện tử kim loại - chất mang (Metal-Support Interaction - MSI): Nghiên cứu chứng minh quy luật biến đổi cấu trúc electron khi đưa cluster lên chất mang: mật độ điện tích trên cluster $Ni_5$ bị tái phân bố do tương tác với các tâm oxy của $MgO$ hoặc mạng carbon sp2 của $AC$, làm giảm năng lượng hấp phụ ban đầu nhưng đồng thời hạ thấp năng lượng hoạt hóa của các bước phân cắt liên kết trung gian.

graph TD
    A["Cluster Ni5 (Lưỡng tháp tam giác biến dạng Jahn-Teller, 8e- độc thân)"] --> B["Hấp phụ H2: Phân ly tự phát thành 2H*"]
    A --> C["Hấp phụ CO2: Chuyển dịch e- d(Ni) -> LUMO 2π*(CO2)"]
    B --> D{"Bề mặt thế năng phản ứng (PES)"}
    C --> D
    D -->|"CO Route (Ưu tiên động học & nhiệt động)"| E["Trung gian CO* + O*"]
    D -.->|"Formate Route (Hàng rào Ea cao)"| F["Trung gian HCOO*"]
    E --> G["Hydrogen hóa từng bước: C* -> CH* -> CH2* -> CH3* -> CH4"]
    G --> H["Khí Methane (CH4) - Độ chọn lọc tối đa"]
    style E stroke:#2ecc71,stroke-width:2px;
    style G stroke:#2ecc71,stroke-width:2px;
    style H stroke:#2ecc71,stroke-width:3px;
    style F stroke:#e74c3c,stroke-dasharray: 5 5;

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

Khung phân tích của luận án là sự tích hợp đa ngành giữa ba trụ cột:

1. Lý thuyết Phiếm hàm Mật độ (Kohn-Sham DFT): Mô tả trạng thái vi mô của hệ đa electron thông qua mật độ electron $\rho(\mathbf{r})$ phụ thuộc 3 tọa độ không gian, khắc phục sự bùng nổ biến số của lý thuyết hàm sóng truyền thống ($3N$ biến).
2. Lý thuyết Trạng thái Chuyển tiếp (Transition State Theory) và CI-NEB: Khảo sát bề mặt thế năng đa chiều (PES), tìm kiếm chính xác điểm yên ngựa bậc một (first-order saddle point) tương ứng với trạng thái chuyển tiếp của từng bước cơ bản.
3. Phương pháp phân tích quần cư electron hiện đại: Sử dụng phân tích điện tích theo biến dạng Voronoi (Voronoi Deformation Density) và bậc liên kết theo Mayer (Mayer Bond Order)—hai kỹ thuật có ưu điểm vượt trội là hầu như không phụ thuộc vào kích thước bộ hàm cơ sở, mang lại giá trị vật lý sát với thực nghiệm.

Điều kiện biên của khung phân tích được xác định rõ ràng: các tính toán được thực hiện tại nhiệt độ $0\text{ K}$ trong pha khí; các lớp nguyên tử chất mang phía dưới được cố định tọa độ (fixed) để bảo toàn cấu trúc mạng tinh thể chuẩn, trong khi toàn bộ cluster kim loại và các tiểu phân chất phản ứng được thả lỏng hoàn toàn (relaxed) theo mọi bậc tự do hình học.

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

Thiết kế nghiên cứu

Thiết kế nghiên cứu tuân theo hệ hình triết học thực chứng tính toán (Computational Positivism), tiếp cận thực tại khách quan thông qua các định luật cơ bản của cơ học lượng tử ab initio:

• Thiết kế đa cấp độ (Multi-level Design): Nghiên cứu được triển khai đồng thời ở 4 cấp độ: (i) Cấp độ hạt cơ bản và orbital (mật độ electron, cấu hình spin, mật độ trạng thái riêng PDOS, HOMO-LUMO); (ii) Cấp độ liên kết nguyên tử (độ dài liên kết, góc liên kết, bậc liên kết Mayer, điện tích Voronoi); (iii) Cấp độ tương tác pha dị thể (giao diện cluster $Ni_5$ - chất mang $MgO/AC$); (iv) Cấp độ động học hệ thống (bề mặt thế năng PES gồm 33 đường phản ứng, 93 trạng thái chuyển tiếp).
• Mô hình tính toán chuẩn xác: Mô hình cụm nguyên tử $Ni_5$; mô hình bề mặt $MgO(200)$ gồm các lớp mạng tinh thể lập phương tâm diện; mô hình than hoạt tính ($AC$) được xây dựng từ các vòng thơm đa nhân đại diện cho cấu trúc graphene xốp; hệ nghiên cứu lớn nhất đạt quy mô 95 nguyên tử.

Quy trình nghiên cứu rigorous

Quy trình tính toán lượng tử được chuẩn hóa nghiêm ngặt trên nền tảng phần mềm QUANTUM tích hợp mã nguồn mở SIESTA (Spanish Initiative for Electronic Simulations with Thousands of Atoms):

• Phiếm hàm tương quan trao đổi: Sử dụng phiếm hàm xấp xỉ gradient suy rộng GGA-PBE (Perdew-Burke-Ernzerhof), đáp ứng độ chính xác cao đối với các hệ chứa kim loại chuyển tiếp và biến thiên mật độ electron phức tạp.
• Bộ hàm cơ sở: Bộ hàm cơ sở số nguyên tử DZP (Double Zeta plus Polarization), bổ sung các hàm phân cực cho phép mô tả chính xác sự biến dạng của đám mây electron khi hình thành liên kết hóa học.
• Thế giả bảo toàn chuẩn: Sử dụng thế giả Kleinman-Bylander dạng Troullier-Martins cho các electron vùng lõi, với ngưỡng cắt hàm sóng tương đương sóng phẳng (mesh cutoff) cực cao đạt $2040,75\text{ eV}$ (~$150\text{ Ry}$).
• Kích thước hộp mô phỏng: Đặt hệ trong hộp siêu mạng tinh thể định kỳ $25\times 25\times 25\text{ \AA}^3$, khoảng không gian chân không đủ lớn để triệt tiêu hoàn toàn tương tác giả tạo giữa các tế bào lân cận.
• Tiêu chuẩn hội tụ: Tối ưu hóa cấu trúc bằng thuật toán Quasi-Newton với tiêu chuẩn hội tụ lực cực đại nhỏ hơn $0,05\text{ eV/\AA}$.
• Xác định trạng thái chuyển tiếp CI-NEB: Sử dụng chuỗi 7 ảnh (images) nội suy tuyến tính dọc theo đường cực tiểu năng lượng (MEP), tiêu chuẩn hội tụ lực trên mỗi ảnh đạt dưới $0,1\text{ eV/\AA}$.
• Trạng thái spin: Kích hoạt tùy chọn phân cực spin (Spin Polarized = True) trong toàn bộ các phép tính, tự động tối ưu hóa số electron độc thân của hệ tại trạng thái năng lượng thấp nhất.

Data và phân tích

• Năng lượng hấp phụ: Xác định theo biểu thức nhiệt động $E_{ads} = E_{AB} - E_A - E_B$.
• Bậc liên kết Mayer ($B_{AB}$): Tính toán trực tiếp từ các ma trận mật độ spin $\alpha$ và spin $\beta$ theo công thức: $$B_{AB} = 2 \sum_{\mu \in A} \sum_{\nu \in B} \left[ (P^\alpha S){\mu\nu} (P^\alpha S){\nu\mu} + (P^\beta S){\mu\nu} (P^\beta S){\nu\mu} \right]$$
• Kiểm chứng độ tin cậy của phương pháp (Robustness & Benchmarking): Luận án tính toán năng lượng ion hóa của cluster $Ni_5$ dạng lưỡng tháp tam giác đạt $IP(Ni_5) = 5,91\text{ eV}$. So sánh với giá trị đo đạc thực nghiệm là $6,17\text{ eV}$, sai số tương đối chỉ là $4,2%$. Độ dài liên kết $Ni-Ni$ tính toán đạt $2,270 - 2,614\text{ \AA}$ (trung bình $2,40\text{ \AA}$), tương thích cao với dữ liệu thực nghiệm và lý thuyết quốc tế của Reddy ($2,27\text{ \AA}$) và Castro ($2,34\text{ \AA}$).

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

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

Luận án mang lại 5 phát hiện mang tính đột phá được chứng minh bằng các dữ liệu định lượng cụ thể:

1. Sự ổn định cấu trúc và tính chất từ của nanocluster $Ni_5$: "Các cấu trúc phẳng đều không bền, do có năng lượng cao hơn đáng kể so với hai đồng phân $Ni_5(7)$ và $Ni_5(8)$... cấu trúc lưỡng tháp tam giác được xác định là bền nhất." Đồng phân $Ni_5(8)$ có $E_{rel} = 0,0\text{ eV}$, $E_b = 3,60\text{ eV}$, tổng bậc liên kết Mayer đạt 5,1 và 8 electron độc thân, cấu trúc không đối xứng với khoảng cách $Ni-Ni$ ngắn nhất là $2,270\text{ \AA}$. Sự phân bố điện tích Voronoi không đồng nhất trên các nguyên tử $Ni_1 \rightarrow Ni_5$ lần lượt là $+0,001; +0,022; +0,011; +0,006; +0,004$, minh chứng thực nghiệm cho hiệu ứng Jahn-Teller.

2. Cơ chế hấp phụ phân ly $H_2$ và hoạt hóa $CO_2$: Trên các tâm xúc tác $Ni_5$, phân tử $H_2$ bị hấp phụ phân ly tự phát không cần vượt qua hàng rào năng lượng lớn, tạo thành 2 nguyên tử hydro liên kết tại các đỉnh hoặc cạnh của cluster. Đối với $CO_2$, quá trình diễn ra là hấp phụ hóa học tỏa nhiệt; liên kết $C=O$ nhận điện tích từ cluster kim loại, làm giảm mật độ điện tử liên kết và khiến góc $O-C-O$ bị uốn cong mạnh, tạo điều kiện thuận lợi cho phản ứng cắt đứt liên kết $C-O$.

3. Xác lập con đường phản ứng tối ưu (CO-route) cho quá trình tạo Methane: "Đã đề xuất và tính toán cơ chế hydrogen hóa $CO_2$ tạo thành $CH_4, CH_3OH, HCHO$ và $HCOOH$ trên ba hệ xúc tác $Ni_5, Ni_5/MgO$ và $Ni_5/AC$ gồm 33 đường phản ứng, 93 bước phản ứng... CH4 là sản phẩm hydrogen hóa $CO_2$ ưu tiên nhất. Đường phản ứng tạo thành $CH_4$ tối ưu là đường phản ứng đi qua sản phẩm trung gian là CO (CO route)." Các con đường tạo methanol ($CH_3OH$), formaldehyde ($HCHO$) và acid formic ($HCOOH$) đều có rào cản thế năng $E_a$ cao hơn hoặc kém thuận lợi về mặt nhiệt động học ($\Delta E$).

4. Định danh các giai đoạn quyết định tốc độ phản ứng (RDS): Trên bề mặt thế năng PES của cả ba hệ xúc tác, luận án định lượng chính xác năng lượng hoạt hóa của từng bước cơ bản. Giai đoạn phân cắt liên kết $C-O$ trong các tiểu phân trung gian chứa carbon và hydro (như $HCOH^* \rightarrow CH^* + OH^$ hoặc $CO^ \rightarrow C^* + O^*$) đòi hỏi năng lượng hoạt hóa cao nhất, đóng vai trò kiểm soát động học toàn bộ chuỗi phản ứng methane hóa.

5. Làm rõ vai trò điều biến điện tử của chất mang $MgO$ và $AC$: Khi đưa $Ni_5$ lên chất mang $MgO$ và $AC$, khả năng hấp phụ $CO_2$ và $H_2$ giảm nhẹ do sự chuyển dịch mật độ điện tích từ cluster kim loại sang chất mang. Tuy nhiên, chất mang giúp phân tán các tâm hoạt động, làm thay đổi cấu trúc electron d của Ni, giảm thiểu hiện tượng thiêu kết cụm nguyên tử và ngăn chặn sự hình thành phức carbonyl độc hại $Ni(CO)_5$, nâng cao độ bền và độ chọn lọc tổng thể của xúc tác.

Implications đa chiều

• Ý nghĩa học thuật và lý thuyết: Khẳng định tính đúng đắn của con đường trung gian $CO$ trên xúc tác nanocluster nickel; cung cấp bộ cơ sở dữ liệu bề mặt thế năng PES định lượng chính xác phục vụ nghiên cứu động học vi mô.
• Đổi mới phương pháp luận: Chuẩn hóa quy trình tính toán DFT/CI-NEB kết hợp phân tích điện tích Voronoi và bậc liên kết Mayer cho các hệ xúc tác dị thể quy mô lớn (>90 nguyên tử).
• Ứng dụng công nghệ và thực tiễn: Cung cấp nguyên lý thiết kế hệ xúc tác $Ni/MgO$ và $Ni/AC$ có độ chọn lọc $CH_4$ cao, chi phí thấp, hoạt động ổn định ở nhiệt độ thấp để thay thế các xúc tác kim loại quý đắt đỏ (Ru, Rh, Pt).
• Hỗ trợ chính sách khí hậu: Cung cấp luận cứ khoa học vững chắc cho việc phát triển công nghệ thu hồi và tái sử dụng carbon (CCUS), đóng góp trực tiếp vào mục tiêu giảm phát thải ròng bằng 0 (Net Zero) theo khuyến nghị của Ủy ban Liên chính phủ về Biến đổi Khí hậu (IPCC) và Cơ quan Năng lượng Quốc tế (IEA).

Limitations và Future Research

Nhằm duy trì tính khách quan và liêm chính học thuật, luận án chỉ rõ các giới hạn phương pháp luận:

1. Giới hạn nhiệt độ tính toán: Toàn bộ các giá trị $\Delta E$ và $E_a$ thu được từ tính toán DFT tĩnh ở trạng thái cơ bản tại $0\text{ K}$, chưa xét đến các hiệu ứng entropy nhiệt động học và năng lượng dao động điểm không (ZPE) ở nhiệt độ phản ứng công nghiệp thực tế (thường từ $250 - 450^\circ\text{C}$).
2. Giới hạn kích thước và mô hình chất mang: Bề mặt chất mang $MgO(200)$ và $AC$ được xây dựng dưới dạng các phiến tinh thể phẳng lý tưởng hóa, chưa mô phỏng hết các khuyết tật mạng phức tạp (defects, bước mạng, góc cạnh) và độ xốp ba chiều bất trật tự trong vật liệu thực tế.
3. Hiệu ứng bao phủ bề mặt và dung môi: Mô hình tính toán khảo sát các phản ứng đơn phân tử cô lập, chưa xét đến mức độ che phủ bề mặt cao (high surface coverage) của các tiểu phân hydro và nước đồng sản phẩm ($H_2O$), những yếu tố có thể làm thay đổi động học bề mặt.

Chương trình nghiên cứu mở rộng trong tương lai bao gồm:

• Mở rộng tính toán mô hình vi động học (microkinetic modeling) kết hợp động lực học phân tử ab initio (AIMD) tại nhiệt độ thực ($500 - 700\text{ K}$) để thu được các thông số vận tốc phản ứng và tần số luân chuyển (TOF).
• Nghiên cứu cơ chế phản ứng trên các cụm cluster kim loại kích thước lớn hơn ($Ni_{10}, Ni_{13}$) hoặc các cluster hợp kim đa kim loại pha tạp dị nguyên tử ($Ni-Fe, Ni-Ce, Ni-Co$) nhằm tối ưu hóa hoạt tính xúc tác.
• Nghiên cứu chi tiết cơ chế kháng tạo muội carbon (coking resistance) và độ bền chống ngộ độc carbonyl trong điều kiện phản ứng dài hạn.

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

Công trình nghiên cứu tạo ra những tác động sâu rộng trên nhiều bình diện:

• Tác động học thuật: Là tài liệu tham khảo chuẩn mực trong lĩnh vực hóa học tính toán và xúc tác hấp phụ tại Việt Nam; mở ra tiềm năng trích dẫn cao trong các nghiên cứu quốc tế về chuyển hóa $CO_2$.
• Chuyển đổi công nghiệp năng lượng: Mở đường cho việc thương mại hóa công nghệ tổng hợp khí thiên nhiên nhân tạo (SNG - Synthetic Natural Gas) từ nguồn phát thải $CO_2$ công nghiệp kết hợp với hydro xanh từ quá trình điện phân nước.
• Giá trị xã hội và môi trường: Góp phần giảm thiểu phát thải hàng triệu tấn khí nhà kính, tận dụng nguồn phụ phẩm nông nghiệp và chất thải công nghiệp để sản xuất than hoạt tính làm chất mang xúc tác giá rẻ.

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

• Nghiên cứu sinh và Học viên cao học: Tiếp cận quy trình tính toán lượng tử hiện đại, phương pháp CI-NEB định vị trạng thái chuyển tiếp và cách xây dựng bề mặt thế năng phản ứng.
• Các nhà khoa học chuyên ngành Hóa lý và Xúc tác: Kế thừa hệ thống dữ liệu nhiệt động và động học của 33 đường phản ứng để phát triển các lý thuyết xúc tác dị thể mới.
• Kỹ sư R&D trong ngành Hóa dầu và Năng lượng xanh: Ứng dụng các thông số tối ưu về tỷ lệ kim loại/chất mang để chế tạo xúc tác công nghiệp $Ni/MgO$ và $Ni/AC$ hiệu năng cao.
• Nhà hoạch định chính sách môi trường: Nắm bắt cơ sở khoa học để xây dựng lộ trình đầu tư công nghệ Power-to-Gas và dự án chuyển dịch năng lượng tuần hoàn.

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?
Trả lời: Đóng góp độc đáo nhất là việc mở rộng định lý Jahn-Teller và lý thuyết liên kết kim loại chuyển tiếp trên cấu trúc nanocluster $Ni_5$. Luận án chứng minh $Ni_5$ lưỡng tháp tam giác bị biến dạng phi đối xứng để giải tỏa sự suy biến orbital d, giảm số electron độc thân từ 10 xuống 8 với năng lượng liên kết $E_b = 3,60\text{ eV}$, tạo ra các tâm hoạt tính có mật độ điện tích không đồng nhất tối ưu cho việc kích hoạt phân tử $CO_2$.

2. Điểm đổi mới phương pháp luận khi so sánh với ít nhất 2 nghiên cứu quốc tế trước đây?
Trả lời: So với nghiên cứu của Grabow & Mavrikakis (2011) trên mặt phẳng $Cu(111)$ và Ren và cộng sự (2015) trên mặt phẳng $Ni(111)$, luận án đã đổi mới phương pháp luận bằng cách mô hình hóa hệ xúc tác nanocluster phân tán đa chiều ($Ni_5$) gắn trên chất mang dị thể phức tạp ($MgO$ và $AC$), sử dụng thuật toán CI-NEB 7 ảnh hội tụ đồng thời kết hợp phân tích điện tích Voronoi và bậc liên kết Mayer không phụ thuộc bộ hàm cơ sở.

3. Phát hiện bất ngờ nhất được hỗ trợ bởi dữ liệu tính toán là gì?
Trả lời: Phát hiện bất ngờ nhất là việc đưa cluster $Ni_5$ lên chất mang $MgO$ và $AC$ làm giảm nhẹ năng lượng hấp phụ ban đầu của $CO_2$ và $H_2$, nhưng lại tối ưu hóa hàng rào năng lượng hoạt hóa của các bước phân cắt liên kết tiếp theo trong chuỗi phản ứng $CO-route$, đồng thời ngăn chặn hoàn toàn xu hướng tạo phức carbonyl độc hại $Ni(CO)_5$.

4. Quy trình tái lập (Replication protocol) có được cung cấp đầy đủ không?
Trả lời: Luận án cung cấp quy trình tái lập hoàn chỉnh: phần mềm QUANTUM/SIESTA; phiếm hàm GGA-PBE; bộ hàm cơ sở DZP; thế giả Kleinman-Bylander Troullier-Martins với ngưỡng cắt sóng phẳng $2040,75\text{ eV}$; hộp mô phỏng $25\times 25\times 25\text{ \AA}^3$; tiêu chuẩn hội tụ lực $0,05\text{ eV/\AA}$ (tối ưu hóa) và $0,1\text{ eV/\AA}$ (CI-NEB).

5. Chương trình nghiên cứu 10 năm được định hình như thế nào từ kết quả này?
Trả lời: Định hình lộ trình 3 giai đoạn: (i) 1-3 năm: Mô hình hóa vi động học ở nhiệt độ cao và khảo sát cụm hợp kim đa kim ($Ni-Fe, Ni-Co$); (ii) 4-6 năm: Tổng hợp thực nghiệm vật liệu nano $Ni_5/AC$ và kiểm chứng hoạt tính trong lò phản ứng dòng liên tục; (iii) 7-10 năm: Tích hợp hệ xúc tác vào trạm thí nghiệm pilot Power-to-Gas quy mô bán công nghiệp.

Kết luận

1. Xác lập cấu hình lưỡng tháp tam giác không đối xứng ($Ni_5(8)$) là cấu trúc bền vững nhất của cluster $Ni_5$ với 8 electron độc thân, $E_b = 3,60\text{ eV}$, minh chứng cho hiệu ứng biến dạng Jahn-Teller.
2. Chứng minh cơ chế hấp phụ phân ly tự phát của $H_2$ và hấp phụ hóa học kích hoạt liên kết $C=O$ của $CO_2$ thông qua sự chuyển dịch điện tích từ orbital d kim loại vào orbital phản liên kết LUMO $2\pi^*$.
3. Khảo sát toàn diện 33 đường phản ứng với 93 bước cơ bản, khẳng định con đường trung gian $CO$ (CO route) là đường phản ứng tối ưu nhất để tạo ra khí methane ($CH_4$) trên cả ba hệ xúc tác $Ni_5$, $Ni_5/MgO$ và $Ni_5/AC$.
4. Định danh chính xác các giai đoạn phân cắt liên kết $C-O$ trong các hợp chất trung gian là bước quyết định tốc độ phản ứng (RDS).
5. Làm sáng tỏ vai trò của chất mang $MgO$ và than hoạt tính $AC$ trong việc điều biến cấu trúc electron d của cụm kim loại, nâng cao độ phân tán và triệt tiêu nguy cơ ngộ độc carbonyl.
6. Mở ra 3 hướng nghiên cứu mới: tính toán vi động học cho xúc tác cluster đa kim, phát triển vật liệu xúc tác công nghiệp giá rẻ cho phản ứng Sabatier, và ứng dụng mô phỏng lượng tử trong thiết kế vật liệu chuyển hóa khí nhà kính.