Tổng quan luận án

Luận án tiến sĩ chuyên ngành Vật liệu điện tử (mã số: 9440123) của nghiên cứu sinh Ngô Thị Lan, được thực hiện tại Học viện Khoa học và Công nghệ thuộc Viện Hàn lâm Khoa học và Công nghệ Việt Nam dưới sự hướng dẫn khoa học của PGS. Nguyễn Văn Đăng và PGS. Nguyễn Thanh Tùng (Hà Nội, 2024), mang tiêu đề: "Nghiên cứu tương tác vật lý giữa điện tử tự do và điện tử định xứ trong các hệ nano cluster hợp kim Au9M2+ (M = Sc-Ni) và AgnCr (n = 2-12) bằng phương pháp phiếm hàm mật độ".

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

Nano cluster (cụm nguyên tử kích thước vài nguyên tử đến vài nanomet) thể hiện các hiệu ứng bề mặt và hiện tượng giam giữ lượng tử, dẫn đến cấu trúc và tính chất hóa lý khác biệt căn bản so với vật liệu dạng khối. Khi kết hợp kim loại quý (Au, Ag) với kim loại chuyển tiếp dãy 3d (Sc đến Ni), hệ nano cluster hợp kim xuất hiện trạng thái đồng tồn tại giữa hai nhóm điện tử:

  • Điện tử hóa trị chuyển động tự do bắt nguồn từ phân lớp $s$ của kim loại quý, hình thành các mức năng lượng lớp vỏ điện tử chung chi phối độ bền cấu trúc của toàn cụm.
  • Điện tử định xứ phân bố cục bộ trên các orbital phân lớp $d$ chưa bão hòa của kim loại chuyển tiếp, đóng vai trò nguồn gốc tạo ra mômen từ spin và các tâm hoạt tính xúc tác.

Khoảng trống nghiên cứu được xác định rõ: các nghiên cứu trước đây (điển hình như luận án của TS. Nguyễn Thị Mai, 2021) đã phân tích cấu trúc điện tử và độ bền của một số hệ cluster như $\text{Si}_n\text{Mn}^{2+}$, $\text{Co}n\text{O}m^+$, $\text{MnCr}$ và $\text{Au}{19}\text{M}$ ($\text{M} = \text{Sc-Ni}$), chứng minh $\text{Au}{19}\text{Cr}$ là một siêu nguyên tử bền vững với 20 điện tử điền đầy và 5 điện tử chưa ghép cặp trên orbital $3d$. Tuy nhiên, các công trình này tồn tại hạn chế là chưa làm rõ khả năng xúc tác và chưa xác định các trạng thái động học trong phản ứng liên kết với hydrogen ($\text{H}2$). Bên cạnh đó, hệ $\text{Au}{10}^{2+}$ có cấu trúc tứ diện đối xứng cao ($T_d$), khe năng lượng $\text{HOMO-LUMO}$ lớn ($3,88\text{ eV}$) và 8 điện tử hóa trị điền đầy ($1S^2 1P^6$) mang tính bền vững nhưng vỏ đóng kín làm hạn chế các tính chất điện tử mong muốn. Việc thay thế 1 nguyên tử Au bằng kim loại chuyển tiếp M để tạo thành $\text{Au}_9\text{M}^{2+}$ nhằm duy trì độ bền vỏ đồng thời tạo các đặc tính điện tử mới vẫn chưa được giải quyết đầy đủ. Đối với hệ $\text{Ag}_n\text{Cr}$ ($n = 2-12$), sự biến đổi của tương tác $s-d$ theo kích thước cụm và số lượng điện tử $s$ khi pha tạp nguyên tử Cr ([$\text{4s}^1\text{3d}^5$]) cũng chưa được làm rõ một cách hệ thống.

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

Luận án xác định hai mục tiêu chính:

  1. Làm rõ tương tác giữa các điện tử tự do và điện tử định xứ trong các hệ nano cluster kim loại quý pha tạp kim loại chuyển tiếp $\text{Au}_9\text{M}^{2+}$ ($\text{M} = \text{Sc-Ni}$) và $\text{Ag}_n\text{Cr}$ ($n = 2-12$).
  2. Xác định ảnh hưởng của cấu trúc điện tử tới sự phát triển cấu trúc hình học bền, độ bền vững (năng lượng liên kết, năng lượng phân ly) và tác động của tương tác $s-d$ tới động học phản ứng hấp phụ/phân ly $\text{H}_2$ của các hệ nano cluster hợp kim nêu trên.

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

  • Đối tượng nghiên cứu:
    1. Hệ nano cluster hợp kim $\text{Au}_9\text{M}^{2+}$ với $\text{M}$ là các nguyên tử kim loại chuyển tiếp dãy $3d$ gồm: $\text{Sc}$, $\text{Ti}$, $\text{V}$, $\text{Cr}$, $\text{Mn}$, $\text{Fe}$, $\text{Co}$ và $\text{Ni}$.
    2. Hệ nano cluster hợp kim $\text{Ag}_n\text{Cr}$ với kích thước $n$ thay đổi từ 2 đến 12.
  • Phạm vi nghiên cứu:
    • Khảo sát tương tác vật lý $s-d$, cấu trúc hình học bền ở trạng thái cơ bản, cấu trúc điện tử, độ bội spin, mật độ trạng thái (DOS, pDOS), giản đồ orbital phân tử (MO).
    • Đánh giá độ bền tĩnh và độ bền nhiệt động lực học thông qua năng lượng liên kết trung bình ($\text{BE}$), chênh lệch năng lượng bậc hai ($\Delta_2\text{E}$), năng lượng phân ly ($\text{DE}$).
    • Động học quá trình hấp phụ phân tử $\text{H}2$ và hấp phụ phân ly thành 2 nguyên tử H trên bề mặt cluster (năng lượng hấp phụ $E{ads}$, rào cản hoạt hóa, trạng thái chuyển tiếp TS, trạng thái trung gian IS, tọa độ phản ứng nội tại IRC).
    • Không gian nghiên cứu: Mô phỏng hóa tính toán lượng tử pha khí ở nhiệt độ $0\text{ K}$ (có hiệu chỉnh năng lượng dao động điểm không ZPE) trên hệ thống máy tính hiệu năng cao.

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

Các hướng nghiên cứu tiền đề

Lịch sử nghiên cứu nano cluster và sự hình thành lớp vỏ điện tử đã được tác giả hệ thống hóa qua các mốc công trình:

  • Năm 1967, Robbins và cộng sự lần đầu tiên tạo ra các cụm nguyên tử trong pha khí.
  • Năm 1984, Knight và cộng sự phát hiện sự biến đổi tính chất theo kích thước trên cụm $\text{Na}_N$ thông qua phổ khối, xác định các kích thước ổn định nổi bật tại $N = 8, 20, 40, 58, 92$, đặt nền tảng cho mô hình vỏ điện tử Jellium hình cầu với sự điền đầy các phân lớp $1S, 1P, 1D, 2S, 1F...$
  • Năm 1985, Clemenger và Nilsson phát triển mô hình vỏ elip (biến dạng nén - oblate và kéo dãn - prolate) nhằm mô tả sự tách mức năng lượng và hiệu ứng Jahn-Teller ở các cụm có lớp vỏ mở.
  • Khái niệm "siêu nguyên tử" (superatoms) hình thành qua các phát hiện: $\text{C}{60}$ (Kroto và cộng sự, 1985, khe $\text{HOMO-LUMO} \approx 1,8\text{ eV}$), cụm $\text{Al}{13}$ (Harms và cộng sự, 1990, cấu trúc nhị thập diện bền tương tự halogen), cụm $\text{Au}_{20}$ (J. Li và cộng sự, 2003, cấu trúc tứ diện $T_d$, khe $\text{HOMO-LUMO} \approx 1,7\text{ eV}$, 20 điện tử điền đầy cấu hình $1S^2 1P^6 2S^2 1D^{10}$).
  • Năm 2003, Kumar và cộng sự chứng minh siêu nguyên tử mang từ tính trên các hệ $\text{X}n\text{M}$ ($\text{X} = \text{Si}, \text{Ge}$; $\text{M} = \text{Be}, \text{Mg}, \text{Zn}, \text{Cd}, \text{Mn}$; $n = 8-14$), điển hình $\text{M@X}{12}$ có mômen từ $\approx 5\mu_B$.
  • Năm 2006, Reveles và cộng sự đưa ra lý thuyết về sự đồng tồn tại giữa lớp vỏ hóa trị điền đầy và lớp vỏ từ tính tạo bởi điện tử định xứ trên hệ $\text{Na}_8\text{V}$.
  • Năm 2005, Li và cộng sự dùng kỹ thuật bay hơi laser và quang phổ quang điện tử (PES) nghiên cứu các cụm $\text{M@Au}_6$ ($\text{M} = \text{Ti}, \text{V}, \text{Cr}$), xác định mômen từ spin từ $2\mu_B$ đến $4\mu_B$. Các nghiên cứu quang phân ly trên $\text{Au}n\text{M}$ chỉ ra $\text{Sc}, \text{Ti}$ đóng góp cả điện tử $4s$ và $3d$ vào vỏ chung, trong khi $\text{Cr}, \text{Mn}, \text{Fe}, \text{Co}, \text{Ni}$ chỉ đóng góp $4s$, giữ lại các điện tử $3d$ định xứ. Các công trình trên $\text{Au}{24}\text{M}$ ghi nhận mômen từ spin tương ứng $3, 6, 5, 4, 3, 2\mu_B$ cho $\text{M} = \text{V}, \text{Cr}, \text{Mn}, \text{Fe}, \text{Co}, \text{Ni}$.
  • Về hấp phụ và xúc tác hydrogen: Nghiên cứu của J. Vanbuel và cộng sự trên $\text{Al}_n\text{V}_2\text{H}2^+$, X. Ma và cộng sự trên $\text{Mg}{55}$ pha tạp $\text{Ti}, \text{Nb}$, Phala, Buckart, Lang trên hệ $\text{Au}_n\text{-H}$, Fang và Kuang trên $\text{Au}_n\text{Pt}$, Mondal, Li, Perera trên $\text{Ag}_n\text{H}$.

Tại Việt Nam, các hướng tính toán lượng tử cấu trúc cluster đã được triển khai bởi các nhóm nghiên cứu: nhóm PGS. Nguyễn Thị Minh Huệ (ĐH Sư phạm Hà Nội - động học, nhiệt phân tử pha khí), PGS. Ngô Tuấn Cường (ĐH Sư phạm Hà Nội - quang và cấu trúc $\text{Ag}, \text{Si}$), PGS. Vũ Thị Ngân (ĐH Quy Nhơn - cụm $\text{Si}_n$), PGS. Phạm Cẩm Nam (ĐH Bách Khoa Đà Nẵng - $\text{Si}$ pha tạp $\text{Au}, \text{Ag}, \text{Cu}$), PGS. Trần Văn Tân (ĐH Đồng Tháp - $\text{Si}, \text{Ge}$ pha tạp kim loại chuyển tiếp), TS. Nguyễn Minh Tâm (ĐH Phan Thiết - cụm Boron pha tạp), TS. Phạm Vũ Nhật (ĐH Cần Thơ - cụm $\text{Au}, \text{Ag}$ và phổ hồng ngoại), cùng nhóm nghiên cứu tại Viện Khoa học vật liệu (VAST) hợp tác với ĐH KU Leuven (Bỉ).

Vị trí và khoảng trống luận án giải quyết

Kế thừa các kết quả về cụm tứ diện $\text{Au}_{19}\text{Cr}$ và các hệ kim loại quý pha tạp từ luận án của TS. Nguyễn Thị Mai (2021), luận án của NCS. Ngô Thị Lan tập trung giải quyết khoảng trống học thuật về:

  1. Cơ chế tương tác chi tiết giữa orbital $s$ của kim loại quý và orbital $d$ của kim loại chuyển tiếp trong các hệ cluster kích thước nhỏ $\text{Au}_9\text{M}^{2+}$ và $\text{Ag}_n\text{Cr}$.
  2. Xác định định lượng đường phản ứng, trạng thái chuyển tiếp, năng lượng rào cản hoạt hóa trong quá trình tương tác, hấp phụ phân tử và phân ly $\text{H}_2$ trên các tâm cluster.

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

Khung lý thuyết sử dụng

  • Lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT): Sử dụng các phép gần đúng gradient tổng quát (GGA), phiếm hàm trao đổi tương quan BP86 và phiếm hàm có hiệu chỉnh tương tác xa LC-BP86 kết hợp với các bộ hàm cơ sở thích hợp (như LanL2DZ, Def2-TZVP...).
  • Lý thuyết cấu trúc vỏ điện tử:
    • Mô hình Jellium: Giải phương trình Schrödinger cho hạt chuyển động trong giếng thế 3 chiều: $$\hat{H} = \frac{\hat{p}^2}{2m} + \frac{1}{2}m\omega_0^2 r^2 - U\hbar\omega_0 \left[ \hat{l}^2 - \frac{n(n+3)}{2} \right]$$ với các mức năng lượng sắp xếp theo trật tự: $1s^2 1p^6 1d^{10} 2s^2 1f^{14} 2p^6 1g^{18} 2d^{10} 1h^{22} 3s^2 2f^{14} 1i^{26} 3p^6...$ Các số kỳ diệu điền đầy lớp vỏ gồm: 2, 8, 20, 34, 40...
    • Mô hình vỏ elip Clemenger-Nilsson: Giải thích sự phân tách suy biến $(2l+1)$ do biến dạng hình học nén (oblate) hoặc kéo dãn (prolate).
  • Lý thuyết vỏ hình học (Geometric Shell Theory): Mô tả các dạng hình học đối xứng cao bền vững độc lập với vỏ điện tử như lưỡng tháp ngũ giác ($N = 7$), nhị thập diện ($N = 13$), lập phương tâm khối BCC ($N = 15$), hai nhị thập diện lồng ($N = 19$).
  • Phương pháp phân tích cấu trúc điện tử: Mật độ trạng thái toàn phần (DOS), mật độ trạng thái từng phần (pDOS), phân tích orbital phân tử tự nhiên (NBO), xác định độ rộng vùng cấm HOMO-LUMO.

Các công thức đánh giá độ bền và năng lượng

  1. Năng lượng liên kết trung bình ($\text{BE}$): $$BE(\text{A}_n\text{B}_m) = \frac{n E(\text{A}) + m E(\text{B}) - E(\text{A}_n\text{B}_m)}{n + m}$$
  2. Chênh lệch năng lượng bậc hai ($\Delta_2\text{E}$): $$\Delta_2 E(\text{A}n) = E(\text{A}{n+1}) + E(\text{A}_{n-1}) - 2E(\text{A}_n)$$
  3. Năng lượng phân ly ($\text{DE}$): Năng lượng tối thiểu cần cung cấp để phá vỡ liên kết tạo thành hai phân mảnh nguyên tử, phân tử hoặc cluster con theo các kênh phân ly cạnh tranh.
  4. Năng lượng hấp phụ ($E_{ads}$): Xác định mức độ liên kết giữa phân tử $\text{H}_2$ hoặc 2 nguyên tử H với khung nano cluster.

Phương pháp tính toán và nguồn dữ liệu

  • Công cụ phần mềm: Gói phần mềm Gaussian 09 kết hợp phần mềm hiển thị GaussView.
  • Quy trình tính toán:
    1. Tối ưu hóa cấu trúc hình học và tìm cấu trúc bền thấp năng lượng nhất ở các trạng thái spin khác nhau ($2S+1$).
    2. Tính toán tần số dao động để đảm bảo cấu trúc không có tần số ảo (trạng thái cực tiểu năng lượng) hoặc có duy nhất 1 tần số ảo (trạng thái chuyển tiếp TS).
    3. Quét tọa độ phản ứng nội tại (IRC) để liên kết trạng thái chuyển tiếp với trạng thái đầu và trạng thái cuối của phản ứng hấp phụ/phân ly $\text{H}_2$.
  • Đánh giá độ tin cậy của phương pháp: Tác giả tiến hành tính toán kiểm chứng độ dài liên kết ($\text{BL}$), độ bội spin và năng lượng phân ly ($\text{DE}$) trên các hệ đime $\text{Au}_2$, $\text{AuM}$, $\text{AuH}$, $\text{MH}$, $\text{Ag}_2$, $\text{Cr}_2$, $\text{AgCr}$ bằng nhiều phiếm hàm và bộ hàm cơ sở khác nhau, so sánh trực tiếp với số liệu thực nghiệm để lựa chọn bộ thông số tối ưu trước khi tính toán cho các cluster phức tạp.
+-----------------------------------------------------------------------------------+
|                            QUY TRÌNH TÍNH TOÁN DFT                                |
|                                                                                   |
|  [Khởi tạo cấu trúc] ---> [Tối ưu hóa hình học & Spin] ---> [Tính tần số dao động]|
|                                                                     |             |
|                                                                     v             |
|  [Động học phản ứng IRC] <-- [Xác định TS/IS/Hóa hấp phụ] <-- [Phân tích DOS/pDOS]|
+-----------------------------------------------------------------------------------+

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

Chương 1: Tổng quan về nano cluster hợp kim

Chương 1 hệ thống hóa cơ sở lý thuyết và các công trình liên quan đến nano cluster kim loại và hợp kim:

  • Trình bày khái niệm nano cluster, hiệu ứng giam giữ lượng tử và bản chất của trạng thái trung gian giữa nguyên tử đơn lẻ và vật liệu khối.
  • Phân tích chi tiết mô hình Jellium vỏ cầu, mô hình giếng thế điều hòa có số hạng hiệu chỉnh $U[l^2 - n(n+3)/2]$, mô hình Clemenger-Nilsson cho vỏ biến dạng elip và lý thuyết vỏ hình học đối với kim loại chuyển tiếp.
  • Tổng kết các đặc tính vật lý quan trọng: Độ bền tĩnh ($\text{BE}, \Delta_2\text{E}$), độ bền nhiệt động ($\text{DE}$ qua các kênh phân ly), từ tính phụ thuộc kích thước/cấu hình (như mômen từ tăng cao ở cụm nhỏ $\text{Fe}, \text{Co}, \text{Ni}, \text{Mn}, \text{Cr}$), và hoạt tính xúc tác (phân tách nước trên $\text{Au}n$, sản xuất $\text{H}2$ trên $\text{Au}{12}\text{Ag}{32}/\text{TiO}_2$).
  • Làm rõ bản chất tương tác $s-d$ trong cluster hợp kim kim loại quý - kim loại chuyển tiếp: điện tử $s$ tạo vỏ tự do, điện tử $d$ định xứ tạo tâm từ và xúc tác. Phân tích hiện trạng nghiên cứu trong và ngoài nước, chỉ rõ cơ sở lựa chọn hai hệ $\text{Au}_9\text{M}^{2+}$ và $\text{Ag}_n\text{Cr}$.

Chương 2: Phương pháp nghiên cứu

Chương 2 trình bày chi tiết các phương pháp tính toán hóa học lượng tử:

  • Trình bày cơ sở toán học của DFT, các gần đúng mật độ cục bộ (LDA, LSDA), gần đúng gradient tổng quát (GGA), các hiệu chỉnh tương tác khoảng cách xa (Long-range correction - LC).
  • Thiết lập quy trình tính toán trên phần mềm Gaussian 09 và GaussView từ bước tối ưu hóa cấu trúc, khảo sát các trạng thái đa bội spin, đến tính toán năng lượng liên kết, phân tích NBO, DOS, pDOS và khảo sát bề mặt thế năng phản ứng (IRC).
  • Thực hiện đánh giá độ tin cậy của các phiếm hàm và bộ hàm cơ sở thông qua so sánh các thông số tính toán ($\text{BL}, \text{DE}$, spin) với số liệu thực nghiệm của các hệ đime: $\text{Au}_2$, $\text{AuM}$, $\text{AuH}$, $\text{MH}$ ($\text{M} = \text{Sc-Ni}$), $\text{Ag}_2$, $\text{Cr}_2$, $\text{AgCr}$. Kết quả khẳng định độ chính xác và tính phù hợp của phiếm hàm BP86 và LC-BP86 trong việc mô tả tương tác kim loại - kim loại và kim loại - hydro.

Chương 3: Tương tác điện tử s-d trong các hệ nano cluster hợp kim

Chương 3 đi sâu phân tích cấu trúc điện tử và bản chất tương tác $s-d$ trên hai hệ đối tượng:

  1. Hệ nano cluster hợp kim $\text{Au}_9\text{M}^{2+}$ ($\text{M} = \text{Sc-Ni}$):
    • Xác định cấu hình điện tử lớp ngoài của nguyên tử M và cấu trúc điện tử của cụm. Phân tích trạng thái spin bền $(2S+1)$ và số lượng điện tử hóa trị chưa ghép cặp định xứ trên orbital $3d$.
    • Giản đồ orbital phân tử (MO) và phân bố mật độ trạng thái (DOS, pDOS) chỉ ra: các nguyên tử đầu dãy $3d$ ($\text{Sc}, \text{Ti}$) có xu hướng đóng góp nhiều điện tử hóa trị vào lớp vỏ điện tử tự do chung của cụm; trong khi các nguyên tử cuối dãy ($\text{Cr}, \text{Mn}, \text{Fe}, \text{Co}, \text{Ni}$) giữ lại các điện tử $3d$ định xứ không tham gia vào vỏ tự do, tạo ra các trạng thái spin cao.
  2. Hệ nano cluster hợp kim $\text{Ag}_n\text{Cr}$ ($n = 2-12$):
    • Khảo sát sự biến đổi cấu trúc điện tử khi tăng dần số nguyên tử Ag từ $n = 2$ đến $n = 12$.
    • Phân tích mức độ lai hóa orbital giữa $5s$ của Ag và $3d$ của Cr, xác định sự biến đổi số điện tử định xứ trên nguyên tử Cr theo kích thước cụm và quy luật đóng/mở lớp vỏ điện tử tự do.
Hệ nano cluster Cấu hình/Đặc trưng tương tác điện tử Bản chất lớp vỏ điện tử & Trạng thái spin
$\text{Au}_9\text{M}^{2+}$ ($\text{M} = \text{Sc, Ti}$) Điện tử $4s$ và một phần $3d$ tham gia tạo vỏ Vỏ điện tử chung điền đầy, số điện tử định xứ $3d$ ít, spin thấp
$\text{Au}_9\text{M}^{2+}$ ($\text{M} = \text{Cr, Mn, Fe, Co, Ni}$) Điện tử $4s$ tham gia tạo vỏ, điện tử $3d$ định xứ Tồn tại đồng thời vỏ điện tử tự do và các điện tử $3d$ định xứ, spin cao
$\text{Ag}_n\text{Cr}$ ($n = 2-12$) Tương tác giữa $n$ điện tử $5s$ (Ag) và phân lớp $4s^1 3d^5$ (Cr) Biến đổi tuần hoàn số điện tử vỏ tự do và số điện tử $3d$ định xứ theo kích thước $n$

Chương 4: Ảnh hưởng của tương tác s-d đến các tính chất của nano cluster hợp kim

Chương 4 khảo sát mối quan hệ giữa cấu trúc điện tử $s-d$ với các đặc tính hình học, độ bền và hoạt tính xúc tác liên kết $\text{H}_2$:

  1. Hệ nano cluster $\text{Au}_9\text{M}^{2+}$ ($\text{M} = \text{Sc-Ni}$):
    • Xác định cấu trúc hình học bền: Sự thay thế nguyên tử M làm biến dạng cấu trúc tứ diện nguyên bản của $\text{Au}_{10}^{2+}$ tùy thuộc vào bán kính nguyên tử và mức độ lai hóa $s-d$.
    • Đánh giá độ bền vững qua năng lượng liên kết trung bình ($\text{BE}$) và năng lượng phân ly ($\text{DE}$) theo các kênh phân ly khác nhau.
    • Khảo sát tương tác với $\text{H}_2$: Xác định cấu trúc bền của $\text{Au}_9\text{M}^{2+}\text{-H}_2$ (hấp phụ phân tử) và $\text{Au}_9\text{M}^{2+}\text{-2H}$ (hấp phụ phân ly). Tính toán đường phản ứng và bề mặt năng lượng tự do cho các hệ $\text{Au}_9\text{Sc}^{2+}$, $\text{Au}_9\text{Ti}^{2+}$, $\text{Au}_9\text{Mn}^{2+}$, $\text{Au}_9\text{Fe}^{2+}$, xác định cơ chế và rào cản năng lượng của quá trình hoạt hóa phân ly liên kết $\text{H-H}$.
  2. Hệ nano cluster $\text{Ag}_n\text{Cr}$ ($n = 2-12$):
    • Sự tiến hóa cấu trúc hình học bền từ $n = 2$ đến $12$.
    • Phân tích các cực đại trên đường cong $\text{BE}$, chênh lệch năng lượng bậc hai $\Delta_2\text{E}$ và $\text{DE}$ để chỉ ra các kích thước cụm có độ ổn định nhiệt động vượt trội.
    • Khảo sát tương tác với $\text{H}2$: Xác định năng lượng hấp phụ $E{ads}$, sự giãn độ dài liên kết $\text{H-H}$ ($d_{\text{H-H}}$), phân bố DOS/pDOS khi có mặt $\text{H}_2$/$\text{2H}$. Tính toán đường phản ứng và năng lượng tự do tương đối cho quá trình hấp phụ và phân ly $\text{H}_2$ trên các cụm $\text{Ag}_n\text{Cr}$ ($n = 2-4, 6-7$).

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

Đóng góp mới về mặt khoa học và thực tiễn

  1. Về cấu trúc và tương tác điện tử: Xác định được cấu trúc hình học bền nhất ở trạng thái cơ bản và cấu trúc điện tử tương ứng của hai dãy nano cluster hợp kim $\text{Au}_9\text{M}^{2+}$ ($\text{M} = \text{Sc-Ni}$) và $\text{Ag}_n\text{Cr}$ ($n = 2-12$). Luận án làm sáng tỏ quy luật tương tác giữa điện tử tự do (từ kim loại quý Au, Ag) và điện tử định xứ (trên phân lớp $3d$ của kim loại chuyển tiếp), chứng minh sự tồn tại đồng thời của lớp vỏ điện tử tự do đóng kín và các điện tử $3d$ định xứ chưa ghép cặp.
  2. Về quy luật biến đổi tính chất: Làm rõ ảnh hưởng trực tiếp của mức độ lai hóa $s-d$ và số lượng điện tử $3d$ định xứ lên độ bền năng lượng ($\text{BE}, \Delta_2\text{E}, \text{DE}$) cũng như sự biến dạng hình học của cụm so với khung đối xứng chuẩn.
  3. Về động học phản ứng xúc tác hydrogen: Làm rõ cơ chế và xây dựng bức tranh đầy đủ về động học phản ứng tương tác, hấp phụ phân tử và phân ly $\text{H}_2$ trên bề mặt các cụm $\text{Au}_9\text{M}^{2+}$ và $\text{Ag}_n\text{Cr}$. Luận án xác định được các trạng thái trung gian (IS), trạng thái chuyển tiếp (TS) và rào cản năng lượng hoạt hóa, cung cấp cơ sở khoa học để giải thích vai trò xúc tác của nguyên tử pha tạp kim loại chuyển tiếp trong việc kích hoạt phá vỡ liên kết $\text{H-H}$.

Kiến nghị và giải pháp đề xuất

  • Đề xuất hoàn thiện và chuẩn hóa quy trình mô phỏng lượng tử kết hợp phiếm hàm BP86/LC-BP86 và phân tích IRC để khảo sát động học phản ứng trên các hệ nano cluster hợp kim pha tạp.
  • Đề xuất tiếp tục mở rộng nghiên cứu tương tác giữa các hệ cluster này với các phân tử khí khác ($\text{CO}, \text{O}_2, \text{CO}_2, \text{NO}_x$) nhằm đánh giá toàn diện tiềm năng xúc tác dị thể.
  • Kiến nghị triển khai các nghiên cứu thực nghiệm phối hợp (sử dụng kỹ thuật chùm phân tử, bẫy ion và phổ laser) để kiểm chứng trực tiếp các trạng thái cấu trúc và phổ hấp thụ đã dự đoán lý thuyết.

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

Hạn chế của công trình

  • Các phép tính toán mô phỏng được thực hiện trên các cụm nano cluster đơn lẻ ở pha khí (chân không) tại nhiệt độ $0\text{ K}$, chưa xét đến ảnh hưởng của môi trường dung môi hoặc hiệu ứng tương tác với chất nền mang (substrate) như $\text{TiO}_2, \text{Al}_2\text{O}_3$ hay vật liệu carbon thường gặp trong điều kiện làm việc thực tế.
  • Khảo sát tương tác xúc tác trong luận án mới chỉ giới hạn trên đối tượng phân tử hydrogen ($\text{H}_2$), chưa mở rộng sang các phân tử phản ứng khác hoặc phản ứng đa phân tử phức tạp.

Hướng nghiên cứu tiếp theo

  • Mở rộng tính toán tương tác $s-d$ trên các hệ nano cluster hợp kim có kích thước lớn hơn hoặc chứa từ hai nguyên tử pha tạp kim loại chuyển tiếp trở lên.
  • Nghiên cứu cơ chế xúc tác và hấp phụ của hệ cluster khi được đặt trên các bề mặt chất mang (supported clusters) để tiệm cận với điều kiện thực nghiệm xúc tác công nghiệp.
  • Khảo sát thêm các phản ứng oxy hóa khử và chuyển hóa năng lượng khác như phản ứng khử $\text{CO}_2$, phân tách nước ($\text{H}_2\text{O}$) hoặc phản ứng oxy hóa $\text{CO}$.

Giá trị tham khảo

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

  • Đối tượng thụ hưởng:
    • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Vật liệu điện tử, Vật lý chất rắn, Hóa học tính toán, Khoa học và Công nghệ nano.
    • Các kỹ sư và nhà nghiên cứu trong lĩnh vực vật liệu lưu trữ năng lượng (hydrogen storage), vật liệu spintronics và xúc tác dị thể kích thước nguyên tử.
  • Nội dung đáng tham khảo:
    • Quy trình và phương pháp luận tính toán DFT xác định cấu trúc cluster, trạng thái spin và phân tích mật độ trạng thái (DOS, pDOS) tại Chương 2 và Chương 3.
    • Số liệu chi tiết về năng lượng phân ly, năng lượng liên kết và tọa độ phản ứng phân ly $\text{H}_2$ trên các tâm kim loại chuyển tiếp dãy $3d$ tại Chương 4.
    • Phương pháp kiểm chứng độ tin cậy của mô hình tính toán lượng tử thông qua việc đối chiếu với số liệu thực nghiệm trên các hệ đime kim loại.

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

1. Bản chất của tương tác $s-d$ trong các hệ nano cluster kim loại quý pha tạp kim loại chuyển tiếp được luận án mô tả như thế nào?

Trong các hệ này, các nguyên tử kim loại quý (Au, Ag) đóng góp các điện tử hóa trị phân lớp $s$ có độ linh động cao để tạo thành đám mây điện tử tự do chung, điền đầy các mức năng lượng lớp vỏ Jellium và quyết định độ bền khung hình học. Trong khi đó, các nguyên tử kim loại chuyển tiếp dãy $3d$ sở hữu các điện tử phân lớp $d$ có tính định xứ cao. Tùy thuộc vào bản chất nguyên tố pha tạp, các điện tử $3d$ không tham gia hoặc chỉ tham gia một phần vào vỏ chung, phần còn lại tồn tại dưới dạng các điện tử chưa ghép cặp định xứ, tạo nên từ tính (mômen từ spin cao) và đóng vai trò tâm hoạt tính xúc tác.

2. Tại sao tác giả lại lựa chọn nghiên cứu hệ $\text{Au}9\text{M}^{2+}$ dựa trên khung cấu trúc của $\text{Au}{10}^{2+}$?

Cụm $\text{Au}{10}^{2+}$ nguyên bản có cấu trúc tứ diện đối xứng cao ($T_d$), khe năng lượng $\text{HOMO-LUMO}$ rất lớn ($3,88\text{ eV}$) và 8 điện tử hóa trị điền đầy cấu hình vỏ đóng kín $1S^2 1P^6$, mang độ bền vững tương tự một siêu nguyên tử khí trơ. Tuy nhiên, chính lớp vỏ đóng kín này làm cho $\text{Au}{10}^{2+}$ trơ về mặt hóa học và không có từ tính. Việc thay thế 1 nguyên tử Au bằng 1 nguyên tử kim loại chuyển tiếp M ($\text{Sc-Ni}$) tạo ra hệ $\text{Au}_9\text{M}^{2+}$ nhằm mục đích duy trì độ bền khung cấu trúc đồng thời đưa vào các điện tử $3d$ định xứ chưa bão hòa, tạo ra các tính chất điện tử, từ tính và hoạt tính xúc tác mới.

3. Sự đóng góp điện tử của các kim loại chuyển tiếp dãy $3d$ vào lớp vỏ điện tử tự do của cụm $\text{Au}_9\text{M}^{2+}$ thay đổi như thế nào từ Sc đến Ni?

Theo kết quả phân tích cấu trúc điện tử và orbital phân tử trong luận án, các nguyên tố nhẹ ở đầu dãy $3d$ ($\text{Sc}, \text{Ti}$) có xu hướng đóng góp cả điện tử $4s$ và một phần điện tử $3d$ vào lớp vỏ điện tử tự do chung của cụm, dẫn đến trạng thái spin cơ bản thấp. Ngược lại, các nguyên tố nặng hơn từ giữa đến cuối dãy ($\text{Cr}, \text{Mn}, \text{Fe}, \text{Co}, \text{Ni}$) chủ yếu chỉ đóng góp điện tử $4s$, các điện tử $3d$ còn lại duy trì trạng thái định xứ cục bộ trên orbital nguyên tử, tạo ra các trạng thái spin cao tương ứng với số điện tử chưa ghép cặp.

4. Quá trình tương tác giữa phân tử hydrogen ($\text{H}_2$) và các nano cluster được khảo sát thông qua những trạng thái nhiệt động và động học nào?

Luận án khảo sát quá trình tương tác qua hai giai đoạn chính: hấp phụ phân tử ($\text{Cluster-H}2$) và hấp phụ phân ly ($\text{Cluster-2H}$). Tác giả tính toán năng lượng hấp phụ $E{ads}$, độ giãn liên kết $\text{H-H}$, xác định cấu trúc trạng thái ban đầu, trạng thái trung gian (IS), trạng thái chuyển tiếp (TS) với tần số ảo duy nhất, và quét đường phản ứng nội tại (IRC) để xác định chính xác năng lượng kích hoạt rào cản phân ly $\text{H}_2$ trên từng tâm cluster.

5. Phương pháp tính toán DFT trong luận án được kiểm chứng độ tin cậy bằng cách nào trước khi áp dụng cho các cụm phức tạp?

Tác giả đã tiến hành tính toán kiểm chứng độ dài liên kết ($\text{BL}$), độ bội spin và năng lượng phân ly ($\text{DE}$) trên các hệ phân tử đime mẫu gồm $\text{Au}_2$, $\text{AuM}$, $\text{AuH}$, $\text{MH}$ ($\text{M} = \text{Sc-Ni}$), $\text{Ag}_2$, $\text{Cr}_2$, $\text{AgCr}$ bằng nhiều phiếm hàm (như BP86, LC-BP86) kết hợp các bộ hàm cơ sở khác nhau. Các kết quả tính toán được so sánh trực tiếp với dữ liệu thực nghiệm sẵn có để chứng minh độ chính xác và tính nhất quán của phương pháp trước khi triển khai trên hệ $\text{Au}_9\text{M}^{2+}$ và $\text{Ag}_n\text{Cr}$.

Kết luận

Luận án tiến sĩ của tác giả Ngô Thị Lan là công trình nghiên cứu hệ thống bằng phương pháp phiếm hàm mật độ (DFT) về tương tác vật lý $s-d$ trong các hệ nano cluster hợp kim $\text{Au}_9\text{M}^{2+}$ ($\text{M} = \text{Sc-Ni}$) và $\text{Ag}_n\text{Cr}$ ($n = 2-12$). Tác giả đã làm sáng tỏ cơ chế đồng tồn tại giữa lớp vỏ điện tử tự do điền đầy và các điện tử $3d$ định xứ, giải thích mối quan hệ giữa cấu trúc điện tử với độ bền tĩnh, độ bền nhiệt động và sự biến đổi hình học của cụm. Công trình đã làm rõ động học phản ứng, xác định rào cản năng lượng của quá trình hấp phụ và phân ly phân tử $\text{H}_2$, cung cấp cơ sở dữ liệu lý thuyết đáng tin cậy cho định hướng thiết kế các vật liệu siêu nguyên tử ứng dụng trong xúc tác và lưu trữ năng lượng.