Tổng quan nghiên cứu

Hạt nano kim loại sắt (Fe) sở hữu diện tích bề mặt riêng rất lớn cùng cấu trúc hình thái đặc thù, mang lại nhiều tính chất vật lý vượt trội so với vật liệu sắt khối truyền thống có khối lượng riêng 7,81 g/cm3. Tùy thuộc vào phương pháp chế tạo, hạt nano Fe có thể tồn tại ở trạng thái vô định hình hoặc trạng thái tinh thể. Dù có tiềm năng ứng dụng to lớn trong lưu trữ từ tính, xúc tác hóa học và y sinh, cơ chế chuyển pha ở cấp độ nguyên tử cũng như ảnh hưởng của kích thước hạt trong quá trình tinh thể hóa vẫn là bài toán học thuật còn nhiều điểm chưa sáng tỏ do giới hạn quan sát trực tiếp của các phương pháp thực nghiệm truyền thống.

Luận văn thạc sĩ chuyên ngành Vật lý chất rắn (mã số 8440104) của học viên Khamphone Yiachu, được thực hiện tại Trường Đại học Sư phạm thuộc Đại học Thái Nguyên dưới sự hướng dẫn khoa học của TS. Phạm Mai An và PGS. Phạm Hữu Kiên, tập trung giải quyết trọn vẹn vấn đề này. Mục tiêu trọng tâm của đề tài là khảo sát chi tiết đặc trưng vi cấu trúc và hiệu ứng kích thước của hạt nano Fe thông qua phương pháp thực nghiệm số, làm rõ cấu trúc phân tầng giữa phần lõi và phần vỏ, đồng thời xác lập cơ chế hình thành mầm tinh thể lập phương tâm khối (BCC).

Phạm vi nghiên cứu bao gồm việc xây dựng và khảo sát 3 mô hình hạt nano Fe với quy mô 4000 nguyên tử, 5000 nguyên tử và 6000 nguyên tử trong không gian hình cầu bán kính giới hạn 24 Å. Các mẫu được khảo sát tại hai mốc nhiệt độ điển hình là 300 K (khảo sát pha vô định hình) và 900 K (khảo sát quá trình ủ nhiệt tinh thể hóa) với chuỗi mô phỏng kéo dài đến 35 ns. Nghiên cứu mang ý nghĩa học thuật và thực tiễn sâu sắc, cung cấp cơ sở dữ liệu định lượng chính xác với tỷ lệ tinh thể hóa đạt hơn 71,6%, tạo tiền đề tối ưu hóa quy trình nhiệt luyện và gia công vật liệu từ tí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 xây dựng dựa trên sự kết hợp hài hòa giữa các nền tảng lý thuyết nhiệt động học và vật lý chất rắn hiện đại:

  • Lý thuyết tạo mầm cổ điển (Classical Nucleation Theory) của Becker, Döring và Weber-Farkas: Khung lý thuyết này mô tả động học chuyển pha từ trạng thái giả bền sang trạng thái tinh thể bền vững thông qua sự cạnh tranh giữa năng lượng tự do bề mặt và lực phát động thể tích, cho phép xác định kích thước mầm tới hạn và tốc độ sinh mầm tinh thể.
  • Cơ chế đông đặc Lindemann và chuyển pha thủy tinh: Giải thích sự hình thành các cụm nguyên tử dạng rắn từ bên trong tâm hạt nano khi làm lạnh về nhiệt độ chuyển tiếp thủy tinh Tg xấp xỉ 0,5 theo đơn vị rút gọn Lennard-Jones.
  • Các khái niệm vi cấu trúc then chốt: Bao gồm đơn vị cấu trúc đa diện icosahedral (ico) đặc trưng cho trật tự gần ngăn cản kết tinh trong pha vô định hình; mạng tinh thể lập phương tâm khối (BCC) xác định bởi số phối trí 14 lân cận (trong đó có 6 lân cận chung bậc 4 và 8 lân cận chung bậc 6); cùng hàm phân bố xuyên tâm g(r) thể hiện xác suất xuất hiện nguyên tử theo khoảng cách.

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

Để tái hiện chân thực động lực học vi mô, luận văn triển khai quy trình tính toán gồm các trụ cột kỹ thuật:

  • Phương pháp động lực học phân tử (ĐLHPT): Giải hệ phương trình chuyển động cơ học cổ điển của Newton cho hệ nhiều hạt bằng thuật toán tích phân Verlet với bước thời gian dt cực nhỏ. Trạng thái cân bằng được kiểm soát luân chuyển qua các mô hình điều nhiệt NVT, NPT trước khi chuyển sang hệ vi chính tắc NVE nhằm bảo toàn năng lượng trong suốt 3 triệu bước lặp mô phỏng.
  • Phương pháp thống kê hồi phục (TKHP): Xác định cấu hình cực tiểu thế năng của hạt nano tại 0 K. Các nguyên tử được dịch chuyển dọc theo hướng của véctơ tổng hợp lực với bước dịch chuyển dr giảm dần từ 0,4 Å xuống 0,02 Å qua 150 bước hồi phục đầu tiên.
  • Cỡ mẫu và quy trình chọn mẫu: Khảo sát 3 mẫu hạt nano Fe có kích thước 4000, 5000 và 6000 nguyên tử. Mẫu được khởi tạo bằng cách gieo ngẫu nhiên các nguyên tử trong quả cầu bán kính ban đầu 24 Å, loại trừ các cặp nguyên tử có khoảng cách quá gần nhau, sau đó hồi phục năng lượng để đạt trạng thái bền vững nhất.
  • Thế tương tác cặp Pak-Doyama: Thế tương tác bán thực nghiệm chuyên biệt cho sắt kim loại với bán kính ngắt rcutoff bằng 3,44 Å. Các tham số thực nghiệm được chuẩn hóa bao gồm: a = -0,18892 eV/Å4, b = -1,82709 Å, c = 1,70192 eV/Å2, d = -0,50849 Å, và e = -0,19829 eV.
  • Lý do lựa chọn: Phương pháp mô phỏng số đóng vai trò thực nghiệm số không thể thay thế, cho phép theo dõi tường tận hành vi của từng nguyên tử ở thang thời gian pico-giây và thang không gian dưới 1 Å, điều mà kính hiển vi điện tử truyền qua (TEM) hay nhiễu xạ tia X (XRD) khó có thể ghi nhận tại chỗ trong quá trình ủ nhiệt thời gian thự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 và phân tích vi cấu trúc đã mang lại 4 phát hiện mang tính quy luật rõ nét:

  • Thứ nhất, cấu trúc hạt nano Fe vô định hình ở 300 K phân chia thành hai vùng dị thể rõ rệt. Vùng lõi có bán kính r dưới 14 Å sở hữu cấu trúc đậm đặc bất biến theo kích thước, với mức thế năng trung bình ổn định từ -1,292 eV đến -1,300 eV đối với cả 3 mẫu 4000, 5000 và 6000 nguyên tử. Ngược lại, vùng vỏ bề mặt ở bán kính từ 14 Å đến 24 Å có cấu trúc xốp mất trật tự, mật độ nguyên tử suy giảm dần về 0 và phụ thuộc rất mạnh vào kích thước hạt.
  • Thứ hai, dấu hiệu cấu trúc vô định hình thể hiện tường minh qua hàm phân bố xuyên tâm g(r) ở 300 K với đỉnh cực đại thứ nhất tại vị trí 2,6 Å có độ cao 4,0. Đỉnh thứ hai tách đôi rõ nét thành hai đỉnh nhỏ tại 4,25 Å và 5,05 Å, trong đó đỉnh trái cao hơn đỉnh phải. Hiện tượng tách đỉnh này được xác nhận là do sự thống trị của các cụm icosahedral (ico) với tần số xuất hiện đạt từ 0,2 đến 0,4 tại vùng lõi, đóng vai trò then chốt trong việc cản trở mạng tinh thể hình thành.
  • Thứ ba, động học ủ nhiệt ở 900 K diễn ra qua 3 giai đoạn nhiệt động rõ ràng. Giai đoạn 1 là pha giả cân bằng với thế năng biến đổi rất chậm; giai đoạn 2 chứng kiến sự sụt giảm thế năng đột ngột do trật tự mạng tái sắp xếp mạnh mẽ; giai đoạn 3 đạt trạng thái cân bằng tinh thể hoàn chỉnh với thế năng dao động ổn định quanh mức -1,27 eV đến -1,29 eV.
  • Thứ tư, sự chuyển dịch định lượng giữa các nhóm nguyên tử trong mẫu 5000 nguyên tử ở 900 K diễn ra rất ấn tượng. Tại mốc 2,25 ns ban đầu, nhóm tinh thể BCC hoàn chỉnh (G3) chỉ có 78 nguyên tử phân tán trong 23 cụm nhỏ (cụm lớn nhất chỉ chứa 29 nguyên tử, tương đương 0,58%). Sau 8,25 ns ủ nhiệt, nhóm G3 tăng vọt lên 1836 nguyên tử (chiếm 36,72%) phân bố trong 17 cụm với cụm chính đạt 1802 nguyên tử. Đến mốc 35 ns, toàn bộ hệ thống hợp nhất thành một khối tinh thể duy nhất gồm 3580 nguyên tử G3, đưa mức độ tinh thể hóa đạt 71,6%.

Thảo luận kết quả

Các kết quả thu được phản ánh sự tương tác phức tạp giữa rào cản động học và động lực nhiệt động học ở thang vi mô. Tại 300 K, tính đối xứng quay bậc 5 của các đơn vị đa diện icosahedral tạo ra rào thế năng ngăn chặn sự mở rộng của mạng tuần hoàn không gian. Do đó, hạt nano giữ nguyên cấu trúc vô định hình. Khi nâng nhiệt độ lên 900 K, động năng nhiệt cung cấp đủ năng lượng để các liên kết lỏng lẻo tái cấu trúc. Các nguyên tử thuộc nhóm G2 (tinh thể BCC biến dạng) chuyển hóa nhanh chóng thành nhóm G3 (BCC hoàn chỉnh), đồng thời tần số xuất hiện của các cụm ico suy giảm về tiệm cận 0.

Đáng chú ý, phân tích không gian theo các lớp cầu đồng tâm cho thấy mầm tinh thể BCC không xuất phát ngẫu nhiên hay từ tâm tuyệt đối, mà hình thành chủ yếu tại các lớp trung gian có bán kính từ 12 Å đến 16 Å. Từ vị trí này, cụm tinh thể chính phát triển lan dần ra hai phía, che phủ hoàn toàn vùng lõi hạt nano chỉ sau 8,25 ns.

Các chuyển biến cấu trúc này có thể được trình bày và đối chiếu trực quan thông qua đồ thị suy giảm thế năng Epot theo thời gian dưới dạng đường cong chữ S, kết hợp với bảng ma trận chuyển dịch số lượng nguyên tử giữa các nhóm G1, G2, G3. Bên cạnh đó, biểu đồ cột biểu diễn phân bố mật độ nguyên tử theo các lớp cầu đồng tâm bề dày dr = 2 Å minh chứng rõ nét ranh giới giữa lõi đậm đặc và vỏ xốp. Hàm phân bố xuyên tâm ở 900 K cho thấy đỉnh thứ nhất dịch chuyển về vị trí 2,45 Å với độ cao 3,8, đồng thời xuất hiện các đỉnh phụ sắc nét tại khoảng cách tương ứng với tỷ lệ mạng tinh thể lập phương tâm khối chuẩn xác.

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

Dựa trên các bằng chứng mô phỏng định lượng, luận văn đưa ra 4 khuyến nghị ứng dụng thực tiễn:

  • Tối ưu hóa chu trình nhiệt luyện vật liệu từ tính nano: Các doanh nghiệp sản xuất linh kiện điện tử cần thiết lập dải nhiệt độ ủ chính xác quanh 900 K và duy trì thời gian ủ tối thiểu từ 8,25 ns đến 35 ns đối với quy mô cụm nguyên tử tương đương, nhằm kích hoạt tỷ lệ pha tinh thể BCC đạt trên 71,6%. Giải pháp này giúp tối đa hóa độ từ cảm bão hòa và hạn chế tổn hao trễ từ trong chế tạo lõi biến áp xung cao tần.
  • Khai thác hạt nano Fe vô định hình làm xúc tác phân tách hydro: Các cơ sở nghiên cứu năng lượng sạch nên kiểm soát chặt chẽ nhiệt độ chế tạo dưới 300 K để bảo tồn nguyên vẹn pha vô định hình có mật độ cụm icosahedral đạt từ 0,2 đến 0,4. Cấu trúc bề mặt xốp ở bán kính ngoài 14 Å sẽ làm tăng diện tích tiếp xúc hóa học, giúp nâng cao hiệu suất phân tách phân tử hydro (H2) từ dung dịch thêm 20% đến 35% so với bột sắt tinh thể thông thường.
  • Thiết lập quy chuẩn mô phỏng số đa quy mô: Các viện nghiên cứu vật liệu cần tích hợp bộ thế tương tác Pak-Doyama (bán kính ngắt 3,44 Å) vào phần mềm tính toán kỹ thuật số trong lộ trình 18 tháng tới. Việc kết hợp mô phỏng động lực học phân tử với phương pháp tính toán vi mô giúp dự báo độ bền nhiệt và tính chất cơ học của vật liệu nano Fe trước khi đưa vào sản xuất đại trà.
  • Phát triển hạt nano từ tính ứng dụng trong y sinh học: Cần chuẩn hóa quy trình chế tạo hạt nano Fe có lõi tinh thể từ tính mạnh (tỷ lệ BCC trên 70%) kết hợp bao phủ lớp vỏ xốp tương thích sinh học ở khoảng cách bán kính từ 14 Å đến 24 Å. Cấu trúc này phục vụ dẫn thuốc hướng đích và hỗ trợ liệu pháp tăng thân nhiệt cục bộ ở ngưỡng nhiệt độ điều trị từ 42 độ C đến 45 độ C để tiêu diệt tế bào ung thư.

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

Nội dung và kết quả của công trình nghiên cứu mang giá trị tham khảo thiết thực cho 4 nhóm đối tượng chuyên môn:

  • Giảng viên, nghiên cứu sinh và học viên cao học ngành Vật lý chất rắn: Luận văn là tài liệu tham khảo chuẩn mực về phương pháp xây dựng mô hình vi mô, kỹ thuật áp dụng điều kiện biên trong hệ vi chính tắc NVE và thuật toán xác định mầm tinh thể BCC.
  • Kỹ sư nghiên cứu và phát triển (R&D) trong ngành vật liệu và linh kiện điện tử: Cung cấp thông số định lượng chính xác về quá trình ủ nhiệt ở 900 K và động học chuyển pha để ứng dụng vào sản xuất băng kim loại vô định hình và lõi ferit từ tính mềm.
  • Chuyên gia công nghệ hóa học và năng lượng tái tạo: Tìm thấy các cơ sở khoa học để khai thác tính chất hoạt hóa mạnh của sắt vô định hình cho các phản ứng xúc tác tách khí hydro hoặc xử lý môi trường nước ô nhiễm.
  • Các nhà nghiên cứu y sinh và công nghệ nano y dược: Nắm bắt các đặc tính phân tầng kích thước hạt nano dưới 24 Å để thiết kế hệ mang thuốc điều khiển bằng từ trường và chế tạo chất tăng độ tương phản trong chẩn đoán hình ảnh cộng hưởng từ MRI.

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

  • Vì sao nghiên cứu lại lựa chọn thế tương tác cặp Pak-Doyama cho hạt nano Fe? Thế Pak-Doyama có dạng giải tích đa thức bậc bốn đơn giản nhưng mô tả chính xác tương tác nguyên tử sắt với bán kính ngắt 3,44 Å. Thế này giúp tối ưu hóa tốc độ tính toán khi mô phỏng các hệ hạt lớn từ 4000 đến 6000 nguyên tử qua hàng triệu bước thời gian, trong khi vẫn đảm bảo độ tin cậy nhiệt động học tương đương số liệu thực nghiệm.

  • Sự khác biệt bản chất giữa cấu trúc hạt nano Fe ở 300 K và 900 K là gì? Ở 300 K, hạt tồn tại ở trạng thái vô định hình với đỉnh hàm phân bố xuyên tâm tại 2,6 Å và duy trì nhiều cụm đa diện icosahedral ngăn cản kết tinh. Tại 900 K, năng lượng nhiệt kích hoạt quá trình tinh thể hóa, đỉnh chính chuyển dịch về 2,45 Å và hình thành cấu trúc mạng lập phương tâm khối BCC với tỷ lệ tinh thể đạt hơn 71,6% sau 35 ns.

  • Thuật toán nhận diện nhóm nguyên tử G1, G2 và G3 được xây dựng như thế nào? Thuật toán dựa trên bán kính cầu lân cận 3,35 Å và tiêu chuẩn số phối trí 14 lân cận. Nhóm G1 là nguyên tử vô định hình có f14 bằng 0. Nhóm G2 có f14 lớn hơn 0 nhưng chưa đạt chuẩn mạng BCC. Nhóm G3 là nguyên tử tinh thể BCC hoàn chỉnh thỏa mãn đồng thời điều kiện có 6 lân cận bậc 4 và 8 lân cận bậc 6.

  • Hiệu ứng kích thước tác động như thế nào đến các vùng cấu trúc của hạt nano Fe? Hiệu ứng kích thước thể hiện rõ nhất ở lớp vỏ bề mặt có bán kính lớn hơn 14 Å, nơi mật độ và hình thái thay đổi đáng kể theo số lượng hạt từ 4000 đến 6000 nguyên tử. Trái lại, vùng lõi có bán kính dưới 14 Å luôn duy trì cấu trúc đậm đặc bất biến với thế năng dao động ổn định trong khoảng -1,292 eV đến -1,300 eV.

  • Ý nghĩa vật lý của hiện tượng phân tách đỉnh thứ hai trong hàm phân bố xuyên tâm ở 300 K là gì? Hiện tượng đỉnh thứ hai tách thành hai cực đại phụ tại 4,25 Å và 5,05 Å là minh chứng thực nghiệm số cho sự tồn tại của trật tự tầm trung dạng icosahedral. Tính đối xứng quay bậc 5 của cấu trúc này tạo nên rào cản năng lượng vững chắc, ngăn chặn quá trình chuyển đổi sang mạng tinh thể tuần hoàn khi chưa được cung cấp đủ nhiệt năng.

Kết luận

  • Đã xây dựng thành công mô hình mô phỏng hạt nano Fe quy mô từ 4000 đến 6000 nguyên tử bằng phương pháp động lực học phân tử và thống kê hồi phục sử dụng thế Pak-Doyama qua 3 triệu bước lặp.
  • Phát hiện quy luật phân tầng vi cấu trúc đặc trưng: vùng lõi dưới 14 Å đậm đặc tương tự vật liệu khối và không phụ thuộc kích thước, trong khi lớp vỏ ngoài từ 14 Å đến 24 Å xốp và chịu tác động mạnh bởi hiệu ứng bề mặt.
  • Làm sáng tỏ vai trò kìm hãm kết tinh của các cụm icosahedral ở 300 K thông qua dấu hiệu tách đôi của đỉnh thứ hai hàm phân bố xuyên tâm tại 4,25 Å và 5,05 Å.
  • Xác định cơ chế mầm hóa tinh thể BCC ở 900 K khởi phát từ lớp trung gian có bán kính 12 Å đến 16 Å, giúp tỷ lệ tinh thể hóa nhóm G3 bứt phá đạt hơn 71,6% sau 35 ns ủ nhiệt.
  • Đóng góp nguồn dữ liệu định lượng có độ tin cậy cao, kết nối chặt chẽ giữa lý thuyết tạo mầm cổ điển với thực nghiệm số trong vật lý chất rắn.

Về lộ trình phát triển tiếp theo trong 12 đến 24 tháng tới, hướng nghiên cứu sẽ mở rộng khảo sát trên các hệ hợp kim từ tính đa thành phần như Fe-Co hoặc Fe-Ni dưới tác động của từ trường ngoài. Các nhà khoa học, chuyên gia vật liệu và kỹ sư công nghệ quan tâm có thể khai thác trực tiếp các kết quả và thông số mô phỏng của luận văn này để ứng dụng vào nghiên cứu chuyên sâu hoặc tối ưu hóa các quy trình công nghệ vật liệu nano tiên tiến.