Silicát là nhóm vật liệu được ứng dụng rộng rãi trong nhiều lĩnh
vực công nghiệp (điện tử, quang học, siêu dẫn, cơ khí.) cũng như
đời sống (gốm, men, thủy tinh.). Ví dụ như nhôm-silicát được ứng
dụng trong công nghiệp hóa học, công nghiệp cao su, sản xuất ra các
vật liệu chuyên dụng như gạch men, da giày nhân tạo v.v., đặc biệt
trong một số ứng dụng công nghệ cao như chì-silicát dùng che chắn
phóng xạ.
Hiện tại cấu trúc và động học không đồng nhất của các hệ vật
liệu silicát vẫn đang là vấn đề mang tính chất thời sự, được nhiều nhà
khoa học quan tâm nghiên cứu. Các nghiên cứu trước đó đã chỉ ra
rằng các hệ vật liệu silicát có cấu trúc mạng ngẫu nhiên liên tục, bao
gồm các đơn vị cấu trúc SiOx. Trong đó sự chuyển pha mật độ được
cho có liên quan đến sự thay đổi tỷ phần các đơn vị cấu trúc. Ngoài
ra, cấu trúc không đồng nhất được tạo ra là do sự phân bố khác nhau
của các đơn vị cấu trúc trong mô hình, từ đó hình thành nên vùng
giàu cation. Giữa các đơn vị cấu trúc được kết nối với nhau thông
qua nguyên tử O, mức độ polymer hóa (DOP) cũng được đánh giá
qua các nguyên tử Oxy cầu Qn. Tuy nhiên, thông tin chi tiết về sự
hình thành các cụm với kích thước bao nhiêu, các đơn vị cấu trúc
phân bố đồng đều hay tách ra thành các cụm riêng biệt vẫn chưa
được làm rõ. Đặc biệt hơn trong hệ silicát ba nguyên, với sự có mặt
của hai loại cation thì sự kết nối giữa các đơn vị cấu trúc của hai loại
như thế nào? có thực sự tồn tại các vùng mà ở đó các cation tách
nhau ra không? Trong phạm vi luận án này chúng tôi sẽ làm rõ hơn,
chi tiết hơn về vấn đề này. Ngoài ra, để làm sáng tỏ hiện tượng
không đồng nhất cấu trúc, chúng tôi đưa ra hai phương pháp mới để
phân tích là simplex và shell-core (SC)
27 trang |
Chia sẻ: tranhieu.10 | Lượt xem: 1515 | Lượt tải: 0
Bạn đang xem trước 20 trang tài liệu Luận án Nghiên cứu cấu trúc và sự không đồng nhất động học trong vật liệu silicát ba nguyên PbO.SiO2, Al2O3.2SiO2 và Na2o.2SiO2 ở trạng thái lỏng và vô định hình, để xem tài liệu hoàn chỉnh bạn click vào nút DOWNLOAD ở trên
BỘ GIÁO DỤC VÀ ĐÀO TẠO
TRƯỜNG ĐẠI HỌC BÁCH KHOA HÀ NỘI
NGUYỄN VĂN YÊN
NGHIÊN CỨU CẤU TRÚC VÀ SỰ KHÔNG ĐỒNG NHẤT
ĐỘNG HỌC TRONG VẬT LIỆU SILICÁT BA NGUYÊN
PbO.SiO2, Al2O3.2SiO2 và Na2O.2SiO2 Ở TRẠNG THÁI
LỎNG VÀ VÔ ĐỊNH HÌNH
Chuyên ngành : VẬT LÝ KỸ THUẬT
Mã số: 62520401
TÓM TẮT LUẬN ÁN TIẾN SĨ VẬT LÝ KỸ THUẬT
HÀ NỘI – 2017
Công trình được hoàn thành tại:
Trường Đại học Bách khoa Hà Nội
Người hướng dẫn khoa học: 1. PGS.TS. NGUYỄN VĂN HỒNG
2. PGS.TS. LÊ THẾ VINH
Phản biện 1:
Phản biện 2:
Phản biện 3:
Luận án được bảo vệ trước Hội đồng đánh giá luận án tiến sĩ cấp
Trường họp tại Trường Đại học Bách khoa Hà Nội
Vào hồi .. giờ, ngày .. tháng .. năm
Có thể tìm hiểu luận án tại thư viện:
1. Thư viện Tạ Quang Bửu - Trường ĐHBK Hà Nội
2. Thư viện Quốc gia Việt Nam
1
MỞ ĐẦU
1. Lý do chọn đề tài
Silicát là nhóm vật liệu được ứng dụng rộng rãi trong nhiều lĩnh
vực công nghiệp (điện tử, quang học, siêu dẫn, cơ khí...) cũng như
đời sống (gốm, men, thủy tinh...). Ví dụ như nhôm-silicát được ứng
dụng trong công nghiệp hóa học, công nghiệp cao su, sản xuất ra các
vật liệu chuyên dụng như gạch men, da giày nhân tạo v.v..., đặc biệt
trong một số ứng dụng công nghệ cao như chì-silicát dùng che chắn
phóng xạ.
Hiện tại cấu trúc và động học không đồng nhất của các hệ vật
liệu silicát vẫn đang là vấn đề mang tính chất thời sự, được nhiều nhà
khoa học quan tâm nghiên cứu. Các nghiên cứu trước đó đã chỉ ra
rằng các hệ vật liệu silicát có cấu trúc mạng ngẫu nhiên liên tục, bao
gồm các đơn vị cấu trúc SiOx. Trong đó sự chuyển pha mật độ được
cho có liên quan đến sự thay đổi tỷ phần các đơn vị cấu trúc. Ngoài
ra, cấu trúc không đồng nhất được tạo ra là do sự phân bố khác nhau
của các đơn vị cấu trúc trong mô hình, từ đó hình thành nên vùng
giàu cation. Giữa các đơn vị cấu trúc được kết nối với nhau thông
qua nguyên tử O, mức độ polymer hóa (DOP) cũng được đánh giá
qua các nguyên tử Oxy cầu Qn. Tuy nhiên, thông tin chi tiết về sự
hình thành các cụm với kích thước bao nhiêu, các đơn vị cấu trúc
phân bố đồng đều hay tách ra thành các cụm riêng biệt vẫn chưa
được làm rõ. Đặc biệt hơn trong hệ silicát ba nguyên, với sự có mặt
của hai loại cation thì sự kết nối giữa các đơn vị cấu trúc của hai loại
như thế nào? có thực sự tồn tại các vùng mà ở đó các cation tách
nhau ra không? Trong phạm vi luận án này chúng tôi sẽ làm rõ hơn,
chi tiết hơn về vấn đề này. Ngoài ra, để làm sáng tỏ hiện tượng
không đồng nhất cấu trúc, chúng tôi đưa ra hai phương pháp mới để
phân tích là simplex và shell-core (SC).
Vấn đề động học không đồng nhất (DH) trong các hệ chất lỏng
đã được ghi nhận, trong đó quan tâm đến vùng chuyển động nhanh
và các vùng chuyển động chậm. DH được nghiên cứu khá chi tiết
2
trong các hệ keo qua việc phân tích các hạt chuyển động nhanh và
chuyển động chậm. Theo tác giả Antonio M. Puertas, những hạt
chuyển động nhanh thì lân cận của chúng cũng chuyển động nhanh,
hạt chuyển động chậm thì lân cận của chúng cũng chuyển động
chậm. Điều này tạo nên vùng nhanh chậm tách nhau ra, hơn nữa mật
độ vùng nhanh và chậm là rất khác nhau. Ngoài ra, phân bố không
gian của các vùng nhanh và chậm cũng được Claudio Donati quan
tâm nghiên cứu, trong đó tác giả cho rằng phân bố không gian của
những hạt chuyển động nhanh thì lớn hơn những hạt chuyển động
chậm. Điều này dẫn đến vùng mật độ cao là vùng chứa các hạt
chuyển động chậm và mật độ thấp chứa các hạt chuyển động nhanh.
Ngoài các hệ keo thì các kết quả tương tự cũng được quan sát thấy
trên các chất lỏng nguội nhanh, tuy nhiên với các hệ ôxít lỏng đặc
biệt là các hệ silicát thì số lượng các công trình nghiên cứu về DH là
rất ít. Cho đến nay chỉ mới ghi nhận tác giả K.D. Vargheese và các
cộng sự của ông nghiên cứu về DH trên hệ nhôm-silicát năm 2010.
Tuy nhiên tác giả chỉ mới dừng lại ở mức độ định tính, theo tác giả
trong mô hình có tồn tại các vùng chuyển động nhanh và chuyển
động chậm, và các vùng chuyển động chậm là vùng giàu Si và O, các
vùng chuyển động nhanh là vùng giàu Al và Ca. Ngoài ra mật độ
vùng chuyển động nhanh như thế nào? liệu rằng các hạt chuyển động
nhanh có sự kết cụm hay không? Các liên kết của chúng giảm như
thế nào theo thời gian? Tất cả những vấn đề này vẫn chưa được làm
rõ, hơn nữa có tồn tại hay không sự tương quan giữa động học và cấu
trúc trong hệ vật liệu silicát và cụ thể ra sao? Đây đều là những câu
hỏi còn bỏ ngỏ, trên cơ sở đó chúng tôi chọn đề tài “Nghiên cứu cấu
trúc và sự không đồng nhất động học trong vật liệu Silicát ba
nguyên PbO.SiO2, Al2O3.2SiO2 và Na2O.2SiO2 ở trạng thái lỏng và
vô định hình” nhằm đưa đến cái nhìn rõ ràng hơn về cấu trúc cũng
như động học của hệ silicát ba nguyên.
2. Mục đích đối tượng và phạm vi nghiên cứu
3
Đối tượng nghiên cứu của luận án là các hệ silicát ba nguyên
gồm PbO.SiO2 (lỏng), Al2O3.2SiO2 (lỏng và vô định hình) và
Na2O.2SiO2 (lỏng). Trong đó phạm vi nghiên cứu của đề tài chủ yếu
là cấu trúc và động học với nội dung như sau:
+ Nghiên cứu vi cấu trúc của các hệ vật liệu PbO.SiO2 (lỏng),
Al2O3.2SiO2 (lỏng và vô định hình) dưới ảnh hưởng của áp suất.
Trong đó đánh giá vi cấu trúc thông qua các thông số cơ bản như
hàm phân bố xuyên tâm, số phối trí và phân bố góc.
+ Nghiên cứu cấu trúc không đồng nhất của hệ vật liệu Al2O3.2SiO2
và Na2O.2SiO2 lỏng qua việc phân tích các quả cầu simplex và các
hạt shell-core (SC).
+ Nghiên cứu về DH của các hệ vật liệu Al2O3.2SiO2 và Na2O.2SiO2
lỏng qua việc phân tích chuyển động nhanh và chậm của một nhóm
các nguyên tử.
+ Nghiên cứu mối tương quan giữa động học và cấu trúc của hệ vật
liệu Al2O3.2SiO2 và Na2O.2SiO2 lỏng.
3. Phương pháp nghiên cứu
+ Phương pháp mô phỏng động lực học phân tử
+ Phương pháp phân tích cấu trúc địa phương
+ Phương pháp simplex và SC phân tích cấu trúc không đồng nhất.
+ Phương pháp phân tích cấu trúc bằng kỹ thuật trực quan hóa.
+ Phương pháp phân tích DH qua việc đánh giá chuyển động nhanh
và chuyển động chậm của một nhóm các nguyên tử.
4. Ý nghĩa khoa học và thực tiễn của đề tài
Các kết quả của Luận án cho một cái nhìn tổng quan về mặt cấu
trúc của các hệ silicát ba nguyên PbO.SiO2 (lỏng), Al2O3.2SiO2 (lỏng
và vô định hình) và Na2O.2SiO2 (lỏng), đây sẽ là những thông tin
hữu ích cho những nghiên cứu tiếp theo. Bên cạnh đó, những nghiên
cứu của chúng tôi cũng cho thấy sự ảnh hưởng của áp suất hay số hạt
lên vi cấu trúc của các hệ silicát ba nguyên.
4
Luận án cũng cho thấy một số cấu trúc đặc biệt tồn tại trong các
mô hình vật liệu silicát ba nguyên như, giữa hai đơn vị liền kề tồn tại
các liên kết góc, liên kết cạnh và liên kết mặt. Ngoài ra trong mô
hình nhôm-silicát lỏng tồn tại các subnet Si-O, Al-O. Trong đó,
subnet Si-O có kích thước rất lớn, subnet Al-O có kích thước nhỏ
hơn nhưng số lượng thì lớn hơn.
Luận án chỉ ra trong các mô hình Na2O.2SiO2, Al2O3.2SiO2 lỏng
tồn tại DH, với vùng chuyển động nhanh và vùng chuyển động chậm
của các nguyên tử O và các cation, trong đó mật độ vùng chuyển
động nhanh và chuyển động chậm là khác nhau.
5. Những đóng góp mới của luận án
Luận án lần đầu tiên chỉ ra được cấu trúc không đồng nhất của
hệ Al2O3.2SiO2 và Na2O.2SiO2 lỏng bằng phương pháp phân tích
mới là simplex và Shell-core (SC). Phương pháp simplex cho phép
tìm ra các quả cầu đi qua bốn nguyên tử bất kỳ mà trong đó xảy ra
các trường hợp như không có nguyên tử nào bên trong, hoặc chỉ chứa
nguyên tử O, hoặc chỉ chứa cation. Với phương pháp phân tích cấu
trúc bằng SC thì cho phép tìm ra các quả cầu với phần lõi chứa
cation phần vỏ chứa O. Các kết quả phân tích cho thấy các nguyên tử
phân bố không đồng đều về không gian cũng như thành phần hóa
học, trong đó hình thành các cụm giàu cation hay O. Ngoài ra, với hệ
nhôm-silicát còn tồn tại các subnet Si-O lớn và các subnet Al-O nhỏ
hơn.
Luận án cũng lần đầu tiên xác định DH cho hệ nhôm và natri-
silicát thông qua phân tích chuyển động của từng loại nguyên tử, và
đánh giá DH thông qua sự kết cụm của các nguyên tử chuyển động
nhanh và chuyển động chậm. Trong đó đối với mô hình nhôm-silicát
lỏng thì vùng chuyển động nhanh là vùng có mật độ cao các nguyên
tử O, vùng chuyển động chậm là vùng có mật độ thấp các nguyên tử
O, với hệ natri-silicát thì ngược lại.
Sự tương quan giữa cấu trúc và động học trong hệ nhôm và natri-
silicát lỏng cũng được nghiên cứu, trong mô hình tồn tại các hạt (siêu
5
phân tử SM và SLP) có kích thước lớn gồm nhiều nguyên tử chuyển
động cùng nhau theo thời gian. Trong đó thời gian sống các hạt này
có liên quan đến các liên kết bền Si-O và các liên kết kém bền Al-O.
6. Cấu trúc của luận án
Ngoài phần mở đầu và kết luận, thì nội dung của luận án được
chia là 5 chương cụ thể như sau:
Chương 1 tổng quan: nội dung chính của chương này là trình bày
những hiểu biết chung về cấu trúc cũng như động học của ba hệ
silicát là nhôm-silicát, chì-silicát và natri-silicát. Chương 2 phương
pháp nghiên cứu: nội dung của chương này trình bày về cách xây
dựng các mô hình PbO.SiO2 lỏng, Al2O3.2SiO2 lỏng và vô định hình
và Na2O.2SiO2 lỏng. Phương pháp phân tích cấu trúc simplex và SC,
phương pháp phân tích DH. Chương 3 vi cấu trúc của hệ chì-silicát
và nhôm-silicát trong đó phân tích cấu trúc của hai hệ thông qua các
thông số cơ bản như hàm phân bố xuyên tâm, phân bố góc và phân
bố số phối trí. Chương 4 phân tích cấu trúc không đồng nhất bằng
phương pháp simplex và shell-core, trong đó chỉ ra các cấu trúc
không đồng nhất của hệ nhôm và natri-silicát lỏng bằng cách phân
tích simplex cũng như shell-core. Chương 5 động học không đồng
nhất của hệ natri và nhôm silicát trong đó chỉ ra sự không đồng nhất
bằng các chuyển động nhanh nhất và chậm nhất xét cho từng loại
loại nguyên tử.
Chương 1 TỔNG QUAN
Vật liệu silica và silicát có vai trò quan trọng trong các lĩnh vực
kỹ thuật cũng như đời sống. Vì vậy để tối ưu hóa các quá trình công
nghệ thì các đặc trưng của chúng vẫn không ngừng được nghiên cứu
bằng cả thực nghiệm và mô phỏng. Trong đó đặc trưng về cấu trúc
và động học đang là vấn đề mang tính thời sự, nhận được nhiều sự
quan tâm từ các nhà khoa học trong và ngoài nước. Vật liệu silicát
chính là sự pha trộn của các ôxít với với silica, các ôxít khác nhau
khi pha trộn cho những tính chất và đặc tính rất khác nhau. Ngoài ra
6
tỷ lệ pha trộn cũng ảnh hưởng không nhỏ đến đặc tính của vật liệu,
hơn nữa điện tích các nguyên tố là khác nhau điều này cũng tạo ra sự
cân bằng điện tích khác nhau dẫn đến cấu trúc mật độ địa phương và
động học cũng rất khác nhau.
Cấu trúc và động học của hệ silicát ba thành phần được nghiên
cứu trong thời gian rất dài, đã cho những kết quả rất quan trọng, đặt
nền móng cho những nghiên cứu tiếp theo. Một trong những kết quả
đó, là sự tồn tại cấu trúc mạng trong cấu trúc mạng silicát, trong đó
tồn tại các đơn vị cấu trúc kết nối với nhau thông qua các nguyên tử
O cầu. Với hệ nhôm-silicát bằng phương pháp thực nghiệm cộng
hưởng từ hạt nhân tác giả S.H. Risbud lần đầu tiên phát hiện ra sự
tồn tại các đơn vị cấu trúc AlO4, AlO5 và AlO6 năm 1987. Những kết
quả này sau đó được S.Sen và R.K. Sato khẳng định lại, ngoài ra với
tỷ lệ ôxít nhôm tăng lên thì số phối trí của O quanh Al cũng có xu thế
tăng lên và vuợt trội so với Si (J. Am. Ceram. Soc. 70, C10 (1987), J.
Chem. Phys. 95 4483(1991), J. Phys.Chem. B 108, 7557 (2004)). Sự
tách pha vi mô được tìm có liên quan đến vùng giàu Al và giàu Si,
mạng cấu trúc Al-O được hình thành các cụm lớn và một phần lồng
vào mạng của Si-O. Với hệ chì-silicát, một trong những nghiên cứu
sớm phải kể đến Bair 1936 bằng phương pháp thực nghiệm X-ray,
cấu trúc của nó được cho là tương tự như sôđa silicát. Nhiều nghiên
cứu sau đó cho thấy chì-silicát có cấu trúc tồn tại các mạng nguyên
tử, trong đó ôxít PbO có cấu trúc tạo thành các chuỗi polymer (J.
Am. Ceram. Soc 19(1–12), 339–347(1936), Condens. Matter, 13(43),
9781–9797(2001), J. Am. Ceram. Soc, 88(6), 1591–1596 (2005),
Physical Review B, 82(13), 134209 (2010)). Với hệ natri-silicát, Na
đóng vai trò là nguyên tố thay đổi cấu trúc mạng silica, nghiên cứu
đã chỉ ra rằng cấu trúc của natri-silicát không tồn tại các đơn phân tử
Na2O, SiO2, Na2Si2O5 và Na2SiO3 một cách rời rạc tách biệt. Kỹ
thuật tán xạ Nơtron cũng đã được sử dụng để thăm dò cấu trúc của
các hệ silicát của kiềm. Tác giả Misawa và các cộng sự (J. Non-
Cryst. Solids 37, 85-97(198)) đã khảo sát cấu trúc trật tự gần của hệ
7
kiềm silicát là natri-silicát và liti-silicát và thấy rằng, nó gồm các
mạng tứ diện của Si. Song song với việc nghiên cứu cấu trúc của hệ
silicát là vấn đề động học, và giữa cấu trúc và động học có một sự
tương quan, tức là một thay đổi của cấu trúc cũng ảnh hưởng lên quá
trình động học. Các quá trình động học cũng được cho là liên quan
chặt chẽ đến cấu trúc của nguyên tử O cầu, với sự khảo sát về mức
độ polymer hóa (DOP). Nó được đo bằng tỷ số (Q3 + Q4)/(Q1 + Q2),
khi tỷ số này tăng lên cho thấy DOP tăng lên, điều này dẫn đến làm
tăng hệ số khuếch tán và giảm độ nhớt mô hình ( ISIJ International,
Vol. 52 No. 3, pp. 342(2012), Materials Sciences and Applications,
5, 73 (2014)). Ngoài sự ảnh hưởng của nguyên tử O cầu, thì nguyên
tử O không cầu và O tricluter cũng là yếu tố gây nên những ảnh
hưởng đến quá trình động học. Ở đây cơ chế chủ yếu của khuếch tán
là do sự đứt gãy và hình thành liên tục theo chuỗi các liên kết các
NBO và O tricluter (J. Non-Cryst. Solids 357, 1780 (2011)).
Chương 2 PHƯƠNG PHÁP TÍNH TOÁN
Trong luận án chúng tôi sử dụng phương pháp mô phỏng động
lực học phân tử với thế tương tác cặp Born – Meyer và thế tương tác
ba thành phần (J.Phys.Chem.C, 113(49) 20773-20784(2009), Z.
Phys. A: Hadrons Nucl, 75(1) 1-18(1932)). Với mô hình chì và
nhôm-silicát chúng tôi sử dụng thế tương tác cặp Born – Meyer có
phương trình:
U(rij) = qiqj
ị
+ Aijexp(-Bijrij) (2.2)
Trong đó các giá trị qi và qj tương ứng là các điện tích của nguyên tử
thứ i và thứ j, Aij và Bij tương ứng là các thông số thế.
Với hệ Na2O.2SiO2 được xây dựng với thế tương tác hai và ba
thành phần, với thế hai thành phần có phương trình như sau:
8
Uij = qiqj
ị
+ fo(bi + bj)exp(
ị ) +
ị
+ D1ijexp(-rij) + D2ijexp(-
rij) (2.3)
Trong đó Uij là thế năng tương tác, qi và qj tương ứng là điện tích
của nguyên tử i và j, rij là khoảng cách giữa nguyên tử thứ i và
nguyên tử thứ j, fo là hằng số có giá trị là 41.865 kJ nm-1mol-1. Số
hạng thứ nhất là lực tương tác tĩnh điện, số hạng thứ hai là các tương
tác gần, thành phần thứ ba mô tả tương tác Van der Waals, thành
phần thứ tư và thứ năm liên quan đến ảnh hưởng của bán kính các
liên kết cộng hóa rị. Ngoài ra hệ Na2O.2SiO2 còn được xây dựng với
thế tương tác ba thành phần theo như phương trình:
Uijk= -f[cos{2(ị-)}-1]ị (2.4)
Kij =
[(ị
)]
(2.5)
Với f là hằng số lực; ị là góc giữa ba nguyên tử i-j-k; , gr, rm là
các thông số điều chỉnh các liên kết hóa trị.
Từ các thế tương tác cặp và thế tương tác ba thành phần như ở
trên, qua các bước chạy thống kê hồi phục với NPT và NVE, chúng
tôi xây dựng được bốn mẫu mô hình, gồm chì-silicát lỏng, nhôm-
silicát lỏng và vô định hình và Natri-silicát lỏng. Riêng với hệ natri-
silicát thì được xây dựng tại viện Riken nhật bản ở nhiệt độ 1873K
với 7992 nguyên tử trong đó có 1776 nguyên tử Na và Si và 4440
nguyên tử O, mẫu được xây dựng với hai mô hình ở áp suất 0 GPa và
8 GPa.
Chương 3 VI CẤU TRÚC CỦA HỆ NHÔM-SILICÁT VÀ CHÌ-
SILICÁT
Trong chương này chúng tôi sẽ nghiên cứu cấu trúc vi mô của
các hệ ba nguyên silicát của nhôm và chì. Cấu trúc vi mô của các hệ
sẽ được nghiên cứu chi tiết thông qua các đơn vị cấu trúc, phân bố
góc và hàm phân bố xuyên tâm. Trật tự gần được khảo sát trong một
9
đơn vị cấu trúc TOx (trong đó T là các cation, x là số nguyên tử Oxy
lân cận nguyên tử T), trật tự khoảng trung được khảo sát giữa hai
đơn vị cấu trúc liên kề. Qua đó vi cấu trúc của hệ được khảo sát một
cách chi tiết như sự phụ thuộc của vi cấu trúc vào sự thay đổi áp suất
hoặc thay đổi thành phần hóa học.
Bảng 3.1. So sánh thực nghiệm của các đặc trưng cấu trúc của hệ
PbO.SiO2 lỏng.
rLK và θT-O-T Luận án Thực nghiệm [68]
rSi-O(Å) 1.64 1.60
rPb-O(Å) 2.34 2.28
rO-O(Å) 2.64 2.58
rSi-Si(Å) 3.24 2.88
rPb-Si(Å) 3.60 3.10
rPb-Pb(Å) 3.64 3.20
θO-Si-O(độ) 1050 109.50
θO-Pb-O(độ) 600 620
Để khẳng định các mô phỏng của các hệ ba nguyên silicát có độ
tin cậy cao, chúng tôi đã so sánh các đặc trưng cấu trúc của hệ với
các kết quả thực nghiệm cũng như các kết quả mô phỏng khác. Thật
vậy các thông số về đặc trưng cấu trúc cơ bản của các hệ chì- silicát
lỏng và nhôm-silicát vô định hình và lỏng được tổng hợp như trong
bảng 3.1, bảng 3.2 và bảng 3.3. Trong đó khoảng cách liên kết rLK là
Bảng 3.2. So sánh thực nghiệm và mô phỏng các thông số cấu trúc của
hệ Al2O3.2SiO2 vô định hình.
rLK và θT-O-T Luận án Mô phỏng [6, 126] Thực nghiệm
[132]
rSi-O(Å) 1.60 1.65 1.59
rAl-O(Å) 1.66 1.67 1.70
rO-O(Å) 2.60 2.68 -
rSi-Si(Å) 3.18 3.15 -
rAl-Si(Å) 3.18 3.23 -
rAl-Al(Å) 3.16 3.15 -
θO-Si-O(độ) 980 1010 -
θO-Al-O(độ) 950 92.50 -
10
vị trí đỉnh đầu tiên của hàm phân bố xuyên tâm, θT-O-T là các góc liên
kết. Các kết quả trên cho thấy phù hợp với thực nghiệm và mô
phỏng, tức là các hệ chì-silicát và nhôm-silicát mà chúng tôi xây
dựng có độ tin cậy cao. Từ những kết quả này chúng tôi tiến hành
các nghiên cứu khác về cấu trúc cũng như động học không đồng
nhất, phần này sẽ được trình bày chi tiết hơn trong các chương sau.
Bảng 3.3. So sánh thực nghiệm và mô phỏng các thông số cấu trúc của hệ
Al2O3.2SiO2 lỏng.
rLK và θT-O-T Luận án Mô phỏng [6, 126] Thực nghiệm [51]
rSi-O(Å) 1.58 1.63 1.61
rAl-O(Å) 1.64 1.64 1.69
rO-O(Å) 2.64 2.68 -
rSi-Si(Å) 3.16 3.10 -
rAl-Si(Å) 3.16 3.16 -
rAl-Al(Å) 3.16 3.10 -
θO-Si-O(độ) 920 940 -
θO-Al-O(độ) 900 870 -
Chương 4 PHÂN TÍCH CẤU TRÚC BẰNG PHƯƠNG PHÁP
SIMPLEX VÀ SHELL-CORE
4.1. Phương pháp simplex
4.1.1. Mô hình Na2O.2SiO2 lỏng
Voi- simplex (VS) là quả cầu đi qua bốn nguyên tử bất kỳ mà ở
bên trong là trống rỗng ký hiệu các loại VS là abc, với a, b và c
tương ứng là nguyên tử O, Si và Na. Cation-simplex(CS) và oxy-
simplex(OS) là các vòng tròn đi qua bốn nguyên tử chỉ chứa cation
hoặc Oxy bên trong, ký hiệu a1b1c1 trong đó a1, b1 và c1 tương ứng là
số lượng các nguyên tử O, Si và Na. Trong bảng 4.1 là đặc trưng các
VS của hệ Na2O.2SiO2 lỏng tại 1873K với áp suất 0.1 MPa và 8
GPa, trong đó mVS là tỷ lệ VS trung bình tính trên mỗi nguyên tử,
RVS là bán kính VS tương ứng.
Từ kết quả cho trong bảng 4.1 có thể thấy tỷ lệ VS cho loại 211
là lớn nhất, nhưng chúng có bán kính nhỏ nhất, tỷ lệ VS giảm theo
thứ tự từ 211 đến 013. Khi nén thì bán kính các simplex đều giảm,
11
đồng thời các simplex lớn bị biến mất, phân bố bán kính simplex có
giá trị trong khoảng 1.6 Å đến 2.8 Å. Trong bảng 4.2 và bảng 4.3 là
số OS và CS, trong đó số nguyên tử oxy có trong các OS là từ 3 đến
9 và tỷ lệ loại 300 là lớn nhất. Khi số nguyên tử oxy bên trong tăng
lên, thì bán kính các OS tăng nhưng số lượng các OS giảm dần. Các
CS có số lượng từ 1 đến 5 nguyên tử, khi nén các CS có kích thước
lớn hơn 3 Å bị mất đi. Ngoài ra, những CS có kích thước lớn thì số
lượng các nguyên tử cation bên trong cũng tăng lên, điều này đúng
cho cả CS ở áp suất cao và áp suất thấp. Có thể thấy trong mô hình
tồn tại các vùng có điện tích âm và vùng điện tích dương, và tạo nên
sự không đồng nhất về thành phần hóa học. Sự dịch chuyển của các
nguyên tử O bên trong các OS có kích thước lớn thì chậm hơn các
nguyên tử còn lại, điều này cho thấy các vùng có điện tích âm là
vùng chuyển độn