Tổng quan nghiên cứu

Vật liệu nano từ tính với kích thước hạt dưới 100 nm đang trở thành trọng tâm nghiên cứu của khoa học vật liệu hiện đại nhờ các tính chất nhiệt, quang, điện và từ vượt trội so với vật liệu khối truyền thống. Trong nhóm hợp kim kim loại chuyển tiếp và á kim, hợp kim sắt - bo (FeB) đóng vai trò then chốt trong chế tạo lõi biến áp điện, thiết bị lưu trữ thông tin và chất xúc tác hiệu năng cao. Tuy nhiên, sự hiểu biết về cơ chế hình thành pha tinh thể và pha thủy tinh ở cấp độ nguyên tử trong các hạt nano kim loại vẫn còn nhiều hạn chế do diện tích bề mặt lớn và cấu trúc vỏ bất đối xứng gây khó khăn cho các phương pháp thực nghiệm truyền thống.

Nghiên cứu được thực hiện tại Trường Đại học Sư phạm thuộc Đại học Thái Nguyên vào năm 2017 nhằm giải quyết bài toán trên thông qua việc phân tích toàn diện cấu trúc vi mô của hạt nano Fe100-xBx (với x nhận các giá trị 0, 5, 10). Mục tiêu trọng tâm của đề tài là khảo sát sự tiến hóa cấu trúc của hạt nano chứa 5.000 nguyên tử trong dải nhiệt độ rộng từ 200K đến 1500K ở áp suất phòng 0 GPa, đồng thời làm sáng tỏ cơ chế tinh thể hóa và động học phục hồi cấu trúc khi ủ nhiệt đẳng nhiệt tại nhiệt độ 900K.

Ý nghĩa khoa học của nghiên cứu được thể hiện qua các chỉ số định lượng cụ thể về sự phân vùng cấu trúc vi mô. Kết quả chỉ ra rằng mật độ nguyên tử ở vùng lõi hạt nano đạt mức ổn định 0,0825 nguyên tử trên mỗi angstrom khối, trong khi mật độ ở lớp vỏ bề mặt có độ dày 6,0 Å chỉ đạt 0,0304 nguyên tử trên mỗi angstrom khối. Những thông số này cung cấp nền tảng vật lý vững chắc để định hướng công nghệ xử lý nhiệt luyện và thiết kế các loại vật liệu từ mềm chất lượng cao trong công nghiệp.

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 xây dựng trên hai nền tảng lý thuyết cốt lõi của vật lý chất rắn và vật lý thống kê:

Thứ nhất là Lý thuyết mầm và sự phát triển mầm cổ điển (Classical Nucleation Theory). Lý thuyết này giải thích sự biến đổi pha từ trạng thái vô định hình sang tinh thể dựa trên sự cạnh tranh giữa năng lượng tự do thể tích và năng lượng tự do bề mặt của các đám nguyên tử. Khi độ siêu bão hòa tăng lên, bán kính mầm tới hạn giảm xuống, tạo điều kiện thuận lợi cho mầm tinh thể tự phát triển vượt qua rào cản năng lượng.

Thứ hai là Cơ học thống kê và nhiệt động lực học của vật chất ngưng tụ phi tinh thể. Khung lý thuyết này cho phép mô tả trạng thái cân bằng và cấu trúc trật tự tầm gần thông qua các hàm phân bố mật độ xác suất không gian.

Các khái niệm cơ bản được áp dụng xuyên suốt bao gồm:

  • Hàm phân bố xuyên tâm (Radial Distribution Function - RDF): Đại lượng biểu diễn xác suất tìm thấy nguyên tử ở khoảng cách r so với một nguyên tử gốc, dùng để nhận diện trật tự tầm gần và tầm xa.
  • Số phối trí (Coordination Number - SPT): Số lượng nguyên tử lân cận gần nhất bao quanh một nguyên tử trung tâm trong bán kính ngắt xác định.
  • Mầm tinh thể cơ bản: Cấu hình gồm một nguyên tử trung tâm có 14 lân cận gần nhất thỏa mãn cấu trúc mạng lập phương tâm khối (BCC).
  • Trạng thái vô định hình và thủy tinh: Trạng thái rắn không có trật tự tầm xa nhưng duy trì trật tự tầm gần, hình thành khi làm nguội nhanh từ pha lỏng.

Mô hình nghiên cứu sử dụng thế tương tác cặp bán thực nghiệm Pak-Doyama. Thế này mô tả chính xác tương tác giữa các cặp nguyên tử Fe-Fe với bán kính ngắt 3,44 Å, Fe-B với bán kính ngắt 3,09 Å và B-B với bán kính ngắt 3,79 Å.

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

Nguồn dữ liệu của luận văn được thu thập từ mô hình tính toán số hạt nano gồm 5.000 nguyên tử được phân bố ban đầu trong không gian hình cầu bán kính 27,5 Å đến 28,0 Å với điều kiện biên tự do.

Phương pháp nghiên cứu kết hợp hai kỹ thuật mô phỏng tính toán chủ đạo:

  • Phương pháp Thống kê hồi phục (TKHP): Áp dụng ở nhiệt độ 0K nhằm tối thiểu hóa thế năng của hệ. Từng nguyên tử được dịch chuyển theo hướng lực tổng hợp với bước dịch chuyển ban đầu 0,4 Å trong 150 bước đầu, sau đó giảm dần xuống 0,02 Å cho đến khi hệ đạt trạng thái cân bằng năng lượng sau khoảng 100.000 bước tính toán.
  • Phương pháp Động lực học phân tử (ĐLHPT): Sử dụng thuật toán Verlet để giải hệ phương trình vi phân chuyển động Newton theo bước thời gian tích phân. Phương pháp này theo dõi chính xác quỹ đạo, vận tốc và năng lượng của từng nguyên tử theo thời gian ủ nhiệt đẳng nhiệt từ 20.000.000 đến 70.000.000 bước tại 900K trong hệ tập hợp vi chính tắc (NVE) và tập hợp đẳng nhiệt - đẳng tích (NVT).

Phương pháp mô phỏng số được lựa chọn vì nó cho phép quan sát chi tiết sự dịch chuyển của từng nguyên tử và các mầm kết tinh ở kích thước vài angstrom - điều mà các thiết bị phân tích thực nghiệm như nhiễu xạ tia X hay kính hiển vi điện tử truyền qua rất khó ghi nhận liên tục theo thời gian thực. Toàn bộ thuật toán thống kê số phối trí và nhận dạng mầm tinh thể BCC được lập trình và xử lý tự động bằng ngôn ngữ C++.

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 cấu trúc và cơ chế chuyển pha của hạt nano Fe100-xBx đã mang lại các phát hiện khoa học quan trọng:

  • Sự chuyển pha từ vô định hình sang tinh thể lập phương tâm khối (BCC): Khi ủ nhiệt tại 900K, hạt nano sắt và Fe95B5 sau 100.000 bước ĐLHPT vẫn giữ cấu trúc vô định hình. Tuy nhiên, khi kéo dài thời gian ủ nhiệt lên 70.000.000 bước ĐLHPT, mạng tinh thể BCC đã hình thành hoàn chỉnh trong vùng lõi của hạt nano.
  • Đặc trưng tiến hóa của hàm phân bố xuyên tâm: Đối với mẫu Fe thuần ở 900K sau 100.000 bước, đỉnh thứ nhất của hàm phân bố xuyên tâm xuất hiện ở vị trí 2,55 Å với độ cao 3,60. Đỉnh thứ hai tách thành hai đỉnh phụ ở vị trí 4,25 Å (độ cao 1,63) và 5,10 Å (độ cao 1,23). Sau 70.000.000 bước, đỉnh phụ bên phải tăng cao vượt trội so với đỉnh bên trái và xuất hiện các cực đại mới ở khoảng cách lớn trên 10 Å, khẳng định sự thiết lập của trật tự tầm xa dạng tinh thể.
  • Sự biến thiên năng lượng qua 3 giai đoạn rõ rệt: Đồ thị năng lượng của hệ hạt nano Fe theo thời gian ủ nhiệt tại 900K cho thấy năng lượng sụt giảm mạnh khoảng 0,032 eV trên mỗi nguyên tử trong giai đoạn chuyển pha cấu trúc chính. Năng lượng trung bình của mẫu Fe95B5 ổn định quanh mức -1,2435 eV đến -1,2573 eV, còn mẫu Fe90B10 dao động quanh mức -1,2363 eV đến -1,2585 eV.
  • Sự dịch chuyển số phối trí trung bình: Số phối trí của hệ tăng từ 13 (đặc trưng cho trạng thái vô định hình) lên 14 (đặc trưng chuẩn của cấu trúc mạng BCC), đồng thời khoảng cách giữa các nguyên tử lân cận có xu hướng giãn nhẹ khi tăng nhiệt độ.

Thảo luận kết quả

Cơ chế tinh thể hóa của hạt nano sắt và FeB bắt nguồn từ sự hình thành và phát triển của các mầm tinh thể cơ bản trong vùng lõi. Một mầm tinh thể BCC chuẩn được xác định khi nguyên tử trung tâm có đủ 14 nguyên tử lân cận trong bán kính ngắt 3,35 Å, trong đó 6 nguyên tử có 4 lân cận chung và 8 nguyên tử có 6 lân cận chung. Sự sụt giảm 0,032 eV năng lượng hệ chính là động lực nhiệt động thúc đẩy sự định hướng lại liên kết giữa các nguyên tử sắt.

Khi so sánh với mẫu sắt khối vô định hình có độ cao đỉnh thứ nhất là 3,51 tại vị trí 2,55 Å, hạt nano sắt có độ cao đỉnh thứ nhất đạt 3,60, phản ánh mức độ trật tự cục bộ cao hơn ở vùng lõi. Việc bổ sung 5% đến 10% nguyên tử Bo làm giảm độ cao đỉnh thứ nhất xuống mức 3,21 tại vị trí 2,51 Å. Nguyên nhân là do bán kính nguyên tử Bo nhỏ hơn sắt, khi xen kẽ vào mạng tinh thể sẽ làm biến dạng trường thế cục bộ và làm chậm tốc độ phát triển của mầm tinh thể.

Dữ liệu mô phỏng được biểu diễn trực quan và sinh động thông qua hệ thống đồ thị hàm phân bố xuyên tâm g(r), biểu đồ phân bố số phối trí theo nhiệt độ, đường cong suy giảm năng lượng tự do và các ảnh chụp không gian 3 chiều phân tách rõ rệt giữa nguyên tử vùng lõi và lớp vỏ xốp bề mặt.

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

Dựa trên các kết quả mô phỏng vi cấu trúc và động học chuyển pha, luận văn đưa ra 4 khuyến nghị ứng dụng thực tiễn:

  • Tối ưu hóa quy trình nhiệt luyện công nghiệp cho hợp kim từ mềm: Các cơ sở sản xuất và viện nghiên cứu vật liệu nên thiết lập chế độ ủ nhiệt hạt nano FeB trong dải nhiệt độ 850K đến 950K với thời gian ủ tương đương trên 50.000.000 chu kỳ dao động mạng. Giải pháp này giúp kiểm soát tỷ lệ kết tinh mạng BCC đạt trên 85%, tối ưu hóa độ từ thẩm và giảm tổn hao trễ từ cho lõi biến áp điện cao tần trong vòng 12 đến 18 tháng tới.
  • Mở rộng quy mô mô hình mô phỏng lên kích thước lớn: Nhóm nghiên cứu tính toán vật lý chất rắn cần mở rộng mô hình lên quy mô từ 500.000 đến 1.000.000 nguyên tử bằng các thuật toán tính toán song song. Mục tiêu là khảo sát chi tiết hơn sự ảnh hưởng của độ dày lớp vỏ xốp (hiện chiếm mật độ 0,0304 nguyên tử trên mỗi angstrom khối) đối với các hạt nano có đường kính trên 20 nm trong giai đoạn 6 đến 12 tháng.
  • Chuẩn hóa bộ thế tương tác nhiều hạt cho hợp kim sắt - á kim: Các phòng thí nghiệm mô phỏng tính toán cần nâng cấp từ thế tương tác cặp Pak-Doyama sang các thế tương tác nhúng nhiều hạt (Embedded Atom Method - EAM). Việc này nhằm nâng cao độ chính xác dự báo năng lượng tự do bề mặt lên trên 95% trong lộ trình 24 tháng.
  • Kiểm soát nghiêm ngặt nồng độ pha tạp Bo ở ngưỡng 5%: Các đơn vị luyện kim và chế tạo vật liệu nano nên duy trì tỷ lệ Bo trong hợp kim ở mức 5% để vừa giữ được khả năng tạo thủy tinh khi làm nguội nhanh, vừa đảm bảo tốc độ tạo mầm tinh thể đồng đều khi gia nhiệt xử lý.

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

Tài liệu luận văn mang lại giá trị học thuật và ứng dụng thực tiễn cao cho 4 nhóm đối tượng chuyên môn:

  • Các nhà nghiên cứu khoa học vật liệu và vật lý chất rắn: Nắm bắt bức tranh chi tiết về cơ chế chuyển pha vô định hình - tinh thể ở cấp độ nguyên tử, sử dụng các số liệu hàm phân bố xuyên tâm và số phối trí làm dữ liệu đối chuẩn cho các nghiên cứu thực nghiệm tán xạ tia X và neutron.
  • Kỹ sư công nghệ vật liệu từ tính và luyện kim bột: Khai thác thông số nhiệt độ ủ 900K và tỷ lệ hợp kim Fe95B5 để thiết kế dây chuyền công nghệ nhiệt luyện bột nano từ mềm dùng trong chế tạo linh kiện điện tử viễn thông và xe điện.
  • Giảng viên, học viên cao học và sinh viên ngành Vật lý: Sử dụng công trình như một tài liệu tham khảo chuẩn mực về phương pháp mô phỏng động lực học phân tử, kỹ thuật lập trình nhận diện mầm tinh thể bằng C++ và phương pháp xử lý dữ liệu cấu trúc vi mô.
  • Chuyên gia phát triển phần mềm mô phỏng vật liệu: Ứng dụng quy trình tích hợp giữa phương pháp Thống kê hồi phục ở 0K và mô hình Động lực học phân tử Verlet để xây dựng các gói phần mềm mô phỏng vật liệu nano kim loại mới.

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

  • Phương pháp mô phỏng động lực học phân tử mang lại ưu thế gì vượt trội so với thực nghiệm? Phương pháp cho phép theo dõi chính xác vị trí và vận tốc của từng nguyên tử trong 5.000 hạt theo từng bước thời gian siêu nhỏ. Điều này giúp quan sát trực tiếp sự hình thành mầm tinh thể BCC ở kích thước 3,35 Å mà các thiết bị thực nghiệm hiện nay chưa thể ghi nhận liên tục theo thời gian thực.

  • Vì sao hàm phân bố xuyên tâm của vật liệu vô định hình lại có hiện tượng tách đôi ở cực đại thứ hai? Hiện tượng tách đỉnh thứ hai thành hai đỉnh phụ ở 4,25 Å và 5,10 Å là dấu hiệu đặc trưng của trạng thái vô định hình đặc sít. Nó phản ánh sự sắp xếp bất đối xứng của các đa diện phối trí tầm ngắn trước khi hệ chuyển sang trật tự tầm xa của mạng tinh thể.

  • Việc pha tạp nguyên tử Bo từ 5% đến 10% ảnh hưởng như thế nào đến hạt nano sắt? Nguyên tử Bo xen kẽ vào mạng làm giảm độ cao đỉnh thứ nhất của hàm phân bố xuyên tâm từ 3,60 xuống 3,21. Nồng độ Bo càng cao làm tăng độ ổn định của pha vô định hình và cản trở tốc độ phát triển mầm tinh thể trong quá trình ủ nhiệt.

  • Tiêu chuẩn nào được sử dụng để xác định một mầm tinh thể cơ bản trong mô hình? Một mầm tinh thể BCC được xác định khi nguyên tử trung tâm có đúng 14 nguyên tử lân cận trong bán kính ngắt 3,35 Å. Trong số 14 lân cận này, bắt buộc phải có 6 nguyên tử có 4 lân cận chung và 8 nguyên tử có 6 lân cận chung.

  • Sự khác biệt về mật độ giữa vùng lõi và vùng bề mặt hạt nano là bao nhiêu? Vùng lõi hạt nano có cấu trúc đặc sít với mật độ đạt 0,0825 nguyên tử trên mỗi angstrom khối, tương đương vật liệu khối. Ngược lại, lớp vỏ bề mặt có cấu trúc xốp với mật độ chỉ đạt 0,0304 nguyên tử trên mỗi angstrom khối, thấp hơn vùng lõi khoảng 63,15%.

Kết luận

  • Luận văn đã xây dựng thành công mô hình mô phỏng hạt nano Fe, Fe95B5 và Fe90B10 gồm 5.000 nguyên tử bằng sự kết hợp giữa phương pháp Thống kê hồi phục và Động lực học phân tử với thế Pak-Doyama.
  • Làm sáng tỏ cấu trúc phân vùng của hạt nano hình cầu bán kính 27,5 Å, chỉ ra sự chênh lệch mật độ rõ rệt giữa vùng lõi (0,0825 nguyên tử/Å3) và lớp vỏ xốp bề mặt (0,0304 nguyên tử/Å3).
  • Xác định chính xác động học chuyển pha vô định hình sang mạng tinh thể lập phương tâm khối (BCC) ở nhiệt độ 900K qua 70.000.000 bước mô phỏng, đi kèm mức giảm năng lượng 0,032 eV trên mỗi nguyên tử.
  • Thiết lập tiêu chuẩn và thuật toán nhận dạng mầm tinh thể BCC với số phối trí 14, giải thích định lượng vai trò cản trở kết tinh của nguyên tử Bo ở nồng độ 5% và 10%.
  • Cung cấp cơ sở lý thuyết và dữ liệu số quan trọng để định hướng công nghệ xử lý nhiệt luyện hạt nano từ mềm trong các ứng dụng biến áp điện và lưu trữ từ tính.

Trong giai đoạn tiếp theo từ 6 đến 12 tháng, các nghiên cứu mở rộng nên tập trung vào việc nâng cấp số lượng nguyên tử lên trên 100.000 hạt và áp dụng thế tương tác nhúng nhiều hạt để hoàn thiện bức tranh chuyển pha. Quý độc giả và các nhà nghiên cứu quan tâm có thể khai thác các thông số cấu trúc vi mô trong luận văn để ứng dụng vào công tác nghiên cứu và giảng dạy chuyên sâu.