Nghiên Cứu Hóa Học Lượng Tử Về Cấu Trúc Và Tính Aromatic Của Một Số Cụm Boron

Tài liệu nghiên cứu A quantum chemical research of structure and aromaticity of some boron clusters, tổng hợp lý thuyết và thực hành, cung cấp kiến thức chuyên sâu về .

• Tác giả: Duong Van Long
• Người hướng dẫn: Prof. Nguyen Minh Tho
• Trường học: Quy Nhon University
• Chuyên ngành: Theoretical and Physical Chemistry
• Đề tài: Nghiên Cứu Hóa Học Lượng Tử Về Cấu Trúc Và Tính Aromatic Của Một Số Cụm Boron
• Loại tài liệu: Doctoral Dissertation
• Năm xuất bản: 2023
• Địa điểm: Binh Dinh

Mục lục chi tiết

• DECLARATION
• ACKNOWLEDGEMENTS
• Table of Contents
• List of symbols and notations
• List of Figures
• List of Tables
• 1. GENERAL INTRODUCTION
• 1.1. Objectives and scope of the research
• 1.2. Novelty and scientific significance
• 1.3. Overview of the research
• 1.4. Objectives of the research
• 1.5. Search for lower-lying isomers
• 1.6. ELF – The electron localization function
• 1.7. Ring current maps
• 1.8. Bond order and net atomic charge
• 1.9. The Hückel and Baird rules
• 1.10. Hollow cylinder model
• 2. THEORETICAL BACKGROUNDS AND COMPUTATIONAL METHODS
• 2.1. Schrödinger equation
• 2.2. The Born–Oppenheimer Approximation
• 2.3. Ab initio computational method
• 2.4. The Hartree-Fock Method
• 2.5. Density Functional Theory
• 2.6. Benchmarking the functional and basis set in DFT
• 2.7. Post-Hartree-Fock methods
• 3. RESULTS AND DISCUSSION
• 3.1. The Hückel rule and the ribbon model: The cases of B2Si3q and B3Si2p clusters. Motivation for the study
• 3.2. The benchmarking tests
• 3.3. Ribbon aromaticity model versus the Hückel electron count
• 3.4. The disk aromaticity on the quasi-planar boron cluster B700/2-. Motivation of the study
• 3.5. The quasi-planar B700/2-
• 3.6. Disk model and electron count rule
• 3.7. Binary boron lithium clusters B12Lin with n = 1–14: the disk-cone model for the B12Li4 cluster
• 3.8. Motivation of the study
• 3.9. The growth pattern of B12Lin with n = 0 – 14
• 3.10. Relative stabilities of clusters
• 3.11. A mixed cone-disk model
• 3.12. B14FeLi2 and the hollow cylinder model
• 3.13. Motivation of the study
• 3.14. Stability of B14FeLi2 and its potential applications
• 4. GENERAL CONCLUSIONS AND FUTURE DIRECTIONS
• LIST OF PUBLICATIONS CONTRIBUTING TO THE DISSERTATION

Tóm tắt

I. Nghiên Cứu Hóa Học Lượng Tử Tổng Quan Cụm Boron Aromatic

Nghiên cứu về cụm Boron đã phát triển mạnh mẽ trong sáu thập kỷ qua, kể từ khi khái niệm "cụm nguyên tử kim loại" được định nghĩa. Hóa học lượng tử đóng vai trò then chốt trong việc tìm kiếm cấu trúc ổn định, giải thích tính ổn định và làm sáng tỏ các tính chất hóa học đặc trưng của chúng. Khái niệm tính thơm (aromaticity) nổi lên như một chủ đề liên quan mật thiết đến khoa học cụm, mặc dù đây là một khái niệm không hoàn toàn rõ ràng. Các mô hình định tính và định lượng khác nhau hỗ trợ lẫn nhau, tạo nên sự phức tạp trong việc hiểu và giải thích hiện tượng này. Bài viết này tập trung vào việc thiết lập các mô hình thích hợp cho tính thơm dựa trên hình dạng hình học, sử dụng các phương pháp toán học chặt chẽ để giải thích sự ổn định. Lý thuyết hàm mật độ (DFT) và các phương pháp tính toán lượng tử khác được sử dụng rộng rãi. Trích dẫn: 'Gần sáu thập kỷ đã trôi qua kể từ khi định nghĩa về thuật ngữ cluster được hình thành chính thức khi F. Cotton, vào năm 1964, lần đầu tiên tổng quát hóa một định nghĩa làm việc về “cluster nguyên tử kim loại” là “một nhóm hữu hạn các nguyên tử kim loại được giữ lại với nhau chủ yếu, hoặc ít nhất là ở một mức độ đáng kể, bởi các liên kết trực tiếp giữa các nguyên tử kim loại, ngay cả khi một số nguyên tử phi kim loại cũng có thể liên kết chặt chẽ với cluster”'.

1.1. Ứng Dụng Hóa Học Lượng Tử Trong Nghiên Cứu Cụm Boron

Các nghiên cứu thực nghiệm và lý thuyết về cluster đã phát triển mạnh mẽ, biến cluster nguyên tử thành một lĩnh vực khoa học đa ngành. Khoa học cluster chủ yếu liên quan đến việc tìm kiếm các cấu trúc phân tử ổn định của một nhóm nguyên tử. Tiếp theo đó là giải thích sự ổn định của các cấu trúc đó, cùng với việc làm sáng tỏ các tính chất hóa học đặc trưng của chúng và các ứng dụng tiềm năng. Các phương pháp tính toán lượng tử cung cấp một công cụ mạnh mẽ để dự đoán và giải thích các tính chất này.

1.2. Giới Thiệu Về Tính Thơm Aromaticity Trong Hóa Học Lượng Tử

Tính thơm (Aromaticity), mặc dù là một khái niệm cơ bản trong hóa học hiện đại, lại không có một định nghĩa rõ ràng. Điều này xuất phát từ sự tồn tại của nhiều mô hình định tính và thậm chí định lượng, vừa hỗ trợ vừa mâu thuẫn lẫn nhau. Mô hình Hückel là một trong những mô hình nổi tiếng nhất, ban đầu được hình thành cho các hydrocarbon phẳng. Tuy nhiên, việc áp dụng một cách máy móc quy tắc (4n+2) mà bỏ qua việc giải phương trình thế tục cho từng cấu trúc cụ thể dẫn đến nhiều sai lệch.

II. Thách Thức Mô Hình Aromatic Phù Hợp Cho Cụm Boron

Việc áp dụng mô hình Hückel một cách tùy tiện cho các hệ thống không phẳng và ba chiều đã dẫn đến những diễn giải sai lệch. Sự tiện lợi của quy tắc đếm electron (4n + 2) đã khiến nó bị lạm dụng đến mức bản chất và nguồn gốc của quy tắc thường bị lãng quên. Do đó, việc thiết lập các mô hình phù hợp cho tính thơm (aromaticity) dựa trên hình dạng hình học, sử dụng các phương pháp toán học chặt chẽ là rất quan trọng. Điều này bao gồm việc sử dụng các mô hình đĩa tròn, mô hình dải và mô hình hình trụ rỗng để nhấn mạnh sự khác biệt và tương đồng của các quy tắc đếm electron khi hình dạng hình học của các loài được xem xét khác biệt đáng kể so với vòng tròn phẳng của hydrocarbon hữu cơ. Cần có sự so sánh và đối chiếu giữa các phương pháp tính toán khác nhau để đảm bảo tính chính xác của kết quả. Trích dẫn: 'Sự tiện lợi của quy tắc đếm electron (4n + 2) đã khiến nó bị lạm dụng đến mức bản chất và nguồn gốc của quy tắc thường bị lãng quên và dẫn đến 1 diễn giải sai sót.'

2.1. Hạn Chế Của Mô Hình Hückel Trong Cụm Boron Không Phẳng

Mô hình Hückel, ban đầu được phát triển cho các phân tử hydrocarbon phẳng, thường được áp dụng không chính xác cho các cụm Boron có cấu trúc không phẳng hoặc ba chiều. Việc chỉ tập trung vào quy tắc (4n+2) mà bỏ qua các yếu tố hình học và cấu trúc phân tử đặc thù có thể dẫn đến những kết luận sai lệch về tính thơm (aromaticity).

2.2. Yêu Cầu Về Phương Pháp Toán Học Chặt Chẽ Hơn

Để đánh giá tính thơm (aromaticity) của cụm Boron một cách chính xác, cần thiết phải sử dụng các phương pháp toán học chặt chẽ hơn, dựa trên hình dạng hình học thực tế của cụm. Các mô hình như mô hình đĩa tròn, mô hình dải và mô hình hình trụ rỗng cung cấp các cách tiếp cận phù hợp hơn để mô tả sự phân bố electron và các tính chất hóa học của các cụm này.

2.3. Vai Trò của Phương Pháp Tính Toán Lượng Tử Chính Xác

Sử dụng các phương pháp tính toán lượng tử tiên tiến, như lý thuyết hàm mật độ (DFT) với các hàm trao đổi tương quan phù hợp và các bộ cơ sở đủ lớn, là rất quan trọng để đảm bảo tính chính xác của các kết quả nghiên cứu về cấu trúc phân tử và tính chất hóa học của cụm Boron.

III. Phương Pháp Nghiên Cứu DFT và Các Mô Hình Aromatic Đặc Thù

Luận án này sử dụng phương pháp tính toán lượng tử để xác định cấu trúc hình học, cấu hình electron và tính ổn định nhiệt động của một số cụm Boron và cụm Boron pha tạp. Các cụm Boron được khảo sát bao gồm B2Si3q và B3Si2p ở các trạng thái tích điện khác nhau, B700/2- trung tính và dianionic, B12Lin với n = 0 – 14 và B14FeLi2. Mô hình dải kết hợp với mô hình Hückel để giải thích các tính chất liên quan đến các cụm B2Si3q và B3Si2p. Tính ổn định của đồng phân quasi-planar của B700/2- và B12Li4 hình nón được hiểu rõ thông qua mô hình đĩa. Mô hình hình trụ rỗng đóng góp vào việc làm sáng tỏ các tính chất của B14FeLi2. Phần mềm tính toán hóa học chuyên dụng được sử dụng để thực hiện các tính toán này. Trích dẫn: 'Các cụm boron được khảo sát trong luận án bao gồm B2Si3q và B3Si2p ở các trạng thái tích điện khác nhau, B700/2- trung tính và dianionic, B12Lin với n = 0 – 14 và B14FeLi2.'

3.1. Áp Dụng Mô Hình Ribbon và Hückel cho B2Si3q và B3Si2p

Mô hình ribbon (dải) được sử dụng kết hợp với mô hình Hückel để giải thích các tính chất liên quan đến cụm Boron B2Si3q và B3Si2p. Mô hình này đặc biệt hữu ích trong việc mô tả sự liên kết hóa học và tính thơm (aromaticity) trong các cấu trúc có hình dạng dải hoặc chuỗi.

3.2. Mô Hình Đĩa Cho Tính Ổn Định Của B700 2 và B12Li4

Tính ổn định của đồng phân quasi-planar của B700/2- và B12Li4 hình nón được hiểu rõ thông qua mô hình đĩa. Mô hình này phù hợp để mô tả sự liên kết hóa học và tính thơm (aromaticity) trong các cấu trúc gần phẳng hoặc có dạng đĩa.

3.3. Sử Dụng Mô Hình Hình Trụ Rỗng Cho B14FeLi2

Mô hình hình trụ rỗng đóng góp vào việc làm sáng tỏ các tính chất của cụm Boron B14FeLi2. Mô hình này phù hợp để mô tả sự liên kết hóa học và các tính chất hóa học trong các cấu trúc có dạng ống hoặc hình trụ.

IV. Kết Quả Cấu Trúc và Tính Aromatic Độc Đáo Của Cụm Boron

Nghiên cứu này khám phá sự đa dạng hình học của các cụm Boron và các cụm Boron pha tạp, xác định các hình học ổn định của chúng. Các kết quả thu được tiết lộ sự đa dạng hình học và tính chất độc đáo. Mô hình orbital phân tử và năng lượng liên kết được sử dụng để phân tích sự ổn định của các cấu trúc phân tử. Kết quả nghiên cứu này cung cấp thông tin chi tiết về sự liên kết, tính thơm (aromaticity) và các tính chất hóa học của các cụm này, góp phần vào sự hiểu biết cơ bản về hóa học cụm và có thể mở đường cho các ứng dụng tiềm năng trong tương lai. Trích dẫn: 'Luận án này tiến hành khảo sát các cấu trúc của các cụm boron và các cụm boron pha tạp khác nhau và xác định hình học ổn định của chúng. Các kết quả thu được cho thấy sự đa dạng về hình học...'

4.1. Nghiên Cứu Chi Tiết Về Cấu Trúc Phân Tử B2Si3q và B3Si2p

Nghiên cứu chi tiết về cấu trúc phân tử của B2Si3q và B3Si2p, bao gồm việc xác định các đồng phân ổn định và phân tích sự phân bố electron. Các phương pháp tính toán lượng tử giúp xác định hình học phân tử và các tính chất hóa học của các cụm này.

4.2. Khám Phá Tính Aromatic Trong B700 2 B12Lin và B14FeLi2

Khám phá tính thơm (aromaticity) trong các cụm Boron B700/2-, B12Lin và B14FeLi2, sử dụng các mô hình đĩa tròn, mô hình dải và mô hình hình trụ rỗng để giải thích sự ổn định và sự phân bố electron. Phân tích orbital phân tử đóng vai trò quan trọng trong việc hiểu rõ các tính chất này.

4.3. Liên Hệ Giữa Cấu Trúc và Tính Chất Hóa Học

Phân tích mối liên hệ giữa cấu trúc phân tử và tính chất hóa học của các cụm Boron được nghiên cứu. Nghiên cứu này cung cấp thông tin chi tiết về cách cấu trúc hình học ảnh hưởng đến sự liên kết hóa học, tính ổn định và các tính chất hóa học khác của các cụm này.

V. Ứng Dụng Cụm Boron Tiềm Năng Trong Vật Liệu và Hóa Học

Các nghiên cứu về cụm Boron có tiềm năng ứng dụng rộng rãi trong nhiều lĩnh vực. Ví dụ, sự hiểu biết về tính chất hóa học và cấu trúc phân tử của cụm Boron có thể dẫn đến việc phát triển các vật liệu mới với các tính chất độc đáo. Các ứng dụng tiềm năng bao gồm chất xúc tác, vật liệu nano, và các thiết bị điện tử. Việc sử dụng phần mềm tính toán hóa học ngày càng trở nên quan trọng trong việc thiết kế và dự đoán các tính chất của vật liệu. Các kết quả từ các phương pháp tính toán lượng tử có thể giúp tối ưu hóa các tính chất hóa học của các cụm Boron cho các ứng dụng cụ thể. Trích dẫn: 'Các nghiên cứu về cụm Boron có tiềm năng ứng dụng rộng rãi trong nhiều lĩnh vực. Ví dụ, sự hiểu biết về tính chất hóa học và cấu trúc phân tử của cụm Boron có thể dẫn đến việc phát triển các vật liệu mới với các tính chất độc đáo.'

5.1. Cụm Boron Như Chất Xúc Tác Tiềm Năng

Cụm Boron có thể được sử dụng làm chất xúc tác trong nhiều phản ứng hóa học khác nhau do tính chất hóa học độc đáo của chúng. Các phương pháp tính toán lượng tử có thể giúp thiết kế các cụm Boron có hoạt tính xúc tác cao.

5.2. Ứng Dụng Trong Vật Liệu Nano

Các cấu trúc phân tử đặc biệt của cụm Boron làm cho chúng trở thành vật liệu nano đầy hứa hẹn. Các ứng dụng tiềm năng bao gồm cảm biến, vật liệu gia cường và các thiết bị lưu trữ năng lượng.

5.3. Cụm Boron Trong Thiết Bị Điện Tử

Cụm Boron có thể được sử dụng trong các thiết bị điện tử do khả năng dẫn điện và các tính chất hóa học độc đáo của chúng. Các nghiên cứu về cấu trúc phân tử và tính chất điện tử của cụm Boron là rất quan trọng cho ứng dụng này.

VI. Kết Luận Hướng Nghiên Cứu Tương Lai Về Cụm Boron Aromatic

Nghiên cứu này đã cung cấp những hiểu biết sâu sắc về cấu trúc phân tử, tính thơm (aromaticity) và tính chất hóa học của một số cụm Boron. Các kết quả này góp phần vào sự hiểu biết cơ bản về hóa học cụm và mở đường cho các nghiên cứu trong tương lai. Các hướng nghiên cứu tiếp theo có thể tập trung vào việc khám phá các loại cụm Boron mới, phát triển các mô hình lý thuyết tiên tiến hơn và khám phá các ứng dụng tiềm năng trong nhiều lĩnh vực. Sự kết hợp giữa hóa học lượng tử, phương pháp tính toán lượng tử, và các thí nghiệm thực nghiệm sẽ tiếp tục thúc đẩy lĩnh vực nghiên cứu cụm Boron phát triển. Trích dẫn: 'Nghiên cứu này đã cung cấp những hiểu biết sâu sắc về cấu trúc phân tử, tính thơm (aromaticity) và tính chất hóa học của một số cụm Boron. Các kết quả này góp phần vào sự hiểu biết cơ bản về hóa học cụm và mở đường cho các nghiên cứu trong tương lai.'

6.1. Phát Triển Các Mô Hình Lý Thuyết Tiên Tiến Hơn

Việc phát triển các mô hình lý thuyết tiên tiến hơn để mô tả sự liên kết hóa học và tính thơm (aromaticity) trong cụm Boron là một hướng nghiên cứu quan trọng. Các mô hình này cần phải có khả năng mô tả chính xác sự tương tác giữa các electron và các hạt nhân, cũng như ảnh hưởng của các yếu tố môi trường.

6.2. Khám Phá Ứng Dụng Mới Của Cụm Boron

Việc khám phá các ứng dụng mới của cụm Boron trong các lĩnh vực như vật liệu nano, xúc tác và điện tử là một hướng nghiên cứu đầy hứa hẹn. Các ứng dụng này có thể mang lại những lợi ích đáng kể cho xã hội.

6.3. Kết Hợp Tính Toán Lượng Tử và Thực Nghiệm

Sự kết hợp giữa hóa học lượng tử, các phương pháp tính toán lượng tử và các thí nghiệm thực nghiệm là rất quan trọng để thúc đẩy lĩnh vực nghiên cứu cụm Boron phát triển. Các kết quả từ các phương pháp tính toán lượng tử có thể giúp giải thích các kết quả thực nghiệm và hướng dẫn các thí nghiệm trong tương lai.

Nội dung chính

Tổng quan luận án

Khoa học về cụm nguyên tử (cluster science) tập trung vào việc xác định các cấu trúc bền nhiệt động của một nhóm hữu hạn nguyên tử, giải thích nguồn gốc độ bền và khảo sát các tính chất hóa lý cùng tiềm năng ứng dụng. Trong lĩnh vực này, khái niệm thơm (aromaticity) giữ vai trò then chốt nhưng thường bị lạm dụng khi quy tắc đếm electron cổ điển Hückel $(4n + 2)$ vốn xây dựng cho hydrocacbon phẳng mạch vòng bị áp dụng cơ học cho các hệ cấu trúc không phẳng, dạng dải băng, đĩa phẳng lớn hoặc dạng ống ba chiều, dẫn đến các diễn giải sai lệch về bản chất liên kết.

Nghiên cứu của luận án tiến sĩ chuyên ngành Hóa lý thuyết và Hóa lý (Mã số: 9440119) do nghiên cứu sinh Dương Văn Long 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 nhằm giải quyết các khoảng trống lý thuyết trên thông qua các mục tiêu cụ thể:

1. Khảo sát cấu trúc hình học, cấu hình electron và độ bền nhiệt động của một số hệ cụm boron nguyên chất và cụm boron pha tạp.
2. Thiết lập và kiểm chứng sự phù hợp của các mô hình aromaticity hiện đại (mô hình dải băng - ribbon model, mô hình đĩa tròn - disk model, mô hình hình trụ rỗng - hollow cylinder model) dựa trên việc giải phương trình sóng tương ứng với hình học cấu trúc.
3. Tổng quát hóa quy tắc đếm electron $(4n + 2m)$ cho cấu trúc dạng đĩa và đề xuất mô hình lai disk-cone cho các cụm dạng nón.
4. Đánh giá tiềm năng ứng dụng thực tế của các cụm boron bền, bao gồm khả năng lưu trữ hydro của cụm $B_{12}Li_8$ và tính chất quang điện tử của vật liệu dây nano phát triển từ đơn vị $B_{14}FeLi_2$.

Đối tượng nghiên cứu của luận án bao gồm các hệ cụm boron và boron pha tạp: cụm $B_2Si_3^q$ ($q = -2$ đến $+2$), $B_3Si_2^p$ ($p = -3$ đến $+1$), cụm boron bán phẳng trung hòa và dianion $B_{70}^{0/2-}$, dãy cụm nhị nguyên boron-lithium $B_{12}Li_n$ ($n = 0 - 14$) và cụm pha tạp kim loại chuyển tiếp $B_{14}FeLi_2$.

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 chính thức định nghĩa từ năm 1964. Tiếp đó, các nghiên cứu cấu trúc boron phát triển mạnh mẽ do nguyên tử boron có đặc tính thiếu hụt electron (electron-deficient), tạo nên nhiều dạng hình học đặc biệt từ phẳng, bán phẳng (quasi-planar), lồng icosahedral, fullerene đến ống nano (tubular).

Về mặt lý thuyết aromaticity:

• Quy tắc Wade-Mingos dự đoán độ bền của borane dựa trên số cặp electron khung $(2n + 2)$.
• Quy tắc Hückel $(4n + 2)$ cho hệ đơn vòng liên hợp đóng vỏ ở trạng thái singlet và quy tắc Baird $(4n)$ cho hệ mở vỏ ở trạng thái triplet được thiết lập cho các hydrocacbon thơm cổ điển.
• Năm 2012, Tai và cộng sự lần đầu tiên áp dụng mô hình hạt chuyển động trong hộp thế đĩa tròn 2D để giải thích độ bền cho các cụm phẳng; mô hình này tiếp tục được kiểm chứng trên các cụm $B_{19}^-$, $B_{20}^{2-}$, $B_{30}$, $B_{50}$.
• Đối với các cấu trúc ống hai vòng (double-ring - DR) như $B_{20}$ trung hòa hay $B_{14}^{2+}$, quy tắc $(4n + 2)$ vẫn được ghi nhận; tuy nhiên với các cấu trúc ba vòng (triplet-ring - TR) như $B_{27}^+$ và $B_{42}$, quy tắc $(4n + 2)$ không còn phù hợp mà chuyển sang quy tắc $(4n + 2m)$, dẫn đến sự hình thành mô hình hình trụ rỗng (HCM).
• Đối với cấu trúc dạng dải băng hai chuỗi (double-chain ribbon) kéo dài một chiều như $B_{10}H_2^{2-}$, $B_{14}H_2^{2-}$, $B_{10}Si_2^{2-}$ và $B_{12}Si_2^{2-}$, các nghiên cứu trước đây (2017) đã đề xuất mô hình hộp thế một chiều (1D box) giải thích cấu hình bền đặc trưng $\pi^{2(n+1)}\sigma^{2n}$.

Vị trí của luận án là mở rộng và hoàn thiện các mô hình aromaticity phi cổ điển nói trên, làm rõ ranh giới giữa mô hình Hückel và mô hình ribbon, mở rộng quy tắc đếm electron cho đĩa phẳng $B_{70}^{0/2-}$, đề xuất mô hình đĩa-nón (disk-cone) cho $B_{12}Li_4$, và phân tích tương tác dập tắt từ tính của Fe trong cụm hình trụ $B_{14}FeLi_2$.

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

Luận án kết hợp các phương pháp hóa học lượng tử tính toán từ cơ bản (ab initio) đến phiếm hàm mật độ (DFT) và các công cụ phân tích hàm sóng topo:

Phương pháp tìm kiếm cấu trúc đáy và tối ưu hóa hình học

• Thuật toán tìm kiếm ngẫu nhiên (stochastic search) kết hợp giữa phương pháp dịch chuyển ngẫu nhiên (random kick) và thuật toán di truyền (Genetic Algorithm - GA / MEGA) để tạo ra lượng lớn các cấu trúc ban đầu.
• Tối ưu hóa sơ bộ ở mức DFT/TPSSh/6-31G(d).
• Tối ưu hóa chuẩn xác và tính toán tần số dao động điều hòa (harmonic vibrational frequencies) ở mức TPSSh/6-311+G(d) hoặc def2-TZVP bằng phần mềm Gaussian 16 nhằm đảm bảo cấu trúc thu được là điểm cực tiểu thực trên mặt thế năng (không có tần số ảo, $N_{imag} = 0$).
• Tính toán đơn điểm năng lượng cao cấp bằng phương pháp Coupled Cluster hiệu chỉnh kích thích ba lần: (U)CCSD(T)/cc-PVTZ, aug-cc-pVTZ hoặc ngoại suy tập hàm cơ sở đầy đủ (CBS).

Các công cụ phân tích liên kết và tính thơm

• Hàm định vị electron (ELF): Phân tích hàm định vị electron $\eta(\mathbf{r})$ dựa trên mật độ động năng Pauli ($D_P$) và mật độ động năng Thomas-Fermi ($D_h$): $$\eta(\mathbf{r}) = \frac{1}{1 + (D_P / D_h)^2}$$ Phân tách ELF thành các thành phần định vị cục bộ $\sigma$ ($\text{ELF}{\sigma l}$), giải tỏa $\sigma$ ($\text{ELF}{\sigma d}$) và giải tỏa $\pi$ ($\text{ELF}_\pi$) bằng chương trình Dgrid-5.
• Bản đồ dòng điện cảm ứng từ (Ring current maps): Sử dụng phương pháp CTOCD-DZ2 thông qua phần mềm SYSMOIC và phương pháp dị hướng của mật độ dòng cảm ứng (ACID) để xác định tính diatropic (thơm, dòng thuận chiều kim đồng hồ) hoặc paratropic (chống thơm, dòng ngược chiều kim đồng hồ).
• Phân tích điện tích nguyên tử thực và bậc liên kết: Sử dụng phương pháp phân chia mật độ DDEC6 để tính điện tích nguyên tử thực (NAC), bậc liên kết (BO) và tổng bậc liên kết trung bình (ASBO).

Phương pháp / Phần mềm	Mục đích sử dụng trong luận án
Gaussian 16	Tối ưu hóa cấu trúc hình học, tính tần số dao động, tính năng lượng DFT và CCSD(T).
GA / Random Kick	Tạo lập và tìm kiếm quét toàn không gian các cấu trúc đồng phân năng lượng thấp.
Dgrid-5	Tính toán và biểu diễn mặt đẳng trị của hàm ELF toàn phần và các thành phần $\text{ELF}\sigma$, $\text{ELF}\pi$.
SYSMOIC / ACID	Tính toán tensor mật độ dòng điện cảm ứng từ, vẽ bản đồ dòng vòng (ring current).
DDEC6	Phân tích mật độ electron để tính điện tích nguyên tử thực (NAC) và bậc liên kết (BO).

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

Chương 1: Giới thiệu chung và các mô hình aromaticity (General Introduction)

Chương này thiết lập nền tảng toán học cho các mô hình aromaticity dựa trên phương trình Schrödinger:

1. Mô hình dải băng (Ribbon model): Xem electron giải tỏa chuyển động tự do trong hộp thế một chiều chiều dài $l$: $$\psi_n(x) = \sqrt{\frac{2}{l}} \sin\left(\frac{n\pi x}{l}\right), \quad E_n = \frac{n^2 h^2}{8 m l^2}$$ Mối liên hệ năng lượng obitan phân tử chuẩn tắc (CMO) từ DFT và mô hình 1D được tối thiểu hóa qua tổng bình phương độ lệch (RSS): $$RSS_\sigma = \sqrt{\sum_{i=1}^{N_\sigma} (E_{DFT, n\sigma} - E_{n\sigma})^2}$$ Độ dài hiệu dụng $l_\pi$ và $l_\sigma$ tìm được phù hợp chặt chẽ với kích thước hình học của cụm (ví dụ $B_{14}H_2^{2-}$ có $l_\pi = 12.3\text{ \AA}$, $l_\sigma = 11.2\text{ \AA}$ so với khoảng cách hai nguyên tử B xa nhất là $10.6\text{ \AA}$). Hiện tượng tự khóa (self-locking phenomenon): Các liên kết ở hai đầu dải băng (như B1–B14 và B7–B8 trong $B_{14}H_2^{2-}$) bị co ngắn đáng kể (~$1.58\text{ \AA}$), tạo ra hàng rào thế năng Coulomb cao giam giữ các electron giải tỏa chuyển động tuần hoàn dọc theo khung dải băng.

2. Mô hình đĩa tròn (Disk model): Giải phương trình vi phân trong tọa độ cực $(r, \theta)$ cho hạt chuyển động trên mặt phẳng giới hạn bởi bán kính $R$:

• Hàm góc: $\Theta_m(\theta) = \frac{1}{\sqrt{2\pi}} e^{i m \theta}$ với $m = 0, \pm 1, \pm 2, \dots$ (trạng thái $m=0$ không suy biến, $m \ge 1$ suy biến bậc hai).
• Hàm xuyên tâm: $\mathcal{R}_m(r) = J_m(\alpha r)$ (hàm Bessel loại 1 bậc $m$).
• Điều kiện biên $\mathcal{R}m(R) = J_m(\alpha R) = 0$ xác định các mức năng lượng thông qua nghiệm $x{m,n} = \alpha_{m,n} R$.

3. Mô hình hình trụ rỗng (Hollow Cylinder Model - HCM): Áp dụng cho các cụm dạng ống có bán kính trong $R_0 = R - r$, bán kính ngoài $R_1 = R + r$ và chiều cao $L$.

Chương 2: Cơ sở lý thuyết và phương pháp tính toán (Theoretical Backgrounds and Computational Methods)

Trình bày chi tiết lý thuyết lượng tử: phương trình Schrödinger, phép gần đúng Born-Oppenheimer, phương pháp Hartree-Fock (HF), lý thuyết phiếm hàm mật độ (DFT), và các phương pháp hậu Hartree-Fock (MPn, CI, CASSCF, CC, CCSD(T)). Kết quả kiểm chuẩn (benchmark) cho thấy:

• Phiếm hàm meta-GGA TPSSh kết hợp tập hàm cơ sở 6-311+G(d) cho độ chính xác cao nhất về mặt năng lượng và tối ưu hóa hình học đối với các hệ chứa đồng thời nguyên tử B và Si khi so sánh với tính toán mức cao (U)CCSD(T)/CBS.
• Phiếm hàm B3LYP mô phỏng phổ quang electron tách hạt (PES - VDE và ADE) và phổ quang phổ ion hóa hai màu hồng ngoại - tử ngoại (IR-UV2CI) sát nhất với dữ liệu thực nghiệm đối chứng của cụm $B_2Si_3^-$.

Chương 3: Kết quả và thảo luận (Results and Discussion)

1. Quy tắc Hückel và mô hình nanoribbon: Các cụm B2Si3^q và B3Si2^p

• Khảo sát các trạng thái điện tích $q = -2, -1, 0, +1, +2$ cho $B_2Si_3^q$ và $p = -3, -2, -1, 0, +1$ cho $B_3Si_2^p$.
• Xác định sự cạnh tranh giữa hai dạng đồng phân chính: dạng dải băng phẳng (Ribbon - R) và dạng vòng phẳng Hückel (Hückel - H).
• Đối với cấu hình electron $\pi^{2(n+1)}\sigma^{2n}$, cụm thể hiện tính thơm dải băng bền vững (ví dụ dianion $B_{14}H_2^{2-}$ với cấu hình $\pi^6\sigma^4$). Ngược lại, cấu hình $\pi^{2n}\sigma^{2n}$ dẫn đến trạng thái chống thơm dải băng (ribbon antiaromaticity) do sự đứt gãy các hốc ELF giữa các khối hình thoi (rhombus unit), làm xuất hiện tần số dao động ảo (như trạng thái singlet của $B_{12}H_2^{2-}$ có tần số $-178\text{ cm}^{-1}$).
• Trạng thái triplet của $B_{12}H_2^{2-}$ với cấu hình $\pi^{2n+1}\sigma^{2n-1}$ có các hốc $\text{ELF}\alpha$ và $\text{ELF}\beta$ bù trừ liên tục, phục hồi độ bền và tạo nên tính thơm triplet dải băng (triplet ribbon aromaticity) với tần số dao động dương nhỏ nhất $68\text{ cm}^{-1}$.

2. Tính chất disk aromaticity trên cụm boron bán phẳng B70^(0/2-)

• Cấu trúc bán phẳng của cụm 70 nguyên tử boron ($B_{70}$) được tạo lập thành công bằng thuật toán bước cóc topo (topological leapfrog algorithm) từ cấu trúc ban đầu $B_{16}$ có 13 đỉnh.
• Khảo sát phân bố năng lượng chỉ ra trạng thái đáy của cấu trúc bán phẳng $B_{70}$ trung hòa là trạng thái triplet.
• Cụm dianion $B_{70}^{2-}$ có cấu hình vỏ electron đóng đạt độ bền nhiệt động vượt trội.
• Sự tương ứng giữa các obitan phân tử $\pi$-MO tính toán từ DFT của $B_{70}^{2-}$ và các mức năng lượng giải từ hàm Bessel của mô hình đĩa tròn xác nhận cấu trúc tuân theo quy tắc đếm electron thơm dạng đĩa tổng quát $(4n + 2m)$ với $m$ là số mức năng lượng không suy biến ($m=0$). Bản đồ mật độ dòng từ xác nhận dòng thuận từ mạnh mẽ (diatropic current) phân bố trên bề mặt.

3. Các cụm nhị nguyên boron-lithium B12Li_n (n = 1–14) và mô hình disk-cone cho cụm B12Li4

• Khảo sát quá trình pha tạp Li vào khung $B_{12}$ từ $n = 1$ đến $14$ ở mức (U)CCSD(T)/cc-PVTZ + ZPE//TPSSh/6-311+G(d).
• Phân tích năng lượng liên kết trên từng nguyên tử, khoảng cách HOMO-LUMO, năng lượng phân ly và sai phân năng lượng bậc hai xác định các cấu trúc bền đặc biệt: đồng phân 4C ($B_{12}Li_4$), 8A ($B_{12}Li_8$) và 14C-Th ($B_{12}Li_{14}$).
• Cụm $B_{12}Li_4$ có dạng nón (cone-shaped) đẳng electron với $B_{13}Li$. Luận án đề xuất mô hình thơm lai đĩa-nón (disk-cone model) để giải thích sự phân bố obitan hóa trị và dòng vòng cảm ứng từ bên trong khung nón.
• Đánh giá khả năng hấp phụ hydro: Cấu trúc 8A ($B_{12}Li_8$) có khả năng hấp phụ liên tiếp các phân tử $H_2$ tạo thành các phức $B_{12}Li_8\text{--}8H_2$, $B_{12}Li_8\text{--}16H_2$, $B_{12}Li_8\text{--}24H_2$, $B_{12}Li_8\text{--}32H_2$ và tối đa $B_{12}Li_8\text{--}40H_2$ (ở mức tính toán $\omega\text{B97XD/6-311++G(2d,2p)}$), cho thấy tiềm năng làm vật liệu lưu trữ hydro.

4. Cụm B14FeLi2 và mô hình hình trụ rỗng (Hollow Cylinder Model)

• Cụm $B_{14}FeLi_2$ có cấu trúc đối xứng dạng con quay (teetotum structure) với trục Li–Fe–Li nằm dọc theo tâm khung ống hình trụ $B_{14}$.
• Phân tích tương tác giữa obitan $d$ của nguyên tử Fe với các obitan $\sigma$ và $\pi$ của khung trụ rỗng $B_{14}$ giải thích nguyên nhân từ tính của Fe bị dập tắt hoàn toàn (closed-shell singlet).
• Thiết kế mô hình dây dẫn nano từ đơn vị cơ sở: chuỗi $[Li\text{-}B_{14}Fe\text{-}Li]\text{-}[B_{14}Fe]\text{-}[Li\text{-}B_{14}Fe\text{-}Li]$ và các dây nano $B_{28}Fe_2Li_2Mg$, $B_{42}Fe_3Li_2Mg_2$, $B_{56}Fe_4Li_2Mg_3$, $B_{70}Fe_4Mg_4Li_2$ dạng lăng trụ và đối lăng trụ (prism / antiprism forms) thể hiện dải hấp thụ quang học UV-Vis phù hợp cho ứng dụng vật liệu quang điện.

Chương 4: Kết luận chung và định hướng phát triển (General Conclusions and Future Directions)

Tổng kết các kết quả tính toán cơ bản và hệ thống hóa các quy tắc đếm electron cho từng họ cấu trúc cụm boron.

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

1. Phân biệt và mở rộng mô hình dải băng (Ribbon model): Luận án đã làm rõ sự khác biệt giữa mô hình Hückel cổ điển và mô hình ribbon; chỉ ra rằng cấu hình electron $\pi^{2(n+1)}\sigma^{2n}$ quyết định tính thơm và độ bền của dải băng, giải thích hiện tượng tự khóa (self-locking) ở hai đầu mút và chứng minh trạng thái triplet thơm dải băng với cấu hình $\pi^{2n+1}\sigma^{2n-1}$.
2. Kiểm chuẩn phương pháp tính toán: Xác lập cơ sở lựa chọn phiếm hàm TPSSh cho việc tối ưu hóa hình học và xác định năng lượng tương đối của các cụm boron pha tạp silic, và phiếm hàm B3LYP cho việc mô phỏng phổ quang electron (PES) và phổ dao động ion hóa (IR-UV2CI).
3. Phát hiện trạng thái đáy triplet của B70 bán phẳng và tổng quát hóa quy tắc đĩa: Chứng minh cụm bán phẳng $B_{70}$ có trạng thái đáy triplet và dianion $B_{70}^{2-}$ có độ bền nhiệt động cao; tổng quát hóa quy tắc đếm electron $(4n + 2m)$ cho các cấu trúc đĩa tròn thông qua phương trình sóng Bessel.
4. Đề xuất mô hình thơm disk-cone và đánh giá vật liệu lưu trữ hydro: Xây dựng mô hình lai disk-cone giải thích cấu trúc nón $B_{12}Li_4$; phát hiện cụm $B_{12}Li_8$ có khả năng hấp phụ tối đa 40 phân tử $H_2$ ($B_{12}Li_8\text{--}40H_2$).
5. Ứng dụng mô hình hình trụ rỗng (HCM) trên B14FeLi2: Giải thích bản chất dập tắt từ tính của Fe khi nội hấp trong lồng trụ $B_{14}$ và thiết kế thành công các chuỗi dây nano đơn chiều ứng dụng trong linh kiện quang điện tử.

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

• Hạn chế: Các kết quả khảo sát dừng lại ở quy mô tính toán mô phỏng lượng tử trên hệ cô lập trong chân không (khí trơ); các dự đoán về khả năng hấp phụ hydro và cấu trúc dây nano chưa có điều kiện kiểm chứng trực tiếp bằng tổng hợp thực nghiệm hóa học do độ phức tạp trong việc kiểm soát phản ứng của cụm boron.
• Hướng nghiên cứu tiếp: Tiếp tục mở rộng khảo sát pha tạp các kim loại chuyển tiếp khác vào khung ống boron kích thước lớn hơn; nghiên cứu động học phản ứng giải hấp hydro của cụm $B_{12}Li_8$ ở các điều kiện nhiệt độ và áp suất khác nhau; tối ưu hóa tính chất truyền dẫn điện tử của các cấu trúc dây nano phát triển từ $B_{14}FeLi_2$.

Giá trị tham khảo

Luận án là tài liệu tham khảo chuyên sâu đối với:

• Nghiên cứu sinh, học viên cao học và giảng viên ngành Hóa lý thuyết, Hóa lý, Vật lý chất rắn và Khoa học vật liệu tính toán, đặc biệt trong việc ứng dụng các phương pháp phân tích hàm sóng topo (ELF, ACID, DDEC6) và giải các bài toán hạt trong hộp thế hình học đặc thù.
• Các nhóm nghiên cứu thực nghiệm trong lĩnh vực vật liệu lưu trữ năng lượng (pin lithium, vật liệu hấp phụ khí hydro) và công nghệ nano quang điện tử nhằm định hướng tổng hợp các cụm nguyên tử bền.

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

1. Vì sao quy tắc Hückel cổ điển không thể giải thích thỏa đáng độ bền của dải băng boron?

Quy tắc Hückel $(4n + 2)$ được xây dựng dựa trên việc giải phương trình thế vòng tròn đối xứng cho các obitan $p_z$ đơn lẻ trong hệ hydrocacbon phẳng mạch vòng liên hợp. Ngược lại, cấu trúc dải băng boron là hệ hai chuỗi kéo dài một chiều, trong đó cả hệ electron giải tỏa $\pi$ và $\sigma$ đều tham gia quyết định độ bền với cấu hình đóng vỏ tối ưu là $\pi^{2(n+1)}\sigma^{2n}$, tương ứng với nghiệm của phương trình sóng cho hạt chuyển động trong hộp thế một chiều (1D box model).

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

Hiện tượng tự khóa là sự co ngắn đáng kể của các liên kết boron–boron ở hai đầu tận cùng của dải băng (ví dụ liên kết B1–B14 và B7–B8 trong $B_{14}H_2^{2-}$ chỉ dài khoảng $1.58\text{ \AA}$, ngắn hơn nhiều so với các liên kết phía trong). Sự co ngắn này làm tăng lực đẩy Coulomb giữa các hạt nhân ở biên, tạo ra một hàng rào thế năng cao phản xạ electron, giam giữ các electron $\pi$ và $\sigma$ giải tỏa chuyển động ổn định trong một hộp thế một chiều dọc theo thân dải băng.

3. Quy tắc đếm electron $(4n + 2m)$ trong mô hình đĩa tròn (Disk model) có ý nghĩa gì?

Trong mô hình đĩa tròn 2D, phương trình sóng Bessel cho các mức năng lượng có hai nhóm: nhóm mức không suy biến ứng với số lượng tử góc $m = 0$ (mỗi mức chứa tối đa 2 electron) và nhóm mức suy biến bậc hai ứng với $m \ge 1$ (mỗi cặp mức chứa tối đa 4 electron). Do đó, tổng số electron để điền đầy hoàn toàn các vỏ electron tạo nên hệ thơm đóng vỏ bền vững tuân theo quy tắc $(4n + 2m)$, trong đó $n$ là số cặp mức suy biến đã điền đầy và $m$ là số mức không suy biến đã điền đầy.

4. Cụm nhị nguyên B12Li_n nào có triển vọng ứng dụng lưu trữ hydro cao nhất và dung lượng đạt được là bao nhiêu?

Cụm $B_{12}Li_8$ (đồng phân 8A) là ứng viên triển vọng nhất trong dãy $B_{12}Li_n$ ($n = 1 - 14$). Kết quả tối ưu hóa hình học ở mức $\omega\text{B97XD/6-311++G(2d,2p)}$ cho thấy cụm $B_{12}Li_8$ có khả năng hấp phụ ổn định tối đa 40 phân tử $H_2$ tạo thành cấu trúc phức bền $B_{12}Li_8\text{--}40H_2$.

Kết luận

Luận án của nghiên cứu sinh Dương Văn Long đã nghiên cứu hệ thống 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 ($B_2Si_3^q$, $B_3Si_2^p$, $B_{70}^{0/2-}$, $B_{12}Li_n$, $B_{14}FeLi_2$) bằng các phương pháp hóa học lượng tử chuẩn xác cao. Công trình đã phân biệt rõ ràng phạm vi áp dụng của mô hình Hückel và mô hình dải băng, tổng quát hóa quy tắc đếm electron $(4n + 2m)$ cho mô hình đĩa tròn trên hệ bán phẳng $B_{70}^{2-}$, đề xuất mô hình lai disk-cone và kiểm chứng mô hình hình trụ rỗng. 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 $B_{12}Li_8$ trong công nghệ lưu trữ hydro và các cấu trúc dây nano phát triển từ $B_{14}FeLi_2$ trong lĩnh vực vật liệu quang điện tử.