Mô Phỏng Quá Trình Tạo Màng Fe Phẳng Hai Chiều Từ Trạng Thái Lỏng

Luận văn thạc sĩ kỹ thuật nghiên cứu vật lý kỹ thuật mô phỏng quá trình tạo màng fe phẳng hai chiều từ trạng thái lỏng, khảo sát thực trạng, phân tích nguyên nhân, đề xuất giải

• Tác giả: Vương Phú Tài
• Người hướng dẫn: GS. Võ Văn Hoàng
• Trường học: Đại học Quốc gia TP. Hồ Chí Minh
• Chuyên ngành: Vật lý kỹ thuật
• Đề tài: Mô Phỏng Quá Trình Tạo Màng Sắt Hai Chiều Từ Trạng Thái Lỏng
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2015
• Địa điểm: TP. Hồ Chí Minh

Mục lục chi tiết

• LỜI CAM ĐOAN
• 1. CHƯƠNG 1: TỔNG QUAN TÀI LIỆU
• 1.1. Tổng quan về vật liệu hai chiều
• 1.2. Tổng quan nghiên cứu thực nghiệm và mô phỏng của màng sắt hai chiều
• 1.2.1. Nghiên cứu thực nghiệm của màng sắt hai chiều
• 1.2.2. Nghiên cứu mô phỏng của màng sắt hai chiều
• 2. CHƯƠNG 2: PHƯƠNG PHÁP NGHIÊN CỨU
• 2.1. Sơ lược về mô phỏng động lực học phân tử
• 2.2. Cấu trúc của một chương trình mô phỏng động lực học phân tử
• 2.3. Mô phỏng động lực học phân tử cho màng sắt phẳng hai chiều
• 2.4. Các chi tiết về mô hình vật liệu được mô phỏng
• 2.5. Thế tương tác EAM. Các điều kiện biên
• 2.6. Các phần mềm sử dụng trong luận văn
• 3. CHƯƠNG 3: KẾT QUẢ VÀ THẢO LUẬN
• 3.1. Các tính chất nhiệt động lực học của quá trình đông đặc
• 3.2. Sự biến đổi năng lượng toàn phần theo nhiệt độ cho mỗi nguyên tử
• 3.3. Sự biến đổi của nhiệt dung riêng theo nhiệt độ trên mỗi nguyên tử
• 3.4. Sự thay đổi cấu trúc mô hình trong quá trình đông đặc
• 3.5. Sự thay đổi của hàm phân bố xuyên tâm
• 3.6. Trật tự định hướng các liên kết nguyên tử
• 3.7. Sự thay đổi số phối vị
• 3.8. Sự thay đổi độ dài liên kết giữa các nguyên tử
• 3.9. Sự thay đổi phân bố góc
• 3.10. Cấu trúc chi tiết của mô hình ở nhiệt độ 300 K
• 3.11. Phân bố số phối vị
• 3.12. Phân bố độ dài liên kết giữa các nguyên tử
• 3.13. Phân bố góc liên kết
• 3.14. Hình ảnh nhiễu xạ của mô hình đạt được ở 300 K
• 3.15. Các dạng khuyết tật
• 3.16. Cơ chế nguyên tử của quá trình chuyển pha
• 3.17. Sự phụ thuộc vào nhiệt độ của các nguyên tử có số phối vị Z = 6
• 3.18. Thể hiện trực quan sự chuyển pha của mô hình
• 4. CHƯƠNG 4: KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN CỦA ĐỀ TÀI
• 4.1. Hướng phát triển của đề tài
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng Quan Mô Phỏng Quá Trình Tạo Màng Sắt 2D Từ Lỏng

Nghiên cứu mô phỏng quá trình tạo màng sắt 2D từ trạng thái lỏng ngày càng trở nên quan trọng. Việc tạo màng ironene (màng sắt đơn lớp) bằng thực nghiệm gặp nhiều khó khăn do đặc điểm cấu trúc vi mô hai chiều. Xác định nhiệt độ đông đặc và khảo sát hiện tượng đông đặc bằng thực nghiệm là thách thức lớn. Mô phỏng là giải pháp tối ưu, cung cấp thông tin chi tiết về cơ chế nguyên tử trong quá trình đông đặc. Các nghiên cứu trước đây đã khảo sát hằng số mạng và moment từ, nhưng quá trình chuyển pha khi làm lạnh từ trạng thái lỏng vẫn chưa được nghiên cứu kỹ lưỡng. Nghiên cứu này tập trung vào các tính chất nhiệt động lực học và sự thay đổi cơ chế nguyên tử của quá trình đông đặc. Học viên quyết định thực hiện nghiên cứu “Mô phỏng quá trình tạo màng Fe phẳng hai chiều từ trạng thái lỏng”.

1.1. Giới thiệu về Màng Sắt Hai Chiều Ironene

Màng sắt hai chiều, hay còn gọi là ironene, là vật liệu tiềm năng cho các ứng dụng trong tương lai. Cấu trúc đơn lớp nguyên tử của màng sắt 2D mang lại những tính chất độc đáo so với vật liệu sắt khối. Khảo sát quá trình hình thành màng sắt từ trạng thái lỏng giúp hiểu rõ hơn về cấu trúc và tính chất của vật liệu này. Các nghiên cứu mô phỏng đóng vai trò quan trọng trong việc khám phá những đặc tính của màng sắt 2D trước khi thực hiện các thí nghiệm thực tế.

1.2. Vai trò của Mô Phỏng trong Nghiên cứu Vật Liệu 2D

Do tính phức tạp của việc nghiên cứu màng sắt hai chiều bằng thực nghiệm, mô phỏng trở thành công cụ không thể thiếu. Mô phỏng giúp dự đoán và giải thích các hiện tượng xảy ra trong quá trình tạo màng sắt. Các phương pháp mô phỏng động lực học phân tử (MD) và mô phỏng Monte Carlo được sử dụng rộng rãi để nghiên cứu màng sắt 2D. Kết quả mô phỏng cung cấp những thông tin quan trọng để định hướng các nghiên cứu thực nghiệm và phát triển ứng dụng.

II. Thách Thức Nghiên Cứu Trạng Thái Lỏng Trong Môi Trường 2D

Nghiên cứu trạng thái lỏng và quá trình đông đặc trong không gian hai chiều đặt ra nhiều thách thức. Các lý thuyết và mô hình vật lý được phát triển cho vật liệu khối ba chiều có thể không còn áp dụng được cho màng sắt 2D. Sự dao động nhiệt của các nguyên tử trong mạng 2D có thể dẫn đến sự tan chảy mạng tinh thể. Khảo sát sự thay đổi cấu trúc và tính chất của màng sắt trong quá trình chuyển pha từ lỏng sang rắn đòi hỏi các phương pháp mô phỏng tiên tiến và khả năng tính toán mạnh mẽ. Theo Peierls và Landau, sự tồn tại của vật liệu 2D mâu thuẫn với lý thuyết bền vững về mặt nhiệt động lực học.

2.1. Ảnh hưởng của Nhiệt Độ Màng Đến Quá Trình Đông Đặc

Nhiệt độ màng là yếu tố quan trọng ảnh hưởng đến quá trình đông đặc của màng sắt 2D. Khi nhiệt độ màng giảm, các nguyên tử mất dần động năng và bắt đầu hình thành các liên kết. Quá trình kết tinh xảy ra khi nhiệt độ màng đạt đến một giá trị nhất định. Mô phỏng giúp xác định chính xác nhiệt độ màng tại đó quá trình chuyển pha xảy ra.

2.2. Vai trò của Năng Lượng Bề Mặt trong Cấu Trúc Màng

Năng lượng bề mặt đóng vai trò quan trọng trong việc xác định cấu trúc của màng sắt 2D. Các nguyên tử trên bề mặt có số liên kết ít hơn so với các nguyên tử bên trong, dẫn đến năng lượng bề mặt cao hơn. Xu hướng giảm thiểu năng lượng bề mặt thúc đẩy quá trình hình thành các cấu trúc ổn định. Mô phỏng giúp tính toán năng lượng bề mặt và dự đoán cấu trúc của màng sắt 2D.

2.3. Tác động của Áp Suất Hơi và Tốc Độ Bay Hơi

Áp suất hơi và tốc độ bay hơi ảnh hưởng đến quá trình tạo màng sắt 2D từ trạng thái lỏng. Áp suất hơi cao có thể dẫn đến sự ngưng tụ nhanh chóng, trong khi tốc độ bay hơi quá cao có thể làm cho màng không đồng đều. Điều chỉnh các thông số này là rất quan trọng để tạo ra các màng sắt 2D chất lượng cao. Mô phỏng giúp tối ưu hóa các thông số này để đạt được kết quả tốt nhất.

III. Phương Pháp Mô Phỏng Động Lực Học Phân Tử Cho Màng Sắt

Mô phỏng động lực học phân tử (MD) là phương pháp mạnh mẽ để nghiên cứu quá trình tạo màng sắt 2D từ trạng thái lỏng. Phương pháp này mô phỏng chuyển động của các nguyên tử dựa trên các định luật vật lý cổ điển và các thế tương tác giữa các nguyên tử. Kết quả mô phỏng MD cung cấp thông tin chi tiết về vị trí, vận tốc và năng lượng của các nguyên tử theo thời gian. Các nghiên cứu trước đây đã khảo sát và tính toán được hằng số mạng, độ lớn moment từ.

3.1. Giới thiệu về Thế Tương Tác EAM

Thế tương tác EAM (Embedded Atom Method) là một trong những thế tương tác phổ biến được sử dụng trong mô phỏng MD cho kim loại. Thế EAM mô tả năng lượng của mỗi nguyên tử dựa trên mật độ điện tử xung quanh nguyên tử đó. Thế EAM cho phép mô phỏng chính xác các tính chất của kim loại, bao gồm cả màng sắt. Các chi tiết về mô hình vật liệu được mô phỏng.

3.2. Chi tiết về Điều Kiện Biên trong Mô Phỏng

Điều kiện biên đóng vai trò quan trọng trong mô phỏng MD. Điều kiện biên tuần hoàn thường được sử dụng để mô phỏng hệ vô hạn với số lượng nguyên tử hữu hạn. Điều kiện biên giúp giảm thiểu ảnh hưởng của các bề mặt và cạnh của hệ thống mô phỏng. Minh họa điều kiện biên tuần hoàn.

3.3. Các Phần Mềm Sử Dụng trong Mô Phỏng Vật Liệu

Có nhiều phần mềm khác nhau được sử dụng trong mô phỏng vật liệu, bao gồm cả màng sắt 2D. Các phần mềm này cung cấp các công cụ để xây dựng mô hình, thực hiện mô phỏng, và phân tích kết quả. Lựa chọn phần mềm phù hợp phụ thuộc vào yêu cầu cụ thể của từng nghiên cứu.

IV. Kết Quả Tính Chất Nhiệt Động Lực Học Của Quá Trình Đông Đặc

Các kết quả mô phỏng MD cho thấy sự thay đổi rõ rệt trong các tính chất nhiệt động lực học của màng sắt trong quá trình đông đặc. Năng lượng toàn phần giảm khi nhiệt độ giảm, cho thấy sự hình thành các liên kết ổn định. Nhiệt dung riêng có một đỉnh tại nhiệt độ đông đặc, cho thấy sự chuyển pha từ lỏng sang rắn. Các hàm phân bố xuyên tâm (RDF) cho thấy sự thay đổi trong trật tự nguyên tử khi màng sắt đông đặc. Nhiệt độ nóng chảy của ironene tinh thể trong không gian 2D được xác định vào khoảng 2640 K.

4.1. Sự Biến Đổi Năng Lượng Toàn Phần Theo Nhiệt Độ Màng

Sự biến đổi năng lượng toàn phần theo nhiệt độ màng cung cấp thông tin quan trọng về quá trình đông đặc. Khi nhiệt độ màng giảm, năng lượng toàn phần giảm do sự hình thành các liên kết giữa các nguyên tử. Đồ thị năng lượng toàn phần theo nhiệt độ màng có thể được sử dụng để xác định nhiệt độ đông đặc.

4.2. Sự Biến Đổi Nhiệt Dung Riêng Theo Nhiệt Độ Màng

Sự biến đổi nhiệt dung riêng theo nhiệt độ màng cũng cung cấp thông tin về quá trình chuyển pha. Nhiệt dung riêng có một đỉnh tại nhiệt độ đông đặc, cho thấy sự hấp thụ hoặc giải phóng nhiệt trong quá trình chuyển pha. Đỉnh nhiệt dung riêng có thể được sử dụng để xác định nhiệt độ đông đặc chính xác hơn.

4.3. Phân Tích Hàm Phân Bố Xuyên Tâm RDF

Hàm phân bố xuyên tâm (RDF) mô tả sự phân bố không gian của các nguyên tử trong màng sắt. Khi màng sắt chuyển từ trạng thái lỏng sang rắn, các đỉnh trong RDF trở nên sắc nét hơn, cho thấy sự hình thành trật tự nguyên tử. Phân tích RDF giúp xác định cấu trúc tinh thể của màng sắt sau khi đông đặc.

V. Ứng Dụng Tiềm Năng của Màng Sắt 2D Trong Vật Liệu Nano

Màng sắt 2D có tiềm năng ứng dụng to lớn trong lĩnh vực vật liệu nano. Với độ dày chỉ một lớp nguyên tử, màng sắt 2D có diện tích bề mặt lớn, thích hợp cho các ứng dụng xúc tác và cảm biến. Tính chất từ của màng sắt 2D cũng có thể được khai thác trong các thiết bị lưu trữ từ tính mật độ cao. Nghiên cứu về ứng dụng màng sắt trong các thiết bị nano vẫn đang được tiến hành.

5.1. Ứng Dụng Màng Sắt 2D Trong Xúc Tác

Diện tích bề mặt lớn của màng sắt 2D làm cho nó trở thành vật liệu xúc tác tiềm năng. Các nguyên tử trên bề mặt có hoạt tính cao, có thể thúc đẩy các phản ứng hóa học. Màng sắt 2D có thể được sử dụng trong các ứng dụng như xử lý khí thải và sản xuất năng lượng sạch.

5.2. Ứng Dụng Màng Sắt 2D Trong Cảm Biến

Màng sắt 2D có thể được sử dụng để chế tạo các cảm biến nhạy bén. Sự thay đổi trong tính chất điện hoặc từ của màng sắt 2D khi tiếp xúc với các chất khác có thể được sử dụng để phát hiện các chất đó. Màng sắt 2D có thể được sử dụng trong các ứng dụng như phát hiện khí độc và giám sát môi trường.

5.3. Sử dụng Màng Sắt 2D Trong Lưu Trữ Dữ Liệu Từ Tính

Tính chất từ tính của màng sắt 2D mở ra khả năng sử dụng nó trong các thiết bị lưu trữ dữ liệu từ tính. Kích thước nhỏ của màng sắt 2D cho phép tăng mật độ lưu trữ dữ liệu. Nghiên cứu về các thiết bị lưu trữ từ tính dựa trên màng sắt 2D vẫn đang được tiến hành.

VI. Kết Luận Tiềm Năng Và Hướng Phát Triển Của Nghiên Cứu Màng Sắt

Nghiên cứu mô phỏng quá trình tạo màng sắt 2D từ trạng thái lỏng đã cung cấp những hiểu biết sâu sắc về cấu trúc và tính chất của vật liệu này. Các kết quả mô phỏng cho thấy tiềm năng ứng dụng to lớn của màng sắt 2D trong nhiều lĩnh vực. Các nghiên cứu trong tương lai nên tập trung vào việc khám phá các phương pháp tổng hợp màng sắt 2D chất lượng cao và phát triển các ứng dụng cụ thể của vật liệu này. Sau quá trình làm lạnh, đề tài thu được màng sắt phẳng hai chiều có cấu trúc lưới tam giác với hằng số mạng là 2,455 Å, chứa khuyết tật mất nguyên tử (vacancy).

6.1. Hướng Phát Triển Của Đề Tài

Đề tài có thể được phát triển theo nhiều hướng khác nhau. Một hướng là nghiên cứu ảnh hưởng của các tạp chất và khuyết tật đến tính chất của màng sắt 2D. Một hướng khác là nghiên cứu các phương pháp điều chỉnh tính chất của màng sắt 2D để đáp ứng các yêu cầu cụ thể của từng ứng dụng.

6.2. Nghiên cứu các Phương pháp Tổng Hợp Màng Sắt 2D

Phát triển các phương pháp tổng hợp màng sắt 2D chất lượng cao là rất quan trọng để khai thác tiềm năng của vật liệu này. Các phương pháp như lắng đọng hơi hóa học (CVD) và phún xạ có thể được sử dụng để tổng hợp màng sắt 2D. Nghiên cứu nên tập trung vào việc tối ưu hóa các điều kiện tổng hợp để đạt được các màng sắt 2D có cấu trúc hoàn hảo và tính chất mong muốn.

6.3. Phân tích cấu trúc và khuyết tật của Màng Sắt 2D

Phân tích cấu trúc chi tiết của màng sắt 2D ở nhiệt độ 300K, phân bố số phối vị, phân bố độ dài liên kết giữa các nguyên tử, phân bố góc liên kết, và ảnh nhiễu xạ của mô hình đạt được ở 300 K. Các dạng khuyết tật ảnh hưởng đến tính chất của màng.

Nội dung chính

Tổng quan nghiên cứu

Vật liệu hai chiều (2D) đã trở thành chủ đề nghiên cứu nổi bật trong khoa học vật liệu hiện đại, đặc biệt sau thành công của graphene vào năm 2004. Trong đó, màng sắt đơn lớp (ironene) với cấu trúc phẳng hai chiều được xem là vật liệu từ tiềm năng, có thể ứng dụng trong các thiết bị lưu trữ với mật độ cao và kích thước nhỏ gọn. Tuy nhiên, việc nghiên cứu quá trình tạo màng Fe 2D từ trạng thái lỏng gặp nhiều khó khăn do đặc tính cấu trúc và tính chất nhiệt động học phức tạp của kim loại. Nhiệt độ đông đặc và cơ chế chuyển pha của màng sắt 2D chưa được khảo sát kỹ lưỡng, đặc biệt là trong không gian hai chiều.

Mục tiêu của luận văn là mô phỏng quá trình tạo màng Fe phẳng hai chiều từ trạng thái lỏng, khảo sát sự thay đổi cấu trúc và các tính chất nhiệt động học trong quá trình chuyển pha từ lỏng sang rắn. Nghiên cứu được thực hiện trong phạm vi mô hình ironene gồm 6400 nguyên tử, sử dụng phương pháp mô phỏng động lực học phân tử (MD) cổ điển với thế tương tác EAM, trong khoảng thời gian làm lạnh từ 4300 K xuống 300 K. Kết quả nghiên cứu cung cấp thông tin chi tiết về nhiệt độ chuyển pha, cấu trúc tinh thể, các khuyết tật và cơ chế nguyên tử trong quá trình đông đặc màng Fe 2D.

Nghiên cứu có ý nghĩa quan trọng trong việc hiểu rõ cơ chế chuyển pha và cấu trúc của vật liệu 2D kim loại, góp phần mở rộng kiến thức về vật liệu từ 2D và hỗ trợ phát triển các ứng dụng công nghệ nano, vi điện tử và lưu trữ dữ liệu.

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

Khung lý thuyết áp dụng

Luận văn dựa trên hai lý thuyết và mô hình chính:

1. Lý thuyết động lực học phân tử (Molecular Dynamics - MD): Phương pháp này mô phỏng chuyển động của các nguyên tử trong hệ dựa trên các định luật cơ học cổ điển Newton. Phương trình chuyển động được giải bằng thuật toán Verlet với bước thời gian khoảng 10⁻¹⁵ giây, cho phép theo dõi sự biến đổi vị trí và vận tốc của từng nguyên tử trong quá trình làm lạnh.

2. Thế tương tác nhúng nguyên tử (Embedded Atom Method - EAM): Đây là thế tương tác bán thực nghiệm phù hợp cho các kim loại chuyển tiếp như Fe, cho phép mô phỏng chính xác các tương tác nhiều nguyên tử và các hiệu ứng bề mặt, khuyết tật. Tổng năng lượng hệ được biểu diễn gồm năng lượng tương tác cặp và năng lượng nhúng của nguyên tử vào mật độ electron xung quanh.

Các khái niệm chính được sử dụng bao gồm:

• Hàm phân bố xuyên tâm (Radial Distribution Function - RDF): Mô tả sự phân bố khoảng cách giữa các nguyên tử, phản ánh cấu trúc trật tự trong hệ.
• Số phối vị (Coordination Number): Số nguyên tử lân cận trong mặt cầu phối vị đầu tiên, dùng để đánh giá mức độ kết tinh và khuyết tật.
• Trật tự định hướng liên kết (Bond Orientation Order - ψ₆): Đại lượng đo độ đối xứng và trật tự của các liên kết nguyên tử trong mạng tinh thể 2D.
• Nhiệt dung riêng (Specific Heat Capacity): Đại lượng nhiệt động học phản ánh sự thay đổi năng lượng toàn phần theo nhiệt độ, dùng để xác định nhiệt độ chuyển pha và bậc chuyển pha.

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

Nguồn dữ liệu chính là mô hình mô phỏng ironene gồm 6400 nguyên tử sắt, được xây dựng dựa trên các thông số cấu trúc thu được từ tính toán ab initio. Mô hình được nung nóng đến 4300 K để đạt trạng thái lỏng, sau đó làm lạnh tuyến tính với tốc độ 2×10¹⁰ K/s xuống 300 K.

Phương pháp phân tích bao gồm:

• Tính toán năng lượng toàn phần và nhiệt dung riêng theo nhiệt độ để xác định nhiệt độ chuyển pha.
• Phân tích hàm phân bố xuyên tâm g(r), số phối vị, phân bố góc và độ dài liên kết để khảo sát sự thay đổi cấu trúc.
• Đánh giá trật tự định hướng liên kết ψ₆ để xác định mức độ kết tinh.
• Quan sát các khuyết tật cấu trúc như khuyết tật nút khuyết (vacancy).
• Trực quan hóa sự chuyển pha qua phân bố nguyên tử có số phối vị Z=6.

Các phần mềm sử dụng gồm LAMMPS cho mô phỏng MD, Fortran để tạo mô hình và xử lý dữ liệu, Origin để biểu diễn đồ thị, VESTA để minh họa cấu trúc. Điều kiện biên áp dụng là biên tuần hoàn theo trục x và biên phản xạ đàn hồi theo trục y.

Thời gian nghiên cứu từ tháng 7 đến tháng 12 năm 2015, với cỡ mẫu 6400 nguyên tử được chọn nhằm đảm bảo tính đại diện và khả năng tính toán hiệu quả.

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

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

1. Nhiệt độ chuyển pha và bậc chuyển pha:
   Nhiệt độ đông đặc của màng Fe 2D được xác định khoảng 2640 K dựa trên sự thay đổi đột ngột của năng lượng toàn phần và đỉnh riêng biệt của nhiệt dung riêng. Quá trình chuyển pha thuộc loại I, không có pha trung gian, thể hiện qua đồ thị năng lượng và nhiệt dung riêng.

2. Sự thay đổi cấu trúc qua hàm phân bố xuyên tâm (RDF):
   Ở 4300 K, RDF mịn với đỉnh thấp, đặc trưng chất lỏng. Khi nhiệt độ giảm gần 2600 K, các đỉnh phụ xuất hiện rõ, độ cao đỉnh chính tăng, cho thấy sự hình thành cấu trúc tinh thể. Ở 300 K, RDF thể hiện cấu trúc kết tinh hoàn chỉnh với các đỉnh rõ nét.

3. Trật tự định hướng liên kết ψ₆ và số phối vị:
   Giá trị ψ₆ tăng từ gần 0 ở trạng thái lỏng lên 0,943 ở 300 K, chứng tỏ sự hình thành mạng lưới tam giác tinh thể. Số phối vị trung bình tăng từ 4,796 ở trạng thái lỏng lên 5,891 ở 300 K, gần với giá trị lý tưởng 6 của lưới tam giác. Tại nhiệt độ chuyển pha, khoảng 91% nguyên tử có số phối vị Z=6, còn lại là các khuyết tật.

4. Khuyết tật cấu trúc:
   Mô hình chứa khoảng 5% nguyên tử có số phối vị khác 6, chủ yếu là khuyết tật nút khuyết (vacancy). Đây là dạng sai hỏng điểm phổ biến trong kim loại, ảnh hưởng đến tính chất từ và cơ học của màng Fe 2D.

Thảo luận kết quả

Sự xác định nhiệt độ chuyển pha ở khoảng 2640 K phù hợp với các nghiên cứu lý thuyết và thực nghiệm trước đây, đồng thời khẳng định tính khả thi của mô phỏng MD với thế EAM trong việc mô phỏng quá trình đông đặc vật liệu 2D kim loại. Đồ thị năng lượng và nhiệt dung riêng minh họa rõ ràng quá trình chuyển pha loại I, không có pha trung gian, phù hợp với đặc tính của màng Fe 2D.

Phân tích RDF, ψ₆ và số phối vị cho thấy quá trình kết tinh diễn ra dần dần khi làm lạnh, với sự hình thành mạng lưới tam giác đặc trưng. Sự tồn tại của khuyết tật cấu trúc là điều không tránh khỏi do ảnh hưởng của điều kiện biên và tốc độ làm lạnh cao, tương tự như các kết quả thực nghiệm về màng sắt đơn lớp trên nền graphene.

So sánh với các nghiên cứu thực nghiệm, mô hình mô phỏng không quan sát được sự chuyển về trạng thái vô định hình sau khi kết tinh do quá trình làm lạnh liên tục và nhanh, trong khi thực nghiệm có thể xảy ra do tác động của bức xạ điện tử và điều kiện giữ mẫu. Kết quả này cung cấp cái nhìn sâu sắc về cơ chế nguyên tử trong quá trình đông đặc, hỗ trợ phát triển các vật liệu 2D kim loại với tính chất điều khiển được.

Dữ liệu có thể được trình bày qua các biểu đồ năng lượng toàn phần, nhiệt dung riêng, đồ thị RDF, biểu đồ phân bố số phối vị và ψ₆ theo nhiệt độ, cùng hình ảnh trực quan 2D thể hiện sự xuất hiện của nguyên tử có số phối vị Z=6.

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

1. Tối ưu hóa tốc độ làm lạnh trong mô phỏng:
   Giảm tốc độ làm lạnh xuống khoảng 10⁹ K/s để quan sát rõ hơn pha trung gian và sự chuyển đổi giữa trạng thái tinh thể và vô định hình, giúp mô phỏng sát thực tế hơn. Chủ thể thực hiện: nhóm nghiên cứu vật lý tính toán, thời gian 6 tháng.

2. Mở rộng mô hình với kích thước lớn hơn và điều kiện biên đa dạng:
   Tăng số nguyên tử lên khoảng 10.000-20.000 để giảm ảnh hưởng biên và khảo sát các dạng khuyết tật phức tạp hơn. Áp dụng điều kiện biên khác như biên tuần hoàn cả hai chiều. Chủ thể: phòng thí nghiệm vật liệu nano, thời gian 1 năm.

3. Kết hợp mô phỏng với các phương pháp tính toán ab initio:
   Sử dụng DFT để tính toán các đặc tính điện tử và từ tính của màng Fe 2D sau khi mô phỏng cấu trúc, nhằm đánh giá ảnh hưởng của khuyết tật đến moment từ. Chủ thể: nhóm nghiên cứu lý thuyết vật liệu, thời gian 9 tháng.

4. Thí nghiệm xác nhận mô hình mô phỏng:
   Hợp tác với các phòng thí nghiệm thực nghiệm để chế tạo màng Fe 2D trên nền graphene hoặc các vật liệu 2D khác, sử dụng kính hiển vi điện tử truyền qua (TEM) và kỹ thuật quang phổ để so sánh cấu trúc và tính chất với mô hình mô phỏng. Chủ thể: viện nghiên cứu vật liệu, thời gian 1-2 năm.

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

1. Nhà nghiên cứu vật liệu 2D và nano:
   Luận văn cung cấp dữ liệu chi tiết về quá trình chuyển pha và cấu trúc màng Fe 2D, hỗ trợ phát triển vật liệu mới với tính chất từ và cơ học đặc biệt.

2. Chuyên gia mô phỏng động lực học phân tử:
   Tham khảo phương pháp mô phỏng, lựa chọn thế tương tác EAM và kỹ thuật phân tích dữ liệu trong mô hình kim loại 2D.

3. Kỹ sư phát triển thiết bị lưu trữ và vi điện tử:
   Hiểu rõ cơ chế tạo màng Fe 2D giúp thiết kế vật liệu từ có mật độ lưu trữ cao, kích thước nhỏ, ứng dụng trong công nghệ nano.

4. Sinh viên và học viên cao học ngành vật lý kỹ thuật, khoa học vật liệu:
   Tài liệu tham khảo quý giá về mô phỏng quá trình chuyển pha, phân tích cấu trúc và tính chất nhiệt động học của vật liệu 2D kim loại.

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

1. Phương pháp mô phỏng động lực học phân tử có ưu điểm gì trong nghiên cứu vật liệu 2D?
   MD cho phép mô phỏng chi tiết chuyển động nguyên tử theo thời gian thực, giúp quan sát cơ chế chuyển pha và thay đổi cấu trúc ở cấp độ nguyên tử, phù hợp với vật liệu 2D có tính chất phức tạp.

2. Tại sao chọn thế tương tác EAM cho mô phỏng màng Fe 2D?
   EAM phản ánh chính xác tương tác nhiều nguyên tử trong kim loại chuyển tiếp, bao gồm hiệu ứng bề mặt và khuyết tật, vượt trội hơn các thế cặp đơn giản như Morse hay Lennard-Jones.

3. Nhiệt độ chuyển pha của màng Fe 2D là bao nhiêu và ý nghĩa của nó?
   Khoảng 2640 K, đây là nhiệt độ mà màng Fe chuyển từ trạng thái lỏng sang rắn, xác định tính ổn định cấu trúc và tính chất vật lý của vật liệu 2D.

4. Khuyết tật cấu trúc ảnh hưởng thế nào đến tính chất của màng Fe 2D?
   Khuyết tật như vacancy làm biến dạng mạng tinh thể, ảnh hưởng đến moment từ, tính dẫn điện và cơ học, cần được kiểm soát để ứng dụng hiệu quả.

5. Làm thế nào để mô phỏng sát thực tế hơn quá trình tạo màng Fe 2D?
   Giảm tốc độ làm lạnh, tăng kích thước mô hình, áp dụng điều kiện biên đa dạng và kết hợp với tính toán ab initio, đồng thời phối hợp với thí nghiệm thực tế để hiệu chỉnh mô hình.

Kết luận

• Nhiệt độ chuyển pha của màng Fe phẳng hai chiều được xác định khoảng 2640 K, thuộc chuyển pha loại I không có pha trung gian.
• Quá trình làm lạnh từ trạng thái lỏng đến rắn tạo thành cấu trúc lưới tam giác với hằng số mạng 2,455 Å, gần với giá trị lý thuyết.
• Mô hình chứa khoảng 5% khuyết tật nút khuyết, ảnh hưởng đến cấu trúc và tính chất vật liệu.
• Phương pháp mô phỏng động lực học phân tử kết hợp thế EAM là công cụ hiệu quả để nghiên cứu quá trình chuyển pha và cấu trúc vật liệu 2D kim loại.
• Hướng phát triển tiếp theo bao gồm tối ưu mô hình, mở rộng kích thước, kết hợp tính toán ab initio và thí nghiệm xác nhận.

Đề nghị các nhà nghiên cứu và kỹ sư trong lĩnh vực vật liệu 2D tiếp tục khai thác kết quả này để phát triển vật liệu từ mới, nâng cao hiệu suất thiết bị nano và vi điện tử.

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

Chương 1 “Tổng quan tài liệu”, trình bày sự ra đời và vai trò của vật liệu hai chiều, bằng chứng thực nghiệm và các kết quả mô phỏng sự tồn tại của màng sắt một lớp nguyên tử. Chương 2 “Phương pháp nghiên cứu”, giới thiệu về phương pháp MD, thế tương tác và các kỹ thuật tính toán của mô hình. Chương 3 “Kết quả và thảo luận”, trình bày các kết quả về quá trình chuyển pha, các tính chất liên quan đến cấu trúc trong quá trình làm lạnh và cơ chế nguyên tử trong quá trình đông đặc. Chương 4 “Kết luận và hướng phát triển của đề tài”, chương này trình bày về các kết luận thu được và hướng phát triển của đề tài.

Thông qua việc tính toán và phân tích các đại lượng, đề tài đã thu được một số kết quả nổi bật. Nhiệt độ nóng chảy của ironene tinh thể trong không gian 2D được xác định vào khoảng 2640 K. Sau quá trình làm lạnh, đề tài thu được màng sắt phẳng hai chiều có cấu trúc lưới tam giác với hằng số mạng là 2,455 Å, chứa khuyết tật mất nguyên tử (vacancy). Dựa trên các kết quả thu được từ đề tài, học viên dưới sự hướng dẫn của GS.

Võ Văn Hoàng, đã tiến hành phân tích chi tiết hơn và đã có hai báo cáo tham gia hội nghị: - V. Hoàng, “Mô phỏng quá trình hình thành màng sắt đơn lớp 2D từ trạng thái lỏng”, đã tham gia Hội nghị Khoa học và Công nghệ lần thứ 14 tại Đại học Bách khoa thành phố Hồ Chí Minh (30/10/2015). Hằng, “Ironene – Vật liệu 2D mới”, đã tham gia Hội nghị quốc tế lần thứ 5 về Công nghệ nano và ứng dụng tại thành phố Vũng Tàu (11-14/11/2015). TỔNG QUAN TÀI LIỆU 1.

Tổng quan về vật liệu hai chiều Vật liệu khối ba chiều (3D) đã có nhiều ứng dụng thực tế. Vấn đề đặt ra ở đây là có hay không sự tồn tại của vật liệu hai chiều (2D) và chúng có ứng dụng như thế nào trong thực tế? Điều này đã dẫn đến các nghiên cứu lý thuyết và thực nghiệm về vật liệu 2D vô cùng thú vị từ những năm 30 của thế kỷ 20. Hai nhà vật lý nổi tiếng Peierls và Landau với những nghiên cứu của họ [1,2] đã cho rằng việc tồn tại vật liệu 2D mâu thuẫn với lý thuyết bền vững về mặt nhiệt động lực học. Các nguyên tử trong mạng 2D khi dao động nhiệt sẽ gây ra sự tan chảy mạng tinh thể khiến cho sự tồn tại của vật liệu 2D là bất khả thi và việc tìm kiếm là vô vọng.

Kết quả thực nghiệm [3] cũng chứng minh được độ bền vững nhiệt động học của màng mỏng giảm dần theo sự giảm bề dày của màng vật liệu 2D. Khi màng đạt đến độ dày tương ứng với vài chục lớp nguyên tử thì vật liệu 2D sẽ suy thoái và co cụm lại thành vật liệu 3D. Như vậy, bằng lý thuyết và thực nghiệm lúc bấy giờ đều đưa đến kết luận thống nhất là vật liệu 2D không thể tồn tại. Năm 2004, sự thành công của nhóm nhà khoa học Novoselov và các cộng sự [4] trong việc chế tạo vật liệu 2D graphene đã thật sự làm thay đổi những giả thuyết trong vật lý lý thuyết trước đây.

Tuy nhiên, việc khám phá ra graphene của Novoselov và các cộng sự không bác bỏ những nghiên cứu của Peierls và Landau. Lý thuyết không chấp nhận màng 2D tồn tại trên một mặt phẳng tuyệt đối, tức là chiều cao của mặt phẳng bằng không, nhưng điều đó không đồng nghĩa với việc màng 2D tựa trên nền là vật liệu khối 3D hay tự thân hiện hữu trong không gian 3D là không thể tồn tại. Điều này hoàn toàn phù hợp với quan sát của Novoselov và các cộng sự [5]. Quá trình nghiên cứu graphene không những được chú ý đến trong khoảng thời gian gần đây mà còn được nghiên cứu từ thế kỷ trước.

Lịch sử nghiên cứu vật liệu đã ghi nhận không ít nỗ lực để tìm ra vật liệu 2D đầu tiên với bề dày một lớp nguyên tử. Năm 1859, khi cho graphite vào axit mạnh, nhà hóa học Benjamin Brodie [6] quan sát được một dạng vật liệu mà ông gọi là “graphon”. Kỹ thuật ngày nay đã 4 khẳng định “graphon” chính là graphene oxit được cấu tạo từ các tinh thể nhỏ graphene có mật độ cao được bao phủ bởi nhóm hydroxyl và epoxide [7]. Trên lĩnh vực lý thuyết, Phil Wallace dựa trên nghiên cứu về cấu trúc dải của graphene [8] đã đề cập đến sự tồn tại của vật liệu 2D này vào năm 1947.

Lý thuyết này dường như bị lãng quên đến khi graphene thực sự ra đời bằng thực nghiệm thì nó nhanh chóng đóng vai trò soi đường cho các nghiên cứu sau này. Vogt sử dụng kính hiển vi điện tử truyền qua (TEM) quan sát mẫu sau khi làm khô một giọt graphene oxide, kết quả cho thấy một số mảnh có độ dày vài nanomet [9]. Năm 1962, Ulrich Hofmann và Hanns-Peter Boehm đã cố gắng tìm kiếm những mảnh vỡ mỏng nhất có thể có của các mảnh graphene oxide và xác định một số trong số chúng gần như có cấu trúc đơn lớp [10]. Sự thành công của Novoselov và các cộng sự trong việc chế tạo graphene có ý nghĩa vô cùng quan trọng cả lý thuyết và ứng dụng.

Về mặt lý thuyết, graphene đã làm thay đổi những quan điểm trước đó về việc có hay không tồn tại vật liệu 2D. Về mặt ứng dụng, graphene được xem là một ứng cử viên vô cùng sáng giá cho các ứng dụng trong vi điện tử, pin năng lượng mặt trời, các thiết bị quang điện tử,… Sau sự thành công đầy hứa hẹn của graphene, các nghiên cứu thực nghiệm về vật liệu 2D lại thu hút sự quan tâm to lớn như hBN, MoS2, MoSe2 [11,12], silicene [13,14]. Tương tự như graphene, silicene đã được tiên đoán sự tồn tại của nó trên lý thuyết từ khá lâu [13] nhưng những bằng chứng thực nghiệm của vật liệu này mới xuất hiện gần đây [15-23]. Tuy xuất hiện trễ hơn graphene nhưng silicene có khả năng là vật liệu bán dẫn nên cũng hứa hẹn có những đóng góp to lớn trong ngành công nghiệp bán dẫn và điện tử [15-18].

Tổng quan nghiên cứu thực nghiệm và mô phỏng của màng sắt hai chiều 1. Nghiên cứu thực nghiệm của màng sắt hai chiều Các vật liệu 2D xuất hiện ngày càng nhiều trong thực tế đã khuyến khích các nhà khoa học đặt ra một vấn đề nữa là liệu rằng các vật liệu 2D có tính chất từ, mà đại diện là sắt có thể tồn tại hay không. Vật liệu từ 2D nếu tồn tại sẽ đóng vai trò quan 5 trọng trong các thiết bị lưu trữ đặc trưng bởi mật độ lưu trữ cao và kích thước được thu nhỏ tối đa [24]. Sự quan tâm của cộng đồng khoa học với vật liệu từ 2D trên lý thuyết và thực nghiệm thật ra không mới và đã bắt đầu từ sự khám phá ra khả năng tăng cường moment từ của vật liệu bằng cách giảm số phối vị của nguyên tử [25,26].

Nghiên cứu đã chỉ ra rằng, cụm sắt từ một lớp nguyên tử có moment từ trên mỗi nguyên tử lớn hơn nhiều lần so với vật liệu từ 3D [24,27]. Tuy vậy, liên kết kim loại và cấu trúc 3D khiến kim loại không thể tồn tại ở dạng một lớp nguyên tử. Trước đây đã có một số báo cáo của màng mỏng sắt như spin của sắt trên nền là tinh thể Cu (100) [28] hay tính chất cơ học trên tinh thể MgO (001) [29]. Tuy nhiên, các kết quả trên thu được đối với màng sắt đa lớp gồm nhiều đơn lớp sắt.

Ngoài ra, các màng sắt đa lớp này tương tác với nền, ảnh hưởng đến các đại lượng như moment từ của màng sắt. Các nghiên cứu lý thuyết và thực nghiệm hiện nay tập trung vào việc hình thành màng sắt một lớp nguyên tử trên vật liệu 2D khác ví dụ như graphene [30-32] hay cụm nguyên tử sắt trên graphene [33]. Trong đó, phải kể đến kết quả của Jiong Zhao và các cộng sự [34] đã lần đầu tiên tạo thành công màng sắt đơn lớp và đồng thời gần như mở ra phương pháp mới để chế tạo màng sắt đơn lớp hai chiều. Trong thí nghiệm của Jiong Zhao và các cộng sự, màng sắt đơn lớp được tạo ra như sau: đầu tiên, màng graphene được tạo ra trước bằng phương pháp CVD (chemical vapor deposition) trên nền Ni/Mo [35].

Màng graphene đơn lớp được đưa lên lưới carbon sau khi sử dụng dung môi FeCl3 để tách màng graphene khỏi nền Ni/Mo [36]. Sau quá trình chuyển, các nguyên tử sắt còn sót lại từ dung môi FeCl3 tạo thành tinh thể sắt có kích thước nano trên bề mặt của graphene như BCC (Body Centered Cubic), HCP (hexagonal close-packed) và màng sắt đơn lớp dưới chiếu xạ điện tử của kính hiển vi điện tử truyền qua quang cầu sai được điều chỉnh ở điện thế thấp (LVACTEM).1 biểu diễn màng graphene và các cấu trúc tinh thể sắt hình thành trên màng graphene. Hình chụp màng sắt đơn lớp hai chiều bằng kính hiển vi điện tử truyền qua cường độ lớn. (A) cho thấy cấu trúc của màng sắt đơn lớp hai chiều bằng thực nghiệm.

(B, D) hình chụp sau khi xử lý bằng FFT để lọc nhiễu. (C, E, F) biểu diễn màng sắt và hằng số mạng theo mô phỏng và thực nghiệm. cho ta thấy rõ hơn màng sắt đơn lớp hai chiều bằng kính hiển vi điện tử truyền qua cường độ lớn. Kết quả cho thấy, màng sắt đơn lớp hai chiều có dạng lưới ô vuông với hằng số mạng là 2,65 ± 0,05 Å theo thực nghiệm.

Jiong Zhao và các cộng sự cũng đã tiến hành tính toán DFT (density functional theory) dựa trên xấp xỉ gradient tổng quát (GGA) với hàm Perdew-Burke-Ernzerhof [37]. Tính toán chỉ ra rằng, hằng số mạng bền vững nhất cho lưới ô vuông của màng sắt đơn lớp là xấp xỉ 2,35 Å. Tuy nhiên, tính toán năng lượng cho hai kết quả này lệch nhau không lớn (0,2 eV trên mỗi nguyên tử) nên có thể chấp nhận được. Sở dĩ có sự sai khác này do trong tính toán mô phỏng ta đã bỏ qua một vài yếu tố ảnh hưởng như: spin quỹ đạo hay các tính chất cơ học giữa graphene và màng sắt đơn lớp hai chiều.

Màng sắt dạng lưới tam giác (a-d) [38] 8 Bằng khảo sát thực nghiệm cũng như sử dụng bức xạ điện tử, Wei Yang và Hongtao Wang phát hiện ra màng sắt đơn lớp còn tồn tại dưới dạng lưới tam giác [38] bền vững hơn lưới ô vuông trong thí nghiệm của Jiong Zhao và các cộng sự. cho thấy màng sắt dưới dạng lưới tam giác trên nền graphene. biểu diễn các dạng liên kết giữa màng sắt và biên graphene thu được bằng mô phỏng và thực nghiệm.