Tổng quan nghiên cứu

Vật liệu kim loại chuyển tiếp đóng vai trò then chốt trong khoa học vật liệu hiện đại, trong đó Niken nguyên chất sở hữu mật độ thực 7,72 g/cm³ và được ứng dụng rộng rãi từ các hợp kim chịu nhiệt đến linh kiện điện tử nano. Tuy nhiên, việc nắm bắt chi tiết vi cấu trúc và cơ chế khuếch tán nguyên tử ở cấp độ vi mô vẫn là một thách thức lớn khi các kỹ thuật thực nghiệm như nhiễu xạ tia X hay hiển vi điện tử truyền qua gặp nhiều rào cản đo lường ở trạng thái nhiệt độ cực hạn. Vấn đề cốt lõi đặt ra là làm sáng tỏ bản chất chuyển pha và động học khuếch tán giữa trạng thái kim loại lỏng, vùng quá lạnh và trạng thái rắn vô định hình.

Luận văn thạc sĩ chuyên ngành Vật lý chất rắn với đề tài nghiên cứu về 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 được thực hiện tại Trường Đại học Sư phạm thuộc Đại học Thái Nguyên. Mục tiêu cụ thể của công trình là xây dựng thành công hai mô hình mẫu Niken chứa 4.394 nguyên tử và 10.000 nguyên tử, sau đó tiến hành ủ nhiệt và khảo sát toàn diện trong dải nhiệt độ rộng từ 1.380 K xuống đến 118 K. Nghiên cứu tập trung giải mã các biến thiên cấu trúc thông qua hàm phân bố xuyên tâm, phân bố số phối trí, thống kê đơn vị cấu trúc và tính toán dịch chuyển bình phương trung bình. Kết quả nghiên cứu mang ý nghĩa học thuật và ứng dụng sâu sắc, cung cấp dữ liệu định lượng chính xác với độ tin cậy vượt trội, tái hiện chính xác tỷ lệ các đỉnh cực đại cấu trúc với sai số so với thực nghiệm dưới 4%, mở ra nền tảng lý thuyết vững chắc cho công nghệ luyện kim và chế tạo vật liệu vô định hình tiên tiến.

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 được thiết lập vững chắc trên nền tảng cơ học thống kê hiện đại và động lực học cổ điển Newton. Khung lý thuyết trọng tâm xoay quanh lý thuyết trật tự tầm gần trong vật liệu vô định hình và lý thuyết mầm kết tinh cổ điển. Trong trạng thái vô định hình, tính tuần hoàn không gian bị phá vỡ hoàn toàn, biến các nguyên tử thành các tâm tán xạ ngẫu nhiên nhưng vẫn duy trì cấu trúc sắp xếp trật tự tầm gần đặc trưng bởi các khối đa diện 20 mặt. Ba khái niệm nền tảng được sử dụng xuyên suốt bao gồm hàm phân bố xuyên tâm g(r) biểu thị xác suất tìm thấy nguyên tử ở khoảng cách r; thừa số cấu trúc S(k) mô tả cường độ tán xạ đàn hồi qua biến đổi Fourier không gian ba chiều; và chỉ số cặp liên kết Honeycutt-Andersen phân loại chi tiết các liên kết cục bộ như cấu hình 1551 đại diện cho đối xứng ngũ giác và 1421 đại diện cho trật tự tinh thể.

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

Nguồn dữ liệu của luận văn được kiến tạo thông qua sự kết hợp giữa phương pháp động lực học phân tử cổ điển và phương pháp thống kê hồi phục ở 0 K. Cỡ mẫu nghiên cứu được thiết lập chuẩn hóa ở hai mức quy mô gồm 4.394 nguyên tử và 10.000 nguyên tử, được gieo ngẫu nhiên trong hộp mô phỏng lập phương có áp dụng điều kiện biên tuần hoàn nhằm triệt tiêu hoàn toàn hiệu ứng kích thước bề mặt. Lý do lựa chọn phương pháp động lực học phân tử kết hợp thuật toán tích phân Verlet với bước thời gian mô phỏng 0,001 ps là khả năng giải chính xác hệ phương trình vi phân chuyển động đa hạt, theo dõi trực tiếp quỹ đạo không gian và vận tốc của từng nguyên tử.

Tương tác giữa các nguyên tử Niken được mô tả thông qua thế tương tác cặp Pak-Doyama với bán kính cắt r_c = 3,44 Å cùng hệ số thế tương tác a = -0,12929 eV/Å⁴, b = -1,82709 Å, c = 1,16473 eV/Ų, d = -2,50849 Å và e = -0,13705 eV. Quy trình tính toán được thực hiện theo tiến trình thời gian chặt chẽ: ban đầu mẫu được khử trạng thái nhớ bằng cách ủ nhiệt tại 3.000 K trong hơn 5.000.000 bước mô phỏng, sau đó làm lạnh tuần tự qua 10 mốc nhiệt độ từ 1.400 K xuống 200 K và khảo sát chi tiết từ 1.380 K đến 118 K trong các hệ điều hòa vi chính tắc NVE và đẳng nhiệt đẳng tích NVT. Phương pháp thống kê hồi phục sử dụng bước dịch chuyển ban đầu 0,4 Å và giảm dần về 0,02 Å giúp mẫu đạt trạng thái cân bằng thế năng ổn định tuyệt đối.

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

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

Thứ nhất, nghiên cứu phát hiện sự biến đổi rõ nét của hàm phân bố xuyên tâm g(r) khi nhiệt độ giảm từ 1.380 K xuống 118 K. Tại trạng thái vô định hình 300 K và 118 K, đỉnh thứ hai của hàm g(r) bị tách đôi thành hai đỉnh phụ riêng biệt với đỉnh bên trái luôn cao hơn đỉnh bên phải, khẳng định sự hình thành của cấu trúc đa diện 20 mặt. Tỷ lệ vị trí các đỉnh phụ so với đỉnh thứ nhất đạt các giá trị chuẩn xác gồm r2/r1 = 1,71, r3/r1 = 1,93, r4/r1 = 2,55 và r5/r1 = 3,38, tương đồng trên 96% so với dữ liệu thực nghiệm nhiễu xạ tia X của Niken vô định hình. Ngược lại, ở pha lỏng từ 900 K đến 1.380 K, đỉnh thứ hai hoàn toàn trơn nhẵn và không xuất hiện sự phân tách.

Thứ hai, phân bố số phối trí và cấu trúc địa phương thể hiện sự chuyển dịch mạnh mẽ theo mức độ hồi phục. Trong pha lỏng ở 1.000 K, số phối trí trải rộng từ 10 đến 16 với giá trị trung bình tập trung quanh mức 12. Khi chuyển dần sang pha tinh thể ở 873 K và pha vô định hình dưới 642 K, số nguyên tử có số phối trí 14 tăng vọt, đồng thời các cặp liên kết 1421 tăng trưởng mạnh mẽ trong khi cặp 1551 giảm dần trong suốt tiến trình hồi phục từ 1.000.000 bước đến 7.000.000 bước.

Thứ ba, đồ thị dịch chuyển bình phương trung bình theo thời gian đã vạch rõ ranh giới động học khuếch tán. Ở pha lỏng tại 995 K, 1.167 K và 1.380 K, dịch chuyển bình phương trung bình tăng tuyến tính theo thời gian mô phỏng, minh chứng cho cơ chế khuếch tán liên tục của các nguyên tử. Ở các mức nhiệt độ thấp 118 K, 336 K, 516 K và 642 K thuộc pha vô định hình, đường cong dịch chuyển bình phương trung bình xuất hiện vùng nằm ngang kéo dài, phản ánh hiện tượng đóng băng khuếch tán và nguyên tử bị giam hãm trong các lồng cấu trúc địa phương.

Thảo luận kết quả

Nguyên nhân sâu xa của hiện tượng tách đỉnh thứ hai trong hàm g(r) bắt nguồn từ sự hình thành trật tự tầm gần dạng khối đa diện 20 mặt khi động năng nguyên tử suy giảm, ngăn cản quá trình tái sắp xếp mạng tinh thể chu kỳ dài. Kết quả mô phỏng phù hợp cao với các nghiên cứu động lực học phân tử của các nhà khoa học trên thế giới, xác nhận nhiệt độ chuyển pha thủy tinh của Niken nằm trong khoảng 670 K đến 800 K tùy thuộc vào tốc độ làm nguội từ 1,9×10¹² K/s đến 5,7×10¹² K/s. Dữ liệu nghiên cứu được biểu diễn trực quan và sinh động thông qua các đồ thị hàm phân bố xuyên tâm hai chiều, biểu đồ cột biểu diễn phổ phân bố số phối trí, và đặc biệt là hệ thống hình ảnh mô phỏng không gian ba chiều hiển thị rõ nét quá trình phôi tinh thể liên kết thành các đám mầm lớn trước khi chuyển hóa cấu trúc hoàn toàn.

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

Thứ nhất, xây dựng và chuẩn hóa quy trình tính toán động lực học phân tử cho các hệ kim loại chuyển tiếp và hợp kim vô định hình đa cấu tử. Mục tiêu cụ thể là nâng cao độ chính xác dự báo các thông số nhiệt động học lên 15% trong vòng 12 tháng, do các nhóm nghiên cứu vật lý tính toán tại các trường đại học và viện nghiên cứu chuyên ngành chủ trì thực hiện.

Thứ hai, tích hợp các mô hình thế tương tác nhiều hạt như thế nhúng EAM kết hợp thuật toán học máy vào mô phỏng vật liệu nano. Định hướng mục tiêu là mở rộng quy mô mẫu mô phỏng lên 100.000 nguyên tử và rút ngắn thời gian tính toán 30% trước quý 4 năm 2027, do các chuyên gia công nghệ vật liệu và lập trình mô phỏng phối hợp triển khai.

Thứ ba, ứng dụng cơ chế kiểm soát tốc độ làm lạnh từ 10¹² K/s đến 10¹³ K/s vào thực tiễn luyện kim bột và công nghệ in 3D kim loại tiên tiến. Giải pháp này hướng tới mục tiêu tăng cường độ cứng cơ học và độ bền chống ăn mòn của hợp kim Niken vô định hình lên ít nhất 20% trong giai đoạn 2027 - 2028, do các doanh nghiệp chế tạo vật liệu chịu nhiệt liên kết thử nghiệm.

Thứ tư, mở rộng khảo sát cơ chế khuếch tán nguyên tử tại các bề mặt tiếp xúc và ranh giới hạt dưới tác động của gradient nhiệt độ và áp suất cao. Đề xuất hoàn thiện cơ sở dữ liệu hệ số khuếch tán cho 5 hệ vật liệu kim loại mới trong thời hạn 24 tháng tới, do các phòng thí nghiệm trọng điểm quốc gia về vật lý chất rắn đảm nhiệm.

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

Thứ nhất, các nhà nghiên cứu và nghiên cứu sinh chuyên ngành Vật lý chất rắn, Khoa học vật liệu và Vật lý tính toán. Luận văn cung cấp phương pháp luận chuẩn mực về mô phỏng động lực học phân tử, thuật toán Verlet, cách thiết lập thế tương tác Pak-Doyama và phương pháp thống kê hồi phục ở 0 K.

Thứ hai, kỹ sư công nghệ luyện kim, vật liệu nano và chế tạo máy cao cấp. Nhóm đối tượng này có thể khai thác các quy luật biến thiên cấu trúc theo nhiệt độ và các đặc trưng khuếch tán để tối ưu hóa công nghệ làm nguội nhanh và kiểm soát sự hình thành pha tinh thể trong sản xuất thực tế.

Thứ ba, các lập trình viên và chuyên gia phát triển phần mềm mô phỏng động học phân tử. Tài liệu mang lại bộ thuật toán chi tiết về điều khiển nhiệt độ theo mô hình NVT, hiệu chỉnh áp suất NPT và kỹ thuật trực quan hóa dữ liệu không gian ba chiều cho hệ hàng chục nghìn nguyên tử.

Thứ tư, giảng viên đại học và học viên cao học ngành Vật lý học. Luận văn là tài liệu tham khảo giảng dạy mẫu mực về cấu trúc kim loại lỏng, vật liệu vô định hình, hàm phân bố xuyên tâm và cơ sở lý thuyết vi mô của vật lý chất ngưng tụ.

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

Hàm phân bố xuyên tâm g(r) phản ánh điều gì trong vật liệu Niken vô định hình? Hàm phân bố xuyên tâm g(r) định lượng xác suất tìm thấy nguyên tử ở khoảng cách r quanh nguyên tử trung tâm. Điểm đặc trưng lớn nhất ở trạng thái vô định hình tại 300 K là đỉnh thứ hai tách thành hai đỉnh phụ với tỷ lệ r2/r1 đạt 1,71, chứng minh sự tồn tại của trật tự tầm gần dạng đa diện 20 mặt.

Thế tương tác Pak-Doyama có ưu điểm gì trong nghiên cứu kim loại Niken? Thế cặp Pak-Doyama được xây dựng dựa trên dữ liệu thực nghiệm chính xác với bán kính cắt 3,44 Å. Thế năng này mô tả chuẩn xác lực hút và lực đẩy giữa các cặp nguyên tử Niken, giúp tái hiện cấu trúc vi mô với chi phí tính toán tối ưu hơn so với các thế đa hạt phức tạp.

Phương pháp thống kê hồi phục đóng vai trò gì trong mô hình mô phỏng? Phương pháp thống kê hồi phục dịch chuyển từng nguyên tử theo chiều hợp lực với bước dịch chuyển giảm dần từ 0,4 Å về 0,02 Å. Phương pháp này giúp hệ triệt tiêu động năng về 0 K, tối thiểu hóa thế năng và đưa mẫu về trạng thái cân bằng cấu trúc bền vững nhất.

Cơ chế khuếch tán của Niken biến đổi ra sao giữa trạng thái lỏng và rắn? Ở trạng thái lỏng từ 995 K đến 1.380 K, các nguyên tử khuếch tán tự do liên tục thể hiện qua đường dịch chuyển bình phương trung bình dốc tuyến tính. Khi chuyển sang trạng thái rắn vô định hình dưới 642 K, đồ thị xuất hiện các thềm phẳng nằm ngang, báo hiệu sự đóng băng khuếch tán.

Làm thế nào để xác định chính xác nhiệt độ chuyển pha thủy tinh của Niken? Nhiệt độ chuyển pha thủy tinh được xác định thông qua sự thay đổi độ dốc đột ngột của đường phụ thuộc nội năng theo nhiệt độ, sự xuất hiện của vùng nằm ngang thứ hai trên đồ thị dịch chuyển bình phương trung bình và sự phân tách đỉnh thứ hai trên đường cong hàm phân bố xuyên tâm.

Kết luận

  • Luận văn đã xây dựng thành công mô hình mô phỏng Niken quy mô 4.394 và 10.000 nguyên tử bằng phương pháp động lực học phân tử kết hợp thế Pak-Doyama trong dải nhiệt độ từ 1.380 K đến 118 K.
  • Làm sáng tỏ đặc trưng cấu trúc vô định hình với hiện tượng tách đỉnh thứ hai của hàm phân bố xuyên tâm và xác lập các tỷ lệ đỉnh thực nghiệm chính xác gồm r2/r1 = 1,71 và r3/r1 = 1,93.
  • Phân tích chi tiết sự tiến hóa của số phối trí và các cặp liên kết Honeycutt-Andersen 1421 và 1551 trong suốt quá trình hồi phục nhiệt động từ 1.000.000 đến 7.000.000 bước.
  • Chứng minh rõ nét cơ chế chuyển đổi từ khuếch tán tự do ở pha lỏng sang trạng thái đóng băng chuyển động nguyên tử ở pha vô định hình thông qua đại lượng dịch chuyển bình phương trung bình.
  • Khẳng định lộ trình mở rộng nghiên cứu sang các hệ hợp kim đa nguyên tố trong giai đoạn 2026 - 2028, đóng góp nguồn dữ liệu lý thuyết nền tảng phục vụ đắc lực cho công nghệ chế tạo vật liệu tiên tiến.