Tổng quan luận án

Nghiên cứu về vật liệu nano dạng cluster (cụm nguyên tử) là một lĩnh vực thuộc khoa học vật liệu và hóa học lượng tử. Cluster là tập hợp từ một vài đến hàng ngàn nguyên tử ở kích thước nanomet hoặc nhỏ hơn, sở hữu các đặc tính vật lý và hóa học khác biệt so với vật liệu khối tương ứng. Xu hướng thu nhỏ kích thước của các linh kiện điện tử đặt ra yêu cầu tìm kiếm các vật liệu bán dẫn kích cỡ nano mới dựa trên các khối đơn vị cấu trúc silicon. Tuy nhiên, khác với nguyên tử carbon lai hóa $sp^2$ tạo nên cấu trúc lồng fullerene $C_{60}$ bền vững, nguyên tử silicon có xu hướng tạo liên kết lai hóa $sp^3$ dạng tứ diện, khiến cấu trúc lồng của cluster silicon thuần khiết ($Si_n$) kém ổn định về mặt nhiệt động và khó ứng dụng thực tế. Nhằm khắc phục hạn chế này, hướng nghiên cứu pha tạp các nguyên tố kim loại (đặc biệt là kim loại chuyển tiếp dãy 3d) vào khung silicon được triển khai nhằm tìm kiếm các cấu trúc bền vững hơn và có đặc tính electron thích hợp.

Khoảng trống nghiên cứu được tác giả xác định xuất phát từ hai khía cạnh:

  • Về mặt thực nghiệm: Các kỹ thuật đo đạc hiện nay (như quang phổ photoelectron, quang phổ phân ly đa hồng ngoại IR-MPD, phổ khối lượng) chưa thể cung cấp thông tin trực tiếp và đầy đủ về cấu trúc hình học chi tiết cũng như bản chất liên kết bên trong cluster.
  • Về mặt lý thuyết tính toán: Cluster có cấu trúc rất linh động với nhiều đồng phân có mức năng lượng xấp xỉ nhau tồn tại đồng thời. Các nghiên cứu trước đây phần lớn tập trung vào quy luật hình thành hình học hoặc ảnh hưởng của điện tích ở các kích thước riêng lẻ mà chưa phân tích sâu bản chất liên kết hóa học (sự xen phủ orbital, tương tác Si–kim loại và kim loại–kim loại) cũng như sự cạnh tranh khi pha tạp đồng thời hai nguyên tử kim loại khác nhau.

Mục tiêu nghiên cứu của luận án bao gồm:

  1. Xác định cấu trúc hình học của các đồng phân bền nhất, nhóm điểm đối xứng và trạng thái electron của các dãy cluster silicon pha tạp kim loại ở các trạng thái điện tích trung hòa, cation và anion.
  2. Khảo sát các thông số năng lượng định lượng: năng lượng liên kết trung bình ($E_b$), năng lượng phân ly liên kết ($D$), biến thiên năng lượng bậc hai ($\Delta^2 E$), năng lượng ion hóa đoạn nhiệt ($AIE$), ái lực electron ($EA$) và độ rộng vùng cấm năng lượng ($E_{gap} = E_{LUMO} - E_{HOMO}$).
  3. Làm sáng tỏ quy luật hình thành, phát triển cấu trúc và bản chất liên kết hóa học thông qua phân tích orbital phân tử (MO), phân tích orbital liên kết tự nhiên (NBO), cấu hình electron tự nhiên (NEC) và lý thuyết cộng hưởng tự nhiên (NRT).

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

  • Đối tượng nghiên cứu: Các hệ cluster silicon thuần khiết và cluster silicon pha tạp một hoặc hai nguyên tử kim loại chuyển tiếp dãy 3d (Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn) và một số kim loại nhóm chính (Li, Na, K).
  • Phạm vi khảo sát:
    • Dãy cluster silicon pha tạp một nguyên tử kim loại thay đổi kích thước: $Si_nSc^{-/0/+}$ ($n = 2-11$), $Si_nFe^{-/0/+}$ ($n = 8-12$).
    • Dãy cluster silicon pha tạp một nguyên tử kim loại thay đổi nguyên tố: $Si_2M$ ($M = Li, Na, K, Cu, Cr$), $Si_3M$ ($M = Sc-Zn$).
    • Dãy cluster silicon pha tạp hai nguyên tử kim loại giống nhau: $Si_nTi_2$ ($n = 1-8$), $Si_2M_2$ ($M = Sc-Zn$), $Si_3M_2$ ($M = Sc-Zn$).
    • Dãy cluster silicon pha tạp hai nguyên tử kim loại khác nhau thay đổi kích thước: $Si_nScTi^{0/+/-}$ ($n = 1-9$).
  • Phương pháp luận: Hóa học tính toán dựa trên thuyết phiếm hàm mật độ (DFT).

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

Luận án đã hệ thống hóa lịch sử phát triển của khoa học cluster cùng các công trình nghiên cứu thực nghiệm và lý thuyết tiêu biểu:

  • Về lịch sử nghiên cứu cluster: Hợp chất calomel $Hg_2Cl_2$ được biết đến từ thế kỷ 12 tại Ấn Độ và được xác định liên kết thủy ngân vào năm 1930; ứng dụng cluster hạt vàng tạo màu cho kính; Rundle và Dahl phát hiện liên kết Mn-Mn trong phân tử $Mn_2(CO){10}$; Paolo Chini (thập niên 1970) tổng hợp các cluster kim loại lớn họ bạch kim như $[Rh{13}(CO){24}H_3]^{2-}$ qua nhiễu xạ tia X đơn tinh thể; Linus Pauling nghiên cứu cụm bát diện $Mo_6$ trong "$MoCl_2$"; Albert Cotton xác định cấu trúc $Re_3Cl_9$ và liên kết bậc bốn trong $K_2Re_2Cl_8$; Alfred Stock phát hiện các boron hydrides; các cluster sinh học $Fe_4S_4$ và nitrogenase $MoFe_7S_9$; các hợp chất Zintl anion ($[Bi_3]^{3-}$, $[Sn_9]^{4-}$, $[Pb_9]^{4-}$, $[Sb_7]^{3-}$, $Ni@Pb{10}^{2-}$, $[Ir@Sn_{12}]^{3-}$); Walter Knight và cộng sự (đầu những năm 1980) điều chế cụm kim loại kiềm lên đến 100 nguyên tử; các cluster kim loại có lõi ổn định $Al_{77}$, $Ga_{84}$, $Pd_{145}$; sự phát hiện fullerene $C_{60}, C_{70}, C_{84}$.
  • Về cluster silicon thuần khiết: Honig (1954) phát hiện các cluster silicon trung hòa từ $Si_2$ đến $Si_7$ bằng phương pháp hóa hơi silicon qua phổ khối lượng; Tsong (1984) dùng kỹ thuật giải hấp photostimulation xác định các cation $Si_n^+$ ($n = 3-16$); Gun'ko nghiên cứu silica gel Si-60 và pyrocarbon/silica gel bằng TEM, $^{1}H\text{-NMR}$ và hấp phụ khử nitơ ở 77,35 K; J. Nordlund (2003) mô phỏng sự lắng đọng cụm $Si_{20}$ trên bề mặt Si(001) ở 300 K và 1000 K tại các mức năng lượng 1 eV, 5 eV, 10 eV; Holthausen tổng hợp hợp chất khối mười hai $(Si_{20})^{12-}$ và silafullerane chứa ion chloride nội tiếp qua phân tích tinh thể tia X, quang phổ $^{29}Si\text{ NMR}$ và nhiễu xạ bột.
  • Về cluster silicon pha tạp: Tác giả điểm lại các công trình pha tạp kim loại kiềm ($Li, Na$), kim loại kiềm thổ ($Be, Mg, Ba$), phi kim ($B, C, N$), kim loại chuyển tiếp dãy 3d ($Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn$), dãy 4d ($Y, Zr, Nb, Mo, Ru, Rh, Pd, Ag$), dãy 5d ($Hf, Ta, W, Re, Os, Au$). Tiêu biểu gồm các nghiên cứu về $Si_nNi$ ($n = 1-17$), $Si_{12}M$ ($M = Sc-Zn$), $Mo_2Si_n$ ($n = 9-16$), $Si_nAu_2$ ($n = 1-10$), $T_2Si_n$ ($T = Cr, Mn, Fe, Co, Ni$), $Zr_2Si_n$, $Sc_2Si_n^-$, $V_2Si_n^-$, $V_3Si_n^-$, $Si_nAl_m$, $Si_nMg_m$.

Vị trí của luận án: Kế thừa các kết quả về hình học của các hệ cluster đơn lẻ, luận án mở rộng khảo sát có hệ thống từ cluster silicon pha tạp 1 nguyên tử sang 2 nguyên tử kim loại (cùng loại và khác loại), tập trung giải quyết mối liên hệ định lượng giữa cấu trúc hình học, độ bền nhiệt động, sự phân bố electron và bản chất liên kết hóa học.

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

Luận án trình bày hệ thống cơ sở lý thuyết hóa học lượng tử và phương pháp tính toán:

  • Các phương pháp gần đúng orbital phân tử (MO):
    • Phương pháp Hartree-Fock (HF): Dựa trên phương pháp trường tự hợp (SCF), sử dụng định thức Slater đơn mô tả hệ $N$ electron và phương trình Roothaan (MO-LCAO-SCF) cho hệ phân tử.
    • Phương pháp nhiễu loạn Moller-Plesset (MP2, MP3, MP4): Bổ sung hiệu chỉnh thế tương quan electron bậc hai, bậc ba và bậc bốn dựa trên lý thuyết nhiễu loạn Rayleigh-Schrodinger.
    • Phương pháp tương tác cấu hình (CI, CID, CISD, MRCI): Sử dụng tổ hợp tuyến tính các hàm trạng thái cấu hình (CSFs) để tính toán năng lượng tương quan electron và các trạng thái kích thích.
    • Phương pháp chùm tương tác (Coupled Cluster - CC / CPMET): Phát triển bởi Fritz Coester, Hermann Kümmel, Jirí Cizek và Josef Paldus, sử dụng toán tử cụm hàm mũ để tính toán chính xác năng lượng tương quan.
  • Thuyết phiếm hàm mật độ (DFT):
    • Mô hình Thomas-Fermi (1927) biểu diễn năng lượng hệ theo mật độ electron $\rho(r)$.
    • Hai định lý Hohenberg-Kohn (1964) chứng minh mật độ electron xác định trạng thái cơ bản và năng lượng hệ có tính chất biến phân.
    • Phương pháp Kohn-Sham (1965) quy bài toán nhiều hạt về hệ phương trình một electron với thế tương quan trao đổi $V_{XC}$.
    • Các phiếm hàm trao đổi (LDA, LSDA của Slater, B88, PW91) và phiếm hàm tương quan (LYP của Lee-Yang-Parr, PW91 của Perdew-Wang, VWN của Vosko-Wilk-Nusair).
    • Phương pháp DFT hỗn hợp: B3LYP và B3P86. Phương pháp B3P86 kết hợp phiếm hàm trao đổi 3 thông số Becke (B3), phiếm hàm tương quan P86 và 20% năng lượng trao đổi Hartree-Fock.
  • Phương pháp phân tích orbital liên kết tự nhiên (NBO) và lý thuyết cộng hưởng (NRT):
    • Phân tích mật độ electron dựa trên orbital tự nhiên (NO), orbital nguyên tử tự nhiên (NAO: lõi, hóa trị, Rydberg) và orbital liên kết tự nhiên (NBO: một lõi CR, cặp electron riêng LP, liên kết hai tâm BD, phản liên kết BD*, không liên kết LP*).
    • Thuyết cộng hưởng tự nhiên (NRT) định lượng sự dịch chuyển mật độ electron và đóng góp cấu trúc Lewis.
  • Thuyết lớp vỏ electron (Phenomenological Shell Model - PSM): Mô hình điền electron vào các lớp vỏ $1S, 1P, 1D, 2S, 1F, 2P, 1G, 2D...$ đối với các cluster dạng cầu.
  • Mức lý thuyết và phần mềm tính toán thực tế:
    • Mức lý thuyết: B3P86/6-311+G(d). Năng lượng tổng của các cluster được hiệu chỉnh năng lượng điểm không (ZPE).
    • Kiểm chứng độ tin cậy với thực nghiệm trên các phân tử đime:
      • Đime $Si_2$: Tính toán đạt độ dài liên kết $r = 2,27\text{ \AA}$, tần số dao động $\omega_e = 495,53\text{ cm}^{-1}$, năng lượng liên kết $E_b = 3,17\text{ eV}$, phù hợp với số liệu thực nghiệm ($r = 2,25\text{ \AA}$, $\omega_e = 509 \pm 10\text{ cm}^{-1}$, $E_b = 3,21\text{ eV}$).
      • Phân tử $SiSc$: Khoảng cách $r = 2,52\text{ \AA}$, bền nhất ở trạng thái spin quartet ($2S+1 = 4$).
      • Phân tử $SiTi$: Khoảng cách $r = 2,41\text{ \AA}$, bền nhất ở trạng thái spin quintet ($2S+1 = 5$).
    • Phần mềm sử dụng: Gaussian 09 (phiên bản E.01) để tối ưu hóa cấu trúc và tính toán năng lượng; Gaussview 05 hỗ trợ dựng hình và hiển thị; PyMolyze, gOpenMol, Dgrid-4 để phân tích cấu trúc electron và orbital; Molden, JMol, Corel Drawx3 để vẽ cấu trúc; NBO 5.G để phân tích liên kết; Excel và Origin để xử lý dữ liệu.

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

Cluster silicon pha tạp một nguyên tử kim loại

Nội dung khảo sát tập trung vào các hệ cluster $Si_nSc^{-/0/+}$ ($n = 2-11$), $Si_nFe^{-/0/+}$ ($n = 8-12$), $Si_2M$ ($M = Li, Na, K, Cu, Cr$) và $Si_3M$ ($M = Sc-Zn$).

Đối với hệ $Si_nSc^{-/0/+}$ ($n = 2-11$), quy luật hình thành cấu trúc được phân hóa theo trạng thái điện tích:

  • Dạng trung hòa ($Si_nSc$) và cation ($Si_nSc^+$): Ưu tiên trạng thái spin thấp (singlet, doublet). Khi $n \le 6$, cấu trúc hình thành theo quy luật thế một nguyên tử Si trong khung $Si_{n+1}$ bằng nguyên tử Sc, tương ứng với dạng tam giác ($Si_3$), hình thoi ($Si_4$), lưỡng tháp tam giác ($Si_5$), bát diện lệch ($Si_6$) và lưỡng tháp ngũ giác ($Si_7$). Ở kích thước $n = 7$, cấu trúc vẫn giữ dạng lưỡng tháp ngũ giác; $n = 8$ có dạng lưỡng tháp ngũ giác gắn thêm 2 nguyên tử Si đối diện ($C_{2v}$); $n = 9, 10$ chuyển dần sang kiểu trụ tam giác; $n = 11$ hình thành bằng cách thế Sc vào $Si_{12}$ (chưa xuất hiện cấu trúc lồng).
  • Dạng anion ($Si_nSc^-$): Khi $n < 6$, cấu trúc hình thành theo quy luật thế ($Si_3Sc^-$ hình thoi, $Si_4Sc^-$ lưỡng tháp tam giác, $Si_5Sc^-$ bát diện, $Si_6Sc^-$ lưỡng tháp ngũ giác với số phối trí của Sc đạt 5). Khi $n \ge 6$, anion chuyển sang phát triển theo quy luật cộng (gắn thêm nguyên tử Si vào cấu trúc anion nhỏ hơn), làm tăng dần số lưỡng tháp ngũ giác lồng vào nhau ($Si_6Sc^-, Si_7Sc^-, Si_8Sc^-$ có 1 lưỡng tháp ngũ giác; $Si_9Sc^-, Si_{10}Sc^-$ có 2 lưỡng tháp ngũ giác lồng nhau; $Si_{11}Sc^-$ có 3 lưỡng tháp ngũ giác lồng nhau), tạo tiền đề cho cấu trúc Frank-Kasper $Si_{16}Sc^-$.

Số liệu tính toán năng lượng liên kết trung bình, chỉ số liên kết Wiberg và năng lượng phân ly liên kết của $Si_nSc^{-/0/+}$ ($n = 2-11$) được tổng hợp trong các bảng dưới đây:

Bảng 1: Năng lượng liên kết trung bình ($E_b$, eV) của cluster $Si_nSc^{-/0/+}$ ($n = 2-11$)

Kích thước ($n$) $Si_nSc^-$ $Si_nSc$ $Si_nSc^+$
2 2,75 2,28 2,22
3 3,17 2,61 2,69
4 3,45 2,89 2,90
5 3,62 3,18 3,17
6 3,60 3,25 3,33
7 3,66 3,30 3,29
8 3,71 3,37 3,43
9 3,76 3,40 3,43
10 3,71 3,42 3,46
11 3,74 3,47 3,48

Bảng 2: Chỉ số liên kết Wiberg tổng của nguyên tử Sc trong cluster $Si_nSc^{-/0/+}$ ($n = 2-11$)

Kích thước ($n$) $Si_nSc^-$ $Si_nSc$ $Si_nSc^+$
2 2,49 2,52 2,76
3 2,80 2,09 2,66
4 3,83 2,87 2,78
5 3,45 3,02 3,04
6 3,71 2,17 2,99
7 3,59 3,06 3,00
8 3,86 2,97 3,39
9 4,02 3,51 3,43
10 4,27 3,17 3,27
11 4,35 3,15 3,34

Bảng 3: Năng lượng phân ly liên kết $D_1$ (tách Sc) và $D_2$ (tách Si) của $Si_nSc^{-/0/+}$ (eV)

$n$ $D_1(Si_nSc)$ $D_1(Si_nSc^-)$ $D_1(Si_nSc^+)$ $D_2(Si_nSc)$ $D_2(Si_nSc^-)$ $D_2(Si_nSc^+)$
2 4,50 1,85 6,17 5,93 4,99 5,65
3 3,81 2,74 5,69 4,86 6,34 5,37
4 3,18 4,98 4,97 5,40 5,83 5,00
5 5,56 4,61 6,61 5,80 5,76 5,81
6 3,70 4,88 6,05 4,90 4,72 5,58
7 5,13 4,88 5,29 4,90 5,39 4,17
8 5,32 5,83 6,95 5,28 5,37 5,89
9 6,64 6,93 6,65 4,92 5,44 4,67
10 5,65 7,23 5,91 4,93 4,52 5,09
11 5,75 7,06 6,20 5,23 5,34 4,99

Đối với các hệ pha tạp 1 nguyên tử kim loại khác:

  • Hệ $Si_nFe^{0/+/-}$ ($n = 8-12$): Đồng phân bền nhất của $Si_8Fe$ có cấu trúc lưỡng tháp ngũ giác kép đối xứng $C_{2v}$ (trạng thái $^{3}B_1$); $Si_9Fe$ có cấu trúc hở đối xứng $C_s$; $Si_{10}Fe$ có cấu trúc lồng; $Si_{11}Fe$ có dạng lưỡng lăng trụ tứ giác; $Si_{12}Fe$ có dạng lăng trụ lục giác đối xứng $D_{3h}$ với Fe ở tâm lồng.
  • Hệ $Si_2M$ ($M = Li, Na, K, Cu, Cr$) và $Si_3M$ ($M = Sc-Zn$): Phân tích bậc liên kết, độ dài liên kết và cấu hình electron hóa trị cho thấy nguyên tử kim loại chuyển dịch mật độ electron sang khung silicon, tạo liên kết mang đặc tính phân cực.

Cluster silicon pha tạp hai nguyên tử kim loại giống nhau

Khảo sát được thực hiện trên các hệ $Si_nTi_2$ ($n = 1-8$), $Si_2M_2$ ($M = Sc-Zn$) và $Si_3M_2$ ($M = Sc-Zn$):

  • Hệ $Si_nTi_2$ ($n = 1-8$): Khảo sát sự phụ thuộc của năng lượng liên kết trung bình $E_b$, năng lượng phân ly $D_1, D_2$, biến thiên năng lượng bậc hai $\Delta^2 E$ và vùng cấm HOMO-LUMO vào kích thước cluster so với dãy $Si_n$ và $Si_{n+1}Ti$.
  • Hệ $Si_2M_2$ và $Si_3M_2$ ($M = Sc-Zn$): Khảo sát 4 dạng đồng phân cấu trúc chính, xác định bậc liên kết $N_{Si-M}$, $N_{M-M}$, độ dài liên kết $d_{Si-M}, d_{M-M}$, điện tích nguyên tử và sự chuyển mật độ electron ($EDT$). Phân tích cấu hình electron tự nhiên (NEC) chỉ ra độ biến thiên mật độ electron trung bình trên phân lớp 3d ($\Delta d$) và 4s ($\Delta s$) của nguyên tử kim loại $M$ khi tham gia liên kết trong khung silicon.

Cluster silicon pha tạp hai nguyên tử kim loại khác nhau

Khảo sát hệ cluster $Si_nScTi^{0/+/-}$ ($n = 1-9$) ở ba trạng thái điện tích:

  • Xác định cấu trúc hình học bền nhất và trạng thái electron của dãy $Si_nScTi^{+/0/-}$.
  • Tính toán năng lượng phân ly liên kết $D_1, D_2, D_3$, biến thiên năng lượng bậc hai $\Delta^2 E$, vùng cấm năng lượng $\Delta E_{HOMO-LUMO}$ và các thông số ion hóa ($AIE$), ái lực electron ($EA$).
  • Phân tích điện tích tự nhiên của Sc và Ti, độ biến thiên mật độ electron phân lớp d và s ($\Delta d_c, \Delta s_c; \Delta d_a, \Delta s_a; \Delta d_n, \Delta s_n$) và xây dựng giản đồ orbital phân tử (MO) cho cấu trúc $Si_4ScTi$.

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

Các kết quả và đóng góp của luận án được thể hiện qua các điểm chính:

  1. Về mặt lý thuyết và quy luật cấu trúc:
    • Xác định được quy luật phát triển hình học của cluster silicon pha tạp kim loại chuyển tiếp: dạng trung hòa và cation $Si_nSc^{0/+}$ tuân theo quy luật thế nguyên tử trong khung silicon kích thước $n+1$, trong khi dạng anion $Si_nSc^-$ ($n \ge 6$) phát triển theo quy luật cộng nguyên tử Si, hình thành các cấu trúc lưỡng tháp ngũ giác lồng nhau đặc khít.
    • Làm rõ ảnh hưởng của điện tích: Việc nhận thêm 1 electron (tạo anion) làm tăng độ bền nhiệt động rõ rệt (năng lượng liên kết trung bình $E_b$ tăng khoảng 0,5 eV, chỉ số liên kết Wiberg của Sc tăng từ ~2,49 lên 4,35), trong khi việc mất 1 electron (tạo cation) ít làm thay đổi cấu trúc và độ bền so với dạng trung hòa.
    • Giải thích bản chất liên kết hóa học qua phân tích NBO/NRT và MO: Định lượng sự chuyển dịch mật độ electron ($EDT$) từ kim loại sang silicon, xác định vai trò tham gia tạo liên kết của các phân lớp 3d và 4s thông qua biến thiên mật độ electron $\Delta d$ và $\Delta s$.
  2. Về mặt thực tiễn:
    • Cung cấp cơ sở dữ liệu về cấu trúc hình học, trạng thái spin, năng lượng ion hóa đoạn nhiệt ($AIE$), ái lực electron ($EA$), năng lượng vùng cấm và tần số dao động điều hòa của các hệ cluster silicon pha tạp kim loại dãy 3d.
    • Cung cấp các thông số tham chiếu phục vụ cho việc định hướng và đối chứng với các nghiên cứu thực nghiệm đo phổ khối lượng, phổ photoelectron và phổ IR-MPD.

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

  • Hạn chế:
    • Các phương pháp tính toán hóa học lượng tử dựa trên mô hình gần đúng DFT luôn tồn tại sai số nhất định đối với các hệ chứa kim loại chuyển tiếp có phân lớp d phức tạp.
    • Đặc tính linh động của cluster với sự tồn tại đồng thời của nhiều đồng phân có mức năng lượng gần nhau dẫn đến khó khăn trong việc khẳng định cấu trúc duy nhất trong điều kiện thực tế nếu không có thực nghiệm đối chứng trực tiếp.
  • Hướng nghiên cứu tiếp:
    • Mở rộng nghiên cứu cấu trúc electron trên các hệ cluster pha tạp khác để khái quát hóa quy tắc đếm electron cho cluster silicon.
    • Khảo sát giới hạn cuối về độ bền nhiệt động có thể đạt được của các cluster silicon kích thước nhỏ.
    • Tìm kiếm và phân loại các nhóm nguyên tử pha tạp tối ưu nhằm nâng cao độ bền và định hướng tính chất của vật liệu nano silicon.

Giá trị tham khảo

Luận án là tài liệu tham khảo cho:

  • Nghiên cứu sinh, học viên cao học, giảng viên và các nhà nghiên cứu chuyên ngành Hóa vô cơ, Hóa học lý thuyết và Hóa học tính toán, Vật lý chất rắn và Khoa học vật liệu nano.
  • Nội dung tham khảo trọng tâm:
    • Quy trình và mức lý thuyết tính toán DFT B3P86/6-311+G(d) cho các hệ kim loại chuyển tiếp liên kết với bán dẫn.
    • Phương pháp phân tích liên kết hóa học NBO/NRT và giản đồ MO trên các hệ cụm nguyên tử.
    • Tập hợp số liệu về cấu trúc, độ dài liên kết, bậc liên kết Wiberg và các đại lượng năng lượng ($E_b, D, \Delta^2 E, E_{gap}, AIE, EA$) của các hệ cluster $Si_nSc$, $Si_nFe$, $Si_3M$, $Si_nTi_2$, $Si_3M_2$, $Si_nScTi$.

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

1. Tại sao mức lý thuyết B3P86/6-311+G(d) được lựa chọn trong luận án?
Mức lý thuyết B3P86/6-311+G(d) kết hợp phiếm hàm trao đổi Becke B3, phiếm hàm tương quan Perdew 1986 và 20% năng lượng trao đổi Hartree-Fock cùng bộ hàm cơ sở hóa trị tách ba có bổ sung hàm phân cực d và hàm khuếch tán. Khi kiểm chứng trên đime $Si_2$, phương pháp cho độ dài liên kết $r = 2,27\text{ \AA}$, tần số dao động $\omega_e = 495,53\text{ cm}^{-1}$ và năng lượng liên kết $E_b = 3,17\text{ eV}$, phù hợp với số liệu thực nghiệm ($r = 2,25\text{ \AA}$, $\omega_e = 509 \pm 10\text{ cm}^{-1}$, $E_b = 3,21\text{ eV}$).

2. Quy luật hình thành cấu trúc của $Si_nSc^{-/0/+}$ ($n = 2-11$) khác nhau như thế nào giữa các trạng thái điện tích?
Dạng trung hòa và cation ($Si_nSc^{0/+}$) phát triển chủ yếu theo quy luật thế nguyên tử Si trong khung $Si_{n+1}$ bằng nguyên tử Sc. Ngược lại, dạng anion ($Si_nSc^-$) ở kích thước $n \ge 6$ phát triển theo quy luật cộng nguyên tử Si vào cấu trúc nhỏ hơn, làm tăng dần số lưỡng tháp ngũ giác lồng nhau và tăng số phối trí của nguyên tử Sc.

3. Sự thay đổi điện tích ảnh hưởng như thế nào đến độ bền của cluster silicon pha tạp Scandium?
Việc thêm 1 electron (tạo anion $Si_nSc^-$) làm tăng độ bền đáng kể, nâng năng lượng liên kết trung bình $E_b$ cao hơn dạng trung hòa và cation khoảng 0,5 eV và tăng chỉ số liên kết Wiberg của Sc (đạt 4,35 ở $n = 11$). Việc bớt 1 electron (tạo cation $Si_nSc^+$) làm thay đổi độ bền không đáng kể so với dạng trung hòa.

4. Vì sao cluster silicon thuần khiết kém bền hơn cluster carbon và cần phải pha tạp kim loại?
Khác với nguyên tử carbon lai hóa $sp^2$ tạo nên cấu trúc lồng khép kín bền vững ($C_{60}$), nguyên tử silicon có xu hướng tạo liên kết lai hóa $sp^3$ tứ diện, khiến cấu trúc lồng thuần khiết kém bền về mặt nhiệt động. Việc pha tạp kim loại giúp ổn định khung cấu trúc, tăng năng lượng liên kết trung bình và điều chỉnh các đặc tính electron của cluster.

5. Những phần mềm nào được sử dụng để khảo sát cấu trúc và liên kết trong luận án?
Luận án sử dụng phần mềm Gaussian 09 (bản E.01) để tính toán hóa học lượng tử; Gaussview 05 để hiển thị cấu trúc; NBO 5.G để phân tích liên kết tự nhiên; PyMolyze, gOpenMol, Dgrid-4 để khảo sát cấu trúc electron và orbital; Molden, JMol, Corel Drawx3 để hiển thị hình học; Excel và Origin để xử lý dữ liệu.

Kết luận

Luận án của tác giả Phạm Ngọc Thạch đã sử dụng phương pháp hóa học tính toán phiếm hàm mật độ ở mức lý thuyết B3P86/6-311+G(d) kết hợp phân tích NBO/NRT để khảo sát hệ thống cấu trúc, độ bền và tính chất của các dãy cluster silicon pha tạp một và hai nguyên tử kim loại. Kết quả nghiên cứu đã xác định hình học của các đồng phân bền nhất, quy luật phát triển cấu trúc theo trạng thái điện tích, ảnh hưởng của sự chuyển dịch mật độ electron và bản chất liên kết giữa kim loại với khung silicon. Công trình cung cấp tập dữ liệu lý thuyết phục vụ cho việc nghiên cứu vật liệu nano bán dẫn trên cơ sở cluster silicon.