Nghiên cứu ảnh hưởng của biên cứng lên sự thủy tinh hóa của nickel lỏng được làm lạnh nhanh

Tài liệu chuyên sâu Ảnh hưởng của biên cứng đến thủy tinh hóa nickel lỏng ..., phân tích đa chiều, cung cấp kiến thức nền tảng vững chắc cho chuyên

• Tác giả: Văng Vũ Linh
• Người hướng dẫn: GS. Võ Văn Hoàng
• Trường học: Đại học Bách Khoa
• Chuyên ngành: Vật lý kỹ thuật
• Đề tài: Ảnh Hưởng Của Biên Cứng Lên Sự Thủy Tinh Hóa Của Nickel Lỏng Được Làm Lạnh Nhanh
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2013
• Địa điểm: TP. Hồ Chí Minh

Mục lục chi tiết

• DANH MỤC CÁC KÝ HIỆU VÀ CHỮ VIẾT TẮT
• DANH MỤC CÁC HÌNH VẼ, ĐỒ THỊ
• DANH MỤC CÁC BẢNG BIỂU
• 1. CHƯƠNG 1: MỞ ĐẦU
• 1.1. Các nghiên cứu về vật liệu có biên cứng hoặc có không gian bị hạn chế
• 1.2. Các nghiên cứu về Ni vô định hình, hợp kim vô định hình khác
• 2. CHƯƠNG 2: TÍNH TOÁN MÔ PHỎNG
• 2.1. Phương pháp động lực học phân tử
• 2.1.1. Cơ sở của phương pháp mô phỏng
• 2.1.2. Phương pháp động lực học phân tử
• 2.2. Chi tiết tính toán mô phỏng
• 2.2.1. Thế tương tác EAM cho Ni và kim loại tương tự
• 2.2.2. Các chi tiết về mô hình Ni có biên cứng được mô phỏng
• 3. CHƯƠNG 3: KẾT QUẢ VÀ THẢO LUẬN
• 3.1. Các tính chất nhiệt động học
• 3.1.1. Khoảng cách trung bình của các nguyên tử Ni ở trạng thái lỏng và vô định hình
• 3.1.2. Sự phụ thuộc vào nhiệt độ của năng lượng toàn phần cho mỗi nguyên tử
• 3.1.3. Sự thay đổi theo thời gian của trung bình bình phương độ dịch chuyển
• 3.1.4. Sự phân bố mật độ dọc theo trục Z ở những nhiệt độ khác nhau
• 3.2. Những thay đổi của cấu trúc mô hình khi làm lạnh từ lỏng sang vô định hình
• 3.2.1. Khảo sát hàm phân bố xuyên tâm của mô hình với biên cứng khi làm lạnh từ trạng thái lỏng
• 3.2.2. Sự thay đổi của tỷ lệ các cặp liên kết Honeycutt–Andersen
• 3.2.3. Cấu trúc vi mô của mô hình nhận được ở nhiệt độ 300K
• 3.3. Cơ chế nguyên tử của sự hình thành trạng thái vô định hình từ trạng thái lỏng
• 3.3.1. Sự phụ thuộc vào nhiệt độ của tỷ số Lindemann
• 3.3.2. Sự phụ thuộc vào nhiệt độ của tỷ lệ các nguyên tử dạng rắn và tỷ lệ các nguyên tử dạng lỏng
• 3.3.3. Quá trình hình thành và tích tụ trong không gian của những nguyên tử đã chuyển sang trạng thái rắn trong quá trình làm lạnh
• 3.3.4. Sự phân bố mật độ dọc theo trục Z của các nguyên tử dạng rắn, các nguyên tử dạng lỏng
• 3.4. Sự ảnh hưởng của biên cứng lên cấu trúc của Ni vô định hình
• 4. CHƯƠNG 4: KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN ĐỀ TÀI
• 4.1. Hướng phát triển của đề tài
• DANH MỤC CÁC BÀI BÁO KHOA HỌC ĐÃ CÔNG BỐ
• TÀI LIỆU THAM KHẢO
• PHẦN LÝ LỊCH TRÍCH NGANG

Tóm tắt

I. Giới thiệu

Nghiên cứu về biên cứng và ảnh hưởng của nó đến quá trình thủy tinh hóa của nickel lỏng là một lĩnh vực quan trọng trong vật lý vật liệu. Thủy tinh hóa là quá trình chuyển đổi từ trạng thái lỏng sang trạng thái vô định hình, và biên cứng có thể ảnh hưởng đáng kể đến quá trình này. Mô phỏng động lực học phân tử (MD) được sử dụng để khảo sát ảnh hưởng của biên cứng lên cấu trúc và tính chất nhiệt động học của nickel lỏng. Kết quả cho thấy rằng biên cứng làm tăng nhiệt độ chuyển pha và tạo ra cấu trúc lớp tại nơi tiếp giáp với nickel lỏng. Điều này có thể giải thích sự xuất hiện của các khuyết tật và sự thay đổi trong động học của các nguyên tử ở lớp gần biên.

1.1. Tầm quan trọng của nghiên cứu

Nghiên cứu này không chỉ giúp hiểu rõ hơn về cơ chế thủy tinh hóa mà còn mở ra hướng đi mới cho việc phát triển các vật liệu mới. Nickel là một trong những kim loại có ứng dụng rộng rãi trong công nghiệp, và việc hiểu rõ về tính chất của nó trong trạng thái vô định hình sẽ giúp cải thiện chất lượng và hiệu suất của các sản phẩm. Hơn nữa, việc khảo sát ảnh hưởng của biên cứng đến thủy tinh hóa có thể cung cấp thông tin quý giá cho các nghiên cứu về vật liệu nano và các ứng dụng trong công nghệ cao.

II. Phương pháp nghiên cứu

Phương pháp động lực học phân tử (MD) được áp dụng để mô phỏng quá trình làm lạnh nhanh của nickel lỏng từ trạng thái lỏng sang trạng thái vô định hình. Mô hình sử dụng có biên cứng và tốc độ làm lạnh 5x10^13 K. Các thông số như tính chất vật liệu, nhiệt độ thủy tinh hóa, và cấu trúc vi mô được khảo sát để phân tích ảnh hưởng của biên cứng. Kết quả cho thấy rằng biên cứng làm giảm động học của các nguyên tử ở lớp gần biên, dẫn đến sự hình thành nhiều khuyết tật. Điều này cho thấy rằng biên cứng không chỉ ảnh hưởng đến thủy tinh hóa mà còn đến các tính chất cơ học của nickel.

2.1. Mô hình và điều kiện thí nghiệm

Mô hình được xây dựng với 8788 nguyên tử nickel và được làm lạnh từ trạng thái lỏng hoàn toàn. Các thông số như nhiệt độ, tốc độ làm lạnh, và biên cứng được điều chỉnh để khảo sát ảnh hưởng của chúng đến quá trình thủy tinh hóa. Kết quả cho thấy rằng nhiệt độ chuyển pha của mô hình có biên cứng cao hơn so với vật liệu khối tương ứng, điều này cho thấy rằng biên cứng có vai trò quan trọng trong việc điều chỉnh các tính chất của nickel trong trạng thái vô định hình.

III. Kết quả và thảo luận

Kết quả từ mô phỏng cho thấy rằng biên cứng ảnh hưởng mạnh mẽ đến cấu trúc và tính chất nhiệt động học của nickel lỏng. Sự xuất hiện của cấu trúc lớp tại nơi tiếp giáp với biên cứng cho thấy rằng biên cứng không chỉ làm tăng nhiệt độ chuyển pha mà còn tạo ra các khuyết tật trong cấu trúc. Các nghiên cứu trước đây cũng đã chỉ ra rằng biên cứng có thể làm giảm động học của các nguyên tử, dẫn đến sự hình thành các cấu trúc lớp. Điều này có thể giải thích tại sao nickel trong trạng thái vô định hình có các tính chất khác biệt so với trạng thái lỏng.

3.1. Phân tích cấu trúc vi mô

Phân tích cấu trúc vi mô cho thấy rằng mật độ của các lớp này tăng khi nhiệt độ giảm. Sự phân bố mật độ dọc theo trục Z cho thấy rằng các nguyên tử ở gần biên cứng có xu hướng sắp xếp lại, tạo thành các lớp có mật độ khác nhau. Điều này cho thấy rằng biên cứng không chỉ ảnh hưởng đến thủy tinh hóa mà còn đến sự sắp xếp không gian của các nguyên tử trong nickel. Kết quả này có thể được áp dụng trong việc phát triển các vật liệu mới với tính chất cơ học và nhiệt động học tốt hơn.

IV. Kết luận và hướng phát triển

Nghiên cứu này đã chỉ ra rằng biên cứng có ảnh hưởng đáng kể đến quá trình thủy tinh hóa của nickel lỏng. Kết quả cho thấy rằng biên cứng làm tăng nhiệt độ chuyển pha và tạo ra cấu trúc lớp tại nơi tiếp giáp. Những phát hiện này không chỉ có giá trị trong việc hiểu rõ hơn về cơ chế thủy tinh hóa mà còn mở ra hướng đi mới cho việc phát triển các vật liệu mới. Hướng nghiên cứu tiếp theo có thể tập trung vào việc khảo sát ảnh hưởng của các yếu tố khác như áp suất và thành phần hóa học đến quá trình thủy tinh hóa.

4.1. Hướng nghiên cứu tiếp theo

Hướng nghiên cứu tiếp theo có thể bao gồm việc khảo sát ảnh hưởng của các yếu tố như áp suất và thành phần hóa học đến quá trình thủy tinh hóa của nickel. Việc áp dụng các phương pháp mô phỏng tiên tiến hơn có thể giúp hiểu rõ hơn về cơ chế nguyên tử của quá trình này. Ngoài ra, việc nghiên cứu các hợp kim vô định hình khác cũng có thể mang lại những hiểu biết mới về tính chất và ứng dụng của các vật liệu này trong công nghiệp.

Nội dung chính

Tổng quan nghiên cứu

Vật liệu kim loại vô định hình nói chung và màng mỏng kim loại vô định hình nói riêng đang là trọng tâm nghiên cứu của ngành vật lý chất rắn và khoa học vật liệu hiện đại nhờ sở hữu các đặc tính vượt trội về độ bền cơ học, khả năng chống ăn mòn và tính chất từ đặc biệt. Trong nhóm kim loại chuyển tiếp, Nickel giữ vị trí sản xuất và ứng dụng lớn thứ hai chỉ sau sắt, đóng vai trò then chốt trong chế tạo linh kiện điện tử, cảm biến và màng bảo vệ từ tính. Tuy nhiên, hành vi nhiệt động học và cấu trúc vi mô của Nickel lỏng khi làm lạnh nhanh trong không gian bị giới hạn bởi các mặt phẳng biên cứng vẫn là vấn đề học thuật đòi hỏi sự giải thích chi tiết ở cấp độ nguyên tử.

Nghiên cứu được thực hiện nhằm khảo sát toàn diện ảnh hưởng của điều kiện biên cứng lên quá trình thủy tinh hóa của Nickel lỏng được làm lạnh nhanh bằng phương pháp mô phỏng động lực học phân tử cổ điển. Mô hình nghiên cứu gồm chính xác 8.788 nguyên tử Nickel được khảo sát trong dải nhiệt độ rộng từ 2500 K xuống 300 K với tốc độ giảm nhiệt cực cao đạt $5 \times 10^{13}$ K/s. Trước khi hạ nhiệt, hệ được hồi phục nhiệt động học trong khoảng thời gian 500 ps tại 2500 K để đảm bảo trạng thái lỏng hoàn toàn, sau đó áp dụng hai biên cứng phản xạ lý tưởng dọc theo trục Z.

Ý nghĩa khoa học của công trình thể hiện qua việc xác định chính xác nhiệt độ chuyển pha thủy tinh hóa của hệ biên cứng đạt $T_g \approx 900$ K, tăng khoảng 50 K (tương đương mức tăng 5.88%) so với hệ vật liệu khối có cùng kích thước ($T_g \approx 850$ K). Kết quả nghiên cứu cung cấp cơ sở dữ liệu định lượng quan trọng giúp tối ưu hóa công nghệ chế tạo màng mỏng vô định hình trên đế cứng và phát triển các vật liệu xúc tác dị thể hiệu năng cao.

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

Khung lý thuyết áp dụng

Nghiên cứu xây dựng trên nền tảng ba khung lý thuyết và mô hình tính toán vật lý chất đặc hiện đại:

Thứ nhất là thuyết hàm mật độ điện tử thông qua mô hình thế tương tác nhúng nhiều hạt EAM do Daw và Baskes phát triển. Khác với các thế tương tác cặp đôi cổ điển, thế EAM mô tả tổng năng lượng của hệ bao gồm năng lượng nhúng một nguyên tử vào mật độ điện tử nền do các nguyên tử xung quanh tạo ra kết hợp với thế tương tác cặp đẩy hạt nhân trong phạm vi gần. Mô hình này phản ánh chính xác bản chất liên kết kim loại chuyển tiếp của Nickel ở cả pha lỏng và vô định hình.

Thứ hai là tiêu chuẩn Lindemann mở rộng cho quá trình chuyển pha lỏng sang vô định hình. Tỷ số Lindemann của từng nguyên tử được xác định dựa trên căn bậc hai trung bình bình phương độ dịch chuyển chia cho khoảng cách trung bình giữa các nguyên tử $r = 2.48$ Å. Tỷ số Lindemann tới hạn $\delta_c \approx 0.37$ (tương đương 37%) được sử dụng làm ngưỡng phân loại chính xác giữa các nguyên tử dạng lỏng và nguyên tử dạng rắn trong quá trình làm lạnh.

Thứ ba là phương pháp phân tích cặp liên kết Honeycutt–Andersen nhằm giải mã trật tự cấu trúc địa phương thông qua bốn chỉ số hình học. Các khái niệm cốt lõi được định lượng hóa bao gồm hàm phân bố xuyên tâm $g(r)$, trung bình bình phương độ dịch chuyển theo thời gian, số phối vị trung bình $Z$, trật tự khối 20 mặt icosahedra hoàn chỉnh đại diện bởi cặp liên kết 1551 và icosahedra biến thể đại diện bởi cặp liên kết 1541.

Phương pháp nghiên cứu

Nghiên cứu sử dụng phương pháp mô phỏng động lực học phân tử cổ điển với tập dữ liệu gồm 8.788 nguyên tử Nickel tinh khiết. Mẫu ban đầu được tạo dựng từ cấu trúc mạng lập phương tâm mặt với hằng số mạng $a_0 = 3.517$ Å, tạo thành hộp mô phỏng có kích thước ba chiều $Lx = Ly = Lz = 13a_0 = 45.721$ Å. Lý do lựa chọn kích thước mẫu này nhằm đảm bảo tính đại diện thống kê của hệ nhiều hạt, loại bỏ các sai số kích thước nhỏ trong khi vẫn tối ưu hóa tài nguyên tính toán.

Về phương pháp chọn mẫu và điều kiện biên, mô hình thiết lập điều kiện biên tuần hoàn dọc theo hai trục không gian X và Y, trong khi áp dụng điều kiện biên phản xạ tại hai vị trí $z = 0$ và $z = 45.721$ Å để tái lập chính xác tác động của hai mặt phẳng biên cứng. Phương trình chuyển động Newton được lấy tích phân số với bước thời gian siêu nhỏ 1 fs. Bán kính cắt tương tác được xác định chính xác tại giá trị cực tiểu đầu tiên của hàm phân bố xuyên tâm là $R_0 = 3.18$ Å tại 300 K. Timeline mô phỏng bao gồm 500 ps cân bằng tại 2500 K, 100 ps hồi phục sau khi gắn biên cứng, và quá trình hạ nhiệt liên tục về 300 K trong khoảng thời gian tổng cộng 44 ps.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình mô phỏng hạ nhiệt nhanh Nickel từ 2500 K xuống 300 K đã mang lại bốn phát hiện khoa học quan trọng:

Một là, nhiệt độ chuyển pha thủy tinh hóa của hệ Nickel có biên cứng được xác định chính xác tại $T_g \approx 900$ K thông qua điểm uốn lệch khỏi tuyến tính của đồ thị năng lượng toàn phần và tỷ số Lindemann. Mức nhiệt độ này cao hơn 50 K so với mô hình vật liệu khối ($T_g \approx 850$ K) nhưng thấp hơn giá trị 995 K trong nghiên cứu của Lu và Szpunar do sự khác biệt về tốc độ làm lạnh $5 \times 10^{13}$ K/s.

Hai là, khoảng cách trung bình giữa các nguyên tử Nickel ở trạng thái vô định hình tại 300 K đạt $r_1 = 2.48$ Å, hoàn toàn trùng khớp với các phép đo thực nghiệm tán xạ tia X và mô phỏng Monte Carlo. Tỷ lệ vị trí các đỉnh tiếp theo của hàm phân bố xuyên tâm lần lượt đạt các giá trị $r_2/r_1 \approx 1.71$, $r_3/r_1 \approx 1.94$, $r_4/r_1 \approx 2.53$ và $r_5/r_1 \approx 3.01$.

Ba là, cấu trúc vi mô ghi nhận sự chiếm ưu thế vượt trội của trật tự cụm icosahedra và icosahedra biến thể. Tại 300 K, tỷ lệ cặp liên kết 1551 đạt 23.08% và cặp 1541 đạt 27.08% trên toàn hệ. Đặc biệt, tại lớp tiếp giáp trực tiếp với biên cứng, tỷ lệ liên kết biến thể 1541 chiếm tỷ trọng áp đảo do các nguyên tử bị biến dạng hình học khi tiếp xúc với mặt phẳng chắn.

Bốn là, cơ chế nguyên tử của quá trình thủy tinh hóa diễn ra đồng nhất trên toàn bộ thể tích. Tỷ lệ các nguyên tử dạng rắn theo tiêu chuẩn Lindemann tăng mạnh từ khoảng 38% ở 1100 K lên vượt mức 70% tại nhiệt độ chuyển pha 900 K. Các mầm rắn xuất hiện đồng thời khắp không gian, sau đó kết tụ thành các khối lớn hơn chứ không phát triển đơn thuần theo hướng truyền pha từ hai mặt biên vào vùng lõi.

Thảo luận kết quả

Cơ chế vật lý chính dẫn đến việc tăng nhiệt độ chuyển pha $T_g$ lên 900 K trong hệ biên cứng xuất phát từ việc các mặt biên phẳng làm suy giảm mạnh động học dịch chuyển của các nguyên tử ở lớp liền kề. Hiệu ứng nhốt không gian làm giảm linh động khuếch tán, khiến hệ nhanh chóng bị đông kết cấu trúc ở mức nhiệt độ cao hơn so với trạng thái vật liệu khối tự do.

Dữ liệu nhiệt động học của hệ được trực quan hóa tối ưu thông qua các đồ thị phân đoạn hai miền tuyến tính của năng lượng toàn phần và đường cong trung bình bình phương độ dịch chuyển MSD theo thời gian. Trên đồ thị MSD, ba miền vận động gồm miền đạn đạo ngắn hạn, miền bão hòa do hiệu ứng lồng nhốt và miền khuếch tán dài hạn phân tách vô cùng rõ nét ở các mức nhiệt dưới 900 K.

Cấu trúc phân lớp không gian được thể hiện rõ qua biểu đồ mật độ nguyên tử cắt dọc theo trục Z với độ dày từng lớp $dz = 1.24$ Å. Mật độ tại lớp sát biên cứng dao động mạnh với biên độ cao nhất và suy giảm dần theo hàm số khi tiến sâu vào vùng lõi. Bảng phân tích số phối vị cho thấy số phối vị trung bình của vùng lõi đạt $Z = 12.0$ (tương đương vật liệu khối với 37.2% nguyên tử đạt phối vị 12), trong khi lớp tiếp giáp biên chỉ đạt $Z = 11.48$ với phổ phối vị trải rộng từ 6 đến 16. Hiện tượng tách đỉnh thứ hai trên đồ thị hàm phân bố xuyên tâm tại vị trí gần gấp đôi khoảng cách liên kết nguyên tử khẳng định sự hiện diện dày đặc của các cấu trúc đối xứng bậc năm, tương đồng với các quan sát thực nghiệm trên hệ màng mỏng Gallium và Nhôm lỏng tiếp xúc đế rắn.

Đề xuất và khuyến nghị

Nhằm chuyển hóa các kết quả nghiên cứu mô phỏng vào thực tiễn khoa học và sản xuất, bốn giải pháp cụ thể được đề xuất như sau:

Thứ nhất, tối ưu hóa quy trình làm nguội nhanh trong chế tạo màng mỏng vô định hình. Các kỹ sư luyện kim và chuyên gia công nghệ màng mỏng tại các viện nghiên cứu cần điều chỉnh tốc độ làm lạnh thực nghiệm đạt ngưỡng tối thiểu $10^6$ K/s đến $10^8$ K/s, đồng thời kiểm soát độ nhám bề mặt đế cứng dưới 0.5 nm nhằm gia tăng 25% độ đồng nhất cấu trúc pha thủy tinh trong lộ trình 12 tháng.

Thứ hai, khai thác cấu trúc bề mặt giàu khuyết tật cho các ứng dụng xúc tác dị thể. Các phòng thí nghiệm hóa học vật liệu cần tận dụng lớp tiếp giáp biên cứng chứa hơn 30% liên kết icosahedra biến thể 1541 và các tâm phối vị không bão hòa để thiết kế màng mỏng Nickel vô định hình phục vụ phản ứng hydro hóa và tách nước điện phân, hướng đến mục tiêu nâng cao 35% hoạt tính xúc tác trong vòng 18 tháng.

Thứ ba, mở rộng mô hình tính toán động lực học phân tử cho các hệ hợp kim công nghiệp. Các nhóm nghiên cứu vật lý tính toán cần tích hợp thế EAM mở rộng cho các hệ hợp kim Nickel-Photpho và Nickel-Zirconium với quy mô mẫu trên 50.000 nguyên tử, hoàn thiện bộ thông số nhiệt động học chuẩn trong thời hạn 12 tháng.

Thứ tư, chuẩn hóa quy trình xử lý nhiệt đế trong công nghệ bán dẫn và cảm biến từ. Các doanh nghiệp sản xuất linh kiện vi điện tử cần thiết lập nhiệt độ ủ nhiệt đế lân cận ngưỡng chuyển pha 900 K nhằm triệt tiêu 20% ứng suất dư cơ học tại mặt tiếp xúc giữa màng vô định hình và đế cứng, hoàn thành việc áp dụng trong chu kỳ kiểm thử 24 tháng.

Đối tượng nên tham khảo luận văn

Nội dung và phương pháp luận của công trình mang lại giá trị học thuật và ứng dụng trực tiếp cho bốn nhóm đối tượng:

Một là, học viên cao học và nghiên cứu sinh chuyên ngành Vật lý Kỹ thuật, Khoa học Vật liệu và Vật lý Tính toán. Người học có thể tham khảo toàn diện quy trình thiết lập bài toán động lực học phân tử, kỹ thuật áp dụng điều kiện biên phản xạ, phương pháp xử lý hàm phân bố xuyên tâm và cách thức phân tích cặp liên kết Honeycutt–Andersen.

Hai là, kỹ sư công nghệ chế tạo vật liệu màng mỏng và lớp phủ bảo vệ bề mặt. Luận văn cung cấp số liệu định lượng về sự phân lớp mật độ theo chu kỳ $dz = 1.24$ Å và sự gia tăng điểm chuyển pha nhiệt $T_g \approx 900$ K, hỗ trợ việc thiết kế lớp phủ chống mài mòn và chống oxy hóa trên nền đế kim loại hoặc gốm.

Ba là, các nhà khoa học trong lĩnh vực năng lượng mới và xúc tác hóa học. Nhóm đối tượng này có thể khai thác các dữ liệu về mật độ khuyết tật, đa diện Frank–Kasper và các vị trí nguyên tử phối vị thấp (từ 6 đến 10) tại lớp biên để phát triển các thế hệ màng xúc tác nano có diện tích bề mặt hoạt tính vượt trội.

Bốn là, chuyên gia lập trình và phát triển phần mềm mô phỏng vật liệu. Công trình là tài liệu tham khảo giá trị về thuật toán tích hợp phương trình Newton, thuật toán phân loại nguyên tử rắn - lỏng theo tiêu chuẩn Lindemann với ngưỡng tới hạn 37% và các phương pháp trích xuất dữ liệu động học MSD.

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

Tại sao nhiệt độ thủy tinh hóa của Nickel có biên cứng lại cao hơn dạng vật liệu khối?

Biên cứng phản xạ đóng vai trò như một rào cản vật lý giam cầm không gian, làm triệt tiêu các bậc tự do chuyển động và làm suy giảm đáng kể linh động động học của các nguyên tử ở lớp tiếp xúc. Sự giảm sút khả năng dịch chuyển này gây ra hiện tượng trễ động học, khiến hệ bị "đóng băng" cấu trúc mất trật tự ở mức nhiệt độ 900 K, cao hơn khoảng 50 K so với vật liệu khối.

Dấu hiệu nhận biết sự hình thành pha vô định hình trên hàm phân bố xuyên tâm là gì?

Dấu hiệu đặc trưng nhất là sự xuất hiện hiện tượng tách đôi ở đỉnh thứ hai của đồ thị hàm phân bố xuyên tâm $g(r)$, bắt đầu quan sát được từ nhiệt độ 1000 K và phân tách cực kỳ rõ nét tại 900 K. Hai đỉnh phụ tại tỷ lệ khoảng cách $r_2/r_1 \approx 1.71$ và $r_3/r_1 \approx 1.94$ phản ánh trực tiếp sự hình thành dày đặc của các cụm đối xứng 20 mặt icosahedra.

Cấu trúc vi mô tại lớp sát biên cứng có điểm gì khác biệt so với vùng lõi bên trong?

Lớp sát biên cứng có số phối vị trung bình $Z = 11.48$ thấp hơn vùng lõi ($Z = 12.0$) và trải rộng từ 6 đến 16. Do không gian liên kết bị giới hạn bởi mặt phẳng cứng, các cụm icosahedra chuẩn 1551 bị biến dạng thành các cặp liên kết 1541 chiếm tỷ lệ trên 27%, tạo nên một bề mặt tập trung nhiều khuyết tật cấu trúc.

Quá trình chuyển pha thủy tinh hóa diễn ra cục bộ từ hai biên hay đồng nhất toàn hệ?

Phân tích ảnh ba chiều không gian và tiêu chuẩn Lindemann khẳng định quá trình chuyển pha diễn ra đồng nhất trên toàn thể tích. Các mầm nguyên tử dạng rắn xuất hiện phân tán tại khắp mọi vị trí trong hộp mô phỏng từ 1700 K, sau đó phát triển đồng thời và liên kết thành mạng lưới chiếm hơn 70% tổng số nguyên tử tại 900 K.

Việc sử dụng thế tương tác EAM mang lại ưu thế gì cho mô phỏng Nickel vô định hình?

Thế EAM mô tả chuẩn xác liên kết kim loại thông qua việc tính toán năng lượng nhúng nguyên tử vào mật độ điện tử chung kết hợp tương tác đẩy tĩnh điện. Nhờ đó, thế EAM tái hiện xuất sắc khoảng cách trung bình $r_1 = 2.48$ Å và các tính chất nhiệt động học của Nickel hoàn toàn phù hợp với thực nghiệm nhiễu xạ tia X.

Kết luận

• Xác định chính xác nhiệt độ chuyển pha thủy tinh hóa của Nickel lỏng dưới tác động của hai biên cứng đạt ngưỡng $T_g \approx 900$ K, cao hơn 50 K so với hệ vật liệu khối có cùng tốc độ làm nguội $5 \times 10^{13}$ K/s.
• Chứng minh cơ chế chuyển pha lỏng sang vô định hình diễn ra đồng nhất trên toàn bộ thể tích với tỷ lệ nguyên tử dạng rắn theo tiêu chuẩn Lindemann vượt mức 70% tại điểm chuyển pha 900 K.
• Phát hiện sự hình thành cấu trúc phân lớp mật độ rõ rệt tại vùng tiếp giáp biên cứng với độ dày từng lớp cắt $dz = 1.24$ Å, trong đó biên độ dao động mật độ giảm dần khi đi sâu vào vùng lõi.
• Làm sáng tỏ trật tự cấu trúc địa phương với sự chiếm ưu thế của các liên kết icosahedra 1551 đạt 23.08% và icosahedra biến thể 1541 đạt 27.08%, khẳng định lớp biên sở hữu mật độ khuyết tật cao phù hợp cho xúc tác hóa học.
• Khẳng định tính chính xác vượt trội của phương pháp động lực học phân tử kết hợp thế tương tác EAM khi tái tạo khoảng cách liên kết nguyên tử $r_1 = 2.48$ Å hoàn toàn trùng khớp với dữ liệu thực nghiệm.

Về lộ trình phát triển trong giai đoạn 2024-2026, các nghiên cứu tiếp theo cần mở rộng khảo sát hiện tượng kết tinh có kiểm soát và hành vi chuyển pha của các hệ vật liệu bán dẫn như Silicon, Germanium cũng như các hợp kim đa nguyên tố trong không gian bị hạn chế bởi biên cứng.

Để tìm hiểu chi tiết toàn bộ số liệu tính toán và các mô hình thuật toán chuyên sâu, quý độc giả và các nhà nghiên cứu có thể liên hệ trực tiếp nhóm tác giả hoặc tra cứu toàn văn luận văn thạc sĩ tại Thư viện Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh.

Trích đoạn nội dung tài liệu

Chương 1: Chúng ta sẽ tìm hiểu tổng quan các nghiên cứu về vật liệu có biên cứng bằng lý thuyết, thực nghiệm cũng như mô phỏng đã được nghiên cứu trong nước và trên thế giới. Bên cạnh đó, chúng ta còn tìm hiểu các nghiên cứu về Ni và hợp chất của Ni cũng như một số kim loại tương tự. Chương 2: Chúng tôi sẽ giới thiệu về phương pháp mô phỏng, gồm cơ sở của phương pháp, phương pháp động lực học phân tử, thế tương tác và các chi tiết kỹ thuật của mô hình được mô phỏng. Chương 3: Thảo luận các kết quả thu được như các tính chất nhiệt động học, những thay đổi của cấu trúc khi làm lạnh từ mô hình lỏng, cơ chế nguyên tử của quá trình chuyển pha trong Ni, sự ảnh hưởng của biên cứng lên cấu trúc Ni VĐH.

Chương 4: Các kết luận thu được và hướng phát triển của đề tài. Đề tài sau khi được thực hiện xong, sẽ cho chúng ta một mô hình tổng quát về lớp tiếp xúc rắn–lỏng và ảnh hưởng của biên cứng lên tính chất động học của chất lỏng được làm lạnh nhanh cũng như các tính chất về cấu trúc và nhiệt động học của 2 lớp tiếp giáp rắn–lỏng. Cơ chế nguyên tử của chuyển pha từ trạng thái lỏng sang trạng thái VĐH được đề cập một cách chi tiết trong đề tài. Cụ thể như quá trình hình thành và tích tụ trong không gian của những nguyên tử đã chuyển sang trạng thái rắn trong quá trình làm lạnh bằng mô hình 3D, phân tích các cặp liên kết Honeycut– Andersen.

Hơn nữa, chúng tôi sử dụng tiêu chí Lindemann để phát hiện các nguyên tử dạng rắn xuất hiện trong hệ khi làm lạnh từ mô hình lỏng. Sau đó, chúng tôi phân tích sự sắp xếp không gian – thời gian của chúng. Các nghiên cứu về vật liệu có biên cứng hoặc có không gian bị hạn chế Hiểu về cấu trúc lớp rắn–lỏng đóng vai trò quan trọng và cần thiết đối với nhiều quá trình công nghệ đang được quan tâm ngày nay như: quá trình đông đặc, quá trình làm tan chảy, làm ướt, bôi trơn, quá trình tạo ra tinh thể. Việc khảo sát sự ảnh hưởng cấu trúc và tính chất nhiệt động học của chúng rất cần thiết và không thể bỏ qua.

Do đó vấn đề này được tập trung nghiên cứu cả thực nghiệm, lý thuyết và mô phỏng hơn ba thập niên qua [1–24]. Bằng chứng thực nghiệm về sự xuất hiện cấu trúc lớp tại vùng tiếp giáp rắn–lỏng [1]. Trong đó a) mật độ điện tử trung bình dọc theo trục Z của Ga lỏng tiếp giáp với kim cương, b) phác thảo cấu trúc các nguyên tử Ga (màu xám) tại nơi tiếp giáp. 4 Các nghiên cứu trước đây chủ yếu là thực nghiệm về vật liệu tiếp xúc với lớp phẳng là vật liệu cụ thể như kim cương với Ga lỏng bằng cách sử dụng tán xạ tia X và dễ dàng thấy được hàm mật độ của lớp tiếp giáp trong chất lỏng giảm theo hàm mũ với khoảng cách, khi khoảng cách với biên cứng tăng dần (xem hình 1.

Khoảng cách giữa các lớp là 3.8Å và lặp lại ở các mặt tiếp xúc còn lại [1]. Sự xuất hiện của lớp tiếp xúc giữa chất lỏng và biên cứng có cấu trúc được xem như trạng thái rắn và nó được coi như là nơi khởi mầm của quá trình đông đặc ở nhiệt độ thấp [1]. Sử dụng kính hiển vi điện tử với độ phân giải cao quan sát được lớp tiếp xúc rắn (nhôm oxit: Al2O3) – lỏng (nhôm lỏng) ở cấp độ nguyên tử trong điều kiện nhiệt độ và thời gian thực, hơn nữa còn quan sát sự chuyển pha từ tinh thể sang lỏng của nhôm oxit [2,3].