Tổng quan luận án

Luận án tiến sĩ Hóa học chuyên ngành Hóa lý thuyết và Hóa lý (Mã số: 9440119) của tác giả Dương Văn Long, được thực hiện tại Trường Đại học Quy Nhơn dưới sự hướng dẫn của GS. Nguyễn Minh Thơ và PGS. TS. Nguyễn Phi Hùng, với tiêu đề: "A Quantum Chemical Research of Structure and Aromaticity of Some Boron Clusters" (Nghiên cứu hóa học lượng tử về cấu trúc và tính thơm của một số cụm boron).

                      HỘI ĐỒNG VÀ THÔNG TIN ĐÀO TẠO

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

Khoa học về cụm nguyên tử (cluster science) tập trung vào việc tìm kiếm các cấu trúc bền vững của các nhóm nguyên tử, giải thích độ bền nhiệt động học, làm sáng tỏ tính chất hóa lý và khám phá ứng dụng thực tiễn. Trong số đó, các cụm boron nhận được sự quan tâm lớn do boron là nguyên tố thiếu hụt electron (electron-deficient), tạo ra nhiều dạng cấu trúc hình học đa dạng từ phẳng (planar), gần phẳng/bán phẳng (quasi-planar), lồng icosahedral, lồng fullerene đến dạng ống (tubular).

Độ bền của cụm boron gắn liền chặt chẽ với khái niệm tính thơm (aromaticity). Tuy nhiên, các quy tắc kinh điển như quy tắc đếm electron Hückel $(4n + 2)$ ban đầu chỉ được xây dựng cho các hydrocarbon thơm dạng vòng phẳng. Khi áp dụng quy tắc $(4n + 2)$ cho các cụm phi hydrocarbon, các cấu trúc 2D phức tạp hoặc cấu trúc 3D, việc lạm dụng quy tắc mà bỏ qua giải phương trình thế đã dẫn đến những diễn giải sai lệch bản chất. Do đó, yêu cầu cấp thiết đặt ra là phải kiểm chứng, cải tiến và xây dựng các mô hình tính thơm mới dựa trên giải phương trình sóng lượng tử chuẩn xác cho từng dạng hình học đặc thù của cụm boron.

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

  1. Khảo sát một loạt các cụm boron nguyên chất và cụm boron pha tạp nhằm xác định các cấu trúc bền nhiệt động, cấu hình electron và xác định tiềm năng ứng dụng của chúng.
  2. Khảo sát và làm rõ tính đa dạng về hình học của các cụm boron; phân tích việc sử dụng các mô hình tính thơm tương ứng cho từng dạng cấu trúc hình học cụ thể.
  3. Kiểm chứng tính phù hợp của các mô hình tính thơm hiện đại (mô hình đĩa tròn - disk model, mô hình dải ruy-băng - ribbon model, mô hình hình trụ rỗng - hollow cylinder model), tổng quát hóa quy tắc đếm electron $(4n + 2m)$ cho cấu trúc dạng đĩa tròn và đề xuất mô hình tính thơm đĩa - nón (disk-cone model).

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

  • Đối tượng nghiên cứu: Các hệ cụm boron và boron pha tạp cụ thể:
    • Cụm pha tạp silic: $\text{B}_2\text{Si}_3^q$ ($q = -2, -1, 0, +1, +2$) và $\text{B}_3\text{Si}_2^p$ ($p = -3, -2, -1, 0, +1$).
    • Cụm boron kích thước lớn: $\text{B}{70}$ ở trạng thái trung hòa và dianion ($\text{B}{70}^{0/2-}$).
    • Hệ cụm nhị phân boron - liti: $\text{B}_{12}\text{Li}n$ với $n = 0 - 14$ (đặc biệt là $\text{B}{12}\text{Li}4$ và $\text{B}{12}\text{Li}_8$).
    • Cụm pha tạp kim loại chuyển tiếp: $\text{B}_{14}\text{FeLi}_2$ và các cấu trúc dây nano dẫn xuất.
  • Phạm vi phương pháp: Phương pháp hóa học lượng tử ab initio và lý thuyết hàm mật độ (DFT), kết hợp thuật toán tìm kiếm cấu trúc ngẫu nhiên, giải phương trình sóng giải tích và các phương pháp phân tích liên kết hóa học.

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

Khái niệm cụm nguyên tử kim loại được F. Cotton đề xuất lần đầu vào năm 1964. Khái niệm tính thơm khởi nguồn từ phát hiện benzen ($\text{C}_6\text{H}_6$) của Faraday năm 1825, sau đó mô hình Hückel (1931) ra đời mô tả hệ liên hợp $\pi$ của annulene vòng phẳng với quy tắc $(4n + 2)$ electron $\pi$. Năm 1972, Baird đề xuất quy tắc cho trạng thái triplet, chỉ ra hệ vòng liên hợp có $4n$ electron $\pi$ mang tính thơm triplet. Trong hóa học borane, quy tắc Wade-Mingos sử dụng số cặp electron khung $(2n + 2)$ để dự đoán tính thơm của các cấu trúc lồng đa diện. Quy tắc Möbius cũng được mở rộng cho một số phức chất kim loại chuyển tiếp.

Tác giả / Mô hình Đối tượng / Bản chất mô hình Phạm vi và quy tắc đếm electron
Hückel (1931) Hydrocarbon vòng phẳng liên hợp $(4n + 2)$ electron $\pi$ cho trạng thái singlet đóng vỏ; $4n$ là phản thơm
Baird (1972) Vòng phẳng liên hợp trạng thái triplet $4n$ electron $\pi$ là thơm triplet; $(4n + 2)$ là phản thơm triplet
Wade-Mingos (1971, 1984) Cụm borane đa diện (polyhedral boranes) $(2n + 2)$ electron khung delocalized
Tai và cộng sự (2012) Mô hình đĩa tròn (Disk Model) Hạt chuyển động tự do trên mặt phẳng tròn có thế tường vô hạn
Nhóm nghiên cứu (2017) Mô hình dải ruy-băng (Ribbon Model) Hạt chuyển động trong hộp thế 1D, cấu hình chuẩn $\pi^{2(n+1)}\sigma^{2n}$
Mô hình trụ rỗng (HCM) Cụm dạng ống nhiều lớp (TR, nanotubes) Quy tắc $(4n + 2m)$ với $m$ mức năng lượng không suy biến

Các nghiên cứu trước đây chỉ ra rằng cụm boron nhỏ như $\text{B}3^-$ và $\text{B}{12}$ phẳng thỏa mãn quy tắc Hückel $(4n + 2)$ cho cả electron $\pi$ và $\sigma$; $\text{B}7^-$ trạng thái triplet tuân theo quy tắc Baird. Tuy nhiên, khi kích thước cụm tăng lên, các cấu trúc như $\text{B}{19}^-$ hay $\text{B}{20}^{2-}$ chứa 12 electron $\pi$ vi phạm quy tắc Hückel nhưng lại phù hợp với mô hình đĩa tròn (Disk Model). Cụm ống hai vòng (double-ring - DR) nhỏ nhất $\text{B}{20}$ và cation $\text{B}{14}^{2+}$ tuân theo $(4n + 2)$, nhưng các cấu trúc ba vòng (triple-ring - TR) như $\text{B}{27}^+$ và $\text{B}{42}$ đòi hỏi quy tắc $(4n + 2m)$ thông qua mô hình hình trụ rỗng (HCM). Cụm dải chuỗi kép $\text{B}{22}\text{H}_2^{2-}$ có thể kéo dài tới 17.0 Å nhờ cấu hình electron đặc trưng $\pi^{2(n+1)}\sigma^{2n}$.

Khoảng trống nghiên cứu luận án lựa chọn giải quyết:

  1. Phân định ranh giới và mối quan hệ giữa mô hình Hückel kinh điển và mô hình dải ruy-băng trong việc giải thích cấu trúc pha tạp $\text{B}_2\text{Si}_3^q$ và $\text{B}_3\text{Si}_2^p$.
  2. Kiểm chứng và giải thích trạng thái cơ bản triplet cùng tính thơm của cụm boron bán phẳng lớn $\text{B}_{70}^{0/2-}$ bằng mô hình đĩa tròn và đề xuất quy tắc đếm electron $(4n + 2m)$ cho hệ đĩa.
  3. Làm sáng tỏ quy luật phát triển cấu trúc khi pha tạp liti vào cụm $\text{B}{12}$ ($\text{B}{12}\text{Li}n$, $n = 0 - 14$), đề xuất mô hình đĩa - nón (disk-cone model) cho cấu trúc nón $\text{B}{12}\text{Li}4$, và đánh giá khả năng lưu trữ hydro của $\text{B}{12}\text{Li}_8$.
  4. Làm rõ cơ chế dập tắt từ tính của Fe và độ bền của cấu trúc dạng con quay $\text{B}_{14}\text{FeLi}_2$ bằng mô hình hình trụ rỗng (HCM) và thiết kế vật liệu dây nano 1D.

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

Cơ sở lý thuyết và các mô hình lượng tử

  • Phương trình Schrödinger và Gần đúng Born–Oppenheimer: Nền tảng phân tách chuyển động của hạt nhân và electron.
  • Phương pháp trường tự phối Hartree–Fock (HF) và Hậu Hartree–Fock: Bao gồm lý thuyết nhiễu loạn Møller–Plesset ($MPn$), Tương tác cấu hình (CI, MRCI), Trường tự phối đa cấu hình không gian hoạt động hoàn toàn (CASSCF), Cụm liên kết (Coupled Cluster: CC, CCSD, CCSD(T)) và ngoại suy về tập hàm cơ sở đầy đủ (CBS).
  • Lý thuyết hàm mật độ (DFT): Sử dụng các phiếm hàm LDA, GGA, meta-GGA và hybrid meta-GGA (TPSSh, B3LYP, wB97XD).
  • Mô hình dải ruy-băng (Ribbon Model): Mô tả electron $\pi$ và $\sigma$ delocalized như các hạt chuyển động tự do trong hộp thế một chiều (1D box) có chiều dài tương ứng $l_\pi$ và $l_\sigma$: $$E_n = \frac{n^2 h^2}{8m l^2}$$ Hàm sóng chuẩn hóa có dạng: $\psi_n(x) = \sqrt{\frac{2}{l}} \sin\left(\frac{n\pi x}{l}\right)$. Tối thiểu hóa phi tuyến tổng bình phương phần dư (RSS) để xác định các hằng số $C_1, C_2$ và độ dài hộp thế $l_\pi, l_\sigma$.
  • Mô hình đĩa tròn (Disk Model): Giải phương trình Schrödinger 2D trong hệ tọa độ cực $(r, \theta)$ với thế thành vô hạn tại $r = R$: Hàm góc: $\Theta_m(\theta) = \frac{1}{\sqrt{2\pi}} e^{im\theta}$ với $m = 0, \pm 1, \pm 2, \dots$ Hàm bán kính dẫn về phương trình Bessel bậc $m$: $\mathcal{R}m(r) = J_m(\alpha r)$, với điều kiện biên $J_m(\alpha R) = 0$. Mức năng lượng xác định thông qua các nghiệm không âm $\alpha{m,n} R$ của hàm Bessel $J_m(x)$.
  • Mô hình hình trụ rỗng (Hollow Cylinder Model - HCM): Mô tả hạt chuyển động trong không gian giới hạn giữa bán kính trong $R_0 = R - r$ và bán kính ngoài $R_1 = R + r$, chiều cao $L$.

Công cụ và kỹ thuật phân tích hóa học lượng tử

  • Thuật toán tìm kiếm cấu trúc đáy (Stochastic search): Kết hợp kỹ thuật tạo dịch chuyển ngẫu nhiên (random kick) và thuật toán di truyền (Genetic Algorithm - GA / MEGA) để quét bề mặt thế năng (PES).
  • Tối ưu hóa hình học và tần số dao động: Sử dụng phần mềm Gaussian 16 và gói mô phỏng VASP. Các cấu trúc trong khoảng chênh lệch năng lượng $\sim 2\text{ eV}$ được tái tối ưu hóa ở mức lý thuyết TPSSh/6-311+G(d) và tính toán tần số dao động điều hòa để khẳng định điểm cực tiểu thực (không có tần số ảo, $N_{\text{imag}} = 0$).
  • Hàm định xứ electron (Electron Localization Function - ELF): Sử dụng phần mềm Dgrid-5 để phân tách ELF tổng thành các thành phần định xứ $\sigma$ ($\text{ELF}{\sigma l}$), giải tỏa $\sigma$ ($\text{ELF}{\sigma d}$) và giải tỏa $\pi$ ($\text{ELF}_\pi$).
  • Bản đồ dòng vòng cảm ứng từ (Ring current maps):
    • Phương pháp CTOCD-DZ2 thực thi trong chương trình SYSMOIC, tính toán tenxơ mật độ dòng cảm ứng.
    • Phương pháp dị hướng của mật độ dòng cảm ứng (ACID) để xác định tính thơm (dòng diatropic cùng chiều kim đồng hồ) và phản thơm (dòng paratropic ngược chiều kim đồng hồ).
  • Điện tích nguyên tử thực (NAC) và Bậc liên kết (Bond Order): Sử dụng phương pháp phân chia mật độ electron DDEC6.
  • Dữ liệu thực nghiệm đối sánh: Dữ liệu phổ quang electron (Photoelectron Spectroscopy - PES) và phổ quang ion hóa hai màu hồng ngoại - tử ngoại cộng hưởng (IR-UV2CI).

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

Chương 1 & 2: Cơ sở lý thuyết, Phương pháp tính toán và Khảo sát chuẩn hóa (Benchmarking)

  • Luận án thiết lập quy trình tính toán chuẩn mực cho các cụm chứa boron và silic. Kết quả tính toán chuẩn hóa (benchmark) so sánh độ lệch năng lượng đơn điểm $\text{(U)CCSD(T)/CBS}$ cho thấy phiếm hàm TPSSh là lựa chọn phù hợp nhất cho việc tối ưu hóa hình học cấu trúc cụm hỗn hợp $\text{B-Si}$.
  • Đối với việc mô phỏng các giá trị năng lượng tách electron thẳng đứng (VDE), năng lượng tách electron đoạn nhiệt (ADE) trong phổ PES của $\text{B}_2\text{Si}_3^-$ và phổ dao động điều hòa so với phổ IR-UV2CI thực nghiệm của $\text{B}_2\text{Si}_3$, phiếm hàm B3LYP/6-311+G(d) cho kết quả tiệm cận tốt nhất với số liệu thực nghiệm.

Chương 3: Quy tắc Hückel và mô hình dải ruy-băng: Trường hợp các cụm $\text{B}_2\text{Si}_3^q$ và $\text{B}_3\text{Si}_2^p$

  • Khảo sát cấu trúc: Tối ưu hóa cấu trúc các đồng phân năng lượng thấp của $\text{B}_2\text{Si}_3^q$ ($q = -2$ đến $+2$) và $\text{B}_3\text{Si}_2^p$ ($p = -3$ đến $+1$) ở mức TPSSh/6-311+G(d).
  • Phát hiện cơ chế tiến hóa cấu trúc: Sự thay thế đẳng điện tử một đơn vị $\text{B}^-$ bằng nguyên tử $\text{Si}$ trong khung dải ruy-băng trianion $\text{B}_3\text{Si}_2^{3-}$ tại hai vị trí khác nhau dẫn đến hai nhánh đồng phân: dạng dải ruy-băng (Ribbon - R) và dạng vòng Hückel (Hückel - H).
  • Hiện tượng tự khóa (Self-locking phenomenon): Phân tích ELF cho thấy sự co ngắn mạnh của các liên kết biên tận cùng (terminal bonds) tạo nên hàng rào lực đẩy Coulomb liên nguyên tử cao, đóng vai trò như hai thành vô hạn của hộp thế 1D giam giữ các electron $\pi$ và $\sigma$ delocalized.
  • Quy luật cấu hình electron:
    • Cấu hình $\pi^{2(n+1)}\sigma^{2n}$ (như $\text{B}{14}\text{H}2^{2-}$) mang lại độ bền vượt trội nhờ sự xen kẽ liên tục và hỗ trợ lẫn nhau của các bồn $\text{ELF}\pi$ và $\text{ELF}{\sigma d}$.
    • Cấu hình $\pi^{2n}\sigma^{2n}$ (như singlet $\text{B}_{12}\text{H}2^{2-}$) dẫn đến sự gián đoạn các bồn ELF, tạo nên các mắt xích yếu hình thoi ($\text{B}2\text{–B}3\text{–B}{10}\text{–B}{11}$) và mang tính phản thơm dạng ruy-băng ($N{\text{imag}} = 1$ tại $-178\text{ cm}^{-1}$).
    • Trạng thái triplet $\text{B}_{12}\text{H}_2^{2-}$ với cấu hình $\pi^{2n+1}\sigma^{2n-1}$ loại bỏ được mắt xích yếu và thể hiện tính thơm triplet dạng ruy-băng.
                   CẤU HÌNH ELECTRON VÀ TÍNH THƠM DẢI RUY-BĂNG

Chương 4: Tính thơm đĩa trên cụm boron bán phẳng $\text{B}_{70}^{0/2-}$

  • Cấu trúc và trạng thái cơ bản: Tạo cấu trúc bán phẳng gồm 70 nguyên tử boron bằng thuật toán bước cóc tô pô (topological leapfrog algorithm) từ cấu trúc ban đầu $\text{B}_{16}$ (13 đỉnh).
  • Phát hiện quan trọng: Cụm $\text{B}{70}$ trung hòa dạng bán phẳng có trạng thái cơ bản là trạng thái triplet. Khi nhận thêm 2 electron để tạo dianion $\text{B}{70}^{2-}$, cấu trúc đạt độ bền nhiệt động học rất cao.
  • Áp dụng mô hình đĩa tròn và quy tắc đếm electron: Các orbital phân tử $\pi$ (CMO) của $\text{B}_{70}^{2-}$ tương ứng chính xác với các mức năng lượng giải tích từ các nghiệm của hàm Bessel $J_m(x)$.
  • Tổng quát hóa quy tắc đếm electron: Luận án đề xuất quy tắc đếm $(4n + 2m)$ cho các cấu trúc dạng đĩa tròn, trong đó $m$ là số mức năng lượng không suy biến (ứng với $m=0$). Bản đồ mật độ dòng từ tổng của $\text{B}{70}^{2-}$ và triplet $\text{B}{70}$ khẳng định dòng cảm ứng diatropic mạnh mẽ đặc trưng cho tính thơm.

Chương 5: Cụm nhị phân boron - liti $\text{B}_{12}\text{Li}n$ ($n = 1 - 14$) và mô hình đĩa - nón cho $\text{B}{12}\text{Li}_4$

  • Quy luật phát triển cấu trúc: Khảo sát toàn diện chuỗi $\text{B}_{12}\text{Li}_n$ ($n = 0 - 14$) ở mức $\text{(U)CCSD(T)/cc-PVTZ + ZPE // TPSSh/6-311+G(d)}$. Phân tích các thông số: năng lượng liên kết trên mỗi nguyên tử, giá trị trung bình tổng bậc liên kết (ASBO), khe năng lượng HOMO-LUMO (HLG), năng lượng phân ly và sai phân năng lượng bậc hai ($\Delta^2 E$).
  • Ứng dụng tích trữ Hydro của $\text{B}_{12}\text{Li}_8$: Kết quả chỉ ra đồng phân $8A$ của $\text{B}_{12}\text{Li}_8$ có độ bền vượt trội. Khảo sát hấp phụ phân tử hydro từ $8A\text{-}8\text{H}_2$ đến $8A\text{-}40\text{H}2$ ở mức lý thuyết wB97XD/6-311++G(2d,2p) cho thấy $\text{B}{12}\text{Li}_8$ là ứng viên tiềm năng cho vật liệu lưu trữ hydro và pin ion liti.
  • Mô hình đĩa - nón (Disk-Cone Model) cho $\text{B}_{12}\text{Li}_4$: $\text{B}_{12}\text{Li}4$ có cấu trúc dạng nón bền vững, đẳng điện tử với $\text{B}{13}\text{Li}$. Bản đồ dòng dòng mật độ trên mặt nón chỉ ra dòng chảy diatropic mạnh. Luận án đề xuất mô hình đĩa - nón giải thích sự chuyển dịch mức orbital hóa trị khi cấu trúc phẳng uốn cong thành dạng nón.

Chương 6: $\text{B}_{14}\text{FeLi}_2$ và mô hình hình trụ rỗng (HCM)

  • Cấu trúc dạng con quay (Teetotum): $\text{B}_{14}\text{FeLi}2$ có cấu trúc bền nhiệt động với trục $\text{Li-Fe-Li}$ nằm dọc theo tâm của khung ống boron $\text{B}{14}$ hai vòng.
  • Giải thích bằng mô hình HCM: Sự kết hợp giữa các orbital phân tử của khung $\text{B}_{14}$ (được phân loại theo mô hình HCM gồm các tập orbital định xứ, tiếp tuyến - tangential, và xuyên tâm - radial) với các orbital nguyên tử $d$ của Fe dẫn đến việc ghép đôi hoàn toàn các electron, dập tắt hoàn toàn từ tính của nguyên tử sắt (quenched magnetic moment).
  • Phát triển vật liệu nano: Tính toán quang phổ hấp thụ electron (UV-Vis) và tối ưu hóa cấu trúc các dây nano 1D được thiết kế từ đơn vị cơ sở $\text{B}{14}\text{FeLi}2$, bao gồm dây nano $[\text{Li-B}{14}\text{Fe-Li}]\text{-}[\text{B}{14}\text{Fe}]\text{-}[\text{Li-B}{14}\text{Fe-Li}]$ và chuỗi dây nano lăng trụ/phản lăng trụ $\text{B}{28}\text{Fe}_2\text{Li}2\text{Mg}$, $\text{B}{42}\text{Fe}_3\text{Li}_2\text{Mg}2$, $\text{B}{56}\text{Fe}_4\text{Li}_2\text{Mg}3$, $\text{B}{70}\text{Fe}_4\text{Mg}_4\text{Li}_2$ cho các thiết bị quang điện (photovoltaic devices).

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

  1. Đóng góp về phương pháp luận và mô hình lý thuyết:
    • Xác lập rõ ràng ranh giới phân biệt giữa mô hình Hückel kinh điển và các mô hình hiện đại (mô hình dải ruy-băng, mô hình đĩa tròn, mô hình hình trụ rỗng), ngăn chặn việc lạm dụng quy tắc $(4n + 2)$ cho các cấu trúc không phù hợp.
    • Làm rõ hiện tượng tự khóa (self-locking phenomenon) ở các liên kết biên của dải ruy-băng boron thông qua phân tích hàm $\text{ELF}{\sigma l}, \text{ELF}{\sigma d}$ và $\text{ELF}_\pi$.
    • Tổng quát hóa quy tắc đếm electron $(4n + 2m)$ cho các cấu trúc dạng đĩa tròn (Disk Model) dựa trên nghiệm hàm Bessel và xây dựng thành công mô hình lai đĩa - nón (Disk-Cone Model).
    • Đưa ra cơ chế tương tác MO trong mô hình HCM giải thích hiện tượng dập tắt từ tính của nguyên tử kim loại Fe trong lồng boron hình ống $\text{B}_{14}\text{FeLi}_2$.
  2. Đóng góp về phát hiện cấu trúc và tính chất vật liệu:
    • Phát hiện trạng thái cơ bản triplet của cụm boron bán phẳng lớn $\text{B}{70}$ và độ bền nhiệt động học cao của dianion $\text{B}{70}^{2-}$.
    • Xác định $\text{B}_{12}\text{Li}8$ là cấu trúc bền vững nhất trong chuỗi $\text{B}{12}\text{Li}_n$, có khả năng hấp phụ tối đa tới 40 phân tử $\text{H}_2$, định hướng ứng dụng làm vật liệu lưu trữ hydro.
    • Thiết kế các mô hình dây nano đơn chiều (1D nanowires) từ cấu trúc $\text{B}_{14}\text{FeLi}_2$ kết hợp kim loại $\text{Mg}$ và $\text{Li}$, mở ra hướng ứng dụng cho vật liệu quang điện tử.

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

  • Hạn chế:
    • Nghiên cứu tập trung chủ yếu vào các tính toán mô phỏng lý thuyết lượng tử pha khí ở nhiệt độ $0\text{ K}$, chưa khảo sát đầy đủ ảnh hưởng của hiệu ứng dung môi, chất nền hoặc nhiệt độ hữu hạn lên động học chuyển hóa cấu trúc.
    • Các cấu trúc phức tạp như $\text{B}_{70}$ và chuỗi dây nano mới chỉ dừng lại ở tính toán ab initio/DFT, đòi hỏi các điều kiện thực nghiệm chuyên biệt (như phổ khối, bẫy ion, kỹ thuật chùm laser) để tổng hợp và kiểm chứng trực tiếp.
  • Hướng nghiên cứu tiếp:
    • Mở rộng mô hình đĩa - nón và mô hình hình trụ rỗng (HCM) cho các hệ cụm boron pha tạp kim loại chuyển tiếp khác ($d$ và $f$-block elements).
    • Nghiên cứu cơ chế động học hấp phụ và giải hấp hydro trên vật liệu $\text{B}_{12}\text{Li}_8$ ở điều kiện nhiệt độ và áp suất thực tế.
    • Khảo sát sâu hơn các đặc tính dẫn điện và hiệu suất quang điện của các chuỗi dây nano đa kim loại phát triển từ $\text{B}_{14}\text{FeLi}_2$.

Giá trị tham khảo

  • Nghiên cứu sinh và Giảng viên chuyên ngành Hóa lý thuyết, Hóa lý, Hóa học tính toán:
    • Tài liệu tham khảo về quy trình giải tích cơ học lượng tử áp dụng cho các hạt chuyển động trong các hộp thế có hình học đặc thù (hộp 1D, đĩa tròn 2D, ống trụ rỗng 3D).
    • Phương pháp luận phân tích liên kết hiện đại kết hợp giữa hàm định xứ electron phân tách ($\text{ELF}\sigma, \text{ELF}\pi$), tenxơ mật độ dòng cảm ứng (SYSMOIC, ACID) và phân tích điện tích DDEC6.
  • Nhà khoa học trong lĩnh vực Khoa học Vật liệu và Năng lượng mới:
    • Dữ liệu cấu trúc, năng lượng liên kết và dung lượng hấp phụ hydro của cụm liti-boron phục vụ phát triển vật liệu tích trữ năng lượng.
    • Thiết kế cấu trúc các khối cơ bản (building blocks) dạng dây nano kim loại - boron cho linh kiện quang điện tử và cảm biến.

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

1. Tại sao quy tắc Hückel $(4n + 2)$ kinh điển không còn phù hợp để giải thích độ bền cho các cụm boron lớn và cụm dạng dải?

Quy tắc Hückel $(4n + 2)$ được thiết lập dựa trên việc giải phương trình thế cho các phân tử hydrocarbon vòng phẳng đơn lớp. Đối với cụm boron, do tính chất thiếu hụt electron và sự phân bố đa hướng của cả electron $\sigma$ và $\pi$ delocalized, hình học của cụm biến đổi sang dạng dải dài (ribbon), dạng đĩa (disk) hoặc ống trụ (hollow cylinder). Việc áp dụng quy tắc $(4n + 2)$ bỏ qua bản chất hàm sóng hình học dẫn đến dự đoán sai lệch (ví dụ cụm dải có $4n$ electron vẫn bền, hoặc cụm đĩa $\text{B}{19}^-$, $\text{B}{20}^{2-}$ chứa 12 electron $\pi$ nhưng vẫn mang tính thơm).

2. Hiện tượng tự khóa (self-locking phenomenon) trong mô hình dải ruy-băng là gì?

Hiện tượng tự khóa là sự co ngắn bất thường của các liên kết boron–boron ở hai đầu tận cùng của khung dải ruy-băng (như liên kết $\text{B}1\text{–B}{14}$ và $\text{B}_7\text{–B}8$ trong $\text{B}{14}\text{H}_2^{2-}$). Sự co ngắn này làm tăng lực đẩy Coulomb liên hạt nhân, tạo ra một hàng rào thế năng cao ở hai biên, giam giữ các electron $\pi$ và $\sigma$ giải tỏa chuyển động tự do bên trong như một hạt trong hộp thế một chiều (1D box).

3. Cấu hình electron nào đảm bảo độ bền tối ưu cho cấu trúc dải ruy-băng boron?

Cấu hình electron tối ưu là $\pi^{2(n+1)}\sigma^{2n}$. Cấu hình này tạo ra sự phân bố xen kẽ liên tục và hỗ trợ đồng thời giữa các bồn $\text{ELF}\pi$ và $\text{ELF}{\sigma d}$, triệt tiêu các điểm yếu liên kết trong cấu trúc. Ngược lại, cấu hình $\pi^{2n}\sigma^{2n}$ làm ngắt quãng các bồn ELF, tạo ra các mắt xích hình thoi yếu và khiến cấu trúc rơi vào trạng thái phản thơm.

4. Quy tắc đếm electron trong mô hình đĩa tròn (Disk Model) được biểu diễn như thế nào?

Dựa trên việc giải phương trình Schrödinger trong hệ tọa độ cực với điều kiện biên hàm Bessel $J_m(\alpha R) = 0$, mức năng lượng của hệ đĩa phụ thuộc vào các nghiệm $\alpha_{m,n} R$. Do các trạng thái với $m \neq 0$ có tính suy biến bậc hai, còn trạng thái $m = 0$ không suy biến, quy tắc đếm electron tổng quát cho hệ đĩa tròn đóng vỏ bền vững có dạng $(4n + 2m)$, trong đó $m$ là số mức năng lượng không suy biến được điền đầy.

5. Nguyên nhân nào khiến từ tính của nguyên tử sắt (Fe) bị dập tắt trong cụm $\text{B}_{14}\text{FeLi}_2$?

Trong cấu trúc con quay $\text{B}{14}\text{FeLi}2$, nguyên tử Fe nằm ở tâm của ống boron $\text{B}{14}$. Các orbital $d$ của Fe tương tác xen phủ mạnh mẽ với các orbital phân tử hóa trị của khung trụ $\text{B}{14}$ (được phân loại theo mô hình HCM thành các nhóm orbital định xứ, tiếp tuyến và xuyên tâm), dẫn đến việc toàn bộ các electron hóa trị đều được ghép đôi thành các vỏ đóng kín, dập tắt mômen từ tự do của Fe.

Kết luận

Luận án của NCS. Dương Văn Long đã nghiên cứu toàn diện cấu trúc hình học, cấu hình electron và tính thơm của các cụm boron nguyên chất và pha tạp ($\text{B-Si}$, $\text{B-Li}$, $\text{B-Fe}$, $\text{B}{70}$) bằng các phương pháp hóa học lượng tử tiên tiến. Công trình phân định ranh giới giữa quy tắc Hückel kinh điển và các mô hình hiện đại, làm sáng tỏ hiện tượng tự khóa ở cấu trúc dải ruy-băng, tổng quát hóa quy tắc đếm electron $(4n + 2m)$ cho mô hình đĩa tròn và đề xuất mô hình lai đĩa - nón. Các kết quả tính toán cũng chỉ ra tiềm năng ứng dụng thực tiễn của cụm $\text{B}{12}\text{Li}8$ trong lưu trữ hydro và các dây nano dẫn xuất từ $\text{B}{14}\text{FeLi}_2$ trong linh kiện quang điện tử.