Mô phỏng đặc trưng cấu trúc và cơ chế khuếch tán trong vật liệu niken lỏng và rắn

ĐẠI HỌC THÁI NGUYÊN
TRƯỜNG ĐẠI HỌC SƯ PHẠM

VŨ THỊ THANH HƯƠNG

MÔ PHỎNG ĐẶC TRƯNG CẤU TRÚC
VÀ CƠ CHẾ KHUẾCH TÁN TRONG VẬT LIỆU
NIKEN LỎNG VÀ RẮN

LUẬN VĂN THẠC SĨ VẬT LÝ

THÁI NGUYÊN - 2019


ĐẠI HỌC THÁI NGUYÊN
TRƯỜNG ĐẠI HỌC SƯ PHẠM

VŨ THỊ THANH HƯƠNG

MÔ PHỎNG ĐẶC TRƯNG CẤU TRÚC
VÀ CƠ CHẾ KHUẾCH TÁN TRONG VẬT LIỆU
NIKEN LỎNG VÀ RẮN
Ngành: Vật lý chất rắn
Mã ngành: 8440104

LUẬN VĂN THẠC SĨ VẬT LÝ

Người hướng dẫn khoa học: TS. Phạm Mai An

THÁI NGUYÊN - 2019

LỜI CAM ĐOAN
Tôi xin cam đoan đây là đề tài riêng của tôi, do chính tôi thực hiện dưới sự
hướng dẫn của TS. Phạm Mai An và trên cơ sở nghiên cứu các tài liệu tham
khảo. Đề tài này không trùng với kết quả của tác giả khác đã công bố. Nếu sai
tôi hoàn toàn chịu trách nhiệm trước hội đồng.
Thái Nguyên, tháng 6 năm 2019
Học viên

Vũ Thị Thanh Hương


LỜI CẢM ƠN
Luận văn được hoàn thành dưới sự hướng dẫn và chỉ bảo tận tình của
TS. Phạm Mai An. Tôi xin bày tỏ lòng biết ơn sâu sắc đến Thầy. Thầy đã tận
tình hướng dẫn, giúp đỡ tôi trong suốt quá trình học tập, giảng giải cho tôi các
vấn đề liên quan và nghiên cứu để hoàn thành luận văn.
Tôi xin gửi lời cảm ơn tới Ban giám hiệu nhà trường, Ban chủ nhiệm khoa
và thầy cô giáo trong khoa Vật lý, Trường Đại học Sư phạm - Đại học Thái
Nguyên, đã tận tình giảng dạy và tạo điều kiện thuận lợi cho tôi trong suốt quá
trình học tập và nghiên cứu đề tài luận văn.
Cuối cùng xin được bày tỏ lòng biết ơn tới gia đình, bạn bè và các bạn học
viên lớp cao học Vật lý Chất rắn K25 đã dành tình cảm, luôn động viên và giúp
đỡ tôi để vượt qua những khó khăn để hoàn thành luận văn này.

Thái Nguyên, tháng 6 năm 2019
Học viên

Vũ Thị Thanh Hương

MỤC LỤC
Danh mục các ký hiệu và chữ viết tắt ............................................................... ii
Danh mục các bảng biểu................................................................................... iii
Danh mục các hình vẽ và đồ thị .........................................................................v
MỞ ĐẦU............................................................................................................1
Chương 1. TỔNG QUAN VẬT LIỆU KIM LOẠI .......................................3
1.1. Cấu trúc và tính chất vật liệu kim loại vô định hình……………… ...........3
1.2. Cấu trúc tinh thể kim loại lỏng……………………....................................4
1.3. Mô phỏng cấu trúc và tính chất động học trong kim loại Ni ......................7
Chương 2. PHƯƠNG PHÁP NGHIÊN CỨU ..............................................12
2.1. Phương pháp mô phỏng vật liệu................................................................12
2.2. Phương pháp dựng mẫu Ni........................................................................18
2.3. Phương pháp xác định hàm phân bố xuyên tâm........................................19
2.4. Phương pháp xác định cấu trúc đơn giản ..................................................23
2.5. Phương pháp xác cấu trúc tinh thể ............................................................25
Chương 3. KẾT QUẢ VÀ THẢO LUẬN .....................................................30
3.1. Khảo sát cấu trúc Ni theo nhiệt độ và mức độ hồi phục ...........................30
3.2. Cơ chế tạo pha thủy tinh (rắn), tạo pha tinh thể và tạo pha lỏng ..............37
3.3. Khảo sát cơ chế khuếch tán trong kim loại Ni ..........................................47
KẾT LUẬN......................................................................................................50
TÀI LIỆU THAM KHẢO..............................................................................52


DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT
Chữ viết tắt

Chữ viết đầy đủ

CST

Tứ diện cầu

DCBPTB

Dịch chuyển bình phương trung bình

ĐLHPT

Động lực học phân tử

DSR

Phân bố bán kính (Distribution of simplex radius)

ĐVCT

Đơn vị cấu trúc

EAM

Thế tương tác nhúng

FCC

Lập phương tâm mặt

HPBXT

Hàm phân bố xuyên tâm

SPT

Số phối trí

TKHP

Thống kê hồi phục

TSCT

Thừa số cấu trúc

TSKG

Thừa số không Gauss

VĐH

Vô định hình

Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên




DANH MỤC CÁC BẢNG BIỂU
Trang
Bảng 1.1.

Vị trí và tỉ lệ độ cao các đỉnh trong HPBXT cặp, g(r) của một số

4

kim loại VĐH.
Bảng 1.2.

Khối lượng riêng và một số dữ liệu của hệ số cấu trúc a(K) của

6

một số kim loại, trong đó K1 và K2 là vị trí đỉnh cao thứ nhất và
thứ hai a(K)
Bảng 2.1.

Hệ số của thế tương tác cặp nguyên tử Pak-Doyama.

19

Bảng 2.2.

Phương pháp CNA nhận diện cấu trúc tinh thể. Ví dụ: nguyên

26

tử được phối trí bởi mạng hcp có 6 liên kết loại (421) và 6 liên
kết loại (422). Tức là bất kì hai nguyên tử lân cận nào trong tinh
thể hcp cũng có 4 lân cận chung mà giữa chúng có 2 liên kết
trong đó có thể là 1 (421) hoặc liên kết 2 (422).
Bảng 2.3.

Quan hệ giữa các loại biểu đồ CNA trong các cấu trúc tinh thể

29

fcc, bcc và hcp.
Bảng 3.1.

Đặc trưng cấu trúc của Ni vô định hình ở 300 K được so sánh

30

với số liệu thực nghiệm trong các công trình [1], [5]. Ở đây ri
and gi là vị trí và độ cao của đỉnh thứ i trong HPBXT, Z là số
phối trí trung bình.

Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên




DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ
Trang
Hình 1.1.

Thừa số cấu trúc a(K) của một số kim loại VĐH theo dữ liệu

3

thực nghiệm của Leung P.K.

Hình 1.2.

Hệ số cấu trúc của Ga lỏng ở nhiệt độ khác nhau: a) Trường

5

hợp 323 K, b) Trường hợp 333 K .
Hình 1.3.

HPBXT của Fe lỏng và VĐH, đường 1 là Fe VĐH, đường 2

6

là Fe lỏng ở nhiệt độ 1560 K .
Hình 1.4.

Ảnh 3D biểu diễn các mầm và sự phát triển của các pha tinh

9

thể: a) t= 5 ps; b) t= 10 ps; c) t= 110 ps; d) t= 120 ps; e) t=
130 ps; f) t= 150 ps. Từ a) đến c) sự tạo thành mầm bởi sự
sắp xếp lại do khuếch tán, d) và e) sự phát triển tinh thể bởi
chuyển động tập thể của nguyên tử, f) liên kết các đám khác
nhau thành tinh thể với các hướng tinh thể khác nhau.
Hình 2.1.

Sơ đồ khối mô tả phương pháp ĐLHPT

15

Hình 2.2.

Sơ đồ khối mô tả phương pháp TKHP.

18

Hình 2.3.

Các loại ĐVCT: a) 4- ĐVCT, b) 5- ĐVCT, c) 6- ĐVCT, d) 7-

24

ĐVCT, e) 8– ĐVCT
Hình 2.4.

Minh hoạ cấu trúc địa phương của cặp nguyên tử được đặc

27

trưng bởi ba thông số kmn. Các cặp nguyên tử màu xám, các
nguyên tử lân cận chung màu trắng .
Hình 2.5.

Minh hoạ cấu trúc địa phương của cặp nguyên tử được đặc

28

trưng bởi các thông số hkmn. Các cặp nguyên tử màu
xám,các nguyên tử lân cận chung màu trắng

Hình 3.1

Hàm phân bố xuyên tâm đối với các mẫu 300 K; 500 K ở các

31

bước hồi phục 106, 3×106, 7×106 bước và số liệu thực nghiệm.
Hình 3.2.

Hàm phân bố xuyên tâm đối với các mẫu 650 K; 750 K ở các

32

bước hồi phục 106, 2×106, 3×106, 4×106 , 5×106 , 6×106 bước.
Hình 3.3.

Hàm phân bố xuyên tâm đối với các mẫu 900 K; 1000 K ở

34

các bước hồi phục 106, 3×106, 7×106 bước.
Hình 3.4.

Phân bố số phối trí (SPT) đối với các mẫu VĐH 300 K, 500

34

K tại các bước hồi phục 106, 3×106, 7×106 bước.
Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên

Hình 3.5.

Phân bố số phối trí (SPT) đối với các mẫu 650 K, 750 K tại

36

106, 2×106, 3×106, 4×106 , 5×106 , 6×106 bước.
Hình 3.6.

Phân bố số phối trí (SPT) đối với các mẫu lỏng 900 K, 1000

36

Hình 3.7.

K tại các bước hồi phục 106, 3×106, 7×106 bước.
Sự phụ thuộc của năng lượng theo thời gian đối với các mẫu

37

300 K, 1000 K.
Hình 3.8.

Sự phụ thuộc của năng lượng theo thời gian đối với các mẫu

38

650 K, 750 K.
Hình 3.9.

Trực quan hóa 3 chiều về sự sắp xếp các nguyên tử trong

40

mẫu 750 K: (A), (B) và (C) vẽ mẫu ở các thời gian hồi phục
106, 3×106 và 7×106 bước
Hình 3.10.
Hình 3.11.

Hàm phân bố xuyên tâm của các mẫu Ni tại các nhiệt độ khác
nhau.
Phân bố số phối trí (SPT) của các mẫu Ni tại các nhiệt độ

41
41

khác nhau.
Hình 3.12.

Sự phụ thuộc nhiệt của năng lượng đối với các mẫu trong

42

khoảng nhiệt độ 118 K đến 1380 K.
Hình 3.13.

Sự phụ thuộc nhiệt độ của tỉ số gmin/gmax đối với các mẫu

42

Hình 3.14.

trong khoảng nhiệt độ 118 K đến 1380 K.
Phân bố bán kính của ĐVCT đối với các pha VĐH, pha tinh

43

thể và pha lỏng..
Hình 3.15.

Phân bố khoảng cách cạnh, aS của ĐVCT đối với các pha

44

Hình 3.16.

VĐH, pha tinh thể và pha lỏng
Trực quan hóa 3 chiều về sự sắp xếp nguyên tử trong các

46

mẫu 1380 K, 995 K và 873 K: (A), (B) và (C) vẽ mẫu với
SPT khác 14.
Hình 3.17.

Trực quan hóa 3 chiều về sự sắp xếp nguyên tử trong các

46

mẫu 1380 K, 995 K và 873 K: (D), (E) và (F) vẽ mẫu với
nguyên tử có SPT bằng 14
Hình 3.18.

Khoảng cách dịch chuyển bình phương trung bình theo thời

48

gian mô phỏng đối với các pha VĐH (118 K, 336 K, 516 K,
642 K) và pha tinh thể (873 K).
Hình 3.19.

Khoảng cách dịch chuyển bình phương trung bình theo thời
gian mô phỏng đối với pha lỏng (995 K, 1167 K, 1380 K).

Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên

48




MỞ ĐẦU
1. Lý do chọn đề tài
Niken (Ni) là vật liệu kim loại có nhiều đóng góp quan trọng trong lĩnh
vực khoa học ứng dụng và khoa học vật liệu. Nhiều nghiên cứu về cấu trúc và
sự biến đổi cấu trúc đã được thực hiện bằng cả thực nghiệm, lý thuyết và mô
phỏng. Các kết quả nghiên cứu thực nghiệm bằng cách sử dụng phương pháp

tán xạ raman, nhiễu xạ tia X, nhiễu xạ neutron, hiển vi điện tử truyền qua,…
đã cung cấp nhiều dữ liệu quan trọng về cấu trúc của các vật liệu kim loại
thông qua phân tích hàm phân bố xuyên tâm (HPBXT) hoặc thừa số cấu trúc.
Những kết quả nghiên cứu thực nghiệm đã chỉ ra rằng, ở trạng thái vô định
hình (VĐH), đỉnh thứ hai của HPBXT tách thành hai đỉnh nhỏ, sự tách đỉnh
này thường được cho là liên quan đến sự tồn tại của các khối đa diện 20 mặt.
Bằng phương pháp mô phỏng, HPBXT của kim loại nhôm bắt đầu có sự tách
đỉnh ở cực đại thứ hai ở dưới nhiệt độ 500 K. Trong suốt quá trình làm lạnh từ
từ kim loại Nhôm (Al) từ 700 xuống 300 K thì thấy xuất hiện ba trạng thái:
lỏng, tinh thể yếu và trạng thái VĐH. Tuy nhiên, trong các nghiên cứu trình
bày ở trên vẫn chưa chỉ ra được cơ chế khuếch tán đối với kim loại lỏng và
VĐH. Do đó cơ chế động học trong các kim loại vẫn đang là vấn đề còn nhiều
tranh luận và cần có nghiên cứu mới. Lý do trên khiến chúng tôi chọn đề tài
“Mô phỏng đặc trưng cấu trúc và cơ chế khuếch tán trong vật liệu Ni lỏng và
rắn”. Ở đây, chúng tôi thực hiện nghiên cứu cấu trúc và cơ chế khuếch tán đối
với kim loại Ni trong khoảng nhiệt độ từ 1380 K- 118 K bởi thế tương tác cặp
Pak- Doyama.
2. Mục tiêu đề tài
Tạo các mẫu Ni chứa 4394 và 10000 nguyên tử với mật độ thực của Ni
trong khoảng nhiệt độ từ 1380 K đến 118 K. Phân tích vi cấu trúc của Ni trên
các
Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên




mẫu tạo được thông qua các đại lượng như HPBXT, phân bố số phối trí (SPT),
thống kê đơn vị cấu trúc (ĐVCT) và sử dụng phương pháp trực quan hóa.
Đánh giá nhiệt độ chuyển pha từ pha lỏng sang pha VĐH. Qua đó phân
tích cơ chế khuếch tán của Ni ở trạng thái lỏng và rắn.

3. Phương pháp nghiên cứu
Đề tài sử dụng phương pháp thống kê hồi phục, phương pháp động lực
học phân tử và các phương pháp phân tích cấu trúc vật liệu.
4. Cấu trúc đề tài
Đề tài gồm phần mở đầu, chương 1, 2, 3, kết luận và tài liệu tham khảo.
Trong đó, chương 1 trình bày tổng quan về vật liệu kim loại như: cấu trúc, tính
chất, ứng dụng của kim loại vô định hình, kim loại lỏng, các kết quả mô phỏng
cấu trúc và tính chất của Ni. Chương 2 trình bày phương pháp mô phỏng,
phương pháp phân tích cấu trúc kim loại. Chương 3 trình bày kết quả và những
thảo luận của đề tài.


Chương 1
TỔNG QUAN VẬT LIỆU KIM LOẠI
Nội dung chương này được trình bày như sau: đầu tiên là tổng quan về đặc
trưng cấu trúc và tính chất của vật liệu kim loại trạng thái vô định hình và lỏng.
Sau đó là phương pháp chế tạo và ứng dụng của các vật liệu kim loại. Cuối
cùng, chúng tôi trình bày các kết quả mô phỏng về kim loại Ni.
1.1. Cấu trúc và tính chất vật liệu kim loại vô định hình
Cấu trúc: Hình 1.1 mô tả thừa số cấu trúc (TSCT) của một số màng mỏng
kim loại VĐH từ số liệu thực nghiệm của P.K. Leung [8,10]. Từ hình ảnh phổ
trên hình 1.1 có thể thấy, vị trí và độ cao cực đại của TSCT đặt trong khoảng 29
đến 33 (nm-1) và từ 3,6 đến 4,2. Đỉnh thứ hai (cực đại thứ hai) trong TSCT bị
tách thành hai đỉnh nhỏ, cụ thể là độ cao đỉnh nhỏ bên trái cao hơn độ cao đỉnh
nhỏ bên phải.

Hình 1.1. Thừa số cấu trúc a(K) của một số kim loại
VĐH theo dữ liệu thực nghiệm của Leung P.K. [15]

Tương tự đối với TSCT, HPBXT g(r) của kim loại VĐH, đỉnh thứ hai của
cũng bị tách thành hai đỉnh nhỏ, độ cao đỉnh nhỏ bên trái cao hơn độ cao đỉnh
nhỏ bên phải. Hình dạng TSCT cũng như HPBXT của các kim loại có dạng khá
giống nhau. Vị trí các cực đại trong HPBXT được liệt kê trong bảng 1.1. Từ
bảng 1.1 có thể thấy, các tỉ số ri/r1 (i = 2, 3, 4 và 5) đối với các kim loại VĐH
khác nhau là bằng nhau.
Bảng 1.1. Vị trí và tỉ lệ độ cao các đỉnh trong HPBXT cặp, g(r) của một số
kim loại VĐH công bố trong các công trình [1,2,6,15,20].
Kim loại

r2/r1

r3/r1

r4/r1

r5/r1

Cr

1,66

1,91

2,49

3,36

Fe

1,67

1,96

2,51

3,38

Co

1,96

1,93

2,49

3,35

Ni

1,71

1,93

2,55

3,38

Tính chất: Như chúng ta biết, kim loại ở trạng thái VĐH có nhiều đặc
điểm khác với trạng thái tinh thể. Ở trạng thái tinh thể, số lượng tâm tán xạ tăng

khi nồng độ khuyết tật tăng. Trái lại, với kim loại VĐH, tính tuần hoàn bị mất,
nên mọi nguyên tử đều trở thành tâm tán xạ, do đó các tính chất dẫn điện, dẫn
nhiệt và các hiệu ứng trong kim loại VĐH cũng bị thay đổi đáng kể. Các tính
chất khác như lực kháng từ, tính chống ăn mòn, điện trở suất, độ bền cơ học, độ
cứng của kim loại VĐH đều có giá trị cao hơn hẳn so với kim loại này khi ở
trạng thái tinh thể.
1.2. Cấu trúc và tính chất của kim loại lỏng
Như đã biết, khi phân tích cường độ chùm tia tán xạ qua kim loại lỏng đơn
nguyên tử người ta có thể xác định được đặc trưng cấu trúc của kim loại. Ở hình
1.2 cho thấy, TSCT của kim loại Ga ở nhiệt độ 323 K đến 333 K. Nhận thấy
đường cong TSCT của kim loại Ga lỏng thay đổi đáng kể khi nhiệt độ thay đổi
từ 323 K đến 333 K. Cụ thể, độ cao đỉnh thứ nhất tăng từ 2,45 đến 2,75, hình


dạng của đỉnh trồi ra sau cực đại chính là khá rõ ràng. Hình 1.3 cho thấy
HPBXT ở trạng thái lỏng và VĐH của kim loại Fe là khác nhau đáng kể. Độ
cao cực đại thứ nhất của HPBXT của Fe lỏng là 2,5 thấp hơn nhiều so với cực
đại thứ nhất HPBXT của Fe VĐH là 3,20. Đối với Fe VĐH, cực đại thứ hai bị
tách thành hai đỉnh nhỏ, độ cao đỉnh nhỏ bên trái cao hơn độ cao đỉnh nhỏ bên
phải. Trong khi đó, với Fe lỏng thì không xuất hiện sự tách đỉnh ở cực đại thứ
hai.
Kim loại lỏng không có trật tự xa trong phân bố nguyên tử như trong trạng
thái tinh thể. Do vậy, kim loại lỏng thuộc nhóm vật liệu có cấu trúc không trật
tự. Tuy nhiên trong giới hạn của mặt cầu phối vị thứ nhất thì phân bố nguyên tử
trong kim loại lỏng cũng có một trật tự đặc trưng cho từng chất khác nhau và
được gọi là trật tự gần.
Bởi vì kim loại lỏng có cấu trúc không trật tự nên kim loại lỏng có các
tính chất rất khác thường về độ dẫn điện, từ tính và một số tính chất cơ học
khác. Ví dụ, ở trạng thái lỏng các kim loại đất hiếm Eu, Er, Gd có độ dẫn điện
thấp hơn đáng kể so với các kim loại khác. Giá trị điện trở suất của Er lớn hơn

gấp mười lần so với giá trị điện trở suất của kim loại Cu. Nhìn chung, điện trở
suất của kim loại lỏng lớn hơn rất nhiều so với kim loại tương ứng khi chúng ở
trạng thái rắn [1].

Hình 1.2. Hệ số cấu trúc của Ga lỏng ở nhiệt độ khác
nhau: a) Trường hợp 323 K, b) Trường hợp 333 K [3].


Các kết quả thực nghiệm cho thấy TSCT của nhiều kim loại lỏng là rất
giống nhau như có thể thấy trong bảng 1.2. Điều này có nghĩa là cấu trúc của
kim loại lỏng nói chung là khá tương tự nhau. Nếu như dựng đồ thị của hệ số
cấu trúc theo đơn vị K/K1, với phần lớn kim loại lỏng như Hg, Al, Pb, Fe và Co thì
độ cao đỉnh đầu tên gần như trùng khít và độ cao đỉnh thứ hai thì rất giống
nhau.

Hình 1.3. HPBXT của Fe lỏng và VĐH, đường 1 là Fe
VĐH, đường 2 là Fe lỏng ở nhiệt độ 1560 K [9].
Bảng 1.2. Khối lượng riêng và một số dữ liệu của hệ số cấu trúc a(K) của một
số kim loại , trong đó K1 và K2 là vị trí đỉnh cao thứ nhất và thứ hai a(K)
Kim loại

T, 0C

d, g/cm3

K1, nm-1

Li

180

0,504

24,9

Đỉnh phụ,
nm-1
-

Na

105

0,928

20,3

Al

670

2,37

Ge

980

Pb

K2, nm-1

K2/K1

45,5

1,83

-

37,5

1,85

26,8

-

49,6

1,85

5,56

25,6

32,4

51,1

2,00

340

10,66

22,8

-

42,3

1,86

Fe

1550

7,01

29,8

-

54,6

1,83

Co

1550

7,70

30,2

-

56,0

1,85


1.3. Mô phỏng cấu trúc và tính chất động học trong kim loại Ni
Trong mục này, chúng tôi trình bày chi tiết một số kết quả mô phỏng
ĐLHPT về cấu trúc, các tính chất động học trên trong kim loại Ni thông qua
một số công trình đã được công bố, cụ thể như sau:
Zhang Yanning và các cộng sự [21] đã nghiên cứu tính chất nhiệt động
học, động học và sự hồi phục cấu trúc trong Ni ở trạng thái thủy tnh và trạng
thái siêu lỏng dưới nhiệt độ tới hạn bằng phương pháp mô phỏng ĐLHPT dựa
trên thế tương tác nhúng với mẫu 4000 nguyên tử. Họ đã sử dụng các phương
pháp: (1) Phân tích cấu trúc thông qua hàm tương quan cặp; (2) Xác định ranh
giới giữa vùng siêu lỏng và vùng vô định hình thông qua tham số WendtAbraham; (3) Sử dụng các chỉ số Honeycutt- Anderson (HA) index; (4) Sử
dụng các hàm tương quan thời gian. Kết quả mô phỏng của Zhang Yanning và
các cộng sự chỉ ra rằng, nhiệt độ chuyển thủy tinh của Ni là Tg = 800 K và nhiệt
độ tới hạn được xác định bằng 840 K, ở đây tốc độ làm lạnh 2,86×1013 K/s. Họ
phát hiện ra ba vùng đặc biệt trong suốt quá trình hồi phục với các tính chất
nhiệt động, động học và cấu trúc rất khác biệt. Trong vùng nhiệt độ 840 K1000 K sự tinh thể trong trạng thái siêu lỏng xảy ra với thời giản ủ nhiệt là τ; ở
trong vùng nhiệt độ thấp hơn 810 K- 840 K, giá trị nội năng cao hơn đáng kể
trong thời gian ủ dài và trạng thái cuối cùng của hệ được đặc trưng bởi cấu
trúc khá phức tạp, kết quả này được cho là do tồn tại sự pha trộn giữa các

phần tnh thể và thủy tnh. Trong khi đó, khi nhiệt độ dưới 810 K, sự tinh
thể hóa của thủy tnh Ni tăng và thời gian τ kéo dài nhanh chóng với sự giảm
nhiệt độ. Tác giả và cộng sự chỉ ra rằng dáng điệu động học xảy ra trong vùng
tinh thể trong trạng thái siêu lỏng có thể so sánh với vùng tinh thể hóa của
thủy tnh nhận được và nó được đặc trưng bởi sự xuất hiện khửu (đường
nằm ngang) thứ hai trong dịch chuyển bình phương trung bình (DCBPTB) và
đường ngang (khửu) không bằng không của thừa số không Gauss (TSKG) với
thời gian ủ lớn. Điều


này có nghĩa là, quá trình tnh thể từ trạng thái siêu lỏng hoặc thủy tinh thông
qua sự mầm hóa của các nhân tnh thể. Sự hồi phục ở nhiệt độ gần Tg chịu ảnh
hưởng rất phức tạp theo sự biến đổi động học. Sự hồi phục nhiệt động học và
động học kéo theo sự sắp xếp lại cấu trúc nguyên tử. Cặp liên kết 1551
đặc trưng đối xứng 5 địa phương giảm trong quá trình hồi phục; tương ứng,
cặp liên kết 1421 mô tả trật tự cấu trúc tăng. Thời gian ủ có thể được xác
định khi số các cặp liên kết thay đổi hình dạng mạnh, trong khi “vùng phức
tạp” chỉ biến đổi nhỏ trong cấu trúc hệ được tìm thấy.
Bằng phương pháp mô phỏng, L.Wang và cộng sự [19] cho thấy mối liên
hệ giữa sự mầm hóa đến cấu trúc và động học không đồng nhất trong kim loại
Ni trạng thái siêu lạnh. Tính chất động học và cấu trúc không đồng nhất
trong Ni trạng thái siêu lạnh đã được xác định với mô phỏng ĐLHPT với thế
tương tác nhiều hạt (thế nhúng- EAM). Đặc trưng về thời gian hồi phục đã
được đánh giá bởi các hàm tương quan thời gian, DCBPTB và TSKG. Họ đã
phát hiện, sau khi trải qua đường ngang ra đầu tiên, đường cong DCBPTB
cho thấy đường ngang thứ hai xuất hiện tại nhiệt độ 1090 K. Ở khoảng thời
gian đầu sự hồi phục cấu trúc, phôi xuất hiện bởi sự sắp xếp lại do khuếch tán
của các nguyên tử. Khi thời gian tăng, mầm tnh thể được tạo thành trong
trạng thái lỏng và chúng lớn dần bởi chuyển động tập thể của các nguyên tử
để tạo thành các đám lớn, kết quả này nhận được dựa trên đặc trưng của sự

xuất hiện điểm ngang khác không trong đường cong TSKG và hàm tương
quan thời gian. Ở thời gian ủ dài hơn, đám tnh thể lớn với hướng tinh thể
khác nhau liên kết với nhau, sau đó tinh thể tạo thành.
Sử dụng phương pháp phân tích cặp và phương pháp kiểm tra trực quan
L.Wang và cộng sự đã phân tích quá trình mầm trong Ni trạng thái siêu lỏng.
Họ cho rằng, các đám mầm với kích thước nhỏ kết nối với nhau và trở
Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên




thành một đám lớn, tinh thể tạo thành cuối cùng có cấu trúc tinh thể . Trên
hình 1.4, tác giả và cộng sự cho thấy ảnh chụp 3D (3 chiều) về sự tạo thành
mầm và phát

Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên




triển mầm tnh thể. Họ chỉ ra rằng, ở trạng thái lỏng trong các số hạng của lý
thuyết mầm cổ điển, có nhiều đám gọi là phôi (embryo) như thấy trong
hình
1.4a và 1.4b, kích thước của chúng là khá nhỏ, các nguyên tử trong phôi
tập hợp (liên kết) với nhau vào cùng một thời điểm bằng cách sắp xếp
lại do khuếch tán. Tuy nhiên, không phải tất cả các phôi đều phát triển
thành mầm, mà chỉ có phôi nào mà kích thước của chúng đạt kích thước tới
hạn mới có thể trở thành mầm như thấy trong hình 1.4c. Một mầm xuất hiện
trong trạng thái lỏng, chúng phát triển lên là do sự chuyển động tập thể của
các nguyên tử và sau đó liên kết với nhau trở thành một tinh thể lớn. Cuối

cùng, quá trình này được minh họa bởi hình 1.4d, 1.4e và 1.4f.

Hình 1.4. Ảnh 3D biểu diễn các mầm và sự phát triển của các pha tinh thể: a)
t= 5 ps; b) t= 10 ps; c) t= 110 ps; d) t= 120 ps; e) t= 130 ps; f) t= 150 ps. Từ a)
đến c) sự tạo thành mầm bởi sự sắp xếp lại do khuếch tán, d) và e) sự
phát
triển tinh thể bởi chuyển động tập thể của nguyên tử, f) liên kết các đám
khác


Số hóa bởi Trung tâm Học liệu và Công nghệ thông tn – Đại học Thái Nguyên

nhau thành tinh thể với các hướng tinh thể khác nhau
[19].

htp://lrc.tnu.edu.vn


Theo tác giả và cộng sự, đám nguyên tử có các nhân với cấu trúc tinh thể
tạo bởi 12 cặp liên kết 1421 đóng vai trò quan trọng trong quá trình mầm
hóa và phát triển thành tinh thể trong quá trình làm lạnh, sau quá trình
phát triển này, cuối cùng đám tnh thể thống trị cấu trúc của hệ. Mô phỏng
của họ đã cung cấp một sự hiểu biết tương đối đầy đủ về động học và cấu trúc
trong trạng thái siêu lỏng phát triển thành tinh thể.
Y.Zhang và cộng sự [22] sử dụng mô phỏng ĐLHPT dựa trên thế tương
tác nhiều hạt (thế nhúng - EAM) để nghiên cứu sự khác nhau về sự không đồng
nhất cấu trúc trong quá trình tnh thể hóa và quá trình chuyển thành pha thủy
tinh. Quá trình hồi phục động học đối với Ni siêu lỏng theo các quá trình giảm
nhiệt khác nhau đã được thực hiện. Khảo sát các kết quả về thế năng
và HPBXT nhận được từ dữ liệu mô phỏng, họ đã khẳng định rằng hệ

chuyển thành trạng thái tinh thể khi tốc độ làm lạnh bằng 1,9×1012 K/s, ứng
với nhiệt độ tinh thể là Ts = 990 K. Trong khi đó, hệ chuyển thành trạng thái
thủy tnh xảy ra khi tốc độ làm lạnh tăng đến 5,7×1012 K/s ứng với nhiệt độ
chuyển pha thủy tinh Tg = 670 K. Theo Y.Zhang và cộng sự, trong quá trình
tnh thể và chuyển thủy tnh, sự khác nhau về các hàm tương quan thời gian
DCBPTB và TSKG rất phù hợp với nhau và cho thấy sự không đồng nhất động
học cũng rất khác nhau.
Theo Y.Zhang và cộng sự trong suốt quá trình hồi phục cho thấy, quá
trình tinh thể hóa và quá trình chuyển thành pha thủy tnh có mối liên hệ chặt
chẽ với nhau. Đối với quá trình tnh thể hóa từ trạng thái nóng chảy, Y.Zhang
và cộng sự cho rằng sự tạo thành mầm là do có sự sắp xếp lại nhưng không
khuếch tán của các nguyên tử, mà chúng được đặc trưng bởi sự xuất hiện
điểm ngang thứ hai trong DCBPTB. Tuy nhiên, trong suốt quá trình chuyển
Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên




sang pha thủy tnh, chuyển động của nguyên tử được cho là do có sự khuếch
tán ngay cả trong trạng thái thủy tnh. Thêm nữa, điểm ngang thứ hai
trong DCBPTB

Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên




không quan sát thấy và đặc điểm không Gauss trong đồ thị TSKG trở lên rất rõ
nét khi ở nhiệt độ thấp.
Các kết quả phân tích trên cho thấy đặc trưng cấu trúc, cơ chế khuếch tán

trong kim loại Ni trạng thái lỏng, tinh thể và VĐH vẫn còn nhiều vấn đề cần
được làm rõ thêm và đây cũng chính là nội dung và mục têu mà đề tài đặt
ra cần được nghiên cứu.

Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên




Chương 2
PHƯƠNG PHÁP NGHIÊN CỨU
Nội dung trong chương này, chúng tôi trình bày phương pháp động
lực học phân tử, phương pháp thống kê hồi phục. Tiếp theo, chúng tôi trình
bày phương pháp xây dựng mẫu Ni. Sau đó, chúng tôi trình bày phương pháp
xác định hàm phân bố xuyên tâm. Cuối cùng chúng tôi trình bày cách xác định
các đơn vị cấu trúc đơn giản trong mẫu mô phỏng.
2.1. Phương pháp mô phỏng vật liệu
Phương pháp động lực học phân tử: Phương pháp ĐLHPT là một công
cụ cho phép chúng ta xây dựng mô hình vật liệu nano dựa trên hệ phương
trình chuyển động của Newton. Phương trình chuyển động được khảo sát với
vận tốc chuyển động của hạt tính bằng thuật toán Verlet theo bước thời gian
dt [11].
Xét một hệ gồm N nguyên tử được gieo trong khối hình lập phương
cạnh L. Tọa độ ban đầu của nguyên tử có thể lấy ngẫu nhiên nhưng phải thỏa
mãn điều kiện không có bất kỳ hai nguyên tử nào quá gần nhau. Dưới tác dụng
của lực tương tác, các nguyên tử sẽ dịch chuyển dần đến vị trí cân bằng. Trạng
thái cân bằng của mô hình được xác định bởi nhiệt độ và áp suất. Chuyển
động của các nguyên tử trong mô hình tuân theo định luật cơ học cổ điển
Newton. Đối với hệ gồm N hạt, phương trình chuyển động của định luật II
Newton có thể viết như sau:

Phương pháp động lực học phân tử cổ điển dựa trên phương trình chuyển
động Newton:

Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên




Fi  mi a i .

(2.1)

d2
mi ai  mi r 2 i  F i r1...., rN  .
dt

(2.2)

Số hóa bởi Trung tâm Học liệu và Công nghệ thông tin – Đại học Thái Nguyên