Tổng quan luận án

Nghiên cứu về các cụm nguyên tử (atomic clusters) có kích thước dưới nanomet là một lĩnh vực trọng tâm trong khoa học vật liệu hiện đại. Ở kích thước từ vài đến vài chục nguyên tử, tính chất vật lý và hóa học của vật liệu không tuân theo các quy luật của vật liệu khối do sự chi phối của hiệu ứng giam giữ lượng tử và tỷ lệ diện tích bề mặt trên thể tích lớn. Cụm nguyên tử của các nguyên tố không từ tính ở dạng khối (Ti, V, Cr) có thể mang từ tính; cụm nguyên tử của nguyên tố bán dẫn (Si, Ge) có thể mang đặc tính kim loại; cụm kim loại (Ga) thể hiện tính siêu dẫn; và cụm kim loại quý (Au) thể hiện hoạt tính xúc tác cao. Việc thêm hoặc bớt chỉ một nguyên tử có thể làm biến đổi cấu trúc hình học và cấu trúc điện tử của toàn bộ hệ.

Sự tham gia của các nguyên tử kim loại chuyển tiếp (với các điện tử chưa ghép cặp trên phân lớp d) vào cấu trúc cụm nguyên tử đóng vai trò quan trọng trong việc điều chỉnh độ bền, mô men từ, đặc tính hấp thụ quang và hoạt tính xúc tác. Tuy nhiên, việc nghiên cứu các hệ này bằng các phương pháp tính toán lượng tử gặp nhiều khó khăn do sự xuất hiện của các điện tử phân lớp d dẫn đến số lượng lớn đồng phân cấu trúc và cấu hình spin có mức năng lượng gần nhau. Quy luật tương tác giữa các điện tử định xứ d của kim loại chuyển tiếp và các điện tử tự do thuộc lớp vỏ điện tử của cụm nguyên tử, cũng như cơ chế phân ly và sự dập tắt từ tính ở cấp độ nguyên tử, vẫn còn nhiều vấn đề cần làm sáng tỏ.

Luận án tiến sĩ của tác giả Nguyễn Thị Mai với đề tài "Nghiên cứu cấu trúc hình học và điện tử của cụm nguyên tử chứa kim loại chuyển tiếp bằng phương pháp tính toán lượng tử định hướng ứng dụng cho vật liệu xúc tác" (chuyên ngành Vật liệu điện tử, mã số 9440123) được thực hiện nhằm giải quyết các khoảng trống nghiên cứu nêu trên.

Mục tiêu nghiên cứu của luận án được xác định gồm:

  1. Làm rõ cấu trúc hình học bền, xây dựng cấu trúc điện tử, cơ chế phân ly, quy luật tồn tại và biến đổi điện tử hóa trị của các cụm nguyên tử chứa kim loại chuyển tiếp thuộc ba hệ đại diện: cụm silicon - kim loại chuyển tiếp, cụm oxide kim loại chuyển tiếp, và cụm kim loại quý - kim loại chuyển tiếp.
  2. Dựa trên các kết quả cấu trúc và điện tử, thiết kế và đề xuất các cụm nguyên tử bền vững mang tính chất đặc biệt như những "siêu nguyên tử" (superatoms) phục vụ định hướng phát triển vật liệu nano.

Các câu hỏi và giả thuyết nghiên cứu chính được đặt ra:

  • Quy luật tiến hóa cấu trúc hình học và trật tự sắp xếp không gian của cụm nguyên tử thay đổi như thế nào khi số lượng nguyên tử $n$ tăng dần?
  • Sự tương tác và lai hóa giữa các điện tử định xứ 3d của kim loại chuyển tiếp với các điện tử hóa trị tự do diễn ra theo cơ chế nào, và ảnh hưởng ra sao đến tính bền vững, vùng cấm năng lượng HOMO-LUMO và mô men từ spin?
  • Kênh phân ly nhiệt động lực học nào chiếm ưu thế khi các cụm nguyên tử bị kích thích bởi bức xạ quang phân ly?
  • Có hay không sự hình thành và điều kiện dập tắt từ tính cục bộ (mô hình hiệu ứng Kondo ở cấp độ phân tử) trong các cụm kim loại quý pha tạp kim loại chuyển tiếp?

Đối tượng và phạm vi nghiên cứu của luận án bao gồm ba nhóm hệ cụm nguyên tử cô lập trong pha khí:

  • Cụm nguyên tử của kim loại chuyển tiếp và silicon: $\text{Si}_n\text{Mn}^{2+}$ ($n = 1-10$).
  • Cụm nguyên tử oxide của kim loại chuyển tiếp: $\text{Co}_n\text{O}_m^+$ ($n = 5-9, n-2 \le m \le n+2$).
  • Cụm nguyên tử của kim loại chuyển tiếp và kim loại quý: $\text{Au}_{19}\text{TM}$ ($\text{TM} = \text{Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu}$) và $\text{M}_n\text{Cr}$ ($\text{M} = \text{Cu, Ag, Au}$; $n = 2-20$).

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

Lịch sử nghiên cứu cụm nguyên tử bắt đầu từ các thí nghiệm tạo cụm trong pha khí của Robbins và cộng sự vào năm 1967. Đến năm 1984, nhóm nghiên cứu của Knight và cộng sự đã quan sát thấy sự biến đổi tuần hoàn về độ bền của cụm kim loại kiềm $\text{Na}_n$ tại các kích thước $n = 8, 20, 40, 58$, tương ứng với các lớp vỏ điện tử điền đầy theo mô hình Jellium và mô hình Clemenger-Nilsson. Đối với kim loại chuyển tiếp (như Fe, Ti), độ bền của cụm nguyên tử được chứng minh tuân theo quy luật vỏ hình học xếp chặt đối xứng (pentagonal bipyramid tại $n = 7$, icosahedron tại $n = 13$, BCC tại $n = 15$, icosahedron kép tại $n = 19$).

Trong lĩnh vực xúc tác, các công trình từ những năm 1980 về hạt nano vàng và cụm $\text{Au}_{16-35}$ của Haruta và các nhóm nghiên cứu đã chứng minh hoạt tính xúc tác cao trong phản ứng oxi hóa CO. Lei và cộng sự khảo sát động học quá trình epoxy hóa propylene ở nhiệt độ dưới 200°C bằng xúc tác cụm bạc $\text{Ag}_n$, ghi nhận hoạt tính nổi trội tại $\text{Ag}_3$. Các nghiên cứu tương tác pha khí của $\text{Pt}n^+$ và $\text{NbPt}{n-1}^+$ với CO giải thích hiện tượng giảm năng lượng hấp phụ CO khi pha tạp Nb tại $n = 18-22$. Đối với quá trình hấp phụ và phân ly $\text{H}_2$, các cụm $\text{Al}_n\text{V}^+$ ($n = 2, 4, 10$) và $\text{Al}_n\text{Rh}^{2+}$ ($n = 10-13$) được xác định rõ về rào thế kích hoạt và động lực học phản ứng.

Về từ tính và quang tính, các nghiên cứu lý thuyết và thực nghiệm xác định tính chất từ phụ thuộc vào cấu trúc hình học (ví dụ cụm phẳng $\text{Li}4$ có mô men từ bằng 0 trong khi đồng phân tứ diện có mô men từ $4,\mu\text{B}$); các nguyên tố không từ tính như V, Rh, Pd ở dạng cụm thể hiện từ tính rõ rệt; mô men từ trên mỗi nguyên tử trong cụm $\text{Co}_N\text{Mn}_M$ tăng theo nồng độ Mn; mô men từ của cụm Fe, Ni, Co tăng khi kích thước cụm giảm. Sự pha tạp Pd vào cụm $\text{Au}_n^+$ làm thay đổi cấu trúc từ dạng phẳng sang 3D và gây dập tắt hấp thụ quang học do giảm mật độ trạng thái s gần mức Fermi.

Đối với cụm silicon pha tạp kim loại chuyển tiếp, các kỹ thuật thực nghiệm như phổ hồng ngoại phân ly đa photon (IRMPD), quang phổ phân ly chọn lọc theo khối lượng (MSPFS), quang phổ lưỡng sắc từ trường tia X (XMCD) và phổ quang điện tử (PES) kết hợp với tính toán DFT đã được áp dụng. Bin Yang và cộng sự nghiên cứu $\text{Cr}_2\text{Si}_n^-$ ($n = 1-12$), chỉ ra cấu trúc ngoại diện tại $n = 1-7$ và cấu trúc lồng lăng kính phản lục giác tại $n = 6, 12$. Cụm $\text{Si}_6\text{M}^+$ ($\text{M} = \text{Cu, Ag, Cr}$) cho thấy độ bền cao trong phổ quang phân ly laser; cụm $\text{Si}_n\text{Co}^{2+}$, $\text{Si}_n\text{Cu}^+$, $\text{Si}_n\text{V}^+$, $\text{Si}_n\text{Nb}^+$ được xác định vị trí chiếm chỗ ưu tiên của nguyên tử pha tạp trên khung silicon.

Tại Việt Nam, các hướng nghiên cứu tính toán mô phỏng cụm nguyên tử được phát triển bởi nhóm của Ngô Tuấn Cường (cụm bạc và tính chất quang), Vũ Thị Ngân (tính chất phản ứng cụm silicon), nhóm tại Đại học Tôn Đức Thắng và Viện Khoa học và Công nghệ tính toán TP.HCM (cụm boron, nhôm), Trần Văn Tân (cụm Si và Ge pha tạp kim loại chuyển tiếp), Phạm Vũ Nhật (cụm vàng và phổ hồng ngoại). Hướng nghiên cứu cụm kim loại chuyển tiếp kết hợp thực nghiệm đo độ lệch từ Stern-Gerlach và quang phân ly tại Đại học KU Leuven (Bỉ) được triển khai tại Viện Khoa học vật liệu từ năm 2014.

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

  • Thiết lập quy trình tính toán phiếm hàm mật độ chuẩn xác cho các hệ cụm nhị nguyên chứa kim loại chuyển tiếp có độ suy biến spin cao.
  • Làm rõ bản chất tương tác giữa điện tử định xứ 3d và điện tử tự do trong các cấu trúc cụm silicon hóa trị cao ($\text{Si}_n\text{Mn}^{2+}$), cụm oxide kim loại chuyển tiếp ($\text{Co}_n\text{O}m^+$) và cụm kim loại quý pha tạp ($\text{Au}{19}\text{TM}$, $\text{M}_n\text{Cr}$).
  • Dự đoán và xác lập mô hình siêu nguyên tử mang từ tính có cấu hình vỏ điện tử đóng kín.

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

Luận án sử dụng phương pháp tính toán hóa học lượng tử trên nền tảng cơ học lượng tử và lý thuyết phiếm hàm mật độ (Density Functional Theory - DFT).

Khung lý thuyết và mô hình tính toán

  • Phương trình Schrödinger trạng thái dừng và gần đúng Born-Oppenheimer: Tách biệt chuyển động của hạt nhân và điện tử, cho phép giải phương trình vi phân cho hệ nhiều điện tử trong trường thế hạt nhân cố định: $$\hat{H}\Psi = E\Psi$$
  • Phương pháp tính toán Ab-initio và DFT: Sử dụng gần đúng Hartree-Fock (HF), tương tác cấu hình (CI), và lý thuyết phiếm hàm mật độ với các gần đúng mật độ cục bộ (LDA, LSDA) và gần đúng gradient tổng quát (GGA).
  • Mô hình vỏ điện tử và vỏ hình học:
    • Mô hình Jellium: Coi các lõi ion dương tạo thành một giếng thế liên tục, các điện tử hóa trị chuyển động tự do điền đầy các mức năng lượng theo thứ tự: $1\text{S}^2 1\text{P}^6 1\text{D}^{10} 2\text{S}^2 1\text{F}^{14} 2\text{P}^6 1\text{G}^{18} 2\text{D}^{10} 1\text{H}^{22} 3\text{S}^2 2\text{F}^{14} 1\text{I}^{26} 3\text{P}^6 \dots$ Các cụm có số điện tử 2, 8, 20, 40, 58... đạt trạng thái vỏ đóng kín ("siêu nguyên tử").
    • Mô hình Clemenger-Nilsson: Áp dụng toán tử biến dạng elipsoid cho các cụm có lớp vỏ chưa điền đầy nhằm giải thích sự giảm tính suy biến năng lượng theo định lý Jahn-Teller.
    • Mô hình vỏ hình học đối xứng: Áp dụng cho các cụm kim loại chuyển tiếp với các cấu trúc hình học bền vững xếp chặt ($n = 7, 13, 15, 19$).

Các thông số đánh giá độ bền và tính chất điện tử

  1. Năng lượng liên kết trung bình trên mỗi nguyên tử ($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 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 ($DE$): Tính bằng hiệu số giữa tổng năng lượng các tiểu phân sản phẩm và năng lượng của cụm mẹ, có hiệu chỉnh năng lượng dao động điểm không (ZPE).
  4. Tham số nhiệt động lực học: Biến thiên năng lượng tự do Gibbs ($\Delta G$), biến thiên entropy ($\Delta S$), biến thiên enthalpy ($\Delta H$) ở nhiệt độ xác định (300 K).
  5. Cấu trúc điện tử và từ tính: Bề rộng vùng cấm HOMO-LUMO ($E_g$), khoảng cách năng lượng $\delta_1$ (LUMO $\beta$ - HOMO $\alpha$) và $\delta_2$ (LUMO $\alpha$ - HOMO $\beta$), mật độ trạng thái điện tử (DOS), phân tích orbital phân tử (MO), phân bố mật độ spin và điện tích tự nhiên.

Nguồn dữ liệu và quy trình tính toán

  • Tạo tập hợp các cấu trúc đồng phân khởi tạo ban đầu thông qua các thuật toán tối ưu hóa hình học.
  • Khảo sát các bộ phiếm hàm (B3LYP, PBE, BP86, TPSSh, M06, v.v.) và hệ hàm cơ sở (cc-pVTZ, def2-TZVP, LanL2DZ, v.v.) kết hợp thế giả định (pseudopotential) cho các nguyên tử kim loại nặng.
  • Đối sánh kết quả tính toán năng lượng phân ly và tần số dao động với dữ liệu thực nghiệm quang phân ly chọn lọc khối lượng (MSPFS), phổ hồng ngoại và phổ khối thời gian bay (TOF-MS).
  • Tính toán được thực hiện trên các hệ thống máy tính hiệu năng cao tại Trung tâm Tin học và Tính toán (Viện Hàn lâm KHCN Việt Nam), Viện Khoa học vật liệu, Trường Đại học Sư phạm Hà Nội và Trường Đại học Tôn Đức Thắng.
Hệ nghiên cứu Đối tượng cụ thể Phương pháp tính toán chính Dữ liệu thực nghiệm đối sánh
Silicon - Kim loại chuyển tiếp $\text{Si}_n\text{Mn}^{2+}$ ($n = 1-10$) DFT / Tối ưu hóa đa cấu hình spin Phổ khối TOF, số liệu quang phân ly
Oxide kim loại chuyển tiếp $\text{Co}_n\text{O}_m^+$ ($n = 5-9$) DFT / Nhiệt động lực học ($\Delta G, \Delta S$) Phổ quang phân ly MSPFS (laser 355 nm)
Kim loại quý - Kim loại chuyển tiếp $\text{Au}_{19}\text{TM}$, $\text{M}_n\text{Cr}$ ($n = 2-20$) DFT / Phân tích MO, DOS, Spin Hiệu ứng Kondo phân tử, phổ PES

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

Chương 1. Tổng quan về cụm nguyên tử kim loại chuyển tiếp

Chương 1 hệ thống hóa cơ sở lý thuyết và các nghiên cứu về cụm nguyên tử. Tác giả phân tích sự biến đổi tính chất theo kích thước từ nguyên tử đơn lẻ đến vật liệu khối; bản chất mô hình vỏ điện tử Jellium và mô hình biến dạng Clemenger-Nilsson đối với cụm kim loại kiềm và kim loại quý; mô hình vỏ hình học đối với cụm kim loại chuyển tiếp.

Chương này tổng kết các nghiên cứu về hoạt tính xúc tác của cụm kim loại (Au, Ag, Pt, Pd, Al, Rh), vai trò của các điện tử chưa ghép cặp trong phản ứng oxi hóa CO, phân ly $\text{H}_2$ và epoxy hóa propylene. Đồng thời, chương 1 phân tích các đặc trưng từ tính (siêu thuận từ, tăng mô men từ ở kích thước nhỏ, sự xuất hiện từ tính ở cụm phi từ) và tính chất quang của cụm kim loại quý pha tạp. Cuối cùng, tác giả tổng quan tình hình nghiên cứu các hệ cụm nhị nguyên gồm silicon - kim loại chuyển tiếp, oxide kim loại chuyển tiếp và kim loại quý - kim loại chuyển tiếp trên thế giới và tại Việt Nam.

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

Chương 2 trình bày chi tiết cơ sở toán học và vật lý của các phương pháp tính toán lượng tử. Nội dung bao gồm phương trình Schrödinger, phép gần đúng Born-Oppenheimer, phương pháp Hartree-Fock và tương tác cấu hình (CI), lý thuyết phiếm hàm mật độ (DFT) với các hàm gần đúng LDA, LSDA và GGA.

Quy trình tính toán tối ưu hóa năng lượng của các cụm nguyên tử được thiết lập gồm các bước: xây dựng không gian cấu trúc ban đầu, quét trạng thái spin, tối ưu hóa hình học, kiểm tra tần số dao động để đảm bảo cấu trúc điểm cực tiểu năng lượng (không có tần số ảo), tính toán năng lượng hiệu chỉnh ZPE, năng lượng liên kết $BE$, chênh lệch bậc hai $\Delta^2 E$, năng lượng phân ly $DE$, và các hàm trạng thái nhiệt động lực học ($\Delta H, \Delta S, \Delta G$).

Chương 3. Cụm nguyên tử của kim loại chuyển tiếp và silicon SinMn2+

Chương 3 trình bày các kết quả nghiên cứu hệ cụm nguyên tử $\text{Si}_n\text{Mn}^{2+}$ ($n = 1-10$):

  • Khảo sát phiếm hàm và bộ hàm cơ sở: Đánh giá độ bội, độ dài liên kết và năng lượng phân ly của các cụm $\text{SiTM}$ ($\text{TM} = \text{Sc, V, Cr, Mn, Fe, Cu}$) ở các mức lý thuyết khác nhau nhằm chọn lựa phiếm hàm phù hợp.
  • Tiến hóa cấu trúc hình học: Với $n = 1-5$, cụm nguyên tử có xu hướng tạo thành các cấu trúc mở dạng chóp phẳng hoặc bán lồng; từ $n = 6-10$, cấu trúc chuyển dịch sang dạng lồng không gian 3D bọc lấy nguyên tử Mn ở tâm hoặc phối trí mặt có số liên kết cực đại.
  • Cấu trúc điện tử và phân bố điện tích: Phân tích giản đồ orbital phân tử cho thấy sự chuyển giao điện tích từ khung silicon sang nguyên tử mangan và trạng thái ion hóa $2+$; xác định mật độ spin và mô men từ của cụm.
  • Độ bền và năng lượng phân ly: Đồ thị năng lượng liên kết $BE$ tăng dần theo kích thước cụm; chênh lệch năng lượng bậc hai $\Delta^2 E$ chỉ ra các kích thước có độ bền tương đối vượt trội; tính toán định lượng 8 kênh phân ly khả dĩ (1)-(8), chỉ rõ các kênh giải phóng nguyên tử hoặc ion trung hòa chiếm ưu thế về mặt năng lượng.

Chương 4. Cụm nguyên tử oxide của kim loại chuyển tiếp ConOm+

Chương 4 tập trung vào nghiên cứu thực nghiệm quang phân ly kết hợp tính toán lượng tử cho hệ cụm $\text{Co}_n\text{O}_m^+$ ($n = 5-9, n-2 \le m \le n+2$):

  • Kết quả thực nghiệm quang phân ly: Sử dụng phổ quang phân ly chọn lọc khối lượng (MSPFS) với nguồn laser bước sóng 355 nm, mật độ năng lượng $80,\text{mJ/cm}^2$. Phổ khối ghi nhận các kênh phân ly chính của cụm $\text{Co}_n\text{O}_m^+$, đặc biệt là các cụm thiếu oxy dạng $\text{Co}n\text{O}{n-2}^+$.
  • Cấu trúc hình học và cấu hình spin: Khảo sát các mức lý thuyết trên cụm $\text{Co}_5\text{O}_4^+$ và $\text{Co}_6\text{O}_4^+$, xác định các đồng phân có năng lượng thấp nhất, độ bội spin và vùng cấm HOMO-LUMO.
  • Năng lượng phân ly và biến thiên nhiệt động lực học: Tính toán năng lượng phân ly $DE$, biến thiên entropy $\Delta S$ trên một đơn vị nhiệt độ ($\text{eV/K}$) và biến thiên năng lượng tự do Gibbs $\Delta G$ tại 300 K cho các quá trình bay hơi một nguyên tử O, phân tử $\text{O}_2$, nguyên tử Co và phân tử $\text{CoO}$. Kết quả tính toán giải thích sự phù hợp về mặt nhiệt động học của các kênh phân ly bay hơi $\text{CoO}$ và $\text{O}_2$ được quan sát trực tiếp trong thực nghiệm.

Chương 5. Cụm nguyên tử của kim loại chuyển tiếp và kim loại quý Au19TM, MnCr

Chương 5 khảo sát hai hệ cụm nguyên tử kim loại chuyển tiếp kết hợp kim loại quý:

  • Hệ cụm $\text{Au}_{19}\text{TM}$ ($\text{TM} = \text{Sc-Cu}$):
    • Khảo sát cấu trúc hình học bền vững nhất của khung $\text{Au}_{19}$ khi pha tạp một nguyên tử kim loại chuyển tiếp 3d.
    • Phân tích cấu trúc điện tử, giản đồ MO, mật độ trạng thái (DOS) và phân bố mật độ spin.
    • Xác định cụm $\text{Au}_{19}\text{Cr}$ có cấu trúc hình học đối xứng cao, đạt trạng thái vỏ điện tử đóng kín với 20 điện tử hóa trị ($1\text{S}^2 1\text{P}^6 2\text{S}^2 1\text{D}^{10}$), đồng thời duy trì 5 điện tử $3\text{d}$ chưa ghép cặp định xứ trên nguyên tử Cr. Cụm này thể hiện đặc tính của một "siêu nguyên tử" mang từ tính cao và có độ bền năng lượng nổi trội.
  • Hệ cụm $\text{M}_n\text{Cr}$ ($\text{M} = \text{Cu, Ag, Au}$; $n = 2-20$):
    • Khảo sát sự tiến hóa cấu trúc hình học, sự thay đổi độ dài liên kết $\text{Cr-M}$ và năng lượng phân ly của các cụm $\text{CrM}$ ở các mức lý thuyết khác nhau.
    • Xác định cấu hình điện tử trên orbital $3\text{d-Cr}$, độ bội spin và giản đồ orbital phân tử qua từng kích thước ($n = 2-12$ và $n = 13-20$).
    • Phân tích giá trị vùng cấm HOMO-LUMO ($E_g$) cùng các thông số $\delta_1$ và $\delta_2$, năng lượng liên kết $BE$, chênh lệch năng lượng bậc hai $\Delta^2 E$, và năng lượng phân ly theo các kênh. Kết quả làm sáng tỏ điều kiện xuất hiện và cơ chế dập tắt từ tính của tạp chất từ Cr trong đám mây điện tử tự do của kim loại quý.

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

Luận án đạt được các kết quả và đóng góp mới về lý luận, phương pháp và thực tiễn:

Đóng góp về lý luận và khoa học

  1. Xác định một cách hệ thống cấu trúc hình học trạng thái cơ bản, cấu trúc điện tử, phân bố mật độ spin và quy luật tiến hóa theo kích thước của ba nhóm cụm nguyên tử: $\text{Si}_n\text{Mn}^{2+}$ ($n = 1-10$), $\text{Co}_n\text{O}m^+$ ($n = 5-9$), $\text{Au}{19}\text{TM}$ ($\text{TM} = \text{Sc-Cu}$) và $\text{M}_n\text{Cr}$ ($\text{M} = \text{Cu, Ag, Au}$; $n = 2-20$).
  2. Làm rõ bản chất tương tác điện tử giữa các điện tử hóa trị tự do thuộc lớp vỏ chung của cụm và các điện tử định xứ chưa ghép cặp thuộc phân lớp $3\text{d}$ của nguyên tử kim loại chuyển tiếp.
  3. Xác lập mối tương quan giữa cấu trúc hình học đối xứng xếp chặt, số lượng điện tử hóa trị và quy luật điền đầy lớp vỏ điện tử (Jellium) trong các hệ cụm nhị nguyên chứa kim loại chuyển tiếp.
  4. Chứng minh cụm nguyên tử $\text{Au}_{19}\text{Cr}$ là một "siêu nguyên tử" bền vững, có cấu trúc đối xứng cao, sở hữu lớp vỏ điện tử đóng kín 20 điện tử ($1\text{S}^2 1\text{P}^6 2\text{S}^2 1\text{D}^{10}$) và đồng tồn tại 5 điện tử $3\text{d}$ chưa ghép cặp định xứ trên nguyên tử Cr tạo ra mô men từ lớn.
  5. Làm sáng tỏ cơ chế phân ly nhiệt động lực học của các cụm oxide $\text{Co}_n\text{O}_m^+$, chứng minh sự tương hợp giữa kết quả tính toán năng lượng phân ly, biến thiên năng lượng tự do Gibbs $\Delta G$ và số liệu quang phân ly thực nghiệm (MSPFS).

Đóng góp về thực tiễn và ứng dụng

  1. Thiết lập quy trình tính toán mô phỏng lượng tử tin cậy đối với các hệ cụm nguyên tử kim loại chuyển tiếp phức tạp có độ suy biến spin cao.
  2. Cung cấp dữ liệu cơ bản về cấu trúc, năng lượng liên kết, vùng cấm HOMO-LUMO và các kênh phân ly ưu tiên, làm cơ sở dự đoán tính chất hóa lý trước khi tiến hành tổng hợp thực nghiệm.
  3. Đề xuất các cấu trúc cụm nguyên tử bền vững có hoạt tính từ và xúc tác tiềm năng, định hướng làm các đơn vị cấu trúc nền tảng (building blocks) để chế tạo vật liệu nano tiên tiến.

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

Hạn chế

  • Các mô hình tính toán và dữ liệu thực nghiệm trong luận án chủ yếu thực hiện trên các cụm nguyên tử cô lập trong pha khí ở điều kiện chân không, chưa tính đến ảnh hưởng tương tác của môi trường xung quanh như chất mang (substrate) rắn hoặc dung môi phân tán.
  • Chưa khảo sát chi tiết đường phản ứng động học hoàn chỉnh (năng lượng trạng thái chuyển tiếp) của các chu trình xúc tác cụ thể trên bề mặt các cụm nguyên tử được đề xuất.

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

  1. Mở rộng tính toán mô phỏng tương tác và hoạt tính xúc tác trực tiếp của các "siêu nguyên tử" (như $\text{Au}_{19}\text{Cr}$) với các phân tử phản ứng nhỏ ($\text{CO}, \text{O}_2, \text{H}_2, \text{CO}_2$).
  2. Nghiên cứu ảnh hưởng của các bề mặt giá mang (như graphene, $\text{TiO}_2, \text{Al}_2\text{O}_3$) lên cấu trúc hình học, cấu trúc điện tử và độ bền của các cụm nguyên tử.
  3. Phối hợp các phương pháp tổng hợp hóa học ướt hoặc lắng đọng chùm cụm nguyên tử (cluster beam deposition) để chế tạo và kiểm chứng tính chất của các cấu trúc cụm bền vững trong điều kiện thực tế.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên khảo cho các đối tượng:

  • Nghiên cứu sinh, học viên cao học và sinh viên thuộc các chuyên ngành Khoa học vật liệu, Vật liệu điện tử, Hóa học lý thuyết và Hóa lý, Vật lý chất rắn, Công nghệ nano.
  • Giảng viên và nhà nghiên cứu trong lĩnh vực tính toán lượng tử cấu trúc vật liệu, mô phỏng DFT cho hệ có cấu hình spin phức tạp, xúc tác dị thể và vật liệu từ tính dưới nanomet.

Các phần có giá trị tham khảo chuyên sâu:

  • Chương 2: Quy trình thiết lập và tối ưu hóa đa cấu hình spin cho các hệ chứa kim loại chuyển tiếp d.
  • Chương 4: Phương pháp tính toán các hàm trạng thái nhiệt động lực học ($\Delta G, \Delta S, DE$) để đối sánh trực tiếp với dữ liệu khối phổ quang phân ly MSPFS.
  • Chương 5: Phương pháp phân tích orbital phân tử (MO) và mật độ trạng thái (DOS) để xác định cấu trúc "siêu nguyên tử" và tương tác giữa điện tử tự do với tạp chất từ.

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

1. Luận án nghiên cứu những nhóm đối tượng cụ thể nào?
Luận án tập trung vào ba nhóm cụm nguyên tử chứa kim loại chuyển tiếp trong pha khí:

  • Cụm silicon và kim loại chuyển tiếp: $\text{Si}_n\text{Mn}^{2+}$ ($n = 1-10$).
  • Cụm oxide kim loại chuyển tiếp: $\text{Co}_n\text{O}_m^+$ ($n = 5-9, n-2 \le m \le n+2$).
  • Cụm kim loại chuyển tiếp và kim loại quý: $\text{Au}_{19}\text{TM}$ ($\text{TM} = \text{Sc-Cu}$) và $\text{M}_n\text{Cr}$ ($\text{M} = \text{Cu, Ag, Au}$; $n = 2-20$).

2. Tại sao cụm nguyên tử $\text{Au}_{19}\text{Cr}$ lại được xem là một "siêu nguyên tử" bền vững?
Cụm $\text{Au}{19}\text{Cr}$ có cấu trúc hình học đối xứng cao, đạt cấu hình vỏ điện tử đóng kín với 20 điện tử hóa trị điền đầy theo mô hình Jellium ($1\text{S}^2 1\text{P}^6 2\text{S}^2 1\text{D}^{10}$), mang lại độ bền hóa học và năng lượng liên kết cao. Đồng thời, cụm duy trì 5 điện tử $3\text{d}$ chưa ghép cặp định xứ trên nguyên tử Cr, tạo ra mô men từ lớn ($5,\mu\text{B}$) và hoạt tính tương tác bề mặt mạnh.

3. Phương pháp thực nghiệm nào được tác giả sử dụng để đối sánh với kết quả tính toán DFT đối với cụm cobalt oxide?
Tác giả sử dụng kỹ thuật phổ quang phân ly chọn lọc khối lượng (Mass-Selected Photofragmentation Spectroscopy - MSPFS) kết hợp nguồn laser xung bước sóng 355 nm (mật độ năng lượng $80,\text{mJ/cm}^2$) và khối phổ kế thời gian bay (TOF-MS) để xác định các mảnh ion phân ly và kênh bay hơi ưu tiên của cụm $\text{Co}_n\text{O}_m^+$.

4. Quy luật tiến hóa cấu trúc hình học của cụm $\text{Si}_n\text{Mn}^{2+}$ biến đổi ra sao theo kích thước $n = 1-10$?
Ở các kích thước nhỏ ($n = 1-5$), cụm có cấu trúc dạng chóp hoặc bán phẳng. Khi số lượng nguyên tử silicon tăng lên ($n = 6-10$), cấu trúc phát triển thành dạng lồng không gian ba chiều bao quanh hoặc phối trí tối đa xung quanh nguyên tử Mn ở tâm, giúp tối ưu hóa số lượng liên kết $\text{Si-Si}$ và $\text{Si-Mn}$.

5. Những chỉ số nào được dùng để xác định độ bền của cụm nguyên tử trong luận án?
Độ bền tĩnh được đánh giá qua năng lượng liên kết trung bình trên mỗi nguyên tử ($BE$), chênh lệch năng lượng bậc hai ($\Delta^2 E$), và bề rộng vùng cấm HOMO-LUMO. Độ bền nhiệt động lực học được đánh giá thông qua năng lượng phân ly ($DE$), biến thiên entropy ($\Delta S$), biến thiên enthalpy ($\Delta H$) và biến thiên năng lượng tự do Gibbs ($\Delta G$) ở 300 K.

Kết luận

Luận án tiến sĩ của tác giả Nguyễn Thị Mai đã hoàn thành việc nghiên cứu có hệ thống cấu trúc hình học, cấu trúc điện tử, độ bền nhiệt động và cơ chế phân ly của ba hệ cụm nguyên tử chứa kim loại chuyển tiếp ($\text{Si}_n\text{Mn}^{2+}$, $\text{Co}_n\text{O}m^+$, $\text{Au}{19}\text{TM}$ và $\text{M}n\text{Cr}$) bằng lý thuyết phiếm hàm mật độ kết hợp đối sánh dữ liệu thực nghiệm quang phân ly. Nghiên cứu đã làm sáng tỏ cơ chế tương tác giữa các điện tử định xứ $3\text{d}$ và lớp vỏ điện tử tự do, đồng thời thiết kế thành công mô hình "siêu nguyên tử" mang từ tính $\text{Au}{19}\text{Cr}$. Các kết quả đạt được cung cấp cơ sở dữ liệu và quy trình mô phỏng tin cậy, đóng góp vào sự phát triển của lĩnh vực vật liệu điện tử và xúc tác cấu trúc nano.