Nghiên cứu khả năng hấp thụ một số phân tử hữu cơ chứa vòng benzen lên bề mặt vật liệu tio2

BỘ GIÁO DỤC VÀ ĐÀO TẠO
TRƯỜNG ĐẠI HỌC QUY NHƠN

HUỲNH THỊ MỸ PHÚC

NGHIÊN CỨU KHẢ NĂNG HẤP PHỤ MỘT SỐ PHÂN
TỬ HỮU CƠ CHỨA VÒNG BENZEN LÊN BỀ MẶT
VẬT LIỆU TiO2

LUẬN VĂN THẠC SĨ HĨA HỌC

Bình Định - Năm 2019


BỘ GIÁO DỤC VÀ ĐÀO TẠO
TRƯỜNG ĐẠI HỌC QUY NHƠN

HUỲNH THỊ MỸ PHÚC

NGHIÊN CỨU KHẢ NĂNG HẤP PHỤ MỘT SỐ PHÂN
TỬ HỮU CƠ CHỨA VÒNG BENZEN LÊN BỀ MẶT
VẬT LIỆU TiO2

Chun ngành :
Mã số

:

Hóa vơ cơ
8440113

Người hướng dẫn: PGS.TS. NGUYỄN TIẾN TRUNG


LỜI CAM ĐOAN
Tôi xin cam đoan luận văn này là cơng trình nghiên cứu của riêng tơi.
Các số liệu và kết quả trong luận văn là trung thực và chưa từng được cơng bố
trong các cơng trình khoa học khác. Việc tham khảo các nguồn tài liệu được
trích dẫn và ghi nguồn theo đúng quy định.

Tác giả
Huỳnh Thị Mỹ Phúc


LỜI CẢM ƠN

Lời đầu tiên cho em xin được bày tỏ lòng biết ơn sâu sắc đến Thầy
PGS.TS. Nguyễn Tiến Trung, người đã ln tận tình hướng dẫn, chỉ bảo, động
viên em trong suốt q trình nghiên cứu để hồn thành luận văn này.
Em xin trân trọng cảm ơn quý Thầy Cơ giáo trong Khoa Hóa, Trường
Đại học Quy Nhơn đã trang bị cho em những kiến thức khoa học bổ ích để em
có thể hồn thành khóa luận tốt nghiệp.
Em cũng xin chân thành cảm ơn NCS. Nguyễn Ngọc Trí đã ln
nhiệt tình giúp đỡ, hướng dẫn tận tình trong suốt q trình thực hiện đề tài
này.
Ngồi ra, tơi xin chân thành cảm ơn các anh, chị, em học viên cao học
K20 và các bạn sinh viên đang làm NCKH tại Phịng Thí nghiệm Hóa học
tính tốn và Mơ phỏng (LCCM), Trường Đại học Quy Nhơn luôn luôn giúp
đỡ, hướng dẫn và đóng góp ý kiến để tơi hồn thành tốt luận văn.
Cuối cùng, tôi xin dành lời cảm ơn gia đình ln ủng hộ, động viên,
giúp đỡ và tạo điều kiện để tơi có thêm động lực và thời gian để hồn thành

luận văn này.
Bình Định, ngày 07 tháng 8 năm
2019
Học viên
Huỳnh Thị Mỹ Phúc


1

MỞ ĐẦU
1. Lý do chọn đề tài
Các dẫn xuất chứa vòng benzen như benzaldehyde (C 6H5CHO), benzoic
acid (C6H5COOH), aniline (C6H5NH2), phenol (C6H5OH), benzenesulfonic
acid (C6H5SO3H) cùng nhiều các hợp chất hữu cơ khác là một trong những
thành phần quan trọng trong cơng nghiệp như: ngành hóa chất, phân bón, khai
thác và chế biến dầu mỏ, ngành dệt nhuộm, giấy, nhựa, chất tẩy rửa, ngành
chế biến thực phẩm,... ngành cơng nghiệp hóa chất khác [5], [6], [7], [8], [9].
Tuy nhiên, dư lượng lớn các loại hóa chất hữu cơ trên cùng với một lượng lớn
các chất gây ô nhiễm môi trường khác như nước thải axit (hoặc kiềm), nhiều
chất lơ lửng và dư lượng hóa chất hữu cơ gây hại có trong nước thải (các hợp
chất dẫn xuất của benzen có hàm lượng cao) gây nên các tác động xấu đến đời
sống của các động vật và con người, gây ra những biến chứng khôn lường [9],
[10], [11]. Điều đáng quan tâm nhất là con người vừa là nhân tố gây nên sự ô
nhiễm môi trường vừa nạn nhân chủ yếu của sự ơ nhiễm đó.
Đáng chú ý, các nghiên cứu về thực nghiệm và mô phỏng cho thấy bề
mặt vật liệu Titanium dioxide (TiO2) có nhiều ưu điểm như năng lượng vùng
cấm rộng, độ ổn định nhiệt cao, tương đối trơ về mặt hóa học. Do đó, TiO 2 đã
và đang được nghiên cứu, ứng dụng rộng rãi trong các lĩnh vực khoa học và
công nghệ [12], [13], [14], [15]. Trong số các dạng tồn tại của TiO 2 (rutile,
anatase, brookite), rutile là pha bền nhất của nó và bền ở nhiệt độ cao. Pha

anatase mặc dù kém bền hơn nhưng có độ hoạt động và khả năng xúc tác
quang cao hơn pha rutile, trong dạng hạt nano, cấu trúc anatase hoạt động hơn
cấu trúc rutile và brookite [16]. Khi xem xét các bề mặt hoạt động của hai pha
rutile và anatase thấy rằng bề mặt rutile-110 và bề mặt anatase-101 là các bề
mặt bền nhất và được nghiên cứu nhiều nhất [13], [14], [17], [18], [19]. Các


2

nghiên cứu chỉ ra hai bề mặt này có khả năng tương tác tốt với chất hữu cơ
chứa các nhóm chức –COOH [18], [20], [21], [22], [23] nhóm –NO 2 [18],
nhóm –CHO [24], [25] ... Tuy nhiên bản chất các tương tác giữa các phân tử
hữu cơ và bề mặt vật liệu TiO2 chưa được quan tâm đúng mức và vai trò của
chúng đến độ bền các phức thu được quá trình hấp phụ trên các bề mặt vật
liệu này cịn chưa được làm rõ.
Thêm vào đó, để tìm kiếm các vật liệu có khả năng hấp phụ tốt các chất
hữu cơ gây ô nhiễm môi trường chúng ta cần phải hiểu được sự tương tác giữa
các phân tử đó với bề mặt vật liệu và sự bền hóa của cấu trúc hình thành. Vì
vậy, việc hiểu biết về bản chất, vai trò của các tương tác giữa các phân tử hữu
cơ gây ơ nhiễm mơi trường chứa vịng benzen như benzaldehyde benzoic
acid, aniline, phenol, benzenesulfonic acid với bề mặt vật liệu TiO 2 là cơ sở
cho việc hiểu rõ hơn tương tác của các phân tử có nhóm chức khác nhau với
bề mặt các vật liệu nói chung. Đặc biệt, tính tốn về hóa học lượng tử đã và
đang giúp các nhà khoa học hiểu rõ hơn về các hiện tượng trong tự nhiên, đặc
biệt là sự hấp phụ các phân tử lên các bề mặt vật liệu [26], [27], [28]. Và theo
hiểu biết của chúng tôi, vẫn chưa có cơng trình nào nghiên cứu về khả năng
hấp phụ của các dẫn xuất vòng benzen lên các bề mặt vật liệu TiO 2 ở cấp độ
phân tử.
Từ những lý do nêu trên, chúng tôi chọn đề tài: “Nghiên cứu khả năng
hấp phụ một số phân tử hữu cơ chứa vòng benzen lên bề mặt vật liệu

TiO2”. Các kết quả thu được từ việc nghiên cứu này hi vọng sẽ đưa ra các
nhận định về khả năng hấp phụ của các hợp chất có nhóm chức khác nhau
trên bề mặt vật liệu, và từ đó là cơ sở để tìm kiếm vật liệu phù hợp để hấp phụ
những dòng hợp chất hữu cơ khác nhau, góp phần thiết thực để loại bỏ các
hợp chất hữu cơ độc hại bằng phương pháp hấp phụ.


3

2. Tổng quan tài liệu và tình hình nghiên cứu đề tài
Những năm gần đây, việc nghiên cứu thực nghiệm và lý thuyết để xử lý
các dư lượng các chất hữu cơ gây ô trong môi trường nước cũng như môi
trường đất và trong cá thể sinh vật được tiến hành với các phương pháp khác
nhau [29], [30]. Kết quả cho thấy các phương pháp oxi hóa, lọc cơng nghệ
cao, hấp phụ trên các vật liệu than hoạt tính, than sinh học, vật liệu khoáng, sử
dụng xúc tác quang, … mang lại những thành cơng nhất định, góp phần loại
bỏ các chất gây ô nhiễm môi trường, làm sạch nguồn nước, chuyển hóa các
chất ơ nhiễm thành các chất thân thiện môi trường [29], [30], [31], [32]. Trong
số các phương pháp đó, phương pháp hấp phụ bằng vật liệu cũng đang được
nhiều nhà khoa học quan tâm nghiên cứu [19], [22], [24], [25]. Vật liệu TiO 2
là một oxit kim loại bán dẫn, có tính chất vật lý và hóa học đặc biệt như diện
tích bề mặt lớn, hoạt tính cao, ổn định, khơng độc hại, chi phí sản xuất thấp ...
[33]. Do nhiều ưu điểm vượt trội nên vật liệu TiO 2 đang được ứng dụng nhiều
trong xúc tác quang, hấp phụ xử lý loại bỏ các chất gây ô nhiễm, làm sạch
nguồn nước, chuyển hóa các chất khó phân hủy thành các chất thân thiện môi
trường [34], [35], [36], [37].
Trong nhiều năm qua, nhiều cơng trình đã được công bố về khả năng hấp
phụ của TiO2 với các chất hữu cơ: benzoic acid, acetaldehyde, acetic acid,
methylacetate,... Năm 2006, Maria L. Sushko cùng các cộng sự [38] đã sử
dụng một mơ hình cluster đơn giản để nghiên cứu sự hấp phụ của một số

nhóm methyl, benzyl và cacboxylic trên bề mặt rutile-TiO 2 (110). Tuy nhiên,
kết quả thu được chưa mơ tả sự hình thành, cắt đứt các liên cộng hóa trị, sự
truyền tải điện tích giữa phân tử và bề mặt trong quá trình hấp phụ. Năm
2012, Grinter và các cộng sự [22] thực hiện khảo sát sự hấp phụ và khả năng
phản ứng của acetic acid lên bề mặt anatase-TiO 2 (101) bằng phương pháp
kính hiển vi quét điện tử xuyên hầm (STM). Kết quả cho thấy phân tử acetic


4

acid có khả năng tương tác tương đối mạnh với bề mặt anatase-TiO 2 (101).
Năm 2014, Xiang Zhao và các cộng sự [39] đã nghiên cứu sự hấp phụ của
benzoic acid trên bề mặt rutile-TiO2 (110) bằng phương pháp DFT trong đó có
xem xét đến sự ảnh hưởng của liên kết hydro và lực Van der Waals đến độ bền
tương đối của các cấu trúc siêu tầng khác nhau. Năm 2016, Maria Buchholz
và các cộng sự [23] đã nghiên cứu sự hấp phụ của benzoic acid và terephthalic
acid lên bề mặt rutile-TiO2 (110) thấy rằng, sự hình thành các cấu hình hấp
phụ khác nhau dựa vào cấu trúc của các acid ban đầu. Kết quả cho thấy cả hai
phân tử hấp phụ phân li trên TiO2 (110) thông qua tách proton hình thành liên
kết O-H tại nguyên tử O cầu nối (O br) và các nhóm carboxylate liên kết với
các vị trí nguyển tử Ti phối trí 5 (Ti5c).
Gần đây, năm 2017, Martin Setvincùng các cộng sự [25] đã nghiên cứu
lý thuyết và thực nghiệm sự hấp phụ formaldehyde trên bề mặt anatase-TiO2,
(101) thấy rằng quá trình hấp phụ phức tạp và các chuỗi paraformaldehyde là
cấu hình hấp phụ thuận lợi nhất. Cùng năm đó, nhóm của Meijing Yao [40]
nghiên cứu cơ chế hấp phụ các chất hữu cơ dễ bay hơi có chứa carbonyl:
acetaldehyde, acetone và methylacetate trên bề mặt chất xúc tác anatase-TiO2
(001) bằng cách sử dụng các tính tốn DFT. Kết quả chỉ ra rằng các tương tác
chủ yếu thơng qua sự tương tác giữa nhóm carbonyl với các vị trí Ti phối trí
năm của bề mặt TiO2. Khi nguyên tử H của nhóm aldehyde được thay thế bởi

các nhóm khác, kích thước của nhóm thế thay thế ít tác động đến năng lượng
hấp phụ. Năm 2018, Guofei Wu và các cộng sự [41] nghiên cứu về sự tương
tác của phân tử formaldehyde (HCHO) với bề mặt anatase-TiO 2 (101) khơng
và có phân tử nước và oxy được nghiên cứu bằng các tính tốn DFT cho thấy
phân tử HCHO hấp phụ trên bề mặt anatase-TiO2 (101) có năng lượng hấp
phụ ít âm hơn và khơng bị phân tách khi xảy ra hấp phụ, cịn khi có mặt của
phân tử oxi và nước năng lượng hấp phụ âm hơn, tất cả các phân tử đều bị


5

thay đổi liên kết. Sự phá vỡ liên kết, hình thành các tương tác bền trên bề mặt
TiO2 được tìm thấy trong nghiên cứu gần đây của nhóm Trung và cộng sự
[18] về sự hấp phụ axit (formic, acetic, benzoic), phenol và nitrobenzene trên
bề mặt rutile-TiO2 (110) sử dụng các tính tốn DFT. Kết quả cho thấy, các q
trình hấp phụ trên được đặc trưng như là sự hấp phụ hóa học và vai trị của
tương tác O···Ti và O/C-H···O đóng góp vào độ bền của các phức cũng được
phân tích một cách chi tiết.
Nhìn chung, các nghiên cứu gần đây chủ yếu tập trung vào các hướng: i)
nghiên cứu quá trình hấp phụ trên cơ sở nghiên cứu các tương tác hình thành
khi hấp phụ lên mặt vật liệu TiO 2 đơn tinh thể; ii) thay đổi điều kiện, kĩ thuật
để nghiên cứu để khảo sát quá trình hấp phụ; iii) dần chuyển sang hướng
nghiên cứu sự thay đổi liên kết tại nhóm chức của các phân tử hữu cơ và các
liên kết của bề mặt vật liệu TiO2 trong quá trình hấp phụ. Tuy nhiên, khi
nghiên cứu quá trình hấp phụ các chất hữu cơ lên các bề mặt vật liệu TiO 2 đã
cho thấy ngoài tương tác giữa nhóm chức của phân tử hữu cơ với các nguyên
tử trên bề mặt vật liệu TiO 2 thì quá trình hấp phụ xuất hiện một vài tương tác
yếu như: liên kết hydro, tương tác Van der Waals, tương tác tĩnh điện… giữa
các nguyên tử trên vòng benzen của các phân tử bị hấp phụ-chất hấp phụ. Các
tương tác này góp phần làm tăng khả năng hấp phụ các phân tử hữu cơ lên bề

mặt vật liệu. Do đó, việc nghiên cứu lý thuyết bằng phương pháp hóa học
lượng tử với các cơng cụ hiện đại có độ chính xác cao sẽ giúp chúng ta có cái
nhìn sâu sắc hơn về khả năng hấp phụ các chất hữu cơ chứa vòng benzen lên
bề mặt vật liệu TiO2 ở cấp độ phân tử.
3. Mục đích nghiên cứu
- Xác định cấu trúc hình học bền của các phức hình thành khi hấp phụ
các phân tử hợp chất hữu cơ chứa vòng benzen lên trên các bề mặt vật liệu
TiO2 khác nhau.


6

-

So sánh và đánh giá khả năng hấp phụ các phân tử hữu cơ lên các bề

mặt vật liệu của TiO2 khác nhau.
-

Làm rõ vai trò và bản chất của các loại tương tác đóng góp vào độ bền

của các phức do quá trình hấp phụ.
4.

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

4.1. Đối tượng nghiên cứu
-

Các hợp chất hữu cơ chứa vòng benzen như: benzaldehyde

(C6H5CHO), benzoic acid (C6H5COOH), aniline (C6H5NH2), phenol
(C6H5OH), benzenesulfonic acid (C6H5SO3H).
- Các bề mặt vật liệu TiO2 như rutile (110), anatase (101).
4.2. Phạm vi nghiên cứu
Nghiên cứu cấu trúc, độ bền các phức, vai trò của các tương tác và khả
năng hấp phụ một số hợp chất hữu chứa vòng benzen lên các bề mặt vật liệu
TiO2 bằng phương pháp hóa học tính tốn.
5.
-

Phương pháp nghiên cứu
Sử dụng phần mềm Material Studio 2010 để xây dựng các bề mặt vật

liệu TiO2 như rutile (110), anatase (101) [42].
- Sử dụng chương trình mơ phỏng vật liệu VASP với phiếm hàm PBE
để
tối ưu hóa các cấu trúc hình học và tính các thơng số đặc trưng về năng lượng.
-

Sử dụng các phần mềm AIM 2000, NBO 5.G [43], [44] tại mức lý

thuyết B3LYP/6-31+G(d,p) để xác định sự có mặt các tương tác yếu như liên
kết hydro, tương tác acid-base Lewis, các tương tác yếu khác và độ bền tương
quan của chúng trong việc làm bền phức.
-

Sử dụng phần mềm đồ họa hỗ trợ như: Gausview 5.0, VESTA, paint

3D để xác định hình học topo của các monomer, phức[45].

7

6. Ý nghĩa khoa học và thực tiễn của đề tài
Nghiên cứu lý thuyết bằng phương pháp hóa học tính toán, đặc biệt là
nghiên cứu về khả năng hấp phụ các chất hữu cơ chứa vòng benzen lên bề mặt
vật liệu TiO2 là một hướng nghiên cứu mới, lý thú ở Việt Nam và thế giới. Kết
quả đạt được của luận văn sẽ giúp hiểu rõ hơn về bản chất quá trình hấp phụ
trên các bề mặt vật liệu TiO2, giúp các nhà nghiên cứu chọn vật liệu phù hợp
cho việc loại bỏ các chất hữu cơ độc hại có nhóm chức khác nhau. Đáng chú
ý, kết quả đề tài sẽ góp phần giúp các nhà khoa học có sự đánh giá về diễn
biến các quá trình diễn ra trên bề mặt vật liệu trong quá trình hấp phụ, đồng
thời giúp các nhà thực nghiệm giảm chi phí thử nghiệm và thời gian thực hiện
các nghiên cứu.
7. Cấu trúc luận văn
Luận văn gồm 74 trang, trong đó có 10 bảng và 8 hình. Phần mở đầu 7
trang, phần nội dung gồm: Chương 1: 26 trang, Chương 2: 8 trang, Chương 3:
31 trang, Kết luận và kiến nghị 2 trang.


8

CHƯƠNG 1. CƠ SỞ LÍ THUYẾT HĨA HỌC LƯỢNG TỬ
1.1. Phương trình Schrưdinger
Việc xác định trạng thái của hệ lượng tử tại thời điểm t từ trạng thái hiện
tại t0 cần một hàm sóng thể hiện sự thay đổi theo thời gian. Chính vì vậy
người ta đưa ra một hàm gọi là hàm sóng hay hàm trạng thái (Ψ). Đối với hệ
một hạt chuyển động trong không gian một chiều phương trình Schrưdinger
phụ thuộc thời gian có dạng đơn giản nhất: [3] [46]

(x, t) t

−

+ V(x, t)

i

(x, t) (1.1)

Trong đó: ћ =

khối lượng của hạt, i =
của hệ phụ thuộc vào cả biến tọa độ x và biến thời gian t. Hàm sóng Ψ(x,t) là
hàm liên tục, hữu hạn, đơn trị, khả vi, nói chung là phức và chuẩn hóa.
Năm 1926, phương trình Schrưdinger được xây dựng bởi nhà vật lý
người Áo là Erwin Schrödinger trên cơ sở hợp nhất giữa thuyết lượng tử
Planck và thuyết lưỡng tính sóng hạt của Louis De Broglie. Phương trình này
cho biết hàm sóng tại thời điểm t nếu biết hàm sóng tại thời điểm t0. Khi biết
hàm sóng ta sẽ xác định được xác suất tìm thấy hạt trong vùng tọa độ giữa x
và x+dx tại thời điểm t.
Tuy nhiên, trong hóa học hầu hết các hệ lượng tử đều được khảo sát ở
trạng thái dừng – trạng thái mà mật độ xác suất tìm thấy hệ khơng biến đổi
theo thời gian mà chỉ biến đổi theo tọa độ (|Ψ(x,t)|2 = |Ψ(x)|2). Do đó, phương
trình Schrưdinger khơng phụ thuộc thời gian của hệ một hạt, một chiều là:
−

ħ2 d2Ψ(x)

2

+ V(x)Ψ(x) = EΨ(x) (1.2)

2m dx

Trong đó: (x) là hàm sóng chỉ phụ thuộc tọa độ không gian.


Hoặc viết đơn giản dưới dạng:


9

Trong đó, Ψ là hàm riêng của tốn tử Hamilton Ĥ, E là trị riêng năng
lượng của toán tử Ĥ. Phươngtrình (1.2) và (1.3) là phương trình Schrưdinger
trạng thái dừng.
Khi giải phương trình hàm riêng – trị riêng (1.3) sẽ thu được nghiệm
gồm Ψ và E. Hàm riêng Ψ mô tả trạng thái của hệ lượng tử và cho phép rút ra
được tất cả các thông tin về hệ lượng tử. Xét hệ gồm M hạt nhân và N
electron. Toán tử Hamilton Ĥ tổng quát được xác định theo biểu thức:

Trong đó:
e và

hạt nhân

các e
các hạt nhân
U nn : là thế năng tương tác đẩy tĩnh điện

giữa Được biểu diễn đầy đủ theo phương trình sau:

H

ˆ
=

N1

M

2

(1.5)

2

i=1
i

A=1

Trong đó:

A, B: kí hiệu cho hạt nhân A và B
MA: tỉ số khối lượng của hạt nhân A đối với khối lượng của 1 e
i, j : kí hiệu cho electron trong hệ
ZA, ZB: số đơn vị điện tích các hạt nhân A và B tương ứng
rij: khoảng cách giữa electron thứ i và thứ j
riA: khoảng cách giữa electron thứ i và hạt nhân A

RAB: khoảng cách giữa hai hạt nhân A và B


2

là tốn tử Laplace có dạng:


10

Trên thực tế, chỉ có thể giải chính xác phương trình Schrưdinger đối với
hệ 1 electron và 1 hạt nhân (bài toán nguyên tử H và những ion giống H), cịn
các bài tốn khác chỉ có thể giải gần đúng vì ta khơng thể xác định chính xác
thành phần tương tác đẩy giữa hai electron. Sự gần đúng Born−Oppenheimer
coi hạt nhân đứng yên so với electron là sự gần đúng đầu tiên cho phép tách
chuyển động của hạt nhân và electron. Khi coi hạt nhân đứng yên so với
electron, sự chuyển động của electron không phụ thuộc vào vận tốc hạt nhân
mà phụ thuộc vào vị trí của hạt nhân.
Áp dụng cho phương trình (1.4): số hạng thứ hai

ˆ = ; và số hạng Tn 0

cuối cùng ˆ = const(C). Toán tử Hamilton của cả hệ trở thành toán tử U nn
Hamilton của các electron ứng với năng lượng electron toàn phần Eel.

ˆ =−N
1

i

(1.6)

2

Heli=1 2

Khi xét sự chuyển động của hạt nhân trong trường trung bình của các
electron. Tốn tử hạt nhân có dạng:

ˆ
H
nucl

Như vậy, với sự gần đúng Born−Oppenheimer thì hàm sóng đầy đủ cho
hệ N electron, M hạt nhân sẽ là :

(r

i

, RA

)= (r
el

i

, RA

)

nucl

(R )
A

(1.8)

1.2. Nguyên lí khơng phân biệt các hạt đồng nhất, mơ hình các hạt độc
lập
Đối với hệ nhiều electron, phương trình Schrưdinger khơng thể giải
chính xác vì tương tác giữa các electron i và j không thể xác định do không
thể xác định tường minh vị trí của chúng trong khơng gian. Ta chỉ có thể xác
định được mật độ xác suất ở một vị trí xác định của hạt là bao nhiêu. Như


11

vậy, “các hạt trong một hệ đồng nhất là không thể phân biệt được”. Đây là nội
dung của nguyên lý không biệt các hạt đồng nhất.
Trong nguyên tử nhiều electron, mỗi electron chuyển động độc lập với
các electron khác (bỏ qua tương tác đẩy giữa các electron) trong một trường
trung bình có đối xứng cầu (trường xun tâm) tạo bởi hạt nhân nguyên tử và
các electron còn lại.
1.3. Nguyên lý phản đối xứng (nguyên lý loại trừ Pauli) [3] [46]
Trên thực tế, quang phổ phát xạ của nguyên tử hydro bị tách thành hai
phần tương tự nhau dưới tác dụng của từ trường ngoài. Điều này cho thấy bộ
3
số lượng tử n, l, ml được phát hiện trước đây chưa mô tả đầy đủ trạng
thái

electron. Năm 1928, Dirac giải phương trình Schrưdinger theo thuyết tương
đối của Einstein đối với ngun tử hydro và làm xuất hiện số lượng tử spin
(ms = ±1/2).
Mặt khác, tốn tử Hamilton Ĥ trong phương trình (1.3) chỉ phụ thuộc
vào tọa độ không gian của electron, để mơ tả đầy đủ thuộc tính của một
electron ta cần chỉ rõ yếu tố spin, do đó việc đưa spin của electron vào phần
khơng gian của hàm sóng là cần thiết. Gọi α(ω) và β(ω) là hai hàm spin ứng
với spin−up (↑) và spin−down (↓). Hai hàm này trực chuẩn:
∫α*(ω)α(ω)dω = ∫β*(ω)β(ω)dω = 1 hay
∫α*(ω)β(ω)dω = ∫β*(ω)α(ω)dω = 0 hay
Như vậy, một electron lúc này không chỉ được mô tả bởi tọa độ khơng
gian r mà cịn tọa độ spin ω, ký hiệu x={r,ω}. Khi đó, hàm sóng N electron
được viết: Ψ(x1,x2,…, xN). Bởi vì tốn tử Hamilton khơng có yếu tố spin, vì
vậy việc làm cho hàm sóng phụ thuộc vào spin không mang lại bất kỳ cải
thiện nào. Để thỏa mãn điều này cần phải thêm điều kiện đối với một hàm
sóng. Sự thêm này được phát biểu: một hàm sóng nhiều electron phải phản


12

đối xứng với sự trao đổi (hoán đổi) tọa độ x (gồm không gian và spin) của bất
kỳ 2 electron: [3] [46]
Ψ(x1,…,xi,…,xj,…xN) = -Ψ((x1,…,xj,…,xi,…xN) (1.11) Đây chính là
nguyên lý phản đối xứng hay nguyên lý loại trừ Pauli.
Như vậy, hàm sóng mơ tả trạng thái của electron khơng chỉ thỏa mãn phương
trình Schrưdinger mà cịn phản đối xứng.
1.4. Hàm sóng của hệ nhiều electron
Trong sự gần đúng Born−Oppenheimer và mô hình hạt độc lập, Hamilton Ĥ elcủa hệ có thể được xem gần đúng
bằng tổng các Hamilton 1ee và hằng số C [24], [26]. Do đó, hàm sóng Ψ el của hệ có thể là tích Hartree của các hàm sóng
orbital−spin 1ee : [3] [46]

el

Trong đó:

i
i

(x1 , x2 , x3 ,..., xN )= 1 (x1 ) 2 (x2 ) 3 (x3 )... N (xN

)

(1.12)

= i ( r ). ( )
là hàm orbital-spin 1ee thứ i;

xi là tọa độ khái quát của electron i;
i

( ) là hàm spin, có thể là α hoặc β.

Hàm sóng dạng tích Hartree khơng thỏa mãn nguyên lý không phân
biệt các hạt đồng nhất và nguyên lý phản đối xứng. Do đó ta phải tổ hợp tuyến
tính một cách thích hợp các hàm tích này, tổ hợp đó cịn được viết dưới dạng
định thức Slater. Đối với hệ vỏ kín có số chẵn electron, hàm sóng toàn phần
của hệ là một định thức Slater:

el

= (N!)−1 / 2

(x1) 1 (x2)... 1 (xN)
(x1) 2 (x2)... 2 (xN)

1
2

................................
N

(x1) N (x2)...

N

(xN)


13

Với

−1/ 2

(N!)

hàm sóng. Hàng của định thức (1.13) được ký hiệu tương ứng với hàm
orbital−spin, và cột ứng với electron. Khi đổi tọa độ hai electron tương ứng
với sự hoán vị 2 cột sẽ làm thay đổi dấu của định thức. Nếu có hai electron có
cùng spin chiếm cùng một obital thì tương ứng với hai hàng của định thức

bằng nhau và định thức sẽ bằng zero. Vì vậy, chỉ có tối đa 1 electron có thể
chiếm một orbital−spin (nguyên lý loại trừ Pauli). Hàm sóng Ψ el có thể được
biểu diễn ngắn gọn dưới dạng đường chéo chính của định thức Slater, với quy
ước đã có mặt thừa số chuẩn hóa (N!)
(1.14)
Việc đối xứng hóa một tích Hartree để thu được một định thức Slater
mang lại những hiệu ứng trao đổi vì điều kiện hàm sóng

2

bất biến khi đổi chỗ bất

kỳ hai electron. Hàm định thức Slater gắn với tương quan trao đổi, có nghĩa rằng có kể đến sự chuyển động tương quan của hai electron với spin song song. Tuy nhiên, chuyển động
của những electron với spin đối song chưa được xét đến, mặc dù đây là sự đóng góp chính đến năng lượng tương quan. Vì vậy thơng thường ta gọi hàm sóng định thức đơn Slater là
hàm sóng khơng tương quan. Đối với hệ kín (N=2n), hàm sóng một định thức Slater có thể mơ tả tốt trạng thái của hệ, còn đối với hệ mở có số lẻ electron (N=2n+1) hàm sóng phải là
tổ hợp tuyến tính của nhiều định thức Slater.

1.5. Cấu hình electron
Cấu hình electron cho biết sự phân bố electron trong hệ lượng tử vào
các orbital và có thể được phân loại như sau: [3] [46]
-

Cấu hình vỏ đóng (Closed shell): Hệ có 2n electron chiếm n orbital

khơng gian.


14

2n

Cấu hình vỏ mở (Opened shell): Hệ có 2n+1 electron, trong đó

electron chiếm n orbital khơng gian, một electron cịn lại chiếm orbital thứ n+1. Cấu hình hạn chế (Restricted): là cấu hình mà một hàm khơng gian

được sử dụng chung cho hai hàm spin α và hàm spin βkhi có 2 electron ghép
đơi. Các electron độc thân (nếu có) thuộc về các hàm không gian khác nhau.
Như vậy, chỉ các MO bị chiếm bởi các electron không ghép đôi mới được xử
lý riêng biệt. Phương pháp Hartree−Fock với cấu hình này có thể là RHF (với
cấu hình vỏ đóng) hoặc ROHF (với cấu hình vỏ mở hạn chế cho trạng thái
kích thích).
Cấu hình khơng hạn chế (Unrestricted): là cấu hình mà các hàm
spin
và thuộc hai hàm khơng gian khác nhau, nghĩa là không suy biến năng
lượng. Tất cả các electron đều được xử lý ở các orbital không gian riêng biệt.
Cấu hình này có thể sử dụng với hệ có 2n hoặc 2n+1 electron ở các dạng vỏ
đóng, mở và trạng thái kích thích. Phương pháp HF với cấu hình này được gọi
là là UHF, phương pháp này cho kết quả tốt trong trường hợp hệ là gốc, ion.
1.6. Bộ hàm cơ sở
1.6.1. Orbital kiểu Slater và Gaussian
Để giải phương trình Schrưdinger cho phân tử, người ta dùng hàm MO
là tổ hợp tuyến tính các AO (MO−LCAO): [3] [46]
Ψi =

C νi Φi
ν =1n

(1.15)

Cνi là các hệ số tổ hợp; Φi là các AO cơ sở. Tập hợp các hàm Φi được

gọi là bộ cơ sở. Mỗi AO cơ sở gồm phần bán kính và phần góc:
Φ(r,θ,φ) = R(r).Y(θ,φ)
Theo cách biểu diễn tốn học khác nhau của phần bán kính, có 2 kiểu
hàm cơ sở thường được sử dụng:


AO kiểu Slater (Slater Type Orbital, STO):


15

STO

Φ
AO kiểu Gaussian (Gaussian Type Orbital, GTO):
r-R A

ΦGTO = CG.

Trong đó:

r: toạ độ orbital

2

e

RA: toạ độ hạt nhân nguyên tử A
CS, CG: các hệ số (bao gồm phần góc)
η, α: thừa số mũ của các hàm STO và GTO tương ứng, là

các số dương xác định kích thước của orbital[24], [26].
Nhược điểm của hàm Slater là khi dùng để tính những tích phân 2
electron nhiều tâm (nhiều ngun tử) thì rất khó hội tụ. Những phép tính này
được đơn giản rất nhiều khi dùng hàm Gaussian. Tuy nhiên, hàm Gaussian mô
tả không tốt trạng thái ở gần nhân và cả ở những khoảng cách lớn xa nhân
vì dΦ
dr

= nhưng cực trị của hàm sóng khơng phải đạt được tại hạt
0

nhân, và khi r → ∞ thì ΦGTO giảm quá nhanh. Để có bộ hàm cơ sở tốt hơn, có
thể làm theo 2 cách:
Tổ hợp tuyến tính n hàm GTO thành 1 hàm STO thu được các
bộ hàm
cơ sở STO−nG.
Ví dụ: bộ hàm STO−3G, STO−4G, …
- Tổ hợp tuyến tính một số hàm GTO gọi là các hàm Gaussian ban đầu
(PGTO, Primitive Gaussian−type Orbital) thu được hàm Gaussian rút gọn, kí
hiệu là CGF (Contracted Gaussian Functions):
Φ

CGF

k

a i Ψ iGTO (1.16)

=
i

Trong đó: ai: hệ số rút gọn (contraction coefficient) được chọn sao cho


ΦCGF giống hàm STO nhất; k: bậc rút gọn.


16

1.6.2. Một số khái niệm về bộ hàm cơ sở
Một bộ hàm cơ sở tối thiểu (minimal basis sets) bao gồm tất cả các
orbital vỏ trong và orbital vỏ hóa trị.
Ví dụ:

Nguyên tử H: 1s;
Nguyên tử C: 1s, 2s, 2px, 2py, 2pz

Bộ hàm cơ sở hóa trị tách (split valence basis set): tách đôi, tách ba, …
số hàm cơ sở cho mỗi AO hóa trị, làm thay đổi kích thước của orbital.
Ví dụ:

Nguyên tử H: 1s, 1s’;
Nguyên tử C: 1s, 2s, 2s’, 2px, 2py, 2pz , 2px’, 2py’, 2pz’

Bộ hàm cơ sở mở rộng (extended basis sets): là bộ cơ sở tối thiểu bao
gồm thêm các orbital của lớp vỏ bên ngồi vỏ hóa trị (orbital ảo – virtual
orbital)
Hàm phân cực: thêm vào các AO có momen góc lớn hơn cho nguyên tử
nặng hoặc nguyên tử H vào bộ cơ sở hóa trị tách, có thể làm biến đổi hình
dạng các orbital.

Hàm khuyếch tán: là những hàm s, p có kích thước lớn, mơ tả các
orbital có khơng gian lớn hơn. Bộ cơ sở có hàm khuyếch tán quan trọng với
các hệ có electron ở xa hạt nhân như các phân tử có cặp electron riêng, các
anion, các hệ ở trạng thái kích thích, hệ có thế ion hóa thấp, hệ tương tác
yếu…
1.6.3. Phân loại bộ hàm cơ sở
*

Bộ cơ sở kiểu Pople:
Bộ cơ sở STO−nG: tổ hợp n GTO thành STO, với n = 2 6.
Thực tế n

> 3, kết quả rất ít thay đổi so với n = 3, do đó bộ hàm STO−3G được sử dụng
rộng rãi nhất và cũng là bộ cơ sở cực tiểu.
Bộ cơ sở k−nlmG: với k là số hàm GTO dùng làm orbital lõi,
bộ số