Hành vi biến dạng và cơ chế hư hỏng của 3C-SiC trong quá trình thụt nano +
PAM-XIAMEN cung cấp tấm bán dẫn 3C-SiC cho mục đích nghiên cứu. Để biết chi tiết, vui lòng liên hệ:[email protected]
Cacbua silic khối (3C-SiC) thể hiện tiềm năng lớn trong các lĩnh vực tiên tiến như hệ thống vi cơ điện tử, điện toán lượng tử và thiết bị quang điện tử, nhờ dải rộng, độ dẫn nhiệt cao, cường độ trường phân hủy tới hạn cao và độ trơ hóa học tuyệt vời. Tuy nhiên, độ cứng cực cao và độ giòn nội tại của 3C-SiC khiến việc gia công siêu chính xác của nó trở nên vô cùng khó khăn và hư hỏng bề mặt/dưới bề mặt do gia công gây ra đã trở thành nút thắt chính hạn chế các ứng dụng cao cấp của nó. Do đó, sự hiểu biết sâu sắc về hành vi biến dạng và cơ chế hư hỏng của 3C-SiC khi tải tiếp xúc có ý nghĩa lớn về mặt lý thuyết và kỹ thuật để tối ưu hóa các kỹ thuật xử lý của nó.
Trong những năm gần đây, các nhà nghiên cứu đã khám phá rộng rãi phản ứng cơ học của 3C-SiC bằng cách kết hợp các thí nghiệm thụt nano với mô phỏng động lực phân tử. Zhu và cộng sự. đã nghiên cứu một cách có hệ thống hành vi thụt nano của 3C-SiC trên các mặt phẳng tinh thể (010), (110) và (111) bằng cách sử dụng mô phỏng động lực phân tử (MD), cho thấy các đặc tính dị hướng của biến dạng đàn hồi và tiến hóa trật khớp, đồng thời giới thiệu khối tứ diện Thompson để giải thích cơ chế nội tại của biến dạng dẻo. Trong khi đó, nhiệt độ, yếu tố chính trong môi trường sử dụng thực tế, ảnh hưởng đáng kể đến tính chất cơ học của 3C-SiC. Kourbani và cộng sự. được phát hiện thông qua mô phỏng MD rằng việc tăng nhiệt độ sẽ làm giảm độ cứng và mô đun Young của vật liệu và thúc đẩy quá trình vô định hình. Mặt khác, các phương pháp cơ học liên tục truyền thống có những hạn chế trong việc xử lý các vấn đề không liên tục như sự hình thành và lan truyền vết nứt. Để giải quyết vấn đề này, Yang và Zhu đã sử dụng động lực học để mô phỏng quá trình phát triển hư hỏng của 3C-SiC dưới tác dụng thụt đầu dòng và trầy xước, tái tạo thành công hiện tượng bật vào và hình thành các vết nứt xuyên tâm. Những nghiên cứu này cung cấp tài liệu tham khảo đa quy mô để hiểu được đặc tính cơ học và lý thuyết gia công của 3C-SiC.
1. Biến dạng đàn hồi dị hướng và tính chất dẻo ban đầu của 3C-SiC
3C-SiC có cấu trúc hỗn hợp kẽm và phản ứng cơ học của nó phụ thuộc rất nhiều vào hướng tinh thể. Thông qua mô phỏng MD, Zhu et al. nhận thấy rằng trong quá trình thụt nano trên các mặt phẳng (010), (110) và (111), các vùng biến dạng đàn hồi thể hiện sự đối xứng khác nhau rõ rệt: mặt phẳng (010) hiển thị tính đối xứng bốn lần, với vùng biến dạng hình chóp mở rộng dọc theo hướng <110>; mặt phẳng (110) thể hiện tính đối xứng trung tâm, với vùng biến dạng chia thành hai phần; mặt phẳng (111) thể hiện sự đối xứng bậc ba, với ba vùng biến dạng đàn hồi cũng mở rộng dọc theo hướng <110>. Kết quả này chỉ ra rằng hướng <110> là đường tối ưu để truyền ứng suất đàn hồi.
Hình 1 Hình thái vùng biến dạng đàn hồi của 3C-SiC ở độ sâu vết lõm 1,1nm: mặt phẳng (a,d) (010); (b,e) (110) mặt phẳng; (c,f) (111) mặt phẳng; mũi tên đỏ chỉ hướng truyền biến dạng đàn hồi
Từ các đường cong tải-độ sâu, chiều dài của giai đoạn đàn hồi thay đổi đáng kể giữa các mặt phẳng tinh thể khác nhau. Đường cong Ph cho mặt phẳng (010) rất khớp với mô hình tiếp xúc Hertz với h <2.0nm, trong khi các giai đoạn đàn hồi của mặt phẳng (110) và (111) ngắn hơn đáng kể, cho thấy năng suất nhựa sớm hơn. Đáng chú ý, các sự kiện bật lên xuất hiện trên các đường cong của mặt phẳng (010) và (110), nhưng không được quan sát thấy đối với mặt phẳng (111). Thông qua phân tích chuyển vị nguyên tử, Zhu et al. xác nhận rằng sự kiện bật lên tương ứng với chuyển động có hướng của một số lượng lớn nguyên tử dọc theo hướng <110>, dẫn đến sự giải phóng đột ngột năng lượng đàn hồi dự trữ. Ví dụ: khi độ sâu vết lõm trên mặt phẳng (110) tăng từ 2,2nm lên 2,3nm thì số lượng nguyên tử có độ dịch chuyển lớn hơn 0,1nm tăng từ 9,558 lên 13,767; sự di chuyển nguyên tử tập thể này là cơ chế vi mô trực tiếp của pop-in.
Hình 2 (ac) Tải các đường cong độ sâu của vết lõm nano trên các bề mặt 3C-SiC (010), (110) và (111); (d) So sánh mô hình tiếp xúc đàn hồi ba đường Ph và mô hình tiếp xúc đàn hồi Hertz
Để mô tả một cách định lượng ảnh hưởng của sự phân bố ứng suất đến quá trình khởi tạo dẻo, Zhu et al. tính toán phân bố ứng suất cắt von Mises. Ở độ sâu thụt 2,0nm, độ dốc ứng suất giảm trong vùng biến dạng đàn hồi; khi độ sâu đạt tới 2,1nm và các điểm trật khớp bắt đầu tạo mầm, ứng suất ở vùng trật khớp tăng lên 8–10GPa, trong khi ứng suất giữa các vùng bị biến dạng và không bị biến dạng giảm xuống 5–8GPa. Điều này chỉ ra rằng việc mở rộng trật khớp sẽ làm tiêu tan ứng suất cao cục bộ một cách hiệu quả, cho phép một phần cấu trúc bị biến dạng trở lại cấu trúc kim cương hoàn hảo.
2. Phân tích sự tạo mầm, nhân giống và phân tích tứ diện Thompson
Biến dạng dẻo của 3C-SiC chủ yếu bị chi phối bởi sự trượt lệch vị trí. Sử dụng thuật toán trích xuất sai vị trí, Zhu et al. xác định rằng trật khớp ban đầu là một vòng lặp trật khớp hoàn hảo với vectơ Burgers b = 1/2[1-10], tạo mầm từ vùng ứng suất cắt cao bên dưới đầu lõm. Trượt lệch vị trí biến đổi trình tự sắp xếp nguyên tử từ ABCABC ban đầu thành ABABC, trong đó các nguyên tử silicon lệch khỏi vị trí ban đầu và hình thành liên kết mới với các nguyên tử carbon lân cận dọc theo hướng <110>. Quá trình này đánh dấu sự chuyển đổi không thể đảo ngược từ tính đàn hồi sang tính dẻo.
Hình 3. Mặt cắt ngang của vùng trật khớp trong quá trình thụt nano: (a) chế độ xem cục bộ; (b) mặt cắt ngang; mũi tên nét đứt màu đỏ biểu thị sự trượt của các nguyên tử silicon từ vị trí liên kết ban đầu sang vị trí mới
Khi độ sâu vết lõm tăng lên, sự lệch vị trí của phôi dần dần phát triển thành các dải băng lỗi xếp chồng lên nhau và các vòng trật khớp hình lăng trụ. Lấy mặt phẳng (010) làm ví dụ: ở độ sâu 2,1nm, hai đứt gãy xếp chồng (SF1 và SF2) nằm trên mặt phẳng (111) giao nhau, tạo thành trật khớp thanh cầu thang 1/2[1-10] (SRD), sau đó kết nối thành một vòng trật khớp dải băng; khi độ sâu tăng lên 2,6nm, sự lệch vị trí “tách ra”, tạo thành các vòng lệch vị trí hình lăng trụ bị cô lập. Tương tự, trên các mặt phẳng (110) và (111), trình tự tiến hóa từ trật khớp phôi → trật khớp ruy băng → trật khớp lăng trụ cũng được quan sát thấy, nhưng tốc độ tiến hóa và hình thái hình học khác nhau do tính dị hướng.
Hình 4. Sự phát triển của trật khớp với độ sâu vết lõm trên mặt phẳng (010): (a) h = 2,1nm; (b) h = 2,2nm; (c) h = 2,5nm; (d) h = 2,6nm
Để giải thích một cách có hệ thống tính bất đẳng hướng này, Zhu et al. giới thiệu mô hình tứ diện Thompson. Đối với vết lõm trên mặt phẳng {100}, xếp chồng các lỗi trên bốn mặt phẳng trượt {111} (ACE, AEG, ACG, CEG) giao nhau, mỗi mặt phẳng tạo ra sự lệch vị trí của thanh cầu thang dọc theo hướng <110>, với phản ứng được biểu thị là δβ = δE + Eβ, v.v. Đối với mặt phẳng {110}, một số mặt phẳng trượt vuông góc với bề mặt vết lõm, làm cho các vòng lệch vị trí hình lăng trụ mở rộng vuông góc với bề mặt; đối với mặt phẳng {111}, ba vị trí lệch dải băng thể hiện sự đối xứng ba lần, cuối cùng tạo thành ba vòng lệch vị trí hình lăng trụ. Tứ diện Thompson cung cấp một khuôn khổ hình học ngắn gọn nhưng mạnh mẽ để hiểu tính dị hướng dẻo của gốm có cấu trúc kim cương lập phương.
3. Ảnh hưởng của nhiệt độ đến tính chất cơ học và biến dạng
Trong các ứng dụng thực tế, thiết bị 3C-SiC có thể hoạt động trong phạm vi nhiệt độ rộng, do đó không thể bỏ qua ảnh hưởng của nhiệt độ. Kourbani và cộng sự. đã nghiên cứu phản ứng thụt nano của 3C-SiC trên mặt phẳng (010) trong phạm vi nhiệt độ từ 50K đến 1000K bằng mô phỏng MD. Họ phát hiện ra rằng khi nhiệt độ tăng lên, tải trọng thụt tối đa giảm từ 2555,4nN (50K) xuống 2197,4nN (1000K). Bằng cách sử dụng phương pháp Oliver–Pharr, mô đun Young được phát hiện giảm tuyến tính từ 328GPa (50K) xuống 293,5GPa (1000K), với giá trị ở nhiệt độ phòng (300K) xấp xỉ 318GPa, gần với giá trị thực nghiệm là 314GPa; độ cứng tương tự giảm từ 45,3GPa xuống 37,9GPa. Hiệu ứng làm mềm này phát sinh từ sự gia tăng khoảng cách giữa các nguyên tử và giảm năng lượng liên kết cộng hóa trị Si-C khi nhiệt độ tăng.
Nhiệt độ cũng ảnh hưởng đáng kể đến hành vi trật khớp. Kourbani và cộng sự. phân tích hình thái trật khớp ở các nhiệt độ khác nhau bằng DXA. Trật khớp hoàn hảo với b = 1/2<110> và trật khớp một phần Shockley với b = 1/6<112> được quan sát thấy ở mọi nhiệt độ, với trật khớp hoàn hảo chiếm ưu thế. Khi nhiệt độ tăng lên, tổng chiều dài trật khớp tăng lên rõ rệt và sự hình thành các vòng trật khớp hình lăng trụ kiểu “lasso” trở nên tích cực hơn. Điều này chỉ ra rằng kích hoạt nhiệt thúc đẩy quá trình tạo mầm và trượt trật khớp, tăng cường khả năng biến dạng dẻo.
Hình 5 Phân tích DXA và IDS: cấu trúc lệch vị trí của 3C-SiC sau khi thụt tới độ sâu tối đa (h = 2,5nm) ở các nhiệt độ khác nhau (50K, 100K, 300K, 500K, 800K, 1000K)
Đặc biệt quan tâm là vai trò của nhiệt độ trong việc điều chỉnh sự chuyển pha và vô định hình. Sử dụng phương pháp Xác định cấu trúc kim cương, Kourbani et al. phát hiện ra rằng các nguyên tử có cấu trúc kim cương lục giác (tức là cấu trúc wurtzite B4) xuất hiện trong vùng thụt lề, như trong Hình 6. Khi nhiệt độ tăng, số lượng nguyên tử “khác” (đại diện cho trạng thái vô định hình) giảm, trong khi số lượng nguyên tử HD (cấu trúc wurtzite) cũng giảm ở mức 1000K. Phân tích hàm phân bố xuyên tâm cho thấy rằng sau khi thụt đầu dòng, các đỉnh trật tự tầm xa biến mất, xác nhận rằng sự vô định hình thực sự xảy ra và trở nên rõ rệt hơn ở nhiệt độ cao hơn. Phân tích số phối trí còn tiết lộ thêm rằng khi nhiệt độ tăng, tỷ lệ nguyên tử có CN = 4 giảm, trong khi số nguyên tử có độ phối trí thấp (CN = 1, 2, 3) tăng lên, cho thấy nhiệt độ thúc đẩy quá trình chuyển đổi của các nguyên tử từ mạng hoàn hảo sang cấu trúc rối loạn. Những phát hiện này cho thấy tác động kép của nhiệt độ đến độ ổn định cấu trúc của 3C-SiC: một mặt, kích hoạt nhiệt thúc đẩy trượt trật khớp; mặt khác, nó làm trầm trọng thêm thiệt hại vô định hình.
Hình 6. Sự phân bố các nguyên tử bị biến dạng trong 3C-SiC sau khi thụt đầu dòng ở các nhiệt độ khác nhau (loại trừ các nguyên tử cấu trúc kim cương lập phương): HD, HD (hàng xóm thứ 1), HD (hàng xóm thứ 2) đại diện cho các nguyên tử liên quan đến cấu trúc wurtzite; “khác” đại diện cho các nguyên tử vô định hình
4. Mô phỏng hư hỏng và gãy xương dựa trên Peridynamics
Mặc dù mô phỏng MD có thể giải quyết biến dạng dẻo ở quy mô nguyên tử, quy mô tính toán của chúng thường bị giới hạn ở phạm vi nanomet (hàng chục nanomet), khiến việc so sánh trực tiếp với các thí nghiệm ở quy mô micromet trở nên khó khăn. Hơn nữa, MD không hiệu quả trong việc xử lý sự lan truyền vết nứt trong quá trình gãy giòn. Để giải quyết những vấn đề này, Yang và Zhu đã mô phỏng quá trình thụt đầu dòng và làm xước của 3C-SiC bằng cách sử dụng phương pháp cận động lực học (PD). Peridynamics là một lý thuyết phi tiêu điểm thay thế các phương trình vi phân từng phần bằng các phương trình tích phân, làm cho nó phù hợp một cách tự nhiên để xử lý các vấn đề không liên tục như sự hình thành và lan truyền vết nứt.
Trong các mô phỏng vết lõm, Yang và Zhu đã quan sát thấy hiện tượng bật lên rõ ràng trên đường cong tải trọng-độ sâu, rất phù hợp với hành vi nổ nhựa được quan sát trong mô phỏng MD. Các đường cong tải ở các tốc độ thụt khác nhau (1–20m/s) phần lớn trùng nhau, nhưng trong quá trình dỡ tải, tốc độ cao hơn dẫn đến giảm tải nhanh hơn. Quan trọng hơn, các mô phỏng PD đã tái tạo thành công hình thái vết nứt do vết lõm gây ra: khi độ sâu vết lõm đạt 15μm với tốc độ 5m/s, có thể quan sát rõ ràng các vết nứt hướng tâm, vết nứt bên và vết nứt ở giữa, phù hợp với hình thái vết nứt vết lõm thực nghiệm. Điều này chứng tỏ rằng PD có thể nắm bắt được các đặc tính đứt gãy của vật liệu giòn ở quy mô liên tục.
Hình 7 Hình thái hư hỏng vết lõm của 3C-SiC ở tốc độ vết lõm 5m/s: (a–c) góc nhìn từ trên xuống; (d–i) mặt cắt ngang
Trong các mô phỏng vết xước, Yang và Zhu đã kiểm tra sự thay đổi của lực ma sát và lực pháp tuyến ở các độ sâu vết xước khác nhau (1–4μm). Ở độ sâu 1μm và 2μm, không xuất hiện vết nứt nào ở hai bên rãnh xước và việc loại bỏ vật liệu chủ yếu là nhựa; khi độ sâu tăng lên 3μm và 4μm, các vết nứt hướng tâm xuất hiện ở các cạnh rãnh và cả lực ma sát và lực pháp tuyến đều biểu hiện sự dao động định kỳ. Năng lượng cắt cụ thể (năng lượng tiêu thụ trên một đơn vị thể tích vật liệu bị loại bỏ) tăng phi tuyến tính khi độ sâu vết xước giảm, đạt tới 22GPa ở 1μm, cho thấy hiệu ứng kích thước đáng kể.
Yang và Zhu đã mô phỏng thêm quá trình cào xước hai vết lõm để mô phỏng sự tương tác của nhiều vật liệu mài mòn trong quá trình mài thực tế. Khi khoảng cách giữa hai vết lõm là 50μm, chiều rộng vết xước xấp xỉ 92μm, gấp hơn hai lần chiều rộng của vết xước một vết lõm và khối lượng loại bỏ vật liệu là lớn nhất; ở khoảng cách 60μm, hai vết xước độc lập và khối lượng loại bỏ chỉ gấp đôi so với một vết lõm đơn lẻ; ở khoảng cách 40μm, các vùng thiệt hại chồng lên nhau, nhưng tổng khối lượng loại bỏ giảm. Hiện tượng này được cho là do hiệu ứng chồng chất của trường dịch chuyển: ở một khoảng cách thích hợp, các vết nứt kết nối giữa hai mũi ấn, dẫn đến việc loại bỏ vật liệu hợp tác; khoảng cách quá nhỏ sẽ dẫn đến sự triệt tiêu lẫn nhau của các trường chuyển vị, trong khi khoảng cách quá lớn sẽ không tạo ra sự tương tác. Phát hiện này cung cấp cơ sở lý thuyết để tối ưu hóa các thông số của quá trình nghiền.
Hình 8 Hình thái hư hỏng của vết xước kép (độ sâu vết lõm 3μm): (a,d) khoảng cách 40μm, vết xước không đều; (b, e) khoảng cách 50μm, chiều rộng vết xước tăng đáng kể, khối lượng loại bỏ vật liệu tối đa; (c,f) khoảng cách 60μm, hai vết xước độc lập
Hình 9 Mặt cắt ngang của vết xước kép: (a) khoảng cách 40μm; (b) khoảng cách 50μm; (c) khoảng cách 60μm
Thông qua sự kết hợp giữa mô phỏng MD và PD, các nhà nghiên cứu đã tiết lộ biến dạng dị hướng, tiến hóa trật khớp, hiệu ứng nhiệt độ và cơ chế hư hỏng/gãy gãy của 3C-SiC dưới quá trình thụt nano. Nghiên cứu phát hiện ra rằng các trạng thái đàn hồi và dẻo phụ thuộc rất nhiều vào hướng tinh thể, với các sai lệch tạo mầm dọc theo hệ thống trượt {111<110> và phát triển thành các vòng lăng trụ; tăng nhiệt độ làm giảm độ cứng và mô đun đồng thời thúc đẩy quá trình vô định hình. Phương pháp PD tái tạo thành công hiện tượng lồi lõm trong vết lõm và sự lan truyền vết nứt khi trầy xước, đồng thời xác định khoảng cách giữa hai vết lõm tối ưu để nâng cao hiệu quả loại bỏ vật liệu. Công việc trong tương lai có thể phát triển hơn nữa các mô hình ghép nối đa quy mô và các hàm tiềm năng được học bằng máy để đạt được các dự đoán định lượng từ cơ chế nguyên tử đến quy trình gia công.
Cho dù bạn cần wafer 3C-SiC để nghiên cứu hay cho các ứng dụng công nghiệp, vui lòng liên hệ với chúng tôi qua email tại[email protected] và [email protected].
Người giới thiệu:
- Zhu, B., Zhao, D., Zhang, Z., Niu, Y., Zhang, Z., Zhao, J., … & Zhao, H. (2024). Nghiên cứu nguyên tử về hành vi biến dạng và hiệu ứng dị hướng của 3C-SiC dưới phương pháp thụt nano. Tạp chí Nghiên cứu và Công nghệ Vật liệu, 28, 2636-2647.
- Kourbani, H., Bogtob, S., Samiri, A., Hassani, A., & Hasnaoui, A. (2025). Thụt nano nhiệt của silicon cacbua (3C) 3C-SiC sử dụng mô phỏng động lực phân tử. Truyền thông thể rắn, 403, 116013.
- Xu, Y., & Zhu, P. (2022). Mô phỏng chu động của hư hỏng khi thụt và trầy xước 3C-SiC. Tạp chí Nghiên cứu Vật liệu, 37(23), 4381-4391.